高斯噪聲的產生
 function[gsrv1,gsrv2]=gngauss(m,sgma)
%[gsrv1,gsrv2]=gngauss(m,sgma)
%[gsrv1,gsrv2]=gngauss(sgma)
%[gsrv1,gsrv2]=gngauss
%gngauss.m程序產生兩個均值為m,標準偏差為sgma的相互獨立的高斯隨機變量.
%如果缺少其中一個輸入參數,則取均值為0,如果既沒有均值也沒有給出方差,本
%函數將產生兩個標準高斯隨機變量.

  if nargin==0;
    m=0;sgma=1;
  else if 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);

  end