CN102564454B - 基于日地月方位信息的自主导航***的测量数据模拟方法 - Google Patents
基于日地月方位信息的自主导航***的测量数据模拟方法 Download PDFInfo
- Publication number
- CN102564454B CN102564454B CN201110441935.XA CN201110441935A CN102564454B CN 102564454 B CN102564454 B CN 102564454B CN 201110441935 A CN201110441935 A CN 201110441935A CN 102564454 B CN102564454 B CN 102564454B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- moon
- pulse
- earth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 19
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims abstract description 95
- 238000005259 measurement Methods 0.000 claims abstract description 27
- 230000009466 transformation Effects 0.000 claims description 7
- 230000010354 integration Effects 0.000 claims description 6
- 239000000523 sample Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000004576 sand Substances 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims description 2
- 238000012795 verification Methods 0.000 abstract description 4
- 230000008859 change Effects 0.000 abstract description 3
- 230000006835 compression Effects 0.000 abstract 1
- 238000007906 compression Methods 0.000 abstract 1
- 230000003287 optical effect Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 5
- 230000000007 visual effect Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 241001272720 Medialuna californiensis Species 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 238000010408 sweeping Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 241000959638 Geastrum saccatum Species 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004898 kneading Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Landscapes
- Navigation (AREA)
Abstract
基于日地月模拟测量信息的自主导航***的测量数据模拟方法,分别包括参考脉冲、地球脉冲、月球脉冲和太阳脉冲的模拟。模拟测量数据精确和全面地反映了各种复杂因素对日地月一体化敏感器测量的影响,这包括地球扁率、月球形状变化的影响。本发明实现了日地月一体化敏感器脉冲信息的模拟,可以有效地通过数值仿真的形式对敏感器的***误差进行定量的评估和对导航***的性能进行验证。
Description
技术领域
本发明涉及一种自主导航仿真试验方法,特别是一种基于日地月一体化敏感器测量信息的数值模拟方法,属于自主导航技术领域。
背景技术
自主导航技术是指卫星在不依赖地面***支持的情况下,仅依靠星载测量设备在轨实时地确定卫星的位置和速度,也称自主轨道确定。对于卫星***来讲,自主导航有利于降低卫星对地面的依赖程度,提高***生存能力,例如战时,当地面测控站遭到敌方的破坏和干扰时,仍能完成轨道的确定和保持,这对军事卫星来讲具有非常重要的意义。此外,自主导航还可以有效减轻地面测控站的负担,降低地面支持成本,从而降低整个航天计划的研制费用。自主导航是卫星实现自主控制的基本前提和基础,也是构造星座、天基组网的关键技术之一。
美国Microcosm公司研制的MANS自主导航***,利用专用的一体化敏感器根据日、地、月的在轨测量数据实时确定航天器的轨道和三轴姿态,是完全意义上的自主导航***(Anthony J.Autonomous Space NavigationExperiment[C].AIAA 92-1710,1992.和Collins J T and Conger R E.MANS:Autonomous Navigation and Orbit Control for Communication Satellites[C].AIAA 94-1127-CP,1994.以及Hosken R W,Wertz J R.Microcosmautonomous navigation system on-orbit operation[C].Keystone,Colorado:The 18th Annual AAS Guidance and Control Conference,AAS 95-074,1995.)。由于太阳、地球、月亮是自然天体,其运动规律清楚且容易识别,导航资源不受限制。因此,研究该自主导航技术具有重要的工程应用价值。早在1994年3月美国空军发射的“空间试验平台-D”航天器就搭载该***进行了在轨飞行试验(Wertz J R.Implementing autonomous orbit control[C].Breckenridge,Colorado:Proceedings of the Annual AAS Guidance andControl Conference,AAS 96-004,1996.)。北京控制工程研究所对日地月自主导航***导航敏感器的特点及***的导航算法进行了理论分析和研究(李捷,陈义庆.航天器自主导航技术的新进展[J].控制工程,1997,(1):76~81.)。
但是,由于传统的基于日地月测量的自主导航算法存在着导航精度低的问题,大大影响和制约了基于日地月测量的自主导航***的发展和应用(TaiFrank,Noerdlinger Peter D.A Low Cost Autonomous Navigation System[C].AAS 89-001,1989.和Daniele Mortari.Moon-Sun Attitude Sensor[J].Journalof Spacecraft and Rockets,1997,34(3):360~364.以及李季陆,陈义庆,孙承启.一种MANS自主导航***的导航精度分析[J].控制工程,1998,(5):1~6.)。对于这种导航***来说,
对于基于日地月自主导航***来说,制约其导航精度的主要因素来源于地心方向和月心方向的确定。而影响地心方向确定的一个重要因素为地球扁率,影响月心方向确定的因素为月相。而这些因素的影响都直接体现在敏感器的脉冲数据上。因此有必要建立日地月导航的仿真***。但目前的半物理仿真试验验证***还无法模拟地球扁率对日地月一体化敏感器测量的影响;也还无法实现对月亮形状变化的高精度模拟,因此,有必要研究导航***的测量原理,对敏感器的脉冲进行模拟,并对导航方法进行仿真验证,建立基于测量信息的自主导航仿真验证***。
发明内容
本发明的技术解决问题是:克服现有不足,提供了一种日地月自主导航***的测量数据模拟方法,能够用于研究导航***误差的修正方案,提高日地月自主导航***的导航性能。
本发明的技术解决方案是:基于日地月方位信息的自主导航***的测量数据模拟方法,实现如下:
(1)日地月一体化敏感器参考脉冲,参考脉冲tref由下式给出:
tref=BR/ωrot+kTs (1)
其中:BR为对称面相对于测量坐标系X轴的滞后角,ωrot为扫描转速,Ts为扫描周期,k为扫描圈数;
(2)日地月一体化敏感器地球脉冲,步骤如下:
第一步:由轨道动力学模型积分得到飞行器的位置信息r,投影在惯性下为其中上标表示转置,上标I代表该变量投影在惯性系下,并定义姿态模式,给出飞行器的姿态RISE;
第二步:在一个扫描周期内,取任意时刻t,结合姿态信息,计算该时刻扫描视线的单位向量在惯性坐标系中的表示
其中:γi为扫描锥的半锥角,BRi为扫描探头相对于敏感器参考线的滞后角,i=1,2分别对应小椎和大锥;
第三步:定义Fun:=b2-4ac并计算Fun的大小,其中:
上式中:ee为地球扁率,RE为地球参考椭球的赤道半径;
第四步:若Fun=0,则t时刻即为脉冲时刻。否则的话,则取另外一个时刻,重复第一步到第四步,计算对应的Fun;由两个时刻的Fun值,利用二分法找到Fun=0对应的时间t,也即脉冲时刻;
(3)日地月一体化敏感器月球脉冲模拟,步骤如下:
第一步:由当前日期利用星历判别月相;
第二步:由轨道动力学模型积分得到飞行器的位置信息r,投影在惯性下为由星历得到惯性系下的月球位置,并定义姿态模式,给出飞行器的姿态RISE;
第三步:利用坐标变换计算星月矢量在敏感器测量坐标系下的表示,从而可以得到方位角φm和高度角δm;
第四步:计算不考虑月相影响的月球脉冲t1moon和t2moon:
其中:σm=arcsin(tanδm/cosβs),βs为敏感器狭缝相对于扫描球经度线的倾斜角,θ为狭缝底部的球心角度;
第五步:若当前的月相为满月,则第四步得到的脉冲即为敏感器的输出月球脉冲,若为弦月,则需要在第四步输出的脉冲的基础上添加延迟δt作为输出脉冲,其中δt为下式给出:
其中:θx为扫描轴与xm即白道平面内与日月方向垂直的轴的夹角,δtmax为最大延迟量。上弦月的时候δt<0,下弦月的时候δt>0;
(4)日地月一体化敏感器太阳脉冲模拟,步骤如下:
第一步:由轨道动力学模型积分得到飞行器的位置信息r,投影在惯性下为由星历得到惯性系下的月球位置,并定义姿态模式,给出飞行器的姿态RISE;
第二步:利用坐标变换计算星日矢量在敏感器测量坐标系下的表示,从而可以得到太阳的方位角φs和高度角δs;
第三步:计算太阳脉冲t1sun和t2sun:
其中:σs=arcsin(tanδs/cosβs),βs为敏感器狭缝相对于扫描球经度线的倾斜角,θ为狭缝底部的球心角度;
通过上述计算得到日地月一体化敏感器的参考脉冲、地球脉冲、月球脉冲和太阳脉冲,从而完成了日地月一体化敏感器脉冲信息的模拟。
本发明与现有技术相比的优点在于:从发表的文献和专利可以看出,目前对基于日地月信息自主导航方法的数学仿真采用直接模拟月心、地心和太阳方向,没有模拟敏感器的直接输出脉冲,因此无法考虑地球扁率和月相变化的影响。本发明直接模拟了敏感器的输出脉冲数据,且考虑了地球扁率和月相的影响。模拟测量数据能够精确和全面地反映复杂因素(地球扁率、月相)对日地月一体化敏感器测量的影响。
附图说明
图1为导航敏感器的结构;
图2为导航敏感器的安装示意图;
图3为导航敏感器的视场;
图4为本发明的基准脉冲示意图;
图5为地平扫描矢量几何关系图;
图6为本发明的红外地平脉冲模拟流程图;
图7为本发明的月相所对应的隶属度函数;
图8为敏感器敏感月球示意图;
图9为月球光线在敏感器测量坐标系的投影;
图10为本发明的上弦月时扫描轴与月球方位示意图;
图11为本发明的扫描转轴方向与脉冲延迟示意图;
图12为本发明的月球脉冲模拟流程图;
图13为敏感器敏感太阳示意图;
图14为太阳光线在敏感器测量坐标系的投影;
图15为本发明的太阳脉冲模拟流程图;
图16为本发明的半锥角38°的扫描锥扫描到的地球脉冲(不考虑地球扁率);
图17为本发明的半锥角为73°的扫描锥扫到的地球脉冲(不考虑地球扁率);
图18为本发明的太阳脉冲(无数据的地方表明卫星处于阴影区);
图19为本发明的月球脉冲(无数据的地方表明卫星处于阴影区);
图20为本发明的没有扁率的月地夹角的误差角;
图21为本发明的没有扁率的日地夹角的误差角;
图22为本发明的半锥角38°的扫描锥扫描到的地球脉冲(考虑地球扁率);
图23为本发明的半锥角73°的扫描锥扫描到的地球脉冲(考虑地球扁率);
图24为本发明的存在扁率的月地夹角误差角;
图25为本发明的存在扁率的日地夹角误差角;
图26为本发明的上弦月时的脉冲(带延迟);
图27为本发明的上弦月时的月地误差角;
图28为本发明实现流程图。
具体实施方式
首先,对日地月一体化敏感器的测量几何关系进行分析,分析敏感器扫描一圈得到的脉冲个数;测量几何的分析结果可以很直观的得到参考脉冲的模拟方法;根据一体化敏感器测量几何关系,利用地球椭球方程,建立对地球脉冲的模拟,给出敏感器扫进、扫出地球的脉冲时刻;然后建立月相模型,根据轨道、星历及姿态信息,利用敏感器测量几何关系确定敏感器狭缝视场扫到月球脉冲;最后根据太阳方向,利用敏感器测量几何关系确定敏感器狭缝视场扫到太阳脉冲。
一、敏感器测量几何分析
导航敏感器由双圆锥扫描式红外地球敏感器和两个扇形狭缝视场的扫描式日、月敏感器构成。其结构图如图1所示。
圆锥扫描式红外地球敏感器具有单一的光学扫描头部。利用反射镜结构得到两个红外视场,扫描后红外视场的轨迹是两个共轴的圆锥。光学头部扫描一圈,热电检测器最多可以检测到四个地平穿越信号,由信号出现的时刻可以确定地心方向矢量相对于卫星的方位,并可求得卫星到地心的距离。在双圆锥扫描式红外敏感器的基础上增加了两个可见光敏感器。在光学头部的扫描过程中,扇形视场扫过空间与扫描转轴夹角在23°和87°之间的球带区域。利用硅光二极管检测器可以敏感到太阳和月球,根据太阳、月球在扇形视场中出现的时刻可以求得其方向矢量相对于卫星的方位。检测器具有多个光强阈值,可以辨别太阳和月球信号,并且可以剔除地球信号。
自主导航***采用两个导航敏感器,其相对于卫星的安装结构如图2所示。导航敏感器视场关系的几何关系见图3。如图所示,将两个圆锥红外扫描视场分别记为红外扫描锥1和红外扫描锥2,相应的半锥角分别记为γ1,γ2,这里的γ1=38°,γ2=73°。将两个扇形狭缝式可见光敏感器视场分别记为狭缝视场1和狭缝视场2,将狭缝视场1和狭缝视场2的对称面记为M1-M2面。狭缝视场1在敏感器赤道平面上超前于狭缝视场2的焦距记为θ,狭缝视场1和狭缝视场2相对于扫描转轴的倾斜角度记为β,这里θ=4°,β=16°,设两个直径为2.5°的红外视场1和红外视场2在敏感器赤道平面上沿扫描方向相对于M1-M2面滞后的角距分别为BR1,BR2。
在导航敏感器扫描一圈的过程中,当红外扫描1和红外扫描2扫入、扫出地球,狭缝视场1和狭缝视场2扫到太阳、月球,以及狭缝视场1和狭缝视场2的对称面通过与敏感器固连的基准点时,导航敏感器均给出相关的脉冲信号,由此可以得到9个脉冲时刻。
本发明实现过程如图28所示,本发明具体实现如下:
二、参考脉冲的模拟
该脉冲指的是对称面M1-M2通过基准点X轴的脉冲时刻,记为tref,如图4所示,可以由下式计算:
tref=BR/ωrot+kTs (6)
其中:BR为对称面相对于XSE的滞后角,ωrot为扫描转速,Ts为扫描周期,k为扫描圈数。
三、地球敏感器脉冲模拟
如图5所示,为红外扫描探头的单位矢量,对于t时刻,其在第一敏感器测量坐标系下的分量可以表示为:
是地心指向地球等效椭球面的矢量,满足椭球方程:
其中:ee为地球扁率。
令红外扫描探头的矢量与地球等效椭球面交线长度为l,由图5的几何关系可以得到:
其中,为敏感器测量坐标系到惯性坐标系的旋转矩阵,RSEI为惯性坐标系到敏感器测量坐标系的旋转矩阵,在脉冲模拟的时候,姿态是事先给定的。因此是已知的。是从地心指向卫星的矢量在惯性系下的表示,可以由轨道动力学积分得到。
令将(9)代入到椭球方程(8)中有:
方程(10)是关于l的二次方程,对其化简可得:
al2+bl+c=0 (11)
其中:
当红外扫描矢量扫描不到地平圈的时候,(11)没有实根;恰好扫地平圈边缘的时候,(11)有两个相等的实根;扫到地球椭球内部的时候,(11)有两个不等的实根。
由二次方程的根的性质可知,对于扫进、扫出地球的脉冲时刻,(11)式有:
Fun=b2-4ac=0
对于扫入地球的脉冲时刻,Fun以递增的形式穿越0点;对于扫出地球的脉冲时刻,Fun以递减的形式穿越0点。
地球红外脉冲模拟流程图见图6:
四、月球脉冲模拟
(1)月相模型
对于月球的脉冲来说,仅考虑四种月相:新月、上弦月、满月和下弦月,当月相为新月的时候,是观测不到月球的,不会输出月球脉冲;弦月的时候,观测到的月球光心和月球质心是不重合的,因此需要在原理的基础上加一个微小的延迟(上弦月延迟为负,下弦月延迟为正);满月的时候不用加延迟。月相严格上分为八个相,令p为月球照亮部分的面积与总面积的百分比,这八个月相分别为:新月(new moon p=0)、峨眉月(waxing crescent moon 0<p<0.5)、上弦月(first quarter half moon p=0.5)、渐盈凸月(waxing gibbousmoon 0.5<p<1)、满月(full moon p=1)、渐亏凸月(waning gibbousmoon 0.5<p<1)、下弦月(last quarter half moon p=0.5)、残月(waning cresentmoon 0<p<0.5)。由于月球脉冲的模拟仅考虑4个相,这就需要利用模糊数学的原理将其他月相模糊近似为这四种月相。将月相看作是逻辑变量,它的逻辑值有:新月、弦月和满月。将p作为月相的测度,p是连续的,需要对其进行模糊化(fuzzification)。如图7所示,新月的隶属度函数(Z型)μN(p)为:
弦月的隶属度函数(∏型)μQ(p)为:
满月的隶属度函数(Z型)μF(p)为:
由模糊集合的并运算可知,脉冲模拟时的月相判别如下:
其中pi的含义如图7所示,表示μN(p)、μQ(p)、μF(p)的边界点。
(2)满月期间的月球脉冲模拟
满月时候,月球的光心与质心重合,月球可以近似为光点。如图8所示,Q1Q2和Q3Q4为两个狭缝敏感器用来敏感太阳和月球信号,两个狭缝对称分布,相对于扫描球经度线的倾斜角为β=16°。狭缝底部的球心角度为θ=4°。在t时刻,月球方向矢量在测量坐标系下的坐标可以通过星历和坐标变换求得,从而两个方位角δm,φm是已知的(见图9)。当GQ12与GS1重合的时候,该时刻就是狭缝视场1扫描到月球的脉冲时刻t1moon;当GQ34与GS1重合的时候,该时刻就是狭缝视场2扫描到太阳的脉冲时刻t2moon。
操作步骤如下:
i).计算t时刻的地月矢量,和地星矢量在惯性系下的表示;
ii).利用坐标变换计算星月矢量在敏感器测量坐标系下的表示,从而可以到两个方位角δm,φm;
iii).计算σm=arcsin(tanδm/cosβs);
iv).由以下两个方程得到脉冲时刻:
v).输出脉冲时刻t1moon,t2moon。
(3)上下弦期间的月球脉冲模拟
这种情况与满月情形下的脉冲模拟类似,但与满月又有所区别。上下弦月的情形下,月球的质心与扫描的光心不重合。模拟机制产生脉冲的时候,应在标准脉冲的基础上加上延迟。下面详细说明延迟的计算方法。
如图10所示,。当扫描轴ZSE与Oxm平行的时候,敏感器由上至下对半月进行扫描,光心与质心是重合的,脉冲没有延迟。当ZSE在o-ymzm平面的时候,敏感器从左至右对半月进行扫描,光心pm与质心om不重合,需要在标准脉冲上加一个延迟δt(上弦月的时候δt<0,下弦月的时候δt>0)。当ZSE平行于ozm或oym得时候,此时的δt最大记为δtmax,下面给出δtmax的近似计算方法。如图10所示,ompm≈Rm/2,Rm为月球半径,oom=384400km,从而:
所以:
由于扫描轴方向对月球上下弦的影响是存在着很复杂的关系,很难用精确的数学模型来描述。只能根据扫描几何图形,在一定程度上对其进行逼近真实的扫描机制。设θx为扫描轴与xm的夹角,根据扫描图形,将脉冲延迟模型近似如下(见图11):zse与xm(或-xm)重合的时候,δt=0(θx=0(π));zse在ymzm平面的时候,δt=δtmax(θx=π/2);脉冲延迟与θx呈线性关系。
因此脉冲的延迟可以用一个简单的线性模型来表示:
(上弦月的时候δt<0,下弦月的时候δt>0)
月球脉冲模拟流程图见图12。
五、太阳脉冲模拟
如图13所示,在t时刻,太阳方向矢量在测量坐标系下的坐标可以通过星历和坐标变换求得,从而两个方位角δs,φs是已知的(见图14)。当GQ12与GS1重合的时候,该时刻就是狭缝视场1扫描到太阳的脉冲时刻t1sun;当GQ34与GS1重合的时候,该时刻就是狭缝视场2扫描到太阳的脉冲时刻t2sun。
操作步骤如下:
i).计算t时刻的地日矢量,和地星矢量在惯性系下的表示;
ii).利用坐标变换计算星日矢量在敏感器测量坐标系下的表示,从而可以得到两个方位角δs,φs;
iii).计算σs=arcsincot(δs)/cosβ
iv).由以下两个方程得到脉冲时刻:
v).输出脉冲时刻t1sun,t2sun
太阳脉冲模拟流程图见图15。
六、仿真及误差评估
日地月一体化敏感器的相关物理参数说明如下:日地月一体化敏感器扫描转速为ωrot=240r/min。圆锥地平敏感器红外地平视场圆锥的两个半锥角分别为γ1=38°和γ2=73°。扫描探头相对于M1M2的滞后角度为BR=0°,日地月一体化敏感器测量坐标系相对于体系的安装角度为γI=π/6,βI=π/6,αI=π/6。扇形日月敏感器相对于与扫描转轴的倾斜角为βs=16°,狭缝视场1在日地月一体化敏感器赤道平面内超前于狭缝视场2的角距θ=4°。
真实轨道的初始六要素为:长半轴为6878km,轨道偏心率为0,轨道倾角为92°,升交点赤经为π/6,近地点幅角为0°,真近角为0°。
(1)脉冲精度评估
轨道的初始时刻取为2012年5月6日0时0分0秒。这时候的月相为满月,因此月球脉冲没有延迟。如果不考虑地球扁率,得到的地球脉冲见图16和图17。太阳脉冲和月球脉冲分别见图18和图19。
为了验证脉冲产生的精度,由于地球方位确定时没有考虑地球扁率,所以将不考虑地球扁率的地球脉冲和日、月脉冲计算日地夹角、月地夹角与实际的日地夹角、月地夹角相比,得到的误差角见图20和图21。从图中可以看出来,误差角的精度达到了10-10度,从而验证了脉冲的精度。
(2)扁率对地心方位确定的影响
轨道的初始时刻取为2012年5月6日0时0分0秒。如果考虑地球扁率,得到的地球脉冲见图22和图23,太阳脉冲和月球脉冲还是如图18和图19。
如果由脉冲计算月星地、日星地方位角的时候不进行扁率修正,则其与真实方位角之间的误差结果见图24和图25。由图中可以看出,如果不进行扁率修正,月地夹角的最大误差为0.13日地夹角的最大误差为0.12°。
(3)月相对月心方位确定的影响
上弦月和下弦月对脉冲的延迟的原理是相似的,只是存在正负的关系。这里仅给出上弦月时的仿真结果。仿真的初始时刻取为2012年5月28日0时0分0秒。月球脉冲如图26所示。脉冲得到的月地夹角与真实的月地夹角的误差角曲线见图27,由图中可知月地夹角的最大误差角为0.12度,且计算的月地夹角要大于真实的月地夹角。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (1)
1.基于日地月方位信息的自主导航***的测量数据模拟方法,其特征在于实现如下:
(1)日地月一体化敏感器参考脉冲,参考脉冲tref由下式给出:
tref=BR/ωrot+kTs (1)
其中:BR为对称面相对于测量坐标系X轴的滞后角,ωrot为扫描转速,Ts为扫描周期,k为扫描圈数;
(2)日地月一体化敏感器地球脉冲,步骤如下:
第一步:由轨道动力学模型积分得到飞行器的位置信息r,投影在惯性下为其中上标表示转置,上标I代表该变量投影在惯性系下,并定义姿态模式,给出飞行器的姿态RISE;
第二步:在一个扫描周期内,取任意时刻t,结合姿态信息,计算该时刻扫描视线的单位向量在惯性坐标系中的表示
其中:γi为扫描锥的半锥角,BRi为扫描探头相对于敏感器参考线的滞后角,i=1,2分别对应小椎和大锥;
第三步:定义Fun=b2-4ac并计算Fun的大小,其中:
上式中:ee为地球扁率,RE为地球参考椭球的赤道半径;
第四步:若Fun=0,则t时刻即为脉冲时刻;否则的话,则取另外一个时刻,重复第一步到第四步,计算对应的Fun;由两个时刻的Fun值,利用二分法找到Fun=0对应的时间t,也即脉冲时刻;
(3)日地月一体化敏感器月球脉冲模拟,步骤如下:
第一步:由当前日期利用星历判别月相;
第二步:由轨道动力学模型积分得到飞行器的位置信息r,投影在惯性下为由星历得到惯性系下的月球位置,并定义姿态模式,给出飞行器的姿态RISE;
第三步:利用坐标变换计算星月矢量在敏感器测量坐标系下的表示,从而可以得到方位角φm和高度角δm;
第四步:计算不考虑月相影响的月球脉冲t1moon和t2moon:
其中:σm=arcsin(tanδm/cosβs),βs为敏感器狭缝相对于扫描球经度线的倾斜角,为狭缝底部的球心角度;
第五步:若当前的月相为满月,则第四步得到的脉冲即为敏感器的输出月球脉冲,若为弦月,则需要在第四步输出的脉冲的基础上添加延迟δt作为输出脉冲,其中δt为下式给出:
其中:θx为扫描轴与xm即白道平面内与日月方向垂直的轴的夹角,δtmax为最大延迟量;上弦月的时候δt<0,下弦月的时候δt>0;
(4)日地月一体化敏感器太阳脉冲模拟,步骤如下:
第一步:由轨道动力学模型积分得到飞行器的位置信息r,投影在惯性下为由星历得到惯性系下的月球位置,并定义姿态模式,给出飞行器的姿态RISE;
第二步:利用坐标变换计算星日矢量在敏感器测量坐标系下的表示,从而可以得到太阳的方位角φs和高度角δs;
第三步:计算太阳脉冲t1sun和t2sun:
其中:σs=arcsin(tanδs/cosβs),βs为敏感器狭缝相对于扫描球经度线的倾斜角,为狭缝底部的球心角度;
通过上述计算得到日地月一体化敏感器的参考脉冲、地球脉冲、月球脉冲和太阳脉冲,从而完成了日地月一体化敏感器脉冲信息的模拟。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110441935.XA CN102564454B (zh) | 2011-12-23 | 2011-12-23 | 基于日地月方位信息的自主导航***的测量数据模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110441935.XA CN102564454B (zh) | 2011-12-23 | 2011-12-23 | 基于日地月方位信息的自主导航***的测量数据模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102564454A CN102564454A (zh) | 2012-07-11 |
CN102564454B true CN102564454B (zh) | 2015-02-11 |
Family
ID=46410501
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110441935.XA Active CN102564454B (zh) | 2011-12-23 | 2011-12-23 | 基于日地月方位信息的自主导航***的测量数据模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102564454B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112179334B (zh) * | 2020-09-15 | 2023-03-14 | 中国科学院微小卫星创新研究院 | 基于两步Kalman滤波的星光导航方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4688092A (en) * | 1986-05-06 | 1987-08-18 | Ford Aerospace & Communications Corporation | Satellite camera image navigation |
US5109346A (en) * | 1990-02-01 | 1992-04-28 | Microcosm, Inc. | Autonomous spacecraft navigation system |
CN101178312A (zh) * | 2007-12-12 | 2008-05-14 | 南京航空航天大学 | 基于多信息融合的航天器组合导航方法 |
CN102175259A (zh) * | 2010-12-31 | 2011-09-07 | 北京控制工程研究所 | 基于地日月一体化敏感器的自主导航仿真试验*** |
-
2011
- 2011-12-23 CN CN201110441935.XA patent/CN102564454B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4688092A (en) * | 1986-05-06 | 1987-08-18 | Ford Aerospace & Communications Corporation | Satellite camera image navigation |
US5109346A (en) * | 1990-02-01 | 1992-04-28 | Microcosm, Inc. | Autonomous spacecraft navigation system |
CN101178312A (zh) * | 2007-12-12 | 2008-05-14 | 南京航空航天大学 | 基于多信息融合的航天器组合导航方法 |
CN102175259A (zh) * | 2010-12-31 | 2011-09-07 | 北京控制工程研究所 | 基于地日月一体化敏感器的自主导航仿真试验*** |
Non-Patent Citations (2)
Title |
---|
基于"日-地-月"信息的卫星自主导航技术研究;黄翔宇等;《哈尔滨工业大学学报》;20021031;第34卷(第5期);第643-646,683页 * |
基于日地月方位信息的自主导航***的脉冲仿真和方位确定;荆武兴等;《航天控制》;20110630;第29卷(第3期);第69-74,79页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102564454A (zh) | 2012-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nimmo et al. | Mean radius and shape of Pluto and Charon from New Horizons images | |
CN104573251B (zh) | 一种星载光学遥感器全视场表观光谱辐亮度确定方法 | |
CN101344391B (zh) | 基于全功能太阳罗盘的月球车位姿自主确定方法 | |
Konopliv et al. | The Vesta gravity field, spin pole and rotation period, landmark positions, and ephemeris from the Dawn tracking and optical data | |
Winn et al. | Measurement of spin-orbit alignment in an extrasolar planetary system | |
CN102162731B (zh) | 基于日地月集成敏感器脉冲数据的卫星高精度自主导航方法 | |
Li et al. | The global image of the Moon obtained by the Chang’E-1: Data processing and lunar cartography | |
CN103900576B (zh) | 一种深空探测自主导航的信息融合方法 | |
CN106871932A (zh) | 基于金字塔搜索地形匹配的星载激光在轨指向检校方法 | |
CN103605887B (zh) | 一种在轨卫星成像波浪水面太阳耀光时空分布仿真模型 | |
CN106124170A (zh) | 一种基于高精度姿态信息的相机光轴指向计算方法 | |
CN102426025B (zh) | 遥感卫星姿态机动时偏流修正角的仿真分析方法 | |
CN102175241A (zh) | 一种火星探测器巡航段自主天文导航方法 | |
CN102288201B (zh) | 用于星敏感器的精度测量方法 | |
CN106840212A (zh) | 基于地面激光光斑质心位置的星载激光在轨几何检校方法 | |
CN102538819A (zh) | 基于双圆锥红外和星敏感器的自主导航半物理仿真试验*** | |
CN105004354A (zh) | 大斜视角下无人机可见光和红外图像目标定位方法 | |
CN106643741A (zh) | 一种卫星相对小行星视觉自主导航方法 | |
CN102288200B (zh) | 用于星敏感器的精度测量*** | |
CN105444778A (zh) | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 | |
CN105446346A (zh) | 遥感卫星对月相对定标姿态调整方法 | |
CN106679653A (zh) | 一种基于星敏感器和星间链路的heo卫星群相对测量方法 | |
CN111505608B (zh) | 一种基于星载激光单片足印影像的激光指向在轨检校方法 | |
CN105486315A (zh) | 遥感卫星对月绝对定标姿态调整方法 | |
RU2304549C2 (ru) | Автономная бортовая система управления космического аппарата гасад-2а |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |