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

    基于χ2檢驗(yàn)的慣性/衛(wèi)星緊耦合系統(tǒng)故障檢測(cè)與隔離方法*

    2016-03-24 08:09:33趙子陽李醒飛張家川
    火力與指揮控制 2016年2期
    關(guān)鍵詞:偽距慣導(dǎo)慣性

    楊 濤,趙子陽,李醒飛,張家川

    (1.天津大學(xué)精密儀器與光電子工程學(xué)院,天津 300072;2.天津航海儀器研究所,天津 300131)

    ?

    基于χ2檢驗(yàn)的慣性/衛(wèi)星緊耦合系統(tǒng)故障檢測(cè)與隔離方法*

    楊濤1,2,趙子陽2,李醒飛1,張家川2

    (1.天津大學(xué)精密儀器與光電子工程學(xué)院,天津300072;2.天津航海儀器研究所,天津300131)

    摘要:研究了一種基于χ2檢驗(yàn)的慣性/衛(wèi)星緊耦合系統(tǒng)故障檢測(cè)與隔離方法,該方法可利用故障檢測(cè)函數(shù)自動(dòng)識(shí)別單個(gè)衛(wèi)導(dǎo)故障,并實(shí)時(shí)剔除故障星重構(gòu)觀測(cè)信息序列,從而避免了對(duì)緊耦合系統(tǒng)的影響。構(gòu)建數(shù)字化仿真環(huán)境對(duì)所研究的故障檢測(cè)與隔離方法進(jìn)行了仿真,仿真結(jié)果表明,該方法可在衛(wèi)星發(fā)生故障時(shí)及時(shí)檢測(cè)并隔離故障星,并自動(dòng)實(shí)現(xiàn)觀測(cè)信息重構(gòu),保證了緊耦合系統(tǒng)的精度和完好性。

    關(guān)鍵詞:殘差χ2檢驗(yàn),慣性/衛(wèi)星,緊耦合,故障檢測(cè)與隔離

    0 引言

    慣性/衛(wèi)星緊耦合系統(tǒng)直接利用衛(wèi)星偽距、偽距率等原始觀測(cè)量信息估計(jì)慣導(dǎo)誤差和接收機(jī)鐘差,大大提高了組合系統(tǒng)的精度和連續(xù)性[1 -4]。然而,在工程實(shí)際中發(fā)現(xiàn),在未進(jìn)行完好性監(jiān)測(cè)時(shí),緊耦合濾波器的穩(wěn)定性易受錯(cuò)誤的星歷參數(shù)、接收機(jī)原始觀測(cè)量跳數(shù)等因素的影響。衛(wèi)星導(dǎo)航系統(tǒng)曾多次出現(xiàn)過星歷參數(shù)故障,其中最近一次發(fā)生在2014年4月2日,由于地面控制站注入了錯(cuò)誤的星歷參數(shù),GLONASS在軌衛(wèi)星星歷參數(shù)誤差驟然增大,導(dǎo)致定位誤差一度達(dá)數(shù)十千米[5]。此外,當(dāng)衛(wèi)導(dǎo)接收機(jī)發(fā)生故障時(shí),也會(huì)導(dǎo)致其輸出的偽距和偽距率等測(cè)量參數(shù)出現(xiàn)野值,這些都會(huì)嚴(yán)重干擾緊耦合系統(tǒng)的正常運(yùn)行。為此,開發(fā)一種可實(shí)時(shí)運(yùn)行的系統(tǒng)級(jí)故障檢測(cè)與隔離方法顯得十分必要。

    國(guó)內(nèi)外已有許多學(xué)者開展了接收機(jī)自主完好性監(jiān)測(cè)(RAIM)方法及慣導(dǎo)輔助的完好性監(jiān)測(cè)方法研究[6 -9],傳統(tǒng)的RAIM方法對(duì)星座分布等外部條件的要求較嚴(yán)格,使用受到一定限制,且絕大部分慣導(dǎo)輔助的衛(wèi)導(dǎo)完好性監(jiān)測(cè)研究成果都側(cè)重于理論建模與仿真,未考慮故障檢測(cè)與隔離方法的實(shí)時(shí)運(yùn)行。本文擬首先建立慣性/衛(wèi)星緊耦合系統(tǒng)的數(shù)學(xué)模型,之后研究?jī)煞N基于χ2檢驗(yàn)的故障檢測(cè)與隔離方法,利用模擬的慣性和衛(wèi)星系統(tǒng)開展數(shù)字仿真,對(duì)所研究的故障檢測(cè)與隔離方法進(jìn)行驗(yàn)證,最后得出結(jié)論。

    1 慣性/衛(wèi)星緊耦合系統(tǒng)數(shù)學(xué)模型

    1.1系統(tǒng)狀態(tài)方程

    對(duì)于艦船等載體,可選取慣導(dǎo)及衛(wèi)導(dǎo)接收機(jī)誤差參數(shù)等作為誤差狀態(tài)[10]:

    XI=[△Ve,△Vn,△準(zhǔn),△λ,△φx,△φy,△φz,△εx,△εy,εφz,△ax,△ay,△tu,△tur]

    其中:△Ve,△Vn分別為慣導(dǎo)東向速度誤差、北向速度誤差;△準(zhǔn),△λ分別為慣導(dǎo)緯度誤差、經(jīng)度誤差;△φx,△φy,△φz分別為慣導(dǎo)在載體坐標(biāo)系x、y、z方向的姿態(tài)角誤差;△εx,△εy,εφz分別為慣導(dǎo)載體坐標(biāo)系x、y、z方向陀螺漂移誤差;△ax,△ay分別為慣導(dǎo)載體坐標(biāo)系x、y方向加速度計(jì)零位誤差;△tu,△tur分別為衛(wèi)導(dǎo)接收機(jī)鐘差及鐘差變化率。

    系統(tǒng)狀態(tài)方程為,

    X=FX+GW(1)

    式中,X為狀態(tài)變量,F(xiàn)為狀態(tài)轉(zhuǎn)移矩陣,G為噪聲驅(qū)動(dòng)陣,W為系統(tǒng)噪聲。

    1.2系統(tǒng)量測(cè)方程

    采用偽距量測(cè)時(shí),慣性/衛(wèi)星緊耦合系統(tǒng)的量測(cè)信息為接收機(jī)輸出的多星座測(cè)量偽距與預(yù)測(cè)偽距之間的差值,預(yù)測(cè)偽距一般利用慣導(dǎo)位置信息及衛(wèi)星星歷計(jì)算得到。系統(tǒng)量測(cè)方程為:

    Z=HX+V(2)

    式中,Z為系統(tǒng)觀測(cè)量,H為量測(cè)矩陣,V為量測(cè)噪聲。

    2 基于χ2檢驗(yàn)的緊耦合系統(tǒng)故障檢測(cè)與隔離方法

    2.1總體思路

    慣性/衛(wèi)星緊耦合系統(tǒng)故障檢測(cè)與隔離方法的總體思路是:利用慣導(dǎo)設(shè)備提供的短時(shí)間精密導(dǎo)航信息,為衛(wèi)導(dǎo)系統(tǒng)提供序貫新息,利用慣導(dǎo)輸出的用戶位置及衛(wèi)星星歷計(jì)算預(yù)測(cè)偽距,通過各顆可見星測(cè)量偽距與預(yù)測(cè)偽距的一致性校驗(yàn)來判斷GNSS故障。即使在可見星數(shù)少于5顆的條件下,也可利用該方法對(duì)緊耦合系統(tǒng)進(jìn)行完好性監(jiān)測(cè),并對(duì)故障星進(jìn)行檢測(cè)與隔離。

    在緊組合卡爾曼濾波器中,完好性監(jiān)測(cè)在濾波器的時(shí)間更新和量測(cè)更新之間進(jìn)行,步驟如下:

    ①對(duì)濾波器的被估計(jì)狀態(tài)(包括慣導(dǎo)誤差和衛(wèi)導(dǎo)鐘差等狀態(tài))進(jìn)行時(shí)間更新,計(jì)算狀態(tài)估計(jì)預(yù)測(cè)值和均方誤差估計(jì)預(yù)測(cè)值;

    ②利用預(yù)測(cè)的誤差狀態(tài)信息計(jì)算量測(cè)預(yù)報(bào)值及預(yù)報(bào)偽距差,進(jìn)而計(jì)算殘差rk和殘差方差A(yù)k,并進(jìn)行衛(wèi)導(dǎo)故障檢測(cè);

    ③當(dāng)故障檢測(cè)函數(shù)超過告警門限時(shí),利用m-1法進(jìn)行m次故障檢測(cè)排除故障星,最后利用無故障的可見星信息重組觀測(cè)量Zk進(jìn)行量測(cè)更新,完成一個(gè)濾波循環(huán)。

    整個(gè)故障檢測(cè)與隔離的流程圖如圖1所示。

    圖1慣性/衛(wèi)星緊耦合系統(tǒng)故障檢測(cè)與隔離流程圖

    2.2殘差χ2檢測(cè)與隔離方法

    卡爾曼濾波器中,殘差可記為:

    在慣導(dǎo)輔助衛(wèi)導(dǎo)接收機(jī)完好性監(jiān)測(cè)時(shí),Hkk/k-1=為一步狀態(tài)估計(jì)量測(cè)預(yù)報(bào)值,是由k時(shí)刻慣導(dǎo)誤差傳播方程狀態(tài)更新推出的誤差狀態(tài)估計(jì)值乘以量測(cè)陣Hk得到的量測(cè)預(yù)報(bào)值,即偽距差的預(yù)報(bào)值。由于卡爾曼濾波屬于線性最小方差估計(jì),因此,殘差rk是白噪聲序列,其方差為

    定義檢測(cè)函數(shù)為:

    顯然λk符合χ2分布。由于贊k/ k-1由{Zk-1,Zk-2,…,Z1}計(jì)算得到,可以得到rk各元素不相關(guān),因此,λk的自由度等于rk的維數(shù)m(即緊耦合系統(tǒng)可見衛(wèi)星數(shù)),此時(shí)λk~χ2(m)。確定故障判定準(zhǔn)則為:

    式中,故障判定門限TD可由可見衛(wèi)星數(shù)m和虛警率確定,當(dāng)故障檢測(cè)周期為1 s,虛警率為0.002/h時(shí),判定門限與可見衛(wèi)星數(shù)關(guān)系如表1所示。

    表1檢測(cè)門限與可見星數(shù)的關(guān)系

    當(dāng)檢測(cè)函數(shù)超過門限時(shí),可以判斷存在故障,然后利用m-1組合法進(jìn)行m次殘差檢驗(yàn),即一次排除一顆可見星,利用余下的m-1顆可見星重組偽距差量測(cè)量,并進(jìn)行殘差檢測(cè),以識(shí)別故障星。

    2.3狀態(tài)χ2檢測(cè)與隔離方法

    由于卡爾曼濾波采用遞推方式實(shí)現(xiàn)狀態(tài)估計(jì),當(dāng)量測(cè)量Zk存在慢變故障時(shí)(如衛(wèi)導(dǎo)時(shí)鐘漂移故障),濾波器將有可能將定位結(jié)果的偏差視為用戶狀態(tài)的正常變化而修正Pk,導(dǎo)致濾波器不再正確跟蹤用戶真正狀態(tài),即出現(xiàn)“故障跟蹤”現(xiàn)象。此時(shí),在同樣故障偏差條件下,慢變故障對(duì)應(yīng)的檢測(cè)函數(shù)小于快變故障對(duì)應(yīng)的λk,使得殘差檢驗(yàn)法不易檢測(cè)出故障,采用狀態(tài)外推檢驗(yàn)法可以解決這一問題。

    無關(guān),僅僅是通過慣導(dǎo)誤差傳播方程外推出的偽距差量測(cè)預(yù)報(bào)值。因此,Z贊sk不會(huì)受s時(shí)刻到k時(shí)刻區(qū)間內(nèi)帶故障的量測(cè)更新影響,所以不會(huì)出現(xiàn)殘差χ2檢測(cè)與隔離方法的“故障跟蹤”現(xiàn)象。

    遞推周期m=k-s應(yīng)根據(jù)衛(wèi)導(dǎo)慢變故障檢測(cè)的設(shè)計(jì)要求進(jìn)行選擇。遞推周期太短,故障偏差導(dǎo)致的濾波器量測(cè)量殘差在短的時(shí)間內(nèi)無法發(fā)展到可以被檢測(cè)出的程度,因而易產(chǎn)生漏檢;遞推周期太長(zhǎng),一方面會(huì)導(dǎo)致狀態(tài)外推的精度降低,使得濾波狀態(tài)在遞推中會(huì)越來越偏離真實(shí)值,造成無故障時(shí)殘差超差,另一方面,遞推周期過長(zhǎng)也嚴(yán)重影響了故障檢測(cè)的反應(yīng)速度。重新定義殘差與檢測(cè)函數(shù)為:

    3 仿真驗(yàn)證與分析

    仿真系統(tǒng)主要包括用戶真實(shí)軌跡模擬、慣導(dǎo)模擬、衛(wèi)導(dǎo)模擬和故障模擬4個(gè)模塊,仿真系統(tǒng)的流程圖如圖2所示。

    圖2故障檢測(cè)與隔離算法仿真流程圖

    仿真參數(shù)的設(shè)置如下:仿真總時(shí)間3 600 s;初始位置為東經(jīng)117.143 23°,北緯39.182 33°,海拔0 m;衛(wèi)導(dǎo)的輸出頻率為1 Hz,遮蔽角為5°;慣導(dǎo)的輸出頻率為100 Hz;采用全部可見星濾波,濾波器周期為1 Hz;在第1 100 s~1 500 s模擬2號(hào)星突變故障,故障大小為100 m。

    仿真的主要步驟如下:

    ①進(jìn)行用戶真實(shí)軌跡模擬,并在真實(shí)軌跡中加入慣導(dǎo)器件誤差模擬慣導(dǎo)輸出,結(jié)合衛(wèi)導(dǎo)星座模擬模塊生成的衛(wèi)星位置模擬測(cè)量偽距和預(yù)測(cè)偽距。

    ②在選定的故障衛(wèi)星量測(cè)偽距上疊加一個(gè)故障值,然后將其送入慣性/衛(wèi)星緊耦合濾波器進(jìn)行濾波,同時(shí)在濾波器中加入上述故障檢測(cè)與隔離算法,以檢測(cè)和隔離故障星。

    ③利用濾波器誤差狀態(tài)的估計(jì)值來校正慣導(dǎo)輸出,將全部仿真時(shí)間內(nèi)校正后的慣導(dǎo)輸出和相應(yīng)的用戶真實(shí)位置的數(shù)據(jù)一一記錄并進(jìn)行比較,繪制誤差曲線圖;記錄每一次濾波的故障檢測(cè)函數(shù)值,與相應(yīng)時(shí)刻的檢測(cè)門限值做比較,并繪制故障檢測(cè)函數(shù)曲線圖。

    圖3未進(jìn)行故障檢測(cè)與隔離前的經(jīng)度誤差

    圖4未進(jìn)行故障檢測(cè)與隔離時(shí)的緯度誤差

    圖3和圖4給出了衛(wèi)導(dǎo)故障時(shí),未進(jìn)行故障檢測(cè)與隔離前的組合導(dǎo)航定位結(jié)果,圖5給出了對(duì)應(yīng)的檢測(cè)函數(shù)及門限曲線圖。

    圖5未進(jìn)行故障檢測(cè)與隔離前的檢測(cè)函數(shù)及門限值

    圖6進(jìn)行故障檢測(cè)與隔離后的經(jīng)度誤差

    圖6和圖7給出了衛(wèi)導(dǎo)故障時(shí),進(jìn)行故障檢測(cè)與隔離后的組合導(dǎo)航定位結(jié)果,圖8給出了對(duì)應(yīng)的檢測(cè)函數(shù)及門限曲線圖。

    圖7進(jìn)行故障檢測(cè)與隔離后的緯度誤差

    圖8進(jìn)行故障檢測(cè)與隔離后的檢測(cè)函數(shù)及門限值

    從仿真結(jié)果可看出,未進(jìn)行故障檢測(cè)與隔離前,當(dāng)可見星出現(xiàn)故障時(shí)會(huì)造成緊耦合系統(tǒng)定位精度嚴(yán)重下降,最大位置誤差可達(dá)70 m,此時(shí),故障函數(shù)遠(yuǎn)遠(yuǎn)超出門限值;加入故障檢測(cè)與隔離模塊后,當(dāng)可見星發(fā)生故障時(shí),殘差χ2檢測(cè)與隔離方法能及時(shí)檢測(cè)并隔離故障星,進(jìn)行系統(tǒng)重構(gòu)后,故障函數(shù)值遠(yuǎn)遠(yuǎn)低于門限值,系統(tǒng)定位精度也未受到影響。

    4 結(jié)論

    針對(duì)傳統(tǒng)基于衛(wèi)導(dǎo)自身信息的RAIM方法只有在可見衛(wèi)星數(shù)大于或等于5顆時(shí)可用的問題,提出一種基于χ2檢驗(yàn)的慣性/衛(wèi)星緊耦合系統(tǒng)故障檢測(cè)與隔離方法,該方法取消了對(duì)可見衛(wèi)星數(shù)的限制,當(dāng)導(dǎo)航衛(wèi)星發(fā)生故障時(shí)可實(shí)時(shí)檢測(cè)并隔離障星,并利用剩余無故障可見星重構(gòu)量測(cè)值序列,從而保證了組合導(dǎo)航系統(tǒng)的定位信息的精度完好性。

    參考文獻(xiàn):

    [1]XU X S,ZHANG G L. Implementation of tightly coupled GPS/INS navigation algorithm on DSP[C]//Computer Design and Applications,2010.

    [2]WANG X L,LI Y F. An innovative scheme for SINS/GPS ultra-tight integration system with low-grade IMU[J]. Aerospace Science and Technology,2012,23(1):452-460.

    [3]鄭辛,付夢(mèng)印.SINS/GPS緊耦合組合導(dǎo)航[J].中國(guó)慣性技術(shù)學(xué)報(bào),2011,19(1):33-37.

    [4]閆捷,徐曉蘇,張濤,等.艦載小型化SINS/GNSS緊組合導(dǎo)航系統(tǒng)設(shè)計(jì)[J].中國(guó)慣性技術(shù)學(xué)報(bào),2013,21(6):775-780.

    [5]Bhatti U I,Ochieng W Y. Performance of rate detector algorithms for an integrated GPS/INS systems in the presence of slowly growing error[J]. GPS Solutions,2012,16(3):293-301.

    [6]劉海穎,岳亞洲,楊毅鈞.基于多解分離的GNSS/Inertial組合系統(tǒng)完好性監(jiān)測(cè)[J].中國(guó)慣性技術(shù)學(xué)報(bào),2012,20 (1):63-68.

    [7]劉海穎,馮成濤,王惠南.一種慣性輔助衛(wèi)星導(dǎo)航系統(tǒng)及其完好性檢測(cè)方法[J].宇航學(xué)報(bào),2011,32(4):775-780.

    [8]鐘麗娜,劉建業(yè),李榮冰,等.基于LSSVM的慣性/衛(wèi)星緊組合導(dǎo)航系統(tǒng)緩變故障檢測(cè)方法[J].中國(guó)慣性技術(shù)學(xué)報(bào),2012,20(6):754-758.

    [9]陳坡,孫付平,景曉鵬,等.基于慣導(dǎo)輔助的GNSS完好性檢測(cè)方法研究[J].大地測(cè)量與地球動(dòng)力學(xué),2013,33(2):101-104.

    [10]王惠南. GPS導(dǎo)航原理與應(yīng)用[M].北京:科學(xué)出版社,2003.

    聲明

    近期有不法機(jī)構(gòu)冒用《火力與指揮控制》編輯部名義,以非法手段騙取作者的稿件及財(cái)物,因此,我刊鄭重聲明:

    從2005年至今,我刊一直使用以下兩種投稿方式:

    ①電子投稿。投稿郵箱為:HLYZ@chinajoumal.net.cn;hlyz207@126.com,在線投稿尚未啟用,啟用后予以告知。

    ②信函投稿。郵寄地址為:山西省太原市193信箱《火力與指揮控制》編輯部。郵編:030006。

    除此之外無其他任何方式向作者索稿。投稿時(shí)可電話咨詢:0351-8725026或0351-8725316。

    特此聲明!

    《火力與指揮控制》編輯部

    2016年2月25日

    Approach for Fault Detection and Isolation in INS/GNSS Tightly- coupled System Based on Chi- square Test

    YANG Tao1,2,ZHAO Zi-yang2,LI Xing-fei1,ZHANG Jia-chuan2
    (1. Tianjin University School of Precision Instrument and Opto-electronics Engineering,Tianjin 300072,China;2. Tianjin Navigation Instruments Research Institute,Tianjin 300131,China)

    Abstract:An approach for fault detection and isolation of INS/GPS tightly-coupled system was studied,which could recognize single navigation satellite error automatically by the utilization of fault detection function,could eliminate the fault satellite and reconstruct measurements array,and could avoid the influence of the fault on system. The digital simulation environment was constructed,then simulations of method for fault detection and isolation were conducted,results of which indicated that the proposed method could detect and isolate the fault satellite in time,and could reconstruct measurements information,which ensured the accuracy and integrity of tightly-coupled system.

    Key words:chi-square test of filter residual,INS/GNSS,tightly-coupled,fault detection and isolation

    作者簡(jiǎn)介:楊濤(1979-),男,河北南宮人,博士后,高級(jí)工程師。研究方向:慣性/衛(wèi)星組合導(dǎo)航技術(shù)。

    *基金項(xiàng)目:海軍重點(diǎn)預(yù)研項(xiàng)目(401010301);天津市企業(yè)博士后創(chuàng)新基金資助項(xiàng)目(2014-010)

    收稿日期:2014-12-20修回日期:2015-02-14

    文章編號(hào):1002-0640(2016)02-0169-04

    中圖分類號(hào):V249.32

    文獻(xiàn)標(biāo)識(shí)碼:A

    猜你喜歡
    偽距慣導(dǎo)慣性
    你真的了解慣性嗎
    沖破『慣性』 看慣性
    自適應(yīng)模糊多環(huán)控制在慣導(dǎo)平臺(tái)穩(wěn)定回路中的應(yīng)用
    無人機(jī)室內(nèi)視覺/慣導(dǎo)組合導(dǎo)航方法
    無處不在的慣性
    北斗偽距觀測(cè)值精度分析
    基于Bagging模型的慣導(dǎo)系統(tǒng)誤差抑制方法
    GNSS偽距粗差的開窗探測(cè)及修復(fù)
    普遍存在的慣性
    基于多線程的慣導(dǎo)邏輯仿真器設(shè)計(jì)
    每晚都被弄得嗷嗷叫到高潮| 国产成人欧美| bbb黄色大片| 国内精品久久久久久久电影| 国产av不卡久久| 成人免费观看视频高清| 日本免费一区二区三区高清不卡| 亚洲va日本ⅴa欧美va伊人久久| 岛国在线观看网站| 日本在线视频免费播放| 搞女人的毛片| 日韩成人在线观看一区二区三区| or卡值多少钱| 中文字幕精品亚洲无线码一区 | 老司机深夜福利视频在线观看| 久久精品人妻少妇| 伊人久久大香线蕉亚洲五| 日本黄色视频三级网站网址| 午夜免费成人在线视频| 女人高潮潮喷娇喘18禁视频| 久久久精品欧美日韩精品| 欧美日韩一级在线毛片| 亚洲 欧美一区二区三区| 亚洲成人免费电影在线观看| 他把我摸到了高潮在线观看| 中国美女看黄片| 亚洲九九香蕉| www国产在线视频色| 亚洲国产日韩欧美精品在线观看 | 香蕉国产在线看| √禁漫天堂资源中文www| 天天一区二区日本电影三级| 搡老岳熟女国产| 叶爱在线成人免费视频播放| 国产成人一区二区三区免费视频网站| 免费观看精品视频网站| 9191精品国产免费久久| 人人妻,人人澡人人爽秒播| 国产一区二区激情短视频| 最近在线观看免费完整版| 精品第一国产精品| 999精品在线视频| 日本精品一区二区三区蜜桃| 国产欧美日韩一区二区精品| 精品一区二区三区av网在线观看| 亚洲电影在线观看av| 丰满的人妻完整版| 99riav亚洲国产免费| 正在播放国产对白刺激| 一区二区三区激情视频| 国产成年人精品一区二区| 午夜福利视频1000在线观看| 黄色a级毛片大全视频| 一区二区三区激情视频| 午夜影院日韩av| 男女床上黄色一级片免费看| 人成视频在线观看免费观看| 中文字幕人妻熟女乱码| 亚洲av中文字字幕乱码综合 | www国产在线视频色| 超碰成人久久| 国产又色又爽无遮挡免费看| 亚洲精华国产精华精| 欧洲精品卡2卡3卡4卡5卡区| 一二三四在线观看免费中文在| 午夜福利一区二区在线看| 国产精品 国内视频| 午夜福利成人在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 色综合亚洲欧美另类图片| 午夜福利在线观看吧| √禁漫天堂资源中文www| 一二三四在线观看免费中文在| 免费av毛片视频| 欧美成人一区二区免费高清观看 | 色综合欧美亚洲国产小说| 欧美 亚洲 国产 日韩一| 精华霜和精华液先用哪个| 亚洲精品中文字幕一二三四区| 精品久久蜜臀av无| 777久久人妻少妇嫩草av网站| 久久亚洲精品不卡| a在线观看视频网站| 日本一本二区三区精品| 日韩欧美 国产精品| 琪琪午夜伦伦电影理论片6080| 久久午夜亚洲精品久久| 成人精品一区二区免费| 日韩成人在线观看一区二区三区| 母亲3免费完整高清在线观看| 精品欧美一区二区三区在线| 亚洲专区国产一区二区| 欧美激情 高清一区二区三区| 桃红色精品国产亚洲av| 丁香六月欧美| 91成人精品电影| 日本免费a在线| 免费看日本二区| 亚洲午夜精品一区,二区,三区| 男女做爰动态图高潮gif福利片| 日韩大码丰满熟妇| 日韩大尺度精品在线看网址| 男人舔奶头视频| 俄罗斯特黄特色一大片| 美女扒开内裤让男人捅视频| 熟女电影av网| 欧美激情极品国产一区二区三区| av片东京热男人的天堂| 在线十欧美十亚洲十日本专区| 日韩欧美一区视频在线观看| 高清毛片免费观看视频网站| 久久 成人 亚洲| 日本五十路高清| 日日干狠狠操夜夜爽| 国产久久久一区二区三区| 成人特级黄色片久久久久久久| 日本精品一区二区三区蜜桃| 精品一区二区三区四区五区乱码| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品98久久久久久宅男小说| 国产精品一区二区三区四区久久 | 日本黄色视频三级网站网址| 成人亚洲精品一区在线观看| 人妻久久中文字幕网| 变态另类成人亚洲欧美熟女| 久久精品影院6| 亚洲人成伊人成综合网2020| 曰老女人黄片| 国产91精品成人一区二区三区| 午夜精品在线福利| 搡老熟女国产l中国老女人| 亚洲第一av免费看| 青草久久国产| 91九色精品人成在线观看| 搡老岳熟女国产| 欧美成人一区二区免费高清观看 | 国产精品永久免费网站| 国产成人欧美| 搡老妇女老女人老熟妇| 级片在线观看| 嫩草影视91久久| 免费女性裸体啪啪无遮挡网站| 午夜久久久在线观看| 精品国产超薄肉色丝袜足j| 琪琪午夜伦伦电影理论片6080| 亚洲成国产人片在线观看| 变态另类成人亚洲欧美熟女| 黄色视频不卡| 老司机午夜福利在线观看视频| 变态另类成人亚洲欧美熟女| 色综合站精品国产| 老司机福利观看| 日本成人三级电影网站| 久久香蕉国产精品| 亚洲欧美精品综合久久99| 国产精品亚洲一级av第二区| 天天添夜夜摸| 正在播放国产对白刺激| 欧美又色又爽又黄视频| 一级片免费观看大全| 精品久久久久久久久久免费视频| av中文乱码字幕在线| 免费人成视频x8x8入口观看| 精品国产乱子伦一区二区三区| av天堂在线播放| 亚洲国产精品成人综合色| 久久热在线av| 2021天堂中文幕一二区在线观 | 亚洲成人久久性| 久久中文字幕一级| 欧美久久黑人一区二区| 2021天堂中文幕一二区在线观 | 日韩欧美 国产精品| 久久中文字幕一级| 日本一本二区三区精品| 亚洲自拍偷在线| www.www免费av| 热re99久久国产66热| 午夜a级毛片| 麻豆国产av国片精品| 中文字幕精品免费在线观看视频| 国产精品二区激情视频| 可以免费在线观看a视频的电影网站| 国产精品免费视频内射| 欧美性猛交黑人性爽| 人妻久久中文字幕网| 欧美亚洲日本最大视频资源| АⅤ资源中文在线天堂| www.www免费av| 久久国产精品人妻蜜桃| 国产99白浆流出| 少妇被粗大的猛进出69影院| 亚洲aⅴ乱码一区二区在线播放 | 欧美一区二区精品小视频在线| 性色av乱码一区二区三区2| 成人18禁高潮啪啪吃奶动态图| 好男人电影高清在线观看| 满18在线观看网站| 日日夜夜操网爽| 午夜福利欧美成人| 国产成人av教育| 国产乱人伦免费视频| 久久久水蜜桃国产精品网| 成年女人毛片免费观看观看9| 亚洲三区欧美一区| 免费观看人在逋| 亚洲人成网站高清观看| 美女大奶头视频| 国产精华一区二区三区| 免费在线观看视频国产中文字幕亚洲| 日韩三级视频一区二区三区| 亚洲精品国产一区二区精华液| 色哟哟哟哟哟哟| 搡老妇女老女人老熟妇| 啦啦啦观看免费观看视频高清| 无遮挡黄片免费观看| 午夜免费激情av| 欧美日本视频| 亚洲国产欧洲综合997久久, | 99精品欧美一区二区三区四区| 国产亚洲欧美98| 精品国产亚洲在线| 大型av网站在线播放| 窝窝影院91人妻| av欧美777| 黄色片一级片一级黄色片| 一级a爱视频在线免费观看| 变态另类丝袜制服| 国产精品 欧美亚洲| 好男人电影高清在线观看| 中亚洲国语对白在线视频| 99国产精品99久久久久| 亚洲av美国av| 国产久久久一区二区三区| 欧美又色又爽又黄视频| 精品欧美国产一区二区三| 熟妇人妻久久中文字幕3abv| 琪琪午夜伦伦电影理论片6080| 国产精品永久免费网站| 美女大奶头视频| 亚洲色图 男人天堂 中文字幕| 免费在线观看日本一区| 91九色精品人成在线观看| 午夜日韩欧美国产| 此物有八面人人有两片| 久久久久久国产a免费观看| 久久这里只有精品19| 亚洲自偷自拍图片 自拍| 国产一区在线观看成人免费| 久久精品国产亚洲av高清一级| 天堂影院成人在线观看| 黑人巨大精品欧美一区二区mp4| 国产成年人精品一区二区| 桃色一区二区三区在线观看| 777久久人妻少妇嫩草av网站| 在线免费观看的www视频| 国产亚洲av嫩草精品影院| 麻豆一二三区av精品| 国产亚洲精品一区二区www| 国产真人三级小视频在线观看| 黄片大片在线免费观看| 中文字幕久久专区| 妹子高潮喷水视频| 免费女性裸体啪啪无遮挡网站| 高清在线国产一区| 国产99白浆流出| 国产av不卡久久| 99热6这里只有精品| netflix在线观看网站| 在线观看一区二区三区| 成人手机av| 国产精品一区二区免费欧美| 好男人电影高清在线观看| 久久久久久久精品吃奶| xxx96com| 一区二区三区激情视频| 成人免费观看视频高清| 十分钟在线观看高清视频www| 精品少妇一区二区三区视频日本电影| 久久中文看片网| 无人区码免费观看不卡| 激情在线观看视频在线高清| 亚洲成av片中文字幕在线观看| 久久这里只有精品19| 白带黄色成豆腐渣| 97人妻精品一区二区三区麻豆 | 又大又爽又粗| 精品第一国产精品| 国产成人欧美| 欧美性猛交黑人性爽| xxxwww97欧美| 亚洲一区二区三区不卡视频| 丁香欧美五月| 人成视频在线观看免费观看| 男人操女人黄网站| 日日爽夜夜爽网站| 国产野战对白在线观看| 制服丝袜大香蕉在线| 亚洲免费av在线视频| 成年版毛片免费区| 久久精品国产清高在天天线| 两性夫妻黄色片| 久久精品国产综合久久久| 美女大奶头视频| 欧美成人一区二区免费高清观看 | 久久香蕉国产精品| 亚洲第一欧美日韩一区二区三区| 欧美不卡视频在线免费观看 | 欧美性猛交╳xxx乱大交人| 午夜福利成人在线免费观看| 国产一级毛片七仙女欲春2 | 色综合婷婷激情| 亚洲片人在线观看| 午夜福利视频1000在线观看| 亚洲av中文字字幕乱码综合 | 桃红色精品国产亚洲av| 可以在线观看毛片的网站| 久久国产精品男人的天堂亚洲| www.www免费av| 日日摸夜夜添夜夜添小说| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品久久视频播放| 国产精品电影一区二区三区| 在线视频色国产色| 亚洲精品久久国产高清桃花| 看免费av毛片| 天天躁狠狠躁夜夜躁狠狠躁| 精品欧美国产一区二区三| 国产久久久一区二区三区| 日本一本二区三区精品| 久久国产乱子伦精品免费另类| www.www免费av| 一区二区三区国产精品乱码| 久久精品aⅴ一区二区三区四区| 99久久精品国产亚洲精品| 国产精品亚洲一级av第二区| 成在线人永久免费视频| 精品久久久久久久久久免费视频| 免费看a级黄色片| 免费看十八禁软件| 午夜日韩欧美国产| 1024香蕉在线观看| 亚洲欧美激情综合另类| 亚洲电影在线观看av| 两个人视频免费观看高清| 色综合婷婷激情| АⅤ资源中文在线天堂| 国产久久久一区二区三区| 欧美日本亚洲视频在线播放| 午夜精品在线福利| 国产精品国产高清国产av| 国产午夜精品久久久久久| 久久久久国产一级毛片高清牌| 亚洲欧美精品综合一区二区三区| 亚洲性夜色夜夜综合| 在线播放国产精品三级| 村上凉子中文字幕在线| 久久九九热精品免费| 夜夜躁狠狠躁天天躁| avwww免费| 欧美丝袜亚洲另类 | 亚洲欧美日韩无卡精品| 99国产精品一区二区蜜桃av| 91麻豆av在线| 黄网站色视频无遮挡免费观看| 成人免费观看视频高清| 高清在线国产一区| 久久人人精品亚洲av| 国产精品久久久久久人妻精品电影| 欧美一级a爱片免费观看看 | 欧美成人午夜精品| 听说在线观看完整版免费高清| 日本一本二区三区精品| 日韩 欧美 亚洲 中文字幕| 欧美乱妇无乱码| 久久久久免费精品人妻一区二区 | 欧美乱妇无乱码| www日本在线高清视频| 哪里可以看免费的av片| 成人三级做爰电影| 日韩精品青青久久久久久| 国产成人影院久久av| 又黄又粗又硬又大视频| 国产av不卡久久| 久久伊人香网站| 精品人妻1区二区| 国产亚洲av嫩草精品影院| 在线十欧美十亚洲十日本专区| 波多野结衣高清作品| 欧美一级毛片孕妇| 精品电影一区二区在线| 久久精品91无色码中文字幕| 一本综合久久免费| 国产成人精品无人区| 手机成人av网站| 亚洲七黄色美女视频| 又黄又爽又免费观看的视频| 中亚洲国语对白在线视频| 亚洲五月色婷婷综合| 男人的好看免费观看在线视频 | 亚洲久久久国产精品| 看免费av毛片| 黑丝袜美女国产一区| 波多野结衣巨乳人妻| 成在线人永久免费视频| 岛国在线观看网站| 日韩高清综合在线| 色哟哟哟哟哟哟| 欧美日韩一级在线毛片| 欧美日本亚洲视频在线播放| 欧美另类亚洲清纯唯美| 免费电影在线观看免费观看| 亚洲欧美日韩高清在线视频| 视频区欧美日本亚洲| 亚洲国产日韩欧美精品在线观看 | 精品国产亚洲在线| 18禁黄网站禁片午夜丰满| 午夜福利在线在线| 人妻久久中文字幕网| 久久婷婷人人爽人人干人人爱| 中文字幕最新亚洲高清| 99在线人妻在线中文字幕| 中文字幕高清在线视频| 久久中文字幕一级| 国产黄片美女视频| 久久午夜亚洲精品久久| 神马国产精品三级电影在线观看 | 色婷婷久久久亚洲欧美| www.www免费av| 丝袜在线中文字幕| 在线看三级毛片| 国产又爽黄色视频| 91大片在线观看| 99在线人妻在线中文字幕| 午夜日韩欧美国产| 欧美成人性av电影在线观看| 精品少妇一区二区三区视频日本电影| 国产熟女午夜一区二区三区| 亚洲精品在线观看二区| 美女扒开内裤让男人捅视频| 岛国视频午夜一区免费看| 成人欧美大片| 少妇裸体淫交视频免费看高清 | 午夜视频精品福利| xxxwww97欧美| 欧美黑人精品巨大| 中文字幕久久专区| 久久中文字幕人妻熟女| 丝袜在线中文字幕| 在线观看日韩欧美| 日韩成人在线观看一区二区三区| 色尼玛亚洲综合影院| 天堂动漫精品| 久久久久免费精品人妻一区二区 | 精品国内亚洲2022精品成人| videosex国产| 女性被躁到高潮视频| 精品久久久久久久末码| 欧美亚洲日本最大视频资源| 亚洲免费av在线视频| 日本 av在线| 欧美久久黑人一区二区| 国产一区二区三区在线臀色熟女| 欧美激情久久久久久爽电影| 成人手机av| 波多野结衣高清无吗| 久久精品国产99精品国产亚洲性色| 久久精品aⅴ一区二区三区四区| 人人澡人人妻人| 成人免费观看视频高清| 亚洲一区中文字幕在线| 后天国语完整版免费观看| 国产av在哪里看| 午夜久久久在线观看| www.999成人在线观看| 中文字幕高清在线视频| 色播亚洲综合网| 国产一区二区三区视频了| 精品一区二区三区四区五区乱码| 久久精品夜夜夜夜夜久久蜜豆 | 中文资源天堂在线| 一级片免费观看大全| 黄色 视频免费看| 97超级碰碰碰精品色视频在线观看| 巨乳人妻的诱惑在线观看| 亚洲天堂国产精品一区在线| videosex国产| 国产黄片美女视频| 亚洲成人久久爱视频| 久久国产精品男人的天堂亚洲| 久久热在线av| 极品教师在线免费播放| 在线观看免费视频日本深夜| 精品一区二区三区四区五区乱码| 日韩欧美一区视频在线观看| 黄色毛片三级朝国网站| 日韩av在线大香蕉| 999精品在线视频| 99久久精品国产亚洲精品| 日本 av在线| 日日夜夜操网爽| 嫩草影院精品99| 欧美日本视频| 精品久久久久久久毛片微露脸| 亚洲一区中文字幕在线| 欧美乱妇无乱码| 国产一区二区三区视频了| 亚洲色图 男人天堂 中文字幕| 99国产极品粉嫩在线观看| 日韩欧美在线二视频| 亚洲精品粉嫩美女一区| 亚洲最大成人中文| 欧美色视频一区免费| 亚洲国产精品sss在线观看| 观看免费一级毛片| 久久人妻福利社区极品人妻图片| 一级黄色大片毛片| 老司机午夜十八禁免费视频| 欧美黑人精品巨大| 一区二区三区国产精品乱码| 色精品久久人妻99蜜桃| 国产精品久久久久久亚洲av鲁大| 成人18禁在线播放| 91成年电影在线观看| 一二三四在线观看免费中文在| 黄色片一级片一级黄色片| 日本 av在线| 久久国产亚洲av麻豆专区| 亚洲黑人精品在线| 国产精品综合久久久久久久免费| 色av中文字幕| 女生性感内裤真人,穿戴方法视频| 两性午夜刺激爽爽歪歪视频在线观看 | a在线观看视频网站| 人人妻,人人澡人人爽秒播| 精品久久久久久久久久久久久 | 国产成人av激情在线播放| 黄频高清免费视频| 色综合亚洲欧美另类图片| 国产又黄又爽又无遮挡在线| 午夜a级毛片| 制服丝袜大香蕉在线| 国产亚洲精品综合一区在线观看 | 日韩三级视频一区二区三区| 久久精品国产综合久久久| av视频在线观看入口| 精品少妇一区二区三区视频日本电影| 1024手机看黄色片| 韩国精品一区二区三区| 女人被狂操c到高潮| 国产又色又爽无遮挡免费看| 午夜成年电影在线免费观看| 久久精品国产清高在天天线| 高潮久久久久久久久久久不卡| 国产精品九九99| 亚洲国产欧美网| 欧美激情高清一区二区三区| 亚洲色图av天堂| 一边摸一边抽搐一进一小说| 久久婷婷成人综合色麻豆| 亚洲成人久久爱视频| 18禁观看日本| 亚洲欧美激情综合另类| 美女免费视频网站| 国内久久婷婷六月综合欲色啪| 岛国视频午夜一区免费看| 丁香六月欧美| 久久人妻福利社区极品人妻图片| 国语自产精品视频在线第100页| 黄片大片在线免费观看| 男人舔女人的私密视频| 免费在线观看亚洲国产| 亚洲真实伦在线观看| 亚洲男人的天堂狠狠| 国产精品免费一区二区三区在线| 热re99久久国产66热| 亚洲av熟女| 欧美亚洲日本最大视频资源| 日韩一卡2卡3卡4卡2021年| 日韩精品免费视频一区二区三区| 午夜久久久在线观看| 国产在线观看jvid| 国产又黄又爽又无遮挡在线| 国产精品爽爽va在线观看网站 | 久久人人精品亚洲av| 国产久久久一区二区三区| 岛国视频午夜一区免费看| 99久久99久久久精品蜜桃| 母亲3免费完整高清在线观看| 亚洲熟女毛片儿| 国产国语露脸激情在线看| www日本黄色视频网| xxx96com| 国产精品日韩av在线免费观看| 超碰成人久久| 一区二区三区精品91| 大型黄色视频在线免费观看| 久久草成人影院| 日韩欧美一区二区三区在线观看| 亚洲欧美日韩高清在线视频| 国产一区二区三区在线臀色熟女| 国语自产精品视频在线第100页| 男人的好看免费观看在线视频 | 亚洲av成人不卡在线观看播放网| 性欧美人与动物交配| 一区二区日韩欧美中文字幕| 午夜免费鲁丝| 久久中文字幕一级| 又黄又粗又硬又大视频| 母亲3免费完整高清在线观看| 色哟哟哟哟哟哟| 国内精品久久久久久久电影| 亚洲第一av免费看| 久久人妻av系列| 淫秽高清视频在线观看|