数字信号处理主要知识点整理复习总结课件

,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,数字信号处理课程 知识点概要,数字信号处理课程 知识点概要,1,第1章 数字信号处理概念知识点,1、掌握连续信号、模拟信号、离散时间信号、数字信号的特点及相互关系（时间和幅度的连续性考量）,2、数字信号的产生；,3、典型数字信号处理系统的主要构成。,量化、编码,采样,模拟信号,离散时间信号,数字信号,第1章 数字信号处理概念知识点1、掌握连续信号、模拟信号、离,2,A / D,变换器,通用或,专用,计算机,采样,保持器,D/ A,变换器,模拟,低通,滤波器,模拟,信号,数字信号,模拟,信号,连续时间,信号,连续时间,信号,数字信号处理系统,A / D通用或采样D,3,1. 周期序列的判断与周期T的求取。,基本概念题（填空、判断、选择）。,本章典型题型与习题讲解,：,2. 判断系统是否是线性非时变系统。,Linear system : 齐次性与叠加性,即,y,1,(n)=T,x,1,(n) ，,y,2,(n)=T,x,2,(n),y,(n)=Ta,x,1,(n) b,x,2,(n) = a,y,1,(n) b,y,2,(n),*加权信号和的响应=响应的加权和。,Time-invariant: 时不变特性,即,y,(n-n,0,)=T,x,(n-n,0,),1. 周期序列的判断与周期T的求取。基本概念题（填空、判断,4,习题1. 判断下面的序列是否是周期的，若是周期的，确定其周期。,（1）,（3）,解,：,（1）,（2）,这是无理数，因此是非周期序列。,A是常数；,这是有理数，因此是周期序列，周期是T=14；,习题1. 判断下面的序列是否是周期的，若是周期的，确定其周期,5,4. 线性卷积的计算。,5.模拟信号数字处理的方法与过程；采样、恢复的概念；采样定理及采样后产生的影响；预滤波、平滑滤波的作用；,4. 线性卷积的计算。5.模拟信号数字处理的方法与过程；采样,6,数字信号处理主要知识点整理复习总结课件,7,数字信号处理主要知识点整理复习总结课件,8,数字信号处理主要知识点整理复习总结课件,9,第二部分 离散时间系统,1、线性时不变系统的判定,2、线性卷积,3、系统稳定性与因果性的判定,4、线性时不变离散时间系统的表示方法,5、,系统分类及两种分类之间的关系,第二部分 离散时间系统,10,1、线性系统：对于任何线性组合信号的,响应等于系统对各个分量的响应的线性组合。,线性系统,判别准则,若,则,2、时不变系统：系统的参数不随时间而变化，不管输入信号作用时间的先后，输出信号的响应的形状均相同，仅是出现时间的不同,若,则,时不变系统,判别准则,1、线性系统：对于任何线性组合信号的响应等于系统对各个分量的,11,3、线性卷积,y(n)的长度L,x,L,h,1,两个序列中只要有一个是无限长序列，则卷积之后是无限长序列,卷积是线性运算，长序列可以分成短序列再进行卷积，但必须看清起点在哪里,3、线性卷积y(n)的长度LxLh1,12,系统,时域充要条件,Z域充要条件,因果,h(n)0 (n0),ROC: R1 ,Z,稳定,h(n),n=-,ROC: 包含单位圆,4、系统的稳定性与因果性,系统时域充要条件Z域充要条件因果h(n)0 (n0)RO,13,5、,差分方程描述系统输入输出之间的运算关系,N阶线性常系数差分方程的一般形式：,其中,a,i,、,b,i,都是常数。,离散系统差分方程表示法有两个主要用途：, 求解系统的瞬态响应；,由差分方程得到系统结构；,5、差分方程描述系统输入输出之间的运算关系,14,6、线性时不变离散时间系统的表示方法,线性常系数差分方程,单位脉冲响应 h(n),系统函数 H(z),频率响应 H(e,jw,),零极点图（几何方法）,7、系统的分类,IIR和FIR,递归和非递归,6、线性时不变离散时间系统的表示方法,15,例1. 判断下列系统是否为线性系统。,解,：(a),故为线性系统。,例1. 判断下列系统是否为线性系统。解：(a)故为线性系统。,16,（b),故为线性系统。,（b)故为线性系统。,17,故不是线性系统。,(c),可见：,故不是线性系统。(c)可见：,18,(d),故不是线性系统。,可见：,(d)故不是线性系统。可见：,19,例2 判断系统 是否是移不变系统。,其中a和b均为常数,解：,故为移不变系统。,例2 判断系统,20,例3 判断系统 是否是移不变系统。,解：,故不是移不变系统。,又：,显然,例3 判断系统,21,例4.,判断下列系统是否为移不变系统。,解：,故不是移不变系统。,又：,显然,（a）,例4. 判断下列系统是否为移不变系统。解：故不是移不变系统。,22,故是移不变系统。,又：,显然,（b),故是移不变系统。又：显然（b),23,一个常系数线性差分方程是否表征一个线性移不变系统，这,完全由边界条件决定。,例如：差分方程,（c) 边界条件 时，既不是线性的也不是移不变的。,（a) 边界条件 时，是线性的但不是移不变的。,（b) 边界条件 时，是线性移不变的。,一个常系数线性差分方程是否表征一个线性移不变系统，这（c),24,令,.,所以：,令.所以：,25,.,所以：,可见 是移一位的关系， 亦是移一位的关系。因此是移不变系统。,.所以：可见,26,代入差分方程，得：,代入差分方程，得：,27,.,所以：,因此为线性系统。,.所以：因此为线性系统。,28,3. 判断系统是否是因果稳定系统。,Causal and Noncausal System（因果系统）,causal system: (1) 响应不出现于激励之前,(2) h(n)=0, n0 （线性、时不变系统）,Stable System （稳定系统）,(1) 有界输入导致有界输出,(2) （线性、时不变系统）,(3) H(z)的极点均位于Z平面单位圆内（因果系统）,*实际系统一般是因果系统；,* y(n)=,x,(-n)是非因果系统,因n0时的输入;,3. 判断系统是否是因果稳定系统。 Causal and,29,（b）由于 领先于 ，故为非因果系统。,例5 判断下列系统是否为因果系统。,（a) 为因果系统，由定义可知。,解：,（b）由于 领先于,30,由于 由目前和过去的输入所决定，故为,因果系统。,由于,n,=-1时，有y（-1）=,x,(1)；,也就是 领先于 ，故为非因果系统。,由于 由目前和过去的输入所,31,第2章回顾要点与难点,1、Z变换,Z变换的定义、零极点、收敛域,逆Z变换（部分分式法）,Z变换的性质及Parseval定理,2、离散时间傅里叶变换,DTFT的定义、性质,DTFT与Z变换的关系,DTFT存在的条件,3、DFT,DFT定义，与Z变换的关系，DFT性质,4、FFT,5、DFT的应用,第2章回顾要点与难点,32,2.1节知识点,1、DTFT的定义：,正变换：,反变换：,基本性质。,常见变换对；,离散时间信号的频域（频谱）为周期函数；,2.1节知识点 正变换：反变换：基本性质。常见变换对；离散,33,Condition:,(DTFT)序列傅立叶变换,(IDTFT)序列傅立叶反变换,注,：,周期序列不满足该绝对可和的条件，因此它的DTFT不存在。,1. DTFT的计算及其性质。,方法1：根据定义式求解,Condition:(DTFT)序列傅立叶变换(IDTFT),34,一般序列,共轭对称序列,共轭反对称序列,一般实序列,偶序列,奇序列,方法2：根据DTFT的性质求解（特别是对称性）,一般序列共轭对称序列共轭反对称序列一般实序列偶序列奇序列方法,35,（a）序列分成实部与虚部时:,其中,序列分成实部与虚部两部分，实部对应的FT具有共轭对称性，虚部和j一起对应的FT具有共轭反对称性。,（a）序列分成实部与虚部时:其中 序列分成实,36,其中,（b）序列分成共轭对称 与共轭反对称 时:,序列的共轭对称部分,x,e,(n)对应着FT的实部,X,R,(e,j,)，而序列的共轭反对称部分x,o,(n)对应着FT的虚部j,X,I,(e,j,),。,其中（b）序列分成共轭对称 与共轭反对称,37,例1：若序列h(n)是实因果序列，其DTFT的实部如下式： H,R,(e,j,)1+cos,求序列h(n)及其傅里叶变换H(e,j,).,解：,例1：若序列h(n)是实因果序列，其DTFT的实部如下式：,38,数字信号处理主要知识点整理复习总结课件,39,数字信号处理主要知识点整理复习总结课件,40,2、Z 变换表示法：,1) 级数形式（定义）,2) 解析表达式,（根据常见公式）,（注意:表示收敛域上的函数，同时注明收敛域）,3、Z 变换收敛域的特点：,1) 收敛域是一个圆环，有时可向内收缩到原点，有时可向外扩展到，只有,x,(n)=(n)的收敛域是整个Z 平面,2) 在收敛域内没有极点，,X,(z)在收敛域内每一点上都是解析函数。,2、Z 变换表示法：,41,4、几类序列Z变换的收敛域,(1) 有限长序列:X(z)=, x(n)z,-n, (n,1, n n,2,), 0, n,1, n n,2,0|z|,展开式出现z的负幂,n,1, n n,2, 0,0,|z|,展开式出现z的正幂,n,1, 0 0|z|,R,x,n,1, 0, n,2,= , R,x,|z|,展开式出现z的正幂,Z 变换的收敛域包括 点是因果序列的特征。,4、几类序列Z变换的收敛域 Z 变换的收敛域包括 点是因,42,(3) 左边序列,X(z)=, x(n)z,-n, (n,1, n n,2, n,1,=-),n,1,= -, n,2, 0, |z| 0, 0|z| R,x, R,x,|z| R,x,R,x, R,x, 空集,(3) 左边序列,43,5、部分分式法进行逆Z变换,求极点,将X(z)分解成部分分式形式,通过查表，对每个分式分别进行逆Z变换,注：左边序列、右边序列对应不同收敛域,将部分分式逆Z变换结果相加得到完整的x(n)序列,6、Z变换的性质,移位、反向、,乘指数序列、卷积,5、部分分式法进行逆Z变换,44,常用序列z变换（可直接使用）,常用序列z变换（可直接使用）,45,7、DTFT与Z变换的关系,采样序列在单位圆上的Z变换等于该序列的DTFT,序列频谱存在的条件Z变换的收敛域包含单位圆,8、Parseval定理重要应用计算序列能量：,即时域中对序列求能量与频域中求能量是一致,7、DTFT与Z变换的关系采样序列在单位圆上的Z变换等于该,46,分析计算题（计算证明、分析问答）。,本章典型题型与习题讲解,：,分析计算题（计算证明、分析问答）。 本章典型题型与习题讲解：,47,数字信号处理主要知识点整理复习总结课件,48,方法2. 幂级数法(长除法),左边序列：将X(z)的分子、分母按Z的升幂排列,右边序列：将X(z)的分子、分母按Z的降幂排列,对于大多数单阶极点的序列，常常用这种部分分式展开法求逆Z变换。,方法3. 部分分式展开法,3.逆Z变换的计算。,方法1. 用留数定理求逆Z变换,求逆z变换时特别需要注意收敛域的范围，收敛域不同，逆z变换的结果是不同的。如果没有明确告诉收敛域的范围，则求逆z变换时需要讨论。,方法2. 幂级数法(长除法)左边序列：将X(z)的分子、分,49,16. 已知:,求出对应,的各种可能的序列的表达式。,解：,有两个极点，因为收敛域总是以极点为界，因此收敛域有以下三种情况：,三种收敛域对应三种不同的原序列。,时，,（1）当收敛域,令,，因为c内无极点，x(n)=0；,，C内有极点0，但z=0是一个n阶极点，改为求圆外极点留数，圆外极点有,16. 已知:求出对应的各种可能的序列的表达式。解：时，（1,50,那么,（2）当收敛域,时，,，C内有极点0.5；,那么（2）当收敛域时，C内有极点0.5；,51,，C内有极点0.5，0，但0是一个n阶极点，改成求c外极点留数，c外极点只有一个，即2，,最后得到,（3）当收敛域,，C内有极点0.5，0，但0是一个,52,n0，由收敛域判断，这是一个因果序列，因此x(n)=0。,，C内有极点0.5，2；,或者这样分析，C内有极点0.5，2，0，但0是一个n阶极点，改成求c外极点留数，c外无极点，所以x(n)=0。,最后得到,nN，则L=M）。较短的一个需要补0至L（两个序列的长度要求相等）。循环卷积可以用DFT(FFT)实现；,用循环卷积实现线性卷积：L,M+N-1,若不满足这个条件，则只在N-1,n M-1范围内两者相等。,3.循环卷积的计算方法，循环卷积与线性卷积的关系，用DFT计,73,典型题型与习题讲解,：,分析计算题（计算证明、分析问答、判断）。,典型题型与习题讲解：分析计算题（计算证明、分析问答、判断）。,74,数字信号处理主要知识点整理复习总结课件,75,数字信号处理主要知识点整理复习总结课件,76,数字信号处理主要知识点整理复习总结课件,77,数字信号处理主要知识点整理复习总结课件,78,数字信号处理主要知识点整理复习总结课件,79,2.4频域采样定理,如果x（n）的长度为M，则,只有当频,域采样点数N,M时,，才有,可由频域采样 恢复原序列x（n），,否则将产生时域混叠现象。,在z平面的单位圆上的N个等角点上，对z变换进行取样，将导致相应的时间序列周期延拓，延拓周期为N。,2.4频域采样定理 在z平面的单位圆上的N个等角,80,DFS,DFT,线性,线性,序列移位,循环移位,共轭对称性,共轭对称性,周期卷积,循环卷积,DFSDFT线性线性序列移位循环移位共轭对称性共轭对称性周期,81,DFT,选频性,DFT与Z变换,DFT与DTFT,DFT形式下的,Parseval定理,DFT选频性DFT与Z变换DFT与DTFTDFT形式下的Pa,82,重新构造两个长度为L的序列x(n)和y(n), 方法：末尾补零,对x(n)和y(n)进行圆周卷积：,首先对两个序列进行周期延拓,对延拓后的周期序列进行周期卷积,对周期卷积的结果取主值区间,使圆周卷积等于线性卷积而不产生混淆的必要条件是LN+M-1；,步骤如下：,重新构造两个长度为L的序列x(n)和y(n), 方法：末尾补,83,圆 周 卷 积 与 线 性 卷 积 的 性 质 对 比,圆周卷积,线性卷积,针对FFT引出的,一种,表示方法,信号通过线性系统时，信号输出等于,输入与系统单位冲激响应的卷积,两序列长度必须,相等,，,不等时按要求,补足零值点,两序列长度可以,不等,如x1(n)为 N1点，,x2(n)为 N2点,卷积结果长度,与两信号长度相等皆为,N,卷积结果长度为,N=N1+N2-1,圆 周 卷 积 与 线 性 卷 积 的 性 质 对 比圆周卷,84,变量,周期,分辨率,数字频域,模拟频域,离散频域,变量周期分辨率数字频域模拟频域离散频域,85,时域/频域同时采样,对有限时宽的信号x,a,(t)的时域波形和频域波形同时进行取样，其结果是时域波形和频域的都变成了离散的、周期性的波形；,时域内的离散周期信号为 ，频域内离散周期信号为 ，它们之间形成DFS变换对；,分别取它们的一个周期，得到x(n)与X(k)，它们之间形成DFT变换对。,n,N,0,k,0,N,-N,1/T,-N,时域/频域同时采样nN0k0N-N1/T-N,86,第二部分 快速傅里叶变换FFT,1、FFT计算原理。,2、基2时间抽取算法和频率抽取算法。,3、DFT、R-2 FFT算法的运算量比较。,4、实数序列的FFT高效算法。,5、FFT的应用。,第二部分 快速傅里叶变换FFT1、FFT计算原理。,87,主要要求掌握的内容,：,1、FFT、IFFT的计算方法、特点，DIT、DIF的运算流图。,2、FFT应用于频谱分析和快速卷积。,3、DFT、FFT的运算量计算。,4、FFT减少运算量的途径。,本章典型题型与习题讲解,：,作图题（作图、计算）。,主要要求掌握的内容：1、FFT、IFFT的计算方法、特点，D,88,N,点的,FFT,的运算量为,复乘:,C,M,=（N/2）M=（N/,2）,log,2,N,复加:,C,A,=N M=N log,2,N,1. 画出N点（例如8点、16点）FFT的运算流图,2. FFT的特点，FFT减少运算量的途径。,DIT DIF,3. FFT的运算量的计算，与DFT运算量的比较。,FFT算法的基本思想、特点、编程方法,N,点的,DFT,的运算量为,复乘:,C,M,=N,2,复加:,C,A,=N（N-1）,N点的FFT的运算量为1. 画出N点（例如8点、16点）FF,89,例1：如果通用计算机的速度为平均每次复数乘需要5s,每次复数加需要1s,用来计算,N,1024点DFT，问直接计算需要多少时间。用FFT计算呢?照这样计算，用FFT进行快速卷积对信号进行处理时，估计可实现实时处理的信号最高频率。,解：N=1024=2,10,直接计算DFT的运算量：,复乘:,C,M,=N,2,1024,2,2,20,次,复加:,C,A,=N（N-1）102410231047552,直接计算DFT所用的时间为：,例1：如果通用计算机的速度为平均每次复数乘需要5s,每次复,90,用FFT计算DFT,的运算量为,复乘:,C,M,=（N/2）M=（N/,2）,log,2,N1024/2105120,复加:,C,A,=N M=N log,2,N10241010240,用FFT计算DFT所用的时间为：,快速卷积时，要计算一次N点FFT（H(k)已经计算好存入ROM中了，不需用FFT计算出H(k)）；N次频域复数乘法（H(k)*X(k)）；一次N点IFFT（也是用FFT实现的）。所以，计算1024点快速卷积的计算时间约为,用FFT计算DFT的运算量为用FFT计算DFT所用的时间为：,91,所以，每秒种处理的采样点数（即采样速率）为,.,3.实数序列的FFT高效算法。,由采样定理可知，可实时处理的信号最高频率为,实际实现时，fmax要比这个小一些。,所以，每秒种处理的采样点数（即采样速率）为.3.实数序列的F,92,3. 已知,和,是两个N点实序列,和,的DFT，若要从,和,求,和,，为提高运算效率，试设计用一次N点IFFT来完成。,解：因为,和,均为实序列，所以，,和,为共轭对称序列，j,为共轭反对称序列。可令,和j,分别作为复序列,分量和共轭反,对称分量，即,计算一次N点IFFT得到,3. 已知和是两个N点实序列和的DFT，若要从和求和，为提高,93,由DFT的共轭对称性可知，,故,由DFT的共轭对称性可知，故,94,2.6节知识点,连续信号的频谱分析 (利用DFT的选频性),过程：采样截短DFT,效应：混叠原因：采样、频谱泄漏,泄漏原因：截短,栅栏效应原因：DFT,DFT的分辨率,2.6节知识点连续信号的频谱分析 (利用DFT的选频性),95,DFT的应用（频谱分析、分段卷积）。,频谱分析：DFT代替频谱分析引起的误差（混叠现象、栅栏效应、截断效应频谱泄漏、谱间干扰）；提高谱分辨率的方法；分段卷积（重叠相加法、重叠保留法）,15. 用微处理机对实数序列作谱分析,要求谱分辨率,，信号最高频率为1kHZ，试确定以下各参数：,（3）最少采样点数,（4）在频带宽度不变的情况下，将频率分辨率提高一倍的N值。,；,（1）最小记录时间,（2）最大取样间隔,；,；,DFT的应用（频谱分析、分段卷积）。频谱分析：DFT代替频谱,96,第3章回顾要点与难点,(1)数字滤波器频响应能模仿模拟滤波器频响,(2)因果稳定的模拟系统变换为数字系统仍为因果稳定的,S到Z平面的映射关系满足条件,第3章回顾要点与难点(1)数字滤波器频响应能模仿模拟滤波,97,主要内容,：,1、 数字滤波器的分类及特性。,2、 数字信号系统的信号流图。,3、 IIR滤波器的结构和信号流图：直接型；级联型；并联型。,4、 FIR数字滤波器的结构和信号流图：直接型；快速卷积型、频率采样型。,3.1 数字滤波器的结构,主要内容：1、 数字滤波器的分类及特性。3.1 数字滤,98,本章主要要求掌握的内容,：,1、数字信号系统的信号流图描述方法。,2、 IIR滤波器的信号流图：直接型；级联型；并联型。,3、 FIR数字滤波器的实现流图：直接型；级联型；线性相位型。,1.画出滤波器的实现结构（实现流图）。,本章主要要求掌握的内容：1、数字信号系统的信号流图描述方法。,99,IIR数字滤波器的直接I型结构,IIR数字滤波器的直接I型结构,100,两条延时链中对应的延时单元内容完全相同,可合并,得,两条延时链中对应的延时单元内容完全相同,可合并,得,101,数字信号处理主要知识点整理复习总结课件,102,数字信号处理主要知识点整理复习总结课件,103,2、FIR数字滤波器：非递归结构，无反馈，但在频率采样结构等某些结构中也包含有反馈的递归部分。,（1）直接型（卷积型、横截型）,（2）级联型,（3）线性相位型,（4）频率采样型,2、FIR数字滤波器：非递归结构，无反馈，但在频率采样结构等,104,105,直接型的转置,：,直接型的转置：,106,数字信号处理主要知识点整理复习总结课件,107,FIR数字滤波器要点与难点,1、线性相位：系统的相频特性是频率的线性函数,群时延:,偶对称,奇对称,FIR数字滤波器要点与难点1、线性相位：系统的相频特性是,108,2、四种线性相位FIR滤波器,2、四种线性相位FIR滤波器,109,四种线性相位FIR DF特性,第一类 ，h(n)偶、N奇，四种滤波器都可设计。,第二类,， h(n)偶、 N偶，可设计低、带通滤波器,不能设计高通和带阻。,第三类,， h(n)奇、 N奇，只能设计带通滤波器，,其它滤波器都不能设计。,第四类,， h(n)奇、 N偶，可设计高通、带通滤波,器，不能设计低通和带阻。,四种线性相位FIR DF特性第二类 ， h(n)偶、 N偶，,110,小结,1、相位特性只取决于h(n)的对称性，而与h(n)的值无关。,2、幅度特性取决于h(n)。,3、设计FIR数字滤波器时，在保证h(n)对称的条件下，只要完成幅度特性的逼近即可。,注意：当H()用H()表示时，当H()为奇对称时，其相频特性中还应加一个固定相移,小结1、相位特性只取决于h(n)的对称性，而与h(n)的值无,111,3、线性相位FIR滤波器的零点特性,零点必须是互为倒数的共轭对,3、线性相位FIR滤波器的零点特性零点必须是互为倒数的共轭,112,数字信号处理主要知识点整理复习总结课件,113,数字信号处理主要知识点整理复习总结课件,114,数字信号处理主要知识点整理复习总结课件,115,作图题,典型题型与习题讲解,：,作图题典型题型与习题讲解：,116,1. 设系统用下面的差分方程描述：,试画出系统的直接型、级联型和并联型结构。,解：,将上式进行Z变换,1. 设系统用下面的差分方程描述：试画出系统的直接型、级联型,117,（,1,）按照系统函数,，画出直接型结构如,图（一）,所示。,（1）按照系统函数，画出直接型结构如图（一）所示。,118,（2）将,的分母进行因式分解,按照上式可以有两种级联型结构：,(a),(b),画出级联型结构如,图（二）（b）,所示,画出级联型结构如,图（二）（a）,所示,（2）将的分母进行因式分解 按照上式可以有两种级联型,119,级联型结构,图（二）（a）,级联型结构,图（二）（b）,级联型结构图（二）（a）级联型结构图（二）（b）,120,（3）将,进行部分分式展开,（3）将进行部分分式展开,121,根据上式画出并联型结构如图（三）所示。,根据上式画出并联型结构如图（三）所示。,122,第2部分要点与难点,(1)数字滤波器频响应能模仿模拟滤波器频响,(2)因果稳定的模拟系统变换为数字系统仍为因果稳定的,S到Z平面的映射关系满足条件,第2部分要点与难点(1)数字滤波器频响应能模仿模拟滤波器,123,主要内容,：,1、数字滤波器的设计方法：IIR的设计方法分类。,2、理想滤波器的特性及逼近方法：理想滤波器的特性；连续函数逼近方法。,3、模拟滤波器设计：几种逼近函数及特点；模拟滤波器逼近函数设计方法。,4、模拟滤波器的数字仿真：冲激响应不变法；双线性变换法。,5、数字滤波器的频率变换。,IIR数字滤波器的设计,主要内容：1、数字滤波器的设计方法：IIR的设计方法分类。,124,主要要求掌握的内容,：,1、数字滤波器的概念、技术指标、设计过程、设计方法。,2、,IIR数字滤波器的设计与模拟滤波器设计的关系；转换方法：冲激响应不变法；双线性变换法；,3、 Butterworth数字低通滤波器的设计。,4、 IIR数字滤波器频带变换方法（由低通，设计高通、带通、带阻滤波器）,5、IIR滤波器的特点。,综合设计题（计算）。,本章典型题型与习题讲解,：,主要要求掌握的内容：1、数字滤波器的概念、技术指标、设计过程,125,思路：,脉冲响应不变法,思路：脉冲响应不变法,126,脉冲响应不变法的,映射关系,脉冲响应不变法的映射关系,127,S 平面,Z 平面,脉冲响应不变法满足变换的映射条件，但映射关系,不是一一对应,的,。,S 平面Z 平面 脉冲响应不变法满足变换的映射条件，但映射,128,脉冲响应不变法,优点：,时域脉冲响应的模仿性能好,频率坐标的变换是线性的，与是线性关系。,脉冲响应不变法,缺点：,有频谱周期延拓效应. 只能用于带限的频响特性，如衰减特性很好的低通或带通；,脉冲响应不变法优点：,129,S,1,平面,Z,平面,S,平面,一一对应,双线性变换法,S1平面Z平面S平面一一对应双线性变换法,130,优点：S平面与Z平面是单值的一一对应关系,与成,非线性,关系,缺点,:,不会产生混叠现象；,映射关系,优点：S平面与Z平面是单值的一一对应关系与成非线性关系,131,畸变：,经,双线性变换,后，频率发生了非线性变化， 相应地，数字滤波器的幅频特性,在临界频率点会发生非线性变化。这种频率点的畸变可以通过,预畸,来加以校正。,注意：预畸不能在整个频率段消除非线性畸变，只能消除模拟和数字滤波器在特征频率点的畸变。,畸变：经双线性变换后，频率发生了非线性变化， 相应地，数字,132,设计步骤：,三：通过变量代换求H(z),设计步骤：三：通过变量代换求H(z),133,置换过程:,频响,：,置换过程:频响：,134,1. IIR滤波器的设计与实现。,1. IIR滤波器的设计与实现。,135,数字信号处理主要知识点整理复习总结课件,136,冲激不变法（或称为脉冲响应不变法）步骤：,（1）将模拟滤波器的传递函数Ha（s）展开成部分分式的形式：,（2）将由第（1）步所得到的s,k,代入到下式中：,（3）设一个T值，并将T值和ze,j,代入到上式中即可得到数字滤波器的频率响应。,T的选取应按照滤波器最高截止频率的2倍以上选取（T过大时，频率混叠现象严重。）,3. IIR模拟滤波器到数字滤波器的转换方法,冲激不变法（或称为脉冲响应不变法）步骤：（1）,137,（1）确定数字低通技术指标：,通带截止频率 、通带衰减 、,阻带截止频率 、阻带衰减 ；,（2）将数字低通指标转换成模拟低通指标：,（ 和 不变 ）,边界频率的变换关系：,频率预畸变,双线性变换法步骤：,（1）确定数字低通技术指标：频率预畸变双线性变换法步骤：,138,（3）设计模拟低通滤波器；,（4）转换成数字低通滤波器：,这里的采样间隔T可任意选取,通常取 T=1或T=2,（3）设计模拟低通滤波器；,139,4. IIR模拟滤波器到数字滤波器转换特性与对应关系,脉冲响应不变法,H,a,(s)的极点s,i,映射到z平面，其极点变为e,SiT,稳定条件：,产生频率混叠现象，不适合高通、带阻滤波器的设计。,（S Z）,4. IIR模拟滤波器到数字滤波器转换特性与对应关系脉冲响应,140,例：,.,设,h,(,t,)表示一模拟滤波器的单位冲激响应，,用脉冲响应不变法，将此模拟滤波器转换成数字滤波器（,h(n),表示单位取样响应，即,h(n),h,a,(nT),）。确定系统函数,H(z),，并把,T,作为参数，证明：,T,为任何值时，数字滤波器是稳定的，并说明数字滤波器近似为低通滤波器还是高通滤波器。,H,a,(s)的极点s,1,0.9，数字滤波器系统函数应为,例：.设h(t)表示一模拟滤波器的单位冲激响应，Ha(s),141,H（z)的极点为,画出T=0.5和T=1时的幅频响应，由图可以看出数字滤波器近似是低通滤波器。,H（z)的极点为画出T=0.5和T=1时的幅频响应，由图可以,142,（S Z）,双线性变换法,稳定条件：,消除了频率混叠，但产生了频率畸变现象，需要预畸变处理。,（S Z）双线性变换法稳定条件： 消除了,143,5. 已知模拟滤波器的传输函数为：,（2）,试用脉冲响应不变法和双线性变换法分别将其转换为数字滤波器，设T=2s。,（1）,解：,（1）用脉冲响应不变法,方法1 直接按脉冲响应不变法设计公式，,的极点为：,5. 已知模拟滤波器的传输函数为：（2）试用脉冲响应不变法和,144,代入T=2s,代入T=2s,145,方法2 直接套用4题(2)所得公式，为了套用公式，先对,为一常数，,的分母配方，将 化成4题中的标准形式：,由于,所以,方法2 直接套用4题(2)所得公式，为了套用公式，先对为,146,对比可知，,套用公式得,对比可知，套用公式得,147,或通分合并两项得,148,（2）用双线性变换法,（2）用双线性变换法,149,150,3.3.4节要点,1. 从模拟滤波器低通原型到各种数字滤波器的频率变换,了解设计IIR数字滤波器的两种变换法,其中第二种要求会低通变换和高通变换,2.从数字滤波器低通原型到各种数字滤波器的频率变换,已知 ，会利用表，求,3.3.4节要点1. 从模拟滤波器低通原型到各种数字滤波器的,151,主要内容,：,1、FIR滤波器的设计方法分类。,2、FIR滤波器的线性相位特性：线性相位特性；实现 FIR滤波器的线性相位特性的条件。,3、 FIR,滤波器的窗函数截取方法：理想滤波特性的傅立叶级数逼近；窗函数截取的吉布斯效应和解决方法；常用的窗函数。,4、 FIR 滤波器的窗函数设计法设计步骤。,5、 FIR滤波器的频率取样设计法。,第 3部分 FIR数字滤波器的设计,主要内容：1、FIR滤波器的设计方法分类。第 3部分 FI,152,第3部分FIR数字滤波器要点与难点,1、线性相位：系统的相频特性是频率的线性函数,群时延:,偶对称,奇对称,第3部分FIR数字滤波器要点与难点1、线性相位：系统的相,153,3.4.2节要点,1、窗口设计法步骤；,2 、线性相位理想低通FIR DF 的设计 (会求h(n)；,3 、窗口函数对理想特性的影响；,(过渡带, 肩峰, Gibbs效应, 窗函数的要求, 常用窗函数的名称 ),3.4.2节要点1、窗口设计法步骤；,154,4、窗口法设计原理：,卷积关系,4、窗口法设计原理：卷积关系,155,数字信号处理主要知识点整理复习总结课件,156,窗口函数对理想特性的影响：,改变了理想频响的边沿特性，形成过渡带，宽为,，等于W,R,()的主瓣宽度。（决定于窗长）,过渡带两旁产生肩峰和余振（带内、带外起伏），取决于 W,R,()的旁瓣，旁瓣多，余振多；旁瓣相对值大，肩峰强 ，与 N无关。（决定于窗口形状）,N增加,过渡带宽减小,肩峰值不变。当N增加时，幅值变大，频率轴变密，而最大肩峰永远为8.95%，这种现象称为吉布斯（Gibbs）效应。,窗口函数对理想特性的影响：,157,窗函数的要求：,窗谱主瓣宽度要窄，以获得较陡的过渡带；,相对于主瓣幅度，旁瓣要尽可能小，使能量尽量集中在主瓣中，这样就 可以减小肩峰和余振，以提高阻带衰减和通带平稳性。,但实际上这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑制。,肩峰值的大小决定了滤波器通带内的平稳程度和阻带内的衰减，所以对滤波器的性能有很大的影响。,窗函数的要求：窗谱主瓣宽度要窄，以获得较陡的过渡带；,158,3.4.3节要点与难点,基本思想,在某些离散频率点上的值准确地等于所需滤波器在这些频率点处的值，其它频率处的特性则有较好的逼近,线性相位FIR DF的约束条件,线性相位低通FIR DF设计,会求各采样点的H(k),增大阻带衰减的两种方法,内插公式,3.4.3节要点与难点基本思想在某些离散频率点上的值准确,159,增大阻带衰减的两种方法：,1）加宽过渡带宽，以牺牲过渡带换取阻带衰减的增加,2）如果要进一步增加阻带衰减，但又不允许再增加过渡带宽，可增加采样点数N。,增大阻带衰减的两种方法：,160,主要要求掌握的内容,：,1、 FIR滤波器的线性相位特性和实现条件。四种基本类型的FIR滤波器。,2、窗函数截取的吉布斯效应和解决方法。,3、各种窗函数； FIR 滤波器的窗函数设计法。,4、频率采样法设计FIR滤波器。,5、FIR与 IIR数字滤波器的比较。,典型题型与习题讲解,：,综合设计题（计算）。,主要要求掌握的内容：1、 FIR滤波器的线性相位特性和实现条,161,1. FIR滤波器的设计与实现。,1. FIR滤波器的设计与实现。,162,数字信号处理主要知识点整理复习总结课件,163,3. 设FIR滤波器的系统函数为,求出该滤波器的单位取样响应,，,判断是否具有线性相位，求出其幅度特性和相位特性，并画出其直接型结构和线性相位型结构和线性相位型结构。,解：对FIR数字滤波器，其系统函数为,所以，其单位脉冲响应为,3. 设FIR滤波器的系统函数为求出该滤波器的单位取样响应,164,由 的取值可知 满足,所以，该FIR滤波器具有第一类线性相位特性。设其频率响应函数为,由 的取值可知 满足,165,直接型结构,线性相位型结构,直接型结构线性相位型结构,166,幅度特性函数为,相位特性函数为,由,画出直接型结构和线性相位型结构分别如,图（一）,和,图（二）,所示。,幅度特性函数为相位特性函数为由 画出直接型,167,4. 用矩形窗设计线性相位低通滤波，逼近滤波器传输函数,为,（1）求出相应于理想低通的单位脉冲响应,（2）求出矩形窗设计 的表达式，确定a与N之间的关系；,（3）N取奇数或偶数时对滤波特性有什么影响？,4. 用矩形窗设计线性相位低通滤波，逼近滤波器传输函数为（1,168,解：（1）,（2）为了满足线性相位条件，要求,为矩形窗函数长度。加矩形窗函数得到,解：（1）（2）为了满足线性相位条件，要求为矩形窗函数长度。,169,（3）N取奇数时，幅度特性函数 关于 三点偶对称，可实现各类幅频特性；N取偶数时， 关于 奇对称，所以不能实现高通、带阻和点阻滤波特性。,（3）N取奇数时，幅度特性函数 关于,170,IIR与FIR数字滤器的比较,FIR,IIR,设计方法,一般无解析的设计公式，要借助计算机程序完成,利用AF的成果，可简单、有效地完成设计,设计结果,可得到幅频特性（可以多带）和线性相位（最大优点）,只能得到幅频特性，相频特性未知（一大缺点），如需要线性相位，须用全通网络校准，但增加滤波器阶数和复杂性,稳定性,极点全部在原点（永远稳定）无稳定性问题,有稳定性问题,阶数,高,结构,非递归,递归系统,运算误差,一般无反馈，运算误差小,有反馈，由于运算中的四舍五入会产生极限环,快速算法,可用FFT实现，减少运算量,无快速运算方法,低,IIR与FIR数字滤器的比较FIRIIR设计方法一般无解,171,