协同能量检测matlab仿真

单节点与协作能量检测性能的ROC曲线，出两个图，一个是虚警概率和漏检概率的比较，一个是检测概率和漏检概率的比较。

clc;
clear all;
u=20;
r=10;
n1=2;
n2=3;
n3=4;
qf=0.0001:0.01:1;
%qf=1e-19:1e-19:1e-17;
vt=gaminv（（1-qf）u2）;
qd=marcumq（sqrt（2*r）sqrt（vt）u）;
%qd=qf+（（r+1）/r）^（u-1）*exp（-vt/（2*（r+1）））.*gammainc（（vt*r）/（2*（r+1））u）;


Qf1=1-（1-qf）.^n1;
Qd1=1-（1-qd）.^n1;
Qf2=1-（1-qf）.^n2;
Qd2=1-（1-qd）.^n2;
Qf3=1-（1-qf）.^n3;
Qd3=1-（1-qd）.^n3;
qm=1-qd;
Qm1=1-Qd1;
Qm2=1-Qd2;
Qm3=1-Qd3;



figure（1）;
semilogx（qfqd‘--‘Qf1Qd1‘-*‘Qf2Qd2‘-+‘Qf3Qd3‘-o‘）;
axis（[1e-3 1 0 1]）
%axis（[0 1e-17 0 1]）
grid;
xlabel（‘Probability of False Alarm Pf‘）;
ylabel（‘Probability of Detection Pd‘）;
h=legend（‘单节点‘‘n=2‘‘n=3‘‘n=4‘4）;
figure（2）
loglog（qfqm‘-*‘Qf1Qm1‘-*‘Qf2Qm2‘-+‘Qf3Qm3‘-o‘）;
grid;
axis（[1e-3 1 1e-4 1]）;
h=legend（‘单节点‘‘n=2‘‘n=3‘‘n=4‘3）;
xlabel（‘Probability of False Alarm Pf‘）;
ylabel（‘Probability of Miss Pm‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        863  2011-12-29 13:05  xietong.m

----------- ---------  ---------- -----  ----

                  863                    1