CN106091816B - 一种基于滑模变结构理论的半捷联空空导弹制导方法 - Google Patents
一种基于滑模变结构理论的半捷联空空导弹制导方法 Download PDFInfo
- Publication number
- CN106091816B CN106091816B CN201610366190.8A CN201610366190A CN106091816B CN 106091816 B CN106091816 B CN 106091816B CN 201610366190 A CN201610366190 A CN 201610366190A CN 106091816 B CN106091816 B CN 106091816B
- Authority
- CN
- China
- Prior art keywords
- guidance
- missile
- target
- line
- sight
- 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
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000001133 acceleration Effects 0.000 claims description 17
- 230000000007 visual effect Effects 0.000 claims description 9
- 230000003247 decreasing effect Effects 0.000 claims description 5
- 230000004438 eyesight Effects 0.000 claims description 5
- 230000010355 oscillation Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 abstract description 3
- 238000001914 filtration Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 230000003190 augmentative effect Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000003416 augmentation Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 238000003331 infrared imaging Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
本发明公开了一种基于滑模变结构理论的半捷联空空导弹制导方法,包括以下几个步骤:步骤1:获取视线坐标系下弹目视线转率微分模型;步骤2:获取运动跟踪滑模制导律;步骤3:获取期望视线转率xd;本发明将sigmoid函数应用于运动跟踪滑模制导律的滑模面设计中,在制导初期增大弹目视线角速率,提高对目标运动状态的估计,在末制导阶段令弹目视线角速率为零,降低脱靶量,提高导弹命中精度。
Description
技术领域
本发明属于武器技术、制导方法领域,涉及空空导弹制导方法研究,具体涉及一种半捷联红外成像近距格斗空空导弹制导方法。
背景技术
随着未来空战中目标机动能力的不断增强,对导弹的制导能力提出了进一步的要求。为快速发现并跟踪目标,实现对目标的快速捕获和精确跟踪,将导弹的导引头由传统的稳定平台式改为半捷联方式,由于导引头的改变使得导弹在获取弹目相对运动状态方式上发生了改变,传统的制导律已不能满足制导要求,需要设计新型制导律。
传统的制导律如比例制导律和改进型增广比例制导律以弹目相对角速率作为制导信息进行制导律生成。由于导弹制导精度高要求弹目相对运动变化小,而半捷联体制下弹目相对角速率的滤波估计精度高又要求弹目相对运动变化大,二者相互制约。须提出一种既增强系统可观测性,改善系统滤波精度又提高导弹命中精度的制导律设计方法。
发明内容
本发明的目的是为了解决上述问题,提出一种基于滑模变结构理论的半捷联空空导弹制导方法。
本发明的一种基于滑模变结构理论的半捷联空空导弹制导方法,包括以下几个步骤:
步骤1:获取视线坐标系下弹目视线转率微分模型;
步骤2:获取运动跟踪滑模制导律;
步骤3:获取期望视线转率xd;
本发明的优点在于:
(1)将sigmoid函数应用于运动跟踪滑模制导律的滑模面设计中,在制导初期增大弹目视线角速率,提高对目标运动状态的估计,在末制导阶段令弹目视线角速率为零,降低脱靶量,提高导弹命中精度;
(2)采用运动跟踪滑模制导律均能保证制导信息滤波器对制导信息的良好估计,尤其是对弹目相对距离、相对速度、视线角速率以及目标法向机动加速度的估计精度满足制导需求,从而使得导弹制导精度得到保证;
(3)制导律通过其表达式中的变结构项来消除目标机动和参数估计误差所带来的影响。这样即便弹目相对距离、相对速度和目标法向加速度估计不准确,导弹仍有很高的制导精度。
附图说明
图1为本发明的方法流程图;
图2为弹目相对运动关系三维图;
图3为目标法向加速度估计误差;
图4为弹目视线角估计误差;
图5弹目视线角速率估计误差;
图6弹目相对速度估计误差。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是一种基于滑模变结构理论的半捷联空空导弹制导方法,流程如图1所示,包括以下几个步骤:
步骤1:获取视线坐标系下弹目视线转率微分模型;
三维空间内,建立导弹与目标间的相对运动关系如图2所示,其中:oxGyGzG为惯性空间坐标系,oxSySzS为弹目视线坐标系,拦截弹与坐标系原点重合,下标M和T分别表示导弹和目标,oxS为弹目视线LOS(line-of-sight)。
在视线坐标系下,根据理论力学中的哥氏定理,建立弹目相对运动关系模型:
其中,
as=[asx,asy,asz]T=[astx-asmx,asty-asmy,astz-asmz]T;
rs=[R,0,0]T;
式中:R为导弹与目标间的相对距离,为相对距离对时间的导数,即接近速率,astx,asty,astz分别为目标在视线坐标系下的加速度分量,asmx,asmy,asmz分别为导弹在视线坐标系下的加速度分量,asx,asy,asz分别为导弹与目标在视线坐标系下的相对加速度分量,ωx,ωy,ωz分别为视线角速率在视线坐标系下的分量。
将各量代入式(1)中并展开,得到视线角速率的微分方程如下:
考虑到ωx在拦截过程中很小,忽略ωxωz和ωxωz这两项,上式后可以简化为
比较其中只有关于导弹和目标相对加速度这一项的符号不同,为了方便后续推导,设弹目视线转率在视线坐标系下的微分方程为:
其中,当x取ωy时,am=-asmz,at=-astz,当x取ωz时,am=asmy,at=asty。
步骤2:根据滑模变结构理论,获取运动跟踪滑模制导律。
根据滑模变结构理论,为弹目视线转率x设定期望视线转率参考指令xd,则实际运动与预定运动的偏差为选取切换函数为
令滑模趋近律为
式中,k和ε均为大于零的常数,sgn是符号函数。
把上面两式结合起来得到
由的定义和弹目视线转率微分方程可得
其中,am和at分别为导弹和目标的机动加速度。
注意到综合得到运动跟踪滑模制导律(MTSMG)为:
步骤3:获取预定运动xd的规律。
由步骤2可知,通过设计预定运动xd的规律,利用运动跟踪滑模制导律可使弹目视线转率按照给定的规律变化,从而可在显著提高系统可观测性的同时,提高导弹的命中精度。在被动制导的情况下,对于运动跟踪滑模制导律需要提供弹目相对距离R、相对速度和目标法向加速度at的估计。制导律通过其表达式中的变结构项εsgne来消除目标机动和参数估计误差所带来的影响。这样即便弹目相对距离R、相对速度和目标法向加速度at估计不准确,导弹仍有很高的制导精度。
现在分析预定运动xd的规律,从提高系统可观测性的角度来看,希望弹目视线角转率不为零,且具有一定的幅值;从提高导弹末制导命中精度的要求来看,又希望弹目视线角转率在末制导时段为零。为此,引入sigmoid函数:
其中:e表示指数;
将频率与幅值随时间依次递减的正弦函数叠加至sigmoid函数中,构造如下xd变化规律:
其中:α表示幅值系数,其控制xd变化的幅值,h表示衰减系数,其控制xd幅值递减时的速度,β表示振荡频率,控制xd振荡频率。tgo0为制导过程中从导弹发射开始到减小弹目视线角转率的时间,用来控制xd在末制导到来之前幅值递减的时刻。
若记导弹从发射至命中目标历时为tgo,xd从幅值急剧减小至命中目标的时间为t0,则有tgo0=tgo-t0。
如此,通过调节参数构造合适的xd变化规律即可在满足系统可观测性要求的同时提高导弹命中精度。
从步骤1~3完成了整个制导方法的推导与设计过程。步骤2中的式(9)给出一种运动跟踪滑模制导律,该制导律能保证弹目视线转率跟踪上期望视线转率参考指令xd,在后续的步骤3中又给出了该参考指令xd的一种设计方法如式(11),该参考指令xd的变化规律不仅能满足系统可观测性要求,同时提高导弹末段的命中精度。结合公式(9)和(11)就得到了本发明提出的“一种基于滑模变结构理论的半捷联空空导弹制导方法”。
步骤4:仿真验证。
仿真的条件:假设弹体初始姿态角为θ0=ψ0=φ0=0deg,弹体俯仰偏航滚转3个方向的角速度扰动为幅值[10,30,30]deg,频率为1Hz的正弦波,光轴初始位置与视线轴一致,即初始失调角均为0deg,设制导律为标准比例导引法,导航比K=5,选取导弹初始位置为惯性坐标系原点,弹目视线所在铅垂面为惯性坐标系的纵向平面,弹目初始距离3000m,导弹速度为600m/s,目标速度为300m/s,失调角的测量噪声为0.1deg,在此条件上,研究多种态势下的滤波效果。首先假设拦截过程中,在目标在弹道坐标系下,目标采用如下三种机动方式以规避导弹的攻击。
方式1:目标采取常值形式的机动
方式2:目标采取斜坡形式的机动
方式3:目标采取正弦形式的机动
上式中,aTx,aTy,aTz为目标在弹道坐标系下3个方向的加速度,amax为目标最大加速度,t0为目标开始机动时间,T为正弦运动周期,g表示重力常数,目标最大过载为12g,假定交战时段内,目标纵向机动很小,法向和侧向机动概率一样,则认为目标纵向加速度为0,法向和侧向机动加速度amax最大约为8g。
仿真弹目初始态势如下表所示,选择3组初始态势,针对不同初始态势下的大机动目标进行拦截过程仿真。仿真终止条件为相对距离r≤40m(盲区距离rb=40m)。仿真算法采用定步长四阶龙格库塔法,仿真步长取为1ms。
弹目初始态势列表
通过上面的分析,为比较不同制导律的制导效果,选定扩展卡尔曼滤波器,在距离不可观的情况下,目标机动方式为常值机动,初始态势为目标弹道偏角为-30°,目标弹道倾角为105°下分别应用标准比例制导律、改进型增广比例制导律和运动跟踪滑模制导律并调整制导律参数进行仿真分析,得到滤波模型各个状态估计误差,如下图3、图4、图5、图6所示。
图3~5中的前两条曲线分别为比例导引律(PNG)和修正比例导引律(MAPNG),后三条曲线均为本发明提出的运动跟踪滑模制导律(MTSMG),只是设置参数略有不同。对仿真结果进行分析,可以看出,在距离不可观测的情况下,曲线1(PNG)得到的拦截弹道,会使得滤波器得到的弹目相对距离、相对速度、以及弹目视线转率和目标机动加速度等估计信息发散,错误的视线转率和目标机动加速度会导致最终脱靶量较高,如后文中脱靶量表所示。曲线2(MAPNG)得到的拦截弹道,会使滤波器发散速度有所降低,但是滤波器收敛速度仍然较慢,收敛精度较低。而采用MTSMG(曲线3,4,5)得到的拦截弹道,不仅能抑制滤波器的发散,而且滤波器收敛速度更快。同时从曲线3,4,5可以看出,采用MTSMG制导方法,通过选择不同的参数,滤波的收敛速度和估计精度有略微差异,但滤波精度都能满足制导系统要求。
为了验证不同情况下各制导律的命中精度,假定弹目相对距离不可观,滤波模型采用扩展卡尔曼滤波方法,分别在目标取3种不同机动方式、3种不同弹目初始态势下分别对标准比例制导律、改进型增广比例制导律和运动跟踪滑模制导律进行仿真,并计算制导结束后的脱靶量信息。
分别得到如下三个表格如下:
目标常值机动方式下不同初始态势不同制导律的脱靶量
目标斜坡机动方式下不同初始态势不同制导律的脱靶量
目标正弦机动方式下不同初始态势不同制导律的脱靶量
结果表示:三种目标机动方式和三种弹目初始态势下,运动跟踪滑模制导律的制导精度要明显好于改进型增广比例制导律,而改进型增广比例制导律的制导精度则要明显好于标准比例制导律。同时可以看出目标常值机动、弹目初始态势为态势1和态势3的情况下,标准比例制导律的脱靶量分别是22.83m和29.79m,在目标斜坡机动,弹目初始态势为态势3的情况下,标准比例制导律的脱靶量为52.81m。实际上在这些情况下由于滤波模型中对于标准比例制导律所需的弹目视线角速率估计误差过大,滤波值发散,最终导弹命中目标失败。而上述这些情况下,采用运动跟踪滑模制导律均能对滤波状态值进行良好的估计,尤其是对弹目相对距离、相对速度、视线角速率以及目标法向机动加速度的估计精度满足制导需求,从而使得导弹制导精度得到保证。
Claims (1)
1.一种基于滑模变结构理论的半捷联空空导弹制导方法,包括以下几个步骤:
步骤1:获取视线坐标系下弹目视线转率微分模型;
弹目视线转率在视线坐标系下的微分方程为:
其中,当x取ωy时,am=-asmz,at=-astz,当x取ωz时,am=asmy,at=asty,R为导弹与目标间的相对距离,为相对距离对时间的导数,astx,asty,astz分别为目标在视线坐标系下的加速度分量,asmx,asmy,asmz分别为导弹在视线坐标系下的加速度分量,ωx,ωy,ωz分别为视线角速率在视线坐标系下的分量;
步骤2:获取运动跟踪滑模制导律;
运动跟踪滑模制导律为:
其中:am和at分别为导弹和目标的机动加速度,k和ε均为大于零的常数,sgn是符号函数,xd为期望视线转率,为实际运动与预定运动的偏差,
步骤3:获取期望视线转率xd;
期望视线转率xd为:
其中:e表示指数,t表示时间,α表示幅值系数,其控制xd变化的幅值,h表示衰减系数,其控制xd幅值递减时的速度,β表示振荡频率,控制xd振荡频率;tgo0为制导过程中从导弹发射开始到减小弹目视线角转率的时间,用来控制xd在末制导到来之前幅值递减的时刻,tgo0=tgo-t0,tgo为导弹从发射至命中目标历时,t0为xd从幅值急剧减小至命中目标的时间;
通过调节式(11)中参数,获取符合要求的期望视线转率,完成半捷联空空导弹制导。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610366190.8A CN106091816B (zh) | 2016-05-27 | 2016-05-27 | 一种基于滑模变结构理论的半捷联空空导弹制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610366190.8A CN106091816B (zh) | 2016-05-27 | 2016-05-27 | 一种基于滑模变结构理论的半捷联空空导弹制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106091816A CN106091816A (zh) | 2016-11-09 |
CN106091816B true CN106091816B (zh) | 2018-01-02 |
Family
ID=57230232
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610366190.8A Active CN106091816B (zh) | 2016-05-27 | 2016-05-27 | 一种基于滑模变结构理论的半捷联空空导弹制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106091816B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107179021B (zh) * | 2017-06-14 | 2018-08-28 | 北京理工大学 | 一种驾束制导体制下多弹协同高精度制导控制方法 |
CN108052112B (zh) * | 2017-12-01 | 2020-10-02 | 哈尔滨工业大学 | 基于pn制导律辨识的多飞行器威胁度获取方法 |
CN110220416B (zh) * | 2019-05-15 | 2021-12-10 | 南京理工大学 | 一种自适应快速弹道跟踪制导方法 |
CN110186328B (zh) * | 2019-06-03 | 2020-04-03 | 北京航空航天大学 | 一种多目标拦截的导弹分配方法及系统 |
CN113917841B (zh) * | 2021-09-28 | 2024-02-09 | 中国星网网络应用有限公司 | 基于二阶滑模的前向拦截制导方法及系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB8825367D0 (en) * | 1988-10-29 | 1989-04-19 | British Aerospace | Semi-strapdown gyromagnetic compass |
CN103090728B (zh) * | 2013-01-07 | 2014-07-23 | 北京理工大学 | 一种基于滑模控制的带末角约束制导方法 |
CN103512426B (zh) * | 2013-09-06 | 2015-05-06 | 北京理工大学 | 一种次优的带末角约束制导方法 |
CN104266546B (zh) * | 2014-09-22 | 2015-10-07 | 哈尔滨工业大学 | 一种基于视线的有限时间收敛主动防御制导控制方法 |
-
2016
- 2016-05-27 CN CN201610366190.8A patent/CN106091816B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106091816A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106091816B (zh) | 一种基于滑模变结构理论的半捷联空空导弹制导方法 | |
CN106681348B (zh) | 考虑全捷联导引头视场约束的制导控制一体化设计方法 | |
CN110008502B (zh) | 考虑全捷联导引头视场约束的三维制导控制一体化设计方法 | |
CN103090728B (zh) | 一种基于滑模控制的带末角约束制导方法 | |
CN111551080B (zh) | 一种攻击时间控制的制导方法 | |
RU2530705C2 (ru) | Способ и система для оценки траектории движущегося тела | |
CN105486307B (zh) | 针对机动目标的视线角速率估计方法 | |
CN111591470B (zh) | 一种适应推力可调模式的飞行器精确软着陆闭环制导方法 | |
CN111692919B (zh) | 超近射程的飞行器精确制导控制方法 | |
CN105486308B (zh) | 估计弹目视线角速率的快收敛Kalman滤波器的设计方法 | |
Kim et al. | Guidance and control for missiles with a strapdown seeker | |
CN112648886B (zh) | 一种组合制导目标拦截方法及系统 | |
CN112525003B (zh) | 一种带落角约束的扩展比例导引方法 | |
CN110764523B (zh) | 基于抗饱和平滑变换的比例积分前置导引攻击目标的方法 | |
CN114815888B (zh) | 一种仿射形式的制导控制一体化控制方法 | |
CN109976380A (zh) | 基于卡尔曼滤波估计的隔离度辨识校正方法及系统 | |
CN110703793B (zh) | 采用姿态角测量的飞行器积分比例导引攻击机动目标的方法 | |
CN109827478A (zh) | 一种带落角约束与过载约束的制导方法 | |
CN115033024A (zh) | 基于攻击时间和角度共同约束的多导弹三维协同制导方法 | |
CN116974208B (zh) | 基于捷联导引头的旋翼无人机目标打击控制方法及系统 | |
CN114489101A (zh) | 一种无人机的末端制导控制方法及系统 | |
CN107422744B (zh) | 一种基于径向速度控制的交会时间控制方法 | |
CN110471283A (zh) | 一种带碰撞角约束的三维鲁棒制导律构建方法 | |
Wang et al. | A guidance and control design with reduced information for a dual-spin stabilized projectile | |
CN114216368B (zh) | 一种变速条件下弹群协同制导方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |