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

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