gngauss.m

来自「功能如名称主要就是介绍蒙特卡罗方法的使用」· M 代码 · 共 19 行

function[gsrv1,gsrv2]=gngauss(m,sgma)
% [gsrv1,gsrv2]=gngauss(m,sgma)
% [gsrv1,gsrv2]=gngauss(sgma)
% [gsrv1,gsrv2]=gngauss
% gngauss带有（均值m和标准偏差sgma生成两个独立标准高斯随机变量,
% 如果输入参数丢失，则以均值为0，标准偏差为给定的参数，如果没有
% 给定均值和变量，则产生两个标准高斯随机变量。

if nargin==0,
    m=0;sgma=1;
elseif nargin==1,
    sgma=m;m=0;
end;
u=rand;   %(0,1)一致随机变量
z=sgma*(sqrt(2*log(1/(1-u))));  %瑞利分布随机变量
u=rand; %（0，1）另一个一致性随机变量
gsrv1=m+z*cos(2*pi*u);
gsrv2=m+z*sin(2*pi*u);