CN108828938B - Finite time control method of four-rotor aircraft based on inverse proportion function enhanced index approach law and fast terminal sliding mode surface - Google Patents

Finite time control method of four-rotor aircraft based on inverse proportion function enhanced index approach law and fast terminal sliding mode surface Download PDF

Info

Publication number
CN108828938B
CN108828938B CN201810519823.3A CN201810519823A CN108828938B CN 108828938 B CN108828938 B CN 108828938B CN 201810519823 A CN201810519823 A CN 201810519823A CN 108828938 B CN108828938 B CN 108828938B
Authority
CN
China
Prior art keywords
sliding mode
formula
mode surface
aircraft
coordinate system
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
CN201810519823.3A
Other languages
Chinese (zh)
Other versions
CN108828938A (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201810519823.3A priority Critical patent/CN108828938B/en
Publication of CN108828938A publication Critical patent/CN108828938A/en
Application granted granted Critical
Publication of CN108828938B publication Critical patent/CN108828938B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Abstract

A finite time control method of a four-rotor aircraft based on an inverse proportion function enhanced exponential approach law and a fast terminal sliding mode surface comprises the following steps: step 1, determining a transfer matrix from a body coordinate system based on a four-rotor aircraft to an inertial coordinate system based on the earth; step 2, analyzing a four-rotor aircraft dynamic model according to a Newton Euler formula; and 3, calculating a tracking error, and designing a controller according to the fast terminal sliding mode surface and the first derivative thereof. Aiming at a four-rotor aircraft system, the method combines sliding mode control based on inverse proportion function enhanced index approach law and rapid terminal sliding mode control, can increase approach speed when the system is far away from a sliding mode surface, can reduce buffeting, improves rapidity and robustness of the system, realizes rapid and stable control, can realize limited time control of tracking errors, and solves the problem that the tracking errors tend to 0 only when time tends to be infinite in the traditional sliding mode surface.

Description

Finite time control method of four-rotor aircraft based on inverse proportion function enhanced index approach law and fast terminal sliding mode surface
Technical Field
The invention relates to a finite time control method of a four-rotor aircraft based on an inverse proportion function enhanced exponential approximation law and a fast terminal sliding mode surface.
Background
The four-rotor aircraft has attracted wide attention of domestic and foreign scholars and scientific research institutions due to the characteristics of simple structure, strong maneuverability and unique flight mode, and is rapidly one of the hotspots of international research at present. Compared with a fixed-wing aircraft, the rotary-wing aircraft can vertically lift, has low requirement on the environment, does not need a runway, reduces the cost and has great commercial value. The development of aircrafts makes many dangerous high-altitude operations easy and safe, so as to cause deterrence to other countries in the military aspect and greatly increase the working efficiency in the civil aspect. The four-rotor aircraft has strong flexibility, can realize rapid transition of motion and hovering at any time, and can be competent for more challenging flight tasks with less damage risk. In the field of scientific research, because a four-rotor aircraft has the dynamic characteristics of nonlinearity, under-actuation and strong coupling, researchers often use the four-rotor aircraft as an experimental carrier for theoretical research and method verification. An aircraft flight control system is built by relying on a small four-rotor aircraft to carry out high-performance motion control research on the aircraft, and the method is a hot research field of the current academic world.
The approach law sliding mode control has the characteristics that discontinuous control can be realized, the sliding mode is programmable and is not related to system parameters and disturbance. The approach law sliding mode not only can reasonably design the speed of reaching the sliding mode surface, reduce the time of the approach stage, improve the robustness of the system, but also can effectively weaken the buffeting problem in the sliding mode control. Currently, in the field of four-rotor control, approach law sliding mode control is less used. The enhanced approach law further accelerates the approach speed of the system to the sliding mode surface and simultaneously enables the buffeting to be smaller on the basis of the traditional approach law.
Disclosure of Invention
In order to solve the problems that the traditional sliding mode surface can not realize limited time control, further accelerate the approaching speed of an approaching law and reduce buffeting, the invention adopts the rapid terminal sliding mode control and the enhanced index approaching law based on an inverse proportion function, avoids the singularity problem by the switching control idea, accelerates the approaching speed of a system to the sliding mode surface, reduces buffeting and realizes limited time control.
The technical scheme proposed for solving the technical problems is as follows:
a finite time control method of a four-rotor aircraft based on an inverse proportion function enhanced exponential approach law and a fast terminal sliding mode surface comprises the following steps:
step 1, determining a transfer matrix from a body coordinate system based on a four-rotor aircraft to an inertial coordinate system based on the earth;
Figure GDA0002969686990000021
wherein psi, theta and phi are respectively the yaw angle, pitch angle and roll angle of the aircraft, and represent the angle of the aircraft sequentially rotating around each axis of the inertial coordinate system, and TψTransition matrix, T, representing psiθA transition matrix, T, representing thetaφA transition matrix representing phi;
step 2, analyzing a four-rotor aircraft dynamic model according to a Newton Euler formula, wherein the process is as follows:
2.1, the translation process comprises the following steps:
Figure GDA0002969686990000022
wherein x, y and z respectively represent the position of the four rotors under an inertial coordinate system, m represents the mass of the aircraft, g represents the gravity acceleration, mg represents the gravity borne by the four rotors, and the resultant force U generated by the four rotorsr
2.2, the rotation process comprises the following steps:
Figure GDA0002969686990000023
wherein tau isx、τy、τzRespectively representing the axial moment components, I, in the coordinate system of the machine bodyxx、Iyy、IzzRespectively representing the rotational inertia component of each axis on the coordinate system of the machine body, x represents cross product, wp、wq、wrRespectively representing the attitude angular velocity components of each axis on the coordinate system of the body,
Figure GDA0002969686990000024
respectively representing the attitude angular acceleration components of all axes on the coordinate system of the machine body;
considering that the change of the attitude angle is small when the aircraft is in a low-speed flight or hovering state, the change is considered to be
Figure GDA0002969686990000039
Then the formula (3) is represented as the formula (4) in the rotation process
Figure GDA0002969686990000031
2.3, connecting the vertical type (1), (2) and (4), and obtaining the dynamic model of the aircraft as shown in the formula (5)
Figure GDA0002969686990000032
Wherein
Figure GDA0002969686990000033
Figure GDA0002969686990000034
Figure GDA0002969686990000035
Ux、Uy、UzThe input quantities of the three position controllers are respectively;
according to the formula (5), decoupling calculation is carried out on the position and posture relation, and the result is as follows:
Figure GDA0002969686990000036
wherein phidIs the desired signal value of phi, thetadDesired signal value of theta, psidFor desired signal values of ψ, the arcsin function is an arcsine function and the arctan function is an arctangent function;
equation (5) can also be written in matrix form as follows:
Figure GDA0002969686990000037
wherein X1=[x,y,z,φ,θ,ψ]T,
Figure GDA0002969686990000038
B(X)=diag(1,1,1,b1,b2,b3),U=[Ux,Uy,Uzxyz]T
Step 3, calculating a tracking error, and designing a controller according to the fast terminal sliding mode surface and the first derivative thereof, wherein the process is as follows:
3.1, defining the tracking error and its first and second differentials:
e=X1-Xd (8)
Figure GDA0002969686990000041
Figure GDA0002969686990000042
wherein, Xd=[xd,yd,zdddd]T,xd,yd,zddddConductive desired signals of x, y, z, phi, theta, psi, respectively;
3.2, designing a quick terminal sliding mode surface:
Figure GDA0002969686990000043
wherein, sigα(x)=|x|α·sign(x),α1>α2>1,λ1>0,λ2>0;
Derivation of equation (11) yields:
Figure GDA0002969686990000044
order to
Figure GDA0002969686990000045
Formula (12) is simplified to formula (13)
Figure GDA0002969686990000046
But because of
Figure GDA0002969686990000047
In existence of
Figure GDA0002969686990000048
When α (e) is 0 and β (e) is not equal to 0, the negative power term of (a) causes a singularity problem;
consider the method of handover control:
Figure GDA0002969686990000049
wherein q isi(e),αi(e),βi(e) The i-th element, i ═ 1,2,3,4,5,6, q (e), α (e), β (e), respectively;
combining formula (13) and formula (14) to obtain:
Figure GDA0002969686990000051
conjunctive formula (7), formula (10) and formula (15) yields:
Figure GDA0002969686990000052
3.3 design enhanced approach law
Figure GDA0002969686990000053
Wherein
Figure GDA0002969686990000054
N-1(X) is the inverse of N (X), k1>0,k2More than 0, more than 0 and less than 1, more than 0, more than 1, and p is a positive integer;
3.4, combined vertical (16) and formula (17), to obtain a controller
Figure GDA0002969686990000055
Wherein B is-1(X) is the inverse of B (X).
Further, the control method further includes the steps of:
step 4, property specification, the process is as follows:
4.1, proving accessibility of sliding forms:
designing Lyapunov functions
Figure GDA0002969686990000056
The derivation is performed on both sides of the function to obtain:
Figure GDA0002969686990000057
because of the scalar quantity
Figure GDA0002969686990000058
The constant is larger than 0, so the formula (18) is constantly smaller than 0, the accessibility of the sliding mode is met, and the system can reach the sliding mode surface;
4.2, enhanced effect description:
when the system moves away from the sliding mode, | s | is large, N(s) approaches δ,
Figure GDA0002969686990000059
the approach speed of the system is accelerated; when the system approaches the sliding mode, | s | approaches 0, N(s) approaches μ,
Figure GDA00029696869900000510
the buffeting of the system is reduced.
The technical conception of the invention is as follows: aiming at a four-rotor aircraft system, by combining index approach law sliding mode control and rapid terminal sliding mode control, a four-rotor aircraft finite time control method based on inverse proportion function enhanced index approach law and rapid terminal sliding mode surface is designed. The quick terminal sliding mode surface can realize the limited time control of the tracking error, and solves the problems that the time tends to be infinite and the error tends to be 0 in the traditional sliding mode surface. Based on the inverse proportion function enhanced approach law, the approach speed can be increased when the sliding mode face is far away, buffeting can be reduced, the rapidness and robustness of the system are improved, and rapid and stable control is achieved.
The invention has the beneficial effects that: compared with the traditional index approach law sliding mode control, the method can increase the approach speed when the system is far away from the sliding mode surface, reduce buffeting and shorten the arrival time of the sliding mode, thereby enabling the system to realize stable convergence more quickly. In addition, the invention utilizes the quick terminal sliding mode, solves the problems that the time tends to be infinite and the error tends to be 0 in the traditional sliding mode surface, and realizes the limited time control.
Drawings
Fig. 1 is a schematic diagram of the position tracking effect of a four-rotor aircraft, in which a dotted line represents the conventional exponential approach law control, and a dotted line represents the finite time control of the four-rotor aircraft based on the inverse proportional function enhanced exponential approach law and the fast terminal sliding mode surface.
Fig. 2 is a schematic diagram of position tracking error of a quadrotor, wherein a dotted line represents conventional exponential approximation law control, and a dotted line represents finite time control of the quadrotor based on an inverse proportional function enhanced exponential approximation law and a fast terminal sliding mode surface.
Fig. 3 is a schematic diagram of the attitude angle tracking effect of a quadrotor aircraft, wherein a dotted line represents the control of the conventional exponential approximation law, and a dotted line represents the finite time control of the quadrotor aircraft based on the inverse proportional function enhanced exponential approximation law and the fast terminal sliding mode surface.
Fig. 4 is a schematic diagram of attitude angle tracking error of a quadrotor aircraft, wherein a dotted line represents traditional exponential approach law control, and a dotted line represents limited time control of the quadrotor aircraft based on an inverse proportional function enhanced exponential approach law and a fast terminal sliding mode surface.
FIG. 5 is a schematic diagram of position controller input under finite time control of a quadrotor aircraft based on an inverse proportional function enhanced exponential approach law and a fast terminal sliding mode surface.
Fig. 6 is a schematic diagram of the position controller inputs under conventional exponential approach law control for a four-rotor aircraft.
FIG. 7 is an input schematic diagram of an attitude angle controller under finite-time control of a quadrotor aircraft based on an inverse proportional function enhanced exponential approach law and a fast terminal sliding mode surface.
FIG. 8 is a schematic diagram of attitude angle controller inputs under conventional exponential approximation law control for a quad-rotor aircraft.
FIG. 9 is a control flow diagram of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Referring to fig. 1-9, a finite-time control method of a four-rotor aircraft based on an inverse proportional function enhanced exponential approach law and a fast terminal sliding mode surface includes the following steps:
step 1, determining a transfer matrix from a body coordinate system based on a four-rotor aircraft to an inertial coordinate system based on the earth;
Figure GDA0002969686990000071
wherein psi, theta and phi are respectively the yaw angle, pitch angle and roll angle of the aircraft, and represent the angle of the aircraft sequentially rotating around each axis of the inertial coordinate system, and TψTransition matrix, T, representing psiθA transition matrix, T, representing thetaφA transition matrix representing phi;
step 2, analyzing a four-rotor aircraft dynamic model according to a Newton Euler formula, wherein the process is as follows:
2.1, the translation process comprises the following steps:
Figure GDA0002969686990000072
wherein x, y and z respectively represent the position of the four rotors under an inertial coordinate system, m represents the mass of the aircraft, g represents the gravity acceleration, mg represents the gravity borne by the four rotors, and the resultant force U generated by the four rotorsr
2.2, the rotation process comprises the following steps:
Figure GDA0002969686990000081
wherein tau isx、τy、τzRespectively representing the axial moment components, I, in the coordinate system of the machine bodyxx、Iyy、IzzRespectively representing the rotational inertia component of each axis on the coordinate system of the machine body, x represents cross product, wp、wq、wrRespectively representing the attitude angular velocity components of each axis on the coordinate system of the body,
Figure GDA0002969686990000082
respectively representing the attitude angular acceleration components of all axes on the coordinate system of the machine body;
considering that the change of the attitude angle is small when the aircraft is in a low-speed flight or hovering state, the change is considered to be
Figure GDA0002969686990000083
Then the rotation process is Chinese (3 can be expressed as formula (4)
Figure GDA0002969686990000084
2.3, connecting the vertical type (1), (2) and (4), and obtaining the dynamic model of the aircraft as shown in the formula (5)
Figure GDA0002969686990000085
Wherein
Figure GDA0002969686990000086
Figure GDA0002969686990000087
Figure GDA0002969686990000088
Ux、Uy、UzThe input quantities of the three position controllers are respectively;
according to the formula (5), decoupling calculation is carried out on the position and posture relation, and the result is as follows:
Figure GDA0002969686990000091
wherein phidIs the desired signal value of phi, thetadDesired signal value of theta, psidFor desired signal values of ψ, the arcsin function is an arcsine function and the arctan function is an arctangent function;
equation (5) can also be written in matrix form as follows:
Figure GDA0002969686990000092
wherein
Figure GDA0002969686990000093
B(X)=diag(1,1,1,b1,b2,b3),U=[Ux,Uy,Uzxyz]T
Step 3, calculating a tracking error, and designing a controller according to the fast terminal sliding mode surface and the first derivative thereof, wherein the process is as follows:
3.1, defining the tracking error and its first and second differentials:
e=X1-Xd (8)
Figure GDA0002969686990000094
Figure GDA0002969686990000095
wherein, Xd=[xd,yd,zdddd]T,xd,yd,zddddConductive desired signals of x, y, z, phi, theta, psi, respectively;
3.2, designing a quick terminal sliding mode surface:
Figure GDA0002969686990000096
wherein, sigα(x)=|x|α·sign(x),α1>α2>1,λ1>0,λ2>0;
Derivation of equation (11) yields:
Figure GDA0002969686990000097
order to
Figure GDA0002969686990000101
Formula (12) is simplified to formula (13)
Figure GDA0002969686990000102
But because of
Figure GDA0002969686990000103
In existence of
Figure GDA0002969686990000104
When α (e) is 0 and β (e) is not equal to 0, the negative power term of (a) causes a singularity problem;
consider the method of handover control:
Figure GDA0002969686990000105
wherein q isi(e),αi(e),βi(e) The i-th element, i ═ 1,2,3,4,5,6, q (e), α (e), β (e), respectively;
combining formula (13) and formula (14) to obtain:
Figure GDA0002969686990000106
conjunctive formula (7), formula (10) and formula (15) yields:
Figure GDA0002969686990000107
3.3 design enhanced approach law
Figure GDA0002969686990000108
Wherein
Figure GDA0002969686990000109
N-1(X) is the inverse of N (X), k1>0,k2More than 0, more than 0 and less than 1, more than 0, more than 1, and p is a positive integer;
3.4, combined vertical (16) and formula (17), to obtain a controller
Figure GDA00029696869900001010
Wherein B is-1(X) is the inverse of B (X);
step 4, property specification, the process is as follows:
4.1, proving accessibility of sliding forms:
designing Lyapunov functions
Figure GDA0002969686990000111
The derivation is performed on both sides of the function to obtain:
Figure GDA0002969686990000112
because of the scalar quantity
Figure GDA0002969686990000113
The constant is larger than 0, so the formula (18) is constantly smaller than 0, the accessibility of the sliding mode is met, and the system can reach the sliding mode surface;
4.2, enhanced effect description:
when the system moves away from the sliding mode, | s | is large, N(s) approaches δ,
Figure GDA0002969686990000114
the approach speed of the system is accelerated; when the system approaches the sliding mode, | s | approaches 0, N(s) approaches μ,
Figure GDA0002969686990000115
the buffeting of the system is reduced.
In order to verify the effectiveness of the method, the invention provides a contrast between an enhanced index approximation law sliding mode control method based on an inverse proportion function and a traditional index approximation law sliding mode control method:
for more efficient comparison, all parameters of the system are consistent, i.e. xd=yd=zd=20、ψd0.5, slip form surface parameters: lambda [ alpha ]1=0.2、λ2=0.7、α1=2、α21.1, 0.1, and the approach law parameter: k is a radical of1=0.6、k20.8, δ -0.5, p-1, γ -1, μ -2, quad-rotor aircraft parameters: m 0.625, L0.1275, Ixx=2.3×10-3、Iyy=2.4×10-3、Izz=2.6×10-3G ═ 10, sampling parameters: t is ts=0.007,N=5000。
1-4, the finite-time control of the quadrotor aircraft based on the inverse proportion function enhanced exponential approach law and the fast terminal sliding mode surface can reach the expected position more quickly; with reference to fig. 5-8, the limited-time control of the quadrotor based on the inverse proportional function enhanced exponential approximation law and the fast terminal sliding mode surface has smaller buffeting.
In conclusion, the finite time control of the four-rotor aircraft based on the inverse proportion function enhanced index approach law and the fast terminal sliding mode surface can reduce the buffeting and the tracking time at the same time, improve the tracking performance and enable the system to enter stable convergence more quickly.
While the foregoing has described a preferred embodiment of the invention, it will be appreciated that the invention is not limited to the embodiment described, but is capable of numerous modifications without departing from the basic spirit and scope of the invention as set out in the appended claims.

Claims (2)

1. A finite time control method of a four-rotor aircraft based on an inverse proportion function enhanced exponential approach law and a fast terminal sliding mode surface is characterized by comprising the following steps:
step 1, determining a transfer matrix from a body coordinate system based on a four-rotor aircraft to an inertial coordinate system based on the earth;
Figure FDA0002969686980000011
wherein psi, theta and phi are respectively the yaw angle, pitch angle and roll angle of the aircraft, and represent the angle of the aircraft sequentially rotating around each axis of the inertial coordinate system, and TψTransition matrix, T, representing psiθA transition matrix, T, representing thetaφA transition matrix representing phi;
step 2, analyzing a four-rotor aircraft dynamic model according to a Newton Euler formula, wherein the process is as follows:
2.1, the translation process comprises the following steps:
Figure FDA0002969686980000012
wherein x, y and z respectively represent the position of the four rotors under an inertial coordinate system, m represents the mass of the aircraft, g represents the gravity acceleration, mg represents the gravity borne by the four rotors, and the resultant force U generated by the four rotorsr
2.2, the rotation process comprises the following steps:
Figure FDA0002969686980000013
wherein tau isx、τy、τzRespectively representing the axial moment components, I, in the coordinate system of the machine bodyxx、Iyy、IzzRespectively representing the rotational inertia component of each axis on the coordinate system of the machine body, x represents cross product, wp、wq、wrRespectively representing the attitude angular velocity components of each axis on the coordinate system of the body,
Figure FDA0002969686980000014
respectively representing the attitude angular acceleration components of all axes on the coordinate system of the machine body;
considering that the aircraft is in a low-speed flight or hovering state, consider
Figure FDA0002969686980000021
Then the formula (3) is represented as the formula (4) in the rotation process
Figure FDA0002969686980000022
2.3, connecting the vertical type (1), (2) and (4), and obtaining the dynamic model of the aircraft as shown in the formula (5)
Figure FDA0002969686980000023
Wherein
Figure FDA0002969686980000024
Figure FDA0002969686980000025
Figure FDA0002969686980000026
Ux、Uy、UzThe input quantities of the three position controllers are respectively;
according to the formula (5), decoupling calculation is carried out on the position and posture relation, and the result is as follows:
Figure FDA0002969686980000027
wherein phidIs the desired signal value of phi, thetadDesired signal value of theta, psidFor desired signal values of ψ, the arcsin function is an arcsine function and the arctan function is an arctangent function;
equation (5) can also be written in matrix form as follows:
Figure FDA0002969686980000028
wherein X1=[x,y,z,φ,θ,ψ]T,
Figure FDA0002969686980000031
B(X)=diag(1,1,1,b1,b2,b3),U=[Ux,Uy,Uzxyz]T
Step 3, calculating a tracking error, and designing a controller according to the fast terminal sliding mode surface and the first derivative thereof, wherein the process is as follows:
3.1, defining the tracking error and its first and second differentials:
e=X1-Xd (8)
Figure FDA0002969686980000032
Figure FDA0002969686980000033
wherein, Xd=[xd,yd,zdddd]T,xd,yd,zddddConductive desired signals of x, y, z, phi, theta, psi, respectively;
3.2, designing a quick terminal sliding mode surface:
Figure FDA0002969686980000034
wherein, sigα(x)=|x|α·sign(x),α1>α2>1,λ1>0,λ2>0;
Derivation of equation (11) yields:
Figure FDA0002969686980000035
order to
Figure FDA0002969686980000036
Formula (12) is simplified to formula (13)
Figure FDA0002969686980000037
But because of
Figure FDA0002969686980000038
In existence of
Figure FDA0002969686980000039
When α (e) is 0 and β (e) is not equal to 0, the negative power term of (a) causes a singularity problem;
consider the method of handover control:
Figure FDA0002969686980000041
wherein q isi(e),αi(e),βi(e) The i-th element, i ═ 1,2,3,4,5,6, q (e), α (e), β (e), respectively;
combining formula (13) and formula (14) to obtain:
Figure FDA0002969686980000042
conjunctive formula (7), formula (10) and formula (15) yields:
Figure FDA0002969686980000043
3.3 design enhanced approach law
Figure FDA0002969686980000044
Wherein N(s) is defined as
Figure FDA0002969686980000045
N-1(X) is the inverse of N (X), k1>0,k2More than 0, more than 0 and less than 1, more than 0, more than 1, and p is a positive integer;
3.4, combined vertical (16) and formula (17), to obtain a controller
Figure FDA0002969686980000046
Wherein B is-1(X) is the inverse of B (X).
2. The limited-time control method for a quadrotor aircraft based on an inverse proportional function enhanced exponential approach law and a fast terminal sliding-mode surface according to claim 1, characterized in that the control method further comprises the following steps:
step 4, property specification, the process is as follows:
4.1, proving accessibility of sliding forms:
designing Lyapunov functions
Figure FDA0002969686980000047
The derivation is performed on both sides of the function to obtain:
Figure FDA0002969686980000048
because of the fact that
Figure FDA0002969686980000051
The constant is larger than 0, so the formula (18) is constantly smaller than 0, the accessibility of the sliding mode is met, and the system can reach the sliding mode surface;
4.2, enhanced effect description:
when the system moves away from the sliding mode, | s | is large, N(s) approaches δ,
Figure FDA0002969686980000052
the approach speed of the system is accelerated; when the system approaches the sliding mode, | s | approaches 0, N(s) approaches μ,
Figure FDA0002969686980000053
the buffeting of the system is reduced.
CN201810519823.3A 2018-05-28 2018-05-28 Finite time control method of four-rotor aircraft based on inverse proportion function enhanced index approach law and fast terminal sliding mode surface Active CN108828938B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810519823.3A CN108828938B (en) 2018-05-28 2018-05-28 Finite time control method of four-rotor aircraft based on inverse proportion function enhanced index approach law and fast terminal sliding mode surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810519823.3A CN108828938B (en) 2018-05-28 2018-05-28 Finite time control method of four-rotor aircraft based on inverse proportion function enhanced index approach law and fast terminal sliding mode surface

Publications (2)

Publication Number Publication Date
CN108828938A CN108828938A (en) 2018-11-16
CN108828938B true CN108828938B (en) 2021-06-18

Family

ID=64145794

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810519823.3A Active CN108828938B (en) 2018-05-28 2018-05-28 Finite time control method of four-rotor aircraft based on inverse proportion function enhanced index approach law and fast terminal sliding mode surface

Country Status (1)

Country Link
CN (1) CN108828938B (en)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102436671A (en) * 2011-08-16 2012-05-02 上海交通大学 Virtual viewpoint drawing method based on depth value non-linear transformation
WO2013021453A1 (en) * 2011-08-09 2013-02-14 トヨタ自動車株式会社 Sliding mode controller, and internal combustion engine system control apparatus
CN103197551A (en) * 2013-03-12 2013-07-10 南京航空航天大学 One-way sliding-mode control method of single rotor wing unmanned aerial vehicle
CN103217981A (en) * 2013-04-18 2013-07-24 东南大学 Four-rotor aircraft speed control method based on integral variable structure control
CN105676641A (en) * 2016-01-25 2016-06-15 南京航空航天大学 Nonlinear robust controller design method based on back-stepping and sliding mode control technologies and aimed at nonlinear model of quad-rotor unmanned plane
CN105759832A (en) * 2016-05-20 2016-07-13 武汉科技大学 Four-rotor aircraft sliding mode variable structure control method based on inversion method
CN105786021A (en) * 2016-05-27 2016-07-20 天津工业大学 Quadrotor unmanned aerial vehicle control method
CN105911866A (en) * 2016-06-15 2016-08-31 浙江工业大学 Finite-time full-order sliding mode control method of quadrotor unmanned aircraft
CN106886149A (en) * 2017-02-23 2017-06-23 哈尔滨工业大学 A kind of spacecraft robust finite time saturation Attitude tracking control method
CN107490966A (en) * 2017-08-22 2017-12-19 浙江工业大学 A kind of aircraft finite time Adaptive Attitude control method based on improvement power Reaching Law
CN107688295A (en) * 2017-08-29 2018-02-13 浙江工业大学 A kind of quadrotor finite time self-adaptation control method based on fast terminal sliding formwork
CN107817818A (en) * 2017-12-07 2018-03-20 中国人民解放军国防科技大学 Finite time control method for flight path tracking of uncertain model airship
CN107977011A (en) * 2017-12-26 2018-05-01 电子科技大学 Quadrotor UAV Flight Control method based on Fractional Control Algorithm

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7453227B2 (en) * 2005-12-20 2008-11-18 Intuitive Surgical, Inc. Medical robotic system with sliding mode control

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013021453A1 (en) * 2011-08-09 2013-02-14 トヨタ自動車株式会社 Sliding mode controller, and internal combustion engine system control apparatus
CN102436671A (en) * 2011-08-16 2012-05-02 上海交通大学 Virtual viewpoint drawing method based on depth value non-linear transformation
CN103197551A (en) * 2013-03-12 2013-07-10 南京航空航天大学 One-way sliding-mode control method of single rotor wing unmanned aerial vehicle
CN103217981A (en) * 2013-04-18 2013-07-24 东南大学 Four-rotor aircraft speed control method based on integral variable structure control
CN105676641A (en) * 2016-01-25 2016-06-15 南京航空航天大学 Nonlinear robust controller design method based on back-stepping and sliding mode control technologies and aimed at nonlinear model of quad-rotor unmanned plane
CN105759832A (en) * 2016-05-20 2016-07-13 武汉科技大学 Four-rotor aircraft sliding mode variable structure control method based on inversion method
CN105786021A (en) * 2016-05-27 2016-07-20 天津工业大学 Quadrotor unmanned aerial vehicle control method
CN105911866A (en) * 2016-06-15 2016-08-31 浙江工业大学 Finite-time full-order sliding mode control method of quadrotor unmanned aircraft
CN106886149A (en) * 2017-02-23 2017-06-23 哈尔滨工业大学 A kind of spacecraft robust finite time saturation Attitude tracking control method
CN107490966A (en) * 2017-08-22 2017-12-19 浙江工业大学 A kind of aircraft finite time Adaptive Attitude control method based on improvement power Reaching Law
CN107688295A (en) * 2017-08-29 2018-02-13 浙江工业大学 A kind of quadrotor finite time self-adaptation control method based on fast terminal sliding formwork
CN107817818A (en) * 2017-12-07 2018-03-20 中国人民解放军国防科技大学 Finite time control method for flight path tracking of uncertain model airship
CN107977011A (en) * 2017-12-26 2018-05-01 电子科技大学 Quadrotor UAV Flight Control method based on Fractional Control Algorithm

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A novel reaching law for smooth sliding mode control using inverse hyperbolic function;Muhammad Asad等;《2012 International Conference on Emerging Technologies》;20121231;1-6 *
Li Xing等.Model-free adaptive sliding mode control method for uncertain robot system.《Journal of Nanjing University of Science and Technology》.2015,第39卷(第6期),655-60. *
基于扩张状态观测器的永磁同步电机滑模变结构位置伺服控制;陈强等;《新型工业化》;20151231;第5卷(第8期);17-25 *
小型四旋翼无人飞行器非线性控制方法研究;宋占魁;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20150715(第07期);C031-42 *
近空间飞行器鲁棒自适应滑模控制;蒲明;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20131015(第10期);C031-12 *

Also Published As

Publication number Publication date
CN108828938A (en) 2018-11-16

Similar Documents

Publication Publication Date Title
CN107976902B (en) Enhanced constant-speed approach law sliding mode control method of quad-rotor unmanned aerial vehicle system
CN107976903B (en) Enhanced double-power-order approach law sliding mode control method of quad-rotor unmanned aerial vehicle system
CN107957682B (en) Enhanced fast power-order approach law sliding mode control method of quad-rotor unmanned aerial vehicle system
CN108828937B (en) Finite time control method of four-rotor aircraft based on exponential enhancement type exponential approaching law and fast terminal sliding mode surface
CN108845497B (en) Finite time control method of four-rotor aircraft based on hyperbolic tangent enhanced index approach law and fast terminal sliding mode surface
CN108153148B (en) Enhanced index approach law sliding mode control method of quad-rotor unmanned aerial vehicle system
CN108646773B (en) Self-adaptive control method of four-rotor aircraft based on exponential enhancement type double-power approach law and fast terminal sliding mode surface
CN108563128B (en) Self-adaptive control method of four-rotor aircraft based on exponential enhancement type rapid power approximation law and rapid terminal sliding mode surface
CN108536018B (en) Four-rotor aircraft self-adaptive control method based on inverse proportion function enhanced double-power approach law and fast terminal sliding mode surface
CN108549401B (en) Finite time control method of four-rotor aircraft based on hyperbolic sine enhanced index approach law and fast terminal sliding mode surface
CN108829128B (en) Four-rotor aircraft finite time control method based on logarithm enhancement type exponential approaching law and fast terminal sliding mode surface
CN108762076B (en) Finite time control method of four-rotor aircraft based on inverse proportion function enhanced constant speed approach law and rapid terminal sliding mode surface
CN108845586B (en) Finite time control method of four-rotor aircraft based on hyperbolic sine enhanced constant-speed approach law and fast terminal sliding mode surface
CN108563126B (en) Self-adaptive control method of four-rotor aircraft based on hyperbolic sine enhanced power approximation law and fast terminal sliding mode surface
CN108563125B (en) Self-adaptive control method of four-rotor aircraft based on exponential enhancement type power approach law and fast terminal sliding mode surface
CN108803638B (en) Self-adaptive control method of four-rotor aircraft based on hyperbolic tangent enhanced rapid power approach law and rapid terminal sliding mode surface
CN108829117B (en) Self-adaptive control method of four-rotor aircraft based on logarithm enhancement type power approach law and fast terminal sliding mode surface
CN108549241B (en) Self-adaptive control method of four-rotor aircraft based on arc tangent enhanced double-power approach law and fast terminal sliding mode surface
CN108803319B (en) Self-adaptive control method of four-rotor aircraft based on logarithm enhancement type fast power approach law and fast terminal sliding mode surface
CN108549400B (en) Self-adaptive control method of four-rotor aircraft based on logarithm enhanced double-power approach law and fast terminal sliding mode surface
CN108828938B (en) Finite time control method of four-rotor aircraft based on inverse proportion function enhanced index approach law and fast terminal sliding mode surface
CN108628332B (en) Finite time control method of four-rotor aircraft based on arc tangent enhanced index approach law and fast terminal sliding mode surface
CN108803320B (en) Finite time control method of four-rotor aircraft based on exponential enhancement type constant velocity approach law and rapid terminal sliding mode surface
CN108762075B (en) Four-rotor aircraft finite time control method based on logarithm enhanced constant-speed approach law and fast terminal sliding mode surface
CN108828936B (en) Finite time control method of four-rotor aircraft based on arc tangent enhanced constant velocity approach law and fast terminal sliding mode surface

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20181116

Assignee: Foshan chopsticks Technology Co.,Ltd.

Assignor: JIANG University OF TECHNOLOGY

Contract record no.: X2024980000084

Denomination of invention: Finite time control method for quadcopter based on inverse proportional function enhanced exponential approaching law and fast terminal sliding mode surface

Granted publication date: 20210618

License type: Common License

Record date: 20240104

EE01 Entry into force of recordation of patent licensing contract