gaussn.m

来自「通信系统仿真与分析中经常用到的高斯白噪声的产生方法」· M 代码 · 共 18 行

function un=gaussn(num1)
a=1;
M=2^36;
b=20061054;
amd=5^13;
delt=0.5;
xn=a;
yn=a;
un=zeros(1,num1);
for n=1:num1
    xn=mod(amd*xn,M);
    en=xn/M;
    yn=mod((amd*yn+b),M);
    vn=yn/M;
    cn=((-2)*log(en))^0.5*cos(2*pi*vn);
    un(n)=delt*cn;
end