当量正态化的matlab程序

案例：已知非线性状态方程 567fr-0.5H2 =0, f 服从正态分布， =0.6 ，变异系数 =0.131 ， r 服从正态分布， =2.18 ， =0.03 ； H 服从对数正态分布， =32