999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

離散切換系統觀測器存在性討論及降維觀測器設計

2018-01-08 02:59:56朱芳來蔡明郭勝輝
自動化學報 2017年12期
關鍵詞:系統設計

朱芳來 蔡明 郭勝輝,2

離散切換系統觀測器存在性討論及降維觀測器設計

朱芳來1蔡明1郭勝輝1,2

對具有未知輸入的離散切換系統討論了未知輸入觀測器(Unknown input observer,UIO)設計方法.首先,對一般離散系統的未知輸入觀測器匹配條件的Lyapunov-type表示,進行了等價性論證;然后,基于不具有未知輸入的離散切換系統的穩定性理論,對具有未知輸入的離散切換系統提出了一種切換降維觀測器設計方法.通過矩陣分塊確定出的觀測器增益矩陣,使得降維觀測器能直接消去未知輸入的影響;然后,在此基礎上提出了一種未知輸入代數重構方法;最后,通過仿真驗證了方法的有效性和正確性.

離散系統,未知輸入估計,觀測器匹配條件,切換系統

現代控制理論中,未知輸入觀測器(Unknown input observer,UIO)設計具有重要意義,自上世紀60年代提出以來,一直受到學者們的廣泛的關注和深入研究[1?2].對于連續受控系統,文獻[3]研究了未知輸入觀測器匹配條件,并將該條件等價于一個Lyapunov-type矩陣代數方程.之后,在這一框架下,未知輸入觀測器設計得到了廣泛的研究[4?10].如文獻[7]在匹配條件滿足的前提下,通過系統擴展并結合奇異系統設計手段,提出了一種滑模觀測器設計方法,以此實現了故障和輸出噪聲的估計.由于未知輸入觀測器匹配條件較為苛刻,很多實際的控制系統并不滿足,因而如何突破該匹配條件的限制,成為近期UIO設計研究熱點之一[11?14],如文獻[12]基于構造輔助輸出的方式,提出了一種解決方案.文獻[13]同樣通過構造輔助輸出,提出了降維觀測器結合高階滑模觀測器的觀測器匹配條件突破方法,并考慮了未知輸入重構問題.文獻[14]將未知輸入分解為代數抑制未知輸入和模型未知輸入,以此來達到滿足觀測器匹配條件的目的.然而,值得強調的是,以上這些成果都是基于連續系統討論的.隨著數字系統廣泛的應用,離散系統下的相關研究受到了關注[15?16].基于這樣的關注,本文首先對離散系統未知輸入觀測器匹配條件進行了討論,并將其轉換為一個Lyapunov-type代數方程組的等價形式.

切換系統通常含有若干個子系統,同時具有一個規定子系統之間切換方式的切換律.切換系統可以描述一大類實際的系統,如脈沖寬度調制(Pulse width modulation,PWM)電路、直流斬波電路、網絡控制系統和飛行器姿態等[17?25].早期針對切換系統的研究主要是考慮其穩定性,如文獻[17?18]針對離散切換系統,討論了平均駐留時間(Average dwell time,ADT)下的切換系統穩定問題,而文獻[19]對切換系統的穩定性研究成果進行了更為深入的總結.針對切換系統的狀態估計的研究相對較少,如文獻[20]針對離散切換系統,設計了異步濾波器來處理一類具有電力電子現實背景的問題,并討論了其穩定性.文獻[21]同時針對連續系統和離散系統設計了延遲觀測器,同時達到了異步切換的效果.對于具有未知輸入的切換系統,文獻[22]針對廣義離散切換系統,設計了未知輸入觀測器來估計系統狀態,并基于線性矩陣不等式證明了系統的穩定性.文獻[23]則在未知輸入匹配條件滿足的前提下,基于系統強可檢測性的假設,針對連續切換系統,設計了未知輸入切換觀測器.文獻[24]在駐留時間假設下研究了切換系統的指數鎮定問題,并消除了切換超調的不利影響.總之,目前針對切換系統的觀測器設計的成果不多,降維未知輸入切換觀測器的成果還鮮有報道.

