MTD with FFT and FIR

% 用FFT和FIR两种方法实现MTD

close all;clear all;clc;
fz = 10e3;          % 帧频率（雷达信号的重频）
tz = 1/fz;          % 帧周期（雷达信号的重复周期）
fs = 8e6;           % 采样频率
ts = 1/fs;        
f_doppler = 4e3;  % 多目标的多普勒频率
N = round（tz/ts）;

% 产生雷达回波信号
s_pc_1 = [zeros（1200）11zeros（1N-202）];
s_pc = repmat（s_pc_1116）;
n = 1:16*N;
s_doppler = cos（n*f_doppler/fs*2*pi）;
s_pc = s_pc.*s_doppler;
s_noise = 0.1*randn（1N*16）;
s_pc = s_pc+s_noise;
figure;plot（0:ts:（16*N-1）*tss_pc‘.‘）;xlabel（‘t单位：S‘）title（‘回波信号‘）;
s_pc = reshape（s_pcN16）;


% 采用FFT法进行MTD处理
M =16;
s_mtd =zeros（N16）;
for i = 1:N
    s_temp = s_pc（i:）;
    s_mtd（i:） = fft（s_tempM）;
end
s_mtd1=abs（s_mtd）
s_mtd2=s_mtd1‘
for i = 1:N
    fir_accumulation_result（i） = sum（s_mtd2（:i））/16;
end
figure;
plot（[0:ts:（N-1）*ts]fir_accumulation_result）;xlabel（‘t‘）;title（‘非相参积累处理结果1‘）
for z=1:1:16
s_MTD_result=s_mtd2（z:）
N = 800;
cfar_result（z1） = s_MTD_result（11）/（sqrt（2）/pi*mean（s_MTD_result（12:17）））;
for i = 2:16
    noise_mean = sqrt（2）/pi*（mean（s_MTD_result（11:i-1））+mean（s_MTD_result（1i+1:i+16）））/2;
    cfar_result（zi） = s_MTD_result（1i）/noise_mean;
end
for i = 17:N-17
    noise_mean = sqrt（2）/pi*max（mean（s_MTD_result（1i-16:i-1））mean（s_MTD_result（1i+1:i+16）））;
    cfar_result（zi） = s_MTD_result（1i）/noise_mean;
end
for k = N-16:N-1
    noise_mean = sqrt（2）/pi*（mean（s_MTD_result（1k-16:k-1））+mean（s_MTD_result（1k+1:N）））/2;
    cfar_result（zk） = s_MTD_result（1k）/noise_mean;
end
cfar_result（zN） = s_MTD_result（1N）/（sqrt（2）/pi*mean（s_MTD_result（1N-16:N-1）））;
end
for i = 1:N
    fir_accumulation_result（i） = sum（cfar_result（:i））/16;
end
figure;
plot（[0:ts:（N-1）*ts]fir_accumulation_result）;xlabel（‘t‘）;title（‘非相参积累处理结果‘）

figure;
plot（0:ts:tz-tscfar_result）
xlabel（‘t（单位:S）‘）;title（‘采用恒虚警率处理后的结果‘）;
figure;
t = 0:ts:tz-ts;
subplot（221）;plot（tabs（s_mtd（:1）））;title（‘多普勒滤波器组（第1通道）‘）
subplot（222）;plot（tabs（s_mtd（:2）））;title（‘多普勒滤波器组（第2通道）‘）
subplot（223）;plot（tabs（s_mtd（:3）））;title（‘多普勒滤波器组（第3通道）‘）
subplot（224）;plot（tabs（s_mtd（:4）））;title（‘多普勒滤波器组（第4通道）‘）
figure;
subplot（221）;plot（tabs（s_mtd（:5）））;title（‘多普勒滤波器组（第5通道）‘）
subplot（222）;plot（tabs（s_mtd（:6）））;title（‘多普勒滤波器组（第6通道）‘）
subplot（223）;plot（tabs（s_mtd（:7）））;title（‘多普勒滤波器组（第7通道）‘）
subplot（224）;plot（tabs（s_mtd（:8）））;title（‘多普勒滤波器组（第8通道）‘）
figure;
subplot（221）;plot（tabs（s_mtd（:9）））;title（‘多普勒滤波器组（第9通道）‘）
subplot（222）;plot（tabs（s_mtd（:10）））;title（‘多普勒滤波器组（第10通道）‘）
subplot（223）;plot（tabs（s_mtd（:11）））;title（‘多普勒滤波器组（第11通道）‘）
subplot（224）;plot（tabs（s_mtd（:12）））;title（‘多普勒滤波器组（第12通道）‘）
figure;
subplot（221）;plot（tabs（s_mtd（:13）））;title（‘多普勒滤波器组（第13通道）‘）
subplot（222）;plot（tabs（s_mtd（:14）））;title（‘多普勒滤波器组（第14通道）‘）
subplot（223）;pl