CN109813301B - 一种最佳导航星方位快速确定方法 - Google Patents
一种最佳导航星方位快速确定方法 Download PDFInfo
- Publication number
- CN109813301B CN109813301B CN201910083726.9A CN201910083726A CN109813301B CN 109813301 B CN109813301 B CN 109813301B CN 201910083726 A CN201910083726 A CN 201910083726A CN 109813301 B CN109813301 B CN 109813301B
- Authority
- CN
- China
- Prior art keywords
- gradient
- current position
- determining
- star
- optimal navigation
- 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
Landscapes
- Navigation (AREA)
Abstract
本发明公开了一种最佳导航星方位快速确定方法,首先在给出初始点、步长、误差的基础上,计算出当前位置损失函数的梯度,然后利用步长和损失函数的梯度计算出当前位置下降的距离,对位置进行更新并判断更新位置与当前位置差值的绝对值是否小于误差值,若小于则更新的位置即为最佳导航星方位,否则继续进行迭代。本发明利用梯度下降法确定最佳导航星,迭代次数少,减小了弹上计算量,加快了运算速度。
Description
技术领域
本发明涉及制导控制技术领域,可应用于远程弹道导弹、运载火箭的主动段星光/惯性复合制导,尤其适用于复合制导理论最佳星的确定。
背景技术
现有的远程弹道导弹、运载火箭的主动段一般采用惯性制导技术,即通过加速度计、陀螺等惯性器件测量载体的视加速度、角速度等信息,经导航积分可获得载体的实时运行状态,并传递给制导系统。制导系统计算得到指导指令,控制载体的运动状态达到期望值。在惯性制导过程中,不可避免地存在定位定向误差、初始对准误差、惯性器件漂移误差等误差因素的影响,使惯性制导的精度降低,为此提出星光/惯性复合制导方法。
星光/惯性复合制导,是一种星光制导和惯性制导相结合的复合制导方式,更准确地说,是在纯惯性制导的基础上辅以星光制导。它利用在弹上测量恒星所获得的准确的惯性空间方位基准,来估计惯性导航平台(物理平台或数学平台)与发射惯性坐标系之间的失准角,并根据失准角估计值来修正惯性导航误差造成的落点偏差,以达到综合利用惯性导航信息和星光信息来提高导弹制导精度的目的。
根据观测矢量个数的不同,星光/惯性制导可以分为单星方案和双星方案,但是单星方案需要测量最佳方位的恒星才能与双星方案达到同样的精度。在单星调平台方案的实际应用中,一方面要保证导弹命中精度,同时提高导弹的机动快速发射能力,尽量缩短确定最佳星的时间;另一方面,单星调平台方案落点精度指标是关于测星方位的复杂非线性函数,难以从解析或半解析的方法入手进行解决,难以保证时间指标的要求。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种最佳导航星方位快速确定方法,在保证精度的前提下,减少迭代次数,降低弹上的计算量,加快计算速度。
为解决上述技术问题,本发明所采用的技术方案是:一种最佳星方位快速确定方法,具体操作过程如下:
Step1:给定初始点X0=[e0 σ0]、步长α、误差ε;
Step2:确定出当前位置Xi=[ei σi](i=0,1,2,…)损失函数(即圆概率偏差)CEP=f(ei,σi)的梯度,对于ei、σi,其梯度的表达式如下:
Step3:判断更新后位置的梯度是否小于当前位置的梯度,若小于,转Step4,若不小于,迭代结束(当i=0时,跳过此步骤)。
Step4:用步长α乘以损失函数的梯度,得到当前位置下降的距离:
Step5:更新位置Xi+1=[ei+1 σi+1],其更新的表达式如下:
并判断下式是否成立:
若成立则停止迭代(若只成立其中一种情况,则该种情况停止迭代,继续迭代另一种情况),否则转Step2。
与现有技术相比,本发明所具有的有益效果为:本发明利用梯度下降法,计算出当前位置损失函数的梯度以及下降的距离,并利用下降的距离更新位置,从而确定出最佳星的方位。其优点在于:在保证精度的前提下,迭代次数少,显著降低了弹上的计算量,加快了计算速度。
附图说明
图1利用梯度下降法进行最佳星方位确定的流程图。
具体实施方式
以某远程弹道导弹星光/惯性复合制导系统为例,说明本发明实施的方式。
假设发射时间为2018年1月1日0时0分0秒,发射点的纬度为10度,经度为10度,发射点高程为0m,发射方位角为45度,星敏感器的安装角为[2010]度。
S1.给定初始点X0=[20 10]°、α=0.0001、ε=0.1°;
S2.计算出当前位置损失函数的梯度:
S3.用步长α乘以损失函数的梯度,得到当前位置下降的距离:
S4.更新ei、σi:
并判断下式是否成立:
下式不成立,转S2。
将利用梯度下降法与单纯形调优法确定出的最佳导航星方位进行比较,得到的结果如下表所示:
表1梯度下降法与单纯形的比较结果
Claims (3)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910083726.9A CN109813301B (zh) | 2019-01-29 | 2019-01-29 | 一种最佳导航星方位快速确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910083726.9A CN109813301B (zh) | 2019-01-29 | 2019-01-29 | 一种最佳导航星方位快速确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109813301A CN109813301A (zh) | 2019-05-28 |
CN109813301B true CN109813301B (zh) | 2021-04-02 |
Family
ID=66605539
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910083726.9A Active CN109813301B (zh) | 2019-01-29 | 2019-01-29 | 一种最佳导航星方位快速确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109813301B (zh) |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6470270B1 (en) * | 2000-11-07 | 2002-10-22 | The Boeing Company | Lost-in-space stellar attitude acquisition using multiple star trackers |
CN101763087B (zh) * | 2009-12-29 | 2011-08-31 | 浙江大学 | 一种基于非线性共轭梯度法的工业过程动态优化系统及方法 |
CN103217157B (zh) * | 2012-01-18 | 2016-02-03 | 北京自动化控制设备研究所 | 一种惯导/里程计自主组合导航方法 |
CN106441280B (zh) * | 2012-09-17 | 2019-05-14 | 常州工学院 | 用于星敏感器的筛选导航星的方法 |
CN102865865B (zh) * | 2012-09-17 | 2015-07-22 | 常州工学院 | 一种采用十字链表的星敏感器筛选导航星的方法 |
CN103591956B (zh) * | 2013-11-27 | 2016-02-03 | 北京理工大学 | 一种基于可观测性分析的深空探测器自主导航方法 |
CN106154291B (zh) * | 2016-07-08 | 2019-03-05 | 北京航空航天大学 | 基于卫星几何结构的sbas可用性预测方法与装置 |
CN108594274B (zh) * | 2018-02-01 | 2021-09-03 | 沈阳航空航天大学 | 一种自适应混沌粒子群优化的组合导航选星方法 |
-
2019
- 2019-01-29 CN CN201910083726.9A patent/CN109813301B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109813301A (zh) | 2019-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101793523B (zh) | 一种组合导航和光电探测一体化系统 | |
CN107618678B (zh) | 卫星姿态角度偏差下的姿控信息联合估计方法 | |
EP3213034B1 (en) | Inertial navigation system with compensation of roll scale factor error | |
CN109724624B (zh) | 一种适用于机翼挠曲变形的机载自适应传递对准方法 | |
CN110672128B (zh) | 一种星光/惯性组合导航及误差在线标定方法 | |
CN107478110B (zh) | 一种基于状态观测器的旋转弹姿态角计算方法 | |
CN109269504B (zh) | 一种具有末端约束的姿态机动路径规划方法 | |
US5193064A (en) | Method and apparatus of integrating Global Positioning System and Inertial Navigation System without using accelerometers | |
CN102168978B (zh) | 一种船用惯性导航系统摇摆基座开环对准方法 | |
CN113503892B (zh) | 一种基于里程计和回溯导航的惯导系统动基座初始对准方法 | |
CN111623770B (zh) | 一种基于速度误差开环修正的提高惯性制导精度的方法 | |
CN113551668A (zh) | 一种航天器惯性/恒星星光矢量/星光折射组合导航方法 | |
JP5338464B2 (ja) | 慣性航法装置、飛翔体及び航法データ算出方法 | |
CN109813301B (zh) | 一种最佳导航星方位快速确定方法 | |
CN111637883B (zh) | 基于速度和位置误差闭环反馈的提高惯性制导精度的方法 | |
Bin et al. | Attitude dynamics aiding for line-of-sight angular rate reconstruction of strap-down seeker | |
Zhang et al. | Strapdown stellar-inertial guidance system for launch vehicle | |
CN109813302B (zh) | 最佳可用导航星快速确定方法 | |
CN114111805B (zh) | 运载火箭多源多类测量数据位置基准高精度对齐方法 | |
US20210116247A1 (en) | Method for harmonising two inertial measurement units with one another and navigation system implementing this method | |
CN110765638B (zh) | 扰动引力场对星光/惯性复合制导影响的计算方法 | |
Zhang et al. | Deep data fusion method for missile-borne inertial/celestial system | |
CN113124903B (zh) | 传递对准下基于姿态匹配的最小二乘陀螺零偏快速估计 | |
CN114608579B (zh) | 基于惯性星光组合的星矢量气动光学零位误差校正方法 | |
CN113884704B (zh) | 运载火箭多源多类测量数据速度基准高精度对齐方法 |
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 |