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

    多星座組合RAIM技術(shù)在列車定位中的應(yīng)用

    2022-06-10 03:49:10鐵凱博米根鎖王彥快
    導(dǎo)航定位學(xué)報(bào) 2022年3期
    關(guān)鍵詞:偽距可用性星座

    鐵凱博,米根鎖,王彥快

    (蘭州交通大學(xué), 自動化與電氣工程學(xué)院, 蘭州 730070)

    0 引言

    列車定位精度直接影響列車的運(yùn)行安全,提高列車定位的精度和可靠性已經(jīng)成為不可忽視的重要目標(biāo)[1]。全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system, GNSS)在鐵路領(lǐng)域已廣泛應(yīng)用,但在列車實(shí)際運(yùn)行過程中,衛(wèi)星出現(xiàn)故障,會導(dǎo)致定位結(jié)果不準(zhǔn)確[2],這就要求導(dǎo)航系統(tǒng)能對所用衛(wèi)星進(jìn)行實(shí)時(shí)檢測,剔除出現(xiàn)故障的衛(wèi)星。因此,有學(xué)者提出將接收機(jī)自主完好性監(jiān)測(receiver autonomous integrity monitoring,RAIM)技術(shù)[3]應(yīng)用于列車定位。文獻(xiàn)[4]提出了RAIM技術(shù)在列車定位中的可用性預(yù)測方法,并對該技術(shù)在鐵路領(lǐng)域內(nèi)的可用性做出了判斷;文獻(xiàn)[5]將RAIM技術(shù)應(yīng)用于北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDou navigation satellite system,BDS)的列車定位中,提出沿股道水平保護(hù)級(along track protect level,ATPL)算法,提高了RAIM技術(shù)在列車定位中的可用性;文獻(xiàn)[6]分析了BDS完好性問題在列車定位中的風(fēng)險(xiǎn),并設(shè)計(jì)了BDS RAIM可用性預(yù)測系統(tǒng),保障了RAIM預(yù)測結(jié)果的可靠性;文獻(xiàn)[7]在故障衛(wèi)星檢測環(huán)節(jié)中引入加權(quán)因子,降低檢測過程中的漏警率,提高了列車定位精度。文獻(xiàn)[8]解決了GNSS性能指標(biāo)與鐵路可用性指標(biāo)的映射問題,為基于GNSS的 RAIM 算法在鐵路領(lǐng)域中的應(yīng)用提供了參考。通過分析文獻(xiàn)可知:目前主要針對單星座RAIM 在列車定位中的應(yīng)用進(jìn)行了研究,單星座RAIM應(yīng)于列車定位時(shí),常會出現(xiàn)因可見星不足或幾何分布較差而降低算法可用性的情況,進(jìn)而影響算法的故障檢測率和故障識別率。由于可見星數(shù)目的限制,單星座 RAIM 在某些高海拔鐵路定位精度并不高[9],甚至不能使用,RAIM算法在鐵路領(lǐng)域內(nèi)的研究還處于起步階段。

    鑒于此,本文提出了將多星座RAIM技術(shù)應(yīng)用于列車定位,通過增加可見星數(shù)目提高RAIM算法的可用性。在構(gòu)造多星座系統(tǒng)觀測方程的基礎(chǔ)上,采用奇偶矢量法對導(dǎo)航衛(wèi)星進(jìn)行故障檢測。通過青藏鐵路列車定位數(shù)據(jù)對全球定位系統(tǒng)(global positioning system, GPS)、伽利略衛(wèi)星導(dǎo)航系統(tǒng)(Galileo satellite navigation system, Galileo)和 BDS多星座組合RAIM和GPS單星座RAIM進(jìn)行對比仿真,分析多星座RAIM在列車定位中的檢測性能。

    1 多星座系統(tǒng)列車定位原理

    1.1 時(shí)空基準(zhǔn)統(tǒng)一

    為構(gòu)成基于列車定位的多星座導(dǎo)航系統(tǒng),需將不同星座的時(shí)間系統(tǒng)進(jìn)行統(tǒng)一,坐標(biāo)系統(tǒng)進(jìn)行轉(zhuǎn)換,采用相同的時(shí)間和坐標(biāo)標(biāo)準(zhǔn),實(shí)現(xiàn)多星座系統(tǒng)的數(shù)據(jù)融合[10-11]。本文將多星座時(shí)間、坐標(biāo)系統(tǒng)統(tǒng)一轉(zhuǎn)換為以GPS為基準(zhǔn)。轉(zhuǎn)換過程如下:

    1)將Galileo系統(tǒng)時(shí)(Galileo system time, GST)轉(zhuǎn)化為GPS時(shí)(GPS time, GPST)的計(jì)算公式為

    式中:TGPST為 GPS時(shí);TGST為 Galileo系統(tǒng)時(shí);TGGTO為GPST與GST間的偏差;A0G、A1G、t0G和SWNOG分別為時(shí)間偏差的常數(shù)項(xiàng)、變化率、參考時(shí)間和參考星期個(gè)數(shù);TTOW為星期內(nèi)的時(shí)間;T是常數(shù),為604 800;SWN為Galileo的星期個(gè)數(shù);mod為求余函數(shù)。

    2)將北斗時(shí)(BDS time, BDT)轉(zhuǎn)化為GPST的計(jì)算公式為

    式中,TBDT為北斗系統(tǒng)時(shí)。

    3)將格洛納斯衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system, GLONASS)時(shí)(GLONASS time,GLONASST)轉(zhuǎn)化為GPST的計(jì)算公式為

    式中:TGLONASST為GLONASS時(shí);n為協(xié)調(diào)世界時(shí)(coordinated universal time, UTC)與國際原子時(shí)(international atomic time, TAI)之間不斷調(diào)整的參數(shù),由國際地球自轉(zhuǎn)服務(wù)組織(International Earth Rotation Service, IERS)提供。

    GPS采用1984世界大地坐標(biāo)系(world geodetic coordinate system 1984, WGS84);Galileo采用伽利略地球參考框架(Galileo terrestrial reference frame,GTRF);BDS采用北斗坐標(biāo)系(BeiDou coordinate system,BDCS);GLONASS采用俄羅斯大地坐標(biāo)框架(Parametry Zelmy1990, PZ90)。由于列車偽距定位的精度為米級,BDCS坐標(biāo)和GTRF坐標(biāo)轉(zhuǎn)換為WGS84坐標(biāo)產(chǎn)生的誤差可以忽略不計(jì)[12],只需將PZ90坐標(biāo)轉(zhuǎn)換為WGS84坐標(biāo)即可,采用布爾薩-沃爾夫(Bursa-Wolf)模型[13]進(jìn)行PZ90和WGS84間的坐標(biāo)轉(zhuǎn)換,Bursa-Wolf模型亦稱為黑爾默特(Helmert)7參數(shù)模型,其計(jì)算公式為

    1.2 列車定位原理

    在列車運(yùn)行過程中,可見星連續(xù)不斷地向車載接收機(jī)發(fā)送星歷參數(shù)和時(shí)間信息,車載計(jì)算機(jī)計(jì)算出與N顆可見星之間的偽距(ρ1,ρ2,ρ3,…,ρN),通過聯(lián)立方程進(jìn)行列車位置的定位解算,即

    式中:ρi(i=1,2,...,N)為第i號衛(wèi)星到車載接收機(jī)的偽距,N為可見星數(shù)目;(xi,yi,zi)為第i號衛(wèi)星的空間位置;(xu,yu,zu)為列車位置;c為光速;tu為車載接收機(jī)與衛(wèi)星間的鐘差,其定位原理如圖1所示。

    圖1 列車偽距定位原理

    若參與解算的可見星發(fā)生故障,會造成ρi的測量值出現(xiàn)偏差,導(dǎo)致列車定位解算結(jié)果不準(zhǔn)確。因此,需在列車定位解算之前,采用RAIM技術(shù)對衛(wèi)星的健康狀況進(jìn)行判斷,剔除出現(xiàn)故障的衛(wèi)星。

    2 基于列車定位的多星座RAIM算法

    在列車定位中,RAIM算法采用冗余技術(shù)對車載接收機(jī)接收到的導(dǎo)航信息進(jìn)行一致性校驗(yàn),其本質(zhì)是通過增加可見星數(shù)目來增加信息冗余,實(shí)現(xiàn)對衛(wèi)星的故障檢測[14]。當(dāng)列車導(dǎo)航系統(tǒng)至少監(jiān)測到4顆衛(wèi)星時(shí),可進(jìn)行用戶位置的定位解算;至少監(jiān)測到5顆衛(wèi)星時(shí),可實(shí)現(xiàn)故障衛(wèi)星的檢測;至少監(jiān)測到6顆衛(wèi)星時(shí),可對故障衛(wèi)星進(jìn)行識別。因此,可見星數(shù)目對 RAIM 算法的可用性有很大影響,通過采用多星座系統(tǒng)增加可見星數(shù)目進(jìn)而改善幾何分布情況對 RAIM 算法的可用性會有相應(yīng)提升[15]。

    在鐵路沿線中,海拔較高時(shí),會使部分衛(wèi)星的檢測效果變差[16],影響算法的可用性。多星座組合系統(tǒng)通過增加線路中可見衛(wèi)星數(shù)目增加冗余觀測信息,進(jìn)而提高 RAIM 算法在列車定位中的可用性。當(dāng)車載接收機(jī)監(jiān)測到的可見衛(wèi)星數(shù)N>n時(shí)(n=4+m,m為星座數(shù)),可進(jìn)行故障衛(wèi)星的檢測與識別。對于單星座系統(tǒng)至少監(jiān)測到n=4+1=5顆衛(wèi)星;對于雙星座系統(tǒng)至少監(jiān)測到n=4+2=6顆衛(wèi)星;對于三星座系統(tǒng)至少監(jiān)測到n=4+3=7顆衛(wèi)星;以此類推。圖2為多星座 RAIM算法的程序流程圖。圖2中n為故障檢測所需的最少可見星數(shù)。

    圖2 多星座RAIM算法程序流程

    2.1 構(gòu)造多星座系統(tǒng)觀測方程

    假設(shè)在某一時(shí)刻共有N顆可見星,對式(5)進(jìn)行泰勒展開,并考慮觀測誤差,線性化后的多星座系統(tǒng)觀測方程為

    式中:y?RN×1為觀測偽距與計(jì)算偽距的差值;x?Rl×1為列車的三維位置信息和m個(gè)星座對應(yīng)的時(shí)鐘誤差參數(shù),l=m+3;ε?RN×1為觀測偽距噪聲矢量,服從均值為 0,方差為的正態(tài)分布。H?RN?l為N ?l維的觀測矩陣,其表達(dá)式為

    式中:A為列車與衛(wèi)星之間的方向余弦值;B為對應(yīng)的各星座車載接收機(jī)鐘差;Aij(i=1,…,m,j=1,…,3)為對應(yīng)第i個(gè)星座觀測矩陣Ai的第j列元素,Ai可表示為

    2.2 基于奇偶矢量法的故障衛(wèi)星檢測識別

    假設(shè)偽距誤差ε=O,則系統(tǒng)的線性化觀測方程為

    將系數(shù)矩陣H分解成

    式中:Q是N?N階的正交矩陣;R為N?l階的上三角矩陣;O為零矩陣。

    將式(10)代入式(9)并左乘QT,可得

    式中,QT和R可分別表示為

    式中:Qx為QT的前l(fā)行,Q為QT的后N-l行;Rx為R的前l(fā)行,O是R的后N-l行,為零矩陣。

    將式(12)代入式(11)中可得

    可得最小二乘解

    若考慮到觀測誤差ε,即y=Hx+ε ,則有

    式中:p為奇偶空間上的奇偶矢量;Q為奇偶空間矩陣;奇偶空間矢量p是由ε在奇偶空間矩陣Q上投影所得,可反映故障衛(wèi)星偏差量等信息。因此,可基于奇偶矢量p構(gòu)造檢測統(tǒng)計(jì)量,進(jìn)行故障衛(wèi)星的檢測識別。

    衛(wèi)星偽距殘差向量平方和為v=pTp。根據(jù)統(tǒng)計(jì)分布理論,當(dāng)偽距噪聲矢量ε中的各個(gè)分量相互獨(dú)立,且服從均值為 0,方差為的正態(tài)分布時(shí),v/服從自由度為N-l的中心X2分布;均值不為0時(shí),v/服從自由度為N?l的非中心χ2分布,非中心參數(shù)為λ,即

    構(gòu)造檢驗(yàn)統(tǒng)計(jì)量ri為

    式中:ri為第i顆衛(wèi)星的檢驗(yàn)統(tǒng)計(jì)量;Qi為第i顆衛(wèi)星的奇偶空間矩陣。

    已知誤警率PFA,則有

    式中:P(ri>Tr)為檢測統(tǒng)計(jì)量大于檢測門限的概率,為標(biāo)準(zhǔn)正態(tài)分布的概率密度函數(shù)。

    根據(jù)式(18)可求出檢測門限Tr,將每個(gè)檢測統(tǒng)計(jì)量ri與Tr進(jìn)行比較,若ri>Tr,說明第i顆衛(wèi)星出現(xiàn)故障。

    2.3 多星座RAIM算法可用性判斷

    為保障列車定位結(jié)果的可靠性,在可見星數(shù)目達(dá)到要求后,還需對沿線可見星的空間幾何分布進(jìn)行預(yù)測,對RAIM算法的可用性進(jìn)行判斷。

    本文采用水平保護(hù)級(horizontal protect level,HPL)THPL來判斷RAIM可用性。根據(jù)車載接收機(jī)接收到的列車定位數(shù)據(jù)計(jì)算THPL值,當(dāng)THPL小于水平告警閾值(horizontal alert limit,HAL)THAL時(shí),算法可用,反之則無效,車載需重新接收衛(wèi)星信號進(jìn)行可用性判斷。

    利用最大特征斜率Lmax乘以臨界偏差Pbias計(jì)算出THPL的值。首先,計(jì)算每顆衛(wèi)星的特征斜率為

    式中:Li為第i顆衛(wèi)星的特征斜率;G1i、G2i分別為矩陣G中元素,G= (HTH)?1HT;S為矩陣ii S的對角元素,S=I? (HTH)?1HT。

    然后根據(jù)檢測門限Tr,基于漏檢概率PMD計(jì)算出臨界偏差Pbias=σ0., λ的計(jì)算方法為

    水平保護(hù)限值可由式(21)計(jì)算出,即

    式中,SHPL為計(jì)算出的水平保護(hù)限值。

    將SHPL與HAL進(jìn)行比較,若SHPL

    3 多星座RAIM算法仿真分析

    利用矩陣實(shí)驗(yàn)室(matrix laboratory, MATLAB)軟件對算法進(jìn)行仿真分析,本次實(shí)驗(yàn)選用 GPS、Galileo和BDS三星座系統(tǒng)。首先根據(jù)3個(gè)星座系統(tǒng)的衛(wèi)星及軌道信息參數(shù),對星座進(jìn)行模擬,模擬中選擇24顆GPS衛(wèi)星,24顆Galileo衛(wèi)星和14顆BDS衛(wèi)星,選擇列車用戶位置在青藏鐵路某段線路的固定位置處,用戶觀測衛(wèi)星的高度截至角為15°。

    3.1 多星座RAIM算法可用性驗(yàn)證

    首先對多星座系統(tǒng)和單星座系統(tǒng)下的可見星數(shù)量進(jìn)行仿真對比。圖3為GPS/Galileo/BDS多星座系統(tǒng)與GPS單星座系統(tǒng)在某天同時(shí)刻下可見星數(shù)量的對比情況。由圖3可知,在24 h內(nèi)的任意時(shí)刻,多星座系統(tǒng)的可見星數(shù)量均高于單星座系統(tǒng),由于多星座系統(tǒng)增加了衛(wèi)星數(shù)量,所以在相同路況下,可見星的數(shù)量也隨之增多。隨著可見星數(shù)量增多,更容易滿足 RAIM 算法對冗余觀測信息的需求,在列車運(yùn)行中的某些特殊路段,如城市、山川,多星座 RAIM 能夠更好地保障導(dǎo)航信息的連續(xù)性和完好性。

    圖3 可見星數(shù)量對比

    圖4為 GPS/Galileo/BDS多星座系統(tǒng)和 GPS單星座系統(tǒng)在 24 h內(nèi)對應(yīng)的位置精度因子(position dilution of precision, PDOP)曲線。

    圖4 位置精度因子曲線對比

    由圖4可知,多星座系統(tǒng)各時(shí)刻的位置精度因子均小于GPS單星座系統(tǒng),說明多星座系統(tǒng)能夠有效改善可見星幾何分布的情況,進(jìn)而提高算法的可用性。

    圖5表示青藏鐵路某段中密度線路在GPS/Galileo/BDS多星座系統(tǒng)和 GPS單星座系統(tǒng)下的SHPL值對比情況。

    圖5 SHPL對比

    根據(jù)鐵路領(lǐng)域完好性水平要求,中密度線路的水平告警閾值為20 m[17]。由圖5可知,GPS單星座系統(tǒng)SHPL值的波動范圍為11.2~27.0 m,多星座系統(tǒng)SHPL值的波動范圍為6.0~18.7 m,在34 min之后,GPS單星座系統(tǒng)的SHPL開始出現(xiàn)大于水平告警閾值的情況,發(fā)生虛警現(xiàn)象。多星座系統(tǒng)的SHPL值始終小于水平告警閾值,未出現(xiàn)虛警,且在同一時(shí)刻,多星座系統(tǒng)的SHPL值均小于 GPS單星座系統(tǒng),說明列車在相同運(yùn)行路況下,多星座組合RAIM算法的可用性比GPS單星座系統(tǒng)更好。

    3.2 多星座RAIM算法故障檢測識別

    本次實(shí)驗(yàn)所用數(shù)據(jù)為無故障衛(wèi)星時(shí)的列車定位數(shù)據(jù),需通過人為給衛(wèi)星注入偽距偏差來模擬衛(wèi)星出現(xiàn)故障,進(jìn)而對算法進(jìn)行驗(yàn)證。選用線路數(shù)據(jù)中幾何構(gòu)型較好的21顆無故障可見星,依次編號為1號至21號,其中:1號至8號為GPS衛(wèi)星;9號至16號為Galileo衛(wèi)星;17號至21號為BDS衛(wèi)星。設(shè)定誤警率P= 1 × 1 0?6,漏警率P= 1 × 1 0?6,F(xiàn)AMD假設(shè)各導(dǎo)航系統(tǒng)的偽距觀測噪聲水平相同,取σ0=5 m 。

    在第80—120秒的時(shí)間間隔內(nèi),給9號衛(wèi)星注入60 m的持續(xù)故障,其結(jié)果如圖6所示。

    圖6 持續(xù)故障檢測

    由圖6可知,在第80—120秒的時(shí)間間隔內(nèi),檢測統(tǒng)計(jì)量超出了預(yù)設(shè)閾值,算法檢測出了衛(wèi)星出現(xiàn)持續(xù)偽距偏差。

    在第20—80秒的時(shí)間間隔內(nèi),給第17號衛(wèi)星注入0~60 m的慢變斜坡故障,其結(jié)果如圖7所示。

    圖7 慢變斜坡故障檢測

    由圖7可知,在第32秒時(shí),檢測統(tǒng)計(jì)量超出閾值,檢測到有故障發(fā)生。隨著偽距偏差的增大,檢測統(tǒng)計(jì)量也不斷增大,算法檢測出了衛(wèi)星出現(xiàn)慢變斜坡偽距偏差。

    通過故障識別算法得到兩種故障情況下各衛(wèi)星的檢測統(tǒng)計(jì)量如表1所示,檢測統(tǒng)計(jì)量最大的衛(wèi)星即為故障星。

    表1 各星座系統(tǒng)衛(wèi)星檢測統(tǒng)計(jì)量的值

    由表1可知,在持續(xù)故障情況下,9號衛(wèi)星的檢測統(tǒng)計(jì)量最大,為算法識別出的故障星。在慢變斜坡故障情況下,17號衛(wèi)星的檢測統(tǒng)計(jì)量最大,為算法識別出的故障星。算法有效識別出了在持續(xù)故障和慢變斜坡故障兩種情況下預(yù)設(shè)的故障星。

    3.3 多星座RAIM算法檢測性能分析

    為檢驗(yàn)多星座 RAIM 算法的檢測性能,利用3.2節(jié)列車定位數(shù)據(jù)對算法的故障檢測率和故障識別率進(jìn)行仿真分析。

    給某顆衛(wèi)星注入 0~200 m的偽距偏差,采樣時(shí)間為24 h,采樣間隔1 min,仿真結(jié)果如圖8和圖9所示。

    圖8 故障檢測率對比

    圖9 故障識別率對比

    由仿真結(jié)果可知,當(dāng)偽距偏差達(dá)到100 m時(shí),多星座系統(tǒng)的故障檢測率和故障識別率均可達(dá) 100%,而對于GPS單星座系統(tǒng)當(dāng)偽距偏差達(dá)到136 m時(shí),故障檢測率才達(dá)到100%,偽距偏差達(dá)到160 m時(shí),故障識別率才達(dá)到94.1%。表明多星座組合RAIM的故障檢測性能優(yōu)于GPS單星座RAIM。

    4 結(jié)束語

    本文利用多星座組合RAIM技術(shù),實(shí)時(shí)檢測接收機(jī)自身以及衛(wèi)星故障,進(jìn)而提高列車運(yùn)行過程中定位的精度。在列車定位過程中,采用多星座系統(tǒng)可在相同路況下獲得更多的可見星,從而提高RAIM算法的可用性及信息冗余度。在故障檢測部分,利用奇偶矢量法計(jì)算量少、檢測速度快的優(yōu)勢,可實(shí)現(xiàn)故障衛(wèi)星的快速檢測。結(jié)合青藏鐵路列車定位數(shù)據(jù)分別從算法可用性、故障檢測識別和檢測性能三個(gè)方面進(jìn)行了仿真分析。結(jié)果表明,多星座 RAIM 算法的可用性較單星座系統(tǒng)更好,且能迅速地檢測出定位過程中出現(xiàn)的故障衛(wèi)星,并予以識別。由于信息冗余度的增加,多星座系統(tǒng)的故障檢測率和故障識別率較單星座系統(tǒng)均有所提高,位置精度因子較單星座系統(tǒng)有明顯的降低。多星座RAIM技術(shù)可有效提高列車定位精度,保證鐵路運(yùn)輸行車安全。

    猜你喜歡
    偽距可用性星座
    基于文獻(xiàn)計(jì)量學(xué)的界面設(shè)計(jì)可用性中外對比研究
    包裝工程(2023年24期)2023-12-27 09:18:26
    基于輻射傳輸模型的GOCI晨昏時(shí)段數(shù)據(jù)的可用性分析
    北斗偽距觀測值精度分析
    星座
    健康女性(2017年3期)2017-04-27 21:34:31
    12星座之我愛洗澡
    GNSS偽距粗差的開窗探測及修復(fù)
    星座
    都市麗人(2016年12期)2016-12-14 09:54:15
    星座
    都市麗人(2016年7期)2016-08-03 10:06:05
    聯(lián)合碼偽距和載波寬巷組合的相對定位技術(shù)研究
    空客A320模擬機(jī)FD1+2可用性的討論
    河南科技(2015年7期)2015-03-11 16:23:13
    亚洲精品粉嫩美女一区| 亚洲av第一区精品v没综合| 757午夜福利合集在线观看| 在线观看免费日韩欧美大片| 男女下面进入的视频免费午夜 | 免费观看人在逋| 欧美老熟妇乱子伦牲交| 动漫黄色视频在线观看| 亚洲情色 制服丝袜| 亚洲精品国产色婷婷电影| 9191精品国产免费久久| 亚洲av五月六月丁香网| 嫩草影视91久久| 99re在线观看精品视频| 国产一区二区在线av高清观看| 亚洲精品美女久久av网站| 亚洲欧美一区二区三区黑人| 18禁黄网站禁片午夜丰满| 亚洲国产精品久久男人天堂| 黄色成人免费大全| 午夜福利欧美成人| 亚洲人成77777在线视频| a在线观看视频网站| 美女 人体艺术 gogo| 99香蕉大伊视频| 麻豆一二三区av精品| 丝袜在线中文字幕| av视频免费观看在线观看| 国产成人精品在线电影| 这个男人来自地球电影免费观看| 午夜福利成人在线免费观看| 麻豆一二三区av精品| 在线观看www视频免费| 美女免费视频网站| 午夜亚洲福利在线播放| 老司机午夜十八禁免费视频| 精品国产乱子伦一区二区三区| 亚洲av五月六月丁香网| 国产不卡一卡二| 国产精品日韩av在线免费观看 | 日本黄色视频三级网站网址| 男女做爰动态图高潮gif福利片 | 真人做人爱边吃奶动态| ponron亚洲| 久久精品国产清高在天天线| 亚洲av第一区精品v没综合| 亚洲人成伊人成综合网2020| 精品久久久久久久人妻蜜臀av | 人成视频在线观看免费观看| 欧美成人免费av一区二区三区| 中文字幕高清在线视频| 一区二区三区高清视频在线| 19禁男女啪啪无遮挡网站| 18美女黄网站色大片免费观看| 亚洲成国产人片在线观看| 中文字幕av电影在线播放| 又紧又爽又黄一区二区| 国产一卡二卡三卡精品| 最近最新免费中文字幕在线| 色av中文字幕| 9热在线视频观看99| 日韩成人在线观看一区二区三区| 免费看美女性在线毛片视频| 国产区一区二久久| 成人精品一区二区免费| 99久久久亚洲精品蜜臀av| 久热这里只有精品99| 制服丝袜大香蕉在线| 日韩有码中文字幕| 久久 成人 亚洲| 亚洲五月色婷婷综合| 在线观看一区二区三区| 日韩av在线大香蕉| 少妇的丰满在线观看| 精品高清国产在线一区| 嫩草影院精品99| 麻豆av在线久日| 亚洲国产中文字幕在线视频| 老司机在亚洲福利影院| 亚洲精品久久成人aⅴ小说| av在线天堂中文字幕| 欧美日韩一级在线毛片| 色播在线永久视频| 日韩成人在线观看一区二区三区| 日韩av在线大香蕉| 久久久久精品国产欧美久久久| 欧美性长视频在线观看| av在线播放免费不卡| 国产精品 国内视频| 亚洲国产欧美网| av片东京热男人的天堂| 欧美人与性动交α欧美精品济南到| 久久久精品国产亚洲av高清涩受| 免费在线观看黄色视频的| 最近最新免费中文字幕在线| 成人手机av| 欧美色欧美亚洲另类二区 | 美国免费a级毛片| 精品不卡国产一区二区三区| 国产精品免费视频内射| 国产欧美日韩一区二区三区在线| 国产精品久久久av美女十八| 日韩 欧美 亚洲 中文字幕| 人人妻人人爽人人添夜夜欢视频| 国产亚洲欧美在线一区二区| 日本免费a在线| 久久精品成人免费网站| 三级毛片av免费| 色综合欧美亚洲国产小说| 中文字幕久久专区| 亚洲欧美精品综合一区二区三区| 亚洲成人久久性| 国产精品香港三级国产av潘金莲| 国产亚洲av高清不卡| 在线视频色国产色| 成在线人永久免费视频| 色尼玛亚洲综合影院| 一二三四在线观看免费中文在| 午夜福利一区二区在线看| 一级,二级,三级黄色视频| √禁漫天堂资源中文www| 国产精品乱码一区二三区的特点 | 99香蕉大伊视频| 一级毛片高清免费大全| 美女午夜性视频免费| 久久亚洲真实| 免费一级毛片在线播放高清视频 | 亚洲avbb在线观看| 黑人巨大精品欧美一区二区mp4| 在线免费观看的www视频| 在线观看午夜福利视频| 高潮久久久久久久久久久不卡| 麻豆久久精品国产亚洲av| 亚洲,欧美精品.| 午夜精品在线福利| 国产成+人综合+亚洲专区| 亚洲第一av免费看| 91国产中文字幕| 精品国内亚洲2022精品成人| 成人三级黄色视频| 日韩av在线大香蕉| 99国产精品99久久久久| 国产高清激情床上av| 少妇裸体淫交视频免费看高清 | 国产私拍福利视频在线观看| 男人的好看免费观看在线视频 | 一边摸一边抽搐一进一小说| 国内精品久久久久精免费| www日本在线高清视频| 女人高潮潮喷娇喘18禁视频| 亚洲黑人精品在线| 日韩大尺度精品在线看网址 | www.精华液| 女同久久另类99精品国产91| 亚洲三区欧美一区| 看黄色毛片网站| 午夜福利高清视频| 老鸭窝网址在线观看| 中文字幕最新亚洲高清| 97人妻天天添夜夜摸| 国产精品久久久人人做人人爽| 如日韩欧美国产精品一区二区三区| 国产日韩一区二区三区精品不卡| 老熟妇仑乱视频hdxx| 久久久久久久午夜电影| 少妇的丰满在线观看| 久热爱精品视频在线9| 亚洲成av片中文字幕在线观看| 好男人在线观看高清免费视频 | 精品国产超薄肉色丝袜足j| 国产精品日韩av在线免费观看 | 美女高潮喷水抽搐中文字幕| 国产区一区二久久| 黄网站色视频无遮挡免费观看| 亚洲片人在线观看| 精品久久久精品久久久| 久久狼人影院| 国产成人精品久久二区二区91| 久久伊人香网站| 国产成人欧美| 91av网站免费观看| 国产精品 欧美亚洲| 女人被躁到高潮嗷嗷叫费观| 精品免费久久久久久久清纯| 夜夜躁狠狠躁天天躁| 日本免费一区二区三区高清不卡 | 中文字幕av电影在线播放| 黄色 视频免费看| 国产精品日韩av在线免费观看 | 国产精品一区二区三区四区久久 | 亚洲男人的天堂狠狠| 黄网站色视频无遮挡免费观看| 亚洲久久久国产精品| e午夜精品久久久久久久| 美国免费a级毛片| 国产蜜桃级精品一区二区三区| 国产av又大| 9191精品国产免费久久| 亚洲五月天丁香| 一夜夜www| 久久国产精品人妻蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 香蕉国产在线看| 亚洲国产欧美网| 国产主播在线观看一区二区| 久9热在线精品视频| 国产熟女xx| 亚洲第一av免费看| www.www免费av| 久久人人爽av亚洲精品天堂| 在线视频色国产色| 欧美日韩福利视频一区二区| 国产乱人伦免费视频| 99在线人妻在线中文字幕| 91麻豆av在线| 不卡一级毛片| 国产精品日韩av在线免费观看 | 久久精品国产清高在天天线| 人人妻人人澡人人看| 亚洲色图av天堂| 成人免费观看视频高清| 最新美女视频免费是黄的| 97碰自拍视频| 久久精品成人免费网站| 不卡av一区二区三区| 欧美午夜高清在线| 老司机午夜十八禁免费视频| 日本黄色视频三级网站网址| 69av精品久久久久久| 久久久久精品国产欧美久久久| www.自偷自拍.com| 母亲3免费完整高清在线观看| 热re99久久国产66热| www.自偷自拍.com| 一级,二级,三级黄色视频| 日韩大尺度精品在线看网址 | 熟女少妇亚洲综合色aaa.| 精品国产亚洲在线| 国产三级黄色录像| 久久热在线av| 午夜福利欧美成人| 午夜精品久久久久久毛片777| 欧美乱色亚洲激情| 国产午夜精品久久久久久| 欧美一级毛片孕妇| 国产一区二区三区视频了| 色综合站精品国产| 国产日韩一区二区三区精品不卡| 午夜a级毛片| av福利片在线| 国产精品1区2区在线观看.| 国产精品自产拍在线观看55亚洲| 精品久久久久久久久久免费视频| 亚洲国产毛片av蜜桃av| 欧美日本中文国产一区发布| 多毛熟女@视频| 搡老熟女国产l中国老女人| 国产精品久久久人人做人人爽| e午夜精品久久久久久久| 国产一区二区三区视频了| 国产精品久久久久久人妻精品电影| 欧洲精品卡2卡3卡4卡5卡区| 日韩成人在线观看一区二区三区| 欧美日韩一级在线毛片| 999精品在线视频| 18禁美女被吸乳视频| 欧美大码av| 97人妻精品一区二区三区麻豆 | 免费女性裸体啪啪无遮挡网站| 国产精品日韩av在线免费观看 | 久久久水蜜桃国产精品网| 97超级碰碰碰精品色视频在线观看| 亚洲欧美精品综合一区二区三区| 国产精品精品国产色婷婷| 多毛熟女@视频| 国产三级在线视频| 欧美日韩福利视频一区二区| 动漫黄色视频在线观看| 亚洲欧美一区二区三区黑人| 热re99久久国产66热| 中文字幕人妻丝袜一区二区| 国产成人啪精品午夜网站| 丝袜在线中文字幕| 国产精品久久视频播放| 老司机深夜福利视频在线观看| 欧美老熟妇乱子伦牲交| 丰满人妻熟妇乱又伦精品不卡| 婷婷精品国产亚洲av在线| 桃红色精品国产亚洲av| 一边摸一边做爽爽视频免费| 电影成人av| 国产伦一二天堂av在线观看| 国产精品日韩av在线免费观看 | 欧美日韩亚洲综合一区二区三区_| www.精华液| 激情视频va一区二区三区| 91成年电影在线观看| 叶爱在线成人免费视频播放| 精品人妻1区二区| 色播亚洲综合网| 黄色丝袜av网址大全| 丝袜美足系列| 女人被躁到高潮嗷嗷叫费观| 亚洲狠狠婷婷综合久久图片| 国产精品亚洲av一区麻豆| 在线观看免费视频日本深夜| 久久香蕉国产精品| 午夜老司机福利片| 久久精品亚洲熟妇少妇任你| 搡老岳熟女国产| 欧美成人一区二区免费高清观看 | 亚洲人成网站在线播放欧美日韩| 久久香蕉国产精品| 黄色片一级片一级黄色片| 无限看片的www在线观看| 免费看a级黄色片| 欧美成人性av电影在线观看| 操美女的视频在线观看| 久久国产精品人妻蜜桃| 99久久国产精品久久久| 精品国产乱子伦一区二区三区| 国产99白浆流出| 在线天堂中文资源库| 一进一出抽搐动态| 熟妇人妻久久中文字幕3abv| 一级a爱片免费观看的视频| 一级,二级,三级黄色视频| 国产人伦9x9x在线观看| 黄色成人免费大全| 国产真人三级小视频在线观看| 男女下面插进去视频免费观看| 精品无人区乱码1区二区| 757午夜福利合集在线观看| 日韩 欧美 亚洲 中文字幕| 久久久久九九精品影院| 欧美最黄视频在线播放免费| 国产精品久久久久久亚洲av鲁大| 久久香蕉激情| 身体一侧抽搐| 一边摸一边抽搐一进一小说| 亚洲成人免费电影在线观看| 18美女黄网站色大片免费观看| 757午夜福利合集在线观看| 黄色视频,在线免费观看| 日日爽夜夜爽网站| av有码第一页| 波多野结衣巨乳人妻| 亚洲aⅴ乱码一区二区在线播放 | 可以在线观看的亚洲视频| 黄网站色视频无遮挡免费观看| a级毛片在线看网站| 欧美绝顶高潮抽搐喷水| 成在线人永久免费视频| 搡老熟女国产l中国老女人| 亚洲av电影在线进入| 黄色毛片三级朝国网站| 亚洲中文av在线| 叶爱在线成人免费视频播放| 亚洲国产欧美网| 欧美大码av| 免费在线观看亚洲国产| 一进一出抽搐gif免费好疼| 欧美色视频一区免费| 日韩中文字幕欧美一区二区| 18禁裸乳无遮挡免费网站照片 | 亚洲国产精品合色在线| а√天堂www在线а√下载| 中文字幕精品免费在线观看视频| 我的亚洲天堂| www国产在线视频色| 老熟妇仑乱视频hdxx| 真人一进一出gif抽搐免费| 在线观看www视频免费| www.www免费av| 少妇 在线观看| 黄频高清免费视频| 国产精品99久久99久久久不卡| cao死你这个sao货| 一二三四在线观看免费中文在| 免费搜索国产男女视频| 少妇 在线观看| 日韩高清综合在线| 一二三四在线观看免费中文在| av网站免费在线观看视频| 9热在线视频观看99| 人成视频在线观看免费观看| 激情在线观看视频在线高清| 人成视频在线观看免费观看| 亚洲欧美激情在线| 男人的好看免费观看在线视频 | 欧美性长视频在线观看| cao死你这个sao货| 一级a爱视频在线免费观看| 久久亚洲真实| 亚洲人成网站在线播放欧美日韩| 丝袜人妻中文字幕| 精品不卡国产一区二区三区| 青草久久国产| 一区二区三区高清视频在线| 在线观看免费午夜福利视频| 欧美一区二区精品小视频在线| 亚洲人成电影免费在线| 97碰自拍视频| 亚洲第一欧美日韩一区二区三区| 欧美av亚洲av综合av国产av| 亚洲精品国产一区二区精华液| 在线观看一区二区三区| 正在播放国产对白刺激| 国产一区二区三区综合在线观看| 91精品国产国语对白视频| 欧美日本中文国产一区发布| 丝袜美足系列| 午夜亚洲福利在线播放| 9热在线视频观看99| 一级,二级,三级黄色视频| 亚洲五月色婷婷综合| 中国美女看黄片| 亚洲中文字幕日韩| 桃红色精品国产亚洲av| 午夜精品国产一区二区电影| 热99re8久久精品国产| 97碰自拍视频| 中文字幕色久视频| 99re在线观看精品视频| 最近最新中文字幕大全电影3 | 91av网站免费观看| 激情在线观看视频在线高清| 大香蕉久久成人网| 亚洲午夜理论影院| 亚洲精品国产区一区二| 久久香蕉激情| av中文乱码字幕在线| 在线av久久热| av欧美777| 国产精品久久电影中文字幕| 日韩欧美一区二区三区在线观看| 欧美大码av| 久久精品影院6| 91精品国产国语对白视频| 亚洲人成77777在线视频| 不卡av一区二区三区| 老鸭窝网址在线观看| 亚洲人成网站在线播放欧美日韩| 日韩有码中文字幕| 日本五十路高清| 夜夜夜夜夜久久久久| 国产在线观看jvid| 亚洲第一电影网av| 深夜精品福利| 一级毛片高清免费大全| 午夜亚洲福利在线播放| 一级a爱视频在线免费观看| 午夜精品在线福利| 嫩草影视91久久| 欧美中文日本在线观看视频| 久热这里只有精品99| 午夜久久久久精精品| 99热只有精品国产| 可以在线观看毛片的网站| 国产精品日韩av在线免费观看 | 人成视频在线观看免费观看| 亚洲欧美激情在线| 18禁黄网站禁片午夜丰满| 免费看a级黄色片| 亚洲精品国产一区二区精华液| 久久精品国产清高在天天线| 可以在线观看毛片的网站| 国内毛片毛片毛片毛片毛片| 日日摸夜夜添夜夜添小说| 久久中文看片网| 一卡2卡三卡四卡精品乱码亚洲| 乱人伦中国视频| 在线观看免费日韩欧美大片| 亚洲va日本ⅴa欧美va伊人久久| АⅤ资源中文在线天堂| 久久久国产成人精品二区| 妹子高潮喷水视频| 视频在线观看一区二区三区| 国产精品亚洲美女久久久| 色综合站精品国产| 久久青草综合色| 色播在线永久视频| 88av欧美| 亚洲精品国产色婷婷电影| 国产一卡二卡三卡精品| 美女扒开内裤让男人捅视频| 欧美av亚洲av综合av国产av| 久久草成人影院| 婷婷精品国产亚洲av在线| x7x7x7水蜜桃| 色在线成人网| 欧美丝袜亚洲另类 | 国产精品98久久久久久宅男小说| 日本 欧美在线| 欧美精品亚洲一区二区| 啪啪无遮挡十八禁网站| 757午夜福利合集在线观看| av视频在线观看入口| 亚洲国产精品合色在线| 亚洲av熟女| 亚洲一区中文字幕在线| 亚洲欧美激情在线| av网站免费在线观看视频| 老熟妇乱子伦视频在线观看| 日韩三级视频一区二区三区| 宅男免费午夜| 女人精品久久久久毛片| 神马国产精品三级电影在线观看 | 最近最新中文字幕大全电影3 | 99国产精品一区二区蜜桃av| 成人欧美大片| 十八禁网站免费在线| 国产人伦9x9x在线观看| 中文字幕人妻丝袜一区二区| 一级毛片女人18水好多| videosex国产| 咕卡用的链子| 国产av精品麻豆| 欧美久久黑人一区二区| av在线播放免费不卡| 狂野欧美激情性xxxx| 涩涩av久久男人的天堂| 国产亚洲欧美98| 一区二区三区高清视频在线| 12—13女人毛片做爰片一| 国产欧美日韩精品亚洲av| 国产日韩一区二区三区精品不卡| 国产蜜桃级精品一区二区三区| 欧美久久黑人一区二区| 久久中文看片网| 级片在线观看| 久久婷婷人人爽人人干人人爱 | 亚洲国产毛片av蜜桃av| 国产精品久久电影中文字幕| 精品乱码久久久久久99久播| 我的亚洲天堂| 亚洲人成77777在线视频| 欧美乱色亚洲激情| 久久亚洲真实| 亚洲国产高清在线一区二区三 | 亚洲欧美一区二区三区黑人| 欧美日韩精品网址| 国产成人影院久久av| 麻豆成人av在线观看| 中国美女看黄片| 亚洲成人免费电影在线观看| 97人妻精品一区二区三区麻豆 | 亚洲中文字幕一区二区三区有码在线看 | bbb黄色大片| 两个人看的免费小视频| 国产午夜精品久久久久久| netflix在线观看网站| 久久人妻福利社区极品人妻图片| 热re99久久国产66热| 国产精品av久久久久免费| 欧美日本视频| 亚洲精品中文字幕一二三四区| 久久国产亚洲av麻豆专区| 久久久久久大精品| 91在线观看av| 亚洲成av片中文字幕在线观看| 校园春色视频在线观看| 欧美一级毛片孕妇| 亚洲第一av免费看| 亚洲人成电影观看| 亚洲欧美日韩另类电影网站| 母亲3免费完整高清在线观看| 午夜福利影视在线免费观看| 国产精品99久久99久久久不卡| 欧美激情久久久久久爽电影 | 少妇裸体淫交视频免费看高清 | 国产乱人伦免费视频| 宅男免费午夜| 人妻丰满熟妇av一区二区三区| 1024香蕉在线观看| 久热这里只有精品99| 日本一区二区免费在线视频| 亚洲欧美精品综合一区二区三区| 中亚洲国语对白在线视频| 手机成人av网站| 91成年电影在线观看| 国产成+人综合+亚洲专区| 亚洲自拍偷在线| 欧美成狂野欧美在线观看| 波多野结衣高清无吗| aaaaa片日本免费| 涩涩av久久男人的天堂| 高清在线国产一区| 天堂√8在线中文| 一级毛片女人18水好多| 亚洲国产看品久久| 国产aⅴ精品一区二区三区波| 一级黄色大片毛片| 好男人电影高清在线观看| 亚洲aⅴ乱码一区二区在线播放 | а√天堂www在线а√下载| 最近最新中文字幕大全电影3 | 欧美精品亚洲一区二区| 大香蕉久久成人网| 啦啦啦观看免费观看视频高清 | 久久久国产欧美日韩av| 亚洲一码二码三码区别大吗| 亚洲人成网站在线播放欧美日韩| 国产成年人精品一区二区| 国产精品野战在线观看| 精品国产乱码久久久久久男人| 级片在线观看| 在线天堂中文资源库| 亚洲av电影不卡..在线观看| 国产成年人精品一区二区| 免费人成视频x8x8入口观看| 精品国内亚洲2022精品成人| 长腿黑丝高跟| 国产亚洲欧美精品永久| 在线观看一区二区三区|