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

    濾波辨識(3):方程誤差自回歸系統(tǒng)的濾波遞推廣義參數(shù)估計

    2022-12-30 15:52:18劉海波劉喜梅
    關(guān)鍵詞:新息廣義梯度

    丁 鋒,劉海波,劉喜梅

    (1.青島科技大學(xué) 自動化與電子工程學(xué)院,山東 青島 266061;2.江南大學(xué) 物聯(lián)網(wǎng)工程學(xué)院,江蘇 無錫 214122)

    梯度方法和牛頓方法是研究優(yōu)化問題的基本方法。用于系統(tǒng)辨識就得到梯度辨識方法和牛頓辨識方法。對于線性參數(shù)系統(tǒng),牛頓辨識方法退化為最小二乘辨識方法。

    梯度方法和牛頓方法與輔助模型辨識思想、多新息辨識理論、遞階辨識原理、耦合辨識概念等相結(jié)合[1-7],誕生出許多辨識方法,如輔助模型隨機梯度方法、輔助模型最小二乘方法、輔助模型牛頓辨識方法、多新息隨機梯度辨識方法、多新息最小二乘辨識方法、多新息牛頓辨識方法、遞階隨機梯度辨識方法、遞階最小二乘辨識方法、遞階牛頓辨識方法、耦合隨機梯度辨識方法、耦合最小二乘辨識方法、耦合牛頓辨識方法等。

    梯度方法、最小二乘方法、牛頓方法與濾波辨識理念相結(jié)合,誕生出濾波梯度辨識方法、濾波最小二乘辨識方法、濾波牛頓辨識方法等?!肚鄭u科技大學(xué)學(xué)報(自然科學(xué)版)》的連載論文研究了線性回歸系統(tǒng)、方程誤差系統(tǒng)和輸出誤差系統(tǒng)的遞階遞推辨識方法和遞階迭代辨識方法[8-16]。

    最近的連載論文研究了有限脈沖響應(yīng)滑動平均系統(tǒng) 的遞階遞推增廣參數(shù)辨識方法[9]和遞階增廣迭代參數(shù)辨識方法[10],以及濾波增廣遞推參數(shù)辨識方法[17]和濾波增廣迭代辨識方法[18]。本研究利用濾波辨識理念,研究方程誤差自回歸系統(tǒng),即受控自回歸自回歸系統(tǒng)的濾波遞推廣義辨識方法。

    基于數(shù)據(jù)濾波的辨識方法一般是針對有色噪聲干擾的系統(tǒng)而言。其基本思想是,在系統(tǒng)方程兩邊同乘以一個濾波器(通常取為噪聲模型傳遞函數(shù)的逆),對輸入輸出數(shù)據(jù)進行濾波,濾波后的系統(tǒng)是一個白噪聲干擾的系統(tǒng),而原系統(tǒng)的傳遞函數(shù)不變,然后使用濾波后的輸入輸出數(shù)據(jù)等進行辨識。因為濾波器(噪聲模型的傳遞函數(shù))是未知的,濾波后的輸入輸出數(shù)據(jù)也是未知的,實際中通常采用噪聲模型傳遞函數(shù)的估計進行濾波,所以基于數(shù)據(jù)濾波的辨識方法只能以遞推方式或迭代方式實現(xiàn)。

    如果知道系統(tǒng)先驗知識,濾波器也可以采用參數(shù)已知的固定模型,否則可采用動態(tài)變化模型,如即將討論的用噪聲傳遞函數(shù)的估計作為濾波器。

    辨識中數(shù)據(jù)濾波與通信和信號處理中的濾波有本質(zhì)不同。后者濾波是剔除被污染信號中的噪聲,使得噪聲在信號中的比重降低,如低通濾波器和高通濾波器。而濾波辨識方法只是改變系統(tǒng)模型的結(jié)構(gòu),不改變系統(tǒng)的輸入-輸出關(guān)系??傊?濾波辨識是將有色噪聲干擾的系統(tǒng)模型變換為白噪聲干擾的模型,即模型白色化(不是數(shù)據(jù)白色化),它能夠提高參數(shù)辨識精度。

    1 系統(tǒng)描述與濾波辨識模型

    方程誤差自回歸模型(EEAR模型)又稱為受控自回歸自回歸模型。考慮下列受控自回歸自回歸模型(controlled autoregressive autoregressive model,CARAR模型)描述的動態(tài)隨機系統(tǒng):

    其中{u(t)}和{y(t)}分別是系統(tǒng)的輸入和輸出序列,{v(t)}是零均值、方差為σ2的隨機白噪聲序列,A(z),B(z)和C(z)是單位后移算子z-1的多項式:

    設(shè)階次na,nb和nc已知。記n:=na+nb+nc。且設(shè)t≤0時,所有變量的初值均為零,如y(t)=0,u(t)=0,v(t)=0。

    置有關(guān)參數(shù)向量和信息向量如下:

    定義中間變量:

    取濾波器L(z):=C(z)。定義濾波輸入uf(t),濾波輸出yf(t)分別為

    定義濾波輸出信息向量φaf(t)和濾波輸入信息向量φbf(t)如下:

    式(1)兩邊乘以L(z)得到

    這個濾波后的模型是一個方程誤差模型(即CAR模型)。將它寫為向量形式為

    將式(6)代入式(10)得到

    其中

    在辨識模型(11)中,參數(shù)向量θ包含了系統(tǒng)的所有參數(shù)ai,bi和ci。目標(biāo)就是利用系統(tǒng)的輸入輸出數(shù)據(jù)u(t)和y(t),研究和提出估計參數(shù)向量θ的辨識方法。

    對于辨識模型(11),由于多項式C(z)是未知的(即L(z)是未知的),故uf(t),yf(t),以及φaf(t)和φbf(t)都是未知的,且由它們構(gòu)成的信息向量φ(t)也是未知的,故需要采用估計值代替未知變量來推導(dǎo)基于數(shù)據(jù)濾波的辨識方法。

    2 濾波廣義隨機梯度辨識算法

    對于EEAR系統(tǒng)的辨識模型(11),定義梯度準(zhǔn)則函數(shù):

    這些遞推關(guān)系式無法實現(xiàn),因為式(13)~(14)右邊信息向量φ(t)是未知的,所以需要先構(gòu)造它們的估計。由噪聲模型參數(shù)估計

    構(gòu)造多項式C(z)的估計:

    根據(jù)定義式(7)~(8)的結(jié)構(gòu),用yf(t-i)和uf(ti)的估計和構(gòu)造濾波輸出信息向量φaf(t)和濾波輸入信息向量φbf(t)的估計:

    將式(13)~(14)中右邊未知信息向量φ(t)用其估計代替,得到式(20)~(22),聯(lián)立式(19),(17)~(18),(2)~(3),(15)~(16),可以得到估計EEAR系統(tǒng)參數(shù)向量θ的基于濾波的廣義隨機梯度辨識算法(filtering-based generalized stochastic gradient identification algorithm),簡稱為濾波廣義隨機梯度算法(filtered generalized stochastic gradient algorithm,F-GSG算法):

    F-GSG算法(20)~(30)計算參數(shù)估計向量^θ(t)的步驟如下。

    1)初始化:令t=1。置初值,i)=1/p0,i=1,2,…,max[na,nb,nc],p0是一個大正數(shù),如p0=106。給定數(shù)據(jù)長度Le。

    2)采集輸入輸出數(shù)據(jù)u(t)和y(t),用式(24)~(27)構(gòu)造信息向量和φc(t),用式(23)構(gòu)造信息向量。

    3)用式(22)計算r(t),用式(21)計算新息e(t)。

    6)如果t<Le,t就增加1,轉(zhuǎn)到步驟2);否則輸出參數(shù)估計,終止遞推計算過程。

    在F-GSG辨識算法(20)~(30)中,e(t)=為辨識新息。定義辨識殘差。

    3 濾波多新息廣義隨機梯度算法

    為改進F-GSG算法的收斂速度,一種方法是借助于多新息辨識理論,通過擴展新息維數(shù),推導(dǎo)出基于濾波的多新息廣義隨機梯度辨識算法。

    設(shè)正整數(shù)p表示新息長度?;贔-GSG算法(20)~(30),將系統(tǒng)輸出y(t)和信息向量擴展為堆積輸出向量Y(p,t)和堆積信息矩陣(p,t):

    將式(20)中標(biāo)量新息e(t)=y(tǒng)(t)-^φT(t)^θ(t-1)∈?擴展為新息向量

    這是新息長度p=1的基于濾波的多新息廣義隨機梯度算法。將上式中和E(1,t)中的“1”換為p,得到式(34),聯(lián)立式(31)~(33)和(22)~(30),就得到新息長度為p的估計EEAR系統(tǒng)參數(shù)向量θ的基于濾波的多新息廣義隨機梯度辨識算法(filtering-based multi-innovation generalized stochastic gradient identification algorithm),簡稱為濾波多新息廣義隨機梯度算法(filtered multi-innovation generalized stochastic gradient algorithm,FMI-GSG算法):

    1)初始化:令t=1,給定新息長度p。置初值,p0是一個大正數(shù),如p0=106。給定數(shù)據(jù)長度Le。

    2)采集輸入輸出數(shù)據(jù)u(t)和y(t),用式(40)~(43)構(gòu)造信息向量和φc(t),用式(39)構(gòu)造信息向量。

    3)用式(37)和(38)構(gòu)造堆積輸出向量Y(p,t)和堆積信息矩陣。

    4)用式(35)計算新息向量E(p,t),用式(36)計算r(t)。

    6)如果t<Le,t就增加1,轉(zhuǎn)到步驟2);否則輸出參數(shù)估計,終止遞推計算過程。

    引理1對于F-MI-GSG辨識算法(34)~(46),新息向量E(p,t)與殘差向量

    滿足關(guān)系:

    4 濾波多新息廣義投影辨識方法

    下面利用最速下降法推導(dǎo)濾波多新息廣義投影辨識方法??紤]從j=t-p+1到j(luò)=t的數(shù)據(jù)窗里的p組數(shù)據(jù)。定義堆積輸出向量Y(p,t)和堆積信息矩陣Φ(p,t)如下:

    基于EEAR系統(tǒng)的濾波辨識模型(11),定義滑動數(shù)據(jù)窗準(zhǔn)則函數(shù):

    假設(shè)步長(step size)為μ(t)。使用最速下降法(負(fù)梯度搜索),極小化準(zhǔn)則函數(shù)J2(θ),可以得到梯度遞推關(guān)系:

    其中

    下面求最佳步長(best step-size)μ(t)。將代入準(zhǔn)則函數(shù)J2(θ)中,可得

    由此可求得最佳步長為

    由于堆積信息矩陣Φ(p,t)是未知的,所以算法(48)~(51)無法實現(xiàn)。解決的方法是用φ(t-i)的估計構(gòu)造堆積信息矩陣Φ(p,t)的估計

    式(49)~(51)中未知Φ(p,t)用其估計代替,得到式(53)~(55),聯(lián)立式(48),(52)和(39)~(46),就得到辨識EEAR系統(tǒng)參數(shù)向量θ的基于濾波的多新息廣義投影辨識算法(filtering-based multi-innovation generalized projection identification algorithm),簡稱為濾波多新息廣義投影算法(filtered multi-innovation generalized projection algorithm,F-MI-GProj算法):

    因為收斂因子的計算式比較復(fù)雜,故對其進行簡化。由于對于任意實向量x和非負(fù)定對稱矩陣Q,不等式成立,其中為矩陣Q的最大特征值,所以

    于是,收斂因子可以保守取為

    因為計算矩陣的跡(trace)比計算特征值簡單,所以收斂因子可以更保守取為

    因此,將式(53)修改為

    則式(54)和(56)~(67)構(gòu)成了簡化的濾波多新息廣義投影算法(F-MI-GProj算法)。

    5 濾波遞推廣義梯度辨識算法

    定義堆積輸出向量Y(t)和堆積輸入信息矩陣Φ(t)如下:

    根據(jù)式(11),定義二次準(zhǔn)則函數(shù):

    定義遞推關(guān)系:

    遞推關(guān)系式(68)~(71)不可實現(xiàn),因為右邊包含了未知信息向量φ(t)。解決方案是用其估計代替,得到式(72)~(75),聯(lián)立式(23)~(30),便得到辨識EEAR系統(tǒng)參數(shù)向量θ的基于濾波的遞推廣義梯度辨識算法(filtering-based recursive generalized gradient identification algorithm),簡稱為濾波遞推廣義梯度算法(filtered recursive generalized gradient algorithm,F-RGG算法):

    1)初始化:令t=1,置初值,,…,max[na,nb,nc],p0是一個大正數(shù),如p0=106。給定參數(shù)估計精度ε。

    2)采集輸入輸出數(shù)據(jù)u(t)和y(t),用式(77)~(80)構(gòu)造信息向量和φc(t),用式(76)構(gòu)造信息向量。

    3)用式(73)計算r(t),用式(74)計算向量ξ(t),用式(75)計算矩陣R(t)。用式(72)刷新參數(shù)估計向量)。

    6 濾波多新息遞推廣義梯度算法

    設(shè)正整數(shù)p表示新息長度。定義堆積輸出向量Y(p,t)和堆積輸入信息矩陣Φ(p,t)如下:

    基于F-RGG算法(72)~(83),將式(73)~(75)中系統(tǒng)輸出y(t)和輸入信息向量擴展為堆積輸出向量Y(p,t)和堆積信息矩陣),得到式(87)~(89),聯(lián)立式(72),(84)~(85)和(76)~(83),便得到辨識EEAR系統(tǒng)參數(shù)向量θ的基于濾波的多新息遞推廣義梯度辨識算法(filtering-based multi-innovation recursive generalized gradient identification algorithm),簡稱為濾波多新息遞推廣義梯度算法(filtered multi-innovation recursive generalized gradient algorithm,F-MI-RGG算法):

    1)初始化:令t=1, 給定新息長度p和精度指標(biāo)ε。置初值,,p0是一個大正數(shù),如p0=106。

    2)采集輸入輸出數(shù)據(jù)u(t)和y(t),用式(93)~(96)構(gòu)造信息向量和φc(t),用式(92)構(gòu)造信息向量。

    3)用式(90)和(91)構(gòu)造堆積輸出向量Y(p,t)和堆積信息矩陣。

    4)用式(87)計算r(t),用式(88)計算向量ξ(t),用式(89)計算矩陣R(t)。用式(86)刷新參數(shù)估計向量。

    7 濾波遞推廣義最小二乘算法

    參照遞推最小二乘算法的推導(dǎo),極小化J3(θ),可以得到下列最小二乘遞推關(guān)系:

    這些遞推關(guān)系無法實現(xiàn),因為式(100)~(102)右邊的信息向量φ(t)是未知的,在辨識算法中使用它的估計代替,得到式(103)~(106),聯(lián)立式(23)~(30),便得到辨識EEAR系統(tǒng)參數(shù)向量θ的基于濾波的遞推廣義最小二乘辨識算法(filteringbased recursive generalized least squares identification algorithm),簡稱為濾波遞推廣義最小二乘算法(filtered recursive generalized least squares algorithm,F-RGLS算法):

    1)初始化:令t=1。置初值,,,p0是一個大正數(shù),如p0=106。給定參數(shù)估計精度ε。

    2)采集輸入輸出數(shù)據(jù)u(t)和y(t),用式(108)~(111)構(gòu)造信息向量和φc(t),用式(107)構(gòu)造信息向量。

    3)用式(104)計算新息e(t),用式(105)和(106)計算增益向量L(t)和協(xié)方差陣P(t)。根據(jù)式(103)刷新參數(shù)估計向量。

    8 濾波多新息廣義最小二乘算法

    基于F-RGLS算法(103)~(114),將系統(tǒng)輸出y(t)和信息向量分別擴展為堆積輸出向量Y(p,t)和堆積信息矩陣:

    得到式(118)~(121),聯(lián)立式(115)~(117)和(107)~(114),便得到辨識EEAR系統(tǒng)參數(shù)向量θ的基于濾波的多新息遞推廣義最小二乘辨識算法(filtering-based multi-innovation recursive generalized least squares identification algorithm),簡稱為濾波多新息廣義最小二乘算法(filtered multi-innovation generalized least squares algorithm,F-MIGLS算法):

    F-MI-GLS算法(118)~(131)計算參數(shù)估計向量的步驟如下。

    1)初始化:令t=1,給定新息長度p和精度指標(biāo)ε。置初值,…,max[na,nb,nc],p0是一個大正數(shù),如p0=106。

    2)采集輸入輸出數(shù)據(jù)u(t)和y(t),用式(125)~(128)構(gòu)造信息向量和φc(t),用式(124)構(gòu)造信息向量。

    3)用式(122)和(123)構(gòu)造堆積輸出向量Y(p,t)和堆積信息矩陣。

    4)用式(119)計算新息向量E(p,t),用式(120)和(121)計算增益矩陣L(t)和協(xié)方差陣P(t)。根據(jù)式(118)刷新參數(shù)估計向量。

    多新息辨識方法包含多新息遞推辨識方法和多新息迭代辨識方法。對于多新息遞推辨識方法,常常省略遞推兩個字,將多新息遞推最小二乘辨識方法簡稱為多新息最小二乘辨識方法。如,濾波多新息遞推廣義最小二乘辨識算法簡稱為濾波多新息廣義最小二乘算法,濾波多新息遞推廣義最小二乘辨識算法簡稱為濾波多新息廣義最小二乘算法。

    9 結(jié) 語

    針對方程誤差自回歸(EEAR)系統(tǒng),研究和提出了濾波(多新息)廣義隨機梯度辨識方法、濾波多新息廣義投影辨識方法、濾波(多新息)遞推廣義梯度辨識方法、濾波(多新息)廣義最小二乘辨識方法。盡管這些濾波廣義辨識方法是針對自回歸噪聲干擾下的方程誤差自回歸隨機系統(tǒng)提出的,但是其思想可以推廣到有色噪聲干擾下的線性和非線性多變量隨機系統(tǒng)中[19-31]。

    猜你喜歡
    新息廣義梯度
    Rn中的廣義逆Bonnesen型不等式
    一個改進的WYL型三項共軛梯度法
    傳遞函數(shù)辨識(21):線性回歸系統(tǒng)的遞階遞推參數(shù)估計
    一種自適應(yīng)Dai-Liao共軛梯度法
    從廣義心腎不交論治慢性心力衰竭
    一類扭積形式的梯度近Ricci孤立子
    M估計的強跟蹤SVD-UKF算法在組合導(dǎo)航中的應(yīng)用
    電子科技(2018年7期)2018-07-23 05:30:32
    有限群的廣義交換度
    自適應(yīng)卡爾曼濾波在航空重力異常解算的應(yīng)用研究
    基于新息正交性自適應(yīng)濾波的慣性/地磁組合導(dǎo)航方法
    亚洲欧美一区二区三区国产| 人人妻人人澡人人爽人人夜夜| 国产黄色视频一区二区在线观看| 99久久精品热视频| 成年人午夜在线观看视频| 亚洲av男天堂| 国产探花极品一区二区| 国产黄片视频在线免费观看| 永久网站在线| 青春草国产在线视频| 少妇裸体淫交视频免费看高清| 精品少妇黑人巨大在线播放| 国产片特级美女逼逼视频| 欧美日韩视频精品一区| 国产精品麻豆人妻色哟哟久久| 久久人人爽av亚洲精品天堂 | 成人特级av手机在线观看| 久久久久久久午夜电影| 日韩,欧美,国产一区二区三区| 韩国av在线不卡| 日本猛色少妇xxxxx猛交久久| 黄片无遮挡物在线观看| 国产爱豆传媒在线观看| 国产伦理片在线播放av一区| 国产在线一区二区三区精| 日韩欧美精品免费久久| 久久久久久久久大av| 亚洲av欧美aⅴ国产| 国产精品不卡视频一区二区| 少妇 在线观看| 丰满乱子伦码专区| 亚洲欧美精品专区久久| 80岁老熟妇乱子伦牲交| 极品教师在线视频| 激情五月婷婷亚洲| 成人黄色视频免费在线看| 97超碰精品成人国产| 欧美日韩视频高清一区二区三区二| 亚洲一区二区三区欧美精品 | 国产精品成人在线| 天堂网av新在线| 国产欧美另类精品又又久久亚洲欧美| 亚洲激情五月婷婷啪啪| 18+在线观看网站| av线在线观看网站| 99久久精品国产国产毛片| 亚洲av男天堂| 久久99蜜桃精品久久| 在线看a的网站| 日韩成人伦理影院| 亚洲精品456在线播放app| 午夜精品国产一区二区电影 | 国内揄拍国产精品人妻在线| 91精品一卡2卡3卡4卡| av女优亚洲男人天堂| 国产精品女同一区二区软件| 97超视频在线观看视频| 国产一区二区三区av在线| 亚洲欧美一区二区三区黑人 | 在线 av 中文字幕| 高清av免费在线| 哪个播放器可以免费观看大片| 国产男女超爽视频在线观看| 欧美97在线视频| 美女视频免费永久观看网站| 最后的刺客免费高清国语| 国产av不卡久久| 亚洲成人精品中文字幕电影| 精品酒店卫生间| 亚洲三级黄色毛片| 欧美极品一区二区三区四区| 国产黄色视频一区二区在线观看| 91在线精品国自产拍蜜月| 国产精品av视频在线免费观看| 亚洲精品日韩在线中文字幕| 亚洲成人av在线免费| 美女被艹到高潮喷水动态| 亚洲精品中文字幕在线视频 | 欧美成人一区二区免费高清观看| 国产视频首页在线观看| 日本色播在线视频| 免费观看av网站的网址| 看黄色毛片网站| 亚洲天堂国产精品一区在线| 成人鲁丝片一二三区免费| 禁无遮挡网站| 精品少妇久久久久久888优播| av免费在线看不卡| 国产精品无大码| 国产成人精品久久久久久| 日韩,欧美,国产一区二区三区| 成人综合一区亚洲| 亚洲精品456在线播放app| 久久久色成人| 国产综合精华液| 亚洲av.av天堂| 午夜福利在线在线| 亚洲电影在线观看av| 一个人看视频在线观看www免费| 精品亚洲乱码少妇综合久久| a级毛色黄片| 国产老妇伦熟女老妇高清| 欧美xxxx黑人xx丫x性爽| 超碰av人人做人人爽久久| 七月丁香在线播放| 美女主播在线视频| 永久免费av网站大全| 亚洲成色77777| 人人妻人人爽人人添夜夜欢视频 | 精品人妻一区二区三区麻豆| 亚洲国产色片| 亚洲综合精品二区| 国产精品久久久久久av不卡| 超碰av人人做人人爽久久| 亚洲自拍偷在线| 国产精品.久久久| 欧美97在线视频| 内地一区二区视频在线| 高清毛片免费看| 国产成人精品一,二区| 美女视频免费永久观看网站| 久久99精品国语久久久| 插逼视频在线观看| 久久久久久久久久人人人人人人| 国产精品.久久久| 在线天堂最新版资源| 在线天堂最新版资源| 亚洲经典国产精华液单| 夫妻性生交免费视频一级片| 国产白丝娇喘喷水9色精品| .国产精品久久| 国产毛片在线视频| 欧美高清成人免费视频www| 夫妻性生交免费视频一级片| 成人黄色视频免费在线看| 精品国产三级普通话版| 国产一区二区在线观看日韩| 免费观看a级毛片全部| 深夜a级毛片| 国产精品成人在线| 免费高清在线观看视频在线观看| 最新中文字幕久久久久| 国产精品秋霞免费鲁丝片| 亚洲精品自拍成人| 天天躁夜夜躁狠狠久久av| 精品国产露脸久久av麻豆| 亚洲电影在线观看av| 国产av码专区亚洲av| 夜夜看夜夜爽夜夜摸| 欧美精品人与动牲交sv欧美| 婷婷色麻豆天堂久久| 日韩av不卡免费在线播放| 久久久久久九九精品二区国产| 午夜福利在线在线| 中文欧美无线码| 久久99热这里只有精品18| 免费不卡的大黄色大毛片视频在线观看| 激情 狠狠 欧美| 久久6这里有精品| 嫩草影院入口| 精品久久久精品久久久| 日本一二三区视频观看| 久久久精品欧美日韩精品| 高清欧美精品videossex| 久久久久网色| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩一本色道免费dvd| 精品一区二区三区视频在线| 综合色av麻豆| 99久久中文字幕三级久久日本| 别揉我奶头 嗯啊视频| 国产亚洲精品久久久com| 18禁在线播放成人免费| 免费av毛片视频| 国产欧美日韩精品一区二区| 免费在线观看成人毛片| 亚洲av一区综合| 高清日韩中文字幕在线| 国产黄色免费在线视频| 午夜免费观看性视频| 亚洲va在线va天堂va国产| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲国产精品成人综合色| 久久99蜜桃精品久久| 搞女人的毛片| 久久久久久久亚洲中文字幕| 秋霞在线观看毛片| av国产久精品久网站免费入址| 亚洲国产欧美在线一区| av免费观看日本| 天天躁夜夜躁狠狠久久av| 国产熟女欧美一区二区| 青春草视频在线免费观看| 中文字幕制服av| av女优亚洲男人天堂| 亚洲va在线va天堂va国产| 在线观看免费高清a一片| 麻豆国产97在线/欧美| 美女主播在线视频| 99精国产麻豆久久婷婷| 欧美激情国产日韩精品一区| 老女人水多毛片| 99热这里只有精品一区| 九九爱精品视频在线观看| 成人亚洲精品av一区二区| 国产精品伦人一区二区| av黄色大香蕉| .国产精品久久| 成年人午夜在线观看视频| 少妇猛男粗大的猛烈进出视频 | 久久久久久国产a免费观看| 国产精品成人在线| 欧美精品人与动牲交sv欧美| 国产av码专区亚洲av| 亚洲经典国产精华液单| 亚洲国产日韩一区二区| 亚洲性久久影院| 午夜免费观看性视频| 直男gayav资源| 免费在线观看成人毛片| av卡一久久| 日本av手机在线免费观看| 男女国产视频网站| 国产精品无大码| 婷婷色麻豆天堂久久| 国产一区二区在线观看日韩| 国产午夜精品久久久久久一区二区三区| 男女啪啪激烈高潮av片| 亚洲欧美日韩无卡精品| 亚洲成人精品中文字幕电影| 亚洲激情五月婷婷啪啪| 日韩欧美 国产精品| 国产欧美亚洲国产| 久久久久久久久久成人| 国产成人午夜福利电影在线观看| 天堂网av新在线| 亚洲天堂国产精品一区在线| 蜜桃亚洲精品一区二区三区| 最近的中文字幕免费完整| 大香蕉97超碰在线| 美女内射精品一级片tv| 久久综合国产亚洲精品| 日韩强制内射视频| 97超视频在线观看视频| 又粗又硬又长又爽又黄的视频| 亚洲伊人久久精品综合| 国产精品一及| 国产精品秋霞免费鲁丝片| 97人妻精品一区二区三区麻豆| 亚洲成人av在线免费| 免费看光身美女| 搞女人的毛片| 精品久久久久久久人妻蜜臀av| 新久久久久国产一级毛片| 伊人久久国产一区二区| 国产精品伦人一区二区| 一区二区三区乱码不卡18| 狂野欧美激情性xxxx在线观看| 在线观看一区二区三区| 两个人的视频大全免费| 天堂中文最新版在线下载 | 一本一本综合久久| 中文天堂在线官网| 欧美+日韩+精品| 99九九线精品视频在线观看视频| 色哟哟·www| 亚洲国产av新网站| 国产成人91sexporn| 午夜福利高清视频| 成人高潮视频无遮挡免费网站| 国产伦理片在线播放av一区| 欧美变态另类bdsm刘玥| av线在线观看网站| 极品教师在线视频| 国产 精品1| 精品一区在线观看国产| 亚洲国产欧美在线一区| 亚洲自偷自拍三级| 亚洲精品乱久久久久久| 网址你懂的国产日韩在线| 亚洲在久久综合| 国产成人freesex在线| 人妻一区二区av| 真实男女啪啪啪动态图| 久久精品国产亚洲网站| 最近2019中文字幕mv第一页| 天堂俺去俺来也www色官网| 午夜福利视频1000在线观看| 日韩av免费高清视频| www.色视频.com| 日韩欧美 国产精品| 我的老师免费观看完整版| 97超碰精品成人国产| 国产午夜精品一二区理论片| a级一级毛片免费在线观看| 国产精品不卡视频一区二区| 97超碰精品成人国产| 看黄色毛片网站| 搞女人的毛片| 亚洲av在线观看美女高潮| 中文在线观看免费www的网站| 日本三级黄在线观看| 成年av动漫网址| 九九爱精品视频在线观看| 久久这里有精品视频免费| 一级爰片在线观看| 美女被艹到高潮喷水动态| 可以在线观看毛片的网站| 大又大粗又爽又黄少妇毛片口| 亚洲av中文字字幕乱码综合| 亚洲欧美成人精品一区二区| 熟女av电影| 成年女人在线观看亚洲视频 | 99热6这里只有精品| 欧美xxxx性猛交bbbb| 青春草国产在线视频| 欧美 日韩 精品 国产| 99热这里只有是精品50| 国产毛片a区久久久久| 国产熟女欧美一区二区| 亚洲成色77777| 久久久久性生活片| 国产精品熟女久久久久浪| 中文精品一卡2卡3卡4更新| 欧美精品国产亚洲| 2022亚洲国产成人精品| 午夜福利在线在线| 亚洲精品乱码久久久久久按摩| 亚洲精品一二三| 波多野结衣巨乳人妻| 女人久久www免费人成看片| 特级一级黄色大片| 亚洲成人av在线免费| 亚洲最大成人中文| 国产精品人妻久久久久久| 成人一区二区视频在线观看| 日韩视频在线欧美| 可以在线观看毛片的网站| 精品人妻一区二区三区麻豆| 久久久a久久爽久久v久久| 成人综合一区亚洲| 欧美日本视频| 国产欧美亚洲国产| 亚洲欧美日韩东京热| 国产高清不卡午夜福利| 亚洲av中文字字幕乱码综合| 国产精品人妻久久久久久| 天天一区二区日本电影三级| 欧美最新免费一区二区三区| 日韩欧美 国产精品| 亚洲在久久综合| 大话2 男鬼变身卡| av网站免费在线观看视频| 亚洲欧美精品专区久久| 免费观看在线日韩| 交换朋友夫妻互换小说| 日韩中字成人| 高清日韩中文字幕在线| 国产一区二区亚洲精品在线观看| av在线天堂中文字幕| 午夜福利视频精品| 2021少妇久久久久久久久久久| 国产亚洲av嫩草精品影院| 国产精品99久久99久久久不卡 | 我的女老师完整版在线观看| 777米奇影视久久| 久久久精品94久久精品| 亚洲内射少妇av| 舔av片在线| 好男人视频免费观看在线| 亚洲av免费高清在线观看| 特大巨黑吊av在线直播| 亚洲天堂国产精品一区在线| 丝瓜视频免费看黄片| av专区在线播放| 秋霞伦理黄片| 国产成人freesex在线| 国产精品一区二区三区四区免费观看| 午夜精品国产一区二区电影 | 欧美xxxx黑人xx丫x性爽| 日韩亚洲欧美综合| 成年免费大片在线观看| 成年av动漫网址| 18禁动态无遮挡网站| 国产探花在线观看一区二区| 国产欧美亚洲国产| 大话2 男鬼变身卡| 欧美xxⅹ黑人| 日本午夜av视频| 中国国产av一级| 一区二区三区乱码不卡18| 亚洲国产色片| 国产 一区 欧美 日韩| 人人妻人人看人人澡| 成人毛片a级毛片在线播放| 3wmmmm亚洲av在线观看| 国产亚洲午夜精品一区二区久久 | 免费观看无遮挡的男女| 亚洲四区av| 久久久久久久国产电影| 亚洲丝袜综合中文字幕| 免费看光身美女| 91久久精品电影网| 亚洲欧美成人综合另类久久久| 少妇被粗大猛烈的视频| 国产一区亚洲一区在线观看| 亚洲美女视频黄频| 国产视频首页在线观看| 国语对白做爰xxxⅹ性视频网站| 精品少妇久久久久久888优播| 亚洲综合精品二区| 亚洲美女搞黄在线观看| 特大巨黑吊av在线直播| 亚洲国产精品国产精品| 欧美激情久久久久久爽电影| 免费观看av网站的网址| 男女下面进入的视频免费午夜| 亚洲精品色激情综合| 激情 狠狠 欧美| 日韩 亚洲 欧美在线| 久久久久久久大尺度免费视频| 成人欧美大片| 一个人看的www免费观看视频| 国产永久视频网站| 99久久精品热视频| 午夜免费鲁丝| 天堂俺去俺来也www色官网| 国产大屁股一区二区在线视频| 熟女av电影| 麻豆久久精品国产亚洲av| 国产一区二区在线观看日韩| 欧美一区二区亚洲| 一级毛片我不卡| av国产久精品久网站免费入址| 亚洲精品456在线播放app| 在线亚洲精品国产二区图片欧美 | 国内精品美女久久久久久| 国产 一区精品| 欧美另类一区| 午夜爱爱视频在线播放| 日日啪夜夜爽| 国产欧美日韩精品一区二区| 免费观看在线日韩| 亚洲欧美成人综合另类久久久| 国产黄a三级三级三级人| 免费大片黄手机在线观看| 97超碰精品成人国产| 黄色日韩在线| 最近最新中文字幕免费大全7| 亚洲国产精品成人久久小说| 免费av观看视频| 久久99蜜桃精品久久| 成人漫画全彩无遮挡| 久久99热6这里只有精品| 80岁老熟妇乱子伦牲交| 91午夜精品亚洲一区二区三区| 久久ye,这里只有精品| 蜜桃亚洲精品一区二区三区| 老司机影院毛片| 丰满人妻一区二区三区视频av| 日日摸夜夜添夜夜爱| 蜜桃久久精品国产亚洲av| 色网站视频免费| 一边亲一边摸免费视频| 国产乱来视频区| 搞女人的毛片| av免费观看日本| 免费观看性生交大片5| 亚洲精品中文字幕在线视频 | 日本色播在线视频| 日韩三级伦理在线观看| 卡戴珊不雅视频在线播放| 国产视频内射| 国产精品人妻久久久久久| 白带黄色成豆腐渣| 欧美bdsm另类| 久久97久久精品| 免费黄网站久久成人精品| 国产日韩欧美在线精品| 免费观看在线日韩| 亚洲av.av天堂| 国产成年人精品一区二区| av在线播放精品| 亚洲久久久久久中文字幕| 天天一区二区日本电影三级| 日韩欧美精品v在线| 97超视频在线观看视频| 一二三四中文在线观看免费高清| 99热这里只有精品一区| 人人妻人人看人人澡| 黄色一级大片看看| 亚洲一区二区三区欧美精品 | 日本猛色少妇xxxxx猛交久久| 欧美性感艳星| 午夜爱爱视频在线播放| 国产伦理片在线播放av一区| 免费人成在线观看视频色| 在线a可以看的网站| 亚洲综合色惰| 国产成人免费观看mmmm| 天美传媒精品一区二区| 一区二区av电影网| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 91精品伊人久久大香线蕉| 亚洲精品国产av蜜桃| 欧美三级亚洲精品| av免费在线看不卡| 丰满人妻一区二区三区视频av| 色播亚洲综合网| 国产一级毛片在线| 97在线人人人人妻| 国产一区二区三区综合在线观看 | 观看美女的网站| 日韩欧美精品免费久久| 在线观看一区二区三区激情| 九九爱精品视频在线观看| 人妻 亚洲 视频| 欧美激情国产日韩精品一区| 在线播放无遮挡| 国产有黄有色有爽视频| 哪个播放器可以免费观看大片| 国产高清国产精品国产三级 | 男插女下体视频免费在线播放| videos熟女内射| 久久精品国产亚洲av天美| 人妻 亚洲 视频| 日产精品乱码卡一卡2卡三| a级毛片免费高清观看在线播放| 高清毛片免费看| av免费观看日本| 麻豆精品久久久久久蜜桃| 欧美精品人与动牲交sv欧美| 国产 精品1| 亚洲精品亚洲一区二区| 99热这里只有是精品在线观看| 高清午夜精品一区二区三区| 毛片一级片免费看久久久久| 国产成人一区二区在线| 爱豆传媒免费全集在线观看| 国产精品国产三级专区第一集| 亚洲国产欧美人成| 久久6这里有精品| 久久精品国产a三级三级三级| av黄色大香蕉| 内射极品少妇av片p| 寂寞人妻少妇视频99o| 欧美亚洲 丝袜 人妻 在线| videossex国产| 欧美激情久久久久久爽电影| 午夜亚洲福利在线播放| 亚洲国产精品999| 97在线视频观看| 男人添女人高潮全过程视频| 午夜福利网站1000一区二区三区| 久久精品综合一区二区三区| 成人亚洲精品av一区二区| 波野结衣二区三区在线| 大片电影免费在线观看免费| 99热全是精品| 人人妻人人澡人人爽人人夜夜| 欧美bdsm另类| 国产亚洲精品久久久com| 最近最新中文字幕免费大全7| 老司机影院成人| 亚洲国产最新在线播放| 真实男女啪啪啪动态图| 日韩一本色道免费dvd| 在线观看免费高清a一片| 老师上课跳d突然被开到最大视频| 身体一侧抽搐| av在线天堂中文字幕| 免费黄网站久久成人精品| 身体一侧抽搐| 黄色怎么调成土黄色| 国产在线一区二区三区精| 十八禁网站网址无遮挡 | 别揉我奶头 嗯啊视频| 国产免费一级a男人的天堂| 欧美一区二区亚洲| av线在线观看网站| 夜夜爽夜夜爽视频| 街头女战士在线观看网站| av国产精品久久久久影院| 国产成人精品一,二区| 久久99精品国语久久久| 国产成人免费无遮挡视频| 夜夜爽夜夜爽视频| 黄片无遮挡物在线观看| 成人无遮挡网站| 精品久久久久久久人妻蜜臀av| 国产精品偷伦视频观看了| 少妇猛男粗大的猛烈进出视频 | 国产亚洲最大av| 亚洲人成网站在线播| 欧美区成人在线视频| 国产高潮美女av| 色播亚洲综合网| 尤物成人国产欧美一区二区三区| 免费观看在线日韩| 麻豆国产97在线/欧美| 国产人妻一区二区三区在| 日日摸夜夜添夜夜添av毛片| 亚洲自拍偷在线| 精品久久久久久久末码| 可以在线观看毛片的网站| 天堂中文最新版在线下载 | 亚洲av成人精品一区久久| 午夜日本视频在线| 少妇猛男粗大的猛烈进出视频 | 国产成人a∨麻豆精品| 亚洲精品成人av观看孕妇| 亚洲国产日韩一区二区| 九色成人免费人妻av| 国产一级毛片在线| 精品久久久久久久人妻蜜臀av|