結合上述背景,本文首先對離散系統未知輸入觀測器的存在性進行討論,然后將其應用于離散切換系統,設計了降維未知輸入切換觀測器,并提出了未知輸入代數重構方法.本文的貢獻在于:1)將連續系統未知輸入觀測器的秩條件和Lyapunov-type代數方程組兩者的等價性結論,推廣到了離散系統;2)基于等價性條件,給出了切換系統未知輸入觀測器存在性的前提條件;3)提出了一種切換降維觀測器設計方法,該方法可以直接消除未知輸入的影響,并提出了一種未知輸入代數重構方法.論文的結構如下:第1節對問題進行了描述,對一般系統,討論了觀測器匹配條件與一個Lyapunov-type代數方程的等價性;第2節提出了降維未知輸入切換系統觀測器設計方法;第3節提出了未知輸入代數重構;第4節給出了仿真分析;第5節對全文進行了總結.

1 問題描述

考慮一類具有未知輸入的離散線性系統

其中,x(k)∈Rn,u(k)∈Rm是已知控制輸入,y(k)∈Rp是可測輸出,η(k)∈Rq是未知輸入.系數矩陣 A ∈ Rn×n,B ∈ Rn×m,C ∈ Rp×n和D ∈Rn×p為已知常數矩陣.不失一般性,假設(A,C)可觀測,B,D為列滿秩,C為行滿秩,并有n≥p≥q.

定義1.定義為系統(1)的Rosen-brock矩陣.當復數z0滿足秩條件rank(R(z0))<n+q時,稱z0為系統(1)的不變零點.

定義2.若系統(1)的所有不變零點均在 z平面單位圓內,則稱系統(1)為最小相系統.其等價描述為:對任何滿足|z|≥1的z∈C,有rank(R(z))=n+q.

假設1.系統(1)是最小相系統.

假設2.rank(D)=rank(CD)=q.

注1.無論是連續系統還是離散系統,假設1和2都是未知輸入觀測器設計的基本前提條件.針對連續系統,文獻[3]給出了一個與假設1和2等價的Lyapunov-type描述,而這樣的等價描述,為未知輸入觀測器設計帶來了極大的便利.接下來我們在引理1中將該等價描述推廣到離散系統,并給出嚴格的證明和算法.

引理1.假設1和2成立的充分必要條件是:對于任意對稱正定矩陣Q∈Rn×n,如下Lyapunov矩陣代數方程組

關于矩陣L∈Rn×p、G∈Rq×p和對稱正定矩陣P∈Rn×n有解.

證明.詳見附錄.

2 切換降維觀測器設計

本節將上述結論應用于切換系統觀測器設計中,結合切換系統的穩定性理論,提出一種切換降維觀測器設計方法.為此,考慮如下具有未知輸入的離散切換系統

其中,x(k)∈Rn,u(k)∈Rm是已知控制輸入,y(k)∈Rp是可測輸出,η(k)∈Rq是未知輸入. 系數矩陣Aσ(k)∈Rn×n,Bσ(k)∈Rn×m和Dσ(k)∈Rn×q為常數切換矩陣,C ∈Rp×n為已知常數矩陣.不失一般性,假設.切換信號 σ(k):N+→ ?={1,2,···,N},其中N 為子系統數.設采樣時刻 k1,k2,···,kv,kv+1,···為切換系統(4)的切換時刻,且不妨假設當k∈[kv,kv+1),第i個子系統被激活,即假設當k∈[kv,kv+1)時,σ(k)=i.

假設3.系統(4)的任意子系統均為最小相系統,而且存在常數0<?<1,對任何滿足|z|≥?的z∈C,有其中

假設 4.

定義3.對任意0≤kl≤ks,以Nσ(k)(kl,ks)表示區間[kl,ks)上的切換次數.若存在常數τa>0使得

成立,則稱常數τa為切換信號σ(k)的平均駐留時間(Average dwell time,ADT),稱N0為抖振界.

引理2[17].考慮離散切換系統

若存在常數0<α<1,β>1,并存在C1函數Vσ(k):Rn→ R,σ(k)∈ ?,和兩類K∞函數 ε1和ε2使得 ?σ(k)=i∈ ?,均有:

且 ?(σ(kv)=i,σ(kv?1)=j)∈ ?×?,i/=j,

則平均駐留時間(ADT)滿足

的切換系統(5)是全局漸近穩定的.

則對應于k∈[kv,kv+1)上的第i個子系統有

引理3.假設3和4的Lyapunov-type等價描述為:存在0<γ1<1,使得對任意的對稱正定矩陣Qi∈Rn×n,如下的矩陣方程組

關于對稱正定矩陣Pi∈Rn×n和矩陣Li∈Rn×p和 Gi∈ Rq×p有解.

