%3.1通过计算机产生的仿真数据,对似然比检测方法的检测概率、虚警概率、漏警概率进行仿真计算
%产生检测信号
clear all;close all;clc;
t1=0:49;
s1=t1/50;
t2=50:149;
s2=-t2/50+2;
t3=150:199;
s3=t3/50-4;
s=[s1,s2,s3]; %产生信号
%4.1通过改变P(H1)和P(H0)来改变判决的门限(风险系数C_10 和C_01不变),观察似然比检测方法的检测概率、虚警概率、漏警概率的变化
p0=0.4;%改变此处的值完成第4.1问
p1=0.6;
for t=1:5000 %5000次仿真
%5.1改变信号的方差,观察似然比检测方法的检测概率、虚警概率、漏警概率的变化
n=5.*randn(1,200); %产生噪声 改变此处的值完成第5.1问 改变标准差的值
x=s+n;%产生含有噪声的信号
%利用似然比检测方法(最小错误概率准则),对信号是否到达进行检测
x1=x.*s;
x2=s.*s;
if 2*sum(x1)-sum(x2)>2*25*log(p0/p1)
tally1(t)=1; %检测到信号
' 未检测到信号 '
else
tally1(t)=0; %未检测到信号
' 检测到信号 '
end;
end;
n1=sum(tally1);%检测到信号的次数
n0=5000-n1;
Pd=n1/5000;%检测概率
Pm=n0/5000;%漏警概率
%计算虚警概率程序
for t=1:5000 %5000次仿真
n=5.*randn(1,200); %产生噪声 改变此处的值完成第5.1问
x=n;%产生含有噪声的信号
x1=x.*s;
x2=s.*s;
if 2*sum(x1)-sum(x2)>2*25*log(p0/p1)% 虚警概率 由于噪声总是客观存在的,当噪声信号的幅度超过检测门限时,
%雷达(或其他检测系统)就会被误认为发现目标,这种错误称为"虚警",它的发生概率称为虚警概率。
tally1(t)=1; %检测到信号
' 未检测到信号 ';
else
tally1(t)=0; %未检测到信号
' 检测到信号 ';
end;
end;
n2=sum(tally1);
Pf=n2/5000;
P=[Pd Pm Pf];
'检测概率Pd, 漏警概率Pm,虚警概率Pf='
P