测试技术-3-2.ppt

1、1,第三章 数字信号处理技术,1、数字信号处理的主要研究内容,数字信号处理主要研究用数字序列来表示测试信号，并用数学公式和运算来对这些数字序列进行处理。内容包括数字波形分析、幅值分析、频谱分析和数字滤波等。,3.1 数字信号处理概述,2,2、测试信号数字化处理的基本步骤,3,3、数字信号处理的优势,1)用数学计算和计算机显示代替复杂的电路 和机械结构,4,2)计算机软硬件技术发展的有力推动,a)多种多样的工业用计算机。,5,b)灵活、方便的计算机虚拟仪器开发系统,6,3.2 模数(A/D)和数模(D/A),采样利用采样脉冲序列，从信号中抽取一系列离散值，使之成为采样信号x(nTs)的过程。,编

2、码将经过量化的值变为二进制数字的过程。,量化把采样信号经过舍入变为只有有限个有效数字的数，这一过程称为量化。,1、A/D转换,7,8,2) A/D转换器的技术指标,(1) 分辨率; 用输出二进制数码的位数表示。位数越多，量化误差越小，分辨力越高。常用有8位、10位、12位、16位等。,为了缩短量化一个采样值所需的时间，一般均把采样值经过舍入变为只有有限个有效数字的整数，这个整数最大为255时，用二进制表示有8bit，最大为4095时，用二进制表示有12bit，分别表示量化能达到的精度。采样后的数据总是离散的及有限的，这就引出了离散傅里叶变换的问题，即希望电子计算机算出的离散傅里叶变换与原信号的

3、连续傅里叶变换是一致的。,9,(3) 模拟信号的输入范围; 如，5V， +/-5V，10V，+/-10V等。,(2) 转换速度; 指完成一次转换所用的时间，如:1ms(1KHz)； 10us(100kHz) 时间间隔t必须大于量化一个采样值所需的时间才能保证采样能连续进行下去，f =1/t称为采样频率，各种A/D转换器件都规定了它能工作的最高采样频率。,10,2、D/A转换过程和原理,D/A转换器是把数字信号转换为电压或电流信号的装置。,11,3.3 采样定理,采样是将采样脉冲序列p(t)与信号x(t)相乘，取离散点x(nt)的值的过程。,12,X(0), X(1), X(2), , X(n)

4、,13,每周期应该有多少采样点 ？,14,15,频域解释,16,采样定理,为保证采样后信号能真实地保留原始模拟信号信息，信号采样频率必须至少为原信号中最高频率成分的2倍。这是采样的基本法则，称为采样定理。,Fs 2 Fmax,17,需注意，满足采样定理，只保证不发生频率混叠，而不能保证此时的采样信号能真实地反映原信号x(t)。工程实际中采样频率通常大于信号中最高频率成分的3到5倍。,18,频混计算：,混迭频率=Fs-信号频率,19,A/D采样前的抗混迭滤波：,20,3.4 信号的截断、能量泄漏,为便于数学处理，对截断信号做周期延拓，得到虚拟的无限长信号。,用计算机进行测试信号处理时，不可能对无

5、限长的信号进行测量和运算，而是取其有限的时间片段进行分析，这个过程称信号截断。,21,周期延拓后的信号与真实信号是不同的，下面我们就从数学的角度来看这种处理带来的误差情况。,设有余弦信号x(t), 用矩形窗函数w(t)与其相乘，得到截断信号: y(t) =x(t)w(t),将截断信号谱 XT()与原始信号谱X()相比较可知，它已不是原来的两条谱线，而是两段振荡的连续谱. 原来集中在f0处的能量被分散到两个较宽的频带中去了，这种现象称之为频谱能量泄漏。,22,周期延拓信号与真实信号是不同的：,能量泄漏误差,23,克服方法之一：信号整周期截断,24,3.5 DFT与FFT,由于计算机的应用与普及，

