CN112286059A - 一种具有攻击角和视场角约束的二阶前置角重塑制导方法 - Google Patents
一种具有攻击角和视场角约束的二阶前置角重塑制导方法 Download PDFInfo
- Publication number
- CN112286059A CN112286059A CN202011243991.8A CN202011243991A CN112286059A CN 112286059 A CN112286059 A CN 112286059A CN 202011243991 A CN202011243991 A CN 202011243991A CN 112286059 A CN112286059 A CN 112286059A
- Authority
- CN
- China
- Prior art keywords
- missile
- angle
- guidance
- formula
- attack
- 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
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明属于导弹制导系统设计领域,涉及一种具有攻击角和视场角约束的二阶前置角重塑制导方法,包括:建立导弹拦截目标的制导方程;将制导方程转化为以导弹与目标相对距离为自变量的微分方程组,得到制导律表达式,同时设计基于相对距离的二次函数,作为导弹前置角变化曲线;基于边界条件和攻击角度约束求解导弹前置角变化曲线中的未知参数,得到解析式制导律;通过限制导弹前置角曲线的最大幅值,实现导引头视场角约束,并得到攻击角度可达域。本发明采用基于前置角重塑的解析方法,不需要对模型进行线性化,也不需要对制导参数进行反复试凑,适用于具有攻击角度约束和导引头视场受限的制导任务,具有设计使用简单、可靠性高、适用性强等特点。
Description
技术领域
本发明属于导弹制导系统设计领域,具体来说,涉及一种具有攻击角和导引头视场角约束的二阶前置角重塑制导方法。
背景技术
制导系统是导弹的核心技术之一,当导弹的气动外形及动力系统确定后,制导系统在很大程度上决定了导弹的作战性能。制导系统的基本任务是产生合适的制导指令使得导弹获得最小的末端脱靶量。这一类任务指标单一,设计较为简单,方法也已经成熟。然而,现代战争赋予导弹越来越复杂的作战任务,终端脱靶量已不仅仅是衡量完成军事任务的唯一指标。针对于采用了复合装甲的新型坦克,配备了先进防空系统的军舰等防御能力强的目标,不仅希望导弹能够实现脱靶量最小,还希望能够以一定的攻击角度实现对目标的打击,从而达到更高的作战效能。因此针对这一问题,需要开展具有攻击角度约束的导弹制导方法的研究。
导弹在末段导引飞行过程中,为了防止丢失目标,需要保证导引头能一直锁定目标,这就需要限制导弹的前置角在导引头的视场范围以内,这种约束称为导引头的视场约束。而且,随着全捷联雷达导引头的工程化应用,该类导引头的视场范围较为狭小,且制导系统需要的部分制导信息无法通过该类导引头直接测得,因此,需要针对这类导弹武器系统设计新型的制导方法。目前,同时满足攻击角度约束、导引头视场约束的制导方法是国内外的研究热点。已有的大部分具有攻击角度约束和视场约束的制导方法,往往需要将导弹和目标相对运动的数学模型进行线性化,或者需要假设导弹的飞行速度是个常值,这些给所设计的制导方法带来了一定的局限性,降低了方法的适用性。
因此,有必要设计一种具有攻击角度和导引头视场约束的解析制导律,在克服已有保守性的同时,提高制导系统的适用性、可靠性、简化设计方法等方面。
发明内容
为了克服上述现有方法的不足和缺点,本发明设计了一种具有攻击角度和视场角约束的二阶前置角重塑制导方法,首先,建立导弹在纵向平面内拦截静止目标的制导方程;其次,将制导方程转化为以导弹与目标之间相对距离为自变量的微分方程组,得到制导律的表达式,同时,设计一个基于导弹与目标之间相对距离的二次函数,以此作为导弹前置角的变化曲线;然后,基于边界条件和攻击角度约束求解导弹前置角变化曲线中的未知参数,得到具有解析形式的制导律;最后,通过限制导弹前置角曲线的最大幅值,实现导引头视场角约束,并进一步得到攻击角度可达域。与已有方法不同,本发明采用的解析制导方法,在克服已有保守性、提高制导系统的适用性、可靠性的同时,不仅可以实现攻击角度和导引头视场角约束,而且不需要对导弹和目标相对运动的数学模型进行线性化、不需要做导弹常值速度假设、不需要对制导参数的进行反复试凑,适用于具有攻击角度约束和捷联导引头视场受限的制导任务,具有设计方法简单、适用性强等特点。
为实现上述目的,本发明提供了一种具有攻击角度和视场角约束的二阶前置角重塑制导方法,包括如下步骤:
S1:考虑导弹与静止目标所在的纵向平面为攻击平面,建立导弹拦截静止目标的制导方程;
S2:将步骤S1中建立的制导方程转化为以导弹与静止目标间相对距离为自变量的微分方程组,得到制导律的表达式,同时设计基于所述相对距离的二次函数,得到导弹前置角变化曲线;
S3:基于边界条件和攻击角度约束求解步骤S2中得到的导弹前置角变化曲线中的未知参数,得到具有解析形式的制导律;
S4:通过限制步骤S2中得到的导弹前置角变化曲线的最大幅值,实现导引头视场角约束,进而得到攻击角度可达域。
进一步,步骤S1中所建立的导弹拦截静止目标的制导方程为:
其中,R表示导弹和静止目标的相对距离,VM为导弹的飞行速度,σM为导弹的前置角,λ为导弹相对于静止目标的视线角,γM为导弹的弹道倾角,g表示重力加速度,aM表示导弹的法向加速度,也是需要进行设计的制导律;上标·表示一阶求导;m表示导弹的质量,T为导弹发动机的推力,D为导弹所受的空气阻力。
进一步,步骤S2具体过程为:
将制导方程(1)转化基于导弹与静止目标的相对距离R的微分方程组:
由式(3),得到制导律表达式为:
设计以相对距离R为自变量的二次函数:
将式(5)作为导弹前置角变化曲线,其中,相对距离R的定义域为R∈[0,R0],R0是导弹与静止目标的初始相对距离,a,b,c是需要设计的参数。
进一步,步骤S3具体过程为:
导弹前置角变化曲线的边界条件如下:
结合式(5)、式(6),得到:
求式(7)关于相对距离R的一阶导数,得
定义期望攻击角度为λd,对式(2)从R=R0到R=0进行积分,得到:
将tanσM项在σM=0处展开成泰勒级数,忽略三次以上的高阶项后,将其代入式(9),并结合式(7),得
其中,λd,λ0,σM0均为已知量,则参数a直接通过求解式(10)得到,将式(8)代入式(4),得到解析形式的制导律为:
进一步,步骤S4具体过程为:
结合式(7)和式(12),计算得到参数a的范围
Γa={a|amin≤a≤amax} (13)
将式(7)代入式(9),得
求式(14)关于参数a的一阶导数,得
则结合式(10)、式(13)和式(15)得到期望攻击角度λd的可达域如下:
λmax≤λd≤λmax (16)
其中,λmin和λmax分别为,
本发明的有益效果:
1)本发明中所设计的解析形式的制导律,设计方法简单、工程实用性较强,不仅可以实现攻击角度的约束,同时还可以保证导引头视场约束;
2)本发明采用基于前置角重塑的解析方法,不需要对导弹和目标相对运动的数学模型进行线性化,不需要做导弹常值速度假设,不需要对制导参数进行反复试凑,具有设计使用方法简单、适用性强等特点。
附图说明
图1为本发明实施例的具有攻击角和导引头视场角约束的二阶前置角重塑制导方法流程图;
图2为本发明实施例的导弹拦截地面静止目标的纵向平面示意图;
图3为本发明是实施例的制导律在不同攻击角度约束下的仿真示意图。
具体实施方式
下面将结合附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。
如图1所示,本实施例的具有攻击角和导引头视场角约束的二阶前置角重塑制导方法,包括如下步骤:
S1:考虑导弹与静止目标所在的纵向平面为攻击平面,建立导弹拦截静止目标的制导方程;
如图2所示,M与T分别表示导弹与静止目标,根据二维平面内导弹和静止目标之间的运动学和动力学关系,得到导弹在攻击平面内拦截静止目标的制导方程为:
其中,R表示导弹和目标的相对距离,VM为导弹的飞行速度,σM为所述导弹的前置角,λ表示导弹相对于目标的视线角,γM为导弹的弹道倾角,aM表示导弹的法向加速度,也是需要进行设计的制导律,g表示重力加速度,m表示导弹的质量,T为导弹发动机的推力,D表示导弹所受的空气阻力。
S2:将制导方程(1)转化为以导弹与目标的相对距离为自变量的微分方程组,得到制导律表达式,同时,设计基于导弹与目标相对距离的二次函数,作为期望的导弹前置角变化曲线。具体过程如下:
将制导方程(1)转化基于导弹与目标的相对距离R的微分方程组:
由式(3),可以得到制导律表达式为:
同时,设计以相对距离R为自变量的二次函数,作为导弹前置角的变化曲线,如下所示:
其中,自变量R的范围是R∈[0,R0],R0是导弹与目标的初始相对距离,a,b,c是需要设计的参数。
S3:基于边界条件和攻击角度约束求解导弹前置角变化曲线中的待定参数,得到具有解析形式的制导律。具体过程如下:
考虑导弹前置角变化曲线(式(5))的边界条件如下:
通过式(6)可以得到:
结合式(5)、式(6-1),可以得到:
求式(7)关于相对距离R的一阶导数,可得
考虑期望的攻击角度为λd,对式(2)从R=R0到R=0进行积分,得到:
将tanσM项在σM=0处展开成泰勒级数,忽略三次以上的高阶项,将其代入式(9),并结合式(7),可以得到
式(9-1)整理可得
其中,λd,λ0,σM0均为已知量,则参数a可以直接通过求解式(10)得到。将式(8)代入式(4),得到解析形式的制导律为:
S4:通过限制导弹前置角曲线的最大幅值,实现导引头视场角约束,并得到前置角曲线的参数a的取值范围,结合参数a与攻击角度的关系,得到攻击角度可达域。具体过程如下:
结合式(12)和式(12-1),可以解得参数a的一个安全取值范围是:
结合式(12)和式(12-3),可以解得参数a的一个安全取值范围是
综上,可以得到参数a的安全取值范围为
将式(7)代入式(9),可以得到
求λd关于参数a的一阶导数,可得
即,期望攻击角度λd是关于参数a的单调增函数,则结合式(10)、式(13)、式(15)可以得到期望攻击角度λd的可达域如下:
λmax≤λd≤λmax (16)
其中,
基于以上分析,当选择期望攻击角度约束在式(16)的范围内时,采用式(11)所示的解析制导律时,则可以实现对静止目标的精确打击,并同时满足攻击角度和导引头视场约束。
应当理解,本文中所有变量上标“·”都是该变量的导数,除非该变量的导数有实际物理含义。
下面以某反坦克导弹拦截静止目标为例,说明本发明所提出的方法的有效性。定义导弹的初始速度为200m/s,导弹初始位置为(0,0)km,目标位于(10,0)km处,导弹最大法向加速度限制为100m/s2,导弹导引头的最大视场角限定为45度(deg),需要设置导弹初始的前置角小于45度以便使得导弹的前置角在整个制导过程中满足视场约束,本实施例选取初始前置角为30度。针对该地面静止目标的拦截仿真示意图如图3所示,其中攻击角度约束选取为:0deg,-30deg,-60deg,-90deg。图3中的图a)-d)分别表示导弹的飞行轨迹、导弹的法向加速度(制导律)的变化趋势,导弹视线角的变化趋势、导弹前置角的变化趋势。从图3可以看出,本发明所提出的制导方法可以以期望的攻击角度成功地拦截静止目标,并保证了导弹前置角在制导过程中始终小于所设定的导弹导引头的最大视场角,从而实现了导引头视场约束。
综上,本发明所提出的解析制导律有效地实现了对攻击角度和导引头视场约束的要求,不需要对模型进行线性化,不需要做导弹常值速度的假设,不需要对制导参数进行反复试凑,具有较好的实用性、较强适应性、较高可靠性。
对于本领域的普通技术人员来说,在不脱离本发明创造构思的前提下,还可以对本发明的实施例做出若干变型和改进,这些都属于本发明的保护范围。
Claims (5)
1.一种具有攻击角度和视场角约束的二阶前置角重塑制导方法,其特征在于,包括如下步骤:
S1:考虑导弹与静止目标所在的纵向平面为攻击平面,建立导弹拦截静止目标的制导方程;
S2:将步骤S1中建立的制导方程转化为以导弹与静止目标间相对距离为自变量的微分方程组,得到制导律的表达式,同时设计基于所述相对距离的二次函数,得到导弹前置角变化曲线;
S3:基于边界条件和攻击角度约束求解步骤S2中得到的导弹前置角变化曲线中的未知参数,得到具有解析形式的制导律;
S4:通过限制步骤S2中得到的导弹前置角变化曲线的最大幅值,实现导引头视场角约束,进而得到攻击角度可达域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011243991.8A CN112286059B (zh) | 2020-11-10 | 2020-11-10 | 一种具有攻击角和视场角约束的二阶前置角重塑制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011243991.8A CN112286059B (zh) | 2020-11-10 | 2020-11-10 | 一种具有攻击角和视场角约束的二阶前置角重塑制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112286059A true CN112286059A (zh) | 2021-01-29 |
CN112286059B CN112286059B (zh) | 2022-01-28 |
Family
ID=74351121
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011243991.8A Active CN112286059B (zh) | 2020-11-10 | 2020-11-10 | 一种具有攻击角和视场角约束的二阶前置角重塑制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112286059B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116502420A (zh) * | 2023-04-12 | 2023-07-28 | 北京理工大学 | 一种基于序列二次规划方法的框架角约束变权制导方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105021092A (zh) * | 2015-06-30 | 2015-11-04 | 北京航天长征飞行器研究所 | 一种捷联寻的导引头的制导信息提取方法 |
CN106681348A (zh) * | 2017-01-13 | 2017-05-17 | 西北工业大学 | 考虑全捷联导引头视场约束的制导控制一体化设计方法 |
CN109614756A (zh) * | 2019-01-08 | 2019-04-12 | 北京航空航天大学 | 一种具有攻击时间和导引头视场约束的制导律的解析方法 |
CN110414159A (zh) * | 2019-08-01 | 2019-11-05 | 北京航空航天大学 | 一种基于圆的渐开线的攻击角度约束制导方法 |
CN110645844A (zh) * | 2019-09-04 | 2020-01-03 | 南京理工大学 | 具有攻击角度约束的高速拦截制导方法 |
CN110879604A (zh) * | 2019-12-25 | 2020-03-13 | 中国人民解放军海军潜艇学院 | 一种带落角控制的飞行器航向导引方法 |
-
2020
- 2020-11-10 CN CN202011243991.8A patent/CN112286059B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105021092A (zh) * | 2015-06-30 | 2015-11-04 | 北京航天长征飞行器研究所 | 一种捷联寻的导引头的制导信息提取方法 |
CN106681348A (zh) * | 2017-01-13 | 2017-05-17 | 西北工业大学 | 考虑全捷联导引头视场约束的制导控制一体化设计方法 |
CN109614756A (zh) * | 2019-01-08 | 2019-04-12 | 北京航空航天大学 | 一种具有攻击时间和导引头视场约束的制导律的解析方法 |
CN110414159A (zh) * | 2019-08-01 | 2019-11-05 | 北京航空航天大学 | 一种基于圆的渐开线的攻击角度约束制导方法 |
CN110645844A (zh) * | 2019-09-04 | 2020-01-03 | 南京理工大学 | 具有攻击角度约束的高速拦截制导方法 |
CN110879604A (zh) * | 2019-12-25 | 2020-03-13 | 中国人民解放军海军潜艇学院 | 一种带落角控制的飞行器航向导引方法 |
Non-Patent Citations (4)
Title |
---|
HUAN JIANG 等: "Cooperative guidance with multiple constraints using convex", 《AEROSPACE SCIENCE AND TECHNOLOGY》 * |
李宏宇: "多约束下的制导炮弹最优末制导律研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
李晓宝等: "考虑攻击角度和视场角约束的自适应终端滑模制导律", 《控制与决策》 * |
王宁宇: "考虑落角/视场角约束高速空对地导弹", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116502420A (zh) * | 2023-04-12 | 2023-07-28 | 北京理工大学 | 一种基于序列二次规划方法的框架角约束变权制导方法 |
CN116502420B (zh) * | 2023-04-12 | 2024-01-30 | 北京理工大学 | 一种基于序列二次规划方法的框架角约束变权制导方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112286059B (zh) | 2022-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109597423B (zh) | 一种基于参考视线角信号的多约束末制导律的设计方法 | |
CN109506517B (zh) | 一种带约束的中制导弹道优化方法 | |
US20080206718A1 (en) | Apparatus, method and computer program product for weapon flyout modeling and target damage assessment | |
CA2585501C (en) | Method of determining a fire guidance solution | |
CN112069605B (zh) | 一种带有攻击时间约束的比例导引律设计方法 | |
CN108073742B (zh) | 基于改进粒子滤波算法的拦截导弹末段飞行状态估计方法 | |
CN112286059B (zh) | 一种具有攻击角和视场角约束的二阶前置角重塑制导方法 | |
CN114035616B (zh) | 一种飞行器对移动目标打击控制方法及系统 | |
Wan et al. | Suboptimal midcourse guidance with terminal-angle constraint for hypersonic target interception | |
CN114153143B (zh) | 一种导弹非奇异固定时间滑模制导律的设计方法 | |
CN112464451B (zh) | 基于作战仿真系统的防空导弹武器命中概率修正方法 | |
CN113639586B (zh) | 一种抗雷达关机的制导方法、系统及介质 | |
CN111832193B (zh) | 一种超视距空战中空中目标的威胁评估方法 | |
CN113359819A (zh) | 一种带有碰撞角约束和加速度限制的最优制导律 | |
CN106643298A (zh) | 一种基于预置命中点的大气层内反导拦截器中制导方法 | |
Shi et al. | Design of bias proportional navigation guidance law for motion target with impact angle constraint | |
Zuoe et al. | Study on vertical attack penetration probability of anti-ship missile | |
Yu et al. | Flight Trajectory Correction of Target Missile and Design of Tracking Guidance | |
Dutta et al. | Estimation of drag coefficient from radar-tracked flight data of a cargo shell | |
CN116125808A (zh) | 一种基于纯方位测量的多约束机动突防滑模制导方法 | |
Zhu | Survey of Interception Strategies for Near-space hypersonic vehicles | |
Baranowski et al. | The analysis of the 35 mm artillery projectile's motion model parameters' identification based on the recorded flight trajectory | |
CN112729015B (zh) | 一种基于拦截弹导引律辨识的机动突防方法 | |
Tan | A Novel Sliding Mode Guidance Law with Impact Angle and Time Constraints | |
Zhigang et al. | Study on Optimal Guidance Method for Intercepting Hypersonic Weapons |
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 |