第五章 数字滤波器结构 DF (Digital Filter)

Slides:



Advertisements
Similar presentations
3 的倍数的特征 的倍数有 : 。 5 的倍数有 : 。 既是 2 的倍数又是 5 的倍数有 : 。 12 , 18 , 20 , 48 , 60 , 72 , , 25 , 60 ,
Advertisements

一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
第八章 第四节 机动 目录 上页 下页 返回 结束 一个方程所确定的隐函数 及其导数 隐函数的微分法.
2.6 隐函数微分法 第二章 第二章 二、高阶导数 一、隐式定义的函数 三、可微函数的有理幂. 一、隐函数的导数 若由方程 可确定 y 是 x 的函数, 由 表示的函数, 称为显函数. 例如, 可确定显函数 可确定 y 是 x 的函数, 但此隐函数不能显化. 函数为隐函数. 则称此 隐函数求导方法.
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
信号与系统 第三章 傅里叶变换 东北大学 2017/2/27.
§3.4 空间直线的方程.
3.4 空间直线的方程.
圆的一般方程 (x-a)2 +(y-b)2=r2 x2+y2+Dx+Ey+F=0 Ax2+Bxy+Cy2+Dx+Ey+ F=0.
第五章 二次型. 第五章 二次型 知识点1---二次型及其矩阵表示 二次型的基本概念 1. 线性变换与合同矩阵 2.
第三章 函数逼近 — 最佳平方逼近.
第四章 数字滤波器基础 本章要点 数字滤波器 Z变换 数字滤波器的组成 数字滤波器的类型 差分方程的传递函数 Z平面的零-极点分布图
分式的乘除.
第2章 时域离散信号和系统的频域分析 教学内容包括: 序列的傅立叶变换定义及性质 Z变换的定义与收敛域 利用z变换分析信号和系统的频域特性.
第一章 绪论.
第一章 行列式 第五节 Cramer定理 设含有n 个未知量的n个方程构成的线性方程组为 (Ⅰ) 由未知数的系数组成的n阶行列式
内容提要 数字滤波器基本原理 IIR 数字滤波器 FIR 数字滤波器 数字滤波器的实现问题 数字滤波器的结构及有限字长效应
第四节 对数留数与辐角原理 一、对数留数 二、辐角原理 三、路西定理 四、小结与思考.
恰当方程(全微分方程) 一、概念 二、全微分方程的解法.
一、原函数与不定积分 二、不定积分的几何意义 三、基本积分公式及积分法则 四、牛顿—莱布尼兹公式 五、小结
第四章 函数的积分学 第六节 微积分的基本公式 一、变上限定积分 二、微积分的基本公式.
第5章 定积分及其应用 基本要求 5.1 定积分的概念与性质 5.2 微积分基本公式 5.3 定积分的换元积分法与分部积分法
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
第二章 导数与微分 第二节 函数的微分法 一、导数的四则运算 二、复合函数的微分法.
2-7、函数的微分 教学要求 教学要点.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
实验二 FFT与DFT计算时间的比较及圆周卷积代替线性卷积的有效性实验
1.4 离散时间系统与差分方程 y(n)= T[x(n)]
Signals and Systems Lecture 28
第5章 §5.3 定积分的积分法 换元积分法 不定积分 分部积分法 换元积分法 定积分 分部积分法.
第2章 Z变换 Z变换的定义与收敛域 Z反变换 系统的稳定性和H(z) 系统函数.
数字滤波器结构 DF (Digital Filter)
Matlab 中IIR数字滤波器设计相关函数
第二章 矩阵(matrix) 第8次课.
元素替换法 ——行列式按行(列)展开(推论)
实验四 滤波器传输函数的零点和极点 对滤波特性的影响
§2 求导法则 2.1 求导数的四则运算法则 下面分三部分加以证明, 并同时给出相应的推论和例题 .
第8章 数字滤波系统的网络结构与分析 课程名称:数字信号处理 任课教师:张培珍 授课班级:信计
第一章 函数与极限.
第四章习题.
晶体管及其小信号放大 -单管共射电路的频率特性.
概 率 统 计 主讲教师 叶宏 山东大学数学院.
Three stability circuits analysis with TINA-TI
线 性 代 数 厦门大学线性代数教学组 2019年4月24日6时8分 / 45.
语音信号的短时分析技术 对语音信号采用分段(或称分帧处理),称为短时分析。 短时平均能量 短时能量分析 窗口形状的选择 窗口的长度
晶体管及其小信号放大 -单管共射电路的频率特性.
第四章 数字滤波器的基本结构(3).
复习.
成绩是怎么算出来的? 16级第一学期半期考试成绩 班级 姓名 语文 数学 英语 政治 历史 地理 物理 化学 生物 总分 1 张三1 115
2019/5/4 实验三 离散傅立叶变换的性质及应用 06:11:49.
§6.7 子空间的直和 一、直和的定义 二、直和的判定 三、多个子空间的直和.
2019/5/11 实验四 FIR滤波器的特性及应用 05:31:12.
2019/5/11 实验三 线性相位FIR滤波器的特性 05:31:30.
第三章 函数的微分学 第二节 导数的四则运算法则 一、导数的四则运算 二、偏导数的求法.
第六节 用频率特性法分析系统性能举例 一、单闭环有静差调速系统的性能分析 二、单闭环无静差调速系统的性能分析
学习任务三 偏导数 结合一元函数的导数学习二元函数的偏导数是非常有用的. 要求了解二元函数的偏导数的定义, 掌握二元函数偏导数的计算.
第15讲 特征值与特征向量的性质 主要内容:特征值与特征向量的性质.
线 性 代 数 厦门大学线性代数教学组 2019年5月12日4时19分 / 45.
2019/5/21 实验一 离散傅立叶变换的性质及应用 实验报告上传到“作业提交”。 11:21:44.
§2 方阵的特征值与特征向量.
§7.3 离散时间系统的数学 模型—差分方程 线性时不变离散系统 由微分方程导出差分方程 由系统框图写差分方程 差分方程的特点.
欢迎大家来到我们的课堂 §3.1.1两角差的余弦公式 广州市西关外国语学校 高一(5)班 教师:王琦.
信号发生电路 -非正弦波发生电路.
定义5 把矩阵 A 的行换成同序数的列得到的矩阵,
第四节 向量的乘积 一、两向量的数量积 二、两向量的向量积.
第四章 函数的 积分学 第七节 定积分的换元积分法     与分部积分法 一、定积分的换元积分法 二、定积分的分部积分法.
第三节 数量积 向量积 混合积 一、向量的数量积 二、向量的向量积 三、向量的混合积 四、小结 思考题.
§4.5 最大公因式的矩阵求法( Ⅱ ).
Presentation transcript:

