现代信号处理信号的域分析

现代信号处理信号的域分析第1页，共57页。引言以时间为自变量描述物理量的变化是信号最基本、最直观的表达形式。在时域内对信号进行滤波、放大、统计特征计算、相关性分析等处理，统称为信号的时域分析。通过时域分析方法，可以有效提高信噪比，求取信号波形在不同时刻的相似性和关联性，获得反映机械设备运行状态的特征参数，为机械系统动态分析和故障诊断提供有效的信息。第2页，共57页。第二章信号的时域分析2.1信号的预处理2.2信号的采样2.3时域统计分析2.4相关分析及应用第3页，共57页。2.1信号的预处理传感器获取的信号往往比较微弱，并伴随着各种噪声。不同类型的传感器，其输出信号的形式也不尽相同。为了抑制信号中的噪声，提高检测信号的信噪比，便于信息提取，须对传感器检测到的信号进行预处理。所谓信号预处理，是指在对信号进行变换、提取、识别或评估之前，对检测信号进行的转换、滤波、放大等处理。第4页，共57页。2.1信号的预处理常用的信号预处理方法信号类型转换应变测力传感器、热电阻传感器输出的信号均为电阻信号，为了便于后续处理常用电桥将电阻信号转变为电压信号信号放大常用的信号放大器包括：测量放大器、隔离放大器、可编程增益放大器等信号滤波（本节重点介绍）去除均值在计算信号的标准差等统计量时，需要去除信号均值去除趋势项常用的趋势项消除方法有滤波法、多项式拟合法第5页，共57页。2.1.1信号的滤波处理

信号滤波处理是消除或减弱干扰噪声，保留有用信号的过程。把实现滤波功能的系统称之为滤波器。滤波器可分为两大类，即经典滤波器和现代滤波器。第6页，共57页。1.经典滤波器定义：当噪声和有用信号处于不同的频带时，噪声通过滤波器将被衰减或消除，而有用信号得以保留分类根据幅频特性的不同，滤波器分为低通滤波器、高通滤波器、带通滤波器、带阻滤波器等类型。根据处理信号类型的不同，滤波器可分为模拟滤波器和数字滤波器。对于数字滤波器来说，根据滤波器的单位脉冲响应序列长度的无限和有限，数字滤波器可进一步分为无限冲击响应滤波器(IIR)和有限冲击响应滤波器(FIR)两类第7页，共57页。1)经典滤波器原理经典滤波概念和方法建立在频域分析基础上滤波后的信号为滤波器的传递函数或滤波器的频率响应函数在噪声频带和有用成分频带分离的情况下，通过设计如下的滤波器函数滤波器的单位脉冲响应函数(2.1.1)第8页，共57页。2)理想模拟滤波器理想模拟滤波器是一个理想化的模型，对其讨论有助于进一步了解和改进实际滤波器的性能理想模拟滤波器的幅频特性曲线图2.1.1第9页，共57页。理想模拟滤波器由于理想低通滤波器具有矩形幅频特性和线性相位特性。同时，理想高通、带通和带阻滤波器均可以由理想低通滤波器串联得到。因此，以后均以理想低通滤波器为例来说明理想低通滤波器的矩形幅频、相频特性(2.1.4)第10页，共57页。理想模拟滤波器理想低通滤波器的单位脉冲响应函数(2.1.5)第11页，共57页。3)实际滤波器及其基本参数实际的滤波器为了物理上可实现，通常在通带和阻带之间设置过渡带通带边缘频率阻带边缘频率截止频率(滤波器的半功率点)第12页，共57页。实际滤波器及其基本参数实际滤波器的幅频特性幅值在通带和阻带内一般不严格为1和0。它们分别允许的波动量分别为波动的大小分别用通带和阻带内的衰减实际滤波器的参数还有：波纹幅度、带宽、品质因数和倍频程选择性等第13页，共57页。4)数字滤波器的设计数字滤波器有无限冲击响应IIR型滤波器和有限冲击响应FIR型滤波器之分IIR型数字滤波器的传递函数是FIR型数字滤波器的传递函数是滤波器的单位脉冲响应函数第14页，共57页。数字滤波器的设计方法简述FIR型数字滤波器的设计方法主要是建立在对理想滤波器频率特性作某种近似的基础上的。这些近似方法有窗函数法、频率抽样法等；IIR型数字滤波器的设计属于间接设计法。IIR型数字滤波器目前最通用的设计方法是利用已经很成熟的模拟滤波器的设计方法来进行设计。而模拟滤波器的设计方法又有巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)和椭圆滤波器等不同的设计方法。第15页，共57页。数字滤波器的设计方法简述不论是FIR型数字滤波器还是IIR型数字滤波器的设计都包括三个步骤Matlab数字滤波器设计演示(FDATool.fda,SPTool.fda)第16页，共57页。2．现代滤波器当噪声频带和有用信号频带相互重叠时，经典滤波器就无法实现滤波功能现代滤波器也称统计滤波器，从统计的概念出发对信号在时域进行估计，在统计指标最优的意义下，用估计值去逼近有用信号，相应的噪声也在统计最优的意义下得以减弱或消除常用的统计滤波器有维纳滤波器和卡尔曼滤波器两类第17页，共57页。1）维纳滤波器20世纪40年代第二次世界大战期间，由于军事上的需要，Wiener提出并解决了平稳过程的最佳线性滤波问题采用线性最小均方误差估计准则，设计的最佳滤波器称为维纳(Wiener)滤波器h(t):脉冲响应函数，按最小均方误差准则确定(2.1.8)(2.1.9)第18页，共57页。维纳滤波器维纳滤波器可根据t时刻，及t以前时刻的观测值x(t)，实现以下三个方面的应用第19页，共57页。2）卡尔曼滤波器维纳滤波器由于计算量大，难以作实时处理，故不能广泛应用，同时它对非平稳信号的滤波也无能为力。60年代初由于航天事业发展的需要，卡尔曼(Kalman)和布西(Bucy)在解决非平稳、多输入输出随机序列的估计问题中引入了状态变量，在克服维纳滤波某些局限的基础上，提出了被后人称为卡尔曼滤波的新滤波方法。该方法在雷达、通信、控制、生物和勘探等领域得到了广泛的应用第20页，共57页。卡尔曼滤波器卡尔曼滤波是线性最小均方误差滤波器的另一种处理方法。卡尔曼滤波建立在已知随机信号模型的基础上，包括模型的阶次、模型的参数和激励白噪声的统计特性。其原理是用信号前一时刻的估计值与测量值的误差项的加权平均作为当前时刻的估计值。(2.1.10)信号前一时刻估计值测量值的误差项常数a和c分别表示参数模型和测量模型的参数，用Gn表示时刻的加权系数。第21页，共57页。维纳滤波器与卡尔曼滤波器对比卡尔曼滤波的特点是把信号的先验知识用信号的模型表达出来；在时域上引入状态变量法进行处理；采用递推型的线性最小均方误差算法。卡尔曼滤波和维纳滤波都是在应用随机信号和观测噪声的前二阶矩的统计特性，以线性最小均方估计解决随机信号的滤波问题。维纳滤波需要给出随机信号和噪声的有理谱形式；卡尔曼滤波则要求把随机信号规定为白噪声驱动的线性系统的输出。维纳滤波理论适应于平稳随机过程；卡尔曼滤波适用于有限初始时间的非平稳随机过程第22页，共57页。第二章信号的时域分析2.1信号的预处理2.2信号的采样2.3时域统计分析2.4相关分析及应用第23页，共57页。将连续信号转换成离散的数字序列过程就是信号的采样，它包含了离散和量化两个主要步骤。本节主要介绍采样过程中采样与混频、量化与误差、采样长度与分辨率及窗函数与泄露四方面的内容。第24页，共57页。2.2.1采样与混频设模拟信号为x(t)，间距为

