%例程4－4    基于巴特沃斯模擬濾波器設計數字帶通濾波器

Wp=[0.3*pi,0.4*pi];
Ws=[0.2*pi,0.5*pi];
Ap=3;
As=18;
[N,wn]=buttord(Wp/pi,Ws/pi,Ap,As); %計算巴特沃斯濾波器階次和截止頻率
[b,a]=butter(N,wn,'bandpass');      %頻率變換法設計巴特沃斯帶通濾波器
[db,mag,pha,grd,w]=freqz_m(b,a);   %數字濾波器響應
Plot(w/pi,mag);
Title('數字濾波器幅頻響應|H(ejOmega)|')