gngauss.m

来自「QAM的蒙特卡洛仿真QAM的蒙特卡洛仿真QAM的蒙特卡洛仿真.」· M 代码 · 共 20 行

function[gsrv1,gsrv2]=gngauss(m,sgma) 
              [gsrv1,gsrv2]=gngauss(m,sgma) 
              [gsrv1,gsrv2]=gngauss(sgma) 
              [gsrv1,gsrv2]=gngauss 
%             GNGAUSS产生两个均值为m、标准差sgma的、互相独立的 
%             高斯随机变量。如果缺少其中一个输入参数，则取均值 
%             为0。如果既没有给出均值也没有给出方差，本函数将 
%             产生两个标准高斯随机变量。 
if nargin==0, 
	m=0;sgma=1; 
elseif nargin==1, 
	sgma=m;m=0; 
end; 
 
u=rand;                %在区间(0,1)内的均匀随机变量 
z=sgma*(sqrt(2*log(1/(1-u))));   %瑞利分布随机变量 
u=rand;                %在区间(0,1)内的均匀随机变量 
gsrv1=m+z*cos(2*pi*u); 
gsrv2=m+z*sin(2*pi*u);