CA-CFAR代码实现

clear all;clc;close all;
%%%% 仿真数据 %%%%
T=0.25;M=100;tc=T/M;fc=1/tc;%%% 一个扫频周期T 为250ms ； M为每个周期采样点100 ； fc采样频率
N=256;n=1:N;%% 一个相关积累时间64
%%%
fw=10;%% 雷达工作频率10MHz
fb=0.102*sqrt(fw);%% bragg 频率
fs=1;%% 信号频率
a_db=-10;a=10^(a_db/20);%% 正Bragg峰幅值
b_db=-12;b=10^(b_db/20);%% 负Bragg峰幅值
c_db=-35;c=10^(c_db/20);%% 噪声基底

s_db=-30;s=10^(s_db/20);%% 信号幅度
bragg=a.*exp(j*2*pi*fb*T.*n)+b.*exp(-j*2*pi*fb*T*n);%%% 正弦序列产生 Bragg 峰
signal=s.*exp(-j*2*pi*fs*T*n);
nosic=c.*random(‘Normal’,0,1,1,N);
x=bragg+nosic+signal;
%%%
Fx=abs(fftshift(fft(x)));
freq=(-1/(2*T):1/(T*N):(1/(2*T)-1/(T*N)));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% CA-CFAR %%%%%%%%%%%%%%%%%%%%%%%%%
M=N;P_fa=10.^(-6);
R=16;%参考单元数
n=R/2;
L_slipper=R+1;%