證明.首先,在假設3 下,由矩陣C?=(1/?)C和所確定的未知輸入切換系統的任何第i個系統均是最小相系統.事實上,?|z|≥1,由于|z|≥1且0<?<1?|?z|=?|z|≥?,故在假設3下有:

關于對稱正定矩陣Pi∈Rn×n和矩陣Li∈Rn×p和Gi∈Rq×p有解. 如果取γ1=1? ?2,則有0<γ1<1,這時式(8)即為式(7"). □

其中,Ai,11∈Rp×p,Bi,11∈Rp×m,Di,11∈Rp×q,Pi,11∈ Rp×p.

對系統(4)進行線性變換θ(k)=Tσ(k)x(k),

則當系統處于第i個子系統時,一方面,由式(7)的第二式知

注意到KiDi,1+Di,2=0,因而有:

由于切換系統在切換時刻,新子系統的狀態初值為前一采樣時刻的狀態值,即有x(kv)=x(kv?1),所以有:

于是有:在假設3和4下,由引理3知,存在常數0<γ1<1,使得(7)關于對稱正定矩陣Pi∈Rn×n和矩陣Li∈Rn×p和Gi∈Rq×p有解.如果對對稱正定矩陣Pi解作進一步的假設,則提出如下切換降維觀測器設計方法:

定理2.在假設3和4下,如果還存在常數γ2>1使得對于?(σ(k)=i,σ(k?1)=j)∈?×?,i/=j,滿足:

且假設系統(4)切換信號σ(k)的平均駐留時間(ADT)滿足:

則系統

是系統(4)的切換降維未知輸入觀測器使得limk→∞(x(k)??x(k))=0.

證明.定義誤差方程為選取Lyapunov函數為

當系統處于第i個子系統的時刻時,即當k∈[kv,kv+1),由式(14)可得:

且λ1=mini∈‘{λmin(Pi,22)},λ2=maxi∈‘{λmax(Pi,22)},λmin(Pi,22) 表示Pi,22最小的特征值,λmax(Pi,22) 表示Pi,22(i=1,2,···,N)最大的特征值.

由式(7)第一式,利用Schur引理,取Wi=PiLi,可得:

分解上式,有:

再經初等行、列變換,可得:

故有:

由Schur補定理,可得:

也即

對于第i個子系統,式(14)即為

從式(9)和(13)中的第一式,可得:

則有:

由式(16)和(18)可得:

用式(10)減去式(13)中第二式,我們有:

對于?(σ(kv)=i,σ(kv?1)=j)∈?×?,i/=j,從式(11)可得:

也即

因此有

進一步地

也即

注意到式(15)、(18)、(20)和式(12),根據引理2推知觀測器誤差方程(17)是漸近穩定的,這表明式(13)是式(4)的切換降維觀測器. □

3 未知輸入代數重構

本節對未知輸入提出一種代數延遲重構方法,式(4)兩邊同時左乘矩陣C,有:

其中

定理3.在假設3和4的前提下,

是未知輸入η(k)的漸近收斂延遲重構.

證明.由式(21)和(22)可得,未知輸入重構誤差

4 仿真

考慮離散切換系統(4)的系數矩陣為

假設輸入為u(k)=5sin(2k),未知輸入為η(k)=sin(50k),γ1=0.6,γ2=3.8.可以算出:

以及τ?a=1.457.切換序列如圖1所示,圖2給出狀態估計和未知輸入重構,由圖2可以看出,所提方法是正確有效的.

圖1 切換序列Fig.1 Switched sequence

圖2 系統狀態估計和未知輸入重構Fig.2 Estimated states and unknown input reconstruction

5 總結

本文對具有未知輸入的離散切換系統,提出了一種能直接消去未知輸入影響的切換降維觀測器設計方法.為了達到設計目的,事先將文獻[3]中有關連續系統觀測器匹配條件的等價性條件的相關結論,推廣到一般的離散系統,針對離散系統給出了類似于連續系統中觀測器匹配條件的Lyapunov-type等價表示形式,并給出了嚴格的推導過程.然后,將所得結論應用于離散切換系統,提出了降維切換觀測器設計方法,并以此討論了切換觀測器存在的前提條件.在同樣的前提下,提出了一種未知輸入代數重構方法.對同時具有未知輸入和測量噪聲的切換系統,設計切換觀測器并討論存在性,是值得進一步研討的議題.

附錄

為了證明引理1的結論,首先引入下面引理.

引理A1[3].假設2成立的充分必要條件是:存在可逆矩陣T∈Rn×n,S∈Rp×p,使得:

引理A2.對于任意對稱正定矩陣Q∈Rn×n,矩陣代數方程組(3)關于L∈Rn×p,G∈Rq×p和P∈Rn×n有解的充分必要是:=TTPT,=T?1LS,=GS和=TTQT是如下的代數矩陣方程

的解.

證明.直接驗證即可得知.引理1的證明:(必要性:假設1和2成立?式(3)有解).定義

可知

其中,?>0是任意的正常數,為任意矩陣,顯然,是對稱正定矩陣.在假設1和2的前提下,由引理4可知,可檢測,因此,存在矩陣使得如下Lyapunov代數方程

充分性:式(3)有解?假設1和2成立).由于P∈Rn×n為對稱正定矩陣,故矩陣D和DTPD有相同的零空間,即rank?DTPD¢=rank(D),由DTP=GC,得DTPD=GCD,即

