%?一維距離像仿真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=23;????%壓縮后信噪比
sigma_noise=Amp_max*10^（-SNR/20）;?????%噪聲σ，信號功率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;????%產生額定信噪比的噪聲
data_r=data_r+noise;
figure
plot（ab