第五章 数字滤波器结构 DF (Digital Filter)

第一节 引 言

一、什么是数字滤波器 顾名思义:其作用是对输入信号起到滤波的作用;即DF是由差分方程描述的一类特殊的离散时间系统。 它的功能:把输入序列通过一定的运算变换成输出序列。不同的运算处理方法决定了滤波器的实现结构的不同。

二、数字滤波器的工作原理 x(n) y(n) h(n) 则LTI系统的输出为:

三、数字滤波器表示方法 有两 种表示方法:方框图表示法;流图表示法. 数字滤波器中,信号只有延时,乘以常数和相加三种运算。 所以DF结构中有三个基本运算单元:加法器,单位延时,乘常数的乘法器。

1、方框图、流图表示法 方框图表示法: 信号流图表示法: 单位延时 系数乘 相加 a a Z-1 Z-1 把上述三个基本单元互联,可构成不同数字网络或运算结构,也有方框图表示法和流图表示法。

2.例子 例:二阶数字滤波器: 其方框图及流图结构如下: b0 b0 y(n) x(n) y(n) x(n) a1 Z-1 a1 Z-1 看出:可通过流图或方框图看出系统的运算步骤和运算结构。 以后我们用流图来分析数字滤波器结构。DF网络结构或DF运算结构二个术语有微小的差别,但大抵一样,可以混用。

四、数字滤波器的分类 滤波器的种类很多,分类方法也不同。 1.从功能上分;低、带、高、带阻。 2.从实现方法上分:FIR、IIR 3.从设计方法上来分:Chebyshev(切比雪夫),Butterworth(巴特沃斯) 4.从处理信号分:经典滤波器、现代滤波器 等等。

1、经典滤波器 假定输入信号x(n)中的有用成分和希望去除的成分,各自占有不同的频带。当x(n)经过一个线性系统(即滤波器)后即可将欲去除的成分有效地去除。但如果信号和噪声的频谱相互重叠,那么经典滤波器将无能为力。 |H(ejw)| |Y(ejw)| |X(ejw)| 无用 有用 wc w wc wc w

2.现代滤波器 它主要研究内容是从含有噪声的数据记录(又称时间序列)中估计出信号的某些特征或信号本身。一旦信号被估计出,那么估计出的信号将比原信号会有高的信噪比。 现代滤波器把信号和噪声都视为随机信号,利用它们的统计特征(如自相关函数、功率谱等)导出一套最佳估值算法,然后用硬件或软件予以实现。 现代滤波器理论源于维纳在40年代及其以后的工作,这一类滤波器的代表为:维纳滤波器,此外,还有卡尔曼滤波器、线性预测器、自适应滤波器。 本课程主要讲经典滤波器,外带一点自适应滤波器.

