第五章信号处理初步

1、 第五章第五章 重点：相关分析及其应用、功率谱分 析及其应用； 难点：信号数字化中出现的问题：采 样、混叠、截断、泄漏和窗函 数等概念。 测试工作的目的是获得有用的信息，但测试信号中既含有有用信息，也含有大量干扰噪声。因此，信号处理的任务：修正测试系统的某些误差；滤除干扰噪声，提取有用信息；分离信、噪，提高信噪比。 研究信号的构成和特征值（对信号结构无影响）。 把信号经过必要的变换，以获得所需信息的过程（有可能改变信号本身结构）。二者没明确的界限，有时作同义语。 信号分析信号分析：信号处理信号处理： 信号处理可用模拟信号处理系统和数字信号系统来实现。 模拟信号处理系统：由实现模拟运算功能的 电

2、路组成。 数字信号处理系统：由微型计算机和相关软 件组成（也可用专用信 号处理机）。 第一节 随着微电子技术和信号处理技术的发展，在工程测试中，数字信号处理方法得到广泛的应用。从传感器获取的信号大多为模拟信号，进行数字信号处理之前，一般先要对信号作预处理和数字化处理。而数字传感器可直接通过接口与计算机联接，将数字信号送给计算机（或数字信号处理器）进行处理。基本步骤：数字信号处理系统如下图：数字信号处理系统如下图：物理信号物理信号对象对象传传感感器器电信号电信号预预处处理理电信号电信号A/D转换转换数字信号数字信号计计算算机机显显示示D/A转换转换电信号电信号控制控制物理信号物理信号 在数字处理

3、之前，对信号用模拟方法进行的处理。 如对输入信号的幅值处理，使幅值与A / D转换器的动态范围相适应；衰减信号中不感兴趣的高频部分，减少频混的影响；隔离被分析信号中的直流分量，消除趋势项及直流分量。1 1、预处理：、预处理： A / D转换转换：将预处理以后的模拟信号变为数 字信号，存入指定的地方。 分析计算分析计算：对采集到的数字信号进行分析、 计算（可用数字运算器件组成的 信号处理器，也可用通用计算机）。4. 4. 结果显示结果显示：一般采用数据和图形显示结果。 数字信号处理首先把一个连续变化的模拟信号转化成数字信号，然后由计算机处理，从中提取有用的信息。 众所周知，数字计算机不能对一个连

4、续的模拟信号计算，其原因为： 必须经过A / D转换器把一个模拟信号转变成数字量存放到一定的内存单元，然后进行计算； 2. 计算机内存容量总是有限的，它不可能存放无限 多的采样数据。因此“数值离散”和“点数有限”是 使用数字计算机进行傅立叶变换的两大特点。为 了区别常见的傅立叶变换，称之为有限离散傅立 叶变换，简称离散傅立叶变换（DFT）。 下面从图形上认识这一过程及出现的问题:某一连续信号x（t）及傅立叶变换X（f）； 为了用计算机计算，必须将x（t）变换成有限长的离散序列。采样：采样：如图53为采样函数s（t）及频谱S（f）。 由图可知，时域函数的离散时域函数的离散导致频域图形的频域图形的

5、周期化周期化。这是离散傅立叶变换的第一次误差。时域截断时域截断：取出N个有限点，这在数学上可理解为把采样后的信号乘以一个矩形窗函数（图55），窗宽为T 。那么 。 sTNT 如图56所示。可见时域函数的截断截断导致频域函数出现皱波（泄漏）皱波（泄漏），这是离散傅立叶变换引入的第二次误差 。 图56中的频域函数仍不是计算机能接受的离散信号。 至此，得到了如图58所示的离散傅立叶变换对，它在时域上和频域上都是用离散值表示的，且在时域、频域都周期化。但计算机中仅存储N个采样点，分别表示时域、频域波形的一个周期，可用此近似原来的连续傅立叶变换。那么，必须要做好该信号处理过程中采样、截断、DFT计算，否

