• <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)航方法
    搞女人的毛片| 又爽又黄无遮挡网站| 深爱激情五月婷婷| 欧美av亚洲av综合av国产av| 久久香蕉精品热| 久久久久免费精品人妻一区二区| 搞女人的毛片| 韩国av一区二区三区四区| 亚洲不卡免费看| 69人妻影院| 久久久久久久久大av| 免费看a级黄色片| 国产伦在线观看视频一区| 婷婷丁香在线五月| 欧美成人一区二区免费高清观看| 淫秽高清视频在线观看| 国产私拍福利视频在线观看| 亚洲片人在线观看| 无遮挡黄片免费观看| 变态另类成人亚洲欧美熟女| 久久精品国产亚洲av涩爱 | 亚洲成a人片在线一区二区| 成人午夜高清在线视频| 国内久久婷婷六月综合欲色啪| 亚洲欧美一区二区三区黑人| 国产高清视频在线观看网站| 美女大奶头视频| 老熟妇仑乱视频hdxx| 欧美最新免费一区二区三区 | 欧美在线一区亚洲| 国内揄拍国产精品人妻在线| 老汉色∧v一级毛片| 激情在线观看视频在线高清| 日本一本二区三区精品| 色尼玛亚洲综合影院| 桃色一区二区三区在线观看| 精品熟女少妇八av免费久了| 精品欧美国产一区二区三| 一本一本综合久久| 成人午夜高清在线视频| 日本a在线网址| 国产精品三级大全| 欧美黄色片欧美黄色片| 一区福利在线观看| 欧美av亚洲av综合av国产av| 久久欧美精品欧美久久欧美| 1000部很黄的大片| 色哟哟哟哟哟哟| 亚洲人成网站在线播放欧美日韩| av天堂中文字幕网| 在线观看66精品国产| 亚洲国产精品合色在线| 91在线观看av| 久久久久久久久中文| 国产成人aa在线观看| 国产伦人伦偷精品视频| 天堂动漫精品| 男女午夜视频在线观看| 国产极品精品免费视频能看的| 日韩国内少妇激情av| netflix在线观看网站| 美女cb高潮喷水在线观看| 男女视频在线观看网站免费| 精品免费久久久久久久清纯| 天美传媒精品一区二区| 狠狠狠狠99中文字幕| 精品久久久久久久毛片微露脸| 18禁黄网站禁片免费观看直播| 亚洲成人久久性| 日本在线视频免费播放| 国产高清videossex| 成人鲁丝片一二三区免费| 国产精品久久视频播放| av在线蜜桃| 国产乱人伦免费视频| 又黄又粗又硬又大视频| 美女被艹到高潮喷水动态| 91在线观看av| 国产乱人伦免费视频| 激情在线观看视频在线高清| 欧美区成人在线视频| 全区人妻精品视频| 亚洲电影在线观看av| 国内精品一区二区在线观看| 老司机深夜福利视频在线观看| 老司机在亚洲福利影院| 精品不卡国产一区二区三区| 网址你懂的国产日韩在线| 在线a可以看的网站| 99国产精品一区二区三区| 久久久久免费精品人妻一区二区| 小蜜桃在线观看免费完整版高清| 国产精品99久久久久久久久| 欧美日韩瑟瑟在线播放| 国产乱人伦免费视频| 全区人妻精品视频| 久久久久久人人人人人| 久久精品人妻少妇| 日韩精品青青久久久久久| 国内少妇人妻偷人精品xxx网站| 午夜福利高清视频| 久久久久亚洲av毛片大全| 国产真实伦视频高清在线观看 | 99热这里只有是精品50| 欧美乱码精品一区二区三区| 国产精品久久久久久精品电影| 五月伊人婷婷丁香| 99视频精品全部免费 在线| 久久精品国产亚洲av涩爱 | 麻豆国产97在线/欧美| 日本一本二区三区精品| 欧美3d第一页| 色精品久久人妻99蜜桃| 欧美在线一区亚洲| 欧美性猛交╳xxx乱大交人| 91av网一区二区| 成人高潮视频无遮挡免费网站| 国产综合懂色| 亚洲男人的天堂狠狠| 日韩亚洲欧美综合| 中文字幕人成人乱码亚洲影| 美女免费视频网站| 嫩草影院入口| 亚洲欧美日韩高清在线视频| 午夜福利高清视频| 国产精品精品国产色婷婷| aaaaa片日本免费| 中文字幕高清在线视频| 亚洲欧美日韩无卡精品| 国产精品 国内视频| www.熟女人妻精品国产| 麻豆久久精品国产亚洲av| 久久久久久久久久黄片| 午夜精品久久久久久毛片777| 日日夜夜操网爽| 日韩免费av在线播放| 丰满乱子伦码专区| 成年女人毛片免费观看观看9| 欧美日本视频| 在线观看一区二区三区| 99久久无色码亚洲精品果冻| 亚洲人成网站高清观看| 精品电影一区二区在线| 一夜夜www| 亚洲人成网站在线播| 精品福利观看| 九色成人免费人妻av| 国产精品美女特级片免费视频播放器| 无限看片的www在线观看| 人人妻人人澡欧美一区二区| 久久亚洲真实| 啦啦啦韩国在线观看视频| 日韩欧美国产一区二区入口| 欧美zozozo另类| svipshipincom国产片| 熟女电影av网| 国产亚洲欧美在线一区二区| 中出人妻视频一区二区| 一区二区三区激情视频| 亚洲av日韩精品久久久久久密| aaaaa片日本免费| 欧美+日韩+精品| 国产真人三级小视频在线观看| 岛国在线免费视频观看| 国产成人av激情在线播放| 国内少妇人妻偷人精品xxx网站| 成人特级av手机在线观看| 午夜激情福利司机影院| 久久精品亚洲精品国产色婷小说| 国产免费av片在线观看野外av| www.www免费av| 男女那种视频在线观看| 人人妻人人澡欧美一区二区| 俺也久久电影网| 最新美女视频免费是黄的| 日本黄大片高清| 久久久国产精品麻豆| 国产97色在线日韩免费| 熟女电影av网| 国产精品女同一区二区软件 | 国产亚洲精品久久久com| av专区在线播放| 男人舔奶头视频| 我要搜黄色片| 黑人欧美特级aaaaaa片| 村上凉子中文字幕在线| 禁无遮挡网站| 又黄又粗又硬又大视频| 亚洲av免费在线观看| 女警被强在线播放| 欧美日韩精品网址| 国产亚洲欧美在线一区二区| 日本在线视频免费播放| 国产男靠女视频免费网站| 成人午夜高清在线视频| 一级毛片女人18水好多| 亚洲精品成人久久久久久| 免费人成视频x8x8入口观看| 可以在线观看毛片的网站| 国产三级在线视频| 男插女下体视频免费在线播放| 一区二区三区国产精品乱码| 久久久久久久久大av| 麻豆成人午夜福利视频| 又黄又爽又免费观看的视频| 哪里可以看免费的av片| 欧美极品一区二区三区四区| 成人av一区二区三区在线看| 女人十人毛片免费观看3o分钟| 久久久久精品国产欧美久久久| 99久久久亚洲精品蜜臀av| 免费大片18禁| 久久人妻av系列| 亚洲精品亚洲一区二区| 日韩欧美精品免费久久 | 亚洲av五月六月丁香网| 国产一区二区在线av高清观看| 国产69精品久久久久777片| 一个人看的www免费观看视频| 五月玫瑰六月丁香| 亚洲成a人片在线一区二区| 国产乱人伦免费视频| 麻豆国产97在线/欧美| 99久国产av精品| 国产老妇女一区| 色综合站精品国产| 国产在视频线在精品| 午夜福利在线在线| 国产精品久久视频播放| 国产伦一二天堂av在线观看| 亚洲18禁久久av| 国产亚洲欧美98| 三级男女做爰猛烈吃奶摸视频| 国产一区二区在线av高清观看| 午夜福利成人在线免费观看| 97碰自拍视频| 桃红色精品国产亚洲av| 少妇高潮的动态图| 国产精品爽爽va在线观看网站| 内地一区二区视频在线| 午夜福利在线在线| 亚洲最大成人手机在线| 国产精品国产高清国产av| 波多野结衣高清作品| 少妇的丰满在线观看| 亚洲精品在线美女| 久久久成人免费电影| 国产真实乱freesex| av片东京热男人的天堂| 成人18禁在线播放| 亚洲国产欧美网| 麻豆一二三区av精品| 在线观看66精品国产| 亚洲精品456在线播放app | 成人欧美大片| 国产探花极品一区二区| 91在线观看av| 看免费av毛片| 久久精品国产综合久久久| 亚洲片人在线观看| 在线看三级毛片| 亚洲国产高清在线一区二区三| 偷拍熟女少妇极品色| 久久人妻av系列| 搡老岳熟女国产| 熟女人妻精品中文字幕| 在线天堂最新版资源| 女生性感内裤真人,穿戴方法视频| 欧美日韩亚洲国产一区二区在线观看| 欧美日本亚洲视频在线播放| 成人av在线播放网站| 欧美不卡视频在线免费观看| 黄色视频,在线免费观看| 制服人妻中文乱码| 日韩 欧美 亚洲 中文字幕| 成人性生交大片免费视频hd| 在线观看免费视频日本深夜| 亚洲av熟女| 国产精品电影一区二区三区| 久久久久久久久久黄片| 高清在线国产一区| av国产免费在线观看| 狂野欧美白嫩少妇大欣赏| 欧美最新免费一区二区三区 | 亚洲欧美日韩卡通动漫| 欧美+日韩+精品| 有码 亚洲区| 久久久久精品国产欧美久久久| 色在线成人网| 狂野欧美白嫩少妇大欣赏| 国产亚洲精品久久久久久毛片| 香蕉丝袜av| 99久久精品一区二区三区| 久久99热这里只有精品18| av片东京热男人的天堂| 亚洲av成人不卡在线观看播放网| 好男人在线观看高清免费视频| 婷婷亚洲欧美| 在线观看美女被高潮喷水网站 | 国产免费男女视频| 淫秽高清视频在线观看| 免费人成在线观看视频色| 男女做爰动态图高潮gif福利片| x7x7x7水蜜桃| ponron亚洲| 国模一区二区三区四区视频| 国产黄片美女视频| 12—13女人毛片做爰片一| 国产乱人视频| 最近视频中文字幕2019在线8| 欧美色欧美亚洲另类二区| 99精品在免费线老司机午夜| 精品不卡国产一区二区三区| 亚洲成人久久性| 看片在线看免费视频| 高潮久久久久久久久久久不卡| a级一级毛片免费在线观看| 亚洲国产高清在线一区二区三| 99久久精品热视频| 岛国在线免费视频观看| 蜜桃久久精品国产亚洲av| 成人av一区二区三区在线看| 少妇裸体淫交视频免费看高清| 精品久久久久久久久久久久久| 在线观看av片永久免费下载| 亚洲欧美一区二区三区黑人| 两个人看的免费小视频| 国产三级中文精品| 久久6这里有精品| 国产色婷婷99| 2021天堂中文幕一二区在线观| 国产精品久久视频播放| 91在线观看av| 亚洲熟妇中文字幕五十中出| 午夜精品一区二区三区免费看| 国产精品一及| 国内久久婷婷六月综合欲色啪| 五月伊人婷婷丁香| 欧美日韩亚洲国产一区二区在线观看| 一a级毛片在线观看| 欧美黄色淫秽网站| 亚洲美女黄片视频| 成人性生交大片免费视频hd| 欧美激情久久久久久爽电影| 九色成人免费人妻av| av片东京热男人的天堂| 青草久久国产| 中文字幕av在线有码专区| 国产精品久久久人人做人人爽| 少妇人妻精品综合一区二区 | 国产成年人精品一区二区| 99久久精品国产亚洲精品| 欧美日韩精品网址| 久久精品综合一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 国产精品综合久久久久久久免费| 日本免费a在线| 国产 一区 欧美 日韩| 男人舔奶头视频| 中亚洲国语对白在线视频| 波多野结衣巨乳人妻| 黄色女人牲交| 51国产日韩欧美| 欧美不卡视频在线免费观看| 亚洲精品粉嫩美女一区| 欧美成人免费av一区二区三区| 全区人妻精品视频| 亚洲国产欧美人成| 蜜桃亚洲精品一区二区三区| 国产精品98久久久久久宅男小说| 国产一区二区亚洲精品在线观看| 婷婷丁香在线五月| 一进一出抽搐gif免费好疼| 亚洲av第一区精品v没综合| 日本一本二区三区精品| 精华霜和精华液先用哪个| tocl精华| 伊人久久精品亚洲午夜| 国产高清视频在线观看网站| 亚洲欧美精品综合久久99| 老汉色∧v一级毛片| 亚洲欧美日韩卡通动漫| 亚洲av中文字字幕乱码综合| 久久精品国产清高在天天线| 中出人妻视频一区二区| 久久这里只有精品中国| 99久久成人亚洲精品观看| 九九久久精品国产亚洲av麻豆| 蜜桃久久精品国产亚洲av| 婷婷六月久久综合丁香| 尤物成人国产欧美一区二区三区| 免费av毛片视频| 国产av一区在线观看免费| www.色视频.com| 99久久成人亚洲精品观看| 丰满人妻一区二区三区视频av | 亚洲成人中文字幕在线播放| 99久久综合精品五月天人人| av中文乱码字幕在线| 女人被狂操c到高潮| 亚洲自拍偷在线| 久久精品人妻少妇| 国内精品久久久久久久电影| 村上凉子中文字幕在线| 亚洲久久久久久中文字幕| 在线观看一区二区三区| 国产aⅴ精品一区二区三区波| 欧美黄色片欧美黄色片| 国产成人av激情在线播放| av欧美777| 欧美中文综合在线视频| 久久久久久久久久黄片| 99久久精品一区二区三区| 90打野战视频偷拍视频| 亚洲va日本ⅴa欧美va伊人久久| 狠狠狠狠99中文字幕| 岛国在线观看网站| 91久久精品国产一区二区成人 | 搡老妇女老女人老熟妇| 手机成人av网站| 久久久久国内视频| 国产午夜福利久久久久久| 成人av在线播放网站| 久久欧美精品欧美久久欧美| 亚洲第一欧美日韩一区二区三区| 神马国产精品三级电影在线观看| 久久久国产成人精品二区| 欧美av亚洲av综合av国产av| aaaaa片日本免费| 搡老岳熟女国产| 亚洲不卡免费看| 久久中文看片网| 精品久久久久久久久久免费视频| 国产成人福利小说| 亚洲中文日韩欧美视频| 国产高清有码在线观看视频| h日本视频在线播放| 亚洲五月婷婷丁香| 我的老师免费观看完整版| 51国产日韩欧美| 无遮挡黄片免费观看| 色综合欧美亚洲国产小说| www.色视频.com| 国产色婷婷99| 国产精品精品国产色婷婷| 欧美精品啪啪一区二区三区| 午夜亚洲福利在线播放| 一区二区三区国产精品乱码| 亚洲精品美女久久久久99蜜臀| 国产av一区在线观看免费| 亚洲欧美日韩无卡精品| 特级一级黄色大片| 国产高清视频在线观看网站| 亚洲欧美精品综合久久99| 成人国产一区最新在线观看| 最近最新免费中文字幕在线| 一个人看视频在线观看www免费 | 午夜视频国产福利| 一进一出好大好爽视频| 国产在线精品亚洲第一网站| 琪琪午夜伦伦电影理论片6080| 又黄又粗又硬又大视频| 亚洲精品久久国产高清桃花| 一区二区三区高清视频在线| 欧美丝袜亚洲另类 | 在线观看av片永久免费下载| 日本撒尿小便嘘嘘汇集6| 国产成人系列免费观看| 亚洲第一电影网av| 亚洲18禁久久av| netflix在线观看网站| 日本a在线网址| 亚洲专区国产一区二区| 99精品久久久久人妻精品| 神马国产精品三级电影在线观看| 欧美成人免费av一区二区三区| 欧美在线黄色| 最近最新中文字幕大全电影3| 长腿黑丝高跟| 国产伦一二天堂av在线观看| www.www免费av| 国产午夜福利久久久久久| 青草久久国产| 日本黄大片高清| 欧美日韩福利视频一区二区| 国产爱豆传媒在线观看| 久久久久性生活片| 村上凉子中文字幕在线| 成人鲁丝片一二三区免费| 国产精品久久久久久亚洲av鲁大| 热99re8久久精品国产| 一进一出抽搐gif免费好疼| 男人和女人高潮做爰伦理| 九九在线视频观看精品| 嫩草影院入口| 国语自产精品视频在线第100页| 观看美女的网站| 久久精品国产综合久久久| 国产精品精品国产色婷婷| 久久精品国产综合久久久| 国产精品爽爽va在线观看网站| 少妇的逼水好多| a级毛片a级免费在线| 真人做人爱边吃奶动态| 亚洲熟妇熟女久久| 怎么达到女性高潮| АⅤ资源中文在线天堂| 国产成人a区在线观看| 五月玫瑰六月丁香| 啦啦啦观看免费观看视频高清| 男女床上黄色一级片免费看| 亚洲在线观看片| 哪里可以看免费的av片| 美女被艹到高潮喷水动态| 国产av不卡久久| 老司机深夜福利视频在线观看| 嫁个100分男人电影在线观看| 日韩人妻高清精品专区| 婷婷丁香在线五月| 长腿黑丝高跟| 亚洲国产精品合色在线| 国产淫片久久久久久久久 | 国产精品免费一区二区三区在线| 国产成人av激情在线播放| 制服丝袜大香蕉在线| 一区二区三区激情视频| 性欧美人与动物交配| 日韩欧美一区二区三区在线观看| 国产精品野战在线观看| 免费观看人在逋| 亚洲精品一卡2卡三卡4卡5卡| 18禁黄网站禁片午夜丰满| 国内揄拍国产精品人妻在线| 天堂av国产一区二区熟女人妻| 免费人成视频x8x8入口观看| 国产爱豆传媒在线观看| 蜜桃久久精品国产亚洲av| 亚洲国产欧美人成| 亚洲av二区三区四区| 成人性生交大片免费视频hd| 亚洲片人在线观看| 亚洲中文日韩欧美视频| 欧美成人一区二区免费高清观看| 人人妻,人人澡人人爽秒播| 国产久久久一区二区三区| 丰满乱子伦码专区| 麻豆成人av在线观看| 在线看三级毛片| 免费电影在线观看免费观看| 亚洲成av人片在线播放无| 女人十人毛片免费观看3o分钟| 日韩欧美 国产精品| 欧美另类亚洲清纯唯美| 久久精品综合一区二区三区| 成年版毛片免费区| 中文字幕精品亚洲无线码一区| 免费av毛片视频| 免费av观看视频| 国模一区二区三区四区视频| 午夜福利在线在线| 色尼玛亚洲综合影院| АⅤ资源中文在线天堂| 亚洲国产精品合色在线| 一进一出抽搐动态| 国产高清激情床上av| 在线国产一区二区在线| 欧美成人一区二区免费高清观看| 男女午夜视频在线观看| 精品久久久久久久久久久久久| 国产精品精品国产色婷婷| 欧美区成人在线视频| 欧美色欧美亚洲另类二区| 国产精品一及| 精品久久久久久久人妻蜜臀av| 欧美成人一区二区免费高清观看| 国产三级在线视频| 午夜精品一区二区三区免费看| 久久久久久久精品吃奶| 国产三级中文精品| 国产高清三级在线| 日本免费a在线| 12—13女人毛片做爰片一| 亚洲不卡免费看| 乱人视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 欧美一级毛片孕妇| 欧美日韩中文字幕国产精品一区二区三区| 国语自产精品视频在线第100页| 很黄的视频免费| avwww免费| 99热这里只有是精品50| a级毛片a级免费在线| 人妻久久中文字幕网| 久久国产精品影院| tocl精华| 亚洲真实伦在线观看| 免费看日本二区| 久久婷婷人人爽人人干人人爱| 真实男女啪啪啪动态图| 十八禁网站免费在线| 老司机深夜福利视频在线观看| 一边摸一边抽搐一进一小说| 中文字幕av成人在线电影| 久久精品国产99精品国产亚洲性色| 男女做爰动态图高潮gif福利片| 99热精品在线国产| 日本免费a在线| 一级黄色大片毛片| 日韩高清综合在线| 精品一区二区三区人妻视频| 亚洲男人的天堂狠狠| 美女黄网站色视频| 久久久久久久久中文| 99国产精品一区二区蜜桃av| 在线观看免费视频日本深夜|