CN112666960B - L1 augmentation self-adaption-based control method for rotary wing aircraft - Google Patents

L1 augmentation self-adaption-based control method for rotary wing aircraft Download PDF

Info

Publication number
CN112666960B
CN112666960B CN202011466461.XA CN202011466461A CN112666960B CN 112666960 B CN112666960 B CN 112666960B CN 202011466461 A CN202011466461 A CN 202011466461A CN 112666960 B CN112666960 B CN 112666960B
Authority
CN
China
Prior art keywords
control
representing
rate
loop
self
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
CN202011466461.XA
Other languages
Chinese (zh)
Other versions
CN112666960A (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.)
Xi'an Innno Aviation Technology Co ltd
Northwestern Polytechnical University
Original Assignee
Xi'an Innno Aviation Technology Co ltd
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xi'an Innno Aviation Technology Co ltd, Northwestern Polytechnical University filed Critical Xi'an Innno Aviation Technology Co ltd
Priority to CN202011466461.XA priority Critical patent/CN112666960B/en
Publication of CN112666960A publication Critical patent/CN112666960A/en
Application granted granted Critical
Publication of CN112666960B publication Critical patent/CN112666960B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Feedback Control In General (AREA)

Abstract

The invention relates to a rotor craft control method based on L1 augmentation self-adaption, which comprises the steps of firstly establishing a rotor craft kinematics and dynamics model, designing cascade PID control rate and L1 augmentation self-adaption control rate, and obtaining total control input of the rotor craft. The method can estimate interference on line in real time, and can effectively realize stable control of the rotor craft and large disturbance by matching with the traditional cascade PID control. Practical engineering application shows that the L1-based augmentation self-adaption combined with the traditional cascade PID control can be easily realized in engineering, and the mode has better robustness, stronger disturbance rejection capability and better self-adaption effect than the traditional cascade PID control.

Description

