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

    基于新息速率抗差估計的INS/GNSS組合導航系統(tǒng)欺騙檢測算法

    2021-10-17 23:52:02呂志偉張倫東高揚駿
    中國慣性技術學報 2021年3期
    關鍵詞:新息抗差檢測時間

    張 超,呂志偉,2,張倫東,高揚駿,2

    (1. 信息工程大學,鄭州 450000;2. 地理信息工程國家重點實驗室,西安 710054)

    慣性導航系統(tǒng)(Inertial Navigation System, INS)和全球衛(wèi)星導航系統(tǒng)(Global Navigation Satellite System, GNSS)的組合導航系統(tǒng)在軍用和民用領域得到了廣泛而深入的應用。但由于GNSS 信號微弱且民碼結構公開,使其容易受到欺騙干擾的影響。欺騙干擾是指干擾源播發(fā)虛假衛(wèi)星信號,誘導目標接收機進行捕獲與跟蹤,進而輸出錯誤信息的技術[1]。GNSS欺騙技術的分類方式有很多,按攻擊階段分類可分為捕獲階段欺騙干擾和跟蹤階段欺騙干擾,而跟蹤階段欺騙干擾按欺騙信號的碼相位和載波頻率與真實信號是否保持同步,又可分為同步式欺騙干擾和異步式欺騙干擾[2]。在欺騙干擾環(huán)境下,組合導航系統(tǒng)的GNSS 模塊鎖定欺騙信號并輸出錯誤的偽距和偽距率,進而影響測量更新階段的狀態(tài)誤差估計值并輸出錯誤的導航結果。因此,欺騙檢測對于提高組合導航系統(tǒng)的完好性和可用性顯得十分重要。

    INS/GNSS 組合導航系統(tǒng)的欺騙檢測算法可分為非耦合和耦合兩類[3]。非耦合欺騙檢測是指GNSS 與INS 獨立工作,將GNSS 輸出的位置、速度轉換為加速度和角速度,與INS 輸出的加速度和角速度進行一致性檢驗。耦合欺騙檢測是指GNSS 與INS 輸出的信息通過Kalman 濾波等方式融合,通常使用殘差或新息序列作為檢驗統(tǒng)計量,檢測方式采用假設檢驗。這類方法的缺點是在欺騙干擾影響下組合導航系統(tǒng)的反饋校正機制會使殘差或新息序列偏離正常值導致誤警率和漏警率升高[4]。耦合欺騙檢測又可分為“快照法”和“連續(xù)法”。“快照法”是以當前時刻的新息序列構成檢驗統(tǒng)計量,適合檢測GNSS 測量值中階躍類型的誤差,對緩變類型的誤差檢測時間較長,例如新息欺騙檢測、殘差卡方檢測屬于該類算法,在單路衛(wèi)星信號受到0.1 m/s 緩變類型的干擾時,新息欺騙檢測算法的檢測時間為360 s[5];“連續(xù)法”是將檢測窗口內的新息序列構成檢驗統(tǒng)計量,適合檢測緩變類型的誤差,但檢測窗口長度的設置需要在告警頻率和檢測成功率之間權衡,例如自主完好性檢測外推法屬于該類算法,在單顆衛(wèi)星受到0.1 m/s 緩變類型的干擾時,檢測窗口設置為150 s 的自主完好性檢測外推法的檢測時間為162 s[5]。

    針對GNSS 測量值中存在緩變類型的干擾,Bhatti提出了一種新息速率檢測算法[6],將新息的變化率作為檢測量,在單顆衛(wèi)星受到0.1 m/s 緩變類型的干擾時檢測時間為110 s,但該方法需要與自主完好性檢測外推法結合使用,無法針對單個GNSS 測量值檢測。針對組合導航系統(tǒng)反饋校正對耦合欺騙檢測的影響,Shizhuang Wang 將抗差估計與基于新息的檢測算法結合,單顆衛(wèi)星受到0.1 m/s 緩變類型的干擾時檢測時間縮短至65 s[7];Chuang Zhang 提出的檢測算法在單顆衛(wèi)星受到0.5 m/s 緩變類型的干擾時漏警率為26%[8],但這兩種算法需要設置檢測窗口;許睿等提出了MEDLL 輔助的GNSS/INS 系統(tǒng)欺騙信號檢測方法[9],成功檢測出牽引速率為2 m/s 的欺騙信號。近幾年有學者將人工智能(例如神經網絡[10]、支持向量機[4])、芯片級原子鐘[11]等應用到欺騙檢測算法中,但對計算能力要求較高,與現(xiàn)有系統(tǒng)兼容性不強。本文主要針對同步式欺騙干擾,其造成的影響可視為緩變類型的誤差[12]。通過將基于新息速率的檢測算法與抗差估計結合,在無需設置檢測窗口的情況下有效提高了檢測時間,降低了漏警率和誤警率。

    本文首先根據(jù)同步式欺騙原理建立了基于GNSS測量值的欺騙仿真模型,并介紹了INS/GNSS 緊耦合的新息欺騙檢測算法;然后針對新息欺騙檢測算法存在的缺陷,提出了一種基于新息速率抗差估計的INS/GNSS 組合導航系統(tǒng)欺騙檢測算法;最后通過仿真驗證了新算法的有效性。

    1 GNSS 欺騙仿真模型

    對于INS/GNSS 組合導航系統(tǒng)的欺騙仿真分析,基于GNSS 測量值的欺騙模型可以滿足仿真需求,避免了復雜的欺騙信號仿真,容易實現(xiàn)。本文模擬同步式欺騙信號。

    第i 顆衛(wèi)星在t 時刻的真實偽距測量值可表示為:

    其中τ(i)為信號傳播延遲,c 為光速, δtu和δt(i)為接收機鐘差和衛(wèi)星鐘差。

    異步式欺騙信號的碼相位和載波頻率與真實信號無法保持同步,所以其偽距測量值是任意的。同步式欺騙信號初期與真實信號保持同步,當接收機跟蹤環(huán)路鎖定欺騙信號后,為防止接收機自主完好性檢測,采取逐步誘導的方式使偽距測量值逐漸偏離真實值。

    綜上所述,基于GNSS 測量值的欺騙模型可表示為:

    tInit為欺騙信號鎖定接收機跟蹤環(huán)路的時刻,a為誘導速率,b 為真實偽距與欺騙偽距的偏差,且a≠ 0。當b= 0時, s ( t )代表同步式欺騙;當b≠0時, s ( t )代表異步式欺騙。

    2 新息欺騙檢測算法

    新息序列是實際觀測向量與一步預測狀態(tài)向量計算的觀測向量的差值,由系統(tǒng)模型誤差和觀測量誤差共同決定,當系統(tǒng)模型誤差較小時,新息序列可以反映觀測量誤差。

    設觀測矩陣為Hk,觀測向量為 Zk,一步預測狀態(tài)向量協(xié)方差矩陣為Pk|k-1,觀測噪聲協(xié)方差矩陣為Rk,則新息序列 rk與其協(xié)方差矩陣 Vk的表達式如下:

    歸一化后的新息可以表示為:

    其中rki為k 時刻新息序列的第i 個值(i=1…n,n為可見衛(wèi)星個數(shù)),Vkii為k 時刻rki的方差。 ωi代表k時刻歸一化后的第i 個新息值,可以反映第i 個GNSS測量值誤差。

    如果接收到的第i 路衛(wèi)星信號是真實信號,在無其他故障的情況下,第i 個測量值不存在粗差, ωi服從標準正態(tài)分布:

    如果接收到的第i 路衛(wèi)星信號是欺騙信號,則第i個測量值存在粗差, ωi服從非中心正態(tài)分布:

    根據(jù)實際應用對于導航系統(tǒng)完好性的需求,設定誤警率 Pfa,則對應的第i 個測量值的誤警率為[9]:

    構造假設H0:ωi~N(0,1)代表第i 路衛(wèi)星信號狀態(tài)正常,H1: ωi~N(δ ,1)代表第i 路衛(wèi)星信號狀態(tài)異常,根據(jù)標準正態(tài)分布概率密度函數(shù),由 α0可計算出門限值 Td,公式如下:

    根據(jù)假設檢驗原理,可由下式判斷第i 路衛(wèi)星信號是否存在異常:

    3 基于新息速率抗差估計的INS/GNSS 組合導航系統(tǒng)欺騙檢測算法

    3.1 基于新息速率的檢測算法

    新息速率欺騙檢測算法通過歸一化新息 ωi的變化速率來判斷GNSS 測量值是否存在異常??紤]到測量噪聲的影響,通常采用Kalman 濾波實時估計新息ωi的速率。

    Kalman 濾波估計的狀態(tài)向量定義為:

    系統(tǒng)模型定義為:

    觀測模型定義為:

    其中, ωi作為測量值輸入, ni為觀測噪聲。

    在無其他干擾情況下,新息速率 vi服從高斯分布,給定誤警率Pfd,由式(9)計算出單個測量值對應的誤警率 α0,則 vi對應的檢測門限為:

    其中,Q-1是高斯分布的逆,Pvi是狀態(tài)協(xié)方差陣中 vi的方差。

    同步欺騙導致GNSS 測量值誤差隨時間積累,如式(3)所示,a 越大,誤差積累越快。新息欺騙檢測算法存在最小可檢測誤差(Minimum Detectable Blunders, MDB)[8],只有當誤差積累到MDB 大小時才能檢測出同步欺騙。由于誤差累積會導致新息呈增長或減少趨勢,新息速率欺騙檢測算法通過檢測新息增長或減少的速率判定是否存在欺騙干擾,無需等誤差累積到一定程度才能被檢測。所以,新息速率欺騙檢測算法比新息欺騙檢測算法的檢測時間更短。但計算新息速率的濾波器在收斂之前無法進行檢測。

    3.2 基于新息速率抗差估計的檢測算法

    組合導航系統(tǒng)的反饋校正機制會使前一時刻GNSS 測量值誤差傳遞至下一時刻系統(tǒng)更新階段的狀態(tài)向量估計值中,進而使下一時刻原本狀態(tài)正常的新息偏離正常值,對基于新息序列的檢測算法造成干擾。欺騙干擾對新息序列的影響可表達為[4]:

    其中 Φk+1,k為狀態(tài)轉移矩陣,Kk為濾波增益矩陣,Δk=(Δ1…Δi…Δn)為 k 時刻欺騙干擾對GNSS 測量值的影響,Δr~k+1為k+1 時刻新息序列的偏移量。從式(16)可以看出,k 時刻欺騙干擾的影響會使k+1 時刻的新息序列發(fā)生偏移。

    油菜葉片在不同生育時期呈現(xiàn)不同的形態(tài),主要包括長柄葉、短柄葉和無柄葉3種類型。選用3種類型油菜葉片在室內和戶外環(huán)境下的圖像,如圖2(a)4幅圖像所示。利用本文方法所提取的結果(為便于觀察,將葉脈圖像作求補運算)如圖2(b)4幅圖像所示。

    通過引入抗差估計可以很好地削弱這種影響。選取IGG-3 等價權函數(shù)[13],使用新息速率計算等價權,具體公式如下:

    vi為第i 個GNSS 測量值對應的新息速率, wi為第i 個測量值的等價權。將 k1設定為新息速率的門限值 vD, k0設為0.5 倍的 k1。

    等價權矩陣可表示為:

    文獻[7]對增益矩陣Kk進行調整,具體公式如下:

    將KR代替式(16)中的Kk,可得:

    從式(20)可以看出,Δk中的元素Δi增大或減少,導致k 時刻新息 ωi增大或減少,進而新息速率vi會增大,等價權 wi降低。因此無論 Δi增大或減小其對應的等價權均會降低,從而降低k 時刻的欺騙干擾對k+1 時刻新息序列的影響。

    基于新息速率抗差估計的INS/GNSS 組合導航系統(tǒng)欺騙檢測算法流程如圖1 所示。組合導航濾波器輸出各路歸一化新息,計算對應的新息速率和等價權,等價權組成等價權矩陣用于計算KR,將KR替代原增益矩陣Kk輸入到組合導航濾波器。

    圖1 檢測算法流程圖Fig.1 Flow chart of detection algorithm

    3.3 多路衛(wèi)星信號受欺騙時的算法改進

    同步欺騙影響下的GNSS 測量值誤差隨時間累積的快慢與誘導速率有關。誘導速率越小,誤差累積的越慢,對應的新息速率 vi變化的就越慢, vi到達門限值 vD的時間就越久,對應的等價權 wi減至0 的時間變久,最終導致抗差估計的抑制效果降低。尤其多路衛(wèi)星信號受到低誘導速率的同步欺騙時,抗差估計的抑制效果會顯著降低。為此,改進等價權函數(shù)中參數(shù)k1,令 k1=κvD,0< κ≤1 。通過調節(jié)參數(shù)κ ,可以達到縮短等價權 wi減至0 的時間,從而提高抗差估計的抑制效果。

    4 仿真試驗

    為簡化表達,M1表示新息欺騙檢測算法,M2表示新息速率欺騙檢測算法,M3 表示基于新息速率抗差估計的INS/GNSS 組合導航系統(tǒng)欺騙檢測算法。

    針對INS/GNSS 緊耦合模式進行仿真實驗,設計了3 個算例。1)在1 顆衛(wèi)星受到同步欺騙時比較M1和M2 的檢測能力;2)在2 顆衛(wèi)星受到同步欺騙時比較M1、M2 和M3 的檢測能力;3)在3 顆衛(wèi)星受到同步欺騙時比較等價權函數(shù)中參數(shù) k1的取值對M3 檢測能力的影響。

    4.1 仿真條件

    設30 顆GNSS 衛(wèi)星均勻分布在6 個等圓軌道上。GNSS 模塊采用雙頻數(shù)據(jù)消除電離層影響。接收機誤差設為高斯白噪聲序列。誤警率設置為10-5,M1 的門限值為 4.95 m,M2 的門限值為 0.00276 m/s。GNSS 模塊參數(shù)見表1;IMU 模塊參數(shù)配置見表2;Kalman 濾波器配置參數(shù)見表3。

    表1 GNSS 模塊參數(shù)Tab.1 Parameters of GNSS model

    表2 IMU 模塊參數(shù)Tab.2 Parameters of IMU model

    表3 Kalman 濾波器參數(shù)Tab.3 Parameters of Kalman filter

    模擬飛機飛行軌跡,時長746 s,并據(jù)此產生GNSS 和IMU 測量值。飛行軌跡如圖2 所示。

    圖2 飛行軌跡Fig.2 Flight trajectory

    欺騙場景均采用同步式欺騙干擾,設置如表4 所示。

    表4 欺騙場景設置Tab.4 Spoofing scene setup

    4.2 仿真結果與分析

    1)實驗1

    通過實驗1 驗證M2 比M1 對同步欺騙的檢測時間短。M1 和M2 的仿真結果如圖3(a)和3(b)所示。

    圖3 M1 與M2 仿真結果比較Fig.3 Comparison of simulation results between M1 and M2

    由圖3(a)可知,M1 在300~500 s 內檢測失敗,將欺騙時間延長至 600 s,得到 M1 的檢測時間為210 s;由圖3(b)可知,M2 的檢測時間為80 s。所以,在1 路衛(wèi)星信號受到同步欺騙時,M2 比M1的檢測時間縮短了61.9%。

    2)實驗2

    通過實驗2 驗證抗差估計抑制新息序列偏離正常值的效果。M1、M2和M3的仿真結果如圖4(a)、4(b)和4(c)所示。由圖4(a)可知,M1 在300~500 s 內未檢測出1、2 號衛(wèi)星所受欺騙,若將欺騙時間延長至600 s,得到M1 的檢測時間為225 s;由圖4(b)可知,M2 的檢測時間為80 s,成功檢測出了1、2 號衛(wèi)星,但狀態(tài)正常的3、4 號衛(wèi)星的新息速率卻超過了門限;由圖4(c)可知,M3的檢測時間為76 s,其余衛(wèi)星的新息速率正常。

    圖4 M1、M2 和M3 仿真結果比較Fig.4 Comparison of simulation results among M1, M2 and M3

    為更好說明抗差估計的抑制效果,對實驗2 進行100 循環(huán)的蒙特卡洛仿真,得到三種算法的漏警率和誤警率如表5 所示。

    表5 實驗2 蒙特卡洛仿真結果Tab.5 Monte Carlo simulation results of Exp.2

    實驗2 表明:在2 路衛(wèi)星信號受到同步欺騙時,M2 相比M1 的檢測時間縮短了64.4%,漏警率降低至0%,但誤警率卻顯著上升。M3 不僅漏警率保持0%,誤警率也降低至0%,同時檢測時間在M2 的基礎上又縮短了5%。

    3)實驗3

    為驗證參數(shù) k1的取值對M3 檢測能力的影響,通過公式 k1=κvD調節(jié)參數(shù) k1的大小,κ =1 時的檢測結果如圖5(a)所示,κ =0.6 的檢測結果如圖5(b)所示。由圖5(a)可知,1 號、2 號和3 號衛(wèi)星的檢測時間分別為101 s、82 s 和108 s,但狀態(tài)正常的4 號衛(wèi)星卻超過了門限值;由圖5(b)可知,1 號、2 號和3 號衛(wèi)星的檢測時間分別為88 s、75 s 和88 s,其余狀態(tài)正常的衛(wèi)星新息速率均正常。

    圖5 κ =1 與κ =0.6 的仿真結果比較Fig.5 Comparison of simulation results between κ =1 and κ =0.6

    為更好說明參數(shù) k1的取值對M3 檢測能力的影響,對實驗3 進行100 循環(huán)的蒙特卡洛仿真,得到κ =1 和κ =0.6 的漏警率和誤警率如表6 所示。

    表6 實驗3 蒙特卡洛仿真結果Tab.6 Monte Carlo simulation results of Exp. 3

    實驗3 表明,κ =0.6 的檢測時間在κ =1 的基礎上進一步縮短了12%左右,同時誤警率降低了81%。但是漏警率有所升高。

    5 結 論

    本文針對基于新息的欺騙檢測算法對同步式欺騙干擾的檢測時間較長、漏警率和誤警率較高的問題,提出了一種基于抗差估計的INS/GNSS 導航系統(tǒng)欺騙檢測算法。該算法將新息的變化速率作為檢測量,并采用抗差估計的方法削弱了欺騙干擾對新息序列產生的影響,通過調參數(shù) k1,進一步提高了檢測能力。仿真實驗驗證了新算法的可行性,相比新息欺騙檢測算法,新算法對同步欺騙的檢測時間縮短了60%以上,漏警率和誤警率維持在4%以內。

    未來研究工作展望:1)新算法的檢測時間、誤警率和漏警率隨欺騙施加時長、誘導速率的變化規(guī)律有待進一步研究;2)新算法在不同精度等級的慣導系統(tǒng)下檢測能力的變化有待進一步研究。

    猜你喜歡
    新息抗差檢測時間
    傳遞函數(shù)辨識(21):線性回歸系統(tǒng)的遞階遞推參數(shù)估計
    對兩種細菌鑒定法在血液檢驗中的應用效果進行分析
    新型溶血素與傳統(tǒng)溶血素在臨床血常規(guī)檢驗中的應用研究
    M估計的強跟蹤SVD-UKF算法在組合導航中的應用
    電子科技(2018年7期)2018-07-23 05:30:32
    ABL90血氣分析儀在急診科的應用研究
    改善單頻PPP參數(shù)收斂速度的抗差估計方法
    不同檢測時長對粉煤灰砌塊放射性檢測結果的影響
    自適應卡爾曼濾波在航空重力異常解算的應用研究
    基于新息正交性自適應濾波的慣性/地磁組合導航方法
    地形簡化對DEM不確定性的抗差性研究
    精品国产亚洲在线| 美女高潮喷水抽搐中文字幕| 岛国在线免费视频观看| 精品日产1卡2卡| 精品国内亚洲2022精品成人| 亚洲成人免费电影在线观看| 亚洲美女视频黄频| 久久久久久久精品吃奶| 欧美黄色淫秽网站| 久久精品国产综合久久久| 两个人的视频大全免费| 国产国拍精品亚洲av在线观看 | 欧美高清成人免费视频www| 黄色成人免费大全| 亚洲国产精品999在线| 久久久国产精品麻豆| 亚洲最大成人手机在线| 久久久精品大字幕| 悠悠久久av| 国产精品,欧美在线| 91久久精品国产一区二区成人 | 99久久久亚洲精品蜜臀av| 欧美午夜高清在线| 亚洲精华国产精华精| 老司机福利观看| 好看av亚洲va欧美ⅴa在| 在线观看av片永久免费下载| 午夜福利免费观看在线| 可以在线观看的亚洲视频| 久久国产精品人妻蜜桃| 午夜两性在线视频| 欧美午夜高清在线| 亚洲欧美日韩高清专用| 黑人欧美特级aaaaaa片| 在线观看66精品国产| 欧美在线黄色| 夜夜爽天天搞| 观看免费一级毛片| 日韩中文字幕欧美一区二区| 国产91精品成人一区二区三区| 国产精品精品国产色婷婷| 九九热线精品视视频播放| 久久久色成人| 真人做人爱边吃奶动态| 91字幕亚洲| 麻豆国产97在线/欧美| 亚洲av免费高清在线观看| 真实男女啪啪啪动态图| 12—13女人毛片做爰片一| 一进一出好大好爽视频| 久久精品91无色码中文字幕| 男女那种视频在线观看| 91久久精品电影网| 亚洲精品粉嫩美女一区| 国产日本99.免费观看| 脱女人内裤的视频| 国产三级中文精品| 在线观看美女被高潮喷水网站 | 国内久久婷婷六月综合欲色啪| 手机成人av网站| 啦啦啦免费观看视频1| 超碰av人人做人人爽久久 | 久久久久久人人人人人| 亚洲精品美女久久久久99蜜臀| 一本久久中文字幕| 观看免费一级毛片| 中文字幕av成人在线电影| 19禁男女啪啪无遮挡网站| 精品一区二区三区视频在线 | 超碰av人人做人人爽久久 | 91字幕亚洲| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美激情综合另类| 熟妇人妻久久中文字幕3abv| av欧美777| 狠狠狠狠99中文字幕| 国产精品久久久久久久久免 | 国产精品久久久久久久电影 | 免费av观看视频| 久久久久亚洲av毛片大全| 露出奶头的视频| 午夜激情欧美在线| 亚洲国产欧美网| 精品人妻1区二区| 亚洲av电影在线进入| 在线免费观看不下载黄p国产 | 精品人妻一区二区三区麻豆 | 国产伦人伦偷精品视频| 亚洲国产中文字幕在线视频| 黄色成人免费大全| 麻豆国产97在线/欧美| 国产麻豆成人av免费视频| 午夜日韩欧美国产| 男女下面进入的视频免费午夜| 午夜亚洲福利在线播放| www.熟女人妻精品国产| 久久精品91蜜桃| 亚洲av美国av| 精品欧美国产一区二区三| 成人无遮挡网站| 亚洲人成网站在线播放欧美日韩| 女人十人毛片免费观看3o分钟| 悠悠久久av| 国产精品久久久久久久电影 | 亚洲精华国产精华精| 国产亚洲欧美98| 亚洲第一电影网av| 啦啦啦韩国在线观看视频| 少妇熟女aⅴ在线视频| 99视频精品全部免费 在线| 久久精品人妻少妇| 久久亚洲精品不卡| 久久久国产精品麻豆| 亚洲五月天丁香| 成人性生交大片免费视频hd| 精华霜和精华液先用哪个| 女人高潮潮喷娇喘18禁视频| 国产私拍福利视频在线观看| 国产精品女同一区二区软件 | 亚洲无线观看免费| 亚洲国产色片| 18禁美女被吸乳视频| 国产成人影院久久av| 美女 人体艺术 gogo| 国产久久久一区二区三区| 欧美乱码精品一区二区三区| av天堂中文字幕网| 中文字幕av成人在线电影| 国产伦精品一区二区三区视频9 | 国产日本99.免费观看| 国内精品久久久久精免费| 日韩精品中文字幕看吧| 精品电影一区二区在线| 久久这里只有精品中国| 一个人免费在线观看的高清视频| 欧美成人性av电影在线观看| 亚洲成人久久爱视频| 亚洲欧美精品综合久久99| 久久人妻av系列| 成人一区二区视频在线观看| 99精品在免费线老司机午夜| 亚洲狠狠婷婷综合久久图片| 久久久久久大精品| 亚洲精品久久国产高清桃花| 三级毛片av免费| 欧美色欧美亚洲另类二区| 免费看十八禁软件| 一区二区三区免费毛片| 国产一区二区激情短视频| 又粗又爽又猛毛片免费看| 悠悠久久av| 人妻夜夜爽99麻豆av| 国语自产精品视频在线第100页| 人人妻人人澡欧美一区二区| 国产野战对白在线观看| 欧美日本亚洲视频在线播放| 女同久久另类99精品国产91| 精品久久久久久久人妻蜜臀av| 亚洲无线观看免费| 1000部很黄的大片| 国产一区二区在线观看日韩 | 亚洲国产精品成人综合色| 中文字幕精品亚洲无线码一区| 亚洲第一电影网av| 两性午夜刺激爽爽歪歪视频在线观看| 1000部很黄的大片| 最新中文字幕久久久久| 色综合站精品国产| 欧美成人a在线观看| 女人十人毛片免费观看3o分钟| 婷婷丁香在线五月| 成人精品一区二区免费| 色综合亚洲欧美另类图片| 久久久色成人| 手机成人av网站| 长腿黑丝高跟| 国产精品久久久久久久久免 | 亚洲美女黄片视频| 欧美xxxx黑人xx丫x性爽| 中文在线观看免费www的网站| 午夜福利高清视频| 搡老妇女老女人老熟妇| 一级毛片女人18水好多| 亚洲精品成人久久久久久| 亚洲五月婷婷丁香| 色老头精品视频在线观看| 亚洲精品美女久久久久99蜜臀| 欧美3d第一页| 亚洲成人久久性| 欧美日韩乱码在线| 亚洲在线观看片| 国产三级在线视频| 真人一进一出gif抽搐免费| 国产午夜精品论理片| 宅男免费午夜| 亚洲第一欧美日韩一区二区三区| 在线观看av片永久免费下载| 桃红色精品国产亚洲av| 日本熟妇午夜| 嫩草影院精品99| av黄色大香蕉| 丁香欧美五月| 热99re8久久精品国产| 国内精品一区二区在线观看| 有码 亚洲区| 亚洲久久久久久中文字幕| 一进一出抽搐gif免费好疼| 18禁美女被吸乳视频| 淫妇啪啪啪对白视频| 久久国产乱子伦精品免费另类| 高清毛片免费观看视频网站| 成人国产综合亚洲| 亚洲性夜色夜夜综合| 久久伊人香网站| 日韩免费av在线播放| 久久久久久久精品吃奶| 国产亚洲精品一区二区www| 高清日韩中文字幕在线| 校园春色视频在线观看| 亚洲真实伦在线观看| 午夜福利在线观看吧| 午夜免费男女啪啪视频观看 | 波多野结衣巨乳人妻| 国产高清激情床上av| 99热6这里只有精品| 欧美最新免费一区二区三区 | 久久亚洲真实| 久久精品亚洲精品国产色婷小说| 国产成人福利小说| 最新美女视频免费是黄的| 村上凉子中文字幕在线| 精品熟女少妇八av免费久了| 九九在线视频观看精品| 欧美成人性av电影在线观看| 国产熟女xx| 岛国视频午夜一区免费看| 国产69精品久久久久777片| 日本黄色视频三级网站网址| 夜夜爽天天搞| 伊人久久精品亚洲午夜| 嫩草影院精品99| 男女那种视频在线观看| 五月伊人婷婷丁香| 亚洲精品一区av在线观看| 久久久久久久精品吃奶| 亚洲午夜理论影院| 亚洲国产欧洲综合997久久,| 可以在线观看的亚洲视频| 亚洲一区二区三区不卡视频| 深夜精品福利| 国产成人aa在线观看| 无限看片的www在线观看| 亚洲av成人av| 男女那种视频在线观看| 国产主播在线观看一区二区| 无人区码免费观看不卡| 亚洲精品在线美女| 精品久久久久久久末码| 婷婷亚洲欧美| а√天堂www在线а√下载| 两个人看的免费小视频| 丁香欧美五月| 日本撒尿小便嘘嘘汇集6| 午夜福利在线观看免费完整高清在 | 欧美性猛交黑人性爽| 成人18禁在线播放| 国产一区在线观看成人免费| 亚洲成a人片在线一区二区| 国产午夜福利久久久久久| 51午夜福利影视在线观看| 成人一区二区视频在线观看| 国产91精品成人一区二区三区| 老司机午夜福利在线观看视频| 精品一区二区三区视频在线观看免费| 欧美成人a在线观看| 国产午夜精品论理片| 国产黄a三级三级三级人| 国产精品永久免费网站| av欧美777| 在线看三级毛片| 欧美日韩黄片免| 男女视频在线观看网站免费| 国产精品乱码一区二三区的特点| 舔av片在线| 亚洲天堂国产精品一区在线| 91av网一区二区| 成人性生交大片免费视频hd| 欧美日韩黄片免| 色精品久久人妻99蜜桃| 网址你懂的国产日韩在线| 亚洲性夜色夜夜综合| 精品一区二区三区av网在线观看| 亚洲美女视频黄频| 欧美一级毛片孕妇| 午夜影院日韩av| 操出白浆在线播放| 久久国产精品影院| 99国产精品一区二区蜜桃av| 亚洲国产精品合色在线| 狠狠狠狠99中文字幕| 国产亚洲精品综合一区在线观看| 国产私拍福利视频在线观看| 蜜桃久久精品国产亚洲av| 久久欧美精品欧美久久欧美| 亚洲人成电影免费在线| 亚洲 欧美 日韩 在线 免费| 国产成人影院久久av| 白带黄色成豆腐渣| 最近最新中文字幕大全免费视频| 老司机深夜福利视频在线观看| 国产精品三级大全| 国产激情偷乱视频一区二区| 亚洲精品美女久久久久99蜜臀| 丰满人妻熟妇乱又伦精品不卡| 久久人人精品亚洲av| 女警被强在线播放| 国产精品综合久久久久久久免费| 麻豆国产av国片精品| 国产精品亚洲一级av第二区| 久久久久国内视频| 免费av不卡在线播放| 亚洲欧美日韩卡通动漫| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕熟女人妻在线| 一区二区三区国产精品乱码| 亚洲av日韩精品久久久久久密| 一级毛片女人18水好多| 少妇的逼好多水| 制服人妻中文乱码| 午夜免费观看网址| 亚洲天堂国产精品一区在线| 丰满人妻一区二区三区视频av | 亚洲真实伦在线观看| 女警被强在线播放| 嫩草影视91久久| 99精品久久久久人妻精品| 全区人妻精品视频| 青草久久国产| 三级国产精品欧美在线观看| 久久久国产成人免费| 午夜福利在线观看免费完整高清在 | 久久午夜亚洲精品久久| 欧美日韩黄片免| 国产精华一区二区三区| 色哟哟哟哟哟哟| 久久九九热精品免费| 在线观看66精品国产| 91久久精品国产一区二区成人 | or卡值多少钱| 中文字幕高清在线视频| 黄色视频,在线免费观看| 国产伦精品一区二区三区四那| 国产av不卡久久| 在线播放国产精品三级| 香蕉av资源在线| 尤物成人国产欧美一区二区三区| 欧美午夜高清在线| 黑人欧美特级aaaaaa片| 国产av不卡久久| 啦啦啦韩国在线观看视频| 一二三四社区在线视频社区8| 亚洲在线自拍视频| 老汉色av国产亚洲站长工具| av专区在线播放| 中文字幕人妻丝袜一区二区| 午夜精品一区二区三区免费看| 日本精品一区二区三区蜜桃| 日日夜夜操网爽| 亚洲av成人av| 婷婷精品国产亚洲av| 欧美在线一区亚洲| 小说图片视频综合网站| 一本久久中文字幕| av在线蜜桃| 欧美乱妇无乱码| 国内精品美女久久久久久| 久久午夜亚洲精品久久| 日韩免费av在线播放| 国产美女午夜福利| 欧美日本亚洲视频在线播放| 精品午夜福利视频在线观看一区| 日本 av在线| 校园春色视频在线观看| 欧美午夜高清在线| 天堂网av新在线| 99在线人妻在线中文字幕| 窝窝影院91人妻| 亚洲内射少妇av| 丁香六月欧美| 叶爱在线成人免费视频播放| 国产探花极品一区二区| 女人十人毛片免费观看3o分钟| 少妇高潮的动态图| 欧美成人一区二区免费高清观看| 在线天堂最新版资源| 男人的好看免费观看在线视频| 99视频精品全部免费 在线| 禁无遮挡网站| 两个人的视频大全免费| 国内精品久久久久久久电影| 97人妻精品一区二区三区麻豆| 精品日产1卡2卡| 搡女人真爽免费视频火全软件 | 搡老岳熟女国产| 老司机福利观看| 日本 av在线| 国产高清videossex| 亚洲,欧美精品.| 啦啦啦观看免费观看视频高清| АⅤ资源中文在线天堂| 日韩欧美在线乱码| а√天堂www在线а√下载| 婷婷亚洲欧美| 美女cb高潮喷水在线观看| 国产午夜福利久久久久久| 国产真实乱freesex| 97碰自拍视频| 美女高潮的动态| eeuss影院久久| 成人鲁丝片一二三区免费| 欧美av亚洲av综合av国产av| 国产精品 欧美亚洲| 日韩欧美一区二区三区在线观看| 亚洲欧美日韩高清在线视频| 亚洲av电影在线进入| 又黄又粗又硬又大视频| 精品久久久久久久末码| 悠悠久久av| 很黄的视频免费| 日韩欧美精品v在线| 婷婷丁香在线五月| 国产精品一及| 免费人成在线观看视频色| 一个人免费在线观看的高清视频| 在线十欧美十亚洲十日本专区| 国产高清有码在线观看视频| 在线观看66精品国产| 男人和女人高潮做爰伦理| 一进一出抽搐gif免费好疼| 亚洲欧美激情综合另类| 国产精品香港三级国产av潘金莲| 国产高清三级在线| 国产精品爽爽va在线观看网站| 好看av亚洲va欧美ⅴa在| 最近最新中文字幕大全免费视频| 中文亚洲av片在线观看爽| 国产淫片久久久久久久久 | 男女做爰动态图高潮gif福利片| 久久久久精品国产欧美久久久| 深夜精品福利| 亚洲国产日韩欧美精品在线观看 | 亚洲成人免费电影在线观看| 亚洲成人久久爱视频| 亚洲美女黄片视频| 黄色丝袜av网址大全| 国产精品久久久久久人妻精品电影| 老司机在亚洲福利影院| 国产精品嫩草影院av在线观看 | 免费在线观看日本一区| 嫩草影院精品99| 色老头精品视频在线观看| 高清毛片免费观看视频网站| 亚洲精品亚洲一区二区| 久久性视频一级片| 国内少妇人妻偷人精品xxx网站| 亚洲欧美激情综合另类| 国产主播在线观看一区二区| 免费观看精品视频网站| 欧美一区二区亚洲| 欧美午夜高清在线| 日本精品一区二区三区蜜桃| 亚洲在线观看片| 精品久久久久久久毛片微露脸| 成人高潮视频无遮挡免费网站| 国产精品久久视频播放| 日韩免费av在线播放| www国产在线视频色| 一进一出好大好爽视频| 一区二区三区国产精品乱码| 亚洲国产精品成人综合色| 午夜福利在线观看免费完整高清在 | 国产中年淑女户外野战色| 亚洲精品成人久久久久久| 最近最新中文字幕大全电影3| 99精品欧美一区二区三区四区| 午夜视频国产福利| 99国产精品一区二区蜜桃av| 三级国产精品欧美在线观看| 欧美中文日本在线观看视频| 国产精品一及| netflix在线观看网站| 成人午夜高清在线视频| 欧美三级亚洲精品| aaaaa片日本免费| 国产不卡一卡二| xxx96com| 中文字幕熟女人妻在线| 三级男女做爰猛烈吃奶摸视频| 成人性生交大片免费视频hd| www日本黄色视频网| 国产在线精品亚洲第一网站| 偷拍熟女少妇极品色| 一级黄色大片毛片| 在线观看日韩欧美| 国产不卡一卡二| 日韩大尺度精品在线看网址| 精品免费久久久久久久清纯| 999久久久精品免费观看国产| 国产亚洲精品久久久com| 九九在线视频观看精品| 日韩有码中文字幕| 午夜a级毛片| 成人一区二区视频在线观看| 亚洲av熟女| 999久久久精品免费观看国产| 国产69精品久久久久777片| 亚洲激情在线av| 亚洲午夜理论影院| 日韩高清综合在线| 人人妻,人人澡人人爽秒播| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av日韩精品久久久久久密| 免费无遮挡裸体视频| 国产欧美日韩一区二区精品| 国产一区在线观看成人免费| 高潮久久久久久久久久久不卡| 男女之事视频高清在线观看| 日本黄色视频三级网站网址| 国产精品亚洲av一区麻豆| 99热只有精品国产| 一个人看视频在线观看www免费 | 无遮挡黄片免费观看| 国产aⅴ精品一区二区三区波| 综合色av麻豆| 在线观看免费视频日本深夜| xxxwww97欧美| 岛国视频午夜一区免费看| 又黄又粗又硬又大视频| 久久久色成人| 色噜噜av男人的天堂激情| 成人欧美大片| 亚洲黑人精品在线| 此物有八面人人有两片| 日本三级黄在线观看| 99精品欧美一区二区三区四区| 老司机深夜福利视频在线观看| 亚洲av免费高清在线观看| 美女免费视频网站| 一级毛片女人18水好多| 日韩av在线大香蕉| 日韩精品青青久久久久久| 亚洲不卡免费看| 久久久久久国产a免费观看| 国产精品1区2区在线观看.| 亚洲人成伊人成综合网2020| 国产精品,欧美在线| 熟妇人妻久久中文字幕3abv| 麻豆国产av国片精品| 最新中文字幕久久久久| 色综合婷婷激情| www.www免费av| 男女午夜视频在线观看| 观看免费一级毛片| а√天堂www在线а√下载| 91久久精品电影网| 日本在线视频免费播放| 欧美性感艳星| 999久久久精品免费观看国产| 精品午夜福利视频在线观看一区| 久久香蕉国产精品| 亚洲专区国产一区二区| 免费在线观看亚洲国产| 日韩欧美三级三区| 伊人久久精品亚洲午夜| 操出白浆在线播放| 午夜免费观看网址| 伊人久久精品亚洲午夜| 女警被强在线播放| 国产私拍福利视频在线观看| 一个人看视频在线观看www免费 | 午夜福利欧美成人| 午夜日韩欧美国产| 久久性视频一级片| 男人舔女人下体高潮全视频| 国产激情偷乱视频一区二区| 99riav亚洲国产免费| 国产成人系列免费观看| 国产激情偷乱视频一区二区| 国产又黄又爽又无遮挡在线| 黄色视频,在线免费观看| 国内精品美女久久久久久| 亚洲avbb在线观看| 乱人视频在线观看| 老汉色∧v一级毛片| 亚洲久久久久久中文字幕| 听说在线观看完整版免费高清| 熟女电影av网| 亚洲av不卡在线观看| 一个人免费在线观看的高清视频| 亚洲 欧美 日韩 在线 免费| 欧美极品一区二区三区四区| 国产一区二区三区在线臀色熟女| 国产99白浆流出| 国产久久久一区二区三区| 超碰av人人做人人爽久久 | 男插女下体视频免费在线播放| 国产真人三级小视频在线观看| 亚洲av中文字字幕乱码综合| 男人的好看免费观看在线视频| 国产不卡一卡二| 91九色精品人成在线观看| 国产精品嫩草影院av在线观看 | 午夜两性在线视频| 成人亚洲精品av一区二区| 舔av片在线|