%Samp7_6
wp=[0.35 0.65];N=48;    %通帶邊界頻率（歸一化頻率）和濾波器階數
Fs=50;
b=fir1(N,wp);           %設計FIR帶通濾波器
figure(1)
[H,f]=freqz(b,1,512,Fs);        %以50Hz為采樣頻率求出濾波器頻率響應
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=1;f2=10;f3=20;         %輸入信號的三種頻率成分
t=0:1/50:3;               %時間序列
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t)+0.5*sin(2*pi*f3*t);%輸入信號
%y=filter(b,1,x);          %可以采用過濾器進行濾波
y=fftfilt(b,x);            %采用 fftfilt 對輸入信號濾波
figure(2)
subplot(2,1,1), plot(t,x),title('輸入信號')%繪出輸入信號波形
subplot(2,1,2),plot(t,y)                  %繪出輸出信號波形
hold on;plot(N/2*0.02*ones(1,2),ylim, 'r')    %繪制延遲到的時刻
title('輸出信號'),xlabel('時間/s')