Path: news.mathworks.com!not-for-mail
From: "Sadik " <sadik.hava@gmail.com>
Newsgroups: comp.soft-sys.matlab
Subject: Re: nakagami-m MGF average probability  bpsk
Date: Thu, 15 Jan 2009 20:42:02 +0000 (UTC)
Organization: The MathWorks, Inc.
Lines: 45
Message-ID: <gko72q$l0a$1@fred.mathworks.com>
References: <5473317.1232049736743.JavaMail.jakarta@nitrogen.mathforum.org>
Reply-To: "Sadik " <sadik.hava@gmail.com>
NNTP-Posting-Host: webapp-05-blr.mathworks.com
Content-Type: text/plain; charset="ISO-8859-1"
Content-Transfer-Encoding: 8bit
X-Trace: fred.mathworks.com 1232052122 21514 172.30.248.35 (15 Jan 2009 20:42:02 GMT)
X-Complaints-To: news@mathworks.com
NNTP-Posting-Date: Thu, 15 Jan 2009 20:42:02 +0000 (UTC)
X-Newsreader: MATLAB Central Newsreader 1666517
Xref: news.mathworks.com comp.soft-sys.matlab:511860


edin <edo_heineken@hotmail.com> wrote in message <5473317.1232049736743.JavaMail.jakarta@nitrogen.mathforum.org>...
> 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

Hello edin,

What problem do you have? Please indicate the line at which the program gives error and also give a copy of the error message.

Thanks.