CN108519109B - 一种空间非合作目标相对导航在轨演示安装参数确定方法 - Google Patents
一种空间非合作目标相对导航在轨演示安装参数确定方法 Download PDFInfo
- Publication number
- CN108519109B CN108519109B CN201810359564.2A CN201810359564A CN108519109B CN 108519109 B CN108519109 B CN 108519109B CN 201810359564 A CN201810359564 A CN 201810359564A CN 108519109 B CN108519109 B CN 108519109B
- Authority
- CN
- China
- Prior art keywords
- satellite
- observation
- orbit
- installation
- target
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种空间非合作目标相对导航在轨演示安装参数确定方法,用于在轨演示验证中以星载目标为观测对象,根据观测相机参数和对非合作目标的感知测量与相对导航的观测要求,设计星载目标的安装方位、分离脉冲以及观测相机安装方位,形成可以满足在轨演示验证任务要求的星载目标相对于观测卫星的相对运行轨道。该方法不仅得到观测时间的最优解,又能保证相对导航任务的可观测度、观测卫星和星载目标的安全性。
Description
技术领域
本发明涉及一种空间非合作目标相对导航在轨演示安装参数确定方法,属于卫星轨道设计领域。
背景技术
空间非合作目标,其特征就是无特定标识、无应答装置、无法执行动作配合,而且往往处于失效滚动状态。因此,确定这类非合作目标的位置和状态是对其开展在轨服务的前提。非合作目标感知、测量与辨识在轨在轨演示验证试验中,在轨感知辨识***主要由立体视觉***、智能感知与辨识处理单元、星载目标与目标施放机构等子***组成。观测卫星释放星载目标后,对星载目标进行感知,重构星载目标外形轮廓,辨识位置、速度、角速度等参数。
非合作目标感知测量与相对导航的在轨演示验证任务中,在观测相机的视宽、测距距离等指标条件下,构建合理的相对运动轨道,既为测距、测速、测角等测量任务提供足够的观测时间和良好的观测度,又满足观测卫星与星载目标不相撞的安全性要求,具有重要的工程应用价值。
目前基于合作目标的相对导航技术比较成熟,由于合作目标的相对运动轨道已知,在合作目标存在机动能力的条件下还可以相互配合运动,观测条件好,观测时间充裕,这些特点均与非合作目标的相对导航不同。关于非合作的目标的相对导航,尚无在轨试验关键参数的确定方法。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种空间非合作目标相对导航在轨演示安装参数确定方法,在观测相机视场、观测距离、观测最短时间、相对导航可观测度等约束条件下,要求既满足验证任务的功能性要求,又满足安全性要求。
本发明的技术解决方案是:一种空间非合作目标相对导航在轨演示安装参数确定方法,所述安装分离参数确定方法包括星载目标的安装方向偏置角ΔαT、星载目标的分离初始速度Δv、观测相机安装偏角Δβ,所述方法包括如下步骤:
(1)、以观测卫星质心为坐标原点O,地心指向观测卫星质心方向为X轴正方向,观测卫星飞行轨道的法向为Z轴正方向,按照与X轴、Z轴满足右手定则的原则确定Y轴,建立观测卫星质心轨道坐标系;
(2)、确定星载目标的初始安装位置,使其分离速度方向在观测卫星质心轨道坐标系下的YOZ平面内,定义星载目标的安装方向偏置ΔαT为星载目标分离速度方向与观测卫星飞行轨道的法向之间的夹角;确定观测相机的初始安装位置,使其视轴方向在XOZ平面内,定义观测相机安装偏角Δβ为观测相机视轴方向与观测卫星飞行轨道的法向之间的夹角;
(3)、建立分离后的星载目标在观测卫星质心轨道坐标系下的相对运动模型;
(4)、根据星载目标在观测卫星质心轨道坐标系下的相对运动模型,模拟分离过程,得到最长观测时间Jmax与安装分离参数的对应关系;
(5)、以最长观测时间Jmax为优化目标函数,以所述安装分离参数为设计参数,下述条件为约束函数,根据最长观测时间Jmax与安装分离参数的对应关系,采用遗传算法计算最长观测时间Jmax对应的安装分离参数;所述约束函数为:
Δv≥Δvmin
0≤ΔαT≤θ1/2
其中,Δvmin为观测卫星分离机构所能提供的最小分离速度,θ1为观测相机长边方向视场宽度角,θ2为观测相机短边方向视场宽度角,θE为地球对观测相机的半张角。
所述星载目标在观测卫星质心轨道坐标系下的相对运动模型为:
所述步骤(2)的具体实现为:
(4.1)、改变星载目标的安装位置,得到星载目标在观测卫星质心轨道坐标系下的初始位置(x0,y0,z0);
(4.2)、改变星载目标的安装方向偏置ΔαT、星载目标的分离初始速度Δv,计算星载目标在观测卫星质心轨道坐标系下的初始速度:
(4.3)、根据分离后的星载目标在观测卫星质心轨道坐标系与观测卫星的相对运动模型,实时计算星载目标相对于观测相机的飞行轨迹,并将其转换为球坐标,求得星载目标的方位角α、俯仰角β,具体的转换公式如下:
式中,ρ为星载目标与观测卫星的距离;
(4.4)、改变安装分离参数,重复步骤(4.1)~步骤(4.4),判断星载目标的方位角α、俯仰角β大小,将或或ρ>dmax时所对应的时刻确定为观测相机对最长观测时间Jmax,从而得到最长观测时间Jmax所对应的安装分离参数的对应关系,其中,dmax为观测相机的最大可观测距离。
所述步骤(4)中星载目标的分离初始速度Δv的初值确定为:
Δv=4dmax/T
式中,T为观测卫星飞行轨道周期。
所述步骤(4)中星载目标的安装方向偏置ΔαT确定为:
ΔαT=θ1/4。
所述步骤(4)中观测相机安装偏角Δβ的初值确定为0。
所述步骤(3)采用C-W方程建立分离后的星载目标在观测卫星质心轨道坐标系下的相对运动模型。
本发明与现有技术相比的优点在于:
(1)、本发明明确了非合作目标相对导航的安装分离关键参数,建立了观测时间与安装分离关键参数的关系,给出了相关约束条件;
(2)、本发明相对运行模型采用C-W方程解析进行计算,对于10米内的相对运动计算精度既满足任务要求,又满足优化计算的计算速度要求;
(3)、本发明给出了观测相机对星载目标观测时间最优的计算过程,该计算方法流程初值确定简单,约束条件清晰,有利于最优观测时间的求取;
(4)、本发明星载目标和观测相机初始安装位置均在卫星飞行轨道的法向,在此基础上对星载目标和观测相机的安装方位进行优化设计,不仅得到观测时间的最优解,满又能保证相对导航任务的可观测度、观测卫星和星载目标的安全性。
附图说明
图1为本发明空间非合作目标相对导航在轨演示安装分离参数确定方法流程图;
图2为观测卫星质心轨道坐标系和及其相应的球坐标系的定义。
具体实施方式
以下结合附图和具体实施例对本发明进行详细说明。
如图1所示,本发明提供的一种空间非合作目标相对导航在轨演示安装参数确定方法,所述安装分离参数确定方法包括星载目标的安装方向偏置ΔαT、星载目标的分离初始速度Δv、观测相机安装偏角Δβ,该方法包括如下步骤:
(1)、以观测卫星质心为坐标原点O,地心指向观测卫星质心方向为X轴正方向,观测卫星飞行轨道的法向为Z轴正方向,按照与X轴、Z轴满足右手定则的原则确定Y轴,建立观测卫星质心轨道坐标系;
(2)、确定星载目标的初始安装位置,使其分离速度方向在观测卫星质心轨道坐标系下的YOZ平面内,定义星载目标的安装方向偏置ΔαT为星载目标分离速度方向与观测卫星飞行轨道的法向之间的夹角;确定观测相机的初始安装位置,使其视轴方向在XOZ平面内,定义观测相机安装偏角Δβ为观测相机视轴方向与观测卫星飞行轨道的法向之间的夹角;
观测相机的初始安装位置,观测相机视场长边与观测卫星质心轨道坐标系Y轴平行,初始观测安装偏角Δβ为0。在此基础上对星载目标和观测相机的安装方位进行优化设计,不仅得到观测时间的最优解,满又能保证相对导航任务的可观测度、观测卫星和星载目标的安全性。
(3)、建立分离后的星载目标在观测卫星质心轨道坐标系下的相对运动模型;
本发明的具体实施例采用C-W方程建立分离后的星载目标在观测卫星质心轨道坐标系下的相对运动模型。参考轨道为圆时相对运动,根据C-W方程,星载目标在无外力作用下相对于观测卫星的相对运动解析解在观测卫星质心轨道坐标系中表达式,即所述星载目标在观测卫星质心轨道坐标系下的相对运动模型为:
其中,t为星载目标分离后的飞行时间,x0、y0和z0为星载目标在观测卫星质心轨道坐标系的位置,即星载目标与观测卫星相对位置,和为星载目标在观测卫星质心轨道坐标系下的运动速度,即星载目标与观测卫星相对运动速度;n表示观测卫星的轨道角速度,利用观测卫星轨道半长轴a0和地球引力常数μ,n可以写为:
C-W方程准确度比较高,10米内的相对运动计算精度既满足任务要求,又满足迭代计算的计算速度要求。
(4)、根据星载目标在观测卫星质心轨道坐标系下的相对运动模型,模拟分离过程,得到最长观测时间Jmax与安装分离参数的对应关系;具体实现为:
(4.1)、改变星载目标的安装位置,得到星载目标在观测卫星质心轨道坐标系下的初始位置(x0,y0,z0);
(4.2)、改变星载目标的安装方向偏置ΔαT、星载目标的分离初始速度Δv,计算星载目标在观测卫星质心轨道坐标系下的初始速度:
根据观测卫星的参考轨道的半长轴a0,可得观测卫星飞行轨道周期为T
由于观测相机观测到星载目标的时间最长为观测卫星飞行轨道周期T的1/2,折中考虑取观测相机观测到星载目标的时间为观测卫星飞行轨道周期T的1/4,根据观测相机的最大可观测距离dmax,计算得到星载目标的分离初始速度为:
Δv=4dmax/T
式中,dmax为观测相机的最大可观测距离,T为观测卫星飞行轨道周期。
为了让分离后的星载目标位于观测相机视场内,所述星载目标的安装方向偏置ΔαT确定为:
ΔαT=θ1/4。
(4.3)、根据分离后的星载目标在观测卫星质心轨道坐标系与观测卫星的相对运动模型,实时计算星载目标相对于观测相机的飞行轨迹,并将其转换为球坐标,求得星载目标的方位角α、俯仰角β,具体的转换公式如下:
式中,ρ为星载目标与观测卫星的距离;
观测卫星质心轨道坐标系相应的球坐标定义如图2所示。
(4.4)、改变安装分离参数,重复步骤(4.1)~步骤(4.4),判断星载目标的方位角α、俯仰角β大小,将或或ρ>dmax时所对应的时刻确定为观测相机对最长观测时间Jmax,从而得到最长观测时间Jmax所对应的安装分离参数的对应关系,其中,dmax为观测相机的最大可观测距离。
(5)、以最长观测时间Jmax为优化目标函数,以所述安装分离参数为设计参数,下述条件为约束函数,根据最长观测时间Jmax与安装分离参数的对应关系,采用遗传算法计算最长观测时间Jmax对应的安装分离参数;所述约束函数为:
Δv≥Δvmin
0≤ΔαT≤θ1/2
其中,Δvmin为观测卫星分离机构所能提供的最小分离速度,θ1为观测相机长边方向视场宽度角,θ2为观测相机短边方向视场宽度角,θE为地球对观测相机的半张角:
θE=arcsin(a0/RE);
式中,RE为地球半径。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (7)
1.一种空间非合作目标相对导航在轨演示安装分离参数确定方法,所述安装分离参数包括星载目标的安装方向偏置角ΔαT、星载目标的分离初始速度Δv、观测相机安装偏角Δβ,其步骤在于:
(1)、以观测卫星质心为坐标原点O,地心指向观测卫星质心方向为X轴正方向,观测卫星飞行轨道的法向为Z轴正方向,按照与X轴、Z轴满足右手定则的原则确定Y轴,建立观测卫星质心轨道坐标系;
(2)、确定星载目标的初始安装位置,使其分离速度方向在观测卫星质心轨道坐标系下的YOZ平面内,定义星载目标的安装方向偏置角ΔαT为星载目标分离速度方向与观测卫星飞行轨道的法向之间的夹角;确定观测相机的初始安装位置,使其视轴方向在XOZ平面内,定义观测相机安装偏角Δβ为观测相机视轴方向与观测卫星飞行轨道的法向之间的夹角;
(3)、建立分离后的星载目标在观测卫星质心轨道坐标系下的相对运动模型;
(4)、根据星载目标在观测卫星质心轨道坐标系下的相对运动模型,模拟分离过程,得到最长观测时间Jmax与安装分离参数的对应关系;
(5)、以最长观测时间Jmax为优化目标函数,以所述安装分离参数为设计参数,下述条件为约束函数,根据最长观测时间Jmax与安装分离参数的对应关系,采用遗传算法计算最长观测时间Jmax对应的安装分离参数;所述约束函数为:
Δv≥Δvmin
0≤ΔαT≤θ1/2
其中,Δvmin为观测卫星分离机构所能提供的最小分离速度,θ1为观测相机长边方向视场宽度角,θ2为观测相机短边方向视场宽度角,θE为地球对观测相机的半张角。
3.根据权利要求2所述的一种空间非合作目标相对导航在轨演示安装分离参数确定方法,其特征在于所述步骤(4)的具体实现为:
(4.1)、改变星载目标的安装位置,得到星载目标在观测卫星质心轨道坐标系下的初始位置(x0,y0,z0);
(4.2)、改变星载目标的安装方向偏置角ΔαT、星载目标的分离初始速度Δv,计算星载目标在观测卫星质心轨道坐标系下的初始速度:
(4.3)、根据分离后的星载目标在观测卫星质心轨道坐标系与观测卫星的相对运动模型,实时计算星载目标相对于观测相机的飞行轨迹,并将其转换为球坐标,求得星载目标的方位角α、俯仰角β,具体的转换公式如下:
式中,ρ为星载目标与观测卫星的距离;
4.根据权利要求3所述的一种空间非合作目标相对导航在轨演示安装分离参数确定方法,其特征在于所述步骤(4.2)中星载目标的分离初始速度Δv的初值确定为:
Δv=4dmax/T
式中,T为观测卫星飞行轨道周期。
5.根据权利要求3所述的一种空间非合作目标相对导航在轨演示安装分离参数确定方法,其特征在于所述步骤(4.2)中星载目标的安装方向偏置角ΔαT确定为:
ΔαT=θ1/4。
6.根据权利要求3所述的一种空间非合作目标相对导航在轨演示安装分离参数确定方法,其特征在于所述步骤(4.4)中观测相机安装偏角Δβ的初值确定为0。
7.根据权利要求1所述的一种空间非合作目标相对导航在轨演示安装分离参数确定方法,其特征在于所述步骤(3)采用C-W方程建立分离后的星载目标在观测卫星质心轨道坐标系下的相对运动模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810359564.2A CN108519109B (zh) | 2018-04-20 | 2018-04-20 | 一种空间非合作目标相对导航在轨演示安装参数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810359564.2A CN108519109B (zh) | 2018-04-20 | 2018-04-20 | 一种空间非合作目标相对导航在轨演示安装参数确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108519109A CN108519109A (zh) | 2018-09-11 |
CN108519109B true CN108519109B (zh) | 2020-04-10 |
Family
ID=63429887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810359564.2A Active CN108519109B (zh) | 2018-04-20 | 2018-04-20 | 一种空间非合作目标相对导航在轨演示安装参数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108519109B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118246252A (zh) * | 2024-05-28 | 2024-06-25 | 南京航空航天大学 | 一种非合作目标卫星的脉冲机动检测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101423121A (zh) * | 2008-12-08 | 2009-05-06 | 北京航空航天大学 | 应用双视线测量的非合作目标航天器相对导航制导方法 |
CN103438888A (zh) * | 2013-07-24 | 2013-12-11 | 西北工业大学 | 一种对空间非合作目标自主交会的相对导航方法 |
CN104765373A (zh) * | 2015-01-30 | 2015-07-08 | 北京控制工程研究所 | 一种星上相对运动状态获取方法 |
CN107529375B (zh) * | 2011-04-20 | 2015-07-15 | 上海航天控制工程研究所 | 一种用星上视线角度信息的相对导航方法 |
CN107621266A (zh) * | 2017-08-14 | 2018-01-23 | 上海宇航***工程研究所 | 基于特征点跟踪的空间非合作目标相对导航方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10571224B2 (en) * | 2015-05-04 | 2020-02-25 | Propagation Research Associates, Inc. | Systems, methods and computer-readable media for improving platform guidance or navigation using uniquely coded signals |
-
2018
- 2018-04-20 CN CN201810359564.2A patent/CN108519109B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101423121A (zh) * | 2008-12-08 | 2009-05-06 | 北京航空航天大学 | 应用双视线测量的非合作目标航天器相对导航制导方法 |
CN107529375B (zh) * | 2011-04-20 | 2015-07-15 | 上海航天控制工程研究所 | 一种用星上视线角度信息的相对导航方法 |
CN103438888A (zh) * | 2013-07-24 | 2013-12-11 | 西北工业大学 | 一种对空间非合作目标自主交会的相对导航方法 |
CN104765373A (zh) * | 2015-01-30 | 2015-07-08 | 北京控制工程研究所 | 一种星上相对运动状态获取方法 |
CN107621266A (zh) * | 2017-08-14 | 2018-01-23 | 上海宇航***工程研究所 | 基于特征点跟踪的空间非合作目标相对导航方法 |
Non-Patent Citations (4)
Title |
---|
"Relative navigation to non-cooperative targets in LEO: achievable accuracy from radar tracking measurements";Kahle, R. 等,;《International Journal of Space Science and Engineering》;20141231;第2卷(第1期);1-13页 * |
"基于双目光学图像信息的空间非合作目标自主相对导航方法";王大轶 等,;《中国科学: 物理学 力学 天文学》;20190103;第49卷(第2期);1-8页 * |
"空间非合作目标相对导航技术研究";李文跃 等,;《航空兵器》;20120630(第3期);14-17,50页 * |
"非合作航天器自主相对导航研究综述";王大轶 等,;《控制理论与应用》;20181031;第35卷(第10期);1392-1404页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108519109A (zh) | 2018-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108519110B (zh) | 基于图像的空间非合作目标自主导航在轨验证*** | |
Accardo et al. | Flight test of a radar-based tracking system for UAS sense and avoid | |
CN104316060B (zh) | 空间非合作目标的交会对接方法与装置 | |
US8970401B2 (en) | Using image sensor and tracking filter time-to-go to avoid mid-air collisions | |
CN106595674B (zh) | 基于星敏感器和星间链路的heo卫星编队飞行自主导航方法 | |
CN103913181B (zh) | 一种基于参数辨识的机载分布式pos传递对准方法 | |
CN105184002B (zh) | 一种数传天线指向角度的仿真分析方法 | |
Ahn et al. | Fast alignment using rotation vector and adaptive Kalman filter | |
US10209159B2 (en) | Displacement detecting device for structural object, sharing system of displacement of structural object, and method and program of detecting displacement of structural object | |
CN106595673B (zh) | 面对地球静止轨道目标操作的空间多机器人自主导航方法 | |
CN111102981B (zh) | 一种基于ukf的高精度卫星相对导航方法 | |
CN101762273A (zh) | 一种深空探测器软着陆自主光学导航方法 | |
CN104537898A (zh) | 一种空地协同的无人机感知规避***及其规避方法 | |
CN110412868B (zh) | 一种使用星间光学图像的非合作航天器轨道确定方法 | |
CN103697883B (zh) | 一种基于天际线成像的飞行器水平姿态确定方法 | |
CN106767768A (zh) | 一种双星编队的自主导航方法 | |
Xu et al. | A pose measurement method of a non-cooperative GEO spacecraft based on stereo vision | |
CN108519109B (zh) | 一种空间非合作目标相对导航在轨演示安装参数确定方法 | |
Albrektsen et al. | Robust and secure UAV navigation using GNSS, phased-array radio system and inertial sensor fusion | |
CN110146092B (zh) | 基于导航信息评价的双体小行星探测轨迹优化方法 | |
EP3166849B1 (en) | System for locating the barycenter of at least one object orbiting in space and related process of physical and mechanical characterization of the identified object | |
Nordlund et al. | Probabilistic noncooperative near mid-air collision avoidance | |
CN109341685A (zh) | 一种基于单应变换的固定翼飞机视觉辅助着陆导航方法 | |
CN106643742B (zh) | 一种卫星自主连续观测小行星的方法 | |
KR100918484B1 (ko) | 별과 지구의 영상 인식을 이용한 인공 위성의 위치 및 궤도결정 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |