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

    基于RMSprop的高精度GNSS-R鏡點快速預(yù)測算法

    2021-10-15 04:06:22盛志超顧斯祺孟婉婷
    無線電工程 2021年10期
    關(guān)鍵詞:鏡面反射二分法線段

    盛志超,周 勃,秦 瑾,顧斯祺,孟婉婷

    (上海航天電子技術(shù)研究所,上海 201109)

    0 引言

    全球?qū)Ш叫l(wèi)星系統(tǒng)反射信號(Global Navigation Satellite System Reflectometry,GNSS-R)技術(shù)作為外源助動式新興遙感技術(shù),以其寬刈幅、高時空分辨率、低功耗、低成本和低重量等特點受到了國內(nèi)外遙感領(lǐng)域的廣泛關(guān)注[1]。GNSS-R技術(shù)能夠進(jìn)行中尺度海洋環(huán)境要素的探測,海面測高是其主要應(yīng)用方向之一。如何實現(xiàn)5cm的海面測高精度,一直是GNSS-R研究的重點和難點。

    GNSS-R本質(zhì)上是雙基雷達(dá),GNSS衛(wèi)星發(fā)射的GNSS信號經(jīng)過地球表面反射后,被GNSS-R接收機(jī)接收和處理[2]。但GNSS信號經(jīng)過大氣、電離層和反射面的衰減和散射之后,信號及其微弱,為了能夠更好的捕獲反射信號,GNSS-R接收機(jī)往往只接收來自以鏡面反射點為中心的閃耀區(qū)的信號[3]。鏡面反射點預(yù)測的精確與否對于高精度海面測高應(yīng)用尤為重要,一方面,鏡面反射點是整個測高反演過程中的坐標(biāo)基準(zhǔn);另一方面,鏡面反射點的精度直接影響了反射信號的信噪比,進(jìn)而影響海面高度的反演精度。

    目前,常用的鏡面反射點預(yù)測算法有:Gleason算法[4]、S.C.Wu算法[5]、線段二分法[6]。其中Gleason算法利用向量共線的思想,算法效果不理想,精度差,迭代次數(shù)多;S.C.Wu算法精度高,但搜索完成之后需要進(jìn)行復(fù)雜的校正;線段二分法模型簡單,迭代次數(shù)少,但沒有解決徑向和法向不統(tǒng)一問題[7]。

    本文針對鏡面反射點預(yù)測問題,分析利用梯度下降法和AdaGrad算法預(yù)測鏡點存在的不足,提出了一種基于RMSprop的高精度鏡面反射點預(yù)測算法,該算法相比于現(xiàn)有的算法,通過平滑梯度平方,獲得更高的預(yù)測精度,適合對精度要求很高的GNSS-R海面測高應(yīng)用中。

    1 鏡面反射點定義

    GNSS-R技術(shù)中GNSS衛(wèi)星、地球表面和GNSS接收機(jī)的空間幾何位置關(guān)系如圖 1所示。在理想情況下,鏡面反射點滿足以下條件:① 鏡面反射點在地球表面;② 發(fā)射機(jī)-鏡面反射點-接收機(jī)是所有反射路徑中的最短路徑;③ 滿足斯涅耳定律——入射角等于反射角;④ 鏡面反射點與接收機(jī)、發(fā)射機(jī)和地心共面[8]。

    圖1 GNSS-R空間幾何位置關(guān)系Fig.1 Spatial geometric relations of GNSS-R

    圖2是GNSS空間幾何位置關(guān)系的二維視圖。

    圖2 幾何關(guān)系二維視圖Fig.2 2D view of geometric relations

    圖中,R表示接收機(jī),T表示發(fā)射機(jī),S表示鏡面反射點,O表示地心,M表示OC所在直線與RT所在直線的交點,U表示R關(guān)于OM的對稱點,C表示M關(guān)于RU的對稱點,αr、αt分別表示信號的入射角和反射角。

    根據(jù)上面對鏡面反射點條件的描述,鏡面反射點預(yù)測問題就是一個帶非線性約束條件的非線性優(yōu)化問題,數(shù)學(xué)表達(dá)式可以歸納為:

    (1)

    式中,f(S)是鏡面反射點預(yù)測的目標(biāo)函數(shù);R表示接收機(jī);T表示發(fā)射機(jī);S表示地球表面的鏡面反射點;Sx、Sy、Sz為S在地心坐標(biāo)系下對應(yīng)的三軸坐標(biāo);Rx、Ry、Rz為R在地心坐標(biāo)系下對應(yīng)的三軸坐標(biāo);Tx、Ty、Tz為T在地心坐標(biāo)系下對應(yīng)的三軸坐標(biāo);Ox、Oy、Oz為地心坐標(biāo)系原點的三軸坐標(biāo);αt、αr分別表示信號的入射角和反射角。

    2 基于RMSprop的鏡面反射點預(yù)測算法原理

    2.1 梯度下降法及其問題

    梯度下降法就是利用梯度的方向是函數(shù)變化最快的方向這一準(zhǔn)則進(jìn)行最優(yōu)化求解的[9][10]。針對反射路徑長度f(S),也可以用梯度下降法來實現(xiàn)對鏡面反射點的預(yù)測:

    f(S)=|S-T|+|R-S|。

    (2)

    即信號從GNSS衛(wèi)星發(fā)射經(jīng)鏡面反射點反射后到達(dá)GNSS-R接收機(jī)所走過的路徑長度。

    梯度下降法預(yù)測鏡點的具體求解流程如圖 3所示。首先求得一個粗略的鏡面反射點估計值S1,然后求在該點處反射路徑長度函數(shù)f(S)的梯度,記為f(S1),隨后利用f(S1)更新鏡面反射點位置S2,更新方程為:

    S2=S1-K×f(S1),

    (3)

    圖3 梯度下降法示意Fig.3 Schematic diagram of gradient descent method

    式中,K為學(xué)習(xí)率。

    再計算S2處f(S)的梯度f(S2),更新S3的位置……,對于第n次鏡面反射點預(yù)測,Sn更新表達(dá)式為:

    Sn=Sn-1-K×f(Sn)。

    (4)

    梯度下降法會一直計算到前后2個預(yù)測值之間的距離小于某個閾值,則認(rèn)為到達(dá)鏡面反射點,算法結(jié)束。

    在梯度下降迭代的過程中,學(xué)習(xí)率決定了每一步沿梯度負(fù)方向前進(jìn)的長度。梯度下降法采用了固定的學(xué)習(xí)率,并不會隨著預(yù)測的深入進(jìn)行自適應(yīng)調(diào)節(jié)。所以梯度下降法非常依賴學(xué)習(xí)率K的選擇,學(xué)習(xí)率太小,收斂慢,學(xué)習(xí)率太大,會導(dǎo)致算法在鏡面反射點附近發(fā)生震蕩,無法收斂[11]。

    2.2 Adagrad算法及其問題

    為了解決梯度下降法學(xué)習(xí)率無法調(diào)節(jié)的問題,引入了Adagrad算法[12]。Adagrad算法利用累積梯度平方和對學(xué)習(xí)率進(jìn)行約束,即將式(4)的鏡面反射點更新表達(dá)式改為:

    (5)

    隨著算法遞推的深入,梯度平方和越來越大,每次移動的步長越來越小,最后到達(dá)預(yù)測值。但是,在Adagrad算法執(zhí)行的中后期,分母上梯度平方和的累加將會越來越大,使得每次更新步長趨近于0,預(yù)測會提前結(jié)束。

    2.3 基于RMSprop的鏡面反射點預(yù)測算法

    為了解決Adagrad算法中學(xué)習(xí)率下降過快的問題,需要降低式(5)分母的平方和項,RMSprop通過將平方和變?yōu)橐苿蛹訖?quán)平均平方和來解決,算法鏡面反射點的更新表達(dá)式為:

    (6)

    式中,τn=τn-1×β+(1-β)×f(Sn)2;τn是加權(quán)平均梯度和;β是加權(quán)平均指數(shù)。

    式(6)中,τn完成的是對梯度的平方的一次平滑處理。在更新鏡面反射點Sn+1時,先用梯度f(Sn)除以的值,相當(dāng)于對梯度做了一次歸一化。在鏡面反射點預(yù)測的過程中,豎直方向的梯度比水平方向上的梯度大很多,此時豎直方向的τn也很大,當(dāng)f(Sn)除以之后,歸一化梯度變小,豎直方向的步長變小了;而水平方向梯度小,相應(yīng)的τn也小,水平方向上歸一化的梯度就大了,步長也變大了。

    基于RMSprop的鏡面反射點預(yù)測算法作為梯度下降法的優(yōu)化算法,一方面能夠解決Adagrad算法遞推過程學(xué)習(xí)率下降過快的問題,另一方面可以調(diào)整不同維度上的步長,加快收斂速度,最終實現(xiàn)對鏡面反射點的預(yù)測。

    3 基于RMSprop的鏡面反射點預(yù)測算法流程

    根據(jù)GNSS信號中攜帶的星歷數(shù)據(jù),可以解算出在WGS-84坐標(biāo)系統(tǒng)中接收機(jī)RWGS-84和發(fā)射機(jī)TWGS-84的位置,作為已知條件利用如下步驟進(jìn)行鏡面反射點的預(yù)測。

    基于RMSprop的高精度GNSS-R鏡面反射點預(yù)測算法流程如圖 4所示,其處理流程主要包括準(zhǔn)球映射、鏡面反射點粗估和鏡面反射點精估3個內(nèi)容。

    圖4 算法總流程Fig.4 SOverall flow of algorithm

    3.1 準(zhǔn)球映射

    將WGS-84坐標(biāo)系統(tǒng)映射成單位圓,WGS-84坐標(biāo)系統(tǒng)下的接收機(jī)位置RWGS-84和發(fā)射機(jī)位置TWGS-84將分別被映射成R和T。

    3.2 鏡面反射點的粗估

    根據(jù)如圖 2所示的空間幾何位置關(guān)系,可以得到如下數(shù)學(xué)關(guān)系:

    |RM|=|MU|=|UC|,

    (7)

    |RM|/|RT|=|RM|/(|RM|+|MT|)=|CU|/(|CU|+|MT|)=

    |SR|/(|SR|+|ST|)≈|OR|/(|OR|+|OT|)。

    (8)

    M點的估計值為:

    M=R+RM≈R+[|OR|/(|OR|+|OT|)]×RT。

    (9)

    S1是M的星下點,則:

    S1=(OM/|OM|)×Re,

    (10)

    式中,Re為地球半徑,因為進(jìn)行了準(zhǔn)球映射,此時地球是單位圓,即Re=1。

    3.3 RMSprop鏡面反射點精估

    設(shè)置全局學(xué)習(xí)率ε,鏡面反射點的初始值S1,小常數(shù)δ(典型值10-8,用于被小數(shù)除時的數(shù)值穩(wěn)定),加權(quán)平均系數(shù)β。

    ① 初始化梯度累積變量τ1=0;

    f(Sn)=(Sn-T)/|T-Sn|+(Sn-R)/|R-Sn|,

    (11)

    ③ 計算累積梯度平方和τn:

    τn=τn-1×β+(1-β)×f(Sn)2,

    (12)

    (13)

    ④ 計算更新增量ΔS;

    ⑤ 更新沿梯度方向下降后的預(yù)測值S′:

    S′=Sn+ΔS,

    (14)

    ⑥ 取S′的星下點作為當(dāng)前鏡面反射點預(yù)測值Sn+1:

    Sn+1=(OS′/|OS′|)×Re(S′),

    (15)

    ⑦ 計算當(dāng)前鏡面反射點預(yù)測值與之前的歐氏距離作為更新的幅度,如果大于預(yù)設(shè)門限則返回步驟①,反之則退出循環(huán),以當(dāng)前鏡面反射點預(yù)測值Sn+1作為鏡面反射點最優(yōu)值進(jìn)行反射信號的捕獲跟蹤。

    AMP=|Sn+1-Sn|。

    (16)

    4 仿真與分析

    本文利用Matlab和STK仿真工具進(jìn)行如圖 5所示的聯(lián)合仿真。利用Matlab對STK進(jìn)行軌道參數(shù)配置后,STK完成衛(wèi)星可見性分析,將分析結(jié)果導(dǎo)入Matlab,在Matlab中執(zhí)行鏡面反射點計算并分析算法性能。本次仿真實驗中,GNSS星座采用GPS星座和北斗星座共58顆GNSS衛(wèi)星,LEO衛(wèi)星的軌道高度為500 km,接收入射角最大45°的反射信號,仿真了在軌運行24 h(2019年3月17日—18日)的數(shù)據(jù),采樣時間間隔為1 s。利用GPS的PRN1衛(wèi)星與LEO所得到的9段數(shù)據(jù)共12 341個時間點進(jìn)行下面的仿真。

    圖5 聯(lián)合仿真結(jié)構(gòu)框架Fig.5 Co-simulation framework

    4.1 粗估的S1的精度

    利用本文算法預(yù)測每個時刻的鏡面反射點,將最終預(yù)測結(jié)果作為鏡面反射點預(yù)測值,把粗估的鏡面反射點值與預(yù)測值進(jìn)行比較,分析粗估的精度。

    鏡面反射點粗估性能如圖6所示。

    (a) 鏡面反射點粗估誤差

    (b) 各時刻GNSS衛(wèi)星與LEO衛(wèi)星距離圖6 鏡面反射點粗估性能Fig.6 Performance of specular point by coarse estimation

    圖6(a)中,粗估值與預(yù)測值相距最近65 m,最遠(yuǎn)172 320 m,均值68 290 m。圖 6(b)是各個時刻GNSS衛(wèi)星與LEO衛(wèi)星的距離。圖6(a)中曲線的起伏與圖6(b)中曲線的起伏基本一致,說明粗估的精度與GNSS衛(wèi)星和LEO衛(wèi)星所處的幾何位置有關(guān),二者距離越近,衛(wèi)星高度角就越高,粗估的精度就越高;二者距離越遠(yuǎn),衛(wèi)星高度角就越低,粗估的精度就越低。

    4.2 單個數(shù)據(jù)仿真性能比較

    選取了上述12 341個時間點內(nèi)某一時間點GNSS衛(wèi)星和LEO衛(wèi)星的位置,分別利用Gleason算法、S.C.Wu算法、線段二分法和本文算法進(jìn)行鏡面反射點預(yù)測,分析其收斂性能和算法精度。收斂性能的評判依據(jù)是算法的迭代次數(shù),算法精度的依據(jù)是GNSS信號從GNSS衛(wèi)星發(fā)射經(jīng)鏡面反射點反射到達(dá)LEO衛(wèi)星所走過的路徑長度。

    4.2.1 算法精度

    圖 7是4種算法路徑長度的比較。從圖中可以看出,線段二分法、S.C.Wu算法和Gleason算法得到的反射路徑長度非常接近。其中,S.C.Wu算法得到的反射路徑比線段二分法和Gleanson算法得到的反射路徑長度更小。而本文提出的算法得到的反射路徑長度最短。這說明本文提出的算法獲得的鏡面反射點與真實鏡面反射點距離最近,計算精度最高。

    圖7 4種算法得到的路徑長度Fig.7 Path length obtained by four algorithms

    4.2.2 收斂性能

    將仿真中本文算法計算得到的鏡面反射點作為參考,對算法的收斂性能進(jìn)行分析。圖 8為4種鏡面反射點算法收斂性能的比較圖。通過圖 8分析可以得到,線段二分法在大約10次迭代后收斂至最小點,其曲線下降過程比較震蕩;S.C.Wu算法在大約10次迭代后,曲線達(dá)到最小值,隨后曲線迅速上升且無法收斂,該現(xiàn)象是由未經(jīng)過坐標(biāo)校正引起的;Gleason算法在收斂過程很慢,曲線較為平滑,經(jīng)過大約100次迭代后,曲線開始收斂,其迭代次數(shù)是最多的;本文提出的算法在經(jīng)過約10次迭代后達(dá)到預(yù)測點,迭代次數(shù)最少,運算速度最快。

    圖8 4種算法的收斂性能Fig.8 Convergence performance of four algorithms

    4.3 統(tǒng)計數(shù)據(jù)仿真性能比較

    將12 341個時間點每個時刻GNSS衛(wèi)星和LEO衛(wèi)星的位置分別利用線段二分法和基于RMSprop的鏡面反射點預(yù)測算法進(jìn)行預(yù)測,將結(jié)果進(jìn)行統(tǒng)計意義上的收斂性能比較和精度比較。

    4.3.1 收斂性能

    圖 9是本文算法和線段二分法的效率統(tǒng)計圖。其中橫軸表示對單個數(shù)據(jù)點進(jìn)行計算,收斂到預(yù)測值所需的迭代次數(shù)??v軸表示需要相應(yīng)迭代次數(shù)的數(shù)據(jù)點個數(shù)。通過比較可以得到,本文算法中所有數(shù)據(jù)點計算的迭代次數(shù)區(qū)間為9~15,而線段二分法中所有數(shù)據(jù)點計算的迭代次數(shù)區(qū)間為9~19,且數(shù)據(jù)點落在17~18區(qū)間的最多。利用均值和方差對本文提出的算法和線段二分法進(jìn)行統(tǒng)計意義分析,得到本文算法與線段二分法的收斂性能統(tǒng)計表,如表1所示。

    表1 本文算法與線段二分法的收斂性能統(tǒng)計表Tab.1 Convergence performance list of the proposed algorithm and the dichotomy of line segment

    (a) RMSprop算法的迭代次數(shù)統(tǒng)計

    (b) 線段二分法的迭代次數(shù)統(tǒng)計圖9 2種算法的迭代次數(shù)直方圖Fig.9 Iteration times histogram of two algorithms

    分析表 1可以得到,本文提出的算法運算效率略優(yōu)于線段二分法。

    4.3.2 算法精度

    圖10是本文算法和線段二分法的路徑長度比較圖。

    (a) 本文算法和線段二分法路徑長度比較

    (b) 2種算法路徑長度差值圖10 2種算法的路徑長度比較Fig.10 Path length comparison between two algorithms

    圖10(a)中,2種算法的路徑長度曲線基本重合,其中RMSprop算法的路徑長度要小一點,為了進(jìn)行更直觀的分析,本文將線段二分法路徑長度與RMSprop算法路徑長度進(jìn)行差值,得到如圖10(b)所示的路徑長度差值圖,可以看出,RMSprop算法要比線段二分法算出來的路徑短18km左右。經(jīng)統(tǒng)計分析可以得到,在觀察數(shù)據(jù)點樣本中,任何時刻,本文算法的精度都更高。

    綜上所述,無論是算法精度還是收斂性能,本文算法都明顯優(yōu)于傳統(tǒng)的Gleason算法、S.C.Wu算法和線段二分法。

    5 結(jié)束語

    本文通過研究鏡面反射點定義、梯度下降法及其優(yōu)化算法的相關(guān)知識,提出了一種基于RMSprop的鏡面反射點預(yù)測算法。作為機(jī)器學(xué)習(xí)中梯度下降法優(yōu)化算法的一種,該算法解決了梯度下降法的收斂問題,也不會和AdaGrad算法一樣提前收斂,專門用于對鏡面反射點精確度要求很高的GNSS-R海面測高應(yīng)用中。為了驗證算法的正確性,利用Matlab和STK進(jìn)行聯(lián)合仿真,通過與常用算法的比較,結(jié)果表明,本文提出的算法可以對鏡面反射點進(jìn)行正確預(yù)測,該算法不僅精度高,而且收斂性能好,運算效率高,適用于GNSS-R海面測高應(yīng)用。

    猜你喜歡
    鏡面反射二分法線段
    光滑物體表面反射光偏振特征分析及反射光分離技術(shù)*
    基于二進(jìn)制/二分法的ETC狀態(tài)名單查找算法
    “二分法”求解加速度的分析策略
    “二分法”求解加速度的分析策略
    畫出線段圖來比較
    基于最短路徑的GNSS-R鏡面反射點算法
    怎樣畫線段圖
    我們一起數(shù)線段
    數(shù)線段
    估算的妙招——“二分法”
    日韩制服骚丝袜av| 国产午夜精品一二区理论片| 两个人的视频大全免费| 我的女老师完整版在线观看| 好男人视频免费观看在线| 女人被狂操c到高潮| 99热精品在线国产| 亚洲一区二区三区色噜噜| 婷婷六月久久综合丁香| 色播亚洲综合网| 午夜福利视频1000在线观看| 国模一区二区三区四区视频| 国产成人a∨麻豆精品| 欧美成人一区二区免费高清观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲无线观看免费| 亚洲精品日韩在线中文字幕 | 久久久久性生活片| 尾随美女入室| 18禁黄网站禁片免费观看直播| 嫩草影院精品99| 熟女人妻精品中文字幕| 成人午夜精彩视频在线观看| 亚洲av中文av极速乱| 少妇的逼好多水| 免费观看的影片在线观看| 亚洲不卡免费看| 亚洲国产精品合色在线| 男人舔女人下体高潮全视频| 一边摸一边抽搐一进一小说| 国内精品久久久久精免费| 国产高清不卡午夜福利| 久久久久久久久大av| 99久久无色码亚洲精品果冻| 又爽又黄a免费视频| 人妻制服诱惑在线中文字幕| 少妇被粗大猛烈的视频| 欧美成人一区二区免费高清观看| 国产成人精品一,二区 | 91精品国产九色| 尾随美女入室| 日韩欧美精品v在线| 欧美性猛交╳xxx乱大交人| 少妇人妻精品综合一区二区 | 久久精品国产清高在天天线| 亚洲av免费高清在线观看| 国产综合懂色| 日本黄色视频三级网站网址| 亚洲精品粉嫩美女一区| 日本撒尿小便嘘嘘汇集6| 人人妻人人看人人澡| 日本免费一区二区三区高清不卡| 国产一区二区亚洲精品在线观看| 麻豆乱淫一区二区| 色5月婷婷丁香| av在线观看视频网站免费| 一级毛片电影观看 | 中文字幕制服av| 亚洲欧美日韩东京热| 深夜a级毛片| av在线天堂中文字幕| 日本撒尿小便嘘嘘汇集6| 97超碰精品成人国产| 国产精品综合久久久久久久免费| 狂野欧美激情性xxxx在线观看| 亚洲欧美中文字幕日韩二区| 亚洲第一电影网av| 中文字幕av在线有码专区| 国产精品1区2区在线观看.| 亚洲第一区二区三区不卡| 免费电影在线观看免费观看| av女优亚洲男人天堂| 欧美最黄视频在线播放免费| 2022亚洲国产成人精品| 岛国在线免费视频观看| 国产在视频线在精品| 在线a可以看的网站| 91麻豆精品激情在线观看国产| 一边摸一边抽搐一进一小说| 观看免费一级毛片| 又粗又爽又猛毛片免费看| 综合色av麻豆| 欧美丝袜亚洲另类| 亚洲在线观看片| 丰满的人妻完整版| 毛片女人毛片| 在线国产一区二区在线| 国产 一区精品| 91在线精品国自产拍蜜月| 成熟少妇高潮喷水视频| 可以在线观看毛片的网站| 亚洲人成网站在线播放欧美日韩| 一边亲一边摸免费视频| 国产精品一二三区在线看| 欧美+日韩+精品| 少妇人妻一区二区三区视频| 嫩草影院入口| 啦啦啦观看免费观看视频高清| 国产精品综合久久久久久久免费| 国产黄色小视频在线观看| 91精品国产九色| 免费无遮挡裸体视频| 丰满的人妻完整版| 少妇人妻精品综合一区二区 | 国产亚洲91精品色在线| 青春草国产在线视频 | 成年女人永久免费观看视频| 99热全是精品| h日本视频在线播放| 亚洲av男天堂| 久久精品夜色国产| 国产精品99久久久久久久久| 国产精品久久久久久精品电影| 国产成人影院久久av| 啦啦啦观看免费观看视频高清| 日韩av不卡免费在线播放| 成人特级av手机在线观看| 少妇熟女欧美另类| 国产黄片美女视频| 深夜精品福利| 亚洲一区高清亚洲精品| 免费av不卡在线播放| 久久久成人免费电影| 嫩草影院入口| 可以在线观看毛片的网站| 久久99精品国语久久久| 日韩高清综合在线| 18+在线观看网站| 欧美日韩在线观看h| 天美传媒精品一区二区| 级片在线观看| 男人舔奶头视频| 我要搜黄色片| 国产亚洲av片在线观看秒播厂 | 久久国内精品自在自线图片| 国产午夜精品论理片| 国产极品精品免费视频能看的| www.色视频.com| 麻豆国产av国片精品| 欧美最黄视频在线播放免费| 在线天堂最新版资源| 欧美日韩一区二区视频在线观看视频在线 | 国产久久久一区二区三区| 成人永久免费在线观看视频| 精品一区二区三区人妻视频| 久久精品国产亚洲av香蕉五月| 欧美日韩乱码在线| 午夜激情福利司机影院| 国产色婷婷99| 亚洲经典国产精华液单| 91aial.com中文字幕在线观看| 少妇的逼好多水| 亚洲va在线va天堂va国产| 亚洲欧美日韩高清在线视频| 国产午夜精品久久久久久一区二区三区| 99国产精品一区二区蜜桃av| 成人综合一区亚洲| 亚洲电影在线观看av| 日韩一本色道免费dvd| 亚洲人成网站在线播| 亚洲人成网站高清观看| 麻豆成人av视频| 成人毛片60女人毛片免费| 国产男人的电影天堂91| 国产淫片久久久久久久久| 看片在线看免费视频| 99国产极品粉嫩在线观看| 18禁在线播放成人免费| 国产一级毛片在线| 欧美日韩国产亚洲二区| 亚洲电影在线观看av| 国产不卡一卡二| 国产精品久久久久久av不卡| 中出人妻视频一区二区| 久久这里有精品视频免费| 成年女人看的毛片在线观看| 日韩欧美一区二区三区在线观看| 嫩草影院精品99| 热99re8久久精品国产| 国产午夜精品一二区理论片| 丰满人妻一区二区三区视频av| 午夜激情欧美在线| 久久久久久久久中文| 日韩精品有码人妻一区| 亚洲国产色片| 老熟妇乱子伦视频在线观看| 青春草亚洲视频在线观看| 校园春色视频在线观看| 12—13女人毛片做爰片一| 在线观看美女被高潮喷水网站| 老女人水多毛片| 波多野结衣高清无吗| 亚洲国产精品成人综合色| 国产日韩欧美在线精品| 亚洲欧洲国产日韩| 久久精品国产亚洲网站| 天堂中文最新版在线下载 | 成人亚洲精品av一区二区| 亚洲第一区二区三区不卡| 美女xxoo啪啪120秒动态图| 国产精品永久免费网站| 日韩精品青青久久久久久| 卡戴珊不雅视频在线播放| 亚洲三级黄色毛片| 国产蜜桃级精品一区二区三区| 久久精品夜色国产| 国产一级毛片七仙女欲春2| 欧美日韩乱码在线| 成人欧美大片| 国产精品一及| 性欧美人与动物交配| 亚洲国产高清在线一区二区三| 中国国产av一级| 国产 一区 欧美 日韩| 日日啪夜夜撸| 久久久久网色| 日韩一区二区三区影片| 熟女电影av网| 最近视频中文字幕2019在线8| 午夜久久久久精精品| 神马国产精品三级电影在线观看| 久久国内精品自在自线图片| 亚洲av二区三区四区| 九九热线精品视视频播放| 直男gayav资源| 亚洲av中文av极速乱| 婷婷色av中文字幕| 亚洲色图av天堂| 成年女人永久免费观看视频| av在线老鸭窝| 日韩欧美精品v在线| 国产精品嫩草影院av在线观看| 搞女人的毛片| 国产一区二区三区在线臀色熟女| 青春草视频在线免费观看| 色噜噜av男人的天堂激情| 亚洲成人久久爱视频| 在线观看66精品国产| 国产精品爽爽va在线观看网站| 天堂影院成人在线观看| 欧美日韩精品成人综合77777| 日本一二三区视频观看| 白带黄色成豆腐渣| 综合色丁香网| 久久亚洲国产成人精品v| 亚州av有码| 特大巨黑吊av在线直播| 黄片无遮挡物在线观看| 51国产日韩欧美| 变态另类成人亚洲欧美熟女| 一夜夜www| 一级毛片aaaaaa免费看小| av免费观看日本| 给我免费播放毛片高清在线观看| 欧美+日韩+精品| 国产男人的电影天堂91| 又爽又黄a免费视频| 免费看a级黄色片| 国产亚洲av嫩草精品影院| 在线播放无遮挡| 高清日韩中文字幕在线| 精品日产1卡2卡| 国内精品美女久久久久久| 亚洲av第一区精品v没综合| 午夜精品一区二区三区免费看| 日本与韩国留学比较| 男人的好看免费观看在线视频| 久久久久久久久久成人| 国产一区二区亚洲精品在线观看| 国产亚洲av片在线观看秒播厂 | 欧美成人免费av一区二区三区| .国产精品久久| 国产视频首页在线观看| 午夜精品一区二区三区免费看| 欧美日韩精品成人综合77777| 久99久视频精品免费| 看十八女毛片水多多多| 在线观看美女被高潮喷水网站| 国产av一区在线观看免费| 在线观看免费视频日本深夜| 国产精品无大码| 变态另类成人亚洲欧美熟女| 亚洲欧美精品专区久久| 综合色av麻豆| 国产黄色视频一区二区在线观看 | 日本撒尿小便嘘嘘汇集6| 久久亚洲国产成人精品v| 人人妻人人澡欧美一区二区| 青青草视频在线视频观看| 成人鲁丝片一二三区免费| 国产成人一区二区在线| 亚洲国产精品sss在线观看| 校园人妻丝袜中文字幕| 国产v大片淫在线免费观看| 伦理电影大哥的女人| 色综合亚洲欧美另类图片| 久久久久性生活片| 国产精品一区二区性色av| 色播亚洲综合网| 日韩欧美在线乱码| 国产一区亚洲一区在线观看| 亚洲国产精品国产精品| av国产免费在线观看| 久久久久网色| 亚洲人成网站在线播放欧美日韩| 六月丁香七月| 久久久久免费精品人妻一区二区| 亚州av有码| 极品教师在线视频| 波多野结衣高清无吗| 久久久久九九精品影院| 色综合色国产| 欧美激情在线99| 日韩欧美在线乱码| 亚洲精品色激情综合| 欧美另类亚洲清纯唯美| 精品国产三级普通话版| 国产日本99.免费观看| 身体一侧抽搐| 日本一二三区视频观看| 国产极品精品免费视频能看的| 18禁在线无遮挡免费观看视频| 特级一级黄色大片| 超碰av人人做人人爽久久| 国产av麻豆久久久久久久| 日日撸夜夜添| 99久久成人亚洲精品观看| 五月伊人婷婷丁香| 欧美成人精品欧美一级黄| 男女啪啪激烈高潮av片| 国产精品久久久久久精品电影小说 | 婷婷精品国产亚洲av| 综合色丁香网| 国产伦一二天堂av在线观看| 我要看日韩黄色一级片| 97超碰精品成人国产| 欧美精品国产亚洲| 尤物成人国产欧美一区二区三区| 国产精品av视频在线免费观看| 午夜福利成人在线免费观看| 成人三级黄色视频| 夜夜夜夜夜久久久久| ponron亚洲| 在线免费十八禁| 亚洲电影在线观看av| 国产精品久久视频播放| 黄色日韩在线| 一级毛片我不卡| 国产精品一区二区三区四区免费观看| 日韩成人av中文字幕在线观看| 久久久欧美国产精品| 变态另类丝袜制服| 国产日韩欧美在线精品| 国产一区亚洲一区在线观看| 黄色欧美视频在线观看| 国产一区二区激情短视频| 十八禁国产超污无遮挡网站| 99热只有精品国产| 午夜精品一区二区三区免费看| 91麻豆精品激情在线观看国产| 亚洲美女搞黄在线观看| 男人舔女人下体高潮全视频| 久久久久久久久久久丰满| 欧美bdsm另类| 内地一区二区视频在线| 久久精品国产清高在天天线| 真实男女啪啪啪动态图| 日韩三级伦理在线观看| 国产欧美日韩精品一区二区| 色综合站精品国产| av在线亚洲专区| 青春草亚洲视频在线观看| 91精品一卡2卡3卡4卡| 亚洲精品国产成人久久av| 午夜免费男女啪啪视频观看| 淫秽高清视频在线观看| 国产伦在线观看视频一区| 欧美性感艳星| 亚洲在久久综合| 亚洲七黄色美女视频| 内射极品少妇av片p| 青青草视频在线视频观看| 中国美白少妇内射xxxbb| a级毛片a级免费在线| 小蜜桃在线观看免费完整版高清| 激情 狠狠 欧美| 日韩在线高清观看一区二区三区| 亚洲欧洲日产国产| 亚洲人成网站在线观看播放| 少妇裸体淫交视频免费看高清| 91狼人影院| 一级毛片久久久久久久久女| 成人永久免费在线观看视频| av视频在线观看入口| 一个人看视频在线观看www免费| 国产精品女同一区二区软件| 天天一区二区日本电影三级| 欧美在线一区亚洲| 亚洲色图av天堂| 欧美xxxx黑人xx丫x性爽| 看免费成人av毛片| 亚洲国产精品成人综合色| 人人妻人人澡人人爽人人夜夜 | 亚洲一级一片aⅴ在线观看| 国产亚洲av片在线观看秒播厂 | 亚洲中文字幕一区二区三区有码在线看| 午夜福利视频1000在线观看| 国产高潮美女av| 麻豆成人午夜福利视频| 国产乱人视频| 韩国av在线不卡| 亚洲欧美精品专区久久| 亚洲欧美精品综合久久99| 欧美bdsm另类| 国产色爽女视频免费观看| 久久久久久久久久黄片| 成人欧美大片| 一级毛片久久久久久久久女| 久久鲁丝午夜福利片| 69av精品久久久久久| 白带黄色成豆腐渣| 亚洲精品久久国产高清桃花| 亚洲第一区二区三区不卡| 简卡轻食公司| 99国产极品粉嫩在线观看| 麻豆av噜噜一区二区三区| 3wmmmm亚洲av在线观看| 久久精品国产鲁丝片午夜精品| 久久久成人免费电影| 精品久久久久久久久久免费视频| 一边亲一边摸免费视频| 亚洲人成网站在线观看播放| 日日干狠狠操夜夜爽| 久久久久久久久久久免费av| 亚洲欧美清纯卡通| 1000部很黄的大片| 少妇人妻精品综合一区二区 | av女优亚洲男人天堂| 欧美+亚洲+日韩+国产| 麻豆精品久久久久久蜜桃| 精品人妻视频免费看| 中文字幕免费在线视频6| 最新中文字幕久久久久| 99久久成人亚洲精品观看| 亚洲18禁久久av| 可以在线观看毛片的网站| 久久精品综合一区二区三区| 内射极品少妇av片p| 亚洲国产精品成人综合色| 草草在线视频免费看| 男女视频在线观看网站免费| 亚洲无线观看免费| 国产一区二区激情短视频| 欧洲精品卡2卡3卡4卡5卡区| 久久欧美精品欧美久久欧美| 美女cb高潮喷水在线观看| 免费看光身美女| 美女xxoo啪啪120秒动态图| 亚洲在线自拍视频| 久久精品久久久久久噜噜老黄 | 特大巨黑吊av在线直播| av天堂在线播放| 精品欧美国产一区二区三| 看十八女毛片水多多多| а√天堂www在线а√下载| 两性午夜刺激爽爽歪歪视频在线观看| 一个人观看的视频www高清免费观看| 美女脱内裤让男人舔精品视频 | 亚洲中文字幕日韩| 亚洲成人久久爱视频| 久久国产乱子免费精品| 久久精品国产亚洲av天美| 爱豆传媒免费全集在线观看| 赤兔流量卡办理| 国产高清有码在线观看视频| 在线天堂最新版资源| 最近的中文字幕免费完整| 精品人妻视频免费看| 级片在线观看| 97在线视频观看| 不卡视频在线观看欧美| 一边摸一边抽搐一进一小说| 亚洲激情五月婷婷啪啪| 亚洲欧美精品综合久久99| 久久久久久伊人网av| 免费一级毛片在线播放高清视频| 一级毛片久久久久久久久女| 国产熟女欧美一区二区| 三级毛片av免费| 91aial.com中文字幕在线观看| 国产乱人偷精品视频| 99在线人妻在线中文字幕| 亚洲成av人片在线播放无| 亚洲aⅴ乱码一区二区在线播放| 日本成人三级电影网站| 熟女人妻精品中文字幕| 男女下面进入的视频免费午夜| 国产精品麻豆人妻色哟哟久久 | 国产精品久久久久久亚洲av鲁大| 国内精品宾馆在线| 欧美三级亚洲精品| 国产真实乱freesex| 亚洲自拍偷在线| 中文在线观看免费www的网站| 精品久久国产蜜桃| 一个人看视频在线观看www免费| 亚洲欧美清纯卡通| 午夜福利在线在线| 精品久久久久久久久久久久久| 国产成人福利小说| 国产毛片a区久久久久| 亚洲欧美清纯卡通| 蜜臀久久99精品久久宅男| 亚洲图色成人| 国产精品三级大全| 婷婷亚洲欧美| 国产淫片久久久久久久久| 哪个播放器可以免费观看大片| 国产综合懂色| 青青草视频在线视频观看| 小说图片视频综合网站| 99在线视频只有这里精品首页| av免费观看日本| 国产精品麻豆人妻色哟哟久久 | 色哟哟·www| 天天一区二区日本电影三级| 久久九九热精品免费| 国产精品无大码| 国产亚洲av片在线观看秒播厂 | 五月伊人婷婷丁香| 久久久精品大字幕| 校园人妻丝袜中文字幕| 我的老师免费观看完整版| 久久久国产成人精品二区| 亚州av有码| 国产午夜精品一二区理论片| 免费大片18禁| 熟妇人妻久久中文字幕3abv| 嫩草影院入口| 国产精品久久久久久亚洲av鲁大| 久久99热6这里只有精品| av免费观看日本| 亚洲国产精品久久男人天堂| 久久精品综合一区二区三区| av天堂中文字幕网| 在线免费观看不下载黄p国产| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品1区2区在线观看.| 午夜爱爱视频在线播放| 男女下面进入的视频免费午夜| 国产精品野战在线观看| 亚洲七黄色美女视频| 亚洲中文字幕一区二区三区有码在线看| 久久精品国产99精品国产亚洲性色| 丰满人妻一区二区三区视频av| 伦精品一区二区三区| 日本一二三区视频观看| 亚洲精品日韩av片在线观看| 国产在线精品亚洲第一网站| 国产高清有码在线观看视频| 国产高潮美女av| 欧洲精品卡2卡3卡4卡5卡区| 久久久精品94久久精品| 一级黄色大片毛片| 欧美色欧美亚洲另类二区| www.色视频.com| 亚洲成av人片在线播放无| 两个人视频免费观看高清| 国产一级毛片七仙女欲春2| 国产精品蜜桃在线观看 | 观看免费一级毛片| 免费看美女性在线毛片视频| 天堂中文最新版在线下载 | 99国产极品粉嫩在线观看| 天堂影院成人在线观看| 18禁在线播放成人免费| 寂寞人妻少妇视频99o| 国产v大片淫在线免费观看| 97超视频在线观看视频| 久久久久久久久久成人| 亚洲欧美中文字幕日韩二区| 亚洲欧洲国产日韩| av在线老鸭窝| 淫秽高清视频在线观看| 日本一本二区三区精品| 特级一级黄色大片| 校园人妻丝袜中文字幕| 大香蕉久久网| 国产精品久久久久久av不卡| 波野结衣二区三区在线| 国产精品久久视频播放| 亚洲第一区二区三区不卡| 好男人在线观看高清免费视频| 国产女主播在线喷水免费视频网站 | 国产在线精品亚洲第一网站| 国产一区二区激情短视频| 国产伦理片在线播放av一区 | 国产精品永久免费网站| 天堂√8在线中文| 久久久欧美国产精品| 久久6这里有精品| 欧美区成人在线视频| 麻豆国产av国片精品| 午夜福利在线在线| 中文字幕免费在线视频6| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费看a级黄色片| 又爽又黄a免费视频| 在线观看美女被高潮喷水网站| av在线蜜桃| 久久精品国产亚洲av天美| 99九九线精品视频在线观看视频| 成人毛片60女人毛片免费| 久久精品国产亚洲av天美| 观看美女的网站|