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

    基于UKF的GPS/SINS組合導(dǎo)航直接式濾波算法

    2015-03-11 05:32:42林雪原
    導(dǎo)航定位與授時(shí) 2015年3期
    關(guān)鍵詞:導(dǎo)航系統(tǒng)卡爾曼濾波慣性

    王 林,林雪原

    (1.海軍航空工程學(xué)院 研究生管理大隊(duì),山東 煙臺(tái) 264000;2.信息融合研究所,山東 煙臺(tái) 264000)

    ?

    基于UKF的GPS/SINS組合導(dǎo)航直接式濾波算法

    王 林1,林雪原2

    (1.海軍航空工程學(xué)院 研究生管理大隊(duì),山東 煙臺(tái) 264000;2.信息融合研究所,山東 煙臺(tái) 264000)

    基于無(wú)跡卡爾曼濾波(UKF)方法,使用姿態(tài)、速度、位置等9個(gè)導(dǎo)航參數(shù)組成狀態(tài)向量,以GPS系統(tǒng)輸出的速度、位置組成6維觀測(cè)向量,構(gòu)建直接式結(jié)構(gòu)的UKF濾波器。該濾波器能夠直接反映系統(tǒng)導(dǎo)航參數(shù)的動(dòng)態(tài)過(guò)程,準(zhǔn)確顯示運(yùn)動(dòng)狀態(tài)演變。針對(duì)GPS/SINS組合導(dǎo)航系統(tǒng)的特點(diǎn),構(gòu)建了GPS/SINS組合導(dǎo)航直接式卡爾曼濾波仿真驗(yàn)證系統(tǒng),仿真結(jié)果驗(yàn)證了基于UKF的GPS/SINS組合導(dǎo)航直接式濾波算法的有效性,該直接式非線性濾波算法可使慣性組合導(dǎo)航系統(tǒng)的導(dǎo)航精度得到提高。

    慣性組合導(dǎo)航;無(wú)跡卡爾曼濾波;直接式濾波

    0 引言

    卡爾曼濾波是一種基于模型的線性最小方差估計(jì),其標(biāo)準(zhǔn)離散型算法具有遞推計(jì)算、使用計(jì)算機(jī)執(zhí)行、適用于平穩(wěn)或非平穩(wěn)多維隨機(jī)信號(hào)估計(jì)等獨(dú)特優(yōu)點(diǎn)。所以在隨機(jī)信號(hào)處理中,特別是在組合導(dǎo)航設(shè)計(jì)中應(yīng)用廣泛。根據(jù)卡爾曼濾波方法設(shè)計(jì)的組合導(dǎo)航濾波器,成為組合導(dǎo)航系統(tǒng)的核心。標(biāo)準(zhǔn)卡爾曼濾波只適用于系統(tǒng)方程和量測(cè)均為線性時(shí)的估計(jì),但是工程實(shí)踐中所遇到的物理系統(tǒng)數(shù)學(xué)模型往往是非線性的,即系統(tǒng)方程是非線性的,或者系統(tǒng)和量測(cè)都是非線性的。盡管擴(kuò)展卡爾曼濾波(EKF)可解決系統(tǒng)和量測(cè)為非線性時(shí)的估計(jì),但必須用泰勒級(jí)數(shù)的方法對(duì)原系統(tǒng)和量測(cè)作展開(kāi)并舍去非線性項(xiàng),再用基礎(chǔ)卡爾曼濾波算法對(duì)處理后的系統(tǒng)和量測(cè)方程作濾波計(jì)算,可見(jiàn)EKF本質(zhì)上仍然是基礎(chǔ)卡爾曼濾波。由于在線性化處理過(guò)程中舍棄了非線性的高階項(xiàng),EKF只適用于濾波對(duì)象非線性較弱的情況下,被估計(jì)對(duì)象非線性性質(zhì)越強(qiáng),所產(chǎn)生的誤差就會(huì)越大,嚴(yán)重的會(huì)引起濾波發(fā)散[1-2]。

    為了解決強(qiáng)非線性條件下的估計(jì)問(wèn)題,出現(xiàn)了基于UT變換的UKF算法,為以直接式濾波方法解決強(qiáng)非線性條件下的慣性組合導(dǎo)航提供了契機(jī)。

    本文將UKF應(yīng)用于組合導(dǎo)航,建立慣性組合導(dǎo)航系統(tǒng)的直接式濾波模型,以導(dǎo)航參數(shù)作為狀態(tài)變量,以GPS系統(tǒng)的速度和位置組成觀測(cè)向量,實(shí)現(xiàn)直接式濾波關(guān)鍵計(jì)算流程,以微小型飛行器作為對(duì)象進(jìn)行仿真驗(yàn)證。

    1 UKF濾波算法

    UKF和標(biāo)準(zhǔn)卡爾曼濾波都是線性的最小方差估計(jì),都是建立在模型基礎(chǔ)上的算法,兩者的不同點(diǎn)在于計(jì)算最佳增益陣的方法。標(biāo)準(zhǔn)卡爾曼濾波在計(jì)算最佳增益陣時(shí),建立在系統(tǒng)方程和量測(cè)方程均為線性的基礎(chǔ)上。而UKF在最佳增益陣的計(jì)算中,對(duì)系統(tǒng)方程以及量測(cè)方程的性質(zhì)沒(méi)有任何具體的約束,所以UKF算法既可以在線性條件下使用,也可以在非線性條件下使用[2]。

    設(shè)系統(tǒng)方程和量測(cè)方程為:

    (1)

    其中k-1,k表示離散時(shí)刻,Xk為狀態(tài)預(yù)測(cè),Xk-1為上一時(shí)刻的狀態(tài)估計(jì),其方差為Pk-1,uk-1為系統(tǒng)確定性控制項(xiàng),Wk-1為系統(tǒng)噪聲,其方差為Qk-1,Zk為量測(cè)值,Vk為量測(cè)噪聲,其方差為Rk,系統(tǒng)和量測(cè)噪聲均服從高斯分布。UKF濾波要解決的問(wèn)題是在已知狀態(tài)初始值X0,系統(tǒng)確定性控制序列u0,u1,…,uk-1以及量測(cè)值序列Z1,Z2,…,Zk的情況下,估計(jì)系統(tǒng)的狀態(tài)變量Xk。

    UKF采用遞推計(jì)算方式,主要包括兩個(gè)過(guò)程:狀態(tài)預(yù)測(cè)計(jì)算和修正計(jì)算。已知狀態(tài)變量的維數(shù)為n,UKF用UT變換的方式獲取2n+1個(gè)采樣點(diǎn)(sigma點(diǎn)),實(shí)現(xiàn)狀態(tài)及其相關(guān)統(tǒng)計(jì)量的遞推計(jì)算。在k時(shí)刻,2n+1個(gè)對(duì)稱的采樣sigma點(diǎn)及各點(diǎn)對(duì)應(yīng)的權(quán)值為[1-2]:

    (2)

    UKF算法的核心是UT變換,它能夠產(chǎn)生代表狀態(tài)量統(tǒng)計(jì)特性的最少數(shù)量的采樣點(diǎn),具有一些比較好的特性:UKF方法中的狀態(tài)變量傳遞過(guò)程中通過(guò)的是沒(méi)有被近似處理的系統(tǒng),精度較高,模型較精確;UKF在保證精度的條件下,其計(jì)算量相比粒子濾波(PF)要少很多。這些特點(diǎn)使得UKF成為一種重要的非線性濾波技術(shù)[1]。

    2 基于UKF的慣性組合導(dǎo)航直接式濾波模型

    基于標(biāo)準(zhǔn)卡爾曼濾波及擴(kuò)展卡爾曼濾波的間接式組合導(dǎo)航系統(tǒng),使用線性化的慣性組合導(dǎo)航系統(tǒng)的誤差方程作為狀態(tài)方程,以慣性組合導(dǎo)航系統(tǒng)的誤差作為被估狀態(tài)。而基于UKF的GPS/SINS組合導(dǎo)航直接式濾波則以捷聯(lián)慣性導(dǎo)航系統(tǒng)的機(jī)械編排方程作為狀態(tài)方程,直接使用導(dǎo)航參數(shù)做被估狀態(tài),是一種直接式濾波方法,狀態(tài)方程及量測(cè)方程無(wú)需線性化。

    2.1 狀態(tài)方程與量測(cè)方程

    采用“東北天”地理坐標(biāo)系作為導(dǎo)航坐標(biāo)系,慣性導(dǎo)航系統(tǒng)采用指北方位機(jī)械編排,狀態(tài)方程如式(3)所示。

    (3)

    (4)

    濾波過(guò)程中,噪聲進(jìn)入非線性系統(tǒng)中,參與系統(tǒng)狀態(tài)轉(zhuǎn)移,所以如式(4)所示將系統(tǒng)噪聲和量測(cè)噪聲與狀態(tài)向量一起組成增廣向量Xa,增廣向量維數(shù)n為21,增廣的方差陣為Pa。

    狀態(tài)方程依據(jù)慣性導(dǎo)航系統(tǒng)工作原理展開(kāi)式如式(5)所示。

    (5)

    基于UKF的直接式卡爾曼濾波的觀測(cè)方程直接反映觀測(cè)系統(tǒng)的量測(cè)原理,以GPS輸出的速度和位置作為觀測(cè)量,量測(cè)方程為式(6)所示,V為量測(cè)噪聲。

    (6)

    2.3 基于UKF的慣性組合導(dǎo)航直接式濾波流程

    基于UKF的慣性組合導(dǎo)航直接式濾波系統(tǒng)以狀態(tài)向量、系統(tǒng)噪聲與量測(cè)噪聲組成增廣向量,以GPS輸出的速度和位置作為觀測(cè)量,其濾波流程分為三個(gè)部分:初始化及樣本點(diǎn)計(jì)算、時(shí)間更新、狀態(tài)更新。

    2.3.1 初始化及樣本點(diǎn)計(jì)算

    2.3.2 時(shí)間更新

    2.3.3 狀態(tài)更新

    獲得測(cè)量值及預(yù)測(cè)值后,進(jìn)行新息計(jì)算,通過(guò)濾波增益矩陣,得到每個(gè)狀態(tài)對(duì)應(yīng)的誤差修正量,UKF中新息及濾波增益的計(jì)算都使用到對(duì)狀態(tài)預(yù)測(cè)點(diǎn)集、測(cè)量預(yù)測(cè)點(diǎn)集均值和方差的計(jì)算。

    通過(guò)sigma點(diǎn)集得到的狀態(tài)預(yù)測(cè)均值和方差為

    估計(jì)觀測(cè)值的均值和方差為

    計(jì)算狀態(tài)預(yù)測(cè)與量測(cè)估計(jì)的協(xié)方差為

    計(jì)算濾波增益

    通過(guò)增益矩陣,獲得各個(gè)狀態(tài)的誤差修正量,狀態(tài)和方差更新為

    3 基于UKF的慣性組合導(dǎo)航直接式濾波仿真驗(yàn)證

    以微小型飛行器作為使用載體構(gòu)建圖1所示的濾波仿真驗(yàn)證系統(tǒng),采用飛行器動(dòng)態(tài)仿真和真實(shí)的導(dǎo)航傳感器測(cè)量誤差數(shù)據(jù)合并生成動(dòng)態(tài)數(shù)據(jù)的方法,慣性元件和GPS誤差數(shù)據(jù)從實(shí)驗(yàn)室實(shí)際測(cè)試的數(shù)據(jù)中注入。飛行軌跡按微小型飛行器性能特點(diǎn)及使用需求設(shè)定,包括如下過(guò)程:加速—爬升—平直飛行—爬升—平直飛行—轉(zhuǎn)彎—平直飛行—再一次爬升—平直飛行—轉(zhuǎn)彎—下滑—盤旋狀態(tài),如圖2。仿真時(shí)間為3600s。

    圖1 UKF直接式濾波仿真系統(tǒng)結(jié)構(gòu)圖Fig.1 Structure chart of the UKF direct filtering demonstration system

    圖2 飛行仿真的航跡Fig .2 Simulation of the flight trace

    微小型飛行器的初始姿態(tài)定為:飛行器載體水平,航向90°,位置(118°,29°,50m),UKF濾波初始參數(shù):慣性器件噪聲參數(shù),狀態(tài)初值和初始方差陣的值,UKF中的sigma點(diǎn)伸縮因子等。根據(jù)目前在用的慣性器件的性能,設(shè)定方差值,設(shè)陀螺噪聲的均方根為0.5(°)/s,設(shè)加速度計(jì)噪聲的均方根為10-3g,設(shè)GPS測(cè)速誤差為0.5m/s,定位誤差8m。

    基于UKF的GPS/SINS組合導(dǎo)航仿真系統(tǒng)得到的姿態(tài)曲線如圖3所示,姿態(tài)誤差如圖4所示,定位誤差分別如圖5所示,系統(tǒng)的穩(wěn)態(tài)誤差值如表1所示。

    圖3 姿態(tài)曲線Fig.3 The curve of attitude curve

    圖4 姿態(tài)誤差曲線Fig.4 The curve of attitude error

    圖5 位置誤差曲線Fig.5 The curve of position error

    誤差項(xiàng)姿態(tài)角/(°)橫滾俯仰航向最大值0.33410.29575.5703最小值-0.2491-0.1106-3.5843均 值-0.00090.0000-0.0422均方差0.02380.02920.5561誤差項(xiàng)位置/m經(jīng)度緯度高度最大值6.99696.28757.5442最小值-6.1195-6.1117-7.7811均 值0.05430.0562-0.0365均方差1.92251.92301.3565

    由仿真圖像及數(shù)據(jù)表明,基于UKF的慣性組合導(dǎo)航直接式濾波方法穩(wěn)定,姿態(tài)誤差和位置誤差較小,具有較高的導(dǎo)航精度。在微小型飛行器特性的飛行條件下,該直接式濾波方法能夠適應(yīng)于各種飛行狀態(tài)。

    4 結(jié)論

    本文研究了基于UKF的直接式慣性組合導(dǎo)航濾波方法,建立直接式結(jié)構(gòu)的卡爾曼濾波器,以導(dǎo)航參數(shù)作為狀態(tài)變量,以GPS系統(tǒng)的速度和位置組成觀測(cè)向量,實(shí)現(xiàn)了直接式濾波計(jì)算流程。利用方向余弦矩陣和姿態(tài)角間的關(guān)系,實(shí)現(xiàn)了狀態(tài)向量中姿態(tài)角sigma點(diǎn)的生成與各sigma點(diǎn)姿態(tài)角誤差均值及方差的計(jì)算。本文構(gòu)建了GPS/SINS組合導(dǎo)航直接式卡爾曼濾波仿真驗(yàn)證系統(tǒng),仿真結(jié)果驗(yàn)證了該直接式非線性濾波方法的有效性,采用此方法可有效避免將非線性濾波模型進(jìn)行線性化后的誤差,適于解決高機(jī)動(dòng)條件下飛行器濾波及狀態(tài)估計(jì)。

    [1] 劉建業(yè),曾慶華,趙偉,等.導(dǎo)航系統(tǒng)理論與應(yīng)用[M].西安:西北工業(yè)大學(xué)出版社,2010:296-366.

    [2] 秦永元,張洪鉞,汪叔華.卡爾曼濾波與組合導(dǎo)航原理[M].西安:西北工業(yè)大學(xué)出版社,1998:198-221.

    [3] 高社生,何鵬舉,楊波等.組合導(dǎo)航原理及應(yīng)用[M].西安:西北工業(yè)大學(xué)出版社,2012:92-108.

    [4] 林雪原,劉華章,史佩.一種低成本、低精度SINS/GPS組合導(dǎo)航系統(tǒng)級(jí)試驗(yàn)研究[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2006(10).

    [5] 林雪原,林炳榮.基于量測(cè)多尺度預(yù)處理的GPS/SINS組合導(dǎo)航系統(tǒng)[J].兵工自動(dòng)化,2010(7).

    [6] 林雪原.GPS/SINS組合導(dǎo)航系統(tǒng)與實(shí)現(xiàn)[J].計(jì)算機(jī)測(cè)量與控制,2007(12).

    [7] 林雪原,駱卉子,范文強(qiáng).一種改進(jìn)的GPS動(dòng)態(tài)定位濾波方法[J].海軍航空工程學(xué)院學(xué)報(bào),2010(4).

    [8] 林雪原.GPS/SINS組合導(dǎo)航系統(tǒng)的多尺度融合算法研究[J].電子科技大學(xué)學(xué)報(bào),2011(9).

    [9] 陸海勇.捷聯(lián)慣性導(dǎo)航系統(tǒng)中UKF濾波技術(shù)的應(yīng)用研究[D].南京:碩士學(xué)位論文,南京航空航天大學(xué), 2008.

    [10] 周菊華.UKF濾波方法在捷聯(lián)慣導(dǎo)系統(tǒng)中的應(yīng)用研究[D].南京:碩士學(xué)位論文,東南大學(xué),2009.

    [11] 徐佳.INS/GPS組合導(dǎo)航系統(tǒng)濾波算法的研究[D].沈陽(yáng):碩士學(xué)位論文,東北大學(xué),2006.

    [12] 劉錚.UKF算法及其改進(jìn)算法的研究[D].長(zhǎng)沙:碩士學(xué)位論文,中南大學(xué),2009.

    [13] 徐佳鶴,井龍偉.UKF濾波算法在INS/GPS組合導(dǎo)航系統(tǒng)中的研究[C]//.中國(guó)控制與決策學(xué)術(shù)年會(huì)論文集, 2007: 213-216.

    [14] S J Julier, J K Uhlmann. Unscented Filtering and nonlinear Estimation[A]. Proceeding of the IEEE[C]//.vol.92 No3, March 2004.

    Direct Kalman Filtering Algorithm for the GPS/SINS Integrated Navigation Based on UKF

    WANG Lin1,LIN Xue-yuan2

    (1.Naval Aeronautical and Astronautical University Graduate Students’ Brigade,Yantai 264000,China;2.Research Institute of Information Fusion, Yantai 264000,China)

    A Kalman filter with direct structure is proposed based on UKF filter. The direct filter takes 9 parameters as state vectors such as attitude, velocity and position. And the six-dimension observation vectors are made up of velocity and position provided by Global Positioning System. The filter directly describes navigation parameters dynamic process of the system and reflects the evolution of motion state. A direct filtering demonstration system for the GPS/SINS Integrated Navigation is constructed according to the characteristics of GPS/SINS integrated navigation system. The simulation results verifies the effectiveness of the direct Kalman filtering method based on UKF and the filtering method can improve the navigation accuracy.

    Inertial integrated navigation; Unscented Kalman filter; Direct filtering

    2015 - 02 - 05;

    2015-03-12。

    國(guó)家自然科學(xué)基金(60874112)

    王林(1982 - ),男,碩士研究生,主要從事有關(guān)信息傳輸與處理技術(shù)方向的研究。

    E-mail:wl19820416@sohu.com

    TP273

    A

    2095-8110(2015)03-0043-06

    猜你喜歡
    導(dǎo)航系統(tǒng)卡爾曼濾波慣性
    你真的了解慣性嗎
    沖破『慣性』 看慣性
    說(shuō)說(shuō)“北斗導(dǎo)航系統(tǒng)”
    “北斗”導(dǎo)航系統(tǒng)是怎樣煉成的
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    無(wú)處不在的慣性
    一種GNSS/SINS容錯(cuò)深組合導(dǎo)航系統(tǒng)設(shè)計(jì)
    解讀全球第四大導(dǎo)航系統(tǒng)
    普遍存在的慣性
    基于模糊卡爾曼濾波算法的動(dòng)力電池SOC估計(jì)
    97在线视频观看| 蜜桃久久精品国产亚洲av| 秋霞在线观看毛片| 久久精品国产亚洲网站| 日本黄色片子视频| 亚洲欧美精品综合久久99| 极品教师在线视频| 成年av动漫网址| 精品久久久久久久久久久久久| 特大巨黑吊av在线直播| 欧美另类亚洲清纯唯美| 国产av一区在线观看免费| 看片在线看免费视频| 精品少妇黑人巨大在线播放 | 麻豆乱淫一区二区| kizo精华| 国内少妇人妻偷人精品xxx网站| 免费av毛片视频| 毛片一级片免费看久久久久| 欧美另类亚洲清纯唯美| 99视频精品全部免费 在线| 国产一区二区在线观看日韩| 中文字幕熟女人妻在线| 免费看av在线观看网站| kizo精华| 亚洲精品成人久久久久久| 美女xxoo啪啪120秒动态图| 国产麻豆成人av免费视频| 欧美一级a爱片免费观看看| 中文亚洲av片在线观看爽| 精品国产三级普通话版| 国产黄片视频在线免费观看| 亚洲不卡免费看| 国产毛片a区久久久久| 国产三级在线视频| 少妇裸体淫交视频免费看高清| 国产成人aa在线观看| 国内揄拍国产精品人妻在线| 美女国产视频在线观看| 中文乱码字字幕精品一区二区三区 | 日本欧美国产在线视频| 能在线免费观看的黄片| 好男人视频免费观看在线| 搡老妇女老女人老熟妇| 青春草视频在线免费观看| av视频在线观看入口| 久久国内精品自在自线图片| 欧美xxxx性猛交bbbb| 午夜精品国产一区二区电影 | 久久精品人妻少妇| 中文天堂在线官网| 3wmmmm亚洲av在线观看| 男女边吃奶边做爰视频| 亚洲精品影视一区二区三区av| 少妇人妻精品综合一区二区| 欧美+日韩+精品| 久久99蜜桃精品久久| 成人漫画全彩无遮挡| 麻豆精品久久久久久蜜桃| 亚洲精华国产精华液的使用体验| 国产乱人视频| 国产黄片视频在线免费观看| 尤物成人国产欧美一区二区三区| 国内精品美女久久久久久| 亚洲精品色激情综合| 亚洲欧美清纯卡通| 国产在线男女| av国产免费在线观看| 岛国在线免费视频观看| 精品人妻一区二区三区麻豆| 国产一级毛片七仙女欲春2| 又粗又硬又长又爽又黄的视频| 丰满少妇做爰视频| 日韩 亚洲 欧美在线| 国产高清国产精品国产三级 | 久久久精品94久久精品| 91久久精品国产一区二区成人| 精品欧美国产一区二区三| 少妇人妻精品综合一区二区| 国产在线一区二区三区精 | 亚洲成人中文字幕在线播放| 久久久久久伊人网av| 欧美一区二区亚洲| 2021少妇久久久久久久久久久| 日本三级黄在线观看| 欧美不卡视频在线免费观看| 人妻夜夜爽99麻豆av| 亚洲欧洲国产日韩| 日本wwww免费看| 啦啦啦韩国在线观看视频| 最近最新中文字幕大全电影3| 真实男女啪啪啪动态图| 国产午夜精品一二区理论片| 波野结衣二区三区在线| 成人亚洲欧美一区二区av| 日韩人妻高清精品专区| 最新中文字幕久久久久| 亚洲欧洲日产国产| 欧美zozozo另类| 免费在线观看成人毛片| 国产乱来视频区| 久久99热这里只频精品6学生 | 国产亚洲av片在线观看秒播厂 | 蜜桃久久精品国产亚洲av| 女人久久www免费人成看片 | 老女人水多毛片| 在线观看66精品国产| 人人妻人人澡欧美一区二区| ponron亚洲| 丝袜喷水一区| 熟妇人妻久久中文字幕3abv| 中文亚洲av片在线观看爽| 精华霜和精华液先用哪个| 亚洲av成人av| 欧美性感艳星| 午夜福利视频1000在线观看| 国产视频内射| 国产精品野战在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲真实伦在线观看| 亚州av有码| av在线播放精品| 三级毛片av免费| 亚洲av成人精品一区久久| 夜夜看夜夜爽夜夜摸| 国产男人的电影天堂91| 亚洲一区高清亚洲精品| 九九在线视频观看精品| 国产精品精品国产色婷婷| 老司机福利观看| 人体艺术视频欧美日本| 我要看日韩黄色一级片| 亚洲图色成人| 汤姆久久久久久久影院中文字幕 | 久久久久国产网址| 中文字幕制服av| 日韩视频在线欧美| 免费无遮挡裸体视频| av又黄又爽大尺度在线免费看 | 男女视频在线观看网站免费| 晚上一个人看的免费电影| 国产免费男女视频| 婷婷色麻豆天堂久久 | 国产精品不卡视频一区二区| 成人毛片60女人毛片免费| av播播在线观看一区| 亚洲欧洲国产日韩| 国国产精品蜜臀av免费| 亚洲av电影在线观看一区二区三区 | 伦理电影大哥的女人| 亚洲精品日韩av片在线观看| 狠狠狠狠99中文字幕| 超碰97精品在线观看| 国产精品野战在线观看| 精品一区二区三区人妻视频| 国产成人免费观看mmmm| 久久精品久久精品一区二区三区| 久久久国产成人精品二区| 国产伦理片在线播放av一区| 成人毛片60女人毛片免费| 国产一区二区三区av在线| 国产欧美另类精品又又久久亚洲欧美| 在线观看美女被高潮喷水网站| 51国产日韩欧美| 级片在线观看| 九草在线视频观看| 91aial.com中文字幕在线观看| 男女国产视频网站| 在现免费观看毛片| 精品久久久久久久久久久久久| 99久久九九国产精品国产免费| 黄片无遮挡物在线观看| 久久久精品欧美日韩精品| 国产成人aa在线观看| 联通29元200g的流量卡| 欧美三级亚洲精品| 亚洲综合色惰| 日韩欧美精品免费久久| 亚洲精品456在线播放app| 三级男女做爰猛烈吃奶摸视频| 久久人人爽人人爽人人片va| 亚洲不卡免费看| 久久久久久九九精品二区国产| 天美传媒精品一区二区| 亚洲丝袜综合中文字幕| 久久久久精品久久久久真实原创| 国产成人aa在线观看| 成人高潮视频无遮挡免费网站| 七月丁香在线播放| 嘟嘟电影网在线观看| 久久精品91蜜桃| 国产成年人精品一区二区| av卡一久久| 久热久热在线精品观看| 国产精品永久免费网站| 久久鲁丝午夜福利片| 久久精品人妻少妇| 中文精品一卡2卡3卡4更新| 91aial.com中文字幕在线观看| 九草在线视频观看| 永久免费av网站大全| www.色视频.com| 亚洲精品乱久久久久久| 国产精品三级大全| 欧美性感艳星| 免费在线观看成人毛片| 一级毛片aaaaaa免费看小| 亚洲av二区三区四区| 日韩精品有码人妻一区| 久久99热这里只频精品6学生 | 国产乱人偷精品视频| 99视频精品全部免费 在线| 久久精品久久久久久噜噜老黄 | 一个人看视频在线观看www免费| 亚洲激情五月婷婷啪啪| 亚洲18禁久久av| 禁无遮挡网站| 亚洲天堂国产精品一区在线| 99视频精品全部免费 在线| 国产伦精品一区二区三区视频9| 一个人观看的视频www高清免费观看| 欧美极品一区二区三区四区| 搞女人的毛片| 久久久久精品久久久久真实原创| 亚洲婷婷狠狠爱综合网| 亚洲精品乱码久久久v下载方式| 亚洲,欧美,日韩| 永久免费av网站大全| 亚洲18禁久久av| 日本欧美国产在线视频| 亚洲久久久久久中文字幕| 欧美激情久久久久久爽电影| 日韩欧美精品v在线| 免费大片18禁| 久久人人爽人人片av| 亚洲人成网站在线观看播放| 91久久精品电影网| 九色成人免费人妻av| 亚洲国产精品成人久久小说| 成人一区二区视频在线观看| ponron亚洲| 免费看光身美女| 边亲边吃奶的免费视频| 高清av免费在线| 一区二区三区乱码不卡18| 精品国内亚洲2022精品成人| 神马国产精品三级电影在线观看| 亚洲国产欧美人成| 日韩一区二区三区影片| 又爽又黄a免费视频| 一级黄色大片毛片| 久久99蜜桃精品久久| 亚洲国产欧美在线一区| 99久久无色码亚洲精品果冻| 免费在线观看成人毛片| 国产激情偷乱视频一区二区| 九草在线视频观看| 久久国产乱子免费精品| 69人妻影院| 99热全是精品| 亚洲国产日韩欧美精品在线观看| 只有这里有精品99| 国产精品女同一区二区软件| 国产精品一二三区在线看| 国产伦理片在线播放av一区| 97人妻精品一区二区三区麻豆| 听说在线观看完整版免费高清| 最近2019中文字幕mv第一页| 久久久久久久久久成人| 久久精品国产自在天天线| 性插视频无遮挡在线免费观看| 国产精品永久免费网站| av.在线天堂| 97超视频在线观看视频| 日本欧美国产在线视频| 最近最新中文字幕大全电影3| 最近中文字幕2019免费版| 日本wwww免费看| 国模一区二区三区四区视频| 国产色爽女视频免费观看| 边亲边吃奶的免费视频| 久久久久久大精品| 国产一级毛片在线| 黄色日韩在线| 我的女老师完整版在线观看| 听说在线观看完整版免费高清| 免费在线观看成人毛片| 人人妻人人看人人澡| 亚洲av中文av极速乱| 日韩强制内射视频| 久久精品国产亚洲网站| 大又大粗又爽又黄少妇毛片口| 2022亚洲国产成人精品| 一个人免费在线观看电影| 91精品国产九色| 国内精品一区二区在线观看| 美女国产视频在线观看| 午夜精品在线福利| 小说图片视频综合网站| 97热精品久久久久久| 亚洲性久久影院| 国产精品人妻久久久影院| 久久久久久久久久久丰满| 成人三级黄色视频| 亚洲人与动物交配视频| 精品一区二区三区人妻视频| 一边摸一边抽搐一进一小说| 黄色日韩在线| 校园人妻丝袜中文字幕| 国产午夜精品论理片| 男女边吃奶边做爰视频| 成人美女网站在线观看视频| 日韩制服骚丝袜av| 亚洲精品乱码久久久v下载方式| av免费观看日本| 少妇熟女欧美另类| ponron亚洲| 在现免费观看毛片| 日本黄大片高清| 精品久久久久久久末码| www.av在线官网国产| 国产精品三级大全| 日韩一区二区视频免费看| 97超视频在线观看视频| 又黄又爽又刺激的免费视频.| 国产不卡一卡二| 成人鲁丝片一二三区免费| 熟女电影av网| 国产一区二区亚洲精品在线观看| 国产成人福利小说| 亚洲一级一片aⅴ在线观看| 国内精品宾馆在线| 国产精品不卡视频一区二区| 黄片wwwwww| 日日摸夜夜添夜夜爱| 中文字幕制服av| 国产精品久久久久久久电影| 亚洲av福利一区| 麻豆乱淫一区二区| av福利片在线观看| 国产一区二区亚洲精品在线观看| 日韩视频在线欧美| 亚洲国产欧美在线一区| 国产免费又黄又爽又色| 国产男人的电影天堂91| 青春草视频在线免费观看| 国产成人免费观看mmmm| 又爽又黄a免费视频| 久久亚洲精品不卡| 中文字幕亚洲精品专区| av国产免费在线观看| 国产亚洲一区二区精品| 国产三级在线视频| 99久国产av精品| 国产午夜精品久久久久久一区二区三区| 1024手机看黄色片| 99久久九九国产精品国产免费| 国产男人的电影天堂91| 国产淫片久久久久久久久| 国产探花在线观看一区二区| 亚洲中文字幕日韩| 又粗又硬又长又爽又黄的视频| 国产精品伦人一区二区| 免费av不卡在线播放| 国产私拍福利视频在线观看| 午夜精品在线福利| 一级毛片aaaaaa免费看小| 国产乱人视频| 女人十人毛片免费观看3o分钟| 一卡2卡三卡四卡精品乱码亚洲| 全区人妻精品视频| 最近视频中文字幕2019在线8| av线在线观看网站| 人妻夜夜爽99麻豆av| 久久6这里有精品| 成人午夜精彩视频在线观看| av又黄又爽大尺度在线免费看 | 老司机影院毛片| 免费搜索国产男女视频| 少妇熟女aⅴ在线视频| 中文字幕熟女人妻在线| 高清视频免费观看一区二区 | 国产精品久久久久久久久免| 少妇被粗大猛烈的视频| 18禁动态无遮挡网站| 精品久久久久久久久亚洲| 91精品国产九色| 欧美bdsm另类| 九九爱精品视频在线观看| 少妇被粗大猛烈的视频| 国产精品久久电影中文字幕| 日日撸夜夜添| 丝袜喷水一区| 国产高清视频在线观看网站| 精品久久久久久久人妻蜜臀av| 久久久久久久午夜电影| 99热精品在线国产| 日本猛色少妇xxxxx猛交久久| 成人午夜精彩视频在线观看| 国产精品爽爽va在线观看网站| 精品久久久久久久人妻蜜臀av| 麻豆国产97在线/欧美| 99热精品在线国产| 国产精品精品国产色婷婷| 亚洲欧美成人精品一区二区| 一级二级三级毛片免费看| 高清视频免费观看一区二区 | 欧美极品一区二区三区四区| 亚洲欧洲国产日韩| 99热这里只有精品一区| 草草在线视频免费看| 亚洲美女搞黄在线观看| 国产精品一及| 国产 一区 欧美 日韩| av免费在线看不卡| 大又大粗又爽又黄少妇毛片口| 久久久久久国产a免费观看| 亚洲人成网站在线观看播放| 国产精品精品国产色婷婷| 一级黄色大片毛片| 丰满少妇做爰视频| 永久网站在线| 大香蕉97超碰在线| 波多野结衣高清无吗| 韩国高清视频一区二区三区| 国产一区二区三区av在线| 亚洲精品456在线播放app| 国内揄拍国产精品人妻在线| 精品久久久久久久久久久久久| 亚洲一级一片aⅴ在线观看| 国产免费又黄又爽又色| 中国国产av一级| 国产午夜精品论理片| 久久综合国产亚洲精品| 国内精品一区二区在线观看| 国产在视频线精品| 国产极品精品免费视频能看的| 日本欧美国产在线视频| 亚洲激情五月婷婷啪啪| 免费看av在线观看网站| 亚洲精品日韩av片在线观看| 欧美极品一区二区三区四区| 免费观看a级毛片全部| 人妻少妇偷人精品九色| 免费av不卡在线播放| 99久久中文字幕三级久久日本| 欧美xxxx黑人xx丫x性爽| kizo精华| 国产精品一及| 欧美另类亚洲清纯唯美| 日产精品乱码卡一卡2卡三| 卡戴珊不雅视频在线播放| kizo精华| 国产精品,欧美在线| 18禁动态无遮挡网站| 国产三级中文精品| 日韩欧美精品免费久久| 少妇熟女aⅴ在线视频| 18禁在线播放成人免费| 久久这里有精品视频免费| 国产精品一及| 亚洲美女搞黄在线观看| 美女大奶头视频| 国产精品一二三区在线看| 麻豆一二三区av精品| 日日啪夜夜撸| 久久久色成人| 三级男女做爰猛烈吃奶摸视频| 欧美性猛交╳xxx乱大交人| 国产免费又黄又爽又色| 日韩制服骚丝袜av| 欧美丝袜亚洲另类| 高清午夜精品一区二区三区| 少妇裸体淫交视频免费看高清| 啦啦啦观看免费观看视频高清| 久久久成人免费电影| 91精品国产九色| 汤姆久久久久久久影院中文字幕 | 精品久久国产蜜桃| 成人美女网站在线观看视频| 男女边吃奶边做爰视频| 亚洲无线观看免费| 免费搜索国产男女视频| 91久久精品电影网| 亚洲一区高清亚洲精品| 婷婷六月久久综合丁香| 国产淫语在线视频| 嫩草影院精品99| 国产成人freesex在线| 天天躁日日操中文字幕| 日韩欧美 国产精品| 麻豆av噜噜一区二区三区| 91av网一区二区| 日本三级黄在线观看| 国产男人的电影天堂91| 内射极品少妇av片p| 一二三四中文在线观看免费高清| 日本免费在线观看一区| 一级黄片播放器| av免费观看日本| 尤物成人国产欧美一区二区三区| 久久久久久国产a免费观看| 久久久久国产网址| 一本一本综合久久| av在线观看视频网站免费| 久久精品国产亚洲av天美| 亚洲自偷自拍三级| 尾随美女入室| 18禁在线无遮挡免费观看视频| 国产成人精品久久久久久| 久久久久性生活片| 日韩中字成人| 久久久久久久亚洲中文字幕| 视频中文字幕在线观看| av天堂中文字幕网| 久久精品夜色国产| 色综合亚洲欧美另类图片| 男女那种视频在线观看| 久久久久久久午夜电影| 两性午夜刺激爽爽歪歪视频在线观看| 六月丁香七月| 亚洲第一区二区三区不卡| 成人av在线播放网站| 97人妻精品一区二区三区麻豆| 精品99又大又爽又粗少妇毛片| 亚洲怡红院男人天堂| 亚洲av成人精品一区久久| 亚洲内射少妇av| 国产又色又爽无遮挡免| 一级毛片久久久久久久久女| 一级黄色大片毛片| 成人av在线播放网站| 少妇人妻一区二区三区视频| 高清在线视频一区二区三区 | 三级国产精品片| 3wmmmm亚洲av在线观看| 高清在线视频一区二区三区 | 日韩 亚洲 欧美在线| 中文欧美无线码| 国产av一区在线观看免费| 国产成年人精品一区二区| a级一级毛片免费在线观看| 99九九线精品视频在线观看视频| 国产精品一区二区三区四区免费观看| 好男人视频免费观看在线| 欧美精品一区二区大全| 亚洲国产成人一精品久久久| 久久久成人免费电影| 中文在线观看免费www的网站| 精品久久久久久久久av| 欧美日韩国产亚洲二区| 亚洲av熟女| 九九在线视频观看精品| 最近2019中文字幕mv第一页| 亚洲成人中文字幕在线播放| 一级爰片在线观看| a级毛片免费高清观看在线播放| 久久久久久久久久成人| 美女xxoo啪啪120秒动态图| 美女被艹到高潮喷水动态| 高清日韩中文字幕在线| 干丝袜人妻中文字幕| 亚洲精品亚洲一区二区| 毛片女人毛片| 高清午夜精品一区二区三区| 一夜夜www| 免费观看精品视频网站| 精品午夜福利在线看| 如何舔出高潮| 国产一区二区三区av在线| 国产老妇伦熟女老妇高清| 最近中文字幕2019免费版| 18禁在线播放成人免费| 人人妻人人看人人澡| 国产成人精品婷婷| 亚洲国产精品久久男人天堂| 欧美日韩在线观看h| 桃色一区二区三区在线观看| 91午夜精品亚洲一区二区三区| 日本欧美国产在线视频| 精品一区二区三区人妻视频| 亚洲va在线va天堂va国产| 小说图片视频综合网站| 嫩草影院入口| 不卡视频在线观看欧美| 日本av手机在线免费观看| 国产亚洲精品av在线| 国产高清国产精品国产三级 | 国产乱来视频区| 免费观看a级毛片全部| 女的被弄到高潮叫床怎么办| 日产精品乱码卡一卡2卡三| 又爽又黄无遮挡网站| 99久久中文字幕三级久久日本| 亚洲成av人片在线播放无| 最后的刺客免费高清国语| 1000部很黄的大片| 成年免费大片在线观看| 亚洲精品乱久久久久久| 欧美激情国产日韩精品一区| 中文乱码字字幕精品一区二区三区 | 最新中文字幕久久久久| 国产伦精品一区二区三区四那| 欧美xxxx黑人xx丫x性爽| 亚洲欧美精品专区久久| 菩萨蛮人人尽说江南好唐韦庄 | 日韩三级伦理在线观看| 国产探花极品一区二区| 国产淫片久久久久久久久| 国产成人91sexporn| 中文资源天堂在线| 国产精品嫩草影院av在线观看| 国产精品av视频在线免费观看| 国产亚洲精品av在线| 久久久久九九精品影院|