CN116379847A - 一种导弹在线能量管理制导律构建方法 - Google Patents
一种导弹在线能量管理制导律构建方法 Download PDFInfo
- Publication number
- CN116379847A CN116379847A CN202310651183.2A CN202310651183A CN116379847A CN 116379847 A CN116379847 A CN 116379847A CN 202310651183 A CN202310651183 A CN 202310651183A CN 116379847 A CN116379847 A CN 116379847A
- Authority
- CN
- China
- Prior art keywords
- angle
- sideslip angle
- speed
- missile
- iteration
- 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.)
- Pending
Links
- 238000010276 construction Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 239000013598 vector Substances 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 230000002265 prevention Effects 0.000 abstract description 7
- 230000009471 action Effects 0.000 abstract description 2
- 230000009467 reduction Effects 0.000 abstract description 2
- 238000007726 management method Methods 0.000 description 17
- 230000008859 change Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F42—AMMUNITION; BLASTING
- F42B—EXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
- F42B15/00—Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
- F42B15/01—Arrangements thereon for guidance or control
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Combustion & Propulsion (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Aviation & Aerospace Engineering (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Chemical & Material Sciences (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明公开了一种导弹在线能量管理制导律构建方法,首先采用了变步长积分通过并行计算实现了末速在线预测;其次利用高斯牛顿法计算附加侧滑角幅值进行减速;然后设置速度追踪切换通道进行附加侧滑角方向切换;最后根据附加侧滑角输出量进行过载跟踪实现侧向机动突防并进行减速。该方法作用时间长,计算精度高,减速的同时可增加突防概率,工程实践能力强。
Description
技术领域
本发明属于导弹技术领域,具体涉及一种导弹在线能量管理制导律构建方法。
背景技术
已有的导弹能量管理技术是基于H-V剖面的方法,这一类方法根据当前高度对应的速度剖面进行比对得到bang-bang形式的附加侧滑或倾侧角进行减速。该形式能量管理精度低,且对于桑格尔弹道,H-V剖面高度不单调,因此通常用于末制导,减速能力有限,且大幅切换侧滑角无法满足导弹高精度能量管理的需求。故如何对导弹进行更大范围和更高精度的能量管理都有待研究。
发明内容
为了克服现有技术的不足,本发明提供了一种导弹在线能量管理制导律构建方法,首先采用了变步长积分通过并行计算实现了末速在线预测;其次利用高斯牛顿法计算附加侧滑角幅值进行减速;然后设置速度追踪切换通道进行附加侧滑角方向切换;最后根据附加侧滑角输出量进行过载跟踪实现侧向机动突防并进行减速。该方法作用时间长,计算精度高,减速的同时可增加突防概率,工程实践能力强。
本发明解决其技术问题所采用的技术方案包括如下步骤:
步骤1:末速预测;
建立速度预测微分方程:
其中为射程和速度;/>为导弹当前视线偏角, /> , />分别为标称发射系弹道倾角和当地弹道倾角;/>为引力加速度,/>为海平面引力加速度, m为地球半径,/>为高度指令;/>为第iter次迭代时速度系下的气动阻力,/>为声速;/>为第iter次迭代时的阻力系数,是和马赫数/>以及合攻角/>有关的函数;/>为气动修正系数, />,/>为声速,合攻角/>,其中/>为攻角,/>为第iter次迭代中当状态取/>时的侧滑角;
步骤2:构造剩余射程侧滑角;
剩余射程的侧滑角为未知量,为了保证速度控制余量,假设剩余射程中的侧滑角有如下形式:
步骤3:当前状态需用侧滑角求解侧滑角幅值迭代计算;
步骤4:需用侧滑角设计;
因此,经过迭代输出后,需用侧滑角符号由下式确定:
本发明的有益效果如下:
本发明方法保障了导弹中制导和末制导全程进行减速,并能随机换向,增加突防概率。该方法具有广阔的军事应用前景。
附图说明
图1为本发明方法逻辑框图。
图2为本发明实施例末速分布区间。
图3为本发明实施例末速统计分布。
图4为本发明实施例射程-侧向位置变化曲线。
图5为本发明实施例时间-侧滑角变化曲线。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
为解决导弹高精度能量管理的问题,需要研究一种减速能力强、减速精度高,且能快速计算的末速管理方法。本发明公开了一种导弹在线能量管理制导律,通过根据标称弹道作为参考设置一种预测校正方法。该方法创新性的采用在线并行计算机制,根据当前状态对对导弹落速进行预测,并附加侧滑角减速,最终实现导弹能量管理的同时也增加了突防概率。
本发明首先根据优化得到的低能量弹道作为标称弹道作为输入;其次,采用在线并行计算机制,通过变步长积分快速预测导弹末速;然后采用高斯牛顿法迭代计算当前时刻的附加侧滑角幅值进行减速;最后附加侧滑角的方向切换采用预设的速度追踪通道减速策略进行换向,从而在完成导弹高精度能量管理需求的同时增加了突防概率,算法逻辑框图如图1所示。该方法结构简单,易于实现,工程实践能力强,可广泛应用于各类要求导弹末段能量管理的导弹制导控制系统设计中。
步骤一,末速预测;
其中为射程和速度;/>为导弹当前视线偏角, /> , />分别为标称发射系弹道倾角和当地弹道倾角;/>为引力加速度,/>为海平面引力加速度, m为地球半径,/>为高度指令;/>为第iter次迭代时速度系下的气动阻力,/>为声速;/>为第iter次迭代时的阻力系数,是和马赫数/>以及合攻角/>有关的函数;/>为气动修正系数, />,/>为声速,合攻角/>,其中/>为攻角,/>为第iter次迭代中当状态取/>时的侧滑角;
步骤二,构造剩余射程侧滑角;
剩余射程的侧滑角为未知量,为了保证速度控制余量,假设剩余射程中的侧滑角有如下形式:
步骤三,当前状态需用侧滑角求解侧滑角幅值迭代计算;
其中:
结束迭代;
继续迭代;
其他;
结束;
步骤四,需用侧滑角设计;
步骤三通过末速预测值对当前时刻的需用侧滑角幅值进行迭代计算,在减速过程中需要导弹保持在指向目标的弹目连线上,因此设计换向通道使得视线角保持在指向目标的弹幕连线附近,并完成减速。设/>为迭代结束时的/>。
具体实施例:
结合某制导火箭武器系统算例对本发明做进一步描述。制导火箭发射点海拔0m,目标海拔0m,对射程1200km 处目标进行垂直攻击,落速要求750m/s。
一、末速预测;
设置微分方程:
中需要装订的弹道参考变量:,/> , />,/>,/>分别为标称高度,标称弹道倾角,标称当地弹道倾角,标称攻角,气动修正系数,根据装订的标称射程/>插值得到;根据高度/>和装订的大气表计算对应大气密度/>和声速/>;根据当前速度和声速计算马赫数。
二、构造剩余射程侧滑角;
根据剩余射程侧滑角构造公式:
设置步骤一中的积分侧滑角。其中为预设末段侧滑角幅值;根据攻角/>和侧滑角/>计算合攻角/>,结合气动修正系数/>,马赫数/>,计算气动阻力系数/>;根据气动阻力系数,大气密度和速度计算阻力/>带入步骤一中的微分方程中。
三、当前状态需用侧滑角求解侧滑角幅值迭代计算;
四、需用侧滑角设计;
为了验证本发明的可行性,分别进行了本发明提出的能量管理制导律进行200次蒙特卡洛打靶数学仿真,并采用6自由度导弹动力学方程对所求需用侧滑角进行跟踪响应,得到数学仿真结果图2~图5的曲线。由仿真结果可知:本发明设计的导弹能量管理技术通过弹载计算机并行计算对末速进行实时预测,通过迭代计算得到当前需用侧滑角,通过导弹控制内环响应需用侧滑角实现能量管理,在导弹高精度能量管理并在中制导中实现侧向随机机动,增加了导弹突防概率。
Claims (2)
1.一种导弹在线能量管理制导律构建方法,其特征在于,包括如下步骤:
步骤1:末速预测;
建立速度预测微分方程:
其中为射程和速度;/>为导弹当前视线偏角, /> , />分别为标称发射系弹道倾角和当地弹道倾角;/>为引力加速度,/>为海平面引力加速度, m为地球半径,/>为高度指令;/>为第iter次迭代时速度系下的气动阻力,/>为声速;/>为第iter次迭代时的阻力系数,是和马赫数/>以及合攻角/>有关的函数;/>为气动修正系数, />,/>为声速,合攻角/>,其中/>为攻角,/>为第iter次迭代中当状态取/>时的侧滑角;
步骤2:构造剩余射程侧滑角;
剩余射程的侧滑角为未知量,为了保证速度控制余量,假设剩余射程中的侧滑角有如下形式:
步骤3:当前状态需用侧滑角求解侧滑角幅值迭代计算;
步骤4:需用侧滑角设计;
因此,经过迭代输出后,需用侧滑角符号由下式确定:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310651183.2A CN116379847A (zh) | 2023-06-05 | 2023-06-05 | 一种导弹在线能量管理制导律构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310651183.2A CN116379847A (zh) | 2023-06-05 | 2023-06-05 | 一种导弹在线能量管理制导律构建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116379847A true CN116379847A (zh) | 2023-07-04 |
Family
ID=86969742
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310651183.2A Pending CN116379847A (zh) | 2023-06-05 | 2023-06-05 | 一种导弹在线能量管理制导律构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116379847A (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6685143B1 (en) * | 2003-01-03 | 2004-02-03 | Orbital Research Inc. | Aircraft and missile forebody flow control device and method of controlling flow |
US20060169833A1 (en) * | 2003-05-19 | 2006-08-03 | Giat Industries | Process to control the trajectory of a spinning projectile |
US7513455B1 (en) * | 2005-02-18 | 2009-04-07 | Lockhead Martin Corporation | Ballistic missile interceptor guidance by acceleration relative to line-of-sight |
CN105043171A (zh) * | 2015-06-30 | 2015-11-11 | 北京航天长征飞行器研究所 | 一种带倾角约束的火箭弹纵向导引方法 |
CN111209521A (zh) * | 2020-01-23 | 2020-05-29 | 西安现代控制技术研究所 | 一种单平面双通道气动数据的双平面三通道扩展方法 |
CN111397441A (zh) * | 2019-01-03 | 2020-07-10 | 北京理工大学 | 带有捷联导引头的远程制导飞行器的全射程覆盖制导系统 |
CN111442697A (zh) * | 2020-02-07 | 2020-07-24 | 北京航空航天大学 | 一种基于伪谱法修正的过重补制导方法和弹道整形制导方法 |
CN112762776A (zh) * | 2021-01-22 | 2021-05-07 | 北京理工大学 | 一种火箭弹末端速度估计方法 |
CN115270055A (zh) * | 2022-07-26 | 2022-11-01 | 北京航空航天大学 | 一种运载火箭助推段大攻角弹道解析预测方法及装置 |
CN116049993A (zh) * | 2023-02-16 | 2023-05-02 | 北京航空航天大学 | 运载火箭助推段三维解析制导方法、装置及电子设备 |
-
2023
- 2023-06-05 CN CN202310651183.2A patent/CN116379847A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6685143B1 (en) * | 2003-01-03 | 2004-02-03 | Orbital Research Inc. | Aircraft and missile forebody flow control device and method of controlling flow |
US20060169833A1 (en) * | 2003-05-19 | 2006-08-03 | Giat Industries | Process to control the trajectory of a spinning projectile |
US7513455B1 (en) * | 2005-02-18 | 2009-04-07 | Lockhead Martin Corporation | Ballistic missile interceptor guidance by acceleration relative to line-of-sight |
CN105043171A (zh) * | 2015-06-30 | 2015-11-11 | 北京航天长征飞行器研究所 | 一种带倾角约束的火箭弹纵向导引方法 |
CN111397441A (zh) * | 2019-01-03 | 2020-07-10 | 北京理工大学 | 带有捷联导引头的远程制导飞行器的全射程覆盖制导系统 |
CN111209521A (zh) * | 2020-01-23 | 2020-05-29 | 西安现代控制技术研究所 | 一种单平面双通道气动数据的双平面三通道扩展方法 |
CN111442697A (zh) * | 2020-02-07 | 2020-07-24 | 北京航空航天大学 | 一种基于伪谱法修正的过重补制导方法和弹道整形制导方法 |
CN112762776A (zh) * | 2021-01-22 | 2021-05-07 | 北京理工大学 | 一种火箭弹末端速度估计方法 |
CN115270055A (zh) * | 2022-07-26 | 2022-11-01 | 北京航空航天大学 | 一种运载火箭助推段大攻角弹道解析预测方法及装置 |
CN116049993A (zh) * | 2023-02-16 | 2023-05-02 | 北京航空航天大学 | 运载火箭助推段三维解析制导方法、装置及电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109557814B (zh) | 一种有限时间积分滑模末制导律 | |
CN108362171B (zh) | 一种具有攻击时间和攻击角度约束的制导策略 | |
Cho et al. | Nonsingular sliding mode guidance for impact time control | |
CN103090728B (zh) | 一种基于滑模控制的带末角约束制导方法 | |
CN106647287B (zh) | 一种基于自适应动态规划的输入受限微分对策制导方法 | |
JP5755795B2 (ja) | ミサイルの針路変更を最小化する誘導システムおよび方法 | |
CN109737812B (zh) | 空对地制导武器侧向攻击方法和装置 | |
Kumar et al. | Finite-horizon robust suboptimal control-based impact angle guidance | |
CN109084641B (zh) | 导弹制导方法及装置 | |
CN112305919A (zh) | 一种带碰撞角约束的固定时间滑模制导律设计方法 | |
CN112859921B (zh) | 一种三维同时攻击鲁棒协同制导律设计方法 | |
CN117171877A (zh) | 基于时机博弈的高超声速飞行器机动突防策略设计方法 | |
Fonod et al. | Blinding guidance against missiles sharing bearings-only measurements | |
CN108073742A (zh) | 基于改进粒子滤波算法的拦截导弹末段飞行状态估计方法 | |
CN113359819B (zh) | 一种带有碰撞角约束和加速度限制的最优制导律 | |
CN116379847A (zh) | 一种导弹在线能量管理制导律构建方法 | |
CN115686059A (zh) | 一种基于伪谱法的高超声速飞行器禁飞区规避制导方法 | |
Wang et al. | Autopilot and guidance law design considering impact angle and time | |
CN118244804A (zh) | 一种采用附加侧滑角在线快速迭代速度控制方法 | |
CN115268501A (zh) | 多飞行器协同再入轨迹规划方法、系统、电子设备及介质 | |
CN112346474B (zh) | 一种有限时间收敛的微分对策制导律的设计方法 | |
CN114815878A (zh) | 基于实时优化和深度学习的高超声速飞行器协同制导方法 | |
Shi et al. | Predictive guidance strategies for active aircraft defense | |
Saleem et al. | Three dimensional nonlinear guidance law for exact impact time control | |
Khaikov | Estimate of projectile initial velocity as a solution of a two-point boundary value problem |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20230704 |