%Samp7_5
wp=0.5*pi;ws=0.66*pi;       %濾波器的邊界頻率
wdelta=ws-wp;            %過渡帶寬度
N=ceil(8* pi/wdelta);%求解濾波器的最小階數，根據Hanning 窗主瓣寬
Wn=(0.5+0.66)*pi/2;%截止頻率取通帶和阻帶邊界頻率的中點
b=fir1(N,Wn/pi,hanning(N+1));%設計FIR濾波器,注意fir1要求輸入歸一化頻率
[H,f]=freqz(b,1,512,50); %采用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;