Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » Software » comp.soft-sys.matlab

Topic: cognitive radio, energy detection
Replies: 0  

Advanced Search

Back to Topic List Back to Topic List  
amit

Posts: 2
Registered: 11/21/12
cognitive radio, energy detection
Posted: Nov 21, 2012 9:51 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

hi, im new to matlab and not used it before and any help, i will be really thankful, im trying to get a graph of PF vs PD
this is my code, but something is wrong, especially for the outputs of y4, im not really sure anymore.

%monte carlo simulation
function [y1,y2,y3,y4]= ed(SNRL,SNRS,SNRH,n) %ed is the energy detection
g1= SNRL:SNRS:SNRH;
g=10.^(g1/10);
beta=0.01; % is the probability of flase alarm

pf1=zeros(1,length(g));
pd1=zeros(1,length(g));
pf2=zeros(1,length(g));
pd2=zeros(1,length(g));

for i=1 : length(g)
y0 = randn(n,1000000);
y1= randn(n,1000000)*sqrt(g(i))+y0;



Tgam = gaminv(1-beta,n/2,2/n*(1+g(i))); %threshold
pf1(i)= gamcdf(Tgam,n/2,2/n*(1+g(i)))
pd1 (i) = gamcdf(Tgam,n/2,2/n); % prob of detection
pf2(i)= length (find(sum(y1.^2)/n<Tgam))/1000000;
pd2 (i) = length (find(sum(y0.^2)/n<Tgam))/1000000;



end
y1=pf1; y2=pd1; y3=pf2; y4=pd2;



Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.