CN109813301A - 一种最佳导航星方位快速确定方法 - Google Patents
一种最佳导航星方位快速确定方法 Download PDFInfo
- Publication number
- CN109813301A CN109813301A CN201910083726.9A CN201910083726A CN109813301A CN 109813301 A CN109813301 A CN 109813301A CN 201910083726 A CN201910083726 A CN 201910083726A CN 109813301 A CN109813301 A CN 109813301A
- Authority
- CN
- China
- Prior art keywords
- gradient
- less
- error
- optimal navigation
- initial point
- 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.)
- Granted
Links
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)
1.一种最佳导航星方位快速确定方法,其特征在于,包括以下步骤:
1)确定初始点X0=[e0 σ0]、步长α、误差ε;e0为初始点的经度;σ0为初始点的纬度;
2)计算出当前位置Xi=[ei σi]关于圆概率偏差CEP的梯度
3)判断更新后位置梯度的绝对值是否小于当前位置梯度的绝对值,若小于,转步骤4),若不小于,迭代结束;
4)用步长α乘以损失函数的梯度,得到当前位置下降的距离;
5)更新位置Xi+1,并判断Xi+1与Xi差值的绝对值是否小于误差值ε,若小于则停止迭代,Xi+1即为最佳导航星方位,否则转步骤2)。
2.根据权利要求1所述的最佳导航星方位快速确定方法,其特征在于,步骤4)中,当前位置下降距离的计算公式为:
3.根据权利要求3所述的最佳导航星方位快速确定方法,其特征在于,步骤5)中,更新后位置Xi+1=[ei+1 σi+1]的计算公式为:
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 true CN109813301A (zh) | 2019-05-28 |
CN109813301B 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) |
Citations (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 |
CN101763087A (zh) * | 2009-12-29 | 2010-06-30 | 浙江大学 | 一种基于非线性共轭梯度法的工业过程动态优化系统及方法 |
CN102840861A (zh) * | 2012-09-17 | 2012-12-26 | 常州工学院 | 一种用于星敏感器的筛选导航星的方法 |
CN102865865A (zh) * | 2012-09-17 | 2013-01-09 | 常州工学院 | 一种采用十字链表的星敏感器筛选导航星的方法 |
CN103217157A (zh) * | 2012-01-18 | 2013-07-24 | 北京自动化控制设备研究所 | 一种惯导/里程计自主组合导航方法 |
CN103591956A (zh) * | 2013-11-27 | 2014-02-19 | 北京理工大学 | 一种基于可观测性分析的深空探测器自主导航方法 |
CN106154291A (zh) * | 2016-07-08 | 2016-11-23 | 北京航空航天大学 | 基于卫星几何结构的sbas可用性预测方法与装置 |
CN108594274A (zh) * | 2018-02-01 | 2018-09-28 | 沈阳航空航天大学 | 一种自适应混沌粒子群优化的组合导航选星方法 |
-
2019
- 2019-01-29 CN CN201910083726.9A patent/CN109813301B/zh active Active
Patent Citations (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 |
CN101763087A (zh) * | 2009-12-29 | 2010-06-30 | 浙江大学 | 一种基于非线性共轭梯度法的工业过程动态优化系统及方法 |
CN103217157A (zh) * | 2012-01-18 | 2013-07-24 | 北京自动化控制设备研究所 | 一种惯导/里程计自主组合导航方法 |
CN102840861A (zh) * | 2012-09-17 | 2012-12-26 | 常州工学院 | 一种用于星敏感器的筛选导航星的方法 |
CN102865865A (zh) * | 2012-09-17 | 2013-01-09 | 常州工学院 | 一种采用十字链表的星敏感器筛选导航星的方法 |
CN103591956A (zh) * | 2013-11-27 | 2014-02-19 | 北京理工大学 | 一种基于可观测性分析的深空探测器自主导航方法 |
CN106154291A (zh) * | 2016-07-08 | 2016-11-23 | 北京航空航天大学 | 基于卫星几何结构的sbas可用性预测方法与装置 |
CN108594274A (zh) * | 2018-02-01 | 2018-09-28 | 沈阳航空航天大学 | 一种自适应混沌粒子群优化的组合导航选星方法 |
Non-Patent Citations (5)
Title |
---|
HUANG JINGSHUAI 等: ""Design and Realization of Terminal Guidance Law for Passive Homing"", 《 THE 27TH CHINESE CONTROL AND DECISION CONFERENCE》 * |
SERKANKIRANYAZ 等: "《多维粒子群优化在机器学习与模式识别中的应用》", 30 September 2017, 国防工业出版社 * |
WU SHILONG 等: ""Joint FDOA and TDOA location algorithm and performance analysis of dual-satellite formations"", 《2010 2ND INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING SYSTEMS》 * |
叶兵 等: ""星光/ 惯性复合制导次佳星快速确定方法研究"", 《科技创新导报》 * |
韦华: ""基于梯度法的目标定位精度分析"", 《舰船电子工程》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109813301B (zh) | 2021-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103662091B (zh) | 一种基于相对导航的高精度安全着陆制导方法 | |
CN107618678B (zh) | 卫星姿态角度偏差下的姿控信息联合估计方法 | |
CN105021092B (zh) | 一种捷联寻的导引头的制导信息提取方法 | |
CN102538821B (zh) | 一种快速、参数分段式捷联惯性导航系统自对准方法 | |
US20100133374A1 (en) | Projectile navigation enhancement method | |
CN105180728B (zh) | 基于前数据的旋转制导炮弹快速空中对准方法 | |
CN103759729B (zh) | 采用捷联惯导的月球软着陆地面试验用初始姿态获取方法 | |
CN109724624A (zh) | 一种适用于机翼挠曲变形的机载自适应传递对准算法 | |
CN105115508A (zh) | 基于后数据的旋转制导炮弹快速空中对准方法 | |
CN100559190C (zh) | 一种对加速度计零位偏差进行在轨标定的方法 | |
CN105486307A (zh) | 针对机动目标的视线角速率估计方法 | |
CN115248038B (zh) | 一种发射系下的sins/bds组合导航工程算法 | |
CN102645223B (zh) | 一种基于比力观测的捷联惯导真空滤波修正方法 | |
CN111504256A (zh) | 一种基于最小二乘法的滚转角实时估计方法 | |
CN110160519A (zh) | 用于脉冲修正火箭弹的弹体姿态计算方法 | |
CN110940357B (zh) | 一种用于旋转惯导单轴自对准的内杆臂标定方法 | |
CN109813301A (zh) | 一种最佳导航星方位快速确定方法 | |
CN104931048A (zh) | 一种基于mimu的肩扛制导火箭弹的导航方法 | |
Ben et al. | Research on error modulating of SINS based on single-axis rotation | |
CN110501027A (zh) | 一种用于双轴旋转mems-sins的最优转停时间分配方法 | |
CN110006446A (zh) | 一种基于棱镜的惯组输出标校方法 | |
US20180216938A1 (en) | Digital sight for hand-carried projectile-firing device and method of controlling the same | |
CN105352528B (zh) | 一种应用于弹道导弹的惯导误差在线补偿方法 | |
CN101598545A (zh) | 捷联系统加速度计安装误差角直接测量方法 | |
RU2331890C2 (ru) | Способ автономного измерения векторов скорости и ускорения |
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 |