nakagamirnd_gamma.m

QPSK通过瑞利信道的MATLAB仿真

function na=nakagamirnd_gamma(m1,N1,N2)
%这个函数用gamma分布的方法产生服从Nakagami_m的【N1,N2】维随机矩阵
%-----------------------------2005.5.30  ----------------------------------
% Shuai lujun,Information Engineering college,Dalian Maritime Univ.
mm=m1;                               %nakagami_m的m系数
w=1;                                 %nakagami_m的功率，当他等于1时为归一化的
aa=mm;                               %gamma函数的参数a
bb=w/mm;                             %gamma函数的参数b
cc=gamrnd(aa,bb,N1,N2);
z=sqrt(cc);
for i=1:N1
    g(i,:)=z((i-1)*N2+1:i*N2);
end
na=g;