6、快速傅立叶变化（FFT）技术的出现，数字信号分析发展越来越快，已完全替代模拟分析的趋势。进行数字信号分析，可以采用专用的数据处理机，用硬件进行快速傅立叶变换，也可以在电子计算机的控制下进行A/D转换，即采样，再用软件或硬件进行快速傅立叶变换及计算其它函数。,25,1、离散傅立叶变换,离散傅里叶变换（Discrete Fourier Transform）一词是为适应计算机作傅里叶变换运算而引出的一个专用名词。,周期信号xT(t)的傅里叶变换：,26,对周期信号xT(t)采样，得离散序列xT(n),将积分转为集合：,展开，得连续傅立叶变换计算公式：,用计算机编程很容易计算出指定频率点值：,27,f

7、=？ /计算的频率点 Fs=? N=1024 dt=1.0/Fs pi=3.1415926 XR=0 XI=0,For n=0 To N-1 XR=XR+x(n)*cos(2*pi*f*n*dt)*dt XI=XI+x(n)*sin(2*pi*f*n*dt)*dt Next A=sqr(XR*XR+XI*XI) Q=atn(XI/XR),VBScript 样例,28,采样信号频谱是一个连续频谱，不可能计算出所有频率点值，设频率取样间隔为：,f = fs / N,频率取样点为0,f,2f,3f,.，有：,29,下图为离散傅里叶变换的图解分析。图a中，左边是被分析的信号x（t）右边是希望得到的连续

8、傅里叶变换 X（f）。图 b中，左边是采样脉冲序列0（t），右边是其傅里叶变换 0（f） 。图c中，左边是被分析的信号x（t） 与采样脉冲序列 0（t）的乘积，右边是其傅里叶变换，相应于图a中右边信号与图b中右边信号的卷积，由于被分析信号中包含有较高频率成分，而采样频率相对较低，因而频谱的高端出现混叠。图d中，左边是矩形脉冲W（t），右边是其傅里叶变换W（f）。图e中，左边是有限个采样脉冲序列，这是实际情况，相应于图c中左边信号与图d中左边信号的乘积，右边是其傅立叶变换，由于时域中信号被截断，反映在频谱高端产生了皱波。,30,图f中，右边是频域采样脉冲序列 1（f），左边是其逆傅里叶变换 1（

9、t）。图g中，由于经过离散傅里叶变换后只能得到离散的频谱即右边信号X（n），相应于图 e中右边信号与图 f 中右边信号的乘积，其包络确实近似地表示出连续傅里时变换，左边是其逆傅里叶变换x（k），即还原的被分析信号，其包络也近似地表示出原被分析信号，但已将其周期化。 由以上简单分析可知，若不希望产生频率混叠现象，采样频率必须大于最高被分析频率两倍以上，或必须将大于1/2采样频率以上的频率成分滤掉，即进行抗混滤波，这称为香农（shannon）采样定理；若希望减小皱波，则时域中不能突然截断信号，应根据信号的特点，采用不同的窗函数，减小皱波效应。,31,32,2、快速傅立叶变换 (FFT),DFT为离

10、散信号的分析处理从理论上提供了变换工具，而快速付里叶变换FFT（Fast Fourier Transform）则是作为DFT的一种快速算法而发展起来的。我们知道，有限长度序列数字信号处理中占有重要的地位，从时间域到频率域处理这些序列的手段就是DFT。然而，按照DFT的定义式进行计算，尽管是采用计算机，也要花费很长的时间，因此，DFT的普及应用受到一定的限制。1965年美国学者库勒（J.W.Cooley）和图基（J.W.Tukey）提出了离散付里叶变换的一种快速算法，使计算工作量大为减少，从而使把时域问题转换到频域的高效处理成为可能。,33,34,按定义计算DFT，需要N2次复数乘法及N（N1）

11、次复数加法，若N=1024，则实数乘和实数加各为419万次，计算量非常巨大。FFT只需Nlog2 N数量级的运算，计算量比值为1/100，并且序列越长，计算量的减少越显著，因此大大地节约了运算时间。 从三角公式也可看出有大量重复的cos、sin计算，FFT的作用就是用技巧减少cos、sin项重复计算。,35,FFT的具体实现方法有多种，关于这方面的讨论可参阅数字信号处理等有关文献资料，此处只简要介绍其中最基本、同时也是最重要的基 2时间选抽算法的原理。 FFT的思路是：将序列逐次奇偶对分，直到各子序列只含有一个数，求出单项序列的DFT，（为其自身），然后合成两序列的 DFT，由此再合成四项序列

