Path: news.mathworks.com!newsfeed-00.mathworks.com!nlpi057.nbdc.sbc.com!prodigy.net!news.glorb.com!news2.glorb.com!tr22g12.aset.psu.edu!news.mathforum.org!not-for-mail
From: edin <edo_heineken@hotmail.com>
Newsgroups: comp.soft-sys.matlab
Subject: nakagami-m MGF average probability  bpsk
Date: Thu, 15 Jan 2009 15:01:28 EST
Organization: The Math Forum
Lines: 38
Message-ID: <5473317.1232049736743.JavaMail.jakarta@nitrogen.mathforum.org>
NNTP-Posting-Host: nitrogen.mathforum.org
Mime-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
X-Trace: support1.mathforum.org 1232049736 26224 144.118.30.135 (15 Jan 2009 20:02:16 GMT)
X-Complaints-To: news@news.mathforum.org
NNTP-Posting-Date: Thu, 15 Jan 2009 20:02:16 +0000 (UTC)
Xref: news.mathworks.com comp.soft-sys.matlab:511856


I have problem with this file:

clear;
clc;
M=4;
m=4;
gpsk=(sin(pi/M))^2;
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