%?一維距離像仿真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