gngauss.m

来自「viterbi译码以及psk调制解调的源程序」· M 代码 · 共 19 行

function [ gsrv1,gsrv2 ] = gngauss( m,sgma )
%   [ gsrv1,gsrv2 ] = gngauss( m,sgma )
%   [ gsrv1,gsrv2 ] = gngauss( sgma )
%   [ gsrv1,gsrv2 ] = gngauss
%                       GNGAUSS函数生成两个相互独立的高斯分布的变量，m为均值，sgma为标准方差。
%                       如果只有一个参数的输入，此参数被认定为标准方差，此时均值被默认为"0"。
%                       如果没有参数输入，此时生成两个标准正态分布的变量，均值为"0"，标准方差为"1"。
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);