clear;
clc;
close all;
fr=1e3;%脈沖重復(fù)頻率
tr=1/fr;%脈沖重復(fù)周期
fs=1e6;%采樣頻率
ts=1/fs;%采樣周期
%%%%%%%%%%%%%%%%%%%%%%%%%生成測試模塊%%%%%%%%%%%%%%%%%%%%%%
Point=4000;
doppler_vector=2*pi*linspace(0,fr,Point);
t_doppler=(1:8)*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;


%%%%%%%%%%%%%%%%%%%%  產(chǎn)生8個(gè)FIR濾波器的濾波器組  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
a0=[ 7.135916968387e-010,1.568395747901e-008,5.471539757694e-008,9.376656588003e-008,...
  9.376656588003e-008,5.471539757694e-008,1.568395747901e-008,7.135916968387e-010]*1e10;
n=1:8;
f=linspace(0.61,1.388,8)*pi;
fn=f'*n;
fnexp=exp(j*fn);
a08=repmat(a0,1,8);
a08=reshape(a08,8,8)';
b=a08.*fnexp;
%%%%%%%% FIR濾波器的濾波器組的頻率響應(yīng)曲線  %%%%%%%%%%%%
figure('Name','FIR濾波器的濾波器組的頻率響應(yīng)曲線');
for n=1:8
    freqz(b(n,:),1,1024,'whole',fr)
    hold on;
end
axis([0 fr -50 100]);
%%%%%%%% FIR濾波器的濾波器組的頻率響應(yīng)曲線  %%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%  產(chǎn)生16個(gè)FIR濾波器的濾波器組  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

yt=test_matrix;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  采用FIR進(jìn)行MTD處理  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
M=8;
L=length(yt);
for i=1:L
    s_temp=yt(:,i); 
    s_mtd(:,i)=b*s_temp;    
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  采用FIR進(jìn)行MTD處理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

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('Name','實(shí)測FIR濾波器組的頻率響應(yīng)曲線 1');

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

figure('Name','實(shí)測FIR濾波器組的頻率響應(yīng)曲線 2');

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
axis([0,fr,-70,0])

figure('Name','實(shí)測FIR濾波器組的頻率響應(yīng)曲線包絡(luò)');
plot(t,20*log10(max(abs(s_mtd))/max(max(abs(s_mtd)))),'-r');grid on;