• <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.

    侵犯人妻中文字幕一二三四区| 亚洲 欧美一区二区三区| 日日摸夜夜添夜夜爱| 午夜福利视频精品| 狂野欧美激情性bbbbbb| 伊人久久大香线蕉亚洲五| 在线观看人妻少妇| 亚洲欧美成人精品一区二区| 国产在线视频一区二区| 中文天堂在线官网| 日本wwww免费看| 成人三级做爰电影| 中文字幕av电影在线播放| 久久久精品国产亚洲av高清涩受| 免费不卡黄色视频| 国产精品免费大片| 两个人看的免费小视频| 婷婷色综合www| 99热全是精品| 欧美av亚洲av综合av国产av | 美女大奶头黄色视频| av.在线天堂| www.av在线官网国产| 两个人免费观看高清视频| 亚洲av综合色区一区| 婷婷色综合www| 久久久久久久久免费视频了| 国产精品亚洲av一区麻豆 | 亚洲av电影在线进入| 麻豆精品久久久久久蜜桃| 最近中文字幕高清免费大全6| 日韩av在线免费看完整版不卡| 日本91视频免费播放| 国产精品熟女久久久久浪| 99久久精品国产亚洲精品| 永久免费av网站大全| 亚洲人成网站在线观看播放| 亚洲欧美一区二区三区久久| 亚洲国产av新网站| 在线亚洲精品国产二区图片欧美| 国产亚洲最大av| 欧美日韩亚洲综合一区二区三区_| 亚洲国产欧美网| 国产精品成人在线| 赤兔流量卡办理| 在线观看www视频免费| 91老司机精品| 黑丝袜美女国产一区| 亚洲 欧美一区二区三区| 最近最新中文字幕免费大全7| 亚洲av电影在线进入| 国产精品麻豆人妻色哟哟久久| 精品久久久精品久久久| e午夜精品久久久久久久| 亚洲七黄色美女视频| videos熟女内射| 咕卡用的链子| 国产亚洲最大av| 亚洲精品久久久久久婷婷小说| 亚洲精品中文字幕在线视频| 美女福利国产在线| 捣出白浆h1v1| 熟妇人妻不卡中文字幕| 亚洲欧美一区二区三区久久| 午夜激情久久久久久久| 国产精品偷伦视频观看了| avwww免费| 日韩欧美精品免费久久| 精品少妇黑人巨大在线播放| 婷婷色综合www| 亚洲精品,欧美精品| 日韩 欧美 亚洲 中文字幕| 午夜激情久久久久久久| 夫妻性生交免费视频一级片| 在线观看免费视频网站a站| 日本欧美国产在线视频| 久久热在线av| 天天躁夜夜躁狠狠久久av| 亚洲av电影在线进入| 成年av动漫网址| 国产老妇伦熟女老妇高清| 亚洲精品国产色婷婷电影| 国产成人91sexporn| 黑丝袜美女国产一区| 国产日韩欧美在线精品| 9热在线视频观看99| 欧美精品av麻豆av| 大话2 男鬼变身卡| 国产不卡av网站在线观看| 亚洲欧美一区二区三区国产| 搡老乐熟女国产| 国产国语露脸激情在线看| 国产一区二区激情短视频 | 欧美人与性动交α欧美精品济南到| 在线天堂中文资源库| 国产亚洲一区二区精品| 一区二区三区四区激情视频| 日韩一卡2卡3卡4卡2021年| 十八禁人妻一区二区| 国产野战对白在线观看| 亚洲精品国产av蜜桃| 国产亚洲一区二区精品| 国产成人欧美| 亚洲五月色婷婷综合| 一级片'在线观看视频| 国产精品成人在线| 狠狠婷婷综合久久久久久88av| 男女无遮挡免费网站观看| 午夜激情av网站| 黄色怎么调成土黄色| bbb黄色大片| 国产精品一区二区在线观看99| 国产日韩欧美亚洲二区| 日韩 亚洲 欧美在线| 2018国产大陆天天弄谢| av免费观看日本| 午夜日本视频在线| 久久久久久久久久久免费av| 人妻一区二区av| 一区二区日韩欧美中文字幕| 丝袜喷水一区| 国产xxxxx性猛交| 久久精品熟女亚洲av麻豆精品| 涩涩av久久男人的天堂| 999久久久国产精品视频| 18禁动态无遮挡网站| 久久精品国产综合久久久| 伦理电影免费视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产熟女欧美一区二区| 尾随美女入室| 日韩一区二区视频免费看| 日韩电影二区| 久久婷婷青草| 亚洲精品久久成人aⅴ小说| 日本91视频免费播放| 天天躁狠狠躁夜夜躁狠狠躁| videos熟女内射| 天天躁狠狠躁夜夜躁狠狠躁| 精品亚洲成a人片在线观看| 日日啪夜夜爽| 美女午夜性视频免费| 另类亚洲欧美激情| 纵有疾风起免费观看全集完整版| 亚洲在久久综合| 亚洲精品自拍成人| 大片免费播放器 马上看| 国产在线一区二区三区精| 肉色欧美久久久久久久蜜桃| 成人三级做爰电影| 超碰97精品在线观看| 国产精品.久久久| 最近最新中文字幕免费大全7| 曰老女人黄片| 亚洲,欧美精品.| 亚洲精品自拍成人| 精品少妇久久久久久888优播| 国产无遮挡羞羞视频在线观看| 免费黄网站久久成人精品| 黄色视频不卡| 日本wwww免费看| 精品卡一卡二卡四卡免费| 人人妻,人人澡人人爽秒播 | 成年人午夜在线观看视频| 91老司机精品| 黄色怎么调成土黄色| 国产精品一区二区精品视频观看| 精品视频人人做人人爽| 亚洲欧美一区二区三区黑人| 欧美日韩亚洲国产一区二区在线观看 | av免费观看日本| 亚洲国产欧美在线一区| 啦啦啦视频在线资源免费观看| 国产精品免费视频内射| 蜜桃在线观看..| 大香蕉久久成人网| 黑人欧美特级aaaaaa片| 麻豆乱淫一区二区| 久久久欧美国产精品| 汤姆久久久久久久影院中文字幕| 国产一区二区 视频在线| 一边摸一边做爽爽视频免费| 亚洲,欧美精品.| 制服诱惑二区| 日日爽夜夜爽网站| 在线观看人妻少妇| 日韩一本色道免费dvd| 国产人伦9x9x在线观看| 性色av一级| 熟女av电影| 中文欧美无线码| 激情五月婷婷亚洲| 久久久久视频综合| 最近2019中文字幕mv第一页| 日韩欧美精品免费久久| 超碰成人久久| 大片免费播放器 马上看| 国产免费视频播放在线视频| 亚洲情色 制服丝袜| 久久青草综合色| 一二三四中文在线观看免费高清| 亚洲综合精品二区| 中国国产av一级| 日本av手机在线免费观看| 中文字幕高清在线视频| 一本久久精品| 亚洲精品国产av蜜桃| 男女午夜视频在线观看| 欧美xxⅹ黑人| 欧美日韩视频高清一区二区三区二| 国产成人欧美| 欧美最新免费一区二区三区| 一级毛片 在线播放| 欧美人与性动交α欧美软件| 晚上一个人看的免费电影| 乱人伦中国视频| 国产成人一区二区在线| 又大又黄又爽视频免费| 国产精品免费视频内射| 深夜精品福利| 看非洲黑人一级黄片| 精品少妇内射三级| 久久久精品免费免费高清| 欧美日韩成人在线一区二区| 免费观看性生交大片5| 18禁观看日本| 在线天堂中文资源库| 菩萨蛮人人尽说江南好唐韦庄| 亚洲自偷自拍图片 自拍| 美女午夜性视频免费| 亚洲中文av在线| 狂野欧美激情性xxxx| 一级毛片我不卡| 亚洲欧美色中文字幕在线| 亚洲av日韩在线播放| 国产高清国产精品国产三级| 国产成人欧美| 国产福利在线免费观看视频| 久久久久久久久免费视频了| 国产欧美日韩一区二区三区在线| 男女之事视频高清在线观看 | 美国免费a级毛片| 国产午夜精品一二区理论片| 亚洲精品美女久久久久99蜜臀 | 超碰成人久久| 精品第一国产精品| 亚洲综合色网址| 男女高潮啪啪啪动态图| 国语对白做爰xxxⅹ性视频网站| 日韩中文字幕欧美一区二区 | 国精品久久久久久国模美| 日本欧美国产在线视频| 尾随美女入室| 丰满乱子伦码专区| 制服丝袜香蕉在线| 老汉色∧v一级毛片| 亚洲熟女精品中文字幕| 欧美日韩一级在线毛片| 精品一区在线观看国产| 曰老女人黄片| 少妇人妻精品综合一区二区| 中国三级夫妇交换| 性高湖久久久久久久久免费观看| 国产一区二区三区综合在线观看| 麻豆精品久久久久久蜜桃| 高清av免费在线| 制服丝袜香蕉在线| 亚洲av日韩精品久久久久久密 | 中文字幕人妻丝袜一区二区 | 亚洲一级一片aⅴ在线观看| 亚洲欧美一区二区三区久久| 午夜福利视频在线观看免费| 美女国产高潮福利片在线看| 欧美久久黑人一区二区| 最新的欧美精品一区二区| 国产精品av久久久久免费| 国产精品国产av在线观看| 妹子高潮喷水视频| 亚洲成人手机| 欧美中文综合在线视频| 亚洲欧美成人综合另类久久久| 亚洲av电影在线进入| 一区二区三区激情视频| 国产伦人伦偷精品视频| 婷婷色综合大香蕉| 日本一区二区免费在线视频| 国产精品免费大片| 99九九在线精品视频| 欧美亚洲日本最大视频资源| 久久亚洲国产成人精品v| 哪个播放器可以免费观看大片| 久久这里只有精品19| 国产视频首页在线观看| 欧美日韩精品网址| 亚洲国产毛片av蜜桃av| 国产伦人伦偷精品视频| 国产精品人妻久久久影院| 亚洲精品国产av成人精品| 免费少妇av软件| 婷婷色麻豆天堂久久| 在线 av 中文字幕| 天堂中文最新版在线下载| 啦啦啦视频在线资源免费观看| a级片在线免费高清观看视频| xxx大片免费视频| 亚洲图色成人| 人人妻人人澡人人看| 久久久久久免费高清国产稀缺| 一级,二级,三级黄色视频| 我要看黄色一级片免费的| 19禁男女啪啪无遮挡网站| 久久午夜综合久久蜜桃| 国产午夜精品一二区理论片| av又黄又爽大尺度在线免费看| 亚洲av成人不卡在线观看播放网 | 久久精品国产亚洲av涩爱| 多毛熟女@视频| 无遮挡黄片免费观看| 国产极品粉嫩免费观看在线| 丝袜在线中文字幕| 亚洲色图 男人天堂 中文字幕| 免费黄色在线免费观看| 永久免费av网站大全| 欧美乱码精品一区二区三区| 久久久久久久久久久免费av| 一区二区三区四区激情视频| 欧美老熟妇乱子伦牲交| 精品一区在线观看国产| 人妻 亚洲 视频| 亚洲一级一片aⅴ在线观看| 国产在线一区二区三区精| 欧美人与性动交α欧美精品济南到| 精品少妇一区二区三区视频日本电影 | 久久精品aⅴ一区二区三区四区| 成人国语在线视频| 搡老岳熟女国产| 免费不卡黄色视频| 十八禁高潮呻吟视频| 精品国产一区二区久久| 丝瓜视频免费看黄片| 欧美日韩福利视频一区二区| 18禁动态无遮挡网站| 街头女战士在线观看网站| 日韩 亚洲 欧美在线| 成人毛片60女人毛片免费| 最近手机中文字幕大全| 天堂中文最新版在线下载| 国产在视频线精品| 免费观看人在逋| 王馨瑶露胸无遮挡在线观看| 亚洲av日韩精品久久久久久密 | 2021少妇久久久久久久久久久| 免费高清在线观看日韩| 亚洲欧美一区二区三区久久| 美女脱内裤让男人舔精品视频| 精品国产国语对白av| 女人久久www免费人成看片| 日本av手机在线免费观看| 久久亚洲国产成人精品v| 中国国产av一级| av视频免费观看在线观看| 大片电影免费在线观看免费| 日韩视频在线欧美| 一边亲一边摸免费视频| 秋霞伦理黄片| 久久久久精品国产欧美久久久 | 久久久久国产精品人妻一区二区| 日韩av不卡免费在线播放| 亚洲国产精品成人久久小说| 日韩人妻精品一区2区三区| 亚洲一区中文字幕在线| 欧美精品一区二区大全| 日韩 亚洲 欧美在线| 一级,二级,三级黄色视频| 少妇 在线观看| 成人亚洲欧美一区二区av| 青春草亚洲视频在线观看| 久久久久久久大尺度免费视频| 国产熟女欧美一区二区| 久久热在线av| videos熟女内射| 又大又爽又粗| 亚洲精品国产av成人精品| 亚洲国产日韩一区二区| 国产人伦9x9x在线观看| 97在线人人人人妻| 一级,二级,三级黄色视频| 国产精品久久久久成人av| 热99久久久久精品小说推荐| 国产97色在线日韩免费| 国产xxxxx性猛交| 精品午夜福利在线看| 亚洲国产看品久久| 国产亚洲最大av| 午夜免费男女啪啪视频观看| 搡老乐熟女国产| 尾随美女入室| 中文字幕av电影在线播放| av有码第一页| 日韩免费高清中文字幕av| 欧美日韩一区二区视频在线观看视频在线| 男人爽女人下面视频在线观看| 午夜福利网站1000一区二区三区| 日韩 亚洲 欧美在线| 日韩大码丰满熟妇| 国产精品久久久久久精品古装| 一区在线观看完整版| 晚上一个人看的免费电影| 亚洲天堂av无毛| 晚上一个人看的免费电影| 婷婷色综合大香蕉| 亚洲av成人不卡在线观看播放网 | 欧美精品人与动牲交sv欧美| 最近的中文字幕免费完整| 久久人妻熟女aⅴ| 国产精品熟女久久久久浪| 男的添女的下面高潮视频| a 毛片基地| 视频在线观看一区二区三区| 久久久精品区二区三区| 别揉我奶头~嗯~啊~动态视频 | 亚洲欧美激情在线| 最黄视频免费看| 韩国精品一区二区三区| 国产精品人妻久久久影院| 久久婷婷青草| 欧美久久黑人一区二区| videosex国产| 婷婷色综合www| 交换朋友夫妻互换小说| 日韩一区二区三区影片| 免费在线观看黄色视频的| 九九爱精品视频在线观看| 亚洲国产av影院在线观看| 极品少妇高潮喷水抽搐| 久久久久国产一级毛片高清牌| 精品国产一区二区三区四区第35| www.熟女人妻精品国产| 亚洲成av片中文字幕在线观看| 久久99热这里只频精品6学生| 九草在线视频观看| 欧美乱码精品一区二区三区| 美女大奶头黄色视频| 天堂俺去俺来也www色官网| 亚洲熟女精品中文字幕| 精品亚洲成a人片在线观看| 在线观看三级黄色| 久久久欧美国产精品| 久久人人爽av亚洲精品天堂| 精品酒店卫生间| 好男人视频免费观看在线| 国产在线视频一区二区| 美女福利国产在线| 高清黄色对白视频在线免费看| 国产精品一区二区精品视频观看| 亚洲熟女精品中文字幕| 国产野战对白在线观看| 夫妻午夜视频| 久久鲁丝午夜福利片| 欧美日韩一区二区视频在线观看视频在线| 久久久久久久国产电影| 久久久久精品久久久久真实原创| 国产精品麻豆人妻色哟哟久久| 亚洲精品国产一区二区精华液| av在线老鸭窝| 99热国产这里只有精品6| 女性被躁到高潮视频| videosex国产| 人人妻人人添人人爽欧美一区卜| 亚洲第一av免费看| 欧美日韩福利视频一区二区| 久热爱精品视频在线9| 美女大奶头黄色视频| a级毛片黄视频| 亚洲av中文av极速乱| 女人被躁到高潮嗷嗷叫费观| 亚洲美女视频黄频| 在线观看免费视频网站a站| 丰满迷人的少妇在线观看| 国产成人精品福利久久| 最近手机中文字幕大全| 日韩av在线免费看完整版不卡| 日韩 亚洲 欧美在线| 亚洲熟女精品中文字幕| 欧美av亚洲av综合av国产av | 99久久人妻综合| 高清av免费在线| 亚洲,欧美精品.| 99久久综合免费| 亚洲第一区二区三区不卡| 欧美日韩精品网址| 欧美日韩视频精品一区| 国产亚洲av片在线观看秒播厂| 中文精品一卡2卡3卡4更新| 少妇被粗大的猛进出69影院| 亚洲综合色网址| 在线观看人妻少妇| 亚洲国产av新网站| 熟妇人妻不卡中文字幕| 人人妻人人澡人人看| 久久久久久人妻| 亚洲天堂av无毛| 极品少妇高潮喷水抽搐| 一级毛片电影观看| 国产欧美亚洲国产| 精品久久久精品久久久| 一级,二级,三级黄色视频| 亚洲专区中文字幕在线 | 18禁国产床啪视频网站| 好男人视频免费观看在线| 一级片免费观看大全| 国产一级毛片在线| av又黄又爽大尺度在线免费看| 午夜福利,免费看| 精品卡一卡二卡四卡免费| 久久免费观看电影| 悠悠久久av| 国产男人的电影天堂91| 国产精品二区激情视频| 国产在视频线精品| 在线观看免费日韩欧美大片| videos熟女内射| 在线观看国产h片| 亚洲成色77777| 国产精品免费视频内射| 亚洲欧洲国产日韩| 久久婷婷青草| 中文字幕av电影在线播放| 永久免费av网站大全| 18禁动态无遮挡网站| 大片电影免费在线观看免费| 19禁男女啪啪无遮挡网站| a级毛片黄视频| av福利片在线| 午夜免费观看性视频| 啦啦啦视频在线资源免费观看| 欧美精品一区二区大全| 免费黄网站久久成人精品| 熟女少妇亚洲综合色aaa.| 久久久久精品人妻al黑| 99久国产av精品国产电影| 成人手机av| 在线观看一区二区三区激情| 亚洲情色 制服丝袜| av.在线天堂| 国产日韩欧美视频二区| 日本wwww免费看| 久久天堂一区二区三区四区| 一区二区三区乱码不卡18| 日韩一本色道免费dvd| 日韩制服丝袜自拍偷拍| 国产亚洲午夜精品一区二区久久| 美女高潮到喷水免费观看| 老司机在亚洲福利影院| 久久久久精品久久久久真实原创| 久热这里只有精品99| 久久精品久久精品一区二区三区| 亚洲欧美精品自产自拍| 婷婷色av中文字幕| 精品久久蜜臀av无| 亚洲精品美女久久av网站| 我要看黄色一级片免费的| 波多野结衣一区麻豆| 国产淫语在线视频| 91精品伊人久久大香线蕉| 国产免费福利视频在线观看| 另类精品久久| 成人三级做爰电影| 黄色视频不卡| 黄色视频在线播放观看不卡| 免费在线观看黄色视频的| 亚洲人成网站在线观看播放| 97精品久久久久久久久久精品| 日本wwww免费看| 欧美激情 高清一区二区三区| 国产精品久久久久成人av| 69精品国产乱码久久久| 一区二区日韩欧美中文字幕| 王馨瑶露胸无遮挡在线观看| 18在线观看网站| 一个人免费看片子| 黄色怎么调成土黄色| 在线看a的网站| 欧美日韩亚洲高清精品| 欧美少妇被猛烈插入视频| 叶爱在线成人免费视频播放| 999精品在线视频| 啦啦啦在线观看免费高清www| 伦理电影大哥的女人| 久久人人爽人人片av| 女人精品久久久久毛片| 亚洲第一青青草原| 无限看片的www在线观看| 免费观看a级毛片全部| 成人午夜精彩视频在线观看| av电影中文网址| 亚洲美女黄色视频免费看| 免费少妇av软件| 热re99久久国产66热| 青草久久国产| 欧美在线黄色| 大片免费播放器 马上看| 两性夫妻黄色片| 亚洲国产看品久久| 99久久99久久久精品蜜桃| 久久99精品国语久久久| 国产免费又黄又爽又色| 99久久综合免费| 如日韩欧美国产精品一区二区三区| 欧美97在线视频| netflix在线观看网站| 街头女战士在线观看网站| 日韩一区二区视频免费看| 亚洲欧美清纯卡通| 欧美日韩一级在线毛片|