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

    改進型自適應無跡卡爾曼姿態(tài)算法*

    2017-11-03 12:32:42侯江寬賈華宇
    傳感技術學報 2017年10期
    關鍵詞:磁強計無跡陀螺儀

    侯江寬,馬 珺*,賈華宇

    (1.太原理工大學物理與光電工程學院,山西 晉中 030600;2.太原理工大學信息工程學院,山西 晉中 030600)

    項目來源:山西省自然科學基金項目(2015011050)

    2017-02-28修改日期2017-06-18

    改進型自適應無跡卡爾曼姿態(tài)算法*

    侯江寬1,馬 珺1*,賈華宇2

    (1.太原理工大學物理與光電工程學院,山西 晉中 030600;2.太原理工大學信息工程學院,山西 晉中 030600)

    提出了一種改進型自適應無跡卡爾曼濾波姿態(tài)算法,能夠有效的解決MEMS陀螺儀的漂移和噪聲問題,同時減小運動加速度對加速度計的影響。將改進的自適應Sage-Husa算法與無跡卡爾曼濾波器相結合,使量測噪聲統計特性在線更新,提高系統的抗干擾能力,避免擴展卡爾曼濾波的線性化誤差,可以得到精確的全姿態(tài)角。每次迭代只更新3個歐拉角,提高了系統的解算速度。飛行實驗和分析表明:改進算法能夠有效的提高姿態(tài)解算精度,收斂速度快,自適應能力強,穩(wěn)定可靠,具有較強的魯棒性,在干擾消失時能夠快速得到準確姿態(tài)角。

    傳感器;姿態(tài)解算;無跡卡爾曼濾波器;自適應Sage-Husa算法;多旋翼飛行器

    獲取準確的姿態(tài)是無人飛行器實現自主飛行的前提條件[1],由于MEMS(Micro Electro Mechanical Systems)傳感器的高集成度、低成本,在微型飛行器上大量使用。多旋翼飛行器上的慣性測量單元主要由MEMS三軸陀螺儀、三軸加速度計、三軸磁強計組成。但是MEMS陀螺儀存在漂移問題,會導致累積誤差,因此需要對陀螺儀輸出的數據進行有效的處理才能獲得穩(wěn)定可靠的姿態(tài)信息。

    目前常用的姿態(tài)解算方法有Robert Mahony的互補濾波姿態(tài)解算算法[2]及其改進算法[3]、Sebastian Madgwick提出的基于梯度下降法的姿態(tài)解算算法[7]、基于卡爾曼濾波器的姿態(tài)解算算法[8]以及粒子濾波姿態(tài)解算算法[9]?;パa濾波算法是比較實用的一種算法,適用于對姿態(tài)精確度要求不高的小型飛行器?;パa濾波姿態(tài)解算的改進主要集中于濾波器系數的自適應調節(jié)[5],以提高算法對于噪聲的補償,這種自適應調節(jié)主要是結合了模糊算法,精度并不高。Madgwick將互補濾波的狀態(tài)更新方程由DCM(Direction Cosine Matrix)矩陣變?yōu)樗脑獢滴⒎址匠?大大提高了互補濾波算法的解算速度。Kalman濾波是目前較為適合多傳感器姿態(tài)融合的算法,但常規(guī)Kalman濾波算法需要系統噪聲和量測噪聲的統計特性已知[14]。然而在多旋翼飛行器的多傳感器融合中,系統噪聲和量測噪聲的統計特性往往是未知的且時變的,這會使常規(guī)Kalman濾波估計精度降低,失去算法最優(yōu)性,甚至導致濾波發(fā)散。擴展卡爾曼濾波姿態(tài)解算算法是對非線性函數進行線性化,但是在線性化過程中必須求解非線性函數的Jacobi矩陣,這給嵌入式微處理器帶來了大量的計算負載,降低了系統的實時性,尤其對于復雜模型的系統,比較復雜且容易出錯。另外,擴展卡爾曼濾波算法在系統達到穩(wěn)定狀態(tài)時,會喪失對于系統突變的跟蹤能力,將導致狀態(tài)估計不準甚至發(fā)散[11]。針對擴展卡爾曼濾波算法存在的問題,很多文獻引入了高斯QEKF算法、高斯IEKF算法、強跟蹤EKF算法[11]以及各種改進算法[12],但是這些算法都是次優(yōu)濾波。

    本文將自適應算法和無跡卡爾曼濾波相結合,使用確定的Sigma點集來描述非線性函數的統計特性。Sigma點集經過非線性函數傳播之后,后驗均值和方差可以達到2階精度而擴展卡爾曼濾波算法只能達到1階精度。不像擴展卡爾曼濾波算法對于非線性函數的近似,無跡卡爾曼濾波是對狀態(tài)隨機變量的近似,使用了真實非線性模型,由于不需要像二次濾波那樣計算Hession和Jacobian矩陣,提高了解算速度。由于加速度計和羅盤良好的低頻特性,使用加速度以及羅盤得到的姿態(tài)角作為觀測量來修正陀螺儀得到姿態(tài)角,可以有效抑制陀螺儀漂移和量化誤差。得益于自適應無跡卡爾曼濾波器響應速度快的特點,在導航坐標系補償運動加速度對加速度計影響的前提,該算法可以可以進一步減少運動加速度對于加速度計的影響。另外該算法可以在線實時更新系統量測噪聲,所以有更強的魯棒性和更高的解算精度。

    1 非線性姿態(tài)解算系統建模

    本文使用的是NED(北、東、地)右手坐標系。姿態(tài)解算系統離散時間狀態(tài)空間模型如式(1)所示:

    (1)

    DCM矩陣微分方程由于需要求解9個變量,且涉及大量三角運算,會增加處理器運算負載。使用歐拉角微分方程更新姿態(tài)角不需要重正交化,可以簡化運算過程。因此我們使用歐拉角法來預測3個姿態(tài)角,非線性函數f(Xk-1)可以寫為式(2):

    (2)

    式中:w=[wxwywz]T為三軸陀螺儀得到的機體角速度。當等式俯仰角為90°時方程具有奇異性,橫滾角與俯仰角變得不確定。在實際工程中,通常將俯仰角限制在一定的范圍以防止方程退化,如式(3)所示:

    (3)

    觀測方程的量測向量Z可以由MEMS三軸加速度計和磁強計直接計算。式(1)中的觀測方程可以寫為式(4):

    (4)

    由式(5)描述的坐標系變換矩陣以及式(6)可以得到由加速度描述的橫滾角和俯仰角,如式(7)所示:

    (5)

    (6)

    (7)

    式中:Ab=[axayaz]T表示三軸加速度計的輸出值,An=[0 0g]T為重力加速度在參考坐標系下的向量表示。

    偏航角可以由三軸磁強計獲得,如式(8)所示:

    Zψ=arctan(mx/my)+Δψ

    (8)

    式中:Mb=[mxmymz]T為磁強計在x軸、y軸和z軸的輸出值,Δψ為真北與磁北的夾角。機體在轉動過程中由于橫滾角和俯仰角的變化,每次姿態(tài)解算之后需要對mx和my進行修正,如式(9)所示:

    (9)

    2 改進自適應無跡卡爾曼濾波算法

    由上一部分得到的狀態(tài)空間方程,我們可以利用加速度計和磁強計得到的橫滾角、俯仰角以及偏航角來修正由陀螺儀更新后的姿態(tài)角。無跡卡爾曼濾波融合算法可以利用加速度計提供的準確的靜態(tài)角度補償陀螺儀的橫滾角和俯仰角的漂移誤差,通過引入自適應算法在線估計測量噪聲,可以有效降低對于加速度計震動噪聲的影響。下面介紹了完整的改進型無跡卡爾曼濾波算法的解算步驟。

    ①計算Sigma點集,由于狀態(tài)變量為三維,所以取7個Sigma點。

    (10)

    ②計算Sigma點集的權值

    (11)

    式中:下標m表示均值,c表示協方差,上標表示第幾個采樣點。參數λ=a2(n+κ)-n是縮放比例參數,用于降低總的預測誤差。a的選取控制了采樣點的分布狀態(tài),κ為另外一個縮放比例參數,用于保證矩陣(n+λ)P為半正定矩陣。參數β≥0是一個非負的權系數,可以用來合并方程中高階項的動差。

    ③利用式(10)、式(11)獲得Sigma點集及其對應的權值。

    (12)

    ④計算Sigma點集的一步預測。

    (13)

    ⑤計算系統狀態(tài)量的一步預測及協方差矩陣,由Sigma點集的預測值加權求和得到,其中權值由式(11)得到。

    (14)

    (15)

    ⑥根據一步預測值,由UT變換產生新的Sigma點集。

    (16)

    ⑦將⑥式預測的Sigma點集代入觀測方程,得到預測的觀測量。

    (17)

    ⑧將⑦得到的Sigma點集的觀測預測值,通過加權求和得到系統預測的均值及協方差。

    (18)

    (19)

    (20)

    (21)

    (22)

    ⑨計算卡爾曼增益矩陣

    (23)

    ⑩計算系統的狀態(tài)更新和協方差更新

    (24)

    (25)

    βk=βk-1(βk-1+b)-1,β1=I

    (26)

    式中:b=diag([bφbθbψ])為漸消記憶因子0

    由于加速度計容易受到電機震動等的干擾。本文采用了低通FIR濾波器對加速度信號進行了處理,可以有效的抑制震動噪聲?;跐h明窗的線性相位FIR濾波器的傳遞函數在Z域中可以用式(27)表示:

    (27)

    式中:N是濾波器階數,hk是濾波器脈沖響應由截止頻率決定。為了有效的獲取加速度計的截止頻率,本文使用離散傅里葉變換來計算譜振幅,如式(28)所示:

    (28)

    圖1 姿態(tài)解算流程圖

    式中:DFT表示離散傅里葉變換,Macc(f)是頻率(f)的譜振幅,macc(n)表示加速度計的采樣序列。整個姿態(tài)解算流程如圖1所示。

    3 實驗與分析

    為了驗證該姿態(tài)解算算法的有效性,選取了商用姿態(tài)板對算法進行了驗證。該姿態(tài)板集成了MPU-6050芯片包括MEMS三軸加速度計和MEMS三軸角速度計,以及MEMS三軸磁強計HMC5883L芯片和一個氣壓高度計BPM180。姿態(tài)板使用了STM32F103T8作為處理核心,該處理器基于ARM CortexTM-M3內核,通過板子上的8M晶體和PLL,STM32可以運行在最高72M時鐘。

    實驗使用本文提出改進型自適應無跡卡爾曼濾波算法(簡記為:IAUKF)對姿態(tài)板采集的數據進行了姿態(tài)解算并與擴展卡爾曼濾波算法(EKF)以及Xsens 公司的MTi標準姿態(tài)曲線進行了對比。姿態(tài)板采集的9軸數據如圖2所示。3種方法解算得到的橫滾角、俯仰角以及偏航角分別如圖3所示。

    圖2 傳感器采集信號

    從圖3中可以看出,改進型自適應無跡卡爾曼濾波算法相比擴展卡爾曼濾波啊算具有更高的解算精度。為了進一步驗證算法的數據融合效果,使用3D Optical Measurement系統對本文提出的算法(IAUKF)的靜態(tài)(w<5°/s)和動態(tài)(w≥5°/s)RMS(Root-Mean-Square)誤差進行了檢測,實驗結果如表1所示。

    圖3 IAUKF與EKF以及MTi解算效果對比

    表1 靜態(tài)和動態(tài)RMS誤差

    加速度計主要受到電機震動的影響,可以視為白噪聲,為了驗證觀測噪聲的在線更新效果,加入了白噪聲作為干擾項,圖4對比了自適應濾波前后的橫滾角以及俯仰角姿態(tài)解算效果。

    圖4中可以看到,無跡卡爾曼濾波算法和改進型自適應濾波相結合以后在噪聲擾動下更為穩(wěn)定,具有更好的魯棒性。為了驗證本文提出的算法(IAUFK)的收斂速度將初始姿態(tài)與目標姿態(tài)設定偏差并與互補濾波(CF)以及擴展卡爾曼濾波(EKF)做了對比試驗,如圖5所示。

    圖4 噪聲擾動下橫滾角、俯仰角有無自適應對比

    圖5 3種算法的收斂性對比

    從圖中可以看到,擴展卡爾曼濾波算法收斂速度快于互補濾波算法,而本文提出的算法收斂速度比擴展卡爾曼濾波算法更快。圖6展示了互補濾波算法、擴展卡爾曼濾波算法以及本文提出的算法的解算時間。

    圖6 3種算法運行1 000次每次需要的時間

    從圖6中可以看出,傳統的互補濾波算法每次解算需要的時間最短,本文提出的算法和擴展卡爾曼濾波算法需要的時間近似。為了測試IAUKF算法對于陀螺儀漂移的補償效果,將姿態(tài)板初始姿態(tài)置為零態(tài),然后任意旋轉姿態(tài)板,每隔3 min記錄一次零態(tài)姿態(tài)角,補償前后橫滾角、俯仰角對比如圖7所示。

    從圖7可以看出,在算法補償前陀螺儀由于存在常值漂移、速率隨機游走、角度隨機游走以及各種噪聲的影響,漂移嚴重。本文提出的算法可以有效的補償陀螺漂移。

    磁強計容易受到外部磁場干擾,偏航角誤差會隨著外部磁場變化,數據融合算法會對干擾起到一定的補償作用,但是偏差依然存在。實驗中加入一定強度外部磁場,偏航角產生一定誤差,當外部磁場消失以后偏航角迅速得到修正,如圖8所示。

    圖7 算法補償前后橫滾角、俯仰角

    圖8 磁強計收到干擾時的偏航角

    4 結束語

    本文提出的改進型自適應無跡卡爾曼濾波算法可以有效的補償低成本MEMS傳感器存在的姿態(tài)漂移以及噪聲影響,可以有效提高姿態(tài)解算精度。該算法將改進型自適應Sage-Husa算法與無跡卡爾曼濾波算法相結合,可以在線更新狀態(tài)方程的量測誤差,同時避免了使用擴展卡爾曼濾波算法帶來的高階截斷誤差,然后使用加速度計以及磁強計作為量測信息校準陀螺儀解算歐拉角的誤差,提高了基于慣性測量單元解算的準確性與可靠性。試驗與仿真分析表明,該算法收斂速度快、解算速度與擴展卡爾曼濾波算法相當但是可以提供更高的姿態(tài)解算精度,靜態(tài)RMS誤差小于0.75°,動態(tài)RMS誤差小于1.65°,可以滿足對于高精度姿態(tài)控制的要求。但是,相對于傳統互補濾波算法,該算法在解算速度上存在較大差距,接下來實驗室將集中精力對于算法進行優(yōu)化以提高算法的數據融合速度。

    [1] 李迪,陳向堅,續(xù)志軍. 增益自適應滑??刂破髟谖⑿惋w行器飛行姿態(tài)控制中的應用[J]. 光學精密工程,2013,21(5):1183-1191.

    [2] Mahony R,Hamel T,Pflimlin J M. Complementary Filter Design on the Special Orthogonal Group SO(3)[C]//Decision and Control,2005 and 2005 European Control Conference. Cdc-Ecc’05. IEEE Conference on,2006:1477-1484.

    [3] Yanlin Y U,Jiakui Z U,Zeng G. The Research on Attitude Determination Algorithm Based on Complementary Filter for Unmanned Helicopter[J]. Helicopter Technique,2014.

    [4] Wu J,Zhou Z,Chen J,et al. Fast Complementary Filter for Attitude Estimation Using Low-Cost MARG Sensors[J]. IEEE Sensors Journal,2016,16(18):1-1.

    [5] 王立,章政,孫平. 一種自適應互補濾波姿態(tài)估計算法[J]. 控制工程,2015(5):881-886.

    [6] 孫金秋,游有鵬,傅忠云. 基于共軛梯度法和互補濾波相結合的姿態(tài)解算算法[J]. 傳感技術學報,2014(4):524-528.

    [7] Madgwick S O,Harrison A J,Vaidyanathan A. Estimation of IMU and MARG Orientation Using a Gradient Descent Algorithm.[C]//IEEE International Conference on Rehabilitation Robotics. IEEE Int Conf Rehabil Robot,2011:1-7.

    [8] Zhang S,Yu S,Liu C,et al. A Dual-Linear Kalman Filter for Real-Time Orientation Determination System Using Low-Cost MEMS Sensors.[J]. Sensors,2016,16(2):264.

    [9] Wang Y,Wai D,Tomizuka M. Steady-State Marginalized Particle Filter for Attitude Estimation[C]//ASME 2014 Dynamic Systems and Control Conference. American Society of Mechanical Engineers,2014:V002T25A002-V002T25A002.

    [10] Wei G,Li J. Adaptive Kalman Filtering for the Integrated SINS/DVL System[J]. Journal of Computational Information Systems,2013,16(9):6443-6450.

    [11] 宋康寧,叢爽,鄧科,等. 自適應強跟蹤卡爾曼濾波在陀螺穩(wěn)定平臺中的應用[J]. 中國科學技術大學學報,2015(1):17-22.

    [12] Zou Bo,Zhang Hua,Jiang Jun. 多傳感信息融合的改進擴展卡爾曼濾波定姿[J]. 計算機應用研究,2014,31(4):1035-1038.

    [13] 王丁偉,祖家奎,黃海. 基于MAEKF算法的航姿參考系統設計[J]. 傳感技術學報,2017,30(2):319-325.

    [14] Hajiyev C,Soken H E. Robust Adaptive Kalman Filter for Estimation of UAV Dynamics in the Presence of Sensor/Actuator Faults[J]. Aerospace Science and Technology,2013,28(1):376-383.

    [15] Makni A,Fourati H,Kibangou A Y. Adaptive Kalman Filter for MEMS-IMU Based Attitude Estimation Under External Acceleration and Parsimonious use of Gyroscopes[M]. 2014.

    [16] 嚴恭敏,李四海,秦永元. 慣性儀器測試與數據分析[M]. 國防工業(yè)出版社,2012.

    AnImprovedAdaptiveUnscentedKalmanFilterforAttitudeDetermination*

    HOUJiangkuan1,MAJun1*,JIAHuayu2

    (1.College of Physics and Optoelectronics,Taiyuan University of Science and Technology,Jinzhong Shanxi 030600,China; 2.College of Information Engineering,Taiyuan University of Science and Technology,Jinzhong Shanxi 030600,China)

    An improved adaptive unscented Kalman filter for attitude determination is proposed,which can effectively solve the drift and noise problem of gyroscope and reduce the influence of motion acceleration on the accelerometer. Combining the improved adaptive Sage-Husa algorithm with the unscented Kalman filter,the statistic characteristic of the measurement noise is updated on-line,the anti-jamming capability of the system is improved,and the linearization error of the extended Kalman filter can be avoided. And we can get the precise Attitude angle. Only three Euler angles are updated each iteration which can improve the solution speed of the system. The flight experiment and analysis show that the improved algorithm can effectively improve the accuracy of attitude determination,improve convergence speed,strong adaptability,robustness and stability,and can get the accurate attitude angle quickly when the interference disappears.

    sensor;attitude estimation;unscented Kalman filter;adaptive Sage-Husa algorithm;multi-rotor aircraft

    TP27

    A

    1004-1699(2017)10-1518-07

    10.3969/j.issn.1004-1699.2017.10.011

    侯江寬(1991-),男,漢族,河北省石家莊市,碩士研究生,主要研究方向為導航與控制、智能控制,houjiangkuan0126@163.com;

    馬珺(1980-),女,漢族,山西省太原市人,博士研究生,副教授,碩士研究生導師,主要研究方向為新型傳感器技術、智能控制理論,zymajun@126.com;

    賈華宇(1977-),男,漢族,山西臨汾人,博士研究生,副教授,碩士研究生導師,主要研究方向為數?;旌霞呻娐吩O計。

    猜你喜歡
    磁強計無跡陀螺儀
    小小宋慈大智慧·無形無跡的證據
    磁強計陣列測量一致性校正
    基于EMD的MEMS陀螺儀隨機漂移分析方法
    基于矢量磁強計的磁場梯度張量儀誤差校正方法
    組合導航中磁強計干擾估計與補償方法
    基于LabVIEW的微型磁通門磁強計測試系統搭建
    無跡卡爾曼濾波在電線積冰觀測數據處理中的應用
    我國著名陀螺儀專家——林士諤
    基于無跡卡爾曼濾波的行波波頭辨識
    微機械陀螺儀概述和發(fā)展
    久久国产精品人妻蜜桃| 国产69精品久久久久777片| 欧美三级亚洲精品| 干丝袜人妻中文字幕| 欧美精品国产亚洲| 婷婷精品国产亚洲av| 在线免费观看不下载黄p国产 | 日韩精品青青久久久久久| 亚洲成av人片在线播放无| 国产国拍精品亚洲av在线观看| 亚洲av中文字字幕乱码综合| 少妇猛男粗大的猛烈进出视频 | 亚洲五月天丁香| 啦啦啦观看免费观看视频高清| 亚洲av熟女| 此物有八面人人有两片| 亚洲图色成人| 国产黄a三级三级三级人| 女同久久另类99精品国产91| 国产欧美日韩一区二区精品| 久久精品久久久久久噜噜老黄 | 国产单亲对白刺激| 国产精品一区二区三区四区久久| 小蜜桃在线观看免费完整版高清| 日本爱情动作片www.在线观看 | 久久天躁狠狠躁夜夜2o2o| 国内精品美女久久久久久| 欧美中文日本在线观看视频| 琪琪午夜伦伦电影理论片6080| 欧美精品啪啪一区二区三区| 国产单亲对白刺激| 亚洲av日韩精品久久久久久密| 成年女人看的毛片在线观看| 黄色欧美视频在线观看| 亚洲男人的天堂狠狠| 国产高清视频在线观看网站| av.在线天堂| 午夜福利成人在线免费观看| 99久久无色码亚洲精品果冻| 午夜a级毛片| 国产白丝娇喘喷水9色精品| 日韩欧美在线二视频| 夜夜爽天天搞| eeuss影院久久| 99久久精品一区二区三区| x7x7x7水蜜桃| 丝袜美腿在线中文| 国产不卡一卡二| 免费av不卡在线播放| 久久精品国产亚洲av天美| 国产精品久久久久久av不卡| 最新中文字幕久久久久| 99热这里只有精品一区| 久久天躁狠狠躁夜夜2o2o| 亚洲精品成人久久久久久| 久久精品91蜜桃| 欧美潮喷喷水| 精品无人区乱码1区二区| 听说在线观看完整版免费高清| 国产av在哪里看| 老熟妇乱子伦视频在线观看| 欧美国产日韩亚洲一区| 99久久中文字幕三级久久日本| h日本视频在线播放| 亚洲国产精品久久男人天堂| 精品久久久久久久久久久久久| 日本色播在线视频| 赤兔流量卡办理| 亚洲欧美清纯卡通| 国产av不卡久久| 国产午夜精品久久久久久一区二区三区 | 日韩人妻高清精品专区| 中文在线观看免费www的网站| 可以在线观看的亚洲视频| 亚洲国产欧美人成| 精品人妻偷拍中文字幕| 国产免费男女视频| 欧美中文日本在线观看视频| 深爱激情五月婷婷| 欧美又色又爽又黄视频| 看免费成人av毛片| 欧美高清性xxxxhd video| 欧美xxxx黑人xx丫x性爽| 国国产精品蜜臀av免费| 99久久精品热视频| 国产精品一区二区性色av| 国产精品一区二区免费欧美| 午夜精品久久久久久毛片777| 狠狠狠狠99中文字幕| 美女高潮喷水抽搐中文字幕| 国产一区二区三区av在线 | 韩国av在线不卡| 九色成人免费人妻av| 91久久精品电影网| 国产综合懂色| 久99久视频精品免费| 国产精品爽爽va在线观看网站| 国产男靠女视频免费网站| 少妇的逼水好多| 久久久久性生活片| 欧美成人性av电影在线观看| 国产精品av视频在线免费观看| 成人av一区二区三区在线看| 观看免费一级毛片| 嫩草影院入口| 国产男人的电影天堂91| av天堂在线播放| 又紧又爽又黄一区二区| 欧美最黄视频在线播放免费| 麻豆精品久久久久久蜜桃| 日日干狠狠操夜夜爽| 国产精品一及| 国产高潮美女av| 国产成人一区二区在线| 国产av一区在线观看免费| 亚洲欧美激情综合另类| 久久精品久久久久久噜噜老黄 | 窝窝影院91人妻| 国产麻豆成人av免费视频| 久久精品国产亚洲av香蕉五月| 国产极品精品免费视频能看的| 久久久久久久亚洲中文字幕| 欧美色视频一区免费| 欧美日韩精品成人综合77777| 欧美+日韩+精品| 别揉我奶头 嗯啊视频| 国产亚洲精品久久久com| 久久亚洲精品不卡| 国内少妇人妻偷人精品xxx网站| 老司机深夜福利视频在线观看| 精品久久久久久久久久久久久| 中亚洲国语对白在线视频| 一进一出好大好爽视频| 亚洲专区中文字幕在线| 午夜精品一区二区三区免费看| 欧美xxxx性猛交bbbb| 亚洲精华国产精华精| 亚洲精品国产成人久久av| 国产真实伦视频高清在线观看 | 精品一区二区三区av网在线观看| 最新在线观看一区二区三区| 精品午夜福利在线看| 免费电影在线观看免费观看| 亚洲第一区二区三区不卡| 精品国内亚洲2022精品成人| 狂野欧美白嫩少妇大欣赏| 偷拍熟女少妇极品色| 色吧在线观看| av中文乱码字幕在线| 婷婷精品国产亚洲av| av专区在线播放| 一进一出抽搐动态| av视频在线观看入口| 免费高清视频大片| 91久久精品电影网| 成人特级黄色片久久久久久久| 欧美+日韩+精品| 中文字幕av在线有码专区| 尤物成人国产欧美一区二区三区| 久久九九热精品免费| 深爱激情五月婷婷| 搡老岳熟女国产| 精品一区二区三区人妻视频| 欧美精品啪啪一区二区三区| 国产在线精品亚洲第一网站| 欧美成人免费av一区二区三区| 久久九九热精品免费| 一区二区三区四区激情视频 | 熟女人妻精品中文字幕| 欧美zozozo另类| 日本 欧美在线| 国产亚洲av嫩草精品影院| 免费观看人在逋| 久久中文看片网| 搡老岳熟女国产| 九九在线视频观看精品| 国产精品不卡视频一区二区| 又黄又爽又刺激的免费视频.| 99riav亚洲国产免费| 欧美性感艳星| 亚洲内射少妇av| 三级国产精品欧美在线观看| 日韩精品中文字幕看吧| 色尼玛亚洲综合影院| bbb黄色大片| 亚洲av免费高清在线观看| 日韩欧美精品免费久久| 亚洲精品色激情综合| 亚洲成人免费电影在线观看| 成人无遮挡网站| 国产精品久久久久久精品电影| 少妇的逼水好多| 亚洲图色成人| 春色校园在线视频观看| 18禁在线播放成人免费| 国产精品av视频在线免费观看| 午夜免费成人在线视频| 亚洲av五月六月丁香网| 亚洲欧美日韩高清在线视频| 欧美中文日本在线观看视频| 性色avwww在线观看| 97碰自拍视频| 欧美+日韩+精品| 国内精品久久久久精免费| 亚洲成a人片在线一区二区| 国产人妻一区二区三区在| 岛国在线免费视频观看| 久久久精品欧美日韩精品| 色综合站精品国产| 国产蜜桃级精品一区二区三区| 亚洲人成网站在线播放欧美日韩| 亚洲美女搞黄在线观看 | 国内揄拍国产精品人妻在线| 搞女人的毛片| 亚洲av成人精品一区久久| 国产单亲对白刺激| 精品欧美国产一区二区三| 国产成人影院久久av| 日本免费a在线| 国产黄a三级三级三级人| 有码 亚洲区| 国产精品免费一区二区三区在线| 国产精品av视频在线免费观看| 国产色婷婷99| 少妇的逼好多水| 在线观看一区二区三区| 日本五十路高清| 午夜福利18| 国产av在哪里看| 又粗又爽又猛毛片免费看| 成人av一区二区三区在线看| 99精品在免费线老司机午夜| 日本一本二区三区精品| 亚洲第一电影网av| 久久精品国产99精品国产亚洲性色| 非洲黑人性xxxx精品又粗又长| 国产日本99.免费观看| 午夜爱爱视频在线播放| 99精品久久久久人妻精品| 日本欧美国产在线视频| 久久九九热精品免费| 国内精品宾馆在线| 日韩欧美精品免费久久| 午夜老司机福利剧场| 男女视频在线观看网站免费| 超碰av人人做人人爽久久| 国产黄a三级三级三级人| 国产精品av视频在线免费观看| 久久精品久久久久久噜噜老黄 | 人妻少妇偷人精品九色| 国产精品精品国产色婷婷| 天堂√8在线中文| 免费电影在线观看免费观看| 国产主播在线观看一区二区| 亚洲av中文av极速乱 | 不卡视频在线观看欧美| 亚洲一区高清亚洲精品| 国模一区二区三区四区视频| 免费看a级黄色片| 在线观看av片永久免费下载| 国产亚洲精品久久久久久毛片| 国产精品乱码一区二三区的特点| av专区在线播放| 亚洲三级黄色毛片| 久久久久久久亚洲中文字幕| 成人美女网站在线观看视频| 亚洲第一电影网av| 精品日产1卡2卡| 久久精品国产鲁丝片午夜精品 | 免费高清视频大片| 黄片wwwwww| 国产高清激情床上av| 国产精品无大码| 国产三级中文精品| 精品一区二区三区视频在线观看免费| 天天一区二区日本电影三级| xxxwww97欧美| 最新在线观看一区二区三区| 啦啦啦啦在线视频资源| 91麻豆av在线| 亚洲不卡免费看| 无人区码免费观看不卡| 亚洲一区高清亚洲精品| 非洲黑人性xxxx精品又粗又长| 欧美色欧美亚洲另类二区| 热99re8久久精品国产| 性欧美人与动物交配| 精品国内亚洲2022精品成人| 色综合色国产| www日本黄色视频网| 免费大片18禁| 国产69精品久久久久777片| 久久久久久久久中文| 国产高清视频在线播放一区| 国产精品国产三级国产av玫瑰| 亚洲精品久久国产高清桃花| 欧美成人免费av一区二区三区| 在线观看舔阴道视频| 久久精品人妻少妇| 天堂动漫精品| 午夜免费激情av| 毛片女人毛片| 男人的好看免费观看在线视频| 国产欧美日韩一区二区精品| 一级黄片播放器| 免费在线观看成人毛片| 两人在一起打扑克的视频| 日本黄色片子视频| 久久久久久久久久久丰满 | 在线观看一区二区三区| 国产伦精品一区二区三区视频9| 一进一出好大好爽视频| 国产一区二区三区在线臀色熟女| 狠狠狠狠99中文字幕| 最近在线观看免费完整版| 99热这里只有是精品在线观看| 亚洲久久久久久中文字幕| 免费观看精品视频网站| 国产综合懂色| .国产精品久久| 亚洲欧美激情综合另类| 性色avwww在线观看| 中亚洲国语对白在线视频| 麻豆成人午夜福利视频| 久久99热6这里只有精品| 极品教师在线视频| 又爽又黄无遮挡网站| 99久久成人亚洲精品观看| 一级av片app| 欧美高清性xxxxhd video| 精品国内亚洲2022精品成人| 中文字幕高清在线视频| 亚洲精华国产精华液的使用体验 | 精品久久久久久久久亚洲 | www.色视频.com| 亚洲一级一片aⅴ在线观看| 少妇的逼好多水| 变态另类丝袜制服| 校园人妻丝袜中文字幕| 精品欧美国产一区二区三| 中出人妻视频一区二区| 欧美另类亚洲清纯唯美| 十八禁国产超污无遮挡网站| 韩国av一区二区三区四区| 国产探花极品一区二区| 在线观看免费视频日本深夜| 日韩中文字幕欧美一区二区| 免费看a级黄色片| 国产蜜桃级精品一区二区三区| 久久精品影院6| 狠狠狠狠99中文字幕| 国产精品嫩草影院av在线观看 | 久久99热这里只有精品18| 日本黄大片高清| 搡老岳熟女国产| 欧美黑人欧美精品刺激| 久久中文看片网| 乱人视频在线观看| 女人十人毛片免费观看3o分钟| 日本爱情动作片www.在线观看 | 久久午夜福利片| 一区二区三区免费毛片| 国产视频内射| 久久精品国产亚洲网站| 内地一区二区视频在线| 国产精品久久电影中文字幕| 国产精品一区二区免费欧美| 国产精品久久电影中文字幕| 中文字幕久久专区| 国产精品电影一区二区三区| 老司机深夜福利视频在线观看| 国产精品福利在线免费观看| 日韩欧美免费精品| 午夜福利成人在线免费观看| 欧美激情国产日韩精品一区| 久久这里只有精品中国| 免费人成视频x8x8入口观看| 中文在线观看免费www的网站| 此物有八面人人有两片| 国产精品精品国产色婷婷| 成人亚洲精品av一区二区| 白带黄色成豆腐渣| 他把我摸到了高潮在线观看| 春色校园在线视频观看| 亚洲国产精品久久男人天堂| 成人特级黄色片久久久久久久| 国产亚洲91精品色在线| 99热这里只有是精品在线观看| 久久久久国产精品人妻aⅴ院| 男女做爰动态图高潮gif福利片| 国产亚洲精品久久久com| 看黄色毛片网站| 深夜a级毛片| 欧美黑人巨大hd| 精品人妻一区二区三区麻豆 | 美女 人体艺术 gogo| 国产伦一二天堂av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧美日韩高清专用| 国产精品人妻久久久久久| 亚洲成人精品中文字幕电影| 99国产精品一区二区蜜桃av| 蜜桃久久精品国产亚洲av| 亚洲成人免费电影在线观看| 亚洲熟妇熟女久久| 国产69精品久久久久777片| 久久香蕉精品热| 成人特级av手机在线观看| 如何舔出高潮| 免费观看在线日韩| 亚洲av日韩精品久久久久久密| 校园春色视频在线观看| 麻豆久久精品国产亚洲av| 露出奶头的视频| 如何舔出高潮| 日本黄色视频三级网站网址| 国产精品,欧美在线| 成年人黄色毛片网站| 搡老熟女国产l中国老女人| 国产精品一及| 51国产日韩欧美| 不卡视频在线观看欧美| 久久久成人免费电影| 18禁黄网站禁片午夜丰满| 国语自产精品视频在线第100页| 久久久久久久久久成人| 久久精品国产亚洲av香蕉五月| 成人性生交大片免费视频hd| 欧美精品啪啪一区二区三区| 亚洲av免费在线观看| 国产中年淑女户外野战色| 亚洲专区国产一区二区| 丰满乱子伦码专区| 日韩欧美在线二视频| 国产精品亚洲一级av第二区| 色哟哟·www| 国产精品伦人一区二区| 亚洲国产高清在线一区二区三| 亚洲精品粉嫩美女一区| 成人美女网站在线观看视频| 国产 一区 欧美 日韩| 日韩欧美国产一区二区入口| 观看美女的网站| 午夜精品一区二区三区免费看| 精品久久久久久,| av女优亚洲男人天堂| 婷婷亚洲欧美| 久久精品综合一区二区三区| 日韩欧美在线二视频| 婷婷精品国产亚洲av在线| 97热精品久久久久久| 精品久久久久久久人妻蜜臀av| 热99在线观看视频| 国产中年淑女户外野战色| 女同久久另类99精品国产91| 欧美一区二区亚洲| 午夜久久久久精精品| 神马国产精品三级电影在线观看| 亚洲精华国产精华精| 国产在线男女| 日韩精品有码人妻一区| 在现免费观看毛片| 人妻制服诱惑在线中文字幕| 有码 亚洲区| 欧美国产日韩亚洲一区| 国产午夜精品久久久久久一区二区三区 | 精品久久久久久久久久免费视频| 无遮挡黄片免费观看| 日本欧美国产在线视频| 欧美不卡视频在线免费观看| a级毛片免费高清观看在线播放| 国产精品久久久久久久久免| 久久久久久九九精品二区国产| 欧美成人a在线观看| av在线亚洲专区| 国产av不卡久久| 国产av一区在线观看免费| 国产三级在线视频| av女优亚洲男人天堂| 在线观看66精品国产| 午夜亚洲福利在线播放| 国产伦精品一区二区三区四那| 日韩欧美一区二区三区在线观看| 夜夜夜夜夜久久久久| 婷婷精品国产亚洲av在线| 国产麻豆成人av免费视频| 一区二区三区高清视频在线| 在现免费观看毛片| 亚洲第一区二区三区不卡| 国产乱人伦免费视频| 国内揄拍国产精品人妻在线| 麻豆av噜噜一区二区三区| 国产女主播在线喷水免费视频网站 | 尤物成人国产欧美一区二区三区| 免费看a级黄色片| 麻豆国产av国片精品| 国产精品1区2区在线观看.| 免费观看人在逋| 国产成人一区二区在线| 久久人人爽人人爽人人片va| 欧洲精品卡2卡3卡4卡5卡区| 在线观看av片永久免费下载| 久久久久国内视频| 少妇高潮的动态图| 日韩欧美免费精品| 波多野结衣巨乳人妻| 日本精品一区二区三区蜜桃| 九九爱精品视频在线观看| 亚洲乱码一区二区免费版| 99在线视频只有这里精品首页| 女同久久另类99精品国产91| 69av精品久久久久久| 在线播放无遮挡| 精品久久国产蜜桃| 日韩大尺度精品在线看网址| 久久久久久伊人网av| 国产精品伦人一区二区| av天堂在线播放| 别揉我奶头~嗯~啊~动态视频| 精品福利观看| 国产成人福利小说| 国产av一区在线观看免费| 五月伊人婷婷丁香| 中文字幕精品亚洲无线码一区| 国产亚洲精品久久久久久毛片| 成人午夜高清在线视频| 国产成人一区二区在线| 亚洲国产欧美人成| 精品久久久久久久久久久久久| 精品久久久久久,| 狂野欧美白嫩少妇大欣赏| 国产高清不卡午夜福利| 在线观看美女被高潮喷水网站| 在线播放无遮挡| 亚洲狠狠婷婷综合久久图片| 午夜视频国产福利| 99在线人妻在线中文字幕| 亚洲人成网站高清观看| 男女做爰动态图高潮gif福利片| 干丝袜人妻中文字幕| 日韩精品青青久久久久久| 国产真实乱freesex| av在线蜜桃| 欧美色视频一区免费| 亚洲电影在线观看av| 高清毛片免费观看视频网站| h日本视频在线播放| 岛国在线免费视频观看| 一个人看视频在线观看www免费| 欧美性感艳星| 在线播放无遮挡| 亚洲无线在线观看| 成人精品一区二区免费| 有码 亚洲区| 窝窝影院91人妻| 欧美黑人欧美精品刺激| 国产aⅴ精品一区二区三区波| 日韩精品有码人妻一区| 免费黄网站久久成人精品| 日韩,欧美,国产一区二区三区 | 免费观看人在逋| 国产伦一二天堂av在线观看| 最近视频中文字幕2019在线8| 日韩欧美在线乱码| 国产一区二区在线av高清观看| 99热精品在线国产| 有码 亚洲区| 国产精品野战在线观看| 久久久久久国产a免费观看| 日日夜夜操网爽| 久久久精品欧美日韩精品| 免费av观看视频| 天美传媒精品一区二区| 网址你懂的国产日韩在线| 天天躁日日操中文字幕| 精品不卡国产一区二区三区| 18禁裸乳无遮挡免费网站照片| 夜夜看夜夜爽夜夜摸| 日韩亚洲欧美综合| 永久网站在线| 真人一进一出gif抽搐免费| 国产一级毛片七仙女欲春2| 麻豆成人午夜福利视频| 97碰自拍视频| 成人鲁丝片一二三区免费| 亚洲国产精品sss在线观看| 99久国产av精品| 搡女人真爽免费视频火全软件 | 亚洲精品久久国产高清桃花| 亚洲成人免费电影在线观看| 好男人在线观看高清免费视频| 国产高清视频在线播放一区| 午夜激情欧美在线| av在线观看视频网站免费| 免费看日本二区| avwww免费| 精品人妻一区二区三区麻豆 | av.在线天堂| 亚洲国产色片| 午夜福利在线在线| 久久99热6这里只有精品| 五月伊人婷婷丁香| 尤物成人国产欧美一区二区三区| 88av欧美| a级一级毛片免费在线观看| 91久久精品国产一区二区三区| 亚洲国产日韩欧美精品在线观看| 观看美女的网站| 亚洲精品一卡2卡三卡4卡5卡| 禁无遮挡网站| 国产精品综合久久久久久久免费| 精品无人区乱码1区二区| 成人三级黄色视频| 亚洲国产精品久久男人天堂| 亚洲av免费在线观看|