4-3.m

Matlab数字信号处理与应用源代码

%例程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)|')