《合工大《测试技术》课件第五章》由会员分享，可在线阅读，更多相关《合工大《测试技术》课件第五章（95页珍藏版）》请在金锄头文库上搜索。
1、第五章 信号处理初步 测试的基本任务是获取有用的信息。测试信号中既含有有用信息，也含有大量干扰噪声。 信号处理的任务对信号施加适当的加工变换，滤除干扰噪声，提取有用信息。 信号分析研究信号的构成和特征值； 信号处理信号经过必要的加工变换，以期获得有用信息的过程。 信号分析对信号本身的结构没有影响，而信号处理则有可能改变信号本身的结构。 模拟信号处理系统、数字信号处理系统来实现模拟信号处理，系统由实现模拟运算功能的电路组成。 数字信号处理系统由微型计算机和相关软件组成。信号处理内容很丰富，但本章只能介绍其中的二、三个问题。,第一节 数字信号处理的基本步骤 第二节 信号数字化出现的问题 第三节 相
2、关分析及其应用 第四节 功率谱分析及其应用 第五节 现代信号分析方法简介,第一节 数字信号处理的基本步骤,数字信号处理器或计算机对离散的时间序列进行运算处理。计算机只能处理有限长度的数据，所以首先要把长时间的序列截断，对截取的数字序列有时还要人为地进行加权(乘以窗函数)以成为新的有限长的序列。对数据中的奇异点(由于强干扰或信号丢失引起的数据突变)应予以剔除。对温漂、时漂等系统性干扰所引起的趋势项(周期大于记录长度的频率成分)也应予以分离。如有必要，还可以设计专门的程序来进行数字滤波，然后把数据按给定的程序进行运算，完成各种分析。运算结果可以直接显示或打印，若后接D/A，还可得到模拟信号。如有需
3、要可将数字信号处理结果送人后接计算机或通过专门程序再做后续处理。,第二节 信号数字化出现的问题,一、概述,从以上过程看到，原来希望获得模拟信号x(t)的频域函数X(f)，由于输入计算机的数据是序列长为N的离散采样后信号x(t)s(t)w(t)，所以计算机输出的是X(f)p。X(f)p不是X(f)，而是用X(f)p来近似代替X(f)。处理过程中的每一个步骤：采样、截断、DFT计算都会引起失真或误差，必须充分注意。好在工程上不仅关心有无误差，而更重要的是了解误差的具体数值，以及是否能以经济、有效的手段提取足够精确的信息。只要概念清楚，处理得当，就可以利用计算机有效地处理测试信号，完成在模拟信号处理
4、技术中难以完成的工作。,二、时域采样、混叠和采样定理,采样把连续时间信号变成离散时间序列的过程。这一过程相当于在连续时间信号上“摘取”许多离散时刻上的信号瞬时值。在数学处理上，可看作以等时距的单位脉冲序列(称其为采样信号)去乘连续时间信号，各采样点上的瞬时值就变成脉冲序列的强度。以后这些强度值将被量化而成为相应的数值。,TS采样间隔； N序列长度，N=T/TS； fs采样频率， fs =1/TS。,若采样间隔太小(采样频率高)，则对定长的时间记录来说其数字序列就很长，计算工作量迅速增大；如果数字序列长度一定，则只能处理很短的时间历程，可能产生较大的误差。若采样间隔过大(采样频率低)，则可能丢掉
5、有用的信息。,傅立叶变换的卷积定理,注意到原频谱X(f)是f的偶函数，并以f=0为对称轴；现在新频谱X(f)*S(f)又是以fs为周期的周期函数。因此，如有混叠现象出现，从图中可见，混叠必定出现在f=fs/2左右两侧的频率处。有时将fs/2称为折叠频率。可以证明，任何一个大于折叠频率的高频成分f1都将和一个低于折叠频率的低频成分f2相混淆，将高频f1误认为低频f2。相当于以折叠频率f2/2为轴，将f1成分折叠到低频成分f2上，它们之间的关系为：(f1+f2)2=fs/2这也就是称fs/2为折叠频率的由来。,不产生混叠的条件： a)模拟信号x(t)为带限信号 b) 奈魁斯特采样定理 通常fs=（
6、34）fc,二、量化和量化误差,量化用有限个允许值近似地代替精确值。 量化方法：截尾、舍入截尾将二进制数的多余位舍掉。舍入是将二进制数的多余位舍去或舍去后且在最低有效位上加l，这与十进制中的四舍五入法相似。 信号x(t)可能出现的最大值为A，量化单位为 当信号x(t)落在某一小间隔内，经过舍入方法而变为有限值时，将会产生量化误差e(n)量化误差的最大值为2，可以认为量化误差在(-2，2)区间各点出现的概率是相等的，其概率密度为1，均值为零。求得其标准差: s=0.29 显然，量化单位愈大，则量化误差愈大。 对信号采集时，量化增量的大小与AD转换器位数有关。 如：8位的AD转换器最大为AD转换器
7、允许的工作电压幅值的1256。,三、截断、泄漏和窗函数,信号数字化处理时，需要截断原始信号。 从原理上讲，截断就是将无限长的原始信号乘以时域有限宽的窗函数。,根据傅里叶变换关系: 截断后的频谱为余弦信号的频谱与窗函数频谱的卷积； 产生泄漏 泄漏由原来的两条谱线，变为一个两段连续谱。这表明原来信号和由其中截取的信号两者的频谱不同了。原来集中在0处的能量被分散到两个较宽的频带中去了。只要信号一经截断，就不可避免地引起混叠。减少混叠的方法：（1）增大截断长度T；（2）采用其它的窗函数,窗函数的选择：应考虑被分析信号的性质与处理要求如要求精确读出主瓣频率，而不考虑幅值精度可选用主瓣宽度比较窄而便于分辨
8、的矩形窗，例如测量物体的自振频率等；如分析窄带信号，且有较强的干扰噪声应选用旁瓣幅度小的窗函数，如汉宁窗、三角窗等；如随时间按指数衰减的函数可采用指数窗来提高信噪比,四、频域采样、时域周期延拓和栅栏效应,经过时域采样和截断后，信号的频谱在频域内还是连续的。 如果要使之数字化频率离散化，实行频域采样频域采样与时域采样相似，在频域中用脉冲序列D(f)乘信号的频谱函数，在时域里，其结果则是将信号平移至各脉冲坐标位置重新构图，从而相当于在时域中将窗内的信号波形在窗外进行周期延拓。对一函数实行采样，即是“摘取”采样点上对应的函数值。其效果有如透过栅栏的缝隙观看外景一样，只有落在缝隙前的少数景象被看到，其
9、余景象都被栅栏挡住，视为零。这种现象被称为栅栏效应。不管是时域采样还是频域采样，都有相应的栅栏效应。只不过时域采样如满足采样定理要求，栅栏效应不会有什么影响。而频域采样的栅栏效应则影响较大，“挡住”或丢失的频率成分有可能是重要的或具有特征的成分，以致于整个处理失去意义。,五、频率分辨力、整周期截断,频率采样间隙f也是频率分辨力的指标此间隔越小，频率分辨力越高，被“挡住”的频率成分越少在利用DFT(离散傅利叶变换)将有限时间序列变换成相应的频谱序列的情况下，f和分析的时间信号长度T的关系是:f=fs/N=1/T (7-14)这种关系是DFT算法固有的特征。这种关系往往加剧频率分辨力和计算工作量的
10、矛盾。根据采样定理，若所感兴趣的最高频率为fh，最低采样频率fs应大于2fh。提高频率分辨力就必须增加数据点数N，从而急剧地增加计算工作量。,两条途径:(1)在DFT的基础上，采用“频率细化技术(ZOOM)”，其基本思路是在处理过程中只提高感兴趣的局部频段中的频率分辨力，以此来减少计算工作量。(2)改用其他把时域序列变换成频谱序列的方法 在分析简谐信号的场合下，需要了解某特定频率f0的谱值，希望DFT谱线落在f0上。单纯减小f，并不一定会使谱线落在频率f0上。从DFT的原理来看，谱线落在f0处的条件是：f0f=整数考虑到f是分析时长T的倒数，简谐信号的周期T0是其频率f0的倒数，因此只有截取的
11、信号长度T正好等于信号周期的整数倍时，才可能使分析谱线落在简谐信号的频率上，才能获得准确的频谱。显然这个结论适用于所有周期信号。,因此，对周期信号实行整周期截断是获取准确频谱的先决条件。从概念来说，DFT把时窗内信号向外周期延拓。若事先按整周期截断信号，则延拓后的信号将和原信号完全吻合，接合处无任何畸变。反之，延拓后将在t=kT交接处出现间断点，波形和频谱都发生畸变。其中k为某个整数。,第三节 相关分析及其应用,在测试技术领域中，无论分析两个随机变量之间的关系，还是分析两个信号或一个信号在一定时移前后之间的关系，都需要应用相关分析。主要内容： 相关和相关系数 信号的自相关函数 信号的互相关函数
12、 相关函数的估计,1相关和相关系数,相关当两个随机变量之间具有某种内在关系时，随着某一个变量数值的确定，另一变量却可能取许多不同值，但取值有一定的概率统计规律，这时称两个随机变量存在着相关关系。,相关系数衡量两个随机变量之间相关程度大小的量被称为相关系数。,分析表明，xy1当数据点分布愈接近于一条直线时，xy的绝对值愈接近1，x和y的线性相关程度愈好，将这样的数据回归成直线才愈有意义。 xy的正负号则是表示一变量随另一变量的增加而增加或减小。当xy接近于零，则可认为x、y两变量之间完全无关，但仍可能存在着某种非线性的相关关系甚至函数关系。,二、信号的自相关函数,对各态历经随机信号及功率信号可定
13、义自相关函数Rx()为,自相关函数具有的性质,1）由式（514）有,5)周期函数的自相关函数仍为同频率的周期函数，其幅值与原周期信号的幅值有关，而丢失了原信号的相位信息。,例5.1 求正弦函数 的自相关函数，初始相角为一随机变量。 解: 此正弦函数是一个零均值的各态历经随机过程，其各种平均值可以用一个周期内的平均值表示。该正弦函数的自相关函数为,可见正弦函数的自相关函数是一个余弦函数，在=0时具有最大值，但它不随的增加而衰减至零。它保留了原正弦信号的幅值和频率信息，而丢失了初始相位信息。,图是某一机械加工表面粗糙度的波形，经自相关分析后所得到的自相关图呈现出周期性。这表明造成表面粗糙度的原因中
14、包含有某种周期因素。从自相关图能确定该周期因素的频率，从而可以进一步分析其原因。,三、信号的互相关函数,定义：,互相关函数具有的性质,1）当=0时，Rxy()不一定出现最大值，而在大于零的某处互相关函数出现最大值。 2）平稳随机过程的Rxy()是非奇非偶函数。,3）当时，有两种情况。 3-1)当两随机信号不含有相同频率的周期成分时，则两随机信号不相关，这时，Rxy()uxuy。3-2)当两随机信号含有相同频率的周期成分时，则两随机信号具有某种程度的相关性，表现为Rxy()具有周期分量，且该分量的频率与当两随机信号所含的周期分量的频率相同。,4)两个均值为零具有相同频率的周期信号，其互相关函数中
15、保留了原两信号的频率、幅值及相位信息。,5)如果随机信号x(t)和y(t)均为几个分量叠加而成则其互相关函数为组成它们各个分量两两之间互相关之和。推论：如果随机信号x(t)为几个分量叠加而成则其自相关函数为组成它们各个分量的自相关及它们两两之间互相关之和。,6）自相关函数的范围,5)例题 例5-2 设有两个周期信号x(t)、y(t)由此例可见：两个均值为零且具有相同频率的周期信号，其互相关函数中保留了这两个信号的圆频率、对应的幅值以及相位差值的信息。,例5-3 设有两个周期信号x(t)、y(t),例5-2和例5-3是自相关函数的一个重要特性： 简称为：同频相关，不同频不相关相关滤波就是利用这个
16、原理,4)互相关函数的应用互相关函数可在在噪声背景下提取有用信息的一个非常有效的手段。如果我们对一个线性系统(例如某个部件、结构或某台机床)激振，所测得的振动信号中常常含有大量的噪声干扰。根据线性系统的频率保持性，只有和激振频率相同的成分才可能是由激振而引起的响应，其他成分均是干扰。因此只要将激振信号和所测得的响应信号进行互相关(不必用时移，=0)就可以得到由激振而引起的响应信号幅值和相位差，消除了噪声干扰的影响。这种应用相关分析原理来消除信号中噪声干扰、提取有用信息的处理方法叫做相关滤波。它是利用互相关函数同频相关、不同频不相关的性质来达到滤波效果的。,V=d/d,S=1/2vm,对能量有限信号，有,四、相关函数的估计,使模拟信号不失真地沿时轴平移是一件困难的工作。因此，模拟相关处理技术只适用于几种特定信号(如正弦信号)。在数字信号处理中，信号时序的增减就表示它沿时间轴平移，是一件容易做到的事。所以实际上相关处理都是用数字技术来完成的。对于有限个序列点N的数字信号的相关函数估计，有,
