(40)--第5章 离散傅里叶变换应用.ppt

第五章第五章 离散傅里叶变换的应用离散傅里叶变换的应用 用用DFT逼近连续时间信号的频谱逼近连续时间信号的频谱5.1用用FFT计算线卷积和相关运算计算线卷积和相关运算5.2倒频谱分析倒频谱分析5.3系统频谱响应函数分析及确定系统频谱响应函数分析及确定内容提要内容提要v 离散傅里叶变换（离散傅里叶变换（DFTDFT）及其快速算法（）及其快速算法（FFTFFT）的）的重要性不仅在于理论上的严格性，而且还在于工重要性不仅在于理论上的严格性，而且还在于工程上的实用性，凡是可以利用傅里叶变换进行分程上的实用性，凡是可以利用傅里叶变换进行分析、综合和处理的技术问题，都能利用析、综合和处理的技术问题，都能利用FFTFFT有效地有效地解决。本章将详细介绍和分析利用解决。本章将详细介绍和分析利用DFTDFT逼近傅里叶逼近傅里叶变换时存在的问题和解决的方法，如何利用变换时存在的问题和解决的方法，如何利用DFTDFT实实现快速卷积及相关运算，倒频谱的基本概念及应现快速卷积及相关运算，倒频谱的基本概念及应用等内容。用等内容。第一节第一节 用用DFT逼近连续时间信号的频谱逼近连续时间信号的频谱 v工程上所遇到的信号，包括传感器的输出信号，工程上所遇到的信号，包括传感器的输出信号，大多是连续非周期信号，这种信号无论是在时域大多是连续非周期信号，这种信号无论是在时域或频域都是连续的，其波形和频谱如图或频域都是连续的，其波形和频谱如图5-15-1所示。所示。0 0ax(t)t图5-1 连续非周期信号时域波形和频谱（5-1）（5-2）v由式（由式（5-15-1）、式（）、式（5-25-2）和图）和图5-15-1，可以看，可以看出出v1 1）两式中的积分区间均为（两式中的积分区间均为（，）；）；v2 2）和和 都是连续函数。都是连续函数。v显然，上述两点无法满足计算机进行数字显然，上述两点无法满足计算机进行数字信号处理的要求，若要应用信号处理的要求，若要应用FFTFFT进行分析和进行分析和处理，必须在时、频域进行有限化和离散处理，必须在时、频域进行有限化和离散化处理。有限化和离散化处理是在时、频化处理。有限化和离散化处理是在时、频域对被处理的连续信号近似或逼近，是一域对被处理的连续信号近似或逼近，是一种近似处理。种近似处理。主要内容主要内容时域的有限化和离散化时域的有限化和离散化一频域的有限化和离散化频域的有限化和离散化二误差产生原因及解决办法误差产生原因及解决办法三周期信号的数字谱分析周期信号的数字谱分析四谱分析时谱分析时DFT参数的选择参数的选择五频谱细化技术频谱细化技术六一、时域的有限化和离散化一、时域的有限化和离散化 时域的有限化时域的有限化，就是对信号的延续时间沿时间轴，就是对信号的延续时间沿时间轴进行截断，反映在图进行截断，反映在图5-25-2中，是把时间区间由（中，是把时间区间由（，）限定为（）限定为（0 0，）。）。一、时域的有限化和离散化一、时域的有限化和离散化v 时域的离散化时域的离散化，就是对连续信号进行抽样，采样后，就是对连续信号进行抽样，采样后，有有 则则 ，其结果如图，其结果如图5-25-2所示。所示。一、时域的有限化和离散化一、时域的有限化和离散化v那么，原连续信号的频谱离散化后，可近似表示为那么，原连续信号的频谱离散化后，可近似表示为v经有限化，即经有限化，即n n由（由（，）近似为（）近似为（0 0，），上式可表示为上式可表示为v要进行数字谱分析，上式中的还须进行有限化和离散要进行数字谱分析，上式中的还须进行有限化和离散化。化。（5-3）（5-4）二、频域的有限化和离散化二、频域的有限化和离散化 v时域上的变化必然引起频域上的变化，由于在时时域上的变化必然引起频域上的变化，由于在时域上对域上对 进行了抽样，则在频域上将引起频谱进行了抽样，则在频域上将引起频谱的周期化（是原连续信号频谱的周期延拓，延拓的周期化（是原连续信号频谱的周期延拓，延拓周期为周期为 ），如图），如图5-35-3所示。所示。图5-3 时域离散化后的频谱二、频域的有限化和离散化二、频域的有限化和离散化（5-5）与与时时域域一一样样，对对频频域域也也要要进进行行有有限限化化和和离离散散化化处处理理。频频域域的的有有限限化化，是是在在频频域域轴轴上上取取一一个个周周期期的的频频率率区区间间 。频频域域的的离离散散化化，就就是是对对一一个个周周期期内内的的频频谱谱进行抽样，进行抽样，有（有（）,则则 二、频域的有限化和离散化二、频域的有限化和离散化v需要指出需要指出，上式中，上式中，代表信号截断的时间长度，代表信号截断的时间长度，不是信号周期概念，因为原信号是非周期信号；不是信号周期概念，因为原信号是非周期信号；也不是基频的概念，而是频谱离散化后相邻离散也不是基频的概念，而是频谱离散化后相邻离散点的频率间隔。因此为了与周期信号离散谱的符点的频率间隔。因此为了与周期信号离散谱的符号号 相区别，用相区别，用 来表示非周期信号频谱来表示非周期信号频谱离散化后的频谱。其结果分别如图离散化后的频谱。其结果分别如图5-45-4和式（和式（5-5-6 6）所示。）所示。二、频域的有限化和离散化二、频域的有限化和离散化v由上式可知，由上式可知，与与 ，仅相差，仅相差一个系数一个系数 。同。同理可得理可得（5-6）（5-7）图图5-4 5-4 连续信号频谱的有限化和离散化连续信号频谱的有限化和离散化二、频域的有限化和离散化二、频域的有限化和离散化v有了式（有了式（5-65-6）和式（）和式（5-75-7），就可将对连续信号），就可将对连续信号的谱分析用的谱分析来逼近，从而采用的谱分析用的谱分析来逼近，从而采用FFTFFT算法。算法。同时需要说明的是，对某一信号分析时，由于关同时需要说明的是，对某一信号分析时，由于关心的是信号的结构成分，所以只需确定信号中频心的是信号的结构成分，所以只需确定信号中频率的相对量即可，因此频谱计算可以直接采用率的相对量即可，因此频谱计算可以直接采用DFTDFT，式（，式（5-65-6）就是对非周期连续信号进行数字谱分）就是对非周期连续信号进行数字谱分析的基本原理。析的基本原理。三、误差产生原因及解决办法三、误差产生原因及解决办法 v对连续非周期信号的数字谱分析实质上是对连续非周期信号的数字谱分析实质上是用有限长抽样序列的用有限长抽样序列的DFTDFT（离散谱）来近似（离散谱）来近似无限长连续信号的频谱（连续谱），其结无限长连续信号的频谱（连续谱），其结果必然会产生误差，主要的误差包括：果必然会产生误差，主要的误差包括：栅栅栏效应栏效应、混叠效应混叠效应和和频谱泄漏频谱泄漏三种。三种。三、误差产生原因及解决办法三、误差产生原因及解决办法v（一）栅栏效应（一）栅栏效应 非周期信号具有连续谱，但用非周期信号具有连续谱，但用DFTDFT来计算非周期信来计算非周期信号的频谱时，只能观察到有限个（号的频谱时，只能观察到有限个（N N个）离散频谱个）离散频谱值，而频谱间隔中的值就观察不到了，就好像通值，而频谱间隔中的值就观察不到了，就好像通过栅栏观察景物一样，一部分景物被阻挡了，这过栅栏观察景物一样，一部分景物被阻挡了，这种现象称为栅栏效应。将能够感受到的频谱最小种现象称为栅栏效应。将能够感受到的频谱最小间隔值称为频谱分辨率，一般用表示。频谱分辨间隔值称为频谱分辨率，一般用表示。频谱分辨率反映了谱分析算法能将信号中两个靠得很近的率反映了谱分析算法能将信号中两个靠得很近的谱保持分开的能力。若时域抽样周期为，抽样点谱保持分开的能力。若时域抽样周期为，抽样点数为数为N N，则有，则有 三、误差产生原因及解决办法三、误差产生原因及解决办法vN NT T实实际际上上就就是是信信号号在在时时域域上上的的截截断断长长度度 ，分分辨辨率率 与与 成成反反比比。栅栅栏栏效效应应是是由由于于频频域域的的离离散散化化引引起起的的，使使得得在在频频谱谱抽抽样样间间隔隔之之间间的的频频谱谱无无法法反反映映出出来来，因因此此是是不不可可避避免免的的。为为了了改改善善栅栅栏栏效效应应，提提高高频频率率分分辨辨率率，应应当当增增加加信信号号的的有有效效数数据据长长度度 或或N N，也也可可以以采采用用频频谱谱细细化化技技术术，使使谱谱线线变变密密，从从而而 看看 到到 原原 来来 看看 不不 到到 的的“频频 谱谱 景景 象象”。（5-8）三、误差产生原因及解决办法三、误差产生原因及解决办法v（二）混叠效应（二）混叠效应v时域信号的离散化是通过抽样实现的，当采样频时域信号的离散化是通过抽样实现的，当采样频率率 不够高时，采样信号相对原信号就会产不够高时，采样信号相对原信号就会产生频谱的混叠，引起频谱失真。频谱混叠效应是生频谱的混叠，引起频谱失真。频谱混叠效应是由于时域的离散化引起的，克服的办法是提高采由于时域的离散化引起的，克服的办法是提高采样频率，设法满足采样定理，保证样频率，设法满足采样定理，保证 ，其，其中中 是原信号的最高频率。如果时间记录长度为是原信号的最高频率。如果时间记录长度为 ，则在则在 时间内的采样次时间内的采样次N N必须满足必须满足（5-9）三、误差产生原因及解决办法三、误差产生原因及解决办法v（三）频谱泄漏三）频谱泄漏v频谱泄漏又称截断误差，是由于对信号进行截断，频谱泄漏又称截断误差，是由于对信号进行截断，把无限长的信号限定为有限长，即令有限区间外把无限长的信号限定为有限长，即令有限区间外的函数值均为零值，相当于用一个矩形（窗）信的函数值均为零值，相当于用一个矩形（窗）信号乘相应的信号，如图号乘相应的信号，如图5-55-5所示。所示。图图5-5 5-5 用矩形窗截断信号用矩形窗截断信号 三、误差产生原因及解决办法三、误差产生原因及解决办法v由图由图5-55-5得，得，由频域卷积定理有，由频域卷积定理有，信号被截断后的频谱为信号被截断后的频谱为而原信号而原信号 的频谱是的频谱是显然，显然，和和 是不同的。是不同的。三、误差产生原因及解决办法三、误差产生原因及解决办法v例如，设例如，设 ，有，有v画成频谱图，如图画成频谱图，如图5-65-6所示。所示。三、误差产生原因及解决办法三、误差产生原因及解决办法 图图5-6 5-6 频谱泄漏现象频谱泄漏现象三、误差产生原因及解决办法三、误差产生原因及解决办法v余弦信号被矩形窗信号截断后，两根冲激谱线变余弦信号被矩形窗信号截断后，两根冲激谱线变成了以成了以 为中心的为中心的sincsinc形状的连续谱，相当于形状的连续谱，相当于频谱从频谱从 处处“泄漏泄漏”到其他频率处，也就是说，到其他频率处，也就是说，原来一个周期内只有一个频率上有非零值，而现原来一个周期内只有一个频率上有非零值，而现在几乎所有频率上都有非零值，这就是频谱泄漏在几乎所有频率上都有非零值，这就是频谱泄漏现象。现象。v复杂的信号，造成复杂的复杂的信号，造成复杂的“泄漏泄漏”，他们互相叠，他们互相叠加，结果使信号难以分辨。频谱泄漏是由时域信加，结果使信号难以分辨。频谱泄漏是由时域信号的截断引起的，减小频谱泄漏的方法一般有两号的截断引起的，减小频谱泄漏的方法一般有两种。种。v1 1）增加截断长度。）增加截断长度。v2 2）改变窗口形状。）改变窗口形状。四、周期信号的数字谱分析四、周期信号的数字谱分析 v对于周期连续信号对于周期连续信号 ，若其采样序列为，若其采样序列为 ，则由则由DFSDFS与与DFTDFT的关系，周期连续信号的关系，周期连续信号 的频谱的频谱可由下式近似计算可由下式近似计算 v注意注意上式中的上式中的 与推导与推导DFSDFS时的处理方法不同。时的处理方法不同。同理可得同理可得 v上式中的上式中的 是是 的主值序列。的主值序列。（5-10）（5-11）四、周期信号的数字谱分析四、周期信号的数字谱分析v连续周期信号是非时限信号，若要用连续周期信号是非时限信号，若要用FFTFFT作数字谱作数字谱分析，必须分析，必须在时域进行有限化（截断）和离散化在时域进行有限化（截断）和离散化（抽样）处理（抽样）处理。对于一个带限（频谱为有限区间）。对于一个带限（频谱为有限区间）的周期信号，若抽样频率满足抽样条件，并且作的周期信号，若抽样频率满足抽样条件，并且作整周期截断，不会产生频谱的混叠。整周期截断，不会产生频谱的混叠。五、谱分析时五、谱分析时DFT参数的选择参数的选择v应用应用DFT(DFT(或或FFT)FFT)进行信号的频谱分析时，要根据进行信号的频谱分析时，要根据给定的要求，确定给定的要求，确定DFTDFT的参数。一般情况下，已知的参数。一般情况下，已知（或先估计）：信号的最高频率（或先估计）：信号的最高频率 、频谱分辨率、频谱分辨率 、抽样时能够达到的最高抽样频率、抽样时能够达到的最高抽样频率 。需要。需要确定的参数通常包括：截取的信号长度（数据长确定的参数通常包括：截取的信号长度（数据长度）度）、抽样频率、抽样频率 （或采样间隔（或采样间隔 ）、点数）、点数N N及选择什么样的窗口函数等。及选择什么样的窗口函数等。五、谱分析时五、谱分析时DFT参数的选择参数的选择v选择参数的总原则是：选择参数的总原则是：v尽可能减少混叠、频谱泄漏和栅栏效应等项误差，尽可能减少混叠、频谱泄漏和栅栏效应等项误差，保证信号处理的精度和可靠性。在实际分析中，保证信号处理的精度和可靠性。在实际分析中，根据这个原则，通常采用以下的基本步骤来选定根据这个原则，通常采用以下的基本步骤来选定相应的相应的DFTDFT参数。参数。五、谱分析时五、谱分析时DFT参数的选择参数的选择v1 1）估计待分析信号中频率范围和频率上限）估计待分析信号中频率范围和频率上限 。v2 2）选定抽样频率）选定抽样频率 。v3 3）根据分析精度，确定数据有效长度）根据分析精度，确定数据有效长度 。v4 4）确定点数）确定点数N N。v5 5）选窗口。）选窗口。六、频谱细化技术六、频谱细化技术v简单的频谱细化技术之一简单的频谱细化技术之一v加零法加零法:即在原有即在原有N N点有限长序列后面，人为地增点有限长序列后面，人为地增加零点，使总的序列长为加零点，使总的序列长为L(LN)L(LN)，然后进行谱分，然后进行谱分析。析。v注意注意:补零不能提高频率分辨率。补零不能提高频率分辨率。六、频谱细化技术六、频谱细化技术v补零的补零的好处好处：v1 1）可使数据）可使数据N N为为2 2的整数次幂，以便于使用的整数次幂，以便于使用FFTFFT。v2 2）补零起到对）补零起到对DFTDFT的的 做插值的作用，一方面做插值的作用，一方面克服克服“栅栏栅栏”效应，使谱的外观得到平滑；另一效应，使谱的外观得到平滑；另一方面，由于数据截断的过短，频谱泄漏将严重影方面，由于数据截断的过短，频谱泄漏将严重影响对原频谱的辨认，甚至可能在频谱中出现一些响对原频谱的辨认，甚至可能在频谱中出现一些难以确定的峰值，而插值可以在一定程度上克服难以确定的峰值，而插值可以在一定程度上克服这一现象。这一现象。六、频谱细化技术六、频谱细化技术v另一种频谱细化技术是另一种频谱细化技术是vZoomFFTZoomFFT法法:或称选带傅里叶分析法（或称选带傅里叶分析法（Band Select Band Select Fourier AnalysisFourier Analysis），以下简称），以下简称ZFFTZFFT法。它是在法。它是在频率分析范围内任何感兴趣的频率点附近，选择频率分析范围内任何感兴趣的频率点附近，选择一个窄的频带，以高分辨率集中分析这一窄带，一个窄的频带，以高分辨率集中分析这一窄带，从而获得这一段频谱的精细结构。从而获得这一段频谱的精细结构。v如图如图5-75-7，ZFFTZFFT法的原理及基本步骤简述如下：法的原理及基本步骤简述如下：六、频谱细化技术六、频谱细化技术图5-7 ZFFT法的基本步骤六、频谱细化技术六、频谱细化技术v步骤：步骤：v1 1）对被分析信号）对被分析信号 ，按，按 进行采样，采样进行采样，采样点数为点数为 ，则可获得分辨率为，则可获得分辨率为 的频谱的频谱 ，如图，如图5-7a5-7a、b b所示。所示。v2 2）选择感兴趣的中心频率）选择感兴趣的中心频率 及带宽及带宽B B。v3 3）对频谱）对频谱 作数字频移处理，得频移作数字频移处理，得频移 后的后的信号频谱信号频谱 ，如图，如图5-7c5-7c所示。所示。v4 4）对）对 作数字低通滤波，得带宽为作数字低通滤波，得带宽为 的的窄带频谱窄带频谱 ，如图，如图5-7d5-7d所示。所示。六、频谱细化技术六、频谱细化技术v5 5）对）对 进行傅里叶反变换（进行傅里叶反变换（IDFTIDFT），得窄带），得窄带信号信号 。v6 6）对）对 进行重新采样，设采样频率进行重新采样，设采样频率 ，采样点数为采样点数为 ，得，得 ，如图，如图5-7e5-7e所示。则该所示。则该信号的傅里叶变换结果如图信号的傅里叶变换结果如图5-7f5-7f所示。所示。v7 7）对重新采样序列）对重新采样序列 作作FFTFFT，可获得细化频谱，可获得细化频谱 ，如图，如图5-7g5-7g所示。细化后的频率分辨率为所示。细化后的频率分辨率为 。当。当 时，时，表明分辨率提高，表明分辨率提高K K倍。倍。第二节第二节 用用FFT计算线卷积和相关运算计算线卷积和相关运算 v本节将讨论两有限长序列的圆卷积与线卷积的等本节将讨论两有限长序列的圆卷积与线卷积的等价条件，从而将圆卷积的快速算法应用于线性卷价条件，从而将圆卷积的快速算法应用于线性卷积和相关的快速算法中。积和相关的快速算法中。主要内容主要内容圆卷积的计算方法圆卷积的计算方法一圆卷积与线卷积的关系圆卷积与线卷积的关系二用用FFT计算有限长序列的线卷积计算有限长序列的线卷积三分段快速卷积分段快速卷积重叠相加法重叠相加法四离散时间序列的相关运算离散时间序列的相关运算五一、圆卷积的计算方法一、圆卷积的计算方法v两个长度均为两个长度均为 （如长度不等，将短序列补零至（如长度不等，将短序列补零至等长）的有限长序列等长）的有限长序列 、，其圆卷积，其圆卷积 定义为定义为 或或 (5-13)(5-14)其主要计算方法有下面几种。其主要计算方法有下面几种。一、圆卷积的计算方法一、圆卷积的计算方法v1 1公式法（解析法）公式法（解析法）直接利用圆卷积的上述定义式来求解，如下例所直接利用圆卷积的上述定义式来求解，如下例所述。述。例例5-1 5-1 设设 ，计算，计算5 5点点长圆卷积长圆卷积 。解：解：为为4 4点长序列，在其尾部补零使其成为点长序列，在其尾部补零使其成为5 5点长点长序列序列 ，然后进行圆卷积。，然后进行圆卷积。根据卷积定义式，根据卷积定义式，5 5点长序列的圆卷积为点长序列的圆卷积为一、圆卷积的计算方法一、圆卷积的计算方法v则 时v则 时一、圆卷积的计算方法一、圆卷积的计算方法v同理同理 ；。因此，因此，。一、圆卷积的计算方法一、圆卷积的计算方法v2 2图形法图形法v先变量置换，然后保持其中的一个序列不动，将先变量置换，然后保持其中的一个序列不动，将另外一个序列进行反折、圆移位，最后将两个序另外一个序列进行反折、圆移位，最后将两个序列对应的元素相乘、求和。列对应的元素相乘、求和。一、圆卷积的计算方法一、圆卷积的计算方法v 由于圆移位的特点，还可以利用由于圆移位的特点，还可以利用同心圆法同心圆法求圆卷积，如图求圆卷积，如图5-85-8所示。将所示。将 、分别分别分布在两个同心圆上，内圆按顺时针方向分布在两个同心圆上，内圆按顺时针方向刻度刻度 ，外圈按逆时针方向刻度，外圈按逆时针方向刻度 ，并，并 使使 与与 对齐。然后将两个圆上的对对齐。然后将两个圆上的对应值相乘并相加，则得到应值相乘并相加，则得到 。再将外圆。再将外圆按顺时针方向旋转一位，对应值相乘并相按顺时针方向旋转一位，对应值相乘并相加，得到加，得到 ，如此下去，直到求出，如此下去，直到求出 。一、圆卷积的计算方法一、圆卷积的计算方法 图图5-8 5-8 用同心圆作图求圆卷积用同心圆作图求圆卷积 一、圆卷积的计算方法一、圆卷积的计算方法v3 3矩阵法矩阵法v利用矩阵相乘来计算圆卷积。根据圆卷积的定义，利用矩阵相乘来计算圆卷积。根据圆卷积的定义，式（式（5-145-14）可用矩阵表示为）可用矩阵表示为 即即（5-15）一、圆卷积的计算方法一、圆卷积的计算方法v注意注意:式（式（5-155-15）表示序列）表示序列 不动，序列不动，序列 对对 进行了求模运算。为循环矩阵，其元素进行了求模运算。为循环矩阵，其元素排列是有规律的，第一行表示的是原点（排列是有规律的，第一行表示的是原点（）不动时，序列的反折（倒序）；接下来的各行分不动时，序列的反折（倒序）；接下来的各行分别是上一行序列的循环右移。别是上一行序列的循环右移。一、圆卷积的计算方法一、圆卷积的计算方法v4 4时域圆卷积定理法时域圆卷积定理法v利用利用DFTDFT的一个重要性质的一个重要性质时域圆卷积定理来计时域圆卷积定理来计算圆卷积。若算圆卷积。若 ，则则二、二、圆卷积与线卷积的关系圆卷积与线卷积的关系 v线卷积具有明确的物理意义，直接计算比较复杂。线卷积具有明确的物理意义，直接计算比较复杂。对于两个有限长序列求线卷积能否用圆卷积来代对于两个有限长序列求线卷积能否用圆卷积来代替，即采用替，即采用FFTFFT计算线卷积而使两者结果又完全相计算线卷积而使两者结果又完全相同呢？答案是肯定的，但需要满足一个条件：就同呢？答案是肯定的，但需要满足一个条件：就是将进行线卷积的两序列的长度（设两序列的点是将进行线卷积的两序列的长度（设两序列的点数分别为数分别为 ）均通过补零的办法，加长至）均通过补零的办法，加长至 ，然后再进行点的圆卷积，则圆卷积的结果与线，然后再进行点的圆卷积，则圆卷积的结果与线卷积的结果相同。卷积的结果相同。二、二、圆卷积与线卷积的关系圆卷积与线卷积的关系v设设 分别由点通过补零，加长至分别由点通过补零，加长至N N点，其线点，其线卷积为，可表示为卷积为，可表示为 计算结果的长度要多出一些零值，但非零值长度计算结果的长度要多出一些零值，但非零值长度仍为仍为 点。点。其圆卷积为其圆卷积为 ，可表示为，可表示为二、二、圆卷积与线卷积的关系圆卷积与线卷积的关系v而而 可得可得 （5-165-16）式中，下标式中，下标p p表示序列的周期化；表示序列的周期化；是指对线卷是指对线卷积积 进行周期为进行周期为N N的延拓后得到的周期序的延拓后得到的周期序 列；列；两序列的圆卷积的结果，是的主值序列。两序列的圆卷积的结果，是的主值序列。二、二、圆卷积与线卷积的关系圆卷积与线卷积的关系v上述过程说明上述过程说明:加长至加长至N N点长的点长的 、两序列的两序列的圆卷积圆卷积 与与 线卷积作周期延拓所得到的序线卷积作周期延拓所得到的序列列 的主值序列相同。在这个条件下（两序列的主值序列相同。在这个条件下（两序列均加长至均加长至N N点），就可以通过计算序列的圆卷积来点），就可以通过计算序列的圆卷积来求解线卷积。从式（求解线卷积。从式（5-165-16）的推导过程还可以看）的推导过程还可以看出，如果两序列不加长至出，如果两序列不加长至N N，其线卷积的周期延拓，其线卷积的周期延拓序列将发生重叠或混叠现象（因为序列将发生重叠或混叠现象（因为 线卷积长线卷积长度为度为 ），相应计算出的圆卷积也将产生），相应计算出的圆卷积也将产生失真，圆卷积的主值序列和线卷积就不相同。失真，圆卷积的主值序列和线卷积就不相同。三、三、用用FFT计算有限长序列的线卷积计算有限长序列的线卷积v根据上述圆卷积与线卷积的关系，可以得出用根据上述圆卷积与线卷积的关系，可以得出用FFTFFT求解两求解两序列线卷积的原理框图，如图序列线卷积的原理框图，如图5-9 5-9 所示。其计算的所示。其计算的具体步具体步骤如下骤如下v1 1）若两序列）若两序列 、的长度为的长度为N N，将序列加长至，将序列加长至2N2N1 1，并应修正为，并应修正为2 2的幂次（基的幂次（基2 2算法）；算法）；v2 2）计算）计算 、；v3 3）计算）计算 ；v4 4）计算）计算 。图图5-9 5-9 用用FFTFFT求线性卷积求线性卷积三、三、用用FFT计算有限长序列的线卷积计算有限长序列的线卷积v在在MATLABMATLAB中直接实现线卷积计算的函数有中直接实现线卷积计算的函数有 convconv,conv2,conv2,convnconvn。其中。其中conv2conv2和和 convnconvn分别用于分别用于2 2维、维、n n维的卷积运算。维的卷积运算。convconv则用于向量卷积与多项式乘则用于向量卷积与多项式乘的计算，调用的格式为的计算，调用的格式为c cconvconv（a a，b b）。式中，）。式中，a a、b b表示两个序列，表示两个序列，c ca a*b b。在。在MATLABMATLAB中，序列中，序列可用向量来表示，若向量可用向量来表示，若向量a a的长度为的长度为nana，向量，向量b b的的长度为长度为nbnb，则向量，则向量c c的长度为的长度为nananbnb1 1。四、分段快速卷积四、分段快速卷积重叠相加法重叠相加法v分段快速卷积的方法分段快速卷积的方法:将长序列分成若干小段，每将长序列分成若干小段，每小段分别与短序列作卷积运算，然后将所有的分小段分别与短序列作卷积运算，然后将所有的分段卷积结果相叠加，就是线卷积的最后结果，这段卷积结果相叠加，就是线卷积的最后结果，这种方法又称为重叠相加法。种方法又称为重叠相加法。四、分段快速卷积四、分段快速卷积重叠相加法重叠相加法v设设 的长度为的长度为M M，为一长序列，将为一长序列，将 进行分进行分段，每段的长度为段，每段的长度为 ，将每一段分别与，将每一段分别与 进行进行线卷积，然后将结果重叠相加，如图线卷积，然后将结果重叠相加，如图5-105-10所示。所示。图5-10 叠加相加法的分段以及的重叠情况 四、分段快速卷积四、分段快速卷积重叠相加法重叠相加法v设将设将 分为分为 第第 段段 表示为表示为 （5-175-17）v则则 （5-185-18）四、分段快速卷积四、分段快速卷积重叠相加法重叠相加法v由于由于 长度为长度为 ，的长度为的长度为M M，故，故 的长度的长度为为 ，即的范围为，即的范围为 （5-195-19）v将式（将式（5-195-19）与式（）与式（5-175-17）的范围比较，的范围比较，显显然然 比长点，而比长点，而 的范围是的范围是 （5-20）四、分段快速卷积四、分段快速卷积重叠重叠相加法相加法v将式（将式（5-195-19）与式（）与式（5-205-20）比较，可知的）比较，可知的 后后部分与的部分与的 前部分，有前部分，有 个点发生重叠。个点发生重叠。这样，对于在此范围的每一个这样，对于在此范围的每一个 值，原序列值，原序列 和和 的卷积的卷积 之值应为之值应为 （5-215-21）v这就是说，式（这就是说，式（5-185-18）中的求和并不是将各段线）中的求和并不是将各段线卷积的结果简单地拼接在一起，在某些点上是需卷积的结果简单地拼接在一起，在某些点上是需要前后两段的结果重叠相加的。要前后两段的结果重叠相加的。五、离散时间序列的相关运算五、离散时间序列的相关运算v设序列设序列 ，（，（），称下述运算），称下述运算 为序列为序列 和和 的线性相关。的线性相关。v对对 和和 均为实序列的情形，由式（均为实序列的情形，由式（5-225-22）可以得到可以得到 (5-22)（5-23）五、离散时间序列的相关运算五、离散时间序列的相关运算v若序列若序列 和和 是不同的两个序列，则称相关为是不同的两个序列，则称相关为互相关，若互相关，若 ，则称相关为，则称相关为自相关自相关。v注意：注意：相关和卷积是两个不同的概念，它相关和卷积是两个不同的概念，它们的区别是明显的。除了形式上的差异，们的区别是明显的。除了形式上的差异，它们还有一些其他的不同。例如，对相关它们还有一些其他的不同。例如，对相关运算来说，既不具有交换性，也不具有结运算来说，既不具有交换性，也不具有结合性。一般来说，合性。一般来说，第三节第三节 倒频谱分析倒频谱分析倒频谱的定义倒频谱的定义一倒频谱的应用倒频谱的应用对语言信号的分析对语言信号的分析二一、倒频谱的定义一、倒频谱的定义v设时域连续信号设时域连续信号x(tx(t)的傅里叶变换为的傅里叶变换为 其功率谱为其功率谱为一、倒频谱的定义一、倒频谱的定义v定义定义 （5-255-25）为连续信号为连续信号x x(t(t)的的倒频谱倒频谱，它实质是，它实质是“信号对数信号对数功率谱的功率谱功率谱的功率谱”。v实际工程中常用幅值倒频谱，其表达式为实际工程中常用幅值倒频谱，其表达式为 （5-265-26）一、倒频谱的定义一、倒频谱的定义v在实际数字信号处理时，对有限长序列的倒谱计在实际数字信号处理时，对有限长序列的倒谱计算算步骤为步骤为v(1)(1)对时域信号对时域信号 作作DFT DFT ；v(2)(2)对频域信号对频域信号 取对数：取对数：；v(3)(3)求倒谱：求倒谱：。二、倒频谱的应用二、倒频谱的应用对语言信号的分析对语言信号的分析v元音元音“a a”的对数谱和倒频谱表示在图的对数谱和倒频谱表示在图5-115-11上。从上。从图中可以看到有两个特点：一是有大量的谐波分图中可以看到有两个特点：一是有大量的谐波分量，谐波间距等于语音音调；二是有许多共振峰，量，谐波间距等于语音音调；二是有许多共振峰，即所谓的构形成分，它由声道的形状决定，并确即所谓的构形成分，它由声道的形状决定，并确定了特定的元音声。定了特定的元音声。二、倒频谱的应用二、倒频谱的应用对语言信号的分析对语言信号的分析图图5-11 5-11 元音元音“a a”的对数谱和倒频谱分析的对数谱和倒频谱分析二、倒频谱的应用二、倒频谱的应用对语言信号的分析对语言信号的分析v为分析方便，用为分析方便，用 代替代替 ，则，则可用可用 表示原来声道内发出的语音信号的表示原来声道内发出的语音信号的功率谱，用功率谱，用 表示共振嗓音成分的功率谱，两者表示共振嗓音成分的功率谱，两者合成的元音声的功率谱为合成的元音声的功率谱为v若以对数形式表达，上式可改写成若以对数形式表达，上式可改写成二、倒频谱的应用二、倒频谱的应用对语言信号的分析对语言信号的分析v因为傅里叶变换的线性特性，所以在倒频谱中仍因为傅里叶变换的线性特性，所以在倒频谱中仍保持相加的关系。保持相加的关系。v并简写成并简写成二、倒频谱的应用二、倒频谱的应用对语言信号的分析对语言信号的分析v从图从图5-115-11中还可以看出，有声道产生的构形成分中还可以看出，有声道产生的构形成分与嗓音产生的语音特征，在倒频谱中处于完全不与嗓音产生的语音特征，在倒频谱中处于完全不同的地方，可以明显地加以区别。同的地方，可以明显地加以区别。第四节第四节 系统频率响应函数分析及确定系统频率响应函数分析及确定v本节主要介绍一下系统频率响应函数的基本特性、本节主要介绍一下系统频率响应函数的基本特性、基本测定方法，并以传感器的频率特性分析为例，基本测定方法，并以传感器的频率特性分析为例，说明离散傅里叶变换（说明离散傅里叶变换（DFTDFT）的应用。）的应用。主要内容主要内容频率响应函数的基本特性频率响应函数的基本特性一频率响应函数的测定频率响应函数的测定二传感器的频率特性分析传感器的频率特性分析三一、频率响应函数的基本特性一、频率响应函数的基本特性v频率响应函数的频率响应函数的物理意义物理意义：首先假定输入是频率：首先假定输入是频率 为为 的正弦波，输出也是一个相同频率的波。同的正弦波，输出也是一个相同频率的波。同一频率的输出振幅之比，等于系统的模一频率的输出振幅之比，等于系统的模 ；输；输出和输入的相位差等于系统的相角出和输入的相位差等于系统的相角 。频率响。频率响应函数就是系统的幅频和相频特性。应函数就是系统的幅频和相频特性。一、频率响应函数的基本特性一、频率响应函数的基本特性v频率响应函数的频率响应函数的基本特性基本特性:1 1），即幅频特性为正值函数。，即幅频特性为正值函数。2 2）与与 互为傅里叶变换对。互为傅里叶变换对。3 3）如果系统为单输入，并没有明显的噪声输）如果系统为单输入，并没有明显的噪声输 入。入。4 4）频率响应函数的相频特性与互谱的相位特性）频率响应函数的相频特性与互谱的相位特性完全相同。完全相同。一、频率响应函数的基本特性一、频率响应函数的基本特性v5 5）频率响应函数的幅频特性图上，幅比值可以用）频率响应函数的幅频特性图上，幅比值可以用线性坐标，其单位为输出信号的物理单位除以输线性坐标，其单位为输出信号的物理单位除以输入信号的物理单位。入信号的物理单位。二、频率响应函数的测定二、频率响应函数的测定试验测定常用的方法有试验测定常用的方法有3 3种种:v（1 1）正弦和随机扫频激振法）正弦和随机扫频激振法v（2 2）冲击法）冲击法v（3 3）原型工况实测法）原型工况实测法三、传感器的频率特性分析三、传感器的频率特性分析v以压力传感器为例，一般而言，他可以看作是一线性以压力传感器为例，一般而言，他可以看作是一线性时不变系统，如图时不变系统，如图5-12 5-12 所示。设传感器的单位冲激所示。设传感器的单位冲激响应为响应为 ，输入为，输入为 ，则输出响应，则输出响应 ，根据卷积定义有，根据卷积定义有图图5-12 5-12 传感器模型传感器模型三、传感器的频率特性分析三、传感器的频率特性分析v对动态压力传感器频率特性的测定，目前通用的对动态压力传感器频率特性的测定，目前通用的实验方法如图实验方法如图5-135-13所示。所示。压力信号 源被 测传感器抗混迭滤波器采 样A/DFFTD/A打印输出数据图5-13 传感器动态标定系统框图三、传感器的频率特性分析三、传感器的频率特性分析v在图在图5-135-13中，输入的压力信号源可以是各种形式中，输入的压力信号源可以是各种形式的正弦压力信号发生器，或激波管；虚线框内可的正弦压力信号发生器，或激波管；虚线框内可以是专门的测试信号分析仪，也可以是由瞬态记以是专门的测试信号分析仪，也可以是由瞬态记录仪，通用数字计算机等组成的信号分析系统。录仪，通用数字计算机等组成的信号分析系统。对应于后者的一个实用试验系统原理框图如图对应于后者的一个实用试验系统原理框图如图5-5-1414所示。所示。三、传感器的频率特性分析三、传感器的频率特性分析激波管被 测传感器瞬态记录仪（记录、采样,A/D）通用数字计算机FFT及其他处理程序打 印画出的 曲 线输出响应y（t）图5-14 传感器动态标定系统原理框图三、传感器的频率特性分析三、传感器的频率特性分析v在利用图在利用图5-145-14所示信号分析系统对传感器频率特所示信号分析系统对传感器频率特性进行分析的情况下，通常其输入信号性进行分析的情况下，通常其输入信号 （由激（由激波管产生）是单位阶跃信号。只要对传感器的阶波管产生）是单位阶跃信号。只要对传感器的阶跃输出响应跃输出响应 进行傅里叶变换，便可得到传感进行傅里叶变换，便可得到传感器的频率响应特性，即器的频率响应特性，即 ，从而根据，从而根据 可以求得传感器的频率特性参数。可以求得传感器的频率特性参数。三、传感器的频率特性分析三、传感器的频率特性分析v但对传感器而言，其频率响应总是有限的，因此，但对传感器而言，其频率响应总是有限的，因此，在激波管产生的阶跃信号输入情况下，传感器的在激波管产生的阶跃信号输入情况下，传感器的输出响应输出响应 可能如图可能如图5-15 b5-15 b所示。在这种情况下，所示。在这种情况下，是非时限信号，且归一化稳态值为是非时限信号，且归一化稳态值为1 1，不满足绝对，不满足绝对可积的条件。可积的条件。三、传感器的频率特性分析三、传感器的频率特性分析v为了利用为了利用DFTDF