CN115220472A - Fault-tolerant control method for air-ground heterogeneous formation system under switching topology - Google Patents

Fault-tolerant control method for air-ground heterogeneous formation system under switching topology Download PDF

Info

Publication number
CN115220472A
CN115220472A CN202210819674.9A CN202210819674A CN115220472A CN 115220472 A CN115220472 A CN 115220472A CN 202210819674 A CN202210819674 A CN 202210819674A CN 115220472 A CN115220472 A CN 115220472A
Authority
CN
China
Prior art keywords
model
fault
derivative
follower
tolerant
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
CN202210819674.9A
Other languages
Chinese (zh)
Other versions
CN115220472B (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202210819674.9A priority Critical patent/CN115220472B/en
Publication of CN115220472A publication Critical patent/CN115220472A/en
Application granted granted Critical
Publication of CN115220472B publication Critical patent/CN115220472B/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/101Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • G05D1/104Simultaneous control of position or course in three dimensions specially adapted for aircraft involving a plurality of aircrafts, e.g. formation flying

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 discloses a fault-tolerant control method for an air-ground heterogeneous formation system under a switching topology, which comprises the following steps: (1) Designing a distributed topology network structure under a switching topology, designing a leader-follower formation strategy under the structure, and determining a Laplace matrix of a communication topology relation; (2) Considering the influence of a complex uncertain item, establishing a model for describing the dynamic characteristics of the air-ground heterogeneous multi-agent formation system; (3) Respectively designing a fixed observer based on a fixed time theory to carry out fixed time observation on the complex uncertain item of the second model and the complex uncertain item of the third model; (4) And designing a distributed formation fault-tolerant controller of an XY axis and a track fault-tolerant tracking controller of a Z axis which are converged at fixed time by adopting a high-order sliding mode control technology and a backstepping method according to the observed value and the Laplace matrix. The invention solves the multiple requirements of the air-ground heterogeneous formation system on adaptability, maneuvering responsiveness, fault-tolerant reliability, continuity of control quantity and the like of the switching topology network.

Description

Fault-tolerant control method for air-ground heterogeneous formation system under switching topology
Technical Field
The invention relates to the technical field of fault-tolerant formation control of a heterogeneous multi-agent system, in particular to a fault-tolerant control method of an air-ground heterogeneous formation system under a switching topology.
Background
In recent years, with the continuous development of artificial intelligence and control science, the intelligent agent has been developed in the aspects of environmental adaptability, autonomy, control accuracy, high efficiency and the like. However, as industrialization progresses faster, it becomes more and more difficult to achieve a desired complex control task by means of a single agent, and thus a multi-agent system constituted by a plurality of agents cooperating with each other is produced. When the dynamic model of one or several agents present in a multi-agent system is different from other agents, such a multi-agent system is called a heterogeneous multi-agent system. A multi-agent formation system composed of an aerial unmanned aerial vehicle and a ground unmanned vehicle is a typical heterogeneous multi-agent system, has various payload capacities, task configuration capacities, control and data acquisition capacities, and is widely applied to the fields of space exploration, collaborative rescue, resource exploration, path planning and the like.
However, in the air-ground cooperative formation process, due to the high complexity of the unmanned aerial vehicle structure, the constraints of individuals and their interaction characteristics and mechanisms, etc., the heterogeneous multi-agent system is susceptible to actuator faults, and the faults are propagated and spread through the interactive topological network, which further affects the stability and safety of the whole formation system. Therefore, the design of the distributed interactive network architecture also brings new challenges to the design of the fault-tolerant controller while improving the overall cooperativity. Especially, the sensing area of the sensor is limited due to the influence of environment factors such as variable and unreliable information transmission links and obstacles which execute tasks, which can cause the change of a topological network structure, and further cause the traditional fault-tolerant formation control strategy based on fixed topological design to be unusable. Moreover, conventional asymptotic responses and discontinuous control inputs have difficulty meeting anticipated control requirements in the face of critical technology systems such as air-ground unmanned systems. The formation control strategy of the existing air-ground unmanned system can only solve one problem in one way, and cannot solve the problems of quick maneuverability, continuity of control input, fault-tolerant control under a switching topology network and the like of the formation system at the same time, so that a new fault-tolerant control scheme capable of solving the series of problems needs to be researched to further improve the realization of the air-ground heterogeneous multi-agent formation system on complex tasks and the capability of maintaining the performance.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the defects, the invention discloses a fault-tolerant control method of an air-ground heterogeneous formation system under a switching topology, which designs continuous fixed time fault-tolerant tracking controllers of a Z-axis motion space in an XY plane under a switching topology network respectively, simultaneously solves the multiple requirements of the air-ground heterogeneous formation system on the adaptability, the maneuvering responsiveness, the fault-tolerant reliability, the continuity of a control quantity and the like of the switching topology network, and realizes the fault-tolerant formation control target of the air-ground heterogeneous formation system in an XYZ three-dimensional space.
The technical scheme is as follows: in order to solve the above problems, the present invention provides a fault-tolerant control method for a switching topology down-air-ground heterogeneous formation system, which specifically includes the following steps:
(1) Designing a distributed topology network structure under a switching topology aiming at an air-ground heterogeneous multi-agent formation system; designing a leader-follower formation strategy under the structure, selecting a unmanned vehicle as a leader and the other unmanned vehicles and/or unmanned aerial vehicles as followers, and determining a Laplacian matrix for describing a communication topological relation between the leader and the followers;
(2) Considering the influence of a complex uncertain item, establishing a model for describing the dynamic characteristics of the air-ground heterogeneous multi-agent formation system according to the kinematics equation of the unmanned aerial vehicle and the unmanned aerial vehicle, wherein the model comprises the following steps: a first model, a second model, a third model; the first model is a dynamic model of the unmanned vehicle of the leader, the second model is a dynamic model of the unmanned vehicle follower and the unmanned vehicle follower in an XY plane, and the third model is a dynamic model of the unmanned vehicle follower in a Z axis; the complex uncertainty comprises actuator faults, model uncertainty and external interference;
(3) Respectively designing a fixed observer based on a fixed time theory to carry out fixed time observation on the complex uncertainty of the second model and the complex uncertainty of the third model;
(4) Designing a distributed formation fault-tolerant controller of an XY axis converging in continuous fixed time by adopting a high-order sliding mode control technology and a backstepping method according to the first observation value and the Laplace matrix obtained in the step (1); the first observation value is an observation value of a fixed observer designed aiming at the complex uncertainty of the second model; designing a Z-axis track fault-tolerant tracking controller converging in continuous fixed time by adopting a high-order sliding mode control technology and a backstepping method according to a second observed value; the second observation value is an observation value of a fixed observer designed for the complex uncertainty of the third model.
Further, (1.1) defining that the air-ground heterogeneous multi-agent formation system is provided with N +1 agents, wherein the agents comprise unmanned planes and unmanned vehicles; wherein the unmanned vehicle has N 1 +1, unmanned aerial vehicle has N 2 N = N 1 +N 2 N is the number of followers; leader unmanned vehicle is labeled 0, follower unmanned vehicle is labeled i =1, \8230;, N 1 With follower drone marked i = N 1 +1,…,N 1 +N 2
(1.2) communication topology network undirected graph between followers is G = (S, E, A); s = { S = 1 ,s 2 ,…,s N Is a set of N followers,
Figure BDA0003743672020000021
for the set of communication connections between the followers,
Figure BDA0003743672020000022
denotes the ith 1 A follower and the ith 2 Communication connection of the individual follower and success of the connection, i 1 ∈i、i 2 ∈i;
Figure BDA0003743672020000023
A connection matrix representing the communication weight between followers, if
Figure BDA0003743672020000024
Then
Figure BDA0003743672020000025
Otherwise
Figure BDA0003743672020000026
The in-degree matrix of follower i is D = diag { deg. } in (s 1 ),...,deg in (s N )},
Figure BDA0003743672020000031
The laplace matrix between followers is L = D-a;
(1.3) the communication connection relation matrix between the leader and the follower is B = diag { B } 1 ,...,b i ,...,b N };b i Representing a communication connection between the ith follower and the leader; the interaction network of the whole formation system comprising the leader and the follower is
Figure BDA0003743672020000032
Its laplace matrix is H = L + B;
(1.4) since the communication topological relation between the followers is that the switching topological change occurs along with the time, a set containing all the topological structures is established:
Figure BDA0003743672020000033
theta represents the number of switching topologies present, and theta is greater than or equal to 1; the signal defining the switching topology is σ (t) [ [0, + ∞) → = {1, 2., Θ }, and the switching time is denoted as 0=t = 0 <t 1 <...<t P = t, in the time interval [ t p ,t p+1 ) The topology network between followers in the network is represented as a directed fixed topology network
Figure BDA0003743672020000034
Its Laplace matrix is defined as L σ(t) The communication matrix of the whole formation system is H σ(t) =L σ(t) +B σ(t)
Further, the step (2) specifically comprises:
(2.1) the expression of the first model is:
Figure BDA0003743672020000035
in the formula (I), the compound is shown in the specification,
Figure BDA0003743672020000036
for the first derivative of the position coordinate of the leader in the XY plane, v 0 (t) is the position speed of the leader in the XY plane at time t,
Figure BDA0003743672020000037
is v of 0 (t) first derivative; u. u 0 (t) inputting the control quantity of the position of the leader in the XY plane at the time t;
(2.2) the expression of the second model is:
Figure BDA0003743672020000038
Figure BDA0003743672020000039
in the formula, x i (t)=[x i ,y i ] T Is the position coordinate of the i-th follower in the XY plane,
Figure BDA00037436720200000310
is x i (ii) the second derivative of (t); u. of i xy (t) is the control quantity input of the i-th follower in the XY plane;
Figure BDA00037436720200000311
a complex uncertainty term for the second model;
Figure BDA0003743672020000041
an additive fault parameter for the ith follower; f. of i xy (t) values of the nonlinear portion of the system containing the model uncertainty for the ith follower;
Figure BDA0003743672020000042
external interference for the ith follower;
(2.3) the expression of the third model is:
Figure BDA0003743672020000043
Δ iz (t)=τ iz (t)+f iz (t)+d iz (t)-ρ iz u iz (t)/m i -g
in the formula, z i (t) is the position coordinate of the ith unmanned aerial vehicle on the Z axis,
Figure BDA0003743672020000044
is z i (ii) the second derivative of (t); u. of iz (t) inputting the control quantity of the ith unmanned aerial vehicle on the Z axis; delta iz (t) is the complex uncertainty term of the third model; tau is iz (t) additive fault parameters for the ith drone; f. of iz (t) values for the system non-linearity part of the model uncertainty for the ith drone; d iz (t) is the external interference suffered by the ith unmanned aerial vehicle; rho iz Executing an efficiency loss fault parameter for the ith drone; g is the acceleration of gravity; m is i Is the quality of the ith drone.
Further, the step (3) specifically comprises:
(3.1) designing a fixed time observer aiming at the complex uncertainty of the second model, wherein the formula is as follows:
Figure BDA0003743672020000045
in the formula, σ i1 、σ i2 Is a state variable of the fixed-time observer,
Figure BDA0003743672020000046
is composed of i1 The first derivative of (a) is,
Figure BDA0003743672020000047
is σ i2 The first derivative of (a); m is i1 Is a positive constant; wherein the content of the first and second substances,
Figure BDA0003743672020000048
Figure BDA0003743672020000049
Figure BDA00037436720200000410
in the formula o i1 、o i2 、o i3 Parameters of a fixed time observer;
Figure BDA00037436720200000411
a first derivative of a complex uncertainty term for the second model; w 1 Is a complex uncertainty term
Figure BDA00037436720200000412
The upper bound of (c);
Figure BDA00037436720200000413
is the first derivative of the location coordinate of the ith follower in the XY plane;
convergence time T of observer 1i Satisfies the following conditions:
Figure BDA0003743672020000051
in the formula, v i1 Is a positive constant of a setting, and
Figure BDA0003743672020000052
then, the minimum convergence time T is obtained 1i
(3.2) designing a fixed time observer aiming at the complex uncertainty of the third model, wherein the formula is as follows:
Figure BDA0003743672020000053
in the formula, σ i3 、σ i4 Is a state variable of the fixed-time observer,
Figure BDA0003743672020000054
is σ i3 The first derivative of (a) is,
Figure BDA0003743672020000055
is σ i4 The first derivative of (a); m is i2 Is a positive constant; wherein
Figure BDA0003743672020000056
Figure BDA0003743672020000057
o i5 >0,o i6 >4W 2
Figure BDA0003743672020000058
In the formula o i4 、o i5 、o i6 Parameters of a fixed time observer;
Figure BDA0003743672020000059
a first derivative of a complex uncertainty term for the third model; w 2 Is a complex uncertainty term
Figure BDA00037436720200000510
The upper bound of (c);
Figure BDA00037436720200000511
a first derivative of a position coordinate of the ith unmanned aerial vehicle on the Z axis is obtained;
convergence time T of observer 1iz Satisfies the following conditions:
Figure BDA0003743672020000061
in the formula, v i2 Is provided withIs constant, positive constant, and
Figure BDA0003743672020000062
then, the minimum convergence time T is obtained 1iz
Further, the step (4) specifically includes:
(4.1) defining an expected formation vector containing all followers in the XY plane as
Figure BDA0003743672020000063
The predetermined track of all the unmanned aerial vehicle followers on the Z axis is
Figure BDA0003743672020000064
There is a settable time constant T that is not limited to the initial state of the system max So that the following holds:
Figure BDA0003743672020000065
(4.2) based on the observed value σ i2 Laplace matrix H of communication topological relation among leader and follower σ(t) Designing a distributed fault-tolerant formation controller of an XY axis, wherein the expression is as follows:
Figure BDA0003743672020000066
Figure BDA0003743672020000067
the virtual control law is as follows:
Figure BDA0003743672020000068
Figure BDA0003743672020000069
in the formula, chi 1 、δ 1 、χ 2 、δ 2 、χ 3 And delta 3 Are all positive constant, c 1 、c 2 、c 3 Are all control parameters, d 1 、d 2 、d 3 Are constants with a value greater than 1; zeta i1 、ζ i2 、ζ i3 A third order system which is the controller;
Figure BDA00037436720200000610
is gamma i2 The first derivative of (a);
Figure BDA00037436720200000611
is gamma i1 The first derivative of (a);
Figure BDA00037436720200000612
is the first derivative of the formation tracking error;
(4.3) from the observed value σ i4 Designing a Z-axis track fault-tolerant tracking controller, wherein the expression is as follows:
Figure BDA0003743672020000071
Figure BDA0003743672020000072
and the virtual control amount is designed as:
Figure BDA0003743672020000073
Figure BDA0003743672020000074
wherein, mu 1 、ω 1 、μ 2 、ω 2 、μ 3 And ω 3 Is a positive constant, n, set manually 1 、n 2 、n 3 Are all control parameters, m 1 、m 2 、m 3 Are constants with a value greater than 1;
Figure BDA0003743672020000075
a third order system in which the controllers are all the same;
Figure BDA0003743672020000076
is composed of
Figure BDA0003743672020000077
The first derivative of (a);
Figure BDA0003743672020000078
is composed of
Figure BDA0003743672020000079
The first derivative of (a);
Figure BDA00037436720200000710
the first derivative of the trajectory tracking error in the Z-axis.
Further, the step (4.2) specifically comprises:
(4.2.1) designing a fixed time integral sliding mode surface based on the tracking error:
Figure BDA00037436720200000711
in the formula, e i =[e ix ,e iy ] T ,sig l (e i )=[|e ix | l sgn(e ix ),|e iy | l sgn(e iy )] T Alpha and beta are both positive definite constants and satisfy 0 < alpha < 1 and beta > 1 1 、b 2 Are all constants; e.g. of the type ix For the formation tracking error of the x-axis, e iy A formation tracking error for the y-axis;
taking a first derivative and a second derivative of the sliding mode surface to obtain:
Figure BDA00037436720200000712
Figure BDA00037436720200000713
in the formula, λ i1 =diag{αb 1 |e ix | α-1 ,αb 1 |e iy | α-1i2 =diag{βb 2 |e ix | β-1 ,βb 2 |e iy | β-1 };
Obtaining a third-order system of an integral sliding mode surface:
Figure BDA00037436720200000714
Figure BDA00037436720200000715
Figure BDA00037436720200000716
(4.2.2) based on the design scheme of the fault-tolerant controller of the backstepping method, the following transformations are carried out on the three-order system:
ζ i1 =S i1
ζ i2 =S i2i1
ζ i3 =S i3i2
(4.2.3) designing a distributed fault-tolerant formation controller in an XY plane under a switching topology network based on the established new third-order system, wherein the distributed fault-tolerant formation controller realizes fixed time convergence.
Further, the step (4.3) specifically comprises:
(4.3.1) designing a fixed time integral sliding mode surface based on the tracking error:
Figure BDA0003743672020000081
in the formula (I), the compound is shown in the specification,
Figure BDA0003743672020000082
α z and beta z Are all positive definite constants; b 1z >0,b 2z The constants are artificially set when the value is more than 0;
taking a first derivative and a second derivative of the sliding mode surface, and finally obtaining a third-order system of the integral sliding mode surface:
Figure BDA0003743672020000083
Figure BDA0003743672020000084
Figure BDA0003743672020000085
in the formula (I), the compound is shown in the specification,
Figure BDA0003743672020000086
(4.3.2) the fault-tolerant controller design scheme based on the backstepping method is used for carrying out the following transformation on the three-order system:
Figure BDA0003743672020000087
Figure BDA0003743672020000088
Figure BDA0003743672020000089
(4.3.3) designing a track fault-tolerant tracking controller on the Z axis based on the established new third-order system under a switching topology network, wherein the track fault-tolerant tracking controller realizes fixed time convergence.
Furthermore, the present invention also provides a computer device, comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of any of the above methods when executing the computer program. A computer-readable storage medium, having stored thereon a computer program which, when being executed by a processor, carries out the steps of any of the methods described above.
Has the beneficial effects that: compared with the prior art, the fault-tolerant control method for the air-ground heterogeneous formation system under the switching topology has the advantages that: 1. aiming at the difference between an unmanned aerial vehicle model and an unmanned vehicle model in an air-ground heterogeneous formation system, a two-dimensional dynamic model of all followers in an XY axis and a dynamic model of an unmanned aerial vehicle in a Z axis are respectively established, a fixed time observer is respectively established for the two established models to realize effective estimation on a complex uncertain item, observed values and neighbor interaction information are utilized in an XY plane, a distributed continuous fixed time fault-tolerant controller under a switching topology network is designed, the controller is designed to integrate a high-order sliding mode technology and a backstep method, not only is the artificial setting on convergence time realized, but also the output of the controller is ensured to be continuous, the fault tolerance and the robustness of the formation system are effectively improved, and meanwhile, the generalization of the designed controller in practical application is greatly improved; 2. a track fault-tolerant tracking controller based on an observer is designed on a Z axis by utilizing a local tracking error, a fixed time theory, a high-order sliding mode control technology and a backstepping method, and the safe and rapid tracking of an expected track by a follower of the unmanned aerial vehicle on the vertical height is ensured.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is an overall block diagram of the air-ground heterogeneous formation system according to the present invention;
FIG. 3 is a diagram illustrating a switching topology of an air-ground heterogeneous formation system; FIG. 3 (a), FIG. 3 (b), FIG. 3 (c), FIG. 3 (d) are all a topology structure diagram;
FIG. 4 is a diagram showing the formation effect of the XY space air-ground heterogeneous formation system according to the present invention;
FIG. 5 is a diagram showing formation error convergence of the XY-axis air-ground heterogeneous formation system according to the present invention; FIG. 5 (a) is a diagram of the formation error of the follower on the X-axis, and FIG. 5 (b) is a diagram of the formation error of the follower on the Y-axis;
FIG. 6 is a tracking error convergence diagram of the Z-axis air-ground heterogeneous formation system according to the present invention;
FIG. 7 is a graph showing the continuous control output of the air-ground heterogeneous formation system of XY axes according to the present invention;
fig. 7 (a) is a control amount output diagram of the follower on the X axis, and fig. 7 (b) is a control amount output diagram of the follower on the Y axis;
FIG. 8 is a graph showing the continuous control output of the Z-axis air-ground heterogeneous formation system according to the present invention;
Detailed Description
The technical scheme of the invention is further explained by combining the attached drawings.
As shown in fig. 1, the present invention provides a fault-tolerant control method for a switching topology down-to-ground heterogeneous formation system, which specifically includes the following steps:
firstly, designing a distributed topology network structure under a switching topology aiming at an air-ground heterogeneous multi-agent formation system; designing a leader-follower formation strategy under the structure, selecting a unmanned vehicle as a leader and the other unmanned vehicles or unmanned aerial vehicles as followers, and determining a Laplacian matrix describing the communication topological relation between the leader and the followers;
(1) Defining an air-ground heterogeneous multi-agent formation system with N +1 agents (including unmanned aerial vehicles and unmanned vehicles); wherein N is 1 +1 unmanned vehicles, N 2 An unmanned aerial vehicle; leader unmanned vehicle is labeled 0, follower unmanned vehicle is labeled i =1, \8230;, N 1 With follower drone marked i = N 1 +1,…,N 1 +N 2 ;N=N 1 +N 2 N is the number of followers;
(1.2) communication topology network undirected graph between followers G =(S,E,A);S={s 1 ,s 2 ,…,s N Is the set of N followers,
Figure BDA0003743672020000101
for the set of communication connections between the followers,
Figure BDA0003743672020000102
denotes the ith 1 A follower and the ith 2 Communication connection of the individual follower and success of the connection, i 1 ∈i、i 2 ∈i;
Figure BDA0003743672020000103
A connection matrix representing the communication weight between followers, if
Figure BDA0003743672020000104
Then the
Figure BDA0003743672020000105
Otherwise
Figure BDA0003743672020000106
The in-degree matrix for follower i is D = diag { deg. } in (s 1 ),...,deg in (s N )},
Figure BDA0003743672020000107
The laplace matrix between followers is L = D-a;
(3) The communication connection relation between the leader and the follower is expressed by a diagonal matrix as follows: b = diag { B } 1 ,...,b i ,...,b N }; wherein, b i Represents the communication between the ith follower and the leader, if b i > 0 means that the ith follower can obtain information from the leader, otherwise communication between the i followers and the leader is not possible.
The interaction network of the whole formation system comprising the leader and the follower is
Figure BDA0003743672020000108
Its laplace matrix is H = L + B; that is, H describes the communication connection matrix of the entire formation system topology network, and H describes the communication network in which the leader has at least one path that can reach any follower.
(4) For the case that the communication topology relationship between the followers is the switching topology change over time, it is necessary to establish a set including all the topologies and an expression form of the switching topology.
Define the set of all directed graphs as
Figure BDA0003743672020000109
Θ represents the number of handover topologies present; the signal defining the switching topology is σ (t): [0, + ∞) → = {1, 2., Θ } and the switching time is denoted 0=t = 0 <t 1 <...<t P =t,t p For a particular switching instant, p is the switching order, i.e. at t p The P-th switching occurs at the moment, and P is the total switching frequency. In the time interval t p ,t p+1 ) The topological network among followers (P ≦ P) is represented as a directed fixed topological network
Figure BDA00037436720200001010
Its Laplace matrix is defined as L σ(t) The communication matrix of the whole formation system is H σ(t) =L σ(t) +B σ(t)
As shown in fig. 2, the embodiment has 5 agents, 3 unmanned vehicles and 2 unmanned vehicles; one unmanned vehicle (i = 0) is taken as a leader, two unmanned vehicles (i =1,2) and two unmanned vehicles (i =3,4) are additionally selected as followers, and the leader and the followers jointly form the air-ground heterogeneous multi-agent formation system.
A specific switching topology distributed topology network structure is designed for the air-ground heterogeneous multi-agent formation system and is shown in figure 3, each switching topology structure is maintained for 5 seconds, namely the communication topology of the formation system in 0 s-5 s is shown in figure 3 (a), and in the period of time
Figure BDA0003743672020000111
The communication topology between 5s and 10s is shown in FIG. 3 (b), during which time
Figure BDA0003743672020000112
The communication topology between 10s and 15s is shown in FIG. 3 (c), during which time
Figure BDA0003743672020000113
The communication topology between 15s and 20s is shown in FIG. 3 (d), during which time
Figure BDA0003743672020000114
Secondly, establishing a model for describing the dynamic characteristics of the air-ground heterogeneous multi-agent formation system according to the kinematics equation of the unmanned aerial vehicle and the unmanned vehicle;
(1) Defining a dynamic model of the leader unmanned vehicle as a first model, wherein the expression of the first model is as follows:
Figure BDA0003743672020000115
in the formula, x 0 ∈R 2 For the position coordinates, v, of the leader in the XY plane 0 ∈R 2 Is the position speed of the leader in the XY plane, u 0 ∈R 2 Inputting a control quantity for the position of the leader in the XY plane; r represents a real number domain.
The initial position of the leader in this embodiment is x 0 (0)=[0,0] T The initial state of the follower is x 1 (0)=[-0.1,0.2] T ,x 2 (0)=[0.5,-0.3] T ,x 3 (0)=[-0.4,0] T And x 4 (0)=[0,0.6] T The leader's input is u 0 =[0.02,0.02] T
(2) And (3) fully considering the influences of the fault of the intelligent actuator, model uncertainty and external interference, and then establishing a unified two-dimensional dynamic model of the unmanned aerial vehicle and the unmanned aerial vehicle in the XY plane, wherein the model is defined as a second model. Specifically, the ith follower (i ∈ (1,.., N) is added 1 +N 2 ) The model in the XY plane is simplified to the model expression shown below:
Figure BDA0003743672020000121
wherein, χ i (t)=[x i ,y i ] T ∈R 2 And
Figure BDA0003743672020000122
respectively represents the position coordinate of the ith follower in the XY plane and the input of the control quantity,
Figure BDA0003743672020000123
representing the actuator efficiency loss factor for the ith follower,
Figure BDA0003743672020000124
an additive failure parameter, f, representing the ith follower i xy (t)∈R 2 Indicating that the ith follower contains the nonlinear part of the system of the model uncertainty,
Figure BDA0003743672020000125
represents the external interference suffered by the ith follower, i belongs to (1.. The., N) when the ith follower is an unmanned vehicle 1 ) Then A i =I 2 ,I 2 Is a two-dimensional identity matrix. Conversely, when the ith follower is a drone, i.e., i e (N) 1 +1,...,N 1 +N 2 ) Then A i =diag{1/m i ,1/m i In which m is i Is the mass of the ith drone.
In order to facilitate the design of the controller, the second model expression is further simplified into the following form by designing a complex uncertainty term containing actuator faults, model uncertainty and external interference:
Figure BDA0003743672020000126
Figure BDA0003743672020000127
in the formula (I), the compound is shown in the specification,
Figure BDA0003743672020000128
is a complex uncertainty term for the second model.
(3) In a similar way, the influence of the fault of the intelligent actuator, the uncertainty of the model and the external interference is fully considered, the dynamic model of the unmanned aerial vehicle on the Z axis is a third model, and the expression of the third model is as follows:
Figure BDA0003743672020000129
wherein z is i E R and u iz (t) belongs to R and is the position coordinate and the control input quantity of the ith unmanned aerial vehicle on the Z axis respectively, and rho iz E.g. R and tau iz (t) is the ith unmanned aerial vehicle execution efficiency loss fault parameter and the additive fault parameter respectively, f iz (t) epsilon R and d iz (t) is the model nonlinear part of the ith unmanned aerial vehicle and the external interference, and g represents the gravity acceleration.
In order to facilitate the design of the controller, the second model expression is further simplified into the following form by designing complex uncertainties including actuator faults, model uncertainty and external interference:
Figure BDA0003743672020000131
Δ iz (t)=τ iz (t)+f iz (t)+d iz (t)-ρ iz u iz (t)/m i -g
in the formula,. DELTA. iz (t) is the complex uncertainty term of the third model.
In this embodiment, the execution efficiency loss module is set as:
Figure BDA0003743672020000132
Figure BDA0003743672020000133
Figure BDA0003743672020000134
Figure BDA0003743672020000135
and an additive fault module as shown below:
τ 1 =[-0.2cos(t),-sin(0.5t)] T2 =[-0.15sin(1.2t),-0.45sin(2t)] T ,
the external interference suffered by each follower is designed as follows:
d 1 =[0.1sin(t),-0.3cos(0.7t)] T ,d 2 =[0.45cos(t),0.7sin(2t)] T ,
d 3 =[-0.3cos(0.5t),-0.2sin(1.5t),0.27sin(2t)] T ,
d 4 =[0.5sin(0.25t),0.55cos(1.5t),0.3sin(2t)] T
in addition, considering that each drone follower is affected by model uncertainty, there is an uncertainty change in the design model parameters with 20% amplitude.
Step three, respectively designing a fixed observer pair based on a fixed time theory
Figure BDA0003743672020000136
Δ iz (t) carrying out fixed time observation, wherein the convergence time of the observation error is not limited by the initial state of the system;
(1) Design of fixed time observer pair
Figure BDA0003743672020000141
And (3) carrying out fixed time observation, wherein the expression is as follows:
Figure BDA0003743672020000142
in the formula, σ i1 、σ i2 Are state variables of the fixed time observer.
Wherein m is i1 Is a positive constant set by people, and satisfies the following conditions: m is i1 >1,
Figure BDA0003743672020000143
Observer parameter o i1 、o i2 、o i3 Satisfies the following conditions:
Figure BDA0003743672020000144
o i2 >0,o i3 >4W 1
in the formula, W 1 Is a complex uncertainty term
Figure BDA0003743672020000145
Upper bound of the derivative, i.e. satisfies
Figure BDA0003743672020000146
Complex uncertainty term
Figure BDA0003743672020000147
Can be respectively sigma i2 Estimated within a fixed time. As will be described below
Figure BDA0003743672020000148
A fixed time convergence proof is performed for the example.
By pairs of s 1 Differentiation can yield:
Figure BDA0003743672020000149
order:
Figure BDA00037436720200001410
then differentiating it can yield:
Figure BDA00037436720200001411
to this end, complex uncertainties
Figure BDA00037436720200001412
The system of observation errors of (a) can be expressed as:
Figure BDA00037436720200001413
the error system described above can complete its convergence process in a fixed time, i.e. σ i2 Can realize the uncertainty of the complex in a fixed time
Figure BDA00037436720200001414
And the convergence time T of the observer 1i Satisfies the following conditions:
Figure BDA0003743672020000151
wherein, v i1 Is greater than 0, and when
Figure BDA0003743672020000152
A minimum convergence time T can be obtained 1i
(2) Design of fixed time observer pair Δ iz (t) performing fixed time observation, wherein the expression is as follows:
Figure BDA0003743672020000153
in the formula, σ i3 、σ i4 All are fixed time observersIs measured.
Wherein m is i2 >1,
Figure BDA0003743672020000154
Observer parameter o i4 、o i5 And o i6 Satisfies the following conditions:
Figure BDA0003743672020000155
o i5 >0,o i6 >4W 2
in the formula, W 2 Is a complex uncertainty term Δ iz (t) upper bound of derivative, i.e. satisfies
Figure BDA0003743672020000156
Then, the complex uncertainty term Δ iz (t) may be respectively σ i4 Estimated over a fixed time. The complex uncertainty Δ in Z-space can be obtained using similar proofs as in step (3.1) iz (T) may also be at a fixed time T 1iz Accurate estimation is achieved, and T 1iz Satisfies the following conditions:
Figure BDA0003743672020000157
wherein, v i2 > 0, and when
Figure BDA0003743672020000158
A minimum convergence time T can be obtained 1iz
In this embodiment, each parameter is set as: o. o i1 =o i4 =5,o i2 =o i5 =2,o i3 =o i6 =4。
And step four, after fixed-time observation of complex uncertainty is realized, designing a fault-tolerant formation controller based on an observed value and by utilizing a high-order sliding mode control technology and a backstepping method, wherein the design of the controller is divided into an XY-axis distributed fault-tolerant formation control part and a Z-axis distributed track fault-tolerant tracking control part.
(1) Defining an expected formation vector containing all follower agents in the XY plane
Figure BDA0003743672020000161
Predetermined trajectory of all Z-axis drone followers
Figure BDA0003743672020000162
In an XYZ three-dimensional space, an expected formation task is that a settable time constant T which can not be limited to the initial state of a system exists in the condition that an air-ground heterogeneous formation system suffers from actuator faults, model uncertainty and external interference max So that the following holds:
Figure BDA0003743672020000163
in this embodiment, the formation form expected by each follower is a time-varying vector, specifically:
Figure BDA0003743672020000164
and
Figure BDA0003743672020000165
(2) Designing a fixed time integral sliding mode surface by utilizing the formation tracking error in an XY axis and the Z axis track tracking error to ensure that the fixed time convergence can be realized by the tracking error after the sliding mode surface is reached; then, a new three-order error system is constructed by utilizing the fixed time sliding mode surface so as to facilitate the designed controller to realize continuous output; and then designing the XY-axis distributed formation fault-tolerant controller by using a backstepping method and interactive information between agents.
(2.1) designing a distributed fault-tolerant controller of an XY axis:
(2.1.1) first, a fixed time integral sliding mode surface is designed based on the tracking error as follows:
Figure BDA0003743672020000166
wherein e is i =[e ix ,e iy ] T ,sig l (e i )=[|e ix | l sgn(e ix ),|e iy | l sgn(e iy )] T Alpha and beta are positive constants, and satisfy 0 < alpha < 1 and beta > 1. In addition, b 1 >0,b 2 And the constant is artificially set when the value is more than 0. e.g. of a cylinder ix For the formation tracking error of the x-axis, e iy For the formation tracking error of the y-axis, e i Defined as the error in formation as:
Figure BDA0003743672020000167
the fixed time convergence of the slip form surface is given below, and after reaching the slip form surface, the following lyapunov function is designed:
Figure BDA0003743672020000168
taking the derivative thereof, we can get:
Figure BDA0003743672020000171
wherein M is e i The designed integral sliding mode surface can realize the fixed time convergence of the global formation error after reaching the sliding mode surface, and the upper bound of the convergence time meets the following requirements:
Figure BDA0003743672020000172
(2.2.2) taking the first derivative and the second derivative of the sliding mode surface to obtain:
Figure BDA0003743672020000173
Figure BDA0003743672020000174
wherein λ is i1 =diag{αb 1 |e ix | α-1 ,αb 1 |e iy | α-1 And λ i2 =diag{βb 2 |e ix | β-1 ,βb 2 |e iy | β-1 }。
To this end, the following third order system for the integral sliding mode surface can be obtained:
Figure BDA0003743672020000175
Figure BDA0003743672020000176
Figure BDA0003743672020000177
(2.2.3) based on a back-stepping method, carrying out the following transformation on the third-order system:
ζ i1 =S i1
ζ i2 =S i2i1
ζ i3 =S i3i2
wherein, γ i1 And gamma i2 Is a virtual control quantity to be designed.
(2.2.4) based on the established new third-order system, a distributed fault-tolerant formation controller in the XY plane is designed under a switching topology network, and convergence of formation errors in fixed time is realized. The expression of the controller is:
Figure BDA0003743672020000178
Figure BDA0003743672020000179
the virtual control law is as follows:
Figure BDA0003743672020000181
Figure BDA0003743672020000182
wherein, χ 1 、δ 1 、χ 2 、δ 2 、χ 3 And delta 3 Is a positive constant set by human, control parameter c 1 、c 2 、c 3 E (0, 1), and d 1 ,d 2 ,d 3 >1。
In this embodiment: b is a mixture of 1 =1,b 2 =1.1,α=0.5,β=1.1,c 1 =c 2 =c 3 =0.5,d 1 =d 2 =d 3 =1.4,χ 1 =11,δ 1 =10,χ 2 =1,δ 2 =0.8,χ 3 =1,δ 3 =2。
The proof of the fixed time convergence for this designed controller is given below:
first, the following Lyapunov function is designed:
Figure BDA0003743672020000183
deriving and introducing a control law into the method to obtain:
Figure BDA0003743672020000184
after the fixed time observer finishes observing the complex uncertain items, the following can be obtained:
Figure BDA0003743672020000185
therefore, when the time T ≧ T 1 In which T is 1 =max{T 11 ,T 12 ,...,T 1N },ζ i3 Can realize fixed time convergence and can obtain the convergence time upper bound of T 1 +T 3 Wherein:
Figure BDA0003743672020000191
then, the following lyapunov function is further designed:
Figure BDA0003743672020000192
based on the previously constructed third order system, the following can be derived:
Figure BDA0003743672020000193
and because when T is more than or equal to T 1 +T 3 When there is
Figure BDA0003743672020000194
And ζ i3 =0, at this time, for V 3 The derivation may be:
Figure BDA0003743672020000195
thus, it can be obtained that when T > T 1 +T 3 +T 4 Time, ζ i2 Converge to zero, and wherein:
Figure BDA0003743672020000196
then, the following lyapunov function is designed:
Figure BDA0003743672020000197
at T > T 1 +T 3 +T 4 Then, taking the derivative of it, we can get:
Figure BDA0003743672020000201
wherein the content of the first and second substances,
Figure BDA0003743672020000202
λ min (H σ(t) ) Is H σ(t) Is determined by the minimum characteristic value of (c),
Figure BDA0003743672020000203
λ max (H σ(t) ) Is H σ(t) The maximum eigenvalue of (c). Based on the above analysis, it can be found that 1 Is a fixed time convergence, i.e. for the ith (i ∈ (1.,. N.) 1 +N 2 ) For an agent, the designed fixed-time fault-tolerant controller of the XY plane is at T ≧ T 1 +T 3 +T 4 +T 5 Can reach the sliding form surface S i1 . And the sliding mode surface can be T > T 2 Then realizing the convergence of the formation error, thus obtaining that when T is more than T max =T 1 +T 2 +T 3 +T 4 +T 5 The desired formation task may be implemented.
(2.2) designing a Z-axis track fault-tolerant tracking controller:
similar to the design mode of the distributed fault-tolerant formation controller in the XY plane, firstly, a fixed-time sliding mode surface is constructed based on the tracking error, then a three-order system is constructed by utilizing the sliding mode surface, a virtual control quantity is designed to form a new error system, and then the error system is proved to realize fixed-time convergence under the control of the designed fault-tolerant controller based on the observer based on a back stepping method.
(2.2.1) first, a fixed time integral sliding mode surface is designed based on the tracking error of the Z axis as follows:
Figure BDA0003743672020000204
wherein alpha is z And beta z Is a positive constant and satisfies 0 < alpha z < 1 and beta z Is greater than 1. In addition, b 1z >0,b 2z The constant is artificially set when the value is more than 0.
The three-order system about the Z-axis integral sliding mode surface can be arranged as follows:
Figure BDA0003743672020000211
Figure BDA0003743672020000212
Figure BDA0003743672020000213
wherein the content of the first and second substances,
Figure BDA0003743672020000214
and
Figure BDA0003743672020000215
(2.2.2) further performing the following transformation on the third-order system:
Figure BDA0003743672020000216
Figure BDA0003743672020000217
Figure BDA0003743672020000218
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003743672020000219
and
Figure BDA00037436720200002110
is a virtual control quantity to be designed.
(2.2.3) designing a Z-axis track fault-tolerant tracking controller expression as follows:
Figure BDA00037436720200002111
Figure BDA00037436720200002112
and the virtual control amount is designed as follows:
Figure BDA00037436720200002113
Figure BDA00037436720200002114
wherein, mu 1 、ω 1 、μ 2 、ω 2 、μ 3 And ω 3 Is a positive constant set by human, and a control parameter n 1 ,n 2 ,n 3 E (0, 1), and m 1 ,m 2 ,m 3 Is greater than 1. The controller of this design has a proof of fixed time convergence similar to the proof of the controller described above.
Based on the parameter setting in the above embodiment, simulation verification is performed by applying a Simulink module in Matlab. As shown in fig. 4, the two-dimensional queuing process diagram illustrates that the controller is designed to enable the air-ground queuing system to complete a predetermined fault-tolerant queuing control task under a switching topology network. As shown in fig. 5 and fig. 6, the convergence process of the error-tolerant controller designed by the present invention in the XY axis and Z axis formation error is shown respectively, even under the influence of actuator failure, model uncertainty and external interference, the designed controller can still complete the convergence process quickly and reliably within a specified time, and ensure the formation and maintenance performance of the formation system in switching topology network to the geometric formation configuration within a fixed time, so that the robustness and the error tolerance of the whole formation system can be significantly improved by quickly estimating and effectively compensating the complex uncertainty. As shown in fig. 7 and fig. 8, which respectively show the control quantities of the controller designed by the present invention in the XY axis and the Z axis, it can be found that the controller realizes stable and continuous control output, and can autonomously adjust the amplitude of the control quantity for the sudden occurrence of the actuator fault, and the designed controller has relatively strong adaptive capability and fault compensation capability.

Claims (9)

1. A fault-tolerant control method for an air-ground heterogeneous formation system under switching topology is characterized by comprising the following steps:
(1) Designing a distributed topology network structure under a switching topology aiming at an air-ground heterogeneous multi-agent formation system; designing a leader-follower formation strategy under the structure, selecting a unmanned vehicle as a leader and the other unmanned vehicles and/or unmanned aerial vehicles as followers, and determining a Laplacian matrix for describing a communication topological relation between the leader and the followers;
(2) Considering the influence of a complex uncertainty item, building a model for describing the dynamic characteristics of the air-ground heterogeneous multi-agent formation system according to the kinematic equations of the unmanned aerial vehicle and the unmanned vehicle in the air-ground heterogeneous multi-agent formation system, wherein the model comprises the following steps: a first model, a second model, a third model; the first model is a dynamic model of the unmanned vehicle of the leader, the second model is a dynamic model of the unmanned vehicle follower and the unmanned vehicle follower in an XY plane, and the third model is a dynamic model of the unmanned vehicle follower in a Z axis; the complex uncertainty comprises actuator faults, model uncertainty and external interference;
(3) Respectively designing a fixed observer based on a fixed time theory to carry out fixed time observation on the complex uncertain item of the second model and the complex uncertain item of the third model;
(4) Designing a distributed formation fault-tolerant controller of an XY axis converging in continuous fixed time by adopting a high-order sliding mode control technology and a backstepping method according to the first observation value and the Laplace matrix obtained in the step (1); the first observation value is an observation value of a fixed observer designed aiming at a complex uncertainty item of the second model; designing a Z-axis track fault-tolerant tracking controller which converges in continuous fixed time by adopting a high-order sliding mode control technology and a backstepping method according to a second observation value; the second observation value is the observation value of a fixed observer designed aiming at the complex uncertainty of the third model.
2. The fault-tolerant control method for the air-ground heterogeneous formation system under the switching topology according to claim 1, wherein the step (1) specifically comprises:
(1.1) defining an air-ground heterogeneous multi-agent formation system to be provided with N +1 agents, wherein the agents comprise an unmanned aerial vehicle and an unmanned vehicle; wherein the unmanned vehicle has N 1 +1, unmanned aerial vehicle has N 2 N = N 1 +N 2 N is the number of followers; leader unmanned vehicle label 0, follower unmanned vehicle label i =1, \8230, N 1 Follower drone is marked i = N 1 +1,…,N 1 +N 2
(1.2) communication topology network undirected graph between followers is marked as G = (S, E, A); s = { S = 1 ,s 2 ,…,s N Is a set of N followers,
Figure FDA0003743672010000011
for the set of communication connections between the followers,
Figure FDA0003743672010000012
denotes the ith 1 A follower and the ith 2 Communication connection and success of each follower, i 1 ∈i、i 2 ∈i;
Figure FDA0003743672010000013
A connection matrix representing the communication weight between followers, if
Figure FDA0003743672010000021
Then
Figure FDA0003743672010000022
Otherwise
Figure FDA0003743672010000023
The in-degree matrix for follower i is D = diag { deg. } in (s 1 ),...,deg in (s N )},
Figure FDA0003743672010000024
The laplace matrix between followers is L = D-a;
(1.3) the communication connection relation matrix between the leader and the follower is B = diag { B } 1 ,...,b i ,...,b N },b i Representing a communication connection between the ith follower and the leader; the interaction network of the whole formation system comprising the leader and the follower is
Figure FDA0003743672010000025
Its laplace matrix is H = L + B;
(1.4) since the communication topological relation between the followers is that the switching topological change occurs along with the time, a set containing all the topological structures is established:
Figure FDA0003743672010000026
theta represents the number of existing switching topologies, and theta is more than or equal to 1; the signal defining the switching topology is σ (t): [0, + ∞) → = {1, 2., Θ } and the switching time is denoted 0=t = 0 <t 1 <...<t P = t, in time interval [ t p ,t p+1 ) The topology network between followers in the network is represented as a directed fixed topology network
Figure FDA0003743672010000027
Its Laplace matrix is defined as L σ(t) (ii) a The communication matrix of the whole formation system is H σ(t) =L σ(t) +B σ(t)
3. The fault-tolerant control method for the air-ground heterogeneous formation system under the handover topology according to claim 2, wherein the step (2) specifically comprises:
(2.1) the expression of the first model is:
Figure FDA0003743672010000028
in the formula (I), the compound is shown in the specification,
Figure FDA00037436720100000214
for the first derivative, v, of the position coordinate of the leader in the XY plane at time t 0 (t) is the position speed of the leader in the XY plane at time t,
Figure FDA00037436720100000210
is v of 0 (t) first derivative; u. of 0 (t) inputting the control quantity of the position of the leader in the XY plane at the time t;
(2.2) the expression of the second model is:
Figure FDA00037436720100000211
Figure FDA00037436720100000212
wherein, χ i (t) = [ x ] i ,y i ] T Is the position coordinate of the ith follower in the XY plane,
Figure FDA00037436720100000213
is x i (t) second derivative; u. of i xy (t) is the control quantity input of the i-th follower in the XY plane;
Figure FDA0003743672010000031
a complex uncertainty term for the second model;
Figure FDA0003743672010000032
an additive fault parameter for the ith follower; f. of i xy (t) is the value of the nonlinear part of the system that the i-th follower contains the model uncertainty;
Figure FDA0003743672010000033
external interference suffered by the ith follower;
(2.3) the expression of the third model is:
Figure FDA0003743672010000034
Δ iz (t)=τ iz (t)+f iz (t)+d iz (t)-ρ iz u iz (t)/m i -g
in the formula, z i (t) is the position coordinate of the ith unmanned aerial vehicle on the Z axis,
Figure FDA0003743672010000035
is z i (t) second derivative; u. of iz (t) inputting the control quantity of the ith unmanned aerial vehicle on the Z axis; delta of iz (t) is the complex uncertainty term of the third model; tau. iz (t) additive fault parameters for the ith drone; f. of iz (t) values for the system non-linear part of the model uncertainty for the ith drone; d is a radical of iz (t) external interference suffered by the ith unmanned aerial vehicle; rho iz Executing an efficiency loss fault parameter for the ith drone; g is gravity acceleration; m is a unit of i Is the quality of the ith drone.
4. The fault-tolerant control method for the air-ground heterogeneous formation system under the switching topology according to claim 3, wherein the step (3) specifically comprises:
(3.1) designing a fixed time observer aiming at the complex uncertainty of the second model, wherein the formula is as follows:
Figure FDA0003743672010000036
in the formula, σ i1 、σ i2 Is a state variable of the fixed-time observer,
Figure FDA0003743672010000037
is σ i1 The first derivative of (a) is,
Figure FDA0003743672010000038
is σ i2 The first derivative of (a); m is a unit of i1 Is a positive constant; wherein the content of the first and second substances,
Figure FDA0003743672010000039
Figure FDA00037436720100000310
Figure FDA00037436720100000311
in the formula o i1 、o i2 、o i3 Parameters of a fixed time observer;
Figure FDA0003743672010000041
a first derivative of a complex uncertainty term for the second model; w is a group of 1 Is a complex uncertainty term
Figure FDA0003743672010000042
The upper bound of (c);
Figure FDA0003743672010000043
is the first derivative of the position coordinate of the i-th follower in the XY plane;
convergence time T of observer 1i Satisfies the following conditions:
Figure FDA0003743672010000044
in the formula, v i1 Is a positive constant of a setting, and
Figure FDA0003743672010000045
then, the minimum convergence time T is obtained 1i
(3.2) designing a fixed time observer aiming at the complex uncertainty of the third model, wherein the formula is as follows:
Figure FDA0003743672010000046
in the formula, σ i3 、σ i4 Is a state variable of the fixed-time observer,
Figure FDA0003743672010000047
is σ i3 The first derivative of (a) is,
Figure FDA0003743672010000048
is σ i4 The first derivative of (a); m is i2 Is a positive constant; wherein
Figure FDA0003743672010000049
Figure FDA00037436720100000410
Figure FDA00037436720100000411
In the formula o i4 、o i5 、o i6 Parameters of a fixed time observer;
Figure FDA00037436720100000412
a first derivative of a complex uncertainty term for the third model; w 2 Is a complex uncertainty term
Figure FDA00037436720100000413
The upper bound of (c);
Figure FDA00037436720100000414
the first derivative of the position coordinate of the ith unmanned aerial vehicle on the Z axis is obtained;
convergence time T of observer 1iz Satisfies the following conditions:
Figure FDA0003743672010000051
in the formula, v i2 Is a positive constant of a setting, and when
Figure FDA0003743672010000052
Then, the minimum convergence time T is obtained 1iz
5. The fault-tolerant control method for the air-ground heterogeneous formation system under the handover topology according to claim 1, wherein the step (4) specifically comprises:
(4.1) defining an expected formation vector containing all followers in the XY plane as
Figure FDA0003743672010000053
The predetermined track of all the unmanned aerial vehicle followers on the Z axis is
Figure FDA0003743672010000054
There is a settable time constant T that is not limited to the initial state of the system max So that the following holds:
Figure FDA0003743672010000055
(4.2) based on the observed value σ i2 Laplace matrix H of communication topological relation between leader and follower σ(t) Designing a distributed fault-tolerant formation controller of an XY axis, wherein the expression is as follows:
Figure FDA0003743672010000056
Figure FDA0003743672010000057
the virtual control law is as follows:
Figure FDA0003743672010000058
Figure FDA0003743672010000059
in the formula, x 1 、δ 1 、χ 2 、δ 2 、χ 3 And delta 3 Are all positive constant, c 1 、c 2 、c 3 Are all control parameters, d 1 、d 2 、d 3 Are constants with a value greater than 1; ζ represents a unit i1 、ζ i2 、ζ i3 A third order system in which the controllers are all the same;
Figure FDA00037436720100000510
a second derivative of the position coordinates of the leader in the XY plane;
Figure FDA00037436720100000511
is gamma i2 The first derivative of (a);
Figure FDA00037436720100000512
is gamma i1 The first derivative of (a);
Figure FDA0003743672010000061
a first derivative of a convoy tracking error;
(4.3) from the observed value σ i4 Designing a Z-axis track fault-tolerant tracking controller, wherein the expression is as follows:
Figure FDA0003743672010000062
Figure FDA0003743672010000063
and the virtual control amount is designed as follows:
Figure FDA0003743672010000064
Figure FDA0003743672010000065
wherein, mu 1 、ω 1 、μ 2 、ω 2 、μ 3 And ω 3 Is a positive constant, n, set by man 1 、n 2 、n 3 Are all control parameters, m 1 、m 2 、m 3 Are constants with a value greater than 1;
Figure FDA0003743672010000066
a third order system in which the controllers are all the same;
Figure FDA0003743672010000067
is composed of
Figure FDA0003743672010000068
The first derivative of (a);
Figure FDA0003743672010000069
is composed of
Figure FDA00037436720100000610
The first derivative of (a);
Figure FDA00037436720100000611
the first derivative of the trajectory tracking error in the Z-axis.
6. The fault-tolerant control method for the air-ground heterogeneous formation system under the handover topology according to claim 5, wherein the step (4.2) specifically comprises:
(4.2.1) designing a fixed time integral sliding mode surface based on the tracking error:
Figure FDA00037436720100000612
in the formula, e i =[e ix ,e iy ] T ,sig l (e i )=[|e ix | l sgn(e ix ),|e iy | l sgn(e iy )] T Alpha and beta are both positive definite constants and satisfy 0 < alpha < 1 and beta > 1 1 、b 2 Are all constants; e.g. of the type ix For the formation tracking error of the x-axis, e iy A formation tracking error for the y-axis;
taking a first derivative and a second derivative of the sliding mode surface to obtain:
Figure FDA00037436720100000613
Figure FDA00037436720100000614
in the formula, λ i1 =diag{αb 1 |e ix | α-1 ,αb 1 |e iy | α-1i2 =diag{βb 2 |e ix | β-1 ,βb 2 |e iy | β-1 };
Obtaining a third-order system of an integral sliding mode surface:
Figure FDA0003743672010000071
Figure FDA0003743672010000072
Figure FDA0003743672010000073
(4.2.2) based on the design scheme of the fault-tolerant controller of the backstepping method, the following transformations are carried out on the three-order system:
ζ i1 =S i1
ζ i2 =S i2i1
ζ i3 =S i3i2
(4.2.3) designing a distributed fault-tolerant formation controller in an XY plane based on the established new third-order system under a switching topology network, wherein the distributed fault-tolerant formation controller realizes fixed time convergence.
7. The fault-tolerant control method for the air-ground heterogeneous formation system under the handover topology according to claim 5, wherein the step (4.3) specifically comprises:
(4.3.1) designing a fixed time integral sliding mode surface based on the tracking error:
Figure FDA0003743672010000074
in the formula (I), the compound is shown in the specification,
Figure FDA0003743672010000075
α z and beta z Are all positive definite constants; b is a mixture of 1z >0,b 2z The parameters are artificially set constants when the value is more than 0;
taking a first derivative and a second derivative of the sliding mode surface, and finally obtaining a third-order system of the integral sliding mode surface:
Figure FDA0003743672010000076
Figure FDA0003743672010000077
Figure FDA0003743672010000078
in the formula (I), the compound is shown in the specification,
Figure FDA0003743672010000079
(4.3.2) based on the fault-tolerant controller design scheme of the backstepping method, the following transformations are made to the three-order system:
Figure FDA00037436720100000710
Figure FDA00037436720100000711
Figure FDA00037436720100000712
(4.3.3) designing a track fault-tolerant tracking controller in a Z axis under a switching topology network based on the established new third-order system, wherein the track fault-tolerant tracking controller realizes fixed time convergence.
8. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the method of any one of claims 1 to 7 when executing the computer program.
9. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 7.
CN202210819674.9A 2022-07-13 2022-07-13 Fault-tolerant control method of air-ground heterogeneous formation system under switching topology Active CN115220472B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210819674.9A CN115220472B (en) 2022-07-13 2022-07-13 Fault-tolerant control method of air-ground heterogeneous formation system under switching topology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210819674.9A CN115220472B (en) 2022-07-13 2022-07-13 Fault-tolerant control method of air-ground heterogeneous formation system under switching topology

Publications (2)

Publication Number Publication Date
CN115220472A true CN115220472A (en) 2022-10-21
CN115220472B CN115220472B (en) 2023-07-21

Family

ID=83611506

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210819674.9A Active CN115220472B (en) 2022-07-13 2022-07-13 Fault-tolerant control method of air-ground heterogeneous formation system under switching topology

Country Status (1)

Country Link
CN (1) CN115220472B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116047908A (en) * 2023-01-16 2023-05-02 齐齐哈尔大学 Mixed-order heterogeneous multi-intelligent system collaborative optimal formation control method and equipment
CN116700355A (en) * 2023-08-04 2023-09-05 南京航空航天大学 Fixed wing unmanned aerial vehicle fault-tolerant control method for facing tracking of unmanned aerial vehicle

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106054922A (en) * 2016-06-22 2016-10-26 长安大学 Unmanned aerial vehicle (UAV)-unmanned ground vehicle (UGV) combined formation cooperative control method
CN110161847A (en) * 2018-09-18 2019-08-23 南京航空航天大学 A kind of unmanned plane fleet system sensor fault estimation method based on distributed unusual observer
CN111665848A (en) * 2020-06-28 2020-09-15 北京航空航天大学 Heterogeneous cluster formation tracking control method for unmanned aerial vehicle and unmanned aerial vehicle under topological switching
CN112379596A (en) * 2020-11-17 2021-02-19 南京航空航天大学 RBF neural network-based cooperative fault-tolerant control method for heterogeneous multi-agent system
CN112904723A (en) * 2021-01-19 2021-06-04 南京航空航天大学 Air-ground fixed time cooperative fault-tolerant formation control method under non-matching interference
CN113342033A (en) * 2021-05-26 2021-09-03 南京航空航天大学 Fixed-time-technology-based attitude synchronization fault-tolerant control method for fixed-wing unmanned aerial vehicle
CN114237041A (en) * 2021-11-24 2022-03-25 南京航空航天大学 Air-ground cooperative fixed time fault-tolerant control method based on preset performance
CN114594784A (en) * 2022-01-13 2022-06-07 南京航空航天大学 Neural network-based finite time fault-tolerant control method for formation of fixed-wing unmanned aerial vehicles

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106054922A (en) * 2016-06-22 2016-10-26 长安大学 Unmanned aerial vehicle (UAV)-unmanned ground vehicle (UGV) combined formation cooperative control method
CN110161847A (en) * 2018-09-18 2019-08-23 南京航空航天大学 A kind of unmanned plane fleet system sensor fault estimation method based on distributed unusual observer
CN111665848A (en) * 2020-06-28 2020-09-15 北京航空航天大学 Heterogeneous cluster formation tracking control method for unmanned aerial vehicle and unmanned aerial vehicle under topological switching
CN112379596A (en) * 2020-11-17 2021-02-19 南京航空航天大学 RBF neural network-based cooperative fault-tolerant control method for heterogeneous multi-agent system
CN112904723A (en) * 2021-01-19 2021-06-04 南京航空航天大学 Air-ground fixed time cooperative fault-tolerant formation control method under non-matching interference
CN113342033A (en) * 2021-05-26 2021-09-03 南京航空航天大学 Fixed-time-technology-based attitude synchronization fault-tolerant control method for fixed-wing unmanned aerial vehicle
CN114237041A (en) * 2021-11-24 2022-03-25 南京航空航天大学 Air-ground cooperative fixed time fault-tolerant control method based on preset performance
CN114594784A (en) * 2022-01-13 2022-06-07 南京航空航天大学 Neural network-based finite time fault-tolerant control method for formation of fixed-wing unmanned aerial vehicles

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
WANGLEI CHENG等: "Fixed-Time Fault-Tolerant Formation Control for Heterogeneous Multi-Agent Systems With Parameter Uncertainties and Disturbances", vol. 48, no. 5, pages 2121 - 2133, XP011850945, DOI: 10.1109/TCSI.2021.3061386 *
马亚杰 等: "一种无人机-无人车编队系统容错控制方法", pages 1 - 15 *
马小山 等: "考虑拓扑故障的无人机编队容错控制方法研究", vol. 38, no. 5, pages 1084 - 1093 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116047908A (en) * 2023-01-16 2023-05-02 齐齐哈尔大学 Mixed-order heterogeneous multi-intelligent system collaborative optimal formation control method and equipment
CN116700355A (en) * 2023-08-04 2023-09-05 南京航空航天大学 Fixed wing unmanned aerial vehicle fault-tolerant control method for facing tracking of unmanned aerial vehicle
CN116700355B (en) * 2023-08-04 2023-11-07 南京航空航天大学 Fixed wing unmanned aerial vehicle fault-tolerant control method for facing tracking of unmanned aerial vehicle

Also Published As

Publication number Publication date
CN115220472B (en) 2023-07-21

Similar Documents

Publication Publication Date Title
Cheng et al. Fixed-time fault-tolerant formation control for heterogeneous multi-agent systems with parameter uncertainties and disturbances
CN109324636B (en) Multi-four-rotor master-slave type cooperative formation control method based on second-order consistency and active disturbance rejection
Wang et al. Neural adaptive command filtered control for cooperative path following of multiple underactuated autonomous underwater vehicles along one path
Kuo et al. Intelligent leader-following consensus formation control using recurrent neural networks for small-size unmanned helicopters
CN115220472B (en) Fault-tolerant control method of air-ground heterogeneous formation system under switching topology
CN104865829B (en) Multi-robot system distributed self-adaption neutral net continuous tracking control method of electro
CN112904723B (en) Air-ground fixed time cooperative fault-tolerant formation control method under non-matching interference
Zhuang et al. Fixed-time formation control for spacecraft with prescribed performance guarantee under input saturation
CN112947505B (en) Multi-AUV formation distributed control method based on reinforcement learning algorithm and unknown disturbance observer
CN114237041B (en) Space-ground cooperative fixed time fault tolerance control method based on preset performance
Cheng et al. Fixed-time formation tracking for heterogeneous multiagent systems under actuator faults and directed topologies
Rekabi et al. Distributed nonlinear H∞ control algorithm for multi-agent quadrotor formation flying
CN112650290B (en) High-rise fire unmanned aerial vehicle formation optimization control method and system with disturbance compensation
Mofid et al. Adaptive integral-type terminal sliding mode control for unmanned aerial vehicle under model uncertainties and external disturbances
Liu et al. Adaptive distributed finite-time formation control for multi-UAVs under input saturation without collisions
CN112148025A (en) Unmanned aerial vehicle stability control algorithm based on integral compensation reinforcement learning
Gong et al. Distributed adaptive fault-tolerant formation control for heterogeneous multiagent systems with communication link faults
Hu et al. Formation tracking for nonlinear multi-agent systems with input and output quantization via adaptive output feedback control
Duan et al. Distributed cooperative control of multiple UAVs in the presence of actuator faults and input constraints
CN116225055A (en) Unmanned aerial vehicle autonomous flight path planning algorithm based on state decomposition in complex environment
Muñoz et al. Adaptive fuzzy consensus control strategy for UAS-based load transportation tasks
Imran et al. Adaptive control of quadrotor unmanned aerial vehicle with time-varying uncertainties
Xu et al. Distributed planar formation maneuvering of leader-follower networked systems via a barycentric coordinate-based approach
Chen et al. Adaptive optimal consensus control of multiagent systems with unknown dynamics and disturbances via reinforcement learning
Juan et al. Trajectory tracking control of multi-AUVs formation based on virtual leader

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