L1 augmentation self-adaption-based control method for rotary wing aircraft
Technical Field
The invention belongs to the technical field of unmanned aerial vehicle application, relates to a control method of an unmanned aerial vehicle, and particularly relates to a control method of a rotor craft based on L1 augmentation self-adaption.
Background
With the development of science and technology, rotorcraft have been used in many military fields, such as reconnaissance of battlefield environments, and striking of small objects, etc. In the civil field, the rotor craft plays an important role in aerial photography, security, electric power inspection, oil and gas pipeline inspection and the like. In the actual flight process of the rotorcraft, the rotation speed of the rotor motor is controlled to achieve the desired position and attitude, and the control of the rotorcraft is the core of the rotorcraft capable of achieving the flight, so that the rotorcraft is particularly important. The control quality of the rotorcraft is directly related to the completion degree of the operation task, and is the core technical basis of the operation. Because the rotor craft needs to carry out the operation flow under complicated operating mode, need rotor craft can adapt to different operating modes to can all have fine control quality under various different operating modes. The traditional PID control has wide application because of no need of model, but because rotor flight needs to run in complex severe environment, the traditional PID control is difficult to adapt to complex working condition environments, such as windy weather, or sudden release of load, and the like, under the conditions, the traditional PID control is difficult to have very good control quality or one group of controller parameters are difficult to adapt to a plurality of working conditions, and the design structure of the control rate of the traditional PID+L1 augmentation adaptive control is proposed herein, can adapt to complex working conditions, and does not need to set a plurality of groups of parameters, and has very good control quality under each complex working condition. The cascade PID+L1 amplification self-adaptive control structure is provided for enhancing the self-adaptability and the robustness of the PID system so as to adapt to different working conditions, and the structure can adapt to the working conditions wider than a single cascade PID without resetting PID parameters to adapt to different working conditions like a single cascade PID, and can obtain better control quality.
Disclosure of Invention
Technical problem to be solved
In order to avoid the defects of the prior art, the invention provides a rotor craft control method based on L1 augmentation self-adaption.
Technical proposal
A rotor craft control method based on L1 augmentation self-adaption is characterized by comprising the following steps:
step 1, establishing a rotorcraft kinematics and dynamics model:
rotorcraft kinematic model:
Figure BDA0002832663990000021
wherein: m represents the weight of the unmanned aerial vehicle, g represents the gravitational acceleration, U 1 Represent lift, K d Represents the coefficient of resistance and,
Figure BDA0002832663990000022
the speeds in the three directions are indicated,
Figure BDA0002832663990000023
indicating acceleration in three directions. />
Figure BDA0002832663990000024
The machine system rotates the matrix to the navigation system;
Figure BDA0002832663990000025
wherein: ψ denotes a heading angle, θ denotes a pitch angle, γ denotes a roll angle,
Figure BDA0002832663990000026
a rotation matrix representing the machine system to the navigation system;
kinetic model of rotorcraft:
Figure BDA0002832663990000031
Figure BDA0002832663990000032
representing the derivative of the attitude angle. cos denotes the cosine, sin the sine,
Figure BDA0002832663990000033
Figure BDA0002832663990000034
Figure BDA0002832663990000035
Figure BDA0002832663990000036
Figure BDA0002832663990000037
triaxial angular acceleration, +.>
Figure BDA0002832663990000038
Representing the triaxial angular rate. I xx ,I yy ,I zz Representing the triaxial inertia, L representing the distance from the rotor centre to the aircraft centre of mass; u (U) 1 Representing lift force, U 2 Indicating torque in rolling direction, U 3 Indicating torque in pitch direction, U 4 Torque representing heading direction;
step 2, designing cascade PID control rate:
horizontal channel position loop control rate design:
v t =K p ·(p t -p)
wherein: p is p t Is the desired location; p is the current position, v t To a desired speed, K p The ratio gain is the position loop;
horizontal channel speed loop control rate design:
a t =K v ·(v t -v)
wherein: v t Is the desired speed; v is the current speed, a t For the desired acceleration, ψ is the mean heading angle,θ is pitch angle, γ is roll angle, K v Gain for the speed loop ratio;
Figure BDA0002832663990000041
θ t =arctan(-a x /g)
γ t =arctan(a y ·cosθ/g)
wherein: θ t Is a desired pose; θ is the current pose, a y ,a x Representing a horizontal acceleration;
angle ring control rate design:
ω t =K θ ·(α t -α)
ω t is the desired angular rate; omega is the current angular rate, K θ For the attitude loop proportional gain, alpha represents the dip angle, i.e. roll and pitch;
angular rate ring control rate design:
Δω=ω t
Figure BDA0002832663990000042
wherein: k (K) ω For angular rate loop proportional gain, K I For angular rate loop integral gain, K D For the angular rate loop differential gain,
Figure BDA0002832663990000043
representing an integral operation, s representing a differential operation; Δω is the target angular rate ω t Difference from the current angular rate ω +.>
Figure BDA0002832663990000044
S is a differential operation symbol;
step 3, L1 augmentation adaptive control rate:
the L1 augmented adaptive state prediction equation is:
Figure BDA0002832663990000051
wherein: a state matrix A, a control matrix B, a matrix A SP A 3 x 3 matrix with a determinant of negative values, T representing the period of operation of the controller;
Figure BDA0002832663990000052
wherein:
Figure BDA0002832663990000058
representing status error, ++>
Figure BDA0002832663990000054
Representing an estimated state, x (t) representing a true state;
adaptive rate of L1 adaptation:
Figure BDA0002832663990000055
output of the adaptive controller:
Figure BDA0002832663990000056
wherein: the expression of C(s) is
Figure BDA0002832663990000057
Step 4, general control input of the rotorcraft:
control output u of step 2 b (t) output u of step 3 a (t) adding to obtain the total control output u (t)
u(t)=u b (t)+u a (t)。
Advantageous effects
According to the L1 augmentation self-adaptive rotor craft control method, interference can be estimated on line in real time, and the method is matched with the traditional cascade PID control, so that stable control of the rotor craft can be effectively achieved, and stable control of large disturbance can be achieved. Practical engineering application shows that the L1-based augmentation self-adaption combined with the traditional cascade PID control can be easily realized in engineering, and the mode has better robustness, stronger disturbance rejection capability and better self-adaption effect than the traditional cascade PID control.
The invention has the beneficial effects that:
(1) Compared with cascade PID, the anti-interference capability is stronger;
(2) Compared with cascade PID, the method has stronger adaptability;
(3) Compared with cascade PID, the method has stronger robustness;
(4) Compared with other modern control theory, the method is easier to realize in engineering, and does not need to finely identify the system model;
the L1 augmentation self-adaptive algorithm is combined with cascade PID, so that the adaptability and the robustness of the system can be enhanced, and the control performance of the rotorcraft is greatly improved.
Drawings
FIG. 1 is a horizontal channel control structure;
FIG. 2 height control structure;
FIG. 3 L1 augments an adaptive control architecture;
fig. 4 basic control rate+l1 augmented adaptive control architecture.
Detailed Description
The invention will now be further described with reference to examples, figures:
the first step, a dynamic and kinematic model of the rotorcraft is given for the rotorcraft as a controlled object; secondly, aiming at a model of the rotorcraft, a control strategy of a traditional cascade PID controller is provided, and the output of an angular rate ring in the PID controller is taken as a basic controller u b (output of the base controller); thirdly, we provide a specific implementation method of the L1 controller, and record the output of the L1 adaptive controller as u a (adaptive controller)An output); fourth, the control input of the model of the rotorcraft is u total The concrete representation method is as follows
u total =u a +u b (1)
The steps can be written as
(1) Establishing a six-degree-of-freedom model of the rotor craft;
(2) Aiming at a six-degree-of-freedom model of the rotor craft, a traditional cascade PID control strategy is provided;
(3) Aiming at a six-degree-of-freedom model of the rotorcraft, a specific implementation mode of L1 augmentation self-adaptive control is provided.
(4) And calculating the total control input of the rotorcraft by combining the control output of the cascade PID and the L1 augmentation self-adaptive control output.
The individual implementation steps are described in detail below.
Step 1, establishing a rotorcraft kinematics and dynamics model
The method only describes a traditional PID and L1 augmentation self-adaptive specific control method, does not describe a complex and fine model of the rotorcraft system, and only adopts general system description, for example, engineering researchers want to know the system model of the rotorcraft more deeply, and can refer to related literature data or books for further study.
Taking a four-rotor aircraft as an example, before deriving a six-degree-of-freedom model of the rotor aircraft, it is necessary to perform variable description on some variables that are used later. Definitions for variables to be used later
Table 1 shows the results.
TABLE 1 definition of physical variables
Figure BDA0002832663990000071
Figure BDA0002832663990000081
(1) Establishing a rotorcraft kinematics equation
Rotorcraft are subjected mainly to the following forces: gravity; lifting force; resistance.
According to the aerodynamic theory of the rotorcraft, the lift force U 1 Torsion in three directions (U 2 ,U 3 ,U 4 ) The specific calculation method of (1) is related to the rotating speed of the rotorcraft, and the expression is as follows:
Figure BDA0002832663990000082
Figure BDA0002832663990000083
Figure BDA0002832663990000084
Figure BDA0002832663990000085
the meanings of the respective letters in the formulae are shown in Table 1. According to the stress condition of the rotor craft, the kinematic equation of the rotor craft is obtained by applying Newton's second law:
Figure BDA0002832663990000091
the third term on the right of the above medium number is a drag term, which is proportional to the speed of the aircraft and has opposite signs.
Wherein the matrix is rotated
Figure BDA0002832663990000092
Is represented by the expression:
Figure BDA0002832663990000093
the ψ in the formula (7) represents a heading angle, θ represents a pitch angle, and γ represents a roll angle.
(2) Establishing a kinetic equation for a rotorcraft
The conversion relationship between the angular rate of the rotorcraft and the derivative of the euler angle is as follows:
Figure BDA0002832663990000094
the following equation holds according to the law of conservation of angular momentum:
Figure BDA0002832663990000095
the above detailed development is as follows:
Figure BDA0002832663990000096
then the kinetic equation for the rotorcraft can be derived.
Figure BDA0002832663990000101
Equations (7) and (11) describe the translational and rotational movement of the rotorcraft in space for the kinematic and kinetic model equations of the rotorcraft.
Step 2, design of cascade PID control rate
For the dynamics and kinematics model of the rotorcraft, the rotorcraft can be controlled by using a traditional cascade PID control mode, and the following description is divided into the following parts:
(1) horizontal position channel control structure
The control structure of the horizontal position channel is shown in fig. 1, and variables used in the figure are as follows:
(1)p t a desired location; p represents the current position, and the onboard integrated navigation systemCalculating to obtain;
(2)v t a desired speed; v represents the current speed and is calculated by an onboard integrated navigation system;
(3)a t a desired acceleration; a represents the current acceleration, and can be calculated by an organic vehicle-mounted integrated navigation system; see fig. 1;
(4)θ t a desired pose; θ represents the current gesture, which can be calculated by the organic-vehicle integrated navigation system;
(5)ω t a desired angular rate; ω represents the current angular rate, which can be obtained by an on-board integrated navigation system.
Desired position p t The expected speed v is obtained by making a difference from the current position p through a position proportion controller t The method comprises the steps of carrying out a first treatment on the surface of the Desired speed v t The expected acceleration a is obtained by making a difference with the current speed v through a speed proportional controller t . The above position loop and speed loop control rates can be described by the following formulas:
horizontal channel position loop control rate design:
v t =K p ·(p t -p) (12)
wherein K is p Is the position loop proportional gain.
Horizontal channel speed loop control rate design:
a t =K v ·(v t -v) (13)
wherein K is v Proportional gain for speed loop
Desired acceleration a t Theoretical calculation of the desired attitude θ by small disturbance hypothesis t The calculation formula for the conversion from the desired acceleration into the desired attitude is as follows, the condition of the small disturbance assumption being that the rotorcraft has no altitude-directional movement in the equilibrium position, the body has no heading-directional movement and the heading angle ψ≡0 is zero.
Figure BDA0002832663990000111
Is arranged to obtain
Figure BDA0002832663990000112
And (3) finishing the formula (15) again to obtain a calculation formula from acceleration to inclination angle:
Figure BDA0002832663990000113
equation (16) is the equation from acceleration to tilt angle (roll and pitch). According to (16) the desired acceleration a t Converted into the desired inclination angle theta t . Desired inclination angle theta t The error angle is obtained by making a difference with the current inclination angle, and the expected angular rate omega is obtained by the error angle through the gesture proportional controller t . The attitude loop control rate design can be expressed by the following expression:
ω t =K θ ·(α t -α) (17)
wherein K is θ For attitude loop proportional gain, α here represents tilt angle (roll and pitch).
Desired angular rate omega t The angular output of the angular rate controller obtained by the angular rate PID controller after the difference is made with the current measured angular rate is recorded as u b (u b Referred to as the base controller output, which is used in the subsequent steps), the horizontal position control architecture shown in fig. 1 describes how the control architecture thought of the tandem PID is transferred to the attitude control, i.e., the outer loop is the position loop and the inner loop is the attitude loop. Thus the output u of the base controller b The calculation formula of (2) is as follows:
Figure BDA0002832663990000121
wherein K is ω For angular rate loop proportional gain, K I For angular rate loop integral gain, K D For the angular rate loop differential gain,
Figure BDA0002832663990000122
the integral operation is represented, and s represents the differential operation. Δω is the target angular rate ω t Difference from the current angular rate ω +.>
Figure BDA0002832663990000123
S is the integral operation symbol and s is the differential operation symbol.
Step 3, L1 augmentation adaptive control rate design
The cascade PID+L1 amplification self-adaptive control structure is provided by combining the traditional cascade PID control structure and the L1 self-adaptive structure, so that the self-adaptability and the robustness of the PID system are enhanced, different working conditions can be adapted to, and the structure can be adapted to the working conditions wider than a single cascade PID without resetting PID parameters like a single cascade PID in order to adapt to different working conditions.
This control method is described below in terms of a gesture control loop:
the angular velocity PID control shown in fig. 1 is outputted as a basic controller control signal in fig. 3, so that fig. 3 can be changed to that shown in fig. 4. The combined controller is now called the base controller by the cascade PID control. For the attitude inner loop controller, we choose the triaxial organism angular rate as the state variable x, and record the output of the cascade PID angular rate PID controller as u b
u b Output of PID controller at angular rate (19)
I.e.
Figure BDA0002832663990000131
Then the L1 augmented adaptive state prediction equation is
Figure BDA0002832663990000132
The parameter interpretation in the formula (21) can be described in the section L1 augmentation self-adaption, and when the system matrix A is used in actual engineering, the system matrix A needs to be determined by a system identification methodAnd B, matrix A SP The 3×3 matrix with the determinant being negative can be adjusted in the actual use process.
Figure BDA0002832663990000133
x (t) represents a true value, and actual measured gyroscope data of the integrated navigation system may be used instead.
L1 adaptive self-adaption rate calculation:
Figure BDA0002832663990000134
t represents the running period of the controller, and in the practical use process, the higher the running frequency of the self-adaptive control rate is, the more the total uncertainty of the system can be estimated in real time
Figure BDA0002832663990000135
After the uncertainty estimate of the system is obtained, the L1-adaptive control quantity output u can be calculated by the following equation ad The output of the self-adaptive controller is:
Figure BDA0002832663990000136
the authors have discrete forms of C(s) in actual use:
Figure BDA0002832663990000137
the above bilinear transformation is:
Figure BDA0002832663990000138
and (3) making:
Figure BDA0002832663990000139
f c for the filter cut-off frequency, f s Data sampling frequency, f when the author actually uses s =1kHz,f c =5 Hz, different systems may differ slightly.
Redefining and sorting into:
Figure BDA0002832663990000141
Figure BDA0002832663990000142
then the real-time usage pair
Figure BDA0002832663990000143
Performing second order low pass filtering
Figure BDA0002832663990000144
Wherein the method comprises the steps of
Figure BDA0002832663990000145
The intermediate variable is w [ n ], w [ n-1] is the value of w [ n ] of the last period, and w [ n-2] is the value of w [ n ] of the last period.
Step 4, calculating the total control input of the rotorcraft
The controller output calculated in step 2 and step 3 is used as the control input u (t) of the actuator of the rotorcraft, and the calculation formula is as follows
u(t)=u b (t)+u a (t) (32)
Wherein u is b (t) outputting the calculation engineering for the cascade PID angular rate controller as shown in the formula (18) of the step 2, u a (t) outputting it for L1 augmentation adaptive controllerThe calculation process is shown in (24).
When the self-adaptive control system is actually used, the operation period of the angular rate ring is 500HZ, the self-adaptive control self-adaptive rate operation rate of the L1 is 1000Hz, and the higher the operation frequency of the self-adaptive controller is, the faster the estimation of the uncertainty of the system is, and the better the control effect is.
The same applies to the height channel. The authors have greatly improved the wind resistance, attitude stability and overall stability of the rotorcraft after loading by using the control method in horizontal channel inner rings and vertical altitude inner rings.