12、的DFT，-，最后由两个N/2项序列的DFT合成原序列的DFT。,36,37,38,这里相位因子WkN = e-j2nk/N有三个重要性质： （1）周期性Wk+mNN= WkN； （2）对称性Wk+N/2N=- WkN； （3）换底公式： WmNN= WkN/m。 正是利用这些性质，可以避免DFT计算式中很多不 必要的重复计算，减少计算量，加快DFT的运算速 度。,39,40,41,从上图可见，变换后的输出序列X（k）按正序排列，但在输入序列X（n）的排列次序不是原来的自然顺序，而是变成了0，4，2，6，1，5，3，7。这是由于对原序列作了三次奇偶分解后得到的，如下图所示。掌握这一规律可对N为

13、2的任意次幂的序列均能作出正确的抽取顺序。,数据整序方法 （奇偶分解整序）,0,42,3.6 栅栏效应与窗函数,1、栅栏效应,为提高效率,通常采用FFT算法计算信号频谱，设数据点数为N，采样频率为Fs。则计算得到的离散频率点为:,Xs(Fi) ， Fi = i *Fs / N , i = 0，1，2，.，N/2,如果信号中的频率分量与频率取样点不重合，则只能按四舍五入的原则，取相邻的频率取样点谱线值代替。,43,可见在进行DFT的过程中，最后需对信号的频谱作采样，经过这种采样所能显示出来的频谱仅在各采样点上，而不在此类点上的频谱一律显示不出来，即使在其它点上有重要的峰值也会被忽略，这就是栅栏效

14、应。,44,2 能量泄漏与栅栏效应的关系,频谱的离散取样造成了栅栏效应，谱峰越尖锐，产生误差的可能性就越大。,例如，余弦信号的频谱为线谱。当信号频率与频谱离散取样点不等时，栅栏效应的误差为无穷大。,45,实际应用中，由于信号截断的原因，产生了能量泄漏，即使信号频率与频谱离散取样点不相等，也能得到该频率分量的一个近似值。,从这个意义上说，能量泄漏误差不完全是有害的。如果没有信号截断产生的能量泄漏，频谱离散取样造成的栅栏效应误差将是不能接受的。,46,能量泄漏分主瓣泄漏和旁瓣泄漏，主瓣泄漏可以减小因栅栏效应带来的谱峰幅值估计误差，有其好的一面，而旁瓣泄漏则是完全有害的。,47,3 常用的窗函数,1

15、）矩形窗,48,2）三角窗,49,3）汉宁窗,50,常用窗函数,51,3.7 数字信号分析处理中存在的问题 及解决方法,1）量化误差 时间离散、幅值离散。量化误差的最大值为数字编码最后位所代表值的一半。 途径：加大量化编码的位数。十二位、十六位。 2）混叠 措施： 提高采样频率， 降低信号的最高频率fmax。抗混叠滤波。,52,3）泄漏 措施： 增加截断长度。 采用不同的窗函数。目的使主瓣突出，二是使旁瓣尽快衰减。但实际上两者往往不可兼的得。,4）栅栏效应 措施： 整周期截取; 采用不同的窗函数。,53,图示周期信号作整周期截取的DFT,54,图示 周期函数作非周期截取的DFT,55,通过加窗

16、控制能量泄漏，减小栅栏效应误差：,加矩形窗,加汉宁窗,56,说明： 对x（t）采一段样，相当于用一矩形窗去截取原信号，在矩形窗外的信号值都假设为零，在矩形窗内同原信号。然后在作离散付里叶变换时，又相当于强迫原信号成为窗长度为周期的周期信号。当原信号不是周期信号，或者虽是周期信号，但截取长度不等于整周期时，就歪曲了原信号，如 图所示的正弦波（图中（a）当非整周期截断时，等价的周期函数将不再是正弦波（图中（b），这就改变了信号的频率结构。 又如图正弦波整周期截断时，其频谱图上有一根谱线，但若是非整周期截断，就出现了原来频率以外的许多频率成分，即出现了频率泄漏现象。为了克服这种现象因而采用了各种不同

