[发明专利]一种小波包辅助的自适应抗混总体平均经验模式分解方法

说明书

技术领域

本发明涉及一种非平稳微弱信号处理方法，特别是涉及一种自适应抗混总体平均经验模式分解方法。

技术背景

在雷达、故障诊断、无线电通讯、声学、生物学、地球物理等众多工程技术领域的信号处理中，存在这大量的微弱且具有非平稳的信号处理分析，目的在于从微弱信号中提取所要关注的特征信息。比如在故障诊断领域随着诊断对象日趋复杂，对故障微弱信号进行检测与特征提取已成为故障诊断领域中的关键问题之一。故障微弱信号检测与特征提取需要依赖于先进的现代信号处理方法。又如在航天测控领域的遥测信号处理，由于高频无线电信号受到传播介质(对流层、电离层、太阳等离子区)的影响，同时受到信号接收系统误差的影响，造成最终经过基带转换与A/D采样后的遥测信号受到噪声的影响，微弱的遥测信息极易淹没于噪声中，因此在航天遥测信号处理过程中提取微弱的信息也是遥测信号处理中的一个热点与难点，同样需要依赖于先进的现代信号处理方法。经验模式分解(EMD)方法是在工程中行之有效的非平稳信号处理方法之一，但模式混淆与端点效应是传统EMD方法中固有的缺陷。总体平均经验模式分解方法是近年来对传统EMD方法进行实质创新的结果，以此来达到避免在EMD方法中模式混淆的目的。总体平均经验模式分解方法的提出，较好地解决了模式混淆的问题，但其效果依赖于加入高斯白噪声的大小与总体平均的次数，而总体平均经验模式分解方法中这两个参数需要人为按照经验设定，缺乏可靠性，同时一定程度上破坏了总体平均经验模式分解方法的自适应性。因此，在总体平均经验模式分解方法基础上，建立一种能自动得到加入高斯白噪声大小与总体平均次数的自适应总体平均经验模式分解方法，对于各种不同工况下微弱信号特征提取具有重要的价值。

发明内容：

本发明要解决的技术问题是：提出一种针对不同分析信号的自适应抗混总体平均经验模式分解的信号处理方法。

本发明包括如下步骤：

步骤一：设置期望的信号分解相对误差e、分析信号的积分时间T、采样频率f_s，并计算待分析信号在积分时间内的幅值标准差σ_o。

步骤二：选择小波包的母小波函数，对待分析的信号进行小波包等频带分解，计算分解后最高频带信号的幅值标准差σ_h。

步骤三：计算ε表示信号中的高频成分的幅值标准差σ_h与原始信号幅值标准差σ_o的比值系数；根据总体平均经验模式分解中加入高斯高斯白噪声的可依据准则，建立α与ε的取值关系，自动获取加入高斯白噪声的幅值标准差比值系数α的值；由可得到总体平均经验模式分解中加入的高斯高斯白噪声幅值标准差σ_n。

步骤四：依据加入高斯白噪声的比值系数与总体平均次数选择之间的关系式：可得到总体平均经验模式分解方法中的总体平均次数N。

步骤五：对待分析信号N次加入具有幅值标准差为σ_n的高斯白噪声后，分别对其按照经验模式分解的算法进行信号分解，将分解后对应所得的基本模式分量进行总体平均计算，将总体平均后的基本模式分量作为最终输出。

步骤一中期望的信号分解相对误差e可以是一个先验的信号分解期望指标，在应用中根据期望的信号处理效果、实时性要求进行设置，在工程实践中e的取值范围为0～1％。

步骤二中对小波包母小波函数的选取可以采用Daubechies母小波函数，小波包的分解层数为2层或3层，对信号进行等频带小波包分解后，计算分解所得的最高频带信号的幅值标准差σ_h。

步骤三中所述总体平均经验模式分解中加入高斯白噪声的比值系数α的取值范围为建立总体平均经验模式分解中加入高斯高斯白噪声的可依据准则：由α的值与原始信号幅值标准差σ_o，依据可计算针对不同分析信号的σ_n取值大小。

步骤一中期望的信号分解相对误差e可以为0.5％。

在本技术方案中，提出的小波包辅助的自适应抗混总体平均经验模式分解方法(WPAEEMD)，综合利用了小波变换良好的时频多分辨率特征与经验模式分解的灵活性，能有效改善在经验模式分解避免模式混淆的问题，并能够提高总体平均经验模式分解的自适应性。

附图说明

图1是WPAEEMD信号处理流程图

图2是高斯白噪声影响信号极值点间隔分布示意图

图3是待分析的仿真信号的组成图

图4是仿真信号的小波包分解图