CN113110576B - Self-adaptive fixed time convergence cooperative guidance method capable of realizing continuous switching - Google Patents

Self-adaptive fixed time convergence cooperative guidance method capable of realizing continuous switching Download PDF

Info

Publication number
CN113110576B
CN113110576B CN202110402142.0A CN202110402142A CN113110576B CN 113110576 B CN113110576 B CN 113110576B CN 202110402142 A CN202110402142 A CN 202110402142A CN 113110576 B CN113110576 B CN 113110576B
Authority
CN
China
Prior art keywords
missile
missiles
target
formula
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
Application number
CN202110402142.0A
Other languages
Chinese (zh)
Other versions
CN113110576A (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
Shanghai Institute of Electromechanical Engineering
Original Assignee
Northwestern Polytechnical University
Shanghai Institute of Electromechanical Engineering
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, Shanghai Institute of Electromechanical Engineering filed Critical Northwestern Polytechnical University
Priority to CN202110402142.0A priority Critical patent/CN113110576B/en
Publication of CN113110576A publication Critical patent/CN113110576A/en
Application granted granted Critical
Publication of CN113110576B publication Critical patent/CN113110576B/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

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)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

The invention discloses a continuously switchable self-adaptive fixed time convergence cooperative guidance method, which comprises the steps of firstly establishing a motion relation mathematical model of a missile relative to a target, predicting the residual hit time of each missile on the basis, further defining the consistency error variable of the residual hit time, and finally providing a continuously switchable self-adaptive fixed time convergence cooperative guidance law capable of ensuring that a plurality of missiles hit the target simultaneously. The method can enable a plurality of missiles to finally hit the target at the same time, and compared with the existing fixed-time convergence cooperative guidance method, the method designed by the invention can effectively reduce the impact of the initial guidance instruction.

Description