17、的窗函数。,57,58,小结：,信号截断能量泄漏,FFT栅栏效应,从克服栅栏效应误差角度看，能量泄漏是有利的。,59,3.8 时域统计量的离散运算及平均,采样后的信号是离散的，积分相应地要变为求和运算，所以上面的各运算公式均要按求和来运算。这里就不在赘述了。 由于数据处理设备的限制，采样后的数据点数N不可能非常大，常用的点数一般为1024、2048、4096等，只按一个样本的结果来分析信号，尤其对于平稳随机过程会产生较大的估计误差。解决的办法是，根据数据处理设备的特点，可每次采样 N点，进行处理并保存中间结果，再采样 N点，进行处理与前次中间结果平均，就能减少估计误差。这种方法称为集合平均的方

18、法，它是算术平均的一种推广。,60,常用的平均方法有两种： l）线性平均。假定每次处理的结果为an，平均后的中间结果为An，第n-1次算术平均应为：,61,式中M为任意正整数，当nM时，可以证明前次的中间结果对当前平均值的影响随平均次数n的增加而减弱，M越大则减弱越慢。指数平均适合渐变的非平稳信号，如用于计算均值、均方值或方差等，可观察到其随时间渐变的过程。,62,另一种不同的平均方式称为时域同期平均，这种方式需要与有用信号同步的触发信号，例如需要分析齿轮的啮合频率信号，可在齿轮轴上设置光码盘，利用光码盘的信号启动采样，此时采到的有用信号相位可保持完全一致，对时间信号直接进行平均，就可消除那些

19、相位与触发信号不同步的随机信号，平均后再求功率谱密度函数，可以提高信噪比。,63,3.9 Zoom-FFT,信号分析和故障诊断中经常会遇到信号中有频率很密集的谐波成分，用普通的谱分析方法就很难加以辩别。例如在10kHz范围内作1024点FFT，频率分辩率只有25Hz，即相邻两根谱线的间距为25Hz，如果有相邻5Hz的谐波成分就不能分辨它们了，为此要求有高分辨率的谱分析方法。照相技术中用局部放大的办法，或用变焦距镜头来观察景物的细部以提高分辨率。“频率细化”的英文名词Zoom就是从照相技术中借用过来的，其原意为变焦距。,64,65,从理论上讲，提高频率分辨率只有增加信号的采样长度T，如果要提高1