3.模拟滤波器和数字滤波器 经典滤波器从功能上分又可分为: 低通滤波器(LPAF/LPDF):Low pass analog filter 带通滤波器(BPAF/BPDF):Bandpass analog filter 高通滤波器(HPAF/HPDF):High pass analog filter 带阻滤波器(BSAF/BSDF):Bandstop analog filter 即它们每一种又可分为:数字(Digital)和模拟(Analog)滤波器。

4.模拟滤波器的理想幅频特性 LPAF HPAF BPAF BSAF

5.数字滤波器的理想幅频特性 LPDF HPDF BPDF BSDF ……. ……. ……. …….

五、研究DF实现结构意义 1.滤波器的基本特性(如有限长冲激响应FIR与无限长冲激响应IIR)决定了结构上有不同的特点。 2.不同结构所需的存储单元及乘法次数不同,前者影响复杂性,后者影响运算速度。 3.有限精度(有限字长)实现情况下,不同运算结构的误差及稳定性不同。 4.好的滤波器结构应该易于控制滤波器性能,适合于模块化实现,便于时分复用。

六、本章介绍主要的内容 1.介绍IIR滤波器实现的基本结构。 2.介绍FIR滤波器实现的基本结构。 3.介绍一种特殊的滤波器结构实现形式:格型滤波器结构.

第二节 IIR DF的基本结构

一、IIR DF特点 1.单位冲激响应h(n)是无限长的n→∞ 2.系统函数H(z)在有限长Z平面(0<|Z|<∞)有极点存在。 3.结构上存在输出到输入的反馈,也即结构上是递归型的。 4.因果稳定的IIR滤波器其全部极点一定在单位园内。

二、IIR DF基本结构 IIR DF类型有:直接型、级联型、并联型。 直接型结构:直接I型、直接II型(正准型、典范型)。

1、 IIR DF系统函数及差分方程 一个N阶IIR DF有理的系统函数可能表示为: 以下我们讨论M<=N情况。 则这一系统差分方程为:

2、直接I型 (1)直接I型流图 IIR DF的差分方程就代表了一种最直接的计算公式,用流图表现出来的实现结构即为直接I型结构(即由差分方程直接实现。) 方程看出:y(n)由两部分组成: 第一部分 是一个对输入x(n)的M节延时链结构。即每个延时抽头后加权相加,即是一个横向网络。 第二部分 是一个N节延时链结构网络。不过它是对y(n)延时,因而是个反馈网络。 x(n) b0 y(n) Z-1 b1 a1 Z-1 Z-1 b2 a2 Z-1 Z-1 bM a N-1 Z-1 aN Z-1

(2)结构的特点 此结构的特点为: (1)两个网络级联:第一个横向结构M节延时网络实现零点,第二个有反馈的N节延时网络实现极点。 (2)共需(N+M)级延时单元 (3)系数ai,bi不是直接决定单个零极点,因而不能很好地进行滤波器性能控制。 (4)极点对系数的变化过于灵敏,从而使系统频率响应对系统变化过于灵敏,也就是对有限精度(有限字长)运算过于灵敏,容易出现不稳定或产生较大误差。

3、直接II型(正准型/典范型) (1)直接II型原理 从上面直接型结构的两部分看成两个独立的网络(即两个子系统)。 原理:一个线性时不变系统,若交换其级联子系统的次序,系统函数不变。把此原理应用于直接I型结构。即: (1)交换两个级联网络的次序 (2)合并两个具有相同输入的延时支路。 得到另一种结构即直接II型。

(2)直接II型的结构流图过程1--对调 x(n) b0 y(n) x(n) b0 y(n) Z-1 Z-1 b1 a1 a1 b1 Z-1 bM a N-1 a N-1 bM Z-1 Z-1 Z-1 aN aN Z-1 Z-1 第一部分 第二部分 对调

(3)直接II型的结构流图过程2--合并 由于对调后前后两路都有一条内容完全相同的延时链,可以合并为一条即可。 b0 x(n) b0 y(n) y(n) Z-1 a1 b1 Z-1 Z-1 a1 b1 a2 Z-1 b2 b2 Z-1 a2 Z-1 Z-1 a N-1 Z-1 bM bM a N-1 Z-1 aN Z-1 Z-1 aN 这就是直接II型的结构流图。 合并

(4)直接II型特点 直接II型结构特点: (1)两个网络级联。 第一个有反馈的N节延时网络实现极点; 第二个横向结构M节延时网络实现零点。 (2)实现N阶滤波器(一般N>=M)只需N级延时单元,所需延时单元最少。故称典范型。 (3)同直接I型一样,具有直接型实现的一般缺点。

