Re: nakagami fading
Hi I tried to get the corresponding M files but was denied access. Can someone please give me the m files to model lognormal and gamma gamma fading please?? I will be grateful.
i need a matlab code for nakagami fading model. I am considering fading model as the sum of rayleigh and WGN fading, considering recvd signal Rnak(t)=a1exp(.5-m)g(t)+a2(1-exp.5-m)(1-exp1-m)C+Rray(t) where a1=1/(sqrtexp-1) and a2=sqrtexp/(sqrtexp-1), C is correction factor c=10.4exp(-19.06m)+1.477exp(.07531m) the value of m=1,.5,2 etc:sad:
---------- Post added at 10:49 ---------- Previous post was at 10:48 ----------
anyone working in nakagami fading model simulation, plz frwd the matalb code.
---------- Post added at 11:02 ---------- Previous post was at 10:49 ----------
clear;
clc;
M=4;
m=4;
gpsk=(sin(pi/M))^2; %mpsk SEP
gamasr_dB=[1:0.1:20];
for j=1:length(gamasr_dB)
gamasr(j)=10.^(gamasr_dB(j)/10);
alfa(j)=sqrt((gpsk.*gamasr(j)/m)/(1+(gpsk.*gamasr(j)/m)))*cot(pi/M);
K(j)=(1/pi)*sqrt((gpsk.*gamasr(j)/m)/(1+gpsk.*gamasr(j)/m));
suma1=0;
for k=0
m-1)
suma1=suma1 +(factorial(2*k)/(factorial(k)^2))*1/(4*(1+(gpsk.*gamasr(j)/m))^k);
end;
suma2=0;
for k=1
m-1)
for i=1:k
T(i)=(factorial(2*k)/(factorial(k)^2))/(factorial(2*(k-i))/(((factorial(k-i))^2)*(4^i)*(2*(k-i)+1)));
BR(i)=T(i)*(cos(atan(alfa(j))))^(2*(k-i)+1);
end;
suma2=suma2+BR(i)/(1+gpsk.*gamasr(j)/m)^k;
end;
Ps(j)=((M-1)/M)-K(j)*((pi/2+atan(alfa(j)))*suma1+sin(atan(alfa(j))*suma2));
Pb(j)=Ps(j)/log2(M);
end;
semilogy(gamasr_dB,Pb);
grid on
hold on