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

    姿態(tài)解算與外力加速度同步估計(jì)算法

    2016-05-14 11:07:05孟唐宇浦劍濤方建軍梁嵐珍
    計(jì)算機(jī)應(yīng)用 2016年5期

    孟唐宇 浦劍濤 方建軍 梁嵐珍

    摘要:針對(duì)慣性導(dǎo)航應(yīng)用中,姿態(tài)解算與外力加速度估計(jì)互相干擾的問題,提出一種基于四元數(shù)和擴(kuò)展卡爾曼濾波器的姿態(tài)解算與外力加速度同步估計(jì)算法。首先,利用估計(jì)的外力加速度修正傳感器加速度數(shù)據(jù)得到準(zhǔn)確的反向重力加速度,再結(jié)合地磁場(chǎng)向量通過梯度下降算法解算得到旋轉(zhuǎn)四元數(shù)的測(cè)量值;其次,構(gòu)建擴(kuò)展卡爾曼濾波模型,對(duì)旋轉(zhuǎn)四元數(shù)和外力加速度進(jìn)行更新,得到旋轉(zhuǎn)四元數(shù)的預(yù)測(cè)值和外力加速度的預(yù)測(cè)值;最后,用旋轉(zhuǎn)四元數(shù)的測(cè)量值和測(cè)量得到的加速度數(shù)據(jù)對(duì)預(yù)測(cè)值通過擴(kuò)展卡爾曼濾波的方法進(jìn)行校正,最終得到準(zhǔn)確的旋轉(zhuǎn)四元數(shù)和參考坐標(biāo)系下三軸方向上的外力加速度。實(shí)驗(yàn)表明,通過擴(kuò)展卡爾曼濾波同時(shí)對(duì)姿態(tài)和外力加速度估計(jì)的方法,能夠迅速收斂并準(zhǔn)確得機(jī)體姿態(tài)信息以及外力加速度信息,歐拉角誤差為±1.95°,加速度誤差為±0.12m/s2,并且該算法能有效抑制外力加速度對(duì)姿態(tài)解算的影響,準(zhǔn)確估計(jì)外力加速度。

    關(guān)鍵詞:姿態(tài)解算;四元數(shù);擴(kuò)展卡爾曼濾波;梯度下降法;外力加速度

    中圖分類號(hào):TP273+.5 文獻(xiàn)標(biāo)志碼:A

    Abstract:Aiming at the problem of mutual interference between attitude algorithm and external force acceleration estimation in inertial navigation system, a new method based on quaternion and extended Kalman filter was proposed. Firstly, the acceleration data of the sensor was corrected by using the estimated external force acceleration data to obtain the accurate reverse gravity acceleration, combined with geomagnetic field vector and calculated by the gradient descent algorithm, the rotate quaternions were obtained. Secondly, the extended Kalman filter model was constructed to update the rotate quaternions and external force acceleration, the prediction value of rotate quaternions and the external force were obtained. Finally, the measured values of rotate quaternions and the acceleration data were corrected by Kalman filtering method, the accurate rotate quaternions and the external force acceleration of the three axis directions in reference coordinate system were obtained. The experimental results show that the method for the synchronization estimation of attitude and external force acceleration by extended Calman filter can quickly converge and accurately get the information of the attitude and the external force acceleration, its Euler angle error is ±1.95° and acceleration error is ±0.12m/s2. The method can effectively restrain the influence of the external force acceleration on the attitude algorithm, and accurately estimate the external force.

    Key words:attitude algorithm; quaternion; extended Kalman filter; gradient descent algorithm; external force acceleration

    0 引言

    在慣性導(dǎo)航應(yīng)用中,如機(jī)器人、汽車產(chǎn)業(yè)和可穿戴設(shè)備等領(lǐng)域,常要求運(yùn)動(dòng)著的機(jī)體能夠確定其自身當(dāng)前的姿態(tài),以及其外力加速度數(shù)據(jù),來估算出機(jī)體的位置信息。目前常采用的方法是,用慣性測(cè)量傳感器測(cè)量機(jī)體的慣性數(shù)據(jù),通過算法解算出當(dāng)前載體的姿態(tài),再根據(jù)姿態(tài)信息解算得出外力加速度。

    當(dāng)機(jī)體加速移動(dòng)時(shí),外力加速度的存在會(huì)嚴(yán)重影響姿態(tài)解算的結(jié)果,當(dāng)外力加速度越大時(shí),求解的姿態(tài)角與實(shí)際值偏差越大。姿態(tài)解算算法主要有歐拉角法[1]、方向余弦法[2-3]和四元數(shù)法[4]等。歐拉角法是采用對(duì)反向重力向量的反三角函數(shù)求解的方法,雖然計(jì)算量小但存在萬向節(jié)死鎖問題;方向余弦法是通過旋轉(zhuǎn)矩陣對(duì)載體坐標(biāo)系進(jìn)行轉(zhuǎn)換,計(jì)算過程中要對(duì)旋轉(zhuǎn)矩陣中的9個(gè)數(shù)據(jù)和傳感器數(shù)據(jù)進(jìn)行計(jì)算,所以計(jì)算量通常較大;四元數(shù)法相較于方向余弦法,采用四維數(shù)據(jù)進(jìn)行計(jì)算,計(jì)算量大幅下降。在算法的運(yùn)用中Wang等[5]運(yùn)用四階龍格庫塔法和旋轉(zhuǎn)矩陣的方法進(jìn)行計(jì)算,并利用每次的計(jì)算結(jié)果來校正傳感器的誤差參數(shù)。這樣雖然減小了計(jì)算量但是忽略了外力加速度對(duì)姿態(tài)的影響,實(shí)驗(yàn)結(jié)果可能會(huì)受到外力加速度的影響。Madgwick[6]運(yùn)用梯度下降法解算載體姿態(tài),這樣計(jì)算在靜態(tài)時(shí)準(zhǔn)確度較高,但在動(dòng)態(tài)時(shí)會(huì)產(chǎn)生姿態(tài)角的大幅抖動(dòng),并且過度依賴加速度計(jì)數(shù)據(jù),使得外力加速度對(duì)計(jì)算姿態(tài)的影響更加明顯。在算法的運(yùn)用中如成怡等[7]利用非線性濾波,閻世梁等[8]運(yùn)用互補(bǔ)濾波, Bachmann等[9]卡爾曼濾波,Madgwick等[10]運(yùn)用梯度下降法, Marins等[11]利用擴(kuò)展卡爾曼濾波算法。以上算法中為了消除外力加速度對(duì)姿態(tài)的影響,均將加速計(jì)傳感器采集到的數(shù)據(jù)進(jìn)行歸一化,再進(jìn)行姿態(tài)解算,這樣可以減小外力加速度對(duì)反向重力向量的影響,但仍然沒有從根本上消除外力加速度的影響,反向重力向量估計(jì)依然是不準(zhǔn)確的。本文在姿態(tài)解算過程中,利用擴(kuò)展卡爾曼濾波器模型,對(duì)姿態(tài)與外力加速度同步進(jìn)行估計(jì),不僅在姿態(tài)解算時(shí)消除了外力加速度對(duì)姿態(tài)解算精度的影響,同時(shí)可以得到外力加速度信息。

    在外力加速度估計(jì)算法運(yùn)用中張?zhí)鸬萚12]利用全球定位系統(tǒng)(Global Positioning System, GPS)基站方式進(jìn)行加速度估計(jì),如果在室內(nèi)、地下室等GPS信號(hào)弱的環(huán)境下無法測(cè)量出機(jī)體的運(yùn)動(dòng)加速度。李海濤[13]利用電子羅盤修正加速度計(jì)信號(hào)的方法對(duì)機(jī)體外力加速度進(jìn)行估計(jì),然而電子羅盤信號(hào)容易受到外界電磁干擾,所得外力加速度會(huì)受到一定影響,因此姿態(tài)解算也會(huì)產(chǎn)生一定誤差。

    本研究采用基于四元數(shù)和擴(kuò)展卡爾曼濾波器對(duì)機(jī)體姿態(tài)和外力加速度同時(shí)估計(jì)。利用外力加速度的估計(jì)值對(duì)采集到的加速度數(shù)據(jù)進(jìn)行修正,得到修正后的反向重力加速度再結(jié)合磁場(chǎng)向量,采用梯度下降法估計(jì)出四元數(shù)測(cè)量值與加速度測(cè)量值;在此基礎(chǔ)上,用擴(kuò)展卡爾曼濾波器對(duì)預(yù)測(cè)四元數(shù)和預(yù)測(cè)外力加速度值同時(shí)進(jìn)行修正,從而得到當(dāng)前狀態(tài)下機(jī)體準(zhǔn)確的旋轉(zhuǎn)四元數(shù)和參考坐標(biāo)系下的外力加速度。

    1 四元數(shù)介紹及基本運(yùn)算

    四元數(shù)應(yīng)用在姿態(tài)解算中不僅計(jì)算量小還可以消除萬向節(jié)死鎖問題,并具有能夠全角度解析姿態(tài)角等優(yōu)點(diǎn), 因此采用四元數(shù)作為姿態(tài)解算的數(shù)學(xué)推導(dǎo)演算依據(jù)。

    1.1 四元數(shù)表示形式

    對(duì)三角函數(shù)的±180°反解,得到姿態(tài)歐拉角的全角度解析,所得四元數(shù)為當(dāng)前載體坐標(biāo)系旋轉(zhuǎn)到參考坐標(biāo)系的四元數(shù)。從而有效解決歐拉角法產(chǎn)生的萬向節(jié)死鎖問題,同時(shí)減小運(yùn)算量,并且提高運(yùn)算速度。

    2 擴(kuò)展卡爾曼濾波姿態(tài)解算

    為了預(yù)測(cè)非線性四元數(shù)和外力加速度值,引入擴(kuò)展卡爾曼濾波器[11]。擴(kuò)展卡爾曼濾波器能夠預(yù)測(cè)出當(dāng)前載體坐標(biāo)系旋轉(zhuǎn)到參考坐標(biāo)系的旋轉(zhuǎn)四元數(shù)和外力加速度,并通過測(cè)量值進(jìn)行修正從而得到當(dāng)前狀態(tài)下準(zhǔn)確的旋轉(zhuǎn)四元數(shù)和外力加速度。這樣計(jì)算得出的四元數(shù)不受外力加速度影響,計(jì)算出的姿態(tài)角在機(jī)體非勻速運(yùn)動(dòng)時(shí)波動(dòng)小、外力加速度信息準(zhǔn)確,同時(shí)分解出的外力加速度可用于慣性導(dǎo)航的研究。

    擴(kuò)展卡爾曼濾波姿態(tài)解算基本流程如圖1所示。首先,利用陀螺儀數(shù)據(jù)更新旋轉(zhuǎn)四元數(shù),并且更新外力加速度數(shù)據(jù);其次,估計(jì)預(yù)測(cè)協(xié)方差Pk以及擴(kuò)展卡爾曼濾波增益K。再次,根據(jù)更新得到的旋轉(zhuǎn)四元數(shù)計(jì)算載體坐標(biāo)系下的反向重力向量;最后,根據(jù)測(cè)量值得到狀態(tài)變量的最優(yōu)估計(jì)并且更新協(xié)方差。其中xk表示狀態(tài)變量,k表示當(dāng)前狀態(tài),k-1表示上一時(shí)刻狀態(tài);Pk-1表示上一時(shí)刻協(xié)方差的值。

    當(dāng)前所得四元數(shù),即為載體轉(zhuǎn)向參考坐標(biāo)系的旋轉(zhuǎn)四元數(shù),所得加速度為載體坐標(biāo)系中三軸方向的加速度。通過四元數(shù)相乘的方法,把所得外力加速度轉(zhuǎn)換成參考坐標(biāo)系中三軸加速度,為實(shí)現(xiàn)慣性定位提供有效的加速度信息。

    3 梯度下降法姿態(tài)解算

    實(shí)際應(yīng)用表明,僅把加速計(jì)采集到的數(shù)據(jù)作為反向重力向量,這樣數(shù)據(jù)中隱含的外力加速度會(huì)影響旋轉(zhuǎn)四元數(shù)測(cè)量值求解的準(zhǔn)確度。為解決這一問題,先利用外力加速度的估計(jì)值修正加速計(jì)數(shù)據(jù)以得到準(zhǔn)確的反向重力向量,再根據(jù)修正后的反向重力向量與磁場(chǎng)向量解算出測(cè)量四元數(shù),從而減少外力加速度對(duì)旋轉(zhuǎn)四元數(shù)的影響。

    3.1 利用修正后的反向重力向量解算橫滾角和俯仰角

    4 算法融合及實(shí)驗(yàn)結(jié)果分析

    4.1 梯度下降法和擴(kuò)展卡爾曼濾波融合

    融合算法基本流程如圖2所示。

    梯度下降法循環(huán)計(jì)算2~3次,得到的四元數(shù)和加速計(jì)測(cè)值作為擴(kuò)展卡爾曼濾波的測(cè)量值。融合算法結(jié)構(gòu)如圖3所示,其中Mag代表電子羅盤傳感器所測(cè)量的磁場(chǎng)向量,Acc代表加速度傳感器所測(cè)量的加速度向量,Gyr代表陀螺儀傳感器所測(cè)量的角速度向量,acc代表計(jì)算得出的參考坐標(biāo)系下機(jī)體外力加速度向量。

    4.2 實(shí)驗(yàn)結(jié)果分析

    實(shí)驗(yàn)MCU選用MEGA2560,傳感器選用MPU6050、HMC5883L實(shí)驗(yàn)平臺(tái),MEGA2560以應(yīng)答方式通過串口把九軸數(shù)據(jù)傳送給Matlab進(jìn)行實(shí)驗(yàn)仿真。傳感器坐標(biāo)系即為機(jī)體坐標(biāo)系。

    通過三種方式驗(yàn)證算法可行性和準(zhǔn)確度。

    1)機(jī)體處于水平位置。機(jī)體姿態(tài)模型如圖4所示。所建模型和實(shí)物姿態(tài)一致。

    機(jī)體處于水平位置未加擴(kuò)展卡爾曼濾波與加上擴(kuò)展卡爾曼濾波歐拉角對(duì)比曲線如圖5所示,Pitch、Roll、Yaw角度在加入擴(kuò)展卡爾曼濾波前后均穩(wěn)定在0°,波動(dòng)幅度±1.95°,但加入擴(kuò)展卡爾曼濾波算法后角度波動(dòng)幅度比未加入擴(kuò)展卡爾曼濾波算法波動(dòng)幅度有所減小。歐拉角波動(dòng)是因?yàn)閼T性傳感器數(shù)據(jù)存在噪聲引起,經(jīng)過算法過濾后波動(dòng)幅度減少了10%。在靜止?fàn)顟B(tài)下兩種算法解算出的歐拉角差別不大主要是應(yīng)為沒有外力加速度干擾。

    機(jī)體向Y正方向以0.6m/s2勻加速運(yùn)動(dòng),加入擴(kuò)展卡爾曼濾波算法和未加入擴(kuò)展卡爾曼濾波算法歐拉角對(duì)比曲線如圖8所示,當(dāng)機(jī)體做勻加速運(yùn)動(dòng)時(shí),歐拉角曲線與機(jī)體靜止時(shí)相比,加入擴(kuò)展卡爾曼濾波算法的Roll有所減小,但范圍在±0.5°,而沒有加入擴(kuò)展卡爾曼濾波算法的Roll減小了3°。外力加速度對(duì)加入擴(kuò)展卡爾曼濾波算法姿態(tài)解算的結(jié)果影響很小而對(duì)未加入擴(kuò)展卡爾曼濾波算法姿態(tài)解算影響較大。因機(jī)體處于水平位置,只在Y軸有一恒定加速度,Pitch和Yaw角度均穩(wěn)定在0°,Yaw會(huì)隨運(yùn)動(dòng)方向的誤差有1°的變化波動(dòng)。

    3)機(jī)體在自由運(yùn)動(dòng)時(shí),加入擴(kuò)展卡爾曼濾波算法和未加入擴(kuò)展卡爾曼濾波算法歐拉角對(duì)比曲線如圖10,當(dāng)機(jī)體處于自由運(yùn)動(dòng)時(shí),加入擴(kuò)展卡爾曼濾波算法與未加入擴(kuò)展卡爾曼濾波算法姿態(tài)角曲線Pitch、Roll、Yaw角曲線趨勢(shì)基本相同,但未加入濾波器的算法解算出的歐拉角均偏大,說明機(jī)體在自由運(yùn)動(dòng)時(shí)所產(chǎn)生的外力加速度對(duì)未加入濾波器算法的姿態(tài)解算產(chǎn)生較大影響,誤差為±4.3°,而對(duì)加入濾波器算法的姿態(tài)解算影響較小,誤差為±1.94°。

    5 結(jié)語

    本文在姿態(tài)解算與外力加速度同步估計(jì)中,基于四元數(shù)和擴(kuò)展卡爾曼濾波器。通過對(duì)外力加速度的估計(jì),修正反向重力向量,從而解算出旋轉(zhuǎn)四元數(shù)的測(cè)量值,并利用卡爾曼濾波器對(duì)旋轉(zhuǎn)四元數(shù)與外力加速度進(jìn)行更新。這樣減小了外力速度對(duì)姿態(tài)的影響,并且準(zhǔn)確計(jì)算出姿態(tài)四元數(shù)和外力加速度。最后通過傳感器采集數(shù)據(jù)解算當(dāng)前傳感器的姿態(tài)和加速度,驗(yàn)證了該方法的可靠性。在Matlab中姿態(tài)更新時(shí)間間隔在15ms左右,且采用應(yīng)答方式獲取傳感器信息這會(huì)產(chǎn)生一定時(shí)間損耗。將算法移植到嵌入式系統(tǒng)中進(jìn)行運(yùn)算,這樣不僅最大限度地降低數(shù)據(jù)在傳輸過程中產(chǎn)生的不必要的時(shí)間損耗,而且運(yùn)算速度將大幅提升?,F(xiàn)該算法已在室內(nèi)定位小車項(xiàng)目中應(yīng)用。

    目前的研究比較集中于姿態(tài)的解算,尚未將外力加速度用于慣性導(dǎo)航中,這將是后續(xù)工作中要研究的一個(gè)重點(diǎn)工作。

    參考文獻(xiàn):

    [1]SUSHKOV O. Robot localisation using a distributed multimodal Kalman filter, and friends[D]. Sydney: University of New South Wales, 2006.

    [2]PREMERLANI W, BIZARD P. Direction cosine matrix IMU: theory[EB/OL]. [2015-04-10]. http://gentlenav.googlecode.com/files/DCMDraft2.pdf

    [3]周浩,浦劍濤,梁嵐珍,等. 基于體感的仿人機(jī)器人步態(tài)學(xué)習(xí)與控制[J]. 計(jì)算機(jī)應(yīng)用,2015,35(3):787-791.(ZHOU H, PU J T, LIANG L Z, et al. Gait learning and control of humanoid robot based on Kinect[J]. Journal of Computer Applications, 2015, 35(3):787-791.)

    [4]KUIPERS J B. Quaternions and Rotation Sequences[M]. Princeton: Princeton University Press, 1999: 127-142.

    [5]WANG Y, LI N, CHEN X, et al. Design and implementation of an AHRS based on MEMS sensors and complementary filtering[J]. Advances in Mechanical Engineering, 2014, 6(6): 1-11.

    [6]MADGWICK S O H. An efficient orientation filter for inertial and inertial/magnetic sensor arrays[R/OL]. [2015-04-30]. http://www.xio.co.uklres/doc/madgwickinternaIJeport.pdf.

    [7]成怡,金海林,修春波,等.四軸飛行器組合導(dǎo)航非線性濾波算法[J]. 計(jì)算機(jī)應(yīng)用,2014,34(S1):341-344.(CHENG Y, JIN H L, XIU C B, et al. Nonlinear filter algorithm for quadrotor integrated navigation [J]. Journal of Computer Applications, 2014, 34(S1):341-344.)

    [8]閻世梁,王銀玲,張華.基于改進(jìn)互補(bǔ)濾波器的低成本微小飛行器姿態(tài)估計(jì)方法[J]. 計(jì)算機(jī)應(yīng)用,2013,33(7):2078-2082.(YAN S L, WANG Y L, ZHANG H. Improved complementary filter for attitude estimation of micro air vehicles using lowcost inertial measurement units [J]. Journal of Computer Applications,2013, 33(7):2078-2082.)

    [9]BACHMANN E R, DUMAN I, USTA U Y, et al. Orientation tracking for humans and robots using inertial sensors[C]// Proceedings of the 1999 IEEE International Symposium on Computational Intelligence in Robotics and Automation. Piscataway, NJ: IEEE, 1999: 187-194.

    [10]MADGWICK S O H, HARRISON A J L, VAIDYANATHAN R. Estimation of IMU and MARG orientation using a gradient descent algorithm[C]// Proceedings of the 2011 IEEE International Conference on Rehabilitation Robotics. Piscataway, NJ: IEEE, 2011: 1-7.

    [11]MARINS J L, YUN X, BACHMANN E R, et al. An extended Kalman filter for quaternionbased orientation estimation using MARG sensors[C]// Proceedings of the 2001 IEEE/RSJ International Conference on Intelligent Robots and Systems. Piscataway, NJ: IEEE, 2001, 4: 2003-2011.

    [12]張?zhí)穑?吳美平, 張開東, 等. 基于精密單點(diǎn)定位的航空重力儀運(yùn)動(dòng)加速度估計(jì)[J]. 導(dǎo)航與控制, 2014, 13(4):26-31.(ZHANG T, WU M P, ZHANG K D, et al. Airborne gravitometer kinematic acceleration determination based on precise point positioning[J].Navigation and Control, 2014, 13(4): 26-31.)

    [13]李海濤. 加速度計(jì)/磁強(qiáng)計(jì)捷聯(lián)慣導(dǎo)系統(tǒng)姿態(tài)解算方法研究[D].太原:中北大學(xué),2007.(LI H T. Study on attitude algorithms of accelerometers/magnetometers strapdown inertial navigation system[D].Taiyuan: North University of China, 2007.)

    [14]HENDERSON D M. Euler angles, quaternions, and transformation matrices[R]. Washington, DC:NASA JSC, 1977.

    亚洲18禁久久av| 日本精品一区二区三区蜜桃| 日本在线视频免费播放| 国产淫片久久久久久久久| 在线观看一区二区三区| 白带黄色成豆腐渣| 国产欧美日韩精品一区二区| 老司机福利观看| 亚洲四区av| 激情 狠狠 欧美| 九色成人免费人妻av| 成年免费大片在线观看| 亚洲精品色激情综合| 精品不卡国产一区二区三区| 真实男女啪啪啪动态图| 岛国在线免费视频观看| 一区福利在线观看| 中文亚洲av片在线观看爽| 直男gayav资源| 日本一二三区视频观看| 黑人高潮一二区| 国产高潮美女av| 亚洲中文日韩欧美视频| 老司机午夜福利在线观看视频| 日日摸夜夜添夜夜爱| 色综合亚洲欧美另类图片| 永久网站在线| 亚洲av成人精品一区久久| 国产精品嫩草影院av在线观看| 国产一区二区亚洲精品在线观看| 91午夜精品亚洲一区二区三区| 亚洲成av人片在线播放无| 国产精品美女特级片免费视频播放器| 中文字幕熟女人妻在线| av国产免费在线观看| 网址你懂的国产日韩在线| 亚洲av电影不卡..在线观看| 18禁裸乳无遮挡免费网站照片| 国产av不卡久久| 深夜a级毛片| 黄色配什么色好看| 又爽又黄a免费视频| 春色校园在线视频观看| 精品久久国产蜜桃| 国产av在哪里看| 女同久久另类99精品国产91| 啦啦啦韩国在线观看视频| 欧美丝袜亚洲另类| 欧美日韩乱码在线| 少妇丰满av| 亚洲一区二区三区色噜噜| 尾随美女入室| 内地一区二区视频在线| 男插女下体视频免费在线播放| 一级毛片久久久久久久久女| 男女做爰动态图高潮gif福利片| 国产精品亚洲一级av第二区| 99久久精品一区二区三区| 中文字幕精品亚洲无线码一区| АⅤ资源中文在线天堂| 又黄又爽又刺激的免费视频.| 六月丁香七月| 国产精品一区二区三区四区久久| 久久人人精品亚洲av| 一级黄色大片毛片| 欧美最黄视频在线播放免费| 国产69精品久久久久777片| 成年女人毛片免费观看观看9| 麻豆乱淫一区二区| 国产69精品久久久久777片| 中国美女看黄片| 国产一区二区在线av高清观看| 国产三级在线视频| 大香蕉久久网| 国产真实乱freesex| 超碰av人人做人人爽久久| 日本一二三区视频观看| 校园春色视频在线观看| 国产精品99久久久久久久久| 日韩人妻高清精品专区| 一进一出好大好爽视频| 亚洲av一区综合| 一边摸一边抽搐一进一小说| 欧美中文日本在线观看视频| 人妻丰满熟妇av一区二区三区| 看免费成人av毛片| 亚洲av中文av极速乱| 国产精品av视频在线免费观看| 成年女人永久免费观看视频| 日韩高清综合在线| 99久国产av精品| 精品久久久久久久人妻蜜臀av| 夜夜看夜夜爽夜夜摸| 蜜臀久久99精品久久宅男| 日韩中字成人| 九九久久精品国产亚洲av麻豆| 18禁黄网站禁片免费观看直播| 久久久午夜欧美精品| av中文乱码字幕在线| 亚洲精品亚洲一区二区| av中文乱码字幕在线| 免费看a级黄色片| 大型黄色视频在线免费观看| 精品一区二区三区视频在线观看免费| av.在线天堂| 成年女人看的毛片在线观看| 欧美日韩综合久久久久久| 身体一侧抽搐| 免费在线观看影片大全网站| 国产精品综合久久久久久久免费| 国内少妇人妻偷人精品xxx网站| 精品人妻偷拍中文字幕| 国产乱人视频| 亚洲精品在线观看二区| 国产午夜精品久久久久久一区二区三区 | 亚洲专区国产一区二区| 欧美日韩综合久久久久久| 午夜福利在线观看免费完整高清在 | 亚洲成人久久爱视频| 三级毛片av免费| 啦啦啦啦在线视频资源| 蜜桃久久精品国产亚洲av| 丝袜喷水一区| 亚洲av五月六月丁香网| 日本成人三级电影网站| 欧美日韩一区二区视频在线观看视频在线 | 亚洲真实伦在线观看| 草草在线视频免费看| 国产一区二区在线av高清观看| 少妇猛男粗大的猛烈进出视频 | a级毛片免费高清观看在线播放| 日韩欧美精品免费久久| 美女 人体艺术 gogo| 人妻丰满熟妇av一区二区三区| 婷婷精品国产亚洲av| 亚洲乱码一区二区免费版| 亚洲av成人av| 久久韩国三级中文字幕| 97超级碰碰碰精品色视频在线观看| 直男gayav资源| 亚洲成av人片在线播放无| 欧美日韩综合久久久久久| 深夜a级毛片| 欧美精品国产亚洲| 嫩草影院新地址| 99在线视频只有这里精品首页| 亚洲美女搞黄在线观看 | 国产av不卡久久| 国产精品久久久久久av不卡| 欧美丝袜亚洲另类| 国产精品电影一区二区三区| h日本视频在线播放| 亚洲激情五月婷婷啪啪| 99热网站在线观看| а√天堂www在线а√下载| 欧美xxxx性猛交bbbb| 黄色日韩在线| 高清午夜精品一区二区三区 | 12—13女人毛片做爰片一| 日韩制服骚丝袜av| 亚洲国产欧美人成| 日日干狠狠操夜夜爽| 国产麻豆成人av免费视频| av黄色大香蕉| 两性午夜刺激爽爽歪歪视频在线观看| 老熟妇仑乱视频hdxx| 国产黄a三级三级三级人| 色播亚洲综合网| 亚洲国产精品成人久久小说 | 狂野欧美白嫩少妇大欣赏| 国产av一区在线观看免费| 中文字幕av在线有码专区| 国产午夜福利久久久久久| 一区福利在线观看| 亚洲av一区综合| 少妇人妻精品综合一区二区 | 一区二区三区四区激情视频 | 午夜老司机福利剧场| 国产片特级美女逼逼视频| 精品人妻熟女av久视频| 久久99热6这里只有精品| 免费黄网站久久成人精品| 国产aⅴ精品一区二区三区波| 99热精品在线国产| 国产伦在线观看视频一区| 亚洲av五月六月丁香网| 亚洲一区二区三区色噜噜| 神马国产精品三级电影在线观看| 亚洲成人精品中文字幕电影| 国产精品av视频在线免费观看| 精品久久久久久久久av| 99九九线精品视频在线观看视频| av视频在线观看入口| 一边摸一边抽搐一进一小说| 国产人妻一区二区三区在| 干丝袜人妻中文字幕| 中文字幕人妻熟人妻熟丝袜美| 天堂av国产一区二区熟女人妻| 日韩国内少妇激情av| 高清午夜精品一区二区三区 | 欧美三级亚洲精品| 欧美+亚洲+日韩+国产| 香蕉av资源在线| 长腿黑丝高跟| 波多野结衣高清作品| 综合色av麻豆| 国产精品久久久久久av不卡| 国产真实乱freesex| 久久99热6这里只有精品| 女的被弄到高潮叫床怎么办| 一本精品99久久精品77| 国产在线精品亚洲第一网站| 亚洲国产日韩欧美精品在线观看| 欧美绝顶高潮抽搐喷水| 国产黄片美女视频| 亚洲欧美精品自产自拍| 亚洲真实伦在线观看| 丰满的人妻完整版| 免费高清视频大片| 在线观看美女被高潮喷水网站| 免费人成视频x8x8入口观看| 男人和女人高潮做爰伦理| 97超级碰碰碰精品色视频在线观看| 内地一区二区视频在线| 联通29元200g的流量卡| 搞女人的毛片| 国产熟女欧美一区二区| 免费看a级黄色片| 亚洲av免费在线观看| or卡值多少钱| 久久天躁狠狠躁夜夜2o2o| 国产人妻一区二区三区在| 成年av动漫网址| 成人一区二区视频在线观看| 亚洲国产精品合色在线| 人妻夜夜爽99麻豆av| 一区福利在线观看| 亚洲七黄色美女视频| 欧美3d第一页| 小蜜桃在线观看免费完整版高清| 小说图片视频综合网站| 91久久精品电影网| 在线观看一区二区三区| 久久精品综合一区二区三区| 卡戴珊不雅视频在线播放| 亚洲人与动物交配视频| 成熟少妇高潮喷水视频| 看黄色毛片网站| 午夜福利在线观看免费完整高清在 | 高清毛片免费看| 18+在线观看网站| 久久鲁丝午夜福利片| 成人av一区二区三区在线看| 少妇丰满av| 噜噜噜噜噜久久久久久91| 一本久久中文字幕| 99久国产av精品| 国产精品人妻久久久久久| 夜夜看夜夜爽夜夜摸| 日产精品乱码卡一卡2卡三| 久久久久久久久久黄片| 少妇猛男粗大的猛烈进出视频 | 久久亚洲精品不卡| 精品久久久久久久久av| 亚洲婷婷狠狠爱综合网| 国产真实乱freesex| 国产乱人偷精品视频| 日韩欧美在线乱码| 亚洲一区二区三区色噜噜| 欧美精品国产亚洲| 99riav亚洲国产免费| 午夜久久久久精精品| 精品少妇黑人巨大在线播放 | 国产精华一区二区三区| 亚洲精品456在线播放app| 精品久久久久久久人妻蜜臀av| 亚洲人成网站在线播放欧美日韩| 亚洲自拍偷在线| 国产一区二区激情短视频| 国产精品久久视频播放| 久久草成人影院| 美女被艹到高潮喷水动态| 高清午夜精品一区二区三区 | 国产午夜精品久久久久久一区二区三区 | 此物有八面人人有两片| 日韩欧美一区二区三区在线观看| 久久久精品欧美日韩精品| 好男人在线观看高清免费视频| 老师上课跳d突然被开到最大视频| 亚洲,欧美,日韩| 麻豆久久精品国产亚洲av| 99久久精品热视频| 日日干狠狠操夜夜爽| 美女高潮的动态| 亚洲人成网站高清观看| 国产精品国产高清国产av| 亚洲美女黄片视频| 色综合站精品国产| 简卡轻食公司| 在现免费观看毛片| 人妻久久中文字幕网| 国产高清三级在线| 赤兔流量卡办理| 国产极品精品免费视频能看的| 一进一出好大好爽视频| 国产一区二区亚洲精品在线观看| 成人特级av手机在线观看| 日韩欧美免费精品| 男女视频在线观看网站免费| 如何舔出高潮| 91狼人影院| 日本免费一区二区三区高清不卡| 一进一出好大好爽视频| 免费搜索国产男女视频| 欧美丝袜亚洲另类| 亚洲av成人精品一区久久| 国产精品久久久久久久电影| 日本 av在线| 久久精品综合一区二区三区| 精品欧美国产一区二区三| 国产aⅴ精品一区二区三区波| 两性午夜刺激爽爽歪歪视频在线观看| 国产av在哪里看| 偷拍熟女少妇极品色| 韩国av在线不卡| 在线国产一区二区在线| 亚洲精品亚洲一区二区| 床上黄色一级片| 久久久色成人| 又黄又爽又刺激的免费视频.| 欧美日韩一区二区视频在线观看视频在线 | 九九在线视频观看精品| 精品午夜福利在线看| 日韩欧美精品免费久久| 日韩欧美国产在线观看| 全区人妻精品视频| 精品午夜福利在线看| 国内精品宾馆在线| 亚洲在线观看片| 国产免费一级a男人的天堂| 亚洲av免费在线观看| 亚洲欧美中文字幕日韩二区| 黄色欧美视频在线观看| 亚洲欧美精品自产自拍| 欧美精品国产亚洲| 亚洲国产色片| 人妻丰满熟妇av一区二区三区| 精品一区二区三区人妻视频| 女的被弄到高潮叫床怎么办| 在线国产一区二区在线| 亚洲人成网站在线播放欧美日韩| 99久久中文字幕三级久久日本| 欧美xxxx黑人xx丫x性爽| 国产精品永久免费网站| 亚洲专区国产一区二区| 丰满的人妻完整版| 男人的好看免费观看在线视频| 午夜精品一区二区三区免费看| 国内精品美女久久久久久| 国产高清视频在线观看网站| 国产亚洲av嫩草精品影院| 观看免费一级毛片| 寂寞人妻少妇视频99o| 天堂√8在线中文| av女优亚洲男人天堂| 波多野结衣高清无吗| 免费电影在线观看免费观看| 日韩欧美精品免费久久| 男女那种视频在线观看| 国产av在哪里看| 综合色丁香网| 国产v大片淫在线免费观看| 亚洲av五月六月丁香网| 色综合亚洲欧美另类图片| 男人的好看免费观看在线视频| 免费av观看视频| 日韩欧美精品v在线| 久久午夜福利片| 联通29元200g的流量卡| 18禁黄网站禁片免费观看直播| 亚洲五月天丁香| 观看免费一级毛片| 99久久久亚洲精品蜜臀av| 黄色视频,在线免费观看| 美女高潮的动态| 白带黄色成豆腐渣| 黑人高潮一二区| 亚洲人与动物交配视频| 日本熟妇午夜| 色综合亚洲欧美另类图片| 精品一区二区三区视频在线| 成人精品一区二区免费| 美女免费视频网站| videossex国产| 99热这里只有是精品50| 日韩人妻高清精品专区| 18禁在线播放成人免费| 一卡2卡三卡四卡精品乱码亚洲| 日本成人三级电影网站| 美女黄网站色视频| 免费在线观看成人毛片| 老熟妇乱子伦视频在线观看| 亚洲自偷自拍三级| 国产视频内射| 黄色日韩在线| 国产欧美日韩精品亚洲av| 成人鲁丝片一二三区免费| 亚洲色图av天堂| 亚洲18禁久久av| 99热这里只有是精品在线观看| 不卡视频在线观看欧美| 日日撸夜夜添| 51国产日韩欧美| 亚洲最大成人中文| 久久精品国产清高在天天线| av女优亚洲男人天堂| 熟妇人妻久久中文字幕3abv| 亚洲欧美成人综合另类久久久 | 国产黄片美女视频| 深夜精品福利| 在线免费观看不下载黄p国产| 大香蕉久久网| 免费不卡的大黄色大毛片视频在线观看 | 91麻豆精品激情在线观看国产| 亚洲成人久久爱视频| 少妇熟女aⅴ在线视频| 精品乱码久久久久久99久播| 能在线免费观看的黄片| 伦理电影大哥的女人| 免费不卡的大黄色大毛片视频在线观看 | 国产一区二区三区在线臀色熟女| 午夜福利成人在线免费观看| 国产精品久久久久久久电影| 露出奶头的视频| 午夜精品国产一区二区电影 | 午夜老司机福利剧场| 青春草视频在线免费观看| 久久亚洲精品不卡| 午夜精品国产一区二区电影 | 久久久久精品国产欧美久久久| 深夜a级毛片| 免费在线观看成人毛片| 日韩一本色道免费dvd| 久久国内精品自在自线图片| 插逼视频在线观看| 成人午夜高清在线视频| 久久久欧美国产精品| 哪里可以看免费的av片| 丝袜喷水一区| 国内揄拍国产精品人妻在线| 国语自产精品视频在线第100页| 一a级毛片在线观看| 久久久久性生活片| 亚洲欧美日韩卡通动漫| 99久久精品一区二区三区| 99热网站在线观看| 老司机福利观看| 精品人妻偷拍中文字幕| 午夜免费激情av| 看十八女毛片水多多多| 欧美性感艳星| 日韩大尺度精品在线看网址| 我的女老师完整版在线观看| 午夜福利成人在线免费观看| 亚洲精品在线观看二区| 亚洲国产精品成人久久小说 | 国产白丝娇喘喷水9色精品| 亚洲性夜色夜夜综合| 亚洲欧美日韩高清在线视频| 国产精品久久久久久精品电影| 精品午夜福利视频在线观看一区| av中文乱码字幕在线| 亚洲欧美中文字幕日韩二区| 国产一级毛片七仙女欲春2| 六月丁香七月| 久久久久国产精品人妻aⅴ院| 国产精品精品国产色婷婷| 一级毛片aaaaaa免费看小| 免费不卡的大黄色大毛片视频在线观看 | 男插女下体视频免费在线播放| 久久鲁丝午夜福利片| 亚洲精品粉嫩美女一区| 日本一本二区三区精品| 免费av毛片视频| 国产精品久久电影中文字幕| 搡老岳熟女国产| aaaaa片日本免费| 亚洲av中文av极速乱| 国产淫片久久久久久久久| 国内精品美女久久久久久| 国产亚洲av嫩草精品影院| 欧美日韩综合久久久久久| 97热精品久久久久久| 一级黄色大片毛片| 99国产精品一区二区蜜桃av| 免费av观看视频| a级毛片a级免费在线| a级一级毛片免费在线观看| 亚洲av美国av| 久久久久精品国产欧美久久久| 国产伦精品一区二区三区视频9| 久久久国产成人精品二区| 黄色日韩在线| 偷拍熟女少妇极品色| 久久久久精品国产欧美久久久| 欧美bdsm另类| 亚洲av不卡在线观看| 日韩精品青青久久久久久| 国产片特级美女逼逼视频| 国产免费一级a男人的天堂| 久久久色成人| 1000部很黄的大片| 国产欧美日韩一区二区精品| av在线蜜桃| 18禁裸乳无遮挡免费网站照片| 国产高清不卡午夜福利| 欧美成人一区二区免费高清观看| 精品福利观看| 在线观看午夜福利视频| 欧美区成人在线视频| 性色avwww在线观看| 国产视频一区二区在线看| 男插女下体视频免费在线播放| 亚洲国产高清在线一区二区三| 午夜福利视频1000在线观看| 三级经典国产精品| 精品久久久久久久末码| 在线a可以看的网站| 日本一二三区视频观看| 亚洲一区二区三区色噜噜| 欧美色视频一区免费| 亚洲一区二区三区色噜噜| 在现免费观看毛片| 国产女主播在线喷水免费视频网站 | 午夜福利高清视频| 热99re8久久精品国产| 在线观看av片永久免费下载| 欧美色视频一区免费| 1000部很黄的大片| 别揉我奶头~嗯~啊~动态视频| 精品人妻偷拍中文字幕| 中文字幕精品亚洲无线码一区| 日韩精品青青久久久久久| 秋霞在线观看毛片| 99国产极品粉嫩在线观看| 三级毛片av免费| 中文字幕av成人在线电影| 国产美女午夜福利| 身体一侧抽搐| 禁无遮挡网站| av天堂中文字幕网| 亚洲精品影视一区二区三区av| 成年版毛片免费区| 日日摸夜夜添夜夜爱| 精品福利观看| 91久久精品国产一区二区成人| 色综合色国产| 丝袜美腿在线中文| 亚洲av二区三区四区| 一区二区三区四区激情视频 | 亚洲中文日韩欧美视频| 欧美一区二区国产精品久久精品| 欧美最黄视频在线播放免费| 少妇的逼好多水| 国产探花在线观看一区二区| 国产精品一区二区三区四区久久| 一本精品99久久精品77| 我的女老师完整版在线观看| 又爽又黄a免费视频| www.色视频.com| 尤物成人国产欧美一区二区三区| 久久精品国产清高在天天线| 日韩欧美国产在线观看| 高清日韩中文字幕在线| 99久久成人亚洲精品观看| 国产成人91sexporn| 观看免费一级毛片| 国产精品人妻久久久影院| 午夜精品一区二区三区免费看| 国产亚洲91精品色在线| 国产成人a区在线观看| 国产一区二区亚洲精品在线观看| 美女被艹到高潮喷水动态| 亚洲精品国产成人久久av| 国产亚洲精品综合一区在线观看| 亚洲性夜色夜夜综合| 亚洲性久久影院| av在线播放精品| 欧美最新免费一区二区三区| 性插视频无遮挡在线免费观看| 亚洲人成网站在线播放欧美日韩| 免费观看的影片在线观看| 日本-黄色视频高清免费观看| 男人和女人高潮做爰伦理| 免费在线观看影片大全网站| 中文字幕精品亚洲无线码一区| 深夜a级毛片| 久久久久精品国产欧美久久久| 国产成年人精品一区二区| 青春草视频在线免费观看| 欧美性猛交黑人性爽| 色5月婷婷丁香| 日日干狠狠操夜夜爽| 三级男女做爰猛烈吃奶摸视频| 精品久久久久久久末码| 寂寞人妻少妇视频99o| 可以在线观看的亚洲视频| av中文乱码字幕在线| 色av中文字幕| 国产伦精品一区二区三区视频9| 老女人水多毛片| 深夜精品福利| 麻豆一二三区av精品| 美女xxoo啪啪120秒动态图| 黄色配什么色好看|