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

    Distributed output feedback stationary consensus of multi-vehicle systems in unknown environments

    2018-06-04 02:47:36VahidREZAEIMargaretaSTEFANOVIC
    Control Theory and Technology 2018年2期

    Vahid REZAEI,Margareta STEFANOVIC

    Department of Electrical and Computer Engineering,University of Denver,Denver,CO 80208,U.S.A.

    1 Introduction

    Control of vehicular systems has received a significant attention among policy makers and researchers during the past two decades due to the increased demand in transportation systems,advances in wireless communication devices,and embedded sensing and computation technologies(see[1,2])such that the market of autonomous vehicular systems will expectedly hit$42B by 2025 and,shortly after that,$85B by 2030[3].

    Cooperative analysis and control of multi-vehicle sys-tems have been done from both the systems-theoretic[4,5]and graph-theoretic viewpoints.In the first,the multi-vehicle system is usually considered over a standard string or mesh topology.Along with the advances in wireless and embedded technologies,graph-theoretic tools have created a promising alternative viewpoint in which the behavior of multi-vehicle systems can be analyzed over graphs where nodes represent vehicles and edges indicate inter-vehicle communication.This approach allows to consider more complicated topologies than the standard string or mesh of vehicles[6],and design the multi-vehicle cooperative algorithm independently of the vehicle-level controllers.Reference[7]developed a feedback linearization-like scheme to transform moving robot’s nonlinear dynamics to double integrator model with the goal of cooperative formation;[8]used double integrator models for the formation of unmanned vehicles;single integrators were used in[9]to model multi-robot systems,[6]proposed double integrators to study the relationship between communication topology and the stability of coordination algorithm,and[10]designed vehicle-level controllers for linear models of vehicles and a filtering-based cooperative algorithm for multi-vehicle systems.

    Motivated by their wide applications in the cooperation of multi-vehicle and multi-robot systems,significant theoretical research work has been devoted to the distributed control of first-and second-order agents[11–13].Other than multi-agent system of single integrators,it is known that distributed consensus algorithms usually result in a dynamic agreement in which all trajectories evolve during the time.References[14–17]introduced leaderless stationary consensus problem in which agents agree to stop at the same(fixed)position.Nevertheless,multi-vehicle systems are subject to unknown disturbances such as road profile[18,19]or wind[20–22]which may degrade the performance of consensus algorithm or destabilize it.Additionally,proposing a distributed algorithm with fewer absolute measurements than[14–17]may result in a less expensive multi-vehicle system.

    Distributed disturbance rejection algorithms have been investigated in the literature of multiagent systems to achieve conventional leaderless or leader-follower agreement among agents.In particular,motivated by the effectiveness of integral action in control systems theory,[23]and[24]studied integral action-based consensus algorithms for multi agent systems of single-and double-integrator agents under persistent disturbances,respectively.In order to consider a wider class of the multiagent systems,[25]showed the disturbance-free algorithms do not achieve conventional consensus in the presence of unknown eventually constant disturbances,and developed an observer-based leader-follower disturbance rejection consensus strategy for linear time in variant dynamical agents.Reference[26]generalized the result to guarantee consensus in the presence of a wider class of unknown disturbances(e.g.,with sinusoidal waveforms).However,because persistent disturbances continuously excite the uncontrollable and unobservable agreement dynamics,these disturbance rejection algorithms are not able to guarantee leaderless stationary consensus in multi-vehicle systems.

    Inspired by[26],we develop a dynamic output feedback leaderless stationary consensus algorithm based on the relative output information of vehicles and only a few vehicles’absolute measurements.We propose a framework to transform this dynamic output feedback problem into three low-order subproblems for disturbance rejection,consensus,and observer gain design tasks.Independently of the number of vehicles,consensus and observer gains are systematically found through two robust static state feedback formulations for loworder dynamics subject to fictitious modeling uncertainties.We further prove the proposed framework can be used to guarantee leader-follower stationary consensus in multi-vehicle systems(with a leader whose dynamics are not identical to the follower vehicles),and find analytical solutions for the consensus gains based on the design parameters and inter-vehicle communication graph.We verify the feasibility of proposed leaderless and leader-follower stationary consensus approaches in simulation.

    The rest of this paper is organized as follows.In Section 2,we introduce the notation and preliminary concepts.In Section 3,we provide the first contribution of this paper and propose a distributed leaderless stationary consensus algorithm that allows multi-vehicle systems to operate in unknown environments.In Section 4,we further prove the proposed framework for the leaderless scheme can be applied to the leader-follower stationary consensus in multi-vehicle systems.In Section 5,we discuss some challenges that may arise using conventional leaderless stationary algorithms in multivehicle systems,and verify the effectiveness of the proposed approaches through simulation studies.We summarize this paper in Section 6.

    2 Notation

    In this paper,Rn×ndenotes real-valued n × n matrices,Inidentity matrix with dimension n,0 a matrix of all zeros with appropriate dimension,1nan n×1 vector of all ones,diag{·}a(block)diagonal matrix with elements inside{·}as its diagonal terms,col{xi}a column vector of ordered vectors(or scalars)xifor i∈ {1,2,...,n},and ‖x‖Euclidean norm of the vector x∈Rn.Also,A?0(≥0)represents positive(semi)definite matrix A,and therow permutation matrix P ∈ Rn×nis an Inwith switched rows.For any symmetric matrix A ∈ Rn×n,there exists a unitary transformation matrix T ∈ Rn×nsuch that T?1AT= Λ is satisfied with Λ ? diag{λi}where λi∈ R are eigenvalues of A for i∈{1,2,...,n}[27].Moreover,the Rayleigh-Ritz inequality states λmin(A)‖x‖2≤ xTAx ≤ λmax(A)‖x‖2holds for any symmetric matrix A with λmin(A)and λmax(A)as its minimum and maximum eigenvalues.

    In the leaderless scenario,a set V of N nodes is associated to vehicles and the edge set E indicates the inter-vehicle communication topology which is visualized by undirected graph G.This information is abstracted by adjacency matrix A=[aij] ∈ RN×Nin which aij=aji,and aij=1 whenever the i th vehicle receives information from j th vehicle and aij=0 otherwise where i,j∈{1,2,...,N}.This characterization is valid for graphs with no self loop,i.e.,aii=0.The graph Laplacian matrix L ∈ RN×Nis further defined by Lij= ?aijforandThe following fact holds for any undirected graphs[28].

    In the leader-follower scenario,we introduce a set of N+1 nodes over graph Glfwhere N vehicles serve as followers and communicate over an undirected graph G,and only a few vehicles receive updates from the leader over aset of directed edges.The information about these leader-to-follower connections is lumped in B=diagwhere b=col{bi},bi=1 whenever the i th follower receives information from the leader,bi=0 otherwise,and i∈{1,2,...,N}.We partition the leader-follower graph Laplacian matrixas follows:

    where the first row,with all zeros,implies that the leader does not receive information from followers,and H=L+B ∈ RN×Nwhere L is a graph Laplacian matrix that abstracts the undirected communication topology among followers.Having knowledge about H,it is straightforward to find the leader-to-follower connections B and asymmetric positivedefinite L based on the fact that L 1N=0(see Fact 1).In this sense,we name H a reduced-order Laplacian matrix because it completely characterizes the leader-follower communication graph Glf,and use H(instead of Llf)to develop the ideas of this paper.The following fact holds[29].

    Fact 2There duced-order Laplacian matrix Hispositive definite if and only if the leader-follower graph Glfhas a spanning tree with the leader node i=0 as its root.

    3 Leaderless stationary consensus

    In this section,we develop a systematic framework to design output feedback dynamic stationary consensus algorithm for leaderless multi-vehicle systems.Compared to the existing results in the literature,which are disturbance-free problems with all vehicles’access to their absolute velocity measurements[14–17],the proposed strategy is based on the relative output measurements and a few vehicles’absolute output variables(potentially,only one vehicle),and ensures agreement in the presence of various unknown persistent disturbances.For such a multi-vehicle system with limited information,we first show the relative information-based dynamic consensus task can be reformulated as a set of local stability problems using heterogeneous absolute measurements.Then,in order to have a scalable design for multi-vehicle systems with a high-number of vehicles,we further recast it as three design sub-problems to find disturbance cancellation,robust feedback,and robust observer gains.

    3.1 Problem statement

    We consider a group of moving vehicles modeled by the following dynamics:

    where xi∈R denotes position,vi∈R velocity,ui∈R control input,andrelative output measurement of the i th vehicle.Also,g2∈R are two out putgain scalars which,for g2=0,reduces to the partial measurement scenario for multi-vehicle systems.Whenever i th vehicle has access to its absolute output measurement,we set bi=1 and,otherwise,bi=0.

    Fact 3The triple(C,A,B)is controllable and observable where

    The effects of unknown environment on the vehicles’dynamics are modeled by a set of heterogeneous persistent disturbances di∈R:

    where zi∈Rnzstands for the disturbance-state of the i th vehicle,unknown initial value of the disturbancestate,andand Cd∈ R1×nzare two known constant matrices that determine the shape of disturbance.For this multi-vehicle system,the leaderless stationary consensus task is defined as follows(adopted from[15]):

    which should be achieved in the presence of unknown persistent disturbances for all i,j∈{1,2,...,N}.

    Table 1 Examples of persistent disturbances in(2).

    3.2 Main result

    In order to guarantee the leaderless stationary consensus of multi-vehicle system(1)in the presence of heterogeneous persistent disturbances(2),we propose a dynamic distributed stationary consensus algorithm:

    where αv> 0 is a design scalar,and kcx,kcv∈ R and Kcd∈ R1×nzare the control gains to be determined later in this section.The i th vehicle’s estimated positionvelocityand disturbance state variableare found using the distributed observer:

    where the observer gains kov,kox∈ R and Kod∈ Rnz×1which will be designed in this section.Also,Nidenotes the neighboring set of the i th vehicle over an undirected graph G which satisfies the following assumption.

    Assumption 1The undirected communication graph G is connected(see Fact 1).

    We define observer error variablesand find the observer er-ror dynamics based on the relative measurements in vehicles’neighborhoods:

    Now,the augmented multi-vehicle system and observer error dynamics are given by

    where

    Based on the disturbance generator model(2)and Table 1,it is evident that the dimension nzof disturbance state variable zidepends on the shape of persistent disturbance;thus,wedefine z=col{col{zli},...,col{znzi}}∈RNnzand ez=col{col{ez1i},...,col{eznzi}}∈RNnzwhere

    Note that,in the proposed algorithm(4)and(5),depending on g2and bi,only a few vehicles have access to their absolute output measurements.Thus,the closedloop multi-vehicle system cannot be designed following the existing results in the literature of stationary consensus(see[14–17]).On the other hand,we know the definition of adjacency and Laplacian matrices does not admit self loops in the communication graph(see Section 2).Thus,we use the reduced-order(leader-follower)Laplacian matrix H in(6)to analytically handle this situation.

    We now take the diagonal blocks in(6)and find two sets of differential equations corresponding to multivehicle and observer error dynamics.In particular,˙e=A22e results in the following transformed aggregated observer dynamics:

    in whichandis a unitary transformation matrix that completely diagonalizes the symmetric matrix H aswhere μi> 0 are eigenvalues of H?0 for i∈{1,2,...,N}.We notice that these aggregated error dynamics are composed by N “networked”observer error models:

    where μiresult in heterogeneity of the closed-loop networked observers(with kox,kovand Kodin the loop).The word “networked”refers to the effect of inter-vehicle communication topology G on the localized models(7)throughμi.Moreover,using the firstrow˙ζ=A11ζ+Bd1z of(6),the multi-vehicle system is written as follows:

    in whichwhereis a unitary transformation that completely diagonalizes the symmetric Laplacian matrix L such thatThe eigenvalue λ1=0 in the diagonal matrix Λldetermines the null space of graph Laplacian matrix L(see Fact 1),and results in a subsystem that does not satisfy the controllability condition.We apply another transformation ξ=Pζ where ζ =[xT,vT]Tand P ∈ R2N×2Nis a row switching matrix that results inandfor i∈ {2,3,...,N}.The subscriptastands for agreement subspace which is uncontrollable,anddindicates controllable disagreement subspace.Using a matrix representation,the transformed system is written as follows:

    in which the agreement and disagreement dynamics are decoupled from each other.The uncontrollable agreement dynamics show their effects on the stationary consensus value,and will be discussed in Corollary 1 and Lemma 1.However,the controllable disagreement dynamics are made by N?1 heterogeneous networked vehicle models:

    for all i∈{2,3,...,N}where the heterogeneity is due to the nonzero eigenvalues λiof L,and the word “networked”is used with the same interpretation as in(7).We now use(7)and(10),and characterize the relative output feedback stationary consensus problem(3)as output feedback stabilization tasks based on the heterogeneously-scaled absolute measurements(due to μiand λi),and establish three conditions for consensus gains to ensure stationary agreement(3)in multi-vehicle system(1)using algorithm(4)and(5).

    Proposition 1The dynamic distributed algorithm(4)–(5)guarantees leaderless stationary consensus(3)among vehicles(1)in the presence of persistent disturbances(2)whenever Assumption 1 is satisfied and the following vehicle-level conditions are guaranteed:

    1)The disturbance control gain Kcdshould accommodate the effect of unknown heterogeneous disturbances on the networked vehicle dynamics(10).

    3)For an arbitrarily selected αv> 0,the control gain Kc=[kcx,kcv]should stabilize networked vehicle dynamics:

    with N?1 heterogeneous measurements where heterogeneity is the effect of positive eigenvalues λiof L.(The effect of αvon the final position of vehicles will be discussed at the end of this section).

    The proof of this proposition is immediate based on the aforementioned derivations and noting the fact that,satisfying Part 1 of Proposition 1,separation principle holds for designing observer and controller gains(see the structure of(6)).Although the first part of this proposition is independent of the multi-vehicle system’s dimension,we need to examine the observer gains kox,kovand Kodfor N nonzero eigenvalues μiof H,and control gains kcx,kcvfor N ?1 nonzero eigenvalues λiof L.This puts question on the feasibility of using Proposition 1 for multi-vehicle systems with a high-number of vehicles.

    In the literature,various viewpoints have been proposed to overcome such a problem.The dominant approaches are based on introducing a scalar correction factor into the consensus algorithm:references[30]and[31]introduced a correction factor to modify the algebraic Riccati equation-based control gain and guarantee leaderless consensus in multiagent systems(e.g.,based on the formulation of this paper,we need to implement ckcx,ckcvin consensus signal(4)where c is a correction gain that should be designed independent of the actual gains kcxand kcv),and[32]used the correction gain to modify its relative measurements(e.g.,in this viewpoint we need to use the modified relative measurements cyriin(4)and(5)with a modification factor c).In the rest of this section,we follow a recent viewpoint on addressing this challenge by reformulating the static output feedback problems(11)and(12)as state-feedback robust control challenge for the nominal networked vehicles’models subject to fictitious uncertainties(e.g.,see[26]and[33]),and provide a systematic approach to find appropriate control and observer gains that guarantee stationary consensus for the multivehicle systems operating in unknown environments.At first,we investigate the control design problem in Parts 1 and 3 of Proposition 1 and,later,we will discuss the observer gain design problem.

    Design procedure 1The gains Kcd,kcx,and kcvshould be designed as follows:

    1)Disturbance control gain Kcdshould minimize‖(Kcd+Cd)zTi‖where,based on the definition of disturbance generator model(2),Kcd=?Cdaccommodates all persistent disturbances.

    2)Robust state feedback gains kcxand kcvshould be designed to stabilize networked robot dynamics:

    where E(λi)= λi/λ2? 1 ≥ 0 are(communication graph-induced)factitious modeling uncertainties for i∈{2,3,...,N}.

    We emphasize that all eigenvalues λimight be known due to the knowledge about communication topology G,but we only use the algebraic connectivity λ2and consider the rest as the sources of modeling uncertainties in order to propose a one-step design procedure and find consensus gains in(4)(compared to the two-step methods in[30–32]).This is the reason to call λ2E(λi)uTi“fictitious”modeling uncertainties.Now,we define a second-order state space model˙ξi=Aξi+BuTifor the nominal networked vehicle dynamics in(13):

    and letbe a design matrix and rc>0 be a design scalar to respectively weigh state and control input variables.In the next theorem,we systematically derive two static consensus gains kcxand kcvthat stabilize uncertain dynamics in Step 2)of Design procedure 1 for all i∈{2,3,...,N}.

    Theorem 1Assuming uTi=Kcξihas been found that achieves the minimum of quadratic cost function(15)where UTiis the set of all(admissible)stabilizing signals uTi.Then,Kc=[kcx,kcv]is the required gain to stabilize uncertain networked-vehicle dynamics(13).

    Furthermore,implementing uTi=Kcξi,we know any pairs(ξi,uTi)satisfy two fundamental properties of optimal control theory(e.g.,see Section 3.11 in[34]):

    which satisfies V(ξi(0))=J(ξi(0))for any initial conditions,and find its time deviation along the uncertain trajectories in(13):

    which proves exponential stability of the origin for networked vehicle systems(13)(see Thm.3.6 in[35]).□

    Now,we introduce the networked observer’s nominal dynamics:

    where Aois defined in(11)andAlso,letbe a design matrix and ro>0 be a design scalar.In the next theorem,we systematically find the required observer gainto be used in(5).

    Theorem 2Findthat achieves the minimum value of(19)subject to completely known dynamical system where Uoidenotes a set of admissible control signals for uoi.Then,Kois the required observer gains for dynamical system(11)in Step 2)of Proposition1.

    ProofA sketch of this poof can be given by noticing that the dynamical system in(19)is dual to(18).We can similarly find the dual representation for(11).Then,this theorem is proved following the steps in the proof of Theorem 1 for the dual problem and based on ARE:

    where the unique stabilizingexists by verifying observability and controllability offor□

    Corollary 1Assuming perfect state and disturbance measurements,the observer-free consensus algorithm(4)will result in the following stationary agreement values:

    ProofBased on(9),we know agreement dynamics are decoupled from the disagreement dynamics(respectively determined by ξaand ξd).The agreement dynamics are modeled byandand we write the solution of these differential equations asandwhich result in the following limit behavior:

    because αv> 0 andas(this αvis a design scalar and can be tuned to achieve desirable consensus behavior).Note that the superscriptadenotes the“agreement”value as.Letwheretion such that(see Fact 1).If we rewrite this result based on the agreement and disagreement variables and useis the diagonalizing unitary trans for ma-(based on Theorem 1),we find

    in which

    Lemma 1The observer-based output feedback stationary consensus algorithm(4)and(5)results in an agreement on the following point in the presence of unknown persistent disturbances:

    ProofWe begin from the augmented multi-vehicle and observer error dynamics(6),substitute Kcdby?Cd,and find

    in which

    Based on a row switching transformation Pl=diag{P,whereis defined in(9),we write the disagreement dynamics as follows:

    The solution of second equation iswhere,as∞,the integral converges to a constant βI∈ R(because the error variables go to zero).Thus,αv> 0 results in

    is also a constant value. □

    In summary,based on the formulation(12)(also(13)),Theorem 1 guarantees the stability of disagreement dynamics in(9)is achieved for any “arbitrarily”selected αv> 0,and Lemma 1 shows this αvadds a level of flexibility to tune the internal behavior of multi-vehicles’agreement dynamics(although they remain“cooperatively”uncontrollable according to the partitioning in(9)).Additionally,as expected,the last term inindicates that agreement value depends on the average of velocity and disturbance state estimation errors’transient behavior.This will be discussed in the simulations of Section 5.

    4 Leader-follower stationary consensus

    In Section 3,we have developed an output feedback leaderless stationary consensus algorithm for multivehicle systems modeled by double integrators in which,unlike[14–17],only a few vehicles have access to their absolute output measurements and all of them are subject to unknown persistent disturbances.Now,we adopt the result of Section 3 and develop a systematic approach to design a leader-follower stationary algorithm.In particular,we find the closed-form solution for consensus gains purely based on the design matrices and smallest positive eigenvalue of the reduced-order Laplacian matrix corresponding to the leader-follower communication graph.

    4.1 Problem statement

    In this section,we consider a multi-vehicle system where the followers are modeled by second-order dynamics:

    for i∈{1,2,...,N},bi=1 whenever the i th vehicle is aware of its relative distance to the reference position and potentially is aware of its absolute velocity(depending on g2),and bi=0 otherwise.All variablesare defined similar to the leaderless consensus problem in Section 3,and we emphasize that g1≠0.The reference point is commanded by a stationary leader:

    which,unlike followers,is described by a first-order model and its adjustable initial state value x0(t)=x0(0)denotes the desired position.Now,we define the leader follower stationary consensus as follows:

    where t0≥0 denotes the time of change in the reference command.Before proposing the main result of this section,we make an assumption on the leader-follower communication graph topology Glf.

    Assumption 2The leader-follower graph Glfhas a spanning tree with the leader node i=0 as the root(see Fact 2).

    4.2 Main result

    In this section,we propose a dynamic distributed leader-follower stationary consensus algorithm in order to ensure agreement(24)in multi-vehicle systems(22)and(23):

    whererespectively denote position and velocity,andrepresents the disturbance state of the i th vehicle estimated by the following distributed observer:

    We define observer error variablesand find observer error dynamics

    where

    and the sub-matrices are as follows:

    Based on the augmented system(27),we conclude the separation principle holds,and consensus and observer gains can be designed independent of each other.For the observer design purpose,we find the following networked observer dynamics:

    which shows the proposed leader-follower stationary consensus algorithm(25)and(26)has resulted in the same problem as the leaderless consensus scenario in(11).For the control gain design problem,since H is a symmetric positive-definite matrix,we find a completely controllable diagonal representation:

    based onandwhereis defined such thatdiag{μi}for i∈ {1,2,...,N}.The transformed dynamics(29)are in fact composed by N heterogeneous networked leader-follower tracking error systems:

    in which μi> 0 for all i∈ {1,2,...,N}.Therefore,the following proposition holds in this section.

    Proposition 2Suppose Assumption 2 is satisfied by the communication graph Glf.The dynamic distributed algorithm(25)–(26)ensures leader-follower stationary agreement(24)in multi-vehiclesy stems(22)in the presence of unknown disturbances(2)whenever,in addition to Steps 1 and 2 of Proposition 1,the control gain Kc=[kcxkcv]stabilizes the networked leader-follower tracking error dynamics:

    whereμi> 0are the eigenvalues of reduced-order Laplacian matrix H for all i∈{1,2,...,N}.

    This proposition,along with the observer dynamics(28)and disturbance component in(30),shows the distributed stationary leader-follower algorithm(25)–(26)is formulated such that the disturbance control gain can be found following Step 1)in Design procedure 1,and the observer gain can be designed using Theorem 2 in Section 3.In the next design procedure,we propose a systematic framework to find the consensus gains kcxand kcvbased on a robust control formulation for modified networked leader-follower tracking error dynamics with a homogeneous nominal part and heterogeneous fictitious modeling uncertainties.

    Design procedure 2Design state feedback gains kcx,kcv∈R that stabilize the networked vehicle dynamics with homogeneous nominal model and heterogeneous fictitious modeling uncertainties:

    Note that we originally proposed a dynamic output feedback stationary consensus algorithm(25)and(26)using relative measurements;converted it to three subproblems to design consensus,observer,and disturbance gains in Proposition 2 where the consensus gains kcxand kcvwere the stabilizing solutions for N static“output feedback”networked vehicles using N scaled absolute measurements;and,eventually,reformulated the problem as N static“state feedback”robust stabilization tasks using vehicles’absolute state measurements(Design procedure 2).Now,we introduce a secondorder state space realization˙ξi=Aξi+BuTiwhich models the nominal dynamics of(32):

    Furthermore,we defineas the state weighting,and rc>0 as the control input weighting design matrices.In the next theorem,we propose a systematic framework to find the required consensus gains as a single robust state feedback problem.

    Theorem 3The solution uTi=Kcξito the minimization problem(34),where UTidenotes the set of all stabilizing state feedback controllers uTi,stabilizes the heterogeneous networked vehicle dynamics(32)for all i∈{1,2,...,N}.

    ProofWe mention that,although system matrices(A,B)in(34)is different from(15),the minimization problems are structurally the same such that the fundamental properties(17)are still valid for any pairs ξi,uTiof this leader-follower control theorem.Therefore,a detailed proof can be found by following the steps of proof in Theorem 1.This is omitted for brevity. □

    which result in a unique positive-definite stabilizing solution:

    Consequently,we find closed-form solutions for the consensus gains kcxand kcvexplicitly based on the design matrix Qc,scalar rc,and smallest eigenvalue of reduced-order Laplacian matrix μ1of H:

    using the optimal gain formula

    We note that Pc(1,1)>0 is guaranteed based on the observability and stabilizability ofA,B).As is seen,the position consensus gain kcxis not directly affected by the communication network G(or μ1);however,a network-dependent fraction of it appears in the velocity consensus gain kcv(i.e.,seethat is added toThese closed-form solutions can be used in tuning of the weighting matrices Qcand rc,and also for the communication topology design purpose.We mention that the tuning process can be further simplified by letting q12=q21=0 or using Qc=qcI2for a scalar tuning parameter qc>0.

    5 Simulation verification

    In this section,we set up leaderless and leader follower case studies to be used in the simulation verifications,review some major challenges in conventional(non-stationary)leaderless consensus problem through simulation,verify the feasibility of proposed stationary leaderless algorithm of Section 3 and leader-follower approach of Section 4.

    5.1 Problem setup

    In the leaderless problem,we consider a group of five vehicles modeled by(1)with g1=1 and g2=0,and assume nodes 1 and 2 have access to their absolute position information(i.e.,b1=b2=1).Vehicles are at initial conditions x1(0)=[?10,20]T,x2(0)=[15,?15]T,x3(0)=[10,15]T,x4(0)=[?30,20]T,and x5(0)=[20,?30]T(which are unknown to the designer).Moreover,vehicles are subject to heterogeneous constant disturbances d1=2,d2=5,d3=3,d4=9 and d5=4.We let vehicles to exchange information over the leaderless graph in Fig.1.

    Fig.1 Leaderless communication:By removing node v0 and edges originating from it,nodes v1–v5 and the associated edges represent an undirected leaderless communication topology G with graph Laplacian matrix L.We assume agents v1 and v2 have access to their absolute measurements.Leader-Follower communication:Agents v1–v5 communicate over undirected graph G,nodes v1 and v2 are aware of their relative distance to the leader vehicle v0 such that the entire multi-vehicle system v0–v5 builds a leader-follower communication graph G lf with a reduced-order Laplacian matrix H=L+B where B=diag{1,1,0,0,0}.

    In the leader-follower scenario,we add a leader agent v0modeled by(23)where its initial value can be commanded globally(and we will use a square wave input in simulation),and let vehicles to communicate over Glfin Fig.1.Moreover,vehicles are subject to sinusoidal disturbances:d1=7sin(0.5t),d2=5.5sin(0.5t),d3=6sin(0.5t),d4=2sin(0.5t),and d5=4sin(0.5t).

    5.2 Overview of the literature

    We overview the problems that may arise using the existing conventional leaderless consensus algorithms with disturbance rejection capability.In particular,we first upgrade the leaderless consensus algorithm of[26]to be able to include unstable vehicle dynamics,and find the following agreement and final values using an observer-based approach over the leaderless graph G of Fig.1:

    where we assume all observers are at initial rest condition except?x3(0)=[10,?15]T,and vehicles are subject to unknown constant disturbances(see Section 5.1).A simulation result is depicted in Fig.2 which verifies all vehicles reach to the(upgraded)agreement values in(35).

    Fig.2 Simulation results using an upgraded leaderless consensus algorithm of[26](to include unstable systems)subject to unknown constant disturbances.All observers are at initial rest conditions except?x3(0).As is seen,all position and velocity variables,and their estimations are unbounded.However,the estimations differ from the actual variables due to the nullity in observer dynamics.

    However,we notice a practically undesirable collective behavior where vehicles continuously increase their speed because of the nullity in system and disturbance state estimations(which also results in a mismatch between the actual variables and their estimations).In the next simulation,we assume all observers are at initial rest conditions and find that,independent of vehicles’initial conditionsfor all i∈{1,2,...,5}.As is depicted in Fig.3,although the observer variables remain bounded,the nullity(unobservability)of leaderless consensus is still problematic in the presence of persistent disturbances,and results in an uncontrolled increase in vehicles’speed.

    Fig.3 Simulation results for an upgraded leaderless consensus algorithm of[26]under the same scenario as Fig.2 except,now,all observers are at initial rest conditions.As is shown,although estimation variables converge to zero,the vehicles’consensus behavior remains the same as Fig.2.

    We further modify the leader-follower result of[26]to beableto in cludeun stable agent dynamics and,in Fig.4,show that all vehicles v1–v5agree on the “second-order”leader’s velocity and follow its position with constant speed.

    Fig.4 Simulation results using an upgraded leader-follower consensus algorithm[26](to include unstable systems)subject to unknown constant disturbances as in Figs.2 and 3.

    5.3 Leaderless stationary consensus

    In this section,we verify the feasibility of the proposed observer-based leaderless stationary consensus algorithm in Section 3.We show that the proposed strategy overcomes the undesirable increasing speed problem of the conventional leaderless algorithm in Section 5.2.The multi-vehicle system is implemented under the same scenario as in Section 5.1 over the leaderless graph of Fig.1.

    At first,in Fig.5,we consider an observer-free algorithm and verify that vehicles velocity reach to zero and they agree on the positionas expected by Corollary 1.In Fig.6,we use the proposed observer-based algorithm(4)and(5)where all observers are at initial rest condition,and show all vehicles reach to a“stationary”agreement on their position at.Moreover,unlike the conventional leaderless scenario of Fig.3,estimations are the same as actual position and velocity variables of vehicles.We further note that the difference in agreement position values of Figs.5 and 6 is expected based on Lemma 1.

    Fig.5 Observer-free leader less stationary consensus algorithm based on Corollary 1 in the presence of constant disturbances.All vehicles agree on a fixed position

    Fig.6 Leaderless stationary consensus algorithm of Section 3 where all observers are at initial rest condition and the agreement is on.The dashed line show the agreement value of Fig.5(see Lemma 1).

    In Fig.7,we initialize the second observer at[?10,15,0]Tand show the effect of observer error(trajectories)on the consensus value where,unlike Fig.6,the agreement is on.Finally,based on Theorem 2,we redesign observer gains by setting the state weighting matrix to be 103greater than the first design(see Fig.6),and find a new agreement onwhich is shown in Fig.8.In all of these simulation scenarios,disturbances are estimated precisely as is depicted in Fig.9.

    Fig.7 Leaderless stationary consensus algorithm with all observers are at initial rest condition except.Different from Fig.6,the agreement is on.This shows the effect of observer error trajectories on agreement values(see Lemma1).

    Fig.8 Leaderless stationary consensus under the same configuration as Fig.6 with a new set of design matrices.The agreement is on.This verifies the effect of observer dynamics(error trajectories)on the stationary agreement value.

    Fig.9 In all leader less stationary consensus simulations,disturbances are estimated precisely(with differences in the transient behavior).Top to bottom:w1 to w5(black)and their estimations(red).

    5.4 Leader-follower stationary consensus

    In this section,we investigate the effectiveness of leader-follower stationary algorithm(25)and(26)in ensuring a stationary consensus on a adjustable position with limited information about the leader(as discussed in Section 4).We consider a leader-follower setup subject to sinusoidal disturbances,as introduced in Section5.1.For this setup,we use a square wave to command the desired position of vehicles.Although only two vehicles are aware of their relative distances to the desired reference point,Figs.10 and 11 show all vehicles precisely estimate their positions,velocities,and disturbances;and agree on the commanded stationary point.

    Fig.10 Leader-follower stationary consensus:State variables and their estimations.The dashed back waves represent the leader’s command.

    Fig.11 Leader-follower stationary consensus.Top to bottom are w1 to w5(black)and their estimations(red).

    6 Summary

    We design leaderless and leader-follower stationary consensus algorithms which ensure all vehicles’agreement on a fixed point in the presence of unknown persistent disturbances with only a few vehicles’absolute measurements.In both leaderless and leader follower scenarios,we provide a systematic framework that transform the high-order dynamic relative-output feedback stationary consensus challenge to three loworder subproblems to be used in designing disturbance,consensus,and observer gains.We formulate the consensus and observer gain design tasks as two robust static feedback problems for modified vehicle dynamic subject to fictitious modeling uncertainties which are induced by communication graph topology.In simulation,we discuss the challenges of applying non-stationary disturbance rejection algorithms to multi-vehicle systems,and verify the feasibility of using proposed strategies for multi-vehicle systems in unknown environments where vehicles might be subject to road profile and wind disturbances.

    [1]D.Cruz,J.McClintock,B.Perteet,et al.Decentralized cooperative control.IEEE Control Systems Magazine,2007,27(3):58–78.

    [2]L.Bento,R.Parafita,U.Nunes.Inter-vehicle sensor fusion for accurate vehicle localization supported by V2V and V2I communications.International Conference on Intelligent Transportation Systems,Anchorage:IEEE,2012:907–914.

    [3]R.Langari.Autonomous vehicles:a tutorial on research and development issues.American Control Conference,Seattle:IEEE,2017:4018–4022.

    [4]P.Seiler,A.Pant,J.Hedrick.Disturbance propagation in vehicle strings.IEEE Transactions on Automatic Control,2004,49(10):1835–1842.

    [5]Y.Zhao,P.Minero,V.Gupta.On disturbance propagation in leader-follower systems with limited leader information.Automatica,2014,50(2):591–598.

    [6]W.Ren,E.Atkins.Distributed multi-vehicle coordinated control via local information exchange.International Journal of Robust and Nonlinear Control,2007,17(10/11):1002–1033.

    [7]J.Lawton,R.Beard,B.Young.A decentralized approach to formation maneuvers.IEEE Transactions on Robotics and Automation,2003,19(6):933–941.

    [8]R.Xue,G.Cai.Formation flight control of multi-UAV system with communication constraints.Journal of Aerospace Technology Management,2016,8(2):203–210.

    [9]A.Jadbabaie,J.Lin,S.Morse.Coordination of groups of mobile autonomous agents using nearest neighbor rules.IEEE Transactions on Automatic Control,2003,48(6):988–1001.

    [10]T.Yucelen,E.Johnson.Control of multi-vehicle systems in the presence of uncertain dynamics.International Journal of Control,2013,86(9):1540–1553.

    [11]W.Ren,R.Beard,E.Atkins.A survey of consensus problems in multi-agent coordination.American Control Conference,Portland:IEEE,2005:1859–1864.

    [12]W.Yu,G.Chen,C.Ming.Some necessary and sufficient conditions for second-order consensus in multi-agent systems.Automatica,2010,46(6):1089–1095.

    [13]N.Huang,Z.Duan,G.Chen.Some necessary and sufficient conditions for consensus of second-order multi-agent systems with sampled position data.Automatica,2016,63:148–155.

    [14]C.-L.Liu,F.Liu.Stationary consensus of heterogeneous multiagent systems with bounded communication delays.Automatica,2011,47(9):2130–2133.

    [15]J.Qin,C.Yu,S.Hirche.Stationary consensus of asynchronous discrete-time second-order multi-agent systems under switching topology.IEEE Transactions on Industrial Informatics,2012,8(4):986–994.

    [16]Y.Feng,S.Xu,F.Lewis,B.Zhang.Consensus of heterogeneous first-and second-order multi-agent systems with directed communication topologies.International Journal of Robust and Nonlinear Control,2015,25(3):362–375.

    [17]Y.Pei,J.Sun.Necessary and sufficient conditions of stationary average consensus for second-order multi-agent systems.International Journal of Systems Science,2016,47(15):3631–3636.

    [18]C.Gohrle,A.Schindler,A.Wagner,et al.Road profile estimation and preview control for low-bandwidth active suspension systems.IEEE/ASME Transactions on Mechatronics,2015,20(5):2299–2310.

    [19]I.Youn,M.Khan,N.Uddin,et al.Road disturbance estimation for the optimal preview control of an active suspension system based on tracked vehicle model.International Journal of Automotive Technology,2017,18(2):307–316.

    [20]A.Cho,S.Kim,C.Kee.Wind estimation and airspeed calibration using a UAV with a single-antenna GPS receiver and pitot tube.IEEE Transactions on Aerospace and Electronic Systems,2011,47(1):109–117.

    [21]B.Arain,F.Kendoul.Real-time wind speed estimation and compensation for improved flight.IEEE Trans actions on Aerospace and Electronic Systems,2014,50(2):1599–1606.

    [22]H.Shen,N.Li,H.Griffiths,et al.Tracking control of a small unmanned air vehicle with airflow awareness.American Control Conference,Seattle:IEEE,2017:4153–4158.

    [23]T.Yucelen,M.Egerstedt.Control of multiagent systems under persistent disturbance.American Control Conference,Montreal:IEEE,2012:5264–5269.

    [24]M.Anderson,D.Dimarogonas,H.Sand berg,et al.Distributed control of networked dynamical systems:static feedback integral action and consensus.IEEE Transactions on Automatic Control,2014,59(7):1750–1764.

    [25]W.Cao,W.Zhang,W.Ren.Leader-follower consensus of linear multi-agent systems with unknown external disturbance.Systems&Control Letters,2015,82:64–70.

    [26]V.Rezaei,M.Stefanovic.Distributed leaderless and leader follower consensus of linear multiagent systems under persistent disturbances.Mediterranean Conference on Control and Automation,Athens,Greece:IEEE,2016:386–391.

    [27]R.Horn,C.Johnson,Matrix Analysis.Cambridge:Cambridge University Press,2013.

    [28]M.Mesbahi,M.Egerstedt.Graph Theoretic Methods in Multiagent Networks.Princeton:Princeton University Press,2010.

    [29]W.Ni,D.Cheng.Leader-following consensus of multi-agent systems under fixed and switching topologies.Systems&Control Letters,2010,59(3/4):209–217.

    [30]S.-E.Tuna.LQR-based coupling gain for synchronization of linear systems.arXiv,2008:arXiv:0801.3390.

    [31]H.Zhang,F.Lewis,A.Das.Optimal design for synchronization of cooperative systems:state feedback,observer and output feedback.IEEE Transactions on Automatic Control,2011,56(8):1948–1952.

    [32]Z.Li,Z.Duan,G.Chen,et al.Consensus of multiagent systems and synchronization of complex networks:a unified viewpoint.IEEE Transactions on Circuits and Systems I:Regular Papers,2010,57(1):213–224.

    [33]V.Rezaei,M.Stefanovic.Distributed decoupling of partially unknown interconnected linear multiagent systems:state and output feedback approaches.IFAC-PapersOnLine,2017,50(1):1766–1771.

    [34]D.Kirk.Optimal Control Theory:An Introduction.Englewood Cliffs:Prentice Hall,1970.

    [35]H.Khalil.Nonlinear Control.Englewood Cliffs:Prentice Hall,2014.

    欧美三级亚洲精品| 久久久久国产一级毛片高清牌| 麻豆国产av国片精品| 国产熟女午夜一区二区三区| www国产在线视频色| 久久这里只有精品19| 欧美色视频一区免费| 1024手机看黄色片| www国产在线视频色| √禁漫天堂资源中文www| 欧美丝袜亚洲另类 | 黑人巨大精品欧美一区二区mp4| 91老司机精品| 久久久久性生活片| 国模一区二区三区四区视频 | 国产伦人伦偷精品视频| 叶爱在线成人免费视频播放| 19禁男女啪啪无遮挡网站| 老汉色av国产亚洲站长工具| 99久久精品热视频| 一区福利在线观看| 日韩成人在线观看一区二区三区| 亚洲片人在线观看| 久久久久久久久久黄片| 少妇裸体淫交视频免费看高清 | 国产精品亚洲av一区麻豆| 99久久综合精品五月天人人| 日韩欧美国产在线观看| 亚洲av片天天在线观看| 女同久久另类99精品国产91| 亚洲精品一区av在线观看| av国产免费在线观看| 国内揄拍国产精品人妻在线| 精品国产美女av久久久久小说| 少妇被粗大的猛进出69影院| 国产精品av视频在线免费观看| 美女高潮喷水抽搐中文字幕| 久久香蕉国产精品| 久久中文字幕人妻熟女| 精品久久久久久久末码| 最近在线观看免费完整版| 成年免费大片在线观看| 国产精品久久电影中文字幕| 黄色a级毛片大全视频| 亚洲午夜理论影院| 老司机在亚洲福利影院| 少妇粗大呻吟视频| 国产又色又爽无遮挡免费看| 成人午夜高清在线视频| 99热这里只有精品一区 | 一级黄色大片毛片| av欧美777| 人成视频在线观看免费观看| 国产成人啪精品午夜网站| 亚洲色图av天堂| 国产一区在线观看成人免费| 免费人成视频x8x8入口观看| 国产成人啪精品午夜网站| 老熟妇乱子伦视频在线观看| 黄色女人牲交| 禁无遮挡网站| 精品国内亚洲2022精品成人| 国内少妇人妻偷人精品xxx网站 | 久99久视频精品免费| 精品国产乱子伦一区二区三区| 动漫黄色视频在线观看| 变态另类丝袜制服| 精品熟女少妇八av免费久了| 在线国产一区二区在线| 午夜福利在线在线| 国产精品av久久久久免费| 中文字幕精品亚洲无线码一区| 天堂av国产一区二区熟女人妻 | 欧美一区二区国产精品久久精品 | 成年人黄色毛片网站| 久久香蕉激情| 女警被强在线播放| 久久草成人影院| 老司机午夜十八禁免费视频| 亚洲成人久久爱视频| 成人三级做爰电影| 久久精品国产99精品国产亚洲性色| 亚洲欧美精品综合久久99| 一本精品99久久精品77| 国产av麻豆久久久久久久| 国产精品一区二区精品视频观看| 亚洲国产欧美网| 欧美乱色亚洲激情| 麻豆成人av在线观看| 又粗又爽又猛毛片免费看| 宅男免费午夜| 日本a在线网址| 后天国语完整版免费观看| 男插女下体视频免费在线播放| 首页视频小说图片口味搜索| 精品一区二区三区av网在线观看| 国产精品乱码一区二三区的特点| 国产爱豆传媒在线观看 | 国产成人欧美在线观看| xxx96com| 亚洲一区中文字幕在线| 日本黄大片高清| 国产不卡一卡二| 日韩av在线大香蕉| 午夜久久久久精精品| 国产午夜福利久久久久久| 天堂√8在线中文| 操出白浆在线播放| 国产亚洲精品一区二区www| 日韩欧美精品v在线| а√天堂www在线а√下载| 色综合欧美亚洲国产小说| 曰老女人黄片| 岛国视频午夜一区免费看| 后天国语完整版免费观看| 国产激情偷乱视频一区二区| 在线视频色国产色| 99久久无色码亚洲精品果冻| 久久久精品大字幕| 神马国产精品三级电影在线观看 | 国产精品免费视频内射| 亚洲精品美女久久av网站| 久久久精品国产亚洲av高清涩受| 老司机深夜福利视频在线观看| 美女免费视频网站| 又大又爽又粗| 久久 成人 亚洲| av天堂在线播放| 制服丝袜大香蕉在线| 欧美中文综合在线视频| 亚洲中文字幕一区二区三区有码在线看 | av免费在线观看网站| 国产欧美日韩一区二区三| 久久午夜亚洲精品久久| 亚洲国产日韩欧美精品在线观看 | 日韩欧美 国产精品| 久久99热这里只有精品18| 亚洲人成77777在线视频| 久久人妻福利社区极品人妻图片| 国产高清视频在线观看网站| 亚洲性夜色夜夜综合| 久久精品成人免费网站| 18禁黄网站禁片午夜丰满| 成人18禁高潮啪啪吃奶动态图| 亚洲专区字幕在线| 国产精品综合久久久久久久免费| 99久久精品热视频| 国产伦人伦偷精品视频| 国产成人aa在线观看| 国产aⅴ精品一区二区三区波| 老司机深夜福利视频在线观看| 亚洲专区中文字幕在线| 久久天堂一区二区三区四区| 亚洲av第一区精品v没综合| 狂野欧美白嫩少妇大欣赏| 最近最新中文字幕大全电影3| 色哟哟哟哟哟哟| 90打野战视频偷拍视频| 国产久久久一区二区三区| 妹子高潮喷水视频| 免费看美女性在线毛片视频| 中文字幕最新亚洲高清| 此物有八面人人有两片| 夜夜爽天天搞| 国产av麻豆久久久久久久| 精品久久久久久,| 中亚洲国语对白在线视频| 看免费av毛片| 中文字幕人成人乱码亚洲影| 真人一进一出gif抽搐免费| 亚洲国产欧美网| 日韩精品免费视频一区二区三区| 1024手机看黄色片| 又黄又粗又硬又大视频| av天堂在线播放| 久久久久久久久免费视频了| 国产成+人综合+亚洲专区| 999精品在线视频| 中文字幕人成人乱码亚洲影| 19禁男女啪啪无遮挡网站| a级毛片在线看网站| x7x7x7水蜜桃| 99久久99久久久精品蜜桃| 成人特级黄色片久久久久久久| 18禁黄网站禁片免费观看直播| 国产高清视频在线观看网站| 欧美在线一区亚洲| 97超级碰碰碰精品色视频在线观看| 老熟妇乱子伦视频在线观看| 视频区欧美日本亚洲| 人妻久久中文字幕网| 久久久久久免费高清国产稀缺| or卡值多少钱| 国产又黄又爽又无遮挡在线| 婷婷精品国产亚洲av在线| 69av精品久久久久久| 午夜精品一区二区三区免费看| 久久久国产欧美日韩av| a在线观看视频网站| 波多野结衣巨乳人妻| 国产真人三级小视频在线观看| 亚洲成av人片免费观看| 少妇熟女aⅴ在线视频| 丰满人妻一区二区三区视频av | www日本黄色视频网| 一本一本综合久久| 曰老女人黄片| 亚洲av电影在线进入| 亚洲第一电影网av| 精品人妻1区二区| 久久精品成人免费网站| 无限看片的www在线观看| 18禁黄网站禁片免费观看直播| 叶爱在线成人免费视频播放| 非洲黑人性xxxx精品又粗又长| 两个人看的免费小视频| 曰老女人黄片| 99re在线观看精品视频| 18禁美女被吸乳视频| 国产97色在线日韩免费| 国产午夜精品论理片| 亚洲自拍偷在线| 免费搜索国产男女视频| 久久国产精品人妻蜜桃| 久久久久久九九精品二区国产 | 91老司机精品| 国产精品国产高清国产av| 精品电影一区二区在线| 久久婷婷人人爽人人干人人爱| 国产一区二区激情短视频| 国产av又大| 丁香六月欧美| 啦啦啦韩国在线观看视频| 嫩草影视91久久| 国产精品九九99| 亚洲欧美日韩高清专用| 天堂av国产一区二区熟女人妻 | 色噜噜av男人的天堂激情| 日韩av在线大香蕉| 午夜精品一区二区三区免费看| www.自偷自拍.com| 国内毛片毛片毛片毛片毛片| 久久久国产成人精品二区| 一二三四在线观看免费中文在| 久久久久久久久久黄片| 国产成人aa在线观看| netflix在线观看网站| 两个人的视频大全免费| 精品乱码久久久久久99久播| 美女 人体艺术 gogo| 免费一级毛片在线播放高清视频| 欧美一区二区精品小视频在线| 久久中文字幕人妻熟女| √禁漫天堂资源中文www| 一区二区三区激情视频| 国产私拍福利视频在线观看| 99精品久久久久人妻精品| 亚洲国产欧美网| 国产免费av片在线观看野外av| 欧美人与性动交α欧美精品济南到| 黄色a级毛片大全视频| 12—13女人毛片做爰片一| 露出奶头的视频| 一卡2卡三卡四卡精品乱码亚洲| 夜夜看夜夜爽夜夜摸| 成在线人永久免费视频| 最新在线观看一区二区三区| 国产伦一二天堂av在线观看| 男女下面进入的视频免费午夜| 1024手机看黄色片| 黑人欧美特级aaaaaa片| 天天躁夜夜躁狠狠躁躁| or卡值多少钱| 老司机午夜福利在线观看视频| 欧美性长视频在线观看| 国产高清激情床上av| 欧美人与性动交α欧美精品济南到| 不卡av一区二区三区| 亚洲精品久久成人aⅴ小说| 黄频高清免费视频| 亚洲国产中文字幕在线视频| 国产爱豆传媒在线观看 | 中文资源天堂在线| 日韩精品免费视频一区二区三区| 久久天躁狠狠躁夜夜2o2o| 91大片在线观看| 我的老师免费观看完整版| АⅤ资源中文在线天堂| 很黄的视频免费| 亚洲欧美精品综合一区二区三区| 少妇被粗大的猛进出69影院| 一个人观看的视频www高清免费观看 | 亚洲乱码一区二区免费版| 国产伦人伦偷精品视频| 在线播放国产精品三级| 午夜亚洲福利在线播放| tocl精华| 中文在线观看免费www的网站 | 亚洲aⅴ乱码一区二区在线播放 | 色综合亚洲欧美另类图片| 这个男人来自地球电影免费观看| 精品国产乱子伦一区二区三区| 午夜福利在线在线| 国产成人精品无人区| 九色国产91popny在线| 999久久久精品免费观看国产| 一个人观看的视频www高清免费观看 | 91老司机精品| 国产蜜桃级精品一区二区三区| 又爽又黄无遮挡网站| 国内少妇人妻偷人精品xxx网站 | 窝窝影院91人妻| 免费电影在线观看免费观看| 亚洲天堂国产精品一区在线| 香蕉av资源在线| 亚洲男人的天堂狠狠| 午夜福利18| 99热6这里只有精品| 999精品在线视频| 午夜免费激情av| 精品福利观看| 亚洲av熟女| 亚洲一区二区三区不卡视频| 琪琪午夜伦伦电影理论片6080| 日本一本二区三区精品| 亚洲激情在线av| 国产精品 欧美亚洲| 美女午夜性视频免费| 国产亚洲欧美98| 18禁观看日本| 久久久国产成人免费| 日韩中文字幕欧美一区二区| 精品久久久久久久末码| 精品不卡国产一区二区三区| www.www免费av| 国产精品1区2区在线观看.| 亚洲av熟女| 欧美性长视频在线观看| av视频在线观看入口| 无限看片的www在线观看| 国产精品自产拍在线观看55亚洲| 看黄色毛片网站| 亚洲av成人av| 欧美国产日韩亚洲一区| 好男人电影高清在线观看| 精品免费久久久久久久清纯| 校园春色视频在线观看| 免费在线观看视频国产中文字幕亚洲| av天堂在线播放| 我要搜黄色片| 成年免费大片在线观看| 国产精品一区二区三区四区免费观看 | 国产精品香港三级国产av潘金莲| 成年女人毛片免费观看观看9| 国产野战对白在线观看| 可以免费在线观看a视频的电影网站| 国产精品久久久久久亚洲av鲁大| aaaaa片日本免费| 国产高清视频在线播放一区| 亚洲人与动物交配视频| 亚洲精品在线观看二区| 首页视频小说图片口味搜索| 美女大奶头视频| 性欧美人与动物交配| 亚洲乱码一区二区免费版| 夜夜夜夜夜久久久久| 国产精品久久久久久亚洲av鲁大| 国产97色在线日韩免费| 精华霜和精华液先用哪个| 久久久国产成人精品二区| 亚洲成人精品中文字幕电影| 91麻豆av在线| av欧美777| 欧美中文综合在线视频| av欧美777| √禁漫天堂资源中文www| 欧美黑人精品巨大| 国产aⅴ精品一区二区三区波| 91麻豆精品激情在线观看国产| 日韩欧美国产在线观看| 波多野结衣高清作品| 日韩大尺度精品在线看网址| 亚洲国产欧美人成| 欧美乱码精品一区二区三区| 51午夜福利影视在线观看| 超碰成人久久| 一个人观看的视频www高清免费观看 | 中文字幕av在线有码专区| 99在线人妻在线中文字幕| 日日爽夜夜爽网站| 亚洲狠狠婷婷综合久久图片| 免费在线观看视频国产中文字幕亚洲| 亚洲国产精品成人综合色| 国语自产精品视频在线第100页| 国产精品久久久久久人妻精品电影| 亚洲色图av天堂| 中亚洲国语对白在线视频| 一级黄色大片毛片| 一本一本综合久久| 国产三级中文精品| 久久草成人影院| 一a级毛片在线观看| 欧洲精品卡2卡3卡4卡5卡区| 1024香蕉在线观看| 国产99久久九九免费精品| 变态另类丝袜制服| 在线观看舔阴道视频| 精品国产乱码久久久久久男人| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av熟女| 国产精品98久久久久久宅男小说| 久久性视频一级片| 又粗又爽又猛毛片免费看| aaaaa片日本免费| 亚洲欧美日韩无卡精品| 日本三级黄在线观看| 国产精品,欧美在线| 国产成人av激情在线播放| 黄色视频不卡| 成人永久免费在线观看视频| 夜夜爽天天搞| 999久久久精品免费观看国产| 欧美日韩一级在线毛片| 精品人妻1区二区| 狂野欧美激情性xxxx| 午夜a级毛片| 亚洲第一电影网av| 久久精品影院6| 国产成人精品久久二区二区免费| 黄色视频,在线免费观看| 天天一区二区日本电影三级| 这个男人来自地球电影免费观看| 大型av网站在线播放| 欧美激情久久久久久爽电影| 久久久久九九精品影院| 亚洲成人久久性| 非洲黑人性xxxx精品又粗又长| 无遮挡黄片免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国产精品久久久久久精品电影| 正在播放国产对白刺激| 床上黄色一级片| 中文资源天堂在线| 一个人观看的视频www高清免费观看 | 欧美乱色亚洲激情| 中亚洲国语对白在线视频| 国产99白浆流出| 狂野欧美白嫩少妇大欣赏| 久久国产乱子伦精品免费另类| 国产激情久久老熟女| 国产主播在线观看一区二区| 婷婷丁香在线五月| 色在线成人网| АⅤ资源中文在线天堂| 男女视频在线观看网站免费 | 国产免费av片在线观看野外av| 夜夜夜夜夜久久久久| 久久久久亚洲av毛片大全| 午夜福利成人在线免费观看| 色综合婷婷激情| 99在线人妻在线中文字幕| 亚洲第一电影网av| 国产不卡一卡二| 国产精品亚洲av一区麻豆| 久久久国产欧美日韩av| 久久久精品国产亚洲av高清涩受| 91字幕亚洲| 国产成人精品无人区| 91老司机精品| 精品乱码久久久久久99久播| 午夜久久久久精精品| 美女午夜性视频免费| 在线观看免费午夜福利视频| 日韩欧美 国产精品| 午夜精品在线福利| 亚洲国产高清在线一区二区三| 麻豆av在线久日| 曰老女人黄片| 国产高清有码在线观看视频 | 三级国产精品欧美在线观看 | 十八禁人妻一区二区| 亚洲一区中文字幕在线| 狠狠狠狠99中文字幕| 国产伦一二天堂av在线观看| 嫩草影院精品99| 亚洲av成人一区二区三| 哪里可以看免费的av片| 午夜精品一区二区三区免费看| 国产精品免费一区二区三区在线| 久久伊人香网站| avwww免费| 日韩欧美在线二视频| 午夜福利欧美成人| 丰满人妻熟妇乱又伦精品不卡| 国产av不卡久久| 又黄又粗又硬又大视频| 日本一二三区视频观看| av中文乱码字幕在线| 国产精品永久免费网站| 久久香蕉国产精品| 日韩精品中文字幕看吧| 欧美乱妇无乱码| 国产亚洲精品第一综合不卡| 一本精品99久久精品77| 国产精品自产拍在线观看55亚洲| 三级男女做爰猛烈吃奶摸视频| 日韩有码中文字幕| 啦啦啦免费观看视频1| xxx96com| 香蕉丝袜av| 黑人巨大精品欧美一区二区mp4| 狂野欧美白嫩少妇大欣赏| av在线播放免费不卡| 淫妇啪啪啪对白视频| 极品教师在线免费播放| 熟女电影av网| 亚洲人与动物交配视频| 草草在线视频免费看| 精品国产乱子伦一区二区三区| 在线观看美女被高潮喷水网站 | 免费观看人在逋| 中文亚洲av片在线观看爽| 精品一区二区三区四区五区乱码| 婷婷六月久久综合丁香| 亚洲中文日韩欧美视频| 精品久久久久久久末码| 午夜福利在线在线| 一本精品99久久精品77| 亚洲九九香蕉| 在线观看www视频免费| 国产区一区二久久| 亚洲,欧美精品.| 国产午夜精品论理片| 一本久久中文字幕| 久久久国产欧美日韩av| 黑人欧美特级aaaaaa片| 男女之事视频高清在线观看| 三级国产精品欧美在线观看 | www日本在线高清视频| 午夜精品在线福利| 啦啦啦免费观看视频1| 看免费av毛片| 欧美在线一区亚洲| 大型黄色视频在线免费观看| 法律面前人人平等表现在哪些方面| 久久久水蜜桃国产精品网| 一本大道久久a久久精品| 日本一二三区视频观看| 国产精品久久久av美女十八| 免费在线观看日本一区| 99久久无色码亚洲精品果冻| 男人舔奶头视频| 中文字幕最新亚洲高清| 日韩欧美免费精品| 国产一级毛片七仙女欲春2| 国产v大片淫在线免费观看| 亚洲精品美女久久av网站| 亚洲一区中文字幕在线| 99久久国产精品久久久| 他把我摸到了高潮在线观看| 香蕉久久夜色| 亚洲五月婷婷丁香| 久久精品人妻少妇| 女人被狂操c到高潮| 色噜噜av男人的天堂激情| 日韩欧美国产在线观看| 亚洲精品美女久久久久99蜜臀| 欧美日韩福利视频一区二区| 欧美高清成人免费视频www| 久久中文字幕人妻熟女| 黄色成人免费大全| 国产真实乱freesex| 日韩成人在线观看一区二区三区| 国产高清videossex| 久久香蕉国产精品| 免费在线观看视频国产中文字幕亚洲| 亚洲五月天丁香| 成人av一区二区三区在线看| 国产69精品久久久久777片 | 2021天堂中文幕一二区在线观| 老司机福利观看| 国产蜜桃级精品一区二区三区| 99热这里只有是精品50| 亚洲中文av在线| 一本综合久久免费| 色精品久久人妻99蜜桃| 国产爱豆传媒在线观看 | 久久久久亚洲av毛片大全| www国产在线视频色| 久久中文看片网| 午夜日韩欧美国产| 欧美黑人巨大hd| 欧美一区二区精品小视频在线| 精品欧美一区二区三区在线| 99热只有精品国产| 日韩欧美在线乱码| 男女视频在线观看网站免费 | 亚洲人与动物交配视频| av国产免费在线观看| 国产精品一区二区免费欧美| 视频区欧美日本亚洲| 成人国产综合亚洲| 真人一进一出gif抽搐免费| 女人高潮潮喷娇喘18禁视频| 久久伊人香网站| 黄色视频不卡| 757午夜福利合集在线观看| 两性夫妻黄色片| av天堂在线播放| 欧美中文综合在线视频| 亚洲精品国产精品久久久不卡| 欧美乱色亚洲激情| 俺也久久电影网| 日本五十路高清| www.精华液|