5信号处理初步-修改稿

1中北大学机电工程学院2008年3月机械工程测试技术基础熊诗波、黄长艺主编2第五章信号处理初步机械工程测试技术基础，第三版本章主要讲述内容5.1数字信号处理的基本步骤5.2信号数字化出现的问题5.3相关分析及其应用5.4功率谱分析及其应用35.0概述测试工作的目的获取研究对象（被测对象）的状态和特征信息。信号处理的目的①分离信、噪，提高信噪比；②从信号中提取有用的特征信息；③修正测试系统的某些误差，如传感器的线性误差、温度影响等。信号处理可用模拟信号处理系统和数字信号处理系统来实现。45.1数字信号处理的基本步骤5.1.1测试中的数字信号处理系统系统简图，如图5-1所示。主要环节包括：①信号预处理（信号调理）；②模数转换（A/D转换）③数字信号处理（专用的或通用计算机）；④处理结果显示图5－1数字信号处理系统的简图预处理A/D转换预处理预处理数字信号处理器或计算机结果显示x(t)y(t)A/D转换55.1.2信号的预处理1.信号预处理的目的把信号变换成适于数字处理的形式，以减轻数字处理的困难。2.信号预处理环节包括的内容①电压幅值调理以便适宜于采样，总是希望信号电压峰—峰值足够大，以便充分利用A/D转换器的精确度。6②抗混滤波进入A/D采样之前，先滤去信号中的高频噪声，以提高信噪比；③隔直如果所测信号不应有直流分量，隔离信号中的直流分量；④解调如果原信号经过调制，则应先进行解调。75.1.3模数转换模－数（A/D）转换是将模拟信号采样、量化，并转化为二进制数的过程。5.1.4数字信号处理数字信号处理，主要研究用数字序列表示信号，采用通用计算机或专用数字信号处理器，并用数字计算方法对这些离散的时间序列进行运算处理，以便把信号变换成符合某种需要的形式。8数字信号处理的主要内容①首先要把长时间序列截断，因为计算机只能处理有限长度的数据。②识别、剔除数据中的奇异点。由于强干扰或信号丢失，引起的数据突变，称为奇异点。③分离、去除趋势项。由于温漂、时漂等系统性干扰，引起信号缓慢变化的成分，称为趋势项。其周期大于记录长度的频率成分。④如有必要，还可以设计专用程序，进行数字滤波。⑤按给定的程序，完成各种分析。包括相关分析、频谱分析及信号的识别等。95.2信号数字化出现的问题由传感器获得的信号大多数为模拟信号，用计算机处理这些信号时首先要把模拟信号转化成数字信号。信号数字化过程包含着一系列步骤，每一步骤都可能引起信号失真或信息丢失。本节对信号数字化出现的问题进行讨论，提出解决这些问题的方法。105.2.0概述现以计算一个模拟信号的频谱为例来说明有关的问题，设模拟信号的傅里叶变换为（见图5-2），为了利用数字计算机处理，必须使变换成有限长的离散时间序列。为此，必须对x(t)进行采样和截断。)(tx)(fXtx(t)0│X(f)│f0图5－2原模拟信号及其幅频谱)(tx11时域采样图5－3采样函数及其幅频谱1/Ts0f…………s(t)1t0…………|s(f)|Ts1/Ts采样就是用一个等时距的周期脉冲序列s(t)，也称采样函数（图５-３）去乘x(t)。由式（１-６０）可知，s(t)的傅立叶变换S(f)也是周期脉冲序列。根据傅立叶变换的性质，采样后信号频谱应是X(f)和S(f)的卷积。12图5－4采样信号及其幅频谱0fx(t)s(t)t0│X(f)*S(f)│TsTs21Ts21Ts1卷积的结果是将它平移，使其中心落在S(f)脉冲序列的频率点上，如图（５-４）所示。若X(f)的频带大于1/2Ts，平移后的图形会发生交叠，如图（５-４）中虚线所示。采样后信号的频谱是这些平移后图形的叠加，如图中实线所示。13截断│W(f)│fT1/T-1/T0W(t)t10图5－5时窗函数及其幅频谱T由于计算机只能进行有限长序列的运算，所以必须从采样后信号的时间序列截取有限长的一段来计算，其余部分视为零而不予考虑。这等于把采样后信号乘上一个矩形窗函数，窗宽为T。所截取的时间序列数据点数（序列长度）N=T/Ts。窗函数w(t)的傅立叶变换W(f)如图（５-５）所示。14图5－6有限长离散信号及其频谱0fx(t)s(t)ω(t)t0│X(f)*S(f)*W(f)│Ts21Ts21T时域相乘对应着频域卷积，因此进入计算机的信号为x(t)s(t)w(t)，是长度为N的离散信号（见图５-６）。它的频谱函数是[X(f)*S(f)*W(f)]，是一个频域连续函数。在卷积中，W(f)的旁瓣引起新频谱的皱波。15频域采样图5－7频域采样函数及其时域函数10t…………│D(f)│Td(t)f0……Δf……-T计算机按照一定算法，比如离散傅立叶变换（DFT），将N点长的离散时间序列变换成N点的离散频率序列，并输出。注意到，x(t)s(t)w(t)的频谱是连续的频率函数，而计算后的输出则是离散的频率序列。可见对其频谱[X(f)*S(f)*W(f)]实施了频域的采样处理，使其离散化。16计算机的实际输出是X(f)p=[X(f)*S(f)*W(f)]D(f)（５-１）应当注意到频域采样形成的频域函数离散化，相应地把其时域函数周期化了，因而x(t)p是一个周期函数，如图（５-８）所示。17从以上过程看到，原来希望获得模拟信号x(t)的频域函数X(f)，由于输入计算机的数据是序列长为N的离散采样后信号x(t)s(t)w(t)，所以计算机输出的是X(f)p。X(f)p不是X(f)，而是用X(f)p来近似代替X(f)。处理过程中的每一个步骤：采样、截断、计算都会引起失真或误差，必须充分注意。下面讨论信号数字化出现的主要问题。185.2.1时域采样、混叠和采样定理把连续时间信号变成离散时间序列的过程，称为采样；在数学上，相当于用采样脉冲去乘连续时间信号；它相当于在连续时间信号上，“摘取”许多离散时刻上信号的瞬时值；各采样点上的瞬时值，就变成信号脉冲序列的强度；在随后的量化处理中，这些采样点上的瞬时值，将被量化为相应的二进制数值。19长度为T的连续时间信号，可以表示为，从点t=0开始采样，采样得到的离散时间序列为（5-2）)(txSSSSSnTtSSSTffTTNNTtxnTxNnfnxnTxnxS/1/|)()(1,,2,1,0)/()()(采样频率序列长度采样间隔上式中离散时间序列20采样间隔的选择，是一个需要重视的问题。对照比较下表信息量不够仍有误差可能产生较大误差程可能处理较长的时间历程只能处理很短的时间历如果数字序列长度一定如果数字序列长度一定可能丢掉有用信息计算工作量大数字序列短数字序列长对定长的时间纪录对定长的时间纪录采样频率低采样频率高采样间隔太大采样间隔太小21举例：图5-9说明：图5-9a，按图中采样，所得采样点无法分辨三条曲线；图5-9b，采样间隔太大，采样频率太低，将其中高频信号误认为某种相应的低频信号，出现了所谓的混叠现象图5－9混叠现象tx(t)1234ABCa)b)•••••••••••TS22下面具体解释混叠现象及其避免的办法。计算机采样，需要使用采样函数，时域采样函数及其傅里叶变换的数学描述如下；（5-3）SSSSSSSrSSnSTTfTfTTTrfTfSnTtts/1/11/11)()()(脉冲强度：脉冲间距：脉冲强度：采样频率：脉冲频率：采样间隔：频域时域)(ts)(fS23由频域卷积定理可知时域采样，相当于模拟信号与采样函数的乘积，对应的傅里叶变换为，其数学描述如式（5-4）。见图（5-4）SSSSrSSrSSTfTfTrfXTTrfTfXfSfXtstxfStsfXtx1111)()()()()()()()()(脉冲周期：采样频率：频域时域)(tx)(ts)()(tstx)()(fSfX2425不产生混叠的条件如果要求不产生频率混叠（见图5-10），首先应使被采样的模拟信号成为有限带宽的信号。为此，对不满足此要求的信号，在采样前，使其通过模拟低通滤波器滤去高频成分，使其成为带限信号，为满足下面要求创造条件。这种处理称为抗混叠滤波预处理。其次，应使采样频率大于带限信号的最高频率的2倍，即hSSfTf21hfSf（5-5）26一个满足采样定理，不产生频率混叠的例子，如图5-10所示。对于这种没有混叠的频谱，可以通过频域滤波，完整地取出原信号的频谱，也就有可能从离散序列，准确地恢复原模拟信号；图5－10不产生混叠的条件f│X(f)*S(f)│-fh0fhfhfs··········27285.2.3量化和量化误差模拟信号经采样后在时间轴上已离散，但其幅值仍为连续的模拟电压值。将模拟信号采样后的采样信号x(t)s(t)的电压幅值经过舍入或截尾的方法转变为离散的二进制数码的过程，称副值量化（简称量化）。量化过程实际上是把采样信号x(t)s(t)经过舍入或截尾的方法变为只有有限个有效数字的数的过程。29若信号x(t)可能出现的最大值A，A/D转换器的位数为B，则两个量化电平的间隔△X=A/2B-1，△X称为量化增量或量化步长。该过程是将被测信号的变化范围划分为若干区间。每个区间都用同一个整量化数字来代替。显然，只有那些正好位于区间量化电平上的离散值，才能精确地转换为量化的数字值，而那些位于区间内的离散值，只能用舍入的办法，近似到最接近的量化电平上，经过舍入或截尾的方法而变为有限值时，必然产生量化误差，最大量化误差为△X/2。305.2.4截断、泄漏和窗函数实际上，计算机只能对有限长的信号进行处理，所以，必须截断过长的信号时间历程；截断，就是对信号加窗；窗宽为的矩形窗函数，其数学描述如下：矩形窗函数，及其傅里叶变换，如图5-5所示。其它,00,1)(Tttw)(txT)(tw│W(f)│fT1/T-1/T0W(t)t10图5－5时窗函数及其幅频谱T31采用矩形窗函数截断采样信号，就是将采样信号乘以时域有限宽矩形窗函数；其时域和频域数学描述为（5-7）它们的图形如图5-6所示。所截取的时间序列数据点数，也叫时间序列长度。注意：a）频谱为连续的；b）在混叠处附近引起了皱褶。)()()()()()(fWfSfXtwtstx)(tw)()(tstx)(tw图5－6有限长离散信号及其频谱0fx(t)s(t)ω(t)t0│X(f)*S(f)*W(f)│Ts21Ts21TSTTN/N作卷积时窗函数频谱的旁瓣会引起皱波。32是一个无限带宽的函数，所以，即使是带限信号，在截断后也必然成为无限带宽的信号，这种信号的能量在频率轴分布扩展的现象，称为泄漏；为了减小截断带来的影响，常常需要合理地选择窗函数，这也是数字信号处理中的重要课题之一。窄的主瓣可以提高频率分辨能力；小的旁瓣可以减小泄漏。)(fWsinc)(tx335.2.5频域采样、时域周期延拓和栅栏效应经过时域采样和截断后，其频谱在频域成为连续的；如果要用数字描述频谱，必须在频域采样，使频谱离散化。1.连续频谱的离散化过程使用的频域采样函数为频域采样函数及其对应的时域函数如图5-7所示。1()snDffnT)(fD)(td图5－7频域采样函数及其时域函数10t…………│D(f)│Td(t)f0……Δf……-T34频域采样就是在频域中用脉冲序列D(f)乘信号的频谱函数。在时域相当于将信号与一个周期脉冲序列d(t)做卷积，其结果是将时域信号平移至各脉冲坐标位置重新构图，从而相当于在时域中将窗内的信号波形在窗外进行周期延拓。所以，频率离散化，无疑已将时域信号“改造”成周期信号。352.栅栏效应①对函数实行采样，其效果有如透过栅栏的缝隙观看外景一样，只有落在缝隙前的少数景象被看到，其余景象都被栅栏挡住，视为零；这种现象，称为栅栏效应。②不管时域采样还是频域采样，都有相应的栅栏效应。只是时域采样时，只要满足采样定理，栅栏效应不会有什