Self-adaptive fixed time convergence cooperative guidance method capable of realizing continuous switching
Technical Field
The invention belongs to the technical field of guidance, and particularly relates to a collaborative guidance method.
Background
With the continuous progress and improvement of the counterguidance defense system, the traditional one-to-one guidance attack tactics have great difficulty in effectively striking the target equipped with the defense system, and the ideal combat efficiency is difficult to achieve. The high-density simultaneous hit saturation attack tactics of multiple aircrafts become an effective method for hitting the targets due to the characteristics of strong penetration capacity and high hit rate.
The convergence time obviously influences the execution effect of the cooperative guidance algorithm, and the finite time convergence method is widely applied to the design of the cooperative guidance algorithm due to the advantages of fast convergence time, estimable convergence time upper bound and strong robustness. However, the upper bound of the convergence time of the finite time convergence method has strong dependence on the initial state of the system, and it is difficult to estimate the upper bound of the convergence time when the initial state of the system is unknown. Inspired by this problem, researchers have proposed a fixed time convergence method in which the upper bound of the convergence time does not depend on the initial state. Therefore, the research on the cooperative guidance method capable of realizing fixed time convergence is of great significance, and the initial control quantity of the existing fixed time convergence method is considered to be large in impact and is not beneficial to being directly applied to the cooperative guidance algorithm.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a continuously switchable self-adaptive fixed time convergence cooperative guidance method, which comprises the steps of firstly establishing a motion relation mathematical model of a missile relative to a target, predicting the residual hit time of each missile on the basis, further defining the consistency error variable of the residual hit time, and finally providing a continuously switchable self-adaptive fixed time convergence cooperative guidance law capable of ensuring that a plurality of missiles hit the target simultaneously. The method can enable a plurality of missiles to finally hit the target at the same time, and compared with the existing fixed-time convergence cooperative guidance method, the method designed by the invention can effectively reduce the impact of the initial guidance instruction.
The technical scheme adopted by the invention for solving the technical problem comprises the following steps:
step 1: establishing a missile relative target motion relation mathematical model;
assuming that n missiles strike a static target, the communication relation between every two missiles adopts an adjacency matrix A ═ a in graph theory ij ]It shows that if the i (i is 1,2, …, n) th missile can be built with the j (j is 1,2, …, n) th missileCommunication relation rule a ij 1, otherwise, a ij =0;
Aiming at the striking of a plurality of missiles on a fixed target under a two-dimensional plane, the communication topology among the missiles is assumed to be an undirected connected graph, and a relative kinematic model of the ith missile and the target is described by a formula (1):
Figure BDA0003020799260000021
in the formula, r i Indicates the distance, V, between the ith missile and the target i Representing constant flight speed of the ith missile, q i Is the i-th missile line-of-sight angle, gamma i And phi i Respectively is a track angle and a lead angle of the ith missile; a is i The normal acceleration of the ith missile is perpendicular to the speed direction;
step 2: predicting the residual hit time of the missile;
residual hit time of ith missile
Figure BDA0003020799260000022
The following formula was used for prediction:
Figure BDA0003020799260000023
in the formula N s >2 denotes the navigation ratio;
and 3, step 3: variable xi defining residual hit time consistency error i
Figure BDA0003020799260000024
In the formula (I), the compound is shown in the specification,
Figure BDA0003020799260000025
the residual hit time of the jth missile is obtained;
and 4, step 4: designing continuous switching self-adaptive fixed time convergence cooperative guidance law a i
Figure BDA0003020799260000026
In the formula, k i For guidance gain, 0<μ 1 <1,μ 2 >1,s i =(sign(ξ i +δ)-sign(ξ i -δ))/2,δ>0 denotes the switching boundary, ∈ i Is a positive real number; definition of
Figure BDA0003020799260000027
Is phi i The upper bound of (a) is,
Figure BDA0003020799260000028
is composed of
Figure BDA0003020799260000029
Is determined by the estimated value of (c),
Figure BDA00030207992600000210
the values of (a) are determined by the adaptation law:
Figure BDA00030207992600000211
wherein gamma is d And σ is a positive real number.
Preferably, k is i =1.3,μ 1 =0.5,μ 2 =1.5,γ d =2.7,N s =3,δ=1.2,
Figure BDA00030207992600000212
ε i =9,σ=0.1。
The invention has the following beneficial effects:
the invention designs a continuously switchable self-adaptive fixed time convergence cooperative guidance method, and multiple missiles can finally hit a target at the same time under a designed cooperative guidance law. Compared with the existing fixed-time convergence cooperative guidance method, the method designed by the invention can effectively reduce the impact of the initial guidance instruction.
Drawings
FIG. 1 is a schematic diagram of a multi-missile communication topology of the method of the invention.
FIG. 2 is a conventional fixed time convergence collaborative guided flight trajectory plot.
FIG. 3 is a conventional fixed time convergence coordinated guidance projectile distance curve.
FIG. 4 is a conventional fixed time convergent cooperative guidance consistency error curve.
FIG. 5 is a conventional fixed time convergence coordinated guided normal acceleration curve.
FIG. 6 is a continuously switched adaptive fixed time convergence cooperative guidance flight trajectory curve of the present invention.
FIG. 7 is a continuously switched adaptive fixed time convergence cooperative guided missile distance curve of the present invention.
FIG. 8 is a continuously switched adaptive fixed time convergence cooperative guidance consistency error curve of the present invention.
FIG. 9 is a continuously switched adaptive fixed time convergence cooperative guided normal acceleration curve of the present invention.
Detailed Description
The invention is further illustrated with reference to the following figures and examples.
The invention provides a continuously switched self-adaptive fixed time convergence cooperative guidance method aiming at the problem that multiple missiles hit a static target at the same time. The cooperative guidance method has the characteristic of fixed time convergence, meanwhile, the impact of an initial guidance instruction can be effectively reduced through the continuous switching characteristic, and the gain in the cooperative guidance law realizes real-time self-adjustment through a self-adaptive law.
A continuously switchable self-adaptive fixed time convergence cooperative guidance method comprises the following steps:
step 1: establishing a missile relative target motion relation mathematical model;
assuming that n missiles strike a static target, the communication relation between every two missiles adopts an adjacency matrix A ═ a in graph theory ij ]That is, if the ith (i is 1,2, …, n) missile can establish a communication relationship with the jth (j is 1,2, …, n) missile, a is represented ij 1, otherwise, a ij 0; for convenience of reference, the physical quantities associated with the ith missile are indicated by the subscript i. If any two missiles in the communication topology can find at least one communication path, the communication topology map is called a connected map. If the links in the communication topology are all bidirectional, i.e. a ij =a ji The graph is called an undirected graph, and a directed graph if unidirectional links exist.
Aiming at the striking of a plurality of missiles on a fixed target under a two-dimensional plane, the communication topology among the missiles is assumed to be an undirected connected graph, and a relative kinematic model of the ith missile and the target is described by a formula (1):
Figure BDA0003020799260000041
in the formula, r i Indicates the distance between the ith missile and the target, V i Representing constant flight speed of the ith missile, q i Is the i-th missile line-of-sight angle, gamma i And phi i Respectively is a track angle and a lead angle of the ith missile; a is i The normal acceleration of the ith missile is vertical to the speed direction;
and 2, step: predicting the residual hit time of the missile;
residual hit time of ith missile
Figure BDA0003020799260000042
The following formula is used for prediction:
Figure BDA0003020799260000043
in the formula N s >2 denotes a navigation ratio;
and step 3: variable xi defining residual hit time consistency error i
Figure BDA0003020799260000044
In the formula (I), the compound is shown in the specification,
Figure BDA0003020799260000045
the remaining hit time of the jth missile;
and 4, step 4: designing continuous switching self-adaptive fixed time convergence cooperative guidance law a i
Figure BDA0003020799260000046
In the formula, k i For guidance gain, 0<μ 1 <1,μ 2 >1,s i =(sign(ξ i +δ)-sign(ξ i -δ))/2,δ>0 denotes the switching boundary, ∈ i Is a positive real number; definition of
Figure BDA0003020799260000047
Is phi i The upper bound of (a) is,
Figure BDA0003020799260000048
is composed of
Figure BDA0003020799260000049
Is determined by the estimated value of (c),
Figure BDA00030207992600000410
the values of (a) are determined by the adaptation law:
Figure BDA00030207992600000411
wherein gamma is d And σ is a positive real number.
The specific embodiment is as follows:
the Matlab simulation platform is used for verifying the designed method, and the designed method is combined with the traditional fixed time convergence cooperative guidance law in order to verify the effectiveness of the continuously switched self-adaptive fixed time convergence cooperative guidance law
Figure BDA0003020799260000051
Making a comparison, wherein k 1,i And k 2,i Is the guidance gain.
In the embodiment, 3 missiles at different positions are selected for simulation verification with the background of hitting a static target in the air, the topology of undirected strong connectivity among the missiles is shown in fig. 1, the target position (5000m and 5000m), the speeds of three missiles are all selected to be 200m/s, the initial horizontal axis coordinate position and the initial vertical axis coordinate position are respectively (500 ), (2000,500) and (500,2000), and the guidance law parameter is set to be k i =1.3,μ 1 =0.5,μ 2 =1.5,γ d =2.7,N s =3,δ=1.2,
Figure BDA0003020799260000052
ε i With 9, σ 0.1, the maximum overload limit is 30 g.
Simulation results under the two methods are respectively shown in fig. 2-5 and fig. 6-9, and the plane motion trail and the missile-eye distance curve show that although response curves are different, the two methods can ensure that three missiles hit a target simultaneously, and consistency errors can be rapidly converged. Compared with an acceleration instruction curve, the acceleration initial value of the continuously switched self-adaptive fixed time convergence cooperative guidance law is obviously smaller than that of the traditional fixed time convergence cooperative guidance method, and the result verifies that the designed guidance method can effectively reduce the impact of the initial control quantity.

