清华大学测试与检测技术基础-王伯雄-第3章相关分析.ppt

清华大学测试与检测技术基础-王伯雄-第3章相关分析.ppt

《清华大学测试与检测技术基础-王伯雄-第3章相关分析.ppt》由会员分享，可在线阅读，更多相关《清华大学测试与检测技术基础-王伯雄-第3章相关分析.ppt（70页珍藏版）》请在冰豆网上搜索。

2.2.7.3相关分析,一、相关二、互相关函数与自相关函数三、相关函数的工程意义及应用,一、相关（correlation）,相关：

用来描述一个随机过程自身在不同时刻的状态间，或者两个随机过程在某个时刻状态间线性依从关系的数字特征。

图2.52变量x和y的相关性（a）精确相关（b）中等程度相关（c）不相关,评价变量x和y间线性相关程度的经典方法：

协方差xy：

式中，E表示数学期望值；x=Ex为随机变量x的均值；y=Ey为随机变量y的均值；相关函数xy：

式中x、y分别为x、y的标准偏差，而x和y的方差x2和y2则分别为,（2.142）,（2.143）,（2.144）,（2.145）,利用柯西许瓦兹不等式（Cauchy-Schwarzinequality）可知|xy|1。

当xy=1时，所有数据点均落在y-y=m（x-x）的直线上，因此x,y两变量是理想的线性相关。

当xy=0时，（xi-x）与（yi-y）的正积之和等于其负积之和，因而其平均积xy为0，表示x,y之间完全不相关。

（2.146）,二、互相关函数与自相关函数,对于各态历经过程，可定义时间变量x（t）和y（t）的互协方差（cross-covariance）函数为式中称x（t）与y（t）的互相关（cross-correlation）函数，自变量称为时移。

（2.147）,（2.148）,当y（t）x（t）时，得自协方差（auto-covariance）函数其中称为x（t）的自相关（auto-correlation）函数。

周期函数的自相关函数仍为周期函数，且两者的频率相同，但丢掉了相角信息。

同频相关，不同频不相关。

（2.149）,（2.150）,图2.53典型的自相关函数和互相关函数曲线（a）自相关函数（b）互相关函数,例1求正弦函数x（t）=Asin（t+）的自相关函数。

解：

正弦函数x（t）是一个均值为零的各态历经随机过程，其各种平均值可用一个周期内的平均值来表示。

令t+=,则dt=d/,由此得正弦函数的自相关函数是一个与原函数具有相同频率的余弦函数，它保留了原信号的幅值和频率信息，但失去了原信号的相位信息。

自相关函数可用来检测淹没在随机信号中的周期分量。

自相关和互相关函数的估计和具有限个数据点N的相关函数估计的数字处理表达式则为：

（2.160）,（2.161）,（2.162）,（2.163）,三、相关函数的工程意义及应用,不同类别信号的辨识,图2.55典型信号的自相关函数,相关滤波（filteringbycorrelation）,图4.79相关滤波频谱分析仪原理框图,相关测速和测距,图2.56相关法测量声传播距离,图2.57带钢测速系统,测量流速和流量,图2.58相在法测定流量,2.2.7.4功率谱分析,2.2.7.4功率谱（powerspectrum）分析,一、自功率谱密度函数二、巴塞伐尔（Parseval）定理三、互功率谱密度函数四、自谱和互谱的估计五、工程应用,一、自功率谱密度函数,设x（t）为一零均值的随机过程，且x（t）中无周期性分量，则其自相关函数Rx（）在当时有该自相关函数Rx（）满足傅里叶变换的条件。

对作傅里叶变换可得其逆变换为,（2.167）,（2.168）,Sx（f）为x（t）的自功率谱密度函数（autopowerspectrum），简称自谱或功率谱。

功率谱Sx（f）与自相关函数Rx（）之间是傅里叶变换对的关系，亦即式（2.167）和（2.168）称为维纳辛钦（Wiener-Khintchine）公式。

由于Rx（）为实偶函数，因此亦为Sx（f）实偶函数。

图2.59单边功率谱和双边功率谱,当=0时，根据自相关函数Rx（）和自功率谱密度函数Sx（f）的定义，可得Sx（f）曲线下面和频率轴所包围的面积即为信号的平均功率；Sx（f）就是信号的功率谱密度沿频率轴的分布，故也称为功率谱。

