gngauss.m

来自「短波信道抗多音干扰的性能分析及其仿真」· M 代码 · 共 16 行

function[gsrv1,gsrv2]=gngauss(m,sgma)
% GNGAUSS generates two independent Gaussian random variables
% with mean m and deviation sgma

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))));
u=rand;
gsrv1=m+z*cos(2*pi*u);
gsrv2=m+z*sin(2*pi*u);