发明名称 飞行器导航制导与控制地面仿真系统效能评估方法
摘要 本发明涉及一种飞行器导航制导与控制地面仿真系统效能评估方法,先根据系统功能原理选择相似元,然后确定各相似元的数值和权系数,最后计算仿真系统的相似度。本方法将效能评估问题转化为可行度问题,从动力学、运动学、测量以及导航制导与控制等方面定量分析“天地一致性”问题,基于可信度理论,客观量化分析地面仿真系统的可信度。本发明提供的评估方法原理清晰、计算简单,具有比较广泛的适用范围,可以为其他仿真系统的效能分析提供参考借鉴。
申请公布号 CN104077490A 申请公布日期 2014.10.01
申请号 CN201410326380.8 申请日期 2014.07.03
申请人 哈尔滨工业大学 发明人 李莉;夏红伟;刘飞;马广程;王常虹;李秋实
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 代理人
主权项 一种飞行器导航制导与控制地面仿真系统效能评估方法,先根据系统功能原理选择相似元,然后确定各相似元的数值和权系数,最后计算仿真系统的相似度,其特征在于,具体步骤如下:第1步:选择相似元飞行器导航制导与控制地面仿真系统由运动模拟系统、测量系统、动力学仿真系统和GNC控制系统四部分构成,因此,评价地面仿真系统与实际空间飞行器运行状况的相似度,通过这四个分系统与空间飞行器实际工作情况的相似程度来确定,所以选择这四个分系统与实际情况的相似程度作为相似元;第2步:确定相似元的值(2.1)、运动模拟系统运动模拟系统能够完成X、Y、Z三个平动方向的运动和俯仰、偏航、滚转三个转动方向的运动,运动模拟系统与空间飞行器实际工作情况的相似程度通过运动模拟器的运行精度参数确定;仿真缩比值为K,运动模拟器的位置精度为J1、速度精度J2、动态精度J3,实际卫星运动过程也会产生偏差,如果实际卫星运动系统偏差大于地面仿真系统运动模拟器的偏差值,则认为运动模拟器能够完全模拟实际卫星运动,即相似元为1,在不知道实际卫星运动偏差值的情况下,设实际卫星运动完全按理论值运动,则运动模拟器相对实际卫星最大偏差取为J1、J2和J3中的最大值,相似元的值为Q<sub>1</sub>=1‑max{J1、J2、J3},其中max{...}表示取最大值;(2.2)、测量系统对于测量系统,设Wc和Fc分别表示实际飞行器测量系统的误差和仿真系统的测量误差,并且仿真缩比为K,如果满足K×Fc≤Wc,认为地面测量系统能够完全模拟卫星测量系统,则相似元定为1;否则按照公式<img file="FSA0000106180990000011.GIF" wi="439" he="131" />计算相似元的值;(2.3)动力学仿真系统仿真时,动力学仿真系统对空间扰动力矩进行数学建模并进行Matlab/Simulink仿真,在RTW下自动生成C代码,并经过修改后生成供LabVIEW直接调用的动态链接库,在LabVIEW RT实时系统中运行,对于这部分相似元值计算,通过与STK仿真系统仿真得到的数据进行比较得到结果,计算方法为:<img file="FSA0000106180990000021.GIF" wi="432" he="146" />其中X<sub>D</sub>表示动力学计算出来的卫星的位置和速度,X<sub>STK</sub>表示STK计算出来的卫星的位置和速度;(2.3)GNC控制模块仿真时GNC控制模块采用的制导控制率与航天器实际采用的控制率相同,故相似元值直接取为Q<sub>4</sub>=1;第3步:确定权系数根据系统差别对仿真结果的影响程度分析,这四个系统相似性的重要程度依次为测量系统、运动系统、动力学运动学系统、GNC控制系统,所以它们的权重影响程度成单调减关系,根据层次分析法得判断矩阵P:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>2</mn></mtd><mtd><mn>3</mn></mtd><mtd><mn>6</mn></mtd></mtr><mtr><mtd><mn>1</mn><mo>/</mo><mn>2</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>2</mn></mtd><mtd><mn>3</mn></mtd></mtr><mtr><mtd><mn>1</mn><mo>/</mo><mn>3</mn></mtd><mtd><mn>1</mn><mo>/</mo><mn>2</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>2</mn></mtd></mtr><mtr><mtd><mn>1</mn><mo>/</mo><mn>6</mn></mtd><mtd><mn>1</mn><mo>/</mo><mn>3</mn></mtd><mtd><mn>1</mn><mo>/</mo><mn>2</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000106180990000022.GIF" wi="1139" he="304" /></maths>求出最大特征值λ<sub>max</sub>=4.0104;一致性指标CI=0.0035,平均随机一致性指数RI<sub>(4)</sub>=0.90,于是一致性指标CR=CI/RI<sub>(4)</sub>=0.0039<0.1;故认为判断矩阵一致性是能够接受的;将判断矩阵P按列归一化得:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0.5000</mn></mtd><mtd><mn>0.5217</mn></mtd><mtd><mn>0.4615</mn></mtd><mtd><mn>0.5000</mn></mtd></mtr><mtr><mtd><mn>0.2500</mn></mtd><mtd><mn>0.2609</mn></mtd><mtd><mn>0.3077</mn></mtd><mtd><mn>0.2500</mn></mtd></mtr><mtr><mtd><mn>0.1667</mn></mtd><mtd><mn>0.1304</mn></mtd><mtd><mn>0.1538</mn></mtd><mtd><mn>0.1667</mn></mtd></mtr><mtr><mtd><mn>0.0833</mn></mtd><mtd><mn>0.0870</mn></mtd><mtd><mn>0.0769</mn></mtd><mtd><mn>0.0833</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000106180990000023.GIF" wi="1222" he="291" /></maths>将列元素归一化后的矩阵按行相加并将列向量归一化处理,得:β=[β<sub>1</sub> β<sub>2</sub> β<sub>3</sub> β<sub>4</sub>]<sup>T</sup>=[0.4959 0.2671 0.1544 0.0826]<sup>T</sup>      (3‑3)其中β<sub>1</sub>,β<sub>2</sub>,β<sub>3</sub>,β<sub>4</sub>即为各个相似元的权重系数;第4步:计算仿真系统的相似度由于系统间对应的元素看成相似元,则相似度为:<img file="FSA0000106180990000025.GIF" wi="268" he="124" />这样就得到了仿真系统的相似度Q。
地址 150001 黑龙江省哈尔滨市南岗区一匡街2号E2栋