（2.169）,二、巴塞伐尔（Parseval）定理,设有变换对：

按频域卷积定理有令k=0，有又令h（t）=x（t），得,x（t）为实函数，故X（-f）=X*（f），于是有巴塞伐尔定理：

信号在时域中计算的总能量等于它在频域中计算的总能量。

式（2.170）又称信号能量等式。

|X（f）|2称能量谱，它是沿频率轴的能量分布密度。

在整个时间轴上信号的平均功率可计算为自谱密度函数与幅值谱之间的关系为,（2.170）,（2.171）,（2.172）,对于单边功率谱G（f）也应满足巴塞伐尔定理，故有由此规定Gx（f）的图形如图2.59中所示。

（2.173）,图2.59单边功率谱和双边功率谱,根据信号功率（或能量）在频域中的分布情况，将随机过程区分为窄带随机、宽带随机和白噪声等几种类型。

窄带过程的功率谱（或能量）集中于某一中心频率附近，宽带过程的能量则分布在较宽的频率上，而白噪声过程的能量在所分析的频域内呈均匀分布状态。

三、互功率谱密度函数,若互相关函数Rxy（）满足傅里叶变换的条件，则定义Rxy（）的傅里叶变换为信号x（t）和y（t）的互功率谱密度函数，简称互谱密度函数（crosspowerspectrum）或互谱。

根据维纳辛钦关系，互谱与互相关函数也是一个傅里叶变换对，即因此Sxy（f）的傅里叶逆变换为：

（2.175）,（2.176）,定义信号x（t）和y（t）的互功率为因此互谱和幅值谱的关系为正如Ryx（）Rxy（）一样，当x和y的顺序调换时，Syx（）Sxy（）。

但根据Rxy（-）=Ryx（）及维纳辛钦关系式，不难证明：

其中,（2.177）,（2.178）,（2.179）,Sxy（f）也是含正、负频率的双边互谱，实用中也常取只含非负频率的单边互谱Gxy（f），由此规定自谱是f的实函数，而互谱则为f的复函数，实部Cxy（f）称为共谱（cospectrum），虚部Qxy（f）称为重谱（quadspectrum），即写为幅频和相频的形式：

（2.180）,（2.181）,（2.182）,四、自谱和互谱的估计,定义功率谱亦即自谱的估计值互谱的估计为,（2.183）,（2.184）,（2.185）,五、工程应用,求取系统的频响（frequencyresponse）函数线性系统的传递函数H（s）或频响函数H（j）十分重要，在机器故障诊断等多个领域常要用到它。

例1：

机器由于其轴承的缺陷而在机器运行中会造成冲击脉冲信号，此时若用安装在机壳外部的加速度传感器来接收时，必须考虑机壳的传递函数。

例2：

当信号经过一个复杂系统被传输时，系统各环节的传递函数便必须要加以考虑。

一个线性系统的输出y（t）等于其输入x（t）和系统的脉冲响应h（t）的卷积，即根据卷积定理，上式在频域中化为式中H（f）即为系统的频响函数。

（2.192）,（2.193）,通过自谱和互谱来求取H（f）：

对式（2.193）两端乘以各自的复共轭并取期望值有上式反映出输入与输出的功率谱密度和频响函数间的关系；式中没有频响函数的相位信息，因此不可能得到系统的相频特性。

（2.194）,如果在式（2.193）两端乘以x（f）的复共轭并取期望值，则有由于Sx（f）为实偶函数，因此频响函数的相位变化完全取决于互谱密度函数的相位变化。

式（2.195）将输入、输出的相位关系完全保留了下来，且在这里输入的形式并不一定限制为确定性信号，也可以是随机信号。

（2.195）,通常一个测试系统往往受到内部和外部噪声的干扰。

从而输出也会带入干扰。

输入信号与噪声是独立无关的，因此它们的互相关为零。

结论：

在用互谱和自谱求取系统频响函数时不会受到系统干扰的影响。

旋转机械振动特性检测旋转机械的转轴部件从起动、升速到额定转速的过程共经历了全部转速的变化，因此在各个转速下的振动状态可用来对机器的临界转速、固有频率和阻尼比等各参数进行辨识。

起动和停车过程则包含了丰富的信息。

是常规运行状态下所无法获得的。