6、则均会引起失真和误差。 将连续变化的模拟信号离散化的过程。 采样过程可以看作用等间隔的单位脉冲序列去乘模拟信号。这样各采样点上信号的大小就变成脉冲序列的强度，这些强度值被量化而成为相应的数值。 ：X(0), X(1), X(2), , X(n) 长度为T的连续时间信号x（t），从t=0点开始采样，采样得到的离散时间序列x（n）为 1210 N,nf/nxnTxnxss snTts|txnTx 式中Ts采样间隔；N序列长度，N=T/ Ts ；fs 采样频率， fs =1/ Ts ； 采样间隔太小（采样频率 ），数字序列长，计算工作量大，计算速度慢； 若采样间隔过大（采样频率），则可能丢掉有用的信

7、息。1、混叠现象：混叠现象：混叠现象的解释及其避免的办法：混叠现象的解释及其避免的办法： nsnTtts rssTrfTfS 1 tstx rssrssTrfXTTrfTfXfSfX11 采样函数：频域解释频域解释 0t0f0t0ft00f 采样间隔不当引起频率混叠。2、避免混叠的方法：避免混叠的方法：混叠必定出现在f= fs /2（折叠频率）左右两侧的频率处。可以证明，任何一个高频成分f1和低频率成分f2、 fs的关系为： 2221/f/ffs 不产生频率混叠的条件：不产生频率混叠的条件：1、使被采样的模拟信号x（t）成为有限带宽的信号。2、应使采样频率fs大于带限信号的最高频率fh的2倍（

8、采样定理），即hssfTf21 采样定理：采样定理：采样的基本法则 fs 2 fh实际工作中选fc ：滤波器的截止频率。 csff43 若x（t）不是带限信号，则A/D采样前抗混叠滤波处理： 物理信号物理信号对象对象传传感感器器电信号电信号放放大大调调制制电信号电信号A/D转换转换数字信号数字信号展开展开低通滤波低通滤波放大放大 把采样点上信号的瞬时幅值转换成二进制的数字量。即离散信号变成离散数字信号。 在工作中，采用A / D转换器来完成，A / D的位数是一定的，且不论它的位数有多高，模拟信号采样点的幅值电平总落在两个相邻量化电平之间，就要舍入到相近的一个量化电平上，见上图，由于有舍入，就

9、会产生误差。 量化电平与信号实际电平之间的差 ，( ) n其最大值为 。 2x以后的讨论假设 ( )n为0。 信号的历程是无限的，而不可能对无限长的信号进行处理，所以要截断。 将无限长的信号乘以有限 宽的窗函数。 为便于数学处理，通常对截断的信号做周期延拓，得到虚拟的无限长的信号。 设有余弦信号x(t), 用矩形窗函数w(t)与其相乘，得到截断信号：y(t) =x(t)w(t) 将截断信号谱将截断信号谱 X XT T()()与与原始信号谱原始信号谱X()X()相比较相比较可知，它已不是原来的两可知，它已不是原来的两条谱线，而是两段振荡的条谱线，而是两段振荡的连续谱连续谱. . 原来集中在原来集

10、中在f0f0处处的能量被分散到两个较宽的能量被分散到两个较宽的频带中去了，这种现象的频带中去了，这种现象称之为频谱能量泄漏。称之为频谱能量泄漏。 样本截断引起泄漏样本截断引起泄漏。 泄漏与窗函数(Function of the window)频谱的旁瓣有关，窗函数旁瓣，相应的泄漏误差也将减小。所以选择合适的窗函数也是信号处理中的重要问题之一。 经过、截断后，其频谱在频域是连续的（图56），所以必须使频率离散化，实行 时域采样的结果使时域信号平移至各脉冲坐标位置重新构图时域周期延拓。 ：不论时域或频域 采样，有如透过栅栏取值，总会有落 在栅栏之外的。 频率采样间隔 f 是频率分辨力的指标： 1s