t的采样脉冲函数为p(t)采样过程是x(t)和p(t)相乘，得到离散信号x(n

t)设x(t)的傅里叶变换为x(

)。采样脉冲函数p(t)的傅里叶变换为(2.2.1)(2.2.2)第25页，共57页。模拟信号、采样脉冲函数及其频谱第26页，共57页。离散信号及其频谱根据频域卷积定理可知，则式(2.2.1)所示离散信号x(n

t)的傅里叶变换为离散信号及其频谱频谱混叠，不能由频谱准确地恢复原信号(2.2.3)第27页，共57页。采样定理采样定理：为避免混叠，采样频率

s必须不小于信号中最高频率

max的两倍即或实际中采样频率的选取往往留有余地，一般选取采样频率

s为处理信号中最高频率的2.5～4倍。另外，由于测量信号中的高频部分往往是由干扰引起的噪声或我们不感兴趣的频谱，因此采样前须先对信号进行低通滤波(又称抗混滤波)。然后再根据滤波后信号的最高频率

max设定采样间隔。第28页，共57页。2.2.2量化与误差量化是对信号采样点取值进行数字化转换的过程。量化结果以一定位数的数字近似表示信号在采样点的取值。由于模/数转换器的位数有限，模/数转换器的输出只能表达一系列具有一定间隔的电平。当模拟信号在采样点上的取值落在两个相邻电平之间时，就要舍入到相近的一个电平上，我们把这一过程称之为量化。若设模/数转换器的位数为N，采用二进制编码，转换器转换的电压范围为

