How can i generate Nakagami fading channel ?

I am doing some wireless communications simulations using Matlab. There are two famous channel models Rayleigh fading and Nakagami fading. I can easily generate Rayleigh fading h=(randn(1,1)+1i*randn(1,1))/sqrt(2).i want to genrate Nakagami fading. If anyone has experience please share the code with me. Thanks

[Mohammad Ali]

For nakagami channel you have to write two line instead of one. first you have to define the pdf as

pd = makedist('Nakagami','mu',m,'omega',omega);

Then to generate a random smaples from above pdf write as

h  = random(pd,N,Itr);

where (N,Itr) is size of channel.

Answer 1

[Angga Dwiki Wicaksono]

% The channel coefficient with PDF Nakagami-m distributed random variable

pd1 = makedist('nakagami',m_h,o);

h1 = (random(pd1,N,Nt) + j*random(pd1,N,Nt)).*sqrt(0.5); % Tx-IRS % There is no pathloss

h2 = (random(pd1,Nr,N) + j*random(pd1,Nr,N)).*sqrt(0.5); % IRS-Rx % There is no pathloss

%The SI channel between the transmit and receive antennas of Use

pd2 = makedist('nakagami',m_I,oI);

hI = (random(pd2,Nt,Nr) + j*random(pd2,Nt,Nr)).*sqrt(0.5) ;% Rx-Tx