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