主权项 |
1.一种利用矩阵小波变换的行星齿轮箱复合故障诊断方法,其特征在于,包括下述步骤:(1)矩阵小波变换由多个尺度函数生成向量小波函数,该多个尺度函数和所述向量小波函数共同构成矩阵小波函数,其中,矩阵小波函数由如下多分辨分析生成,设函数Φ=[φ<sub>1</sub>(t),φ<sub>2</sub>(t),…,φ<sub>r</sub>(t)]<sup>T</sup>∈L<sup>2</sup>(R)<sup>r</sup>,r∈N,其中φ<sub>1</sub>(t),φ<sub>2</sub>(t),…,φ<sub>r</sub>(t)为多个尺度函数,对j∈Z,定义:V<sub>j</sub>=clos{2<sup>j/2</sup>φ<sub>i</sub>(2<sup>j</sup>t-k):1≤i≤r,k∈Z}式中,尺度函数空间V<sub>j</sub>由闭包函数clos{ }张成,分辨率为2<sup>j</sup>;小波函数空间W<sub>j</sub>=clos{2<sup>j/2</sup>ψ<sub>i</sub>(2<sup>j</sup>t-k):1≤i≤r,k∈Z},是尺度函数空间V<sub>j</sub>在V<sub>j+1</sub>中的补空间,向量小波函数Ψ=[ψ<sub>1</sub>(t),ψ<sub>2</sub>(t),…,ψ<sub>r</sub>(t)]<sup>T</sup>∈L<sup>2</sup>(R)<sup>r</sup>的伸缩和平移张成了W<sub>j</sub>空间矩阵小波函数的两尺度方程为:<maths num="0001"><![CDATA[<math><mrow><mi>Φ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>Σ</mi><mi>k</mi></munder><msub><mi>H</mi><mi>k</mi></msub><mi>Φ</mi><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>∈</mo><mi>Z</mi></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>Ψ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>Σ</mi><mi>k</mi></munder><msub><mi>G</mi><mi>k</mi></msub><mi>Φ</mi><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>∈</mo><mi>Z</mi></mrow></math>]]></maths>式中,H<sub>k</sub>,G<sub>k</sub>分别表示矩阵小波变换的低通与高通滤波器,Ψ(t)是由尺度函数Φ(t)生成的向量小波函数;(2)将下述提升方法引入步骤(1)矩阵小波函数构造中,通过设计不同的提升算子改变原有矩阵小波函数滤波器的特性,得到不同性质的矩阵小波函数,所述提升方法为:给定一个初始矩阵小波滤波器组<img file="FDA00002188853300013.GIF" wi="537" he="90" />其中H(z)为低通滤波器,<img file="FDA00002188853300014.GIF" wi="117" he="70" />为对偶低通滤波器,G(z)为高通滤波器,<img file="FDA00002188853300015.GIF" wi="106" he="70" />为对偶高通滤波器;构造的新矩阵小波滤波器组<img file="FDA00002188853300016.GIF" wi="705" he="88" />为:H<sub>new</sub>(z)=H(z)G<sub>new</sub>(z)=T(z<sup>2</sup>)(G(z)+S(z<sup>2</sup>)H(z))<maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>H</mi><mo>~</mo></mover><mi>new</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>H</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>S</mi><mo>*</mo></msup><mrow><mo>(</mo><msup><mi>z</mi><mn>2</mn></msup><mo>)</mo></mrow><mover><mi>G</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>G</mi><mo>~</mo></mover><mi>new</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>T</mi><mo>*</mo></msup><mrow><mo>(</mo><msup><mi>z</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mover><mi>G</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中H<sub>new</sub>(z)为新的低通滤波器,<img file="FDA00002188853300019.GIF" wi="166" he="71" />为新的对偶低通滤波器,G<sub>new</sub>(z)为新的高通滤波器,<img file="FDA000021888533000110.GIF" wi="155" he="72" />为新的对偶高通滤波器;提升矩阵S(z)和T(z)均为有限阶,T(z)的行列式为单项式;(3)根据齿轮和轴承出现损伤将引起振动冲击以及故障状态确定的物理属性,利用信号处理的峭度最大和熵最小作为优化目标,在构造的矩阵小波函数库中选择具有优良性质的最优矩阵小波函数;(4)使用所构造的最优矩阵小波函数对振动信号进行自适应分解,将多个故障特征分解到不同分支的信号中进行提取与识别,实现行星齿轮箱复合故障特征的一次性分离和诊断。 |