11、ffNTsf 采样频率 N 采样点数 若按整周期截断信号，则延拓后的信号将和原信号完全重合。 选择窗函数应使频谱的主瓣宽度窄、旁瓣幅度小。窄的主瓣可提高频率的分辨力，小的旁瓣可以减小泄漏。 主瓣最窄2/T，旁瓣较高T，泄漏较大。 适合于要获得精确主峰的频率，而对 幅值精度要求不高的场合。 特点：特点：主瓣宽度约为矩形窗的2倍，但旁瓣低且不会出现负值。 20221TtTttTtwT 2sin222sin222fTcTfTfTTfWT 旁瓣明显降低，有抑制泄漏的作用，但主瓣较宽，致使频率分辨能力较差。 在截断随机信号或非整周期截断周期函数时，为了平滑或削弱截取信号的两端，减小泄漏，宜加汉宁窗。 （

12、四）、指数窗（四）、指数窗 无旁瓣，主瓣很宽，其频率分辨力低。 对脉冲响应类信号宜加指数窗，若适 当选择衰减函数，可起到抑制噪声的 作用。 常用窗函数常用窗函数 在测试技术中，无论分析两个随机变量之间的关系，还是分析两个信号或一个信号在一定时移前后之间的关系，都需要应用相关分析。 通常，两个变量之间存在对应的确定关系，则为函数关系。如球的体积与半径之间的依存性。 随着某一变量数值的确定，另一变量却可能取许多不同值，但取值有一定概率统计规律。如农产品的含水率与干燥时间的关系，每一确定干燥时间，含水率不是一恒定值(在一定范围内)，但有规律，随干燥时间增加，含水率下降。 b）有相关性。 a）不相关。

13、 如果所研究的变量x, y是与时间有关的函数，即x(t)与y(t)：x(t)x(t)y(t)y(t)22 2() xxE x22 2() yyE y yxyxxyyxE 222yxyxyExEyxE 1 xy 则对于变量x和y之间的相关程度常用相关系数来表示：x y ，x y 正相关。x y ，x y 负相关。 0， x 、y之间无关。但可能存在非线性相关或函数关系。 xy 当数据点分布愈来愈近一条直线时， 的绝对值接近1，x与y的线性相关程度愈好。 的正负号表示一变量随另一变量的增加而增加或减少。xyxyxy1xyxy1xyxy10 xyxy0 xy 信号的自相关是描述信号在一个时刻的取值与

14、另一个时刻取值的依赖关系；用下式表示：01( )( ) ()limTxTRx t x tdtT （519） 式中，T为样本记录长度（即观测时间） x(t)x(t)时时延延器器 乘乘法法器器 x(t +)x(t)x(t +)积积分分 器器 Rxx()2222( )xxxxxR1 1）即： 时， 2(0)xxR为最大值。 02 2）自相关函数在 时为最大值，并为该信号的均方值 。2x。0自相关函数的性质：自相关函数的性质： 3 3）当 足够大或 时， 20 xxxR 4 4）自相关函数为偶函数 。()( )xxRR5 5）周期函数的自相关函数仍为同频率的周期 函数，其幅值与原周期信号幅值有关，而

15、丢失了原信号的相位信息。 图517 典型信号的自相关函数 案例：机械加工表面粗糙度自相关分析案例：机械加工表面粗糙度自相关分析 提取出回转误差等周期性的故障源。提取出回转误差等周期性的故障源。被测工件被测工件相关分析相关分析案例：案例：自相关分析测量转速自相关分析测量转速理想信号理想信号干扰信号干扰信号实测信号实测信号自相关系数自相关系数提取周期性转速成分。提取周期性转速成分。自相关分析的主要应用：自相关分析的主要应用：用来检测混淆在干扰信号中用来检测混淆在干扰信号中的确定性周期信号成分。的确定性周期信号成分。 两个各态历经的随机信号x（t）和y（t）的互相关函数 ： ( )xyR01( )(

16、 ) ()limTxyTRx t y tdtT (525)x(t)y(t)时时延延器器 乘乘法法器器 y(t +)x(t)y(t +)积积分分 器器 Rxy()互相关函数的性质互相关函数的性质峰值点峰值点互相关函数不是偶函数，其峰值偏离原点的位置反映了两信号在错开多大时，相关程度最高。 Rxy（）Ryx（）3）频率相同的两个周期信号的互相关函数的仍是 周期函数，其周期与原信号相同。保留了两个 信号的频率，幅值和相位的信息（同频相关， 不同频不相关）。 4）若 ( )( )x ty t、完全独立无关，且均值为0，则 ， ( )0 xyR均值不为0，( )xtxyR ；2)( )( )xyyxRR

