随机信号分析与处理基础.ppt
- 文档编号:16837918
- 上传时间:2023-07-18
- 格式:PPT
- 页数:88
- 大小:1.96MB
随机信号分析与处理基础.ppt
《随机信号分析与处理基础.ppt》由会员分享,可在线阅读,更多相关《随机信号分析与处理基础.ppt(88页珍藏版)》请在冰点文库上搜索。
第一节随机信号的描述随机信号及其概率结构随机信号在时域的数值特征随机信号的频域描述第二节随机信号通过线性系统的分析平稳随机信号通过连续系统平稳随机信号通过离散系统过渡过程分析,第六章随机信号分析与处理基础,第一节随机信号的描述,1随机信号及其概率结构
(1)随机信号相关的基本概念:
“样本空间”(或称“集合”):
用来表示随机信号的全部可能观测到的波形记录的集合,用X(t)表示;“样本函数”(或称“实现”):
用于表示“样本空间”中每一个确定波形的函数,用x(t)表示;随机信号的样本集合:
X(t)=xi(t),i=1,2,。
t=t1时,随机信号的状态为X(t1)=xi(t1),i=1,2,,为一个数值集合。
随机信号的理解方法:
将其看作随机变量的时间过程,分为连续时间随机信号和离散时间随机信号。
例,汽车车架垂直加速度时间历程记录曲线,图中每一条曲线xi(t)都是加速度时间历程的一次试验记录。
x1(t),x2(t),xn(t)构成加速度时间历程的集合,称为样本空间,记作X(t)。
每一记录曲线称为一个样本,记作xn(t)。
由图可见,各条曲线互不相同,显然不可能用明确的函数式描述。
在任意时刻t1,加速度量值X(t)是一个随机变量。
全部加速度记录的样本空间是无穷多个随机变量的集合。
这种随机现象的进行过程用随机过程来描述。
概率结构:
随机信号是随时间变化的随机变量,用概率结构来描述。
对于离散型随机变量,用概率描述;对于连续型随机变量,用概率密度描述。
1)一维概率分布函数表示随机信号X(t)在t1时刻的取值不大于x1的概率。
2)一维概率密度函数表示随机信号在t1时刻的取值落入x1,x1+极小区间的平均概率。
3)n维联合概率分布函数:
4)n维联合概率密度函数:
在实际中,往往只考虑一维和二维的概率分布函数和概率密度函数。
随机信号分类平稳随机信号非平稳随机信号,平稳信号统计特性不随时间的推移而变化的随机信号。
任意时刻ti的随机变量X(ti)求集合平均有,若不依赖于采样时刻ti,x(ti)为常值,即x(ti)=x,则这种随机信号为平稳信号。
1.随机信号在时域的数字特征,数字特征包括均值(或数学期望)、均方值、方差、相关函数、协方差函数、功率谱等。
连续时间随机信号的数字特征离散时间随机信号的数字特征,连续时间随机信号的数字特征,均值(或数学期望):
随机信号x(t)的所有样本函数在同一时刻取值的统计平均值。
离散随机信号的均值:
连续随机信号的均值:
为随机变量x(t)各个样本的摆动中心。
对于平稳随机信号:
均值为一个与时间无关的常数,相当于信号的直流分量。
均方值:
用于表示随机信号的平均功率,表达式为:
对于平稳随机信号:
其均方值仍为一个与时间无关的常数。
方差:
用于表明随机信号各可能值对其平均值的偏离程度,是随机信号取值分散性的度量。
表达式为:
其中,均方差;对于平稳随机信号:
可见其方差也为一个与时间无关的常数。
相关函数与协方差函数:
1)自相关函数:
用于反映随机信号在不同时刻的内在联系,表达式为:
当t1=t2=t时,有:
表明:
随机信号的均方值是它的自相关函数在t1=t2时的特例。
对于平稳随机信号:
可见其自相关函数是时间间隔的函数。
2)自协方差函数:
用随机信号X(t)在两个不同时刻t1、t2取值起伏变化的相依程度来描述随机信号不同时刻的关联关系,表示为:
当t1=t2=t时,有:
对于平稳随机信号:
可见其自协方差函数也是时间间隔的函数。
方差,3)互相关函数:
用于研究两个随机信号x(t)和y(t)的相互关系,表达式为:
其中,表示两随机信号的二维联合概率密度函数;表示两随机信号之间的线性依赖关系。
对于平稳随机信号,当满足下面条件时:
有:
4)互协方差函数:
用随机信号X(t)在两个不同时刻t1、t2取值起伏变化的相依程度来描述随机信号不同时刻的关联关系,表示为:
对于平稳随机信号:
可见互协方差函数也是时间间隔的函数。
平稳随机信号的均值、方差、均方值是与时间无关的常量,相关函数及协方差仅是时移的函数,与随机信号的起止时刻t无关。
平稳随机信号最重要的特点是随机信号在不同时刻具有相同的统计特征。
与平稳随机信号相反,非平稳随机信号的统计特性是随着时间的推移而变化的。
平稳随机信号的重要性质,例:
一个随机信号,其中、均为常数,为区间均匀分布的随机变量,求该随机信号x(t)的均值、均方值、方差、自相关函数及自协方差函数。
解:
随机变量在区间均匀分布,它与时间无关,故其一维、二维概率密度均为:
根据式(6-15)求得该随机信号的均值为:
根据式(6-20)求得该随机信号的自相关函数为:
其中,根据式(6-24)求得其自协方差函数为:
根据式(6-25)求得该随机信号的方差为:
根据式(6-26)求得该随机信号的均方值为:
时间平均表征量:
1)随机信号x(t)的时间均值:
2)随机信号x(t)的时间相关函数:
各态遍历性随机信号:
在一定条件下,平稳随机信号的每一个样本都同样地经历了随机信号其它样本的各种可能状态,因而从一个样本的统计特性(时间平均)就能得到全部样本的统计特性(集平均),此类信号称为各态遍历性随机信号。
各态遍历性随机信号的数字特征:
先通过简单的实验方法或数学方法得到一个各态遍历性随机信号的均值、自相关函数,再运用关系式得到其它数字特征量。
各态遍历性随机信号及其数字特征,各态遍历性随机信号的意义,可以用时间充分长的单个样本函数的时间平均统计参数来代替总体的平均统计值,这给试验信号处理带来了极大的方便。
离散时间随机信号:
时间t的取值是离散的随机信号x(t)。
数字特征:
与连续时间情况并无不同,只是时间变量t应变为限取整数的变量n。
如,此时总集均值EX(t)应表示为:
遍历性随机序列:
对于一个平稳随机序列X(n),若其各种时间平均以概率1收敛于相应的集合平均,则称其为遍历性随机序列。
离散时间随机信号的数字特征,
(1)连续时间情况功率谱(或称功率谱密度函数):
设xi(t)是随机信号x(t)的一个样本,不满足傅立叶变换所要求的平方可积条件,故将其截短,形成,即:
将上式进行傅立叶变换,结合帕斯瓦尔得平均功率谱表达式为:
随机信号的功率谱:
随机信号的频域描述,维纳-辛钦(Wiener-Khinchine)定理:
平稳随机信号x(t)的功率谱是它的自相关函数的傅立叶变换;而x(t)的自相关函数是其功率谱的傅立叶反变换。
用公式表述如下:
互谱(或称互功率密度谱):
对于两个随机信号X(t)、Y(t),若它们是平稳相关的,则其互谱为:
互谱与互相关函数的关系:
为一组傅立叶变换对,满足:
例:
求随机相位余弦信号的功率谱及平均功率。
解:
设随机相位余弦信号为,其中、均为常数,为区间均匀分布的随机变量。
由例6-1得自相关函数为:
则根据维纳-辛钦(Wiener-Khinchine)定理得此功率谱为:
平均功率为:
(2)离散时间情况上述维纳-辛钦(Wiener-Khinchine)定理同样适用于离散时间序列,故可以求平稳随机序列的功率密度谱,以及序列的自相关函数。
离散时间信号功率谱的特点:
1)功率谱是周期性的,因此可作傅立叶级数分解;2)反演变换的积分区间是。
第二节随机信号通过线性系统的分析,随机信号通过线性系统时,可能出现的两种情况:
(1)平稳情况(稳态):
如果输入是平稳随机信号,系统是线性时不变且稳定的,则当系统完成过渡过程进入稳态后,输出也应该是平稳的随机信号;
(2)非平稳情况(暂态):
即系统进入稳态前的过渡过程。
1平稳随机信号通过连续系统
(1)系统响应的时域分析设一线性时不变系统的单位冲激响应为h(t),当时,其输出零状态响应Y(t)的表达式为:
1.输出Y(t)的均值:
2.输出的Y(t)自相关函数:
3.输入与输出之间的互相关函数:
或,
(2)系统响应的频域分析当系统的输入、输出均为平稳随机信号时,可以通过维纳-辛钦(Wiener-Khinchine)定理求取功率谱密度函数。
1.系统输出的功率谱密度:
可见,系统的功率谱传输能力仅与系统的幅频特性有关,而与系统的相频特性无关。
2.系统输入与输出之间的互功率谱密度:
可见,互功率谱密度不仅包含有系统幅频特性函数的幅度信息,还包含有相位信息。
2平稳随机信号通过离散系统随机序列通过离散系统的分析,与连续时间随机信号通过连续系统对输出统计特性的分析计算类似,可以采用时域分析和频域分析两种方法。
具体分析过程同6.2.1类似,详见p260262。
6.2.3过渡过程分析
(1)连续时间信号情况1.零输入响应:
由初始状态决定。
若初始状态是确定性量,则所得响应也是确定性的;若初始状态是随机变量,则通过解析表达式,分析所得结果的统计特征。
例6-6:
p262263,2.零状态响应:
通过卷积来求解。
对于因果系统,设x(t)在t=0时刻接入,则系统输出为:
1)输出y(t)的自相关函数:
2)输入输出的互相关函数:
注:
分别分析零输入响应和零状态响应,再将二者结合可求得系统过渡过程的总响应。
(2)离散时间信号情况1.分析方法:
1)根据初始状态是确定性量还是随机变量,确定系统的零输入响应过程的确定性或随机性;2)通过卷积和分析系统的零状态响应。
2.例6-7:
p265,功率谱分析!
相干函数功率谱分析的应用细化分析倒频谱分析,功率有限信号功率谱能量有限信号能量谱,频谱分析可以研究信号的频率结构工程信号往往是复杂的:
包含周期成分包含随机干扰频谱表现为连续尖脉冲不易识别工程信号分析的关键是降低噪声,提高信噪比频谱分析不改变信噪比功率谱分析,功率谱分析,设工程信号x(t)=s(t)+n(t)其中s(t)为有用信号;n(t)为随机干扰则x(t)的傅里叶变换X()为:
X()=S()+N()可见:
傅里叶变换不会提高信噪比。
对x(t)进行自相关分析Rxx()=Rss()+Rnn()+Rsn()+Rns()Rss()(n(t)零均值,)可见:
相关函数可以提高信噪比,但不反映频谱相关函数的傅里叶变换功率谱,可以提高信噪比,又能反映频率结构,能量谱与功率谱,设能量有限信号x(t)、y(t)相关函数Rxx()、Rxy()定义:
自能量谱密度函数,互能量谱密度函数,能量谱从频域提取信号中的周期分量或同频分量相关函数从时域提取信号中的周期或同频分量,设功率有限信号x(t)、y(t)相关函数Rxx()、Rxy()定义:
自功率谱密度函数,互功率谱密度函数,功率谱从频域提取信号中的周期分量或同频分量相关函数从时域提取信号中的周期或同频分量,自能量谱密度函数,互能量谱密度函数,自功率谱密度函数,互功率谱密度函数,能量谱与功率谱统称为功率谱,功率谱的性质,功率谱密度函数的性质自功率谱Sxx(f)是实偶函数;互功率谱Sxy(f)是非奇非偶复函数;双边谱:
f(-,);功率谱与相关函数包含的信息完全等价。
单边功率谱工程上,负频率没有意义定义f0,)的单边功率谱根据能量守恒准则:
功率谱的性质,单边功率谱自功率谱Gxx(f),Sxx(f),f,互功率谱Gxy(f),功率谱的性质,功率谱的物理意义由功率谱的定义:
Gxx(f)下的面积等于信号的总能量故称Gxx(f)为能量有限信号的能量谱密度函数或功率有限信号的功率谱密度函数Gxx(f)任意频段间的面积该频带下信号的能量,功率谱的计算,工程测试中,只能得到有限区间上的信号x(t):
t(-T/2,T/2);x(t)X()x(t)x(n)根据连续信号的帕斯瓦尔公式:
功率谱的计算,双边互功率谱:
双边自功率谱:
单边自功率谱:
单边自功率谱:
功率谱的计算,得到功率谱计算2种方法:
直接根据相关函数求傅里叶变换先求幅值谱,再利用上述公式计算,相干函数,评价因果系统输入输出之间的因果关系输出是否由输入引起?
输出中有多少由输入引起?
系统的线性特性如何?
相干函数与相干分析,相干函数,Rxy()能从延时域上描述输出与输入的相关关系相干函数则从频域上描述输出与输入的相关关系由:
相干函数,Gxy(f)表示y与x的相关程度Gxy(f)越大,y与x的相关性越强,相干函数定义:
相干函数,相干函数的意义,:
x(t)ky(t),完全相干,:
x(t)与y(t)完全不相关,完全不相干,:
x(t)与y(t)部分相关,部分相干,完全相干:
输出100由输入引起,理想系统;,完全不相干:
输出与输入毫无关系,极端情况;,部分相干:
一般系统,表示系统多输入系统,有干扰,非线性,功率谱分析的应用,频谱分析故障诊断信号检测模式识别,频率响应函数的测定,已知系统的输入、输出,求频响函数,
(1)幅值估计法,无干扰线性系统,只能求出幅频特性无相位信息,频率响应函数的测定,
(2)互谱估计法,无干扰线性系统,既能求出幅频特性又能得到相位信特性,频率响应函数的测定,(3)实际物理系统频响的求法,有干扰系统,干扰n,m不相干,n、m与a、b不相干,a、b系统真正输入、输出,x、y实际测得的输入、输出,系统频率响应由a、b确定:
但a、b不可知,可通过以下近似方法求得,频率响应函数的测定,幅值估计法,幅频误差相频丢失,频率响应函数的测定,互谱估计法,幅频测小相频精确,频率响应函数的测定,互谱估计法,幅频测大相频精确,频率响应函数的测定,平均估计法,相频精确,幅频测小,幅频测大,幅频更接近,细化分析,引出:
功率谱分析基本思想是带通滤波FFT:
分析个别频带上的频率成分分辨率确定、采样点确定一些工程信号频带很宽频率谱线密集为了识别谱图的细微结构,必须提高功率谱分析的分辨率分析仪器具有很宽的频率范围,细化分析,信号时域宽度T,采样点数N,采样间隔,信号频带宽度fc=fs/2,谱线数N,频率分辨率,提高频率分辨率的两种途径:
引起混叠,工作量增加;硬件不允许,均不可行!
措施:
保持N不变,设法降低,细化分析(ZOOMFFT)的基本思想,细化分析的基本思想,移频低通滤波重新采样FFT,感兴趣频段,-f10f1,移频滤波后,感兴趣频带的最高频率f1f00.5fm;有可能在满足采样定理的前提下增大采样间隔;提高频率分辨率,x(t)e-j2f0Gxx(f+f0),移频,确定分析频带f0f1,找出中心频率f0;利用移频技术,将f0移至0Hz处:
由,频带左移f0,感兴趣的频带移至低频端。
移频的实质是:
给时域信号x(t)乘以,或给时域信号x(n)乘以,其中:
或,Gxx(f),-f10f1,-f10f1,低通滤波,以截至频率为f0的低通滤波器滤波得到感兴趣的频段最高频率f1低通滤波的目的:
避免重新采样时发生混叠,f/Hz,Gxx(f+f0),-f10f1,D为正整数细化倍数,重新采样,信号原长为T,采样间隔Ts,采样点数N,则,采样频率fs,频率分辨率,重新采样:
采样间隔DTs,采样点数N,则:
信号长度变为T=DTsN=DT,,重新采样频率,频率分辨率,结论:
只要将信号长度增加D倍,采样间隔加大D倍,采样点数不变时,频率分辨率提高D倍。
快速傅里叶变换,对重新采样的数据进行FFTX(k);频谱重构,得到细化D倍的频谱:
ZOOMFFT的一般步骤,倒频谱分析,或倒功率谱分析Cepstrum,功率谱分析(Spectrum)的思想:
FFT+带通滤波:
分析全频域频率结构功率谱分析(Spectrum)的局限性:
仅适应于线性叠加信号的频谱分析两信号频带不交叠时信号的分离不适用于非线性信号处理,x(t)+y(t),x(t)y(t),信号非线性叠加的例子,测量问题:
语音分析问题:
输入x(t),输出y(t),h(t),H(f),Y(f),X(f),H(f),X(f),Y(f),X(f)与H(f)不可分离,声道,声门,声门冲激,发音器官的特征,语音,声门冲激信号,X(f)与H(f)不可分离,解决问题的思路,频域上“乘积”化为“求和”频谱取对数倒频谱/倒功率谱1962年Tukey等认提出用于语音分析、地震分析、回声分析、故障诊断等领域,倒频谱分析原理,设时域信号,离散傅里叶变换,频谱的模取对数,傅里叶反变换,实倒频谱定义:
信号频谱模的对数的傅里叶逆变换,幅频,化积为和,Cx()只反映x(n)的实部,称为实倒频谱,由实倒频谱不能恢复信号的频谱!
指数运算,倒频谱自变量的量纲,时间,频率,倒时间,频率,倒时间,幅度变化,时间,倒频率,幅度变化,具有时间的量纲,但幅度发生变化称为倒频率,C()称为倒频谱,倒频谱分析原理,一般情况下,X(n)为复数:
C()包含X(n)的实部和虚部,可以恢复信号的频谱,C(n)称为信号的复倒频谱,复倒频谱与信号频谱一一对应,一般在不关心相位信息时,采用实倒谱,倒功率谱,设时域信号,功率谱(实谱),倒功率谱,从倒功率谱可以恢复信号的功率谱!
倒功率谱定义:
信号功率谱对数的傅里叶变换,实偶函数,倒功率谱分析的过程,编辑,倒滤波,倒频谱术语,频谱:
spectrum倒频谱:
cepstrum频率:
frequency倒频率:
quefrency相位:
phase倒相位:
saphe滤波:
filter倒滤波:
lifter,倒频谱分析的应用,传递函数的分离回声分析与剔除,
(1)传递函数分离,输入x(t),输出y(t),H(f),Y(f),X(f),H(f),X(f),Y(f),求H(f)?
h(t),Cx(),Ch(),傅里叶变换FFT:
(1)传递函数分离,Cx和Ch在倒频域上占有不同的频段:
Cx在高频段Ch在低频段倒滤波(倒谱编辑)低通滤波得Ch分离传递函数h(t)高通滤波得Cx分离信号x(t),Cx(),Ch(),指数运算exp():
开平方:
傅里叶反变换IFFT:
傅里叶变换FFT:
指数运算exp():
开平方:
傅里叶反变换IFFT:
实例,某检测系统输出y(t)功率谱对数图,输入x(t),系统脉冲响应函数h(t),lnSy(f),f0,ln|H(f)|2,lnSx(f),Cx(),Ch(),对数频域有周期分量,Sy(f)=|H(f)|2SX(f);倒频域必有脉冲成分,对应于Sx(f);f0与0满足:
0f01;滤除高频,可恢复脉冲相应函数h(t);滤除低频,可恢复被测信号x(t),
(2)回声分析与剔除,噪声测量中会引入回声;回声使声源的功率谱产生畸变;影响声源定位或频率识别;精密测量时应予以剔除。
声源信号x(t);声速V;反射系数;声源至反射壁距离S。
则信号通道时差02S/V,反射壁,声源,x(t),y(t),波形分析,时域:
传声器输出原声回声回声相对原声:
幅度减小、延时时域中波形重叠,不能剔除回声,波形分析,频域:
幅频:
原声谱回声谱,波形完全重叠相频:
原声相位与回声相位叠加频域中,不能剔除回声,对策:
倒频谱分析,时域卷积,倒频分离,回声分析,测量的时域表达式:
回声分析,两边取对数,幂级数展开,回声分析,求倒功率谱:
傅里叶反变换,回声分析,结论:
回声的倒谱是一系列周期为0的冲激信号;随|n|的增大,脉冲幅度快速减小.冲激只发生在0的整数倍处(n0);鸡冠状滤波器可以滤除此周期冲激;由Cx()恢复声源真实信号的倒谱;有多个回声实,倒谱中有多个冲激族。
回声剔除,倒谱编辑倒滤波,鸡冠状滤波器剔除回声,声源信号恢复,真实信号得以恢复!
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 随机 信号 分析 处理 基础
![提示](https://static.bingdoc.com/images/bang_tan.gif)