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

    基于高度計(jì)輔助的GNSS欺騙干擾檢測(cè)

    2024-11-22 00:00:00徐奕禹陳長(zhǎng)風(fēng)袁雪林陳正坤周志健
    關(guān)鍵詞:信號(hào)檢測(cè)模型

    摘 要:

    氣壓高度計(jì)可以獲得精確高程數(shù)據(jù),且不受外界電磁干擾,利用氣壓高度計(jì)輔助全球?qū)Ш叫l(wèi)星系統(tǒng)(global navigation satellite system, GNSS)進(jìn)行欺騙檢測(cè)具有顯著優(yōu)勢(shì)。為了評(píng)估氣壓高度計(jì)在GNSS欺騙檢測(cè)中的能力并提高檢測(cè)概率,首先對(duì)系統(tǒng)誤差進(jìn)行建模分析,并利用差分整合移動(dòng)平均自回歸(autoregressive integrated moving average, ARIMA)模型對(duì)氣壓高度計(jì)和GNSS測(cè)量結(jié)果進(jìn)行處理。同時(shí),通過計(jì)算歐氏距離構(gòu)建數(shù)據(jù)空間,提出一種基于移動(dòng)方差的自適應(yīng)門限檢測(cè)方法。試驗(yàn)結(jié)果表明,在虛警率為10-6的情況下,檢測(cè)率達(dá)到了92%,所提方法能夠在高程信息上有效監(jiān)測(cè)潛在的欺騙威脅,增強(qiáng)了GNSS系統(tǒng)的可靠性。

    關(guān)鍵詞:

    氣壓高度計(jì); 衛(wèi)星導(dǎo)航; 高度數(shù)據(jù)融合; 移動(dòng)方差; 欺騙檢測(cè)

    中圖分類號(hào):

    TN 967.1

    文獻(xiàn)標(biāo)志碼: A""" DOI:10.12305/j.issn.1001-506X.2024.05.03

    GNSS spoofing jamming detection based on altimeter assistance

    XU Yiyu1, CHEN Changfeng2, YUAN Xuelin1,*, CHEN Zhengkun1, ZHOU Zhijian1

    (1. School of Electronic and Communication Engineering, Sun Yat-Sen University, Shenzhen 518107, China;

    2. Southwest China Research Institute of Electronic Equipment, Chengdu 610036, China)

    Abstract:

    The barometric altimeter can obtain accurate elevation data without being affected by external electromagnetic interference. The use of barometric altimeters to assist the global navigation satellite system(GNSS) in spoofing detection has significant advantages. In order to evaluate the ability of the barometric altimeter in GNSS spoofing detection and improve detection probability, the system error is first modeled and analyzed, and the differential integrated autoregressive moving average (ARIMA) model is used to process the barometric altimeter and GNSS measurement results. Meanwhile, a self-adaptive threshold detection method based on moving variance is proposed by calculating Euclidean distance to construct a data space. The experimental results show that with a 1 alarm rate of 10-6, the detection rate reaches 92%, and the proposed method can effectively monitor potential spoofing threats on elevation information, which enhances the reliability of the GNSS system.

    Keywords:

    barometric altimeter; satellite navigation; altitude data fusion; moving variance; spoofing detection

    0 引 言

    目前,隨著全球?qū)Ш叫l(wèi)星系統(tǒng)(global navigation sate-llite system, GNSS)的發(fā)展,衛(wèi)星導(dǎo)航技術(shù)已被廣泛應(yīng)用于各行各業(yè),支撐著許多現(xiàn)代系統(tǒng),可以為智能終端、汽車、飛機(jī)、游輪等提供定位和導(dǎo)航服務(wù),也可以為電網(wǎng)和通信網(wǎng)絡(luò)等提供授時(shí)服務(wù)[1]。然而,衛(wèi)星信號(hào)仍然具有許多明顯的缺陷,例如容易受到干擾和欺騙,這就將其脆弱性暴露了出來(lái)。

    欺騙干擾指的是通過向目標(biāo)設(shè)備發(fā)射由欺騙干擾機(jī)生成的虛假信號(hào)或經(jīng)相關(guān)處理的真實(shí)衛(wèi)星信號(hào),誘使接收機(jī)鎖定該欺騙信號(hào),并誤導(dǎo)目標(biāo)得到錯(cuò)誤的導(dǎo)航定位結(jié)果[2-3],以此達(dá)到欺騙目的。

    近年來(lái),欺騙干擾的案例層出不窮。例如,2017年,在黑海上航行的數(shù)十艘船只報(bào)告稱,其全球定位系統(tǒng)(global positioning system, GPS)設(shè)備出現(xiàn)故障,定位設(shè)備顯示這些船只??吭趦?nèi)陸。2019年,芬蘭某些機(jī)場(chǎng)的GPS信號(hào)受到了有意干擾,導(dǎo)致大批民航客機(jī)上的導(dǎo)航設(shè)備受到影響,飛機(jī)滯飛。在無(wú)人駕駛領(lǐng)域,多項(xiàng)研究表明大部分無(wú)人駕駛汽車存在GNSS安全性漏洞。2020年,在特斯拉駕駛測(cè)試過程中,一輛使用自動(dòng)駕駛導(dǎo)航(navigate on autopilot, NOA)系統(tǒng)的特斯拉Model 3由于被欺騙而駛離了預(yù)定路線。以上這些案例都有力地證明了GNSS信號(hào)的脆弱性[4],因此研究GNSS欺騙干擾檢測(cè)方法對(duì)GNSS安全、可靠而言意義重大。

    目前,將輔助無(wú)源器件結(jié)合GNSS進(jìn)行欺騙檢測(cè)已經(jīng)成為研究熱點(diǎn)。Lee等[5]利用加速度計(jì)并通過監(jiān)測(cè)其輸出信息與GPS輸出參數(shù)是否一致,實(shí)現(xiàn)對(duì)欺騙信號(hào)的檢測(cè)。Curran等[6]研究了慣性測(cè)量單元(inertial measurement unit, IMU)在GNSS欺騙檢測(cè)中的應(yīng)用,并使用低成本的IMU輔助進(jìn)行欺騙檢測(cè)。Lo 等[7]則使用了低成本加速度計(jì)來(lái)檢測(cè)慣性傳感器和GNSS位置之間的差異。雖然上述利用外界輔助信息的欺騙檢測(cè)方法在短時(shí)間內(nèi)具有良好的抗欺騙干擾能力,但慣性設(shè)備的精度會(huì)隨時(shí)間的增加而不斷下降,需要定期校準(zhǔn),很難大范圍推廣應(yīng)用。Borio和Lee等[8-9]研究了在智能終端中使用氣壓計(jì)等無(wú)源傳感器進(jìn)行GNSS欺騙檢測(cè),但沒有單獨(dú)驗(yàn)證氣壓計(jì)所起的作用,而本文更深入地研究了氣壓高度計(jì)在GNSS欺騙檢測(cè)中所起的輔助作用。

    對(duì)于測(cè)高而言,GNSS雖然能夠提供精度較高的用戶終端位置信息,但其輸出頻率較低。在峽谷、隧道、森林等復(fù)雜遮擋環(huán)境中,GNSS信號(hào)容易出現(xiàn)失鎖并造成捕獲困難或信號(hào)丟失,在垂直方向上的精度因子(dilution of precision, DOP)值較差[10-11],此時(shí)若欺騙干擾機(jī)發(fā)射欺騙信號(hào),很容易實(shí)現(xiàn)對(duì)用戶位置等信息的欺騙。

    氣壓高度計(jì)具有成本低、結(jié)構(gòu)簡(jiǎn)單、無(wú)累計(jì)誤差、不受無(wú)線電干擾等優(yōu)點(diǎn)[12-14]。根據(jù)博世官方公布的BMP380氣壓計(jì)參數(shù),該氣壓計(jì)的相對(duì)壓力精度為0.06 hPa,相當(dāng)于0.5 m,遠(yuǎn)小于GPS系統(tǒng)的垂直方向誤差[15]。因此,利用氣壓高度計(jì)輔助GNSS進(jìn)行欺騙檢測(cè)便能彌補(bǔ)GNSS信號(hào)在垂直方向上DOP值較差這一缺陷,提高導(dǎo)航服務(wù)的精度、魯棒性以及連續(xù)性。

    針對(duì)越來(lái)越復(fù)雜的欺騙干擾場(chǎng)景,本文提出了一種氣壓高度計(jì)輔助GNSS進(jìn)行欺騙檢測(cè)的方法,利用整合移動(dòng)平均自回歸(autoregressive integrated moving average, ARIMA)模型對(duì)氣壓高度計(jì)和GNSS測(cè)量結(jié)果進(jìn)行降噪,提高數(shù)據(jù)精度與容錯(cuò)性,并提出了基于移動(dòng)方差的自適應(yīng)門限檢測(cè)方法,該門限可以直接作為檢測(cè)參量進(jìn)行GNSS欺騙信號(hào)檢測(cè)。

    1 系統(tǒng)誤差模型及降噪處理

    對(duì)欺騙信號(hào)的檢測(cè)、驗(yàn)證可以視為二元信號(hào)檢測(cè)問題[16]。H0假設(shè)為欺騙信號(hào)不存在,H1假設(shè)為欺騙信號(hào)存在,這兩個(gè)假設(shè)可以描述如下:

    =0,H0:欺騙信號(hào)不存在

    1,H1:欺騙信號(hào)存在(1)

    在該假設(shè)下,有4種可能出現(xiàn)的判決結(jié)果,本文只考慮利用虛警概率和檢測(cè)概率作為評(píng)價(jià)系統(tǒng)檢測(cè)性能的指標(biāo)。其中,虛警概率Pfa是欺騙信號(hào)不存在而系統(tǒng)做出誤判的概率;檢測(cè)概率PD是欺騙信號(hào)存在且系統(tǒng)做出正確判斷的概率。后文將會(huì)詳細(xì)展開討論檢測(cè)概率與門限的推導(dǎo)。

    數(shù)據(jù)處理過程和欺騙檢測(cè)原理如圖1所示。其中,氣壓高度計(jì)的海拔高度信息可以通過氣壓讀數(shù)并使用壓高方程進(jìn)行計(jì)算獲得,而GNSS的海拔高度信息可以通過NMEA 0183協(xié)議中的全球定位系統(tǒng)固定數(shù)據(jù)(global positioning system fix data, GGA)語(yǔ)句獲取。在獲取數(shù)據(jù)后,通過求解兩者海拔高度間的歐氏距離并推算自適應(yīng)門限,最后進(jìn)行門限檢測(cè),即可判斷欺騙信號(hào)是否存在。

    在實(shí)際場(chǎng)景中,直接得到的高度計(jì)數(shù)據(jù)和導(dǎo)航數(shù)據(jù)會(huì)存在噪聲,并對(duì)最終結(jié)果產(chǎn)生影響,為了使最終結(jié)果精度更高,本文采用ARIMA模型對(duì)信號(hào)進(jìn)行降噪處理,首先對(duì)各模型進(jìn)行分析。

    1.1 載體運(yùn)動(dòng)模型

    針對(duì)手機(jī)等接收載體,首先將載體的運(yùn)動(dòng)狀態(tài)考慮進(jìn)來(lái),主要在垂直方向上建立載體運(yùn)動(dòng)狀態(tài)模型,其模型如下:

    hreal=hideal+δa(2)

    載體運(yùn)動(dòng)時(shí)的加速度可以看作一個(gè)零均值的平穩(wěn)一階馬爾可夫過程[17-18],可表示為

    jideal=-1τaaideal+na(3)

    式中:hreal為載體的真實(shí)高度;δa為高度測(cè)量誤差;hideal為載體的理想運(yùn)動(dòng)高度;aideal為載體運(yùn)動(dòng)時(shí)的加速度;τa為一階馬爾可夫過程的相關(guān)時(shí)間;na為白噪聲。

    1.2 氣壓-高度計(jì)算模型

    海拔高度是指地面上某載體的位置高出海平面的垂直距離。用氣壓計(jì)測(cè)高,其基本原理是:在重力場(chǎng)內(nèi)隨著海拔的增加,氣壓值呈指數(shù)減少,而當(dāng)氣壓與溫度分布標(biāo)準(zhǔn)化時(shí),可以認(rèn)為氣壓高度計(jì)測(cè)得的高度即為海拔高度[19]。據(jù)實(shí)測(cè),在地面層中,高度每升高100 m,氣壓平均降低0.93 kPa,在高層則小于此數(shù)值。表1展示了標(biāo)準(zhǔn)大氣下等壓面和高度的關(guān)系。

    由于氣壓隨高度的變化會(huì)受到諸多因素的影響,常見的影響因素有溫度、時(shí)間、緯度等,因此國(guó)際上規(guī)定使用國(guó)際標(biāo)準(zhǔn)大氣ISO2533-1975來(lái)避免這些因素的影響。ISO2533-1975規(guī)定下的大氣滿足理想氣體方程,同時(shí)海拔高度以平均海平面作為參考高度,即零高度。在該條件下,某一點(diǎn)氣溫與高度的關(guān)系[20-21]可以表示為

    Tr=T0+β(Hr-H0)(4)

    式中:H0為標(biāo)準(zhǔn)海平面的高度,規(guī)定為0 m;Hr為待測(cè)點(diǎn)的海拔高度;T0為標(biāo)準(zhǔn)海平面的氣溫,取288.15 K;β為溫度梯度,取-0.006 5 K/m。

    在實(shí)際應(yīng)用中,可以用壓高方程來(lái)表示氣壓與高度之間的變化關(guān)系[22-24]:

    Hr=H0+T0βPrP0-βRg-1(5)

    式中:R為空氣氣體常數(shù),取287.05 m2/(s2·K);g為重力加速度,取9.8 m/s2;P0為H0處的氣壓值,取101.325 kPa;Pr為氣壓高度計(jì)實(shí)時(shí)測(cè)量的氣壓值。

    將各數(shù)值代入式(5)并化簡(jiǎn),可得待測(cè)點(diǎn)海拔高度Hr與待測(cè)點(diǎn)氣壓Pr的關(guān)系滿足:

    Hr=44 331×1-Pr101.3250.19(6)

    當(dāng)實(shí)際生活中有高度測(cè)量的需求時(shí),利用此壓高公式,只需代入氣壓高度計(jì)的氣壓值,便可直接獲得該點(diǎn)的海拔高度。

    1.3 氣壓-高度誤差修正模型

    由于國(guó)際標(biāo)準(zhǔn)大氣條件適用于待測(cè)點(diǎn)高度較低且溫度與參考平面變化幅度不大的情況。而在實(shí)際測(cè)量中,當(dāng)待測(cè)點(diǎn)的實(shí)際大氣條件不滿足國(guó)際標(biāo)準(zhǔn)大氣條件時(shí),氣壓-高度計(jì)算模型便會(huì)存在相應(yīng)的原理性誤差,需要進(jìn)行誤差修正。

    當(dāng)前已知的修正方法為修正海平面氣壓法[25-27],其原理為:將氣壓高度計(jì)輸出的氣壓值,按標(biāo)準(zhǔn)大氣條件修正到當(dāng)前平均海平面下的氣壓和高度,即利用氣壓高度計(jì)輸出的實(shí)際氣壓P1與實(shí)際高度H1對(duì)標(biāo)準(zhǔn)大氣條件下的氣壓P0和海拔H0進(jìn)行修正,得到當(dāng)前平均海平面氣壓近似值P′0及海拔近似值H′0。

    修正過程首先以當(dāng)前平均海平面為基準(zhǔn),將實(shí)際氣壓值P1和高度值H1代入式(5),可得

    H1=H′0+T′0βP1P′0-βRg-1(7)

    式中:T′0為當(dāng)前平均海平面的溫度。再以標(biāo)準(zhǔn)海平面為基準(zhǔn),將平均海平面參數(shù)P′0和H′0代入式(5),可得

    H′0=H0+T0βP′0P0-βRg-1(8)

    同時(shí),利用當(dāng)前平均海平面的溫度T′0修正標(biāo)準(zhǔn)海平面基準(zhǔn)溫度T0,即有:

    T′0=T0+β(H′0-H0)(9)

    聯(lián)立式(7) ~式(9),便可求解得到修正后基準(zhǔn)點(diǎn)的參考?xì)鈮篜′0、高度值H′0與溫度T′0。此時(shí),式(5)可修正為

    Hr=H′0+T′0βPrP′0-βRg-1(10)

    由此,便得到了修正后的氣壓-高度計(jì)算模型,利用修正后新基準(zhǔn)點(diǎn)的參數(shù),可將氣壓高度計(jì)讀取的氣壓值P1代入式(10),重新計(jì)算得到氣壓高度計(jì)當(dāng)前所在的高度H2,即可得到當(dāng)前待測(cè)點(diǎn)的實(shí)際海拔高度。

    1.4 GNSS高程誤差模型

    GPS/北斗衛(wèi)星導(dǎo)航系統(tǒng)(Beidou navigation satellite system, BDS)采用的是WGS(world geodetic system)-84/2000國(guó)家大地坐標(biāo)系(China geodetic coordinate system, CGCS),該坐標(biāo)系下的基準(zhǔn)面為參考橢球面,待測(cè)點(diǎn)的高度為沿著橢球面的法線,該點(diǎn)到橢球面的距離[28],也即大地高。

    但導(dǎo)航定位描述物體運(yùn)動(dòng)姿態(tài)、位置和速度時(shí),是建立在給定參考系的基礎(chǔ)之上的,當(dāng)參考系不同時(shí),物體的運(yùn)動(dòng)速度、位置和姿態(tài)等信息也隨之不同,因此必須進(jìn)行坐標(biāo)轉(zhuǎn)換,不能直接將大地高用于實(shí)際測(cè)量。

    本文采用GNSS海拔高度與氣壓高度計(jì)的海拔高度進(jìn)行試驗(yàn)驗(yàn)證,因此首先進(jìn)行高度轉(zhuǎn)換,將GNSS大地高轉(zhuǎn)換為海拔高度。高度之間的轉(zhuǎn)換關(guān)系如圖2所示,其中Hp為大地高,hideal為海拔高,Δh為橢球面與大地水準(zhǔn)面的高度之差,可近似認(rèn)為Hp=hideal+Δh。

    一般而言,載體的海拔高度可以通過計(jì)算載體坐標(biāo)位置并轉(zhuǎn)換到用戶坐標(biāo)系下得到。此時(shí),可將GNSS的高度測(cè)量方程表示為

    Hw=hideal+ns(11)

    式中:ns為測(cè)量引入的噪聲,視為白噪聲。

    1.5 基于ARIMA的信號(hào)降噪方法

    在GNSS與氣壓高度計(jì)測(cè)量過程中,會(huì)有干擾噪聲夾雜其中,為了減少干擾噪聲的影響,需要對(duì)GNSS與氣壓高度計(jì)測(cè)量的數(shù)據(jù)進(jìn)行降噪處理,而ARIMA降噪算法便可以有效過濾非廣義平穩(wěn)信號(hào)中的隨機(jī)噪聲[29]。

    信號(hào)降噪的具體原理為:一定階數(shù)的ARIMA(p, q)模型可以擬合任意廣義平穩(wěn)的隨機(jī)過程,若有一平穩(wěn)、服從正態(tài)分布且零均值的時(shí)間序列{xt},t=1,2,…,N,則該序列的ARIMA(p, q)模型可以表示為

    xt=φ1xt-1+φ2xt-2+…+φpxt-p+

    εt-θ1εt-1-θ2εt-2-…-θqεt-q(12)

    式中:xt為時(shí)間序列;φ1,φ2,…,φp為模型自回歸參數(shù);θ1,θ2,…,θq為模型滑動(dòng)平均參數(shù);εt為白噪聲,服從標(biāo)準(zhǔn)正態(tài)分布;p, q為ARIMA(p, q)模型的階數(shù);N為{xt}的長(zhǎng)度。

    而在許多情況下,{xt}并不能同時(shí)滿足上述提到的3個(gè)預(yù)設(shè)條件,此時(shí){xt}不符合ARIMA(p, q)模型的適用條件,需要對(duì){xt}進(jìn)行d階差分使其滿足ARIMA模型的預(yù)設(shè)要求,然后建立差分后的ARIMA(p, q)模型,即將原序列轉(zhuǎn)換成ARIMA(p, d, q)模型,表達(dá)式為

    Φ(L)Δdxt=δ+Θ(L)εt(13)

    式中:L是后移算子,Lxt=xt-1;Φ(L)=1-∑pi=1φiLi;Δ=1-L;Θ(L)=1-∑qj=1εjLj。

    式(13)表示的ARIMA(p, d, q)模型適用于非廣義平穩(wěn)序列。本文使用ARIMA降噪算法除噪,具體過程如圖3所示。

    本文將根據(jù)所確定的 ARIMA模型來(lái)抑制GNSS高度數(shù)據(jù)和氣壓高度計(jì)數(shù)據(jù)的隨機(jī)誤差。在數(shù)據(jù)處理方面,可認(rèn)為GNSS以真值和白噪聲作為系統(tǒng)輸入,并輸出高度信息作為觀測(cè)值。GNSS與氣壓高度計(jì)經(jīng) ARIMA算法降噪處理后的輸出數(shù)據(jù)將在試驗(yàn)驗(yàn)證部分詳細(xì)展示。

    2 高度計(jì)輔助的GNSS欺騙干擾檢測(cè)

    本節(jié)主要針對(duì)高度計(jì)輔助的GNSS欺騙干擾檢測(cè)展開,檢測(cè)的具體原理如下:當(dāng)接收機(jī)處于一個(gè)靜態(tài)位置并進(jìn)行定位時(shí),其捕獲真實(shí)衛(wèi)星信號(hào),獲取高度數(shù)據(jù),得到衛(wèi)星的測(cè)量高度較為平穩(wěn)。當(dāng)欺騙信號(hào)加入時(shí),若欺騙信號(hào)成功欺騙了接收機(jī),此時(shí)接收機(jī)的高度將會(huì)產(chǎn)生較大的波動(dòng)或是跳變,而氣壓高度計(jì)的高度不會(huì)受到欺騙干擾,仍然在原高度處平穩(wěn)變化?;谄垓_信號(hào)與真實(shí)信號(hào)的這一特點(diǎn),試驗(yàn)設(shè)置為:在起始階段,接收機(jī)正常接收一段真實(shí)的衛(wèi)星信號(hào),在試驗(yàn)的后半段加入欺騙信號(hào)來(lái)欺騙接收機(jī),并將接收機(jī)的高度欺騙至預(yù)定高度,最后將衛(wèi)星高度數(shù)據(jù)與高度計(jì)實(shí)測(cè)數(shù)據(jù)進(jìn)行處理及對(duì)比分析。

    首先,需要計(jì)算GNSS高度數(shù)據(jù)與氣壓高度計(jì)數(shù)據(jù)之間的歐氏距離。歐氏距離即二維空間中兩點(diǎn)的真實(shí)距離,其定義如下:

    dist(i)=(xi-yi)2, i=1,2,…,n(14)

    在真實(shí)場(chǎng)景中,由于高度數(shù)據(jù)變化穩(wěn)定,dist(i)可以近似認(rèn)為服從高斯分布,即dist(i)~N(μR,σ2R)。構(gòu)造變量mR:

    mR=(w-1)·ησ2R(15)

    式中:η是寬度為w的時(shí)間窗內(nèi)的移動(dòng)方差計(jì)算值,表示為

    η=1w∑nk=n-w+1(dist(k)-μR)2(16)

    可知mR服從卡方分布[30],即有

    mR=(w-1)·ησ2R∶χ2(k)(17)

    式中:k為自由度,且k=w-1。根據(jù)卡方分布的統(tǒng)計(jì)特性,mR的均值為k,方差為2k,且mR的概率密度函數(shù)為

    fχ2(mR)=mk2-1Re-mR2, mR>0

    0, 其他(18)

    由此,得到mR的累積分布函數(shù):

    Fχ2(mR;k)=γk2,mR2Γk2=Pk2,mR2(19)

    式中:γ(·)是下不完全Gamma函數(shù);P(s,t)是正則化Gamma函數(shù)。因此,mR的虛警率為

    Pfa=∫+∞Tmfχ2(mR)dmR=1-Fχ2(Tm;k)(20)

    由式(20)便可求出固定虛警率時(shí)的門限Tm,進(jìn)而求解得到移動(dòng)方差門限Th:

    Th=σ2RTmw-1(21)

    式中:w為時(shí)間窗寬度。

    得到移動(dòng)方差門限后,便可計(jì)算檢測(cè)概率PD,這需要獲取欺騙環(huán)境下的η分布,該分布取決于欺騙攻擊的程度。由于接收機(jī)無(wú)法獲取欺騙信號(hào)的功率變化信息,當(dāng)存在欺騙信號(hào)時(shí),η的分布較為復(fù)雜,無(wú)法通過推導(dǎo)其概率密度函數(shù)來(lái)計(jì)算PD。因此,本文使用檢測(cè)率DR來(lái)替代PD。DR定義為時(shí)間窗內(nèi)η(n)超過門限的次數(shù)與總樣本數(shù)之比,計(jì)算DR的具體方法如下:通過Pfa=∫+∞Tmfχ2(mR)dmR=1-Fχ2(Tm;k)計(jì)算出給定Pfa時(shí)的閾值Tm,并根據(jù)Th=σ2RTm/(w-1)計(jì)算得到移動(dòng)方差門限Th,最后在每個(gè)檢測(cè)窗口處進(jìn)行欺騙檢測(cè),計(jì)算每個(gè)窗口處移動(dòng)方差超過門限的個(gè)數(shù)。即

    DR=Num{(η(n)>Th)}N(22)

    此時(shí),二元信號(hào)檢測(cè)問題進(jìn)一步表述為,假設(shè)H0對(duì)應(yīng)為真實(shí)信號(hào)與欺騙信號(hào)間的移動(dòng)方差且沒有超過預(yù)先設(shè)定的門限值,即欺騙攻擊不存在。假設(shè)H1對(duì)應(yīng)為真實(shí)信號(hào)與欺騙信號(hào)間的移動(dòng)方差且超過預(yù)先設(shè)定的門限值,即欺騙攻擊存在。這兩個(gè)假設(shè)可以描述如下:

    =var(real,spoofing)<Th,H0

    var(real,spoofing)≥Th,H1(23)

    欺騙檢測(cè)流程圖如圖4所示。首先獲取衛(wèi)星導(dǎo)航系統(tǒng)與氣壓高度計(jì)的海拔高度數(shù)據(jù),并利用ARIMA對(duì)兩種數(shù)據(jù)進(jìn)行降噪處理,隨后求解歐式距離,并計(jì)算歐氏距離的移動(dòng)方差與方差門限,最后進(jìn)行欺騙檢測(cè),判斷條件為真實(shí)信號(hào)與欺騙信號(hào)之間的移動(dòng)方差是否達(dá)到門限。若達(dá)到門限,則認(rèn)為存在欺騙攻擊;否則,則認(rèn)為不存在欺騙攻擊。

    3 試驗(yàn)驗(yàn)證

    3.1 高度信號(hào)降噪處理

    本節(jié)對(duì)接收機(jī)采集的GNSS高度數(shù)據(jù)和氣壓高度計(jì)數(shù)據(jù)進(jìn)行ARIMA降噪處理。采用的氣壓高度計(jì)為BMP280,參數(shù)如表2所示,可知BMP280在垂直方向上精度較好。圖5則展示了本文所使用的氣壓高度計(jì)BMP280的實(shí)物圖。

    進(jìn)一步分析ARIMA模型對(duì)GNSS高度數(shù)據(jù)與氣壓高度計(jì)數(shù)據(jù)降噪的效果,如圖6和圖7所示。在降噪前后,GNSS高度數(shù)據(jù)與氣壓高度計(jì)數(shù)據(jù)的曲線整體誤差性質(zhì)無(wú)明顯變化,但經(jīng)ARIMA降噪處理后,噪聲方差明顯減小。在GNSS高度數(shù)據(jù)方面,降噪處理前后的均值分別為66.294 6和66.294 4,均在GNSS精度允許范圍內(nèi)。在氣壓高度計(jì)數(shù)據(jù)方面,降噪處理前后的均值分別為66.022 7和66.022 6,其前后均值幾乎是一樣的。

    另外,從表3也可得知,經(jīng)ARIMA降噪處理后的GNSS高度數(shù)據(jù)和氣壓高度計(jì)數(shù)據(jù)標(biāo)準(zhǔn)差更小,遠(yuǎn)優(yōu)于原始數(shù)據(jù)的標(biāo)準(zhǔn)差。這說(shuō)明該降噪處理方法能夠有效抑制GNSS與氣壓高度計(jì)實(shí)測(cè)數(shù)據(jù)中的隨機(jī)噪聲,提高了數(shù)據(jù)精度和容錯(cuò)性。

    3.2 欺騙檢測(cè)結(jié)果分析

    本節(jié)設(shè)置如下試驗(yàn)來(lái)檢測(cè)接收機(jī)接收到的衛(wèi)星信號(hào)中是否存在欺騙信號(hào):在接收機(jī)接收的前半段時(shí)間內(nèi)播發(fā)真實(shí)衛(wèi)星信號(hào),在后半段時(shí)間內(nèi)加入欺騙衛(wèi)星信號(hào),并將得到的高度數(shù)據(jù)進(jìn)行處理分析。試驗(yàn)過程如圖8所示,前300 s時(shí)間內(nèi)為接收機(jī)接收真實(shí)信號(hào)的情況,可以看到GNSS高度值與氣壓高度計(jì)高度值基本在一個(gè)小范圍內(nèi)平穩(wěn)變化。在300 s前后加入欺騙信號(hào)后,接收機(jī)接收到的GNSS高度值發(fā)生了明顯的跳變,由原來(lái)的66 m附近的定位高度跳變到75 m附近。由此可見,當(dāng)加入欺騙信號(hào)時(shí),衛(wèi)星高度數(shù)據(jù)對(duì)比原真實(shí)衛(wèi)星數(shù)據(jù)將會(huì)產(chǎn)生較大變化。而高度計(jì)數(shù)據(jù)在欺騙信號(hào)加入前后均無(wú)較大起伏波動(dòng),說(shuō)明此時(shí)接收機(jī)的定位位置無(wú)明顯變化,參照GNSS高度數(shù)據(jù)和氣壓高度計(jì)測(cè)高數(shù)據(jù),可以初步判斷有欺騙信號(hào)加入。

    本文采用的是恒虛警率檢測(cè),并設(shè)定虛警率Pfa為10-6,窗口寬度為35 s,并根據(jù)式(21)計(jì)算移動(dòng)方差門限。由圖9所示,在欺騙信號(hào)未加入時(shí),GNSS與氣壓高度計(jì)高度值之間的歐式距離方差值始終未超過門限,在欺騙信號(hào)加入的這段時(shí)間內(nèi),兩者歐氏距離的方差值發(fā)生了明顯激增,由原來(lái)0.03 m2左右的移動(dòng)方差激增到峰值為0.78 m2的移動(dòng)方差,遠(yuǎn)遠(yuǎn)超過了自適應(yīng)門限Th。因此,判定此時(shí)有欺騙攻擊的存在。

    移動(dòng)方差算法檢測(cè)效果如圖10所示,從圖中可知,基于移動(dòng)方差的算法檢測(cè)率在前270 s內(nèi)基本為0,在270 s至320 s時(shí)間范圍內(nèi)的檢測(cè)率幾乎是100%。這是由于欺騙信號(hào)是在300 s前后加入的,此時(shí)GNSS與氣壓高度計(jì)高度值之間歐氏距離的方差發(fā)生了突變。在320 s之后,由于已經(jīng)將接收機(jī)欺騙到預(yù)定高度,因此檢測(cè)率與前270 s的一致,基本為0。

    本文使用接受者操作特征(receiver operating characteristic, ROC)曲線來(lái)衡量欺騙信號(hào)的檢測(cè)效果。ROC曲線描述了真正率(true positive rate, TPR)和假正率(1 positive rate, FPR)之間的關(guān)系。對(duì)于欺騙檢測(cè)而言,可以將其簡(jiǎn)化為是否存在欺騙信號(hào),也即分類問題。TPR和FPR定義分別如下:

    TPR=TPTP+FN(24)

    FPR=FPFP+TN(25)

    式中:TP為真正,其對(duì)應(yīng)真值為1且預(yù)測(cè)也為1的正確預(yù)測(cè);FP為假正,其對(duì)應(yīng)真值為0且預(yù)測(cè)為1的錯(cuò)誤預(yù)測(cè);TN為真負(fù),其對(duì)應(yīng)真值為0且預(yù)測(cè)也為0的正確預(yù)測(cè);FN為假負(fù),其對(duì)應(yīng)真值為1且預(yù)測(cè)為0的錯(cuò)誤預(yù)測(cè)。

    ROC曲線下面積(area under curve, AUC)用于對(duì)模型的性能以及整體準(zhǔn)確性進(jìn)行評(píng)估。AUC值越接近1,模型性能越好[31],AUC定義如下:

    AUC=S0-n0(n0+1)/2n0n1(26)

    式中:n0是正樣本的數(shù)量;n1是負(fù)樣本的數(shù)量,且S0=∑n0i=1ri。

    虛警率從10-6量級(jí)逐漸增大的ROC曲線如圖11所示??梢钥吹剑?dāng)時(shí)間窗寬度從15 s逐步增加到40 s時(shí),AUC區(qū)域在窗口寬度為35 s時(shí)達(dá)到最大,隨后開始減小。在虛警率為10-6、窗口寬度為35 s時(shí),對(duì)應(yīng)的ROC曲線檢測(cè)率達(dá)到了92%

    表4展示了各個(gè)窗口寬度下的AUC值,可以看到,在時(shí)間窗寬度為35 s時(shí),AUC值達(dá)到了0.982,優(yōu)于其他時(shí)間窗寬度。

    為了更直觀地展現(xiàn)檢測(cè)率與時(shí)間窗寬度、虛警率的關(guān)系,圖12展示了在欺騙檢測(cè)整體持續(xù)時(shí)間內(nèi)的檢測(cè)結(jié)果,其中時(shí)間窗口從2 s依次增加至50 s。隨著時(shí)間窗口大小的增加,欺騙信號(hào)檢測(cè)效果也有一定的改善。時(shí)間窗寬度為35 s時(shí),AUC值達(dá)到了0.982,此時(shí)檢測(cè)效果最佳。

    混淆矩陣可以用于描述模型的分類性能,本文將其用于對(duì)檢測(cè)結(jié)果進(jìn)行直觀的性能評(píng)估?;煜仃嚨?個(gè)單元格分別對(duì)應(yīng)真正、真負(fù)、假正和假負(fù),如表5所示。

    實(shí)際上,ROC曲線上的每個(gè)點(diǎn)都對(duì)應(yīng)一個(gè)混淆矩陣。當(dāng)時(shí)間窗口寬度為35 s時(shí),F(xiàn)PR在10% ~ 90%時(shí)效果相近,因此本文選取FPR為10%時(shí)的數(shù)據(jù)點(diǎn)繪制混淆矩陣。圖13則對(duì)應(yīng) ROC 曲線上 FPR(即虛警率為10%時(shí))的數(shù)據(jù)點(diǎn)所對(duì)應(yīng)的混淆矩陣,其中標(biāo)簽的0和1分別代表欺騙信號(hào)不存在和欺騙信號(hào)存在兩種情況??梢?,算法對(duì)欺騙信號(hào)正確識(shí)別的檢測(cè)率達(dá)到了94.12%,效果較好。

    在二元欺騙檢測(cè)任務(wù)中,通常用以下幾個(gè)評(píng)價(jià)指標(biāo)評(píng)價(jià)欺騙信號(hào)的檢測(cè)效果,分別是準(zhǔn)確率Accuracy、精確率Precision、召回率Recall和F1分?jǐn)?shù)[32],具體表示如下:

    Accuracy=TP+TNTP+FN+FP+TN(27)

    Precision=TPTP+FP(28)

    Recall=TPTP+FN(29)

    F1=2·Precision·RecallPrecision+Recall(30)

    式中:準(zhǔn)確率指系統(tǒng)正確預(yù)測(cè)的導(dǎo)航信號(hào)占總導(dǎo)航信號(hào)數(shù)的百分比;精確率指的是在所有被預(yù)測(cè)為欺騙信號(hào)的導(dǎo)航信號(hào)中,實(shí)際為欺騙信號(hào)的概率;召回率指實(shí)際為欺騙信號(hào)的所有導(dǎo)航信號(hào)中,被預(yù)測(cè)為欺騙信號(hào)的概率;F1分?jǐn)?shù)則綜合兼顧了精確率和召回率這兩個(gè)指標(biāo)。

    表6顯示了不同窗口寬度下的欺騙檢測(cè)效果,從表6可知,在窗口寬度為35 s時(shí),欺騙檢測(cè)效果優(yōu)于其他窗口寬度,同時(shí)召回率達(dá)到了94.12%。值得一提的是,由于欺騙場(chǎng)景的總時(shí)長(zhǎng)相比于無(wú)欺騙場(chǎng)景的時(shí)長(zhǎng)更短,精確率總體較低,在這種情況下,即使是在窗口寬度為35 s的情況下,F(xiàn)1分?jǐn)?shù)也只能達(dá)到75%。

    4 結(jié) 論

    本文使用ARIMA模型對(duì)GNSS數(shù)據(jù)與高度計(jì)數(shù)據(jù)進(jìn)行降噪,提高了GNSS高度數(shù)據(jù)與氣壓高度計(jì)數(shù)據(jù)的精度,同時(shí)提出了一種基于移動(dòng)方差的自適應(yīng)門限檢測(cè)方法。試驗(yàn)結(jié)果表明,在移動(dòng)窗口寬度為35 s、虛警率為10-6的情況下,對(duì)欺騙信號(hào)的檢測(cè)率達(dá)到了92%。

    低成本高度計(jì)已經(jīng)獲得成功應(yīng)用,根據(jù)高度計(jì)數(shù)據(jù)的實(shí)時(shí)反饋,可以對(duì)GNSS欺騙信號(hào)進(jìn)行輔助檢測(cè)。本文所提基于氣壓高度計(jì)的 GNSS 欺騙檢測(cè)方法能夠在高程信息上有效地檢測(cè)可能存在的欺騙信號(hào),增強(qiáng)了GNSS系統(tǒng)的可靠性。

    參考文獻(xiàn)

    [1] 邊少鋒, 胡彥逢, 紀(jì)兵. GNSS欺騙防護(hù)技術(shù)國(guó)內(nèi)外研究現(xiàn)狀及展望[J]. 中國(guó)科學(xué): 信息科學(xué), 2017, 47(3): 275-287.

    BIAN S F, HU Y F, JI B. Research status and prospect of GNSS spoofing protection technology at home and abroad[J]. Scientia Sinica Informationis, 2017, 47(3): 275-287.

    [2] PSIAKI M L, HUMPHREYS T E. GNSS spoofing and detection[J]. Proceedings of the IEEE, 2016, 104(6): 1258-1270.

    [3] BROUMANDAN A, LACHAPELLE G. Spoofing detection using GNSS/INS/Odometer coupling for vehicular navigation[J]. Sensors, 2018, 18(5): 1305.

    [4] KHAN A M, IQBAL N, KHAN M F. Synthetic GNSS intermediate spoofing data generation using field recorded signals[J]. MethodsX, 2018, 18(5): 1305.

    [5] LEE L, KWON K C, AN D S, et al. GPS spoofing detection using accelerometers and performance analysis with probability of detection[J]. International Journal of Control, Automation and Systems, 2015, 13(4): 951-959.

    [6] CURRAN J T, BROUMENDAN A. On the use of low cost IMUs for GNSS spoofing detection in vehicular applications[C]∥Proc.of the International Technical Symposium on Navigation and Timing, 2017.

    [7] LO S, CHEN Y H, REID T, et al. The benefits of low cost accelerometers for GNSS anti-spoofing[C]∥Proc.of the ION Pacific PNT Meeting, 2017: 775-796.

    [8] BORIO D. An experimental evaluation of global navigation sate-llite system/inertial navigation system-verification strategies for vehicular applications[J]. IEEE Intelligent Transportation Systems Magazine, 2020, 12(3): 25-35.

    [9] LEE D K, PETIT M, MIRALLES D, et al. Analysis of raw GNSS measurements derived navigation solutions from mobile devices with inertial sensors[C]∥Proc.of the 32nd International Technical Meeting of the Satellite Division of the Institute of Navigation, 2019: 3812-3831.

    [10] LEE D K, MIRALLES D, AKOS D, et al. Detection of GNSS spoofing using NMEA messages[C]∥Proc.of the European Navigation Conference, 2020.

    [11] MIRALLES D, LEVIGNE N, AKOS D M, et al. Android raw GNSS measurements as the new anti-spoofing and anti-jamming solution[C]∥Proc.of the 31st International Technical Meeting of the Satellite Division of the Institute of Navigation, 2018: 334-344.

    [12] 王昕. 導(dǎo)航接收機(jī)中氣壓高度計(jì)的設(shè)計(jì)與實(shí)現(xiàn)[D]. 西安: 中國(guó)科學(xué)院國(guó)家授時(shí)中心, 2015.

    WANG X. Design and implementation of barometric altimeter in navigation receiver[D]. Xi’an: National Time Service Center, Chinese Academy of Sciences, 2015.

    [13] XIA H, WANG X G, QIAO Y Y, et al. Using multiple barometers to detect the floor location of smart phones with built-in barometric sensors for indoor positioning[J]. Sensors, 2015, 15(4): 7857-7877.

    [14] LEE D K, NEDELKOV F, AKOS D, et al. Barometer based GNSS spoofing detection[C]∥Proc.of the 33rd International Technical Meeting of the Satellite Division of the Institute of Navigation, 2020: 3268-3282.

    [15] ENGE P K. The global positioning system: signals, measurements and performance[J]. International Journal of Wireless Information Networks, 1994, 1(2): 83-105.

    [16] 王文益, 龔婧, 王金銘. 基于SCB方差的GNSS欺騙式干擾檢測(cè)算法[J]. 系統(tǒng)工程與電子技術(shù), 2021, 43(8): 2254-2262.

    WANG W Y, GONG J, WANG J M. GNSS spoofing inter-ference detection based on variance of SCB[J]. Systems Engineering and Electronics, 2021, 43(8): 2254-2262.

    [17] 胡遠(yuǎn)遷, 張全, 牛小驥. 基于一階高斯-馬爾可夫模型的IMU零偏相關(guān)時(shí)間對(duì)GNSS/INS組合導(dǎo)航結(jié)果的影響[J]. 科學(xué)技術(shù)與工程, 2018, 18(34): 236-241.

    HU Y Q, ZHANG Q, NIU X J. Influence of IMU zero partial correlation time based on first-order Gauss-Markov model on GNSS/INS integrated navigation results[J].Science Technology and Engineering, 2018, 18(34): 236-241.

    [18] RAVEH D E. CFD-based models of aerodynamic gust response[J]. Journal of Aircraft, 2007, 44(3): 888-897.

    [19] NASSAR S, EL-SHEIMY N. A combined algorithm of improving INS error modeling and sensor measurements for accurate INS/GPS navigation[J]. GPS Solutions, 2006, 10(1): 29-39.

    [20] QIN H T, ZHANG H T, NIU Y J, et al. Altitude data fusion algorithm based on INS and GNSS and barometer[C]∥Proc.of the IEEE 3rd International Conference on Civil Aviation Safety and Information Technology, 2021: 85-90.

    [21] LIU X X, ZHAO X B, PANG C L. Research on SINS/GNSS/CNS/pressure altimeter federated filtering adaptive fault tole-rance method[C]∥Proc.of the China Satellite Navigation Conference, 2022: 444-454.

    [22] YE L Y, YANG Y K, MA J G, et al. A distributed formation joint network navigation and positioning algorithm[J]. Mathematics, 2022, 10(10): 1627.

    [23] DREVELLE V, BONNIFAIT P. A set-membership approach for high integrity height-aided satellite positioning[J]. GPS Solutions, 2011, 15(4): 357-368.

    [24] YUN J, PARK B. A technique for finding indoor rescuer using GNSS and barometer sensors on android smart phones in building disaster situations[C]∥Proc.of the 34th International Technical Meeting of the Satellite Division of the Institute of Navigation, 2021: 1961-1980.

    [25] WU P, FENG L, TONG H B, et al. Approximate position estimation method of weak-signal receiver of global navigation satellite systems assisted by barometric altimeter[J]. Traitement du Signal, 2022, 39(3): 945-950.

    [26] GUO H, WEI D Y, OUYANG G Z, et al. Research on barometric altimeter assisting GNSS navigation for urban environments[C]∥Proc.of the China Satellite Navigation Conference, 2018: 799-811.

    [27] SIMONETTI M, CRESPILLO O G. Robust modeling of geodetic altitude from barometric altimeter and weather data[C]∥Proc.of the 34th International Technical Meeting of the Sate-llite Division of the Institute of Navigation, 2021: 1176-1189.

    [28] 魯郁. 北斗/GPS雙模軟件接收機(jī)原理與實(shí)現(xiàn)技術(shù)[M]. 北京: 電子工業(yè)出版社, 2016.

    LU Y. The principle and realization technology of BDS/GPS dual-mode software receiver[M]. Beijing: Publishing House of Electronics Industry, 2016.

    [29] BENVENUTO D, GIOVANETTI M, VASS-ALLO L, et al. Application of the ARIMA model on the COVID-2019 epidemic dataset[J]. Data in Brief, 2020, 29: 105340.

    [30] SUN C, WAYN C J, DEMPSTER A G, et al. Moving variance-based signal quality monitoring method for spoofing detection[J]. GPS Solutions, 2018, 22(3): 83.

    [31] HAND D J, TILL R J. A simple generalisation of the area under the ROC curve for multiple class classification problems[J]. Machine Learning, 2001, 45(2): 171-186.

    [32] ZHU X F, HUA T, YANG F, et al. Global positioning system spoofing detection based on support vector machines[J]. IET Radar, Sonar Navigation, 2022, 16(2): 224-237.

    作者簡(jiǎn)介

    徐奕禹(1998—),男,碩士研究生,主要研究方向?yàn)閷?dǎo)航安全、衛(wèi)星導(dǎo)航欺騙干擾檢測(cè)。

    陳長(zhǎng)風(fēng)(1991—),男,工程師,碩士,主要研究方向?yàn)殡娮訉?duì)抗技術(shù)、導(dǎo)航技術(shù)。

    袁雪林(1979—),男,副教授,博士,主要研究方向?yàn)閷?dǎo)航安全、UWB精密測(cè)量、室內(nèi)定位。

    陳正坤(1991—),男,博士研究生,主要研究方向?yàn)镚NSS多波束智能感知、GNSS/INS深組合導(dǎo)航、智能可信導(dǎo)航。

    周志?。?998—),男,碩士研究生,主要研究方向?yàn)樾滦蛯?dǎo)航接收機(jī)設(shè)計(jì)。

    猜你喜歡
    信號(hào)檢測(cè)模型
    一半模型
    “不等式”檢測(cè)題
    “一元一次不等式”檢測(cè)題
    “一元一次不等式組”檢測(cè)題
    信號(hào)
    鴨綠江(2021年35期)2021-04-19 12:24:18
    重要模型『一線三等角』
    完形填空二則
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    基于FPGA的多功能信號(hào)發(fā)生器的設(shè)計(jì)
    電子制作(2018年11期)2018-08-04 03:25:42
    3D打印中的模型分割與打包
    国精品久久久久久国模美| 亚洲aⅴ乱码一区二区在线播放| 一区二区三区免费毛片| 中文欧美无线码| 国产极品天堂在线| 自拍偷自拍亚洲精品老妇| 观看av在线不卡| 久久99热这里只有精品18| 高清av免费在线| 国产精品爽爽va在线观看网站| 久久人人爽av亚洲精品天堂 | 亚洲欧美日韩另类电影网站 | 看十八女毛片水多多多| 国产精品福利在线免费观看| 国产黄片视频在线免费观看| 国产精品伦人一区二区| 欧美人与善性xxx| 国产黄色免费在线视频| 国产男人的电影天堂91| 简卡轻食公司| 97超碰精品成人国产| 成人毛片a级毛片在线播放| 中国三级夫妇交换| 97热精品久久久久久| 亚洲人成网站在线播| 国产成人freesex在线| 在线观看免费日韩欧美大片 | 亚洲av福利一区| 亚洲精品国产av成人精品| 日日啪夜夜爽| 高清午夜精品一区二区三区| 国产精品人妻久久久影院| 亚洲精品国产av成人精品| 新久久久久国产一级毛片| 久久女婷五月综合色啪小说| 亚洲欧洲国产日韩| 欧美成人a在线观看| 狠狠精品人妻久久久久久综合| 一个人看视频在线观看www免费| 蜜臀久久99精品久久宅男| 91久久精品国产一区二区三区| 国产美女午夜福利| 人人妻人人爽人人添夜夜欢视频 | 亚洲av成人精品一二三区| 亚洲色图综合在线观看| 欧美极品一区二区三区四区| www.av在线官网国产| 免费观看a级毛片全部| 国产综合精华液| 成年免费大片在线观看| 成人高潮视频无遮挡免费网站| 精品人妻视频免费看| 偷拍熟女少妇极品色| 三级国产精品片| 18+在线观看网站| 亚洲精品乱码久久久久久按摩| 国产精品三级大全| 韩国高清视频一区二区三区| 久久久久国产网址| 26uuu在线亚洲综合色| 国产av国产精品国产| 久久久久久久大尺度免费视频| 一个人看视频在线观看www免费| 亚洲天堂av无毛| 最近中文字幕高清免费大全6| 一级av片app| 国产成人freesex在线| 精品国产乱码久久久久久小说| 久久 成人 亚洲| 国产爽快片一区二区三区| 伦理电影大哥的女人| 日韩不卡一区二区三区视频在线| 能在线免费看毛片的网站| 中文字幕人妻熟人妻熟丝袜美| 插逼视频在线观看| 精品午夜福利在线看| 成人无遮挡网站| 在线观看人妻少妇| 国内少妇人妻偷人精品xxx网站| 日本猛色少妇xxxxx猛交久久| 黄色视频在线播放观看不卡| 99re6热这里在线精品视频| av免费观看日本| 久久青草综合色| 国产爽快片一区二区三区| .国产精品久久| 少妇人妻久久综合中文| 大香蕉97超碰在线| 18+在线观看网站| av播播在线观看一区| 日韩成人av中文字幕在线观看| 晚上一个人看的免费电影| 黄色配什么色好看| 亚洲自偷自拍三级| 国产乱来视频区| 亚洲,欧美,日韩| 少妇被粗大猛烈的视频| 亚洲av免费高清在线观看| 美女中出高潮动态图| 中文字幕av成人在线电影| 国产亚洲最大av| 美女高潮的动态| 日日摸夜夜添夜夜添av毛片| 亚洲国产av新网站| 国产精品一区www在线观看| 久久久欧美国产精品| 国产精品成人在线| 啦啦啦中文免费视频观看日本| 亚洲精品乱码久久久v下载方式| 欧美日韩综合久久久久久| 亚洲自偷自拍三级| 韩国高清视频一区二区三区| 欧美区成人在线视频| 777米奇影视久久| 国产精品99久久99久久久不卡 | 国产伦精品一区二区三区四那| 国产精品国产三级专区第一集| 少妇裸体淫交视频免费看高清| 特大巨黑吊av在线直播| 日本黄大片高清| 欧美另类一区| 欧美成人午夜免费资源| 在线观看av片永久免费下载| 啦啦啦视频在线资源免费观看| 久久精品国产鲁丝片午夜精品| 亚洲国产高清在线一区二区三| 大片电影免费在线观看免费| 国产精品久久久久久精品古装| 亚洲av二区三区四区| 免费不卡的大黄色大毛片视频在线观看| 我的老师免费观看完整版| 亚洲av福利一区| 免费在线观看成人毛片| av天堂中文字幕网| 九色成人免费人妻av| 天堂俺去俺来也www色官网| 亚洲欧美日韩东京热| 九九久久精品国产亚洲av麻豆| 国内精品宾馆在线| 免费少妇av软件| 天天躁夜夜躁狠狠久久av| 观看美女的网站| 免费黄网站久久成人精品| 熟女电影av网| 偷拍熟女少妇极品色| 高清毛片免费看| 制服丝袜香蕉在线| 久久97久久精品| 丝袜喷水一区| 亚洲精品乱码久久久久久按摩| 久久青草综合色| 观看免费一级毛片| 亚洲av中文av极速乱| 亚洲,一卡二卡三卡| 国产大屁股一区二区在线视频| 久久热精品热| 丰满少妇做爰视频| 国产伦理片在线播放av一区| 一级毛片 在线播放| 熟女人妻精品中文字幕| 国产精品一区二区在线观看99| 自拍欧美九色日韩亚洲蝌蚪91 | 国产伦在线观看视频一区| 草草在线视频免费看| 久久午夜福利片| 日韩一区二区视频免费看| 亚洲真实伦在线观看| 91精品伊人久久大香线蕉| 成人美女网站在线观看视频| 观看av在线不卡| 久久国内精品自在自线图片| 久久久欧美国产精品| 国内揄拍国产精品人妻在线| 国产中年淑女户外野战色| 亚洲欧美日韩东京热| 爱豆传媒免费全集在线观看| 性色av一级| 大陆偷拍与自拍| 九色成人免费人妻av| 嘟嘟电影网在线观看| 一级毛片黄色毛片免费观看视频| 免费观看性生交大片5| tube8黄色片| 久久久午夜欧美精品| 国产成人a∨麻豆精品| 联通29元200g的流量卡| 精品久久久久久电影网| 欧美xxxx性猛交bbbb| 色综合色国产| 高清av免费在线| 欧美 日韩 精品 国产| 大话2 男鬼变身卡| 久久久久久久精品精品| 国产av一区二区精品久久 | 成年免费大片在线观看| 边亲边吃奶的免费视频| av国产久精品久网站免费入址| a级一级毛片免费在线观看| 少妇被粗大猛烈的视频| 国产乱来视频区| 国产熟女欧美一区二区| 18禁在线播放成人免费| 一级毛片aaaaaa免费看小| 少妇被粗大猛烈的视频| 麻豆精品久久久久久蜜桃| 精品国产一区二区三区久久久樱花 | 亚洲,一卡二卡三卡| 国产色婷婷99| 一区在线观看完整版| 日韩亚洲欧美综合| 欧美激情极品国产一区二区三区 | 国产探花极品一区二区| 精品久久久久久久久av| 在线观看人妻少妇| 国内揄拍国产精品人妻在线| 国产69精品久久久久777片| 全区人妻精品视频| 久久 成人 亚洲| 日韩欧美一区视频在线观看 | 六月丁香七月| 亚洲在久久综合| 中文精品一卡2卡3卡4更新| 国产精品伦人一区二区| 婷婷色麻豆天堂久久| 性色avwww在线观看| 日日摸夜夜添夜夜添av毛片| 免费不卡的大黄色大毛片视频在线观看| 国产无遮挡羞羞视频在线观看| 精品熟女少妇av免费看| 六月丁香七月| 久久久久国产网址| 偷拍熟女少妇极品色| 免费少妇av软件| 日韩伦理黄色片| 久久精品夜色国产| 午夜福利网站1000一区二区三区| 能在线免费看毛片的网站| 亚洲无线观看免费| 精品一区在线观看国产| 国产日韩欧美在线精品| 久久久久久伊人网av| 国产亚洲一区二区精品| 国产爱豆传媒在线观看| 国产精品嫩草影院av在线观看| 国产在线视频一区二区| 欧美性感艳星| 联通29元200g的流量卡| 亚州av有码| 国产精品蜜桃在线观看| av在线观看视频网站免费| 久久国内精品自在自线图片| 黄色一级大片看看| a级毛片免费高清观看在线播放| 乱码一卡2卡4卡精品| 国产真实伦视频高清在线观看| 久久久色成人| 日本免费在线观看一区| 六月丁香七月| 五月开心婷婷网| 国产高清国产精品国产三级 | 久久影院123| 亚洲第一av免费看| 又大又黄又爽视频免费| 国产精品蜜桃在线观看| 国产真实伦视频高清在线观看| 国产人妻一区二区三区在| 波野结衣二区三区在线| 少妇猛男粗大的猛烈进出视频| 91精品一卡2卡3卡4卡| 亚洲精品乱久久久久久| 午夜精品国产一区二区电影| 国产精品无大码| 日韩 亚洲 欧美在线| 熟女电影av网| 日韩欧美精品免费久久| 人人妻人人澡人人爽人人夜夜| 一级毛片久久久久久久久女| 18+在线观看网站| 又粗又硬又长又爽又黄的视频| 国产深夜福利视频在线观看| 男女下面进入的视频免费午夜| 97精品久久久久久久久久精品| 亚洲电影在线观看av| 国产成人免费无遮挡视频| 中国美白少妇内射xxxbb| 国产男女内射视频| 国产淫语在线视频| 久久久久视频综合| 丝袜喷水一区| 婷婷色综合www| 黄色配什么色好看| 久久青草综合色| 在线播放无遮挡| 日韩av免费高清视频| 在线观看人妻少妇| 色吧在线观看| 老师上课跳d突然被开到最大视频| 美女主播在线视频| 香蕉精品网在线| 久久99精品国语久久久| 免费在线观看成人毛片| 22中文网久久字幕| 中文乱码字字幕精品一区二区三区| 欧美区成人在线视频| 国产免费视频播放在线视频| 亚洲无线观看免费| 国产精品爽爽va在线观看网站| 联通29元200g的流量卡| 少妇的逼好多水| 99热这里只有是精品50| 日韩中文字幕视频在线看片 | 免费av中文字幕在线| 国产成人a区在线观看| 久久久成人免费电影| 纯流量卡能插随身wifi吗| 国产精品成人在线| 亚洲成人一二三区av| 精品人妻偷拍中文字幕| av女优亚洲男人天堂| 日本免费在线观看一区| 高清欧美精品videossex| 亚洲成色77777| 日日摸夜夜添夜夜爱| 亚洲,欧美,日韩| 久久久久人妻精品一区果冻| 日韩在线高清观看一区二区三区| 舔av片在线| 午夜视频国产福利| 成人特级av手机在线观看| 夜夜爽夜夜爽视频| 午夜激情久久久久久久| 美女中出高潮动态图| 国产一区亚洲一区在线观看| 黄色欧美视频在线观看| 亚洲内射少妇av| 日本欧美视频一区| 国产男人的电影天堂91| 免费人成在线观看视频色| 亚洲精品色激情综合| 久久久a久久爽久久v久久| 菩萨蛮人人尽说江南好唐韦庄| 国产欧美日韩一区二区三区在线 | freevideosex欧美| 国产乱来视频区| 这个男人来自地球电影免费观看 | 五月天丁香电影| 天堂中文最新版在线下载| 日本一二三区视频观看| 在线看a的网站| 国产精品熟女久久久久浪| 熟女电影av网| 在现免费观看毛片| 我的老师免费观看完整版| 日本欧美国产在线视频| 欧美精品一区二区大全| 久久精品夜色国产| av在线观看视频网站免费| 久久人人爽人人片av| 男人舔奶头视频| 日本一二三区视频观看| 熟女av电影| 爱豆传媒免费全集在线观看| 久热久热在线精品观看| 91午夜精品亚洲一区二区三区| 男人和女人高潮做爰伦理| 亚洲欧美精品自产自拍| 在线免费十八禁| 中文天堂在线官网| 亚洲精品成人av观看孕妇| 日韩av不卡免费在线播放| 亚洲欧美中文字幕日韩二区| 男的添女的下面高潮视频| av女优亚洲男人天堂| 交换朋友夫妻互换小说| 少妇猛男粗大的猛烈进出视频| 亚洲人成网站在线观看播放| 成人国产麻豆网| 高清av免费在线| 午夜精品国产一区二区电影| 亚洲精品,欧美精品| 插阴视频在线观看视频| 日日摸夜夜添夜夜添av毛片| 日韩大片免费观看网站| 国产无遮挡羞羞视频在线观看| 日本免费在线观看一区| av在线播放精品| 日韩欧美一区视频在线观看 | 国产永久视频网站| 中国国产av一级| 中文字幕制服av| 三级国产精品欧美在线观看| 嫩草影院新地址| 久久毛片免费看一区二区三区| 日韩中字成人| av不卡在线播放| 人妻制服诱惑在线中文字幕| 一边亲一边摸免费视频| 亚洲人成网站在线播| 国产精品久久久久久久电影| 丰满乱子伦码专区| 男人添女人高潮全过程视频| 亚洲美女视频黄频| 一级av片app| 日韩av不卡免费在线播放| 18禁裸乳无遮挡免费网站照片| 99re6热这里在线精品视频| 激情五月婷婷亚洲| 亚洲第一av免费看| 久久99热这里只有精品18| 日本黄色片子视频| 三级经典国产精品| 黄色怎么调成土黄色| 久久亚洲国产成人精品v| 一本色道久久久久久精品综合| 久久久久久久国产电影| 亚洲一区二区三区欧美精品| 亚洲综合精品二区| 日韩,欧美,国产一区二区三区| 国产亚洲最大av| 精品一区二区三卡| 亚洲伊人久久精品综合| 晚上一个人看的免费电影| 国产精品福利在线免费观看| 亚洲国产色片| 夜夜看夜夜爽夜夜摸| 亚洲,一卡二卡三卡| 中文资源天堂在线| 有码 亚洲区| 少妇精品久久久久久久| 搡女人真爽免费视频火全软件| 热re99久久精品国产66热6| 国产永久视频网站| 国产精品熟女久久久久浪| 老司机影院毛片| 91在线精品国自产拍蜜月| 日韩成人伦理影院| 成人特级av手机在线观看| 亚洲三级黄色毛片| 亚洲精品色激情综合| 直男gayav资源| 3wmmmm亚洲av在线观看| 国产亚洲一区二区精品| 久久久亚洲精品成人影院| 亚洲精品日本国产第一区| 欧美三级亚洲精品| 亚洲精品第二区| 国产黄频视频在线观看| 久久精品国产a三级三级三级| 国产欧美另类精品又又久久亚洲欧美| 久久毛片免费看一区二区三区| 人人妻人人添人人爽欧美一区卜 | 日韩强制内射视频| 精品酒店卫生间| 18禁裸乳无遮挡免费网站照片| 成人毛片60女人毛片免费| 亚洲欧美日韩另类电影网站 | 日韩一区二区三区影片| 欧美bdsm另类| 久久国产精品男人的天堂亚洲 | 精品人妻偷拍中文字幕| 性色av一级| 一区二区av电影网| 久久精品熟女亚洲av麻豆精品| 成人漫画全彩无遮挡| 黄色配什么色好看| 国产精品熟女久久久久浪| 久久久午夜欧美精品| 亚洲av不卡在线观看| 成人国产av品久久久| 男人爽女人下面视频在线观看| 久热这里只有精品99| 国产一区二区三区av在线| 51国产日韩欧美| 亚洲国产高清在线一区二区三| 一级毛片电影观看| 国产高清三级在线| 成人影院久久| 国产白丝娇喘喷水9色精品| 丰满少妇做爰视频| 亚洲精品亚洲一区二区| 免费黄网站久久成人精品| 亚洲av福利一区| 色5月婷婷丁香| 王馨瑶露胸无遮挡在线观看| 日韩一本色道免费dvd| 久热久热在线精品观看| 亚洲自偷自拍三级| 伦精品一区二区三区| 欧美日韩国产mv在线观看视频 | 女的被弄到高潮叫床怎么办| 国产欧美另类精品又又久久亚洲欧美| 日韩成人伦理影院| 亚洲国产最新在线播放| 亚洲av国产av综合av卡| 老熟女久久久| 丰满乱子伦码专区| 一个人看的www免费观看视频| 亚洲国产色片| 最黄视频免费看| 又黄又爽又刺激的免费视频.| 激情 狠狠 欧美| 国产精品一区二区在线观看99| 久久久久精品性色| 日韩伦理黄色片| 中文字幕久久专区| 久久青草综合色| 永久网站在线| 婷婷色av中文字幕| 国产精品人妻久久久影院| 91狼人影院| 性色avwww在线观看| 这个男人来自地球电影免费观看 | 亚洲av不卡在线观看| 国产一区亚洲一区在线观看| 99国产精品免费福利视频| 亚洲精品国产av成人精品| freevideosex欧美| 新久久久久国产一级毛片| 日韩强制内射视频| 大陆偷拍与自拍| 秋霞在线观看毛片| 日日撸夜夜添| 亚洲精品一二三| 午夜老司机福利剧场| 这个男人来自地球电影免费观看 | 18禁动态无遮挡网站| 国产男人的电影天堂91| 亚洲国产日韩一区二区| 国产老妇伦熟女老妇高清| 精华霜和精华液先用哪个| 黄色欧美视频在线观看| 国产av国产精品国产| 亚洲精品aⅴ在线观看| 女的被弄到高潮叫床怎么办| 少妇的逼好多水| 少妇高潮的动态图| 水蜜桃什么品种好| 亚洲怡红院男人天堂| 精品一品国产午夜福利视频| 少妇猛男粗大的猛烈进出视频| 国产永久视频网站| av线在线观看网站| 在线观看国产h片| 久久久久久久久久人人人人人人| 久久久久国产精品人妻一区二区| 久久99热这里只有精品18| 国产又色又爽无遮挡免| 亚洲国产最新在线播放| 久久99蜜桃精品久久| 国产精品麻豆人妻色哟哟久久| 美女xxoo啪啪120秒动态图| 人体艺术视频欧美日本| 热re99久久精品国产66热6| 久久久久视频综合| 日本黄色日本黄色录像| av女优亚洲男人天堂| 视频区图区小说| 午夜视频国产福利| 欧美日韩国产mv在线观看视频 | 极品少妇高潮喷水抽搐| 免费大片18禁| 天堂中文最新版在线下载| 国产精品蜜桃在线观看| 国产精品欧美亚洲77777| 夫妻性生交免费视频一级片| 久久久久性生活片| 黑人高潮一二区| 日韩,欧美,国产一区二区三区| 国产精品久久久久久精品电影小说 | 黄色视频在线播放观看不卡| 九九久久精品国产亚洲av麻豆| 在线播放无遮挡| 久久国内精品自在自线图片| 最黄视频免费看| 少妇的逼水好多| 亚洲国产最新在线播放| 精品人妻视频免费看| 联通29元200g的流量卡| 嫩草影院入口| 国产精品福利在线免费观看| 欧美97在线视频| 少妇猛男粗大的猛烈进出视频| 男人舔奶头视频| 国产色爽女视频免费观看| 精品少妇久久久久久888优播| 中文天堂在线官网| 精品人妻一区二区三区麻豆| 视频中文字幕在线观看| 午夜福利网站1000一区二区三区| 亚洲精品乱码久久久v下载方式| 中国国产av一级| 男女免费视频国产| 成年av动漫网址| 在线观看国产h片| 男女免费视频国产| 国产深夜福利视频在线观看| 欧美性感艳星| 亚洲精品aⅴ在线观看| 精品亚洲乱码少妇综合久久| 免费大片18禁| 欧美精品人与动牲交sv欧美| 极品少妇高潮喷水抽搐| 久久99热这里只频精品6学生| 日韩一本色道免费dvd| 观看免费一级毛片| 亚洲aⅴ乱码一区二区在线播放| 日日啪夜夜撸| 成人美女网站在线观看视频| av卡一久久| 国产精品无大码| 久久久精品免费免费高清| 国模一区二区三区四区视频| 一区在线观看完整版| 久久毛片免费看一区二区三区| 大香蕉97超碰在线|