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

    自適應卡爾曼濾波的電離層TEC預測模型改進

    2018-06-07 00:43:52王建敏黃佳鵬劉梓然祝會忠馬天明
    導航定位學報 2018年2期
    關鍵詞:電離層卡爾曼濾波殘差

    王建敏,黃佳鵬,劉梓然,祝會忠,馬天明

    (1.遼寧工程技術大學 測繪與地理科學學院,遼寧 阜新 123000;2.遼寧科技大學 土木工程學院,遼寧 鞍山 114000)

    0 引言

    由于電離層自身的不穩(wěn)定性,對于穿過該層的全球衛(wèi)星導航系統(tǒng)(global navigation satellite system,GNSS)的電波會產(chǎn)生折射、散射等物理現(xiàn)象。電離層總電子數(shù)(total electron content,TEC)是反映電離層狀況的重要屬性,電離層TEC的波動情況將會對人類生活產(chǎn)生極大的影響。近幾年業(yè)界使用GNSS針對電離層TEC進行監(jiān)測,通過探測電離層的異常情況作為研究地震發(fā)生的先驗信息。因此,研究預報模型對電離層TEC的變化情況進行高精度預報,進而揭示隨著電離層TEC急劇變化可能發(fā)生的物理現(xiàn)象及自然災害,具有重要的科研及實用的意義[1-4]。文獻[1]使用國際GNSS服務(International GNSS Service,IGS)組織公布的電離層TEC進行自回歸滑動平均模型(auto regressive and moving average model,ARMA)的建模與預測,以電離層平靜期和活躍期為研究對象,預報6 d的相對精度達到80 %以上[1]。文獻[2]使用自回歸模型(auto regressive model,AR)對電離層TEC進行預報,使用不同的原始數(shù)據(jù)分別向后預測2 d、3 d,通過比較模型的預報精度,證明AR模型針對電離層TEC具有很好的預報精度[2]。文獻[3]利用AR模型與誤差反向傳播神經(jīng)網(wǎng)絡模型(back propagation,BP)等組合模型完成電離層TEC預報,得到組合模型較單一模型可得到更好的預報結果的結論[3]。文獻[4]利用神經(jīng)網(wǎng)絡模型完成1 d的電離層TEC預報,表明訓練好的神經(jīng)網(wǎng)絡模型能夠充分反映出不同季節(jié)電離層TEC的變化。文獻[5]利用小波神經(jīng)網(wǎng)絡完成電離層TEC短期預報,在南北高緯度地區(qū)得到較高精度。以上研究均是直接基于IGS公布的數(shù)據(jù)完成預報,對于預處理方式大多以差分方式使數(shù)據(jù)平穩(wěn)化;但差分方式多是以簡單數(shù)學模型完成,并不完全遵循數(shù)據(jù)的變化規(guī)律,只能完成數(shù)據(jù)平穩(wěn)化處理。這類模型針對具有劇烈變化性的電離層TEC數(shù)據(jù)進行簡單的預處理,沒有考慮到原始數(shù)據(jù)可能有誤差的存在;因此針對電離層TEC預報的預處理十分必要。

    傳統(tǒng)卡爾曼濾波在對電離層TEC數(shù)據(jù)進行處理過程中存在一定的缺陷,可能會導致濾波發(fā)散。本文采用方差補償自適應卡爾曼濾波,用已有的電離層TEC數(shù)據(jù)對動態(tài)噪聲方差矩陣進行同步估計,減弱濾波模型所造成的誤差。神經(jīng)網(wǎng)絡具有良好的學習功能,可根據(jù)原始數(shù)據(jù)情況進行自我修正,小波分析理論可充分突出不同頻率情況的變化趨勢;小波神經(jīng)網(wǎng)絡兼具了二者的優(yōu)點,可充分擬合變化劇烈的電離層TEC變化情況。將方差補償自適應卡爾曼濾波作為數(shù)據(jù)預處理方式引入小波神經(jīng)網(wǎng)絡對電離層TEC進行預報,用以降低可能存在的噪聲對預測模型及預測精度的影響。

    1 模型基本原理

    1.1 卡爾曼濾波

    卡爾曼濾波能減少隨機噪聲對電離層TEC觀測量的影響。它具有最小無偏差性,是當前應用最廣的一種數(shù)據(jù)預處理方法。它處理數(shù)據(jù)不僅會去除突變的數(shù)據(jù),而且會保持數(shù)據(jù)原有的變化趨勢,這樣有利于預報模型的建立。離散線性系統(tǒng)的卡爾曼濾波方程包括狀態(tài)方程和觀測方程[6]為

    (1)

    L(k)=B(k)X(k)+Δ(k)

    (2)

    通過最小二乘原理可推得卡爾曼濾波方程為

    (3)

    (4)

    Dx(k)=(E-J(k)×B(k))×Dx(k-1)

    (5)

    1.2 方差補償自適應卡爾曼濾波

    方差補償自適應卡爾曼濾波就是在數(shù)據(jù)預處理過程中使用過去的電離層TEC數(shù)據(jù)對動態(tài)噪聲方差矩陣進行實時估計,用以減少因濾波模型誤差造成的噪聲。其基本思路為:使用濾波值與真實值求得估計剩余,即

    (6)

    (7)

    設式(7)中B(k)Φ(k)T=O(k,m);而

    (8)

    因為E(VT(k)V(k))=trace[E(V(k)VT(k))]=traceDVV,則有

    VT(k+i)V(k+i)=traceDVV+δk+i

    (9)

    式中δk+i為零均值隨機變量。

    令E(k)=VT(k)V(k)-trace[B(k)Φ(k)DX(k)ΦT(k)BT(k)]-traceDΔ(k),又記

    E=[E(k+1),……,E(k+N)]T,δ=[δk+i,……,δk+N]T。則有

    E=OdiagDΔΔ+δ

    (10)

    式(10)是關于diagDΔΔ的矩陣,當N≥m時,有唯一解;則記diagDΔΔ估計為

    (11)

    通過式(11)完成卡爾曼濾波即式(5)的修正[8-9]。

    使用方差補償自適應卡爾曼濾波作為電離層TEC預處理方法。方差補償自適應卡爾曼濾波是依據(jù)前一刻電離層TEC的狀態(tài)量,結合觀測方程去估計當前時刻TEC的濾波值,每一個電離層TEC數(shù)據(jù)都對應一個濾波后的濾波值,再以濾波值作為小波神經(jīng)網(wǎng)絡模型的輸入值完成改進模型的建立。

    1.3 小波神經(jīng)網(wǎng)絡

    小波神經(jīng)網(wǎng)絡(wavelet neural network,WNN)是以神經(jīng)網(wǎng)絡為大體框架,將小波分析理論作為核心內(nèi)容的神經(jīng)網(wǎng)絡。從結構形式上來看是將小波分解與BP神經(jīng)網(wǎng)絡相融合,用小波函數(shù)的平移因子充當隱層閾值作用,用小波函數(shù)代替隱層核心函數(shù)Sigmoid函數(shù)。由于小波神經(jīng)網(wǎng)絡需要參數(shù)較多,極大增強了樣本數(shù)據(jù)與模型的關聯(lián)度,使得模型比小波分解具有更優(yōu)的逼近能力。不同網(wǎng)絡層包括的神經(jīng)元個數(shù)也各不相同。

    小波神經(jīng)網(wǎng)絡核心部分表達式為

    (12)

    Morlet小波函數(shù)表達式為

    (13)

    (14)

    隱含層與輸出層之間的權值調整公式為

    (15)

    (16)

    2 實驗與結果分析

    在利用預測模型對電離層進行預測時,對于電離層TEC可能存在的誤差很少有人討論;因此筆者使用方差補償自適應卡爾曼濾波對電離層TEC進行預處理,用以降低原始電離層TEC可能存在的誤差,進一步提高電離層TEC的預測水平?;跉W洲定軌中心(Center for Orbit Determination in Europe,CODE)下載的2016年電離層網(wǎng)格點TEC數(shù)據(jù)對該方法進行實驗與分析,因自 2015 年起全球網(wǎng)格點數(shù)據(jù)以1 h為時間間隔,限于篇幅,選用2016年1月前10 d的數(shù)據(jù)。為驗證該方法是否與數(shù)據(jù)的時間和緯度有關,取年積日為第1天至第10天全球不同經(jīng)緯度的格網(wǎng)點數(shù)據(jù)作為實驗數(shù)據(jù),格網(wǎng)點的經(jīng)緯度包括(87.5°N,125°E)、(67.5°N,125°E)、(45°N,125°E)、(22.5°N,125°E)、(0°,125°E)、(22.5°S,125°E)、(45°S,125°E)、(67.5°S,125°E)、(87.5°S,125°E)[15-16]。利用AKN模型、KN模型與AN模型完成對于對上述基礎數(shù)據(jù)的預報,使用前10 d的TEC值分別預測后6 d的數(shù)據(jù),比較不同情況下的預測精度。AKN模型使用方差補償自適應卡爾曼濾波預處理,再利用小波神經(jīng)網(wǎng)絡完成預測;KN模型使用傳統(tǒng)卡爾曼濾波預處理,再利用小波神經(jīng)網(wǎng)絡完成預測;AN模型直接使用小波神經(jīng)網(wǎng)絡完成預測。為驗證模型在磁暴情況下的預測水平,使用2004年第302天至第311天作為原始數(shù)據(jù),完成2004年第312 天至第316天的數(shù)據(jù)預測,選擇格網(wǎng)點為(45°N,125°E)[17]。選用預測平均精度和殘差值分布數(shù)據(jù)作為精度分析標準。

    預測平均精度

    (17)

    式中:pi為不同預測模型i時刻的預測數(shù)據(jù);qi為CODE公布的對應時刻實際TEC數(shù)據(jù)。

    殘差值V表達式為

    V=pi-qi

    (18)

    圖1為2016-01-10—2016-01-16不同緯度的3個預測模型的預測情況比較。圖1中橫坐標為預測歷元,其中歷元間隔為1 h;縱坐標表示TEC值,以總電子含量單位(total electron content units,TECU)的個數(shù)計,1個TECU等于1016個電子每平方米。

    圖1 不同天數(shù)的預測情況比較

    由圖1中可以看出:在不同緯度完成預測不同天數(shù)情況下,3個模型都有很高的擬合程度;利用方差補償自適應卡爾曼濾波完成對于原始數(shù)據(jù)的預處理,相對另外2種模型具有更高的平均精度。

    表1為3種不同模型在各個緯度下使用10 d數(shù)據(jù)完成預測的1~6 d的平均精度。

    表1 不同緯度的相對精度統(tǒng)計

    從表1中可以看出,在南半球AN模型、KN模型、AKN模型能夠取得較好的預測精度,而且在(87.5°S,125°E)區(qū)域預測精度最高,在該區(qū)域經(jīng)過方差補償自適應卡爾曼濾波預處理的數(shù)據(jù)再利用模型進行預測會有更高的預測精度。由于是由小波神經(jīng)網(wǎng)絡完成預測,可能導致預測精度在某一區(qū)域較低,如(0°,125°E)區(qū)域,可能是由于該區(qū)域在這一時段不穩(wěn)定造成的;但是經(jīng)過卡爾曼濾波的預處理,提高了模型的總體預測精度:從表中可得出AKN模型預測精度相對KN模型和AN模型有不同程度的提高。從表中亦可看出緯度越高預測模型預測精度相應提高,(87.5°N,125°E)和(87.5°S,125°E)是各自半球預測精度最高的。高緯度地區(qū)是電離層異常活躍的地區(qū),結合本次實驗預測結果和原始數(shù)據(jù)曲線可以發(fā)現(xiàn),在實驗區(qū)間處2個高緯度區(qū)域的電離層TEC變化雖然不規(guī)律,但是電離層TEC變化范圍相對其他緯度更小,可能是導致預測模型精度更高的原因。

    圖2為2004年磁暴現(xiàn)象預測結果,2004年CODE公布的數(shù)據(jù)歷元為2 h,5 d共有60個TEC數(shù)據(jù),圖中第20個歷元至第28個歷元能明顯看出與周圍數(shù)據(jù)不同,通過比較3種模型擬合情況結合表1中數(shù)據(jù)可以看出,在磁暴預測情況下,AKN模型相對AN模型、KN模型有更好的預測效果。圖3為不同模型不同緯度的平均預測精度。

    表2為不同模型在不同緯度預測殘差小于1個TECU的統(tǒng)計情況(預測1~6 d的平均值)。從表中可知:AN模型有16.64 %的殘差值小于1個TECU;AKN模型有20.37 %的殘差值小于1個TECU;在各個緯度AKN模型預測殘差小于1個TECU所占的比重相對另外2個模型都有不同程度的提高。

    圖2 磁暴現(xiàn)象預測5 d結果

    圖3 不同緯度的相對精度統(tǒng)計圖

    結合表1、表2會發(fā)現(xiàn):在各個緯度殘差分布情況和預測的相對精度趨勢不相同,這是由于預測精度的計算結果在很大程度上取決于基礎數(shù)據(jù)的大小,相同或相近的預報殘差,若基值不同,相對精度會有較大差異;各個緯度原始TEC數(shù)據(jù)有很大差異,如(67.5°N,125°E)的數(shù)量級是5個TECU,(0°,125°E)的數(shù)量級是50個TECU,所以即使殘差較大,也可能具有略高的平均精度。綜合2種精度評價標準來講,使用方差補償自適應卡爾曼濾波預處理可以提高模型的預測精度,是一種比較理想的提高預測精度的方式。

    表2 不同模型預測殘差小于1個TECU所占百分比

    表3給出預測不同天數(shù)平均精度的統(tǒng)計結果。

    表3 預測不同天數(shù)的平均精度統(tǒng)計

    圖4為預測不同天數(shù)各模型平均精度的比較情況,由AKN模型變化趨勢結合表3可知,預測不同天數(shù)情況下的AKN模型預測精度在54 %~69 %范圍內(nèi),預測平均精度為60.27 %,預測平均精度相對直接使用原始數(shù)據(jù)進行預測提高4.93 %,且在預測6 d過程中保持逐步降低的總體趨勢,超過4 d其相對精度將隨時間明顯降低。

    圖4 不同天數(shù)各模型的相對精度比較

    圖5、表4為預測不同天數(shù)各模型殘差小于1個TECU的統(tǒng)計情況,從中可知:根據(jù)AKN模型的變化過程,預報殘差小于1個TECU的占20.37 %,相對直接使用原始數(shù)據(jù)建模提高3.73 %;且根據(jù)變化曲線來看,AKN模型隨著預測天數(shù)的變化,殘差也在略微改變,但是仍然優(yōu)于另外2個模型。

    圖5 不同模型預測殘差小于1個TECU預測天數(shù)統(tǒng)計

    預測時間/dAN模型的百分比/(%)KN模型的百分比/(%)AKN模型的百分比/(%)116.6722.6827.77219.4420.6021.29316.2017.5921.45416.2015.1616.55516.5715.2718.24613.1915.6616.24平均值16.6417.8320.37

    3 結束語

    本文研究表明,使用方差補償自適應卡爾曼濾波對電離層TEC進行預處理,再利用小波神經(jīng)網(wǎng)絡完成預測結果較為理想。無論是高緯度、中緯度,還是低緯度地區(qū),經(jīng)過自適應卡爾曼濾波預處理的電離層數(shù)據(jù),再利用小波神經(jīng)網(wǎng)絡模型完成預測,能夠在不同程度上提高預測精度,降低殘差值的影響。使用該方法預測樣本長度會對模型有很大影響,其大體規(guī)律為隨時間增加,模型預測精度隨之降低。

    [1] 張小紅,任曉東,吳風波,等.自回歸移動平均模型的電離層總電子含量短期預報[J].測繪學報,2014,34(5):118-124.

    [2] 李秀海,郭達志.應用半?yún)?shù)AR模型的電離層TEC建模與預測[J].測繪科學,2011,36(2):149-151.

    [3] 王建敏,黃佳鵬,祝會忠,等.電離層總電子數(shù)預報方法研究[J].測繪科學,2016,41(12):47-53.

    [4] 李淑慧,彭軍還,徐偉超,等.利用神經(jīng)網(wǎng)絡預報短期電離層TEC變化[J].測繪科學,2013,38(1):8-12.

    [5] 劉先冬,宋力杰,楊曉輝,等.基于小波神經(jīng)網(wǎng)絡的電離層TEC短期預報[J].海洋測繪,2010,30(5):49-51.

    [6] 鄭崴,張貴賓.自適應卡爾曼濾波在航空重力異常解算的應用研究[J].地球物理學報,2016,59(4):1275-1283.

    [7] 茹國寶,馮鵬,鄧勇.帶狀態(tài)切換自適應濾波算法及其在GPS/INS中的應用[J].武漢大學學報(工學版),2015,48(4):574-579.

    [8] 李奕.自適應卡爾曼濾波在變形監(jiān)測數(shù)據(jù)處理中的應用研究:以太原萬達廣場A2區(qū)基坑監(jiān)測為例[D].成都:成都理工大學,2012:21-30.

    [9] 王宇譜,呂志平,陳正生,等.一種新的鐘差預處理方法及在WNN鐘差長期預報中的應用[J].武漢大學學報(信息科學版),2016,41(3):373-379.

    [10] 王宇譜,呂志平,陳正生,等.基于灰色模型與小波神經(jīng)網(wǎng)絡的衛(wèi)星鐘差預報方法[J].大地測量與地球動力學,2014,34(3):155-159.

    [11] 任強強,王躍鋼,騰紅磊,等.基于小波神經(jīng)網(wǎng)絡的盲區(qū)重力數(shù)據(jù)預測[J].大地測量與地球動力學,2016,36(4):359-363.

    [12] 葛玉輝,熊永良,陳志勝,等.基于小波神經(jīng)網(wǎng)絡的GPS可降水量預測研究[J].測繪科學,2015,40(9):28-32.

    [13] 陳標.短期電力負荷的小波神經(jīng)網(wǎng)絡預測[D].長沙:湖南大學,2012.

    [14] 鮑亞東,劉長健,柴洪洲.小波分解改進電離層VTEC時間序列預報模型[J].大地測量與地球動力學,2015,35(5):784-787.

    [15] 李磊,張淑芳,胡青,等.ARMA模型在地震電離層TEC異常探測中的應用[J].大地測量與地球動力學,2015,35(1):62-66.

    [16] 時爽爽,黃勁松,馮建迪.基于IGS的白天和夜間電離層TEC變化特性分析[J].大地測量與地球動力學,2015,35(2):293-297.

    [17] 馮健,鄧忠新,甄衛(wèi)民,等.2004年11月強磁暴期間中國地區(qū)電離層TEC擾動特性分析[J].電波科學學報,2016,31(1):157-165.

    猜你喜歡
    電離層卡爾曼濾波殘差
    基于雙向GRU與殘差擬合的車輛跟馳建模
    一種電離層TEC格點預測模型
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    基于殘差學習的自適應無人機目標跟蹤算法
    基于遞歸殘差網(wǎng)絡的圖像超分辨率重建
    自動化學報(2019年6期)2019-07-23 01:18:32
    基于遞推更新卡爾曼濾波的磁偶極子目標跟蹤
    電離層對中高軌SAR影響機理研究
    雷達學報(2017年6期)2017-03-26 07:53:02
    基于模糊卡爾曼濾波算法的動力電池SOC估計
    電源技術(2016年9期)2016-02-27 09:05:39
    基于擴展卡爾曼濾波的PMSM無位置傳感器控制
    電源技術(2015年1期)2015-08-22 11:16:28
    Linux Shell語言在電離層解算中的應用
    一夜夜www| 国产精品一区二区三区四区久久| 久久久久久人人人人人| 黄色成人免费大全| 两性午夜刺激爽爽歪歪视频在线观看| 成人欧美大片| 99精品在免费线老司机午夜| 午夜久久久久精精品| 国产亚洲精品一区二区www| 夜夜看夜夜爽夜夜摸| 久久欧美精品欧美久久欧美| 亚洲美女黄片视频| 色老头精品视频在线观看| 欧美+亚洲+日韩+国产| 少妇的丰满在线观看| 亚洲精品久久国产高清桃花| 91久久精品国产一区二区成人 | 婷婷六月久久综合丁香| 少妇丰满av| 91九色精品人成在线观看| 久久久久久久久中文| 久久久久亚洲av毛片大全| 99热这里只有是精品50| 成人亚洲精品av一区二区| 国产日本99.免费观看| 老司机午夜十八禁免费视频| 少妇人妻精品综合一区二区 | 女人高潮潮喷娇喘18禁视频| 在线视频色国产色| av欧美777| 日韩欧美一区二区三区在线观看| 麻豆成人av在线观看| 51午夜福利影视在线观看| 麻豆国产97在线/欧美| 久久精品影院6| 亚洲av成人av| 久久久国产精品麻豆| 99热只有精品国产| 亚洲激情在线av| 久久亚洲精品不卡| 久久亚洲精品不卡| 丰满人妻一区二区三区视频av | 免费高清视频大片| АⅤ资源中文在线天堂| www.999成人在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲精品美女久久久久99蜜臀| 日本精品一区二区三区蜜桃| 一夜夜www| 久久久久性生活片| 久久久久免费精品人妻一区二区| 午夜福利18| 色综合站精品国产| 中文字幕精品亚洲无线码一区| 日日摸夜夜添夜夜添小说| 午夜两性在线视频| 黄片小视频在线播放| 母亲3免费完整高清在线观看| 免费看日本二区| 久久精品国产综合久久久| 久久精品91蜜桃| 日韩欧美一区二区三区在线观看| 午夜福利视频1000在线观看| 国产探花在线观看一区二区| 国内精品美女久久久久久| 久久久成人免费电影| 日本a在线网址| av片东京热男人的天堂| 国产午夜精品论理片| 又黄又粗又硬又大视频| 我的老师免费观看完整版| 国产91精品成人一区二区三区| 天天躁日日操中文字幕| 国内久久婷婷六月综合欲色啪| 韩国av一区二区三区四区| 国产一级毛片七仙女欲春2| 免费av不卡在线播放| 有码 亚洲区| 国产一级毛片七仙女欲春2| 久久香蕉精品热| 成人国产一区最新在线观看| 精品一区二区三区人妻视频| 动漫黄色视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 成人高潮视频无遮挡免费网站| 老熟妇乱子伦视频在线观看| 欧美日韩综合久久久久久 | 中文字幕av成人在线电影| 精品电影一区二区在线| 欧美又色又爽又黄视频| 少妇人妻一区二区三区视频| 啦啦啦免费观看视频1| 啦啦啦免费观看视频1| 亚洲国产欧洲综合997久久,| 99精品久久久久人妻精品| 久久久久久久亚洲中文字幕 | 男女做爰动态图高潮gif福利片| а√天堂www在线а√下载| 91麻豆av在线| 变态另类成人亚洲欧美熟女| 中文字幕人妻丝袜一区二区| 日韩有码中文字幕| 国产探花极品一区二区| 丁香欧美五月| 69av精品久久久久久| 欧美中文日本在线观看视频| 国语自产精品视频在线第100页| 国产主播在线观看一区二区| 免费观看的影片在线观看| 级片在线观看| 色哟哟哟哟哟哟| 18禁国产床啪视频网站| 99国产精品一区二区蜜桃av| 桃色一区二区三区在线观看| 亚洲无线观看免费| av在线蜜桃| 国产高清视频在线观看网站| 国产成人aa在线观看| 成人精品一区二区免费| 亚洲国产精品sss在线观看| 国产精品自产拍在线观看55亚洲| 亚洲av不卡在线观看| 婷婷精品国产亚洲av| 久久久成人免费电影| 免费观看的影片在线观看| 一级毛片女人18水好多| 波多野结衣高清无吗| 亚洲无线在线观看| 亚洲欧美激情综合另类| 国产免费男女视频| 久久久久久久久大av| 午夜精品在线福利| 韩国av一区二区三区四区| 一级a爱片免费观看的视频| 两个人视频免费观看高清| 高清在线国产一区| 中文亚洲av片在线观看爽| 久久久久久久久久黄片| 在线观看免费午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看| 丰满人妻一区二区三区视频av | 国产精品99久久99久久久不卡| 搡女人真爽免费视频火全软件 | 一个人看视频在线观看www免费 | 久久精品国产综合久久久| 村上凉子中文字幕在线| 最后的刺客免费高清国语| 久久人妻av系列| 精品国产超薄肉色丝袜足j| 久久精品国产清高在天天线| 99在线视频只有这里精品首页| 小蜜桃在线观看免费完整版高清| 少妇人妻一区二区三区视频| 亚洲五月天丁香| 日韩欧美国产在线观看| 一区二区三区高清视频在线| 免费观看人在逋| 91麻豆av在线| 三级国产精品欧美在线观看| 国产精品精品国产色婷婷| 国产精品亚洲av一区麻豆| 国产精品乱码一区二三区的特点| 美女黄网站色视频| 欧美色视频一区免费| 午夜久久久久精精品| 国产精品98久久久久久宅男小说| 亚洲成人精品中文字幕电影| 五月玫瑰六月丁香| 亚洲av熟女| 女人被狂操c到高潮| 欧美高清成人免费视频www| 国产高清视频在线播放一区| 长腿黑丝高跟| 俄罗斯特黄特色一大片| 一级黄片播放器| 国产高清videossex| 成人高潮视频无遮挡免费网站| 成年女人看的毛片在线观看| 亚洲精品亚洲一区二区| 黄片小视频在线播放| 91在线精品国自产拍蜜月 | www日本在线高清视频| 我的老师免费观看完整版| 久久婷婷人人爽人人干人人爱| 一个人免费在线观看的高清视频| 给我免费播放毛片高清在线观看| 一区福利在线观看| 波多野结衣高清作品| 欧美日韩亚洲国产一区二区在线观看| 国产精品野战在线观看| 黄色日韩在线| 中文字幕av成人在线电影| av国产免费在线观看| 天天添夜夜摸| 久久性视频一级片| 色av中文字幕| 少妇熟女aⅴ在线视频| 精品一区二区三区视频在线观看免费| 欧美又色又爽又黄视频| 丰满的人妻完整版| 在线国产一区二区在线| 日日摸夜夜添夜夜添小说| 麻豆久久精品国产亚洲av| 日本 av在线| 一个人免费在线观看电影| bbb黄色大片| 久久久精品大字幕| 婷婷亚洲欧美| 日韩成人在线观看一区二区三区| 国产三级黄色录像| 日本一二三区视频观看| 午夜激情欧美在线| 精品国内亚洲2022精品成人| 无限看片的www在线观看| 亚洲精品在线观看二区| 日韩欧美三级三区| 夜夜夜夜夜久久久久| 老司机在亚洲福利影院| 国产一区二区亚洲精品在线观看| svipshipincom国产片| 身体一侧抽搐| 两个人视频免费观看高清| 日本一本二区三区精品| 18禁黄网站禁片免费观看直播| 久久天躁狠狠躁夜夜2o2o| 日日夜夜操网爽| 国产午夜精品久久久久久一区二区三区 | 亚洲精品成人久久久久久| 午夜福利视频1000在线观看| 国产精品乱码一区二三区的特点| 99在线人妻在线中文字幕| 欧美3d第一页| 国产一级毛片七仙女欲春2| 日本免费a在线| 国产综合懂色| 嫩草影视91久久| 男女视频在线观看网站免费| 亚洲中文日韩欧美视频| 在线观看一区二区三区| 99精品久久久久人妻精品| 在线十欧美十亚洲十日本专区| 免费av毛片视频| 在线观看免费午夜福利视频| av福利片在线观看| 日本免费a在线| 国产色爽女视频免费观看| 国产精品影院久久| 午夜精品在线福利| 成人国产综合亚洲| 嫩草影院精品99| 欧美国产日韩亚洲一区| 亚洲美女黄片视频| 午夜精品久久久久久毛片777| 亚洲最大成人中文| 床上黄色一级片| 成人亚洲精品av一区二区| 亚洲成人久久性| а√天堂www在线а√下载| 亚洲欧美日韩高清在线视频| 在线观看一区二区三区| av片东京热男人的天堂| 大型黄色视频在线免费观看| 69av精品久久久久久| 国产精品野战在线观看| 丰满人妻一区二区三区视频av | ponron亚洲| 免费看美女性在线毛片视频| 亚洲在线观看片| 亚洲精品成人久久久久久| 精品一区二区三区视频在线 | 久9热在线精品视频| www日本在线高清视频| 久久香蕉国产精品| 久久国产精品人妻蜜桃| 国产真人三级小视频在线观看| 男人和女人高潮做爰伦理| 搞女人的毛片| 国产精品,欧美在线| 夜夜夜夜夜久久久久| 999久久久精品免费观看国产| 少妇裸体淫交视频免费看高清| 久久久久久久亚洲中文字幕 | 亚洲av熟女| 久久久久久久久中文| 黄色女人牲交| 狂野欧美激情性xxxx| 成人av一区二区三区在线看| 国模一区二区三区四区视频| 亚洲国产精品合色在线| 免费看美女性在线毛片视频| 午夜视频国产福利| 黄片小视频在线播放| 狂野欧美激情性xxxx| 国产精品久久久久久亚洲av鲁大| 成人性生交大片免费视频hd| 日本黄大片高清| 国产伦人伦偷精品视频| 国产精品一区二区免费欧美| 99国产精品一区二区三区| 99国产精品一区二区三区| 亚洲一区二区三区不卡视频| 亚洲中文字幕日韩| 免费大片18禁| 18禁裸乳无遮挡免费网站照片| 精品国内亚洲2022精品成人| 99久久精品热视频| 国产v大片淫在线免费观看| 国产精品亚洲一级av第二区| 日韩欧美免费精品| 天堂av国产一区二区熟女人妻| 美女免费视频网站| 一级作爱视频免费观看| 少妇熟女aⅴ在线视频| 亚洲18禁久久av| 99久久综合精品五月天人人| 国产高潮美女av| av视频在线观看入口| 小蜜桃在线观看免费完整版高清| 午夜免费观看网址| 欧洲精品卡2卡3卡4卡5卡区| 久久欧美精品欧美久久欧美| 久久久色成人| 三级男女做爰猛烈吃奶摸视频| 国产精品爽爽va在线观看网站| 国产aⅴ精品一区二区三区波| 日本与韩国留学比较| 一本一本综合久久| av天堂中文字幕网| 亚洲国产精品成人综合色| 欧美日本视频| 天堂动漫精品| 精品人妻一区二区三区麻豆 | 亚洲精品456在线播放app | 日本成人三级电影网站| 亚洲av熟女| 国产成+人综合+亚洲专区| 国产91精品成人一区二区三区| 国产精品久久电影中文字幕| 成人三级黄色视频| 最新在线观看一区二区三区| 午夜免费成人在线视频| 国产成人av教育| 亚洲一区二区三区色噜噜| 成人无遮挡网站| 欧美日韩瑟瑟在线播放| 国内毛片毛片毛片毛片毛片| 91在线观看av| 丰满乱子伦码专区| 免费一级毛片在线播放高清视频| 欧美黑人欧美精品刺激| 九色国产91popny在线| 中文字幕人成人乱码亚洲影| 欧美3d第一页| 国产精品影院久久| 日韩有码中文字幕| 身体一侧抽搐| 精品欧美国产一区二区三| 婷婷精品国产亚洲av| 超碰av人人做人人爽久久 | 国产蜜桃级精品一区二区三区| 亚洲一区高清亚洲精品| 神马国产精品三级电影在线观看| 国产淫片久久久久久久久 | 亚洲一区二区三区色噜噜| 久久天躁狠狠躁夜夜2o2o| 国产探花极品一区二区| 午夜精品久久久久久毛片777| 欧美区成人在线视频| 免费人成视频x8x8入口观看| 成人国产综合亚洲| 欧美黑人巨大hd| 成人高潮视频无遮挡免费网站| 俄罗斯特黄特色一大片| 一个人看视频在线观看www免费 | 一级黄色大片毛片| 欧美xxxx黑人xx丫x性爽| 亚洲av五月六月丁香网| 夜夜爽天天搞| 黄色日韩在线| 午夜免费激情av| 日韩中文字幕欧美一区二区| 男人和女人高潮做爰伦理| 午夜两性在线视频| 久久午夜亚洲精品久久| 我的老师免费观看完整版| 1024手机看黄色片| 悠悠久久av| 久久久色成人| 国产高清视频在线播放一区| 黄色片一级片一级黄色片| 99久久久亚洲精品蜜臀av| 久久精品91无色码中文字幕| 色老头精品视频在线观看| 国内毛片毛片毛片毛片毛片| 美女被艹到高潮喷水动态| 天堂√8在线中文| 欧美中文日本在线观看视频| 国产99白浆流出| 色播亚洲综合网| 网址你懂的国产日韩在线| 国产爱豆传媒在线观看| a在线观看视频网站| 高潮久久久久久久久久久不卡| www日本黄色视频网| 一级作爱视频免费观看| 亚洲人成网站在线播放欧美日韩| 美女cb高潮喷水在线观看| 色综合婷婷激情| 九九热线精品视视频播放| 亚洲精品乱码久久久v下载方式 | 动漫黄色视频在线观看| 高潮久久久久久久久久久不卡| 欧美丝袜亚洲另类 | 18禁黄网站禁片免费观看直播| 午夜福利高清视频| av福利片在线观看| 99精品欧美一区二区三区四区| 特大巨黑吊av在线直播| 国产极品精品免费视频能看的| 露出奶头的视频| 真实男女啪啪啪动态图| 欧美在线黄色| 给我免费播放毛片高清在线观看| 老司机在亚洲福利影院| 99国产综合亚洲精品| 两性午夜刺激爽爽歪歪视频在线观看| 变态另类丝袜制服| 91九色精品人成在线观看| 岛国视频午夜一区免费看| 免费av毛片视频| 老司机午夜十八禁免费视频| 999久久久精品免费观看国产| 热99在线观看视频| 国产v大片淫在线免费观看| 1000部很黄的大片| 99riav亚洲国产免费| 搞女人的毛片| 国产精品日韩av在线免费观看| 日韩欧美国产在线观看| 中文字幕人成人乱码亚洲影| 一区福利在线观看| 国内精品久久久久久久电影| 久久国产乱子伦精品免费另类| 精品久久久久久久末码| 国产又黄又爽又无遮挡在线| 18禁黄网站禁片免费观看直播| 岛国在线免费视频观看| 小蜜桃在线观看免费完整版高清| 亚洲熟妇中文字幕五十中出| 性欧美人与动物交配| 99久久久亚洲精品蜜臀av| 色精品久久人妻99蜜桃| 丁香六月欧美| 悠悠久久av| 一夜夜www| xxxwww97欧美| 在线天堂最新版资源| а√天堂www在线а√下载| 精品国内亚洲2022精品成人| 天堂动漫精品| 亚洲人与动物交配视频| 999久久久精品免费观看国产| 欧美成人性av电影在线观看| 观看免费一级毛片| 久久久久久久午夜电影| 免费在线观看日本一区| 国产一区二区在线av高清观看| 天堂影院成人在线观看| 色综合婷婷激情| 在线观看美女被高潮喷水网站 | 国产欧美日韩一区二区三| 99久久综合精品五月天人人| 亚洲 国产 在线| 麻豆成人午夜福利视频| 婷婷精品国产亚洲av| 亚洲av一区综合| 亚洲午夜理论影院| 变态另类丝袜制服| 国产探花极品一区二区| 欧美性猛交╳xxx乱大交人| 成熟少妇高潮喷水视频| 日韩欧美精品v在线| 国产欧美日韩一区二区精品| 一级黄色大片毛片| 夜夜爽天天搞| 欧美不卡视频在线免费观看| 免费搜索国产男女视频| 91九色精品人成在线观看| 级片在线观看| 亚洲片人在线观看| 午夜福利免费观看在线| 天美传媒精品一区二区| 成熟少妇高潮喷水视频| 免费在线观看成人毛片| 又粗又爽又猛毛片免费看| 欧美zozozo另类| 色综合婷婷激情| 精品乱码久久久久久99久播| 亚洲第一电影网av| 午夜免费观看网址| 国产伦精品一区二区三区四那| av国产免费在线观看| 国产v大片淫在线免费观看| 国产三级中文精品| 热99re8久久精品国产| 久久久久久九九精品二区国产| 给我免费播放毛片高清在线观看| 在线天堂最新版资源| 黄色片一级片一级黄色片| 成人高潮视频无遮挡免费网站| 亚洲aⅴ乱码一区二区在线播放| 国产成人av激情在线播放| 热99在线观看视频| 黄色视频,在线免费观看| 特级一级黄色大片| 国产高清激情床上av| 日本免费a在线| 午夜视频国产福利| 床上黄色一级片| 欧美大码av| 12—13女人毛片做爰片一| 51午夜福利影视在线观看| 免费观看精品视频网站| 日本 av在线| 日日摸夜夜添夜夜添小说| 丰满人妻一区二区三区视频av | 男女视频在线观看网站免费| av天堂中文字幕网| 色精品久久人妻99蜜桃| 免费在线观看日本一区| 亚洲av电影不卡..在线观看| 亚洲欧美日韩无卡精品| 成人特级av手机在线观看| 老司机在亚洲福利影院| 少妇丰满av| 怎么达到女性高潮| 97超视频在线观看视频| 亚洲最大成人手机在线| 国产精品综合久久久久久久免费| 操出白浆在线播放| 天堂影院成人在线观看| 国产精品综合久久久久久久免费| 99热6这里只有精品| 黄色成人免费大全| 精品熟女少妇八av免费久了| 久久久精品大字幕| 日本a在线网址| 欧美+日韩+精品| 男女之事视频高清在线观看| 搞女人的毛片| 国产黄色小视频在线观看| 欧美+日韩+精品| 69人妻影院| 99久久久亚洲精品蜜臀av| 欧美日韩中文字幕国产精品一区二区三区| 在线播放无遮挡| 黄色视频,在线免费观看| 波多野结衣巨乳人妻| 长腿黑丝高跟| 欧美日韩精品网址| 国产精品三级大全| 在线十欧美十亚洲十日本专区| 亚洲精华国产精华精| 亚洲av二区三区四区| 美女cb高潮喷水在线观看| 在线天堂最新版资源| 久久久国产成人精品二区| 高清在线国产一区| 欧美一区二区亚洲| 91九色精品人成在线观看| 国产一区二区激情短视频| 可以在线观看毛片的网站| 丰满人妻一区二区三区视频av | 欧美区成人在线视频| 最近最新免费中文字幕在线| 国产爱豆传媒在线观看| 国产精品日韩av在线免费观看| 久久久精品大字幕| 2021天堂中文幕一二区在线观| 别揉我奶头~嗯~啊~动态视频| 欧美乱码精品一区二区三区| 中文字幕人妻丝袜一区二区| 久久久久久久亚洲中文字幕 | 国产精品99久久久久久久久| 成年人黄色毛片网站| 国产精品亚洲一级av第二区| 亚洲成a人片在线一区二区| 国产三级中文精品| 观看免费一级毛片| 香蕉av资源在线| 国产av一区在线观看免费| 国产高清激情床上av| 日韩欧美在线二视频| 99国产综合亚洲精品| 成人av在线播放网站| 日韩人妻高清精品专区| 精品欧美国产一区二区三| 久久久久久九九精品二区国产| 亚洲国产欧美人成| 免费人成在线观看视频色| av在线蜜桃| 国产视频内射| 3wmmmm亚洲av在线观看| 操出白浆在线播放| 美女大奶头视频| 在线观看美女被高潮喷水网站 | 亚洲国产精品999在线| 非洲黑人性xxxx精品又粗又长| 国产三级中文精品| 久久草成人影院| 麻豆国产97在线/欧美| 99热这里只有精品一区| or卡值多少钱| 国产熟女xx| 一级作爱视频免费观看|