4-4.m

来自「数字信号处理相关源代码」· M 代码 · 共 12 行

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