雷达CFAR恒虚警检测一维蒙特卡洛仿真

% 一维距离像仿真CA-CFAR目标检测过程
% 进行蒙特卡洛仿真，看看不同信噪比下的检测概率


clear all;close all;clc;

%% 雷达参数设定
c=3e8;
fc=1e9;
lamda=c/fc;

Tp=10e-6;  %脉冲宽度
B=20e6;    %带宽
kr=B/Tp;   %调频斜率
fs=30e6;   %快时间采样率
ts=1/fs;   %时间采样间隔
pulse_sample_num=floor（Tp*fs）;  %脉冲采样点数

R=10e3;    %目标距离10km
delay=2*R/c;  %时延
t=（delay-pulse_sample_num*ts）:ts:（delay+2*pulse_sample_num*ts）;
Nr=length（t）;
t_start=t（1）;
r_axis_real=t*c/2;     %绝对距离坐标轴
%% 信号产生、脉冲压缩、加噪
receiver_sig= rectpuls（t-Tp/2-delayTp）.*exp（1i*2*pi*fc*（-delay））.*exp（1i*pi*kr*（t-Tp/2-delay）.^2）;  
figure
plot（tabs（receiver_sig））;

st=rectpuls（t-t_start-Tp/2Tp）.*exp（1i*pi*kr*（t-t_start-Tp/2）.^2）;   %参考信号，实质时间是从0开始的，所以压缩完后反映的是真实距离

fft_st=fft（st）;   %参考信号fft
receive_fft=fft（receiver_sig）;   %接收信号fft 
data_r=ifft（receive_fft.*conj（fft_st））;

Amp_max=max（abs（data_r））;   %压缩后信号的最大幅度
SNR=0:25;  %压缩后信噪比
mont_num=1000;   %蒙特卡洛次数

%% CFAR自适应门限检测
num_cankao=20;   %参考单元
num_protect=4;   %保护单元
Pfa=1e-6;   %恒虚警率
x_axis=（num_cankao/2+num_protect/2+1）:（Nr-num_cankao/2-num_protect/2）;
alpha=num_cankao