gngauss.m

来自「对QAM的频谱、散图、眼图以及误码率进行了仿真。使用matlab语言。」· M 代码 · 共 21 行

function n=gngauss(m,sgma)
%输入格式可以为[gsrv1,gsrv2]=gngauss(m,sgma)
%或[gsrv1,gsrv2]=gngauss(m,sgma)
%或[gsrv1,gsrv2]=gngauss
%函数生成两个统计独立的高斯分布的随机数，以m为均值，sgma为方差
%缺省时m=0,sgma=1
if nargin==0
    m=0;
    sgma=1;
elseif  nargin==1
        sgma=m;
        m=0;
 end
    u=rand;%产生一个（0，1）间均匀分布的随机数u
    z=sgma*(sqrt(2*log(1/(1-u))));%利用上面的u产生一个瑞利分布的随机数
    u=rand;%重新产生（0，1）间均匀分布的随机数u
    n(1)=m+z*cos(2*pi*u);
    n(2)=m+z*sin(2*pi*u);