例子 已知IIR DF系统函数,画出直接I型、直接II型的结构流图。 解:为了得到直接I、II型结构,必须将H(z)代为Z-1的有理式; x(n) 8 y(n) x(n) 8 y(n) Z-1 5/4 -4 注意反馈部分系数符号 Z-1 -4 5/4 Z-1 Z-1 -3/4 11 Z-1 11 -3/4 Z-1 Z-1 -2 1/8 Z-1 -2 1/8 Z-1

作业

4、级联型结构 (1)系统函数因式分解 一个N阶系统函数可用它的零、极点来表示即系统函数的分子、分母进行因式分解:

(2)系统函数系数分析

(3)基本二阶节的级联结构

(4)滤波器的基本二阶节 所以,滤波器就可以用若干个二阶网络级联起来构成。这每一个二阶网络也称滤波器的基本二阶节(即滤波器的二阶节)。一个基本二阶节的系统函数的形式为: 一般用直接II型(正准型、典范型表示) y(n) x(n) a1i Z-1 β1i Z-1 β2i a2i

(5)用二阶节级联表示的滤波器系统 整个滤波器则是多个二阶节级联 y(n) x(n) a11 Z-1 a12 Z-1 a1M Z-1 β11 β12 β1M Z-1 Z-1 …... Z-1 β21 β22 β2M a21 a22 a2M

例子 设IIR数字滤波器系统函数为: x(n) y(n) Z-1 1 Z-1 1 1 1 Z-1 1 1

(6)级联结构的特点 从级联结构中看出: 它的每一个基本节只关系到滤波器的某一对极点和一对零点。 调整β1i,β2i,只单独调整滤波器第I对零点,而不影响其它零点。 同样,调整a1i,a2i,……只单独调整滤波器第I对极点,而不影响其它极点。 级联结构特点: (a)每个二阶节系数单独控制一对零点或一对极点,有利于控制频率响应。 (b)分子分母中二阶因子配合成基本二阶节的方式,以及各二阶节的排列次序不同。

作业

5、并联型 (1)系统函数的部分分式展开 将系统函数展成部分分式的形式:用并联的方式实现DF。 “相加”在电路中实现用并联。如果遇到某一系数为复数,那么一定有另一个为共轭复数,将它们合并为二阶实数的部分分式。

(2)基本二阶节的并联结构 A0 其实现结构为: A1 a1 y(n) x(n) AN1 . aN1 a11 . a21 a1N2 a2N2 Z-1 a1 y(n) x(n) AN1 . Z-1 aN1 β01 a11 Z-1 . a21 β11 β0N2 a1N2 a2N2 β1N2

(3)并联型基本二阶节结构 并联型的基本二阶节的形式: 其中:要求分子比分母小一阶 x(n) β0 y(n) a1 Z-1 β1 Z-1

(4)并联型特点 (1)可以单独调整极点位置,但不能象级联那样直接控制零点(因为只为各二阶节网络的零点,并非整个系统函数的零点)。 (2)其误差最小。因为并联型各基本节的误差互不影响,所以比级联误差还少。若某一支路a1误差为1%,但总系统的误差仍可达到少1%。(因为分成a1,a2…...支路). 注意:(1)为什么二阶节是最基本的?因为二阶节是实系数,而一阶节一般为复系数。 (2)统一用二阶节表示,保持结构上的一致性,有利于时分多路复用。 (3)级联结构与并联结构的基本二阶节是不同的。

(5)例子 其并联结构为: 1 6 1 Z-1 x(n) y(n) -6 1 Z-1 4 Z-1 1

作业

第三节 FIR DF的结构 (有限长冲激响应滤波器)

一、FIR DF的特点 (1)系统的单位冲激响应h(n)在有限个n值处不为零。即h(n)是个有限长序列。 (2)系统函数|H(z)|在|z|>0处收敛,极点全部在z=0处(即FIR一定为稳定系统) (3)结构上主要是非递归结构,没有输出到输入反馈。但有些结构中(例如频率抽样结构)也包含有反馈的递归部分。

二、FIR的系统函数及差分方程 长度为N的单位冲激响应h(n)的系统函数为:

三、FIR滤波器实现基本结构 1.FIR的横截型结构(直接型) 2. FIR的级联型结构 3.FIR的频率抽样型结构

1、FIR直接型结构 (卷积型、横截型) (1)流图 x(n) y(n) h(0) x(n) h(1) Z-1 Z-1 Z-1 Z-1 Z-1 倒下 h(2) h(N) Z-1 h(0) h(1) h(N-1) Z-1 h(N-1) y(n) h(N) Z-1

