1、数字信号处理课程数字信号处理课程 知识点概要知识点概要 1ppt课件第第1章章 数字信号处理概念知识点数字信号处理概念知识点1、掌握连续信号、模拟信号、离散时间信号、数字、掌握连续信号、模拟信号、离散时间信号、数字信号的特点及相互关系（时间和幅度的连续性考量）信号的特点及相互关系（时间和幅度的连续性考量）2、数字信号的产生；、数字信号的产生；3、典型数字信号处理系统的主要构成。、典型数字信号处理系统的主要构成。量化、编码量化、编码采样采样模拟信号模拟信号离散时间信号离散时间信号数字信号数字信号2ppt课件 A/D变换器变换器通用或通用或专用专用计算机计算机采样采样保持器保持器D/A变换器变换器
2、模拟模拟低通低通滤波器滤波器模拟模拟信号信号数字信号数字信号模拟模拟信号信号连续时间连续时间信号信号连续时间连续时间信号信号数字信号处理系统数字信号处理系统3ppt课件1.周期序列的判断与周期周期序列的判断与周期T的求取。的求取。基本概念题（填空、判断、选择）。基本概念题（填空、判断、选择）。本章典型题型与习题讲解本章典型题型与习题讲解：02 判 断是 否 为 有 理 数。2.判断系统是否是线性非时变系统。判断系统是否是线性非时变系统。Linear system:齐次性与叠加性齐次性与叠加性即即 y1(n)=Tx1(n)，y2(n)=Tx2(n)y(n)=Tax1(n)bx2(n)=ay1(n
3、)by2(n)*加权信号和的响应=响应的加权和。Time-invariant:时不变特性时不变特性 即即 y(n-n0)=Tx(n-n0)4ppt课件3()cos()78x nAn1()8()jnx ne3214,73ww1 2,168ww习题1.判断下面的序列是否是周期的，若是周期的，确定其周期。（1）（3）解解：（1）（2）这是无理数，因此是非周期序列。A是常数；这是有理数，因此是周期序列，周期是T=14；5ppt课件4.线性卷积的计算。线性卷积的计算。5.模拟信号数字处理的方法与过程；采样、恢模拟信号数字处理的方法与过程；采样、恢复的概念；采样定理及采样后产生的影响；预复的概念；采样定理
4、及采样后产生的影响；预滤波、平滑滤波的作用；滤波、平滑滤波的作用；6ppt课件7ppt课件8ppt课件9ppt课件 第二部分第二部分 离散时间系统离散时间系统 1、线性时不变系统的判定、线性时不变系统的判定 2、线性卷积、线性卷积 3、系统稳定性与因果性的判定、系统稳定性与因果性的判定 4、线性时不变离散时间系统的表示方法、线性时不变离散时间系统的表示方法 5、系统分类及两种分类之间的关系系统分类及两种分类之间的关系10ppt课件1、线性系统：对于任何线性组合信号的、线性系统：对于任何线性组合信号的响应等于响应等于系统对各个分量的响应的线性组合。系统对各个分量的响应的线性组合。线性系统线性系统
5、判别准则判别准则 若若11()()y nT x n 22()()y nT x n 1212()()()()T ax nbx nay nby n 则则2、时不变系统：系统的参数不随时间而变化，不管、时不变系统：系统的参数不随时间而变化，不管输入信号作用时间的先后，输出信号的响应的形状均输入信号作用时间的先后，输出信号的响应的形状均相同，仅是出现时间的不同相同，仅是出现时间的不同若若 ()()y nT x n 则则00()()T x nny nn时不变系统时不变系统判别准则判别准则 11ppt课件()()()()*()()()()*()kky nx k h nkx nh nx nk h kh nx
6、 n 3、线性卷积、线性卷积 y(n)的长度的长度LxLh1 两个序列中只要有一个是无限长序列，则卷两个序列中只要有一个是无限长序列，则卷积之后是无限长序列积之后是无限长序列 卷积是线性运算，长序列可以分成短序列再卷积是线性运算，长序列可以分成短序列再进行卷积，但必须看清起点在哪里进行卷积，但必须看清起点在哪里12ppt课件系统系统时域充要条件时域充要条件Z域充要条件域充要条件因果因果h(n)0(n0)ROC:R1 Z稳定稳定 h(n)n=-ROC:包含单位圆包含单位圆4、系统的稳定性与因果性、系统的稳定性与因果性13ppt课件5、差分方程差分方程描述系统输入输出之间的运算关系描述系统输入输出
7、之间的运算关系 N阶线性常系数差分方程的一般形式：阶线性常系数差分方程的一般形式：其中其中 ai、bi都是常数。都是常数。离散系统差分方程表示法有两个主要用途：离散系统差分方程表示法有两个主要用途：求解系统的瞬态响应；求解系统的瞬态响应；由差分方程得到系统结构；由差分方程得到系统结构；01()()()MNiiiiy na x nib y ni 14ppt课件6、线性时不变离散时间系统的表示方法、线性时不变离散时间系统的表示方法线性常系数差分方程线性常系数差分方程单位脉冲响应单位脉冲响应 h(n)系统函数系统函数 H(z)频率响应频率响应 H(ejw)零极点图（几何方法）零极点图（几何方法）7、
8、系统的分类、系统的分类IIR和和FIR递归和非递归递归和非递归15ppt课件例1.判断下列系统是否为线性系统。5)(3)()();()()();()()();()()(22nxnydnxnycnxnybnnxnya解：(a)()(),()()()()(222111nxTnnxnynxTnnxnynnxny）)()()()()()()()(2211221122112211nxTanxTanyanyannxannxanxanxaT故为线性系统。16ppt课件（b)()(),()()()()(222212112nxTnxnynxTnxnynxny）)()()()()()()()(2211221122
9、22112211nxTanxTanyanyanxanxanxanxaT故为线性系统。17ppt课件)()(),()()()()(222212112nxTnxnynxTnxnynxny）)()(2)()()()()()(212122222121222112211nxnxaanxanxanxanxanxanxaT故不是线性系统。(c)()()()(22112211nxTanxTanxanxaT可见：)()()()(2222112211nxanxanxTanxTa18ppt课件(d)。加即，系统操作为乘）53)(5)(3(),(5)(3)(5)(3)(222111nxTnxnynxTnxnynxny
10、5)()(3)()(22112211nxanxanxanxaT故不是线性系统。)()()()(22112211nxTanxTanxanxaT22211122115)(35)(3)()(anxaanxanxTanxTa可见：19ppt课件例2 判断系统 是否是移不变系统。其中a和b均为常数bnaxny)()(解：)()()()()()(mnybmnaxmnxTnybnaxnxT故为移不变系统。20ppt课件例3 判断系统 是否是移不变系统。()()sin(20.1)y nx nn解：)1.02sin()()()()1.02sin()()(系统操作nmnxmnxTnynnxnxT故不是移不变系统。
11、又：1.0)(2sin)()(函数操作mnmnxmny显然)()(mnymnxT21ppt课件);1()()()();()()(nxnxnybnnxnya例4.判断下列系统是否为移不变系统。解：系统操作)()()()()(mnnxmnxTnynnxnxT故不是移不变系统。又：函数操作)()()(mnxmnmny显然)()(mnymnxT（a）22ppt课件)1()()()()1()()()1()()(mnxmnxmnxTnynxnxnxTnxnxny故是移不变系统。又：)1()()(mnxmnxmny显然)()(mnymnxT（b)23ppt课件一个常系数线性差分方程是否表征一个线性移不变系统
12、，这一个常系数线性差分方程是否表征一个线性移不变系统，这完全由边界条件决定。完全由边界条件决定。例如：差分方程（c)边界条件 时，既不是线性的也不是移不变的。)()1()(nxnayny（a)边界条件 时，是线性的但不是移不变的。0)0(y（b)边界条件 时，是线性移不变的。0)1(y(1)1y 24ppt课件的情况）解：（0)1(yb)()(1nnx令1)0()1()0(11ayyaayy)1()0()1(11211)2()1()2(aayy.nannayny)()1()(11所以：)()(1nuanyn)()1()(nxnayny25ppt课件)1()(2nnx又令0)1()1()0(22
13、ayy则：.所以：)1()(12nuanyn1)0()0()1(22ayyaayy)1()1()2(22122)1()1()(nannayny可见 是移一位的关系，亦是移一位的关系。因此是移不变系统。)()(21nxnx)()(21nyny26ppt课件)()()()()(111nxTnuanynnxn由上述分析可知：)()1()()1()(2122nxTnuanynnxn)1()()(3nnnx又令：代入差分方程，得：1)10()0()1()0(33ayy27ppt课件1)0()1()0()1(33aayyaaayy233)1()2()1()2(2333)2()3()2()3(aaayy.1
14、3)(nnaany所以：)()()()()(2113nxnxTnuanuanynn)()()()()(21213nxTnxTnynyny因此为线性系统。28ppt课件3.判断系统是否是因果稳定系统。判断系统是否是因果稳定系统。Causal and Noncausal System（因果系统）（因果系统）causal system:(1)响应不出现于激励之前响应不出现于激励之前 (2)h(n)=0,n0 （线性、时不变系统）（线性、时不变系统）Stable System（稳定系统）（稳定系统）(1)有界输入导致有界输出有界输入导致有界输出 (2)（线性、时不变系统）（线性、时不变系统）(3)H(
15、z)的极点均位于的极点均位于Z平面单位圆内（因果系统）平面单位圆内（因果系统）nnh|)(|*实际系统一般是因果系统；*y(n)=x(-n)是非因果系统,因n0时的输入;29ppt课件（b）由于 领先于 ，故为非因果系统。例5 判断下列系统是否为因果系统。)2()()()(nxnxnya)1()()1()(nxnxnyb)()()(kxnycnk)()()(nxnyd（a)为因果系统，由定义可知。)1(ny)(nx解：30ppt课件 由于 由目前和过去的输入所决定，故为因果系统。)(ny由于 n=-1时，有y（-1）=x(1)；也就是 领先于 ，故为非因果系统。)(ny)(nx)()()(kx
16、nycnk)()()(nxnyd31ppt课件第第2章回顾章回顾要点与难点要点与难点1、Z变换变换Z变换的定义、零极点、收敛域变换的定义、零极点、收敛域逆逆Z变换（部分分式法）变换（部分分式法）Z变换的性质及变换的性质及Parseval定理定理2、离散时间傅里叶变换、离散时间傅里叶变换DTFT的定义、性质的定义、性质DTFT与与Z变换的关系变换的关系DTFT存在的条件存在的条件3、DFTDFT定义，与定义，与Z变换的关系，变换的关系，DFT性质性质4、FFT5、DFT的应用的应用32ppt课件nnjjenxeX)(deeXnxnjj)(212.1节知识点节知识点 1、DTFT的定义：的定义：正
17、变换：正变换：反变换：反变换：l基本性质。基本性质。l常见变换对；常见变换对；l离散时间信号的频域（频谱）为周期函数；离散时间信号的频域（频谱）为周期函数；33ppt课件1()()2jjnxnXeednjnjenxeX)()(nnx|)(|Condition:(DTFT)序列傅立叶变换(IDTFT)序列傅立叶反变换注注：周期序列不满足该绝对可和的条件，因此它的周期序列不满足该绝对可和的条件，因此它的DTFT不存在。不存在。1.DTFT的计算及其性质。的计算及其性质。方法1：根据定义式求解34ppt课件一般序列)()()(nxnxnxoe共轭对共轭对称序列称序列共轭反对共轭反对称序列称序列)(*
18、)(21)()(*)(21)(nxnxnxnxnxnxoe一般实序列()()()eox nx nxn偶序列偶序列奇序列奇序列1()()()21()()()2eoxnx nxnxnx nxn1()()*()21()()*()2jjjejjjoXeXeXeXeXeXe方法2：根据DTFT的性质求解（特别是对称性）35ppt课件（a）序列分成实部与虚部时:)()()()()()(jojejireXeXeXnjxnxnx其中nnjrrjeenxnxFTeX)()()(nnjiijoenxjnjxFTeX)()()(序列分成实部与虚部两部分，实部对应的序列分成实部与虚部两部分，实部对应的FT具有共轭对称
19、性，虚部和具有共轭对称性，虚部和j一起对应的一起对应的FT具有共轭具有共轭反对称性。反对称性。36ppt课件)()()()()()(jIjRjoeejXeXeXnxnxnx其中)()(21)(jjjReXeXeX（b）序列分成共轭对称 与共轭反对称 时:)(nxe)(nxo)()(21)(jjjIeXeXejX 序列的共轭对称部分序列的共轭对称部分xe(n)对应着对应着FT的实部的实部XR(ej)，而序列的共轭反对称部分，而序列的共轭反对称部分xo(n)对应着对应着FT的虚部的虚部jXI(ej)。37ppt课件11()1cos122 ()()()1 12 1 01 12jjjRejeRHeee
20、FT hnhnIFT Hennn 例1：若序列h(n)是实因果序列，其DTFT的实部如下式：HR(ej)1+cos 求序列h(n)及其傅里叶变换H(ej).解：38ppt课件000()()()()2()0()2(0)(0)0()2()()()2eeennnh nhnh nh nh nnh nhhnh nh nh nh n为实因果序列时，；时，39ppt课件010()2()0()()0()01 12()1 01 121 1()1 00 ()()11 1eeenjjnjjnjnh nhnnh nhnnh nnhnnnnh nnnH eh n eeee 时，时，时，而 有为 其 它 值22cos2j
21、e40ppt课件2、Z 变换表示法：变换表示法：1)级数形式（定义）级数形式（定义）2)解析表达式解析表达式（根据常见公式）（根据常见公式）（注意（注意:表示收敛域上的函数，同时注明收敛域）表示收敛域上的函数，同时注明收敛域）3、Z 变换收敛域的特点：变换收敛域的特点：1)收敛域是一个圆环，有时可向内收缩到原点，有收敛域是一个圆环，有时可向内收缩到原点，有时可向外扩展到时可向外扩展到，只有，只有x(n)=(n)的收敛域是整个的收敛域是整个Z 平面平面2)在收敛域内没有极点，在收敛域内没有极点，X(z)在收敛域内每一点上在收敛域内每一点上都是解析函数。都是解析函数。()()nnX zx n z
22、41ppt课件4、几类序列、几类序列Z变换的收敛域变换的收敛域(1)有限长序列有限长序列:X(z)=x(n)z-n,(n1 n n2)0 n1 n n2 0|z|展开式出现展开式出现z的负幂的负幂 n1 n n2 0 0|z|展开式出现展开式出现z的正幂的正幂 n1 0 0|z|Rx n1 0,n2=,Rx|z|展开式出现展开式出现z的正幂的正幂Z 变换的收敛域包括变换的收敛域包括 点是因果序列的特征。点是因果序列的特征。42ppt课件(3)左边序列左边序列 X(z)=x(n)z-n,(n1 n n2,n1=-)n1=-,n2 0,|z|0,0|z|Rx,Rx|z|Rx Rx Rx,空集空集4
23、3ppt课件5、部分分式法进行逆、部分分式法进行逆Z变换变换1)求极点求极点2)将将X(z)分解成部分分式形式分解成部分分式形式3)通过查表，对每个分式分别进行逆通过查表，对每个分式分别进行逆Z变换变换注：左边序列、右边序列对应不同收敛域注：左边序列、右边序列对应不同收敛域1)将部分分式逆将部分分式逆Z变换结果相加得到完整的变换结果相加得到完整的x(n)序列序列 6、Z变换的性质变换的性质 移位、反向、乘指数序列、卷积移位、反向、乘指数序列、卷积44ppt课件常用序列常用序列z变换（可直接使用）变换（可直接使用）11111()1|1 z11()0|11()|1z1(1)0|1zNNnnzu n
24、zzzRnzzza u nazazaa unzaa 45ppt课件7、DTFT与与Z变换的关系变换的关系采样序列在单位圆上的采样序列在单位圆上的Z变换等于该序列的变换等于该序列的DTFT 序列频谱存在的条件序列频谱存在的条件Z变换的收敛域包含单位圆变换的收敛域包含单位圆()()()jjjnz enX eX zx n e 8、Parseval定理重要应用定理重要应用计算序列能量：计算序列能量：即时域中对序列求能量与频域中求能量是一致即时域中对序列求能量与频域中求能量是一致221()|()|2jnx nX ed46ppt课件分析计算题（计算证明、分析问答）。分析计算题（计算证明、分析问答）。本章典
25、型题型与习题讲解本章典型题型与习题讲解：47ppt课件48ppt课件方法方法2.幂级数法幂级数法(长除法长除法)左边序列：将左边序列：将X(z)的分子、分母按的分子、分母按Z的升幂排列的升幂排列右边序列：将右边序列：将X(z)的分子、分母按的分子、分母按Z的降幂排列的降幂排列 对于大多数单阶极点的序列，常常用这种部分分式展开法求逆Z变换。方法方法3.部分分式展开法部分分式展开法3.逆逆Z变换的计算。变换的计算。方法方法1.用留数定理求逆用留数定理求逆Z变换变换11Res(),()()knnkkz zX z zzzzX z z求逆求逆z z变换时特别需要注变换时特别需要注意收敛域的范围，收敛域意
26、收敛域的范围，收敛域不同，逆不同，逆z z变换的结果是不变换的结果是不同的。如果没有明确告诉同的。如果没有明确告诉收敛域的范围，则求逆收敛域的范围，则求逆z z变变换时需要讨论。换时需要讨论。49ppt课件1132()11 212X zzz()X z16.已知已知:求出对应求出对应的各种可能的序列的表达式。的各种可能的序列的表达式。解：有两个极点，因为收敛域总是以极点为界，因此收敛域有以下三种情况：三种收敛域对应三种不同的原序列。0.5z 时，时，（1）当收敛域11()()2ncx nX Z zdzj 111115757()()(1 0.5)(1 2)(0.5)(2)nnnzzF zX z z
27、zzzzzz令令，因为，因为c内无极点，内无极点，x(n)=0；0n ，C内有极点内有极点0，但，但z=0是一个是一个n阶极点，阶极点，改为求圆外极点留数，圆外极点有改为求圆外极点留数，圆外极点有1n 50ppt课件那么那么0.52()Re (),0.5Re (),2(57)(57)(0.5)(2)(0.5)(2)(0.5)(2)1 3()2 2 (1)2nnzznnx ns F zs F zzzzzzzzzzzun (57)()(0.5)(2)nzzF zzz0.52z（2）当收敛域）当收敛域时，时，1()Re (),0.53()2nx ns F z 0n，C内有极点0.5；51ppt课件(
28、)Re (),22 2(1)nx ns F zun 1()3()()2 2(1)2nnx nu nun 0n ，C内有极点0.5，0，但0是一个n阶极点，改成求c外极点留数，c外极点只有一个，即2，最后得到最后得到(57)()(0.5)(2)nzzF zzz2z（3）当收敛域）当收敛域52ppt课件1()Re (),0.5Re (),23()2 22nnx ns F zs F znN，则L=M）。较短的一个需要补0至L（两个序列的长度要求相等）。循环卷积可以用DFT(FFT)实现；用循环卷积实现线性卷积：LM+N-1 若不满足这个条件，则只在N-1 n M-1范围内两者相等。73ppt课件典型
29、题型与习题讲解典型题型与习题讲解：分析计算题（计算证明、分析问答、判断）。分析计算题（计算证明、分析问答、判断）。74ppt课件75ppt课件76ppt课件77ppt课件78ppt课件79ppt课件2.4频域采样定理 如果x（n）的长度为M，则只有当频只有当频域采样点数域采样点数N M时时，才有可由频域采样 恢复原序列x（n），否则将产生时域混叠现象。()()()NxnIDFT X kx n)(kX 在z平面的单位圆上的N个等角点上，对z变换进行取样，将导致相应的时间序列周期延拓，延拓周期为N。80ppt课件DFSDFT线性线性线性线性序列移位序列移位循环移位循环移位共轭对称性共轭对称性共轭对
30、称性共轭对称性周期卷积周期卷积循环卷积循环卷积 ()()()()mkNnlNDFS x nmwX kIDFS X klw x n DFSDFSxnXkxnXk ()()()()()mkNnlNNDFT f nwX kIDFT XklRkw x n DFT xnXNk 1010()()()()()()()()()NmNmF kX k Y kf nIDFS F kx m y nmy m x nm 1010()()()()()()()()()()()NNNmNNNmF kX k Y kf nIDFT F kx m y nmRny m x nmRn 81ppt课件DFT选频性选频性DFT与与Z变换变换
31、DFT与与DTFTDFT形式下的形式下的Parseval定理定理2()()kNjkkNNz wX kX zzwe 2()2()/()011()10ojqnjq kjq kNx nenNNkqeX kekq 11*001122001()()()()1|()|()|()()NNnkNNnkx n y nX k YkNx nX kx ny nN 2()()Njww kwNX kX ewN 82ppt课件重新构造两个长度为重新构造两个长度为L的序列的序列x(n)和和y(n),方法：方法：末尾补零末尾补零对对x(n)和和y(n)进行圆周卷积：进行圆周卷积：首先对两个序列进行周期延拓首先对两个序列进行周期
32、延拓对延拓后的周期序列进行周期卷积对延拓后的周期序列进行周期卷积对周期卷积的结果取主值区间对周期卷积的结果取主值区间使圆周卷积等于线性卷积而不产生混淆的必要条使圆周卷积等于线性卷积而不产生混淆的必要条件是件是LN+M-1；步骤如下：步骤如下：83ppt课件圆圆 周周 卷卷 积积 与与 线线 性性 卷卷 积积 的的 性性 质质 对对 比比圆周卷积圆周卷积线性卷积线性卷积针对针对FFT引出的引出的一种一种表示方法表示方法信号通过线性系统时，信信号通过线性系统时，信号输出等于号输出等于输入与系统单输入与系统单位冲激响应的卷积位冲激响应的卷积两序列长度必须两序列长度必须相等相等，不等时按要求不等时按要
33、求补足零值点补足零值点两序列长度可以两序列长度可以不等不等如如x1(n)为为 N1点，点，x2(n)为为 N2点点卷积结果长度卷积结果长度与两信号长度相等皆为与两信号长度相等皆为N卷积结果长度为卷积结果长度为N=N1+N2-184ppt课件 2 2N f、ssf 、sfNkN变量变量周期周期分辨率分辨率数字频域数字频域模拟频域模拟频域离散频域离散频域85ppt课件时域时域/频域同时采样频域同时采样对有限时宽的信号对有限时宽的信号xa(t)的时域波形和频域波形的时域波形和频域波形同时进行取样，其结果是时域波形和频域的都同时进行取样，其结果是时域波形和频域的都变成了离散的、周期性的波形；变成了离散
34、的、周期性的波形；时域内的离散周期信号为时域内的离散周期信号为 ，频域内离散周，频域内离散周期信号为期信号为 ，它们之间形成，它们之间形成DFS变换对；变换对；分别取它们的一个周期，得到分别取它们的一个周期，得到x(n)与与X(k)，它，它们之间形成们之间形成DFT变换对。变换对。nN0k0N-N)nT(x)n(x1/T)k(X)k(X1-N()x n()X k 86ppt课件第二部分第二部分 快速傅里叶变换快速傅里叶变换FFT 1、FFT计算原理。2、基2时间抽取算法和频率抽取算法。3、DFT、R-2 FFT算法的运算量比较。4、实数序列的FFT高效算法。5、FFT的应用。87ppt课件主要
35、要求掌握的内容主要要求掌握的内容：1、FFT、IFFT的计算方法、特点，DIT、DIF的运算流图。2、FFT应用于频谱分析和快速卷积。3、DFT、FFT的运算量计算。4、FFT减少运算量的途径。本章典型题型与习题讲解本章典型题型与习题讲解：作图题（作图、计算）。作图题（作图、计算）。88ppt课件N N点的FFTFFT的运算量为复乘:C CM M=（N/2N/2）M=M=（N/N/2）loglog2 2 N N复加:C CA A=N M=N log=N M=N log2 2 N N1.画出画出N点（例如点（例如8点、点、16点）点）FFT的运算流图的运算流图2.FFT的特点，的特点，FFT减少
36、运算量的途径。减少运算量的途径。DIT DIF3.FFT的运算量的计算，与的运算量的计算，与DFT运算量的比较。运算量的比较。FFT算法的基本思想、特点、编程方法N N点的DFTDFT的运算量为复乘:C CM M=N=N2复加:C CA A=N=N（N-1N-1）89ppt课件例1：如果通用计算机的速度为平均每次复数乘需要5s,每次复数加需要1s,用来计算N1024点DFT，问直接计算需要多少时间。用FFT计算呢?照这样计算，用FFT进行快速卷积对信号进行处理时，估计可实现实时处理的信号最高频率。解：N=1024=210直接计算DFT的运算量：复乘:C CM M=N=N21024102422
37、220次次复加:C CA A=N=N（N-1N-1）102410241023102310475521047552直接计算DFT所用的时间为：-62-65 101024+1047552 10=6.290432sDT 90ppt课件用用FFTFFT计算计算DFTDFT的运算量为复乘:C CM M=（N/2N/2）M=M=（N/N/2）loglog2 2 N N1024/21024/2101051205120复加:C CA A=N M=N log=N M=N log2 2 N N1024102410101024010240用FFT计算DFT所用的时间为：-6-65 105120+10240 10=3
38、5.84msFT 快速卷积时，要计算一次N点FFT（H(k)已经计算好存入ROM中了，不需用FFT计算出H(k)）；N次频域复数乘法（H(k)*X(k)）；一次N点IFFT（也是用FFT实现的）。所以，计算1024点快速卷积的计算时间约为-3-621024 2101024 5 10=76800 scFTT 次复数乘所需时间35.8491ppt课件所以，每秒种处理的采样点数（即采样速率）为-61024=13333.376800 10fs 次/秒max13333.3=22sff=6666.7.3.实数序列的实数序列的FFT高效算法。高效算法。由采样定理可知，可实时处理的信号最高频率为 实际实现时，
39、fmax要比这个小一些。92ppt课件()X k()Y k()x n()y n()X k()Y k()x n()y n3.已知已知和和是两个是两个N点实序列点实序列和和的的DFT，若要从，若要从和和求求和和，为提高运算效率，试设计用一次N点IFFT来完成。)()()()()(kFkFkjYkXkFopep)(Im)(Re)()(nfjnfkFIFFTnf()x n()y n()X k()Y k()Y k解：因为解：因为和和均为实序列，所以，均为实序列，所以，和和为共轭对称序列，为共轭对称序列，j为共轭反对称序列。可令为共轭反对称序列。可令()X k()Y k)(kF和和j分别作为复序列分别作为
40、复序列 分量和共轭反分量和共轭反对称分量，即对称分量，即计算一次计算一次N点点IFFT得到得到93ppt课件)()()()(Im)()()()(RenjykjYIDFTkFIDFTnfjnxkXIDFTkFIDFTnfopep)()(21)()()(21)(nfnfjnynfnfnx由由DFT的共轭对称性可知，的共轭对称性可知，故故94ppt课件2.6节知识点节知识点n连续信号的频谱分析连续信号的频谱分析(利用利用DFT的选频性的选频性)过程：采样截短过程：采样截短DFT效应：混叠效应：混叠原因：采样、频谱泄漏原因：采样、频谱泄漏 泄漏泄漏原因：截短原因：截短 栅栏效应栅栏效应原因：原因：DF
41、TpDFT的分辨率的分辨率 95ppt课件DFT的应用（频谱分析、分段卷积）。的应用（频谱分析、分段卷积）。频谱分析：频谱分析：DFT代替频谱分析引起的误差（混叠现象、栅代替频谱分析引起的误差（混叠现象、栅栏效应、截断效应栏效应、截断效应频谱泄漏、谱间干扰频谱泄漏、谱间干扰）；提高谱分辨）；提高谱分辨率的方法；分段卷积（重叠相加法、重叠保留法）率的方法；分段卷积（重叠相加法、重叠保留法）50FHzminN15.用微处理机对实数序列作谱分析用微处理机对实数序列作谱分析,要求谱分辨率要求谱分辨率，信号最高频率为，信号最高频率为1kHZ，试确定以下各参数：，试确定以下各参数：（3）最少采样点数）最少
42、采样点数（4）在频带宽度不变的情况下，将频率分辨率提高一倍的）在频带宽度不变的情况下，将频率分辨率提高一倍的N值。值。minpT；（1）最小记录时间）最小记录时间maxT（2）最大取样间隔）最大取样间隔；96ppt课件第第3章回顾章回顾要点与难点要点与难点(1)数字滤波器频响应能模仿模拟滤波器频响数字滤波器频响应能模仿模拟滤波器频响 jsasH)(互为映射关系互为映射关系 jezzH)(2)因果稳定的模拟系统变换为数字系统仍为因果稳定的因果稳定的模拟系统变换为数字系统仍为因果稳定的0Re)(sasH互为映射关系互为映射关系1)(zzHS到到Z平面的映射关系满足条件平面的映射关系满足条件 97p
43、pt课件主要内容主要内容：1、数字滤波器的分类及特性。2、数字信号系统的信号流图。3、IIR滤波器的结构和信号流图：直接型；级联型；并联型。4、FIR数字滤波器的结构和信号流图：直接型；快速卷积型、频率采样型。3.1 数字滤波器的结构数字滤波器的结构98ppt课件本章主要要求掌握的内容本章主要要求掌握的内容：1、数字信号系统的信号流图描述方法。2、IIR滤波器的信号流图：直接型；级联型；并联型。3、FIR数字滤波器的实现流图：直接型；级联型；线性相位型。1.画出滤波器的实现结构（实现流图）。画出滤波器的实现结构（实现流图）。99ppt课件IIR数字滤波器的直接数字滤波器的直接I型结构型结构10
44、0ppt课件两条延时链中对应的延时单元内容完全相同两条延时链中对应的延时单元内容完全相同,可合并可合并,得得101ppt课件)(1zH)(zHM102ppt课件103ppt课件2、FIR数字滤波器：非递归结构，无反馈，但在频率数字滤波器：非递归结构，无反馈，但在频率采样结构等某些结构中也包含有反馈的递归部分。采样结构等某些结构中也包含有反馈的递归部分。（1）直接型（卷积型、横截型）直接型（卷积型、横截型）（2）级联型）级联型（3）线性相位型）线性相位型（4）频率采样型）频率采样型104ppt课件 105ppt课件 直接型的转置直接型的转置：106ppt课件107ppt课件FIR数字滤波器数字滤
45、波器要点与难点要点与难点1、线性相位：系统的相频特性是频率的线性函数、线性相位：系统的相频特性是频率的线性函数群时延群时延:gd()d 20)1(N 20)5.0(N2 偶对称偶对称)(nh 奇对称奇对称)(nh108ppt课件2、四种线性相位、四种线性相位FIR滤波器滤波器109ppt课件四种线性相位四种线性相位FIR DF特性特性第一类第一类，h(n)偶、偶、N奇，四种滤波器都可设计。奇，四种滤波器都可设计。第二类第二类，h(n)偶、偶、N偶，可设计低、带通滤波器偶，可设计低、带通滤波器 不能设计高通和带阻。不能设计高通和带阻。第三类第三类，h(n)奇、奇、N奇，只能设计带通滤波器，奇，只
46、能设计带通滤波器，其它滤波器都不能设计。其它滤波器都不能设计。第四类第四类，h(n)奇、奇、N偶，可设计高通、带通滤波偶，可设计高通、带通滤波 器，不能设计低通和带阻。器，不能设计低通和带阻。110ppt课件小结小结1、相位特性只取决于、相位特性只取决于h(n)的对称性，而与的对称性，而与h(n)的的值无关。值无关。2、幅度特性取决于、幅度特性取决于h(n)。3、设计、设计FIR数字滤波器时，在保证数字滤波器时，在保证h(n)对称的条对称的条件下，只要完成幅度特性的逼近即可。件下，只要完成幅度特性的逼近即可。注意：当注意：当H()用用H()表示时，当表示时，当H()为奇对为奇对称时，其相频特性
47、中还应加一个固定相移称时，其相频特性中还应加一个固定相移111ppt课件3、线性相位、线性相位FIR滤波器的零点特性滤波器的零点特性 )1()(nNhnh 11 zHzzHN零点必须是互为倒数的共轭对零点必须是互为倒数的共轭对112ppt课件113ppt课件114ppt课件115ppt课件作图题作图题典型题型与习题讲解典型题型与习题讲解：116ppt课件1.设系统用下面的差分方程描述：设系统用下面的差分方程描述：311()(1)(2)()(1)483y ny ny nx nx n试画出系统的直接型、级联型和并联型结构。试画出系统的直接型、级联型和并联型结构。311()(1)(2)()(1)48
48、3y ny ny nx nx n121311()()()()()483Y zY z zY z zX zX z z112113()31148zH zzz解：解：将上式进行将上式进行Z变换变换117ppt课件（1）按照系统函数）按照系统函数()H z，画出直接型结构如图（一）所示。118ppt课件112113()31148zH zzz11111311(1)(1)24zzz()H z（2）将）将的分母进行因式分解的分母进行因式分解 按照上式可以有两种级联型结构：1111113()11(1)(1)24zH zzz1111113()11(1)(1)24zH zzz(a)(b)画出级联型结构如画出级联型结
49、构如图图（二）（二）（b）所示所示画出级联型结构如图（二）（a）所示119ppt课件级联型结构图（二）（a）级联型结构级联型结构图（二）（图（二）（b）120ppt课件111113()11(1)(1)24zH zzz1()31111()()2424zH zABzzzzz11103()11123()()224zAzzzz1173()11143()()424zBzzzz 107()331124H zzzz()H z（3）将）将进行部分分式展开进行部分分式展开121ppt课件111071073333()1111112424zzH zzzzz根据上式画出并联型结构如图（三）所示。122ppt课件第第2
50、部分部分要点与难点要点与难点(1)数字滤波器频响应能模仿模拟滤波器频响数字滤波器频响应能模仿模拟滤波器频响 jsasH)(互为映射关系互为映射关系 jezzH)(2)因果稳定的模拟系统变换为数字系统仍为因果稳定的因果稳定的模拟系统变换为数字系统仍为因果稳定的0Re)(sasH互为映射关系互为映射关系1)(zzHS到到Z平面的映射关系满足条件平面的映射关系满足条件 123ppt课件主要内容主要内容：1、数字滤波器的设计方法：IIR的设计方法分类。2、理想滤波器的特性及逼近方法：理想滤波器的特性；连续函数逼近方法。3、模拟滤波器设计：几种逼近函数及特点；模拟滤波器逼近函数设计方法。4、模拟滤波器的
