CN110687801A - 一种基于前向制导的自适应积分滑模三维制导律设计方法 - Google Patents
一种基于前向制导的自适应积分滑模三维制导律设计方法 Download PDFInfo
- Publication number
- CN110687801A CN110687801A CN201911157906.3A CN201911157906A CN110687801A CN 110687801 A CN110687801 A CN 110687801A CN 201911157906 A CN201911157906 A CN 201911157906A CN 110687801 A CN110687801 A CN 110687801A
- Authority
- CN
- China
- Prior art keywords
- guidance
- sliding mode
- law
- dimensional
- target
- 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
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)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
Description
技术领域
本发明涉及前向制导控制领域,更具体地,涉及一种基于前向制导的自适应积分滑模三维制导律设计非方法。
背景技术
随着战争环境的复杂化、以及空天武器的多样化和高速飞行器技术的快速发展。尤其空天武器向着高速化、大机动及规避能力强等方向发展,为了对拦截这类新型攻击性空天武器,传统的反导拦截技术已经不再适用,需要现代防御导弹必须采用一定的特殊拦截方式,才能对这类高速类空天武器保证拦截精度以及达到更大的毁伤程度。
随着目标的速度越来越快,导弹和目标的相对速度也随之增高,使得末端拦截时间过短甚至只有数秒钟,导致传统杀伤方式很难保证成功拦截目标。目前的防御武器的拦截方式通常要求拦截者的速度大于目标的速度时采用尾追拦截方式,相反的采用迎头拦截。因此,传统的制导方式很难保证拦截概率。为了拦截高速大机动飞行器目标,提出一种拦截高速飞行器目标的前向制导方法。采用该种拦截方式,弹目之间的速度变的相对较小,末端制导时间变长,这给导弹提供了充分的调整时间,并且能够增大攻击区域。
发明内容
本发明利用自适应方法、积分滑模控制理论,提出了一种基于前向制导自适应积分滑模三维制导律,能够有效拦截高速机动目标。
为了实现上述目的,本发明的技术方案为:
一种基于前向制导的自适应积分滑模三维制导律设计方法,包括以下步骤:
S1、建立三维拦截高超声速飞行器的前向制导相对运动学模型;
S2、根据步骤S1建立的模型定义满足前向制导条件的滑模面;
S3、根据步骤S2定义的滑面模,设计自适应律积分三维制导律;
上述步骤S2定义满足前向制导条件的滑模面的具体实现过程为:
定义积分滑模面如下:
其中,k1为正常数,ei(0)为误差变量ei的初始值,i=1、2,s1与s2分别表示俯仰和偏航方向上的滑面模,α1为正常数且满足0<α1<1。
与现有技术相比,本发明的有益效果是:基于前向制导方式,能够实现当导弹速度低于目标速度时,具备对高速目标的拦截能力。利用积分滑模控制理论和自适应方法,使得设计的自适应积分滑模三维制导律在降低了***的稳态误差的同时,能够提高拦截高速机动目标的制导精度。
附图说明
图1为前向制导拦截示意图。
图2为三维空间的弹目相对运动几何示意图。
图3为在制导律(15)作用下的仿真结果示意图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步描述。
为了成功拦截高速大机动飞行器目标,本发明采用前向制导方法,其前向制导拦截示意图如图1所示。
根据图1所示,采用前向制导方式,则三维空间的弹目相对运动几何关系如图2所示,其中,T为目标,M为拦截弹。拦截高速机动目标的相对模型如下:
其中,Vt和Vm分别是目标的速度和拦截弹的速度。θL和φL分别是视线关于参考坐标系的仰角和方位角。θt和φt为目标速度矢量前置角。θm和φm为拦截弹速度矢量前置角。ayt和azt是目标的加速度,aym和azm是拦截弹的加速度。
为了有效拦截高速目标,则制导过程中要求下式成立:
θm=n1θt (8)
φm=n2φt (9)
其中,n1和n2均是大于1的常数。式(8)和(9)保证了θm和φm分别随着θt和φt衰减而衰减。
定义误差变量ei如下:
ei=xi-niyi (12)
其中,
3自适应滑模制导律设计
定义积分滑模面如下:
其中,k1为正常数,ei(0)为误差变量ei的初始值,0<α1<1。
对式(13)求导可得:
其中,niMi为干扰项且满足|niMi|≤dMi,dMi为未知正常数。
采用自适应律对干扰niMi进行在线估计,根据(14)设计自适应律积分三维制导律如下:
其中,k2,k3,γi和pi为正常数,0<β<1。
定理1:针对***(11),当外部扰动上界未知时,在控制器(15)和自适应律(16)的作用下,可得到如下结论:
(I)滑模变量si是实际有限时间稳定的;
(II)误差变量ei在有限时间内收敛到零的任意小区域内。
证明:选取Lyapunov函数:
其中,
对式(17)求导,利用式(15)和式(16)整理可得:
根据
由式(19),则式(18)可整理得:
根据式(21)和式(22),则(20)可重写为:
根据(23)可得***是实际有限时间稳定的。
由式(23)可得:
算例分析
为了验证制导律的有效性。制导***初始化参数为:弹目相对距离为5000m,目标位置为:(0m,0m,0m),拦截弹位置为(4820m,1020m,-870m),视线角值为θL=-10deg和φL=-12deg,导弹前置角为θm(0)=-20deg和φm(0)=-15deg,目标前置角为θt=-20deg和φt=-15deg,导弹的速度为1500m/s,目标的速度为2100m/s,目标的加速度为2g。
为表明本文所制导策略的鲁棒性,与文献(Liu K,Cao Y,Wang S,et al.Terminalsliding mode control for landing on asteroids based on double power reachinglaw[C]2015IEEE International Conference on Information andAutomation.Lijiang,China,2015:2444-2449)中的制导律简称为SMGL进行仿真对比,为简便对比,自适应积分滑模制导律简化为ATSMGL。制导参数选取:k2=0.23、k3=0.55、l1=0.01、γi=0.025、pi=0.01,β=0.65和α1=0.75,其仿真结果如图3所示。3(a)为弹目运动轨迹示意图,3(b)为相对距离曲线示意图,3(c)为θm和θt的变化曲线示意图,3(d)为φm和φt的变化曲线示意图,3(e)为导弹加速度示意图,3(f)为自适应参数示意图。
图3(a)-(b)分别给出了在制导律SMGL和自适应积终端滑模制导律ATSMGL作用下的弹目运动轨迹和相对距离曲线可知,这两种制导律均能成功拦截目标。从图3(c)-(d)给出的是θm、θt、φm和φt的曲线可看出,在ATSMGL作用下,θm与θt快速保持倍数关系且迅速随着θt收敛到零,而SMGL作用下θt收敛速度较慢。图3(e)为导弹加速度曲线,从图可知,ATSMGL比SMGL作用下加速度曲线变化连续光滑,抖振小。从图3(f)给出的是自适应曲线可知,自适应参数在较短时间区域稳态值。
以上所述的本发明的实施方式,并不构成对本发明保护范围的限定。任何在本发明的精神原则之内所作出的修改、等同替换和改进等,均应包含在本发明的权利要求保护范围之内。
Claims (4)
2.根据权利要求1所述的方法,其特征在于,所述步骤S1的实现过程为:
采用前向制导方式,以末制导段目标的初始位置为参考坐标原点建立模型如下:
其中,Vt和Vm分别是目标的速度和拦截弹的速度,θL和φL分别是视线关于参考坐标系的仰角和方位角,θt和φt为目标速度矢量关于视线坐标系的方向角,即目标速度矢量前置角,θm和φm为拦截弹的速度关于视线坐标系的方向角,即拦截弹速度矢量前置角,ayt和azt是目标的加速度,aym和azm是拦截弹的加速度,
为了有效拦截高速目标,则制导过程中要求下式成立:
θm=n1θt (8)
φm=n2φt (9)
其中,n1和n2是前向制导系数,均是大于1的常数,式(8)和(9)保证了θm和φm分别随着θt和φt衰减而衰减,
定义误差变量ei如下:
ei=xi-niyi (12)
4.根据权利要求3所述的方法,其特征在于,针对***(11),当外部扰动上界未知时,在控制器(15)和自适应律(16)的作用下,可得到如下结论:
(I)滑模变量si是实际有限时间稳定的;
(II)误差变量ei在有限时间内收敛到零的任意小区域内。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911157906.3A CN110687801A (zh) | 2019-11-22 | 2019-11-22 | 一种基于前向制导的自适应积分滑模三维制导律设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911157906.3A CN110687801A (zh) | 2019-11-22 | 2019-11-22 | 一种基于前向制导的自适应积分滑模三维制导律设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110687801A true CN110687801A (zh) | 2020-01-14 |
Family
ID=69117377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911157906.3A Pending CN110687801A (zh) | 2019-11-22 | 2019-11-22 | 一种基于前向制导的自适应积分滑模三维制导律设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110687801A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113917841A (zh) * | 2021-09-28 | 2022-01-11 | 东方红卫星移动通信有限公司 | 基于二阶滑模的前向拦截制导方法及*** |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508434A (zh) * | 2011-11-07 | 2012-06-20 | 河海大学常州校区 | 用于微陀螺仪的自适应模糊滑模控制器 |
CN102862686A (zh) * | 2012-09-28 | 2013-01-09 | 北京理工大学 | 再入飞行器的最优积分滑模姿态控制方法及控制器 |
CN104898431A (zh) * | 2015-06-10 | 2015-09-09 | 北京理工大学 | 一种基于扰动观测器的再入飞行器有限时间控制方法 |
CN104950898A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 |
CN106843265A (zh) * | 2016-12-30 | 2017-06-13 | 哈尔滨工业大学 | 有限时间收敛的三维多导弹协同制导方法与*** |
CN106934120A (zh) * | 2017-02-23 | 2017-07-07 | 哈尔滨工业大学 | 基于前向制导的拦截高超声速飞行器的三维制导律设计方法 |
CN108416098A (zh) * | 2018-02-02 | 2018-08-17 | 北京航空航天大学 | 一种拦截机动目标的攻击时间约束制导律设计方法 |
CN109085848A (zh) * | 2018-08-02 | 2018-12-25 | 西北工业大学 | 空空导弹直接力/气动力有限时间抗饱和控制方法 |
CN109557814A (zh) * | 2018-12-04 | 2019-04-02 | 四川航天***工程研究所 | 一种有限时间积分滑模末制导律 |
CN109977613A (zh) * | 2019-04-19 | 2019-07-05 | 哈尔滨工业大学 | 一种可预先设定调整时间的自适应滑模末制导律设计方法 |
-
2019
- 2019-11-22 CN CN201911157906.3A patent/CN110687801A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508434A (zh) * | 2011-11-07 | 2012-06-20 | 河海大学常州校区 | 用于微陀螺仪的自适应模糊滑模控制器 |
CN102862686A (zh) * | 2012-09-28 | 2013-01-09 | 北京理工大学 | 再入飞行器的最优积分滑模姿态控制方法及控制器 |
CN104898431A (zh) * | 2015-06-10 | 2015-09-09 | 北京理工大学 | 一种基于扰动观测器的再入飞行器有限时间控制方法 |
CN104950898A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 |
CN106843265A (zh) * | 2016-12-30 | 2017-06-13 | 哈尔滨工业大学 | 有限时间收敛的三维多导弹协同制导方法与*** |
CN106934120A (zh) * | 2017-02-23 | 2017-07-07 | 哈尔滨工业大学 | 基于前向制导的拦截高超声速飞行器的三维制导律设计方法 |
CN108416098A (zh) * | 2018-02-02 | 2018-08-17 | 北京航空航天大学 | 一种拦截机动目标的攻击时间约束制导律设计方法 |
CN109085848A (zh) * | 2018-08-02 | 2018-12-25 | 西北工业大学 | 空空导弹直接力/气动力有限时间抗饱和控制方法 |
CN109557814A (zh) * | 2018-12-04 | 2019-04-02 | 四川航天***工程研究所 | 一种有限时间积分滑模末制导律 |
CN109977613A (zh) * | 2019-04-19 | 2019-07-05 | 哈尔滨工业大学 | 一种可预先设定调整时间的自适应滑模末制导律设计方法 |
Non-Patent Citations (2)
Title |
---|
彭双春等: "适于高超声速飞行器的三维非线性滑模制导律", 《***工程与电子技术》 * |
耿洁等: "飞行器再入段最优自适应积分滑模姿态控制", 《宇航学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113917841A (zh) * | 2021-09-28 | 2022-01-11 | 东方红卫星移动通信有限公司 | 基于二阶滑模的前向拦截制导方法及*** |
CN113917841B (zh) * | 2021-09-28 | 2024-02-09 | 中国星网网络应用有限公司 | 基于二阶滑模的前向拦截制导方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110220416B (zh) | 一种自适应快速弹道跟踪制导方法 | |
CN106934120B (zh) | 基于前向制导的拦截高超声速飞行器的三维制导律设计方法 | |
CN106681348B (zh) | 考虑全捷联导引头视场约束的制导控制一体化设计方法 | |
CN110187640B (zh) | 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法 | |
CN105653827B (zh) | 高超声速飞行器Terminal滑模控制器设计方法 | |
CN111551080B (zh) | 一种攻击时间控制的制导方法 | |
CN107132765B (zh) | 一种基于轨迹规划的攻击角度与攻击时间控制方法 | |
CN112782984B (zh) | 多飞行器反拦截协同打击的制导方法、装置和飞行器 | |
CN114003050B (zh) | 一种基于微分博弈的三体对抗策略的主动防御制导方法 | |
CN104019701B (zh) | 一种利用直接力气动力复合控制的前向拦截制导方法 | |
CN114020021A (zh) | 一种多导弹分布式协同制导律的设计方法及*** | |
CN110411289B (zh) | 一种抑制强机弹干扰的分离稳定控制方法 | |
CN110618608B (zh) | 一种复合制导跟踪控制方法与装置 | |
Li et al. | Three‐Dimensional Impact Time and Angle Control Guidance Based on MPSP | |
CN110929216B (zh) | 一种具有抗拖曳的自适应反步制导律设计方法 | |
CN110687801A (zh) | 一种基于前向制导的自适应积分滑模三维制导律设计方法 | |
CN112379600A (zh) | 一种基于通信时变延迟的分布式协同制导律构建方法 | |
CN114153143B (zh) | 一种导弹非奇异固定时间滑模制导律的设计方法 | |
CN111176315B (zh) | 一种基于l1自适应控制的变体巡飞弹控制方法 | |
CN115406312A (zh) | 考虑视场角和舵机延时约束的导弹制导控制一体化方法 | |
Liu et al. | Cooperative guidance for active defence based on line-of-sight constraint under a low-speed ratio | |
AU731363B2 (en) | Neural network trajectory command controller | |
CN110471283B (zh) | 一种带碰撞角约束的三维鲁棒制导律构建方法 | |
CN114815878B (zh) | 基于实时优化和深度学习的高超声速飞行器协同制导方法 | |
CN114819055A (zh) | 一种基于lstm网络的导弹拦截点预测方法 |
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: 20200114 |
|
RJ01 | Rejection of invention patent application after publication |