Claims (1)

1. A rotor craft control method based on L1 augmentation self-adaption is characterized by comprising the following steps:
step 1, establishing a rotorcraft kinematics and dynamics model:
rotorcraft kinematic model:
Figure FDA0003943154200000011
wherein: m represents the weight of the unmanned aerial vehicle, g represents the gravitational acceleration, U 1 Represent lift, K d Represents the coefficient of resistance and,
Figure FDA0003943154200000012
the speeds in the three directions are indicated,
Figure FDA0003943154200000013
indicating acceleration in three directions->
Figure FDA0003943154200000014
The machine system rotates the matrix to the navigation system;
Figure FDA0003943154200000015
wherein: ψ denotes a heading angle, θ denotes a pitch angle, γ denotes a roll angle,
Figure FDA0003943154200000016
a rotation matrix representing the machine system to the navigation system;
kinetic model of rotorcraft:
Figure FDA0003943154200000017
Figure FDA0003943154200000021
represents the differentiation of the attitude angle, cos represents the cosine, sin represents the sine,
Figure FDA0003943154200000022
Figure FDA0003943154200000023
Figure FDA0003943154200000024
Figure FDA0003943154200000025
Figure FDA0003943154200000026
triaxial angular acceleration, +.>
Figure FDA0003943154200000027
Representing the triaxial angular velocity, I xx ,I yy ,I zz Representing the triaxial inertia, L representing the distance from the rotor centre to the aircraft centre of mass; u (U) 1 Representing lift force, U 2 Indicating torque in rolling direction, U 3 Indicating torque in pitch direction, U 4 Torque representing heading direction;
step 2, designing cascade PID control rate:
horizontal channel position loop control rate design:
v t =K p ·(p t -p)
wherein: p is p t For the desired position, p is the current position, v t To a desired speed, K p The ratio gain is the position loop;
horizontal channel speed loop control rate design:
a t =K v ·(v t -v)
wherein: v t For the desired speed, v is the current speed, a t For the desired acceleration, ψ is the heading angle, θ is the pitch angle, γ is the roll angle, K v Gain for the speed loop ratio;
Figure FDA0003943154200000031
θ t =arctan(-a x /g)
γ t =arctan(a y ·cosθ/g)
wherein: θ t For the desired pitch angle, θ is the current pitch angle, a x Represents tangential acceleration along the x-axis, a y Represents normal acceleration along the y-axis, gamma t Is the desired roll angle;
angle ring control rate design:
ω t =K θ ·(θ t -θ)
ω t for the desired angular rate, ω is the current angular rate, K θ For the proportional gain of the attitude loop, theta t The desired pitch angle is the θ, the current pitch angle;
angular rate ring control rate design:
Δω=ω t
Figure FDA0003943154200000032
wherein: k (K) ω For angular rate loop proportional gain, K I For angular rate loop integral gain, K D For the angular rate loop differential gain,
Figure FDA0003943154200000033
representing an integral operation, s representing a differential operation; Δω is the target angular rate ω t Difference from the current angular rate ω +.>
Figure FDA0003943154200000034
S is a differential operation symbol;
step 3, L1 augmentation adaptive control rate:
the L1 augmented adaptive state prediction equation is:
Figure FDA0003943154200000035
wherein: a state matrix A, a control matrix B, a matrix A SP A 3 x 3 matrix with a determinant of negative values, T representing the period of operation of the controller;
Figure FDA0003943154200000041
wherein: x% (t) represents the state error,
Figure FDA0003943154200000042
representing an estimated state, x (t) representing a true state;
adaptive rate of L1 adaptation:
Figure FDA0003943154200000043
output of the adaptive controller:
Figure FDA0003943154200000044
wherein: the expression of C(s) is
Figure FDA0003943154200000045
Figure FDA0003943154200000046
Representing a mapping of s-domain to z-domain, where z is defined as z=e τ is the sampling period and C1 is a given constant.
Step 4, general control input of the rotorcraft:
control output u of step 2 b (t) output u of step 3 a (t) adding to obtain the total control output u (t)
u(t)=u b (t)+u a (t)。
CN202011466461.XA 2020-12-12 2020-12-12 L1 augmentation self-adaption-based control method for rotary wing aircraft Active CN112666960B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011466461.XA CN112666960B (en) 2020-12-12 2020-12-12 L1 augmentation self-adaption-based control method for rotary wing aircraft

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011466461.XA CN112666960B (en) 2020-12-12 2020-12-12 L1 augmentation self-adaption-based control method for rotary wing aircraft

