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);