gngauss.m

来自「两个最近通信常用的数字调制技术的仿真,mpsk」· M 代码 · 共 14 行

function[gsrv1,gsrv2]=gngauss(m,sgma)
% m--mean,  sgma--standard deviation 
if nargin==0,
    m=0;sgma=1;
elseif nargin==1,
    sgma=m;m=0;
end;
u=rand;
z=sgma*(sqrt(2*log(1/(1-u))));     %a Rayleigh distributed random variable
u=rand;
gsrv1=m+z*cos(2*pi*u);
gsrv2=m+z*sin(2*pi*u);