[发明专利]基于主成份分析的发动机多参数使用相关载荷谱仿真方法在审

专利信息
申请号: 201610099334.8 申请日: 2016-02-23
公开(公告)号: CN105760610A 公开(公告)日: 2016-07-13
发明(设计)人: 孙志刚;邢广鹏;宋迎东 申请(专利权)人: 南京航空航天大学
主分类号: G06F17/50 分类号: G06F17/50
代理公司: 南京瑞弘专利商标事务所(普通合伙) 32249 代理人: 黄成萍
地址: 210016*** 国省代码: 江苏;32
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明公开了一种基于主成份分析的发动机多参数使用相关载荷谱仿真方法,通过采用主成份分析法建立了航空发动机多参数载荷谱的数学模型,并利用该模型对正态、非正态、实测航空发动机相关参数载荷谱进行仿真与模拟。本发明提出的基于主成份分析的多参数相关载荷谱仿真方法可通过较少的测量子样,模拟出任意时间段多参数使用载荷谱,结果表明仿真谱与原谱的穿级计数累积频次具有较好的一致性,且保证了各载荷参数间的相关性不变。
搜索关键词: 基于 成份 分析 发动机 参数 使用 相关 载荷 仿真 方法
【主权项】:
一种基于主成份分析的发动机多参数使用相关载荷谱仿真方法,其特征在于:设总共有p个载荷参数,n个采样点,将第j个载荷参数在第i个采样点的观测值记为xij,则将原始载荷数据矩阵记为X,<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>该方法包括如下步骤:(1)将原始载荷数据矩阵X转换为标准化载荷数据矩阵<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>(2)根据标准化载荷数据矩阵X′,计算标准化载荷数据矩阵X′列之间的相关系数矩阵<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>(3)计算相关系数矩阵R的特征值和对应每个特征值的特征向量,将所有特征向量组成的特征向量矩阵记为T,<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>(4)计算标准化载荷数据矩阵X′的主成分向量Y,Y=[y1 y2 … yp],yj为标准化载荷数据矩阵X′的第j个主成分,Y=X′T;(5)采用单参数载荷谱的仿真方法对第j个主成分yj进行单独仿真,得到第j个主成分谱y′j,将主成分谱向量记为Y′,Y′=[y′1 y′2 … y′p];(6)根据x″j=T‑1y′j,对主成分谱y′j进行逆变换,得到标准化载荷数据矩阵X′的仿真谱X″,X″=[x″1 x″2 … x″p];(7)根据将标准化载荷数据矩阵X′的仿真谱X″转换为原始载荷数据矩阵X的仿真谱X″′,X″′=[x″′1 x″′2 … x″′p],为原始载荷数据矩阵X的第j列数据的平均值,σj为原始载荷数据矩阵X的第j列数据的标准差,<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>
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于南京航空航天大学,未经南京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201610099334.8/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top