• <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)航方法
    久久99热这里只有精品18| 熟女人妻精品中文字幕| 欧美中文日本在线观看视频| 日本黄色片子视频| 亚洲在线观看片| 国产色爽女视频免费观看| 亚洲欧美日韩高清在线视频| 日韩在线高清观看一区二区三区| 亚洲电影在线观看av| 神马国产精品三级电影在线观看| 亚洲精品成人久久久久久| 91麻豆精品激情在线观看国产| 长腿黑丝高跟| 亚洲七黄色美女视频| 亚洲精品国产av成人精品 | 成人二区视频| 日日摸夜夜添夜夜添小说| 国产精品不卡视频一区二区| 在线天堂最新版资源| 天天躁夜夜躁狠狠久久av| 少妇丰满av| 国产精品人妻久久久影院| or卡值多少钱| 免费观看的影片在线观看| 国产麻豆成人av免费视频| 老司机午夜福利在线观看视频| 精品人妻一区二区三区麻豆 | 国产高清三级在线| 又爽又黄a免费视频| 欧美色视频一区免费| 免费人成视频x8x8入口观看| 我的女老师完整版在线观看| 成年版毛片免费区| 99热精品在线国产| 天堂影院成人在线观看| 蜜桃亚洲精品一区二区三区| 亚洲四区av| 最后的刺客免费高清国语| av在线天堂中文字幕| 一级毛片电影观看 | 两性午夜刺激爽爽歪歪视频在线观看| 一区二区三区高清视频在线| 蜜臀久久99精品久久宅男| 国产高清激情床上av| 哪里可以看免费的av片| 免费人成视频x8x8入口观看| 免费在线观看成人毛片| 日韩欧美国产在线观看| 看免费成人av毛片| 日本-黄色视频高清免费观看| 欧美又色又爽又黄视频| 男女做爰动态图高潮gif福利片| 欧美精品国产亚洲| 日韩高清综合在线| 欧美色欧美亚洲另类二区| 乱码一卡2卡4卡精品| 国产精品不卡视频一区二区| 大香蕉久久网| 日本与韩国留学比较| 免费电影在线观看免费观看| 丝袜美腿在线中文| 99热全是精品| 国产aⅴ精品一区二区三区波| 中文资源天堂在线| 亚洲成av人片在线播放无| 精品一区二区三区视频在线观看免费| 亚洲av美国av| 久久精品夜色国产| 不卡视频在线观看欧美| 99久久精品热视频| 男插女下体视频免费在线播放| 中文字幕人妻熟人妻熟丝袜美| 变态另类成人亚洲欧美熟女| 黑人高潮一二区| 99久久精品国产国产毛片| 国产大屁股一区二区在线视频| 一级黄色大片毛片| АⅤ资源中文在线天堂| 在线看三级毛片| 国产精品不卡视频一区二区| 成人国产麻豆网| 天天躁日日操中文字幕| 国产久久久一区二区三区| 国语自产精品视频在线第100页| 亚洲美女视频黄频| 最近手机中文字幕大全| 国产伦精品一区二区三区四那| 欧美+亚洲+日韩+国产| or卡值多少钱| 老司机午夜福利在线观看视频| 国内精品一区二区在线观看| 色综合亚洲欧美另类图片| 欧美xxxx黑人xx丫x性爽| 国产av一区在线观看免费| 国产老妇女一区| 日日撸夜夜添| 午夜福利视频1000在线观看| 国产精品野战在线观看| eeuss影院久久| 少妇猛男粗大的猛烈进出视频 | 少妇高潮的动态图| 欧美日韩在线观看h| 综合色丁香网| 国产精品永久免费网站| 日韩人妻高清精品专区| 男女之事视频高清在线观看| 99久久成人亚洲精品观看| 一级黄色大片毛片| 国产成人福利小说| 亚洲成a人片在线一区二区| 国产亚洲精品av在线| 色在线成人网| 欧美xxxx黑人xx丫x性爽| 日本一本二区三区精品| 91久久精品电影网| 男女之事视频高清在线观看| 国产色爽女视频免费观看| 99久久久亚洲精品蜜臀av| 色综合站精品国产| 女的被弄到高潮叫床怎么办| 亚洲第一区二区三区不卡| 波多野结衣高清作品| 99在线视频只有这里精品首页| a级一级毛片免费在线观看| 日韩av不卡免费在线播放| 狂野欧美激情性xxxx在线观看| 久久99热6这里只有精品| 精品少妇黑人巨大在线播放 | 日韩人妻高清精品专区| 色播亚洲综合网| 老师上课跳d突然被开到最大视频| 午夜福利在线观看免费完整高清在 | 人人妻人人澡欧美一区二区| 久久天躁狠狠躁夜夜2o2o| 久久精品国产鲁丝片午夜精品| 亚洲中文字幕日韩| 日韩强制内射视频| 国产av一区在线观看免费| 永久网站在线| 国产精品女同一区二区软件| 欧洲精品卡2卡3卡4卡5卡区| 色哟哟哟哟哟哟| 久久精品国产清高在天天线| 少妇人妻精品综合一区二区 | 日日摸夜夜添夜夜添小说| av天堂中文字幕网| 日韩精品青青久久久久久| 日韩精品青青久久久久久| 久久精品综合一区二区三区| 老司机午夜福利在线观看视频| 国语自产精品视频在线第100页| 久久中文看片网| 在现免费观看毛片| 尾随美女入室| 99热网站在线观看| 在线免费观看不下载黄p国产| 亚洲国产日韩欧美精品在线观看| 亚洲欧美成人综合另类久久久 | 日本免费a在线| 色av中文字幕| 成人毛片a级毛片在线播放| 九九热线精品视视频播放| 午夜爱爱视频在线播放| 白带黄色成豆腐渣| 长腿黑丝高跟| 色av中文字幕| 日韩欧美 国产精品| 听说在线观看完整版免费高清| 如何舔出高潮| 99久久久亚洲精品蜜臀av| 美女xxoo啪啪120秒动态图| 成年av动漫网址| 国产成人精品久久久久久| 九九热线精品视视频播放| 久久草成人影院| 国产黄片美女视频| 久久久久久久午夜电影| 国产精品永久免费网站| 日日啪夜夜撸| 午夜老司机福利剧场| 日韩大尺度精品在线看网址| 久久国产乱子免费精品| 亚洲最大成人手机在线| 内射极品少妇av片p| 久久精品国产清高在天天线| 国产视频内射| 国产午夜福利久久久久久| 亚洲无线在线观看| 欧美另类亚洲清纯唯美| 99久久久亚洲精品蜜臀av| 在线观看美女被高潮喷水网站| 久久精品国产亚洲av天美| 欧美日韩一区二区视频在线观看视频在线 | 成人毛片a级毛片在线播放| 亚洲性夜色夜夜综合| 成年av动漫网址| 免费在线观看影片大全网站| 丝袜美腿在线中文| 久久久欧美国产精品| 成人三级黄色视频| 日韩亚洲欧美综合| 婷婷亚洲欧美| 国产精品爽爽va在线观看网站| 亚洲美女搞黄在线观看 | 久久99热6这里只有精品| 少妇高潮的动态图| 国产一区二区在线av高清观看| 特大巨黑吊av在线直播| 日韩av不卡免费在线播放| 一区二区三区四区激情视频 | 精品一区二区三区视频在线观看免费| 综合色丁香网| 免费一级毛片在线播放高清视频| 久久婷婷人人爽人人干人人爱| 亚洲成人久久爱视频| 久久久久国产网址| 亚洲国产欧洲综合997久久,| 高清毛片免费看| 悠悠久久av| 最近中文字幕高清免费大全6| 十八禁国产超污无遮挡网站| 美女大奶头视频| 男人和女人高潮做爰伦理| 亚洲,欧美,日韩| 亚州av有码| 干丝袜人妻中文字幕| 精品不卡国产一区二区三区| 精品免费久久久久久久清纯| 人人妻人人看人人澡| 国产伦精品一区二区三区视频9| 国产美女午夜福利| 中文亚洲av片在线观看爽| a级毛片a级免费在线| 99久久成人亚洲精品观看| 波多野结衣高清无吗| 成人午夜高清在线视频| 日韩在线高清观看一区二区三区| 亚洲国产欧洲综合997久久,| 九九爱精品视频在线观看| 国产高潮美女av| 少妇的逼好多水| 欧美中文日本在线观看视频| 亚洲激情五月婷婷啪啪| 久久久久国内视频| 老熟妇乱子伦视频在线观看| 中文字幕熟女人妻在线| 日韩大尺度精品在线看网址| 日韩制服骚丝袜av| 99久国产av精品| 亚洲成人久久性| 成人午夜高清在线视频| 国产午夜精品论理片| 国产一区亚洲一区在线观看| 内地一区二区视频在线| 国产精品无大码| 在现免费观看毛片| 精品不卡国产一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 日韩一区二区视频免费看| 最新中文字幕久久久久| 夜夜夜夜夜久久久久| 你懂的网址亚洲精品在线观看 | 免费不卡的大黄色大毛片视频在线观看 | 国产熟女欧美一区二区| 成人午夜高清在线视频| 久久久久免费精品人妻一区二区| 久久精品影院6| 老司机福利观看| 国产爱豆传媒在线观看| 日韩高清综合在线| 白带黄色成豆腐渣| 18禁在线播放成人免费| 中文字幕av在线有码专区| 国产私拍福利视频在线观看| 高清毛片免费看| 人人妻人人澡欧美一区二区| 国产美女午夜福利| 精品一区二区三区人妻视频| 日韩高清综合在线| 国产白丝娇喘喷水9色精品| 欧美性猛交黑人性爽| 国产综合懂色| 国产蜜桃级精品一区二区三区| 99热这里只有是精品在线观看| 亚洲第一区二区三区不卡| 女人十人毛片免费观看3o分钟| 亚洲丝袜综合中文字幕| 国产精品人妻久久久影院| 久久草成人影院| 久久久久久九九精品二区国产| 99九九线精品视频在线观看视频| 少妇熟女欧美另类| 久久精品国产亚洲av涩爱 | 亚洲精品在线观看二区| 久久久久国内视频| 人妻制服诱惑在线中文字幕| 久久久精品欧美日韩精品| 观看免费一级毛片| 99久久精品热视频| 欧美另类亚洲清纯唯美| 亚洲成人中文字幕在线播放| 午夜激情福利司机影院| 身体一侧抽搐| 俺也久久电影网| 久久久久国产精品人妻aⅴ院| 久久久精品94久久精品| 99久久精品一区二区三区| 亚洲综合色惰| 国产91av在线免费观看| 日本色播在线视频| 亚洲一区高清亚洲精品| 精品人妻偷拍中文字幕| 简卡轻食公司| 一进一出抽搐gif免费好疼| 美女内射精品一级片tv| 国产探花极品一区二区| 亚洲性夜色夜夜综合| 一级毛片aaaaaa免费看小| 久久精品国产清高在天天线| 亚洲图色成人| 欧美日韩精品成人综合77777| 国产黄色视频一区二区在线观看 | 欧美绝顶高潮抽搐喷水| 九色成人免费人妻av| 天堂影院成人在线观看| 免费av毛片视频| 中文字幕免费在线视频6| 99国产极品粉嫩在线观看| 午夜爱爱视频在线播放| 22中文网久久字幕| 国产亚洲精品久久久com| 赤兔流量卡办理| 久久久久久久久大av| 国产91av在线免费观看| 亚洲欧美日韩无卡精品| 99在线视频只有这里精品首页| 夜夜看夜夜爽夜夜摸| av卡一久久| 特级一级黄色大片| a级毛片a级免费在线| 精品久久久噜噜| 亚洲av.av天堂| 欧美极品一区二区三区四区| 精华霜和精华液先用哪个| 99热6这里只有精品| 国产极品精品免费视频能看的| 丝袜美腿在线中文| 伦理电影大哥的女人| 偷拍熟女少妇极品色| 国产av麻豆久久久久久久| 国产一区亚洲一区在线观看| 精品久久久久久久久久免费视频| 日韩三级伦理在线观看| 真实男女啪啪啪动态图| 亚洲精品色激情综合| av黄色大香蕉| 久久久久久伊人网av| 精品久久久久久成人av| 五月伊人婷婷丁香| 国产精品久久久久久久久免| 亚洲国产精品sss在线观看| 欧美色欧美亚洲另类二区| 亚洲激情五月婷婷啪啪| 亚洲人与动物交配视频| 大型黄色视频在线免费观看| 美女 人体艺术 gogo| 22中文网久久字幕| 欧美xxxx性猛交bbbb| 日韩人妻高清精品专区| 最近的中文字幕免费完整| 国产免费一级a男人的天堂| 午夜精品在线福利| 两个人的视频大全免费| 成人av在线播放网站| av.在线天堂| 在线观看午夜福利视频| 女同久久另类99精品国产91| 最后的刺客免费高清国语| 亚洲av二区三区四区| 丝袜美腿在线中文| 亚洲经典国产精华液单| 成人无遮挡网站| 日韩,欧美,国产一区二区三区 | 哪里可以看免费的av片| 天堂av国产一区二区熟女人妻| 99在线人妻在线中文字幕| 男女那种视频在线观看| 日本精品一区二区三区蜜桃| 精品国内亚洲2022精品成人| 亚洲精品粉嫩美女一区| 国产在线男女| av在线天堂中文字幕| 97超碰精品成人国产| 露出奶头的视频| 成人精品一区二区免费| 黄色视频,在线免费观看| 波多野结衣巨乳人妻| 欧美一区二区国产精品久久精品| 精品国内亚洲2022精品成人| 国产精品日韩av在线免费观看| 麻豆一二三区av精品| 一级毛片电影观看 | 亚洲欧美日韩高清专用| 国产爱豆传媒在线观看| 国产日本99.免费观看| 日本五十路高清| 中文字幕精品亚洲无线码一区| 日韩av不卡免费在线播放| 一进一出抽搐动态| av中文乱码字幕在线| a级毛色黄片| 欧美又色又爽又黄视频| 国产伦精品一区二区三区视频9| 国产精品久久电影中文字幕| 免费av毛片视频| 中文字幕人妻熟人妻熟丝袜美| 国产精品国产三级国产av玫瑰| 又黄又爽又刺激的免费视频.| 国产私拍福利视频在线观看| 又黄又爽又刺激的免费视频.| 日本精品一区二区三区蜜桃| 国产欧美日韩精品一区二区| or卡值多少钱| 国产探花极品一区二区| 六月丁香七月| 亚洲欧美精品自产自拍| 亚洲成人久久爱视频| 亚洲av中文av极速乱| 国产视频内射| 男人舔女人下体高潮全视频| 级片在线观看| 变态另类丝袜制服| 在线播放无遮挡| 我要搜黄色片| 一个人看的www免费观看视频| 久久午夜福利片| 欧美日韩一区二区视频在线观看视频在线 | 亚洲自偷自拍三级| 亚洲三级黄色毛片| or卡值多少钱| 又黄又爽又免费观看的视频| 美女内射精品一级片tv| 日本黄大片高清| 午夜激情福利司机影院| 国产精品无大码| 中文字幕免费在线视频6| 亚洲av.av天堂| 精品一区二区三区av网在线观看| 99在线人妻在线中文字幕| 午夜日韩欧美国产| av视频在线观看入口| 伦理电影大哥的女人| 国产成人影院久久av| 成年女人看的毛片在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲精华国产精华液的使用体验 | 听说在线观看完整版免费高清| 欧美日韩国产亚洲二区| 日韩欧美免费精品| 免费av不卡在线播放| 国产亚洲av嫩草精品影院| ponron亚洲| 女人被狂操c到高潮| 国产不卡一卡二| 久久久久久久久久久丰满| 久久中文看片网| 青春草视频在线免费观看| 久久久久久久午夜电影| 国产精品99久久久久久久久| 女人十人毛片免费观看3o分钟| 三级男女做爰猛烈吃奶摸视频| 午夜免费激情av| 欧美日韩一区二区视频在线观看视频在线 | av在线天堂中文字幕| 天天一区二区日本电影三级| 日韩高清综合在线| 三级国产精品欧美在线观看| 韩国av在线不卡| 久久国产乱子免费精品| 中文字幕人妻熟人妻熟丝袜美| 欧美精品国产亚洲| 国产精品99久久久久久久久| 亚洲五月天丁香| 欧美三级亚洲精品| 三级毛片av免费| 精品人妻视频免费看| 色av中文字幕| 亚洲18禁久久av| 小蜜桃在线观看免费完整版高清| 校园人妻丝袜中文字幕| 欧美最新免费一区二区三区| 好男人在线观看高清免费视频| 国产男靠女视频免费网站| 女的被弄到高潮叫床怎么办| 黑人高潮一二区| 欧美日本视频| 日日摸夜夜添夜夜爱| 精品一区二区三区视频在线观看免费| 非洲黑人性xxxx精品又粗又长| 3wmmmm亚洲av在线观看| 淫妇啪啪啪对白视频| 变态另类丝袜制服| 国产精品乱码一区二三区的特点| 日韩高清综合在线| 91午夜精品亚洲一区二区三区| 如何舔出高潮| 精品午夜福利在线看| 欧美高清成人免费视频www| 色播亚洲综合网| 国产精品久久久久久久电影| 久久久久性生活片| 一a级毛片在线观看| 日本一本二区三区精品| 俄罗斯特黄特色一大片| av福利片在线观看| 在线免费十八禁| 国产三级在线视频| 又爽又黄a免费视频| 国产精品三级大全| 免费看光身美女| 99久久精品国产国产毛片| 嫩草影视91久久| 成年av动漫网址| 国产成人福利小说| 日本熟妇午夜| 久久久国产成人免费| 日本与韩国留学比较| av黄色大香蕉| 激情 狠狠 欧美| 国产白丝娇喘喷水9色精品| 在现免费观看毛片| 日本色播在线视频| 天堂网av新在线| 综合色丁香网| 亚洲专区国产一区二区| 午夜激情福利司机影院| 九九久久精品国产亚洲av麻豆| 大又大粗又爽又黄少妇毛片口| 偷拍熟女少妇极品色| 在线观看66精品国产| 午夜激情福利司机影院| 国产精品av视频在线免费观看| 啦啦啦韩国在线观看视频| 校园人妻丝袜中文字幕| 日本-黄色视频高清免费观看| 99热这里只有是精品在线观看| 亚洲综合色惰| 亚洲国产日韩欧美精品在线观看| 18禁在线无遮挡免费观看视频 | 蜜桃亚洲精品一区二区三区| 99热只有精品国产| 国产午夜精品论理片| 99久久精品国产国产毛片| 天天躁夜夜躁狠狠久久av| 欧美另类亚洲清纯唯美| 欧美一区二区国产精品久久精品| 亚洲经典国产精华液单| 一级黄片播放器| 日韩制服骚丝袜av| 97超碰精品成人国产| 一本久久中文字幕| 国产男靠女视频免费网站| 国产大屁股一区二区在线视频| 国产精品爽爽va在线观看网站| 亚洲精品一卡2卡三卡4卡5卡| 最近2019中文字幕mv第一页| 九九久久精品国产亚洲av麻豆| 午夜视频国产福利| 亚洲婷婷狠狠爱综合网| 国产男靠女视频免费网站| 中文字幕人妻熟人妻熟丝袜美| 99九九线精品视频在线观看视频| 熟女电影av网| 最新在线观看一区二区三区| 久久中文看片网| 精品久久久久久久久久免费视频| 国产高清有码在线观看视频| 三级国产精品欧美在线观看| 午夜精品国产一区二区电影 | 成人特级av手机在线观看| h日本视频在线播放| 日本一二三区视频观看| 日本色播在线视频| 国产av麻豆久久久久久久| 国产精品久久久久久av不卡| 午夜福利在线在线| 国产综合懂色| 亚洲av第一区精品v没综合| 1000部很黄的大片| 美女高潮的动态| 欧美精品国产亚洲| 日韩在线高清观看一区二区三区| 俺也久久电影网| 美女大奶头视频| 99热精品在线国产| 国产成人a∨麻豆精品| 蜜臀久久99精品久久宅男| 又黄又爽又刺激的免费视频.| 69人妻影院| 日韩人妻高清精品专区| 一本精品99久久精品77| 日韩欧美免费精品| 一进一出抽搐gif免费好疼| 久久精品国产亚洲av天美| 嫩草影院精品99| 免费观看在线日韩| 成年女人毛片免费观看观看9| АⅤ资源中文在线天堂| 网址你懂的国产日韩在线| 韩国av在线不卡| 波多野结衣高清作品| 日本熟妇午夜| 中文字幕精品亚洲无线码一区| 免费电影在线观看免费观看| 久久久久国产精品人妻aⅴ院|