主题:关于Fortran中生成正态分布随机数的问题
要产生i个正态分布的随机数
call random_seed()
do i=1,1000
call random_number(t)
a(i)=t
enddo
这样产生的随机数经过相关系数检验,各自的相关性还比较满意,相关系数较小
但是正态分布并不太好,chi平方检验出结果不理想,而且可以算一下,曲线的面积并不为1,是0.91左右;同样的积分方法计算其他图形,误差都在0.1%以下,所以计算不会错。
所以我想找一下Fortran生成正态分布随机数的方法。
call random_seed()
do i=1,1000
call random_number(t)
a(i)=t
enddo
这样产生的随机数经过相关系数检验,各自的相关性还比较满意,相关系数较小
但是正态分布并不太好,chi平方检验出结果不理想,而且可以算一下,曲线的面积并不为1,是0.91左右;同样的积分方法计算其他图形,误差都在0.1%以下,所以计算不会错。
所以我想找一下Fortran生成正态分布随机数的方法。