clc,clear all,close all;
Wp=[0.2 0.3];                   %通帶兩個邊界角頻率
Ws=[0.1 0.4];                   %阻帶兩個邊界角頻率
Rp=-1;                          %通帶允許起伏
Rs=-25;                         %阻帶最小衰減  
[N,Wn]=buttord(Wp,Ws,Rp,Rs);    %用buttord確定低通原型參數
[B,A]=butter(N,Wn);             %用butter生成帶通濾波器
[H,w]=freqz(B,A,2*N);
subplot(2,1,1);
plot(w/pi,abs(H));
xlabel('角頻率（單位pi）');
ylabel('濾波器幅度');
title('幅頻特性');
subplot(2,1,2);
plot(w/pi,angle(H));
xlabel('角頻率（單位pi）');