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

    Input-to-State Stabilization of Nonlinear Impulsive Delayed Systems: An Observer-Based Control Approach

    2022-07-18 06:17:18YuhanWangXiaodiLiandShijiSongSenior
    IEEE/CAA Journal of Automatica Sinica 2022年7期

    Yuhan Wang,,, Xiaodi Li,,, and Shiji Song, Senior,

    Abstract—This paper addresses the problems of input-to-state stabilization and integral input-to-state stabilization for a class of nonlinear impulsive delayed systems subject to exogenous disturbances. Since the information of plant’s states, time delays, and exogenous disturbances is often hard to be obtained, the key design challenge, which we resolve, is the construction of a state observer-based controller. For this purpose, we firstly propose a corresponding observer which is independent of time delays and exogenous disturbances to reconstruct (or estimate) the plant’s states. And then based on the observations, we establish an observer-based control design for the plant to achieve the inputto-state stability (ISS) and integral-ISS (iISS) properties. With the help of the comparison principle and average impulse interval approach, some sufficient conditions are presented, and moreover, two different linear matrix inequalities (LMIs) based criteria are proposed to design the gain matrices. Finally, two numerical examples and their simulations are given to show the effectiveness of our theoretical results.

    I. INTRODUCTION

    AS an important subclass of hybrid systems, impulsive systems consist of three main components: An ordinary differential equation, a difference equation, and an impulsive law [1], [2]. These systems have drawn considerable attention due to their wide applications in many challenging areas such as networked control systems [3], multi-agent systems [4],optimization problems [5], etc. In addition, time delays are ubiquitous in many fields of science and engineering and might bring about undesired dynamics such as oscillation,instability, and some complicated phenomena, see [6]–[12].Consequently, considering the fact that many dynamical systems are affected by impulsive effects and time delays,stability properties of such systems which are also called impulsive delayed systems have been extensively investigated, see [13]–[20].

    In the past decades, the concept of input-to-state stability(ISS), firstly introduced by Sontag in [21] and then extended to delayed systems by Teel in [22], plays an important role in characterizing the influences of external inputs on nonlinear systems. Roughly speaking, the ISS property implies that a system can be asymptotically stable in input-free case and it also has bounded state under bounded input. As a natural generalization of ISS, integral-ISS (iISS) is a weaker but still very meaningful concept, see [23]. Recently, the problems of ISS and iISS were investigated for (delay-free) impulsive systems in [24]–[26] and were subsequently also studied for impulsive delayed systems in [14], [16], [17], [20]. For example, [24] firstly developed the ISS-Lyapunov function method to impulsive systems with external inputs affecting both the continuous dynamics and discrete dynamics, [14]established Razumikhin-type theorems for ISS and iISS of a class of impulsive delayed systems, and [20] investigated the ISS property of time-delay systems subject to delayed impulse effects. Another interesting topic that has been investigated extensively is the input-to-state stabilization problems, that is,designing control laws for plants such that the closed-loop systems achieve the ISS property with respect to exogenous disturbances, see [27], [28]. In particular, [29], [30] studied the input-to-state stabilization of nonlinear systems from impulsive control point of view. However, there is no result for input-to-state stabilization that takes into account the effects of impulsive perturbations and time delays. It is,therefore, important and necessary to design a control scheme for impulsive delayed systems to achieve ISS from impulsive perturbation point of view.

    It is well known that the full state information of practical systems might not be completely measurable due to the limitation of implementation cost or measurement technology,which implies that the state feedback is generally limited in control engineering. To this end, numerous interesting works have considered the observer-based output feedback control approach, that is, a sort of dynamic output feedback control approach that can estimate plant’s states online, see [31]–[33].Recently, such control approach has been widely studied for impulsive systems. For example, [34] proposed an observerbased fault-tolerant control method for a class of hybrid impulsive systems, [35] studied robust stabilization for a class of linear impulsive systems, and [36] solved the observerbased quasi-synchronization problem of delayed dynamical networks with impulsive effects. However, it should be noted that when the impulsive systems involve unmeasurable time delays and unknown exogenous disturbances, the abovementioned results are infeasible. Since such undesirable factors may complicate the problems of stability analysis,state estimation, and controller design, it is natural to ask whether it is possible to design an observer-based controller for impulsive systems involving unmeasurable time delays to achieve the ISS performance with respect to exogenous disturbances. So far, to the authors’ best knowledge, this problem has not been fully investigated yet, and remains to be challenging, which motivates our present work.

    Statement of Contributions:Motivated by the above discussions and practical requirements, our primary interest in this paper is to design observer-based controllers for input-tostate stabilization of a class of nonlinear impulsive systems with unmeasurable time delays and unknown exogenous disturbances. Based on the comparison principle and average impulse interval approach, some stability criteria are derived and the observer-based output feedback controllers are correspondingly designed. The effectiveness of the developed results is illustrated by numerical simulations. The main contributions of this paper are summarized as follows.

    1) Differently from the classical state observers for timedelay systems in [31], [33], [36], [37], in which the construction of the state observers depends on the full information of time delays, the state observer proposed in this paper is delay-independent and can be applied to the case involving unavailable and unmeasurable time delays.

    2) The observer-based control for ISS of impulsive systems has been extensively studied, such as [30]–[35]. However,those results cannot be applied to the case involving both time delays and exogenous disturbances. While in this paper, some ISS results have been obtained, in which an observer-based output feedback controller which does not need the full information of time delays and exogenous disturbances has been designed.

    3) In order to design the gain matrices of observer-based output feedback controllers conveniently, two different design schemes are proposed in terms of linear matrix inequalities(LMIs), which can be implemented in different cases.

    The remainder of this paper is organized as follows. In Section II, we precisely explain the type of systems considered and provide the relevant notations, assumptions, definitions, and lemmas. Main results are proposed in Section III.Two numerical examples are given in Section IV to show the effectiveness of the proposed design. Finally, Section V concludes the paper.

    II. PRELIMINARIES

    Notations:We use the common definition of class K, K∞,and KL functions from [38]. Let R denote the set of real numbers, R+the setofnonnegativereal numbers,Z+thesetof positive integer numbers, RnandRn×mthen-dimensionalandn×m-dimensional real spaces equipped with the Euclidean norm ||·||, respectively, λmax(A) and λmin(A) the maximum and minimum eigenvalues of symmetric matrix A,respectively. A>0 or A<0 denotes that the matrix A is a symmetricandpositiveornegativedefinitematrix. The notations AT, A?1denotethetransposeandtheinverse ofA,respectively.Iis the identity matrix with appropriate dimensions and ? denotes the symmetric block in one symmetricmatrix.ForanyintervalS1?RandanysetS2?Rk,1≤k≤2n, we putPC(S1,S2)?{υ:S1→S2is continuous everywhere except at finite number of pointst, at which υ(t?), υ(t+) exist and υ(t+)=υ(t)}. Fora,b∈R,a0,wedenotePC([t0?τˉ,t0],Rk)and //υ//τˉ?//υ//[t0?τˉ,t0].

    Considerthe following nonlinearimpulsivedelayedsystem:

    wherex(t)∈Rnis the state,x˙(t) is the upper right-hand derivative ofx(t),u(t)∈Rqis the control input, ω(t)∈Rnis the bounded exogenous disturbance, andy(t)∈Rmis the output (or measurement) withm≤n;A0∈Rn×n,A1∈Rn×n,Ad∈Rn×n,Bc∈Rn×q,Bd∈Rn×n,C∈Rm×n, andD∈Rn×nare constant matrices; τ(t) is the time-varying delay satisfying τ(t)∈[0, τˉ] with τˉ>0;f(x(·))=(f1(x1(·)),...,fn(xn(·)))T∈Rnandg(x(·))=(g1(x1(·)),...,gn(xn(·)))T∈Rnare the nonlinear functions satisfying certain conditions which will be given later. To prevent the occurrence of accumulation points, such as Zeno phenomenon, we only consider that the impulse time sequence {tk}k∈Z+satisfies 0 ≤t0

    Assumption 1:For nonlinear functionsfandg, the structure isknown, andtherearepositiveconstantsandsuch that|fj(s1)?fj(s2)|≤|s1?s2|and |gj(s1)?gj(s2)|≤|s1?s2|for alls1,s2∈R,j=1,...,n. In particular,f(0)=g(0)=0.DenoteLf=diag() andLg=diag() for later use.

    In general, the plant’s states might not be fully available due to the physical limitation or the implementation cost [32],[37], which means that it is very hard to apply a full state feedback control for the plant in such case. Thus, a state observer can be constructed to estimate the plant’s states by utilizing the available measurements. In this paper, the assumption and processing of nonlinear part are not new [39],[40], but we consider the case that the information of time delays and exogenous disturbances cannot be completely obtained. To this end, we firstly consider the following state observer which is independent of time delays and exogenous disturbances:

    whereK∈Rq×nis the control gain matrix to be designed.Under control input (3), the plant can be rewritten as

    wheree(t)?x(t)?(t) is the estimation error between the plant and the observer. It then follows from systems (2) and(4) that the dynamics of the error system is derived as follows:

    whereF(e(·))?f(x(·))?f(x?(·)) andG(e(·))?g(x(·))?g(x?(·)).According to systems (2) and (5), we obtain the augmented closed-loop impulsive system as follows:

    with the following vectors:

    and the following constant matrices:

    For further discussion, we introduce some definitions and lemmas.

    Definition 1 [14]:For given impulse time sequence {tk}k∈Z+,system (1) is said to be stabilized to input-to-state stability(ISS) by controller (3), if there exist functions β ∈KL and γ ∈K∞such that for any ? ∈and any exogenous disturbance ω (t), the solutionx(t) of system (4) satisfies

    Definition 2 [14]:For given impulse time sequence {tk}k∈Z+,system (1) is said to be stabilized to integral-ISS (iISS) by controller (3), if there exist functions β ∈KL and α,γ ∈K∞such that for any ? ∈and any exogenous disturbance ω (t),the solutionx(t) of system (4) satisfies

    Remark 1:The above definitions depend on the choice of the sequence{tk}k∈Z+.However,itisoften ofinterestto characterize ISS(iISS)overclassesof sequences {tk}k∈Z+. To this end, thesystem issaidto bestabilized to uniformly ISS(iISS) over a given class F if the system is ISS (iISS) for every sequence in F with functionsβandγthat are independent of the choice of the sequence, see [24].

    Definition 3 [24]:The average impulse interval (AII) of the impulsetime sequence{tk}k∈Z+is said tobe notless than τ?, if thereexistpositive constantsN0and τ?such that

    whereN(T,t) is the number of impulses of the sequence{tk}k∈Z+occurring on the interval (t,T], andN0is the elasticity number. Denote such kind of impulse time sequences by set F+[τ?,N0]for later use.

    Lemma 1[16]:Letfunctionsh1(t),h2(t)∈PC([t0?,+∞),R+), and φ(t)∈PC([t0,+∞),R+).If thereexist∈R,∈R+, and∈R+such that

    and

    it then holds that for allt≥t0,

    Remark 2:The proof of Lemma 2 is similar to that in [16].However, one can see that the main difference is that we only need thatis a real constant, which can lead to wider applications.

    III. MAIN RESULTS

    The purpose of this paper is to derive some sufficient conditions for input-to-state and integral input-to-state stabilization of system (1) via observer-based control (2) and(3). To do this, we firstly find sufficient conditions for uniformly ISS (iISS) of system (6), which means that the system (1) is stabilized to uniformly ISS (iISS) under the given observer-based control scheme. The proof is based on Lemmas 1 and 2.

    Based on (11) and (12), we derive the following comparison system:

    Assume thatvε(t) is the corresponding maximal solution of system (13) for any given ε >0. Based on Lemma 1, it can be deduced thatV(t)≤v(t)≤vε(t) for allt≥t0. It then follows from Lemma 2 that:

    for allt≥t0. Using Definition 3, we have

    where λ3=λ1?lnμ/τ?>0 and η=μN0. Substituting (13) and(15) into (14) yields that

    From (16) and (18), we observe that

    where Then, the following proof is separated into two parts.

    Part I:In this part, we aim to prove that system (1) is stabilized to uniformly ISS over the class F+[τ?,N0]. Using(19), we get

    for allt≥t0, which implies that

    Remark 3:In particular, if we consider the case thatN0=1,then one can see that the minimal impulse interval is used instead of the average impulse interval, which implies that two consecutive impulse instants must be separated by at leastτ?units of time. Then, we can derive the following corollary from Theorem 1.

    Corollary 1:Given matricesK∈Rq×nandL∈Rn×m, if there exist 2n×2nmatricesP>0,R>0, 2n×2ndiagonal matricesQi>0,i=1,2,3,and positiveconstantsτ?,λ1, λ2,and μ>1,s uchthatτ?>lnμ/(λ1?μλ2)>0,Q3<λ2P, (9),and(10), then by the observer-based controller (3) with (2),system (1) is stabilized to uniformly ISS (iISS) over the class Fτ?,where Fτ?denotesthe class of impulse sequences{tk}k∈Z+satisfyingtk?tk?1≥τ?.

    Next, consider the case thatD=Iwith μ>1, we can formulate Theorem 1 as follows.

    Corollary 2:Suppose thatD= √Iwith μ>1. Given matricesK∈Rq×nandL∈Rn×m, if there exist 2n×2nmatricesP>0,R>0, 2n×2ndiagonal matricesQi>0,i=1,2,3, and positive constantsN0, τ?, λ1, and λ2, such that τ?>lnμ/(λ1?μN0λ2)>0 ,LˇT Q3Lˇ<λ2P, and (9), then by the observer-based controller (3) with (2), system (1) is stabilized to uniformly ISS (iISS) over the class F+[τ?,N0].

    Remark 4:In Theorem 1, Corollaries 1, and 2, the matrix inequality (9) is nonlinear, which limits the design of gain matricesKandLvia the MATLAB LMI Toolbox. To overcome this difficulty, two different cases are fully considered and some necessary matrix transformations need to be applied, which leads to the following results.

    Theorem 2:Suppose that the control input matrixBcis full columnrank.ForgivenpositiveconstantsN0,λ1,λ2,andμ>1,if thereexist apositive constantτ?,n×nmatricesP1>0,P2>0 ,R1>0,R2>0,n×ndiagonal matricesQij>0,i=1,2,3,j=1,2,and constantmatricesX∈Rn×n,Y∈Rn×m, such thatτ?>lnμ/(λ1?μN0λ2)>0 ,Q3<λ2P,(10), and

    where

    with

    which means that

    where

    with

    Then, one can see that

    Next, pre-multiplying and post-multiplying both sides of(21) with, it holds that

    Remark 5:There have been various results on input-to-state stabilization of linear (nonlinear) systems involving unmeasurable states in the literature, see [30], [34], [35], [37],[41]. Similarly to those results, the output feedback controller we designed is also based on a state observer. One can see that the possible time delays are excluded in [30], [34], [35]. In the framework of time delays, [37] studies the observer-based input-to-state stabilization problem for nonlinear delayed system. However, the state observer proposed in [37] depends on the full information of time delays. That is, the time delays need to be fully available and measurable. Although the observer designed in [41] does not require the full information of network-induced delays, the state delays of system itself are excluded. In this paper, we propose an observer-based controller which can stabilize nonlinear impulsive systems involving unmeasurable time delays and unknown exogenous disturbances to ISS (iISS), where the information of time delays and exogenous disturbances is not needed. Moreover,the destabilizing effects of impulsive perturbations are also considered in this paper, which leads to wider applications.

    Remark 6:Theorem 2 presents some LMI-based conditions to design observer-based control input u such that system (1)is stabilized to ISS (iISS). When using Theorems 2 and 3, it is necessary to give proper parameters such that inequalities τ?>lnμ/(λ1?μN0λ2)>0,Q3<λ2P, (10), and (20) hold.Then the feasible solutions of those inequalities, that is,X,Y,P1,P2,R1,R2, andQi j>0,i=1,2,3,j=1,2, can be derived by the MATLAB LMI Toolbox. Based on the above analysis,the gain matrices can be designed by such known matrices and the low bound of the AII constant τ?can be derived. From the algorithm point of view, since positive constants λ1and λ2describe the possible convergence rate of the Lyapunov functionV, it is desirable to give a larger λ1and a smaller λ2when solving LMIs in Theorem 2. In addition, the positive constant μ>1 describes the possible jump amplitude of the Lyapunov functionVat impulse instants. Thus, a smaller μ is preferred for the feasibility of LMI in (10). Similarly, the above analysis is suitable to Theorem 3.

    Remark 7:In Theorem 2, we assume that the control input matrixBcis full column rank, which implies thatBcis invertible. Such assumption is crucial for the design of the controller gainK. While in Theorem 3, we assume that the output matrixCis full row rank, which implies thatCCTis invertible. Such assumption is crucial for the design of the observer gainL. One can see that the proposed results in Theorems 2 and 3 can be applied for different cases. Hence,they are different but complementary to each other. Moreover,how to derive some less conservative conditions for the general case is an interesting topic. This issue will be explored in our future work.

    IV. NUMERICAL EXAMPLES

    In this section, we present two numerical examples and relevant simulations to demonstrate the validity of the above designed observer-based control schemes.

    Example 1:Consider the following nonlinear delayed system:

    with

    where the time delay τ(t) is unmeasurable and satisfies τ(t)∈[0 0.2], and the exogenous disturbance ω(t) is unknown but bounded. When there is no control input (i.e.,u(t)=0), by simulation, Fig. 1 shows that the system (23) is not ISS with the choiceofτ(t)=0.2|sint|,ω(t)=[0.8cost1.1sint]T,initial condition [?1(s) ?2(s)]T=[?1.5 1.4]Tfors∈[?0.2,0].In the following, we consider two different cases.

    Fig. 1. Trajectories of x(t) without control input.

    Case I:The control input matrixBcis full column rank. In this case, we will design the observer-based controller (3) with(2) such that system (23) with the following impulsive perturbations:

    achieves the ISS property based on Theorem 2. We firstly assume thatBc=[0.5 0.3]TandC=[0 7]. Then, for given λ1=5.8, λ2=0.28 ,μ=2.56,andN0=2,bycalculatingthe rangeofvaluesforτ? andsolvingtheLMIs, we obtain τ?∈(0.2371,+∞) and the following gain matrices:

    Thus, according to Theorem 2, system (23) and (24) can be stabilized to uniformly ISS over the class F+[τ?,2] with τ?∈(0.2371,+∞). For simulations, we take impulse time sequence {tk}k∈Z+={t3n?2=n?0.9,t3n?1=n?0.8,t3n=n}n∈Z+∈F+[1/3,2] and still consider τ(t)=0.2|sint|,ω(t)=[0.8cost1.1sint]T, initial condition [?1(s) ?2(s)]T=[?1.5 1.4]Tfors∈[?0.2,0]. Moreover, the observer-based controller (3) with(2) is given with the designed gains in (25). One can see that Fig. 2 shows the control signals, Fig. 3 shows the trajectories of plant states, and Fig. 4 shows the trajectories of observation error signals. The corresponding simulations reflect the validity of our results.

    Fig. 2. Trajectory of u (t) in Case I.

    Fig. 3. Trajectories of x(t) under control input in Case I.

    Fig. 4. Trajectories of e(t) with the initial condition []T =[?2.7 2.8]T in Case I.

    Case II:The output matrixCis full row rank. In this case,to study the input-to-state stabilization of system (23) and(24), we suppose thatBc=[0 5.2]TandC=[3 6]. If we chooseparametersλ1=5, λ2=0.32,μ=2.56,andN0=2,t hen byTheorem 3,wecan obtain τ?∈(0.3238,+∞) and

    Thus, using the above calculated solutions, the ISS property of system (4) is achieved. In simulations, τ(t), ω(t), ?(s), and{tk}k∈Z+aresameasthoseinCase I. The corresponding simulationsare shownin Figs. 5–7.

    Fig. 5. Trajectory of u (t) in Case II.

    Fig. 6. Trajectories of x(t) under control input in Case II.

    Fig. 7. Trajectories of e (t) with the initial condition []T =[?2.7 2.8]T in Case II.

    Example 2:To further demonstrate the flexibility and effectiveness of the proposed design scheme, we consider a chemical reactor system with two reactorsAandB[42], as shown in Fig. 8. Based on mass balances, when there is no impulse, the dynamics of the reactants can be modeled approximately as follows:

    Fig. 8. Chemical reactor system.

    wherex1(·) andx2(·) are the compositions;y(t) is the output;kAandkBare the reaction constants; θAand θBare the reactor residence times;rAandrBare the recycle flow rates;vAandvBare the reactor volumes; ? is the feed rate; the time delayτ(t)is unmeasurable and satisfies τ (t)∈[0 0.5], and the exogenous disturbances ω1(t) and ω2(t) are unknown but bounded. Some other information on system (26) can be found in [42]. For the chemical reactor system (26), we consider thatkA=1,kB=0.5, θA=θB=2 ,rA=0.5,rB=0.25,vA=vB=0.5, and ?=0.5. Then, we obtain the following system:

    with

    When there is no control input (i.e.,u(t)=0), by simulation,Fig. 9 shows that the system (27) is not ISS with the choice of τ(t)=0.5|cost|,ω(t)=[0.2sint0.25cost]T,initial condition[?1(s)?2(s)]T=[0.8?0.3]Tfors∈[?0.5,0].

    In practice, many chemical reactions are affected by discontinuous perturbations at certain instants, which may lead to sudden changes in the states of the chemical reactor system. To describe it, we consider such sudden changes as the impulsive perturbations in the form of (24). Moreover, due to the complexity of the chemical reactor system, the full information of system states is hard to be obtained. Thus, it is more realistic to design the observer-based controller (3) with(2) for the input-to-state stabilization of system (27) with impulses (24). If we choose parameters λ1=5.7 , λ2=1.6,μ=1.69, andN0=2,thenbyTheorem2,wecan obtain τ?∈(0.4643,+∞)andthe following gainmatrices:

    Fig. 9. Trajectories of x(t) without control input.

    Thus, by using above solutions, the observer-based control inputu(t) can be obtained, which can be used to stabilize system (24)–(27) to ISS. For simulations, we take impulse time sequence{tk}k∈Z+={t3n?2=1.5n?1.4,t3n?1=1.5n?1.1,t3n=1.5n}n∈Z+∈F+[0.5,2] andstill consider τ(t)=0.5|cost|,ω(t)=[0.2sint0.25cost]T,initial condition[?1(s)?2(s)]T=[0.8?0.3]Tfors∈[?0.5,0]. The corresponding control inputu(t)can be found in Fig. 10 . With such control input, the simulation results are shown in Figs. 11 and 12. One can see that the designed controller can render the resulting closedloop system as ISS, which illustrates the effectiveness of the observer-based output feedback control approach.

    Fig. 10. Trajectory of u (t).

    Fig. 11. Trajectories of x (t) under control input.

    V. CONCLUSIONS

    Fig. 12. Trajectories of e(t) with the initial condition[]T =[0.5?0.5]T.

    The paper’s key design idea is how to develop a novel control scheme to stabilize a class of impulsive delayed systems to ISS and iISS when the information of the plant’s states cannot be adequately obtained. We have constructed an observer to estimate the plant’s states and designed an observer-based control scheme for input-to-state stabilization.Using the Lyapunov method and comparison principle, some criteria have been presented for the ISS property of the closedloop system and moreover, the corresponding observer-based controller has been designed. Finally, the advantage and effectiveness of our results have been demonstrated by simulation results. Since practical systems are often affected by the random noises [43], another interesting direction is to develop observer-based control for input-to-state stabilization of impulsive stochastic nonlinear systems.

    日本欧美国产在线视频| 日本vs欧美在线观看视频| 久久精品国产亚洲av天美| 曰老女人黄片| 亚洲国产精品一区三区| 国产成人免费无遮挡视频| 国产精品国产av在线观看| 久久久久久人妻| www日本在线高清视频| 国产成人免费观看mmmm| 捣出白浆h1v1| 免费人成在线观看视频色| 午夜91福利影院| 久久人人爽av亚洲精品天堂| xxx大片免费视频| 亚洲美女黄色视频免费看| 国产成人精品无人区| 亚洲中文av在线| 欧美变态另类bdsm刘玥| 亚洲精品国产av蜜桃| 国内精品宾馆在线| 久久午夜福利片| 看免费av毛片| 国产成人精品婷婷| 中文欧美无线码| 日韩人妻精品一区2区三区| 中文字幕另类日韩欧美亚洲嫩草| www.av在线官网国产| 欧美老熟妇乱子伦牲交| 两性夫妻黄色片 | 国产成人91sexporn| 亚洲一区二区三区欧美精品| 免费av不卡在线播放| 十八禁网站网址无遮挡| kizo精华| 午夜视频国产福利| av在线老鸭窝| 精品亚洲乱码少妇综合久久| 欧美成人午夜免费资源| 人妻少妇偷人精品九色| 国产精品一区二区在线不卡| 精品人妻偷拍中文字幕| 制服丝袜香蕉在线| 涩涩av久久男人的天堂| 18禁动态无遮挡网站| 99热6这里只有精品| 国产精品女同一区二区软件| 国产一区二区在线观看av| 成人亚洲欧美一区二区av| 天堂俺去俺来也www色官网| 精品午夜福利在线看| 国产黄频视频在线观看| av有码第一页| 欧美激情国产日韩精品一区| 国产成人午夜福利电影在线观看| 高清在线视频一区二区三区| 99热这里只有是精品在线观看| 好男人视频免费观看在线| 中文字幕av电影在线播放| 看免费av毛片| 免费av中文字幕在线| 99热国产这里只有精品6| 伊人亚洲综合成人网| 少妇高潮的动态图| 欧美97在线视频| 建设人人有责人人尽责人人享有的| 久久久久久久精品精品| 大片电影免费在线观看免费| 这个男人来自地球电影免费观看 | 中国美白少妇内射xxxbb| 狠狠精品人妻久久久久久综合| 人妻系列 视频| 精品国产国语对白av| 插逼视频在线观看| 日韩熟女老妇一区二区性免费视频| www.av在线官网国产| 久久国产精品大桥未久av| 久久久国产精品麻豆| 久久女婷五月综合色啪小说| 大陆偷拍与自拍| 免费av中文字幕在线| 国产精品人妻久久久影院| 亚洲精品国产色婷婷电影| 亚洲av.av天堂| 久久精品国产鲁丝片午夜精品| 国产成人午夜福利电影在线观看| 一二三四在线观看免费中文在 | 午夜福利,免费看| 免费看av在线观看网站| 女的被弄到高潮叫床怎么办| 18禁观看日本| 国产精品熟女久久久久浪| 中文精品一卡2卡3卡4更新| 日韩不卡一区二区三区视频在线| 久久99一区二区三区| 高清视频免费观看一区二区| 捣出白浆h1v1| 大陆偷拍与自拍| 性高湖久久久久久久久免费观看| 国产成人一区二区在线| 美女国产高潮福利片在线看| 最黄视频免费看| 九九爱精品视频在线观看| av国产精品久久久久影院| 日韩一本色道免费dvd| 精品国产露脸久久av麻豆| 侵犯人妻中文字幕一二三四区| 男人操女人黄网站| 国产一区二区在线观看日韩| 好男人视频免费观看在线| 老熟女久久久| 男的添女的下面高潮视频| 九色成人免费人妻av| 欧美日韩av久久| 日韩在线高清观看一区二区三区| 最黄视频免费看| 久久 成人 亚洲| 九九爱精品视频在线观看| 99国产综合亚洲精品| 在线观看三级黄色| 建设人人有责人人尽责人人享有的| 91精品三级在线观看| 香蕉国产在线看| 在线观看美女被高潮喷水网站| 国产亚洲一区二区精品| 久久精品国产亚洲av涩爱| 少妇被粗大猛烈的视频| 日韩免费高清中文字幕av| 午夜影院在线不卡| 中文字幕人妻丝袜制服| 多毛熟女@视频| 亚洲欧美一区二区三区国产| 久久狼人影院| av在线观看视频网站免费| 久久99热这里只频精品6学生| 久久久久久人妻| 国产精品久久久久久久电影| 亚洲国产欧美在线一区| 性色av一级| 免费看光身美女| 少妇人妻 视频| 国语对白做爰xxxⅹ性视频网站| 夜夜爽夜夜爽视频| 黄片无遮挡物在线观看| 一区二区三区乱码不卡18| 国产精品麻豆人妻色哟哟久久| 最近最新中文字幕大全免费视频 | 一本久久精品| 99re6热这里在线精品视频| 国产免费福利视频在线观看| 久久久国产一区二区| av国产精品久久久久影院| 日本欧美国产在线视频| 一级片免费观看大全| 大片电影免费在线观看免费| 在线观看免费高清a一片| 亚洲欧美中文字幕日韩二区| 日本91视频免费播放| 国产不卡av网站在线观看| 欧美精品一区二区大全| 亚洲美女视频黄频| 日本猛色少妇xxxxx猛交久久| 国产免费视频播放在线视频| 少妇人妻精品综合一区二区| 丝袜喷水一区| 最近中文字幕2019免费版| 亚洲av.av天堂| 婷婷色综合www| 亚洲av综合色区一区| 国产成人午夜福利电影在线观看| 国产精品不卡视频一区二区| 成人毛片a级毛片在线播放| 国产麻豆69| xxxhd国产人妻xxx| 97精品久久久久久久久久精品| 国产成人欧美| a级片在线免费高清观看视频| 亚洲人成网站在线观看播放| 亚洲精品av麻豆狂野| 久久国产亚洲av麻豆专区| 中文字幕人妻熟女乱码| 国产熟女午夜一区二区三区| 一区二区三区精品91| 看免费成人av毛片| 久久99一区二区三区| 国产熟女午夜一区二区三区| 久久人人爽av亚洲精品天堂| 九色亚洲精品在线播放| 国产高清三级在线| 大香蕉久久网| 香蕉国产在线看| 91精品国产国语对白视频| 亚洲精品国产av蜜桃| 九九在线视频观看精品| 亚洲欧美成人综合另类久久久| 国产高清不卡午夜福利| 制服丝袜香蕉在线| 又大又黄又爽视频免费| 多毛熟女@视频| 水蜜桃什么品种好| 天天影视国产精品| 岛国毛片在线播放| 午夜日本视频在线| 亚洲av男天堂| 精品卡一卡二卡四卡免费| 三上悠亚av全集在线观看| 国产深夜福利视频在线观看| 国产色爽女视频免费观看| 久久久久久久大尺度免费视频| 精品国产一区二区三区四区第35| 久久精品久久久久久久性| 老司机影院成人| 丝袜美足系列| 黄色视频在线播放观看不卡| 国产精品一区www在线观看| 日韩欧美精品免费久久| 国产黄色视频一区二区在线观看| 如何舔出高潮| 乱人伦中国视频| 日日啪夜夜爽| 一个人免费看片子| av线在线观看网站| 日韩在线高清观看一区二区三区| 日韩,欧美,国产一区二区三区| 亚洲av欧美aⅴ国产| 亚洲一区二区三区欧美精品| 91国产中文字幕| 亚洲av福利一区| 国产成人av激情在线播放| 高清欧美精品videossex| 欧美 日韩 精品 国产| 精品视频人人做人人爽| 九九在线视频观看精品| av免费观看日本| 国精品久久久久久国模美| av女优亚洲男人天堂| 亚洲国产毛片av蜜桃av| 在线精品无人区一区二区三| 国产精品偷伦视频观看了| 午夜福利视频精品| 亚洲,欧美精品.| 久久这里只有精品19| a 毛片基地| 亚洲精品美女久久久久99蜜臀 | 亚洲四区av| 人人妻人人爽人人添夜夜欢视频| 午夜激情久久久久久久| 亚洲精品一区蜜桃| 看非洲黑人一级黄片| 狂野欧美激情性bbbbbb| 日本黄大片高清| 久久久久久久精品精品| 中国美白少妇内射xxxbb| 国产一区亚洲一区在线观看| 精品99又大又爽又粗少妇毛片| 在线天堂最新版资源| 侵犯人妻中文字幕一二三四区| 人体艺术视频欧美日本| 国产成人a∨麻豆精品| 成年女人在线观看亚洲视频| 亚洲国产精品专区欧美| 内地一区二区视频在线| 香蕉国产在线看| 亚洲欧美中文字幕日韩二区| 欧美xxⅹ黑人| 人妻人人澡人人爽人人| 黄色配什么色好看| 欧美性感艳星| 精品一品国产午夜福利视频| 中国三级夫妇交换| 欧美成人午夜免费资源| 国产精品一国产av| 久久99热6这里只有精品| 满18在线观看网站| 高清不卡的av网站| 老女人水多毛片| 免费观看在线日韩| 欧美另类一区| 亚洲精品中文字幕在线视频| 国产高清不卡午夜福利| 免费久久久久久久精品成人欧美视频 | 国产永久视频网站| 青春草亚洲视频在线观看| 久久人人爽人人片av| 成人综合一区亚洲| 这个男人来自地球电影免费观看 | 最近中文字幕2019免费版| 日日撸夜夜添| 色婷婷久久久亚洲欧美| 免费看不卡的av| 2022亚洲国产成人精品| 午夜福利视频精品| 男人添女人高潮全过程视频| 国产熟女欧美一区二区| 一级,二级,三级黄色视频| 日本黄大片高清| av在线播放精品| 亚洲 欧美一区二区三区| 久久青草综合色| 人成视频在线观看免费观看| 精品视频人人做人人爽| 国产免费又黄又爽又色| 女人久久www免费人成看片| 中文字幕人妻熟女乱码| 国产免费一区二区三区四区乱码| 免费黄网站久久成人精品| 日本免费在线观看一区| 精品久久蜜臀av无| 女性被躁到高潮视频| 久久毛片免费看一区二区三区| 国产黄色免费在线视频| 久久国产亚洲av麻豆专区| av黄色大香蕉| 婷婷色麻豆天堂久久| 久久人人爽av亚洲精品天堂| 九九爱精品视频在线观看| 国产成人精品久久久久久| 亚洲欧洲日产国产| 久久99热6这里只有精品| 亚洲成色77777| 免费不卡的大黄色大毛片视频在线观看| 欧美日韩国产mv在线观看视频| 成年女人在线观看亚洲视频| 亚洲熟女精品中文字幕| 日韩欧美精品免费久久| 18+在线观看网站| 免费人妻精品一区二区三区视频| 黄色视频在线播放观看不卡| 亚洲欧美精品自产自拍| 国产精品熟女久久久久浪| 国产免费现黄频在线看| 黄片播放在线免费| 亚洲人成网站在线观看播放| 黑人猛操日本美女一级片| 成人综合一区亚洲| 国产亚洲欧美精品永久| 国产又爽黄色视频| 色吧在线观看| 中文字幕亚洲精品专区| 免费在线观看黄色视频的| 青春草国产在线视频| 男的添女的下面高潮视频| 欧美 亚洲 国产 日韩一| 亚洲经典国产精华液单| 精品卡一卡二卡四卡免费| 九草在线视频观看| 亚洲情色 制服丝袜| 成人午夜精彩视频在线观看| 国产国拍精品亚洲av在线观看| 男女午夜视频在线观看 | 国产精品久久久久久精品电影小说| 精品国产一区二区三区久久久樱花| 久久韩国三级中文字幕| 亚洲av电影在线进入| 黄色怎么调成土黄色| 亚洲精品一二三| 一边亲一边摸免费视频| 人妻 亚洲 视频| 又粗又硬又长又爽又黄的视频| 精品亚洲成a人片在线观看| 国产国语露脸激情在线看| 美女视频免费永久观看网站| 制服诱惑二区| 国产欧美日韩一区二区三区在线| 欧美日韩精品成人综合77777| 999精品在线视频| 国产黄色免费在线视频| 亚洲精品成人av观看孕妇| 人成视频在线观看免费观看| 狂野欧美激情性bbbbbb| 边亲边吃奶的免费视频| 97人妻天天添夜夜摸| 久久人妻熟女aⅴ| 亚洲综合精品二区| 亚洲av日韩在线播放| 免费日韩欧美在线观看| 满18在线观看网站| a 毛片基地| av网站免费在线观看视频| 一区二区三区乱码不卡18| 欧美日韩综合久久久久久| 欧美日韩国产mv在线观看视频| 国产男女超爽视频在线观看| 欧美日韩国产mv在线观看视频| 免费高清在线观看日韩| 一区二区三区乱码不卡18| 九九在线视频观看精品| 最近手机中文字幕大全| 制服丝袜香蕉在线| 国产精品99久久99久久久不卡 | 国产高清不卡午夜福利| 国语对白做爰xxxⅹ性视频网站| 欧美日韩综合久久久久久| 一级毛片我不卡| 色网站视频免费| 一级毛片我不卡| 插逼视频在线观看| 性色avwww在线观看| 香蕉丝袜av| 久久青草综合色| 九草在线视频观看| 国产成人av激情在线播放| 久久久久久久久久成人| 亚洲综合色网址| 亚洲五月色婷婷综合| 波野结衣二区三区在线| 黑人欧美特级aaaaaa片| 亚洲三级黄色毛片| 精品久久久久久电影网| 晚上一个人看的免费电影| 男的添女的下面高潮视频| 久久午夜福利片| 久久精品熟女亚洲av麻豆精品| av有码第一页| 国产激情久久老熟女| av线在线观看网站| 欧美日韩视频精品一区| 三级国产精品片| 久久国产精品大桥未久av| 亚洲欧洲精品一区二区精品久久久 | 捣出白浆h1v1| 亚洲婷婷狠狠爱综合网| 亚洲av日韩在线播放| 九色亚洲精品在线播放| 日本av免费视频播放| 亚洲精品456在线播放app| 街头女战士在线观看网站| 久久久久久久大尺度免费视频| 久久久久精品久久久久真实原创| 日日啪夜夜爽| 日韩av免费高清视频| 亚洲国产看品久久| 精品一区在线观看国产| 看免费成人av毛片| av卡一久久| 成人免费观看视频高清| 欧美人与性动交α欧美软件 | 丝袜在线中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 欧美精品高潮呻吟av久久| 午夜久久久在线观看| 国产精品不卡视频一区二区| 国产成人免费无遮挡视频| 老司机影院毛片| av电影中文网址| 边亲边吃奶的免费视频| 欧美3d第一页| 日本vs欧美在线观看视频| 国产色婷婷99| 美女xxoo啪啪120秒动态图| 精品酒店卫生间| 久久97久久精品| 91精品三级在线观看| 制服人妻中文乱码| 亚洲三级黄色毛片| 日韩中文字幕视频在线看片| 晚上一个人看的免费电影| av播播在线观看一区| 欧美xxxx性猛交bbbb| 人人妻人人澡人人爽人人夜夜| 高清av免费在线| 国产日韩欧美亚洲二区| 日产精品乱码卡一卡2卡三| 美国免费a级毛片| 国产探花极品一区二区| 亚洲成人一二三区av| 少妇被粗大的猛进出69影院 | 国内精品宾馆在线| 香蕉丝袜av| av电影中文网址| 色网站视频免费| 亚洲国产精品成人久久小说| 久久99精品国语久久久| 国产深夜福利视频在线观看| av又黄又爽大尺度在线免费看| 国产又色又爽无遮挡免| 91在线精品国自产拍蜜月| 99视频精品全部免费 在线| 色婷婷久久久亚洲欧美| 两性夫妻黄色片 | 美女主播在线视频| 十八禁高潮呻吟视频| 亚洲综合色网址| 亚洲精品一二三| 国精品久久久久久国模美| 肉色欧美久久久久久久蜜桃| 一二三四在线观看免费中文在 | 99热全是精品| 精品一区二区三区四区五区乱码 | 久久午夜综合久久蜜桃| 国产精品.久久久| videossex国产| 大话2 男鬼变身卡| 热re99久久国产66热| 国产有黄有色有爽视频| 亚洲精品成人av观看孕妇| 国产一区二区在线观看日韩| 亚洲国产毛片av蜜桃av| 丰满迷人的少妇在线观看| videosex国产| 亚洲欧美中文字幕日韩二区| av免费观看日本| 精品亚洲成国产av| 狂野欧美激情性xxxx在线观看| 国产一区二区激情短视频 | 国产日韩欧美在线精品| 两个人看的免费小视频| 成人免费观看视频高清| 美女国产高潮福利片在线看| 两性夫妻黄色片 | 日韩伦理黄色片| 久久久国产欧美日韩av| 中文字幕最新亚洲高清| 亚洲av电影在线进入| 9热在线视频观看99| 街头女战士在线观看网站| 老熟女久久久| 黑人猛操日本美女一级片| 九草在线视频观看| 国产精品嫩草影院av在线观看| 视频在线观看一区二区三区| 韩国精品一区二区三区 | 亚洲,欧美,日韩| 91精品伊人久久大香线蕉| 亚洲国产精品一区二区三区在线| 热re99久久精品国产66热6| 国产极品粉嫩免费观看在线| 男女无遮挡免费网站观看| 久久ye,这里只有精品| 国产探花极品一区二区| 亚洲四区av| 亚洲美女黄色视频免费看| 91成人精品电影| 亚洲av电影在线观看一区二区三区| 波野结衣二区三区在线| 久久久久国产精品人妻一区二区| 国产极品粉嫩免费观看在线| 亚洲丝袜综合中文字幕| xxxhd国产人妻xxx| 满18在线观看网站| 晚上一个人看的免费电影| 国产男人的电影天堂91| 99re6热这里在线精品视频| 又黄又粗又硬又大视频| 亚洲一区二区三区欧美精品| 内地一区二区视频在线| 丰满乱子伦码专区| 国产精品.久久久| 国产 精品1| 亚洲欧美成人综合另类久久久| 青春草视频在线免费观看| 黄片播放在线免费| 欧美xxⅹ黑人| 国产精品久久久久久久久免| 亚洲精品乱久久久久久| 免费在线观看完整版高清| 国产在视频线精品| 中文欧美无线码| 国产一区有黄有色的免费视频| 久久久精品94久久精品| 成人免费观看视频高清| 欧美xxxx性猛交bbbb| 亚洲国产av影院在线观看| 精品一品国产午夜福利视频| 久久 成人 亚洲| 天天躁夜夜躁狠狠久久av| 晚上一个人看的免费电影| 日韩av免费高清视频| 伦理电影大哥的女人| 狂野欧美激情性bbbbbb| 亚洲伊人久久精品综合| 久久毛片免费看一区二区三区| 国产免费一区二区三区四区乱码| 国产永久视频网站| 亚洲av在线观看美女高潮| 精品久久国产蜜桃| 国产成人精品福利久久| 99久久精品国产国产毛片| 欧美日韩视频精品一区| 亚洲三级黄色毛片| 久久久久人妻精品一区果冻| 国产成人精品福利久久| 老司机亚洲免费影院| 成人免费观看视频高清| 亚洲欧美清纯卡通| 你懂的网址亚洲精品在线观看| 麻豆乱淫一区二区| 国产极品粉嫩免费观看在线| av在线老鸭窝| 亚洲精品国产色婷婷电影| 九草在线视频观看| 在线免费观看不下载黄p国产| 欧美精品一区二区大全| 亚洲av在线观看美女高潮| 国产精品女同一区二区软件| 在线看a的网站| 久久久久久久大尺度免费视频| 免费av中文字幕在线| 精品国产露脸久久av麻豆| 三级国产精品片| 女人精品久久久久毛片| 精品人妻熟女毛片av久久网站| 国产片特级美女逼逼视频| 最新中文字幕久久久久| 一区二区日韩欧美中文字幕 | 免费大片黄手机在线观看| a 毛片基地| av免费在线看不卡| 欧美+日韩+精品| 精品一区二区三区视频在线| 极品人妻少妇av视频| 久久av网站| www日本在线高清视频| 欧美日韩精品成人综合77777| 欧美变态另类bdsm刘玥|