[go: up one dir, main page]

CN114995495A - Time and angle constraint three-dimensional sliding mode guidance law design method under speed uncontrollable condition - Google Patents

Time and angle constraint three-dimensional sliding mode guidance law design method under speed uncontrollable condition Download PDF

Info

Publication number
CN114995495A
CN114995495A CN202210600525.3A CN202210600525A CN114995495A CN 114995495 A CN114995495 A CN 114995495A CN 202210600525 A CN202210600525 A CN 202210600525A CN 114995495 A CN114995495 A CN 114995495A
Authority
CN
China
Prior art keywords
angle
missile
sliding mode
sight
guidance
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
Application number
CN202210600525.3A
Other languages
Chinese (zh)
Other versions
CN114995495B (en
Inventor
张科
游骏
孟源洁
张今
弋可行
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202210600525.3A priority Critical patent/CN114995495B/en
Publication of CN114995495A publication Critical patent/CN114995495A/en
Application granted granted Critical
Publication of CN114995495B publication Critical patent/CN114995495B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/107Simultaneous control of position or course in three dimensions specially adapted for missiles
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

The invention relates to a time and angle constraint three-dimensional sliding mode guidance law design method under the condition of uncontrollable speed, and relates to the technical field of missile guidance. The guidance law designed by the invention is perpendicular to the speed direction, and the expected attack time and attack angle are realized only by changing the speed direction of the missile. In a pitching channel, a rapidly-converged nonsingular sliding mode surface and a sliding mode approach law which changes along with the distance of a missile are provided, and an expected sight elevation angle is realized; in a yaw channel, a special sliding mode surface is designed by analyzing the guidance conditions of the yaw channel, and the missile can be ensured to reach the sliding mode surface by the supercoiled sliding mode theory, so that the expected sight azimuth angle and the expected attack time are realized. Under the condition that the speed of the missile is uncontrollable, the missile strikes a target at the expected attack time and the attack angle in the three-dimensional space, the guidance law is small in miss distance, the expected attack time and the expected sight elevation angle can be accurately achieved, and the sight azimuth angle has small errors.

Description

速度不可控条件下时间和角度约束三维滑模制导律设计方法Design method of three-dimensional sliding mode guidance law with time and angle constraints under uncontrollable speed conditions

技术领域technical field

本发明涉及导弹制导技术领域,尤其涉及在速度不可控条件下控制导弹攻击时间和角度的制导律设计。The invention relates to the technical field of missile guidance, in particular to the design of a guidance law for controlling missile attack time and angle under the condition of uncontrollable speed.

背景技术Background technique

对于现有的三维制导律,主要存在两种制导模型,第一种基于视线坐标系建立,在该模型的所求出的制导指令在转换到弹道坐标系后,存在速度方向的分量,需要飞行器能够控制导弹速度。但在导弹的末制导阶段,导弹通常无推力,不能控制导弹速度。针对这一问题,第二种模型基于导弹前置角坐标系建立,制导指令垂直于导弹速度方向,因此不控制速度。For the existing three-dimensional guidance law, there are mainly two guidance models. The first one is based on the line-of-sight coordinate system. After the guidance command obtained in this model is converted to the ballistic coordinate system, there is a component in the direction of velocity, which requires an aircraft Ability to control missile speed. However, in the terminal guidance stage of the missile, the missile usually has no thrust and cannot control the speed of the missile. In response to this problem, the second model is established based on the missile lead angle coordinate system, and the guidance command is perpendicular to the missile speed direction, so the speed is not controlled.

在速度不可控条件下,导弹的三维制导律可分为时间约束和角度约束。时间约束制导律是以导弹的剩余时间估计表达式为基础,或基于一致性理论,设计协同变量进行推导。角度约束制导律则是通过根据模型中视线角加速度与过载的关系,通过滑模控制理论、最优控制理论、自适应控制理论等方法,推导过载指令的表达式。但是目前,在速度不可控条件下,还没有能够同时实现时间约束和角度约束的制导律。Under the condition of uncontrollable speed, the three-dimensional guidance law of missile can be divided into time constraint and angle constraint. The time-constrained guidance law is derived based on the estimation expression of the remaining time of the missile, or based on the consistency theory, by designing cooperative variables. The angle-constrained guidance law is based on the relationship between the line-of-sight angular acceleration and the overload in the model, through sliding mode control theory, optimal control theory, adaptive control theory and other methods to derive the expression of the overload command. But at present, under the condition of uncontrollable speed, there is no guidance law that can realize time constraint and angle constraint at the same time.

在目前能够同时实现时间约束和角度约束的三维制导律中,均是基于控速的三维模型。而在导弹的末制导阶段,导弹通常无推力,仅通过气动力控制导弹姿态,而在速度不可控条件下,目前的制导律仅能实现单一的时间约束,或者角度约束。In the current three-dimensional guidance laws that can realize time constraints and angle constraints at the same time, they are all based on the three-dimensional model of speed control. In the terminal guidance stage of the missile, the missile usually has no thrust, and only controls the attitude of the missile through aerodynamic force. Under the condition of uncontrollable speed, the current guidance law can only achieve a single time constraint or angle constraint.

综上所述,现有方法在速度不可控条件下,无法实现同时对时间和角度进行约束,为了满足制导需求,亟需一种速度不可控条件下时间和角度约束三维滑模制导律设计方法。To sum up, the existing methods cannot constrain time and angle at the same time under the condition of uncontrollable speed. In order to meet the guidance requirements, there is an urgent need for a time- and angle-constrained three-dimensional sliding mode guidance law design method under the condition of uncontrollable speed. .

发明内容SUMMARY OF THE INVENTION

要解决的技术问题technical problem to be solved

针对现有三维制导律在速度不可控条件下,无法同时约束攻击时间和攻击角度的问题,本发明提出了一种速度不可控条件下时间和角度约束三维滑模制导律设计方法。Aiming at the problem that the existing three-dimensional guidance law cannot constrain the attack time and the attack angle at the same time under the condition of uncontrollable speed, the present invention proposes a time- and angle-constrained three-dimensional sliding mode guidance law design method under the condition of uncontrollable speed.

技术方案Technical solutions

本发明对导弹俯仰和偏航通道的制导律分别进行了设计:The present invention designs the guidance laws of the missile pitch and yaw channels respectively:

在俯仰通道,设计了一种快速收敛的非奇异滑模面以及一种随导弹距离变化的趋近律,所推导的制导律能够实现期望视线高低角。In the pitch channel, a fast-converging non-singular sliding mode surface and a reaching law that changes with the missile distance are designed, and the derived guidance law can achieve the desired line of sight angle.

在偏航通道,通过分析导弹制导条件,设计了一种特殊的滑模面,并通过超螺旋滑模理论,保证导弹能够到达滑模面,从而实现了期望视线方位角以及期望攻击时间。In the yaw channel, by analyzing the missile guidance conditions, a special sliding mode surface is designed, and through the super-helix sliding mode theory, it is ensured that the missile can reach the sliding mode surface, thus achieving the desired line of sight azimuth and expected attack time.

一种速度不可控条件下时间和角度约束三维滑模制导律设计方法,其特征在于:在分析导弹制导条件后,基于滑模控制方法,分别设计amy用于控制导弹达到期望视线倾角,amz用于控制导弹达到期望的视线偏角以及期望攻击时间;A time- and angle-constrained three-dimensional sliding mode guidance law design method under the condition of uncontrollable speed, characterized in that: after analyzing the missile guidance conditions, based on the sliding mode control method, respectively design a my to control the missile to reach a desired line of sight inclination, a mz is used to control the missile to reach the desired declination angle and the desired attack time;

导弹制导模型可表示为:The missile guidance model can be expressed as:

Figure BDA0003668940690000021
Figure BDA0003668940690000021

Figure BDA0003668940690000022
Figure BDA0003668940690000022

Figure BDA0003668940690000023
Figure BDA0003668940690000023

Figure BDA0003668940690000024
Figure BDA0003668940690000024

Figure BDA0003668940690000025
Figure BDA0003668940690000025

Figure BDA0003668940690000026
Figure BDA0003668940690000026

其中,R为弹目相对距离,θ为弹道倾角,

Figure BDA0003668940690000027
为弹道偏角,θL为视线高低角,
Figure BDA0003668940690000028
为视线方位角,σm表示前置角,θm
Figure BDA0003668940690000029
为σm在视线坐标系下俯仰和偏航通道的分解,分别表示前置倾角和前置偏角,γm表示前置滚转角,Vm表示导弹速度;Among them, R is the relative distance of the projectile, θ is the ballistic inclination,
Figure BDA0003668940690000027
is the ballistic declination angle, θ L is the line-of-sight angle,
Figure BDA0003668940690000028
is the line-of-sight azimuth, σ m represents the lead angle, θ m and
Figure BDA0003668940690000029
is the decomposition of the pitch and yaw channels of σ m in the line-of-sight coordinate system, representing the forward inclination angle and the forward declination angle, respectively, γ m the forward roll angle, and V m the missile velocity;

对视线高低角和方位角进一步求导可得Further derivation of the line-of-sight angle and azimuth angle can be obtained

Figure BDA0003668940690000031
Figure BDA0003668940690000031

Figure BDA0003668940690000032
Figure BDA0003668940690000032

制导律设计方法如下:The guidance law design method is as follows:

(1)视线高低角方向的制导律设计(1) Guidance law design for the high and low angle directions of the line of sight

首先,为实现期望的视线高低角,设计快速收敛的非奇异滑模面:First, to achieve the desired line-of-sight angle, a fast-converging non-singular sliding mode surface is designed:

Figure BDA0003668940690000033
Figure BDA0003668940690000033

其中,x1=θLLd

Figure BDA0003668940690000034
λ,β均为大于0的常数,α∈[-1,1];设计随着导弹距离变化的自适应趋近律:where x 1LLd ,
Figure BDA0003668940690000034
λ and β are both constants greater than 0, α∈[-1,1]; the adaptive approach law is designed with the change of missile distance:

Figure BDA0003668940690000035
Figure BDA0003668940690000035

其中,k1,k2为大于0的常数,且k1与k2满足:Among them, k 1 and k 2 are constants greater than 0, and k 1 and k 2 satisfy:

Figure BDA0003668940690000036
Figure BDA0003668940690000036

最终可得视线高低角方向制导律为Finally, the guidance law of the line-of-sight direction can be obtained as

Figure BDA0003668940690000037
Figure BDA0003668940690000037

(2)视线方位角方向的制导律设计(2) Guidance law design for line of sight azimuth direction

令ξ=R/Vm,e1=θLLd

Figure BDA0003668940690000039
φ=td-t-ξ,并作出如下假设:Let ξ=R/V m , e 1LLd ,
Figure BDA0003668940690000039
φ=t d -t-ξ, and make the following assumptions:

假设1:在对e2求导时,cosθL这一项变化较小,将其视为一个常数;Assumption 1: When derivation of e 2 , the change of cosθ L is small, and it is regarded as a constant;

假设2:在制导过程中,时间误差项φ≥0恒成立;Assumption 2: During the guidance process, the time error term φ≥0 is always established;

设计如下滑模面:Design the sliding surface as follows:

Figure BDA0003668940690000038
Figure BDA0003668940690000038

对滑模面s2两端求导,并将式代入可得:Taking the derivative of both ends of the sliding surface s 2 and substituting the formula into:

Figure BDA0003668940690000041
Figure BDA0003668940690000041

基于超螺旋滑模控制方法设计趋近律:The reaching law is designed based on the superhelical sliding mode control method:

Figure BDA0003668940690000042
Figure BDA0003668940690000042

Figure BDA0003668940690000043
Figure BDA0003668940690000043

其中,l1和l2均为大于0的常数,w为超滑模的一个自适应参数,防止滑模面出现抖振现象的;将上式代入滑模面s2中,可得控制量amz为:Among them, l 1 and l 2 are both constants greater than 0, and w is an adaptive parameter of the super sliding mode to prevent chattering on the sliding mode surface. Substitute the above formula into the sliding mode surface s 2 to obtain the control variable a mz is:

Figure BDA0003668940690000044
Figure BDA0003668940690000044

Figure BDA0003668940690000045
Figure BDA0003668940690000045

一种计算机系统,其特征在于包括:一个或多个处理器,计算机可读存储介质,用于存储一个或多个程序,其中,当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现上述的方法。A computer system, characterized by comprising: one or more processors, and a computer-readable storage medium for storing one or more programs, wherein when the one or more programs are processed by the one or more programs When the processor is executed, the one or more processors are caused to implement the above method.

一种计算机可读存储介质,其特征在于存储有计算机可执行指令,所述指令在被执行时用于实现上述的方法。A computer-readable storage medium is characterized in that computer-executable instructions are stored, and the instructions, when executed, are used to implement the above-mentioned method.

一种计算机程序,其特征在于包括计算机可执行指令,所述指令在被执行时用于实现上述的方法。A computer program characterized by comprising computer-executable instructions which, when executed, are used to implement the above-mentioned method.

有益效果beneficial effect

本发明设计的制导律垂直于速度方向,仅通过改变导弹的速度方向实现期望的攻击时间和攻击角度。在俯仰通道,提出了一种快速收敛的非奇异滑模面以及一种随导弹距离变化的滑模趋近律,实现了期望视线高低角;在偏航通道,通过分析其制导条件,设计了一种特殊的滑模面,并通过超螺旋滑模理论,保证导弹能够到达滑模面,从而实现期望视线方位角以及期望攻击时间。The guidance law designed by the invention is perpendicular to the speed direction, and the desired attack time and attack angle can be achieved only by changing the speed direction of the missile. In the pitch channel, a fast-converging non-singular sliding mode surface and a sliding mode approach law that changes with the missile distance are proposed to achieve the desired line of sight angle; in the yaw channel, by analyzing the guidance conditions, a A special sliding mode surface, and through the super-helix sliding mode theory, ensures that the missile can reach the sliding mode surface, so as to achieve the desired line of sight azimuth and the desired attack time.

本发明在导弹速度不可控条件下,实现了在三维空间中导弹以期望攻击时间和攻击角度打击目标,该制导律脱靶量小,能够准确达到期望攻击时间和期望视线高低角,视线方位角仅存在较小误差。Under the condition that the speed of the missile is uncontrollable, the present invention realizes that the missile hits the target with the desired attack time and attack angle in the three-dimensional space. There are minor errors.

本发明的创新点如下:The innovation of the present invention is as follows:

(1)通过滑模控制理论,在速度不可控条件下,根据制导条件设计了滑模面,以此推导出了能够同时以期望时间和角度打击目标的制导律,这是其他方法并未实现的。(1) Through the sliding mode control theory, under the condition of uncontrollable speed, the sliding mode surface is designed according to the guidance conditions, and the guidance law that can hit the target at the desired time and angle at the same time is deduced, which is not achieved by other methods. of.

(2)本发明对导弹的三维制导律进行了合理分工,其中amy保证导弹的视线高低角达到期望值,amz保证导弹的视线方位角和攻击时间达到期望值。(2) The present invention has carried out a reasonable division of labor for the three-dimensional guidance law of the missile, wherein a my ensures that the height angle of the missile's line of sight reaches the desired value, and a mz ensures that the azimuth of the missile's line of sight and attack time reach the desired value.

(3)本发明的所有结果是基于三维非线性耦合动力学模型推导的,并未考虑小角度假设以及线性化,制导律的精度高,适用范围广,能够实现全方位的攻击。(3) All the results of the present invention are derived based on the three-dimensional nonlinear coupled dynamics model, without considering the small angle assumption and linearization, the guidance law has high precision, wide application range, and can achieve all-round attack.

附图说明Description of drawings

附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制,在整个附图中,相同的参考符号表示相同的部件。The drawings are for the purpose of illustrating specific embodiments only and are not to be considered limiting of the invention, and like reference numerals refer to like parts throughout the drawings.

图1:三维制导模型示意图;Figure 1: Schematic diagram of the 3D guidance model;

图2:导弹三维轨迹图;Figure 2: Three-dimensional trajectory map of the missile;

图3:导弹估计剩余时间曲线;Figure 3: Missile estimated remaining time curve;

图4:前置角变化曲线;Figure 4: lead angle change curve;

图5:相对速度变化曲线;Figure 5: Relative velocity curve;

图6:导弹视线高低角变化曲线;Figure 6: The change curve of the height and low angle of the missile line of sight;

图7:导弹视线方位角变化曲线;Figure 7: The change curve of the azimuth angle of the missile line of sight;

图8:视线高低角方向过载变化曲线;Figure 8: Variation curve of overload in the high and low angle directions of the line of sight;

图9:视线方位角方向过载变化曲线。Figure 9: Overload variation curve in line-of-sight azimuth direction.

具体实施方式Detailed ways

为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。此外,下面描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。In order to make the objectives, technical solutions and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention, but not to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

根据图1,M,T分别表示导弹与目标的位置,在Y轴朝上的坐标系下,导弹的三维耦合制导模型可描述为:According to Figure 1, M and T represent the positions of the missile and the target respectively. In the coordinate system with the Y axis pointing upward, the three-dimensional coupled guidance model of the missile can be described as:

Figure BDA0003668940690000061
Figure BDA0003668940690000061

Figure BDA0003668940690000062
Figure BDA0003668940690000062

Figure BDA0003668940690000063
Figure BDA0003668940690000063

Figure BDA0003668940690000064
Figure BDA0003668940690000064

Figure BDA0003668940690000065
Figure BDA0003668940690000065

Figure BDA0003668940690000066
Figure BDA0003668940690000066

其中,R为弹目相对距离,θ为弹道倾角,

Figure BDA0003668940690000067
为弹道偏角,θL为视线高低角,
Figure BDA0003668940690000068
为视线方位角,σm表示前置角,θm
Figure BDA0003668940690000069
为σm在视线坐标系下俯仰和偏航通道的分解,分别表示前置倾角和前置偏角,γm表示前置滚转角,Vm表示导弹速度。Among them, R is the relative distance of the projectile, θ is the ballistic inclination,
Figure BDA0003668940690000067
is the ballistic declination angle, θ L is the line-of-sight angle,
Figure BDA0003668940690000068
is the line-of-sight azimuth, σ m represents the lead angle, θ m and
Figure BDA0003668940690000069
is the decomposition of the pitch and yaw channels of σ m in the line-of-sight coordinate system, representing the forward inclination angle and the forward declination angle, respectively, γ m the forward roll angle, and V m the missile velocity.

(1)视线高低角方向的制导律设计(1) Guidance law design for the high and low angle directions of the line of sight

首先,对(5)求导,并将式(2)代入,可得First, take the derivative of (5) and substitute equation (2), we can get

Figure BDA00036689406900000610
Figure BDA00036689406900000610

为实现期望的视线高低角,设计了一种快速收敛的非奇异滑模面:To achieve the desired line-of-sight angle, a fast-converging non-singular sliding surface is designed:

s1=x2+λx1+βsigα(x1) (8)s 1 =x 2 +λx 1 +βsig α (x 1 ) (8)

其中,x1=θLLd

Figure BDA0003668940690000071
θLd为期望视线高低角,λ,β均为大于0的常数,α∈[-1,1]。为了使视线高低角能快速收敛到期望值且保持稳定,本发明提出了一种随着导弹距离变化的趋近律:where x 1LLd ,
Figure BDA0003668940690000071
θ Ld is the desired line-of-sight angle, λ and β are both constants greater than 0, α∈[-1,1]. In order to make the line-of-sight angle quickly converge to the desired value and keep it stable, the present invention proposes a reaching law that changes with the missile distance:

Figure BDA0003668940690000072
Figure BDA0003668940690000072

其中,k1,k2为大于0的常数,且k1与k2满足:Among them, k 1 and k 2 are constants greater than 0, and k 1 and k 2 satisfy:

Figure BDA0003668940690000073
Figure BDA0003668940690000073

对式(8)求导,结合式(7)、式(9),可得俯仰通道制导律为Taking the derivation of Equation (8) and combining Equation (7) and Equation (9), the pitch channel guidance law can be obtained as

Figure BDA0003668940690000074
Figure BDA0003668940690000074

为证明在该制导律下,导弹能实现期望的视线高低角,假设如下李雅普诺夫函数:In order to prove that the missile can achieve the desired line-of-sight angle under this guidance law, the following Lyapunov function is assumed:

Figure BDA0003668940690000075
Figure BDA0003668940690000075

对该李雅普诺夫函数求导,可得:Taking the derivative of this Lyapunov function, we get:

Figure BDA0003668940690000076
Figure BDA0003668940690000076

对上式化简可得:Simplify the above formula to get:

Figure BDA0003668940690000077
Figure BDA0003668940690000077

对上式积分可得:Integrating the above formula can get:

Figure BDA0003668940690000078
Figure BDA0003668940690000078

其中,s1(0)和R(0)分别为开始时刻的滑模面以及导弹与目标的相对距离。根据(15)可知,存在一个有限时间t*使得滑模面s1=0。Among them, s 1 (0) and R(0) are the sliding surface at the start time and the relative distance between the missile and the target, respectively. According to (15), there is a finite time t * such that the sliding surface s 1 =0.

在到达滑模面后,根据式(8)可得如下关系式After reaching the sliding mode surface, the following relation can be obtained according to formula (8)

x2=-λx1-βsigα(x1) (16)x 2 =-λx 1 -βsig α (x 1 ) (16)

为证明导弹的视线高低角能达到期望值,设计李雅普诺夫函数为:In order to prove that the line-of-sight angle of the missile can reach the expected value, the Lyapunov function is designed as:

Figure BDA0003668940690000081
Figure BDA0003668940690000081

对式(17)求导并将式(16)代入可得:Taking the derivative of equation (17) and substituting it into equation (16), we get:

Figure BDA0003668940690000082
Figure BDA0003668940690000082

根据式(18)可知,x1能在有限时间内收敛到0,收敛时间T1满足According to equation (18), x 1 can converge to 0 in a finite time, and the convergence time T 1 satisfies

Figure BDA0003668940690000083
Figure BDA0003668940690000083

综上所述,视线高低角达到期望值的总时间为到达滑模面的时间t*与到达滑模面后x1的收敛时间T1之和,Tθ=t*+T1To sum up, the total time for the line of sight height angle to reach the desired value is the sum of the time t * reaching the sliding mode surface and the convergence time T 1 of x 1 after reaching the sliding mode surface, Tθ=t * +T 1 .

(2)视线方位角方向的制导律设计(2) Guidance law design for line of sight azimuth direction

令ξ=R/Vm,表示估计剩余飞行时间,e1=θLLd,表示视线高低角误差,

Figure BDA0003668940690000084
表示视线方位角误差,其中
Figure BDA0003668940690000085
为期望视线方位角,φ=td-t-ξ。表示期望攻击时间误差。下面作如下假设:Let ξ=R/V m , denote the estimated remaining flight time, e 1LLd , denote the line-of-sight angle error,
Figure BDA0003668940690000084
is the line-of-sight azimuth error, where
Figure BDA0003668940690000085
is the desired line-of-sight azimuth, φ=t d -t-ξ. Indicates the expected attack time error. The following assumptions are made:

假设1:在对e2求导时,cosθL这一项变化较小,将其视为一个常数。Assumption 1: When differentiating e2 , the term cosθL changes little and treats it as a constant.

假设2:在制导过程中,时间误差项φ≥0恒成立。Assumption 2: During the guidance process, the time error term φ≥0 is always established.

对上式中的ξ,e1,e2,φ求导,并将式(1)-(6)代入,可得:Taking derivatives of ξ, e 1 , e 2 , and φ in the above formula, and substituting formulas (1)-(6), we can get:

Figure BDA0003668940690000086
Figure BDA0003668940690000086

Figure BDA0003668940690000087
Figure BDA0003668940690000087

Figure BDA0003668940690000088
Figure BDA0003668940690000088

Figure BDA0003668940690000089
Figure BDA0003668940690000089

设计一个辅助变量Ω,其表达式为:Design an auxiliary variable Ω whose expression is:

Figure BDA0003668940690000091
Figure BDA0003668940690000091

对Ω两端求导并将式(20)-(23)代入,整理可得:Taking the derivative of both ends of Ω and substituting equations (20)-(23), we can get:

Figure BDA0003668940690000092
Figure BDA0003668940690000092

当amy控制导弹的视线高低角达到期望值后,根据式(5),将满足如下条件:When the line of sight of a my control missile reaches the desired value, according to formula (5), the following conditions will be satisfied:

Figure BDA0003668940690000093
Figure BDA0003668940690000093

此时,式(25)可以改写为In this case, equation (25) can be rewritten as

Figure BDA0003668940690000094
Figure BDA0003668940690000094

根据(27),设计如下滑模面:According to (27), the sliding surface is designed as follows:

Figure BDA0003668940690000095
Figure BDA0003668940690000095

对滑模面s2两端求导,并将式(3)代入可得:Taking the derivation of both ends of the sliding surface s 2 and substituting Equation (3) into:

Figure BDA0003668940690000096
Figure BDA0003668940690000096

为削弱滑模面的抖振现象,并保证滑模面能够快速收敛,本发明基于超螺旋滑模控制方法设计了趋近律,首先将式(29)写为:In order to weaken the chattering phenomenon of the sliding mode surface and ensure that the sliding mode surface can converge quickly, the present invention designs a reaching law based on the superhelical sliding mode control method. First, formula (29) is written as:

Figure BDA0003668940690000097
Figure BDA0003668940690000097

其中,in,

Figure BDA0003668940690000098
Figure BDA0003668940690000098

Figure BDA0003668940690000099
Figure BDA0003668940690000099

根据超螺旋滑模理论,趋近律可表示为:According to the superhelical sliding mode theory, the reaching law can be expressed as:

Figure BDA0003668940690000101
Figure BDA0003668940690000101

其中,l1和l2均为大于0的常数。将式(33)代入到式(30)中,可得控制量amz为:Among them, l 1 and l 2 are both constants greater than 0. Substituting Equation (33) into Equation (30), the control amount a mz can be obtained as:

Figure BDA0003668940690000102
Figure BDA0003668940690000102

首先,根据超螺旋滑模理论可知,滑模面能够在有限时间内达到s2=0。在达到s2=0后,假设此时amy也控制导弹达到期望的视线高低角则根据式(24)和(25)有:First, according to the super-helical sliding mode theory, the sliding mode surface can reach s 2 =0 in a limited time. After reaching s 2 =0, assuming that a my also controls the missile to reach the desired line-of-sight angle, according to equations (24) and (25):

Figure BDA0003668940690000103
Figure BDA0003668940690000103

由于cosθm=1,根据三角公式,式(20)和式(23)可以表示为:Since cosθ m =1, according to the trigonometric formula, equations (20) and (23) can be expressed as:

Figure BDA0003668940690000104
Figure BDA0003668940690000104

Figure BDA0003668940690000105
Figure BDA0003668940690000105

首先,设计如下Lyapunov函数:First, design the following Lyapunov function:

V3=0.5φ2 (38)V 3 =0.5φ 2 (38)

对其两端求导,有Taking the derivative of both ends, we have

Figure BDA0003668940690000106
Figure BDA0003668940690000106

由于φ≥0,因此,

Figure BDA0003668940690000107
当且仅当φ=0时,
Figure BDA0003668940690000108
因此φ能够渐进收敛到0。根据假设2,在达到期望时间之前,φ≥0恒成立,假设φ=0,根据上式可得:Since φ≥0, therefore,
Figure BDA0003668940690000107
If and only if φ=0,
Figure BDA0003668940690000108
Therefore, φ can gradually converge to 0. According to Assumption 2, before reaching the desired time, φ≥0 is always established. Assuming that φ=0, according to the above formula, we can get:

Figure BDA0003668940690000109
Figure BDA0003668940690000109

t+ξ=td (41)t+ξ=t d (41)

若导弹能够击中目标,则在击中目标时,根据定义可知,+ξ=R/Vm=0,根据式(41)可知,此时满足,此时满足t=td。综上所述,只要导弹能够击中目标,则此时导弹是在期望时间击中目标的,且能实现期望的视线方位角。If the missile can hit the target, when hitting the target, according to the definition, +ξ=R/V m =0, according to formula (41), it is satisfied at this time, and t=t d is satisfied at this time. To sum up, as long as the missile can hit the target, then the missile hits the target at the desired time and can achieve the desired line of sight azimuth.

结合amy和amz可知,导弹能够实现如下制导目标:Combining a my and a mz , the missile can achieve the following guidance targets:

Figure BDA0003668940690000111
Figure BDA0003668940690000111

即导弹能够在期望时间击中目标,且能实现期望的视线倾角和视线偏角。That is, the missile can hit the target at the desired time, and can achieve the desired line of sight inclination and line of sight declination.

为论证速度不可控条件下时间和角度约束三维滑模制导律设计方法的可行性和有效性,本发明设计了仿真验证试验。导弹的初始条件如表1所示,四枚导弹的视线高低角约束为θLd=[-35° -10° -15° -20°],视线方位角约束为

Figure BDA0003668940690000112
导弹的期望攻击时间Td均为60s,仿真终止条件设置为r<0.5m,导弹的最大过载为30g,仿真步长为0.001s。仿真结果如表2和图2-图9所示In order to demonstrate the feasibility and effectiveness of the time and angle constrained three-dimensional sliding mode guidance law design method under the condition of uncontrollable speed, the present invention designs a simulation verification test. The initial conditions of the missiles are shown in Table 1. The line-of-sight angle constraints of the four missiles are θ Ld = [-35° -10° -15° -20°], and the line-of-sight azimuth is constrained by
Figure BDA0003668940690000112
The expected attack time T d of the missile is 60s, the simulation termination condition is set as r<0.5m, the maximum overload of the missile is 30g, and the simulation step size is 0.001s. The simulation results are shown in Table 2 and Figure 2-Figure 9

表1导弹初始条件Table 1 Missile initial conditions

Figure BDA0003668940690000113
Figure BDA0003668940690000113

表2仿真结果Table 2 Simulation results

Figure BDA0003668940690000114
Figure BDA0003668940690000114

以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明公开的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。The above are only specific embodiments of the present invention, but the protection scope of the present invention is not limited to this. Any person skilled in the art can easily think of various equivalents within the technical scope disclosed by the present invention. Modifications or substitutions should be included within the protection scope of the present invention.

Claims (4)

1.一种速度不可控条件下时间和角度约束三维滑模制导律设计方法,其特征在于:在分析导弹制导条件后,基于滑模控制方法,分别设计amy用于控制导弹达到期望视线倾角,amz用于控制导弹达到期望的视线偏角以及期望攻击时间;1. a time and angle constraint three-dimensional sliding mode guidance law design method under uncontrollable speed conditions, is characterized in that: after analyzing the missile guidance conditions, based on the sliding mode control method, respectively design amy for controlling the missile to reach the desired line of sight inclination , a mz is used to control the missile to reach the desired declination angle and the desired attack time; 导弹制导模型可表示为:The missile guidance model can be expressed as:
Figure FDA0003668940680000011
Figure FDA0003668940680000011
Figure FDA0003668940680000012
Figure FDA0003668940680000012
Figure FDA0003668940680000013
Figure FDA0003668940680000013
Figure FDA0003668940680000014
Figure FDA0003668940680000014
Figure FDA0003668940680000015
Figure FDA0003668940680000015
Figure FDA0003668940680000016
Figure FDA0003668940680000016
其中,R为弹目相对距离,θ为弹道倾角,
Figure FDA0003668940680000017
为弹道偏角,θL为视线高低角,
Figure FDA0003668940680000018
为视线方位角,σm表示前置角,θm
Figure FDA0003668940680000019
为σm在视线坐标系下俯仰和偏航通道的分解,分别表示前置倾角和前置偏角,γm表示前置滚转角,Vm表示导弹速度;
Among them, R is the relative distance of the projectile, θ is the ballistic inclination,
Figure FDA0003668940680000017
is the ballistic declination angle, θ L is the line of sight height angle,
Figure FDA0003668940680000018
is the line-of-sight azimuth, σ m represents the lead angle, θ m and
Figure FDA0003668940680000019
is the decomposition of the pitch and yaw channels of σ m in the line-of-sight coordinate system, representing the forward inclination angle and the forward declination angle, respectively, γ m the forward roll angle, and V m the missile velocity;
对视线高低角和方位角进一步求导可得Further derivation of the line-of-sight angle and azimuth angle can be obtained
Figure FDA00036689406800000110
Figure FDA00036689406800000110
Figure FDA00036689406800000111
Figure FDA00036689406800000111
制导律设计方法如下:The guidance law design method is as follows: (1)视线高低角方向的制导律设计(1) Guidance law design for the high and low angle directions of the line of sight 首先,为实现期望的视线高低角,设计快速收敛的非奇异滑模面:First, to achieve the desired line-of-sight angle, a fast-converging non-singular sliding mode surface is designed:
Figure FDA00036689406800000112
Figure FDA00036689406800000112
其中,x1=θLLd
Figure FDA00036689406800000113
λ,β均为大于0的常数,α∈[-1,1];设计随着导弹距离变化的自适应趋近律:
where x 1LLd ,
Figure FDA00036689406800000113
λ and β are both constants greater than 0, α∈[-1,1]; the adaptive approach law is designed with the change of missile distance:
Figure FDA0003668940680000021
Figure FDA0003668940680000021
其中,k1,k2为大于0的常数,且k1与k2满足:Among them, k 1 and k 2 are constants greater than 0, and k 1 and k 2 satisfy:
Figure FDA0003668940680000022
Figure FDA0003668940680000022
最终可得视线高低角方向制导律为Finally, the guidance law of the line-of-sight direction can be obtained as
Figure FDA0003668940680000023
Figure FDA0003668940680000023
(2)视线方位角方向的制导律设计(2) Guidance law design for line of sight azimuth direction 令ξ=R/Vm,e1=θLLd
Figure FDA0003668940680000024
φ=td-t-ξ,并作出如下假设:
Let ξ=R/V m , e 1LLd ,
Figure FDA0003668940680000024
φ=t d -t-ξ, and make the following assumptions:
假设1:在对e2求导时,cosθL这一项变化较小,将其视为一个常数;Assumption 1: When derivation of e 2 , the change of cosθ L is small, and it is regarded as a constant; 假设2:在制导过程中,时间误差项φ≥0恒成立;Assumption 2: During the guidance process, the time error term φ≥0 is always established; 设计如下滑模面:Design the sliding surface as follows:
Figure FDA0003668940680000025
Figure FDA0003668940680000025
对滑模面s2两端求导,并将式代入可得:Taking the derivative of both ends of the sliding surface s 2 and substituting the formula into:
Figure FDA0003668940680000026
Figure FDA0003668940680000026
基于超螺旋滑模控制方法设计趋近律:The reaching law is designed based on the superhelical sliding mode control method:
Figure FDA0003668940680000027
Figure FDA0003668940680000027
Figure FDA0003668940680000028
Figure FDA0003668940680000028
其中,l1和l2均为大于0的常数,w为超滑模的一个自适应参数,防止滑模面出现抖振现象的;将上式代入滑模面s2中,可得控制量amz为:Among them, l 1 and l 2 are both constants greater than 0, and w is an adaptive parameter of the super sliding mode to prevent chattering on the sliding mode surface. Substitute the above formula into the sliding mode surface s 2 to obtain the control variable a mz is:
Figure FDA0003668940680000031
Figure FDA0003668940680000031
Figure FDA0003668940680000032
Figure FDA0003668940680000032
2.一种计算机系统,其特征在于包括:一个或多个处理器,计算机可读存储介质,用于存储一个或多个程序,其中,当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现权利要求1所述的方法。2. A computer system, characterized by comprising: one or more processors, and a computer-readable storage medium for storing one or more programs, wherein when the one or more programs are executed by the one or more programs When executed by a plurality of processors, the one or more processors are caused to implement the method of claim 1 . 3.一种计算机可读存储介质,其特征在于存储有计算机可执行指令,所述指令在被执行时用于实现权利要求1所述的方法。3. A computer-readable storage medium, characterized by storing computer-executable instructions that, when executed, are used to implement the method of claim 1 . 4.一种计算机程序,其特征在于包括计算机可执行指令,所述指令在被执行时用于实现权利要求1所述的方法。4. A computer program characterized by comprising computer-executable instructions which, when executed, are used to implement the method of claim 1 .
CN202210600525.3A 2022-05-30 2022-05-30 Design method of three-dimensional sliding mode guidance law with time and angle constraints under uncontrollable speed Active CN114995495B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210600525.3A CN114995495B (en) 2022-05-30 2022-05-30 Design method of three-dimensional sliding mode guidance law with time and angle constraints under uncontrollable speed

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210600525.3A CN114995495B (en) 2022-05-30 2022-05-30 Design method of three-dimensional sliding mode guidance law with time and angle constraints under uncontrollable speed

Publications (2)

Publication Number Publication Date
CN114995495A true CN114995495A (en) 2022-09-02
CN114995495B CN114995495B (en) 2024-11-05

Family

ID=83032093

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210600525.3A Active CN114995495B (en) 2022-05-30 2022-05-30 Design method of three-dimensional sliding mode guidance law with time and angle constraints under uncontrollable speed

Country Status (1)

Country Link
CN (1) CN114995495B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116663263A (en) * 2023-05-16 2023-08-29 哈尔滨工业大学 A sliding mode surface-based method and system for constrained guidance of aircraft terminal falling speed and falling angle
CN117891271A (en) * 2024-03-18 2024-04-16 西北工业大学 Three-dimensional cooperative guidance method for high-speed aircraft considering time and angle constraints

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104266546A (en) * 2014-09-22 2015-01-07 哈尔滨工业大学 A Line-of-Sight Based Convergent Active Defense Guidance Control Method
DE102016004363A1 (en) * 2015-08-27 2017-03-02 Diehl Bgt Defence Gmbh & Co. Kg Flight control device for a missile and method for steering a missile
CN106843265A (en) * 2016-12-30 2017-06-13 哈尔滨工业大学 Three-dimensional many guided missile cooperative guidance method and systems of finite time convergence control
CN106934120A (en) * 2017-02-23 2017-07-07 哈尔滨工业大学 Three-dimensional guidance based on preceding from the interception hypersonic aircraft to guidance restrains method for designing
CN108362171A (en) * 2018-02-02 2018-08-03 北京航空航天大学 A kind of Guidance constrained with attack time and angle-of-attack
US20200141699A1 (en) * 2018-11-02 2020-05-07 Bae Systems Information And Electronic Systems Integration Inc. Multi-mode adaptive nonlinear trajectory shaping (nts) guidance law

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104266546A (en) * 2014-09-22 2015-01-07 哈尔滨工业大学 A Line-of-Sight Based Convergent Active Defense Guidance Control Method
DE102016004363A1 (en) * 2015-08-27 2017-03-02 Diehl Bgt Defence Gmbh & Co. Kg Flight control device for a missile and method for steering a missile
CN106843265A (en) * 2016-12-30 2017-06-13 哈尔滨工业大学 Three-dimensional many guided missile cooperative guidance method and systems of finite time convergence control
CN106934120A (en) * 2017-02-23 2017-07-07 哈尔滨工业大学 Three-dimensional guidance based on preceding from the interception hypersonic aircraft to guidance restrains method for designing
CN108362171A (en) * 2018-02-02 2018-08-03 北京航空航天大学 A kind of Guidance constrained with attack time and angle-of-attack
US20200141699A1 (en) * 2018-11-02 2020-05-07 Bae Systems Information And Electronic Systems Integration Inc. Multi-mode adaptive nonlinear trajectory shaping (nts) guidance law

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴放;常思江;李新华;: "考虑导弹速度变化的攻击时间和攻击角度控制滑模制导律", 弹道学报, no. 04, 15 December 2019 (2019-12-15) *
杨芳;张宽桥;余磊;: "自适应非奇异快速终端二阶滑模制导律", 弹道学报, no. 02, 15 June 2020 (2020-06-15) *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116663263A (en) * 2023-05-16 2023-08-29 哈尔滨工业大学 A sliding mode surface-based method and system for constrained guidance of aircraft terminal falling speed and falling angle
CN117891271A (en) * 2024-03-18 2024-04-16 西北工业大学 Three-dimensional cooperative guidance method for high-speed aircraft considering time and angle constraints
CN117891271B (en) * 2024-03-18 2024-05-31 西北工业大学 Three-dimensional collaborative guidance method for high-speed aircraft in consideration of time and angle constraints

Also Published As

Publication number Publication date
CN114995495B (en) 2024-11-05

Similar Documents

Publication Publication Date Title
CN106681348B (en) Guidance and control integrated design method considering full strapdown seeker field of view constraints
CN110008502B (en) An integrated design method for 3D guidance and control considering the full strapdown seeker field of view constraints
CN103090728B (en) Tail angle restraining guidance method based on sliding mode control
CN112782984B (en) Multi-aircraft anti-interception cooperative attack guidance method and device and aircraft
CN111551080B (en) A Guidance Method for Attack Time Control
CN105184109B (en) Disturb trajectory motors in boost phase penetration state deviation analytic method under graviational interaction
CN109557814B (en) A Finite Time Integral Sliding Mode Terminal Guidance Law
CN113642122B (en) Method and system for acquiring elements for remote interception and launch based on single-sided shooting table
CN114995495A (en) Time and angle constraint three-dimensional sliding mode guidance law design method under speed uncontrollable condition
CN115406312B (en) An Integrated Method of Missile Guidance and Control Considering Field of View and Steering Gear Delay Constraints
CN112305919A (en) A Design Method for Fixed-Time Sliding Mode Guidance Law with Collision Angle Constraint
CN109084641B (en) Missile guidance method and device
CN111473696A (en) Guidance rocket vertical attack guidance method based on drop point estimation
CN111591470A (en) A Closed-loop Guidance Method for Accurate Soft Landing of Aircraft Adapting to Adjustable Thrust Mode
CN112013726B (en) An integrated design method for fully strapdown guidance and control based on a third-order model
CN111336871A (en) A vertical attack guidance method based on detour flight
CN112069605A (en) A Design Method of Proportional Guidance Law with Attack Time Constraint
CN115033024A (en) Multi-missile three-dimensional cooperative guidance method based on attack time and angle common constraints
Zhao et al. New three-dimensional second-order sliding mode guidance law with impact-angle constraints
Yang Analysis of optimal midcourse guidance law
CN114815888A (en) Affine form guidance control integrated control method
CN111176315B (en) Variant cruise missile control method based on L1 adaptive control
Wang et al. A guidance and control design with reduced information for a dual-spin stabilized projectile
CN113341710B (en) A compound control method and application of aircraft agile turning
CN114610057A (en) Design method for maneuvering penetration strategy of high-Mach aircraft

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