本发明属于导航设备测试,更具体地,涉及一种基于大地线递推的近地表轨迹重构方法和系统。
背景技术:
1、由于惯性导航设备在移交应用前,需要进行精度测试以评估其性能是否达标,而地球两极所处的地理位置不便于中国等中低纬度国家前往进行实测试验,且由于惯导极区编排算法和工作模式与中低纬度的不同,不能省略这一阶段的验证。
2、因此,前人研究了虚拟极区技术,将中低纬度地区的实测试验数据通过数学方法转换至极区进行等效的导航试验验证。轨迹重构过程是虚拟极区技术中的关键环节之一,目的是在极区规划一条与实测试验相似的运动轨迹。
3、专利cn2019107309381中公开了一种基于横向地理坐标系的轨迹重构方法,通过采取保持横向地理系下的姿态矩阵、速度和高度信息不变的轨迹重构原则,对补偿椭球校正系数后的横向速度积分重构极区位置轨迹。
4、然而,该方法存在以下缺陷和不足:实际使载体相对于当地水平面的运动特征发生改变,有违设计初衷,也造成了试验的额外误差,为后续的精度评估增加了困难。此外,现有技术也不能对重构轨迹的位置和方向灵活设置,局限性较大。
技术实现思路
1、针对现有技术的缺陷,本发明的目的在于提供一种基于大地线递推的近地表轨迹重构方法和系统,旨在解决如何在保持载体对地运动不变的条件下进行轨迹重构的问题。
2、为实现上述目的,第一方面,本发明提供了一种基于大地线递推的近地表轨迹重构方法,由模拟地区第i点向第i+1点的递推包括:
3、step1:计算实测[i,i+1]时段的平均纬度b=(bi+bi+1)/2,进而计算实测当地子午圈半径m和实测当地卯酉圈半径n;
4、step2:计算该时段航迹的大地线方位角变化dk和相邻点间大地线的投影长度si,i+1:
5、dk=sinb×dl
6、
7、step3:计算cosk,并根据db和dl的符号判断k所处的象限,进而反演出平均方位角k:
8、cosk=mdb/si,i+1
9、step4:计算起始方位角ki,i+1,s、终止方位角ki,i+1,e和模拟轨迹该时段初始方位角
10、
11、
12、
13、step5:计算模拟轨迹相邻点间大地线的投影长度
14、
15、step6:根据i时刻模拟位置的纬度计算模拟地区子午圈半径m*和模拟地区卯酉圈半径n*;
16、step7:计算模拟轨迹该时段纬度增量db*、模拟轨迹该时段经度增量dl*和模拟轨迹该时段的大地线方位角变化dk*:
17、
18、
19、
20、step9:计算载体在i+1时刻的模拟位置和模拟轨迹终止方位角
21、
22、
23、
24、其中,该时段纬度增量db=bi+1-bi,该时段经度增量dl=li+1-li;r表示实测轨迹段所处的地表曲率半径;r*表示模拟轨迹段所处的地表曲率半径;表示i时刻模拟位置的经度,hi+1表示载体在i+1时刻大地高度。
25、优选地,step7和step9之间包括:
26、step8:用分别替换step5~step7中的和重复step5~step7,直至满足精度要求。
27、优选地,实测轨迹段所处的地表曲率半径
28、优选地,模拟轨迹段所处的地表曲率半径
29、为实现上述目的,第二方面,本发明提供了一种基于大地线递推的近地表轨迹重构系统,包括:处理器和存储器;所述存储器,用于存储计算机执行指令;所述处理器,用于执行所述计算机执行指令,使得如第一方面所述的方法被执行。
30、总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
31、本发明提出一种基于大地线递推的近地表轨迹重构方法和系统,通过提取实测定位数据中相邻点间大地线的投影长度(si,i+1)、各定位点处大地线转角(ki,i+1,s-ki-1,i,e)和高程变化信息(hi+1-hi),采用保持对地航向的变化和对地高度不变的原则,由人为设定模拟区域的初始地点(包括h0)和初始方向开始,进行模拟航迹递推,完成实测轨迹在地表任意区域的重构。本发明经过数学证明满足重构轨迹与实测具有相同的对地运动特征,为后续评估结果的认定提供了严谨的理论依据;且本发明可适用于任意地点、任意方向的轨迹重构,增加了试验设计的灵活性。
1.一种基于大地线递推的近地表轨迹重构方法,其特征在于,由模拟地区第i点向第i+1点的递推包括:
2.如权利要求1所述的方法,其特征在于,step7和step9之间包括:
3.如权利要求1或2所述的方法,其特征在于,实测轨迹段所处的地表曲率半径
4.如权利要求1或2所述的方法,其特征在于,模拟轨迹段所处的地表曲率半径
5.一种基于大地线递推的近地表轨迹重构系统,其特征在于,包括:处理器和存储器;