(2)框图 x(n) Z-1 Z-1 Z-1 Z-1 h(N-1) h(0) h(1) h(2) ……. y(n)

2、级联型结构 (1)流图 当需要控制滤波器的传输零点时,可将H(z)系统函数分解成二阶实系数因子的形成: 即可以由多个二阶节级联实现,每个二阶节用横截型结构实现。 β01 β02 y(n) β0N/21 x(n) Z-1 Z-1 Z-1 β11 β12 β1N/2 …... Z-1 Z-1 Z-1 β2N/2 β21 β22

(2)级联型结构特点 由于这种结构所需的系数比直接型多,所需乘法运算也比直接型多,很少用。 由于这种结构的每一节控制一对零点,因而只能在需要控制传输零点时用。

作业

3、频率抽样型结构 (1)频率抽样型结构的导入 若FIR DF 的冲激响应为有限长(N点)序列h(n),则有: Z变换 H(z) h(n) 内插 DFT 单位园上频响 H(k) H(ejw) 取主值序列 N等分抽样 所以,对h(n)可以利用DFT得到H(k),再利用内插公式: 来表示系统函数。

(2)频率抽样型滤波器结构 由: 得到FIR滤波器提供另一种结构:频率抽样型结构。它是由两部分级联而成。 其中:级联中的第一部分为梳状滤波器: 第二部分由N个谐振器组成的谐振柜。

(3)梳状滤波器 (a)零、极点特性 由 看出: 它是一个由N节延时单元所组成的梳状滤波器。它在单位园上有N个等分的零点、无极点。

(b)幅频特性及流图 频率响应为: 梳状滤波器信号流图: |H(ejw)| y(n) 幅频曲线: x(n) 1 …... …... w -Z-N w

(4)谐振器 谐振器:是一个阶网络。 Hk(z) H(k) Z-1 谐振器的零极点:此为一阶网络,有一极点:

(5)谐振柜 谐振柜:它是由N个谐振器并联而成的。 这个谐振柜的极点正好与梳状滤波器的一个零点(i=k)相抵消,从而使这个频率(w=2πk/N)上的频率响应等于H(k). 将两部分级联起来,得到频率抽样结构。

(6)频率抽样型结构流图 H(0) Z-1 x(n) y(n) H(1) Z-1 -Z-N H(2) . Z-1 H(N-1) Z-1

(7)频率抽样型结构特点 (1)它的系数H(k)直接就是滤波器在 处的频率响应。因此,控制滤波器的频率响应是很直接的。 (2)结构有两个主要缺点: (a)所有的相乘系数及H(k)都是复数,应将它们先化成二阶的实数,这样乘起来较复杂,增加乘法次数,存储量。 (b)所有谐振器的极点都是在单位园上,由 决定考虑到系数量化的影响,当系数量化时,极点会移动,有些极点就不能被梳状滤波器的零点所抵消。(零点由延时单元决定,不受量化的影响)系统就不稳定了。

(8)修正的频率抽样结构 (a)产生的原因 为了克服系数量化后可能不稳定的缺点,将频率抽样结构做一点修正。即将所有零极点都移到单位园内某一靠近单位园、半径为r(r≤1)的园上,同时梳状滤波器的零点也移到r园上。(即将频率采样由单位园移到修正半径r的园上)

(b)修正的频率抽样结构的系统函数 为了使系数是实数,可将共 轭根合并,这些共轭根在半径为r的圆周上以实轴成对称分布。

(c)修正的频率抽样结构的系统 极点分布 |z|=r N=8 N=7

(d)修正频率结构的复根部分: 第k和第N-k个谐振器合并为一个实系数的二阶网络 因为h(n)是实数,它的DFT也是圆周共轭对称的。 因此,可以将第k和第N-k个谐振器合并为一个二阶网络。

(e)有限Q的谐振器 第k和第N-k个谐振器合并为一个二阶网络的极点在单位园内,而不是在单位园上,因而从频率响应的几何解释可知,它相当于一个有限Q的谐振器。其谐振频率为:

(f)修正频率抽样结构的谐振器的实根部分 除了共轭复根外,还有实根。 当N=偶数时,有一对实根,它们分别为 两点。 r -r 当N=奇数时,只有一个实根z=r(k=0),即只有H0(z).

(g)修正频率抽样结构流图(N=偶数) r y(n) x(n) -r .

(h)修正频率抽样结构流图(N=奇数) r y(n) x(n) .

