动态信号处理(第四、五、六章).ppt
《动态信号处理(第四、五、六章).ppt》由会员分享,可在线阅读,更多相关《动态信号处理(第四、五、六章).ppt(307页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、目录v第四章第四章频谱分析频谱分析v第五章第五章数字滤波数字滤波 v第六章第六章细化选带频谱分析细化选带频谱分析 12/29/2022 简介简介 卷积和相关定理卷积和相关定理 连续傅立叶变换连续傅立叶变换 几种典型信号的傅立叶变换几种典型信号的傅立叶变换 傅立叶变换的性质傅立叶变换的性质 频谱的作法频谱的作法(离散傅立叶变换离散傅立叶变换DFT)DFT)和基本用途和基本用途 从连续傅立叶变换到离散傅立叶变换从连续傅立叶变换到离散傅立叶变换 快速傅立叶变快速傅立叶变FFTFFT(Fast(Fast Fourier Transform)Fourier Transform)窗函数与泄漏窗函数与泄漏
2、提高频谱分析精度的其他措施提高频谱分析精度的其他措施(多段平均和零均值化等多段平均和零均值化等)第四章第四章 频谱分析频谱分析12/29/2022 频谱分析是现代信号处理技术的最基本和常用的方法之频谱分析是现代信号处理技术的最基本和常用的方法之一,在生产实践和科学研究中获得日益广泛的应用。例如,一,在生产实践和科学研究中获得日益广泛的应用。例如,对汽车、飞机、轮船、汽轮机等各类旋转机械、电机、机床对汽车、飞机、轮船、汽轮机等各类旋转机械、电机、机床等机器的主体或部件进行实际运行状态下的谱分析,可以提等机器的主体或部件进行实际运行状态下的谱分析,可以提供设计数据和检验设计效果,或者寻找振源和诊断
3、故障,保供设计数据和检验设计效果,或者寻找振源和诊断故障,保证设备的安全运行等;在声纳系统中,为了寻找海洋水面船证设备的安全运行等;在声纳系统中,为了寻找海洋水面船只或潜艇,需要对噪声信号进行谱分析,以提供有用信息,只或潜艇,需要对噪声信号进行谱分析,以提供有用信息,判断舰艇运动速度、方向、位置、大小等。判断舰艇运动速度、方向、位置、大小等。谱分析方法研究是信号处理技术中一个活跃课题。谱分析方法研究是信号处理技术中一个活跃课题。在频谱分析中需要考虑两个重要的问题:在频谱分析中需要考虑两个重要的问题:提高分析速度,包括采用提高分析速度,包括采用FFT FFT 等各类算法、减少分析点等各类算法、减
4、少分析点 数等;数等;提高分析精度,这需要采用加窗、频谱修正等技术压低提高分析精度,这需要采用加窗、频谱修正等技术压低 旁瓣、抑制泄漏或修正时域截断的影响。旁瓣、抑制泄漏或修正时域截断的影响。12/29/20224.1 4.1 卷积与相关定理卷积与相关定理 卷积卷积 相关分析相关分析12/29/2022 1.卷积积分的定义卷积积分的定义 两个函数的卷积或褶积两个函数的卷积或褶积(convolution),在不在不同的科学领域中,同的科学领域中,代表着不同的物理概念。在这代表着不同的物理概念。在这里,我们首先仅仅把其作为一种数学运算方法来里,我们首先仅仅把其作为一种数学运算方法来研究。研究。函数
5、函数x(t)与与h(t)的卷积积分定义为:的卷积积分定义为:一一.卷积卷积12/29/2022 2.卷积积分的几种图形表示卷积积分的几种图形表示12/29/2022 卷积的图解表示,可以使人们能够直接观察到许卷积的图解表示,可以使人们能够直接观察到许多抽象关系的结果。多抽象关系的结果。函数函数x(t)与与h(t)的卷积运算过程,包含四个部骤:的卷积运算过程,包含四个部骤:(1).反射反射把把h()相对纵轴折叠过去,得到相对纵轴折叠过去,得到h()的镜象的镜象h(-);(2).平移平移将将h(-)移动移动t1,得到得到h(t1-);移动移动t2,得到得到h(t2-);以下类推。以下类推。(3).
6、相乘相乘将位移后的函数将位移后的函数h(ti-)乘以乘以x();(4).积分积分h(t-)与与x()乘积曲线下所包含面积即乘积曲线下所包含面积即为为t时刻的卷积值,时刻的卷积值,A1 代表代表x()h(t1-)下的面积,是下的面积,是t1时刻的卷积值;时刻的卷积值;A2代表代表x()h(t2-)下的面积,是下的面积,是t2时刻时刻的卷积值。的卷积值。当当h(-)沿整个时间轴经过上述四个步骤后,即可沿整个时间轴经过上述四个步骤后,即可得得x(t)与与h(t)的卷积值。的卷积值。12/29/20223.物理意义物理意义 利用卷积运算,描述线性时不变系统的输出与输入的关系,利用卷积运算,描述线性时不
7、变系统的输出与输入的关系,在物理概念上是十分清楚的,即系统的输出在物理概念上是十分清楚的,即系统的输出y(t)是任意输入是任意输入x(t)与与系统脉冲响应系统脉冲响应h(t)的卷积。这一运算过程可概括为下述三个内容:的卷积。这一运算过程可概括为下述三个内容:(1).信号信号x(t)可以表示为许多宽度为可以表示为许多宽度为t 的窄条面积之和,的窄条面积之和,t=nt时的第时的第n 个窄条宽度为个窄条宽度为x(nt)t,在在t0的情况下,窄的情况下,窄条可以看作是强度等于窄条面积的脉冲。条可以看作是强度等于窄条面积的脉冲。(2).如果已知系统的单位脉冲响应是如果已知系统的单位脉冲响应是h(t),根
8、据线性系统的齐根据线性系统的齐次性与时不变性,在次性与时不变性,在t=nt时刻,窄条脉冲引起的响应是:时刻,窄条脉冲引起的响应是:(3).根据线性系统的叠加性,由所有脉冲分量引起的响应根据线性系统的叠加性,由所有脉冲分量引起的响应是许多窄条脉冲响应之和,即是许多窄条脉冲响应之和,即 12/29/2022当当t0时,离散和变为积分,并令时,离散和变为积分,并令td,nt,则有则有积分式:积分式:用卷积方法求任意输入下系统的响应用卷积方法求任意输入下系统的响应12/29/2022 函数的筛选性质函数的筛选性质 4.含有脉冲函数的卷积含有脉冲函数的卷积 函数的定义函数的定义12/29/2022含有脉
9、冲函数的卷积含有脉冲函数的卷积其它其它12/29/2022计算:计算:根据根据函数的筛选特性,得:函数的筛选特性,得:用作图法,将用作图法,将x(t)反折反折,然后沿时间轴然后沿时间轴(-,)作位移、相乘、积分,作位移、相乘、积分,也可得到同样结果。也可得到同样结果。由此可知,计算函数由此可知,计算函数x(t)与脉冲函数的卷积,与脉冲函数的卷积,就是简单地将就是简单地将x(t)在产生脉冲函数的坐标位置上重在产生脉冲函数的坐标位置上重新构图。即信号波形对单位脉冲信号的卷积:以新构图。即信号波形对单位脉冲信号的卷积:以脉冲为中心复制信号波形。脉冲为中心复制信号波形。12/29/2022 5.卷积定
10、理卷积定理 卷积定理是卷积与傅里叶变换之间关系的重要定理,具有时卷积定理是卷积与傅里叶变换之间关系的重要定理,具有时域卷积定理和频域卷积定理两个定理,是傅里叶变换的一个重要域卷积定理和频域卷积定理两个定理,是傅里叶变换的一个重要特性,它在推导和论述相关定理中有重要地位。特性,它在推导和论述相关定理中有重要地位。(1).时域卷积定理时域卷积定理证明证明:两个时间函数卷积的频谱等于各个时间函数频谱的复数乘两个时间函数卷积的频谱等于各个时间函数频谱的复数乘积,即在时域中两个信号的卷积,等效于在频域中频谱相乘。积,即在时域中两个信号的卷积,等效于在频域中频谱相乘。12/29/2022 下面以几个例子来
11、说明时域卷积定理的应用。下面以几个例子来说明时域卷积定理的应用。例例1 脉冲序列与矩形脉冲的卷积与傅里叶变换之间的关系。脉冲序列与矩形脉冲的卷积与傅里叶变换之间的关系。脉冲序列脉冲序列h(t)的傅里叶变换仍为脉冲序列;的傅里叶变换仍为脉冲序列;矩形脉冲函数的傅里叶变换是矩形脉冲函数的傅里叶变换是sinc(t)型函数。型函数。脉冲序列与矩形脉冲的卷积与傅立叶变换之间的关系脉冲序列与矩形脉冲的卷积与傅立叶变换之间的关系 12/29/2022(2).频域卷积定理频域卷积定理 此即为频域卷积定理,它说明两时间函数的频谱此即为频域卷积定理,它说明两时间函数的频谱的卷积等效于时域两函数的乘积。需要指出的是
12、其中的卷积等效于时域两函数的乘积。需要指出的是其中频域卷积为复频谱卷积,不是直接的谱频域卷积为复频谱卷积,不是直接的谱(自谱、幅值自谱、幅值谱和对数谱)卷积。谱和对数谱)卷积。复数相乘模相乘,幅角相加复数相乘模相乘,幅角相加12/29/2022 例例2:调制信号:调制信号2调制信号调制信号(谐波信号直流信号谐波信号直流信号)的卷积与傅的卷积与傅里叶变换之间的关系。里叶变换之间的关系。调制信号调制信号2调制信号调制信号(谐波信号直流信号谐波信号直流信号)的卷积与傅立叶变换之间的关系的卷积与傅立叶变换之间的关系12/29/2022例例3利用频域卷积定理研究余弦信号截断利用频域卷积定理研究余弦信号截
13、断(加矩形窗加矩形窗)后的频谱,后的频谱,此例表明了无限长余弦信号的频域能量集中此例表明了无限长余弦信号的频域能量集中-1/T与与1/T,而截断后的余弦信号的频域能量则在而截断后的余弦信号的频域能量则在-1/T与与1/T点附近分点附近分散。实质上,任一谐波信号进行离散傅里叶变换时其频谱散。实质上,任一谐波信号进行离散傅里叶变换时其频谱都不再是都不再是函数,而是函数,而是sinc(t)型函数型函数.12/29/2022 6.用频域卷积定理证明巴什瓦定理用频域卷积定理证明巴什瓦定理巴什瓦巴什瓦(Parseval)定理是信号分析中的一个重要定理,它表定理是信号分析中的一个重要定理,它表明了时域中信号
14、明了时域中信号X(t)的总能量等于频域中计算的信号总能量。的总能量等于频域中计算的信号总能量。如果:如果:则有则有证明:根据卷积定理有证明:根据卷积定理有 令令q0,x(t)h(t),则有则有 等式左边表示了时域中信号总能量,关于等式右边,根据傅里叶等式左边表示了时域中信号总能量,关于等式右边,根据傅里叶变换性质可知,变换性质可知,X(f)与与X(-f)共轭,而有共轭,而有 因此可证明因此可证明12/29/2022 7.离散卷积及实现方法离散卷积及实现方法 设有设有N点离散等间隔时间序列点离散等间隔时间序列x(i)和和h(i)(i=0,1,2,N-1),其离散卷积的定义为:其离散卷积的定义为:
15、下面以两个矩形脉冲信号为例说明作离散卷积的步骤:下面以两个矩形脉冲信号为例说明作离散卷积的步骤:两个矩形脉冲信号的离散卷积的步骤两个矩形脉冲信号的离散卷积的步骤 12/29/2022离散频域卷积定理:离散频域卷积定理:离散时域卷积定理:离散时域卷积定理:与连续信号的卷积定理不同的是频域卷积多出一系与连续信号的卷积定理不同的是频域卷积多出一系数数1/N。12/29/2022二、相关分析二、相关分析 相关分析包括自相关函数和相关分析包括自相关函数和互相关函数。相关函数可以描述互相关函数。相关函数可以描述两个信号之间有无关系或其相似两个信号之间有无关系或其相似程度,也可以描述同一信号的现程度,也可以
16、描述同一信号的现在值与过去有无关系,或根据过在值与过去有无关系,或根据过去值、现在值估计未来值。去值、现在值估计未来值。12/29/20221.1.相关函数定义相关函数定义相关函数定义相关函数定义 复函数复函数x(t)的自相关函数的定义为:的自相关函数的定义为:当当x(t)为实函数,其自相关函数可简写为:为实函数,其自相关函数可简写为:复函数复函数x(t)与与y(t)互相关函数的定义为:互相关函数的定义为:当当x(t)与与y(t)为实函数时,其互相关函数可简为实函数时,其互相关函数可简 写为:写为:12/29/20222.相关函数用途相关函数用途 相关测量在工程上有许多应用。相关测量在工程上有
17、许多应用。自相关函数可用于检测淹没在噪声中自相关函数可用于检测淹没在噪声中的周期信号,因为周期函数的自相关的周期信号,因为周期函数的自相关仍然是周期函数,而随机噪声的自相仍然是周期函数,而随机噪声的自相关函数则很快衰减。而互相关函数可关函数则很快衰减。而互相关函数可以用于检测两个同频率信号的相位差以用于检测两个同频率信号的相位差和两个信号的同周期成分。和两个信号的同周期成分。12/29/20223.3.相关与卷积的关系相关与卷积的关系相关与卷积的关系相关与卷积的关系 x(t)与与y(t)的卷积为:的卷积为:为便于比较,把相关函数的变量为便于比较,把相关函数的变量t与与互换,这样互相互换,这样互
18、相关函数可表示为:关函数可表示为:比较两式可知:比较两式可知:故有,故有,上式表明,上式表明,x(t)与与y(t)的互相关函数等于的互相关函数等于x(t)与与y(-t)的的卷积。卷积。12/29/2022 这两种运算过程都包括了位移、相乘、积分三个步这两种运算过程都包括了位移、相乘、积分三个步骤,其差别在于相关运算不要对骤,其差别在于相关运算不要对y(t)进行反折,而卷积进行反折,而卷积运算需要反折。运算需要反折。卷积和相关的运算图形比较卷积和相关的运算图形比较 12/29/20224.4.相关定理相关定理相关定理相关定理证明:证明:同理可证明:同理可证明:12/29/2022 当当x(t)y
19、(t)时,则有时,则有 两信号的互相关函数的傅里叶变换等于两信号的互相关函数的傅里叶变换等于其中一个信号的傅里叶变换乘以另一信号的其中一个信号的傅里叶变换乘以另一信号的傅里叶变换的共轭值,此即相关定理。如果傅里叶变换的共轭值,此即相关定理。如果y(t)为实偶函数,则它的傅里叶变换为实偶函数,则它的傅里叶变换Y()是实是实函数,此时相关定理与卷积定理具有相同结函数,此时相关定理与卷积定理具有相同结果果。相关定理相关定理相关定理相关定理12/29/20225.5.用相关定理证明能量定理用相关定理证明能量定理用相关定理证明能量定理用相关定理证明能量定理(巴什瓦定理巴什瓦定理巴什瓦定理巴什瓦定理)能量
20、定理:时域内能量定理:时域内x2(t)曲线所覆盖的面曲线所覆盖的面积等于频域内积等于频域内 覆盖的面积,覆盖的面积,且等于且等于在原点的自相关函数值。即时域内的总能量在原点的自相关函数值。即时域内的总能量等于频域内的总能量,人们常用这个关系式等于频域内的总能量,人们常用这个关系式来检验谱分析程序是否正确。来检验谱分析程序是否正确。12/29/2022证明:证明:根据相关定理,进行傅里叶变换反变换有根据相关定理,进行傅里叶变换反变换有 当当 时有时有由能量有限信号的自相关函数定义由能量有限信号的自相关函数定义当当 时有:时有:因此有:因此有:12/29/20226.6.离散相关函数的实现方法离散
21、相关函数的实现方法离散相关函数的实现方法离散相关函数的实现方法离散相关函数有两种实现方法,一种是离散相关函数有两种实现方法,一种是直接计算方法,即根据离散相关的数学定义直接计算方法,即根据离散相关的数学定义直接计算的方法,此法最大的缺点是计算速直接计算的方法,此法最大的缺点是计算速度慢,在进行多段平均计算时将需要较大量度慢,在进行多段平均计算时将需要较大量的时间。第二种方法是用的时间。第二种方法是用FFT快速算法计算快速算法计算来实现,此法最大的优点是计算速度快,来实现,此法最大的优点是计算速度快,在在进行多段平均计算时将需要的时间少得多,进行多段平均计算时将需要的时间少得多,但此法要求但此法
22、要求FFT反变换的精度要高,若用软反变换的精度要高,若用软件实现,至少采用实型数来计算。件实现,至少采用实型数来计算。12/29/2022(1)离散自相关函数的计算方法离散自相关函数的计算方法 对于离散化的时间序列对于离散化的时间序列x(i)i=1,2,N,自相关函数自相关函数的直接计算式为:的直接计算式为:式中式中k为时间滞后点数为时间滞后点数,M为总的时间滞后数。为总的时间滞后数。通常采用软件通常采用软件FFT快速算法计算自相关,其步骤为:快速算法计算自相关,其步骤为:a.对离散序列对离散序列x1,x2,x512补零得到补零得到1024点的新序列:点的新序列:x1,x2,x1024 b.对
23、序列对序列xi作作FFT并取其模得自谱:并取其模得自谱:c.对自谱取多段平均得:对自谱取多段平均得:d.对对 作逆作逆FFT变换取实部并乘以修正系数得相关变换取实部并乘以修正系数得相关函数:函数:12/29/2022(2)离散互相关函数的计算方法离散互相关函数的计算方法 对于离散化的时间序列对于离散化的时间序列x(i)和和y(i)i=1,2,N,互互相关函数的直接计算式为:相关函数的直接计算式为:式中式中k为时间滞后点数为时间滞后点数,M为总的时间滞后数。为总的时间滞后数。通常采用软件通常采用软件FFT快速算法计算互相关,其步骤为:快速算法计算互相关,其步骤为:a.对离散序列对离散序列x(i)
24、和和y(i)i=1,2,512补零得到补零得到1024点的新序列:点的新序列:x(i),y(i)i=1,2,1024 b.对序列对序列x(i),y(i)作作FFT得互谱:得互谱:c.对互谱取多段平均得:对互谱取多段平均得:d.对对 作逆作逆FFT变换取实部并乘以修正系数得互相变换取实部并乘以修正系数得互相关函数:关函数:12/29/20224.2 连续傅立叶变换连续傅立叶变换 傅立叶级数傅立叶级数 傅立叶变换傅立叶变换12/29/2022 1.1.周期信号与傅里叶级数周期信号与傅里叶级数周期信号与傅里叶级数周期信号与傅里叶级数 一个周期信号一个周期信号f(t),当满足狄义赫利当满足狄义赫利(D
25、irichlet)条件,即:周期性函数在一个周期内只条件,即:周期性函数在一个周期内只有有限个极值点,只有有限个第一类间断点有有限个极值点,只有有限个第一类间断点(这种间这种间断点处左、右极限存在断点处左、右极限存在),而其余各处都是连续的,而其余各处都是连续的,这时可以用一个傅里叶级数来表示:这时可以用一个傅里叶级数来表示:一一.傅里叶级数傅里叶级数12/29/202212/29/2022 2.2.傅里叶级数的物理意义傅里叶级数的物理意义傅里叶级数的物理意义傅里叶级数的物理意义(1).任意一个周期函数,只要满足一定条件,任意一个周期函数,只要满足一定条件,都可以分解为基频都可以分解为基频f0
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 动态 信号 处理 第四
限制150内