Claims (2)

1. A continuously switchable self-adaptive fixed time convergence cooperative guidance method is characterized by comprising the following steps:
step 1: establishing a missile relative target motion relation mathematical model;
assuming that n missiles strike a static target, the communication relation between every two missiles adopts an adjacency matrix A ═ a in graph theory ij ]That is, if the ith missile can establish a communication relationship with the jth missile, j is 1,2, …, n missiles, i is 1,2, …, a ij 1, otherwise, a ij =0;
Aiming at the striking of a plurality of missiles on a fixed target under a two-dimensional plane, the communication topology among the missiles is assumed to be an undirected connected graph, and a relative kinematic model of the ith missile and the target is described by a formula (1):
Figure FDA0003722253700000011
in the formula, r i Indicates the distance, V, between the ith missile and the target i Representing constant flight speed of the ith missile, q i Is the i-th missile line-of-sight angle, gamma i And phi i Respectively is a track angle and a lead angle of the ith missile; a is i The normal acceleration of the ith missile is vertical to the speed direction;
step 2: predicting the residual hit time of the missile;
residual hit time of ith missile
Figure FDA0003722253700000012
The following formula is used for prediction:
Figure FDA0003722253700000013
in the formula N s >2 denotes the navigation ratio;
and step 3: variable xi defining residual hit time consistency error i
Figure FDA0003722253700000014
In the formula (I), the compound is shown in the specification,
Figure FDA0003722253700000015
the remaining hit time of the jth missile;
and 4, step 4: designing continuous switching self-adaptive fixed time convergence cooperative guidance law a i
Figure FDA0003722253700000016
In the formula, k i For guidance gain, 0<μ 1 <1,μ 2 >1,s i =(sign(ξ i +δ)-sign(ξ i -δ))/2,δ>0 denotes the switching boundary, ∈ i Is a positive real number; definition of
Figure FDA0003722253700000021
Is phi i The upper bound of (a) is,
Figure FDA0003722253700000022
is composed of
Figure FDA0003722253700000023
Is determined by the estimated value of (c),
Figure FDA0003722253700000024
the values of (a) are determined by the adaptation law:
Figure FDA0003722253700000025
wherein gamma is d And σ is a positive real number.
2. The adaptive fixed-time convergent cooperative guidance method for continuous switching according to claim 1, wherein k is i =1.3,μ 1 =0.5,μ 2 =1.5,γ d =2.7,N s =3,δ=1.2,
Figure FDA0003722253700000026
ε i =9,σ=0.1。
CN202110402142.0A 2021-04-14 2021-04-14 Self-adaptive fixed time convergence cooperative guidance method capable of realizing continuous switching Active CN113110576B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110402142.0A CN113110576B (en) 2021-04-14 2021-04-14 Self-adaptive fixed time convergence cooperative guidance method capable of realizing continuous switching

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110402142.0A CN113110576B (en) 2021-04-14 2021-04-14 Self-adaptive fixed time convergence cooperative guidance method capable of realizing continuous switching

