主权项 |
1.一种基于改进的自适应形态滤波的非稳态信号检测方法,在待检测设备的壳体上安装加速度传感器,检测设备的振动加速度信号,作为检测信号,其特征在于实现步骤如下:(1)绘制原始信号的时域图,找出信号各局部极大值的位置,并计算信号局部极大值在时域图中与时间轴围成的轮廓区域面积为参考面积:<maths num="0001"><![CDATA[<math><mrow><msub><mi>S</mi><mi>ref</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>S</mi><mrow><mi>local</mi><mo>_</mo><mi>i</mi></mrow></msub></mrow></math>]]></maths>式中,S<sub>ref</sub>为原始信号极值与时间轴围成的轮廓区域面积,S<sub>local_i</sub>=0.5×length<sub>i</sub>×(max<sub>i</sub>+max<sub>i+1</sub>)为相邻两个局部极值与时间轴围成的轮廓区域面积,k为局部极值个数,i=1,2,...k-1,length<sub>i</sub>为局部极值间的距离,max<sub>i</sub>为第i个局部极值;(2)设定结构元素基本步长step,step为信号采样周期,即采样频率倒数,依次增加扫描各倍率长度,从1×step,2×step,…,m×step分别对原始信号进行形态滤波结果与时间轴围成的轮廓区域面积,并与步骤(1)中所得参考面积对比,找出两者面积最接近情况下采用的结构元素的步长,定为最佳结构元素;m为采样周期倍数;(3)根据步骤(2)中得到的最佳结构元素长度,根据下述计算公式利用闭-开算子对信号进行滤波,并分析处理后信号的频谱,判断故障类别:<img file="FDA00001601471600012.GIF" wi="889" he="83" />上式中,Θ,<img file="FDA00001601471600013.GIF" wi="66" he="82" />ο和·分别表示腐蚀算子、膨胀算子、开算子和闭算子,其中,基本腐蚀算子、膨胀算子、开算子、闭算子分别表达如下:腐蚀算子:(fΘg)(n)=min[f(n+m)-g(m)]膨胀算子:<maths num="0002"><![CDATA[<math><mrow><mrow><mo>(</mo><mi>f</mi><mo>⊕</mo><mi>g</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>max</mi><mo>[</mo><mi>f</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>m</mi><mo>)</mo></mrow><mo>+</mo><mi>g</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>开算子:<img file="FDA00001601471600015.GIF" wi="503" he="97" />闭算子:<maths num="0003"><![CDATA[<math><mrow><mrow><mo>(</mo><mi>f</mi><mo>·</mo><mi>g</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>f</mi><mo>⊕</mo><mi>gΘg</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></math>]]></maths>f(n)为待处理信号,长度为N,g(m)为结构元素,长度为M,m=0,1,2...M-1,n=0,1,2...N-1。 |