clear;
clc;
close all;
fr=1e3;%脈沖重復頻率
tr=1/fr;%脈沖重復周期
fs=1e6;%采樣頻率
ts=1/fs;%采樣周期
%%%%%%%%%%%%%%%%%%%%%%%%%生成測試模塊%%%%%%%%%%%%%%%%%%%%%%
Point=4000;
doppler_vector=2*pi*linspace(0,fr,Point);
t_doppler=(1:18)*tr;
test_matrix=t_doppler'*doppler_vector;
test_matrix_i=cos(test_matrix);
test_matrix_q=sin(test_matrix);
test_matrix=test_matrix_i+test_matrix_q*j;

%%%%%%%%%%%%%%固定雜波對消處理%%%%%%%%%%%%%%%%%%%%
yt=test_matrix;
for i=1:16
       % s_mti(i,:)=yt(i,:)+yt(i+2,:)-2*yt(i+1,:);%三脈沖對消
   s_mti(i,:)=yt(i,:)-yt(i+1,:);%兩脈沖對消
end
%s_mti=yt(1:16,:);
%%%%%%%%%%%%%%%%%%%采用FFT進行MTD處理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
M=16;
L=length(s_mti);
for i=1:L
  % s_temp=s_mti(:,i);
  % s_temp=s_mti(:,i).*chebwin(M);%切比雪夫 加權(quán)FFT
    s_temp=s_mti(:,i).*hamming(M);%Hamming 加窗FFT
    s_mtd(:,i)=fft(s_temp,M);
end;
mtd_abs=abs(s_mtd);
mtd_abs=mtd_abs/max(max(mtd_abs));
mtd_abs=20*log10(mtd_abs+eps);
t=linspace(0,fr,Point);
figure

plot(t,abs(s_mtd(1,:)),'-r'),hold on
plot(t,abs(s_mtd(2,:)),'-g'),hold on
plot(t,abs(s_mtd(3,:)),'-b'),hold on
plot(t,abs(s_mtd(4,:)),'-c'),hold on
plot(t,abs(s_mtd(5,:)),'-m'),hold on;
plot(t,abs(s_mtd(6,:)),'-y'),hold on;
plot(t,abs(s_mtd(7,:)),'-k'),hold on
plot(t,abs(s_mtd(8,:)),'-c'),hold on
plot(t,abs(s_mtd(9,:)),'-r'),hold on
plot(t,abs(s_mtd(10,:)),'-g'),hold on
plot(t,abs(s_mtd(11,:)),'-b'),hold on
plot(t,abs(s_mtd(12,:)),'-c'),hold on
plot(t,abs(s_mtd(13,:)),'-m'),hold on
plot(t,abs(s_mtd(14,:)),'-y'),hold on
plot(t,abs(s_mtd(15,:)),'-k'),hold on
plot(t,abs(s_mtd(16,:)),'-c');grid on;
figure

plot(t,mtd_abs(1,:),'-r'),hold on
plot(t,mtd_abs(2,:),'-g'),hold on
plot(t,mtd_abs(3,:),'-b'),hold on
plot(t,mtd_abs(4,:),'-c'),hold on
plot(t,mtd_abs(5,:),'-m'),hold on;
plot(t,mtd_abs(6,:),'-y'),hold on;
plot(t,mtd_abs(7,:),'-k'),hold on
plot(t,mtd_abs(8,:),'-c'),hold on
plot(t,mtd_abs(9,:),'-r'),hold on
plot(t,mtd_abs(10,:),'-g'),hold on
plot(t,mtd_abs(11,:),'-b'),hold on
plot(t,mtd_abs(12,:),'-c'),hold on
plot(t,mtd_abs(13,:),'-m'),hold on
plot(t,mtd_abs(14,:),'-y'),hold on
plot(t,mtd_abs(15,:),'-k'),hold on
plot(t,mtd_abs(16,:),'-c');grid on;
axis([0,fr,-70,0])

figure
plot(max(abs(s_mtd)))