gngauss.m

完成一个利用矩形信号星座图的M=16的QAM通信系统的MONTECARLO仿真

 function [grsv1 grsv2]=gngauss(m,sgma) 
 %输入格式形式可以为[gsrv1,gsrv2]=gngauss(m,sgma) 
 %或[gsrv1,gsrv2]=gngauss(sgma) 
 %或[gsrv1,gsrv2]=gngauss 
 %函数生成两个统计独立的高斯分布的随机数,以m为均值,sgma为方差 
 %确省时m=0,sgma=1 
 if nargin==0 
     m=0;                 sgma=1; 
 elseif  nargin==1 
     sgma=m;               m=0; 
 end 
 u=rand;                                %产生一个(0,1)间均匀分布的随机数u 
 z=sgma*(sqrt(2*log(1/(1-u))));         %利用上面的产生一个瑞利分布随机数 
 u=rand;                                %重新产生(0,1)间均匀分布的随机数u 
 grsv1=m+z*cos(2*pi*u); 
grsv2=m+z*sin(2*pi*u);