• <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ì)
    久久精品国产鲁丝片午夜精品| 亚洲av.av天堂| 亚洲电影在线观看av| 久久韩国三级中文字幕| 美女内射精品一级片tv| 国国产精品蜜臀av免费| 国产精品女同一区二区软件| 日韩强制内射视频| 久久97久久精品| 一级二级三级毛片免费看| 啦啦啦中文免费视频观看日本| 国产成人免费无遮挡视频| 黄色怎么调成土黄色| 免费观看无遮挡的男女| 性色avwww在线观看| 大香蕉97超碰在线| 亚洲精品第二区| 禁无遮挡网站| 国产亚洲一区二区精品| 最近最新中文字幕大全电影3| 国产久久久一区二区三区| 99热这里只有是精品50| 亚洲精品一二三| 伊人久久精品亚洲午夜| 欧美三级亚洲精品| 久久99精品国语久久久| 美女国产视频在线观看| 嫩草影院精品99| 久久99热这里只频精品6学生| 国精品久久久久久国模美| 欧美成人a在线观看| 久久99热这里只频精品6学生| 亚洲精品色激情综合| 国产精品秋霞免费鲁丝片| 成人毛片60女人毛片免费| 精品国产露脸久久av麻豆| 在线a可以看的网站| 国产午夜福利久久久久久| 国产午夜福利久久久久久| 亚洲精品久久久久久婷婷小说| 婷婷色麻豆天堂久久| 免费不卡的大黄色大毛片视频在线观看| 国产黄a三级三级三级人| 少妇人妻久久综合中文| 国产精品女同一区二区软件| 亚洲欧美一区二区三区国产| tube8黄色片| 人妻 亚洲 视频| 99久久九九国产精品国产免费| 国产成人精品久久久久久| 18+在线观看网站| 亚洲色图av天堂| 我的老师免费观看完整版| 中文欧美无线码| 爱豆传媒免费全集在线观看| 五月玫瑰六月丁香| 国产成人a∨麻豆精品| 午夜老司机福利剧场| 亚洲精品第二区| 国产 一区 欧美 日韩| 日韩在线高清观看一区二区三区| 永久网站在线| 最新中文字幕久久久久| 久久久国产一区二区| 亚洲一级一片aⅴ在线观看| 偷拍熟女少妇极品色| 久久精品久久精品一区二区三区| www.色视频.com| 精品久久久久久久久av| 国产成人免费无遮挡视频| 亚洲国产成人一精品久久久| 可以在线观看毛片的网站| 欧美 日韩 精品 国产| videos熟女内射| 看黄色毛片网站| 国产欧美亚洲国产| 午夜福利高清视频| 视频区图区小说| 激情五月婷婷亚洲| 国产淫片久久久久久久久| 在线免费十八禁| 美女高潮的动态| av黄色大香蕉| 欧美+日韩+精品| 国产av码专区亚洲av| 国内少妇人妻偷人精品xxx网站| 亚洲国产av新网站| 97人妻精品一区二区三区麻豆| 亚洲国产精品专区欧美| 亚洲精品久久午夜乱码| 在线免费观看不下载黄p国产| 欧美日韩亚洲高清精品| 久久久久久久久久成人| 国产精品成人在线| 91久久精品电影网| 色吧在线观看| av卡一久久| 人妻系列 视频| 夜夜看夜夜爽夜夜摸| 国产精品无大码| 2018国产大陆天天弄谢| 婷婷色综合www| 欧美极品一区二区三区四区| 久久久精品免费免费高清| 精品酒店卫生间| 亚洲第一区二区三区不卡| 成人亚洲精品一区在线观看 | 久久精品人妻少妇| 久久久久久久久大av| 少妇裸体淫交视频免费看高清| 亚洲最大成人手机在线| 国产 一区 欧美 日韩| 人妻制服诱惑在线中文字幕| 亚洲欧美一区二区三区黑人 | 亚洲国产精品专区欧美| 亚洲国产精品999| 性插视频无遮挡在线免费观看| 色吧在线观看| 神马国产精品三级电影在线观看| 日本-黄色视频高清免费观看| 九九在线视频观看精品| 日韩电影二区| 观看免费一级毛片| 91在线精品国自产拍蜜月| 男女边吃奶边做爰视频| 一本色道久久久久久精品综合| 亚洲av男天堂| av免费在线看不卡| 最近最新中文字幕免费大全7| 男人舔奶头视频| 国产高潮美女av| 成人毛片60女人毛片免费| 人妻夜夜爽99麻豆av| 久久99精品国语久久久| 亚洲成人精品中文字幕电影| 日本av手机在线免费观看| 午夜老司机福利剧场| 欧美成人精品欧美一级黄| 久久久欧美国产精品| 久久亚洲国产成人精品v| 久久热精品热| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲欧美一区二区三区国产| 国产高清三级在线| 一级爰片在线观看| 色综合色国产| 精品国产一区二区三区久久久樱花 | 国产v大片淫在线免费观看| 美女被艹到高潮喷水动态| 夫妻性生交免费视频一级片| 日本午夜av视频| 综合色av麻豆| 亚洲国产精品成人久久小说| 男人添女人高潮全过程视频| 91午夜精品亚洲一区二区三区| 国产爱豆传媒在线观看| 有码 亚洲区| 国产精品熟女久久久久浪| 久久99热6这里只有精品| av专区在线播放| 中文字幕久久专区| 亚洲最大成人中文| 亚洲精品日韩av片在线观看| 亚洲国产欧美在线一区| 69av精品久久久久久| 在线亚洲精品国产二区图片欧美 | 一区二区三区精品91| 少妇的逼好多水| 国产精品女同一区二区软件| 最近手机中文字幕大全| 国产精品伦人一区二区| 国产成人精品一,二区| 最近的中文字幕免费完整| 中文字幕久久专区| 91午夜精品亚洲一区二区三区| 亚洲人成网站在线播| 内射极品少妇av片p| 汤姆久久久久久久影院中文字幕| av在线播放精品| 少妇 在线观看| 成年女人在线观看亚洲视频 | 色哟哟·www| videossex国产| 精品亚洲乱码少妇综合久久| 成人免费观看视频高清| 又爽又黄a免费视频| 中国美白少妇内射xxxbb| 亚洲欧美日韩东京热| 国产精品久久久久久久电影| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 黄色配什么色好看| 亚洲成人中文字幕在线播放| 国产黄频视频在线观看| 亚洲国产精品成人综合色| 99精国产麻豆久久婷婷| 亚洲精品日韩av片在线观看| 天天躁日日操中文字幕| 真实男女啪啪啪动态图| 国产有黄有色有爽视频| 久久久精品94久久精品| 国产黄色视频一区二区在线观看| 久久精品国产a三级三级三级| .国产精品久久| 午夜日本视频在线| 午夜日本视频在线| 狂野欧美激情性xxxx在线观看| 一级黄片播放器| 午夜免费男女啪啪视频观看| 伊人久久精品亚洲午夜| 人人妻人人澡人人爽人人夜夜| 亚洲第一区二区三区不卡| 大码成人一级视频| 水蜜桃什么品种好| 久久精品国产亚洲av涩爱| 亚洲人与动物交配视频| 中文字幕久久专区| av国产精品久久久久影院| 亚洲内射少妇av| 丝袜脚勾引网站| 2022亚洲国产成人精品| 大码成人一级视频| 99久久九九国产精品国产免费| 最新中文字幕久久久久| 搡老乐熟女国产| 国产精品嫩草影院av在线观看| 狂野欧美激情性bbbbbb| 18禁在线无遮挡免费观看视频| 色婷婷久久久亚洲欧美| 青青草视频在线视频观看| 国内揄拍国产精品人妻在线| av免费在线看不卡| 免费看a级黄色片| 日韩av在线免费看完整版不卡| .国产精品久久| 久热这里只有精品99| 中文资源天堂在线| 国产高清不卡午夜福利| 成人毛片a级毛片在线播放| 国产成人午夜福利电影在线观看| 国产老妇女一区| 永久网站在线| 深夜a级毛片| 成人国产麻豆网| 国产精品精品国产色婷婷| 2021少妇久久久久久久久久久| 热re99久久精品国产66热6| 亚洲av成人精品一区久久| 99久久人妻综合| 熟女电影av网| 深爱激情五月婷婷| 免费黄网站久久成人精品| www.av在线官网国产| 久久久精品94久久精品| 午夜精品国产一区二区电影 | 午夜福利高清视频| 欧美一区二区亚洲| av黄色大香蕉| 国产免费一级a男人的天堂| 又大又黄又爽视频免费| 亚洲精品国产成人久久av| 久久女婷五月综合色啪小说 | 亚州av有码| 蜜臀久久99精品久久宅男| 熟女电影av网| 少妇被粗大猛烈的视频| 又大又黄又爽视频免费| 欧美另类一区| av国产免费在线观看| 一级毛片黄色毛片免费观看视频| 黄片wwwwww| 在线观看国产h片| 亚洲,一卡二卡三卡| 人妻少妇偷人精品九色| 亚洲久久久久久中文字幕| 免费大片黄手机在线观看| 亚洲av国产av综合av卡| 免费观看性生交大片5| av黄色大香蕉| 日本与韩国留学比较| 国产成人a∨麻豆精品| 国产综合懂色| 亚洲精品日韩在线中文字幕| 国产亚洲91精品色在线| 亚洲人与动物交配视频| 大香蕉久久网| 丰满人妻一区二区三区视频av| 亚洲成人久久爱视频| 亚洲四区av| 成人欧美大片| 亚洲欧洲国产日韩| av女优亚洲男人天堂| 99热6这里只有精品| 国产综合懂色| 国产午夜精品久久久久久一区二区三区| 午夜福利视频精品| 在线观看人妻少妇| 99久久精品一区二区三区| 久久99精品国语久久久| av免费观看日本| 久久久精品94久久精品| 青春草国产在线视频| 精品少妇黑人巨大在线播放| 久久久久久久久大av| 亚洲美女视频黄频| 欧美激情国产日韩精品一区| 少妇丰满av| 美女cb高潮喷水在线观看| 亚洲图色成人| 一区二区三区免费毛片| 久久精品国产a三级三级三级| 又爽又黄无遮挡网站| 亚洲av中文字字幕乱码综合| 白带黄色成豆腐渣| 街头女战士在线观看网站| 亚洲欧美成人综合另类久久久| 亚洲国产欧美在线一区| 大香蕉97超碰在线| av在线天堂中文字幕| av国产精品久久久久影院| 成人无遮挡网站| 欧美性猛交╳xxx乱大交人| 一个人观看的视频www高清免费观看| 韩国av在线不卡| 极品少妇高潮喷水抽搐| 啦啦啦中文免费视频观看日本| 久久久久九九精品影院| 99热这里只有是精品在线观看| 亚洲av中文字字幕乱码综合| 欧美精品人与动牲交sv欧美| 肉色欧美久久久久久久蜜桃 | 国产熟女欧美一区二区| 日韩欧美精品v在线| 激情五月婷婷亚洲| videossex国产| 人妻 亚洲 视频| 91久久精品电影网| 国产欧美亚洲国产| 午夜福利在线在线| 亚洲最大成人av| 成年免费大片在线观看| 日韩一区二区视频免费看| 久久久亚洲精品成人影院| 久久国内精品自在自线图片| 日韩一区二区视频免费看| 日韩不卡一区二区三区视频在线| 春色校园在线视频观看| 看免费成人av毛片| av福利片在线观看| 国产欧美另类精品又又久久亚洲欧美| 成人毛片a级毛片在线播放| 激情 狠狠 欧美| 亚洲久久久久久中文字幕| 国产爱豆传媒在线观看| 日日摸夜夜添夜夜爱| 日韩一区二区视频免费看| av在线app专区| 2022亚洲国产成人精品| 国语对白做爰xxxⅹ性视频网站| 极品少妇高潮喷水抽搐| 中文欧美无线码| 五月天丁香电影| 久久精品国产亚洲av涩爱| 自拍欧美九色日韩亚洲蝌蚪91 | 午夜激情福利司机影院| 国产一级毛片在线| 91午夜精品亚洲一区二区三区| 最近中文字幕高清免费大全6| 国产免费一区二区三区四区乱码| 久久亚洲国产成人精品v| 国产日韩欧美在线精品| 边亲边吃奶的免费视频| 校园人妻丝袜中文字幕| 午夜精品一区二区三区免费看| 国产精品国产三级国产av玫瑰| 乱系列少妇在线播放| 国产av码专区亚洲av| 日本爱情动作片www.在线观看| 黄色视频在线播放观看不卡| 国产精品一区二区三区四区免费观看| 国产成人91sexporn| 精品久久久久久久末码| 国产精品久久久久久精品古装| 大片免费播放器 马上看| 91狼人影院| 99久久中文字幕三级久久日本| 新久久久久国产一级毛片| 免费看a级黄色片| 丰满少妇做爰视频| 搞女人的毛片| 亚洲色图av天堂| 夜夜爽夜夜爽视频| 联通29元200g的流量卡| 蜜臀久久99精品久久宅男| 久久韩国三级中文字幕| 狠狠精品人妻久久久久久综合| 亚洲精品日韩av片在线观看| 我的老师免费观看完整版| 国产亚洲91精品色在线| 精品人妻一区二区三区麻豆| 插逼视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 午夜爱爱视频在线播放| 午夜视频国产福利| 99热网站在线观看| 国产一级毛片在线| 2022亚洲国产成人精品| 国产片特级美女逼逼视频| 免费看不卡的av| 我的女老师完整版在线观看| 日本黄色片子视频| 亚洲精品乱码久久久久久按摩| 亚洲欧美日韩另类电影网站 | 国产精品国产av在线观看| 久久99热这里只有精品18| 久久久久九九精品影院| 成人午夜精彩视频在线观看| 欧美激情国产日韩精品一区| 日本一本二区三区精品| 熟女电影av网| 99久久中文字幕三级久久日本| 国产免费视频播放在线视频| 欧美3d第一页| 国产探花极品一区二区| freevideosex欧美| 色哟哟·www| 亚洲精品国产av蜜桃| 大话2 男鬼变身卡| 我要看日韩黄色一级片| 成年人午夜在线观看视频| 波多野结衣巨乳人妻| 人妻制服诱惑在线中文字幕| 亚洲怡红院男人天堂| 免费av不卡在线播放| 国产亚洲午夜精品一区二区久久 | 两个人的视频大全免费| 我的女老师完整版在线观看| 亚洲精品乱久久久久久| av女优亚洲男人天堂| 亚洲欧美成人精品一区二区| 精品熟女少妇av免费看| 五月玫瑰六月丁香| 大香蕉97超碰在线| 国产欧美日韩精品一区二区| 国产免费一级a男人的天堂| 一级二级三级毛片免费看| 成人二区视频| 精品国产露脸久久av麻豆| 欧美另类一区| 久热久热在线精品观看| 国产成人精品久久久久久| 免费不卡的大黄色大毛片视频在线观看| 日韩av在线免费看完整版不卡| a级毛片免费高清观看在线播放| 成人鲁丝片一二三区免费| 色综合色国产| 99久久精品热视频| 成人鲁丝片一二三区免费| 成人二区视频| 人妻 亚洲 视频| 精品酒店卫生间| 欧美日韩视频高清一区二区三区二| 亚洲精品一区蜜桃| 又大又黄又爽视频免费| 色播亚洲综合网| 日本-黄色视频高清免费观看| 精品久久久久久久末码| 免费在线观看成人毛片| 欧美日韩综合久久久久久| av国产精品久久久久影院| av黄色大香蕉| 国产 一区精品| av国产久精品久网站免费入址| videos熟女内射| 男女下面进入的视频免费午夜| 麻豆精品久久久久久蜜桃| 看十八女毛片水多多多| 观看美女的网站| 国产午夜精品一二区理论片| 高清午夜精品一区二区三区| 久久久久九九精品影院| h日本视频在线播放| 色视频在线一区二区三区| 日日摸夜夜添夜夜添av毛片| 91久久精品国产一区二区三区| av黄色大香蕉| 一个人看的www免费观看视频| 欧美一级a爱片免费观看看| 永久网站在线| 一级毛片电影观看| 内地一区二区视频在线| 国产精品国产三级专区第一集| 精品人妻视频免费看| 婷婷色综合www| 亚洲天堂国产精品一区在线| av在线老鸭窝| 一区二区三区精品91| 能在线免费看毛片的网站| 在线播放无遮挡| 久久久欧美国产精品| 亚洲精品亚洲一区二区| 免费电影在线观看免费观看| 美女xxoo啪啪120秒动态图| 国产亚洲一区二区精品| 国产一级毛片在线| 人妻夜夜爽99麻豆av| 午夜激情久久久久久久| 免费av毛片视频| 久久97久久精品| 美女高潮的动态| 亚洲欧美精品自产自拍| 欧美少妇被猛烈插入视频| 成年人午夜在线观看视频| 日韩一区二区三区影片| 看免费成人av毛片| 中文精品一卡2卡3卡4更新| 精品人妻视频免费看| 亚洲最大成人手机在线| 亚洲,欧美,日韩| 蜜臀久久99精品久久宅男| 搡女人真爽免费视频火全软件| 亚洲精品乱久久久久久| 女人被狂操c到高潮| 天堂俺去俺来也www色官网| 日韩一区二区视频免费看| 国产精品一区二区在线观看99| 在线观看三级黄色| 青春草国产在线视频| 久久精品国产亚洲av涩爱| 中文精品一卡2卡3卡4更新| 国产成人福利小说| 亚洲美女搞黄在线观看| 人妻夜夜爽99麻豆av| 欧美成人一区二区免费高清观看| 免费av不卡在线播放| 黄片wwwwww| 高清av免费在线| 久久久久性生活片| 人妻夜夜爽99麻豆av| 成人国产av品久久久| 久久精品国产a三级三级三级| 免费播放大片免费观看视频在线观看| 免费看av在线观看网站| 成人黄色视频免费在线看| 免费黄色在线免费观看| 97在线人人人人妻| 性色av一级| 亚洲欧美中文字幕日韩二区| 中文在线观看免费www的网站| 夫妻性生交免费视频一级片| 亚洲av电影在线观看一区二区三区 | 久久久精品94久久精品| 国产亚洲av片在线观看秒播厂| 欧美人与善性xxx| 精品一区二区免费观看| 在线免费观看不下载黄p国产| 欧美最新免费一区二区三区| 国产伦精品一区二区三区四那| 最新中文字幕久久久久| 少妇人妻一区二区三区视频| 午夜免费鲁丝| 日韩亚洲欧美综合| 丝袜喷水一区| 亚洲最大成人中文| 成人国产av品久久久| 欧美xxxx黑人xx丫x性爽| 人体艺术视频欧美日本| 又粗又硬又长又爽又黄的视频| 欧美潮喷喷水| 全区人妻精品视频| 嘟嘟电影网在线观看| 日韩伦理黄色片| 亚洲精品中文字幕在线视频 | 一级二级三级毛片免费看| 国产乱人偷精品视频| 日韩制服骚丝袜av| 国产成人a区在线观看| 色吧在线观看| 99热这里只有是精品50| 网址你懂的国产日韩在线| 久久久精品免费免费高清| 伦精品一区二区三区| www.色视频.com| 亚洲国产精品成人久久小说| 日本黄大片高清| 久久人人爽人人爽人人片va| 亚洲色图av天堂| 永久免费av网站大全| 色婷婷久久久亚洲欧美| 国内精品美女久久久久久| 在线免费十八禁| 亚洲国产高清在线一区二区三| 春色校园在线视频观看| 国产精品.久久久| 能在线免费看毛片的网站| 久久鲁丝午夜福利片| av播播在线观看一区| 国产91av在线免费观看| 22中文网久久字幕| a级毛色黄片| 午夜福利高清视频| 国产精品久久久久久av不卡| 九九在线视频观看精品| 亚洲国产精品成人综合色| 中国美白少妇内射xxxbb| 国产高清有码在线观看视频| 内射极品少妇av片p| 精品一区二区三区视频在线| 国产精品国产三级专区第一集| 午夜精品一区二区三区免费看| 国产亚洲最大av| 久久久久久久久久人人人人人人| 亚洲国产色片| 边亲边吃奶的免费视频| 欧美成人精品欧美一级黄| 国国产精品蜜臀av免费|