因此有rank(D)≤rank(CD)≤rank(D),即有rank(CD)=rank(D),也即假設2成立.

在假設2成立的前提下,由引理A1可知,存在可逆矩陣T∈Rn×n,S∈Rp×p,使得(3)成立.對于引理A2中定義的矩陣進行如下的分塊

其中,由引理A2可知,式(3)有解意味著式(A1)有解,對式(A1)按分塊展開

其中

也即有

1 Wang S H,Wang E,Dorato P.Observing the states of systems with unmeasurable disturbance.IEEE Transactions on Automatic Control,1975,20(5):716?717

2 Bhattacharyya S.Observer design for linear systems with unknown inputs.IEEE Transactions on Automatic Control,1978,23(3):483?484

3 Corless M,Tu J.State and input estimation for a class of uncertain systems.Automatica,1998,34(6):757?764

4 Charandabi B A,Marquez H J.A novel approach to unknown input fi lter design for discrete-time linear systems.Automatica,2014,50(11):2835?2839

5 Yang Jun-Qi,Chen Yan-Tao,Zhu Fang-Lai,Bu Xu-Hui.State and unknown input estimations for uncertain switched linear systems.Control and Decision,2015,30(12):2247?2252(楊俊起,陳滟濤,朱芳來,卜旭輝.不確定切換線性系統狀態和未知輸入估計方法.控制與決策,2015,30(12):2247?2252)

6 Guo S H,Zhu F L,Xu L Y.Unknown input observer design for Takagi-Sugeno fuzzy stochastic system.International Journal of Control,Automation and Systems,2015,13(4):1003?1009

7 Yang J Q,Zhu F L,Yu K J,Bu X H.Observer-based state estimation and unknown input reconstruction for nonlinear complex dynamical systems.Communications in Nonlinear Science and Numerical Simulation,2015,20(3):927?939

8 Chadli M,Karimi H R.Robust observer design for unknown inputs Takagi-Sugeno models.IEEE Transactions on Fuzzy Systems,2013,21(1):158?164

9 Lee D J,Park Y,Park Y S.RobustH∞sliding mode descriptor observer for fault and output disturbance estimation of uncertain systems.IEEE Transactions on Automatic Control,2012,57(11):2928?2934

10 Gupta M K,Tomar N K,Bhaumik S.Full-and reducedorder observer design for rectangular descriptor systems with unknown inputs.Journal of the Franklin Institute,2015,352(3):1250?1264

11 Kalsi K,Lian J M,Hui S,˙Zak S H.Sliding-mode observers for systems with unknown inputs:a high-gain approach.Automaitca,2010,46(2):347?353

12 Floquet T,Edwards C,Spurgeon S K.On sliding mode observers for systems with unknown inputs.International Journal of Adaptive Control and Signal Processing,2007,21(8?9):638?656

13 Zhu F L.State estimation and unknown input reconstruction via both reduced-order and high-order sliding mode observers.Journal of Process Control,2012,22(1):296?302

14 Park T G.Designing fault detection observers for linear systems with mismatched unknown inputs.Journal of Process Control,2013,23(8):1185?1196

15 Floquet T,Barbot J P.State and unknown input estimation for linear discrete-time systems.Automatica,2006,42(11):1883?1889