(i)修正频率抽样结构的特点 (1)结构有递归型部分谐振柜又有非递归部分--梳状滤波器。 (2)它的零、极点数目只取决于单位抽样响应的长度,因而单位冲激响应长度相同,利用同一梳状滤波器、同一结构而只有加权系数β0k,β1k,H(0),H(N/2)不同的谐振器,就能得到各种不同的滤波器 (3)其结构可以高度模块化,适用于时分复用。

(j) 频率抽样结构的应用范围 (1)如果多数频率特性的采样值H(k)为零,例:窄带低通情况下,这时谐振器中剩下少数几个所需要的谐振器,因而可以比直接型少用乘法器,但存储器还是比直接型多用一些。 (2)可以共同使用多个并列的滤波器。例:信号频谱分析中,要求同时将信号的各种频率分量分别滤出来,这时可采用频率采样结构的滤波器,大家共用一个梳状滤波器及谐振柜,只是将各谐振器的输出适当加权组合就能组成各所需的滤波器。这样结构具有很大的经济性。 (3)常用于窄带滤波,不适于宽带滤波。

作业

4.快速卷积结构 (1)原理 设FIR DF的单位冲激响应h(n)的非零值长度为M,输入x(n)的非零值长度为N。 则输出y(n)=x(n)*h(n),且长度L=N+M-1 若将x(n)补零加长至L,补L-N个零点,将h(n)补零加长至L,补L-M个零点。 这样进行L点圆周卷积,可代替x(n)*h(n)线卷积。 其中: 而由圆卷积可用DFT和IDFT来计算,即可得到FIR的快速卷积结构。

(2)快速 卷积结构框图 x(n) X(k) L点DFT Y(k) L点IDFT h(n) L点DFT H(k) 当N,M中够大时,比直接计算线性卷积快多了。

5、线性相位FIR型结构 (1)定义 所谓线性相位:是指滤波器产生的相移与输入信号频率成线性关系。

(2)线性相位FIR DF具有特性 h(n)是因果的,为实数,且满足对称性。即满足约束条件: h(n)=±h(N-1-n)

(3)h(n)为偶、奇对称,N=偶数时 (a)FIR的线性相位的特性 令n’=N-1-n 代入 用n=n’ 应用线性FIR特性: h(n)=h(N-1-n)

(b) 线性相位FIR的结构流图 x(n) x(n-N/2+1) ……. h(0) h(1) h(3) h(N/2-1) h(2) ……. Z-1 Z-1 Z-1 Z-1 x(n-N/2+1) ……. Z-1 Z-1 Z-1 Z-1 Z-1 Z-1 h(0) h(1) h(3) h(N/2-1) h(2) ……. h(N-1) y(n) 其中h(0)=h(N-1),h(2)=h(N-2)……

(4)h(n)为奇、偶对称,N=奇数时 (a)FIR的线性相位的特性 当N=奇数时,有一中间项h((N-1)/2)无法合并,需提出:

