gngauss.m

来自「数字通信系统仿真」· M 代码 · 共 15 行

function grsv=gngauss(m,sgma,length)
%本程序的功能是产生高斯白噪声
%m    表示均值
%sgma   表示方差 
%length   序列长度
%grsv     返回高斯噪声序列列向量
%==========================================================
for i=1:1:length
    u=rand;
    z=sgma*(sqrt(2*log(1/(1-u))));
    u=rand;
    grsv(1,i)=m+z*cos(2*pi*u);
    %grsv(1,i)=m+z*sin(2*pi*u);
end
%=================结束=====================================