16 Guo S H,Zhu F L.Reduced-Order observer design for discrete-time descriptor systems with unknown inputs:an linear matrix inequality approach.Journal of Dynamic Systems,Measurement,and Control,2015,137(8):Article No.084503

17 Zhang L.and Shi P.l2-l∞model reduction for switched LPV systems with average dwell time.IEEE Transactions on Automatic Control,2008,53(10):2443?2448

18 Zhang L X,Shi P.Stability,l2-gain and asynchronousH∞control of discrete-time switched systems with average dwell time.IEEE Transactions on Automatic Control,2009,54(9):2192?2199

19 Lin H,Antsaklis P J.Stability and stabilizability of switched linear systems:a survey of recent results.IEEE Transactions on Automatic Control,2009,54(2):308?322

20 Zhang L X,Cui N G,Liu M,Zhao Y.Asynchronous fi ltering of discrete-time switched linear systems with average dwell time.IEEE Transactions on Circuits and Systems I:Regular Papers,2011,58(5):1109?1118

21 Zhao X D,Liu H,Zhang J F,Li H Y.Multiple-mode observer design for a class of switched linear systems.IEEE Transactions on Automation Science and Engineering,2015,12(1):272?280

22 Koenig D,Marx B,Jacquet D.Unknown input observers for switched nonlinear discrete time descriptor systems.IEEE Transactions on Automatic Control,2008,53(1):373?379

23 Bejarano F J,Pisano A.Switched observers for switched linear systems with unknown inputs.IEEE Transactions on Automatic Control,2011,56(3):681?686

24 Wang Xing-Ping.Exponential stabilization of switched time-varying linear systems.Acta Automatica Sinica,2016,42(9):1440?1444(王興平.時變線性切換系統的指數鎮定.自動化學報,2016,42(9):1440?1444)

25 Wang Xing-Ping,Song Yan-Rong,Cheng Zhao-Lin.Exponential synchronization of time-varying linear multi-agent systems with switching topology.Acta Automatica Sinica,2015,41(8):1528?1532(王興平,宋艷榮,程兆林.切換網絡下時變線性多智能體系統的指數同步.自動化學報,2015,41(8):1528?1532)

Discussions on Existence of Observers and Reduced-order Observer Design for Discrete-time Switched Systems

ZHU Fang-Lai1CAI Ming1GUO Sheng-Hui1,2

The paper discusses the unknown input observer(UIO)design issue for discrete-time switched systems with unknown inputs.First,a Lyapunov-type equivalent representation of observer matching condition for a general discretetime system is given and proved.Then,based on the stability theory of discrete-time switched system without unknown inputs,a reduced-order observer is developed for a discrete-time switched system with unknown inputs.The reducedorder observer can eliminate the in fl uences of the unknown inputs directly because of the special selection of the observer gain matrix determined by matrix block computation.Meanwhile,an algebraic unknown input reconstruction method is proposed.Finally,a simulation example is given to verify the correctness and eあectiveness of the proposed methods.

Discrete-time system,unknown input estimation,observer matching condition,switched system

Zhu Fang-Lai,Cai Ming,Guo Sheng-Hui.Discussions on existences of observers and reduced-order observer design for discrete-time switched systems.Acta Automatica Sinica,2017,43(12):2091?2099

2016-06-16 錄用日期2016-11-23

June 16,2016;accepted November 23,2016

國家自然科學基金(61573256),上海市科委“創新行動計劃”項目(16111106502),浙江省高端創新載體專項項目浙江省高端創新載體專項項目資助

Supported by National Natural Science Foundation of China(61573256),Shanghai Science and Technology Innovation Fund(16111106502),Zhejiang Province High-end Innovative Carrier Project

本文責任編委孫希明

Recommended by Associate Editor SUN Xi-Ming

1.同濟大學電子與信息工程學院 上海201804 2.蘇州科技大學電子與信息工程學院蘇州215009

1.College of Electronics and Information Engineering,Tongji University,Shanghai 201804 2.College of Electrical Engineering and Automation,Suzhou University of Science and Technology,Suzhou 215009

朱芳來,蔡明,郭勝輝.離散切換系統觀測器存在性討論及降維觀測器設計.自動化學報,2017,43(12):2091?2099

DOI10.16383/j.aas.2017.c160471

朱芳來 同濟大學電子與信息工程學院教授.主要研究方向為非線性系統魯棒控制,觀測器設計,基于模型的故障檢測與隔離.本文通信作者.

