• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Nonlinear adaptive flight control system:Performance enhancement and validation

    2023-05-19 03:40:36YueFENGZonghuSUNLioniWUYongshunWANGBinXIWengKhuenHOYnchengYOU
    CHINESE JOURNAL OF AERONAUTICS 2023年4期

    Yue FENG, Zonghu SUN, Lioni WU,*, Yongshun WANG, Bin XI,Weng Khuen HO, Yncheng YOU

    aSchool of Aerospace Engineering, Xiamen University, Xiamen 361102, China

    bDepartment of Electrical and Computer Engineering, National University of Singapore, Singapore 117576, Singapore

    KEYWORDSAdaptive control systems;Control engineering;Flight control systems;Nonlinear control systems;Robustness stability

    AbstractThe problem of decreasing stability margins in L1adaptive control systems is discussed and an out-of-loop L1adaptive control scheme based on Lyapunov’s stability theorem is proposed.This scheme enhances the effectiveness of the adaptation, which ensures that the system has sufficient stability margins to achieve the desired performance under parametric uncertainty,additional delays,and actuator faults.The stability of the developed control system is demonstrated through a series of simulations.Compared with an existing control scheme,the constant adjustment of the stability margins by the proposed adaptive scheme allows their range to be extended by a factor of 4–5,bringing the stability margin close to that of variable gain PD control with adaptively scheduled gains.The engineered practicability of adaptive technology is verified.A series of flight tests verify the practicability of the designed adaptive technology.The results of these tests demonstrate the enhanced performance of the proposed control scheme with nonlinear parameter estimations under insufficient stability margins and validate its robustness in the event of actuator failures.

    1.Introduction

    Adaptive control1has been a well-researched topic in control theory over several decades.It is a promising method in the field of automatic flight control, because it can improve the performance of control systems in the presence of uncertain modeling parameters, various external disturbances, and strong coupling and nonlinearity over a wide flight envelope.Such issues are often encountered in the field of Unmanned Aerial Vehicles (UAVs).2With the rapid development of UAVs in both military and civilian fields, several advances in adaptive flight control theory have helped to improve performance and robustness so that engineering implementations are now possible.In addition to several well-known modifications,3–5many new adaptive control methods have been developed in the past decade, such as neural network adaptive control,6L1adaptive control,7optimal control modification,8adaptive loop transfer recovery,9retrospective cost optimization,10and combined composite model reference adaptive control.11In spite of extensive research on adaptive control theory, the fact remains that few schemes have been subjected to flight validation experiments12,13and many technical problems related to engineering applications remain unresolved.

    L1adaptive control,which is an improved Model Reference Adaptive Control (MRAC) scheme,14is one method that has been employed in successful flight experiments.The L1filter is designed with trade-off between performance and robustness, and the adaptive law follows Lyapunov’s stability theorem15or the piecewise constant method.16Many simulations using L1adaptive control have been performed,17–19but very few have been considered in realistic scenarios, such as under computational latency,transport delays,and actuator dynamics.In terms of flight validation, a flight experiment using an L1adaptive control system was conducted on the NASA Air-STAR flight test vehicle in 2009,20,21and subsequent experiments have also been conducted.Flight tests based on the piecewise constant method have been implemented on Calspan’s Learjet22and a commercial autopilot,23and those based on Lyapunov’s stability theorem have been implemented on the Rascal UAV,24a small indoor UAV25(where the nonlinear adaptive parameters were omitted), and the SmallWhite UAV,26which was the initial application of nonlinear L1adaptive control.The oscillations caused by an insufficient timedelay margin in the AirSTAR flight test were first reported in Ref.27and have been further discussed in detail, and it is thought that a problem with the L1control system reduces the stability margins.

    An adaptive control system with strong nonlinearity is proposed in this paper.The construction of the L1control system has a design weakness in that the bound of the nonlinearity is limited.Systems based on the piecewise constant method can be approximated as Linear Time-Invariant (LTI) systems in most flight conditions, whereas those based on Lyapunov’s stability theorem have been proved to approximate an implementable LTI controller28despite the nonlinearity introduced through the nonlinear parameters.From a theoretical perspective, nonlinear adaptive control systems can adaptively adjust the control gain to improve the dynamic performance,thereby enhancing the control performance and robustness.However,current L1adaptive control systems can be approximately regarded as linear control systems and cannot ensure consistent performance when the system degenerates.This is the reason for the oscillations reported in Ref.27In the case of UAV systems, the margin requirements may be insufficient because of inaccurate models and changes in the installation center of gravity and load, which could result in system oscillations or even divergence.Therefore, there is an urgent need for a flight control system with strong and effective nonlinearity to satisfy the flight control requirements and ensure flight safety,which is the motivation for this study.

    Although the nonlinear parameters were considered in this previous study, the L1adaptive control system still performed as a linear control system,and so the stability margin could be no longer improved when it became insufficient.Nonlinear L1adaptive control systems make small-scale adjustments to the stability margin on the basis of the baseline controller with relatively weak adaptation, and the stability margin inevitably decreases compared with that of the baseline controller.The main factor that limits the ability of nonlinear L1adaptive control systems is the L1filter.The characteristics of this filter imply that its output has a time delay with respect to its input,which causes a phase difference.The feedback signals in nonlinear adaptation are filtered by the L1filter and sent back to the closed-loop system.The phase lag of the feedback signals significantly impacts the phase and time-delay margins.In addition, the phase differences between the feedback signals of the adaptive control law and those of the baseline control law prevent the nonlinear parameters from directly affecting the baseline controller.Based on these considerations, a new L1adaptive control scheme is proposed,in which the phase difference is eliminated and the reduced stability margin is overcome by redesigning the position of the filter in the feedback loop.The problem of decreased stability margins in L1control systems is solved,and the proposed approach enables the control performance to remain consistent with the expected performance in situations where the system has insufficient stability margins or reduced control performance.

    This paper is a follow-up study on the application of L1adaptive control theory in the field of flight control system.The main contribution is the proposal of a new L1adaptive control system with enhanced performance.First, compared with the original L1control system,the nonlinearity and adaptation are enhanced and the range in which the system stability margin can be adjusted is greatly expanded.Additionally,based on the fact that traditional Proportional–Derivative(PD)controllers are widely used in engineering,the adaptation for PD feedback gains implemented in this study has significant practicability.From the perspective of engineering validation, the proposed adaptive control system has been used in flight test verifications with nonlinear parameters, which increases the confidence in adaptive control as a possible flight control method.

    The remainder of this paper is organized as follows.Taking the SD-25 as the research object,Section 2 presents the design of the nonlinear out-of-loop L1adaptive control system based on a mathematical model for the pitch dynamics.Section 3 demonstrates the stability of the proposed control system.Simulation results and comparative stability analysis are given in Section 4,and the flight test validation is described in Section 5.Finally, Section 6 summarizes the main contributions of this study and states the conclusions.

    2.Design of out-of-loop L1adaptive control system

    2.1.Problem formulation

    UAV motion can be abstracted as six-degree-of-freedom rigid body motion, which can be described in the form of 12-state differential equations.29As the dynamics of UAVs can be decoupled,a simplified mathematical model for the UAV pitch dynamics is used for the design of control system, where the variables of interest are the pitch angle θ and the pitch rate q.Considering straight and level flight,the simplified equations of motion are as follows:

    where ρ is the air density,Vais the speed of the UAV through the surrounding air mass, c is the mean chord of the UAV wing, S is the planform area of the UAV wing, Jyis the moment of inertia of the y axis, δeis the elevator deflection,and the coefficients C0m, Cαm, Cqm, and Cδem are stability derivatives and are functions of state variables.

    The system described in Eq.(1)is time-variant and difficult to obtain an analytical solution for sophisticated control system design.The small-perturbation equation30around the singular points can be obtained based on certain assumptions,26where the time-varying coefficients are approximated as constants, and finally a second-order linear time-invariant is obtained.

    Considering x=[θ q]Tas the system state vector, it follows that the second-order state-space equation with zero initial conditions can be expressed as

    where θ is the pitch angle,q is the pitch rate,u(t)is the control signal,y(t)is the regulated output,A is the known state matrix,B is the known control input matrix, with (A, B) controllable,Λ is the known output matrix, Δ(x) is the uncertain part, and a1,a2,b1are the constants coefficient obtained by the linearization of Eq.(1),which are determined by the variables and stability derivatives at singular points.

    This paper considers systems in the presence of unknown constant parameters.The uncertain part Δ(x), which is a type of structured uncertainty matched by the control input, is defined as a linear combination of N known basis functions related to state variables with N unknown constants and can be described as

    where τ is a constant vector,and Φ(x)is the known regression vector.

    Considering the uncertainty caused by the state variables θ and q and the constant deviation of the control surface δΔ,we choose Φ(x)=[θ(t)q(t)δΔ]Tand τ=[τ1τ2τ3]T.

    2.2.Out-of-loop L1adaptive control architecture

    Fig.1 shows the structure of the out-of-loop L1adaptive control system.The L1adaptive control law15introduces a lowpass L1filter to guarantee robustness in the presence of fast adaptation while reducing the stability margin of the system.The concept of the Out-of-Loop L1Adaptive Control(OLLAC) system is to move the L1filter out of the feedback loop so as to overcome the problem of a decreased system stability margin and guarantee greater robustness.

    Remark 1.The filter in the OLLAC system,hereinafter referred to as the Out-of-Loop(OL)filter,is only effective for parameter estimation and can be separated from the closed-loop control module when analyzing the system stability.Therefore, the robustness of the OLLAC system is ensured at a high identification rate and the L1filter has no impact on the stability margins, especially for a low cut-off frequency.

    2.2.1.Baseline controller

    A feedforward and feedback controller,which is widely used in flight control, is selected as the baseline controller.Neglecting the uncertain component, the control law defined by Eq.(4)yields the closed-loop system with the zero initial conditions in Eq.(5):

    where Am?A-BKm, Bm?kcB, kcis the feedforward gain,Km=[kpkd] is the feedback gain, and r(t) is the reference input signal.According to the desired performances for Amand Bm,the baseline controller gains kp,kd,and kcare readily obtained.

    2.2.2.State predictor

    The closed-loop system defined by Eq.(2) with the feedback control law um(t)=-Kmx(t)under zero initial conditions is

    Fig.1 Structure of OLLAC system.

    where ua(t) is the adaptive term of the control signal.

    We consider the following state predictor under zero initial conditions:

    where τ^(t)=[^τ1(t)^τ2(t)^τ3(t)]Tis the estimate of the ideal unknown parameterτ.

    2.2.3.Adaptation law

    Based on the projection method described in the appendix and introducing the autocorrelation modification26to minimize the impact of noise, the adaptive process is governed by the following adaptation laws:

    2.2.4.Control law

    The control law for the OLLAC system is composed of the baseline feedback component um(t) and the adaptive component ua(t), which have the form of

    where the Laplace transform of η(t)is defined as η(s)?C(s)^τT.Here, C(s) is the OL filter, which is a Bounded-Input/Bounded-Output (BIBO)-stable and strictly proper transfer function with a Direct Current (DC) gain of C(0)=1.

    Remark 2.Nonlinearity is mainly reflected in η(t)Φ(x),which is designed to counteract the uncertain component Δ(x).^τ is the vector of adaptive gains, the first two of which are nonlinear terms and compensate for the feedback gain of the baseline controller so as to adjust the dynamic performance.The third gain is a linear term that compensates for the constant deviation of the control surfaces to realize accurate tracking.

    The original adaptive control law in Ref.15 is ua(s)=-C(s)((η′(s)-kcr(s)),where η′(s)is the Laplace transform of η′(t)?^τT(t)Φ(x).In contrast, the OL filter does not affect the feedback signals in the regression vector Φ (x).

    Finally,the out-of-loop L1adaptive controller is defined by the relationships in Eqs.(7)–(10).

    3.Analysis of OLLAC system

    3.1.Stability of closed-loop reference system

    A closed-loop reference system can be constructed for the OLLAC system described in Section 2.2 as

    Lemma 1.A necessary and sufficient condition for the system in Eq.(11) to be Bounded-Input/Bounded-State (BIBS)-stable with respect to r(t) is that H(s)?(sI-Am)-1B is a proper BIBO-stable transfer function.

    Proof.The Laplace transforms of the state equation and control law in Eq.(11) can be written as

    where xin(s)?(sI-Am)-1x0, x0is the initial value of x(t), and the subscript ‘‘r”represents the variables in the reference system.

    As Amis a Hurwitz matrix and both xin(t)and r(t)are uniformly bounded, xr(t) is also uniformly bounded from Eq.(13).Thus, the BIBS stability of the closed-loop system has been proved.Note that this design condition is not required to prove the L1stability31.

    Remark 3.The influence of the L1filter on the closed-loop stability is considered in the proof of Lemma 1 in Ref.7.However, the object on which C(s) acts in Eq.(11) is the uncertain parameter estimation, which is independent of the dynamics of the closed-loop system.According to the time-scale separation principle,32C(0)=1 in the above proof,and thus the stability of the closed-loop system is not affected by C(s).

    3.2.Steady-state performance of error dynamics

    The following error dynamics can be computed by subtracting Eq.(7) from Eq.(6):

    Lemma 2.The tracking error xein Eq.(14) is uniformly bounded and asymptotically stable with the adaptation laws described in Eq.(8).

    Proof.Consider the global radially unbounded quadratic Lyapunov function candidate of the form

    Consequently, Lemma 2 has been proved.

    4.Simulations and analysis

    4.1.Implementation and time domain simulations

    4.1.1.Configuration for simulation

    A closed-loop OLLAC system was constructed in MATLAB to verify the proposed system performance with the SD-25 Vertical Take-Off and Landing (VTOL) UAV as the design object.The system and controller are discretized with a sampling time of Ts= 0.01 s.Controller parameters are designed based on the linear model and the OLLAC system will be eventually verified using a nonlinear dynamic model.

    We considered the steady-state flight condition with a pitch angle of 3.05°,attack angle of 3.05°,pitch angle rate of 0 rad/s,speed of 30 m/s, and altitude of 100 m as the singular point.The linear small-perturbation equation described in Eq.(2)can then be written as

    According to Section 2.2,the baseline controller is designed as BKm=A-Amand kcB=Bm, and so we have negative feedback gains Km=[0.7 0.15] and feedforward gain kc=1.03.

    Remark 4.The exact values of the feedback gains are Km=[0.68 0.17].We replace these with the approximate values following the conventions of engineering applications.The cut-off frequency of the OL filter, the design of which is related to the trade-off between control performance and robustness,affects the speed at which the adaptation acts on the system and is limited by the frequency of system oscillations.Theoretically,the adaption rate should be as large as possible to ensure better learning efficiency.In practice,there is an upper limit that can be determined through simulations.The identification process will be unable to converge when the adaption rate is above this upper limit.The range of nonlinear parameters is normalized to make it symmetrical about zero, and the projection method is realized on the normalized range.

    4.1.2.Comparative experiment with L1adaptive control system

    The time domain responses of the OLLAC system and the original L1adaptive control structure are compared under the same design parameters.The simulation results are compared separately in systems with sufficient and insufficient stability margins.The step responses of nonlinear dynamic model are shown in Fig.2 and Fig.3.The simulation duration is 30 s,and the power of white noise for θ and q is chosen with 0.0001 and 0.2, respectively.To simulate an insufficient stability state in the aircraft,a pure time-delay module has been added to the actuator with dδ=15Ts.

    According to Fig.2,when the system margins are adequate,both adaptive control systems exhibit very similar performance.However, in terms of the 0.15 s delay of the actuator,the system with the L1control method diverges, whereas the pitch angle of the OLLAC system tracks the instructions and the adaptive parameter estimate converges, as shown in Fig.3.A comparison of the curves in Fig.3 shows that small oscillations occur in the first 5 s,and the trend of the nonlinear parameters is nearly the same.Over the next 5 s, the system with the L1controller diverges, while that controlled by OLLAC remains stable.The approximate time-delay margins of the systems under OLLAC and L1adaptive control can be obtained from the simulations in time domain with different values of dδ.The time-delay margin of the OLLAC system is about twice the value for L1adaptive control, which are 0.21 s and 0.08 s, respectively.

    Fig.2 Comparison of step responses between OLLAC and L1adaptive control with adequate stability margins.

    Fig.3 Comparison of step responses between OLLAC and L1adaptive control with insufficient stability margins.

    In conclusion,when a system is divergent due to an insufficient stability margin, the adaptation of the L1controller cannot make it stable.In contrast,steady tracking and convergent identification can be achieved by the OLLAC system.As a result, the adaptive performance of OLLAC is significantly better than that of the conventional system, and the stability margins and control performance are further optimized.

    Remark 5.Fig.3 provides a more intuitive understanding of the adaptation process.The first 10 s is the period in which the nonlinear parameters are identified; hence, the dynamic performance of the closed-loop system is reconstructed from an unstable state to a stable state.Correspondingly,the pitch angle gradually converges from the oscillating state during this period.The identification gradually converges in the next 10 s, and the system reaches a steady state in which the pitch oscillations disappear and the stability margin is improved.

    4.1.3.Comparative experiment with traditional control system

    In this section,we describe experiments conducted to illustrate the advantages of the OLLAC system versus traditional PD controllers.The baseline controller of the OLLAC system,which is a traditional PD controller, is noted as Controller 1.Controller 2 is also a traditional PD controller with different control gains to Controller 1, and the gains of Controller 2 are obtained by superimposing the converged estimates of nonlinear parameters in Fig.3 onto the gains of Controller 1.The simulations are carried out under the same conditions as described in Section 4.1.2, and a pure time-delay module is added to the actuator with dδ=15Ts.The open-loop Bode diagram and the pole-zero map used to compare Controllers 1 and 2 are presented in Fig.4.The step responses of the systems controlled by Controllers 1 and 2 and the OLLAC system are shown in Fig.5.

    As we can see from Fig.4(a), the 0.15 s delay of the actuator exceeds the margin of Controller 1,but remains within that of Controller 2.This is also illustrated in Fig.5 where the step response curve of Controller 1 diverges while that of Controller 2 indicates smooth tracking.Of more interest is that the pitch angle of the OLLAC system in Fig.5 is consistent with that of Controller i in the first few seconds,and converges in essentially the same way as the curves of Controller 2.The estimations of the OLLAC system enable adjustments toward a more stable state.

    In abstract terms, the closed-loop systems of Controllers 1 and 2 can be regarded as two points in different states,having different but unchanged performance.The OLLAC system can be considered equivalent to its baseline controller (Controller 1)at the initial time,because the nonlinear parameters are initially set to 0.Likewise, at the end of the simulation, the OLLAC system can be regarded as essentially equivalent to Controller 2.Therefore, the OLLAC system can be regarded as operating along a line that has one more time axis than the points of Controllers 1 and 2, and it continuously moves from one state (i.e., from Controller 1 or 2) to the other (i.e.,to Controller 2 or 1).

    In summary, the OLLAC system can be regarded as an adaptive PD control system with the baseline controller as its foundation.The nonlinear parameters compensate for the feedback gains of the baseline controller to create greater stability margins and better control performance in the presence of system degradation or uncertainties.Compared with the traditional control method,the robustness of the OLLAC system is enhanced.

    4.2.Performance analysis and verification

    As previously stated, the OLLAC system can move from one point to another.This section discusses the limiting distances between the initial state and the state that can be achieved by the OLLAC system, L1adaptive control, and the PD controller with variable gains.

    The closed-loop systems with the OLLAC,L1adaptive control and gain-scheduling PD controller are noted as Systems A,B and C, respectively.The same parameters are selected for adaptation, and the control gains for the baseline controller of Systems A and B are the same as the initial PD gains in System C.

    4.2.1.Stability margin and robustness

    The stability margin analysis method34for linear systems is not generally suitable for nonlinear systems.To solve this problem,the nonlinear parameters^τ1and^τ2,which introduce nonlinearity into the system, were placed in a mesh.If the mesh is sufficiently dense, the nonlinear parameters in each cell can be processed as constants.As a result, the nonlinear part of the control system is converted to a constant term,and the nonlinear system can be analyzed based on the stability margin analysis of linear systems.

    Two-dimensional grid coordinates were constructed based on the ranges of ^τ1and ^τ2to form a grid with a resolution of 0.02.The stability of Systems A and B was then analyzed through Bode diagrams.Likewise, the coordinates of System C were constructed and indicated the increment based on the initial PD control gains.The comparison results for the effects of adaptive parameters on the stability margins are shown in Fig.6 and Table 1.

    Fig.4 Frequency domain analysis of Controllers 1 and 2.

    Fig.5 Comparison of time-domain response among Controller 1, Controller 2, and OLLAC.

    Fig.6 Effects of adaptive parameters on stability margins for Systems A, B, and C.

    Table 1 Stability margins analysis of Systems A, B, and C.

    Fig.6 shows how the stability margins of the closed-loop Systems A,B,and C change as the gains of the advanced controller vary in response to the nonlinear parameters or gain scheduling.The ranges on the z-axis of Systems A, B, and C indicate the extent to which they can ‘‘move”the system.The traditional controller will be represented as a plane parallel to the x-y plane and the values of z-axis are its stability margins.Apparently, traditional controllers are unable to reconfigure the system dynamic performance.As we can see in Fig.6,Systems A and C display nearly the same trend,while the curvature of the stability margin for System B is significantly smaller than that of the other two systems.This means that,for the same change rate of^τ1and^τ2,the closed-loop performance of System B changes at the smallest rate.From Table 1, when the adaptive parameters change over the same ranges, the phase, gain, and time-delay margins of System A have ranges that are 4, 4.2, and 5 times those of System B,respectively.By extension, the parameters for System B range over a small interval around the figures for the baseline controller under the adaptation.In particular, the upper bounds of the phase and time-delay margins of System B are less than the figure for the baseline controller,which means that the stability margin is reduced when the adaptation of System B is introduced.This explains why System B diverges in Fig.3 under the same parameter identification as System A.

    In summary,the adaptation in the original L1adaptive control system only works around a fixed state,or‘‘datum point,”and the system is only adjusted slightly with a small amplitude in the vicinity of this datum point,which leads to limited adaptation performance.In contrast,the OLLAC system moves the datum point to adjust the system quickly over a larger range and produces more efficient adaptation.

    4.2.2.Control performance

    The control performance of Systems A,B,and C was analyzed under the nominal model.The grid coordinates were consistent with those constructed in Section 4.2.1,and the step responses of Systems A,B and C were compared in the time domain.We define eAas the difference between the step responses of Systems A and C.Similarly, eBis the difference between the step responses of Systems B and C.The mean and variance of eAand eBin each cell are compared in Fig.7.

    As can be seen from Fig.7,the control performance of System A is consistent with that of System C,whereas that of System B is significantly affected by nonlinear parameters.As the nonlinear parameters deviate significantly from the initial values,the performance of System B shows obvious differences to that of System C.Therefore, with respect to robustness and control performance, System A is highly consistent with an automatic variable gain controller.

    4.3.Effects of matrix Q on control objectives

    We now explore the design of the matrix Q in the algebraic Lyapunov equation.According to the analysis in Section 3.2,any positive-definite symmetric matrix Q can be solved to obtain a matrix P that conforms to the stability condition,but the matrix P calculated from different matrices Q will be different.The influence of matrix P on the control system is analyzed in the following.We define the adaptation error ε(t)?xTe(t)PB, and rewrite Eq.(8) as

    Fig.7 Comparison of mean and variance results.

    The value of λ is determined completely by the matrix P.As P is a symmetric positive-definite matrix,λ>0.Thus,Eq.(25)becomes

    where χ=b1/P12is a constant.

    According to Eq.(28), the value of λ determines how the different tracking tasks should be weighted.For λ > 1, the adaptation is driven primarily by eq(t), and the main object of the controller is the pitch rate.For λ < 1, the adaptation is driven primarily by eθ(t), and the main object of the controller is the pitch.Note that the smooth signal ^q(t)cannot be completely tracked by q(t)due to sensor noise,which makes it difficult to achieve eq(t)≠0.Considering the steady-state system, we have ˙^τ(t)=0, and thus eθ(t)=-λeq(t).As a result,there is always a pitch tracking error in the case of measurement noise.

    Taking the SD-25 VTOL UAV as the object, the matrix Q is designed as a diagonal matrix Q=diag(Q11;Q22).If we define ?Q11/Q22, the relationship between ζ and λ is shown in Fig.8.The simulation results in the time domain for λ = 0.5, 1.0, and 2.0 are shown in Fig.9 and Table 2.

    Fig.8 indicates that ζ is a monotonically decreasing function of λ and the curve is approximately exponential.According to Fig.9 and Table 2,higher values of λ result in a greater pitch tracking error.To achieve better tracking performance,λ < 1 can be ensured through the appropriate design of ζ and the matrix Q.

    Fig.8 Curve of ζ - λ.

    5.Experimental results

    5.1.Flight test overview

    An SD-25 VTOL UAV is used to evaluate the OLLAC system,as shown in Fig.10.Past flight tests on the SD-25 have suffered from the problem of an insufficient stability margin,which could potentially lead to oscillations in attitude angles.To solve this problem, the adaptive control system, which is used to track the attitude angles,is designed and implemented in the inner loop and commands are given by the outer loop.

    The SD-25 VTOL UAV is equipped with four elevon segments (inner/outer) and four rudder segments for the V-tail.Its wingspan is 3.49 m and its takeoff weight is 25 kg.The adaptive control algorithm was written in the C programming language and was run at a frequency of 100 Hz on a Freescale MPC5674 autopilot.The sensors on the autopilot include an attitude and heading reference system,a GPS receiver(NovAtel),and a digital atmospheric data system.The planned route was a circle with a radius of 400 m.

    Under the single-variable method, the flight test was completed in two sorties.Sortie I and Sortie II were designed to test the control performance and stability for the original L1control system in Ref.15 and the OLLAC system,respectively.The robustness of the OLLAC system was verified in Sortie II by reducing the elevator efficiency to 60%to imitate actuator failure (Sortie II - actuator failure).

    Fig.9 Time-domain responses for different λ.

    Table 2 Analysis of tracking error for different λ.

    Fig.10 Photograph of SD-25 VTOL UAV.

    5.2.Data analysis

    The experimental parameters were selected to be the same as the designed parameters in Section 4.1.1.To compare the proposed control system with the existing control methods,conditions such as the flight status, outer-loop control strategy,control parameters, and commands were the same for Sorties I and II.

    Table 3 Analysis of pitch tracking error.

    5.2.1.Data analysis for comparative experiment

    Fig.11 and Table 3 compare the flight results over 100 s between Sortie I and Sortie II.The altitude commands corresponding to the two groups of data are consistent with descending 50 m and climbing 50 m.

    For similar pitch commands, Sortie I exhibits smallamplitude high-frequency pitch oscillations, but these are not observed in Sortie II.The data in Table 3 indicate that the indices of the tracking error in Sortie II are smaller than those in Sortie I, and in particular, the mean pitch tracking error of Sortie I is 14.8 times that of Sortie II.

    The flight data show that the problem of the insufficient stability margin in the SD-25 flight still exists under the original L1control system.However, this problem is solved by the OLLAC system, which achieves better control performance and greater stability margin.

    Remark 6.It can be seen from Figs.11(c) and (f) that, for similar pitch commands, the adaptive control component is different.This problem was further analyzed and found to be caused by differences in the linear parameter estimation, which can be understood as a difference in the trim elevator deflection of the UAV due to changes in the center of gravity or weight between the two flights.

    Fig.11 Aircraft response for original L1adaptive control (Sortie I) and OLLAC (Sortie II) systems.

    Fig.12 Aircraft response under OLLAC system with actuator failure.

    5.2.2.Data analysis for robustness verification

    Fig.12 shows flight data for a period of 80 s from Sortie II in which actuator failure was simulated.The dotted rectangles denote the period (12.0 s ≤t ≤49.8 s) in which the elevator effect was reduced to 60 %.

    At the moments of elevator failure and recovery,a tracking error of approximately 1.0°can be observed in Figs.12(a)and(b).It follows that the parameter identification works and converges within the next 7 s (Figs.12(c)-12(f)), and the tracking error gradually decreases to zero within 10 s due to the enhanced adaptation.Therefore, the OLLAC system is robust and can identify and compensate for the uncertainties to achieve zero-error tracking when actuator failure occurs.

    6.Conclusions

    The design, analysis, simulation, and flight test validation of an OLLAC system have been described in this paper.The problem of the insufficient stability margin15,26has been discussed, and a new structure has been proposed to enhance the nonlinearity and effectiveness of the adaptation process.As the theoretical basis for application, the BIBS stability and uniform boundedness of the tracking error have been proved.A method for stability margin analysis of nonlinear control systems has been established, in which the nonlinear system can be regarded as a linear system within each mesh gridded by the nonlinear parameters.The stability margins of the original L1control system, when adjusted by the adaptive scheme, have obvious boundaries and are smaller than those of the baseline controller.The proposed OLLAC system significantly extends these boundaries and achieves unrestricted stability margins.Comparative experiments focusing on the control performance and robustness have been realized in the frequency and time domains, respectively.Compared with the original control system, the OLLAC system reconstructs the closed-loop system in a wider range and at a faster rate,which leads to a high degree of nonlinearity and enhanced adaptation to compensate for uncertainties.

    A practical flight test validation of the proposed control system was realized on a demonstrator with an insufficient stability margin, and the robustness was further verified by imitating an actuator failure.The pitch oscillations observed under a baseline controller and L1adaptive control system were significantly reduced by the proposed system,which exhibits sufficient robustness to ensure accurate tracking and stability.In addition, the influence of the matrix in the Lyapunov equation on the control performance was analyzed in terms of its physical meaning,and a design method was further discussed.The proposed OLLAC scheme adds new capabilities in terms of enhanced performance and robustness for adaptive flight control systems and promotes the inclusion of nonlinear terms in engineering applications, which further increases the viability of adaptive control for UAVs.

    The engineering practicability of the OLLAC system comes from the fact that the baseline PD controller is already widely used in engineering.If MIMO systems are involved,there may be implementation difficulties because of the number of parameters,complex control systems,and other issues.Future research will investigate how to implement and optimize the OLLAC system under MIMO systems, such as by designing an appropriate parameter setting method or handling the MIMO system.If possible,we will implement the OLLAC system for the three-axis attitude control of a vehicle to solve the problem of modal coupling under large maneuvers.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China(No.U21B6003)and the China Scholarship Council (CSC, No.202006310096).

    Appendix A.A projection method35is introduced in the adaptive law to ensure that the adaptive parameters have a boundary.Considering a parameter vector θ with boundary constraints of the form ‖θ‖≤θmax, we define the continuous projection operator as

    日日干狠狠操夜夜爽| a级毛色黄片| 久久久久久伊人网av| 在线观看一区二区三区| 91av网一区二区| 精品久久久久久久久亚洲| 婷婷六月久久综合丁香| 久久这里只有精品中国| 久久精品夜夜夜夜夜久久蜜豆| 国产精品福利在线免费观看| 久久99蜜桃精品久久| 五月伊人婷婷丁香| 日韩欧美一区二区三区在线观看| 又粗又硬又长又爽又黄的视频 | 日本免费一区二区三区高清不卡| 国内精品美女久久久久久| 国产爱豆传媒在线观看| 亚洲精品亚洲一区二区| 成人三级黄色视频| 色综合亚洲欧美另类图片| 波多野结衣高清作品| 亚洲欧洲日产国产| 99热只有精品国产| 长腿黑丝高跟| 99在线人妻在线中文字幕| 成年女人永久免费观看视频| 亚洲无线在线观看| 最近视频中文字幕2019在线8| 最近中文字幕高清免费大全6| 免费搜索国产男女视频| 美女内射精品一级片tv| 久久久国产成人免费| 午夜激情欧美在线| 亚洲欧美日韩高清专用| 国产黄色视频一区二区在线观看 | 亚洲欧美中文字幕日韩二区| 亚洲国产精品成人久久小说 | 亚洲激情五月婷婷啪啪| 只有这里有精品99| av国产免费在线观看| 国产精品国产三级国产av玫瑰| 国产成人福利小说| 久久久精品欧美日韩精品| 欧美一区二区亚洲| 精品久久久久久久末码| 亚洲四区av| av女优亚洲男人天堂| 久久久久性生活片| 女的被弄到高潮叫床怎么办| 校园春色视频在线观看| 天堂影院成人在线观看| 啦啦啦韩国在线观看视频| 国产欧美日韩精品一区二区| 亚洲国产日韩欧美精品在线观看| 熟妇人妻久久中文字幕3abv| 岛国毛片在线播放| 在线观看66精品国产| av.在线天堂| 蜜臀久久99精品久久宅男| 国产精品久久久久久亚洲av鲁大| 欧美bdsm另类| 久久人人爽人人爽人人片va| 激情 狠狠 欧美| 国产黄a三级三级三级人| 国产一区二区在线av高清观看| 国产精品免费一区二区三区在线| 91久久精品国产一区二区三区| 亚洲欧美日韩卡通动漫| 亚洲四区av| 日韩一区二区三区影片| 日本在线视频免费播放| 亚洲三级黄色毛片| 91精品一卡2卡3卡4卡| 国产精品.久久久| 亚洲精品亚洲一区二区| 久久亚洲精品不卡| 国产一区亚洲一区在线观看| 欧美性感艳星| 日本一二三区视频观看| 国产视频内射| 午夜精品国产一区二区电影 | 老司机福利观看| 亚洲,欧美,日韩| 国产在线精品亚洲第一网站| 成人二区视频| 男人狂女人下面高潮的视频| 日韩亚洲欧美综合| 美女 人体艺术 gogo| 99久久成人亚洲精品观看| 国产欧美日韩精品一区二区| 99久久九九国产精品国产免费| 日本免费一区二区三区高清不卡| 国内揄拍国产精品人妻在线| 国产成人精品婷婷| 欧美性猛交╳xxx乱大交人| 久久6这里有精品| 久久久久网色| 欧美成人免费av一区二区三区| 国产日韩欧美在线精品| 91狼人影院| 在线a可以看的网站| 波多野结衣高清无吗| 搡女人真爽免费视频火全软件| 男女边吃奶边做爰视频| 校园人妻丝袜中文字幕| videossex国产| 简卡轻食公司| 又爽又黄a免费视频| 欧美高清成人免费视频www| 日本-黄色视频高清免费观看| 日本与韩国留学比较| 日韩欧美一区二区三区在线观看| 国产精品日韩av在线免费观看| 一区二区三区免费毛片| 韩国av在线不卡| 精品人妻一区二区三区麻豆| 久久精品国产亚洲av涩爱 | av视频在线观看入口| 国产淫片久久久久久久久| 国产精品蜜桃在线观看 | 深夜a级毛片| 色哟哟·www| 午夜免费激情av| 亚洲av不卡在线观看| 免费人成在线观看视频色| 亚洲色图av天堂| 日韩欧美 国产精品| 男女做爰动态图高潮gif福利片| 九色成人免费人妻av| 丝袜美腿在线中文| 给我免费播放毛片高清在线观看| 国产美女午夜福利| 色尼玛亚洲综合影院| 91精品一卡2卡3卡4卡| 在线天堂最新版资源| 国产真实伦视频高清在线观看| 高清日韩中文字幕在线| 久久久精品欧美日韩精品| 亚洲最大成人手机在线| 国产91av在线免费观看| 伦精品一区二区三区| 12—13女人毛片做爰片一| 久久中文看片网| 日韩成人av中文字幕在线观看| 久久久久网色| 中文在线观看免费www的网站| 国产精品国产高清国产av| 亚洲乱码一区二区免费版| 亚洲熟妇中文字幕五十中出| 色播亚洲综合网| 久久精品久久久久久久性| 午夜免费男女啪啪视频观看| 欧美日韩一区二区视频在线观看视频在线 | 国产精品,欧美在线| 久久久久久久午夜电影| 日本一二三区视频观看| 国产久久久一区二区三区| 亚洲精华国产精华液的使用体验 | 好男人在线观看高清免费视频| 中文字幕精品亚洲无线码一区| 真实男女啪啪啪动态图| 国产精品一区二区三区四区免费观看| 在线观看av片永久免费下载| 最好的美女福利视频网| 免费不卡的大黄色大毛片视频在线观看 | 午夜精品国产一区二区电影 | 一本久久中文字幕| 久久久久性生活片| 九九爱精品视频在线观看| 在线国产一区二区在线| 国产亚洲精品久久久久久毛片| 五月伊人婷婷丁香| 国产亚洲91精品色在线| 美女内射精品一级片tv| 欧美3d第一页| 老师上课跳d突然被开到最大视频| 中文资源天堂在线| 欧美精品一区二区大全| 在线免费观看的www视频| 99国产精品一区二区蜜桃av| 国产视频内射| 91狼人影院| 亚洲乱码一区二区免费版| 久久精品综合一区二区三区| 18禁裸乳无遮挡免费网站照片| 欧洲精品卡2卡3卡4卡5卡区| 国产免费男女视频| 国内少妇人妻偷人精品xxx网站| 国产成人a区在线观看| 天堂av国产一区二区熟女人妻| 黄色配什么色好看| 插逼视频在线观看| 22中文网久久字幕| 久久这里只有精品中国| 亚洲国产精品成人综合色| 久久精品国产亚洲网站| 国产在线精品亚洲第一网站| 不卡一级毛片| 亚洲一级一片aⅴ在线观看| 六月丁香七月| 亚洲美女搞黄在线观看| 赤兔流量卡办理| 一区二区三区高清视频在线| 精品久久国产蜜桃| 高清毛片免费观看视频网站| 色尼玛亚洲综合影院| 熟女电影av网| 国产精品免费一区二区三区在线| 18禁裸乳无遮挡免费网站照片| 老师上课跳d突然被开到最大视频| .国产精品久久| 国产真实伦视频高清在线观看| 免费观看人在逋| 国产精品一二三区在线看| 国产探花在线观看一区二区| 国产精品久久久久久久电影| 在线免费观看的www视频| 亚洲欧美清纯卡通| 日本一本二区三区精品| 寂寞人妻少妇视频99o| 黑人高潮一二区| 一本一本综合久久| 亚洲在线观看片| 不卡一级毛片| 亚洲一级一片aⅴ在线观看| 在线观看午夜福利视频| 级片在线观看| 中文亚洲av片在线观看爽| 免费电影在线观看免费观看| 好男人视频免费观看在线| 免费黄网站久久成人精品| 人人妻人人澡人人爽人人夜夜 | 男人的好看免费观看在线视频| 国产精品久久久久久久久免| 亚洲人成网站在线观看播放| 国产女主播在线喷水免费视频网站 | 夜夜夜夜夜久久久久| 国产伦在线观看视频一区| 2022亚洲国产成人精品| 有码 亚洲区| 秋霞在线观看毛片| 久久久久久久久中文| 亚洲av中文av极速乱| 中文字幕熟女人妻在线| 亚洲国产色片| 亚洲欧美成人精品一区二区| 此物有八面人人有两片| 青春草亚洲视频在线观看| 国产中年淑女户外野战色| 午夜免费男女啪啪视频观看| 国内少妇人妻偷人精品xxx网站| 久久久精品欧美日韩精品| 国产中年淑女户外野战色| 成人永久免费在线观看视频| 五月伊人婷婷丁香| 99久久无色码亚洲精品果冻| 国产91av在线免费观看| 亚洲中文字幕日韩| 九九久久精品国产亚洲av麻豆| 久久韩国三级中文字幕| 亚洲欧美精品综合久久99| 在线观看66精品国产| 国产私拍福利视频在线观看| 高清毛片免费观看视频网站| 日本一本二区三区精品| av在线蜜桃| 天堂√8在线中文| 欧美色视频一区免费| 国产在线精品亚洲第一网站| 九草在线视频观看| www日本黄色视频网| 性插视频无遮挡在线免费观看| 99久久人妻综合| 精品不卡国产一区二区三区| av国产免费在线观看| 亚洲精品国产av成人精品| 国产v大片淫在线免费观看| 久久鲁丝午夜福利片| 黄片wwwwww| 成人午夜高清在线视频| 中文字幕制服av| 国产高清不卡午夜福利| 色哟哟·www| 美女 人体艺术 gogo| 欧美潮喷喷水| 亚洲中文字幕一区二区三区有码在线看| 久久婷婷人人爽人人干人人爱| 成人午夜高清在线视频| 又爽又黄无遮挡网站| 国产黄片美女视频| 高清日韩中文字幕在线| 毛片女人毛片| 中文在线观看免费www的网站| 波野结衣二区三区在线| 久久久久久久久大av| 国产成人a∨麻豆精品| 啦啦啦韩国在线观看视频| 蜜臀久久99精品久久宅男| 日日摸夜夜添夜夜爱| 国产av一区在线观看免费| 寂寞人妻少妇视频99o| 精品久久久久久久久av| 韩国av在线不卡| 久久久久免费精品人妻一区二区| 最好的美女福利视频网| 一级毛片电影观看 | 狂野欧美激情性xxxx在线观看| 国产精品蜜桃在线观看 | 久久精品影院6| 最近的中文字幕免费完整| av专区在线播放| 亚洲成人中文字幕在线播放| 99riav亚洲国产免费| 国产欧美日韩精品一区二区| 最好的美女福利视频网| 日本黄大片高清| 国产成人精品久久久久久| a级毛片免费高清观看在线播放| 国产v大片淫在线免费观看| 变态另类成人亚洲欧美熟女| 精品久久国产蜜桃| 国产69精品久久久久777片| 国产精品.久久久| 99久久中文字幕三级久久日本| 久久久久九九精品影院| 国产精品野战在线观看| 午夜福利在线在线| 欧美一区二区亚洲| 国产精品久久久久久久久免| 三级男女做爰猛烈吃奶摸视频| 深夜精品福利| 日韩 亚洲 欧美在线| 日日摸夜夜添夜夜添av毛片| avwww免费| 午夜爱爱视频在线播放| 内射极品少妇av片p| 菩萨蛮人人尽说江南好唐韦庄 | 黄片无遮挡物在线观看| 日日啪夜夜撸| 亚洲欧美成人精品一区二区| 九色成人免费人妻av| 国产一级毛片七仙女欲春2| 亚洲最大成人手机在线| av在线天堂中文字幕| 一个人观看的视频www高清免费观看| 99久久成人亚洲精品观看| 久久久久久九九精品二区国产| 亚洲av第一区精品v没综合| 久久久精品大字幕| 午夜福利在线在线| 六月丁香七月| 欧美最黄视频在线播放免费| 日韩欧美精品v在线| 国内精品宾馆在线| 国产精品不卡视频一区二区| 亚洲欧美成人综合另类久久久 | 能在线免费看毛片的网站| 午夜精品一区二区三区免费看| 欧美又色又爽又黄视频| 亚洲欧美日韩卡通动漫| 国产成人a区在线观看| 中文字幕av成人在线电影| 丝袜美腿在线中文| 最近最新中文字幕大全电影3| 久久欧美精品欧美久久欧美| 欧美zozozo另类| 久久精品国产亚洲av天美| 成年女人看的毛片在线观看| 亚洲内射少妇av| 99久国产av精品| 小说图片视频综合网站| 色综合站精品国产| 精品日产1卡2卡| 欧美精品一区二区大全| 亚洲aⅴ乱码一区二区在线播放| 日韩一区二区视频免费看| 亚洲综合色惰| 狂野欧美白嫩少妇大欣赏| 亚洲欧美日韩东京热| 中文精品一卡2卡3卡4更新| 美女高潮的动态| 中文亚洲av片在线观看爽| 亚洲欧洲日产国产| 特大巨黑吊av在线直播| 国产成人精品久久久久久| 免费看美女性在线毛片视频| 精品久久久久久久久久久久久| 男的添女的下面高潮视频| 欧美zozozo另类| 日韩视频在线欧美| 亚洲真实伦在线观看| 久久国内精品自在自线图片| 亚洲国产日韩欧美精品在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 女人十人毛片免费观看3o分钟| 天堂影院成人在线观看| 村上凉子中文字幕在线| 男插女下体视频免费在线播放| 精品久久久久久久久久久久久| 亚洲精品亚洲一区二区| 欧美成人一区二区免费高清观看| 亚洲人成网站高清观看| 国产一区二区三区av在线 | 观看免费一级毛片| 深夜a级毛片| 久久久久九九精品影院| 国产成人一区二区在线| 日韩三级伦理在线观看| 麻豆av噜噜一区二区三区| 久久精品91蜜桃| 听说在线观看完整版免费高清| 国产中年淑女户外野战色| 一个人看的www免费观看视频| 91久久精品国产一区二区三区| 美女高潮的动态| 男的添女的下面高潮视频| 久久午夜福利片| 国产午夜精品论理片| 麻豆国产97在线/欧美| 亚洲国产精品国产精品| av免费观看日本| 秋霞在线观看毛片| 九草在线视频观看| 久久精品国产清高在天天线| 又爽又黄无遮挡网站| 全区人妻精品视频| 岛国在线免费视频观看| 国产激情偷乱视频一区二区| 婷婷六月久久综合丁香| 啦啦啦啦在线视频资源| 99久久久亚洲精品蜜臀av| 男插女下体视频免费在线播放| 国产精品乱码一区二三区的特点| 97超视频在线观看视频| 成熟少妇高潮喷水视频| 2022亚洲国产成人精品| 给我免费播放毛片高清在线观看| 黄色日韩在线| 免费观看人在逋| 国内精品久久久久精免费| 免费av毛片视频| 99热这里只有是精品在线观看| 我的老师免费观看完整版| 中国美女看黄片| 亚洲最大成人中文| 欧美一区二区国产精品久久精品| 丝袜喷水一区| 五月伊人婷婷丁香| 免费电影在线观看免费观看| 亚洲电影在线观看av| 色播亚洲综合网| 精品少妇黑人巨大在线播放 | 嫩草影院精品99| 亚洲人成网站高清观看| 色综合亚洲欧美另类图片| 日韩一区二区三区影片| 欧美精品一区二区大全| 春色校园在线视频观看| 欧美不卡视频在线免费观看| 淫秽高清视频在线观看| 成人综合一区亚洲| 亚洲精品色激情综合| 欧美一级a爱片免费观看看| 爱豆传媒免费全集在线观看| 在线播放国产精品三级| 亚洲精品日韩av片在线观看| 亚洲国产色片| 日本黄色视频三级网站网址| 日韩欧美精品免费久久| 黄色日韩在线| 免费人成视频x8x8入口观看| 看免费成人av毛片| 日韩强制内射视频| 精品久久久久久久人妻蜜臀av| 波多野结衣高清作品| 99热这里只有是精品50| 久久中文看片网| 直男gayav资源| 亚洲精品成人久久久久久| 在线免费观看不下载黄p国产| 熟女电影av网| 69av精品久久久久久| 久久鲁丝午夜福利片| 国产精品女同一区二区软件| 国产精品野战在线观看| 九色成人免费人妻av| 99热这里只有是精品50| 国产精品一及| 国产乱人偷精品视频| 黄片wwwwww| 欧美成人精品欧美一级黄| 中文欧美无线码| 最近的中文字幕免费完整| ponron亚洲| 精品欧美国产一区二区三| 国产极品精品免费视频能看的| 日日干狠狠操夜夜爽| kizo精华| 国产精品一区二区三区四区久久| 中文在线观看免费www的网站| 国产精品一及| 国产又黄又爽又无遮挡在线| 99热网站在线观看| 中文字幕免费在线视频6| 久久精品国产自在天天线| 精品日产1卡2卡| www日本黄色视频网| 蜜桃亚洲精品一区二区三区| 中文欧美无线码| 男人狂女人下面高潮的视频| 91aial.com中文字幕在线观看| 亚洲一区高清亚洲精品| 国产精品福利在线免费观看| 久久九九热精品免费| 国产熟女欧美一区二区| 一区二区三区免费毛片| 久久这里有精品视频免费| 免费不卡的大黄色大毛片视频在线观看 | 久久精品国产鲁丝片午夜精品| 男女做爰动态图高潮gif福利片| 国产精品蜜桃在线观看 | 舔av片在线| 国产亚洲精品av在线| 亚洲最大成人av| 免费av观看视频| 我要看日韩黄色一级片| 中文资源天堂在线| 精品久久国产蜜桃| 欧美变态另类bdsm刘玥| 爱豆传媒免费全集在线观看| 国产成人a区在线观看| 国产精品国产三级国产av玫瑰| 国产成人精品一,二区 | 久久久久久久久久久丰满| 少妇被粗大猛烈的视频| 如何舔出高潮| 91av网一区二区| 青春草国产在线视频 | 欧美性感艳星| 国产成人一区二区在线| 天美传媒精品一区二区| 国产午夜精品论理片| 一区二区三区四区激情视频 | 亚洲国产精品成人综合色| 欧美不卡视频在线免费观看| 99久久无色码亚洲精品果冻| 老熟妇乱子伦视频在线观看| 亚洲精品影视一区二区三区av| 在线观看66精品国产| 亚洲精品久久国产高清桃花| 大香蕉久久网| 国产高清有码在线观看视频| 白带黄色成豆腐渣| av卡一久久| 亚洲国产高清在线一区二区三| 亚洲av免费高清在线观看| 一个人看的www免费观看视频| 卡戴珊不雅视频在线播放| 内射极品少妇av片p| 少妇猛男粗大的猛烈进出视频 | avwww免费| 国产在视频线在精品| 亚洲自偷自拍三级| 国产欧美日韩精品一区二区| 乱人视频在线观看| 日韩,欧美,国产一区二区三区 | 99久国产av精品国产电影| 日本黄大片高清| 久久精品夜色国产| 自拍偷自拍亚洲精品老妇| 久久精品国产99精品国产亚洲性色| 久久精品久久久久久噜噜老黄 | 国产日韩欧美在线精品| 国产黄a三级三级三级人| 午夜免费男女啪啪视频观看| 免费大片18禁| 国产精品乱码一区二三区的特点| 亚洲高清免费不卡视频| 日本-黄色视频高清免费观看| 麻豆成人午夜福利视频| 欧美又色又爽又黄视频| 久久精品国产99精品国产亚洲性色| 91av网一区二区| 亚洲美女搞黄在线观看| 丰满人妻一区二区三区视频av| 亚洲三级黄色毛片| 国产一区二区三区在线臀色熟女| 日韩在线高清观看一区二区三区| 日本免费a在线| 此物有八面人人有两片| 久久久久性生活片| 国产一区二区三区av在线 | 久久精品夜夜夜夜夜久久蜜豆| 高清在线视频一区二区三区 | 亚洲精品成人久久久久久| 成熟少妇高潮喷水视频| 色吧在线观看| 看黄色毛片网站| 搡老妇女老女人老熟妇| 国产精品一区二区在线观看99 | 九草在线视频观看| 一本久久中文字幕| 成人毛片a级毛片在线播放| 少妇人妻精品综合一区二区 | 国产黄色小视频在线观看| 大香蕉久久网| 国产亚洲精品久久久com| 亚洲精品自拍成人| 如何舔出高潮| 亚洲av.av天堂| 国产精品野战在线观看| 岛国在线免费视频观看| 不卡一级毛片| 99久久无色码亚洲精品果冻| 在线观看免费视频日本深夜|