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;