振动信号的特征频谱提取方法及机械故障诊断分析方法技术

本发明专利技术公开了一种振动信号的特征频谱提取方法，计算振动信号x(t)的时频分布；在时频域中根据交错定理获得切比雪夫范数意义下的最佳逼近函数；利用最佳逼近函数，确定分数阶傅里叶变换的最佳阶次；计算最佳阶次下x(t)的分数阶傅里叶变换结果X

【技术实现步骤摘要】
振动信号的特征频谱提取方法及机械故障诊断分析方法
本专利技术涉及自动化检测技术，特别设计一种振动信号的特征频谱提取方法及机械故障诊断分析方法。
技术介绍
工业生产和日常生活领域中，旋转机械的应用极为广泛，因此旋转机械的故障是关系到生命和生产安全的重要问题，对旋转机械的故障诊断具有突出的实际价值。利用振动信号可以对不同的机械故障进行分析，根据特征频谱的不同判断机械故障类型。但是采集的振动信号中不可避免地存在噪声，而噪声的频谱存在于整个频域范围内，在故障分析中，旋转类机械结构的故障信号能量往往较低，噪声信号的存在降低了故障信号的信噪比，增加了故障诊断的难度。现有滤波方法中，滤波的中心频率一般是先验性的或通过搜索方式确定，前者主观性影响较为严重，而搜索方式的计算量相当大。
技术实现思路
本专利技术要解决的技术问题是提供一种振动信号的特征频谱提取方法，可以提高振动信号故障诊断的准确性，有助于旋转机械故障诊断分析。为解决上述技术问题，本专利技术提供的振动信号的特征频谱提取方法，其包括以下步骤：一.计算振动信号x(t)的时频分布V(f,t)；二.在时频域中，根据交错定理获得切比雪夫范数意义下的最佳逼近函数m(t)；三.利用最佳逼近函数，确定分数阶傅里叶变换的最佳阶次p；四.计算最佳阶次下的振动信号x(t)的分数阶傅里叶变换结果Xp(u)，其中u为分数阶傅里叶域的横轴坐标，确定振动信号在分数阶域内的聚集中心位置，以聚集中心位置为中心在分数阶域内对振动信号分数阶傅里叶变换结果Xp(u)进行滤波；五.对滤波后的振动信号分数阶傅里叶变换结果Xp(u)进行分数阶傅里叶逆变换，获得时域下的振动信号滤波结果，对其进行包络解调，利用离散傅氏变换的快速算法变换得到振动信号特征频谱。较佳的，步骤一中采用一维小波变换计算振动信号的时频分布，以Morlet小波作为母小波，进行计算得到振动信号x(t)的时频分布结果。较佳的，步骤二中，所述切比雪夫范数意义下的最佳逼近函数m(t)的具体计算过程为：1)选取振动信号x(t)的时频分布V(f,t)幅值最高的三个点作为交错点，分别为V1(f1,t1)、V2(f2,t2)、V3(f3,t3)，其中t为时间变量，ti表示第i个交错点横坐标位置，f为频率变量，fi表示第i个交错点纵坐标位置，i＝1,2,3；2)最佳逼近函数m(t)取一次多项式，m(t)＝c0+c1*t，t为时间变量，c0为常数项系数，c1为一次项系数；3)设切比雪夫范数‖r(t)‖c＝α，以f1-m(t1)＝-α，f2-m(t2)＝α，f3-m(t3)＝-α三个方程组成方程组，计算得到最佳逼近函数m(t)的常数项系数c0及一次项系数c1和α的值，即获得切比雪夫范数意义下的最佳逼近函数m(t)。较佳的，步骤三中，利用最佳逼近函数m(t)的一次项系数c1，确定分数阶傅里叶变换的最佳阶次p，较佳的，步骤四中，将计算得到的最佳阶次p的数值增加2后除以4得到余数Re，根据余数数值确定对应的分数阶傅里叶变换结果。较佳的，当Re＝2时，分数阶傅里叶变换结果为Xp(u)＝x(t)；当Re＝0时，分数阶傅里叶变换结果为Xp(u)＝x(-t)；当Re＝3时，分数阶傅里叶变换为普通傅里叶变换；Re＝1时，分数阶傅里叶变换为普通傅里叶逆变换；Re为其余值时，分数阶傅里叶变换结果Xp(u)为其中，β为旋转角度，j为虚数单位。较佳的，在对振动信号进行分数阶傅里叶变换后，取振动信号分数阶傅里叶变换结果Xp(u)幅值的绝对值，在分数阶域内确定幅值最高点的位置，以该位置为振动信号在分数阶域内的聚集中心位置。较佳的，步骤五中，对滤波后的振动信号分数阶傅里叶变换结果Xp(u)进行分数阶傅里叶逆变换，获得时域下的振动信号滤波结果，随后对时域下的信号滤波结果进行经验模态分解，对分解得到的固有模态分量进行包络解调，然后利用离散傅氏变换的快速算法变换得到振动信号频谱。为解决上述技术问题，本专利技术提供的一种根据上述振动信号的特征频谱提取方法得到的振动信号频谱进行机械故障诊断分析方法，所述机械为旋转机械。本专利技术的振动信号的特征频谱提取方法，基于分数阶傅里叶变换自适应地对振动信号进行滤波，然后通过包络分析得到振动信号特征频谱，可以提高振动信号故障诊断的准确性，改善包络分析方法在旋转结构故障诊断的应用效果，有助于旋转机械故障诊断分析。附图说明为了更清楚地说明本专利技术的技术方案，下面对本专利技术所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1是本专利技术的时域振动信号的特征频谱提取方法一实施例的流程图；图2是仿真信号x(t)的时域波形图；图3是对仿真信号未滤波直接进行经验模态分解的固有模态分量频谱图；图4是图3的固有模态分量包络解调后的频谱图；图5是对仿真信号经过滤波后处理后进行经验模态分解的固有模态分量频谱图；图6是图5的固有模态分量包络解调后的频谱图。具体实施方式下面将结合附图，对本专利技术中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本专利技术的一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其它实施例，都属于本专利技术保护的范围。实施例一如图1所示，振动信号的特征频谱提取方法包括以下步骤：一.计算振动信号x(t)的时频分布V(f,t)；二.在时频域中，根据交错定理获得切比雪夫范数意义下的最佳逼近函数m(t)；三.利用最佳逼近函数，确定分数阶傅里叶变换(FRFT)的最佳阶次p；四.计算最佳阶次下的振动信号x(t)的分数阶傅里叶变换结果Xp(u)，其中u为分数阶傅里叶域的横轴坐标，确定振动信号在分数阶域内的聚集中心位置，以聚集中心位置为中心在分数阶域内对振动信号分数阶傅里叶变换结果Xp(u)进行滤波；五.对滤波后的振动信号分数阶傅里叶变换结果Xp(u)进行分数阶傅里叶逆变换，获得时域下的振动信号滤波结果，对其进行包络解调，利用FFT(离散傅氏变换的快速算法)变换得到振动信号特征频谱。实施例一的振动信号的特征频谱提取方法，基于分数阶傅里叶变换自适应地对振动信号进行滤波，然后通过包络分析得到振动信号特征频谱，可以提高振动信号故障诊断的准确性，改善包络分析方法在旋转结构故障诊断的应用效果，有助于旋转机械故障诊断分析。实施例二基于实施例一的振动信号的特征频谱提取方法，步骤一中采用一维小波变换计算振动信号的时频分布，以Morlet小波作为母小波，进行计算得到振动信号x(t)的时频分布结果。实施例三基于实施例一的振动信号的特征频谱提取方法，步骤二中，所述切比雪夫范数意本文档来自技高网...

