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

    融合光流與MARG 傳感器的三維姿態(tài)測(cè)量研究*

    2023-11-21 13:07:28琦,李
    電子器件 2023年5期
    關(guān)鍵詞:光流陀螺儀加速度計(jì)

    石 琦,李 翔

    (桂林電子科技大學(xué)電子工程與自動(dòng)化學(xué)院,廣西 桂林 541004)

    光流傳感器是一種能從連續(xù)圖像流中獲取運(yùn)動(dòng)信息的傳感器,近年來(lái)被廣泛用于無(wú)人機(jī)中,與加速度計(jì)、陀螺儀、磁強(qiáng)計(jì)(合稱MARG 傳感器,即磁強(qiáng)計(jì)magnetometer、加速度計(jì)accelerometer 及速率陀螺rate gyro 的縮寫[1],亦可解釋為磁magnetic、角速度angular rate 和重力gravity 三種測(cè)量對(duì)象的縮寫[2-4])等傳感器相結(jié)合,以實(shí)現(xiàn)無(wú)人機(jī)的飛行控制。光流傳感器用于無(wú)人機(jī)飛控的一種典型方案是采用豎直向下安裝的單個(gè)光流傳感器,通過(guò)其測(cè)量值獲取無(wú)人機(jī)相對(duì)于地面的水平運(yùn)動(dòng),可提高無(wú)人機(jī)定點(diǎn)懸停的位置精度[5-6]。但在此種方案中,光流傳感器的測(cè)量值實(shí)際上是無(wú)人機(jī)的平移、俯仰和橫滾三種運(yùn)動(dòng)共同作用的結(jié)果。因此,要想從光流測(cè)量值準(zhǔn)確獲取無(wú)人機(jī)的平移運(yùn)動(dòng),需要借助慣導(dǎo)(Inertial Navigation System,INS)提供的姿態(tài)信息或者陀螺儀提供的角速度測(cè)量值,將無(wú)人機(jī)姿態(tài)變化引起的光流加以消除,才能得到準(zhǔn)確的水平速度(或位移)[7-11]。

    如上所述,光流傳感器的測(cè)量值可分為平移分量與轉(zhuǎn)動(dòng)分量?jī)刹糠郑謩e與載體的平移運(yùn)動(dòng)和旋轉(zhuǎn)運(yùn)動(dòng)相關(guān)。然而上述將光流傳感器用于無(wú)人機(jī)的方案中,只用到了平移運(yùn)動(dòng)形成的光流(平動(dòng)光流),而對(duì)于載體轉(zhuǎn)動(dòng)造成的光流(轉(zhuǎn)動(dòng)光流)則借助IMU 陀螺儀提供的角速度將其消去。此類用法的缺點(diǎn)在于:一是未能充分利用轉(zhuǎn)動(dòng)光流包含的姿態(tài)信息;二是在消去轉(zhuǎn)動(dòng)光流時(shí)易將慣導(dǎo)或陀螺儀的誤差代入到平移速度(或位移)的計(jì)算中。

    在航天器控制領(lǐng)域,利用轉(zhuǎn)動(dòng)光流獲取載體姿態(tài)信息的方法較為多見(jiàn),如文獻(xiàn)[12]采用擴(kuò)展卡爾曼濾波器(Extended Kalman Filter,EKF)由光流估計(jì)三軸角速度,文獻(xiàn)[13]和[14]采用多種不同的卡爾曼濾波算法由光流估計(jì)航天器的三軸角速度,文獻(xiàn)[15]對(duì)星空成像進(jìn)行預(yù)處理后再提取光流數(shù)據(jù)以計(jì)算角速度。此類方法的問(wèn)題是成像范圍內(nèi)星體數(shù)目較多時(shí),獲取光流所需計(jì)算量較大,且涉及非線性測(cè)量模型,需采用非線性的卡爾曼濾波算法(如EKF)。在無(wú)人機(jī)領(lǐng)域,文獻(xiàn)[16]嘗試采用光流傳感器測(cè)量無(wú)人機(jī)三維姿態(tài),然而其所用算法也采用非線性測(cè)量模型,同樣存在求解復(fù)雜、計(jì)算量大等問(wèn)題。

    為充分利用光流傳感器測(cè)量值中包含的姿態(tài)信息,提出一種簡(jiǎn)單的線性測(cè)量模型,可從光流測(cè)量值中提取角速度,并進(jìn)而將獲得的角速度與陀螺儀測(cè)量值根據(jù)最小方差準(zhǔn)則進(jìn)行融合以減小角速度測(cè)量誤差,最后將融合后的角速度與加速度計(jì)及磁強(qiáng)計(jì)數(shù)據(jù)利用EKF 算法進(jìn)行融合,以提高三維姿態(tài)的估計(jì)精度。通過(guò)實(shí)驗(yàn)驗(yàn)證了所提出的測(cè)量模型和算法的有效性。

    1 光流傳感器測(cè)量模型

    光流傳感器工作原理如圖1 所示,由于通常情況下物距d遠(yuǎn)大于焦距f,故可認(rèn)為像距近似等于f。由圖1 可見(jiàn),凡是垂直于傳感器主光軸的運(yùn)動(dòng)均可產(chǎn)生光流,而旋轉(zhuǎn)運(yùn)動(dòng)與平移運(yùn)動(dòng)均可切割主光軸,故旋轉(zhuǎn)與平移均可產(chǎn)生相應(yīng)的光流。因此,可將光流傳感器的測(cè)量值視為平移光流與轉(zhuǎn)動(dòng)光流的疊加。

    圖1 光流測(cè)量原理

    首先考慮傳感器載體僅有平移運(yùn)動(dòng)的情況。設(shè)垂直于主光軸的平移速度為v,由平移造成的光流為δv,則由圖1 所示幾何關(guān)系可知δv正比于v,如式(1)所示,其中k為沿主光軸的單位矢量。可見(jiàn)平移光流δv與物距d成反比。

    其次考慮傳感器載體轉(zhuǎn)動(dòng)形成的光流。設(shè)載體角速度為ω,由其造成的光流為δω,經(jīng)分析可知二者關(guān)系如式(2)所示。可見(jiàn)轉(zhuǎn)動(dòng)光流δω與物距d無(wú)關(guān)。

    光流傳感器的測(cè)量值應(yīng)是上述兩部分光流的疊加,即

    2 基于光流的角速度測(cè)量

    由式(3)可知,當(dāng)傳感器載體只進(jìn)行旋轉(zhuǎn)運(yùn)動(dòng)時(shí),光流測(cè)量值僅與角速度有關(guān)。光流與角速度的關(guān)系式即式(2)可改寫為矩陣形式,如式(4)所示。由于角速度ω為三維矢量,而轉(zhuǎn)動(dòng)光流δω是像平面內(nèi)的二維矢量,故式(4)中的Uω是一個(gè)2×3 常系數(shù)矩陣,它反映了光流傳感器在b系中的安裝方向。

    將式(4)推廣到采用多個(gè)不共面光流傳感器進(jìn)行測(cè)量的情形,則得到式(5):

    式中:δωk表示第k個(gè)光流傳感器的測(cè)量值(其中k=1,2,…,N),Uωk為第k個(gè)光流傳感器對(duì)應(yīng)的常系數(shù)矩陣。注意每個(gè)光流傳感器的測(cè)量值(即δωk)都是二維矢量(包含像平面內(nèi)兩個(gè)正交方向的光流測(cè)量數(shù)據(jù)),故Vω為2N×3 的常系數(shù)矩陣,由各光流傳感器在b系中的安裝角度確定。式(5)可視為關(guān)于角速度的超定方程組,可利用最小二乘法進(jìn)行求解,如式(6)所示:

    式中:ωof表示由光流測(cè)量值解算出的角速度。實(shí)際應(yīng)用中,式(6)還可以進(jìn)一步簡(jiǎn)寫為式(7),其中W=是一個(gè)3×2N的常系數(shù)矩陣,只要能獲知角速度的真值,即可直接由光流測(cè)量值和相應(yīng)的角速度真值求出W。

    3 光流與MARG 傳感器的融合

    3.1 算法整體流程

    如上所述,光流傳感器可用于測(cè)量三維角速度,將其與MARG 傳感器進(jìn)行數(shù)據(jù)融合有助于改善姿態(tài)估計(jì)精度。為此,設(shè)計(jì)了如下的數(shù)據(jù)融合與姿態(tài)估計(jì)流程:

    ①?gòu)墓饬鱾鞲衅鳒y(cè)量數(shù)據(jù)中提取三維角速度;

    ②將光流傳感器獲取的角速度與陀螺儀測(cè)量得到的角速度進(jìn)行融合;

    ③將融合后的角速度數(shù)據(jù)與加速度計(jì)及磁強(qiáng)計(jì)測(cè)量值一同代入EKF 估計(jì)姿態(tài)。

    算法整體流程框圖如圖2 所示,其中:ωof表示由光流測(cè)量值提取的角速度,ωg為陀螺儀測(cè)得的角速度,ωs為二者融合后的角速度,g為加速度計(jì)測(cè)得的重力矢量,h為磁強(qiáng)計(jì)測(cè)得的地磁矢量。

    圖2 數(shù)據(jù)融合及姿態(tài)估計(jì)算法流程框圖

    3.2 基于最小方差準(zhǔn)則的角速度融合

    光流傳感器和陀螺儀均能測(cè)量載體的三維角速度,將兩者的測(cè)量值進(jìn)行加權(quán)可提高測(cè)量精度,具體做法是按最小方差準(zhǔn)則進(jìn)行融合,如式(8)所示:

    3.3 基于EKF 的姿態(tài)估計(jì)

    選取姿態(tài)四元數(shù)q=(q1q2q3q4)T作為狀態(tài)向量(q1為標(biāo)量部分,其余三個(gè)分量為矢量部分),利用EKF 進(jìn)行姿態(tài)估計(jì),遞推流程如下。

    EKF 算法第一步(狀態(tài)預(yù)測(cè)):

    式中:表示(k-1)時(shí)刻的驗(yàn)后估計(jì),表示k時(shí)刻的狀態(tài)預(yù)測(cè)值即驗(yàn)前估計(jì),Φk/k-1表示從(k-1)時(shí)刻到k時(shí)刻的狀態(tài)轉(zhuǎn)移矩陣,如式(10)所示,ωsx,k、ωsy,k和ωsz,k分別表示k時(shí)刻ωs的三軸分量,T為采樣周期。

    式中:Pk-1/k-1為(k-1)時(shí)刻的驗(yàn)后方差,Pk/k-1為k時(shí)刻的驗(yàn)前方差,Qk為過(guò)程噪聲協(xié)方差陣。Qk由式(12)計(jì)算,其中Rω為角速度的協(xié)方差陣,Sk是由k時(shí)刻的姿態(tài)四元數(shù)按式(13)構(gòu)造的雅可比矩陣,q1~q4依次為k時(shí)刻姿態(tài)四元數(shù)的四個(gè)分量。

    EKF 算法第三步(狀態(tài)修正):

    由k時(shí)刻的加速度計(jì)測(cè)量值gk與磁強(qiáng)計(jì)測(cè)量值hk構(gòu)成EKF 的觀測(cè)量Zk,如式(14)所示:

    式中:量測(cè)矩陣Hk可由式(16)求得。由于和分別為4 維和6 維向量,故Hk為6×4 矩陣。式(16)中:g0為本地重力加速度,hN及hU分別為地磁場(chǎng)的北向和天向分量。

    接下來(lái)由式(17)求出k時(shí)刻的卡爾曼增益Kk,其中Rk為量測(cè)噪聲協(xié)方差陣,由磁強(qiáng)計(jì)的噪聲協(xié)方差陣Rh及加速度計(jì)的協(xié)方差陣Rg根據(jù)式(18)構(gòu)造而成。

    由卡爾曼增益可得k時(shí)刻的驗(yàn)后估計(jì)(狀態(tài)修正)如式(19)所示:

    4 實(shí)驗(yàn)驗(yàn)證

    4.1 實(shí)驗(yàn)條件

    實(shí)驗(yàn)采用的光流模塊型號(hào)為L(zhǎng)C302,將三個(gè)光流模塊布置于兩兩正交的三個(gè)平面上,如圖3 所示。此外,實(shí)驗(yàn)中還采用型號(hào)為IMU901 的MARG 傳感器模塊與圖3 所示光流模塊進(jìn)行數(shù)據(jù)融合與姿態(tài)解算。

    圖3 實(shí)驗(yàn)所用光流傳感器模塊

    采用型號(hào)為WTGAHRS1 的慣導(dǎo)模塊提供三維姿態(tài)基準(zhǔn),其姿態(tài)角精度可達(dá)靜態(tài)0.05°、動(dòng)態(tài)0.1°,可作為三維姿態(tài)的真值,實(shí)驗(yàn)中將其與光流模塊及航姿模塊互相固定,使得三者能夠同步轉(zhuǎn)動(dòng)。數(shù)據(jù)采集時(shí),手持模塊進(jìn)行無(wú)規(guī)則轉(zhuǎn)動(dòng),使用串口接收數(shù)據(jù)進(jìn)行統(tǒng)一處理,并將數(shù)據(jù)導(dǎo)入MATLAB 中進(jìn)行解算。由WTGAHRS1 模塊測(cè)得的三維姿態(tài)變化如圖4 所示。

    圖4 三維姿態(tài)角變化曲線

    4.2 角速度測(cè)量及融合結(jié)果

    以慣導(dǎo)模塊WTGAHRS1 提供的三維角速度數(shù)據(jù)作為真值,首先求出式(7)中的常系數(shù)矩陣W,如式(21)所示:

    得到矩陣W后,再將其代入式(7),即可從光流傳感器測(cè)量值解算出角速度ωof。

    接下來(lái)按式(8)將ωof與陀螺儀測(cè)量的角速度ωg進(jìn)行融合得到ωs。由實(shí)測(cè)數(shù)據(jù)分別計(jì)算ωof和ωg的均方根誤差(Root Mean Square Error,RMSE),結(jié)果如表1 所示。將各軸角速度的RMSE 取平均后代入式(8)即可得到加權(quán)系數(shù),從而由ωof和ωg加權(quán)后得到ωs,其RMSE 亦列于表1 中。由表1 可見(jiàn),式(7)所給出的線性測(cè)量模型能夠有效地從多個(gè)不共面安裝的光流傳感器測(cè)量值中提取三維角速度信息,且測(cè)量精度與低成本MEMS 陀螺儀相當(dāng)。將光流傳感器獲取的角速度信息與陀螺儀測(cè)量值進(jìn)行融合后,可進(jìn)一步提高精度,使三軸角速度的誤差限制在1°/s 以內(nèi)。

    表1 三維角速度誤差(RMSE) 單位:(°)/s

    4.3 EKF 姿態(tài)估計(jì)結(jié)果

    采用3.3 小節(jié)中的EKF 算法進(jìn)行三維姿態(tài)估計(jì),各傳感器的噪聲協(xié)方差陣取值為:Rω=2.5×10-4I3×3,Rh=4×10-4I3×3,Rg=10-6I3×3,其中角速度的噪聲協(xié)方差參考表1 所列RMSE,重力及地磁矢量的噪聲協(xié)方差則由加速度計(jì)和磁強(qiáng)計(jì)噪聲分別經(jīng)過(guò)矢量單位化處理后得到。此外,對(duì)重力及地磁矢量均作單位化處理后,式(16)所示量測(cè)矩陣Hk中的g0相當(dāng)于變?yōu)?,而hN及hU則直接由本地的地磁傾角γ確定,實(shí)驗(yàn)中測(cè)得γ≈39°,即hN≈cos39°,hU≈-sin39°。

    在EKF 算法中分別采用兩種不同來(lái)源的角速度數(shù)據(jù):陀螺儀測(cè)量值ωg及其與光流傳感器融合后的角速度ωs。將EKF 輸出的姿態(tài)四元數(shù)轉(zhuǎn)化為歐拉角,與慣導(dǎo)模塊WTGAHRS1 提供的姿態(tài)角真值進(jìn)行比較,結(jié)果如表2 所示。可見(jiàn),引入光流解算的角速度后,有助于提高姿態(tài)測(cè)量的精度。

    表2 三維姿態(tài)角誤差(RMSE) 單位:(°)

    5 結(jié)束語(yǔ)

    提出一種利用光流傳感器解算角速度并與陀螺儀進(jìn)行融合的方法,采用線性測(cè)量模型由光流測(cè)量值解算出三維角速度,并將其與MARG 傳感器中的陀螺儀按最小方差準(zhǔn)則進(jìn)行融合,再以融合后的角速度構(gòu)建EKF 的狀態(tài)方程,以加速度計(jì)與磁強(qiáng)計(jì)數(shù)據(jù)構(gòu)建EKF 量測(cè)方程。實(shí)驗(yàn)結(jié)果表明,所提出的方法改善了基于MARG 傳感器的姿態(tài)估計(jì)精度,提供了一種姿態(tài)測(cè)量的新途徑,有助于提升姿態(tài)測(cè)量系統(tǒng)的冗余度。

    猜你喜歡
    光流陀螺儀加速度計(jì)
    利用掩膜和單應(yīng)矩陣提高LK光流追蹤效果
    基于加速度計(jì)的起重機(jī)制動(dòng)下滑量測(cè)量
    基于EMD的MEMS陀螺儀隨機(jī)漂移分析方法
    基于物理學(xué)的改善粒子圖像測(cè)速穩(wěn)健光流方法研究
    我國(guó)著名陀螺儀專家——林士諤
    基于遺傳算法的加速度計(jì)免轉(zhuǎn)臺(tái)標(biāo)定方法
    微機(jī)械陀螺儀概述和發(fā)展
    常溫下硅微諧振加速度計(jì)零偏穩(wěn)定性的提高
    MEMS三軸陀螺儀中不匹配干擾抑制方法
    融合光流速度場(chǎng)與背景差分的自適應(yīng)背景更新方法
    久久久久精品国产欧美久久久| 99久久久亚洲精品蜜臀av| 搡老岳熟女国产| 97超视频在线观看视频| 夜夜爽天天搞| 欧美极品一区二区三区四区| 搡老熟女国产l中国老女人| 国产女主播在线喷水免费视频网站 | 男人和女人高潮做爰伦理| 国产久久久一区二区三区| 日韩在线高清观看一区二区三区| 91午夜精品亚洲一区二区三区| 91麻豆精品激情在线观看国产| 一区二区三区四区激情视频 | 亚洲性久久影院| 天天躁夜夜躁狠狠久久av| 搡女人真爽免费视频火全软件 | 真人做人爱边吃奶动态| 亚洲欧美精品自产自拍| 亚洲国产欧洲综合997久久,| 精品人妻熟女av久视频| 亚洲人成网站在线观看播放| ponron亚洲| 亚洲精品色激情综合| 一边摸一边抽搐一进一小说| 亚洲第一电影网av| 人妻制服诱惑在线中文字幕| 2021天堂中文幕一二区在线观| 精品久久国产蜜桃| 午夜激情欧美在线| 欧美人与善性xxx| 在线观看av片永久免费下载| 狂野欧美激情性xxxx在线观看| 国产人妻一区二区三区在| 中文字幕人妻熟人妻熟丝袜美| 少妇熟女aⅴ在线视频| 婷婷精品国产亚洲av在线| 女生性感内裤真人,穿戴方法视频| 久久天躁狠狠躁夜夜2o2o| 干丝袜人妻中文字幕| 国产在线男女| 91午夜精品亚洲一区二区三区| 波多野结衣高清作品| 高清毛片免费观看视频网站| 人妻丰满熟妇av一区二区三区| 欧美区成人在线视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲性夜色夜夜综合| 午夜免费激情av| 欧美国产日韩亚洲一区| av天堂中文字幕网| 久久久久久久久久黄片| 一本一本综合久久| 99久久精品一区二区三区| 人妻夜夜爽99麻豆av| 国产精品一区二区三区四区久久| 99在线视频只有这里精品首页| av在线老鸭窝| 亚洲国产高清在线一区二区三| 国产精品久久久久久精品电影| 国产精品久久视频播放| 又爽又黄a免费视频| .国产精品久久| www.色视频.com| 亚洲高清免费不卡视频| 日本一本二区三区精品| 精品熟女少妇av免费看| 国产一区二区在线av高清观看| 最近在线观看免费完整版| 精品一区二区三区av网在线观看| 日韩av不卡免费在线播放| 国产一区二区亚洲精品在线观看| 国产麻豆成人av免费视频| 人人妻人人看人人澡| 嫩草影视91久久| 一级毛片电影观看 | 国产精品日韩av在线免费观看| 能在线免费观看的黄片| 国产午夜精品论理片| 亚洲成人精品中文字幕电影| av卡一久久| 国产av一区在线观看免费| 成人亚洲欧美一区二区av| 久久久久久伊人网av| 欧美成人一区二区免费高清观看| 国内精品久久久久精免费| 国内精品美女久久久久久| 十八禁网站免费在线| 国产真实乱freesex| 啦啦啦韩国在线观看视频| 国内少妇人妻偷人精品xxx网站| 亚洲18禁久久av| 1024手机看黄色片| 久久久久免费精品人妻一区二区| eeuss影院久久| 国产成人a∨麻豆精品| 色尼玛亚洲综合影院| 亚洲人成网站在线播放欧美日韩| 久久久久久久午夜电影| 国产精品久久久久久久久免| 全区人妻精品视频| 1024手机看黄色片| 精品一区二区三区av网在线观看| 国产午夜精品久久久久久一区二区三区 | 欧美一区二区国产精品久久精品| 国产单亲对白刺激| 夜夜看夜夜爽夜夜摸| 国产精品一区二区三区四区久久| 三级毛片av免费| 99久久精品一区二区三区| 亚洲精品一区av在线观看| 国产精品精品国产色婷婷| 亚洲一区高清亚洲精品| 人人妻人人澡人人爽人人夜夜 | 男人狂女人下面高潮的视频| 国产精品亚洲一级av第二区| 亚洲国产欧美人成| 日本黄色片子视频| 又黄又爽又免费观看的视频| 国产私拍福利视频在线观看| 别揉我奶头~嗯~啊~动态视频| 午夜视频国产福利| 亚洲精品成人久久久久久| 毛片一级片免费看久久久久| 美女被艹到高潮喷水动态| 欧美又色又爽又黄视频| 久久韩国三级中文字幕| 国产伦精品一区二区三区四那| 久久久久久九九精品二区国产| 成人漫画全彩无遮挡| 欧美3d第一页| 国产精品亚洲一级av第二区| 天堂动漫精品| 久久精品国产99精品国产亚洲性色| 一个人看的www免费观看视频| 国产精品久久久久久精品电影| 日韩亚洲欧美综合| 搞女人的毛片| 免费高清视频大片| 狂野欧美白嫩少妇大欣赏| 欧美日本视频| 国产精品久久久久久精品电影| 久久久国产成人免费| 18+在线观看网站| 成人三级黄色视频| 美女cb高潮喷水在线观看| 成人欧美大片| 国产高清激情床上av| 男插女下体视频免费在线播放| 麻豆一二三区av精品| 国产精品福利在线免费观看| 亚洲av免费高清在线观看| 麻豆一二三区av精品| 搡老岳熟女国产| 亚洲精品久久国产高清桃花| 亚洲精华国产精华液的使用体验 | 91在线精品国自产拍蜜月| 一本精品99久久精品77| av天堂在线播放| 久久韩国三级中文字幕| a级一级毛片免费在线观看| 亚洲国产精品sss在线观看| 六月丁香七月| 日韩成人伦理影院| 国产精品,欧美在线| 51国产日韩欧美| 搞女人的毛片| 国产一区二区激情短视频| 女人被狂操c到高潮| 免费av不卡在线播放| 精品久久久久久久久av| 亚洲人成网站在线播| 亚洲国产高清在线一区二区三| 国产精品爽爽va在线观看网站| 欧美成人免费av一区二区三区| 久久精品91蜜桃| 久久久久久久午夜电影| 老师上课跳d突然被开到最大视频| 一进一出抽搐gif免费好疼| 我要看日韩黄色一级片| 十八禁国产超污无遮挡网站| 久久人人精品亚洲av| videossex国产| 美女被艹到高潮喷水动态| 欧美日韩精品成人综合77777| 99久久精品国产国产毛片| 成人精品一区二区免费| 国产伦精品一区二区三区视频9| 一级毛片aaaaaa免费看小| 日韩一本色道免费dvd| 听说在线观看完整版免费高清| 九九爱精品视频在线观看| 人妻夜夜爽99麻豆av| 校园人妻丝袜中文字幕| 一级毛片aaaaaa免费看小| 婷婷色综合大香蕉| 欧美色视频一区免费| 伦精品一区二区三区| 国产精品一及| 成年免费大片在线观看| 三级国产精品欧美在线观看| 久久6这里有精品| 成人美女网站在线观看视频| 久久人人精品亚洲av| 亚洲精品国产av成人精品 | 欧美性感艳星| 在线播放国产精品三级| 91av网一区二区| 午夜免费激情av| 在线播放无遮挡| 久久99热这里只有精品18| av女优亚洲男人天堂| 一个人看视频在线观看www免费| 在线免费观看的www视频| 亚洲人与动物交配视频| 欧美不卡视频在线免费观看| 床上黄色一级片| 成人二区视频| 1000部很黄的大片| 成人无遮挡网站| 久久亚洲国产成人精品v| 久久精品国产亚洲av香蕉五月| 能在线免费观看的黄片| 欧美日韩综合久久久久久| 狂野欧美激情性xxxx在线观看| 久久久久精品国产欧美久久久| av中文乱码字幕在线| 亚洲精品久久国产高清桃花| 尾随美女入室| 哪里可以看免费的av片| 国内久久婷婷六月综合欲色啪| 又粗又爽又猛毛片免费看| 欧美成人一区二区免费高清观看| 成人高潮视频无遮挡免费网站| 久久精品影院6| 女的被弄到高潮叫床怎么办| 蜜臀久久99精品久久宅男| 亚洲精品国产av成人精品 | 最后的刺客免费高清国语| 尾随美女入室| 丰满乱子伦码专区| 亚洲人成网站高清观看| 一进一出好大好爽视频| 啦啦啦啦在线视频资源| 亚洲激情五月婷婷啪啪| 国内精品宾馆在线| 女同久久另类99精品国产91| 岛国在线免费视频观看| 午夜视频国产福利| 少妇的逼好多水| 久久精品夜色国产| 国产高潮美女av| 毛片一级片免费看久久久久| 久久鲁丝午夜福利片| 国产一区二区亚洲精品在线观看| 夜夜夜夜夜久久久久| 高清午夜精品一区二区三区 | 国产精品久久久久久久久免| 自拍偷自拍亚洲精品老妇| 亚洲欧美日韩高清专用| 欧美区成人在线视频| 国产黄色小视频在线观看| 99久久精品一区二区三区| 高清日韩中文字幕在线| 亚洲人成网站在线播| av天堂在线播放| 老司机影院成人| 国语自产精品视频在线第100页| 中国国产av一级| 人妻丰满熟妇av一区二区三区| 亚洲美女视频黄频| 国产视频内射| 中文字幕av成人在线电影| 国产成人影院久久av| 久久久久国产网址| 日本爱情动作片www.在线观看 | 91久久精品国产一区二区成人| 亚洲精华国产精华液的使用体验 | 美女cb高潮喷水在线观看| 久久鲁丝午夜福利片| 亚洲最大成人av| 精品一区二区三区av网在线观看| 美女黄网站色视频| 色尼玛亚洲综合影院| 亚洲精品日韩av片在线观看| 国产视频内射| 日本五十路高清| 国产真实乱freesex| 国产视频一区二区在线看| av在线播放精品| 97超碰精品成人国产| 国产av不卡久久| 老师上课跳d突然被开到最大视频| 精品久久久久久久久亚洲| 51国产日韩欧美| 寂寞人妻少妇视频99o| 夜夜爽天天搞| 国语自产精品视频在线第100页| 在线观看美女被高潮喷水网站| 真人做人爱边吃奶动态| 久久精品综合一区二区三区| 真人做人爱边吃奶动态| 国产精品日韩av在线免费观看| 亚洲av成人精品一区久久| 日韩制服骚丝袜av| 欧美高清性xxxxhd video| 国产 一区精品| 少妇人妻精品综合一区二区 | 最好的美女福利视频网| 麻豆乱淫一区二区| 亚洲真实伦在线观看| 亚洲成人久久爱视频| 国内精品美女久久久久久| 插阴视频在线观看视频| 夜夜看夜夜爽夜夜摸| 美女内射精品一级片tv| 免费av不卡在线播放| 欧美zozozo另类| 国产精品人妻久久久久久| 天堂动漫精品| 成人综合一区亚洲| 国产欧美日韩精品亚洲av| 久久久久精品国产欧美久久久| 男人舔女人下体高潮全视频| 人妻久久中文字幕网| 免费看a级黄色片| 日韩 亚洲 欧美在线| 午夜福利在线观看免费完整高清在 | 男人和女人高潮做爰伦理| 国产亚洲精品久久久久久毛片| 99久久无色码亚洲精品果冻| 国产精品一及| 晚上一个人看的免费电影| 亚洲美女黄片视频| 亚洲精品456在线播放app| 午夜影院日韩av| 久久久久性生活片| 日本五十路高清| 三级毛片av免费| 国产爱豆传媒在线观看| 三级毛片av免费| 天天躁夜夜躁狠狠久久av| av女优亚洲男人天堂| 国内精品美女久久久久久| 国产精华一区二区三区| 亚洲第一区二区三区不卡| 男人的好看免费观看在线视频| 亚洲专区国产一区二区| 久久天躁狠狠躁夜夜2o2o| 网址你懂的国产日韩在线| 狠狠狠狠99中文字幕| 国产免费男女视频| 欧美zozozo另类| 黄色配什么色好看| 性插视频无遮挡在线免费观看| 少妇高潮的动态图| 久久午夜福利片| av天堂中文字幕网| 美女cb高潮喷水在线观看| 全区人妻精品视频| 久久久色成人| 中文字幕av成人在线电影| 能在线免费观看的黄片| 波多野结衣高清无吗| 日韩 亚洲 欧美在线| 国产乱人视频| 国产精品一区www在线观看| 亚洲欧美精品自产自拍| 成人性生交大片免费视频hd| 中文字幕av成人在线电影| 三级毛片av免费| 此物有八面人人有两片| 免费在线观看成人毛片| 亚洲乱码一区二区免费版| 国产精品一二三区在线看| av天堂在线播放| 桃色一区二区三区在线观看| 免费看a级黄色片| 最新在线观看一区二区三区| 亚洲国产精品sss在线观看| 久久草成人影院| 久久国产乱子免费精品| 免费黄网站久久成人精品| 看非洲黑人一级黄片| 中出人妻视频一区二区| 搡女人真爽免费视频火全软件 | 日本色播在线视频| 晚上一个人看的免费电影| 色噜噜av男人的天堂激情| 特大巨黑吊av在线直播| 男女做爰动态图高潮gif福利片| 午夜精品国产一区二区电影 | 亚洲内射少妇av| 国产欧美日韩精品一区二区| 综合色av麻豆| 特级一级黄色大片| av视频在线观看入口| 99久久九九国产精品国产免费| 毛片一级片免费看久久久久| 欧美日韩一区二区视频在线观看视频在线 | 我要搜黄色片| 成人高潮视频无遮挡免费网站| 91久久精品国产一区二区成人| 国产白丝娇喘喷水9色精品| videossex国产| 亚洲电影在线观看av| 最近在线观看免费完整版| 国产乱人视频| 男人狂女人下面高潮的视频| 男女做爰动态图高潮gif福利片| 天堂影院成人在线观看| 国产国拍精品亚洲av在线观看| 日日摸夜夜添夜夜添小说| 亚洲最大成人av| 欧美+日韩+精品| 少妇裸体淫交视频免费看高清| 欧美绝顶高潮抽搐喷水| 天天躁日日操中文字幕| 国产精品亚洲一级av第二区| 麻豆成人午夜福利视频| 2021天堂中文幕一二区在线观| 蜜桃亚洲精品一区二区三区| 草草在线视频免费看| 国产色爽女视频免费观看| 亚洲自偷自拍三级| 日韩欧美在线乱码| 久久6这里有精品| 一级黄色大片毛片| 国产片特级美女逼逼视频| 久久久久久九九精品二区国产| 韩国av在线不卡| 亚洲国产精品sss在线观看| 国产精品久久电影中文字幕| 久久精品国产清高在天天线| 日日摸夜夜添夜夜添小说| 欧美性猛交╳xxx乱大交人| 久久久久国产精品人妻aⅴ院| 国产精品久久久久久精品电影| 久久久久久久亚洲中文字幕| 淫秽高清视频在线观看| 亚洲不卡免费看| 国产精品久久久久久精品电影| 精品一区二区三区av网在线观看| 欧美性猛交黑人性爽| 久久久久免费精品人妻一区二区| 少妇猛男粗大的猛烈进出视频 | 国产精品国产三级国产av玫瑰| 在线播放国产精品三级| 国产黄a三级三级三级人| 人人妻人人看人人澡| 精品免费久久久久久久清纯| 在线a可以看的网站| 国产片特级美女逼逼视频| 蜜臀久久99精品久久宅男| 人妻久久中文字幕网| 99热网站在线观看| 欧美最黄视频在线播放免费| 久久久精品94久久精品| 日韩精品中文字幕看吧| 久久九九热精品免费| 日韩一本色道免费dvd| 性欧美人与动物交配| 久久久成人免费电影| 午夜精品一区二区三区免费看| 久久久久久久午夜电影| 又黄又爽又免费观看的视频| 给我免费播放毛片高清在线观看| 狂野欧美激情性xxxx在线观看| 久久久久久久久中文| 最近中文字幕高清免费大全6| 一级黄色大片毛片| 狂野欧美白嫩少妇大欣赏| 精华霜和精华液先用哪个| 国产精品永久免费网站| АⅤ资源中文在线天堂| av在线亚洲专区| 91久久精品国产一区二区成人| 九九在线视频观看精品| 1024手机看黄色片| 免费看a级黄色片| 成人亚洲精品av一区二区| 高清日韩中文字幕在线| 人人妻人人看人人澡| 91久久精品国产一区二区成人| 午夜福利在线观看免费完整高清在 | 一级av片app| 菩萨蛮人人尽说江南好唐韦庄 | 毛片女人毛片| 国产av在哪里看| 蜜桃久久精品国产亚洲av| 亚州av有码| 亚洲国产日韩欧美精品在线观看| 久久天躁狠狠躁夜夜2o2o| 国产爱豆传媒在线观看| 女的被弄到高潮叫床怎么办| 国产精品一区二区三区四区久久| 欧美人与善性xxx| 日韩高清综合在线| 婷婷精品国产亚洲av在线| 亚洲不卡免费看| 国产黄a三级三级三级人| 国产伦在线观看视频一区| 毛片一级片免费看久久久久| 国内少妇人妻偷人精品xxx网站| 亚洲欧美精品综合久久99| 亚洲精品色激情综合| 久久久久免费精品人妻一区二区| 美女内射精品一级片tv| 亚洲精品亚洲一区二区| 色av中文字幕| 亚洲熟妇中文字幕五十中出| 不卡视频在线观看欧美| 久久久色成人| 亚洲精品成人久久久久久| 国产成人a区在线观看| 97超级碰碰碰精品色视频在线观看| 国产不卡一卡二| 亚洲在线观看片| av在线天堂中文字幕| 人人妻人人澡人人爽人人夜夜 | 日本-黄色视频高清免费观看| 在线免费观看不下载黄p国产| 国产高清激情床上av| 乱码一卡2卡4卡精品| 91久久精品国产一区二区三区| 国产毛片a区久久久久| 亚洲国产色片| 日韩精品有码人妻一区| 国产欧美日韩精品一区二区| 久久韩国三级中文字幕| 91午夜精品亚洲一区二区三区| 亚洲美女搞黄在线观看 | 国产成人freesex在线 | 美女免费视频网站| 免费观看精品视频网站| 亚洲aⅴ乱码一区二区在线播放| 日本成人三级电影网站| 国产视频内射| av免费在线看不卡| 久久久欧美国产精品| 少妇猛男粗大的猛烈进出视频 | 我要看日韩黄色一级片| 中文亚洲av片在线观看爽| 麻豆av噜噜一区二区三区| 日韩欧美三级三区| 免费看光身美女| 久久午夜亚洲精品久久| 日产精品乱码卡一卡2卡三| av黄色大香蕉| 国产精品女同一区二区软件| 欧美日韩一区二区视频在线观看视频在线 | 久久九九热精品免费| 国产高清视频在线播放一区| 韩国av在线不卡| 久久精品国产自在天天线| 国产av麻豆久久久久久久| 91在线精品国自产拍蜜月| 三级男女做爰猛烈吃奶摸视频| 亚洲四区av| 国产午夜精品久久久久久一区二区三区 | 成人亚洲欧美一区二区av| a级毛色黄片| 亚洲精品色激情综合| 特级一级黄色大片| 在线免费十八禁| 18禁黄网站禁片免费观看直播| 99精品在免费线老司机午夜| 免费搜索国产男女视频| 日韩av不卡免费在线播放| 又爽又黄a免费视频| 日韩欧美免费精品| 国产一级毛片七仙女欲春2| 日本黄色片子视频| 天堂√8在线中文| 日韩精品中文字幕看吧| 国产中年淑女户外野战色| 久久中文看片网| 不卡视频在线观看欧美| 熟女人妻精品中文字幕| 久久久久久久久大av| 日韩欧美在线乱码| 给我免费播放毛片高清在线观看| 欧美bdsm另类| a级一级毛片免费在线观看| 午夜老司机福利剧场| 国产女主播在线喷水免费视频网站 | 三级国产精品欧美在线观看| 精品久久久久久成人av| a级毛片免费高清观看在线播放| av天堂中文字幕网| 久久久国产成人精品二区| 亚洲av成人av| 成年版毛片免费区| 亚洲成a人片在线一区二区| 黑人高潮一二区| 白带黄色成豆腐渣| 久久99热6这里只有精品| www.色视频.com| 日韩成人伦理影院| 免费观看的影片在线观看| 国产不卡一卡二| 在线播放无遮挡| 国产欧美日韩一区二区精品| 99久久精品热视频| 欧美绝顶高潮抽搐喷水| 国产麻豆成人av免费视频| 久久精品久久久久久噜噜老黄 | 日本一本二区三区精品| 成人午夜高清在线视频| 日韩精品青青久久久久久| 淫妇啪啪啪对白视频| 69av精品久久久久久| 此物有八面人人有两片| 十八禁国产超污无遮挡网站| 91久久精品国产一区二区成人|