基于LMS-STFT-SPWVD的二值化融合分析方法研究.pdf

第 卷第 期太原科技大学学报 年 月 文章编号：（）收稿日期：基金项目：国家自然科学基金（）；忻州市重点研发计划（）作者简介：张亚斌（），男，硕士研究生，主要研究方向为振动控制算法；通信作者：李占龙副教授，：基于 的二值化融合分析方法研究张亚斌，李占龙，李虹，张正，武鹏，孙宝，刘志奇（太原科技大学，太原 ）摘要：针对二值法存在干扰信号抑制效果、中心频率能量识别效果较差且阈值参数范围难以确定的缺点，提出了 二值化融合时频分析方法。首先，采用 、算法对双分量 含噪信号测试函数进行分析，获取 、时频分布模型；结合图像处理的方法，将整个二维时频分布模型进行二值化处理，使其转化为黑白二值时频图像；将二者进行点乘运算，获取 二值化黑白融合时频图像，并与经典方法（、及 ）进行了对比研究，讨论了其干扰信号抑制效果与中心频率能量识别效果。结果表明，在阈值范围内，二值法突出了信号频率能量，具有较好的抗干扰能力。关键词：图像二值化；含噪信号；时频分析中图分类号：文献标识码：随着我国工程机械及航空航天领域的迅猛发展，重载工程车辆扮演着越来越重要的角色，其重载装备实测信号由复杂多源信号组合而成，往往为典型的复杂非平稳信号，频谱分析无法清晰识别其局部特征，面临严重的交叉频率干扰信号和信号频率分量丢等问题，给系统动态性能预测、故障识别与诊断等带来巨大挑战 。为了解决这些问题，伪维格纳威尔分布（）、平滑伪维格纳威尔分布（）相继被提出，此类算法均改善了交叉频率干扰信号抑制效果并提高了时频分辨率，但信号中心频率能量发散、交叉频率干扰信号未完全抑制及波形易受噪声信号影响已成为此类算法在诸多应用领域的主要缺点 。近些年，针对 系列算法的不足，许多优秀的时频分析方法被提出以解决上述的问题。张凤萍等 提出了 改善上述算法的缺点，但干扰信号抑制效果及中心频率信号识别精度较低。刘希康等 分别提出了 、融合算法应用于地震信号分析，该算法相比 、改善了时频聚集性、干扰信号抑制效果及中心频率信号识别精度，但信噪比及中心频率识别效果欠佳。为了改善这些缺点，文献 提出将阈值参数应用于 指数相乘法中，该算法相比 、提高了干扰信号抑制效果，但信号能量识别效果较差。同时，、等优秀的算法相机涌现，这些算法虽改善了信噪比及中心信号能量识别效果，但由于计算量过大，抗干扰能力及能量识别效果欠佳。而图像处理技术的发展为时频分析技术提供了新渠道，文献 提出了将 应用于 时频图像特性信号提取，该算法缩小了信号计算量，突出了中心频率能量识别效果，但抗干扰能力仍有待进一步提升。王九龙等 将二值化处理技术应用于螺栓松动信号提取，该方法进一步增强了信号能量识别效果，但未对阈值参数取值范围进行确定。针对现有 系列算法及二值化处理技术存在干扰信号抑制效果较差、信号中心频率能量难以识别及阈值参数范围难以确定的缺点，本文提出了一种 二值化融合时频分析方法，并对其分析路线、对比效果及参数特性等综合时频分析性能进行了研究。图像二值化时频分布模型可以等效为一张拥有 个通道的 彩色时频图像，其中 代表红色，代表绿色，代表蓝色。由于彩色图像拥有丰富的图像信息，直接处理 图像会消耗大量的算力。通常在图像处理中，对其进行阈值分割、二值化，使其转变为黑白二值图像。以降低后续步骤的难度，使其满足计算要求。因此，该方法被广泛应用于交通、机械、通信检测等领域 ，其数学模型如式（）所示。（，），（，），（，）；，（）式中，（，）为二值图像（，）处像素的灰度值，（，）为彩色图像（，）的像素值，为阈值门限。由式（）可知，二值化预处理将彩色图像中每个像素点的数值与设定的阈值门限进行比较，当像素点的灰度值小于阈值 时，将像素点的灰度值变为 ，当像素点的灰度值大于阈值 时，将像素点的灰度值变为 ，从而将图像的灰度级变为只有 和 两种像素值。二值法 二值法由 和 组成 。该算法首先采用 和 对非平稳信号进行分析，获取 和 时频分布模型；之后，将时频分布模型转换为彩色时频图像并进行图像二值化预处理，获取 、时频二值图像。随后，将二者时频图像进行点乘运算，获取 融合二值图像并进行对比分析，验证 二值法的有效性，具体步骤如下：步骤（）：将 对原始含噪信号进行分析，获取 时频分布特性模型，其模型为 （，）?（）（）（），（）（），（）（）式中，（，）为 时频域模型，（）为滤波输出时域信号，（）为短时傅里叶变换时域变换窗，（）为滤波输出时频域信号分量，（）为短时傅里叶变换时频域变换窗。步骤（）：将 对原始含噪信号进行分析，获取 时频分布特性模型。步骤（）：将 、时频分布模型转换为 彩色时频图像，并对其设置阈值门限，并对彩色图像进行阈值分割、二值化，使其转化为二值黑白图像，作为融合前数据样本。步骤（）：将二者时频图像进行点乘运算，获取 融合二值时频图像，其模型为：（，）（，）（，）（）式中，（，）项为 二值法时频域融合图像，（，）项为 算法时频域二值图像，（，）项为 算法时频域二值图像。步骤（）：判断当前时频图像融合效果是否优于 、及 三种二值图像融合效果；若达到，则输出 融合二值时频图像并再次确定阈值参数范围；反之，继续调节图像阈值，直至实现较好融合效果为止。对比分析本文以双分量 含噪信号二值图像为测试对象，采用 二值法对其进行分析，对比 、二值时频图像分析效果，验证 二值法的有效性。双分量 含噪信号时频图像双分量 含噪信号因其频率虽时间的变化而变化，保留了连续信号和脉冲信号的特性，被广泛应用于仿真、机械、雷达、声呐信号分析等领域 ，具有广泛的适用性，其数学模型及时域波形图如式（）和图 所示。（）（）（）（），（）（）（）（）（）（）（）第 卷第 期 张亚斌，等：基于 的二值化融合分析方法研究式中，（）为双分量 含噪时域信号，（）、（）为单分量 信号，（）为噪声信号，（），（）；为信噪比系数，（）为高斯白噪声信号。图 双分量 含噪信号时域波形图 为获取时频模型，本文将 滤波器及信噪比参数分别设置为 、，在信噪比为 高斯白噪声环境下，采用 、四种分析方法及参数对其进行分析，如图 所示。图 双分量 含噪信号时频分布模型图 由图 可知，系列算法时频模型与 系列算法模型数据维度不一致。因此，将二者模型融合已不可能。经前期工作发现，将上述模型进行图像二值化预处理，图像数据维度均为 为验证本文算法的有效性，将 系列时频图像阈值设置为 、系列时频图像为 ，对其进行时频二值图像分析，如图 所示。由图 可知，、虽相比 融合算法具有良好的噪声信号及模态混叠抑制效果，但交叉频率干扰信号抑制效果逊色于 而 克服了 、及 三种算法的缺点，提高了信噪比及中心频率能量识别效果，具有较好的时频分布意义。阈值参数影响为确定 二值图像及 二值图像阈值参数范围，本文通过设置三组阈值参数组合，通过对比 、两种二值太原科技大学学报 年图 二值图像效果对比图 时频图像信号能量识别、图像信息真实性及干扰信号抑制效果，进而确定阈值参数取值范围，为后续抗噪性分析提供理论依据，具体参数设置如表 所示。表 二值图像阈值参数组合设置表 参数设置组合 图像阈值 图像阈值信噪比 本文将 二值图像阈值参数设定为 ，二值图像阈值参数设定为 ，图像融合效果如图 所示。图 阈值参数组合 二值图像融合前后效果对比图 由图 可知，在阈值参数组合 条件下，二值图像信息较完整，但 二值图像仅暴露部分图像信息，融合效果较差。为改进图像融合效果，本文将 阈值参数调整为 并将 阈值参数调整为 ，继续对比图像融合前后效果，如图 所示。图 阈值参数组合 二值图像融合前后效果对比图 由图 可知，在阈值参数组合 条件下，所有融合前后二值图像信息均完整并有效抑制了噪声干扰信号，但干扰信号及中心频率能量识别效果仍有待提高。而 二值融合图像再次克服了其他图像的缺点，突出了信号中心频率能量，具有较好的低频噪声信号、交叉频率干扰信号及模态混叠抑制效果。为进一步探究阈值参数对时频图像融合效果的影响，本文再次将 阈值参数调整为 、阈值参数调整为 ，再次对比图像融合前后效果，如图所示。图 阈值参数组合 二值图像融合前后效果对比图 由图 可知，在阈值参数组合 条件下，当阈值参数再次上调后，各融合前二值时频图像抑噪效果逐渐增强，主频信号纯净度较高。融合算法仍具有较好的时频分析效果，但存在波形间断现象。因此，当 阈值参数大于 且 阈值参数大于 时，融合效果欠佳。综上，当 阈值参数取值范围为（，且 阈值参数取值范围为（，时，该融合算法时频分析效果优于 、，具有较好的时频分析效果。抗噪性分析为进一步评价 二值法的抗噪性能，本文采用双分量 不同信噪比二值图像作为测试样本并选取 图像阈值为 第 卷第 期 张亚斌，等：基于 的二值化融合分析方法研究及 图像阈值 作为测试参数，对比不同信噪比 二值融合图像噪声信号抑制效果，具体参数组合设置表及图像融合效果如表 及图 所示。表 二值图像抗噪性参数组合设置表 参数设置组合 图像阈值 图像阈值信噪比 图 双分量 含噪信号融合时频图像抗噪性效果对比图 由图 可知，当信噪比发生变化时，二值融合图像在阈值参数取值范围内均具有较好的抗噪性，信号纯净度较高。结论本文提出了 二值化分析方法，以双线性 含噪信号时频二值图像为分析对象，对比分析了 、及 干扰信号抑制效果及中心频率能量识别效果，并分析了阈值参数融合效果及其抗噪性。（）与 、及 三种算法相比，二值法克服了三种算法的缺点，提高了信噪比及中心频率能量识别效果，是一种有效的时频分析方法。（）当 阈值参数取值范围为（，且 阈值参数取值范围为（，时，融合算法具有较好的信号中心频率能量识别效果。（）当信噪比发生变化时，二值法在阈值范围内均具有较好的抗干扰能力，信号纯净度较高，但信号波形间断现象、抗噪性及融合效果仍有待进一步完善。参考文献：孔凡，李书进 非线性随机动力系统时频子域上的统计线性化方法 机械工程学报，（）：邵旋，康兴无，王旭平，等 基于多元信息融合的小波降噪质量综合评估指标 兵器装备工程学报，（）：，（）：，（）：张凤萍，刘亚奇，梁影，等 高精度引信的 滤波与 算法实现 制导与引信，（）：刘希康，丁志峰，李媛，等 与 在地震波时频分析中的应用 地震地磁观测与研究，（）：付卫红，胡展 一种混合网台跳频参数盲估计算法 北京邮电大学学报，（）：，：王豪，董广明，陈进，等 线调频小波 用于非平稳信号分析研究 北京信息科技大学学报：自然科学版，（）：纪永祯，张渝悦，朱立华，等 基于 的地震高分辨率时频分析 石油物探，（）：牛双晨，程冰洁 基于 的高分辨河道识别方法 石油地球物理勘探，（）：尚东东，张劲东，杜盈，等 基于 和改进 的复合干扰识别 数据采集与处理，（）：王九龙，盛俊杰，蒋家勇，等 时频特征提取技术的螺栓松动识别方法 噪声与振动控制，（）：李伟林，梁卓凡，方遒 一种多阈值 变换车道线检测算法 厦门理工学院学报，（）：太原科技大学学报 年 石炜，张袁祥 轴承缺陷图像二值化处理及形态学边缘检测 舰船电子工程，（）：胡天赐，罗潇远，王思博，等 遥感图像的分割与提取 中国科技信息，（）：石海杰 基于 与 的时频分析方法 硅谷，（）：（），：胡旭超，谭贤四，曲智国，等 基于 信号的风轮机叶片回波模型与分析 计算机仿真，（）：吴贞宇，夏厚培 基于 的混叠多分量 信号时域参数测量 电子科技，（）：朱晓坡，徐仕源，张卓，等 一种多目标水下定位系统信号处理算法 国外电子测量技术，（）：，（，）：，：，檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲 （上接第 页）：，（），：，第 卷第 期 张亚斌，等：基于 的二值化融合分析方法研究