17、与不同， ()( )xyyxRR关函数的对称性，见上图 ，即互相1 1）检测混淆在噪声中的信号）检测混淆在噪声中的信号 激振信号x（t）,所测定的信号为y（t）+ n（t），n（t）为噪声信号。 那么做 互相关分析，y（t）与 n（t）频率无关，那么仅x（t）和 y（t）的互相关函数 存在，并反映了激振的响应幅值和相位差。消除了噪声。 ( )xyR( )( )( )x ty tn tddv 12Sv 据互相关图，发动机与座椅之间相关性差，而车轮与座椅之间互相关明显，那么振动主要由车轮引起。 注意：注意：对能量有限信号进行相关分析时，应按下面定义来计算： dttytxRdttxtxRxyx 矩形

18、脉冲信号、衰减指数函数等 对于随机信号，可用有限时间的样本记录所求得的相关函数值来作为随机信号相关函数的估计，即： TxyTxdttytxTRdttxtxTR0011 对于周期信号，可用一个周期内的观察值的平均值代表整个过程的平均值。 对于有限个序列点N的数字信号的相关函数估计，可写成： NmrrnynxNrRrnxnxNrRNnxyNnx , 2 , 1 , 0111010 x（t）是零均值的随机信号（若x（t）是非零均值的，可以通过适当处理使其均值为0），并假定x（t）中没有周期分量，则： dRRxx0lim 将自相关函数的傅立叶变换 （534） deRfSfjxx2（535） dfefS

19、Rfjxx 2 定义为 的自功率谱密度函数，简称自功率谱自功率谱或自谱自谱。据傅立叶逆变换，有： ( )x t fSx当=0 时，上式为： (前面已述，均方值为信号的平均功率)。 dffSdttxTlimRxTTx0210 这表示自功率谱 与频率轴所包围的面积就是信号的平均功率。因此， 给出了信号中各频率分量的功率沿频率轴的分布，故称 为自功率谱密度函数。 fSx fSx fSx 为实偶函数。 称为信号x（t）的单边功率谱。如图5-22所示。 fSx fSfGxx2 dffXdttx22 称为能谱，它是沿频率轴的能量分布密度，在整个时间轴上信号平均功率为 2fX dffXTdttxTPTTTa

20、v2021lim1lim巴塞伐尔定理巴塞伐尔定理 dffSdttxTlimRxTTx0210根据式得自功率谱密度函数和幅值谱的关系为： 所以功率谱可通过直接对时域信号作傅里叶变换来计算。 21fXTlimfSTx (5-38)功率谱的估计功率谱的估计 用有限长度T的样本记录来计算样本功率谱，并以此作为信号功率谱的初步估计值： 2221fXTfGfXTfSxx 对于数字信号： 2221kXNkGkXNkSxx 1 1）分析信号的成分和结构）分析信号的成分和结构 与幅值谱 均可反映频率结构，但 反映的是信号幅值的平方，因此其频率结构更为明显（见图523）。 fSx fSx fX2 2）故障的判断和

21、分析）故障的判断和分析 据功率谱的变化来判断发生原因，以便排除。 3 3）通过输入、输出自谱的分析，可得系统的）通过输入、输出自谱的分析，可得系统的 幅频特性幅频特性 fSfHfSxy2 但在计算中丢失了相位信息，因而不能得出系统的相频特性。 4 4）用来检测信号的周期成分）用来检测信号的周期成分 周期成分在实测的自功率谱密度图形中以陡峭有限峰值的形态出现。 如果互相关函数满足傅立叶变换条件： dRxy deRfSfjxyxy2 为x（t）和y（t）的互谱密度函数，简称互谱互谱。则定义 互谱和互相关函数构成一对傅立叶变换对，所以二者包含有相同的信息，都可用来描述信号之间的相关性，不同点是互相关函数在时间域上，而互谱在频率域上。 互谱的傅立叶逆变换为： （546） dfefSRfjxyxy 2 互谱估计的计算式如下： ii*yxii*xyfXfYTfSfYfXT