MATLAB作业

在不同功率信噪比下虚警概率P(H1/H0)与检测概率P(H1/H1)的关系曲线

00.10.20.30.40.5

0.60.70.80.91

P(H1/H0)P (H 1/H 1)

Matlab 程序如下：

clear all,close all

clc

d=[0.01,0.5,1,2]; % d^2为功率信噪比，d^2=NA^2/σ^2 eta=0.01:0.1:50; % eta 为似然比检测器的门限值 for j=1:4

for i=1:500

a(i,j)=log(eta(i))/d(j)+d(j)/2;

b(i,j)=log(eta(i))/d(j)-d(j)/2;

y=sym('(1/(2*pi))^(1/2)*exp(-x^2/2)');

Pf(i)= double(int(y,a(i,j),inf)); % Pf 为虚警概率

Pd(i)= double(int(y,b(i,j),inf)); % Pd 为检测概率

end

plot(Pf,Pd)

hold on

end

xlabel('P(H1/H0)')

ylabel('P(H1/H1)')

grid on

在不同门限值下的检测概率P(H1/H1)与功率信噪比的平方根d 的关系曲线

0246

810121400.1

0.2

0.3

0.4

0.5

0.60.7

0.8

0.9

1

d/d B P (H 1/H 1)

Matlab 程序如下：

clear all,close all

clc

d=1:0.1:4; % d^2为功率信噪比，d^2=NA^2/σ^2 d1=10*log(d);

eta=[1,2,5,10,20]; % eta 为似然比检测器的门限值

for j=1:5

a=log(eta(j))./d-d./2;

y=sym('(1/(2*pi))^(1/2)*exp(-x^2/2)');

for i=1:31

Pd(i)= double(int(y,a(i),inf)); %Pd 为检测概率

end

plot(d1,Pd)

hold on

end

grid on

xlabel('d/d B') ylabel('P(H1/H1)')