4-5.m

matlab数字信号与应用源代码4

%例程4－5    基于巴特沃斯模拟滤波器设计数字带阻滤波器

Wp=[0.4*pi,0.7*pi];
Ws=[0.2*pi,0.8*pi];
Ap=1;
As=30;
[N,wn]=buttord(Wp/pi,Ws/pi,Ap,As); %计算巴特沃斯滤波器阶次和截止频率
[b,a]=butter(N,wn,'stop');      %频率变换法设计巴特沃斯带通滤波器
[db,mag,pha,grd,w]=freqz_m(b,a);   %数字滤波器响应
plot(w/pi,mag);
Title('数字滤波器幅频响|H(ejOmega)|')