clc
clear all;
close all;
%{
本程序实现的功能：
CA-CFAR,GO-CFAR,SO-CFAR,OS-CFAR的检测概率与SNR_dB的关系
程序的参考公式来自于文献：基于杂波图的恒虚警处理技术研究
Detection Loss Due to Interfering Targets in Ordered Statistics CFAR
%}
R=32; %参考单元长度
n=R/2; %半滑窗长度
k=3*R/4;
pfa=10^(-6); %虚警概率
%% CA-CFAR检测概率和SNR的关系
T_ca=pfa.^(-1./R)-1; %阈值系数
SNR_dB=0:30;
SNR=10.^(SNR_dB./10);
pd_ca=(1+(T_ca./(1+SNR))).^(-R); %CA-CFAR检测概率
%% OS-CFAR检测概率和SNR的关系
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);
T_os=T(T==abs(T)); %OS-CFAR检测门限
pd_os=k*nchoosek(R,k)*gamma(R-k+1+T_os./(1+SNR))*gamma(k)./gamma(R+T_os./(1+SNR)+1);
%%
pd=0;
syms T4;
for i=1:n
pd=pd+2.*nchoosek(n+i-2,i-1).*(2+T4).^(-n-i+1);
end
T_SO=solve(pd==pfa,T4);
T_SO=double(T_SO);
T_SO=T_SO(T_SO==abs(T_SO)); %求得SO-CFAR的门限因子。
pfa_go=2.*(1+T4).^(-n)-pd;
T_GO=solve(pfa_go==pfa,T4);
T_GO=double(T_GO);
T_GO=T_GO(T_GO==abs(T_GO)); %求得GO-CFAR的门限因子。
pd_so=0;
pd_go
for i=1:n
pd_so=pd_so+2.*nchoosek(n+i-2,i-1).*(2+T_SO./(1+SNR)).^(-n-i+1);
pd_go=pd_go+2.*nchoosek(n+i-2,i-1).*(2+T_GO./(1+SNR)).^(-n-i+1);
%SO-CFAR的检测概率
end
pd_go=2.*(1+T_GO./(1+SNR)).^(-n)-pd_go; %GO-CFAR的检测概率
plot(SNR_dB,pd_ca,'r');
hold on;
plot(SNR_dB,pd_go,'b');
hold on;
plot(SNR_dB,pd_so,'y');
hold on;
plot(SNR_dB,pd_os,'g');
grid on
