ca-cfar的matlab仿真实现

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=-16;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‘011N）;
x=bragg+nosic+signal;
%%%
Fx=abs（fftshift（fft（x）））;
freq=（-1/（2*T）:1/（