Publications (2)

Publication Number Publication Date
CN112666960A CN112666960A (en) 2021-04-16
CN112666960B true CN112666960B (en) 2023-06-23

Family

ID=75405691

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011466461.XA Active CN112666960B (en) 2020-12-12 2020-12-12 L1 augmentation self-adaption-based control method for rotary wing aircraft

Country Status (1)

Country Link
CN (1) CN112666960B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114355963A (en) * 2021-12-27 2022-04-15 海兴东方新能源发电有限公司 Method, system, terminal and medium for realizing height control based on rotor unmanned aerial vehicle
CN114265435B (en) * 2021-12-27 2024-07-05 海兴东方新能源发电有限公司 Method, system and device for realizing multi-airport accurate landing of rotor unmanned aerial vehicle
CN115236972B (en) * 2022-07-25 2024-07-23 中国安全生产科学研究院 Multi-robot formation control method based on double-closed-loop self-adaptive PID

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106292297A (en) * 2016-10-26 2017-01-04 成都市优艾维机器人科技有限公司 Based on PID controller and the attitude control method of L1 adaptive controller
CN111487868A (en) * 2020-04-23 2020-08-04 中国空气动力研究与发展中心设备设计及测试技术研究所 Adaptive control system and method suitable for integral feedback amplification system L1

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017078823A2 (en) * 2015-08-13 2017-05-11 The Board Of Regents For Oklahoma State University Modular autopilot design and development featuring bayesian non-parametric adaptive control
CN106444826A (en) * 2016-09-07 2017-02-22 广西师范大学 Flight control method of QUAV (Quadrotor Unmanned Aerial Vehicle)
CN106933104B (en) * 2017-04-21 2020-05-19 苏州工业职业技术学院 Hybrid control method for attitude and position of four-rotor aircraft based on DIC-PID
CN110456636B (en) * 2019-07-11 2022-04-01 西北工业大学 Self-adaptive control method of aircraft discrete sliding mode based on uncertainty upper bound estimation

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106292297A (en) * 2016-10-26 2017-01-04 成都市优艾维机器人科技有限公司 Based on PID controller and the attitude control method of L1 adaptive controller
CN111487868A (en) * 2020-04-23 2020-08-04 中国空气动力研究与发展中心设备设计及测试技术研究所 Adaptive control system and method suitable for integral feedback amplification system L1

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
cheng haoyu等.Asynchronously finite-time H∞ control for morphing aircraft.《Transactions of the Institute of Measurement and Control》.2018,全文. *
吕蓉蓉 ; 徐亮 ; 薛辰 ; 陆宇平 ; .阵风扰动下大柔性飞行器姿态跟踪控制设计.电子设计工程.2018,(18),全文. *
王易南 ; 陈康 ; 符文星 ; 闫杰 ; .带有攻角约束的高超声速飞行器航迹倾角跟踪控制方法.固体火箭技术.2016,(01),全文. *

