gngauss.m

来自「认知无线电系统中多载波信号带外干扰的抑制」· M 代码 · 共 21 行

 高斯噪声的产生
 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