[發(fā)明專利]基于主成份分析的發(fā)動(dòng)機(jī)多參數(shù)使用相關(guān)載荷譜仿真方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201610099334.8
• 申請(qǐng)日: 2016-02-23
• 公開（公告）號(hào): CN105760610A
• 公開（公告）日: 2016-07-13
• 發(fā)明（設(shè)計(jì)）人: 孫志剛;邢廣鵬;宋迎東
• 申請(qǐng)（專利權(quán)）人: 南京航空航天大學(xué)
• 主分類號(hào): G06F17/50
• 分類號(hào): G06F17/50
• 代理公司: 南京瑞弘專利商標(biāo)事務(wù)所(普通合伙) 32249
• 代理人: 黃成萍

1.一種基于主成份分析的發(fā)動(dòng)機(jī)多參數(shù)使用相關(guān)載荷譜仿真方法，其特征在于：設(shè)總共有p個(gè)載荷參數(shù)，n個(gè)采樣點(diǎn)，將第j個(gè)載荷參數(shù)在第i個(gè)采樣點(diǎn)的觀測(cè)值記為x_ij，則將原始載荷數(shù)據(jù)矩陣記為X，X=x11x12...x1px21x22...x2p............xn1xn2...xnp;]]>

該方法包括如下步驟：

(1)將原始載荷數(shù)據(jù)矩陣X轉(zhuǎn)換為標(biāo)準(zhǔn)化載荷數(shù)據(jù)矩陣X′=x11′x12′...x1p′x21′x22′...x2p′............xn1′xn2′...xnp′;]]>

(2)根據(jù)標(biāo)準(zhǔn)化載荷數(shù)據(jù)矩陣X′，計(jì)算標(biāo)準(zhǔn)化載荷數(shù)據(jù)矩陣X′列之間的相關(guān)系數(shù)矩陣R=r11r12...r1pr21r22...r2p............rp1rp2...rpp;]]>

(3)計(jì)算相關(guān)系數(shù)矩陣R的特征值和對(duì)應(yīng)每個(gè)特征值的特征向量，將所有特征向量組成的特征向量矩陣記為T，T=u11u12...u1pu21u22...u2p............up1up2...upp;]]>

(4)計(jì)算標(biāo)準(zhǔn)化載荷數(shù)據(jù)矩陣X′的主成分向量Y，Y＝[y₁y₂…y_p]，y_j為標(biāo)準(zhǔn)化載荷數(shù)據(jù)矩陣X′的第j個(gè)主成分，Y＝X′T；

(5)采用單參數(shù)載荷譜的仿真方法對(duì)第j個(gè)主成分y_j進(jìn)行單獨(dú)仿真，得到第j個(gè)主成分譜y′_j，將主成分譜向量記為Y′，Y′＝[y′₁y′₂…y′_p]；

(6)根據(jù)x″_j＝T^-1y′_j，對(duì)主成分譜y′_j進(jìn)行逆變換，得到標(biāo)準(zhǔn)化載荷數(shù)據(jù)矩陣X′的仿真譜X″，X″＝[x″₁x″₂…x″_p]；

(7)根據(jù)將標(biāo)準(zhǔn)化載荷數(shù)據(jù)矩陣X′的仿真譜X″轉(zhuǎn)換為原始載荷數(shù)據(jù)矩陣X的仿真譜X″′，X″′＝[x″′₁x″′₂…x″′_p]，為原始載荷數(shù)據(jù)矩陣X的第j列數(shù)據(jù)的平均值，σ_j為原始載荷數(shù)據(jù)矩陣X的第j列數(shù)據(jù)的標(biāo)準(zhǔn)差，σj=1n-1Σi=1n(xij-xj‾)2.]]>