Also Published As

Publication number Publication date
CN112666960A (en) 2021-04-16

Similar Documents

Publication Publication Date Title
CN112666960B (en) L1 augmentation self-adaption-based control method for rotary wing aircraft
CN111766899B (en) Interference observer-based quad-rotor unmanned aerial vehicle cluster anti-interference formation control method
CN107688295B (en) Four-rotor aircraft finite time self-adaptive control method based on rapid terminal sliding mode
CN112346470A (en) Four-rotor attitude control method based on improved active disturbance rejection control
CN106933104B (en) Hybrid control method for attitude and position of four-rotor aircraft based on DIC-PID
CN106249745B (en) The control method of four axis unmanned planes
CN110119089B (en) Immersion constant flow pattern self-adaptive quad-rotor control method based on integral sliding mode
CN109116860B (en) Nonlinear robust control method for three-rotor unmanned aerial vehicle
CN110058600B (en) Fault-tolerant control method for pose system under three-rotor unmanned aerial vehicle steering engine blocking fault
CN102830622A (en) Auto-disturbance-rejection automatic flight control method for four-rotor aircraft
CN107491083B (en) Four-rotor-wing autonomous landing method based on saturation self-adaptive sliding mode control
CN111650951B (en) Full-loop composite dynamic inverse tracking control method for complex track of quad-rotor unmanned aerial vehicle
CN109885074B (en) Finite time convergence attitude control method for quad-rotor unmanned aerial vehicle
CN112558621A (en) Decoupling control-based flying mechanical arm system
CN114281092A (en) Hypersonic aircraft coordination attitude control method based on sliding mode disturbance observer
CN114489094A (en) Unmanned gyroplane wind-disturbance-resistant control method based on acceleration feedback enhancement
CN110275542B (en) Four-rotor aircraft control method based on self-adaptive finite time control
CN110456781B (en) Space stability analysis method of aircraft control system
CN113031655B (en) Flying and lifting system accompanied with load lifting and control method thereof
CN111857171B (en) Unmanned aerial vehicle control method based on quality self-adaptive multilayer neural network
Salazar-Cruz et al. Real-time control of a small-scale helicopter having three rotors
CN112947522A (en) Hard air refueling attitude control method based on finite time observer
Song et al. Research on attitude control of quadrotor uav based on active disturbance rejection control
CN116679548A (en) Three-degree-of-freedom helicopter robust output feedback control method based on time-varying observer
Alcan et al. Robust hovering control of a quadrotor using acceleration feedback

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