MATLAB程序计实验报告.docVIP

MATLAB程序设计实验报告 学 院：地球物理与石油资源学院 班 级： 物探（基）10901 姓 名： 刘豪 合 作 人： 侯颖 学 号： 200906170 班内编号： 09 指导教师： 陈义群 完成日期： 2012年4月20日 一、实验项目：IIR数字滤波器的设计法 二、时间地点：2012/4/17-2012/4/21 三、实验内容 1、数字滤波器知识 数字滤波器与模拟滤波器相对应，在数字信号处理中，如果一个离散时间系统是用来对输入信号进行滤波处理，那么该系统就可以称为数字滤波器（Digital Filter）。因此，数字滤波器实际上就是离散时间系统。 数字滤波器的基本输入/输出关系为： 其中h(n)为滤波器的单位抽样响应。 数字滤波器的转移函数为: 如果A(z)=1,则称滤波器为（有限冲激响应）FIR滤波器；否则为（无限冲激响应）IIR滤波器。令 ,得到滤波器的频率响应为： 定义为滤波器的幅频响应，而 为滤波器的相频响应。 2、IIR数字滤波器的频率变换设计法 IIR数字滤波器的频率变换设计法的基本思想是，由原型低通滤波器，通过频率变换形成不同形式的数字滤波器。其实现过程分为两种：第一种，首选根据滤波器设计要求，设计模拟原型低通滤波器，然后进行频率变换，将其转换为相应的数字滤波器；第二种，首选根据滤波器的设计要求，设计模拟原型低通滤波器，然后利用冲激响应不变法或双线性变换法，将原型低通滤波器数字化成数字原型低通滤波器，最后进行频率变换，将数字原型低通滤波器转换为相应的数字滤波器（高通、带通等）。 MATLAB的典型设计 为了实现频率变换设计IIR数字滤波器，MATLAB在信号处理具箱中提供了大量的IIR数字滤波器设计的相关函数，其中包括IIR滤波器阶次估计函数、模拟低通滤波器原型设计函数，以及模拟滤波器变换函数。每个函数有多钟调用方法 Buttord 计算Butterworth滤波器的阶次和截止频率 Cheb1lord 计算ChebyshevI型滤波器的阶次 Cheb2lord 计算ChebyshevII型滤波器的阶次 elipord 计算椭圆滤波器最小阶次 besselap Bessel模拟低通滤波器原型设计 buttap Butterworth模拟低通滤波器原型设计 cheb1ap Chebyshev I型模拟低通滤波器原型设计 cheb2ap Chebyshev II型模拟低通滤波器原型设计 elipap 椭圆模拟低通滤波器原型设计 lp2bp 把低通模拟滤波器转换成带通滤波器 lp2bs 把低通模拟滤波器转换成带阻滤波器 lp2hp 把低通模拟滤波器转换成高通滤波器 lp2 改变低通模拟滤波器的截止频率 因此，利用这些函数，IIR数字滤波器的频率变化设计将变得非常简单。利用MATLAB设计IIR数字滤波器可以分为以下几类来实现。 ① 按一定规则将数字滤波器的技术指标转换为模拟低通滤波器的技术指标[a] ② 根据转换后的技术指标使用滤波器阶数函数，确定滤波器的最小阶数N和截止频率 。 ③ 利用最小阶数N产生模拟低通滤波器原型。 ④ 利用截止频率 把模拟低通滤波器原型转换成模拟低通、高通、带通或者带阻滤波器。 ⑤ 利用冲激响应不变法或者双线性变换法把模拟滤波器转换成数字滤波器。 四、实验结论 1、MATLAB的典型设计 （1）、设计一个数字信号处理系统，它的采样率为100Hz，希望在该系统中设计一个Butterworth型高通数字滤波器，使其通带中允许的最大衰减为0.5dB，阻带内的最小衰减为40dB，通带上限临界频率为30Hz，阻带下限临界频率为40Hz。 程序如下： （2）、利用冲激响应不变法设计一个低通ChebyshevI型数字滤波器，其通带上限临界频率为0.3Hz，阻带临界频率为0.4Hz，采样频率为1000Hz，在通带内的最大衰减为