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

    非合作航天器相對(duì)運(yùn)動(dòng)估計(jì)的混合濾波器設(shè)計(jì)方法

    2023-06-23 09:50:06張世源
    宇航學(xué)報(bào) 2023年5期
    關(guān)鍵詞:協(xié)方差航天器姿態(tài)

    盧 山,張世源

    (1. 上海航天控制技術(shù)研究所,上海 201109;2. 上海市空間智能控制技術(shù)重點(diǎn)實(shí)驗(yàn)室,上海 201109)

    0 引 言

    近年來(lái),中國(guó)航天領(lǐng)域的在軌服務(wù)任務(wù)日趨頻繁,對(duì)于空間非合作目標(biāo),如故障或失效的飛行器、空間垃圾等的相對(duì)導(dǎo)航技術(shù)研究是當(dāng)下的熱點(diǎn)。對(duì)非合作目標(biāo)的相對(duì)姿態(tài)、相對(duì)位置等信息進(jìn)行精確的濾波估計(jì),是實(shí)現(xiàn)在軌服務(wù)任務(wù)的基礎(chǔ),直接影響控制和制導(dǎo)的精度[1]。

    為了體現(xiàn)追蹤航天器與目標(biāo)航天器間相對(duì)姿態(tài)與相對(duì)位置耦合的實(shí)時(shí)效應(yīng),提升相對(duì)狀態(tài)估計(jì)的精度,文獻(xiàn)[2-4]建立了非合作航天器相對(duì)運(yùn)動(dòng)的軌道-姿態(tài)的一體化濾波模型。這一類(lèi)模型直接將量測(cè)設(shè)備獲取的原始特征點(diǎn)的光學(xué)信息用于相對(duì)狀態(tài)的估計(jì),使?fàn)顟B(tài)變量維度過(guò)高(超過(guò)30維),在現(xiàn)有星上計(jì)算機(jī)的硬件能力下難以實(shí)現(xiàn)實(shí)時(shí)解算。本文考慮在非合作目標(biāo)質(zhì)心位置未知的情況下,將追蹤航天器上量測(cè)設(shè)備輸出的目標(biāo)航天器形心的信息引入濾波方程,建立了追蹤航天器與目標(biāo)航天器相對(duì)運(yùn)動(dòng)的軌道-姿態(tài)的一體化濾波模型。相較以往研究中提出的模型,該模型的形式更加簡(jiǎn)潔且便于計(jì)算,并具有以下特點(diǎn):狀態(tài)方程由軌道相對(duì)運(yùn)動(dòng)模型與姿態(tài)相對(duì)運(yùn)動(dòng)模型構(gòu)成,軌道相對(duì)運(yùn)動(dòng)模型為線性,姿態(tài)相對(duì)運(yùn)動(dòng)模型為非線性,且軌道與姿態(tài)相互解耦;觀測(cè)方程為非線性,且與軌道相關(guān)的狀態(tài)變量和與姿態(tài)相關(guān)的狀態(tài)變量相互耦合。

    利用上述模型進(jìn)行非合作航天器相對(duì)狀態(tài)量的估計(jì)屬于非線性濾波問(wèn)題。文獻(xiàn)[5]提出了一種混合階球面單形-徑向容積卡爾曼濾波器(Mixed-degree spherical simplex-radial cubature Kalman filter,MSSRCKF)。MSSRCKF結(jié)合了三階球面單形規(guī)則和五階徑向規(guī)則,使得算法能夠以接近容積卡爾曼濾波器(Cubature Kalman filter,CKF)[6]的計(jì)算復(fù)雜度獲得接近高階CKF(High-degree CKF,HCKF)[7]的五階精度,適用于高維度的強(qiáng)非線性系統(tǒng)。但由于本文建立的非合作航天器相對(duì)運(yùn)動(dòng)一體化濾波模型的狀態(tài)方程中包含部分線性項(xiàng),當(dāng)使用非線性濾波器對(duì)模型全局進(jìn)行處理時(shí),部分計(jì)算成本將被浪費(fèi)。文獻(xiàn)[8-11]針對(duì)狀態(tài)方程為非線性、觀測(cè)方程為線性的情況,將非線性濾波器的時(shí)間更新公式與線性濾波器的量測(cè)更新公式相結(jié)合,使濾波器能夠?qū)δP椭械木€性與非線性模塊分別采用相對(duì)應(yīng)的濾波方式,提高了計(jì)算效率。文獻(xiàn)[12]使用Rao-Blackwellized粒子濾波器對(duì)狀態(tài)方程部分狀態(tài)變量為線性、線性項(xiàng)與非線性項(xiàng)相互作用的系統(tǒng)模型進(jìn)行估計(jì),實(shí)現(xiàn)了非線性濾波算法采樣點(diǎn)數(shù)量的降維。但上述研究均無(wú)法用于解決狀態(tài)方程中線性項(xiàng)與非線性項(xiàng)相互解耦、觀測(cè)方程中線性項(xiàng)與非線性項(xiàng)相互耦合這一更為復(fù)雜的情況。針對(duì)上述問(wèn)題,本文通過(guò)在時(shí)間更新公式中加入線性項(xiàng)與非線性項(xiàng)之間互協(xié)方差矩陣的更新,將適用于線性系統(tǒng)的卡爾曼濾波器(Kalman filter,KF)[13]與MSSRCKF進(jìn)行了融合,提出了線性-非線性混合卡爾曼濾波器(Linear-nonlinear hybrid Kalman filter,L-NLKF)。該濾波器能夠?qū)σ粋€(gè)濾波方程中的線性與非線狀態(tài)變量分別采用不同的濾波方式進(jìn)行處理,在保證估計(jì)精度的同時(shí),顯著降低了計(jì)算量。

    1 航天器相對(duì)運(yùn)動(dòng)的軌道-姿態(tài)一體化濾波模型

    坐標(biāo)系定義如下:

    OI-XIYIZI:地心慣性坐標(biāo)系FI。原點(diǎn)OI為地心,XI位于赤道平面內(nèi),指向春分點(diǎn),ZI沿地球自轉(zhuǎn)軸方向,向上為正,YI軸與XI軸和ZI軸構(gòu)成右手直角坐標(biāo)系。

    OO-XOYOZO:第二軌道坐標(biāo)系FO。原點(diǎn)OO位于航天器質(zhì)心,ZO軸指向地心,XO軸在軌道平面內(nèi)垂直于ZO軸指向航天器飛行方向,YO軸與XO軸和ZO軸構(gòu)成右手直角坐標(biāo)系。根據(jù)該定義可以分別得到追蹤星第二軌道坐標(biāo)系FCO與目標(biāo)星第二軌道坐標(biāo)系FTO。

    OB-XBYBZB:本體坐標(biāo)系FB。原點(diǎn)OB位于航天器質(zhì)心,XB軸沿航天器主慣量軸指向前,ZB軸垂直于XB軸指向下,YB軸與XB軸和ZB軸構(gòu)成右手直角坐標(biāo)系。根據(jù)該定義可以分別得到追蹤星本體坐標(biāo)系FCB與目標(biāo)星第二軌道坐標(biāo)系FTB。

    1.1 軌道相對(duì)運(yùn)動(dòng)模型

    由于目標(biāo)航天器為非合作,其軌道角速度未知,故將軌道相對(duì)運(yùn)動(dòng)模型投影在坐標(biāo)系FCO中。假設(shè)追蹤星在近距離捕獲階段軌道半徑不發(fā)生較大改變,其軌道角速度n近似為常數(shù)。記坐標(biāo)系FCB中兩星質(zhì)心的相對(duì)位置矢量為ρ=[x,y,z]T,線性化后的目標(biāo)星質(zhì)心相對(duì)追蹤星質(zhì)心的相對(duì)軌道動(dòng)力學(xué)方程表示為

    (1)

    式中:wx,wy,wz為攝動(dòng)力和模型誤差引起的激勵(lì)噪聲。

    由于非合作航天器的質(zhì)心位置未知,而形心位置可以通過(guò)激光雷達(dá)數(shù)據(jù)建立的三維點(diǎn)云得到,故建立目標(biāo)星形心坐標(biāo)系FTF,其三軸指向與坐標(biāo)系FTB一致。將兩星的相對(duì)位置、相對(duì)速度以及目標(biāo)星質(zhì)心在坐標(biāo)系FTF下的坐標(biāo)(a,b,c)作為狀態(tài)量,建立狀態(tài)方程

    (2)

    (3)

    (4)

    (5)

    經(jīng)離散化后的狀態(tài)方程表示為

    Xk+1=ΦXk+Γwk

    (6)

    式中:Φ∈R9×9為系統(tǒng)的狀態(tài)轉(zhuǎn)移矩陣;Γ∈R9×3為系統(tǒng)的干擾矩陣。

    1.2 姿態(tài)相對(duì)運(yùn)動(dòng)模型

    (7)

    由單位四元數(shù)表示的相對(duì)姿態(tài)運(yùn)動(dòng)學(xué)方程為

    (8)

    ω=ωIC-ωIT

    (9)

    式中:ωIC和ωIT分別為追蹤星和目標(biāo)星相對(duì)于坐標(biāo)系FI的角速度。

    坐標(biāo)系FCB下追蹤星相對(duì)目標(biāo)星的姿態(tài)動(dòng)力學(xué)方程為

    (10)

    1.3 狀態(tài)空間模型

    1) 狀態(tài)方程:

    根據(jù)上述建立的相對(duì)運(yùn)動(dòng)模型,狀態(tài)向量x定義為

    (11)

    得到具有加性噪聲的連續(xù)非線性動(dòng)態(tài)系統(tǒng)

    (12)

    式中:w為零均值高斯白噪聲,協(xié)方差為Q;F(x)為關(guān)于狀態(tài)向量x的非線性函數(shù),可結(jié)合式(2)、(8)、(10)得到。

    易知,此處軌道與姿態(tài)的相關(guān)狀態(tài)量是解耦的。由于連續(xù)非線性動(dòng)態(tài)系統(tǒng)不適用于計(jì)算機(jī)數(shù)值仿真,本文采用四階龍格-庫(kù)塔法[15]對(duì)其進(jìn)行離散化處理。

    2) 觀測(cè)方程:

    對(duì)于非合作目標(biāo)的測(cè)量,星上量測(cè)設(shè)備輸出的是坐標(biāo)系FTF與單機(jī)測(cè)量系之間的相對(duì)位置與相對(duì)姿態(tài)信息。假設(shè)單機(jī)測(cè)量系與坐標(biāo)系FCB重合,用rzx表示目標(biāo)星質(zhì)心在坐標(biāo)系FTF中的坐標(biāo)(a,b,c);rld為激光雷達(dá)的測(cè)量值,由追蹤星質(zhì)心指向目標(biāo)星形心;rCT為兩星質(zhì)心連線(追蹤星指向目標(biāo)星)。上述矢量存在幾何關(guān)系如下

    (13)

    建立觀測(cè)方程

    (14)

    2 混合階球面單形-徑向卡爾曼濾波器

    考慮如下的離散非線性動(dòng)態(tài)系統(tǒng)

    (15)

    為了提高濾波精度以及降低計(jì)算復(fù)雜度,文獻(xiàn)[5]將3階球面單形積分準(zhǔn)則與5階徑向積分準(zhǔn)則相結(jié)合,提出了一種混合階球面單形—徑向容積卡爾曼濾波器(MSSRCKF)。MSSRCKF具有優(yōu)異的非線性濾波性能。在維度較高的非線性系統(tǒng)中,MSSRCKF濾波精度接近于具有五階泰勒精度的HCKF,而計(jì)算時(shí)間僅與CKF相當(dāng),適用于非合作航天器相對(duì)運(yùn)動(dòng)的軌道-姿態(tài)一體化濾波模型。

    MSSRCKF的求積節(jié)點(diǎn)設(shè)置為:

    (16)

    (17)

    對(duì)應(yīng)的權(quán)重系數(shù)為:

    (18)

    1) 時(shí)間更新

    (19)

    2) 量測(cè)更新

    (20)

    3 線性-非線性混合卡爾曼濾波器

    式(12)為本文建立的非合作航天器相對(duì)運(yùn)動(dòng)模型的狀態(tài)方程,其中F(x)由式(2)、(8)、(10)構(gòu)成。式(2)為描述航天器軌道相對(duì)運(yùn)動(dòng)的線性方程,式(8)和(10)為描述航天器姿態(tài)相對(duì)運(yùn)動(dòng)的非線性方程,且式(2)與式(8)、(10)中的狀態(tài)變量相互解耦。使用MSSRCKF進(jìn)行濾波估計(jì)時(shí),在式(19)所示的時(shí)間更新步驟中,該算法將對(duì)式(2)、(8)、(10)中包含的全部狀態(tài)變量按式(16)所示的求積節(jié)點(diǎn)設(shè)置要求進(jìn)行確定性采樣。

    確定性采樣法是對(duì)服從高斯分布的狀態(tài)變量經(jīng)過(guò)非線性變化后的分布情況進(jìn)行近似估計(jì)的方法[16]。由于式(2)為線性,服從高斯分布的狀態(tài)變量通過(guò)該方程后的分布仍服從高斯分布,其均值和方差可直接由線性卡爾曼濾波公式求得。使用MSSRCKF的確定性采樣法對(duì)式(2)進(jìn)行處理時(shí),并不會(huì)提高對(duì)應(yīng)狀態(tài)變量的估計(jì)精度,反而會(huì)對(duì)計(jì)算成本造成浪費(fèi)。

    針對(duì)上述問(wèn)題,本文將應(yīng)用于線性系統(tǒng)的KF與MSSRCKF進(jìn)行了融合,提出了線性-非線性混合卡爾曼濾波器(Linear-Nonlinear Hybrid Kalman Filter,L-NLKF)。

    (21)

    過(guò)程噪聲wk也相應(yīng)劃分為w1,k∈R9與w2,k∈R7,對(duì)應(yīng)的協(xié)方差矩陣分別為Q1,k∈R9×9與Q2,k∈R7×7。根據(jù)狀態(tài)變量的劃分情況以及協(xié)方差矩陣的性質(zhì),對(duì)協(xié)方差矩陣Pk進(jìn)行分塊

    (22)

    式中:P1,k∈R9×9為x1,k的協(xié)方差矩陣;P2,k∈R7×7為x2,k的協(xié)方差矩陣;P12,k∈R9×7為x1,k與x2,k的互協(xié)方差矩陣。

    (23)

    式中:f(·)由式(5)和(7)得到。

    由于L-NLKF主要針對(duì)軌道-姿態(tài)一體化濾波模型中狀態(tài)方程的特殊非線性形式進(jìn)行設(shè)計(jì),只需對(duì)MSSRCKF算法的時(shí)間更新方式進(jìn)行改進(jìn),故本節(jié)僅給出L-NLKF算法時(shí)間更新的計(jì)算步驟:

    1) 利用KF進(jìn)行線性狀態(tài)預(yù)測(cè)

    (24)

    2) 利用MSSRCKF進(jìn)行非線性狀態(tài)預(yù)測(cè)

    (25)

    3) 利用卡爾曼濾波器進(jìn)行x1,k協(xié)方差矩陣更新

    P1,k+1|k=ΦP1,kΦT+Q1,k

    (26)

    4) 利用MSSRCKF進(jìn)行x2,k協(xié)方差矩陣更新

    (27)

    5)x1,k與x2,k的互協(xié)方差矩陣的更新

    互協(xié)方差矩陣P12,k的更新是L-NLKF算法的關(guān)鍵,本文將給出較為詳細(xì)的推導(dǎo)過(guò)程。

    P12,k|k+1=E[(x1,k+1|k-Ex1,k+1|k)(x2,k+1|k-

    (28)

    (29)

    (30)

    (31)

    (32)

    通過(guò)這種方式便可將高斯分布的加權(quán)二重積分問(wèn)題轉(zhuǎn)換為如式(30)所示的高斯分布的加權(quán)一重積分,并能夠使用確定性采樣的數(shù)值方法進(jìn)行近似計(jì)算。式(30)中:ηi,k與ωi分別為求積節(jié)點(diǎn)與對(duì)應(yīng)的權(quán)重。若對(duì)式(30)中的求積節(jié)點(diǎn)ηi,k重新進(jìn)行求解,在增加了計(jì)算量的同時(shí),由于一部分求積節(jié)點(diǎn)將被矩陣Θ1與Θ2歸零,其中包含的分布信息將被浪費(fèi),且余下未被歸零的求積節(jié)點(diǎn)不再能夠完整的描述狀態(tài)量的分布情況,使得濾波精度下降。

    (33)

    對(duì)x1,k進(jìn)行確定性采樣時(shí),求積節(jié)點(diǎn)ζi,k應(yīng)滿足條件如下:

    (34)

    ④ζi,k應(yīng)滿足協(xié)方差等式

    (35)

    (36)

    其中條件③、④為使ζi,k的前二階矩與x1,k一致,以保證求積精度。根據(jù)上述條件,本文設(shè)置求積節(jié)點(diǎn)如式(37)所示。需要注意的是,滿足上述條件的求積節(jié)點(diǎn)不唯一,本文選擇其中的一種情況進(jìn)行設(shè)置。當(dāng)bi的解析表達(dá)式難以直接得到時(shí),可以將bi中元素設(shè)為未知數(shù),結(jié)合上述條件利用計(jì)算機(jī)求解多元非線性方程,得到bi的數(shù)值解。

    (37)

    6) 合并線性項(xiàng)與非線性項(xiàng)

    (38)

    (39)

    4 仿真校驗(yàn)

    表1 航天器參數(shù)設(shè)置Table 1 Spacecraft parameter setting

    濾波時(shí)長(zhǎng)設(shè)置為120 s,步長(zhǎng)為1 s;使用均方根誤差(RMSE)來(lái)各濾波算法的精度。

    (40)

    共進(jìn)行NM=100次蒙特卡洛仿真,并將相對(duì)姿態(tài)四元數(shù)的計(jì)算結(jié)果轉(zhuǎn)換為歐拉角,分別記為偏航角ψ、俯仰角θ、滾轉(zhuǎn)角φ。濾波器估計(jì)誤差對(duì)比如圖1-圖5所示。

    圖1 相對(duì)位置估計(jì)的RMSEFig.1 RMSE of the relative position estimation

    圖2 相對(duì)速度估計(jì)的RMSEFig.2 RMSE of the relative velocity estimation

    圖4 相對(duì)姿態(tài)估計(jì)的RMSEFig.4 RMSE of the relative attitude estimation

    圖5 相對(duì)角速度估計(jì)的RMSEFig.5 RMSE of the relative angular velocity estimation

    由圖可知,L-NLKF與MSSRCKF誤差收斂后的RMSE基本一致。

    為了定量描述濾波算法在總體上的精度,定義航天器相對(duì)位置、相對(duì)速度、目標(biāo)質(zhì)心位置與相對(duì)角速度的平均RMSE(Time-averaged RMSE,TARMSE),如式(41)所示

    (41)

    表2 濾波算法性能比較(情況1)Table 2 Performance comparison of filtering algorithms(case 1)

    由于航天器間相對(duì)角速度的大小對(duì)濾波估計(jì)的影響較大,考慮目標(biāo)快速旋轉(zhuǎn)情況,將相對(duì)角速度初值擴(kuò)大100倍,即ω0=[1.32,1.32,2.64]T×10-1rad/s,其余初值不變,再進(jìn)行一組蒙特卡洛仿真。統(tǒng)計(jì)結(jié)果如表3所示。

    表3 濾波算法性能比較(情況2)Table 3 Performance comparison of filtering algorithms(case 2)

    對(duì)比表2、表3可知,與航天器間相對(duì)角速度接近于零(無(wú)相對(duì)轉(zhuǎn)動(dòng))相比,目標(biāo)快速旋轉(zhuǎn)時(shí),相對(duì)位置與目標(biāo)形心位置的估計(jì)精度顯著提升。這是由于目標(biāo)快速旋轉(zhuǎn)時(shí),目標(biāo)形心與質(zhì)心的偏差矢量也在追蹤航天器的觀測(cè)中發(fā)生明顯變化,即更易對(duì)相關(guān)狀態(tài)變量進(jìn)行測(cè)量。

    定義濾波器單步平均運(yùn)算時(shí)間為

    (42)

    計(jì)算得到L-NLKF的單步平均運(yùn)算時(shí)間為1.2 ms;MSSRCKF的單步平均運(yùn)算時(shí)間為2.2 ms,即L-NLKF的單步計(jì)算時(shí)間僅是MSSRCKF的55%。結(jié)合表2、3中的數(shù)據(jù)可知,在非合作目標(biāo)相對(duì)追蹤航天器低速、高速旋轉(zhuǎn)的兩種情況下,L-NLKF均能在保證估計(jì)精度的基礎(chǔ)上,大幅縮短計(jì)算時(shí)間。

    5 結(jié) 論

    本文提出了L-NLKF算法,進(jìn)行了理論推導(dǎo),提高了計(jì)算效率。主要工作包括以下幾個(gè)方面:

    1) 根據(jù)星上量測(cè)設(shè)備獲取的目標(biāo)形心的相對(duì)位置與相對(duì)姿態(tài)信息,建立了非合作航天器相對(duì)運(yùn)動(dòng)的軌道-姿態(tài)一體化模型。相較以往的研究,該模型維度較低、更加便于計(jì)算。

    2) 針對(duì)非線性濾波過(guò)程中,狀態(tài)方程包含部分線性項(xiàng)的問(wèn)題,將KF與MSRRCKF進(jìn)行了融合,提出了L-NLKF。該算法能夠?qū)顟B(tài)方程中的線性項(xiàng)與非線性項(xiàng)分別用不同的濾波方法進(jìn)行處理,同時(shí)利用分塊的方式降低矩陣運(yùn)算的維度顯著提升了計(jì)算效率。

    3) 航天器相對(duì)運(yùn)動(dòng)的仿真實(shí)驗(yàn)對(duì)比表明,本文提出的L-NLKF在保證了濾波精度的前提下,計(jì)算時(shí)間明顯少于MSSRCKF。

    猜你喜歡
    協(xié)方差航天器姿態(tài)
    2022 年第二季度航天器發(fā)射統(tǒng)計(jì)
    攀爬的姿態(tài)
    2019 年第二季度航天器發(fā)射統(tǒng)計(jì)
    2018 年第三季度航天器發(fā)射統(tǒng)計(jì)
    全新一代宋的新姿態(tài)
    2018年第二季度航天器發(fā)射統(tǒng)計(jì)
    跑與走的姿態(tài)
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報(bào)器
    一種基于廣義協(xié)方差矩陣的欠定盲辨識(shí)方法
    縱向數(shù)據(jù)分析中使用滑動(dòng)平均Cholesky分解對(duì)回歸均值和協(xié)方差矩陣進(jìn)行同時(shí)半?yún)?shù)建模
    久久久久久九九精品二区国产| 天堂√8在线中文| 国产黄色小视频在线观看| 一本一本综合久久| 99热这里只有精品一区| 精品国内亚洲2022精品成人| 夜夜看夜夜爽夜夜摸| 欧美bdsm另类| 搡老妇女老女人老熟妇| 51国产日韩欧美| 日本三级黄在线观看| 精品少妇黑人巨大在线播放 | 大话2 男鬼变身卡| 日韩av不卡免费在线播放| 我的女老师完整版在线观看| 99久久精品国产国产毛片| 女的被弄到高潮叫床怎么办| 日韩,欧美,国产一区二区三区 | 极品教师在线视频| 亚洲无线观看免费| 日韩av在线免费看完整版不卡| 亚洲av男天堂| 免费观看的影片在线观看| 高清毛片免费看| 伦精品一区二区三区| 日韩av在线免费看完整版不卡| 麻豆国产97在线/欧美| 国产亚洲av嫩草精品影院| 内地一区二区视频在线| 九草在线视频观看| av国产免费在线观看| 青青草视频在线视频观看| 亚洲欧美一区二区三区国产| 赤兔流量卡办理| 我的老师免费观看完整版| 色吧在线观看| 在线免费十八禁| 国产成人精品一,二区| 亚洲电影在线观看av| 亚洲伊人久久精品综合 | 直男gayav资源| 成年版毛片免费区| 婷婷六月久久综合丁香| 一个人免费在线观看电影| 国产乱来视频区| 99久久精品一区二区三区| 插逼视频在线观看| 久久久成人免费电影| 狠狠狠狠99中文字幕| 国产v大片淫在线免费观看| 亚洲无线观看免费| 亚洲精品日韩av片在线观看| 18禁在线播放成人免费| 麻豆精品久久久久久蜜桃| 一级毛片aaaaaa免费看小| 精品熟女少妇av免费看| 久久久久久久久久黄片| 一本一本综合久久| 青春草视频在线免费观看| 免费av不卡在线播放| 男人舔女人下体高潮全视频| 午夜老司机福利剧场| 99久久人妻综合| 亚洲国产精品合色在线| 精品国产露脸久久av麻豆 | 日韩在线高清观看一区二区三区| 欧美日韩国产亚洲二区| .国产精品久久| 只有这里有精品99| 欧美日韩在线观看h| 日韩精品有码人妻一区| 美女黄网站色视频| 国产精品爽爽va在线观看网站| 能在线免费看毛片的网站| 乱系列少妇在线播放| 亚洲av不卡在线观看| 久久久久久久亚洲中文字幕| 欧美+日韩+精品| 国内精品美女久久久久久| 亚洲自拍偷在线| 日日撸夜夜添| 国产黄色小视频在线观看| 国产成人aa在线观看| 国产精品电影一区二区三区| 两个人的视频大全免费| 丰满乱子伦码专区| 青春草国产在线视频| 日韩一区二区视频免费看| 中文亚洲av片在线观看爽| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品亚洲一区二区| 最近视频中文字幕2019在线8| 狂野欧美白嫩少妇大欣赏| 天堂网av新在线| 欧美日韩国产亚洲二区| 亚洲精品乱码久久久久久按摩| 男女边吃奶边做爰视频| 99热这里只有是精品在线观看| 日韩在线高清观看一区二区三区| 精品国产三级普通话版| 2021少妇久久久久久久久久久| 久久久精品大字幕| 亚洲精品自拍成人| 国产黄色小视频在线观看| 亚洲精品456在线播放app| 伊人久久精品亚洲午夜| 3wmmmm亚洲av在线观看| 丝袜美腿在线中文| 欧美成人午夜免费资源| 亚洲av二区三区四区| 九九久久精品国产亚洲av麻豆| 女人十人毛片免费观看3o分钟| 久久精品久久久久久久性| 国产精品99久久久久久久久| 日本wwww免费看| 一区二区三区乱码不卡18| 国产亚洲av片在线观看秒播厂 | 成人毛片a级毛片在线播放| 国产成人精品婷婷| 秋霞在线观看毛片| 国产成人freesex在线| 免费播放大片免费观看视频在线观看 | 一区二区三区四区激情视频| 亚洲国产精品国产精品| 永久网站在线| 听说在线观看完整版免费高清| 麻豆乱淫一区二区| 成年版毛片免费区| 欧美精品国产亚洲| 国产视频内射| 老司机影院毛片| av又黄又爽大尺度在线免费看 | 九九在线视频观看精品| 长腿黑丝高跟| 欧美3d第一页| 久久精品91蜜桃| 国产一区二区亚洲精品在线观看| 极品教师在线视频| 夫妻性生交免费视频一级片| 国产av不卡久久| 韩国av在线不卡| 不卡视频在线观看欧美| 菩萨蛮人人尽说江南好唐韦庄 | 欧美最新免费一区二区三区| 中文精品一卡2卡3卡4更新| 好男人在线观看高清免费视频| 又黄又爽又刺激的免费视频.| 永久免费av网站大全| 久久久久久大精品| 欧美成人精品欧美一级黄| 亚洲av二区三区四区| 丰满人妻一区二区三区视频av| 欧美色视频一区免费| 少妇的逼水好多| 少妇被粗大猛烈的视频| 国产亚洲av嫩草精品影院| 国产伦精品一区二区三区四那| 亚洲自拍偷在线| 18禁在线无遮挡免费观看视频| 久久久a久久爽久久v久久| 午夜福利在线观看免费完整高清在| 亚洲av福利一区| 熟女人妻精品中文字幕| 亚洲性久久影院| 亚洲怡红院男人天堂| 老司机影院成人| 久久亚洲国产成人精品v| 日韩中字成人| 日韩 亚洲 欧美在线| 可以在线观看毛片的网站| 成年女人看的毛片在线观看| 国产不卡一卡二| 变态另类丝袜制服| 一区二区三区免费毛片| 欧美成人免费av一区二区三区| 精品国内亚洲2022精品成人| 免费无遮挡裸体视频| 久久99蜜桃精品久久| 18禁在线播放成人免费| 韩国av在线不卡| 看黄色毛片网站| ponron亚洲| 精品午夜福利在线看| 男人的好看免费观看在线视频| 在线播放国产精品三级| 少妇高潮的动态图| 日本色播在线视频| 成人漫画全彩无遮挡| 精品人妻熟女av久视频| 午夜精品一区二区三区免费看| 欧美区成人在线视频| 亚洲av成人精品一二三区| 一区二区三区乱码不卡18| 精品久久久久久成人av| 免费观看在线日韩| 午夜亚洲福利在线播放| 天堂网av新在线| 亚洲国产精品sss在线观看| 午夜免费激情av| 日日摸夜夜添夜夜添av毛片| 国产精品久久久久久精品电影| 久99久视频精品免费| 久久鲁丝午夜福利片| 国产精品一区二区三区四区久久| 中文字幕制服av| 精品一区二区三区人妻视频| 毛片一级片免费看久久久久| 日日摸夜夜添夜夜添av毛片| 国产私拍福利视频在线观看| 亚洲欧美精品综合久久99| 亚洲中文字幕一区二区三区有码在线看| 最近视频中文字幕2019在线8| 天堂影院成人在线观看| 日产精品乱码卡一卡2卡三| 国模一区二区三区四区视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产黄a三级三级三级人| 久久久久久久国产电影| 91午夜精品亚洲一区二区三区| 91精品伊人久久大香线蕉| 亚洲丝袜综合中文字幕| 亚洲综合色惰| 中文字幕免费在线视频6| 日韩av不卡免费在线播放| 久久午夜福利片| 免费av观看视频| 亚洲国产精品专区欧美| 成人鲁丝片一二三区免费| 日本-黄色视频高清免费观看| 天美传媒精品一区二区| 欧美精品国产亚洲| 91精品国产九色| 国产午夜福利久久久久久| 搡老妇女老女人老熟妇| 成人一区二区视频在线观看| 最近最新中文字幕免费大全7| 国产69精品久久久久777片| av福利片在线观看| 亚州av有码| 国产真实乱freesex| 国产日韩欧美在线精品| 人体艺术视频欧美日本| 99久久精品热视频| 性色avwww在线观看| 久久精品熟女亚洲av麻豆精品 | 天堂影院成人在线观看| 欧美日韩综合久久久久久| 日日撸夜夜添| 精品欧美国产一区二区三| 免费人成在线观看视频色| 国内少妇人妻偷人精品xxx网站| 久久精品综合一区二区三区| 国产麻豆成人av免费视频| 亚洲国产日韩欧美精品在线观看| 自拍偷自拍亚洲精品老妇| 视频中文字幕在线观看| 成人鲁丝片一二三区免费| 国产淫语在线视频| 2021天堂中文幕一二区在线观| 搡老妇女老女人老熟妇| 国产伦精品一区二区三区四那| 最新中文字幕久久久久| 建设人人有责人人尽责人人享有的 | 久久婷婷人人爽人人干人人爱| 非洲黑人性xxxx精品又粗又长| 亚洲成人av在线免费| 免费大片18禁| 久久亚洲国产成人精品v| 国产亚洲91精品色在线| 亚洲美女搞黄在线观看| 村上凉子中文字幕在线| 亚洲成人久久爱视频| 99视频精品全部免费 在线| 狠狠狠狠99中文字幕| 人人妻人人看人人澡| av又黄又爽大尺度在线免费看 | 亚洲真实伦在线观看| 99在线视频只有这里精品首页| 一级黄色大片毛片| 一二三四中文在线观看免费高清| 久久精品国产亚洲av涩爱| 精品久久国产蜜桃| 成人美女网站在线观看视频| 美女国产视频在线观看| 欧美xxxx性猛交bbbb| 夜夜看夜夜爽夜夜摸| 久久韩国三级中文字幕| 亚洲在久久综合| 综合色av麻豆| 国产在线男女| 床上黄色一级片| 久久人妻av系列| 色噜噜av男人的天堂激情| 日本黄色片子视频| 日日啪夜夜撸| 午夜激情福利司机影院| 永久免费av网站大全| 亚洲va在线va天堂va国产| 看十八女毛片水多多多| 嫩草影院新地址| 国产成人午夜福利电影在线观看| 男女边吃奶边做爰视频| 欧美一级a爱片免费观看看| 99热网站在线观看| 国产乱人偷精品视频| 26uuu在线亚洲综合色| 建设人人有责人人尽责人人享有的 | 国产国拍精品亚洲av在线观看| 直男gayav资源| eeuss影院久久| 99热这里只有是精品在线观看| 精品人妻偷拍中文字幕| 男人和女人高潮做爰伦理| 国产精品1区2区在线观看.| 搡老妇女老女人老熟妇| 一个人看视频在线观看www免费| 舔av片在线| 国产激情偷乱视频一区二区| 中文精品一卡2卡3卡4更新| 国产亚洲av片在线观看秒播厂 | 国产精品三级大全| 国产淫语在线视频| 麻豆国产97在线/欧美| 久久久精品大字幕| 欧美精品国产亚洲| 亚洲三级黄色毛片| 精品久久久久久久久av| 日本一本二区三区精品| 最近最新中文字幕免费大全7| 国产精品永久免费网站| 热99在线观看视频| 内地一区二区视频在线| 婷婷六月久久综合丁香| 国产av一区在线观看免费| 少妇丰满av| 免费看日本二区| 2021天堂中文幕一二区在线观| 精品无人区乱码1区二区| 日韩一本色道免费dvd| 国产精品熟女久久久久浪| 最近中文字幕2019免费版| 亚洲高清免费不卡视频| 国产高清不卡午夜福利| 国语对白做爰xxxⅹ性视频网站| 国产黄片美女视频| 99久久无色码亚洲精品果冻| 国产高清不卡午夜福利| 久久久久九九精品影院| 欧美三级亚洲精品| 中文乱码字字幕精品一区二区三区 | 国产三级中文精品| 中文字幕精品亚洲无线码一区| 精品人妻偷拍中文字幕| 精品无人区乱码1区二区| 亚洲av中文av极速乱| 国产美女午夜福利| 2022亚洲国产成人精品| 亚洲一区高清亚洲精品| 日本色播在线视频| 精品99又大又爽又粗少妇毛片| 精华霜和精华液先用哪个| 嘟嘟电影网在线观看| 久久国产乱子免费精品| 极品教师在线视频| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕亚洲精品专区| 亚洲怡红院男人天堂| av线在线观看网站| 免费电影在线观看免费观看| 国产成人精品一,二区| 精品午夜福利在线看| 国产免费又黄又爽又色| 麻豆国产97在线/欧美| 久久这里有精品视频免费| 色哟哟·www| a级一级毛片免费在线观看| 国产男人的电影天堂91| 成人综合一区亚洲| 亚洲av成人av| 国产精品野战在线观看| 99九九线精品视频在线观看视频| 天天一区二区日本电影三级| 一二三四中文在线观看免费高清| 亚洲成人av在线免费| 久久久久九九精品影院| 永久免费av网站大全| 精品一区二区三区视频在线| 亚洲激情五月婷婷啪啪| 干丝袜人妻中文字幕| 亚洲自拍偷在线| 极品教师在线视频| 亚洲伊人久久精品综合 | 少妇被粗大猛烈的视频| 国产乱人偷精品视频| 国产精品,欧美在线| 十八禁国产超污无遮挡网站| 长腿黑丝高跟| 亚洲成av人片在线播放无| av在线天堂中文字幕| 欧美日韩国产亚洲二区| 久久精品影院6| 免费一级毛片在线播放高清视频| 国产亚洲最大av| 国产精品不卡视频一区二区| 两个人视频免费观看高清| 看十八女毛片水多多多| 免费观看的影片在线观看| 午夜福利网站1000一区二区三区| 热99re8久久精品国产| 亚洲精品久久久久久婷婷小说 | 成人毛片a级毛片在线播放| 欧美xxxx性猛交bbbb| 蜜臀久久99精品久久宅男| 欧美3d第一页| 亚洲怡红院男人天堂| 久久久a久久爽久久v久久| 又粗又硬又长又爽又黄的视频| 青春草视频在线免费观看| 成人特级av手机在线观看| 国产午夜精品论理片| 日本爱情动作片www.在线观看| 欧美不卡视频在线免费观看| 亚洲综合精品二区| 男的添女的下面高潮视频| 日韩一区二区三区影片| 久久国产乱子免费精品| 又爽又黄a免费视频| 精品午夜福利在线看| a级一级毛片免费在线观看| 一级二级三级毛片免费看| 国产免费男女视频| av国产久精品久网站免费入址| 日本免费在线观看一区| 日韩 亚洲 欧美在线| av播播在线观看一区| 亚洲欧美精品综合久久99| 久久99热这里只有精品18| 国产亚洲av片在线观看秒播厂 | 日本爱情动作片www.在线观看| 国语自产精品视频在线第100页| 一区二区三区乱码不卡18| 午夜福利网站1000一区二区三区| 国产亚洲一区二区精品| 亚洲欧美成人综合另类久久久 | 乱人视频在线观看| 看十八女毛片水多多多| 日韩亚洲欧美综合| 国产极品天堂在线| 高清毛片免费看| 国产成年人精品一区二区| 麻豆一二三区av精品| 日本免费a在线| 久久久久久久久久久免费av| 1024手机看黄色片| 91aial.com中文字幕在线观看| 午夜福利成人在线免费观看| 97热精品久久久久久| 小说图片视频综合网站| 国产美女午夜福利| www日本黄色视频网| 中文在线观看免费www的网站| 又粗又硬又长又爽又黄的视频| 亚洲成av人片在线播放无| 岛国在线免费视频观看| 午夜福利在线观看免费完整高清在| 最近的中文字幕免费完整| 亚洲精品成人久久久久久| 直男gayav资源| 国产国拍精品亚洲av在线观看| 中文资源天堂在线| videossex国产| 中文天堂在线官网| 亚洲一区高清亚洲精品| 亚洲av.av天堂| 美女被艹到高潮喷水动态| 视频中文字幕在线观看| 中文亚洲av片在线观看爽| 最近2019中文字幕mv第一页| 色噜噜av男人的天堂激情| 久久久久久久亚洲中文字幕| av在线天堂中文字幕| 黄色日韩在线| 亚洲欧美一区二区三区国产| eeuss影院久久| 国产成人a区在线观看| 最近视频中文字幕2019在线8| 国产白丝娇喘喷水9色精品| 久久国内精品自在自线图片| 国产极品天堂在线| 99久久人妻综合| 国产成年人精品一区二区| 欧美高清成人免费视频www| 91精品伊人久久大香线蕉| 成人毛片60女人毛片免费| a级毛色黄片| 久久久成人免费电影| 99热精品在线国产| 亚洲av熟女| 国产一级毛片七仙女欲春2| 非洲黑人性xxxx精品又粗又长| 男女下面进入的视频免费午夜| 日韩视频在线欧美| 免费观看性生交大片5| 九九在线视频观看精品| 国产精品国产三级国产专区5o | 国产伦精品一区二区三区四那| 最近中文字幕高清免费大全6| 欧美日韩一区二区视频在线观看视频在线 | 久久国产乱子免费精品| 真实男女啪啪啪动态图| 美女xxoo啪啪120秒动态图| 国产精品一区二区三区四区免费观看| 日本欧美国产在线视频| 在线观看av片永久免费下载| 一级av片app| 国产成年人精品一区二区| 日日摸夜夜添夜夜添av毛片| 一本久久精品| av播播在线观看一区| 国产大屁股一区二区在线视频| 我要看日韩黄色一级片| 国产av码专区亚洲av| a级毛色黄片| av免费在线看不卡| 国产伦一二天堂av在线观看| 最近手机中文字幕大全| 只有这里有精品99| 国产综合懂色| 少妇裸体淫交视频免费看高清| 欧美日韩一区二区视频在线观看视频在线 | 国产不卡一卡二| 精品国内亚洲2022精品成人| 欧美激情在线99| 91久久精品电影网| 男人狂女人下面高潮的视频| 国产精品一区二区三区四区久久| 国产精品乱码一区二三区的特点| 亚洲欧美日韩高清专用| 看非洲黑人一级黄片| 精品久久久噜噜| 校园人妻丝袜中文字幕| 日本黄大片高清| 婷婷色av中文字幕| 天堂√8在线中文| 国产精品国产高清国产av| 人妻少妇偷人精品九色| 亚洲怡红院男人天堂| 日韩欧美在线乱码| 97热精品久久久久久| 久久99精品国语久久久| 丝袜喷水一区| 最近中文字幕高清免费大全6| 午夜福利网站1000一区二区三区| 亚洲久久久久久中文字幕| a级毛片免费高清观看在线播放| 免费无遮挡裸体视频| 午夜激情欧美在线| 亚洲综合色惰| 欧美激情国产日韩精品一区| 久久久久网色| 日韩精品有码人妻一区| 日韩欧美在线乱码| 免费观看的影片在线观看| 亚洲成人中文字幕在线播放| 免费人成在线观看视频色| 91精品伊人久久大香线蕉| 婷婷色麻豆天堂久久 | 欧美性猛交╳xxx乱大交人| 久久精品91蜜桃| 女的被弄到高潮叫床怎么办| 久久精品综合一区二区三区| 最近中文字幕2019免费版| 免费看av在线观看网站| 五月玫瑰六月丁香| 亚洲国产精品合色在线| 最近的中文字幕免费完整| 高清日韩中文字幕在线| 特级一级黄色大片| 简卡轻食公司| 99热网站在线观看| 内射极品少妇av片p| 午夜爱爱视频在线播放| 免费av毛片视频| 搞女人的毛片| 国产成人aa在线观看| 国产私拍福利视频在线观看| 免费看光身美女| 成人av在线播放网站| 亚洲怡红院男人天堂| 国产成人一区二区在线| 久久久久网色| 亚洲伊人久久精品综合 | 成人特级av手机在线观看| 中文亚洲av片在线观看爽| 亚洲欧美精品综合久久99| 九九热线精品视视频播放| 亚洲国产高清在线一区二区三| 亚洲欧美精品综合久久99| 亚洲精品国产成人久久av| 中文亚洲av片在线观看爽| 老司机福利观看| 国产高清国产精品国产三级 | 男人狂女人下面高潮的视频| 男女边吃奶边做爰视频| 国产精品久久久久久久电影| 亚洲av电影不卡..在线观看| 舔av片在线| 人妻夜夜爽99麻豆av| 国产乱人偷精品视频| 最近中文字幕高清免费大全6| 日韩av在线免费看完整版不卡| 日本一本二区三区精品| 91精品伊人久久大香线蕉| 黄片无遮挡物在线观看| 国产伦精品一区二区三区四那|