gaussian.m

来自「雷达系统建模与仿真的matlab程序」· M 代码 · 共 23 行

function [y1,y2]=gaussian(n)
%产生数据量为n的两个相互独立高斯分布y1、y2
%---------------------------------------
%
k=1;
y1=zeros(1,n);
y2=zeros(1,n);
while(k<=n)
    u1=junyun(1);
    u2=junyun(1);
    v1=2*u1-1;
    v2=2*u2-1;
    s=v1^2+v2^2;
    if(s>=1)
        continue;
    elseif(s==0)
        k=k+1;
    else
        y1(k)=v1*sqrt(-2*log(s)/s);
        y2(k)=v2*sqrt(-2*log(s)/s);
        k=k+1;
    end
end