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