%Samp7_8
clf;N=300;Fs=1000;         %數據總數和采樣頻率
Order=100;    %濾波器階數
n=[0:N-1];t=n/Fs;     %時間序列
wn=[100 200]/(Fs/2); %邊界頻率轉換為歸一化頻率，據6-19式
b=fir1(Order,wn,'stop');          % 設計 100 階的帶阻濾波器
figure(1)
[H,f]=freqz(b,1,512,Fs);
subplot(2,1,1),plot(f,20*log10(abs(H)))
xlabel('頻率/Hz');ylabel('振幅/dB');grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel('頻率/Hz');ylabel('相位/^o');grid on;
f1=50;f2=150;  %輸入信號頻率
x=sin(2*pi*f1*t)+sin(2*pi*f2*t);       %輸入信號
y=fftfilt(b,x,256);% 對數據進行濾波
figure(2)
subplot(2,1,1);plot(t,x); grid on;                %繪制輸入信號
title('輸入信號');
subplot(2,1,2);plot(t,y);                   %繪制輸出信號
hold on;plot(Order/2/Fs*ones(1,2),ylim, 'r')    %繪制延遲到的時刻
title('輸出信號');
grid on;  xlabel('時間/s')