Publications (2)

Publication Number Publication Date
CN113110576A CN113110576A (en) 2021-07-13
CN113110576B true CN113110576B (en) 2022-09-13

Family

ID=76717233

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110402142.0A Active CN113110576B (en) 2021-04-14 2021-04-14 Self-adaptive fixed time convergence cooperative guidance method capable of realizing continuous switching

Country Status (1)

Country Link
CN (1) CN113110576B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114415724B (en) * 2022-01-14 2024-06-04 东南大学 Three-dimensional multi-missile cooperative proportion guidance law design method allowing speed attenuation

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106843265A (en) * 2016-12-30 2017-06-13 哈尔滨工业大学 Three-dimensional many guided missile cooperative guidance method and systems of finite time convergence control
CN112113466A (en) * 2020-08-19 2020-12-22 南京理工大学 Multi-missile cooperative guidance method with limited time convergence
CN112129292A (en) * 2020-09-02 2020-12-25 北京航空航天大学 Heterogeneous aircraft cooperative guidance method considering attack time and attack angle constraints
CN112631328A (en) * 2020-12-12 2021-04-09 西北工业大学 Multi-target cooperative terminal guidance rate method

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10956803B2 (en) * 2017-04-03 2021-03-23 Wilcox Industries Corp. Modular rocket rescue system and passenger safety wearable band
CN112577373B (en) * 2020-11-27 2022-05-27 南京理工大学 Guidance law for multiple missiles to attack static target simultaneously based on proportion guidance law

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106843265A (en) * 2016-12-30 2017-06-13 哈尔滨工业大学 Three-dimensional many guided missile cooperative guidance method and systems of finite time convergence control
CN112113466A (en) * 2020-08-19 2020-12-22 南京理工大学 Multi-missile cooperative guidance method with limited time convergence
CN112129292A (en) * 2020-09-02 2020-12-25 北京航空航天大学 Heterogeneous aircraft cooperative guidance method considering attack time and attack angle constraints
CN112631328A (en) * 2020-12-12 2021-04-09 西北工业大学 Multi-target cooperative terminal guidance rate method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Adaptive Sliding Mode Guidance With Impact Time and Angle Constraints;Li Guofei;《IEEE Access》;20200204;第8卷;第26926-26932页 *
带落角约束的鸭/尾舵复合控制导弹制导控制一体化设计;凡永华;《弹箭与制导学报》;20170831;第37卷(第4期);第23-26页 *

