• <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福利一区| 日产精品乱码卡一卡2卡三| 欧美+日韩+精品| 欧美丝袜亚洲另类| 偷拍熟女少妇极品色| 国产高潮美女av| 午夜免费观看性视频| 又爽又黄a免费视频| 亚洲国产色片| 91久久精品国产一区二区成人| 黑人高潮一二区| 国产伦精品一区二区三区视频9| 中文字幕av成人在线电影| 久久鲁丝午夜福利片| 美女被艹到高潮喷水动态| 女人被狂操c到高潮| 黄色视频在线播放观看不卡| 免费观看性生交大片5| 成人特级av手机在线观看| 国产成年人精品一区二区| 五月伊人婷婷丁香| 亚洲av国产av综合av卡| 国产老妇女一区| 国产淫语在线视频| 亚洲久久久久久中文字幕| 麻豆成人午夜福利视频| 亚洲婷婷狠狠爱综合网| 久久久精品欧美日韩精品| 免费电影在线观看免费观看| 亚洲欧美日韩无卡精品| 日韩人妻高清精品专区| 国产高清不卡午夜福利| 国产成人freesex在线| 久久久欧美国产精品| 日本一本二区三区精品| 尤物成人国产欧美一区二区三区| 久久99热6这里只有精品| 成人无遮挡网站| 国产黄片视频在线免费观看| 爱豆传媒免费全集在线观看| 少妇人妻一区二区三区视频| 天天躁夜夜躁狠狠久久av| 天天一区二区日本电影三级| 国产伦理片在线播放av一区| 亚洲精品一区蜜桃| 一级av片app| 秋霞在线观看毛片| 青青草视频在线视频观看| 精华霜和精华液先用哪个| 色综合色国产| 国产淫语在线视频| 直男gayav资源| 大香蕉97超碰在线| a级毛片免费高清观看在线播放| 国产精品秋霞免费鲁丝片| 国产爱豆传媒在线观看| 亚洲电影在线观看av| 人人妻人人看人人澡| 极品教师在线视频| 热99国产精品久久久久久7| 91在线精品国自产拍蜜月| 亚洲自偷自拍三级| 精华霜和精华液先用哪个| 成人亚洲欧美一区二区av| a级一级毛片免费在线观看| 一级爰片在线观看| 午夜精品国产一区二区电影 | 97精品久久久久久久久久精品| 王馨瑶露胸无遮挡在线观看| 久久影院123| 国产成人福利小说| 大片电影免费在线观看免费| 亚洲精品乱久久久久久| eeuss影院久久| 亚洲精品久久久久久婷婷小说| 亚洲在久久综合| 日韩欧美一区视频在线观看 | av专区在线播放| 狂野欧美激情性bbbbbb| 18禁裸乳无遮挡免费网站照片| 国产综合精华液| 亚洲精品国产成人久久av| 免费不卡的大黄色大毛片视频在线观看| 午夜福利在线在线| 国产 一区精品| 九九在线视频观看精品| 少妇的逼水好多| 亚洲av成人精品一区久久| 老女人水多毛片| 真实男女啪啪啪动态图| 九九爱精品视频在线观看| 国产 精品1| 一级毛片我不卡| 自拍偷自拍亚洲精品老妇| 成人亚洲精品一区在线观看 | 日本av手机在线免费观看| 免费av毛片视频| 久久久成人免费电影| 国产一区二区在线观看日韩| 少妇人妻一区二区三区视频| 亚洲欧美精品自产自拍| .国产精品久久| 久久久久久久久久久免费av| 不卡视频在线观看欧美| 人人妻人人看人人澡| 国产淫语在线视频| 免费黄频网站在线观看国产| 国产色爽女视频免费观看| 99热这里只有精品一区| 国产69精品久久久久777片| 一本久久精品| 99久国产av精品国产电影| 亚洲高清免费不卡视频| 大片免费播放器 马上看| 高清欧美精品videossex| 午夜福利视频精品| 91久久精品国产一区二区成人| 18禁在线播放成人免费| 欧美少妇被猛烈插入视频| 777米奇影视久久| 一级片'在线观看视频| 亚洲av福利一区| 免费不卡的大黄色大毛片视频在线观看| 国产伦精品一区二区三区视频9| 看非洲黑人一级黄片| 哪个播放器可以免费观看大片| 91aial.com中文字幕在线观看| 国产精品成人在线| 欧美日韩视频精品一区| 久久久久久九九精品二区国产| 久久久亚洲精品成人影院| 久久99精品国语久久久| 国产精品国产三级专区第一集| 91午夜精品亚洲一区二区三区| 男的添女的下面高潮视频| 午夜福利网站1000一区二区三区| 亚洲av日韩在线播放| 老司机影院毛片| 日本-黄色视频高清免费观看| 99精国产麻豆久久婷婷| 国产免费一级a男人的天堂| 国内精品美女久久久久久| 纵有疾风起免费观看全集完整版| 国产免费福利视频在线观看| 日韩一本色道免费dvd| 男女下面进入的视频免费午夜| 男女国产视频网站| 国产高清有码在线观看视频| 亚洲国产欧美人成| 久久久精品欧美日韩精品| 一级毛片久久久久久久久女| 免费大片18禁| 嫩草影院精品99| 综合色丁香网| 69人妻影院| 日日摸夜夜添夜夜爱| 一个人观看的视频www高清免费观看| 一区二区三区四区激情视频| 在线观看免费高清a一片| 日韩免费高清中文字幕av| 日本欧美国产在线视频| 国产中年淑女户外野战色| 欧美日韩亚洲高清精品| 精品人妻一区二区三区麻豆| 观看免费一级毛片| 成年女人在线观看亚洲视频 | 国产高潮美女av| 一个人看的www免费观看视频| 亚洲经典国产精华液单| 国产伦精品一区二区三区四那| 波多野结衣巨乳人妻| 小蜜桃在线观看免费完整版高清| 你懂的网址亚洲精品在线观看| 联通29元200g的流量卡| 中文资源天堂在线| 国产白丝娇喘喷水9色精品| 99精国产麻豆久久婷婷| 亚洲av欧美aⅴ国产| 天天躁夜夜躁狠狠久久av| 欧美精品人与动牲交sv欧美| 黄片wwwwww| 2018国产大陆天天弄谢| 丰满乱子伦码专区| 丰满少妇做爰视频| 人妻少妇偷人精品九色| 性色av一级| 亚洲国产精品国产精品| 大又大粗又爽又黄少妇毛片口| 久久这里有精品视频免费| 在线精品无人区一区二区三 | 简卡轻食公司| 青青草视频在线视频观看| 欧美激情在线99| 国内精品宾馆在线| 国产精品久久久久久精品电影小说 | 亚洲av不卡在线观看| 亚洲精品成人av观看孕妇| 少妇人妻精品综合一区二区| 亚洲精品中文字幕在线视频 | 欧美变态另类bdsm刘玥| 亚洲在久久综合| 在线a可以看的网站| 边亲边吃奶的免费视频| 亚洲真实伦在线观看| 日本午夜av视频| av免费在线看不卡| 国产精品爽爽va在线观看网站| 国产精品久久久久久av不卡| 日韩欧美精品免费久久| 精品久久久久久久末码| 国产精品一区www在线观看| 女的被弄到高潮叫床怎么办| 成人美女网站在线观看视频| 人妻一区二区av| 欧美日韩综合久久久久久| 亚洲精品日韩在线中文字幕| 国产成人91sexporn| av又黄又爽大尺度在线免费看| 在线观看国产h片| 国产一区二区亚洲精品在线观看| 色视频在线一区二区三区| 亚洲精品亚洲一区二区| 黄色一级大片看看| 国精品久久久久久国模美| 亚洲综合色惰| 日韩大片免费观看网站| 夫妻午夜视频| 18禁在线无遮挡免费观看视频| 午夜爱爱视频在线播放| 国产中年淑女户外野战色| 在线观看人妻少妇| 免费大片黄手机在线观看| 免费大片黄手机在线观看| 在线a可以看的网站| 22中文网久久字幕| 国产精品av视频在线免费观看| 精品午夜福利在线看| 久久久精品94久久精品| 亚洲精品影视一区二区三区av| 久久久久性生活片| 亚洲在线观看片| 黄色一级大片看看| 国产午夜精品久久久久久一区二区三区| 日韩三级伦理在线观看| 久久久精品94久久精品| 一二三四中文在线观看免费高清| 欧美丝袜亚洲另类| 国产高清三级在线| 熟女av电影| 99精国产麻豆久久婷婷| 久久久午夜欧美精品| 久久久国产一区二区| av国产久精品久网站免费入址| 简卡轻食公司| 久久精品国产a三级三级三级| 日韩成人伦理影院| 成人毛片a级毛片在线播放| 免费电影在线观看免费观看| 久久久久久久精品精品| 国产精品嫩草影院av在线观看| 韩国高清视频一区二区三区| 夜夜爽夜夜爽视频| 夜夜看夜夜爽夜夜摸| 一区二区av电影网| 中国国产av一级| 日本欧美国产在线视频| 九九爱精品视频在线观看| 久久鲁丝午夜福利片| 亚洲第一区二区三区不卡| 一区二区三区精品91| 亚洲自偷自拍三级| 天天躁日日操中文字幕| 99热这里只有精品一区| 人人妻人人爽人人添夜夜欢视频 | 亚洲av日韩在线播放| 久久99精品国语久久久| 在线观看三级黄色| 97热精品久久久久久| 韩国高清视频一区二区三区| tube8黄色片| 成年版毛片免费区| 欧美日韩在线观看h| 久久精品久久久久久噜噜老黄| 一级a做视频免费观看| 日本av手机在线免费观看| 亚洲欧美日韩另类电影网站 | 国产一区亚洲一区在线观看| 嫩草影院新地址| 麻豆成人午夜福利视频| 国产一级毛片在线| 精品视频人人做人人爽| 男女无遮挡免费网站观看| 黄色配什么色好看| 亚洲精品乱码久久久久久按摩| 高清日韩中文字幕在线| 国产成人a∨麻豆精品| 如何舔出高潮| av专区在线播放| 五月开心婷婷网| 亚洲熟女精品中文字幕| 精品一区二区免费观看| 91久久精品国产一区二区三区| 大片免费播放器 马上看| 国产69精品久久久久777片| 亚洲无线观看免费| 亚洲一区二区三区欧美精品 | 亚洲人与动物交配视频| 日韩视频在线欧美| 精品国产露脸久久av麻豆| 欧美亚洲 丝袜 人妻 在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 91久久精品国产一区二区三区| 欧美成人一区二区免费高清观看| 亚洲人成网站在线播| 日韩成人伦理影院| 亚洲精品亚洲一区二区| 一级片'在线观看视频| 色婷婷久久久亚洲欧美| 直男gayav资源| 国产精品.久久久| 久久99精品国语久久久| 欧美少妇被猛烈插入视频| 欧美xxxx性猛交bbbb| 国产精品人妻久久久久久| 午夜日本视频在线| 一本色道久久久久久精品综合| 日本欧美国产在线视频| 国产色爽女视频免费观看| 寂寞人妻少妇视频99o| 国产成人福利小说| 丝袜脚勾引网站| 日韩视频在线欧美| 男人舔奶头视频| 国国产精品蜜臀av免费| 国产在视频线精品| 老女人水多毛片| 天堂中文最新版在线下载 | 少妇猛男粗大的猛烈进出视频 | 亚洲欧美一区二区三区黑人 | 精品人妻偷拍中文字幕| 亚洲av.av天堂| 日韩欧美精品免费久久| 极品教师在线视频| 日本熟妇午夜| 亚洲怡红院男人天堂| 日日摸夜夜添夜夜爱| 精品99又大又爽又粗少妇毛片| www.av在线官网国产| 亚洲真实伦在线观看| 日韩av免费高清视频| 亚洲成人久久爱视频| 国产精品不卡视频一区二区| 亚洲精品亚洲一区二区| 日韩强制内射视频| 色视频www国产| 国产毛片a区久久久久| 亚洲熟女精品中文字幕| 欧美一区二区亚洲| 国产久久久一区二区三区| 九草在线视频观看| 国产成人91sexporn| 亚洲精品国产av成人精品| 一级爰片在线观看| 亚洲av国产av综合av卡| 18禁在线播放成人免费| 男人狂女人下面高潮的视频| 尾随美女入室| 自拍偷自拍亚洲精品老妇| 国产精品三级大全| 国产精品一区www在线观看| 亚洲性久久影院| 日韩精品有码人妻一区| 亚洲,一卡二卡三卡| 九九爱精品视频在线观看| 欧美精品人与动牲交sv欧美| 观看美女的网站| 中文字幕免费在线视频6| 黄色一级大片看看| 亚洲av成人精品一二三区| 国产亚洲av嫩草精品影院| 午夜免费男女啪啪视频观看| 午夜日本视频在线| 亚洲色图综合在线观看| 最近最新中文字幕大全电影3| 欧美三级亚洲精品| 国产精品久久久久久精品电影| 亚洲va在线va天堂va国产| 男女无遮挡免费网站观看| 久久99蜜桃精品久久| 日日摸夜夜添夜夜爱| 老师上课跳d突然被开到最大视频| 三级男女做爰猛烈吃奶摸视频| 久久久久精品久久久久真实原创| 美女cb高潮喷水在线观看| 亚洲三级黄色毛片| 国产 一区精品| 美女cb高潮喷水在线观看| 欧美日本视频| 国产精品秋霞免费鲁丝片| 婷婷色综合www| 亚洲欧美成人精品一区二区| 国产精品国产三级专区第一集| 欧美性感艳星| 国产精品蜜桃在线观看| 亚洲成人av在线免费| 免费av毛片视频| 亚洲欧美日韩另类电影网站 | 一本久久精品| av线在线观看网站| 日本黄大片高清| 免费大片18禁| 亚洲国产精品国产精品| 中文字幕av成人在线电影| 色视频在线一区二区三区| 国产亚洲av片在线观看秒播厂| 欧美区成人在线视频| 国产精品人妻久久久影院| 色婷婷久久久亚洲欧美| 熟女av电影| 26uuu在线亚洲综合色| 亚洲国产色片| 久久精品国产亚洲av涩爱| 婷婷色综合大香蕉| 男人爽女人下面视频在线观看| 久久精品国产亚洲av涩爱| 免费不卡的大黄色大毛片视频在线观看| 久久精品夜色国产| 亚洲色图av天堂| 国产亚洲午夜精品一区二区久久 | av国产免费在线观看| 国产精品一区www在线观看| .国产精品久久| 视频中文字幕在线观看| 免费黄频网站在线观看国产| 国产黄频视频在线观看| 亚洲欧美日韩卡通动漫| 日韩视频在线欧美| 国产午夜精品久久久久久一区二区三区| 亚洲精品影视一区二区三区av| 国产伦精品一区二区三区四那| 国产精品一及| 一级二级三级毛片免费看| av在线app专区| 一区二区三区精品91| 久久午夜福利片| 亚洲丝袜综合中文字幕| 国产一区二区在线观看日韩| 色婷婷久久久亚洲欧美| 大香蕉久久网| 日韩av免费高清视频| 91精品国产九色| 日韩av在线免费看完整版不卡| 亚洲欧美日韩卡通动漫| 男女无遮挡免费网站观看| 亚洲av成人精品一二三区| 国产精品国产三级国产专区5o| 国产精品一区二区在线观看99| 一级毛片电影观看| videos熟女内射| 3wmmmm亚洲av在线观看| 高清欧美精品videossex| 高清视频免费观看一区二区| 五月玫瑰六月丁香| 国产乱来视频区| 精品国产一区二区三区久久久樱花 | 色播亚洲综合网| 日韩,欧美,国产一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 国产黄色免费在线视频| 性插视频无遮挡在线免费观看| 少妇的逼水好多| 少妇人妻 视频| 三级国产精品片| 九九在线视频观看精品| 久久韩国三级中文字幕| 国产成人精品一,二区| 在线观看一区二区三区| 久久精品国产亚洲网站| 两个人的视频大全免费| 观看美女的网站| 国产在线男女| 亚洲一级一片aⅴ在线观看| 亚洲精品456在线播放app| 夜夜看夜夜爽夜夜摸| 国精品久久久久久国模美| 免费人成在线观看视频色| 免费电影在线观看免费观看| 极品教师在线视频| 亚洲精品乱码久久久久久按摩| 男插女下体视频免费在线播放| 国产精品99久久久久久久久| 亚洲最大成人av| 男女那种视频在线观看| 国产黄频视频在线观看| 日韩亚洲欧美综合| 欧美区成人在线视频| 国产成人a∨麻豆精品| 国产精品99久久久久久久久| 久久久久久久午夜电影| 天堂网av新在线| 七月丁香在线播放| 99热这里只有精品一区| av在线观看视频网站免费| 简卡轻食公司| 精品熟女少妇av免费看| 91久久精品国产一区二区成人| 欧美日本视频| 国国产精品蜜臀av免费| 亚洲三级黄色毛片| 一级毛片 在线播放| 高清毛片免费看| 久久精品国产自在天天线| 女的被弄到高潮叫床怎么办| 精品人妻视频免费看| 97精品久久久久久久久久精品| 亚洲电影在线观看av| 王馨瑶露胸无遮挡在线观看| 免费看光身美女| 国产精品爽爽va在线观看网站| 国产淫语在线视频| 亚洲熟女精品中文字幕| 人妻系列 视频| www.av在线官网国产| 欧美精品国产亚洲| 久久精品国产a三级三级三级| 亚洲第一区二区三区不卡| 综合色丁香网| 精品少妇黑人巨大在线播放| 久久99热这里只频精品6学生| 午夜免费男女啪啪视频观看| 国产成人精品婷婷| 国产色婷婷99| xxx大片免费视频| 涩涩av久久男人的天堂| 国产精品久久久久久精品电影| av线在线观看网站| 少妇人妻 视频| 亚洲美女搞黄在线观看| 亚洲av成人精品一二三区| 中文字幕制服av| 美女主播在线视频| 秋霞在线观看毛片| 亚洲精品亚洲一区二区| 51国产日韩欧美| 欧美精品人与动牲交sv欧美| 中文乱码字字幕精品一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 少妇的逼好多水| 国产精品麻豆人妻色哟哟久久| 亚洲av.av天堂| 涩涩av久久男人的天堂| 亚洲av.av天堂| videos熟女内射| 丰满人妻一区二区三区视频av| 国产大屁股一区二区在线视频| 精品国产三级普通话版| 99热这里只有精品一区| 日本与韩国留学比较| 亚洲成人一二三区av| 欧美日韩精品成人综合77777| 国产v大片淫在线免费观看| 国产一区二区亚洲精品在线观看| 国产淫语在线视频| 久久人人爽人人爽人人片va| 联通29元200g的流量卡| 国产黄a三级三级三级人| 成人综合一区亚洲| 色哟哟·www| 欧美高清成人免费视频www| 久久人人爽人人爽人人片va| 联通29元200g的流量卡| 久久人人爽人人爽人人片va| 可以在线观看毛片的网站| 亚洲欧洲日产国产| 国产又色又爽无遮挡免| 激情五月婷婷亚洲| 亚洲国产精品999| 日本av手机在线免费观看| 亚洲在久久综合| 视频中文字幕在线观看| 国产欧美亚洲国产| 色播亚洲综合网| 老司机影院毛片| 日韩免费高清中文字幕av| 一级毛片黄色毛片免费观看视频| 日本一二三区视频观看| 我的老师免费观看完整版| 国产淫片久久久久久久久| 男的添女的下面高潮视频| 少妇熟女欧美另类| 91aial.com中文字幕在线观看| 香蕉精品网在线| 九九久久精品国产亚洲av麻豆| 国产91av在线免费观看| 国产精品久久久久久精品电影| 伊人久久国产一区二区| 国产精品国产三级国产av玫瑰| 国产精品麻豆人妻色哟哟久久| 自拍欧美九色日韩亚洲蝌蚪91 | 七月丁香在线播放| 国产成人a区在线观看| 大片电影免费在线观看免费| 看黄色毛片网站| 偷拍熟女少妇极品色| 在线精品无人区一区二区三 | 一二三四中文在线观看免费高清| 网址你懂的国产日韩在线| 一级毛片我不卡|