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)的另一個(gè)均勻隨機(jī)變量
  gsrv1=m+z.*cos(2*pi*u);
  gsrv2=m+z.*sin(2*pi*u);