Also Published As

Publication number Publication date
CN113110576A (en) 2021-07-13

Similar Documents

Publication Publication Date Title
CN108168381B (en) A kind of control method of more pieces of guided missile cooperations
CN110928329B (en) Multi-aircraft track planning method based on deep Q learning algorithm
CN111551080B (en) Guidance method for attack time control
CN106352738B (en) More bullet cooperative guidance methods based on output-consistence
CN112129292B (en) Heterogeneous aircraft cooperative guidance method considering attack time and attack angle constraints
CN108153328B (en) Multi-missile collaborative track planning method based on segmented Bezier curve
CN111336871B (en) Vertical attack guidance method based on circuitous flight
CN110187640A (en) For more guided missile cooperation Design of Guidance Law methods of maneuvering target and permission communication delay
CN114063644B (en) Unmanned fighter plane air combat autonomous decision-making method based on pigeon flock reverse countermeasure learning
CN113110576B (en) Self-adaptive fixed time convergence cooperative guidance method capable of realizing continuous switching
CN112379600B (en) Distributed cooperative guidance law construction method based on communication time-varying delay
CN115903865A (en) Aircraft near-distance air combat maneuver decision implementation method
CN112859921A (en) Three-dimensional simultaneous attack robust cooperative guidance law design method
CN115951695A (en) Dynamic tactical control domain resolving method based on three-party game in air combat simulation environment
CN111176315B (en) Variant cruise missile control method based on L1 adaptive control
CN109543135B (en) Multi-AUV (autonomous Underwater vehicle) cooperative simultaneous guidance method based on information maximization
CN114153143A (en) Design method of guided missile nonsingular fixed time sliding mode guidance law
CN110986948A (en) Multi-unmanned aerial vehicle grouping collaborative judgment method based on reward function optimization
CN113625745B (en) Attack time control guidance method based on switching fixed time convergence theory
CN114815878B (en) Hypersonic aircraft collaborative guidance method based on real-time optimization and deep learning
CN116225048A (en) Method for cooperatively striking tail end of fly-round projectile based on falling angle and time constraint
CN115328206A (en) Segmentation guidance method suitable for cooperative strike of multiple unmanned aerial vehicles
CN112379599B (en) Three-dimensional cooperative guidance method for air maneuvering target
CN114580138B (en) Bessel Newton-based missile multi-constraint terminal guidance law design method
CN116753788B (en) Agile turning composite control method based on linear active disturbance rejection

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