function grsv=gngauss(m,sgma)
%輸入格式可以為【gsrv1，gsrv2】=gngauss（m，sgma）
%或【gsrv1，gsrv2】=gngauss（sgma）
%或【gsrv1，gsrv2】=gngauss
%函數(shù)生成兩個(gè)統(tǒng)計(jì)獨(dú)立的高斯分布的隨機(jī)數(shù)，以m為均值，sgma為方差
%缺省是m=0，sgma=1
if (nargin==0)
       sgma=1;    m=0;        
elseif nargin==1
       sgma=m;   m=0;
end;
u=rand;                                            %產(chǎn)生一個(gè)（0，1）間均勻分布的隨機(jī)數(shù)u
z=sgma*(sqrt(2*log(1/(1-u))));                     %利用上面的u產(chǎn)生一個(gè)瑞利分布的隨機(jī)數(shù)u
u=rand;                                            %重新產(chǎn)生（0，1）間均勻分布的隨機(jī)數(shù)u
grsv(1)=m+z*cos(2*pi*u);   grsv(2)=m+z*sin(2*pi*u);