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

    UKF容錯(cuò)濾波在脈沖星組合導(dǎo)航中的應(yīng)用

    2014-10-21 01:08:48北京理工大學(xué)機(jī)電學(xué)院北京100081中國(guó)科學(xué)院空間科學(xué)與應(yīng)用研究中心北京100190
    關(guān)鍵詞:野值脈沖星質(zhì)心

    (1. 北京理工大學(xué) 機(jī)電學(xué)院,北京 100081;2. 中國(guó)科學(xué)院 空間科學(xué)與應(yīng)用研究中心,北京 100190)

    (1. 北京理工大學(xué) 機(jī)電學(xué)院,北京 100081;2. 中國(guó)科學(xué)院 空間科學(xué)與應(yīng)用研究中心,北京 100190)

    針對(duì)X射線脈沖星組合導(dǎo)航系統(tǒng)中測(cè)量野值導(dǎo)致的導(dǎo)航精度下降問(wèn)題,利用基于殘差正交性的Unscented Kalman Filter (UKF)容錯(cuò)濾波方法進(jìn)行野值修正。在深空巡航段,采用基于X射線脈沖星和太陽(yáng)觀測(cè)的組合導(dǎo)航方式,用X射線探測(cè)器測(cè)量脈沖星光子到達(dá)時(shí)間,利用太陽(yáng)敏感器測(cè)量太陽(yáng)視線方向矢量,并利用聯(lián)邦濾波結(jié)構(gòu)進(jìn)行信息融合。仿真結(jié)果表明,基于 UKF的容錯(cuò)濾波算法在脈沖星組合導(dǎo)航系統(tǒng)的應(yīng)用中,能夠?qū)σ爸颠M(jìn)行實(shí)時(shí)修正,避免了導(dǎo)航精度下降,提高了系統(tǒng)的魯棒性和工程實(shí)用價(jià)值。

    脈沖星;太陽(yáng)敏感器;UKF濾波;野值檢測(cè);容錯(cuò)濾波

    在深空探測(cè)活動(dòng)中,往往存在時(shí)間延遲和精度不高的問(wèn)題,而且地面測(cè)控資源有限,因此深空探測(cè)需要自主導(dǎo)航技術(shù)的支撐。X射線脈沖星自主導(dǎo)航技術(shù)具有高精度、高自主性的特點(diǎn),是一種很有應(yīng)用前景的航天器自主導(dǎo)航技術(shù)。該技術(shù)也將在美國(guó)宇航局NICER任務(wù)[1]中得到驗(yàn)證,我國(guó)也在積極開(kāi)展相關(guān)的研制工作[2]。

    然而,在長(zhǎng)時(shí)間的深空飛行中,由于空間中太陽(yáng)活動(dòng)、宇宙射線、星際塵埃、元件故障等多種未知因素的影響,導(dǎo)航敏感器得到的探測(cè)信息中很容易存在野值,導(dǎo)致濾波結(jié)果變差,甚至發(fā)散,嚴(yán)重影響航天器在深空中的自主飛行安全。因此,需要研究基于野值檢測(cè)和信息修正的容錯(cuò)濾波方法在X射線脈沖星自主導(dǎo)航中的應(yīng)用。

    目前對(duì)X射線脈沖星自主導(dǎo)航容錯(cuò)濾波方法的研究仍相對(duì)較少。文獻(xiàn)[3]利用魯棒濾波修正星表誤差,但沒(méi)有考慮野值因素的影響。文獻(xiàn)[4] 利用脈沖星和捷聯(lián)慣導(dǎo)組成組合導(dǎo)航系統(tǒng),考慮了天體遮擋因素,利用聯(lián)邦EKF濾波進(jìn)行導(dǎo)航,但捷聯(lián)慣導(dǎo)誤差隨時(shí)間漂移,有應(yīng)用的局限性。本文采用X射線脈沖星和太陽(yáng)敏感器組成組合導(dǎo)航系統(tǒng),其誤差不隨時(shí)間漂移,利用基于殘差正交性的 UKF濾波方法(residual orthogonal unscented Kalman filter,ROUKF)進(jìn)行野值檢測(cè)。仿真結(jié)果表明,該方法在X射線脈沖星組合導(dǎo)航系統(tǒng)容錯(cuò)濾波中,能有效進(jìn)行野值診斷和修正,為未來(lái)的工程化應(yīng)用提供了借鑒。

    1 數(shù)學(xué)模型

    1.1 巡航段軌道動(dòng)力學(xué)模型

    在執(zhí)行深空任務(wù)時(shí),由于航天器在星際巡航段的飛行時(shí)間最長(zhǎng),遇到的復(fù)雜情況相對(duì)較多,因此選取深空巡航段的導(dǎo)航任務(wù)作為研究對(duì)象。

    在深空巡航段,常將坐標(biāo)系建立在太陽(yáng)質(zhì)心J2000慣性坐標(biāo)系中,除太陽(yáng)中心引力外,還要考慮大行星的攝動(dòng)引力和太陽(yáng)光壓作用。其軌道動(dòng)力學(xué)方程為[5]:

    式中:r和v分別為航天器在太陽(yáng)質(zhì)心J2000慣性坐標(biāo)系中的位置及速度;μs為太陽(yáng)引力常數(shù);np為攝動(dòng)引力的數(shù)目,考慮了太陽(yáng)系內(nèi)所有大行星的引力; μi為行星的引力常數(shù); rri為第i顆行星相對(duì)于航天器的位置; rpi為第i顆行星在太陽(yáng)質(zhì)心慣性坐標(biāo)系中的位置;A為航天器受到太陽(yáng)輻射壓的面積;G為太陽(yáng)光壓系數(shù);m為航天器質(zhì)量;a為未建模的其他攝動(dòng)力。

    導(dǎo)航系統(tǒng)中的狀態(tài)變量取為航天器在太陽(yáng)質(zhì)心J2000慣性坐標(biāo)系中的位置矢量和速度矢量,即

    1.2 X射線脈沖星觀測(cè)模型

    X射線脈沖星高穩(wěn)定的輻射周期和累積脈沖輪廓是高精度自主導(dǎo)航的基礎(chǔ)。X射線光子到達(dá)時(shí)間可以被記錄下來(lái),通過(guò)時(shí)間轉(zhuǎn)換將該時(shí)間轉(zhuǎn)化為光子到達(dá)太陽(yáng)系質(zhì)心的時(shí)間,通過(guò)歷元折疊可以得到某一歷元的累積脈沖輪廓。時(shí)間轉(zhuǎn)換方程為[6]:

    式中,tSSB為經(jīng)轉(zhuǎn)換后到達(dá)太陽(yáng)系質(zhì)心的時(shí)間,tSC為光子到達(dá)X射線探測(cè)器的時(shí)間,n為太陽(yáng)系質(zhì)心J2000坐標(biāo)系中脈沖星的方向矢量,α和λ為脈沖星的赤經(jīng)和赤緯, rsun為太陽(yáng)在太陽(yáng)系質(zhì)心坐標(biāo)系中的位置,c為光速, D0為太陽(yáng)系質(zhì)心到脈沖星的距離,b為太陽(yáng)系質(zhì)心在太陽(yáng)質(zhì)心J2000慣性坐標(biāo)系中的位置。

    通過(guò)脈沖星數(shù)據(jù)庫(kù),可以預(yù)測(cè)獲得同一歷元在太陽(yáng)系質(zhì)心處的累積脈沖輪廓,通過(guò)長(zhǎng)期觀測(cè),可以得到X射線脈沖星的相位模型:

    經(jīng)輪廓互相關(guān)處理可以得到一個(gè)時(shí)間延遲量,該時(shí)間延遲量即為X射線脈沖星導(dǎo)航的基本觀測(cè)量 Δt。

    X射線脈沖星導(dǎo)航的測(cè)量方程為:

    式中,ni=[nixniyniz]T為第i顆X射線脈沖星在太陽(yáng)系質(zhì)心J2000坐標(biāo)系中的方向矢量,i為導(dǎo)航所使用的脈沖星個(gè)數(shù),w1(k)為測(cè)量噪聲。

    1.3 太陽(yáng)敏感器觀測(cè)模型

    由于在巡航段中,幾乎不存在天體遮擋因素,航天器能夠始終獲得太陽(yáng)光線信息,因此在姿態(tài)已知的情況下,用太陽(yáng)敏感器能夠獲得太陽(yáng)相對(duì)于航天器的方向矢量,其測(cè)量模型為[7]:

    式中, ls為太陽(yáng)敏感器的基本觀測(cè)量——太陽(yáng)質(zhì)心相對(duì)于航天器的方向矢量,r為航天器在太陽(yáng)質(zhì)心慣性坐標(biāo)系中的位置。

    則太陽(yáng)敏感器的測(cè)量方程為:

    式中, w2(k)為測(cè)量噪聲。

    2 組合導(dǎo)航系統(tǒng)設(shè)計(jì)

    無(wú)跡卡爾曼濾波(Unscented Kalman Filter,UKF)在解決非線性問(wèn)題時(shí),不需要進(jìn)行狀態(tài)方程和測(cè)量方程的線性化,因此具有更高的精度。

    設(shè)系統(tǒng)的狀態(tài)方程和測(cè)量方程如下:

    則標(biāo)準(zhǔn)的UKF算法如下[8]:

    在UKF濾波器中,該故障檢測(cè)函數(shù)為[10]:

    則式(18)可以修正為:

    式中, f(r)為加權(quán)因子。當(dāng) γk≤ ε時(shí), f(r)=1;當(dāng)γk> ε時(shí), f(r) = 1/γk。在統(tǒng)計(jì)學(xué)上γk服從自由度為觀測(cè)維數(shù)的 χ2分布,ε為在顯著性水平 0.01的情況下 χ2分布的臨界值。

    組合導(dǎo)航系統(tǒng)分為X射線脈沖星導(dǎo)航子系統(tǒng)和太陽(yáng)敏感器導(dǎo)航子系統(tǒng),每個(gè)子濾波器均采用容錯(cuò)UKF算法進(jìn)行野值檢測(cè)和濾波估計(jì),兩個(gè)子濾波器的濾波結(jié)果最終利用聯(lián)邦濾波方法進(jìn)行信息融合。關(guān)于聯(lián)邦濾波器的實(shí)現(xiàn)可參見(jiàn)相關(guān)資料。

    3 數(shù)學(xué)仿真

    3.1 仿真條件

    1)以美國(guó)“深空一號(hào)”任務(wù)巡航段的軌道驗(yàn)證本文算法。航天器在太陽(yáng)質(zhì)心慣性坐標(biāo)系中的位置初值為(126 966 202.276 870, 71 240 023.368 644, 308 951 743.957 506) km ,速度初值為(-17.012 246 566 619, 25.573 125 836 158, 11.411 474 895 983) km/s 。初始時(shí)刻在各方向上的位置誤差為500 km,各方向上的速度誤差為1.5 m/s;

    3)導(dǎo)航用的脈沖星為B0531+21、B1937+ 21、B1957+20、B1821-24,探測(cè)器面積為1 m2,脈沖信號(hào)累積時(shí)間為500 s,測(cè)量更新時(shí)間為500 s,總仿真時(shí)間為30 d。

    4)為檢驗(yàn)該組合導(dǎo)航系統(tǒng)的容錯(cuò)性能,不但對(duì)離散野值進(jìn)行了考察,而且對(duì)連續(xù)野值進(jìn)行了仿真。在離散野值仿真中,分別在4×105s、8×105s、1.2×106s、1.6×106s處的脈沖星導(dǎo)航測(cè)量值中加入 1000 km野值,在5×105s、1×106s、1.5×106s、2×106s處的太陽(yáng)敏感器測(cè)量值中加入0.5 rad野值。在連續(xù)野值仿真中,分別在400 500~410 000 s、1 000 500~1 010 000 s、 2 000 500~2 010 000 s處的脈沖星導(dǎo)航測(cè)量值中連續(xù)加入1000 km野值,在500 500~507 500 s、1 500 500~1 507 500 s處的太陽(yáng)敏感器測(cè)量值中連續(xù)加入0.5 rad野值。

    3.2 仿真結(jié)果

    圖1給出了無(wú)野值情況下,UKF濾波和基于殘差正交性容錯(cuò)UKF濾波的組合導(dǎo)航位置誤差,兩種濾波均能有效收斂。UKF濾波位置精度為0.184 321 km,速度精度為2.569 87×10-6km/s;容錯(cuò)UKF濾波方法的位置精度為0.184 475 km,速度精度為2.570 16×10-6km/s。說(shuō)明在無(wú)野值情況下,容錯(cuò)UKF濾波器與UKF濾波器估計(jì)結(jié)果相當(dāng),都具有較高的導(dǎo)航精度。

    圖 1 UKF濾波和容錯(cuò)UKF濾波對(duì)比(無(wú)野值)Fig.1 UKF and ROUKF results without outliers

    圖2和圖3分別給出了在離散野值和連續(xù)野值條件下,UKF濾波和容錯(cuò)UKF濾波的組合導(dǎo)航位置誤差。由圖2、圖3可知,在較大野值存在時(shí),UKF濾波的估計(jì)效果較差,并且在連續(xù)野值情況下位置誤差更大,在仿真條件下達(dá)到上千千米的位置誤差,而容錯(cuò)UKF濾波能夠有效檢測(cè)野值,并進(jìn)行修正,使濾波結(jié)果始終保持較高的導(dǎo)航精度,具有較強(qiáng)的容錯(cuò)性能。

    圖 2 UKF濾波和容錯(cuò)UKF濾波對(duì)比(離散野值)Fig.2 UKF and ROUKF results with discrete outliers

    圖3 UKF濾波和容錯(cuò)UKF濾波對(duì)比(連續(xù)野值)Fig.3 UKF and ROUKF results with consecutive outliers

    4 結(jié) 論

    本文針對(duì)X射線脈沖星組合導(dǎo)航中出現(xiàn)的野值問(wèn)題,利用基于殘差正交性的UKF容錯(cuò)算法,進(jìn)行野值檢測(cè)和修正,并通過(guò)仿真實(shí)驗(yàn)進(jìn)行了驗(yàn)證。結(jié)果表明,在X射線脈沖星組合導(dǎo)航系統(tǒng)中應(yīng)用ROUKF算法,能夠有效識(shí)別野值,始終保持較高的定位精度,避免了野值存在時(shí)的導(dǎo)航誤差過(guò)大問(wèn)題。

    (References):

    [1] Keith C G, Zaven A, Takashi O, et al. The Neutron star Interior Composition ExploreR (NICER): an Explorer mission of opportunity for soft x-ray timing spectroscopy [C]//SPIE Proceedings, Bellingham: SPIE, 2012: 13-20.

    [2] Fuchang Zuo, Jianwu Chen, Liansheng Li, et al. Analysis and design of grazing incidence x-ray optics for pulsar navigation[C]//SPIE Proceedings, Bellingham, 2013: 23-30.

    [3] 熊凱,魏春嶺,劉良棟. 魯棒濾波技術(shù)在脈沖星導(dǎo)航中的應(yīng)用[J]. 空間控制技術(shù)與應(yīng)用,2008,34(6):8-12.

    XIONG Kai, WEI Chun-ling, LIU Liang-dong. Application of robust filtering in pulsars based navigation[J]. Aerospace Control and Application, 2008, 34(6): 8-12.

    [4] 蘇哲,許錄平,張華,等. 基于XPNAV和SINS的容錯(cuò)組合導(dǎo)航系統(tǒng)[J]. 華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,39(6):24-28.

    SU Zhe, XU Lu-ping, ZHANG hua, et al. Fault-tolerant integrated navigation system based on XPNAV and SINS[J]. Journal of Huazhong University of Science and Technology(Nature Science), 2011, 39(6): 24-28.

    [5] 楊成偉,鄭建華. 基于脈沖星和小行星的組合導(dǎo)航在深空巡航段的應(yīng)用[J]. 中國(guó)慣性技術(shù)學(xué)報(bào),2012,20(5):583-586.

    YANG Cheng-wei, ZHENG Jian-hua. Integrated navigation based on X-ray pulsars and asteroids during interplanetary cruise[J]. Journal of Chinese Inertial Technology, 2012, 20(5): 583-586.

    [6] Li Jianxun, Ke Xizheng. Study on autonomous navigation based on pulsar timing model[J]. Science in China Series G: Physics, Mechanics, Astronomy, 2009, 52(2): 303-309.

    [7] Ke Han, Hao Wang, Binjie Tu, et al. Pico-satellite autonomous navigation with magnetometer and sun sensor data [J]. Chinese Journal of Aeronautics, 2011, 24(1): 46-54.

    [8] Jwo DahJing, Yang Chifan, Chuang Chihhsun, et al. Performance enhancement for ultra-tight GPS/INS integration using a fuzzy adaptive strong tracking unscented Kalman filter[J]. Nonlinear Dynamics, 2013, 73(1): 377-395.

    [9] 榮思遠(yuǎn),常亞武,崔乃剛. UKF容錯(cuò)濾波方法在自主導(dǎo)航中的應(yīng)用研究[J]. 宇航學(xué)報(bào),2006,27(4):605-608.

    RONG Si-yuan, CHANG Ya-wu, CUI Nai-gang. Research on the fault tolerant UKF applied on autonomous navigation[J]. Journal of Astronautics, 2006, 27(4): 605-608.

    UKF容錯(cuò)濾波在脈沖星組合導(dǎo)航中的應(yīng)用

    楊成偉1,鄭建華2,高 東2

    Fault tolerant UKF application in pulsar-based integrated navigation

    YANG Cheng-wei1, ZHENG Jian-hua2, GAO Dong2
    (1. School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, China; 2. Center for Space Science and Applied Research, Chinese Academy of Sciences, Beijing 100190, China)

    In view of the low navigation accuracy caused by the outliers in the pulsar-based integrated navigation system, a residual orthogonal unscented Kalman filter (ROUKF) was employed to correct the outliers. During interplanetary cruise, an integrated navigation was applied which included X-ray pulsar navigation system and sun observation navigation system. The X-ray detector was used to record the arrival times of pulsar photons. The sun sensor was used to measure the sun line-of-sight vector. The information fusion was conducted by the federated filter. Simulation results demonstrate that the ROUKF in the pulsar-based integrated navigation system can effectively correct the outliers. The proposed method improves the robust capability and the application value of the pulsar-based integrated navigation system.

    pulsar; sun sensor; unscented Kalman filter; outlier detection; fault tolerant filter

    1005-6734(2014)06-0759-04

    10.13695/j.cnki.12-1222/o3.2014.06.011

    V448.224

    A

    2014-06-08;

    2014-10-27

    空間科學(xué)預(yù)先研究項(xiàng)目(XDA04074300)

    楊成偉(1986—),男,博士后,從事自主導(dǎo)航研究。E-mail:yangchengwei2009@126.com

    猜你喜歡
    野值脈沖星質(zhì)心
    “中國(guó)天眼”已發(fā)現(xiàn)740余顆新脈沖星
    軍事文摘(2023年12期)2023-06-12 07:51:00
    發(fā)現(xiàn)脈沖星的女天文學(xué)家——貝爾
    科學(xué)(2022年4期)2022-10-25 02:43:42
    重型半掛汽車(chē)質(zhì)量與質(zhì)心位置估計(jì)
    基于回歸模型的采集數(shù)據(jù)清洗技術(shù)
    電光與控制(2022年4期)2022-04-07 12:47:06
    基于GNSS測(cè)量的天宮二號(hào)質(zhì)心確定
    彈道跟蹤數(shù)據(jù)野值剔除方法性能分析
    上海航天(2018年4期)2018-09-07 03:28:46
    導(dǎo)引頭在線野值剔除與平滑方法研究
    RS422總線數(shù)據(jù)野值剔除方法研究
    基于虛擬觀測(cè)值的X射線單脈沖星星光組合導(dǎo)航
    一種海洋測(cè)高衛(wèi)星質(zhì)心在軌估計(jì)算法
    航天器工程(2014年5期)2014-03-11 16:35:53
    亚洲国产精品sss在线观看| 亚洲国产精品成人综合色| 亚洲av日韩精品久久久久久密| АⅤ资源中文在线天堂| 亚洲av日韩精品久久久久久密| 制服诱惑二区| 欧美一级a爱片免费观看看 | 久久久久久久午夜电影| 欧美成人一区二区免费高清观看 | 国产精品乱码一区二三区的特点| 午夜免费成人在线视频| 久久精品人妻少妇| 波多野结衣高清无吗| 久久这里只有精品中国| 狠狠狠狠99中文字幕| 99久久精品国产亚洲精品| 中文在线观看免费www的网站 | 久9热在线精品视频| 制服人妻中文乱码| 精品国产亚洲在线| 中文字幕人妻丝袜一区二区| 在线观看免费日韩欧美大片| 中文字幕最新亚洲高清| 欧美三级亚洲精品| 波多野结衣高清作品| 亚洲五月婷婷丁香| 亚洲av电影不卡..在线观看| 成人高潮视频无遮挡免费网站| 国产男靠女视频免费网站| 亚洲中文av在线| 最近在线观看免费完整版| 两个人看的免费小视频| 欧美日韩黄片免| 99热只有精品国产| 午夜福利视频1000在线观看| 制服人妻中文乱码| 欧美成人午夜精品| 国产成年人精品一区二区| 99riav亚洲国产免费| 丰满人妻熟妇乱又伦精品不卡| 又黄又爽又免费观看的视频| avwww免费| 岛国在线免费视频观看| 国产久久久一区二区三区| 欧美精品啪啪一区二区三区| 国产在线观看jvid| 免费人成视频x8x8入口观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av五月六月丁香网| 宅男免费午夜| 1024手机看黄色片| 看片在线看免费视频| 午夜视频精品福利| 中文字幕熟女人妻在线| 少妇的丰满在线观看| 欧美乱色亚洲激情| 日韩精品免费视频一区二区三区| 夜夜爽天天搞| 91大片在线观看| 亚洲成人免费电影在线观看| 深夜精品福利| 18禁黄网站禁片午夜丰满| 免费一级毛片在线播放高清视频| 亚洲精品在线观看二区| 一夜夜www| 亚洲欧美精品综合久久99| 欧美 亚洲 国产 日韩一| 少妇人妻一区二区三区视频| 亚洲av第一区精品v没综合| 免费看日本二区| 欧美性长视频在线观看| 两个人免费观看高清视频| 午夜福利在线观看吧| 国产片内射在线| 国语自产精品视频在线第100页| 此物有八面人人有两片| 久久婷婷人人爽人人干人人爱| 亚洲成人国产一区在线观看| АⅤ资源中文在线天堂| 免费看美女性在线毛片视频| www.www免费av| 可以在线观看的亚洲视频| 日韩欧美在线二视频| 麻豆一二三区av精品| 国内少妇人妻偷人精品xxx网站 | 欧美日韩福利视频一区二区| 精品高清国产在线一区| 亚洲专区中文字幕在线| 久久人妻福利社区极品人妻图片| 欧美3d第一页| 婷婷精品国产亚洲av| 少妇被粗大的猛进出69影院| 午夜福利视频1000在线观看| 久久 成人 亚洲| 波多野结衣高清无吗| 精品国产乱码久久久久久男人| 欧美中文日本在线观看视频| 亚洲色图av天堂| 久久99热这里只有精品18| 一级片免费观看大全| 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区视频在线观看免费| 国产精品自产拍在线观看55亚洲| 99久久综合精品五月天人人| 777久久人妻少妇嫩草av网站| 丰满的人妻完整版| 久久久国产成人免费| 国产激情偷乱视频一区二区| 国产亚洲精品第一综合不卡| 69av精品久久久久久| 国产av一区在线观看免费| 亚洲 欧美 日韩 在线 免费| 国产一区二区在线av高清观看| 久久这里只有精品19| 无限看片的www在线观看| 国产不卡一卡二| 免费看十八禁软件| 国产成人精品久久二区二区免费| 天天躁狠狠躁夜夜躁狠狠躁| av超薄肉色丝袜交足视频| 国产在线精品亚洲第一网站| 日本免费一区二区三区高清不卡| 亚洲九九香蕉| 亚洲国产精品999在线| 舔av片在线| 欧美一级毛片孕妇| 久久久久久免费高清国产稀缺| 黄色视频不卡| 亚洲成a人片在线一区二区| 757午夜福利合集在线观看| 国产精品av视频在线免费观看| 日本一二三区视频观看| 黑人操中国人逼视频| 少妇人妻一区二区三区视频| 中亚洲国语对白在线视频| 久久热在线av| 亚洲精品一区av在线观看| 欧美成人性av电影在线观看| 九色国产91popny在线| 12—13女人毛片做爰片一| 嫩草影院精品99| 91麻豆av在线| 国产精品久久久久久久电影 | 欧美午夜高清在线| 99在线视频只有这里精品首页| 午夜日韩欧美国产| 黄色女人牲交| 最近在线观看免费完整版| 99久久99久久久精品蜜桃| 亚洲成人中文字幕在线播放| 在线观看免费日韩欧美大片| 国产亚洲精品av在线| 国产真人三级小视频在线观看| 亚洲黑人精品在线| 又大又爽又粗| 99久久精品国产亚洲精品| 特级一级黄色大片| 国产真人三级小视频在线观看| 精品免费久久久久久久清纯| 成人国语在线视频| 成人精品一区二区免费| 老熟妇仑乱视频hdxx| 波多野结衣高清作品| 国产精品 国内视频| 久久草成人影院| 好看av亚洲va欧美ⅴa在| 日韩欧美 国产精品| 亚洲,欧美精品.| 最近在线观看免费完整版| 欧美极品一区二区三区四区| 免费在线观看日本一区| 欧美不卡视频在线免费观看 | 禁无遮挡网站| √禁漫天堂资源中文www| 国产乱人伦免费视频| 婷婷亚洲欧美| 成人av一区二区三区在线看| 一区二区三区高清视频在线| 波多野结衣高清无吗| 亚洲片人在线观看| 国产午夜精品久久久久久| www.www免费av| 桃色一区二区三区在线观看| 日日摸夜夜添夜夜添小说| 国产亚洲精品久久久久久毛片| 成人三级黄色视频| 嫩草影院精品99| 一进一出抽搐动态| 91老司机精品| 妹子高潮喷水视频| 欧美黄色片欧美黄色片| 日本 欧美在线| bbb黄色大片| 国产精品久久电影中文字幕| 国产成人av激情在线播放| 国产精品国产高清国产av| 女人高潮潮喷娇喘18禁视频| 国产一区二区在线av高清观看| 18禁裸乳无遮挡免费网站照片| 欧美另类亚洲清纯唯美| 老司机深夜福利视频在线观看| 此物有八面人人有两片| 九九热线精品视视频播放| av片东京热男人的天堂| 无限看片的www在线观看| 欧美日本亚洲视频在线播放| 国产亚洲欧美在线一区二区| 精品欧美国产一区二区三| 午夜久久久久精精品| 国产精品一区二区三区四区免费观看 | 亚洲av熟女| 久久久久性生活片| 夜夜躁狠狠躁天天躁| 精品国产亚洲在线| 深夜精品福利| 真人一进一出gif抽搐免费| 亚洲在线自拍视频| 国产一级毛片七仙女欲春2| 欧美黑人欧美精品刺激| 成年人黄色毛片网站| 成人一区二区视频在线观看| 欧美三级亚洲精品| 搡老岳熟女国产| 久久久久久久久久黄片| av视频在线观看入口| 久9热在线精品视频| 精品久久久久久久久久久久久| 国内精品一区二区在线观看| 亚洲欧美日韩无卡精品| 亚洲狠狠婷婷综合久久图片| 可以在线观看的亚洲视频| 真人做人爱边吃奶动态| www.自偷自拍.com| 欧美一区二区国产精品久久精品 | 免费一级毛片在线播放高清视频| 一进一出好大好爽视频| 国产成人欧美在线观看| 伊人久久大香线蕉亚洲五| 精品久久久久久久末码| 色老头精品视频在线观看| videosex国产| 中文字幕高清在线视频| 日韩中文字幕欧美一区二区| 91麻豆精品激情在线观看国产| 久久久水蜜桃国产精品网| 久久天躁狠狠躁夜夜2o2o| 亚洲免费av在线视频| 好看av亚洲va欧美ⅴa在| 国产野战对白在线观看| 九九热线精品视视频播放| 首页视频小说图片口味搜索| 国产精品综合久久久久久久免费| 久久天躁狠狠躁夜夜2o2o| 香蕉久久夜色| a级毛片a级免费在线| 人成视频在线观看免费观看| 国产欧美日韩一区二区三| 久久久久久大精品| 精品熟女少妇八av免费久了| 一本大道久久a久久精品| 国产主播在线观看一区二区| 麻豆成人午夜福利视频| 欧美乱码精品一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 国产真实乱freesex| 一本综合久久免费| 亚洲,欧美精品.| 亚洲成av人片免费观看| 成人永久免费在线观看视频| 宅男免费午夜| 国产日本99.免费观看| 久久精品亚洲精品国产色婷小说| 村上凉子中文字幕在线| 国产伦人伦偷精品视频| 欧美日韩亚洲国产一区二区在线观看| 国产一区二区三区视频了| 国内毛片毛片毛片毛片毛片| 成人18禁高潮啪啪吃奶动态图| 麻豆久久精品国产亚洲av| 麻豆国产97在线/欧美 | 搞女人的毛片| 国产主播在线观看一区二区| 女警被强在线播放| 婷婷精品国产亚洲av| 长腿黑丝高跟| 热99re8久久精品国产| 国产成+人综合+亚洲专区| 久久亚洲精品不卡| 俺也久久电影网| 男男h啪啪无遮挡| 亚洲国产看品久久| 18禁裸乳无遮挡免费网站照片| 丁香六月欧美| 此物有八面人人有两片| 听说在线观看完整版免费高清| 在线看三级毛片| 免费无遮挡裸体视频| 变态另类丝袜制服| 91麻豆av在线| 亚洲中文字幕一区二区三区有码在线看 | 在线a可以看的网站| 一a级毛片在线观看| 在线观看免费午夜福利视频| 亚洲一区中文字幕在线| 十八禁网站免费在线| 99久久综合精品五月天人人| 老汉色∧v一级毛片| 欧美日韩乱码在线| 国产精品 欧美亚洲| 国产精品亚洲一级av第二区| 免费高清视频大片| 桃色一区二区三区在线观看| 免费无遮挡裸体视频| 深夜精品福利| 99久久精品国产亚洲精品| 亚洲一区二区三区不卡视频| 香蕉av资源在线| 久久久久久亚洲精品国产蜜桃av| 亚洲成人免费电影在线观看| av在线天堂中文字幕| 法律面前人人平等表现在哪些方面| 全区人妻精品视频| 哪里可以看免费的av片| 国产日本99.免费观看| 久热爱精品视频在线9| 久久久水蜜桃国产精品网| 国产区一区二久久| 日韩国内少妇激情av| 久久精品国产清高在天天线| 全区人妻精品视频| 日韩欧美三级三区| 午夜激情福利司机影院| 亚洲av五月六月丁香网| 九色成人免费人妻av| 高清在线国产一区| 午夜免费成人在线视频| 毛片女人毛片| 可以在线观看毛片的网站| 2021天堂中文幕一二区在线观| 性色av乱码一区二区三区2| 亚洲黑人精品在线| 美女扒开内裤让男人捅视频| 老司机在亚洲福利影院| 长腿黑丝高跟| 制服诱惑二区| 欧美性长视频在线观看| 亚洲,欧美精品.| 亚洲欧美激情综合另类| 岛国在线观看网站| 久久天堂一区二区三区四区| 亚洲五月婷婷丁香| 一级毛片高清免费大全| av天堂在线播放| 欧美黑人巨大hd| 波多野结衣巨乳人妻| 午夜视频精品福利| 哪里可以看免费的av片| 热99re8久久精品国产| 色老头精品视频在线观看| 国产精品av视频在线免费观看| 免费搜索国产男女视频| 九色国产91popny在线| 此物有八面人人有两片| 精品免费久久久久久久清纯| 正在播放国产对白刺激| 999久久久国产精品视频| 51午夜福利影视在线观看| 国内精品久久久久精免费| 少妇人妻一区二区三区视频| 日韩av在线大香蕉| 久久久精品大字幕| 国产精品久久久久久久电影 | 亚洲国产精品久久男人天堂| 国产一级毛片七仙女欲春2| 欧美性猛交╳xxx乱大交人| 久久精品国产99精品国产亚洲性色| av免费在线观看网站| 国产黄a三级三级三级人| 成年免费大片在线观看| 亚洲 欧美 日韩 在线 免费| 久久精品国产亚洲av高清一级| 精品国内亚洲2022精品成人| 久久这里只有精品中国| 亚洲成a人片在线一区二区| 精品一区二区三区视频在线观看免费| 午夜亚洲福利在线播放| 怎么达到女性高潮| 在线观看一区二区三区| 桃色一区二区三区在线观看| 人妻夜夜爽99麻豆av| 精品久久久久久久末码| 成人精品一区二区免费| 亚洲成人国产一区在线观看| 国产区一区二久久| 国产午夜福利久久久久久| 精品免费久久久久久久清纯| 精品久久久久久久久久免费视频| 丁香欧美五月| 国产成人av教育| 精品国产超薄肉色丝袜足j| 亚洲黑人精品在线| 亚洲午夜理论影院| 香蕉av资源在线| 亚洲男人的天堂狠狠| 无遮挡黄片免费观看| 免费在线观看亚洲国产| 熟女电影av网| 欧美一区二区国产精品久久精品 | 曰老女人黄片| 九色成人免费人妻av| 午夜成年电影在线免费观看| 动漫黄色视频在线观看| 九色成人免费人妻av| 国产高清videossex| 日本五十路高清| 在线观看免费日韩欧美大片| 在线观看舔阴道视频| 久久天堂一区二区三区四区| 别揉我奶头~嗯~啊~动态视频| 国产精品98久久久久久宅男小说| 丝袜人妻中文字幕| 老司机靠b影院| 狂野欧美白嫩少妇大欣赏| 成人国产一区最新在线观看| 久久久久性生活片| 久久久久久九九精品二区国产 | 少妇粗大呻吟视频| 久久这里只有精品19| 久久国产乱子伦精品免费另类| 成人亚洲精品av一区二区| 丰满的人妻完整版| 免费在线观看日本一区| 麻豆国产97在线/欧美 | 亚洲18禁久久av| 国产1区2区3区精品| 叶爱在线成人免费视频播放| 欧美不卡视频在线免费观看 | 亚洲一区二区三区色噜噜| 日本一区二区免费在线视频| 美女大奶头视频| 成熟少妇高潮喷水视频| 欧美黄色片欧美黄色片| 午夜日韩欧美国产| а√天堂www在线а√下载| 757午夜福利合集在线观看| 国产成人精品无人区| 变态另类成人亚洲欧美熟女| 一级作爱视频免费观看| 亚洲国产精品sss在线观看| 日韩欧美在线二视频| av天堂在线播放| 热99re8久久精品国产| 亚洲国产精品999在线| 成人三级做爰电影| 午夜a级毛片| 国产人伦9x9x在线观看| 天堂√8在线中文| 色综合婷婷激情| 丝袜美腿诱惑在线| 婷婷精品国产亚洲av| 一本一本综合久久| 2021天堂中文幕一二区在线观| 亚洲精品国产一区二区精华液| 黄片小视频在线播放| 一级毛片高清免费大全| 国产成人精品无人区| 一区二区三区激情视频| 免费在线观看亚洲国产| 久久精品aⅴ一区二区三区四区| 淫妇啪啪啪对白视频| or卡值多少钱| 成年版毛片免费区| 一级作爱视频免费观看| 一级a爱片免费观看的视频| 国产人伦9x9x在线观看| 九九热线精品视视频播放| 少妇熟女aⅴ在线视频| 久久久久国内视频| 日韩欧美国产一区二区入口| 蜜桃久久精品国产亚洲av| 国产高清视频在线观看网站| 国产午夜精品论理片| 男女视频在线观看网站免费 | 性欧美人与动物交配| 在线永久观看黄色视频| 波多野结衣巨乳人妻| 亚洲一码二码三码区别大吗| 亚洲精品在线美女| 桃色一区二区三区在线观看| 国产欧美日韩一区二区三| 精品国内亚洲2022精品成人| 亚洲av第一区精品v没综合| 国产亚洲av高清不卡| 欧美日韩一级在线毛片| av国产免费在线观看| 黑人欧美特级aaaaaa片| av在线天堂中文字幕| 三级男女做爰猛烈吃奶摸视频| 精品熟女少妇八av免费久了| 国产一区在线观看成人免费| 免费在线观看日本一区| а√天堂www在线а√下载| 久久精品91蜜桃| xxxwww97欧美| 日日撸夜夜添| 午夜免费男女啪啪视频观看| 欧美成人a在线观看| 日本免费一区二区三区高清不卡| 午夜精品国产一区二区电影 | 国产真实伦视频高清在线观看| 欧美成人a在线观看| 99视频精品全部免费 在线| 麻豆国产av国片精品| 精品人妻偷拍中文字幕| 99热这里只有是精品在线观看| 欧美日韩在线观看h| 秋霞在线观看毛片| 亚洲乱码一区二区免费版| 国产精品一二三区在线看| 日本色播在线视频| 亚洲欧美精品综合久久99| 美女大奶头视频| 99热这里只有精品一区| 久久精品夜夜夜夜夜久久蜜豆| 午夜福利成人在线免费观看| 日日撸夜夜添| 精品国内亚洲2022精品成人| 听说在线观看完整版免费高清| 国产精品国产高清国产av| videossex国产| 天堂中文最新版在线下载 | 男人狂女人下面高潮的视频| 婷婷色综合大香蕉| 欧美另类亚洲清纯唯美| 精品一区二区免费观看| 国产色爽女视频免费观看| 亚洲自偷自拍三级| 午夜视频国产福利| 亚洲熟妇中文字幕五十中出| 亚洲欧美成人精品一区二区| 亚洲av熟女| 亚洲中文字幕日韩| 午夜爱爱视频在线播放| 男人舔奶头视频| 国产精品久久久久久久久免| 久久精品久久久久久噜噜老黄 | 亚洲av中文字字幕乱码综合| 乱人视频在线观看| 九草在线视频观看| 美女内射精品一级片tv| 日本黄色视频三级网站网址| 久久国内精品自在自线图片| 日韩欧美 国产精品| 精品无人区乱码1区二区| 亚洲真实伦在线观看| 嘟嘟电影网在线观看| 久久精品久久久久久噜噜老黄 | 99久久成人亚洲精品观看| 久久久久久久久久黄片| 免费人成在线观看视频色| 国产精品爽爽va在线观看网站| 欧美日韩国产亚洲二区| 日日摸夜夜添夜夜爱| 国产精品一区二区在线观看99 | 欧美性猛交黑人性爽| 亚洲成人av在线免费| 青春草国产在线视频 | 日本在线视频免费播放| 色综合站精品国产| 国产亚洲精品久久久久久毛片| 夜夜夜夜夜久久久久| 精品午夜福利在线看| 九草在线视频观看| 午夜激情欧美在线| 亚洲一级一片aⅴ在线观看| 不卡视频在线观看欧美| 精品人妻熟女av久视频| 国产成人91sexporn| 欧美人与善性xxx| 久久久久网色| 久久婷婷人人爽人人干人人爱| 午夜免费男女啪啪视频观看| 春色校园在线视频观看| 欧美日韩精品成人综合77777| 看片在线看免费视频| 一边摸一边抽搐一进一小说| 精品人妻视频免费看| 国产欧美日韩精品一区二区| 成人无遮挡网站| 亚洲乱码一区二区免费版| 免费观看精品视频网站| 国产成人a∨麻豆精品| av在线观看视频网站免费| 亚洲欧美日韩高清在线视频| 亚洲欧美日韩卡通动漫| 神马国产精品三级电影在线观看| 国产高清三级在线| 亚洲国产精品成人综合色| 99国产极品粉嫩在线观看| 亚洲在线自拍视频| 天天躁日日操中文字幕| av天堂在线播放| 日韩视频在线欧美| 美女内射精品一级片tv| 国产一区二区三区av在线 | 嫩草影院新地址| 欧美+亚洲+日韩+国产| 又黄又爽又刺激的免费视频.| 日日摸夜夜添夜夜爱| 日产精品乱码卡一卡2卡三| 国内精品久久久久精免费| 国产精品久久视频播放| 美女大奶头视频| 久久久久久久久中文| 国产精品嫩草影院av在线观看|