%fliter_example2.m
%單位沖擊信號(hào)通過帶阻濾波器
t = 1:40;                                  %信號(hào)的時(shí)間
x =zeros(size(t));                         
x(1) = 1;                                  %產(chǎn)生單位沖擊信號(hào)
[b,a] = butter(10,[0.3 0.7],'stop');               %設(shè)計(jì)帶阻濾波器
y = filter(b,a,x);                            %濾波
hold on;                                  %畫圖濾波前后的時(shí)域數(shù)據(jù)
stem(t,x,'marker','o');
stem(t,y,'marker','.');
xlabel('時(shí)間/ts(s)');
legend('原始數(shù)據(jù)','濾波后數(shù)據(jù)'); 
fx=fft(x);                                 %對(duì)單位沖擊信號(hào)進(jìn)行傅立葉變換  
fx=fftshift(fx);
fy=fft(y);                                 %濾波后的信號(hào)進(jìn)行傅立葉變換
fy=fftshift(fy);
figure;
subplot(2,1,1);                            %畫圖顯示濾波前后的數(shù)據(jù)頻譜
f=(t-20)/20;
plot(f,abs(fx),'b-',f,abs(fy),'b-.');              %畫信號(hào)的幅頻曲線
xlabel('數(shù)字頻率/\pi(rad)');
title('幅頻曲線');
legend('原始數(shù)據(jù)的幅頻曲線譜','濾波后數(shù)據(jù)的幅頻曲線');
subplot(2,1,2);  
plot(f,angle(fx),'b-',f,angle(fy),'b-.');          %畫信號(hào)的相頻曲線
xlabel('數(shù)字頻率/\pi(rad)');
title('相頻曲線');
legend('原始數(shù)據(jù)的相頻曲線','濾波后數(shù)據(jù)的相頻曲線');