【技术保护点】
1.一种振动信号的特征频谱提取方法，其特征在于，包括以下步骤：/n一.计算振动信号x(t)的时频分布V(f，t)；/n二.在时频域中，根据交错定理获得切比雪夫范数意义下的最佳逼近函数m(t)；/n三.利用最佳逼近函数，确定分数阶傅里叶变换的最佳阶次p；/n四.计算最佳阶次下的振动信号x(t)的分数阶傅里叶变换结果X

【技术特征摘要】
1.一种振动信号的特征频谱提取方法，其特征在于，包括以下步骤：
一.计算振动信号x(t)的时频分布V(f，t)；
二.在时频域中，根据交错定理获得切比雪夫范数意义下的最佳逼近函数m(t)；
三.利用最佳逼近函数，确定分数阶傅里叶变换的最佳阶次p；
四.计算最佳阶次下的振动信号x(t)的分数阶傅里叶变换结果Xp(u)，其中u为分数阶傅里叶域的横轴坐标，确定振动信号在分数阶域内的聚集中心位置，以聚集中心位置为中心在分数阶域内对振动信号分数阶傅里叶变换结果Xp(u)进行滤波；
五.对滤波后的振动信号分数阶傅里叶变换结果Xp(u)进行分数阶傅里叶逆变换，获得时域下的振动信号滤波结果，对其进行包络解调，利用离散傅氏变换的快速算法变换得到振动信号特征频谱。


2.根据权利要求1所述的振动信号的特征频谱提取方法，其特征在于，
步骤一中采用一维小波变换计算振动信号的时频分布，以Morlet小波作为母小波，进行计算得到振动信号x(t)的时频分布结果。


3.根据权利要求1所述的振动信号的特征频谱提取方法，其特征在于，
步骤二中，所述切比雪夫范数意义下的最佳逼近函数m(t)的具体计算过程为：
1)选取振动信号x(t)的时频分布V(f，t)幅值最高的三个点作为交错点，分别为V1(f1，t1)、V2(f2，t2)、V3(f3，t3)，其中t为时间变量，ti表示第i个交错点横坐标位置，f为频率变量，fi表示第i个交错点纵坐标位置，i＝1，2，3；
2)最佳逼近函数m(t)取一次多项式，m(t)＝c0+c1*t，t为时间变量，c0为常数项系数，c1为一次项系数；
3)设切比雪夫范数||r(t)||c＝α，以f1-m(t1)＝-α，f2-m(t2)＝α，f3-m(t3)＝-α三个方程组成方程组，计算得到最佳逼近函数m(t)的常数项系数c0及一次项系数c1和α的值...

【专利技术属性】
技术研发人员：郑斌，蔡晓华，马骧越，俞英杰，
申请(专利权)人：上海三菱电梯有限公司，
类型：发明
国别省市：上海;31