CN106556287B - 一种积分比例导引非线性修正方法 - Google Patents
一种积分比例导引非线性修正方法 Download PDFInfo
- Publication number
- CN106556287B CN106556287B CN201610944823.9A CN201610944823A CN106556287B CN 106556287 B CN106556287 B CN 106556287B CN 201610944823 A CN201610944823 A CN 201610944823A CN 106556287 B CN106556287 B CN 106556287B
- Authority
- CN
- China
- Prior art keywords
- guidance
- angle
- bullet
- overload
- time
- 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
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F41—WEAPONS
- F41G—WEAPON SIGHTS; AIMING
- F41G3/00—Aiming or laying means
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种积分比例导引非线性修正方法,该方法主要应用于采用角制导的制导律中。通过将非线性影响的偏差等价到对应的制导律角度初始值零位,并在初始角度上进行修正,从而将对应的指令生成中的过载非线性部分消除。该方法弥补了在高空过载响应过低和响应不补充分带来的弹道超调影响,尤其在飞行过程中空域跨度大、指令切换的过程中带来的非线性,可在导弹攻击稳定移动目标时采用。
Description
技术领域
本发明属于制导控制技术领域,具体涉及一种积分比例导引非线性修正方法,该方法保证了过载指令不会由于响应的非线性出现超调现象,实现对移动目标打击时积分比例导引头的弹道特性更接近比例导引的弹道特性。
背景技术
积分比例导引在原有比例导引角速度制导积分而来,使得角制导方式可应用于打击移动目标。以俯仰方向为例,其基本原理为:比例导引头表征为瞬时的速度矢量转动角速度是弹目实现转动角速度的2~6倍;而积分比例导引则取某一时刻为基准,记录下此时的弹目视线角q0和弹道倾角θ0,随着时间变化,弹目视线角和弹道倾角均出现了变化,则制导律满足弹道倾角增量Δθ是弹目视线角增量Δq的 2~6倍。传统对两种制导律的认识为,比例导引制导律的响应更加迅速,而积分比例导引制导律的则由于增量的变化需要一定时间的积累。但在工程的使用过程中发现,制导律响应的滞后只是其中一个方面,另一个方面是,对于任意一段弹道的Δθ和Δq必须严格服从比例项系数的响应关系,若未响应或响应不充分,则未响应对应的Δq变化会进入q0中形成错误的过载指令(比例导引实时刷新,前一时刻指令响应不充分不影响下一时刻),从而导致弹道的超调现象。
发明内容
针对现有技术的缺陷和迫切需求,本发明提供一种积分比例导引非线性修正方法,其目的在于,将积分比例导引制导指令和响应不充分导致的非线性过载指令影响进行修正,从而使得弹道特性更趋近于比例导引。
一种积分比例导引非线性修正方法,包括以下步骤:
(1)设置制导初值:
采用积分比例导引制导律的初始时刻t0,记录初始时刻发射坐标系下弹目视线角q0以及初始时刻发射坐标系下弹道倾角θ0;随着时间变化,记录当前时刻时间t的发射坐标系下弹目视线角q(t)和弹道倾角θ(t)、弹体当前的飞行速度Vm(t)、弹道系下的过载响应值ab(t);
(2)计算过载指令
式中,μ为速度矢量驾驶仪阻尼、N为比例权系数、Tg为指令响应时间常数;
(3)积分累积计算等效弹目视线角值:
(4)依据等效修正的弹目视线角修正过载指令:
本发明的有益技术效果体现在:
本发明公开了一种适用于积分比例导引非线性影响的修正方法,该方法主要应用于采用角制导的制导律中。通过将非线性影响的偏差等价到对应的制导律角度初始值零位,并在初始角度上进行修正,从而将对应的指令生成中的过载非线性部分消除。该方法弥补了在高空过载响应过低和响应不补充分带来的弹道超调影响,尤其在飞行过程中空域跨度大、指令切换的过程中带来的非线性,可在导弹攻击稳定移动目标时采用。
附图说明
图1为考虑制导回路的速度矢量驾驶仪框图;
图2为气动饱和非线性模型图;
图3为指令平滑非线性模型图;
图4为实例修正前后攻角和侧滑角对比图;
图5为实例弹目视线角及初值补偿前后对比图。
具体实施方式
一、理论推导
以俯仰方向为例,比例导引形式为:
N为比例权系数。
将等式两边的速度项约去,假定初始弹道倾角为θ0,t时刻弹道倾角为θ(t),初始弹目视线角为q0,t时刻弹目视线角为q(t),则:
θ(t)-θ0=N(q(t)-q0)
为使得左右等式成立,需弹道倾角的增量满足上述变化关系,因而,对应设计速度矢量驾驶仪如图1所示。
可以将速度矢量驾驶仪的回路阻尼μ设计为合理的值,在此基础上进一步得出速度矢量驾驶仪频率ω、制导时间常数Tg、弹体速度Vm和前向通道系数Kc的关系如下:
则,利用等效速度矢量驾驶仪一阶与二阶系数的对应关系,可求出通道系数和阻尼系数之间的关系为:
得过载指令为:
在任意瞬间,弹体过载与其产生的弹目视线角变化间的关系为:
qb(t)为弹体过载随时间累积产生的弹目视线角变化值。
当前过载产生的弹道矢量变化即为实际的弹道倾角变化,故:
那么:
则,弹目视线角未响应变化的部分对应的初值qf(t)为:
从而有,修正后的积分比例导引制导律为:
对于上述制导律可以看出,若弹体响应是充分的,则修正项本身值为零,因而该制导律设计的修正部分具备通用性,此外,出现非线性的情况主要为气动力饱和或指令平滑,分别如图2和图3所示,非线性部分越大,持续时间越久,影响越严重,使用上述修正制导律后,指令和响应不一致导致的非线性过载均可得到修正。
二、弹目视线度制导非线性过载指令修正方法
积分比例导引非线性修正方法,包括以下步骤:
(1)设置制导初值
采用积分比例导引制导律时,弹载计算机保存初始时刻t0,保存初始时刻发射坐标系下弹目视线角q0,保存初始时刻发射坐标系下弹道倾角θ0,随着时间变化,分别记录当前时刻时间t、发射坐标系下弹目视线角q(t)和弹道倾角θt。弹载计算机还需计算出弹体当前的飞行速度大小V(t)。利用惯导加速度计和陀螺解算弹道系下的过载响应值ab(t)。
导引头测量的弹目视线角信息相对于弹体系,使用过程中需要解算到发射坐标系。使用的方法为,将弹目视线角信息从弹体系解算到地球惯性坐标系,在从地球惯性坐标系转换到发射坐标系。
(2)过载指令解算
选取合适的速度矢量驾驶仪阻尼μ、比例权系数N,并根据自动驾驶仪设计指标确定指令响应时间常数Tg,结合弹载计算机计算的其它参量,采用下述公式进行过载指令计算:
Vm(t)为弹体速度.
(3)等效初始弹目视线角计算
利用积分器累积计算等效弹目视线角值为:
(4)过载指令修正
将等效修正的弹目视线角初始值,实时引入到过载指令计算中,此时,注意阻尼系数、指令响应时间常数以及比例权系数若为改变值,则保证修正过程中的值和过载解算的值一致。计算公式为:
仿真实例:
以高空再入弹道为例,在非线性气动饱和和指令切换的情况下,引入初值修正方法后,攻角和侧滑角的变化规律更符合比例导引制导律,修正现象明显如图4所示;补偿前后的弹目视线角和初值补偿如图5所示,补偿后的弹目视线角输出更加平顺。
Claims (1)
1.一种积分比例导引非线性修正方法,包括以下步骤:
(1)设置制导初值:
采用积分比例导引制导律的初始时刻t0,记录初始时刻发射坐标系下弹目视线角q0以及初始时刻发射坐标系下弹道倾角θ0;随着时间变化,记录当前时刻时间t的发射坐标系下弹目视线角q(t)和弹道倾角θ(t)、弹体当前的飞行速度Vm(t)、弹道系下的过载响应值ab(t);
(2)计算过载指令式中,μ为速度矢量驾驶仪阻尼,N为比例权系数,Tg为指令响应时间常数;
(3)计算等效初始弹目视线角
(4)依据等效修正的弹目视线角修正过载指令:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610944823.9A CN106556287B (zh) | 2016-11-02 | 2016-11-02 | 一种积分比例导引非线性修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610944823.9A CN106556287B (zh) | 2016-11-02 | 2016-11-02 | 一种积分比例导引非线性修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106556287A CN106556287A (zh) | 2017-04-05 |
CN106556287B true CN106556287B (zh) | 2018-04-03 |
Family
ID=58443662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610944823.9A Active CN106556287B (zh) | 2016-11-02 | 2016-11-02 | 一种积分比例导引非线性修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106556287B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107726921B (zh) * | 2017-08-30 | 2019-12-03 | 湖北航天技术研究院总体设计所 | 一种适用于目标机动条件下的主动雷达角制导方法 |
CN111766776B (zh) * | 2020-07-06 | 2022-02-15 | 中国人民解放军海军航空大学 | 一种采用非线性比例与积分型补偿预估器的前置导引方法 |
CN111708382B (zh) * | 2020-07-06 | 2022-02-15 | 中国人民解放军海军航空大学 | 一种基于非线性比例积分的飞行器导引方法 |
CN113587727B (zh) * | 2021-06-22 | 2023-04-14 | 北京航天飞腾装备技术有限责任公司 | 一种基于弹目视线变系数的比例导引方法、系统及介质 |
CN113587740B (zh) * | 2021-07-22 | 2023-06-30 | 北京航天飞腾装备技术有限责任公司 | 一种基于弹目视线角的被动反辐射导引方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5429322A (en) * | 1994-04-22 | 1995-07-04 | Hughes Missile Systems Company | Advanced homing guidance system and method |
CN102706217A (zh) * | 2012-04-17 | 2012-10-03 | 北京理工大学 | 一种控制多枚导弹攻击角度和攻击时间的方法 |
CN104019701A (zh) * | 2014-05-28 | 2014-09-03 | 中国人民解放军海军航空工程学院 | 一种直接力气动力复合控制方法与前向拦截制导方法 |
CN106020215A (zh) * | 2016-05-09 | 2016-10-12 | 北京航空航天大学 | 一种基于单步预测矩阵博弈的近距空战自主决策方法 |
-
2016
- 2016-11-02 CN CN201610944823.9A patent/CN106556287B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5429322A (en) * | 1994-04-22 | 1995-07-04 | Hughes Missile Systems Company | Advanced homing guidance system and method |
CN102706217A (zh) * | 2012-04-17 | 2012-10-03 | 北京理工大学 | 一种控制多枚导弹攻击角度和攻击时间的方法 |
CN104019701A (zh) * | 2014-05-28 | 2014-09-03 | 中国人民解放军海军航空工程学院 | 一种直接力气动力复合控制方法与前向拦截制导方法 |
CN106020215A (zh) * | 2016-05-09 | 2016-10-12 | 北京航空航天大学 | 一种基于单步预测矩阵博弈的近距空战自主决策方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106556287A (zh) | 2017-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106556287B (zh) | 一种积分比例导引非线性修正方法 | |
CN111692919B (zh) | 超近射程的飞行器精确制导控制方法 | |
CN103090728B (zh) | 一种基于滑模控制的带末角约束制导方法 | |
CN111351401B (zh) | 应用于捷联导引头制导飞行器的防侧偏制导方法 | |
Kim et al. | Time-to-go polynomial guidance with trajectory modulation for observability enhancement | |
US8606435B2 (en) | Method and a system for estimating a trajectory of a moving body | |
CN111591470B (zh) | 一种适应推力可调模式的飞行器精确软着陆闭环制导方法 | |
KR20160127734A (ko) | 관성 항법 시스템 | |
CN110032206B (zh) | 远程制导飞行器大落角攻顶控制方法及控制系统 | |
CN110764523B (zh) | 基于抗饱和平滑变换的比例积分前置导引攻击目标的方法 | |
CN110425943B (zh) | 面向变质心飞行器的工程化再入制导方法 | |
CN111324149A (zh) | 一种基于视线角信息与前置导引信息的复合制导方法 | |
CN103604316A (zh) | 一种用于多弹发射的弹道校正方法 | |
CA2886122A1 (en) | Flight vehicle autopilot | |
CN110895418B (zh) | 补偿舵机动力学滞后的低速旋转飞行器控制方法及系统 | |
CN113587740B (zh) | 一种基于弹目视线角的被动反辐射导引方法及系统 | |
Özkan et al. | Performance comparison of the notable acceleration-and angle-based guidance laws for a short-range air-to-surface missile | |
CN111221348B (zh) | 应用于远程制导飞行器的侧偏修正方法 | |
US9347750B2 (en) | Determination of angle of incidence | |
Viswanath et al. | Disturbance observer based sliding mode control for proportional navigation guidance | |
CN105987652A (zh) | 姿态角速率估算系统及应用其的弹药 | |
CN111273682B (zh) | 基于虚拟目标点的侧偏修正方法 | |
Özkan et al. | Comparison of the notable acceleration-and angle-based guidance laws for a short-range air-to-surface missile | |
CN113064443B (zh) | 增益在线调整方法及使用其的阻尼回路控制方法 | |
Fariz et al. | Missile Initial Engagement Determination and Terminal Phase Guidance |
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 |