主权项 |
1.基于非负张量分解振动特征的钢轨伤损探测方法,其特征是它包括具体步骤如下:步骤一:振动加速度传感器阵列(1)采集到被测钢轨的原始振动信号S(t);步骤二:信号转换电路(2)转换原始振动信号S(t)为振动信号x(t);步骤三:将振动信号x(t)进行希尔伯特-黄变换分解,得到各阶内固模态函数和残差,记为:c<sub>1</sub>,c<sub>2</sub>,...c<sub>n</sub>以及r,然后进行相应的希尔伯特变换,得到频率与幅值之间的幅频信号h(w);步骤四:利用希尔伯特-黄变换分析所得到频率幅值之间的幅频信号h(w),结合相应的测量点和伤损种类建立伤损信号的三维张量;步骤五:利用非负张量分解对构建的三维张量信号x<sup>frequency×point×class</sup>进行分解,提取相应的伤损特征信息,得到每种典型伤损情况相应的特征向量,实测振动信号的特征向量通过与典型伤损情况特征向量的对比,判断出伤损的存在及其位置;步骤五包括的子步骤如下:通过步骤四建立的三维张量信号表示为:X<sup>frequency×point×class</sup>=G×<sub>1</sub>A<sup>(1)</sup>×<sub>2</sub>A<sup>(2)</sup>×<sub>3</sub>A<sup>(3)</sup>其中G,A<sup>(1)</sup>,A<sup>(2)</sup>,A<sup>(3)</sup>≥0,利用非负张量分解对其进行分解;子步骤1、随机初始化A<sup>(n)</sup>,n=1,2,3;算最小化的最小均方误差C<sub>new</sub>=‖X<sub>(n)</sub>-A<sup>(n)</sup>Z<sup>(n)T</sup>‖<sup>2</sup>,其中:X<sub>(n)</sub>为原始信号,<maths num="0001"><![CDATA[<math><mrow><msup><mi>Z</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup><mo>=</mo><msup><mi>A</mi><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></msup><mo>|</mo><mo>⊗</mo><mo>|</mo><msup><mi>A</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>|</mo><mo>⊗</mo><mo>|</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>|</mo><mo>⊗</mo><mo>|</mo><msup><mi>A</mi><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>|</mo><mo>⊗</mo><mo>|</mo><msup><mi>A</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>|</mo><mo>⊗</mo><mo>|</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>|</mo><mo>⊗</mo><mo>|</mo><msup><mi>A</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>,</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>3</mn><mo>;</mo></mrow></math>]]></maths>子步骤2、迭代求解<img file="FDA00002530446900012.GIF" wi="491" he="121" />并计算C<sub>new</sub>=‖X<sub>(n)</sub>-A<sup>(n)</sup>Z<sup>(n)T</sup>‖<sup>2</sup>,n=1,2,3;子步骤3、判断<img file="FDA00002530446900013.GIF" wi="272" he="103" />δ为误差精度要求,成立则C<sub>old</sub>=C<sub>new</sub>,继续执行子步骤2;不成立则执行子步骤4;子步骤4、得到满足要求的A<sup>(1)</sup>,A<sup>(2)</sup>,A<sup>(3)</sup>矩阵;其中:A<sup>(1)</sup>,A<sup>(2)</sup>为新的基底,其物理意义分别代表频率和测量点;A<sup>(3)</sup>为伤损类别的特征系数矩阵,其中每一行是一个特征向量,对应一种伤损情况;通过与典型伤损情况特征向量的对比,判断出伤损的存在及其位置。 |