Continue to Site

Welcome to EDAboard.com

Welcome to our site! EDAboard.com is an international Electronics Discussion Forum focused on EDA software, circuits, schematics, books, theory, papers, asic, pld, 8051, DSP, Network, RF, Analog Design, PCB, Service Manuals... and a whole lot more! To participate you need to register. Registration is free. Click here to register now.

Matlab nakagami Fading Channel Simulation?

Status
Not open for further replies.

gr8_halem

Newbie level 4
Joined
Jan 31, 2007
Messages
6
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,321
nakagami matlab

Matlab nakagami Fading Channel Simulation? help me......:cry:
 

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.
 

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
 

Status
Not open for further replies.

Similar threads

Part and Inventory Search

Welcome to EDABoard.com

Sponsor

Back
Top