现代信号处理第八章基于emd的时频分析方法及其应用.pptVIP

构造式(8.3.1)所示的仿真信号，来模拟单自由度结构前三阶模态的响应信号：(8.3.1)其中表示第个脉冲响应信号：(8.3.2)它们得频率分别为Hz，Hz，Hz；阻尼比分别为，，。冲击发生的时刻为0.05s，N表示幅值为1的白噪声。直接采用Laplace小波相关滤波法的不足直接采用Laplace小波相关滤波法的不足01最终的仿真信号及其组成如图所示。02图仿真信号及其组成直接采用Laplace小波相关滤波法的不足对该仿真信号直接进行Laplace相关滤波提取第二阶模态参数为例，结果如下图所示：图8.3.2仿真信号直接提取第二阶模态结果由图可见，相关系数始终处于较低的水平，频率曲线有较大的波动，这说明无法找到与原始信号相似的Laplace小波，难以直接提取准确的模态参数。基于EMD的Laplace小波模态参数识别方法由于直接利用Laplace小波滤波法识别参数遇到困难，故首先对上述仿真信号进行EMD分解，由于信号中的有用部分（冲击响应波形）处于信号中部，两端各有一段无用的白噪声，故不用考虑EMD的边界效应。1/s2图8.3.3仿真信号及其EMD分解结果3基于EMD的Laplace小波模态参数识别由于EMD分解总是先分解出高频分量，所以第一个IMF()就是第三阶模态对应的响应信号，对应第二阶模态，对应第一阶模态的响应信号。对第二个IMF信号进行Laplace相关滤波提取第二阶模态参数，结果如图所示：图8.3.4第二个分量提取第二阶模态结果仿真信号提取结果表1给出了信号前三阶模态参数的理论值、利用频谱细化方法和本文方法识别的结果。表8.3.1仿真信号模态参数识别结果理论值60.00.020400.00.0101000.00.005频谱细化相对误差(%)59.6820.53399.870.033999.9530.005直接相关滤波相对误差(%)59.6500.5830.02525.0399.300.1750.01330.0999.8500.0150.00740.0本方法相对误差(%)59.9800.0330.0200.0399.900.0250.0100.0999.9900.0010.0050.0无阻尼固有频率60.004399.9401000.015得到结构的阻尼固有频率和阻尼比之后，可由下式计算结构的无阻尼固有频率：应用实例为了验证本文所述方法的正确性，搭建了如下图所示的悬臂梁模态识别实验台。采样频率设为3000Hz，采样长度为3000。图悬臂梁模态识别实验台应用实例左图是实测响应信号及其EMD分解结果。右图表示了它们对应的频谱，可见响应信号中包含了悬臂梁的前三阶固有频率，通过EMD分解，响应信号完全分解成了与前三阶模态一一对应的三个分量。图8.3.6实测信号及其EMD分解结果图8.3.7图中各信号对应的频谱应用实例对分解所得第二个分量进行Laplace相关滤波，提取其第二阶模态参数的结果如下图所示。1图8.3.8第二个分量提取第二阶模态2040301EMD的基本理论和算法基于EMD的Laplace小波结构模态参数识别方法研究EMD实用化技术研究 EMD方法在机械设备故障诊断中的应用 02第八章基于EMD的时频分析方法及其应用EMD方法的基本原理EMD方法的完备性和正交性基于EMD的Hilbert变换(HHT)的基本原理和算法EMD方法的基本概念020304018.1EMD的基本理论和算法时间序列的Hilbert变换为:构造解析函数其中幅值函数相位函数相位函数的导数即为瞬时频率(8.1.2)(8.1.1)(8.1.3)(8.1.5)(8.1.4)(8.1.6)瞬时频率的概念瞬时频率的概念(8.1.7)(8.1.8)(8.1.10)(8.1.9)(8.1.11)然而按上述定义求解的瞬时频率在某些情况下是有问题的,考虑如下信号这个信号是解析的，按式(8.1.3)和(8.1.4)可以求解其相位和幅值，得到假设信号幅值是恒定的,频率是正的，信号的频谱瞬时频率的概念”0405当两个正弦频率取，两个频率时，幅值的取值不同，其瞬时频率亦有很大的不同。如图8.