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

    基于軌道動力學(xué)特性的小天體物理參數(shù)感知方法

    2021-01-06 07:54:24燕,秦同,喬
    無人系統(tǒng)技術(shù) 2020年5期
    關(guān)鍵詞:測量模型

    羅 燕,秦 同,喬 棟

    (1. 北京理工大學(xué)宇航學(xué)院,北京100081;2. 清華大學(xué)精儀系,北京100084)

    1 引 言

    小天體是太陽系中較古老的一類天體,探測小天體,獲取其物理參數(shù),是研究太陽系形成與演化、生命起源與進化,以及抵御外來天體撞擊的重要技術(shù)途徑,是人類探索和開發(fā)宇宙資源的重要技術(shù)手段。同時,也為試驗驗證空間技術(shù),尤其是深空探測技術(shù)提供了平臺。目前,小天體物理參數(shù)的獲取方式主要包括近距離光學(xué)成像、多譜段光譜測量和表面采樣探測。此類方式可以感知小天體的尺寸、形狀、自旋等物理參數(shù),以及小天體表面地形、地表物質(zhì)組成等。但光學(xué)觀測的方式精度較低,且無法獲得小天體內(nèi)部質(zhì)量與密度分布等信息。小天體的引力場以及自旋參數(shù)是小天體探測的基本內(nèi)容之一。感知高精度的小天體引力場以及自旋,特別是不規(guī)則形狀小天體引力場與自旋,有助于探索小天體的內(nèi)部結(jié)構(gòu)、質(zhì)量與密度分布等信息,揭示小天體起源,獲得新的科學(xué)發(fā)現(xiàn)。

    小天體引力場以及自旋信息的獲取通常有兩種方式。其一,從傳統(tǒng)的光學(xué)觀測中提取光度或光學(xué)特征點等有效信息,以此為觀測量反演自旋角速度以及外部形貌[1-2],在均勻密度的假設(shè)下,獲得質(zhì)量分布,從而建立引力場模型[3-5]。由于這類方法的信息源來自小天體表面光學(xué)觀測,因此獲得的自旋角速度信息以及引力場模型較為粗糙,且無法體現(xiàn)小天體內(nèi)部結(jié)構(gòu)與質(zhì)量信息。另一種方法是通過探測器軌道數(shù)據(jù)進行反演。未來的探測器需具有全面的感知能力,包括環(huán)境感知及自身狀態(tài)感知[6-7]。由于探測器的運動由其軌道動力學(xué)環(huán)境,即探測器所處引力場決定,且小天體自旋會影響引力場的空間分布,因此借助精確的軌道數(shù)據(jù),可對小天體引力場與自旋進行精確感知[8]。NEAR 探測器在對Eros 433 小行星的近距離飛越過程中,通過地面的無線電測量與星載光學(xué)測量數(shù)據(jù)不僅實現(xiàn)了高精度定軌[9],同時還對Eros 433 的引力場進行了估計,并在美國國家航空航天局官網(wǎng)公布了引力場球諧系數(shù)模型的前8 階球諧系數(shù),是通過軌道數(shù)據(jù)感知小天體引力場的先例。

    高精度測量與動力學(xué)估計是感知小天體自旋角速度與引力場模型的前提?,F(xiàn)階段,對小天體探測器的測量主要為地面無線電測量,包括測距、測速、甚長基線干涉測量,以及星載光學(xué)相機測視線與激光測距儀測距[10-12]。綜合各類測量信息,結(jié)合較粗糙的先驗引力場模型,可以估計探測器在小天體坐標(biāo)系的軌道。然而,地面測控定軌自主性差、精度隨距離增加而降低、且測量幾何構(gòu)型變化弱,無法實現(xiàn)小天體坐標(biāo)系下的精密定軌。此外,小天體自旋與引力場感知需要長時間連續(xù)測量估計,依靠地面測控成本高昂。光學(xué)測量雖自主性強,但受光學(xué)導(dǎo)航原理(圖像特征提取與匹配)的限制,參數(shù)估計精度較低。因此,傳統(tǒng)的測量定軌方式在定軌精度上只能滿足工程需求,無法滿足高精度引力場與自旋參數(shù)感知的需求。

    本文針對無人系統(tǒng)在小天體科學(xué)探測方面的應(yīng)用,研究基于軌道動力學(xué)特性的小天體引力場與自旋角速度感知問題,即如何通過小天體探測器之間的器間測距信息對小天體引力場模型參數(shù)與自旋角速度進行修正。該問題本質(zhì)上為參數(shù)估計問題,可用的測量信息為器間測距,需估計的參數(shù)為引力場模型參數(shù)、自旋角速度以及探測器軌道。在二體動力學(xué)問題中,由于引力場的對稱性,器間相對測量僅能確定探測器軌道的大小、形狀與尺寸,但星座的旋轉(zhuǎn)方位信息,即各個軌道升交點赤經(jīng)、軌道傾角以及近地點幅角均不可估計[13-14]。原因在于,星座整體旋轉(zhuǎn)時,器間的測距、測速信息不會發(fā)生變化,導(dǎo)致器間測量導(dǎo)航系統(tǒng)不完全可觀。但在小天體的不規(guī)則引力場中,引力場各向差異顯著,不再呈現(xiàn)對稱性,因此星座旋轉(zhuǎn)后,各軌道的大小、形狀尺寸也會隨之改變,從而引起器間測量的變化。如此一來,器間測量對軌道各個元素均有一定的敏感性,整個星座定軌系統(tǒng)完全可觀。因此,小天體引力場的不規(guī)則性雖然增加了小天體附近軌道設(shè)計難度,但也由于其復(fù)雜性,改善了參數(shù)感知系統(tǒng)的可觀性,是參數(shù)感知的有利因素。

    小天體的不規(guī)則引力場是實現(xiàn)參數(shù)感知的關(guān)鍵因素,同時,由于先驗引力場模型以及自旋角速度存在誤差,在進行定軌估計時,除了各軌道狀態(tài),還需要將引力場模型系數(shù)以及自旋角速度作為待估狀態(tài)變量進行估計。因此,基于動力學(xué)的小天體參數(shù)感知實質(zhì)上為動力學(xué)參數(shù)估計。以小天體球諧系數(shù)模型為例,通過三維建模與均勻密度假設(shè)得到的引力場模型中各階球諧系數(shù)誤差較大,將各階球諧系數(shù)以及自旋角速度作為狀態(tài)變量,擴充至動力學(xué)估計系統(tǒng)的狀態(tài)向量中,依靠器間觀測量,通過非線性濾波算法同時估計各階球諧系數(shù)、自旋角速度以及探測器軌道參數(shù)。

    2 小天體附近探測器軌道動力學(xué)

    2.1 小天體引力場模型

    小行星形狀不規(guī)則,內(nèi)部密度的分布也不均勻,這使得中心天體對環(huán)繞飛行器的引力作用也不嚴(yán)格等于等效的質(zhì)心引力,而是等于中心天體各部分質(zhì)量引力作用之和,即引力位函數(shù):

    其中,r為飛行器的位置,s為小天體內(nèi)某一點的位置,上述積分應(yīng)該對整個中心天體體積進行積分,中心天體對環(huán)繞飛行器的引力可表示為:

    由于我們無法得到精確的中心天體的質(zhì)量分布結(jié)果,所以上式并不能用來計算中心天體的引力勢。常用的中心天體引力的計算模型有球諧系數(shù)模型[15]、多面體模型[16]、質(zhì)點群模型等[17]。本文用球諧函數(shù)的級數(shù)展開式表示小行星的引力勢能,能夠更直觀地描述小天體模型精度階次。球諧系數(shù)模型如式(3)所示,小行星引力向量可以通過對勢能求梯度得到。

    其中:

    其中,φ和λ分別為飛行器在小天體固連坐標(biāo)系下的緯度和經(jīng)度,φ′和λ′分別為積分單元的緯度和經(jīng)度,Cnm和Snm為非歸一化的引力場模型諧函數(shù)和為歸一化的引力場模型諧函數(shù)。Pnm為勒讓德多項式,δnm為克羅內(nèi)克符號(即n,m相同為1,不同為0),GM為引力常數(shù),R為天體平均半徑。

    由于Cnm和Snm的計算同樣需要知道質(zhì)量分布,故一般也無法計算。目前而言,諧函數(shù)和一般是發(fā)射衛(wèi)星并通過定軌來獲得,從而通過換算得到Cnm和Snm。

    m=n= 0的項為中心天體的質(zhì)心引力部分,當(dāng)n,m不為0 時,代表中心天體引力位偏離中心質(zhì)點的程度。

    2.2 小天體附近探測器軌道動力學(xué)模型

    在小天體固連坐標(biāo)系下,探測器軌道動力學(xué)模型如式(9)所示:

    式中,xi,yi,zi為探測器i的三軸位置,ω為小天體的自旋角速度,nx,ny為未建模的干擾加速度,V為小天體的引力位勢函數(shù)球諧系數(shù)模型。

    3 引力場與自旋參數(shù)感知系統(tǒng)

    3.1 系統(tǒng)狀態(tài)模型

    在小天體固連坐標(biāo)系下建立引力場與自旋感知系統(tǒng),待估系統(tǒng)變量為兩環(huán)繞探測器的位置、速度、引力場球諧系數(shù)、模型系數(shù),以及自旋角速度,系統(tǒng)狀態(tài)變量如式(10)所示:

    小天體各級引力場參數(shù)以及自旋角速度視為常數(shù),其狀態(tài)模型如式(11)所示:

    由式(10)~(11)可得到引力場反演的系統(tǒng)狀態(tài)模型如式(12)所示:

    f(X)的具體表達式參照式(10)~(11)。

    在兩探測器繞飛小天體的過程中,通過器間特高頻頻段無線電通信,實現(xiàn)相對距離測量,測量模型如式(13)所示:

    式中,wρ為測量噪聲,在此假設(shè)為高斯白噪聲。

    3.2 系統(tǒng)狀態(tài)估計方法

    針對小天體引力場反演這一非線性估計問題,采用無跡卡爾曼濾波算法實現(xiàn)狀態(tài)最優(yōu)估計,給出tk時刻的系統(tǒng)先驗狀態(tài)Xk和狀態(tài)誤差協(xié)方差矩陣Pk,后續(xù)的估計流程如下:

    (1)計算2n+ 1個sigma點及其權(quán)值。

    式中,λ=α2(n+k)-n,α決定sigma 點的散布程度,通常取一個小的正值,k通常取0,β用來描述系統(tǒng)狀態(tài)的分析信息,高斯分布情況下,β最優(yōu)取值為2。

    (2)計算sigma 點通過非線性函數(shù)f(X)的傳播結(jié)果。

    進而可得:

    (3)觀測更新。

    4 仿真結(jié)果及分析

    對目標(biāo)小行星2016HO3 進行引力場與自旋參數(shù)感知仿真分析,假設(shè)目標(biāo)天體等效半徑50m,密度1.26×103kg/m3(標(biāo)準(zhǔn)C 類小行星密度,2016HO3富含鐵質(zhì),實際密度要大),套用Bennu 的不規(guī)則引力場模型,前四階引力場參考模型參數(shù)如表1 所示。假設(shè)先驗引力場各階模型誤差為50%,兩探測器在本體系下的初始位置誤差標(biāo)準(zhǔn)差為10m,器間測距精度為10cm/1cm,探測器繞飛高度1km/500m/100m。

    表1 前八階引力場參考模型參數(shù)Table 1 The first 8th order gravity model parameters

    當(dāng)測量精度為10cm 時,不同軌道高度的引力場反演仿真結(jié)果如圖1 所示。軌道高度1km 時,通過240h(包含不可見弧段時長)的測量,可以反演出二階引力場,軌道500m時,通過120h觀測可以反演出三階引力場,軌道高度100m時,通過240h觀測可以反演出四階引力場。不同軌道高度下探測器位置收斂精度近似,三軸位置精度約為20cm,如圖2所示。

    圖1 測量精度10cm時不同軌道高度下感知的引力場模型參數(shù)誤差Fig.1 The gravity parameter errors with 10cm measurement accuracy in situations with different orbits

    圖2 測量精度10cm時探測器軌道確定精度Fig.2 The orbit determination errors with 10cm measurement accuracy

    考慮到2016HO3 的半徑約為50m,因此100m軌道高度已接近工程可實現(xiàn)的極限最低高度。為反演更高階次的引力場模型,需進一步提高測量精度。當(dāng)測量精度達到1cm 時,100m 高度的探測器器間測距反演引力場結(jié)果如圖3 所示,對應(yīng)的探測器軌道精度如圖4 所示。從仿真結(jié)果可以看出,當(dāng)測量精度提升后,引力場反演精度明顯提高,在100m 軌道高度進行器間測量反演,可修正六階引力場參數(shù)。同時,探測器的定軌精度也從20cm 提高至5cm。通過引力場參數(shù)計算得到的真實引力加速度與感知加速度仿真對比如圖5所示。圖6給出了自旋角速度的誤差。假設(shè)光學(xué)觀測獲取的先驗自旋角速度誤差為1%,通過動力學(xué)感知方法,經(jīng)過170h 連續(xù)觀測后,自旋角速度誤差降低至0.01%。

    圖3 測量精度1cm時100m軌道高度120h觀測感知的引力場模型參數(shù)誤差Fig.3 The gravity parameter errors with 1cm measurement accuracy and 100m-height orbit within 120h

    5 結(jié) 論

    圖4 測量精度1cm時探測器軌道確定精度Fig. 4 The orbit determination errors with 1cm measurement accuracy

    圖5 小天體表面真實引力加速度與感知加速度對比Fig. 5 The comparion of real and perceived acceleration errors

    圖6 自旋角速度感知誤差Fig. 6 The self-rotation speed error

    本文研究了基于軌道動力學(xué)特性的小天體引力場與自旋角速度感知方法,利用兩個小天體環(huán)繞探測器之間的無線電測距信息,可同時估計探測器在小天體固連坐標(biāo)系下的軌道、各階球諧系數(shù)模型參數(shù)以及小天體自旋角速度。影響引力場與自旋感知精度的因素主要包括探測器軌道高度和測量精度兩方面。仿真結(jié)果驗證了軌道高度越低,測量精度越高,可以感知的引力場階次越高,自旋角速度精度越高。精確的引力場信息是實施小天體探測任務(wù)的重要保障,本文的研究對空間環(huán)境感知、小天體探測中的參數(shù)估計有重要參考價值,創(chuàng)新點在于以器間測距信息為輸入對小天體的引力場進行了較為精確的估計,同時提高了探測器的定軌精度。但利用球諧系數(shù)計算小行星非球星引力存在著一定的缺點,小行星形狀越不規(guī)則時誤差越大,本文僅套用Bennu 的引力場模型進行了模擬仿真,后期可利用仿真或真實拍攝的圖片初步估計小行星的形狀,并采用更高階的非線性估計方法,進一步驗證該方法和提高對小行星引力場參數(shù)估計的精度。

    猜你喜歡
    測量模型
    一半模型
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    把握四個“三” 測量變簡單
    滑動摩擦力的測量和計算
    滑動摩擦力的測量與計算
    測量的樂趣
    3D打印中的模型分割與打包
    測量
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    日日撸夜夜添| 国产精品嫩草影院av在线观看| 99热网站在线观看| 汤姆久久久久久久影院中文字幕| 女人久久www免费人成看片| 精品少妇黑人巨大在线播放| 九九爱精品视频在线观看| 麻豆精品久久久久久蜜桃| 在现免费观看毛片| 性色avwww在线观看| 91在线精品国自产拍蜜月| 黑人猛操日本美女一级片| 亚洲国产av新网站| 亚洲av国产av综合av卡| 人妻夜夜爽99麻豆av| 国产成人一区二区在线| 午夜福利视频在线观看免费| 亚洲精品国产av蜜桃| 亚洲欧美日韩另类电影网站| 日韩欧美精品免费久久| 天堂8中文在线网| 又粗又硬又长又爽又黄的视频| 草草在线视频免费看| 大香蕉久久成人网| 人妻夜夜爽99麻豆av| 九九在线视频观看精品| 欧美激情国产日韩精品一区| av专区在线播放| 免费少妇av软件| 丝瓜视频免费看黄片| 久热久热在线精品观看| 在线 av 中文字幕| 黄色怎么调成土黄色| 免费久久久久久久精品成人欧美视频 | 欧美日本中文国产一区发布| 免费人成在线观看视频色| 黄色欧美视频在线观看| 免费黄网站久久成人精品| 欧美3d第一页| 亚洲精品成人av观看孕妇| 99久国产av精品国产电影| 精品卡一卡二卡四卡免费| 国产精品一区二区在线不卡| 爱豆传媒免费全集在线观看| 亚洲美女视频黄频| 91aial.com中文字幕在线观看| 91午夜精品亚洲一区二区三区| 满18在线观看网站| 天堂中文最新版在线下载| 国产av码专区亚洲av| videosex国产| 欧美激情国产日韩精品一区| 中国美白少妇内射xxxbb| 成人综合一区亚洲| 久久久a久久爽久久v久久| 涩涩av久久男人的天堂| 99久久精品国产国产毛片| 久久久久网色| 亚洲精品乱码久久久v下载方式| 啦啦啦视频在线资源免费观看| 老女人水多毛片| 国产精品久久久久久精品电影小说| av卡一久久| 亚洲av电影在线观看一区二区三区| 97在线视频观看| 成人国语在线视频| 日本vs欧美在线观看视频| 看非洲黑人一级黄片| 视频在线观看一区二区三区| 成年美女黄网站色视频大全免费 | 日韩不卡一区二区三区视频在线| 午夜日本视频在线| 国产免费一级a男人的天堂| 国产精品久久久久久久电影| 久久久久久久久久久免费av| 色5月婷婷丁香| 在线播放无遮挡| 99久久精品一区二区三区| 亚洲少妇的诱惑av| 18禁观看日本| 下体分泌物呈黄色| 久久久久国产精品人妻一区二区| 下体分泌物呈黄色| 麻豆乱淫一区二区| 性色av一级| 99热这里只有精品一区| 欧美激情 高清一区二区三区| 欧美精品高潮呻吟av久久| 黄色欧美视频在线观看| 欧美一级a爱片免费观看看| 久久精品熟女亚洲av麻豆精品| 日本欧美视频一区| 久久精品熟女亚洲av麻豆精品| 一区二区日韩欧美中文字幕 | 男女国产视频网站| 日韩强制内射视频| 最近中文字幕高清免费大全6| 亚洲精品av麻豆狂野| 男男h啪啪无遮挡| 国产精品嫩草影院av在线观看| 丝袜脚勾引网站| 建设人人有责人人尽责人人享有的| 国产视频内射| 如何舔出高潮| 中文字幕免费在线视频6| 国产精品久久久久久久久免| 晚上一个人看的免费电影| 中文字幕免费在线视频6| 国产深夜福利视频在线观看| 久久精品国产a三级三级三级| 亚洲国产精品国产精品| 丝袜喷水一区| 久久久国产欧美日韩av| 在线观看免费视频网站a站| 精品久久久噜噜| 午夜免费男女啪啪视频观看| 国产不卡av网站在线观看| 日韩,欧美,国产一区二区三区| 亚洲成色77777| 99视频精品全部免费 在线| 久久99一区二区三区| 91aial.com中文字幕在线观看| 制服诱惑二区| 美女国产视频在线观看| xxx大片免费视频| 国产亚洲最大av| 亚洲精品久久午夜乱码| 亚洲精品国产av成人精品| 51国产日韩欧美| 妹子高潮喷水视频| 久久久久久人妻| 精品一品国产午夜福利视频| 热99久久久久精品小说推荐| 不卡视频在线观看欧美| 久久国产精品大桥未久av| 免费观看在线日韩| 亚洲精品乱码久久久久久按摩| 国产熟女午夜一区二区三区 | 三级国产精品片| 日韩一区二区三区影片| 国产欧美日韩一区二区三区在线 | 成人毛片60女人毛片免费| 日韩av在线免费看完整版不卡| 日产精品乱码卡一卡2卡三| 国产一区二区在线观看日韩| 成人国产av品久久久| 国产精品久久久久久精品电影小说| 人妻人人澡人人爽人人| 亚洲精品日本国产第一区| 亚洲五月色婷婷综合| 高清毛片免费看| 欧美97在线视频| 在线精品无人区一区二区三| 精品国产一区二区三区久久久樱花| 欧美精品一区二区大全| 黄片无遮挡物在线观看| 亚洲色图 男人天堂 中文字幕 | 午夜视频国产福利| 国产在线视频一区二区| 老司机影院毛片| 人成视频在线观看免费观看| 哪个播放器可以免费观看大片| 99热国产这里只有精品6| www.色视频.com| 男女边吃奶边做爰视频| 国产av国产精品国产| 久久午夜福利片| 97精品久久久久久久久久精品| 亚洲精品456在线播放app| 国产亚洲av片在线观看秒播厂| h视频一区二区三区| 99国产精品免费福利视频| 日本黄大片高清| 国产黄频视频在线观看| 久久久久久伊人网av| 97在线视频观看| 美女福利国产在线| 曰老女人黄片| 国产精品一二三区在线看| 欧美成人午夜免费资源| 成人午夜精彩视频在线观看| 99re6热这里在线精品视频| 欧美激情极品国产一区二区三区 | 五月开心婷婷网| 高清毛片免费看| 考比视频在线观看| 成人亚洲精品一区在线观看| 久久久久久久久久久久大奶| av福利片在线| 精品久久久久久久久av| av卡一久久| 熟女人妻精品中文字幕| 日韩电影二区| 欧美97在线视频| 国产伦理片在线播放av一区| 一级毛片黄色毛片免费观看视频| 国产午夜精品久久久久久一区二区三区| 国产一区亚洲一区在线观看| 亚洲,一卡二卡三卡| 亚洲精品国产av成人精品| 亚洲av不卡在线观看| 亚洲欧美色中文字幕在线| 不卡视频在线观看欧美| 亚洲av欧美aⅴ国产| www.色视频.com| 9色porny在线观看| 日本欧美国产在线视频| 内地一区二区视频在线| 久久毛片免费看一区二区三区| 九九在线视频观看精品| 亚洲精品日韩av片在线观看| 久久人人爽人人爽人人片va| 亚洲欧美日韩卡通动漫| 国产精品99久久99久久久不卡 | 欧美日韩亚洲高清精品| 亚洲不卡免费看| 亚洲无线观看免费| 两个人的视频大全免费| 日韩精品有码人妻一区| 在线观看美女被高潮喷水网站| 亚洲国产色片| 国产深夜福利视频在线观看| 麻豆成人av视频| 天天躁夜夜躁狠狠久久av| 人人妻人人澡人人看| 美女xxoo啪啪120秒动态图| 简卡轻食公司| 久久久a久久爽久久v久久| 十八禁网站网址无遮挡| 国产免费一级a男人的天堂| 又大又黄又爽视频免费| 国产高清国产精品国产三级| 日本欧美国产在线视频| 老熟女久久久| 精品人妻熟女毛片av久久网站| 男的添女的下面高潮视频| 日韩制服骚丝袜av| 免费大片18禁| 人妻少妇偷人精品九色| 十八禁高潮呻吟视频| 久久亚洲国产成人精品v| 亚洲精品日韩av片在线观看| 一二三四中文在线观看免费高清| 亚洲美女搞黄在线观看| 狠狠婷婷综合久久久久久88av| 亚洲欧洲日产国产| 超碰97精品在线观看| 久久国产精品大桥未久av| 精品午夜福利在线看| 麻豆精品久久久久久蜜桃| 欧美精品高潮呻吟av久久| 伦理电影大哥的女人| 日韩av在线免费看完整版不卡| 各种免费的搞黄视频| 能在线免费看毛片的网站| 99热6这里只有精品| 国产色婷婷99| 国产欧美亚洲国产| 2022亚洲国产成人精品| 十分钟在线观看高清视频www| 国产日韩一区二区三区精品不卡 | 中文字幕免费在线视频6| 男女无遮挡免费网站观看| 插逼视频在线观看| 亚洲美女黄色视频免费看| 国产国拍精品亚洲av在线观看| 国产av国产精品国产| 超碰97精品在线观看| 国产午夜精品一二区理论片| 九九在线视频观看精品| 国产 一区精品| 久久久精品94久久精品| 少妇猛男粗大的猛烈进出视频| 99热6这里只有精品| 国产一区有黄有色的免费视频| 岛国毛片在线播放| 国产精品国产三级专区第一集| 国产欧美另类精品又又久久亚洲欧美| 18禁在线播放成人免费| 在线 av 中文字幕| 国产在线一区二区三区精| 伊人久久精品亚洲午夜| 国产午夜精品久久久久久一区二区三区| 男女边吃奶边做爰视频| av国产久精品久网站免费入址| 一边摸一边做爽爽视频免费| 久久久a久久爽久久v久久| av有码第一页| 久久精品国产自在天天线| 欧美性感艳星| 熟女人妻精品中文字幕| 男女无遮挡免费网站观看| 制服丝袜香蕉在线| 日韩中文字幕视频在线看片| 色婷婷久久久亚洲欧美| 日韩视频在线欧美| av网站免费在线观看视频| 中文精品一卡2卡3卡4更新| 久久久久网色| 大又大粗又爽又黄少妇毛片口| 肉色欧美久久久久久久蜜桃| 一本色道久久久久久精品综合| 纵有疾风起免费观看全集完整版| 国产成人91sexporn| 亚洲,一卡二卡三卡| 丝袜在线中文字幕| 97超碰精品成人国产| 中文天堂在线官网| 自拍欧美九色日韩亚洲蝌蚪91| 成年人午夜在线观看视频| 国产免费又黄又爽又色| 欧美三级亚洲精品| 啦啦啦在线观看免费高清www| 99热6这里只有精品| 日韩电影二区| 99国产综合亚洲精品| 日本wwww免费看| 嫩草影院入口| 亚洲国产毛片av蜜桃av| 九九在线视频观看精品| 久久久久久伊人网av| 国产在视频线精品| 在线观看一区二区三区激情| 一级毛片黄色毛片免费观看视频| 五月玫瑰六月丁香| 美女国产高潮福利片在线看| 一个人免费看片子| 精品人妻偷拍中文字幕| 免费久久久久久久精品成人欧美视频 | 久久精品久久久久久久性| 人妻 亚洲 视频| 成人毛片60女人毛片免费| 国产精品免费大片| 成人无遮挡网站| 午夜福利视频在线观看免费| 亚洲国产精品一区二区三区在线| 国产精品三级大全| 少妇人妻精品综合一区二区| 日韩一区二区视频免费看| 男的添女的下面高潮视频| 有码 亚洲区| 性色av一级| 国产一区二区三区综合在线观看 | 免费观看a级毛片全部| 91久久精品国产一区二区成人| 国产爽快片一区二区三区| 女人久久www免费人成看片| 草草在线视频免费看| 久久精品久久精品一区二区三区| 国产爽快片一区二区三区| 久久精品人人爽人人爽视色| 一级毛片 在线播放| 99久久精品国产国产毛片| 中文字幕久久专区| 国产精品成人在线| 卡戴珊不雅视频在线播放| 丁香六月天网| 午夜福利网站1000一区二区三区| 午夜久久久在线观看| 999精品在线视频| 国产男女超爽视频在线观看| 亚洲av.av天堂| 下体分泌物呈黄色| kizo精华| 国产精品无大码| 91精品伊人久久大香线蕉| 在线观看一区二区三区激情| 亚洲av日韩在线播放| 男女边吃奶边做爰视频| a级毛色黄片| 伊人久久精品亚洲午夜| 日本wwww免费看| 成人无遮挡网站| 天天影视国产精品| 日本免费在线观看一区| 男女高潮啪啪啪动态图| 久久99热6这里只有精品| 亚洲中文av在线| 啦啦啦在线观看免费高清www| 国产精品99久久99久久久不卡 | 国产av精品麻豆| 亚州av有码| 简卡轻食公司| 色5月婷婷丁香| 国产成人免费观看mmmm| 免费观看的影片在线观看| 精品一区二区三卡| 欧美精品一区二区免费开放| 91久久精品国产一区二区成人| 丁香六月天网| 亚洲精品日韩在线中文字幕| 午夜免费鲁丝| 亚洲少妇的诱惑av| 人成视频在线观看免费观看| 黑人欧美特级aaaaaa片| 欧美成人午夜免费资源| 寂寞人妻少妇视频99o| 成人18禁高潮啪啪吃奶动态图 | 一个人看视频在线观看www免费| 哪个播放器可以免费观看大片| 丰满少妇做爰视频| 成年女人在线观看亚洲视频| av播播在线观看一区| 午夜91福利影院| 国产av一区二区精品久久| 街头女战士在线观看网站| videossex国产| 精品视频人人做人人爽| 免费人妻精品一区二区三区视频| 青春草亚洲视频在线观看| 另类精品久久| 99国产综合亚洲精品| 精品久久久久久久久av| 啦啦啦啦在线视频资源| 久久精品人人爽人人爽视色| 日本猛色少妇xxxxx猛交久久| 999精品在线视频| 你懂的网址亚洲精品在线观看| 免费久久久久久久精品成人欧美视频 | 欧美三级亚洲精品| 国产精品三级大全| 亚洲精品中文字幕在线视频| 国产免费现黄频在线看| 精品久久蜜臀av无| 久久人人爽人人片av| 色94色欧美一区二区| 男女国产视频网站| 成年av动漫网址| 天堂中文最新版在线下载| 免费久久久久久久精品成人欧美视频 | 全区人妻精品视频| 久久久久久伊人网av| 亚洲人与动物交配视频| 日本wwww免费看| 欧美日韩av久久| 看免费成人av毛片| 国产国拍精品亚洲av在线观看| 国产视频内射| 国产高清三级在线| 人人妻人人添人人爽欧美一区卜| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产精品一区二区三区在线| 日韩成人av中文字幕在线观看| 欧美 日韩 精品 国产| 亚洲精品aⅴ在线观看| 国产成人免费无遮挡视频| 赤兔流量卡办理| 亚洲伊人久久精品综合| 亚洲精品成人av观看孕妇| 色网站视频免费| 日本av手机在线免费观看| 国产av一区二区精品久久| 中国三级夫妇交换| 飞空精品影院首页| 亚洲激情五月婷婷啪啪| 最后的刺客免费高清国语| a级片在线免费高清观看视频| 99久久中文字幕三级久久日本| 水蜜桃什么品种好| 亚洲av二区三区四区| 成人漫画全彩无遮挡| 亚洲国产毛片av蜜桃av| 看非洲黑人一级黄片| 少妇人妻精品综合一区二区| 色网站视频免费| 国产av国产精品国产| av专区在线播放| 日韩制服骚丝袜av| 美女视频免费永久观看网站| 人成视频在线观看免费观看| 亚洲国产最新在线播放| 国产在线视频一区二区| 国产片内射在线| 在线观看一区二区三区激情| 久久久久精品久久久久真实原创| 秋霞在线观看毛片| 在线看a的网站| 好男人视频免费观看在线| 国产伦理片在线播放av一区| 亚洲av福利一区| 久久久久久久久久久丰满| 如日韩欧美国产精品一区二区三区 | 亚洲高清免费不卡视频| 精品人妻在线不人妻| 我的女老师完整版在线观看| 晚上一个人看的免费电影| 日本与韩国留学比较| 下体分泌物呈黄色| 一级毛片黄色毛片免费观看视频| 国产黄片视频在线免费观看| 国产av一区二区精品久久| 亚州av有码| 亚洲少妇的诱惑av| av福利片在线| 久久久精品区二区三区| 久久久久久久久久久久大奶| 春色校园在线视频观看| 国产精品一区二区在线观看99| 亚洲情色 制服丝袜| 人妻 亚洲 视频| 亚洲精品乱久久久久久| 久久久久久人妻| 久久久午夜欧美精品| 成年人午夜在线观看视频| 国产免费现黄频在线看| 视频中文字幕在线观看| 两个人免费观看高清视频| 一本色道久久久久久精品综合| 久久久久久久国产电影| 只有这里有精品99| 水蜜桃什么品种好| 免费高清在线观看日韩| 综合色丁香网| 久久久精品免费免费高清| 国产一区二区三区av在线| 王馨瑶露胸无遮挡在线观看| 免费看不卡的av| 亚洲熟女精品中文字幕| 夜夜骑夜夜射夜夜干| 伦理电影大哥的女人| 中国国产av一级| 亚洲精品一二三| 男女国产视频网站| 夜夜爽夜夜爽视频| 一本—道久久a久久精品蜜桃钙片| 久久久久久久大尺度免费视频| 欧美日韩成人在线一区二区| 亚洲国产色片| 99视频精品全部免费 在线| 久久精品久久精品一区二区三区| 午夜视频国产福利| xxx大片免费视频| 十八禁高潮呻吟视频| 日本色播在线视频| 精品人妻熟女av久视频| 精品国产露脸久久av麻豆| 啦啦啦视频在线资源免费观看| 中国美白少妇内射xxxbb| 亚洲国产欧美在线一区| 乱人伦中国视频| 久久99蜜桃精品久久| 99热6这里只有精品| 大话2 男鬼变身卡| 91精品国产国语对白视频| 久久亚洲国产成人精品v| 国产精品免费大片| 黄色欧美视频在线观看| 久久久久久伊人网av| 久久精品国产鲁丝片午夜精品| 国产精品一国产av| 国产精品99久久99久久久不卡 | 国产成人一区二区在线| 男女边吃奶边做爰视频| av在线老鸭窝| 亚洲av.av天堂| 嘟嘟电影网在线观看| 一级,二级,三级黄色视频| 成年人午夜在线观看视频| 日本免费在线观看一区| 一级毛片我不卡| 黑丝袜美女国产一区| 91精品国产国语对白视频| 国产成人freesex在线| 免费观看性生交大片5| 91精品一卡2卡3卡4卡| 狂野欧美激情性xxxx在线观看| 亚洲丝袜综合中文字幕| 97超视频在线观看视频| 日日啪夜夜爽| 久久久久精品久久久久真实原创| 999精品在线视频| 亚洲美女黄色视频免费看| 一边摸一边做爽爽视频免费| 亚洲一区二区三区欧美精品| 欧美精品亚洲一区二区| 久久99一区二区三区| 大香蕉97超碰在线| 最近手机中文字幕大全| 最近中文字幕2019免费版| 一区二区三区乱码不卡18| 精品国产露脸久久av麻豆| 丝袜美足系列| 久久久久久久久久人人人人人人| av国产久精品久网站免费入址| 久久国产精品大桥未久av| videossex国产| 亚洲久久久国产精品| 最近中文字幕高清免费大全6| 久久精品国产亚洲av涩爱| 免费观看无遮挡的男女| 看十八女毛片水多多多| 美女cb高潮喷水在线观看| 寂寞人妻少妇视频99o| 老司机亚洲免费影院| 天堂中文最新版在线下载| 国产精品熟女久久久久浪| 人成视频在线观看免费观看| 国产亚洲精品久久久com| 成人亚洲欧美一区二区av| 最后的刺客免费高清国语| av.在线天堂| 亚洲av免费高清在线观看| 天天躁夜夜躁狠狠久久av| 人妻少妇偷人精品九色| 亚洲精品久久久久久婷婷小说| 妹子高潮喷水视频| 搡女人真爽免费视频火全软件| 中文字幕人妻熟人妻熟丝袜美| 国产欧美另类精品又又久久亚洲欧美| 一级爰片在线观看| 多毛熟女@视频| 91精品一卡2卡3卡4卡| 午夜免费鲁丝| 最近中文字幕2019免费版| 插阴视频在线观看视频|