20、0倍分辨率就要采集10倍长的信号来作离散付里叶变换。但这样计算时间就要为原来的100倍（对 DFT方法，计算量与 N2成正比）或33.2倍(对FFT，计算量与NLog2N成正比），而且对信号分析仪来讲还要受到计算FFT点数的限制，一般固定为512、1024或明或2048点。,66,细化分析方法很多，这里以应用较为广泛的复调制细化方法为例说明其原理和步骤。,一般的FFT分析是一种基带的分析方法，在整个分析带宽内，频率是等分频率的，即,67,68,假设原频率分析范围为fmax采样点数为N，采样频率fS2.56 fmax ，则频率分辨率f2.56 fmax N。现若要在f1- f2区间将频率分辨率提

21、高M倍（Zoom Factor）， f2 - f1 fmaxM，则采样长度应为MN点。其分析步骤为：,69,（1）以采样频率fS采集MN点，得离散数据Xt ，t0，1，2，MN1。 （2）将信号序列Xt乘以单位旋转矢量e-j2f1t得到新的序列yt ，按傅氏变换的调制性质，相当对原信号进行频移，即将f1移到频率坐标的原点。 （3）用截止频率为f2- f1的数字低通滤波器对 yt进行数字滤波，得到Zt ，t0，1，2，MN-1。 （4）对Zt每隔M个抽选一个，组成Wt ，即Wt= ZMt，t0，1，2，N-1。 （5）对Wt进行N点FFT就得到了细化的频谱。即在f2- f1范围内的分辨率为f =

22、 f M。,70,复调制法频率细化过程可用下图所示的框图表示，如果要改变分析范围f2- f1，则要重新按上述过程进行。目前信号分析仪可具有10400倍的频率细化能力。,71,Zoom- FFT 原理图,72,图中可以看出，在输出轴齿轮Z3、Z4的啮合频率处调制出以输出轴的旋转频率n*13.17HZ为间隔的一组边频带。,功率谱细化图(细化中心频率为Z3、Z4的啮合频率368HZ),73,第四章 信号的时序分析简介,FFT是谱估计中使用最为广泛的方法。但这种方法存在几个固有的缺陷，其中最突出的是频率分辨力受到采样长度的限制；其次是数据截取加窗的影响，在频率中表现为能量的 “泄漏”。虽然，选用适当的

23、窗函数，可以减小泄漏，然而又将导致谱分辨力和幅值精度的下降，特别是在短数据记录的情况下更为突出，这是在实际情况下经常遇到的问题。而机械冲击响应信号、机械故障源信号等只有很短的数据可用于分析；另一方面，当信号具有缓变的时变谱时，也只有在采样序列较短时，才可视其谱为时不变的。在这些情况下，基于FFT的传统谱估计方法就显得不适用了。,74,为了改善谱分析的性能，扩大信号处理应用的范围而发展了一种适于短数据序列的分析处理方法，即时序分析方法。与 FFT谱分析相对应，时序谱分析方法称为现代谱分析方法。 1、时间序列与时间序列分析 所谓时间序列，是指按时间先后顺序排列的一组数据，在“时序分析”这一学科的研

24、究范围内，时间序列则是广义地指一切有序的随机数据，包括时间上的先后有序和空间上的前后有序。,75,时间序列分析简称时序分析，它把依某一规律变化的信号（数据）看成是依时间变化而变化的先后有序的数据，在一定的假设前提下，依据某一准则建立起数学模型，以此对原时间序列或对产生这一时间序列的系统进行分析辨识。时间序列分析方法从1927年产生至今，已在诸多工程技术领域获得了极为广泛的应用，已成为动态数据处理的一种极为重要的数学工具。将时间序列分析用于机械设备的故障诊断则是近十多年的事情，通过对有关数学模型的时域特性和频域特性分析，可识别机械设备所处的工况，但其更主要的应用则是对机械设备的剩余寿命或其未来发

25、展趋势的预测。,76,时间序列分析是数理统计学科的一个重要分支，是分析随机过程的一个重要数学工具，当然，它也不是包罗万象的，而是有其应用范畴，能够采用时序方法进行分析处理的动态数据应该满足各态历经性假设。,77,2、 时间序列分析的数学模型2.1模型结构 时间序列分析中两类最基本的数学模型是ARMA模型和AR模型。（随机差分方程）,78,公式的含意是在时刻t的输出xt是此系统前n-1 个输出xt-1，xt-2，-，xt-n和由t到t-m时刻中的m个互相独立的白噪声输入的线性和。因此上式可改写为以下形式： xt=1xt-1+-+ nxt-n+ at+ 1at-1 +-+ mat-m,79,综上，

26、可将at 记为at NID（0，a2），意为均值为零，方差等于a2,80,2.2模型定阶 建模的关键是把模型的阶数和参数全部估计出来，其过程是根据xt （t = l，2，N）和一定的准则，选择 n和 m，估计参数i和 j。很明显，如果阶数不同，参数的个数和数值也就不同。 对于合适的阶数和模型参数，模型残差序列应为白噪声。当模型阶数低于实际阶数时，由于数据的动态结构尚未充分反映到模型中，一般残差2a较大。随着阶数的上升，逐步接近实际模型，2a值将下降。,81,但是，若阶数过高，由于参数估计的误差增大，又会使2a值上升。一般用定阶准则来确定其阶数。 2.3模型的特性 自功率谱密度函数是时序模型的重要特性之一。它与信号的 FF