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

    基于新息估計(jì)和正交投影的閉環(huán)子空間模型辨識(shí)

    2016-12-17 08:23:44侯杰劉濤
    自動(dòng)化學(xué)報(bào) 2016年11期
    關(guān)鍵詞:新息閉環(huán)方差

    侯杰 劉濤

    基于新息估計(jì)和正交投影的閉環(huán)子空間模型辨識(shí)

    侯杰1劉濤1

    針對(duì)閉環(huán)控制系統(tǒng)提出一種基于新息估計(jì)和正交投影的閉環(huán)子空間模型辨識(shí)方法.首先采用最小二乘法對(duì)VARX模型(Vector autoregressive with exogenous inputs model)進(jìn)行計(jì)算得到新息估計(jì)值,然后通過將由觀測(cè)輸入輸出數(shù)據(jù)構(gòu)造的Hankel矩陣正交投影到新息數(shù)據(jù)的正交補(bǔ)空間以消除噪聲影響,從而在無噪聲的輸入輸出數(shù)據(jù)奇偶空間中提取得到擴(kuò)展可觀測(cè)矩陣和下三角形Toeplitz矩陣.最后采用平移變換法得到系統(tǒng)矩陣.對(duì)該算法嚴(yán)格分析和證明了實(shí)現(xiàn)一致估計(jì)的條件.通過仿真實(shí)例驗(yàn)證了本文方法的有效性和優(yōu)越性.

    子空間辨識(shí),閉環(huán)辨識(shí),奇偶空間,新息估計(jì),正交投影,一致估計(jì)

    子空間辨識(shí)算法由于能對(duì)多輸入多輸出系統(tǒng)采用統(tǒng)一框架建立狀態(tài)空間模型,在系統(tǒng)辨識(shí)和控制工程領(lǐng)域受到廣泛關(guān)注[1].有一些子空間算法被提出用于開環(huán)辨識(shí)工業(yè)過程,得到一致估計(jì)結(jié)果[2].但是,由于過程操作安全性和穩(wěn)定性的需要,許多工業(yè)過程限制在閉環(huán)條件下進(jìn)行辨識(shí)實(shí)驗(yàn),由于反饋控制作用的影響,使得過程噪聲和輸入存在相關(guān)性,使得傳統(tǒng)開環(huán)子空間方法產(chǎn)生辨識(shí)偏差[3].閉環(huán)系統(tǒng)辨識(shí)因此在近年來受到很多關(guān)注和探討[4],一些閉環(huán)子空間辨識(shí)算法被相繼提出[5].這些算法可被歸納為三類,第一類方法[6?7]采用輔助變量策略消除噪聲影響,保證估計(jì)結(jié)果的一致性;第二類方法[8]采用最小二乘法對(duì)高階VARX模型(Vector autoregressive with exogenous inputs model)進(jìn)行計(jì)算得到馬爾科夫估計(jì)參數(shù),由于VARX模型只包括當(dāng)前時(shí)刻的不可測(cè)噪聲,該噪聲和VARX模型的過去時(shí)刻輸入無關(guān),從而可保證所得參數(shù)的一致性;第三類算法[9]用噪聲預(yù)估值代替真實(shí)值進(jìn)行計(jì)算保證得到一致估計(jì)結(jié)果.

    基于奇偶空間的閉環(huán)子空間辨識(shí)算法SIMPCAwc[6]采用過去時(shí)刻輸入輸出數(shù)據(jù)作為輔助變量來消除噪聲,以得到無噪聲的輸入輸出數(shù)據(jù),然后從無噪聲影響的輸入輸出數(shù)據(jù)奇偶空間中提取得到擴(kuò)展可觀測(cè)矩陣和下三角形Toeplitz矩陣,從而求得系統(tǒng)矩陣,該方法取得較好辨識(shí)精度.然而,文獻(xiàn)[7]指出當(dāng)閉環(huán)系統(tǒng)設(shè)定點(diǎn)輸入激勵(lì)為不相關(guān)白噪聲序列時(shí),雖然引入輔助變量與噪聲不相關(guān),可以有效地消除噪聲,但由于該輔助變量和系統(tǒng)設(shè)定點(diǎn)輸入激勵(lì)也不相關(guān),導(dǎo)致從無噪聲輸入輸出數(shù)據(jù)奇偶空間中提取參數(shù)可能同時(shí)含有過程模型參數(shù)信息和控制器參數(shù)信息,因而無法對(duì)它們進(jìn)行區(qū)分,從而致使過程模型估計(jì)出現(xiàn)偏差.針對(duì)SIMPCAwc[6]辨識(shí)方法存在的問題,本文通過將輸入輸出數(shù)據(jù)正交投影到新息數(shù)據(jù)的正交補(bǔ)空間來消除噪聲,以得到新的無噪聲數(shù)據(jù)矩陣,進(jìn)而從其對(duì)應(yīng)的奇偶空間中提取得到擴(kuò)展可觀測(cè)矩陣和下三角形Toeplitz矩陣.由于新息數(shù)據(jù)的正交補(bǔ)空間數(shù)據(jù)和噪聲無關(guān),且同時(shí)與系統(tǒng)設(shè)定輸入激勵(lì)相關(guān),確保本文方法從新的無噪聲奇偶空間中提取的參數(shù)只包含過程模型參數(shù),有效地保證估計(jì)結(jié)果的一致無偏性.由于新息數(shù)據(jù)為不可測(cè)量數(shù)據(jù),本文通過模型推導(dǎo),得到和待辨識(shí)狀態(tài)空間模型等價(jià)的VARX模型.在此基礎(chǔ)上,采用最小二乘法對(duì)VARX模型進(jìn)行計(jì)算以得到新息的一致估計(jì)值.采用新息一致估計(jì)值代替真實(shí)值,以完成模型參數(shù)估計(jì).為了論證說明本文方法的有效性,嚴(yán)格分析和給出了本文算法保證一致估計(jì)的條件.

    1 問題描述

    本文研究如下線性離散狀態(tài)空間過程模型:

    其中,x(t)∈Rnx,u(t)∈Rnu,y(t)∈Rny分別為系統(tǒng)狀態(tài)和過程輸入和輸出.v(t)∈Rny和w(t)∈Rnx分別為過程測(cè)量噪聲.A,B,C,D分別為相應(yīng)維數(shù)的系統(tǒng)矩陣.本文研究系統(tǒng)在閉環(huán)工作條件下,利用系統(tǒng)輸入和輸出觀測(cè)數(shù)據(jù),辨識(shí)對(duì)象狀態(tài)空間(亦稱子空間)模型.

    由于閉環(huán)反饋控制作用的影響,使得過程測(cè)量噪聲和輸入存在相關(guān)性,若直接通過模型(1)來辨識(shí)系統(tǒng)矩陣,很難消除噪聲對(duì)辨識(shí)結(jié)果的不利影響.因此,采用卡爾曼濾波原理[10],將系統(tǒng)模型(1)等價(jià)表示為新息形式

    其中,K為卡爾曼濾波增益,新息e(t)為零均值白噪聲,當(dāng)i<j時(shí),新息e(j)和輸入輸出{u(i),y(i)}不相關(guān).

    定義過去和將來水平數(shù)分別為p和f,過去和將來輸入向量分別為up(t)=[u(t?p)T···u(t?2)Tu(t?1)T]T和uf(t)=[u(t)T···u(t+f?2)Tu(t+f?1)T]T,定義過去和將來輸入Hankel矩陣Up=[up(t)T···up(N)T]T和Uf=[uf(t)T··· uf(N)T]T,輸出和新息數(shù)據(jù)做類似定義.

    對(duì)式(3)進(jìn)行迭代可得:

    本文將采用模型(2)和(5)對(duì)系統(tǒng)矩陣進(jìn)行辨識(shí).定義Xp=[x(t?p)···x(t?p+N?1)]和Xf=[x(t)···x(t+N?1)].通過對(duì)式(2)進(jìn)行迭代可得:

    其中,擴(kuò)展可觀測(cè)矩陣為Γf= [CT··· (CAf?1)T]T.下三角形Toeplitz矩陣分別為

    2 閉環(huán)子空間辨識(shí)算法和一致性分析

    2.1 閉環(huán)子空間辨識(shí)算法

    在式(6)的基礎(chǔ)上,通過同時(shí)計(jì)算擴(kuò)展可觀測(cè)矩陣Γf和下三角形Toeplitz矩陣Hf實(shí)現(xiàn)對(duì)系統(tǒng)矩陣的辨識(shí).首先將輸入數(shù)據(jù)移至式(6)的左側(cè),得到:

    為求解式(7)得到Γf和Hf的估計(jì)值,需要消除未知狀態(tài)和新息的影響.通過在式(7)的左右側(cè)同時(shí)引入Γf的正交補(bǔ)向量由于

    通過將式(8)正交投影到Ef的正交補(bǔ)空間來消除新息噪聲,

    采用最小二乘法可得Hf1的一致估計(jì)如下:

    本文基于新息估計(jì)和正交投影的閉環(huán)子空間辨識(shí)方法(Closed-loop subspace identification method using innovation estimation and orthogonal projection,CSIMIEOP)可總結(jié)如下:

    步驟2.通過SVD分解對(duì)式(15)進(jìn)行計(jì)算.

    步驟4.通過式(24),式(25)和式(27)求解系統(tǒng)矩陣

    2.2 閉環(huán)一致條件分析

    由第2.1節(jié)閉環(huán)子空間辨識(shí)算法可知,本文采用新息估計(jì)和正交投影消除噪聲,得到一致估計(jì)結(jié)果.辨識(shí)結(jié)果是否一致取決于是否一致,本文對(duì)它們的一致估計(jì)條件進(jìn)行分析和說明,給出如下定理.

    定理1.若

    由于系統(tǒng)為可控可觀測(cè)系統(tǒng),采用文獻(xiàn)[11]定理2可知Υ為滿秩矩陣.同時(shí)如果式(33)第2項(xiàng)為正定滿秩矩陣,根據(jù)文獻(xiàn)[12]給定的秩條件可得:

    因此,矩陣(32)為滿秩矩陣.

    由式(30)、式(32)和式(35)可知

    2.3 閉環(huán)一致條件合理性分析

    第2.2節(jié)定理2的兩個(gè)假設(shè)條件涉及未知新息信息,以上假設(shè)是否合理,將直接決定本文方法能否取得一致估計(jì)結(jié)果.

    由于模型(1)可表示為模型(2),則新息的協(xié)方差可表示為(具體表述可參考文獻(xiàn)[10]第5章),其中P為半正定矩陣,對(duì)于系統(tǒng)噪聲有則可知

    對(duì)于假設(shè)2,由于

    進(jìn)一步,由式(2)可知

    其中,擴(kuò)展可觀性矩陣分別表示為L(zhǎng)1=[Ap?1B ···AB B],L2=[Ap?1K···AK K],下三角形Toeplitz矩陣為

    3 仿真研究

    考慮文獻(xiàn)[6]中研究的閉環(huán)系統(tǒng)

    其中,反饋控制結(jié)構(gòu)為u(t)=?0.6y(t)+r(t).過程噪聲設(shè)置為方差為0.2的白噪聲序列.

    對(duì)系統(tǒng)設(shè)定輸入激勵(lì)r(k)為單位方差白噪聲序列和相關(guān)序列兩種情況進(jìn)行研究.其中相關(guān)序列設(shè)為r(t)=(1+0.8q?1+0.6q?2)r0(t),r0(t)為單位方差白噪聲序列.過去和未來水平數(shù)均設(shè)置為10.在不同數(shù)據(jù)長(zhǎng)度N∈[200,8000]的情況下進(jìn)行1000次Monte Carlo仿真.本文方法的辨識(shí)結(jié)果與SIMPCAwc算法[6]進(jìn)行比較,當(dāng)r(t)為單位方差白噪聲系列時(shí),系統(tǒng)極點(diǎn)(真實(shí)值為0.9)的估計(jì)平均值如圖1所示,系統(tǒng)極點(diǎn)的估計(jì)標(biāo)準(zhǔn)方差如圖2所示.當(dāng)系統(tǒng)設(shè)定輸入r(t)為相關(guān)系列時(shí),系統(tǒng)極點(diǎn)的平均值如圖3所示,系統(tǒng)極點(diǎn)的估計(jì)標(biāo)準(zhǔn)方差如圖4所示.

    圖1 系統(tǒng)設(shè)定輸入激勵(lì)為白噪聲的系統(tǒng)極點(diǎn)估計(jì)平均值Fig.1 Mean value of estimated poles for white noise setpoint excitation

    圖2 系統(tǒng)設(shè)定輸入激勵(lì)為白噪聲的系統(tǒng)極點(diǎn)估計(jì)標(biāo)準(zhǔn)方差Fig.2 Standard deviation of estimated poles for white noise setpoint excitation

    圖3 系統(tǒng)設(shè)定輸入激勵(lì)為相關(guān)序列的系統(tǒng)極點(diǎn)估計(jì)平均值Fig.3 Mean value of estimated poles for correlated quasi-stationary setpoint excitation

    圖4 系統(tǒng)設(shè)定輸入激勵(lì)為相關(guān)序列的系統(tǒng)極點(diǎn)估計(jì)標(biāo)準(zhǔn)方差Fig.4 Standard deviation of estimated poles for correlated quasi-stationary setpoint excitation

    從以上結(jié)果可以看出,當(dāng)系統(tǒng)設(shè)定輸入為單位方差白噪聲系列時(shí),SIMPCAwc算法得出有偏估計(jì)結(jié)果;只有當(dāng)系統(tǒng)設(shè)定輸入激勵(lì)為相關(guān)序列時(shí),SIMPCAwc算法才能保證無偏估計(jì)結(jié)果.本文CSIMIEOP算法對(duì)系統(tǒng)設(shè)定輸入激勵(lì)為無關(guān)序列和相關(guān)序列時(shí)均可得到一致無偏結(jié)果,并且估計(jì)精度優(yōu)于SIMPCAwc算法.

    4 結(jié)論

    本文提出一種基于新息估計(jì)和正交投影的閉環(huán)子空間辨識(shí)算法,對(duì)系統(tǒng)設(shè)定點(diǎn)輸入激勵(lì)為白噪聲無關(guān)序列和相關(guān)序列的情況均可得到一致無偏估計(jì)結(jié)果,并且相對(duì)于近期有關(guān)文獻(xiàn)給出的方法如SIMPCAwc算法[6],能進(jìn)一步提高辨識(shí)精度.同時(shí),嚴(yán)格分析和證明了本文算法保證一致估計(jì)的條件.最后通過仿真實(shí)例驗(yàn)證了本文方法的有效性和優(yōu)越性.

    1 Qin S J.An overview of subspace identification.Computers Chemical and Engineering,2006,30(10?12):1502?1513

    2 Favoreel W,de Moor B,van Overschee P.Subspace state space system identification for industrial processes.Journal of Process Control,2000,10(2?3):149?155

    3 Yang Hua,Li Shao-Yuan.Closed-loop subspace identification based on augmented input with consistency analysis. Acta Automatica Sinica,2007,33(7):703?708 (楊華,李少遠(yuǎn).基于輸入擴(kuò)張的閉環(huán)系統(tǒng)子空間辨識(shí)及其強(qiáng)一致性分析.自動(dòng)化學(xué)報(bào),2007,33(7):703?708)

    4 Wang Le-Yi,Zhao Wen-Xiao.System identification:new paradigms,challenges,and opportunities.Acta Automatica Sinica,2013,39(7):933?942 (王樂一,趙文九虎.系統(tǒng)辨識(shí):新的模式、挑戰(zhàn)及機(jī)遇.自動(dòng)化學(xué)報(bào), 2013,39(7):933?942)

    5 van der Veen G J,van Wingerden J W,Bergamasco M, Lovera M,Verhaegen M.Closed-loop subspace identification methods:an overview.IET Control Theory and Applications,2013,7(10):1339?1358

    6 Wang J,Qin S J.Closed-loop subspace identification using the parity space.Automatica,2006,42(2):315?320

    7 Huang B,Ding S X,Qin S J.Closed-loop subspace identification:an orthogonal projection approach.Journal of Process Control,2005,15(1):53?66

    8 Chiuso A.The role of vector autoregressive modeling in predictor-based subspace identification.Automatica,2007, 43(6):1034?1048

    9 Qin S J,Ljung L.Closed-loop subspace identification with innovation estimation.In:Proceedings of the 13th IFAC System Identification.Rotterdam,Netherlands: IFAC, 2003.887?892

    10 Katayama T.Subspace Methods for System Identification. London:Springer-Verlag,2005,107?137

    11 Jansson M,Wahlberg B.On consistency of subspace methods for system identification.Automatica,1998,34(12): 1507?1519

    12 Liu T,Huang B,Qin S J.Bias-eliminated subspace model identification under time-varying deterministic type load disturbance.Journal of Process Control,2015,25:41?49

    13 Wang J,Qin S J.A new subspace identification approach based on principal component analysis.Journal of Process Control,2002,12(8):841?855

    侯 杰 大連理工大學(xué)控制科學(xué)與工程學(xué)院博士研究生.主要研究方向?yàn)橄到y(tǒng)辨識(shí). E-mail:jiehou.phd@hotmail.com

    (HOU Jie Ph.D.candidate at the School of Control Science and Engineering,Dalian University of Technology.His main research interest is system identification.)

    劉 濤 大連理工大學(xué)控制科學(xué)與工程學(xué)院教授.2012年入選國(guó)家青年千人計(jì)劃,德國(guó)洪堡基金學(xué)者.主要研究方向?yàn)楣I(yè)過程辨識(shí)建模,魯棒過程控制,批量生產(chǎn)過程控制優(yōu)化,過程質(zhì)量監(jiān)測(cè).本文通信作者.E-mail:liurouter@ieee.org

    (LIU Tao Professor at the School of Control Science and Engineering, Dalian University of Technology.He was a recipient of National Thousand Talents Program of China in 2012,and a Humboldt research fellow of Germany.His research interest covers industrial process identification and modeling, robust process control,batch process control and optimization,and process quality monitoring.Corresponding author of this paper.)

    Closed-loop Subspace Model Identification Using Innovation Estimation and Orthogonal Projection

    HOU Jie1LIU Tao1

    In this paper,a closed-loop subspace model identification method using innovation estimation and orthogonal projection is proposed for closed-loop control systems.Firstly,a least-squares algorithm is adopted to estimate the innovation matrix via the vector autoregressive with exogenous inputs(VARX)model.Then,by performing an orthogonal projection of the observed input-output Hankel matrix onto the orthogonal complement space of innovation Hankel matrix to eliminate the influence from noise,the extended observability matrix and lower triangular block Toeplitz matrix are derived from the parity space of noise-free input-output data.Finally,the system matrices are retrieved by using a shiftinvariant approach.The consistent estimation conditions are analyzed with a strict proof.A simulation example is shown to demonstrate the effectiveness and merit of the proposed method.

    Subspace identification,closed-loop identification,parity space,innovation estimation,orthogonal projection,consistent estimation

    侯杰,劉濤.基于新息估計(jì)和正交投影的閉環(huán)子空間模型辨識(shí).自動(dòng)化學(xué)報(bào),2016,42(11):1657?1663

    Hou Jie,Liu Tao.Closed-loop subspace model identification using innovation estimation and orthogonal projection.Acta Automatica Sinica,2016,42(11):1657?1663

    2016-01-19 錄用日期2016-03-10

    Manuscript received January 19,2016;accepted March 10,2016

    國(guó)家自然科學(xué)基金(61473054),第三批國(guó)家青年千人計(jì)劃,中央高?;究蒲袠I(yè)務(wù)費(fèi)重點(diǎn)培育基金(DUT15ZD108)資助

    Supported by National Natural Science Foundation of China (61473054),National Thousand Talents Program of China,and Fundamental Research Funds for the Central Universities of China(DUT15ZD108)

    本文責(zé)任編委方海濤

    Recommended by Associate Editor FANG Hai-Tao

    1.大連理工大學(xué)控制科學(xué)與工程學(xué)院大連116024

    1.School of Control Science and Engineering,Dalian University of Technology,Dalian 116024

    DOI 10.16383/j.aas.2016.c160050

    猜你喜歡
    新息閉環(huán)方差
    方差怎么算
    傳遞函數(shù)辨識(shí)(21):線性回歸系統(tǒng)的遞階遞推參數(shù)估計(jì)
    概率與統(tǒng)計(jì)(2)——離散型隨機(jī)變量的期望與方差
    計(jì)算方差用哪個(gè)公式
    M估計(jì)的強(qiáng)跟蹤SVD-UKF算法在組合導(dǎo)航中的應(yīng)用
    電子科技(2018年7期)2018-07-23 05:30:32
    方差生活秀
    單周期控制下雙輸入Buck變換器閉環(huán)系統(tǒng)設(shè)計(jì)
    黑龍江電力(2017年1期)2017-05-17 04:25:05
    雙閉環(huán)模糊控制在石化廢水處理中的研究
    自適應(yīng)卡爾曼濾波在航空重力異常解算的應(yīng)用研究
    基于新息正交性自適應(yīng)濾波的慣性/地磁組合導(dǎo)航方法
    婷婷成人精品国产| 最新在线观看一区二区三区 | 国产欧美亚洲国产| 熟女av电影| 伦理电影大哥的女人| 欧美xxⅹ黑人| 欧美日韩亚洲综合一区二区三区_| 国产成人午夜福利电影在线观看| avwww免费| 在线观看人妻少妇| 亚洲av综合色区一区| 777久久人妻少妇嫩草av网站| 日韩欧美精品免费久久| 亚洲精品成人av观看孕妇| 热99国产精品久久久久久7| 超碰成人久久| 曰老女人黄片| 一二三四在线观看免费中文在| 秋霞在线观看毛片| 最近2019中文字幕mv第一页| 下体分泌物呈黄色| 亚洲av综合色区一区| 波多野结衣av一区二区av| 老司机影院毛片| 99香蕉大伊视频| 人人妻人人添人人爽欧美一区卜| 黄色 视频免费看| 午夜福利,免费看| 十八禁人妻一区二区| 亚洲精品国产色婷婷电影| 国产 精品1| 丝袜喷水一区| 中文天堂在线官网| 亚洲国产精品一区三区| 国产乱来视频区| 国产午夜精品一二区理论片| 黄片无遮挡物在线观看| 黄色毛片三级朝国网站| 97在线人人人人妻| 亚洲精品aⅴ在线观看| 大码成人一级视频| 狂野欧美激情性bbbbbb| 国产成人一区二区在线| 操出白浆在线播放| 美国免费a级毛片| 新久久久久国产一级毛片| 亚洲人成77777在线视频| 亚洲欧洲国产日韩| 我要看黄色一级片免费的| 精品亚洲成国产av| 人成视频在线观看免费观看| av一本久久久久| 久热爱精品视频在线9| 精品人妻在线不人妻| 伊人亚洲综合成人网| 国产精品免费大片| 日本黄色日本黄色录像| 在线观看免费午夜福利视频| 午夜日韩欧美国产| 国产深夜福利视频在线观看| 你懂的网址亚洲精品在线观看| 99热网站在线观看| 五月天丁香电影| 亚洲第一av免费看| 亚洲第一区二区三区不卡| 少妇的丰满在线观看| 亚洲自偷自拍图片 自拍| 日韩制服丝袜自拍偷拍| 国产熟女欧美一区二区| √禁漫天堂资源中文www| 97人妻天天添夜夜摸| 成人午夜精彩视频在线观看| 日韩大码丰满熟妇| 国产免费视频播放在线视频| 无遮挡黄片免费观看| 国产伦人伦偷精品视频| 国产视频首页在线观看| 中文字幕制服av| 性高湖久久久久久久久免费观看| 国产精品人妻久久久影院| 免费日韩欧美在线观看| 18禁国产床啪视频网站| 中文字幕人妻熟女乱码| 99久久99久久久精品蜜桃| 高清在线视频一区二区三区| 一级片免费观看大全| 伦理电影免费视频| 久久久久久人妻| 亚洲成色77777| 日韩一卡2卡3卡4卡2021年| 亚洲欧美清纯卡通| 老熟女久久久| 高清在线视频一区二区三区| 国产极品粉嫩免费观看在线| 欧美日韩视频精品一区| 99精品久久久久人妻精品| 欧美日韩成人在线一区二区| 欧美日韩视频精品一区| 国产成人欧美在线观看 | 国产探花极品一区二区| 成人免费观看视频高清| 天堂中文最新版在线下载| 久久久亚洲精品成人影院| 观看av在线不卡| 亚洲一级一片aⅴ在线观看| 精品久久久精品久久久| 久久婷婷青草| 婷婷色综合大香蕉| 好男人视频免费观看在线| 天天躁狠狠躁夜夜躁狠狠躁| 欧美在线一区亚洲| 久久ye,这里只有精品| 性色av一级| 国产成人精品福利久久| 成人三级做爰电影| 老司机靠b影院| 国产精品一区二区在线观看99| 国产一区二区在线观看av| 我要看黄色一级片免费的| 亚洲精品久久成人aⅴ小说| 日韩大码丰满熟妇| 精品国产国语对白av| 久久国产精品男人的天堂亚洲| 一级,二级,三级黄色视频| 久久综合国产亚洲精品| 国产亚洲av片在线观看秒播厂| 永久免费av网站大全| 国产人伦9x9x在线观看| 韩国精品一区二区三区| 久久久久久久久久久久大奶| 黑人欧美特级aaaaaa片| 精品人妻一区二区三区麻豆| 国产黄色免费在线视频| 国产免费一区二区三区四区乱码| 欧美黑人精品巨大| 亚洲av成人精品一二三区| 欧美人与性动交α欧美软件| 亚洲精华国产精华液的使用体验| 性少妇av在线| 99热全是精品| 国产在线免费精品| 夫妻午夜视频| 亚洲精品一区蜜桃| 亚洲精品久久久久久婷婷小说| 看免费av毛片| 青春草视频在线免费观看| 汤姆久久久久久久影院中文字幕| 美女福利国产在线| 在线观看免费视频网站a站| 中文欧美无线码| 色播在线永久视频| 另类亚洲欧美激情| 少妇人妻久久综合中文| 国产日韩欧美视频二区| 欧美日韩亚洲综合一区二区三区_| 国产男女内射视频| 18禁观看日本| 中文天堂在线官网| 日本爱情动作片www.在线观看| 女人高潮潮喷娇喘18禁视频| 麻豆av在线久日| 亚洲,一卡二卡三卡| 久久久国产欧美日韩av| 午夜av观看不卡| 欧美日韩国产mv在线观看视频| 精品一品国产午夜福利视频| 久久这里只有精品19| 国产片特级美女逼逼视频| 午夜福利,免费看| 免费不卡黄色视频| 黄色视频在线播放观看不卡| 波野结衣二区三区在线| 久久综合国产亚洲精品| 精品国产一区二区三区四区第35| 美女主播在线视频| 一级a爱视频在线免费观看| 黄色视频不卡| 久久久国产一区二区| 婷婷色av中文字幕| 纵有疾风起免费观看全集完整版| 中文字幕av电影在线播放| 高清在线视频一区二区三区| 国产精品无大码| 一区二区三区乱码不卡18| 久久久精品免费免费高清| 国产欧美日韩综合在线一区二区| 最近2019中文字幕mv第一页| 人体艺术视频欧美日本| 久久99热这里只频精品6学生| 久久久国产欧美日韩av| 欧美日韩一级在线毛片| 精品少妇一区二区三区视频日本电影 | 亚洲欧美一区二区三区黑人| 久久青草综合色| 9191精品国产免费久久| av视频免费观看在线观看| 人人妻人人澡人人看| 美女国产高潮福利片在线看| 色94色欧美一区二区| 亚洲欧洲日产国产| 久久久久久人妻| 国产精品 国内视频| 久久韩国三级中文字幕| 亚洲国产精品999| 亚洲av成人精品一二三区| 国产精品免费视频内射| 免费观看av网站的网址| 卡戴珊不雅视频在线播放| 久久精品aⅴ一区二区三区四区| 哪个播放器可以免费观看大片| av不卡在线播放| 亚洲精品美女久久久久99蜜臀 | 在现免费观看毛片| 成人亚洲精品一区在线观看| 国产精品女同一区二区软件| 日韩欧美一区视频在线观看| 国产精品秋霞免费鲁丝片| 亚洲欧美精品综合一区二区三区| av.在线天堂| 日韩成人av中文字幕在线观看| 亚洲第一av免费看| 免费不卡黄色视频| 成年av动漫网址| 麻豆av在线久日| 亚洲国产最新在线播放| 啦啦啦 在线观看视频| 欧美中文综合在线视频| 免费黄网站久久成人精品| 久久免费观看电影| 免费在线观看视频国产中文字幕亚洲 | 国产精品国产三级国产专区5o| 国产成人一区二区在线| 久久精品久久久久久噜噜老黄| www.精华液| av卡一久久| 天堂俺去俺来也www色官网| 国产成人午夜福利电影在线观看| 叶爱在线成人免费视频播放| 观看av在线不卡| 黄片播放在线免费| 一级黄片播放器| 美女中出高潮动态图| 最近中文字幕2019免费版| 亚洲天堂av无毛| 精品第一国产精品| 日韩一区二区三区影片| 男女无遮挡免费网站观看| 精品国产露脸久久av麻豆| 欧美日韩福利视频一区二区| xxxhd国产人妻xxx| 纯流量卡能插随身wifi吗| 黄色视频不卡| 在线观看免费午夜福利视频| 人妻人人澡人人爽人人| 亚洲国产欧美网| 日韩成人av中文字幕在线观看| 丰满饥渴人妻一区二区三| 欧美激情极品国产一区二区三区| 这个男人来自地球电影免费观看 | 麻豆精品久久久久久蜜桃| 成年人免费黄色播放视频| 制服诱惑二区| 中文字幕高清在线视频| 久久国产精品男人的天堂亚洲| 女性生殖器流出的白浆| 精品一区二区免费观看| 中文字幕高清在线视频| 九九爱精品视频在线观看| 女性生殖器流出的白浆| 久久久精品国产亚洲av高清涩受| 精品少妇内射三级| 亚洲精品一区蜜桃| 精品国产露脸久久av麻豆| 最新在线观看一区二区三区 | 黄片播放在线免费| 热re99久久精品国产66热6| 人妻 亚洲 视频| 亚洲av中文av极速乱| 国产精品一区二区在线不卡| 青春草国产在线视频| 十八禁网站网址无遮挡| www.熟女人妻精品国产| 两个人免费观看高清视频| 成人国语在线视频| 极品少妇高潮喷水抽搐| 日韩一卡2卡3卡4卡2021年| 亚洲国产欧美日韩在线播放| 国产一区二区三区av在线| 五月天丁香电影| 人体艺术视频欧美日本| 精品久久久精品久久久| 免费高清在线观看日韩| 9191精品国产免费久久| 午夜福利在线免费观看网站| 国产成人免费无遮挡视频| 尾随美女入室| 成人18禁高潮啪啪吃奶动态图| 精品视频人人做人人爽| av网站在线播放免费| 日本91视频免费播放| 日本wwww免费看| 国产精品久久久久久久久免| 日本午夜av视频| 欧美成人午夜精品| 免费人妻精品一区二区三区视频| 欧美成人午夜精品| 欧美日韩成人在线一区二区| 国产精品一二三区在线看| 国产精品 国内视频| 热99国产精品久久久久久7| 国产av精品麻豆| av卡一久久| 欧美激情极品国产一区二区三区| 亚洲精品av麻豆狂野| 一本一本久久a久久精品综合妖精| svipshipincom国产片| 日本欧美视频一区| 亚洲av福利一区| 国产亚洲av片在线观看秒播厂| 午夜老司机福利片| 国产精品蜜桃在线观看| 最黄视频免费看| 超碰97精品在线观看| 汤姆久久久久久久影院中文字幕| 日韩中文字幕欧美一区二区 | 亚洲欧洲精品一区二区精品久久久 | 一级毛片黄色毛片免费观看视频| 丰满少妇做爰视频| 欧美最新免费一区二区三区| 日日摸夜夜添夜夜爱| 亚洲国产日韩一区二区| 啦啦啦中文免费视频观看日本| 国产精品二区激情视频| 日本av手机在线免费观看| 欧美激情 高清一区二区三区| 激情五月婷婷亚洲| 亚洲视频免费观看视频| 中国国产av一级| 国产精品国产av在线观看| 高清av免费在线| 国产色婷婷99| 国产精品免费大片| 久久毛片免费看一区二区三区| 国产亚洲午夜精品一区二区久久| 欧美激情 高清一区二区三区| 日韩av免费高清视频| 青春草视频在线免费观看| 日日撸夜夜添| 天天躁狠狠躁夜夜躁狠狠躁| 自拍欧美九色日韩亚洲蝌蚪91| 丝袜喷水一区| 亚洲国产欧美在线一区| avwww免费| 国产亚洲欧美精品永久| h视频一区二区三区| 亚洲精品国产一区二区精华液| 乱人伦中国视频| 亚洲欧洲国产日韩| 男女午夜视频在线观看| 亚洲成国产人片在线观看| 国产亚洲av片在线观看秒播厂| 女人精品久久久久毛片| 亚洲国产毛片av蜜桃av| 2018国产大陆天天弄谢| 考比视频在线观看| 国产日韩一区二区三区精品不卡| 中文字幕人妻丝袜一区二区 | 精品国产露脸久久av麻豆| 国语对白做爰xxxⅹ性视频网站| 精品国产一区二区三区久久久樱花| 大片免费播放器 马上看| 午夜福利在线免费观看网站| 无限看片的www在线观看| 国产有黄有色有爽视频| 国产老妇伦熟女老妇高清| 男女之事视频高清在线观看 | 久久人人爽人人片av| 亚洲av中文av极速乱| 韩国av在线不卡| 超色免费av| 亚洲欧美日韩另类电影网站| 街头女战士在线观看网站| 国产精品久久久久久精品古装| 国产精品一国产av| www.自偷自拍.com| 亚洲一级一片aⅴ在线观看| 丝袜喷水一区| 国产精品三级大全| 亚洲国产精品999| 午夜免费鲁丝| 1024香蕉在线观看| 成人漫画全彩无遮挡| 欧美人与善性xxx| 久久午夜综合久久蜜桃| 日韩视频在线欧美| 亚洲国产av影院在线观看| 国语对白做爰xxxⅹ性视频网站| 国产成人av激情在线播放| 极品少妇高潮喷水抽搐| 91aial.com中文字幕在线观看| 久久精品亚洲熟妇少妇任你| 99re6热这里在线精品视频| 免费在线观看黄色视频的| 少妇人妻精品综合一区二区| 大陆偷拍与自拍| 啦啦啦在线免费观看视频4| 欧美97在线视频| 最黄视频免费看| 丰满少妇做爰视频| 成年人免费黄色播放视频| 伦理电影免费视频| 免费黄频网站在线观看国产| 亚洲熟女精品中文字幕| 亚洲欧美中文字幕日韩二区| 伦理电影大哥的女人| 我的亚洲天堂| 亚洲成色77777| 高清欧美精品videossex| 亚洲七黄色美女视频| 丝袜脚勾引网站| 国产精品国产三级专区第一集| 另类精品久久| 精品一区二区三区四区五区乱码 | 乱人伦中国视频| 操美女的视频在线观看| 午夜影院在线不卡| 国产亚洲午夜精品一区二区久久| 宅男免费午夜| 黄色毛片三级朝国网站| 你懂的网址亚洲精品在线观看| 在线 av 中文字幕| 国产精品国产三级专区第一集| 亚洲成人一二三区av| 看非洲黑人一级黄片| 国产在线一区二区三区精| 中国三级夫妇交换| 久久国产精品男人的天堂亚洲| 亚洲国产中文字幕在线视频| 免费在线观看完整版高清| 一区二区三区四区激情视频| 国产av精品麻豆| 女性被躁到高潮视频| 老鸭窝网址在线观看| 建设人人有责人人尽责人人享有的| 制服诱惑二区| 女人高潮潮喷娇喘18禁视频| 国产男女内射视频| 在线精品无人区一区二区三| 精品国产国语对白av| 大话2 男鬼变身卡| 日韩中文字幕欧美一区二区 | 国产1区2区3区精品| 成人亚洲欧美一区二区av| 最近手机中文字幕大全| 色精品久久人妻99蜜桃| 国产欧美日韩一区二区三区在线| 一本—道久久a久久精品蜜桃钙片| 精品视频人人做人人爽| 亚洲一区二区三区欧美精品| 久久国产亚洲av麻豆专区| 亚洲av电影在线进入| 国产老妇伦熟女老妇高清| 欧美亚洲 丝袜 人妻 在线| 中文字幕制服av| 亚洲成国产人片在线观看| 在线观看免费视频网站a站| 久久久久久久久免费视频了| 久久久久精品人妻al黑| 久久精品国产综合久久久| 亚洲一码二码三码区别大吗| 亚洲精品第二区| 伊人亚洲综合成人网| 老司机影院成人| 国产在视频线精品| 考比视频在线观看| 日韩制服骚丝袜av| 免费看av在线观看网站| 电影成人av| 女性被躁到高潮视频| 夫妻性生交免费视频一级片| 啦啦啦在线免费观看视频4| 欧美老熟妇乱子伦牲交| 男女边吃奶边做爰视频| 欧美变态另类bdsm刘玥| 欧美中文综合在线视频| 精品亚洲成国产av| 久久免费观看电影| 亚洲欧美中文字幕日韩二区| 亚洲精品美女久久av网站| 欧美日韩亚洲国产一区二区在线观看 | 国精品久久久久久国模美| 亚洲欧美成人精品一区二区| 亚洲欧洲日产国产| 国产亚洲av片在线观看秒播厂| 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产a三级三级三级| 人人妻人人澡人人看| 伊人久久大香线蕉亚洲五| 国产国语露脸激情在线看| www.自偷自拍.com| 精品一区二区三区四区五区乱码 | 九色亚洲精品在线播放| 欧美在线一区亚洲| 曰老女人黄片| 国产精品 国内视频| 少妇 在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品第二区| 欧美亚洲 丝袜 人妻 在线| 久久精品久久久久久久性| 99久久精品国产亚洲精品| 国产精品久久久人人做人人爽| 国产1区2区3区精品| 纵有疾风起免费观看全集完整版| 97人妻天天添夜夜摸| 两性夫妻黄色片| 亚洲成人免费av在线播放| 国产精品一区二区在线不卡| 欧美日韩综合久久久久久| 久久久久国产精品人妻一区二区| 少妇人妻 视频| 人妻 亚洲 视频| 在线观看免费高清a一片| 欧美在线一区亚洲| 狠狠精品人妻久久久久久综合| 一边亲一边摸免费视频| 美女高潮到喷水免费观看| 曰老女人黄片| 视频区图区小说| 一区二区日韩欧美中文字幕| 国产1区2区3区精品| 美女大奶头黄色视频| 成年av动漫网址| svipshipincom国产片| 亚洲av日韩精品久久久久久密 | 老司机靠b影院| 亚洲综合色网址| 成年人免费黄色播放视频| 成年人午夜在线观看视频| 观看av在线不卡| 久久精品熟女亚洲av麻豆精品| 一区二区三区四区激情视频| 男的添女的下面高潮视频| 成年av动漫网址| www.熟女人妻精品国产| 99精国产麻豆久久婷婷| 亚洲精品在线美女| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日本中文国产一区发布| 9热在线视频观看99| 日韩不卡一区二区三区视频在线| 两个人免费观看高清视频| 丝袜美足系列| 婷婷色麻豆天堂久久| 男女午夜视频在线观看| 欧美日韩视频高清一区二区三区二| 久久人妻熟女aⅴ| 美女视频免费永久观看网站| 国产男女内射视频| 久久女婷五月综合色啪小说| 久久人人爽人人片av| 成人国语在线视频| 国产在线视频一区二区| 午夜激情av网站| 国产麻豆69| 亚洲,欧美精品.| 亚洲自偷自拍图片 自拍| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美黑人精品巨大| 欧美日韩国产mv在线观看视频| av在线观看视频网站免费| 波多野结衣av一区二区av| 久久精品久久精品一区二区三区| 伊人亚洲综合成人网| 国产乱人偷精品视频| 久久人人爽av亚洲精品天堂| 高清黄色对白视频在线免费看| 国产爽快片一区二区三区| 成人影院久久| 欧美xxⅹ黑人| 熟妇人妻不卡中文字幕| 久久影院123| 99国产综合亚洲精品| 国产精品熟女久久久久浪| 1024视频免费在线观看| 无限看片的www在线观看| 午夜精品国产一区二区电影| 成人黄色视频免费在线看| 天天躁夜夜躁狠狠躁躁| 亚洲av成人精品一二三区| 精品亚洲成国产av| 亚洲欧洲日产国产| 亚洲精品美女久久久久99蜜臀 | 啦啦啦在线免费观看视频4| 99国产综合亚洲精品| 九色亚洲精品在线播放| 国产一区二区三区av在线| 国产日韩欧美亚洲二区| 国产精品 欧美亚洲| 99re6热这里在线精品视频| 9191精品国产免费久久| 亚洲av欧美aⅴ国产| 在线免费观看不下载黄p国产| 人成视频在线观看免费观看| 女性被躁到高潮视频| 亚洲美女黄色视频免费看| 999久久久国产精品视频| 99精品久久久久人妻精品| 晚上一个人看的免费电影| 最新在线观看一区二区三区 | 久久精品国产a三级三级三级| 欧美少妇被猛烈插入视频| av.在线天堂| 日韩制服丝袜自拍偷拍| 人人妻人人添人人爽欧美一区卜|