gngauss.m

来自「基于MATLAB实现信号的调制解调源代码的简单事例.」· M 代码 · 共 11 行

function[gsrv1,gsrv2]=gngauss(m,sgma)
if nargin==0,
    m=0;sgma=1;
elseif nargin==1;
sgma=5;
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);
end;