“瀑布图（waterfallplot）法”：

在机械振动或停车过程中将不同转速下振动的功率谱图迭加而形成的一种图。

图2.61旋转机械的瀑布图,由图可见机器的回转频率n（r/min）及其各次谐波下谱峰高度，由此来得出机器的临界转速、固有频率及阻尼比等数据。

从图可见，机器临界转速约为4000r/min，机器振动的高次谐波分量很小，主要是回转频率处的谱峰，因此可判断转子存在有较严重的失衡。

此外还可看到图中频率60HZ处有一谱峰值，它不随转速升高而改变，判断为电源的脉动干扰。

2.3数字信号处理,数字信号处理（digitalsignalprocessing）：

利用计算机或专用信号处理设备，以数值计算的方法对信号作采集、变换、综合、估值与识别等处理。

一、离散傅里叶变换（DFT）二、离散傅里叶变换的性质三、采样定理四、泄漏与加窗处理五、栅栏效应六、快速傅里叶变换（FFT）,一、离散傅里叶变换（DFT）,对于一个非周期的连续时间信号x（t）来说，它的傅里叶变换应该是一个连续的频谱X（f），其运算公式根据第二章的内容有,（2.199）,（2.200）,图2.63傅里叶变换的几种类型,对于无限连续信号的傅里叶变换共有四种情况：

对于非周期连续信号X（t），频谱X（f）是连续谱；对于周期连续信号，傅里叶变换转变为傅里叶级数，因而其频谱是离散的；对于非周期离散信号，其傅里叶变换是一个周期性的连续频谱；对于周期离散的时间序列，其频谱也是周期离散的。

结论：

若x（t）是周期的，频域中X（f）必然是离散的，反之亦然。

若x（t）是非周期的，则X（f）一定是连续的，反之亦然。

第四种亦即时域和频域都是离散的信号，且都是周期的，给我们利用计算机实施频谱分析提供了一种可能性。

对这种信号的傅里叶变换，我们只需取其时域上一个周期（N个采样点）和频域一个周期（同样为N个采样点）进行分析，便可了解该信号的全部过程。

DFT的定义：

对有限长度的离散时域或频域信号序列进行傅里叶变换或逆变换，得到同样为有限长度的离散频域或时域信号序列的方法，便称为离散傅里叶变换（digitalFouriertransform,DFT）或其逆变换（IDFT）。

离散傅里叶变换的公式：

式中x（n）和X（k）分别为和的一个周期，此处将t和f0均归一化为1。

（2.205）,（2.206）,离散傅里叶变换意义：

可以对任意连续的时域信号进行采样和截断并对其作离散傅里叶变换的运算，得到离散的频谱，该频谱的包络即是对原连续信号真正频谱的估计。

离散傅里叶变换的过程：

时域采样（samplingint-domain）；时域截断（truncationint-doman）；频域采样（samplinginf-domain）。

图2.64离散傅里叶变换的图解过程

（一）,图2.64离散傅里叶变换的图解过程

（二）,图2.64离散傅里叶变换的图解过程（三）,三、采样定理（samplingtheorem）,混叠（aliasing）：

若采样率过低即采样间隔大，则系列的离散时间序列可能不能真正反映原始信号的波形特征，在频域处理时会出现频率混淆。

图2.65不同采样率对采样信号产生的影响

（一）,图2.65不同采样率对采样信号产生的影响

（二）,采样定理：

为避免混叠产生，要求的采样频率fs必须高于信号频率成分中最高频率fmax的两倍，即乃奎斯特（Nyquist）频率：

在给定的采样频率fs条件下，信号中能被分辨的最高频率。

只有低于乃奎斯特频率的频率成分才能被精确地采样，亦即为避免频率混淆，应使被分析信号的最高频率fmax低于乃奎斯特频率。

（2.222）,（2.223）,图2.66混叠产生的条件,四、泄漏（leakage）与加窗（windowing）,图2.67余弦信号加窗截断造成的泄漏现象,抑制或减小泄漏效应的方法：

选择性能更好的特殊窗来替代矩形窗，亦即加窗处理。

评价窗函数的性能指标：

3dB带宽B：

它是主瓣归一化的幅值下降至-3dB时的带宽。

归一化|W（f）|=20lg|W（f）/W（0）|，带宽B的单位为或f。

旁瓣幅度A