function[gsrv1,gsrv2]=gngauss(m,sgma)
%             [gsrv1,gsrv2]=gngauss(m,sgma)
%             [gsrv1,gsrv2]=gngauss(sgma)
%             [gsrv1,gsrv2]=gngauss
%             GNGAUSS產(chǎn)生兩個(gè)均值為m、標(biāo)準(zhǔn)差sgma的、互相獨(dú)立的
%             高斯隨機(jī)變量。如果缺少其中一個(gè)輸入?yún)?shù)，則取均值
%             為0。如果既沒有給出均值也沒有給出方差，本函數(shù)將
%             產(chǎn)生兩個(gè)標(biāo)準(zhǔn)高斯隨機(jī)變量。
if nargin==0,
	m=0;sgma=1;
elseif nargin==1,
	sgma=m;m=0;
end;

u=rand;                %在區(qū)間(0,1)內(nèi)的均勻隨機(jī)變量
z=sgma*(sqrt(2*log(1/(1-u))));   %瑞利分布隨機(jī)變量
u=rand;                %在區(qū)間(0,1)內(nèi)的均勻隨機(jī)變量
gsrv1=m+z*cos(2*pi*u);
gsrv2=m+z*sin(2*pi*u);