%?一維距離像仿真CA-CFAR目標(biāo)檢測(cè)過(guò)程
%?進(jìn)行蒙特卡洛仿真，看看不同信噪比下的檢測(cè)概率


clear?all;close?all;clc;

%%?雷達(dá)參數(shù)設(shè)定
c=3e8;
fc=1e9;
lamda=c/fc;

Tp=10e-6;??%脈沖寬度
B=20e6;????%帶寬
kr=B/Tp;???%調(diào)頻斜率
fs=30e6;???%快時(shí)間采樣率
ts=1/fs;???%時(shí)間采樣間隔
pulse_sample_num=floor（Tp*fs）;??%脈沖采樣點(diǎn)數(shù)

R=10e3;????%目標(biāo)距離10km
delay=2*R/c;??%時(shí)延
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;?????%絕對(duì)距離坐標(biāo)軸
%%?信號(hào)產(chǎn)生、脈沖壓縮、加噪
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）;???%參考信號(hào)，實(shí)質(zhì)時(shí)間是從0開(kāi)始的，所以壓縮完后反映的是真實(shí)距離

fft_st=fft（st）;???%參考信號(hào)fft
receive_fft=fft（receiver_sig）;???%接收信號(hào)fft?
data_r=ifft（receive_fft.*conj（fft_st））;

Amp_max=max（abs（data_r））;???%壓縮后信號(hào)的最大幅度
SNR=0:25;??%壓縮后信噪比
mont_num=1000;???%蒙特卡洛次數(shù)

%%?CFAR自適應(yīng)門(mén)限檢測(cè)
num_cankao=20;???%參考單元
num_protect=4;???%保護(hù)單元
Pfa=1e-6;???%恒虛警率
x_axis=（num_cankao/2+num_protect/2+1）:（Nr-num_cankao/2-num_protect/2）;
alpha=num_cankao