数字电路综合性实验报告

重庆交通大学

综合性设计性实验报告

班级：

学号：

姓名：

实验项目名称：滤波器设计综合实验应用

实验项目性质：综合性设计

实验所属课程：数字信号处理

实验室(中心)：现代电子实验中心

指导教师：

实验完成时间：2014 年12 月28 日

一、实验目的

1.学生自己运用MATLAB设计IIR数字低通滤波器，方法不限。

2.实现信号的滤波。

3.熟悉用冲激响应不变法设计IIR数字滤波器的原理与方法

4.学会巴特沃什型滤波器的设计，各种参数的计算方式

二、实验主要内容及过程

1.设计一模拟IIR模拟低通滤波器并转换为IIR低通滤波器。

（1）模拟滤波器设计采用巴特沃什型滤波器作为原型。

（2）模拟到数字转换采用冲激响应法。

2.利用实现的滤波器对信号进行滤波。

提示：完成此步骤需要与任务一综合考虑。分析有用和噪音信号的频率，并参考噪音信号的频率利用采样定理等知识选取合理的滤波器截止频率进行任务一的设计。

3.过程：用冲激响应不变法设计Butterworth 低通数字滤波器，要求通带频率为0 ≤ω≤0.2π，通带波纹小于1dB，阻带在0.3π≤ω≤π内，幅度衰减大于15dB，采样周期T=0.01s。

（1）用冲激响应不变法设计该数字滤波器 H(z)。 （2）使用MATLAB 软件对滤波器性能进行分析。

（3）假设一个信号t f t f t X 212cos 5.02sin )(ππ+=，其中1f =5Hz,2f =30Hz 。试将原信号与经过该滤波器的输出信号进行比较。 三、设计方案 （一）题目

IIR 数字滤波器的设计 （二）设计的主要思路

1. 根据给定的性能指标和方法不同，首先对设计性能指标中的频率指标，

如数字边界频率进行变换，将其转为模拟频率并且此频率为模拟滤波器原型设计的性能指标。（Ω=ω/T ）及1Ω=20π，2Ω=30π

设计IIR 滤波器时，给出的性能指标通常分数字指标和模拟指标两种。数字性能指标给出通带截止频率p ω ，阻带起始频率s ω ，通带波纹Rp,阻带衰减Rs 等。数字频率p ω 和s ω 的取值范围为0~π ，单位弧度而MATLAB 工具函数常采用归一化频率，p ω 和s ω 的取值范围为0~1，对应于0～π ，此时需进行转换。模拟性能指标给出通带截止频率p ω，阻带起始频率s ω，通带波纹Rp ，阻带衰减Rs 等。模拟频率

p ω和s ω单位为弧度/秒（rad/s ）。

MATLAB 信号处理工具箱中，设计性能指标的转换应根据不同设计方法进行不同处理。

2. 估计模拟滤波器最小阶数和截止频率，利用MATLAB 工具函数buttord 。

3. 设计模拟低通滤波器原型，利用MATLAB 工具函数buttap 。

4. 由原型低通滤波器经频率变换等模拟滤波器，利用MATLAB 工具函数

lp2lp 。

5. 将模拟滤波器离散化获得IIR 数字滤波器，利用MATLAB 工具函数

impinvar 。

6. 编写代码画出该滤波器的幅度与相位特性曲线，用freqz ，impz ，zplane

函数对滤波器性能进行分析。

7. 编写x （t ）原信号代码，调用所设计的滤波器函数得出滤波后的信号图，

并进行比较。 8. 完成设计性实验报告 （三）设计的计算过程

注：前4个步骤可以通过计算完成，实现巴特沃什低通滤波器的设计。在继续后面4个步骤即可完成实验。

根据公式

（3-1）

计算出滤波器的阶次等于N=6， 根据公式

N

k c 2/11.01)110/(1-Ω=Ω- （3-2）

计算得到c Ω=70.29

)/lg(2)]

110/()110lg[(211.01.021ΩΩ--≥

--k k N

表3-1 各阶巴特沃斯因式分解多项式)(s B n

根据查表3-1 各阶巴特沃斯因式分解多项式得到该阶次归一化（即c Ω=1）的巴特沃斯低通原型滤波器的系统函数为

（3-3）

因为c Ω=70.29时,用s/c Ω对)(s H n 中的s 进行置换即可得到要设计的巴特沃斯低通滤波器。 计算详细过程如下图所示：

)932.11)(21)(517.01(1

)(2226s s s s s s s H ++++++=

图4-1 计算详细过程图

(四)设计程序用到的函数

1.buttord

用来估计模拟滤波器最小阶数和截止频率，

[N,Wc]=buttord(Wp,Ws,rp,rs,'s')

2.Buttap

设计模拟低通滤波器原型，[z,p,k]=buttap(N)

3.lp2lp

由模拟原型低通滤波器经频率变换获得模拟低通滤波器，

[b,a]=lp2lp(Bap,Aap,Wc)

4.Impinvar

MATLAB 中提供的冲激响应不变法设计数字滤波器的函数，调用格式为：

[bz,az]=impinvar(b,a[[,Fs[，Fp])

式中，b,a 为模拟滤波器分子和分母多项式系数向量；Fs 为采样频率（所滤波数据），单位Hz，缺省时为1Hz。Fp 为预畸变频率（Prewarped frequency）,是一个“匹配”频率,在该频率上，频率响应在变换前后和模拟频率可精确匹配。一般设计中可以不考虑。bz,az 分别为数字滤波器分子和分母多项式系数向量。前面已提到过，函数输入变量中的[]表示可添加也可略去的内容。

5.Freqz

函数freqz 用于求数字滤波器的频率响应，其调用格式为：

[[h,w]=]freqz(b,a,n[,’whole’]);

或[h,f]=freqz(b,a,n[,’whole’],Fs);

式中，b,a 为数字滤波器分子和分母多项式的系数，n 为复数频率的响应点数，为整数，最好为2 的幂，缺省时为512；Fs 为采样频率，单位Hz。如果给定该值，则f 位置输出为频率Hz，若没有给定，则按角频率（Angular frequency）给定f 的频率矢量；’whole’表示返回的频率f 或w 值包含z 平面整个单位圆频率矢量，即0~2π；缺省时，频率f 或w 值包含z 平面上半单位圆（0~π）之间等间距n 个点频率矢量。h 为复频率响应；w 为n 点频率向量（单位rad）；f 为n 点频率向量(Hz)。

6.Impz

impz 用于产生数字滤波器的冲激响应。调用格式为：