[發(fā)明專利]基于Landweber迭代法的運行工況傳遞路徑分析方法

在審

• 申請?zhí)枺?/td>: 201810283791.1
• 申請日: 2018-04-02
• 公開（公告）號: CN108536946A
• 公開（公告）日: 2018-09-14
• 發(fā)明（設(shè)計）人: 成瑋;朱巖;陳雪峰;褚亞鵬;倪晶磊;陸建濤;加正正;楊志勃
• 申請（專利權(quán)）人: 西安交通大學
• 主分類號: G06F17/50
• 分類號: G06F17/50
• 代理公司: 西安通大專利代理有限責任公司 61200
• 代理人: 徐文權(quán)

1.基于Landweber迭代法的運行工況傳遞路徑分析方法，其特征在于，包括如下步驟：

步驟1，選取所分析機械系統(tǒng)激勵點的臨近位置作為參考點，選取待分析位置為目標點；

步驟2，根據(jù)所分析機械系統(tǒng)設(shè)計若干種工況試驗，使機械系統(tǒng)在每種試驗工況下運行，采集每種試驗工況下參考點和目標點響應信號，得到參考點響應信號矩陣和目標點響應信號矩陣，并對其進行快速傅里葉變換得到相應的頻響函數(shù)矩陣X和Y；

步驟3，設(shè)置迭代初值X₀＝αX^*，計算迭代系數(shù)矩陣X_k＝X_k-1+X₀(I-XX_k-1)，在傳遞率函數(shù)矩陣T_k＝X_kY滿足Morozov誤差原則前持續(xù)迭代系數(shù)矩陣X_k，直至迭代得到的傳遞率函數(shù)矩陣T_k滿足Morozov誤差原則，停止迭代；

步驟4，根據(jù)所分析機械系統(tǒng)進行實際運行工況實驗，得到一組實際運行工況數(shù)據(jù)，根據(jù)該組實際運行工況數(shù)據(jù)和迭代得到的傳遞率函數(shù)矩陣計算得到估計的目標點總貢獻量和各路徑貢獻量

步驟5，針對各條傳遞路徑對各路徑貢獻量進行排序后計算出各傳遞路徑的貢獻量占比，完成運行工況傳遞路徑分析。

2.根據(jù)權(quán)利要求1所述的基于Landweber迭代法的運行工況傳遞路徑分析方法，其特征在于，步驟2、步驟4中均采用振動加速度傳感器測量響應信號。

3.根據(jù)權(quán)利要求1所述的基于Landweber迭代法的運行工況傳遞路徑分析方法，其特征在于，步驟2中試驗工況數(shù)量大于等于參考點數(shù)量，步驟4中實際運行工況數(shù)量大于等于參考點數(shù)量。

4.根據(jù)權(quán)利要求1所述的基于Landweber迭代法的運行工況傳遞路徑分析方法，其特征在于，步驟2中X(ω)＝[x₁(ω),x₂(ω),...,x_r(ω)]為全部試驗工況下所有參考點響應信號，r表示試驗工況數(shù)量，表示第i種試驗工況下所有參考點響應信號，表示試驗工況i下參考點j的振動響應信號，m表示參考點數(shù)量。

5.根據(jù)權(quán)利要求1所述的基于Landweber迭代法的運行工況傳遞路徑分析方法，其特征在于，步驟2中Y(ω)＝[y₁(ω),y₂(ω),...,y_r(ω)]為全部試驗工況下所有目標點響應信號，r表示試驗工況數(shù)量，表示第i種試驗工況下所有目標點響應信號，表示試驗工況i下目標點g的振動響應信號，n表示目標點數(shù)量。

6.根據(jù)權(quán)利要求1所述的基于Landweber迭代法的運行工況傳遞路徑分析方法，其特征在于，步驟3中迭代初值X₀＝αX^*，α為松弛因子，且0<α<1/||X||²，X^*表示參考點頻率響應函數(shù)X的共軛轉(zhuǎn)置。

7.根據(jù)權(quán)利要求1所述的基于Landweber迭代法的運行工況傳遞路徑分析方法，其特征在于，步驟3中計算迭代系數(shù)矩陣X_k＝X_k-1+X₀(I-XX_k-1)，I為單位矩陣，k表示迭代次數(shù)。

8.根據(jù)權(quán)利要求1所述的基于Landweber迭代法的運行工況傳遞路徑分析方法，其特征在于，步驟3中Morozov誤差原則為：當?shù)玫降膫鬟f率函數(shù)矩陣滿足||X(ω)T_k(ω)-Y(ω)||≤τδ時停止迭代，常數(shù)τ＝2，誤差水平δ＝10^-3。

9.根據(jù)權(quán)利要求1所述的基于Landweber迭代法的運行工況傳遞路徑分析方法，其特征在于，步驟4中表示OTPA識別的參考點j對目標點i的傳遞路徑貢獻量。