• <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)背景更新方法
    乱码一卡2卡4卡精品| 亚洲精品影视一区二区三区av| 一边亲一边摸免费视频| 欧美性感艳星| 亚洲精品日本国产第一区| 国产av国产精品国产| 日韩大片免费观看网站| 国产亚洲av片在线观看秒播厂| 国产真实伦视频高清在线观看| 欧美变态另类bdsm刘玥| 亚洲国产成人一精品久久久| 麻豆成人午夜福利视频| 日本猛色少妇xxxxx猛交久久| 免费观看a级毛片全部| 亚洲欧美精品自产自拍| 男人添女人高潮全过程视频| 美女高潮的动态| 国产高清有码在线观看视频| 大码成人一级视频| 国产老妇女一区| 欧美人与善性xxx| 亚洲人成网站在线播| 成年人午夜在线观看视频| 国产在视频线精品| 亚洲综合精品二区| 午夜激情福利司机影院| 老女人水多毛片| 国产爽快片一区二区三区| 婷婷色av中文字幕| 亚洲精品第二区| 国产免费一区二区三区四区乱码| 美女xxoo啪啪120秒动态图| 日韩中字成人| 国产黄色免费在线视频| 男的添女的下面高潮视频| 免费观看av网站的网址| 国产午夜福利久久久久久| 黄色日韩在线| 成人午夜精彩视频在线观看| 午夜精品国产一区二区电影 | 啦啦啦中文免费视频观看日本| 精品久久久久久久久av| 成人国产av品久久久| 人妻制服诱惑在线中文字幕| 免费av观看视频| 欧美成人午夜免费资源| 少妇熟女欧美另类| 69人妻影院| 精品国产三级普通话版| 男男h啪啪无遮挡| 免费av毛片视频| 成人毛片a级毛片在线播放| 国产伦在线观看视频一区| 国产老妇女一区| 成人一区二区视频在线观看| 成人二区视频| 大香蕉97超碰在线| 69人妻影院| 丰满乱子伦码专区| www.色视频.com| 中文字幕av成人在线电影| 国产精品一区二区在线观看99| 日日摸夜夜添夜夜添av毛片| 国产淫语在线视频| 国产伦理片在线播放av一区| 可以在线观看毛片的网站| 欧美成人午夜免费资源| 久久精品国产亚洲网站| 亚洲va在线va天堂va国产| 亚洲aⅴ乱码一区二区在线播放| 日韩成人av中文字幕在线观看| av卡一久久| 国产免费一区二区三区四区乱码| 啦啦啦中文免费视频观看日本| 亚洲国产欧美人成| 在线观看三级黄色| 久久综合国产亚洲精品| 国产成人a∨麻豆精品| 色婷婷久久久亚洲欧美| av在线播放精品| 日韩三级伦理在线观看| 一级片'在线观看视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费av不卡在线播放| 亚洲av中文字字幕乱码综合| 涩涩av久久男人的天堂| 男人狂女人下面高潮的视频| 高清在线视频一区二区三区| 亚洲精品,欧美精品| 大片免费播放器 马上看| 国产精品av视频在线免费观看| 亚洲一级一片aⅴ在线观看| 亚洲伊人久久精品综合| 国产精品无大码| 最新中文字幕久久久久| 日本色播在线视频| 日本色播在线视频| 亚洲熟女精品中文字幕| 精华霜和精华液先用哪个| 永久免费av网站大全| 亚洲性久久影院| 天天一区二区日本电影三级| 成人美女网站在线观看视频| 久久精品国产鲁丝片午夜精品| 国产日韩欧美亚洲二区| 久热这里只有精品99| 黄色欧美视频在线观看| 日韩一区二区视频免费看| 亚洲国产精品国产精品| 久久精品国产a三级三级三级| 极品教师在线视频| 久久久久久久大尺度免费视频| 国产精品99久久99久久久不卡 | 亚洲四区av| 性色avwww在线观看| 97热精品久久久久久| 久久韩国三级中文字幕| videossex国产| 在线 av 中文字幕| 亚洲精品第二区| 一级毛片我不卡| 美女cb高潮喷水在线观看| 91在线精品国自产拍蜜月| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美少妇被猛烈插入视频| 少妇人妻久久综合中文| 国产亚洲91精品色在线| av一本久久久久| 卡戴珊不雅视频在线播放| 下体分泌物呈黄色| 欧美3d第一页| 秋霞在线观看毛片| 乱码一卡2卡4卡精品| 午夜免费男女啪啪视频观看| 少妇高潮的动态图| 成年女人看的毛片在线观看| 少妇高潮的动态图| 亚洲怡红院男人天堂| 别揉我奶头 嗯啊视频| 男插女下体视频免费在线播放| 一级片'在线观看视频| 国产亚洲最大av| 久久久久久久精品精品| 国产精品嫩草影院av在线观看| 视频中文字幕在线观看| 插阴视频在线观看视频| 欧美日韩国产mv在线观看视频 | 免费不卡的大黄色大毛片视频在线观看| 精品午夜福利在线看| 日本爱情动作片www.在线观看| 国产黄片视频在线免费观看| 国产在视频线精品| 如何舔出高潮| 青春草视频在线免费观看| 成人午夜精彩视频在线观看| 大香蕉久久网| 国产免费福利视频在线观看| 91午夜精品亚洲一区二区三区| 欧美日韩视频精品一区| 日本午夜av视频| 又大又黄又爽视频免费| 高清日韩中文字幕在线| 91久久精品国产一区二区三区| 精品一区二区免费观看| 国产伦在线观看视频一区| 人妻夜夜爽99麻豆av| 亚洲av日韩在线播放| 精品国产一区二区三区久久久樱花 | 看十八女毛片水多多多| 插逼视频在线观看| 人体艺术视频欧美日本| 日本黄大片高清| 日本免费在线观看一区| 午夜福利网站1000一区二区三区| 欧美三级亚洲精品| av.在线天堂| 22中文网久久字幕| freevideosex欧美| 少妇人妻一区二区三区视频| 亚洲精品久久久久久婷婷小说| 成人美女网站在线观看视频| 成人美女网站在线观看视频| 亚洲最大成人中文| 国产一区二区三区综合在线观看 | 看黄色毛片网站| 色综合色国产| 国产免费一级a男人的天堂| 亚洲国产精品专区欧美| 观看美女的网站| 最后的刺客免费高清国语| 热re99久久精品国产66热6| 精品一区二区三卡| 亚洲av二区三区四区| 日本与韩国留学比较| eeuss影院久久| 97在线视频观看| 国精品久久久久久国模美| 国产亚洲av嫩草精品影院| 国产精品蜜桃在线观看| 联通29元200g的流量卡| 全区人妻精品视频| 国产黄色免费在线视频| 亚洲国产欧美在线一区| 亚洲欧美一区二区三区黑人 | 国产探花在线观看一区二区| 精品久久久噜噜| 另类亚洲欧美激情| 高清视频免费观看一区二区| 99久久九九国产精品国产免费| 国产高清不卡午夜福利| 男的添女的下面高潮视频| 91aial.com中文字幕在线观看| 亚洲精品日韩在线中文字幕| 国产伦精品一区二区三区视频9| 99久久中文字幕三级久久日本| xxx大片免费视频| 国产精品爽爽va在线观看网站| 欧美日韩视频高清一区二区三区二| 黄色怎么调成土黄色| 性色avwww在线观看| av女优亚洲男人天堂| 国产一级毛片在线| 国产 一区 欧美 日韩| 国产精品麻豆人妻色哟哟久久| 国产乱人偷精品视频| 男女无遮挡免费网站观看| 欧美激情国产日韩精品一区| .国产精品久久| 欧美日韩在线观看h| 视频区图区小说| 久久久精品94久久精品| 性插视频无遮挡在线免费观看| 亚洲国产精品成人久久小说| 精品99又大又爽又粗少妇毛片| 直男gayav资源| 国产精品不卡视频一区二区| 国产熟女欧美一区二区| 国产永久视频网站| 搡女人真爽免费视频火全软件| 国产av不卡久久| 最近中文字幕高清免费大全6| 国产精品嫩草影院av在线观看| 精品人妻熟女av久视频| 欧美 日韩 精品 国产| 国产色爽女视频免费观看| 亚洲图色成人| 免费av不卡在线播放| 婷婷色麻豆天堂久久| 国产成人精品婷婷| 国产成人a∨麻豆精品| 国产亚洲一区二区精品| 男女无遮挡免费网站观看| 久久久久久久久久人人人人人人| 男女啪啪激烈高潮av片| 亚洲精品色激情综合| 18禁在线播放成人免费| 精品一区二区三区视频在线| 国产91av在线免费观看| 日韩电影二区| 精品国产露脸久久av麻豆| 日本与韩国留学比较| 天美传媒精品一区二区| 狂野欧美白嫩少妇大欣赏| 美女内射精品一级片tv| 欧美日韩在线观看h| 我要看日韩黄色一级片| 在线免费观看不下载黄p国产| 水蜜桃什么品种好| 国产大屁股一区二区在线视频| 国产又色又爽无遮挡免| 91久久精品电影网| 国产成年人精品一区二区| 免费观看无遮挡的男女| 中文字幕制服av| 日韩在线高清观看一区二区三区| 国产午夜福利久久久久久| 女人十人毛片免费观看3o分钟| 最近的中文字幕免费完整| 观看免费一级毛片| 欧美高清成人免费视频www| 亚洲欧美精品自产自拍| 国产精品秋霞免费鲁丝片| 久久久久九九精品影院| 国产免费又黄又爽又色| 黄色日韩在线| 亚洲欧美一区二区三区黑人 | 亚洲欧美精品专区久久| 欧美精品国产亚洲| 亚洲国产精品成人综合色| 99久国产av精品国产电影| 熟妇人妻不卡中文字幕| 在线播放无遮挡| 午夜福利网站1000一区二区三区| av国产免费在线观看| 午夜激情久久久久久久| 卡戴珊不雅视频在线播放| 人妻系列 视频| 午夜精品一区二区三区免费看| 不卡视频在线观看欧美| 亚洲激情五月婷婷啪啪| 亚州av有码| 一区二区三区免费毛片| 久热久热在线精品观看| 成人午夜精彩视频在线观看| 日本免费在线观看一区| 免费播放大片免费观看视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲,欧美,日韩| 国产精品一及| 精品少妇久久久久久888优播| 水蜜桃什么品种好| 午夜视频国产福利| 免费看a级黄色片| 综合色丁香网| 在线亚洲精品国产二区图片欧美 | 色综合色国产| 我的女老师完整版在线观看| 大香蕉97超碰在线| 联通29元200g的流量卡| 亚洲精品视频女| 伦理电影大哥的女人| 国产色爽女视频免费观看| 欧美xxⅹ黑人| 日韩伦理黄色片| 亚洲aⅴ乱码一区二区在线播放| 一本一本综合久久| 亚洲三级黄色毛片| 国产亚洲5aaaaa淫片| 欧美少妇被猛烈插入视频| 亚洲精品第二区| 亚洲性久久影院| 深夜a级毛片| 一个人观看的视频www高清免费观看| 女人被狂操c到高潮| 黄色配什么色好看| 久久精品国产自在天天线| 黄色怎么调成土黄色| 建设人人有责人人尽责人人享有的 | 精品国产一区二区三区久久久樱花 | 国产精品嫩草影院av在线观看| 国产一级毛片在线| 肉色欧美久久久久久久蜜桃 | 乱系列少妇在线播放| 成人欧美大片| 亚洲欧美一区二区三区黑人 | 国产成年人精品一区二区| 丝袜美腿在线中文| 国产91av在线免费观看| 赤兔流量卡办理| 精品人妻偷拍中文字幕| 99久久精品热视频| 大码成人一级视频| 免费大片18禁| 国产精品久久久久久精品电影小说 | 免费高清在线观看视频在线观看| 欧美另类一区| 国产精品99久久99久久久不卡 | 中国国产av一级| 国产精品伦人一区二区| 久久久久久久久久人人人人人人| 男女下面进入的视频免费午夜| 丝袜美腿在线中文| 乱码一卡2卡4卡精品| 国内精品美女久久久久久| 在线观看一区二区三区| 丰满少妇做爰视频| 真实男女啪啪啪动态图| av福利片在线观看| 亚洲国产欧美在线一区| 在线免费观看不下载黄p国产| 波多野结衣巨乳人妻| 国内少妇人妻偷人精品xxx网站| 黑人高潮一二区| 国产精品嫩草影院av在线观看| av专区在线播放| 午夜福利高清视频| 97在线视频观看| 女的被弄到高潮叫床怎么办| 免费观看av网站的网址| 亚洲精品国产成人久久av| 欧美成人午夜免费资源| av.在线天堂| 丰满人妻一区二区三区视频av| 一二三四中文在线观看免费高清| 久热久热在线精品观看| 美女脱内裤让男人舔精品视频| 欧美性感艳星| 草草在线视频免费看| 久久精品久久精品一区二区三区| 国产精品国产三级国产专区5o| 国产午夜精品一二区理论片| av免费观看日本| 一级毛片我不卡| 国产在线男女| 欧美日韩精品成人综合77777| 最近中文字幕高清免费大全6| 美女视频免费永久观看网站| 欧美激情久久久久久爽电影| 美女主播在线视频| 欧美xxxx黑人xx丫x性爽| 久久精品国产a三级三级三级| 国产黄片美女视频| 99久久精品热视频| 天美传媒精品一区二区| 九九爱精品视频在线观看| 老司机影院毛片| 亚洲av免费在线观看| 禁无遮挡网站| 午夜精品一区二区三区免费看| 国产免费福利视频在线观看| 欧美另类一区| 久热久热在线精品观看| 成年版毛片免费区| 国产成年人精品一区二区| 18禁在线无遮挡免费观看视频| 久久精品熟女亚洲av麻豆精品| 日本与韩国留学比较| av在线播放精品| 在线观看一区二区三区激情| av播播在线观看一区| 春色校园在线视频观看| 亚洲国产精品999| 成人国产av品久久久| 亚洲人与动物交配视频| 日韩国内少妇激情av| 国产男女内射视频| 亚洲精品456在线播放app| 日韩,欧美,国产一区二区三区| 国产精品国产三级专区第一集| 久久综合国产亚洲精品| 亚洲天堂国产精品一区在线| 寂寞人妻少妇视频99o| 视频中文字幕在线观看| av黄色大香蕉| 国产成人精品久久久久久| 一级毛片久久久久久久久女| 国产日韩欧美亚洲二区| 亚洲精品视频女| 久久午夜福利片| 亚洲人成网站在线观看播放| 亚洲欧美成人综合另类久久久| 性色avwww在线观看| 亚州av有码| 日本黄色片子视频| .国产精品久久| 97在线视频观看| 日本三级黄在线观看| 99九九线精品视频在线观看视频| 国产中年淑女户外野战色| 五月开心婷婷网| 亚洲内射少妇av| 晚上一个人看的免费电影| 国产午夜福利久久久久久| 夫妻午夜视频| 丝瓜视频免费看黄片| 日韩一本色道免费dvd| 性插视频无遮挡在线免费观看| 国产成人a区在线观看| 国产真实伦视频高清在线观看| 在线观看免费高清a一片| 一级毛片aaaaaa免费看小| 黄片wwwwww| 狂野欧美激情性bbbbbb| 99精国产麻豆久久婷婷| 国产精品精品国产色婷婷| 久久久久久久久久久丰满| 午夜日本视频在线| 欧美xxxx黑人xx丫x性爽| 日韩欧美一区视频在线观看 | 毛片女人毛片| 欧美变态另类bdsm刘玥| 99热网站在线观看| 蜜桃久久精品国产亚洲av| 中文字幕亚洲精品专区| 亚洲欧美一区二区三区黑人 | 欧美xxxx黑人xx丫x性爽| 久久这里有精品视频免费| 国产成人91sexporn| 久久久久久久大尺度免费视频| 1000部很黄的大片| 国产成人免费无遮挡视频| 在线精品无人区一区二区三 | 日日撸夜夜添| 久久久久久久大尺度免费视频| 交换朋友夫妻互换小说| 一本久久精品| 国产毛片a区久久久久| 成人特级av手机在线观看| 亚洲精品一区蜜桃| 国产 一区精品| 国产中年淑女户外野战色| 身体一侧抽搐| 日本欧美国产在线视频| 免费观看av网站的网址| videossex国产| 少妇高潮的动态图| 国产精品女同一区二区软件| 久久久精品免费免费高清| 大码成人一级视频| av在线app专区| 成人国产麻豆网| 亚洲精品456在线播放app| 只有这里有精品99| 下体分泌物呈黄色| 国产精品久久久久久久电影| 亚洲国产精品成人综合色| 欧美变态另类bdsm刘玥| 国产精品三级大全| 精品一区二区三卡| av又黄又爽大尺度在线免费看| 欧美丝袜亚洲另类| 亚洲精品日韩在线中文字幕| 丝袜脚勾引网站| 国产高潮美女av| 亚洲av免费高清在线观看| 亚洲成色77777| 女人十人毛片免费观看3o分钟| 欧美高清成人免费视频www| 国产精品人妻久久久影院| 亚洲国产精品国产精品| 国内揄拍国产精品人妻在线| 亚洲精品第二区| 日韩一本色道免费dvd| 丝袜脚勾引网站| 人妻夜夜爽99麻豆av| 成人无遮挡网站| 久久久久国产网址| 亚洲,一卡二卡三卡| 日日摸夜夜添夜夜添av毛片| 男女边吃奶边做爰视频| 国产一区有黄有色的免费视频| 国产一区二区在线观看日韩| 超碰97精品在线观看| 国产亚洲5aaaaa淫片| 麻豆精品久久久久久蜜桃| 91午夜精品亚洲一区二区三区| 国产亚洲av嫩草精品影院| 国产精品.久久久| 91在线精品国自产拍蜜月| 国产成人一区二区在线| 久久影院123| 国产日韩欧美亚洲二区| 2021天堂中文幕一二区在线观| 日韩免费高清中文字幕av| 女人久久www免费人成看片| 午夜亚洲福利在线播放| 国产伦理片在线播放av一区| 男女国产视频网站| 国产午夜精品一二区理论片| 九色成人免费人妻av| 欧美极品一区二区三区四区| 天天躁夜夜躁狠狠久久av| 日本欧美国产在线视频| 亚洲欧美一区二区三区黑人 | 国产 精品1| 亚洲精品视频女| 99热这里只有是精品在线观看| 日韩av不卡免费在线播放| 国产片特级美女逼逼视频| 午夜福利网站1000一区二区三区| 美女内射精品一级片tv| 三级经典国产精品| 国产成人免费无遮挡视频| 九草在线视频观看| 亚洲,一卡二卡三卡| 伦理电影大哥的女人| 亚洲精品影视一区二区三区av| 内射极品少妇av片p| 91在线精品国自产拍蜜月| 久久鲁丝午夜福利片| 亚洲av一区综合| 欧美+日韩+精品| 超碰av人人做人人爽久久| 麻豆成人av视频| 久久精品久久精品一区二区三区| 欧美老熟妇乱子伦牲交| 久久久精品欧美日韩精品| 伦精品一区二区三区| 18禁在线播放成人免费| 在线精品无人区一区二区三 | 国产黄a三级三级三级人| 日韩一区二区视频免费看| 麻豆国产97在线/欧美| 男人添女人高潮全过程视频| 少妇的逼好多水| a级一级毛片免费在线观看| 综合色av麻豆| 午夜免费观看性视频| 欧美少妇被猛烈插入视频| 亚洲精品一二三| 精品少妇久久久久久888优播| 欧美老熟妇乱子伦牲交| 国产精品精品国产色婷婷| 97超碰精品成人国产| 99九九线精品视频在线观看视频| 赤兔流量卡办理| 日韩强制内射视频| 国产伦理片在线播放av一区| 天美传媒精品一区二区| h日本视频在线播放| 男男h啪啪无遮挡| 夜夜看夜夜爽夜夜摸| 一区二区三区乱码不卡18| 国产高潮美女av| 亚洲欧美日韩另类电影网站 | 丝袜喷水一区| 国产成人精品婷婷| 精品久久久久久久人妻蜜臀av| 夫妻午夜视频| 麻豆成人av视频| 亚洲高清免费不卡视频| 亚洲av欧美aⅴ国产| 国产视频内射| 少妇人妻一区二区三区视频| 午夜精品国产一区二区电影 |