Matlab与信号处理知识.pptVIP

（2） 基于MATLAB的IIR滤波器设计 Butterworth数字滤波器设计 椭圆数字滤波器设计 Chebyshev-I数字滤波器设计 Chebyshev-II数字滤波器设计 * Butterworth数字滤波器设计 MATLAB语句： [b,a]=butter(n,wn) (1) [b,a]=butter(n,wn,’ftype’) (2) [z,p,k]=butter(…) (3) [a,b,c,d]=butter(…) (4) 语句(1)~ (4) 都可以设计一个阶数为n，截止频率为 wn的数字低通滤波器。对于语句(2) ，若参数ftype取‘hi- gh’或‘stop’，则可以来设计相应的高通或带阻滤波器。 返回值a和b为系统函数的分子和分母的系数。 * 系统函数为： 返回值z,p,k分别是滤波器的零点、极点和增益。 返回值A、B、C、D构造滤波器的状态方程。状态方程为： 其中u是输入信号，x是状态变量，y是输出信号。 * 例 ： 设计一个阶数为9阶、截止频率为300Hz的高通 Butterworth滤波器，采样频率为1000。并画出它的频率特性。 程序： n=9; wn=300/500; [b,a]=buffer(n,wn,’high’); freqz(b,a,128,1000); * Buttordworth高通滤波器的频率响应 * n=2 * 例 设计一个10阶的带通Butterworth滤波器，它的通带 范围是100~200Hz,采样频率为1000。并画出它的频率特性和冲击响应。 程序：n=10; wn=[100 200]/500; [b,a]=butter(n,wn); freqz(b,a,128,1000); [y,t]=impz(b,a,101); stem(t,y); * Buttordworth带通滤波器的频率响应 * Buttordworth滤波器的时域冲击响应 * 椭圆数字滤波器设计 MATLAB语法： [b,a]=ellip(n,Rp,Rs,wn) (5) [b,a]=ellip(n, Rp,Rs, wn,’ftype’) (6) [z,p,k]=ellip(…) (7) [a,b,c,d]=ellip(…) (8) 语句(5)~ (8)都可以设计一个阶数为n，截止频率为 wn,通带波纹最大衰减为Rp、阻带波纹最小衰减为Rs的 数字低通滤波器。对于语句(6)，若参数ftype取‘high’或 ‘stop’，则可以来设计相应的高通或带阻滤波器。 返回值a和b为系统函数的分子和分母的系数。返回 * 值z,p,k分别是滤波器的零点、极点和增益。返回值A、B、 C、D可以用来构造滤波器的状态方程。 例 采样频率为1000Hz，请设计一个阶数为6阶、截止频 率为300Hz的低通ellip数字滤波器。其中滤波器在通 带的波纹为3dB，阻带的波纹为50dB。 程序： [b,a]=ellip(6,3,50,300/500); freqz(b,a,512,1000);