赞
踩
有序统计量(OS)-CFAR公式的参数设计;T(检测门限)的求解;检测概率Pd的求解;代码如下:
%%%%%%%%%%%%%%输入参数:Pf,R,k%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc;
clear all;
R=24; %参考单元长度
n=R/2; %半滑窗长度
k=3*R/4;
Pfa=1e-6; %虚警概率
SNR_dB=5:1:30; %设置信噪比(dB)
SNR=10.^(SNR_dB./10); %计算信噪比
figure
%%%%%%%%%%%%%%%%%%%%%%%%%%% %计算OS-CFAR门限系数%%%%%%%%%%%%%%%%%%%%%%%%%%
syms T
g=Pfa-k*nchoosek(R,k)*gamma(R-k+1+T)*gamma(k)/gamma(R+T+1);
x=solve(g);
T=double(x);
T1=T(T==abs(T));
%%%%%%%%%%%%%%%%%%%%%%计算恒虚警概率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Pd_SO=k*nchoosek(R,k)*gamma(R-k+1+T1./(1+SNR))*gamma(k)./gamma(R+T1./(1+SNR)+1);
hold on
plot(SNR_dB,Pd_SO,'g');
grid on
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。