%?用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