(b) 线性相位FIR的结构流图 共有(N-3)/2项 x(n) h(0) h(1) h(3) h(2) ……. h(N-1) y(n) Z-1 Z-1 Z-1 Z-1 Z-1 Z-1 Z-1 Z-1 Z-1 Z-1 Z-1 h(0) h(1) h(3) h(2) ……. h(N-1) y(n) 其中h(0)=h(N-1),h(2)=h(N-2)……,h((N-3)/2)=h((N-1)/2

(5)总结:h(n)为偶对称,N=奇、偶数时FIR的线性相位的特性 同理,当h(n)=偶对称时,即h(n)=h(N-1-n),可求出: N=奇数时,

(6)h(n)为奇对称,N=奇、偶数时FIR的线性相位的特性 同理,当h(n)=奇对称时,即h(n)=-h(N-1-n),可求出: N=奇数时,

作业

第四节 格型滤波器

引言 在数字信号处理中,格型(Lattice)网络起着重要的作用。事实证明: (1)由于它的模块化结构便于实现高速并行处理; (2)一个m阶格型滤波器可以产生从1阶到m阶的m个横向滤波器的输出性能; (3)它对有限字长的舍入误差不灵敏。 由于这些优点,使得它在现代谱估计、语音处理、自适应滤波、线性预测和逆滤波等方面已得到广泛应用。

本节讨论: 1.全零点(FIR)格型滤波器 2.全极点(IIR)格型滤波器 3.零、极点(IIR)的格型滤波器

一、全零点(FIR)格型滤波器 一个M阶的FIR滤波器的系统函数H(z)可写成如下形式:

1.全零点格型滤波器网络结构

2.导出格型结构的参量 从上看出,要分析这一格型结构,先讨论如何由横向结构的参量导出格型结构的参量。或由格型结构的参量如何导出横向结构的参量。

在FIR横向结构中有M个 ,共需M次乘法,M次延迟; 在FIR的格型结构中也有M个参数ki(i=1,2,…M), ki称为反射系数,共需2M次乘法,M次延迟。 此格型结构的信号只有正馈通路,没有反馈通路,所以是一个典型的FIR系统。

3、格型网络单元 由上结构可看出:它们是由M个格型网络单元级联而成。每个网络单元有两个输入端和两个输出端,输入信号x(n)同时送到第一级网络单元的两个输入端,而在输出端仅取最后一级网络单元上面的一个输出端作为整个格型滤波器的输出信号y(n).

4、推导出格型结构网络系数{ki}的递推公式 如上图所示的基本格型单元的输入,输出关系如下式: 且 式中:fm(n)、gm(n)分别为第m个基本单元的上、下端的输出序列; fm-1(n)、 gm-1(n)分别为该单元的上、下端的输入序列;

设Bm(z)、Jm(z)分别表示由输入端x(n)至第m个基本单元的上、下端的输出端 fm(n)、 gm(n)对应的系统函数,即: 对上式分别除以F0(z)和G0(z)再代入Bm(z)、 Jm(z)式,得:

反过来 以上两式给出了格型结构中由低阶到高阶(或由高阶到低阶)系统函数的递推关系。

由于上式中同时包含B(z)和J(z)。实际中只给出Bm(z),所以应找出Bm(z)和 Bm-1(z)之间的递推关系。

5、导出km与滤波器系数bm之间的递推关系

6、实际中具体递推步骤 实际工作中,一般先给出H(z)=B(z)=BM(z),要画出H(z)的格型结构,需求出k1,k2,…kM。

7、例子

H(z)格型结构流图如图所示:

作业

二、全极点(IIR)格型滤波器 IIR滤波器的格型结构受限于全极点系统函数,可以根据FIR格型结构开发。设一个全极点系统函数由下式给定:

1、全极点格型网络单元 全极点IIR系统格型结构的基本单元为: 全极点IIR格型结构 基本单元 全零点FIR格型结构 基本单元

全极点(IIR)滤波器格型结构:

例子

IIR格型结构: 最后说明:一般的IIR滤波器既包含零点,又包含极点,它可用全极点格型作为基本构造模块,用所谓的格型梯形结构实现。

作业

三、零、极点系统(IIR系统)的格型结构 一个在有限z平面(0<|z|<)既有极点又有零点的IIR系统的系统函数H(z)可表示为

系统的格型结构流图

看出: (1)若k1=k2=…=kN=0,即所有乘k(或-k)处的联线全断开,则上图将变成一个N阶的FIR系统的横向结构; (2)若c1=c2=…=cN=0,即含c1…CN的联线都断开,C0=1那么上图将变成全极点IIR格型滤波器结构 (3)因此,图上半部分对应于全极点系统;下半部分对应于全零点系统B(z),且下半部分无任何反馈,故参数k1,k2,…,kN,仍可按全极点系统的方法求出,但上半部分对下半部分有影响,所以这里的Ci和全零点系统的bi不会相同。 任务:想办法求出各个ci,i=0,1,…,N。

推导 由上式,设 是由g0(n)至gm(n)之间的系统函数 是由x(n)至gm(n)之间的系统函数 又因为

则 整个系统的系统函数 应是 分别用 加权后的相加(并联)

求解参数c1,c2,…cN的方法 第一种方法: 以N=2为例,则有 其中

令等式两边的同次幂的系数相等,可得 由上式,可求得c2,c1,c0。一般情况下(任意N时,有)

第二种方法: 又定义 则有 对一般项m来说,可有 由此得出: 起始条件为: 已给定

例子 已知

解: (1)极点部分由全极点模型: 按求全极点模型的方法求出,实际上是用求全零点型的方法求得。即:

由H(z)可知: 利用式子 可求得:

其格型流图结构为:

总结本章主要的内容 1. IIR 滤波器实现的基本结构 2. FIR滤波器实现的基本结构 3. 一种特殊的滤波器结构实现形式:格型滤波器结构.

1.IIR DF基本结构 IIR DF类型有: 直接型直接型结构:直接I型、直接II型(正准型、典范型) 级联型 并联型

直接I型 直接I型流图 IIR DF的差分方程就代表了一种最直接的计算公式,用流图表现出来的实现结构即为直接I型结构(即由差分方程直接实现。) 方程看出:y(n)由两部分组成: 第一部分 是一个对输入x(n)的M节延时链结构。即每个延时抽头后加权相加,即是一个横向网络。 第二部分 是一个N节延时链结构网络。不过它是对y(n)延时,因而是个反馈网络。 x(n) b0 y(n) Z-1 b1 a1 Z-1 Z-1 b2 a2 Z-1 Z-1 bM a N-1 Z-1 aN Z-1

直接II型的结构流图 由于对调后前后两路都有一条内容完全相同的延时链,可以合并为一条即可。 b0 x(n) y(n) Z-1 a1 b1 bM a N-1 Z-1 Z-1 aN 这就是直接II型的结构流图。

级联型

级联型的基本二阶节 所以,滤波器就可以用若干个二阶网络级联起来构成。这每一个二阶网络也称滤波器的基本二阶节(即滤波器的二阶节)。一个基本二阶节的系统函数的形式为: 一般用直接II型(正准型、典范型表示) y(n) x(n) a1i Z-1 β1i Z-1 β2i a2i

级联型二阶节表示的滤波器系统 整个滤波器则是多个二阶节级联 y(n) x(n) a11 Z-1 a12 Z-1 a1M Z-1 β11 β12 β1M Z-1 Z-1 …... Z-1 β21 β22 β2M a21 a22 a2M

并联型 将系统函数展成部分分式的形式:用并联的方式实现DF。 “相加”在电路中实现用并联。如果遇到某一系数为复数,那么一定有另一个为共轭复数,将它们合并为二阶实数的部分分式。

并联型基本二阶节结构 并联型的基本二阶节的形式: 其中:要求分子比分母小一阶 x(n) β0 y(n) a1 Z-1 β1 Z-1 a2

二、FIR滤波器 长度为N的单位冲激响应h(n)的系统函数为:

FIR滤波器实现基本结构 (1)FIR的横截型结构(直接型) (2) FIR的级联型结构 (3)FIR的线性型 结构

1.FIR直接型结构 (卷积型、横截型) x(n) y(n) h(0) x(n) h(1) 倒下 h(2) h(N) h(N-1) y(n) Z-1 Z-1 Z-1 Z-1 Z-1 倒下 h(2) h(N) Z-1 h(0) h(1) h(N-1) Z-1 h(N-1) y(n) h(N) Z-1

2.级联型结构 当需要控制滤波器的传输零点时,可将H(z)系统函数分解成二阶实系数因子的形成: 即可以由多个二阶节级联实现,每个二阶节用横截型结构实现。 β01 β02 y(n) β0N/21 x(n) Z-1 Z-1 Z-1 β11 β12 β1N/2 …... Z-1 Z-1 Z-1 β2N/2 β21 β22

3.线性相位FIR型结构 所谓线性相位:是指滤波器产生的相移与输入信号频率成线性关系。

h(n)为偶数,N=奇、偶数时FIR的线性相位的特性 同理,当h(n)=偶对称时,即h(n)=h(N-1-n),可求出: N=奇数时,

h(n)为奇数,N=奇、偶数时FIR的线性相位的特性 当h(n)=奇对称时,即h(n)=-h(N-1-n),可求出: N=奇数时,

4.快速卷积结构 设FIR DF的单位冲激响应h(n)的非零值长度为M,输入x(n)的非零值长度为N。 则输出y(n)=x(n)*h(n),且长度L=N+M-1 若将x(n)补零加长至L,补L-N个零点,将h(n)补零加长至L,补L-M个零点。 这样进行L点圆周卷积,可代替x(n)*h(n)线卷积。 其中: 而由圆卷积可用DFT和IDFT来计算,即可得到FIR的快速卷积结构。

(2)快速 卷积结构框图 x(n) X(k) L点DFT Y(k) L点DFT h(n) L点DFT H(k) 当N,M中够大时,比直接计算线性卷积快多了。

5、频率抽样型结构 若FIR DF 的冲激响应为有限长(N点)序列h(n),则有: Z变换 H(z) h(n) 内插 DFT 单位园上频响 H(k) H(ejw) 取主值序列 N等分抽样 所以,对h(n)可以利用DFT得到H(k),再利用内插公式: 来表示系统函数。

(3)梳状滤波器 (a)零、极点特性 由 看出: 它是一个由N节延时单元所组成的梳状滤波器。它在单位园上有N个等分的零点、无极点。

(6)频率抽样型结构流图 H(0) W-k H(1) W-k H(2) W-k H(N-1) W-k Z-1 x(n) y(n) Z-1 -Z-N H(2) W-k Z-1 H(N-1) W-k Z-1

一、全零点(IIR)格型滤波器 一个M阶的FIR滤波器的系统函数H(z)可写成如下形式:

全零点格型滤波器网络结构

二、全极点(IIR)格型滤波器 IIR滤波器的格型结构受限于全极点系统函数,可以根据FIR格型结构开发。设一个全极点系统函数由下式给定:

全极点(IIR)滤波器格型结构: