%?一維距離像仿真CA-CFAR目標(biāo)檢測(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???假設(shè)目標(biāo)不動(dòng)，這樣就沒有了多普勒
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開始的，所以壓縮完后反映的是真實(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=23;????%壓縮后信噪比
sigma_noise=Amp_max*10^（-SNR/20）;?????%噪聲σ，信號(hào)功率Amp_max^2;所以信噪比SNR=20log（Amp_max/σ）
real_noise=randn（size（t））;
imag_noise=randn（size（t））;
noise=sigma_noise*（real_noise+1i*imag_noise）/2^0.5;????%產(chǎn)生額定信噪比的噪聲
data_r=data_r+noise;
figure
plot（ab