雷达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=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