发明名称 基于主成份分析的发动机多参数使用相关载荷谱仿真方法
摘要 本发明公开了一种基于主成份分析的发动机多参数使用相关载荷谱仿真方法,通过采用主成份分析法建立了航空发动机多参数载荷谱的数学模型,并利用该模型对正态、非正态、实测航空发动机相关参数载荷谱进行仿真与模拟。本发明提出的基于主成份分析的多参数相关载荷谱仿真方法可通过较少的测量子样,模拟出任意时间段多参数使用载荷谱,结果表明仿真谱与原谱的穿级计数累积频次具有较好的一致性,且保证了各载荷参数间的相关性不变。
申请公布号 CN105760610A 申请公布日期 2016.07.13
申请号 CN201610099334.8 申请日期 2016.02.23
申请人 南京航空航天大学 发明人 孙志刚;邢广鹏;宋迎东
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 黄成萍
主权项 一种基于主成份分析的发动机多参数使用相关载荷谱仿真方法,其特征在于:设总共有p个载荷参数,n个采样点,将第j个载荷参数在第i个采样点的观测值记为x<sub>ij</sub>,则将原始载荷数据矩阵记为X,<maths num="0001"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>x</mi><mn>11</mn></msub></mtd><mtd><msub><mi>x</mi><mn>12</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>x</mi><mrow><mn>1</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>21</mn></msub></mtd><mtd><msub><mi>x</mi><mn>22</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>x</mi><mrow><mn>2</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>x</mi><mrow><mi>n</mi><mi>p</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000927682430000011.GIF" wi="541" he="303" /></maths>该方法包括如下步骤:(1)将原始载荷数据矩阵X转换为标准化载荷数据矩阵<maths num="0002"><math><![CDATA[<mrow><msup><mi>X</mi><mo>&prime;</mo></msup><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>x</mi><mn>11</mn><mo>&prime;</mo></msubsup></mtd><mtd><msubsup><mi>x</mi><mn>12</mn><mo>&prime;</mo></msubsup></mtd><mtd><mn>...</mn></mtd><mtd><msubsup><mi>x</mi><mrow><mn>1</mn><mi>p</mi></mrow><mo>&prime;</mo></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mn>21</mn><mo>&prime;</mo></msubsup></mtd><mtd><msubsup><mi>x</mi><mn>22</mn><mo>&prime;</mo></msubsup></mtd><mtd><mn>...</mn></mtd><mtd><msubsup><mi>x</mi><mrow><mn>2</mn><mi>p</mi></mrow><mo>&prime;</mo></msubsup></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mn>1</mn></mrow><mo>&prime;</mo></msubsup></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mn>2</mn></mrow><mo>&prime;</mo></msubsup></mtd><mtd><mn>...</mn></mtd><mtd><msubsup><mi>x</mi><mrow><mi>n</mi><mi>p</mi></mrow><mo>&prime;</mo></msubsup></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000927682430000012.GIF" wi="547" he="310" /></maths>(2)根据标准化载荷数据矩阵X′,计算标准化载荷数据矩阵X′列之间的相关系数矩阵<maths num="0003"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>r</mi><mn>11</mn></msub></mtd><mtd><msub><mi>r</mi><mn>12</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>r</mi><mrow><mn>1</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>r</mi><mn>21</mn></msub></mtd><mtd><msub><mi>r</mi><mn>22</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>r</mi><mrow><mn>2</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>r</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>r</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>r</mi><mrow><mi>p</mi><mi>p</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000927682430000013.GIF" wi="507" he="311" /></maths>(3)计算相关系数矩阵R的特征值和对应每个特征值的特征向量,将所有特征向量组成的特征向量矩阵记为T,<maths num="0004"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>u</mi><mn>11</mn></msub></mtd><mtd><msub><mi>u</mi><mn>12</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>u</mi><mrow><mn>1</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mn>21</mn></msub></mtd><mtd><msub><mi>u</mi><mn>22</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>u</mi><mrow><mn>2</mn><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>u</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>u</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>u</mi><mrow><mi>p</mi><mi>p</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000927682430000014.GIF" wi="534" he="303" /></maths>(4)计算标准化载荷数据矩阵X′的主成分向量Y,Y=[y<sub>1</sub> y<sub>2</sub> … y<sub>p</sub>],y<sub>j</sub>为标准化载荷数据矩阵X′的第j个主成分,Y=X′T;(5)采用单参数载荷谱的仿真方法对第j个主成分y<sub>j</sub>进行单独仿真,得到第j个主成分谱y′<sub>j</sub>,将主成分谱向量记为Y′,Y′=[y′<sub>1</sub> y′<sub>2</sub> … y′<sub>p</sub>];(6)根据x″<sub>j</sub>=T<sup>‑1</sup>y′<sub>j</sub>,对主成分谱y′<sub>j</sub>进行逆变换,得到标准化载荷数据矩阵X′的仿真谱X″,X″=[x″<sub>1</sub> x″<sub>2</sub> … x″<sub>p</sub>];(7)根据<img file="FDA0000927682430000015.GIF" wi="313" he="78" />将标准化载荷数据矩阵X′的仿真谱X″转换为原始载荷数据矩阵X的仿真谱X″′,X″′=[x″′<sub>1</sub> x″′<sub>2</sub> … x″′<sub>p</sub>],<img file="FDA0000927682430000021.GIF" wi="53" he="78" />为原始载荷数据矩阵X的第j列数据的平均值,<img file="FDA0000927682430000022.GIF" wi="277" he="127" />σ<sub>j</sub>为原始载荷数据矩阵X的第j列数据的标准差,<maths num="0005"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>j</mi></msub><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><mover><msub><mi>x</mi><mi>j</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>.</mo></mrow>]]></math><img file="FDA0000927682430000023.GIF" wi="525" he="150" /></maths>
地址 210016 江苏省南京市御道街29号