gngauss.txt

来自「关于瑞利信道」· 文本 代码 · 共 16 行

function [gsrv]=gngauss(m,sgma)
%  [gsrv1,gsrv2]=gngauss(m,sgma)
%  [gsrv1,gsrv2]=gngauss(sgma)
%  [gsrv1,gsrv2]=gngauss
if nargin==0
   m=0;sgma=sqrt(0.5);
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);
gsrv=gsrv1+j*gsrv2;