gngauss.m

来自「该MATLAB程序评估了viterbi算法在衰落信道下的卷积译码的性能」· M 代码 · 共 11 行

function[g1,g2]=gngauss(m,sgma)
if nargin == 0
     m=0;sgma=1;
elseif nargin == 1
     m=0;sgma=1;
end;
u=rand;
z=sgma*(sqrt(2*log(1/(1-u))));
u=rand;
g1=m+z*cos(2*pi*u);
g1=m+z*sin(2*pi*u);