E-mail:zhufanglai@tongji.edu.cn

(ZHU Fang-Lai Professor at the College of Electronics and Information Engineering,Tongji University.His research interest covers nonlinear robust control,modelbased fault detection and isolation.Corresponding author of this paper.)

蔡 明 同濟大學控制理論與控制工程專業碩士研究生.主要研究方向為基于模型的故障檢測與隔離.

E-mail:basycai@163.com

(CAI Ming Master sutdent at the College of Electronics and Information Engineering,Tongji University.His research interest covers model-based fault detection and isolation.)

郭勝輝 蘇州科技大學電子與信息工程學院講師,博士.主要研究方向為基于模型的故障檢測及觀測器設計.

E-mail:12gsh@tongji.edu.cn

(GUO Sheng-Hui Ph.D,lecturer at the College of Electronics and Information Engineering,Suzhou University of Science and Technology.His research interest covers model-based fault detection and observer design.)

猜你喜歡
系統設計
Smartflower POP 一體式光伏系統
工業設計(2022年8期)2022-09-09 07:43:20
WJ-700無人機系統
ZC系列無人機遙感系統
北京測繪(2020年12期)2020-12-29 01:33:58
何為設計的守護之道?
現代裝飾(2020年7期)2020-07-27 01:27:42
《豐收的喜悅展示設計》
流行色(2020年1期)2020-04-28 11:16:38
基于PowerPC+FPGA顯示系統
半沸制皂系統(下)
瞞天過海——仿生設計萌到家
藝術啟蒙(2018年7期)2018-08-23 09:14:18
連通與提升系統的最后一塊拼圖 Audiolab 傲立 M-DAC mini
設計秀
海峽姐妹(2017年7期)2017-07-31 19:08:17
主站蜘蛛池模板: 国产精品女熟高潮视频| 国产欧美性爱网| 亚洲欧美极品| 性激烈欧美三级在线播放| 这里只有精品在线| 成年女人18毛片毛片免费| 欧美日韩精品综合在线一区| 欧美区日韩区| 狂欢视频在线观看不卡| 日韩精品毛片人妻AV不卡| 在线色国产| 四虎国产在线观看| 久久久久久久久18禁秘| 99热这里只有免费国产精品| 99久久精品免费看国产免费软件 | 毛片在线播放网址| 国产一区二区三区视频| 日韩性网站| 一级高清毛片免费a级高清毛片| 亚洲一区二区成人| 久久久久国产一级毛片高清板| 久久久久久国产精品mv| 国产午夜精品鲁丝片| 亚洲成人动漫在线| 亚洲成人网在线观看| 日本在线欧美在线| 国产午夜小视频| 久久综合婷婷| 国产日韩AV高潮在线| 日韩成人高清无码| 熟妇无码人妻| 久久人体视频| 精品一区国产精品| 午夜国产大片免费观看| 538国产视频| jizz国产视频| 国产成人免费观看在线视频| 欧洲熟妇精品视频| 日韩一级毛一欧美一国产| 国产精品美人久久久久久AV| 精品国产一二三区| 毛片免费网址| 国产成人一区在线播放| 人妻中文久热无码丝袜| 女人一级毛片| 国产永久免费视频m3u8| 欧洲亚洲欧美国产日本高清| 白浆免费视频国产精品视频| 国产玖玖视频| 国产小视频免费| 一区二区无码在线视频| 久久精品国产999大香线焦| 一区二区在线视频免费观看| 久久这里只有精品66| 欧美日韩va| 欧美精品高清| AV不卡国产在线观看| 黄色网站不卡无码| 综合亚洲色图| 在线观看欧美国产| 久久久久无码精品| 亚洲—日韩aV在线| 成人亚洲国产| 999国内精品久久免费视频| 在线99视频| 美女国内精品自产拍在线播放 | 亚洲成人黄色在线观看| 国产精品久久久久久搜索| 日本中文字幕久久网站| 在线看片免费人成视久网下载| 免费欧美一级| 伊人成人在线视频| 亚洲高清中文字幕在线看不卡| 国产麻豆精品手机在线观看| 九九九九热精品视频| 国产成人综合久久| 久久国产精品77777| 91区国产福利在线观看午夜| 亚洲人在线| 成人在线欧美| 国产成人禁片在线观看| 欧美影院久久|