雷达恒虚警算法之OS—CFAR的matlab仿真及代码（蒙特卡洛）

  有序统计量（OS）-CFAR的仿真有三个模块：



有序统计量（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