V，则相邻电平之间的增量为(2.2.4)第29页，共57页。2.2.3窗函数或泄漏理论上信号的长度是无限的，但任何观测信号都是在有限时间段内进行观测的。因此，信号采样过程须使用窗函数，将无限长信号截断成为有限长度的信号。从理论上看，截断过程就是在时域将无限长信号乘以有限时间宽度的窗函数第30页，共57页。窗函数或泄漏最简单的窗函数是矩形窗无限带宽泄漏第31页，共57页。窗函数或泄漏泄漏与截断长度、所使用的窗函数等有关。不进行信号截断就没有泄漏误差另外，使用不同的窗函数泄漏大小也不同。泄漏取决于窗函数频谱的旁瓣。如果窗函数的旁瓣小，相应的泄漏也小第32页，共57页。其它窗函数三角窗第33页，共57页。其它窗函数汉宁窗第34页，共57页。2.2.4采样长度与分辩率数字信号的分辨率包括时间分辨率和频率分辨率。数字信号的时间分辨率即采样间隔

t，它反映了数字信号在时域中取值点之间的细密程度。数字信号的频率分辨率为

=2/T，其中T=N

t为数字信号的时间跨度，N为数字信号的长度。频率分辨率表示了数字信号的频谱在频域中取值点之间的细密程度。因此，当采样频率或采样间隔确定后，增大采样点数就可增加信号的时间长度和频率分辨率。第35页，共57页。第二章信号的时域分析2.1信号的预处理2.2信号的采样2.3时域统计分析2.4相关分析及应用第36页，共57页。2.3时域统计分析信号的时域统计分析是指对信号的各种时域参数、指标的估计或计算。常用的时域参数和指标包括：1)均值；2)均方值；3)均方根值；4)方差；5)标准差；6)概率密度函数；7)概率分布函数；8)联合概率密度函数等。本节先介绍常见参数的概念，然后给出它们的应用。第37页，共57页。2.3.1时域指标参数(1)均值(2)均方值、方差第38页，共57页。时域指标参数(3)概率密度函数随机信号的取值落在区间内的概率可用下式表示概率密度函数定义为概率分布函数的定义为第39页，共57页。时域指标参数(5)有量纲参数指标有量纲参数指标包括方根幅值、平均幅值、均方幅值和峰值四种。若随机信号符合平稳、各态历经条件，且均值为零，概率密度函数为p(x)，则有量纲参数指标的定义如下方根幅值平均幅值均方幅值峰值第40页，共57页。时域指标参数上述的有量纲参数指标也可在时域定义第41页，共57页。时域指标参数(6)无量纲参数指标有量纲参数指标不但与机器的状态有关，且与机器的运动参数如转速、载荷等有关。而无量纲参数指标具有对信号幅值和频率变化均不敏感的特点。这就意味着理论上它们与机器的运动条件无关，只依赖于概率密率函数的形状。所以无量纲参数指标是一种较好的机器状态监测诊断参数。无量纲参数指标包括了波形指标、峰值指标、脉冲指标和裕度指标。定义为(2.3.10)第42页，共57页。无量纲参数指标当式(2.3.10)中的l和m取不同值时，就得到如下指标：第43页，共57页。2.3.2参数指标的应用1)利用概率密度函数和概率分布函数进行产品质量控制，研究材料的强度和控制设备的工作稳定性当图(a)代表一批零件的加工尺寸，则根据图(b)、图(c)可判断加工过程的质量高低，进而可评价或判断机床工具是否应该调整、操作工人的技术熟练程度等。第44页，共57页。参数指标的应用2)利用振幅频次分布研究设备的随机疲劳和载荷谱将信号中上升的峰值A称为振幅峰（图(a)），以振幅峰作为横座标，以振幅峰在观测时间内出现的频次作为纵座标，就得到振幅频次图（图(b)）。振幅频次图给出了动态波形峰值出现的频次分布。对振幅频次图沿横坐标进行累计，就得到累计频次图（图(c)）。根据图(b)和(c)人们便可以了解作用于材料的随机载荷谱。这对环境模拟、进行材料的常规疲劳寿命试验以及强化试验的加载方式都是十分重要的依据。第45页，共57页。参数指标的应用概率密度函数用于机器状态判断新变速箱噪声的概率密度曲线如图2.3.3(a)所示，旧变速箱噪声的概率密度曲线如图2.3.3(b)所示第46页，共57页。参数指标的应用机器状态的时域参数指标判断方法。诊断参数指标一般应满足如下要求：①易于测量和计算，所需计算机存储量小。②能敏锐地反映和预报机器的早期故障。③不受机器运行状态，如负载、转速等变化的影响。④能够指示故障的存在，以便及时排查故障。第47页，共57页。参数指标的应用图2.3.4表示了28只汽车后桥齿轮在不同运行状态下，由振动加速度信号计算得到的无量纲参数指标波形指标K的变化很小，没有足够的诊断能力；而峰值指标C和脉冲指标I可以作为齿轮运行状态的优良诊断指标。第48页，共57页。第二章信号的时域分析2.1信号的预处理2.2信号的采样2.3时域统计分析2.4相关分析及应用第49页，共