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

    一種地理系下基于偽線(xiàn)性模型的捷聯(lián)慣導(dǎo)算法

    2022-10-15 13:24:40楊小康嚴(yán)恭敏李思錦李四海劉明雍
    宇航學(xué)報(bào) 2022年9期
    關(guān)鍵詞:捷聯(lián)慣導(dǎo)角速度

    楊小康,嚴(yán)恭敏,李思錦,李四海,劉明雍

    (1. 西北工業(yè)大學(xué)自動(dòng)化學(xué)院,西安 710129; 2. 陜西省飛行控制與仿真技術(shù)重點(diǎn)實(shí)驗(yàn)室,西安 710129;3. 西安微電子技術(shù)研究所,西安 710065; 4. 西北工業(yè)大學(xué)航海學(xué)院,西安 710072)

    0 引 言

    捷聯(lián)式慣導(dǎo)系統(tǒng)與平臺(tái)式不同,需要利用陀螺儀數(shù)據(jù)完成姿態(tài)解算,構(gòu)建虛擬的數(shù)學(xué)平臺(tái),再利用加速度計(jì)數(shù)據(jù)實(shí)現(xiàn)速度和位置解算。因此捷聯(lián)慣導(dǎo)系統(tǒng)的導(dǎo)航解算算法相較于平臺(tái)系統(tǒng)更加復(fù)雜,也更加重要,捷聯(lián)慣導(dǎo)算法研究也一直是捷聯(lián)慣導(dǎo)系統(tǒng)的研究重點(diǎn)。目前受制于慣導(dǎo)器件精度,慣導(dǎo)算法研究的主要方向是組合導(dǎo)航、冗余慣導(dǎo)系統(tǒng)、長(zhǎng)航時(shí)誤差抑制等。但是原子陀螺慣導(dǎo)系統(tǒng)、空天飛機(jī)、高超音速飛行器等新的應(yīng)用場(chǎng)景對(duì)捷聯(lián)慣導(dǎo)解算算法提出了新的挑戰(zhàn),需要針對(duì)捷聯(lián)慣導(dǎo)解算的算法結(jié)構(gòu)或原理進(jìn)行改進(jìn)。在傳統(tǒng)捷聯(lián)慣導(dǎo)算法發(fā)展過(guò)程中,隨著對(duì)算法誤差機(jī)理研究的深入,算法精度不斷提高。文獻(xiàn)[6]利用等效旋轉(zhuǎn)矢量完成姿態(tài)解算,建立了捷聯(lián)慣導(dǎo)算法的基礎(chǔ)。在此基礎(chǔ)上,文獻(xiàn)[7-8]又進(jìn)一步提出了多子樣捷聯(lián)慣導(dǎo)算法,通過(guò)補(bǔ)償圓錐誤差和劃槳誤差提升算法精度。

    由于構(gòu)成捷聯(lián)慣導(dǎo)系統(tǒng)的陀螺儀和加速度計(jì)精度提升,加之大機(jī)動(dòng)高動(dòng)態(tài)條件下導(dǎo)航需求的出現(xiàn),近幾年的捷聯(lián)慣導(dǎo)算法研究重點(diǎn)是提升捷聯(lián)慣導(dǎo)算法精度。文獻(xiàn)[9]通過(guò)建立姿態(tài)更新不可交換誤差的高階誤差模型提出了基于角增量四次叉乘補(bǔ)償項(xiàng)的高精度捷聯(lián)慣導(dǎo)算法。文獻(xiàn)[10]提出使用數(shù)值仿真算法求解不可交換誤差高階補(bǔ)償系數(shù),實(shí)現(xiàn)此類(lèi)算法中補(bǔ)償系數(shù)的快速計(jì)算。文獻(xiàn)[11-12]則是基于函數(shù)迭代方法,構(gòu)建角速度和比力多項(xiàng)式并完成高精度導(dǎo)航解算。

    無(wú)論是傳統(tǒng)的多子樣捷聯(lián)慣導(dǎo)算法還是新提出的幾種高精度捷聯(lián)慣導(dǎo)算法,都是利用姿態(tài)陣、四元數(shù)、羅德里格參數(shù)等進(jìn)行姿態(tài)解算,再根據(jù)對(duì)應(yīng)的坐標(biāo)轉(zhuǎn)換計(jì)算規(guī)則,將載體坐標(biāo)系下的加速度計(jì)輸出轉(zhuǎn)換到導(dǎo)航坐標(biāo)系下。因?yàn)橐陨媳硎咀藨B(tài)的變量并不是在線(xiàn)性空間中定義的,所以速度更新方程是非線(xiàn)性的,不便于直接利用捷聯(lián)慣導(dǎo)系統(tǒng)模型進(jìn)行系統(tǒng)分析、濾波等操作。文獻(xiàn)[13]介紹了偽線(xiàn)性系統(tǒng)相關(guān)算法基礎(chǔ),偽線(xiàn)性系統(tǒng)是指形似線(xiàn)性系統(tǒng)實(shí)則非線(xiàn)性系統(tǒng)的系統(tǒng)。文獻(xiàn)[14-15]介紹了偽線(xiàn)性系統(tǒng)方法在衛(wèi)星姿態(tài)控制、飛行器制導(dǎo)與控制中的應(yīng)用。如果構(gòu)建出捷聯(lián)慣導(dǎo)系統(tǒng)的偽線(xiàn)性系統(tǒng)模型,就可以直接使用線(xiàn)性系統(tǒng)理論完成捷聯(lián)慣導(dǎo)更新,直接使用捷聯(lián)慣導(dǎo)系統(tǒng)模型構(gòu)建線(xiàn)性Kalman濾波器實(shí)現(xiàn)初始對(duì)準(zhǔn)、組合導(dǎo)航等,無(wú)需考慮模型的非線(xiàn)性問(wèn)題。按照此研究路線(xiàn),文獻(xiàn)[16]提出了基于時(shí)間序列的慣性系下偽線(xiàn)性捷聯(lián)慣導(dǎo)算法,建立了完整的偽線(xiàn)性捷聯(lián)慣導(dǎo)算法,但是此算法建立在慣性系下,不能滿(mǎn)足通常的導(dǎo)航系下定位需求,因此也不能直接將此模型應(yīng)用于初始對(duì)準(zhǔn)和組合導(dǎo)航中。

    本文在偽線(xiàn)性模型相關(guān)理論和文獻(xiàn)[16]中偽線(xiàn)性捷聯(lián)慣導(dǎo)算法基礎(chǔ)上,提出地理系下偽線(xiàn)性高精度捷聯(lián)慣導(dǎo)更新算法。與傳統(tǒng)的捷聯(lián)慣導(dǎo)算法相比,偽線(xiàn)性的算法不再將一個(gè)完整的慣導(dǎo)更新過(guò)程分為姿態(tài)更新、速度更新和位置更新三部分,而是使用新構(gòu)建的導(dǎo)航向量統(tǒng)一表示姿態(tài)、速度和位置信息,并且根據(jù)導(dǎo)航向量微分方程設(shè)計(jì)離散的導(dǎo)航向量更新算法。應(yīng)用偽線(xiàn)性捷聯(lián)慣導(dǎo)算法的慣導(dǎo)系統(tǒng)結(jié)構(gòu)如圖1所示,原有的三個(gè)更新回路被導(dǎo)航向量及其線(xiàn)性更新方程代替,更新后的導(dǎo)航向量再轉(zhuǎn)換成姿態(tài)、速度和位置輸出導(dǎo)航結(jié)果。偽線(xiàn)性捷聯(lián)慣導(dǎo)算法因采用了基于泰勒級(jí)數(shù)展開(kāi)的數(shù)值算法,因此在大機(jī)動(dòng)或圓錐運(yùn)動(dòng)等能激勵(lì)出較大不可交換誤差的情況下,理論算法精度遠(yuǎn)高于傳統(tǒng)算法。除此之外,捷聯(lián)慣導(dǎo)系統(tǒng)偽線(xiàn)性模型的另一優(yōu)勢(shì)是可以使用線(xiàn)性的分析與計(jì)算方法。

    圖1 偽線(xiàn)性捷聯(lián)慣導(dǎo)系統(tǒng)結(jié)構(gòu)Fig.1 Structure of pseudo-linear SINS

    最后,為了驗(yàn)證基于偽線(xiàn)性模型的捷聯(lián)慣導(dǎo)算法與傳統(tǒng)算法的差異,用兩種算法分別完成了大機(jī)動(dòng)條件下導(dǎo)航解算和旋轉(zhuǎn)彈飛行仿真導(dǎo)航解算,并比較兩者的導(dǎo)航誤差。

    1 偽線(xiàn)性SINS更新公式

    和傳統(tǒng)的捷聯(lián)慣導(dǎo)算法類(lèi)似,偽線(xiàn)性算法模型也將按照姿態(tài)、速度和位置三個(gè)更新過(guò)程進(jìn)行分析。在捷聯(lián)慣導(dǎo)更新算法的推導(dǎo)過(guò)程中,表示慣性坐標(biāo)系,表示載體坐標(biāo)系,表示導(dǎo)航坐標(biāo)系,表示地球坐標(biāo)系。

    1.1 姿態(tài)更新

    在捷聯(lián)慣導(dǎo)算法中姿態(tài)四元數(shù)微分方程為

    (1)

    (2)

    (3)

    (4)

    表示式(1)模型中的系統(tǒng)矩陣,即

    (5)

    則導(dǎo)航解算的姿態(tài)更新可以表示為

    (6)

    1.2 速度更新

    比力方程是捷聯(lián)慣導(dǎo)解算的核心之一,求解比力方程即可完成速度更新。捷聯(lián)慣導(dǎo)系統(tǒng)導(dǎo)航系下比力方程為

    (7)

    在式(7)中改用四元數(shù)將系下的比力轉(zhuǎn)換到系下,得到:

    (8)

    由式(8)右邊第一項(xiàng)可以看出來(lái),使用四元數(shù)進(jìn)行導(dǎo)航解算時(shí),速度更新方程是非線(xiàn)性的。

    為了得到線(xiàn)性的速度更新方程,對(duì)速度向量進(jìn)行變換,定義新的速度向量——速度四元數(shù)

    (9)

    對(duì)式(9)兩邊求導(dǎo)

    (10)

    將式(8)寫(xiě)成四元數(shù)形式

    (11)

    將式(2)和式(11)代入式(10),得

    (12)

    (13)

    式(12)可整理成

    (14)

    式中:

    (15)

    顯然,式(14)中的速度微分方程是關(guān)于姿態(tài)四元數(shù)和速度四元數(shù)的線(xiàn)性方程。

    1.3 位置更新

    在地理系下通常計(jì)算載體所在緯度、經(jīng)度和高度作為位置更新結(jié)果,因此定義位置向量

    (16)

    式中:為緯度,為經(jīng)度,表示高度。其微分方程為

    (17)

    式中:

    (18)

    根據(jù)式(18)不能直接得出位置四元數(shù)與速度四元數(shù)之間的線(xiàn)性關(guān)系,所以需要研究新的位置四元數(shù)定義。

    由式(9)得

    (19)

    (20)

    (21)

    式(20)寫(xiě)成矩陣形式為

    (22)

    式中:

    (23)

    通過(guò)構(gòu)建位置四元數(shù)向量,位置更新方程被轉(zhuǎn)換成關(guān)于姿態(tài)四元數(shù)、速度四元數(shù)和位置四元數(shù)的偽線(xiàn)性方程。

    2 導(dǎo)航向量及其更新方案

    四元數(shù)微分方程為線(xiàn)性方程,比力方程雖然是非線(xiàn)性方程,但是經(jīng)過(guò)變換得到的速度四元數(shù)微分方程是線(xiàn)性的,同理通過(guò)定義位置四元數(shù)得到了偽線(xiàn)性的位置更新方程,因此可以定義包含四元數(shù)、速度四元數(shù)和位置四元數(shù)的導(dǎo)航向量

    (24)

    (25)

    (26)

    然后由式(6)、式(14)和式(22)可得導(dǎo)航向量微分方程

    (27)

    ()=++…+-1-1

    (28)

    式(27)是線(xiàn)性時(shí)變系統(tǒng),因此在任意更新周期內(nèi),狀態(tài)轉(zhuǎn)移矩陣滿(mǎn)足:

    (29)

    為了求解式(29),將(,0)表示成關(guān)于時(shí)間的次矩陣多項(xiàng)式,即

    (,0)=+++…+

    (30)

    對(duì)式(30)求導(dǎo)可得

    (31)

    將式(28)、(30)和(31)代入式(29),得

    (32)

    整理后得:

    (33)

    式中:是(,0)的展開(kāi)次數(shù),決定其計(jì)算精度,根據(jù)式(30)構(gòu)建狀態(tài)轉(zhuǎn)移矩陣的精度可以達(dá)到(+1),因此在設(shè)計(jì)算法時(shí)可以根據(jù)更新周期、導(dǎo)航算法精度等指標(biāo)選擇合適的展開(kāi)次數(shù)完成計(jì)算。

    根據(jù)式(28)、(30)和(33)可以計(jì)算出狀態(tài)轉(zhuǎn)移矩陣(, 0),然后完成導(dǎo)航向量的更新

    ()=(, 0)(0)

    (34)

    最后將更新后的導(dǎo)航向量轉(zhuǎn)換成地理系下的姿態(tài)、速度和位置。

    3 算法精度評(píng)估

    為了檢驗(yàn)偽線(xiàn)性模型的導(dǎo)航解算效果,在短時(shí)大角速度大加速度條件和設(shè)計(jì)的旋轉(zhuǎn)彈飛行軌跡上進(jìn)行純慣性導(dǎo)航解算,比較偽線(xiàn)性慣導(dǎo)算法與傳統(tǒng)慣導(dǎo)算法的差異,計(jì)算時(shí)統(tǒng)一采用4子樣算法結(jié)構(gòu)。為便于表述評(píng)估結(jié)果,后文中使用“NLNA(Nonlinear navigation algorithm,非線(xiàn)性導(dǎo)航算法)”表示傳統(tǒng)的捷聯(lián)慣導(dǎo)算法,用“PLNA(Pseudo-linear navigation algorithm,偽線(xiàn)性導(dǎo)航算法)”表示新提出的捷聯(lián)慣導(dǎo)算法。

    偽線(xiàn)性算法的4子樣解算結(jié)構(gòu)如下:

    ① 根據(jù)4子樣角增量和速度增量Δ,…,Δ和Δ,…,Δ,計(jì)算角速度和比力多項(xiàng)式系數(shù),得到角速度和比力函數(shù):

    (35)

    (36)

    ② 根據(jù)式(6)、式(14)和式(22)得到不同次系統(tǒng)矩陣的各個(gè)子矩陣。0次系數(shù)矩陣中各項(xiàng)如下:

    (37)

    1次及以上系數(shù)矩陣各項(xiàng)為:

    (38)

    ③ 根據(jù)式(27)得到不同次系數(shù)矩陣

    (39)

    并根據(jù)式(28)得到導(dǎo)航向量系統(tǒng)矩陣的多項(xiàng)式形式:

    =+++

    (40)

    ④ 將代入式(33)得到狀態(tài)轉(zhuǎn)移矩陣不同次項(xiàng),按照式(30)計(jì)算狀態(tài)轉(zhuǎn)移矩陣

    (,0)=+++

    (41)

    ⑤ 最后根據(jù)式(34)完成導(dǎo)航向量更新,并根據(jù)式(9)和式(11)中的定義計(jì)算出地理系下的速度和位置信息。

    3.1 大機(jī)動(dòng)條件慣導(dǎo)解算仿真

    大機(jī)動(dòng)運(yùn)動(dòng)條件對(duì)捷聯(lián)慣導(dǎo)算法有更高的要求,對(duì)于傳統(tǒng)捷聯(lián)慣導(dǎo)算法,在大機(jī)動(dòng)條件下會(huì)激勵(lì)出較大的算法誤差。捷聯(lián)慣導(dǎo)算法研究中一般都采用多項(xiàng)式形式的角速度和加速度模擬大機(jī)動(dòng)條件,并在此條件下檢驗(yàn)算法精度。按照文獻(xiàn)[18-19]中仿真方法,用角速度多項(xiàng)式模擬大機(jī)動(dòng)角運(yùn)動(dòng),設(shè)置大機(jī)動(dòng)條件角速度和加速度多項(xiàng)式如式(42)~(43)所示,對(duì)應(yīng)的角速度與加速度變化如圖2所示。仿真時(shí)按1 kHz采樣得到慣導(dǎo)仿真數(shù)據(jù),同時(shí)以100 kHz采樣完成2子樣導(dǎo)航解算,得到姿態(tài)和速度的參考結(jié)果。兩種捷聯(lián)慣導(dǎo)算法的姿態(tài)誤差和速度誤差如圖3~4所示,導(dǎo)航解算誤差如表1所列,因?yàn)椴煌惴ǖ淖藨B(tài)和速度誤差相差較大,因此不必全部列出,選取導(dǎo)航結(jié)束時(shí)刻的三個(gè)軸上誤差的模值作為姿態(tài)誤差和速度誤差。

    圖2 大機(jī)動(dòng)條件角速度與加速度Fig.2 Angular rate and acceleration in high-dynamic condition

    圖3 大機(jī)動(dòng)條件姿態(tài)誤差Fig.3 Attitude errors in high-dynamic condition

    圖4 大機(jī)動(dòng)條件速度誤差Fig.4 Velocity errors in high-dynamic condition

    (42)

    (43)

    從圖3~4和表1中的姿態(tài)誤差與速度誤差對(duì)比結(jié)果可以明顯看出偽線(xiàn)性算法的導(dǎo)航誤差更小。

    表1 導(dǎo)航解算仿真誤差Table 1 Errors of navigation solution simulation

    這是因?yàn)榛趥尉€(xiàn)性模型的算法在進(jìn)行姿態(tài)和速度更新時(shí),通過(guò)求解姿態(tài)和速度微分方程的高階導(dǎo)數(shù),計(jì)算出高精度的導(dǎo)航結(jié)果。在角速度和加速度為關(guān)于時(shí)間的一般多項(xiàng)式的條件下,此算法可以通過(guò)增加泰勒展開(kāi)次數(shù)減小算法誤差,這一特點(diǎn)也使得偽線(xiàn)性捷聯(lián)慣導(dǎo)算法更加適合應(yīng)用于大機(jī)動(dòng)條件。

    3.2 旋轉(zhuǎn)彈純慣性導(dǎo)航仿真

    旋轉(zhuǎn)彈飛行過(guò)程中含有較長(zhǎng)的大角速度機(jī)動(dòng),容易激勵(lì)出捷聯(lián)慣導(dǎo)算法誤差,因此和傳統(tǒng)算法相比,偽線(xiàn)性模型捷聯(lián)慣導(dǎo)算法更加適用于旋轉(zhuǎn)彈捷聯(lián)慣導(dǎo)解算。為了驗(yàn)證捷聯(lián)慣導(dǎo)算法性能,設(shè)計(jì)旋轉(zhuǎn)彈飛行過(guò)程仿真軌跡如圖5所示,飛行中彈體角速度從20π rad/s逐漸減速至2π rad/s,慣導(dǎo)數(shù)據(jù)更新頻率為1 kHz。

    圖5 旋轉(zhuǎn)彈仿真軌跡Fig.5 Simulation trajectory of spinning projectile

    使用4子樣算法捷聯(lián)慣導(dǎo)解算仿真并統(tǒng)計(jì)姿態(tài)誤差、速度誤差和位置誤差的均方根以分析導(dǎo)航誤差,評(píng)估兩種算法的精度。導(dǎo)航誤差對(duì)比如圖6~7所示,全程的均方根誤差如表2所列。

    從圖6~7和表2可以得出以下結(jié)論:

    表2 旋轉(zhuǎn)彈仿真算法誤差Table 2 Algorithm errors in spinning projectile simulation

    圖6 旋轉(zhuǎn)彈仿真姿態(tài)誤差Fig.6 Attitude errors of spinning projectile simulation

    (1)基于非線(xiàn)性模型的傳統(tǒng)捷聯(lián)慣導(dǎo)算法在大橫滾角速度條件下算法誤差遠(yuǎn)大于偽線(xiàn)性算法;

    (2)從圖7的位置誤差變化圖上可以看出,傳統(tǒng)算法的位置誤差在導(dǎo)航過(guò)程中一直發(fā)散,而偽線(xiàn)性導(dǎo)航算法的位置誤差則在后半程基本不變,這說(shuō)明偽線(xiàn)性算法在旋轉(zhuǎn)彈導(dǎo)航中有明顯的優(yōu)勢(shì);

    圖7 旋轉(zhuǎn)彈仿真速度與位置誤差Fig.7 Velocity and position errors of spinning projectile simulation

    (3)仿真結(jié)果表明,對(duì)于旋轉(zhuǎn)彈,傳統(tǒng)捷聯(lián)慣導(dǎo)算法的算法誤差已不能被忽略,因此捷聯(lián)慣導(dǎo)算法誤差模型不再準(zhǔn)確,如果要進(jìn)行組合導(dǎo)航,應(yīng)當(dāng)以偽線(xiàn)性捷聯(lián)慣導(dǎo)算法為基礎(chǔ),建立組合導(dǎo)航系統(tǒng)線(xiàn)性狀態(tài)空間模型。

    4 結(jié) 論

    本文在已有的高精度捷聯(lián)慣導(dǎo)算法和偽線(xiàn)性系統(tǒng)理論的基礎(chǔ)上,提出了完整的地理系下基于偽線(xiàn)性模型的捷聯(lián)慣導(dǎo)算法。該算法可以求解線(xiàn)性系統(tǒng)方程實(shí)現(xiàn)捷聯(lián)慣導(dǎo)解算更新,且在大機(jī)動(dòng)條件下,例如有較大橫滾角速度的旋轉(zhuǎn)彈的導(dǎo)航解算中,有較高的算法精度。與傳統(tǒng)的捷聯(lián)慣導(dǎo)算法相比,基于偽線(xiàn)性模型的算法具有精度高、結(jié)構(gòu)簡(jiǎn)單、線(xiàn)性化等優(yōu)點(diǎn),對(duì)于捷聯(lián)慣導(dǎo)算法研究與工程應(yīng)用有一定的參考價(jià)值。

    猜你喜歡
    捷聯(lián)慣導(dǎo)角速度
    自適應(yīng)模糊多環(huán)控制在慣導(dǎo)平臺(tái)穩(wěn)定回路中的應(yīng)用
    無(wú)人機(jī)室內(nèi)視覺(jué)/慣導(dǎo)組合導(dǎo)航方法
    彈道導(dǎo)彈的捷聯(lián)慣性/天文組合導(dǎo)航方法
    圓周運(yùn)動(dòng)角速度測(cè)量方法賞析
    基于Bagging模型的慣導(dǎo)系統(tǒng)誤差抑制方法
    捷聯(lián)慣性/天文/雷達(dá)高度表組合導(dǎo)航
    半捷聯(lián)雷達(dá)導(dǎo)引頭視線(xiàn)角速度提取
    基于多線(xiàn)程的慣導(dǎo)邏輯仿真器設(shè)計(jì)
    基于構(gòu)架點(diǎn)頭角速度的軌道垂向長(zhǎng)波不平順在線(xiàn)檢測(cè)
    一種捷聯(lián)式圖像導(dǎo)引頭的解耦算法
    免费大片黄手机在线观看| 国产精品99久久久久久久久| 久久热精品热| 国产一区二区在线观看av| 久久久久人妻精品一区果冻| 久久久久国产精品人妻一区二区| 少妇被粗大的猛进出69影院 | 天堂中文最新版在线下载| 国产黄频视频在线观看| 久久午夜综合久久蜜桃| 免费看av在线观看网站| 色5月婷婷丁香| 男女啪啪激烈高潮av片| 国国产精品蜜臀av免费| 婷婷色综合www| 午夜精品国产一区二区电影| 国产成人精品福利久久| 日韩伦理黄色片| 久久久久人妻精品一区果冻| 蜜桃国产av成人99| 国产伦理片在线播放av一区| 国产精品无大码| 制服人妻中文乱码| 黄色视频在线播放观看不卡| 久久久欧美国产精品| 中文天堂在线官网| 精品一区在线观看国产| 午夜91福利影院| 亚洲欧美中文字幕日韩二区| 在线免费观看不下载黄p国产| 久久人人爽人人片av| 国产精品久久久久久精品古装| 国产黄色视频一区二区在线观看| 人妻系列 视频| 嘟嘟电影网在线观看| 中文字幕制服av| 乱人伦中国视频| 在线观看人妻少妇| 精品视频人人做人人爽| 99国产精品免费福利视频| 久久久久国产网址| 你懂的网址亚洲精品在线观看| 亚洲av中文av极速乱| 亚洲欧洲国产日韩| 日本av手机在线免费观看| 青春草国产在线视频| 蜜桃在线观看..| 亚洲精品亚洲一区二区| 国产精品一区二区在线不卡| 婷婷色av中文字幕| 亚洲在久久综合| av在线观看视频网站免费| 狂野欧美激情性xxxx在线观看| 丝袜脚勾引网站| 黄色欧美视频在线观看| 亚洲一区二区三区欧美精品| 人人妻人人添人人爽欧美一区卜| 久久久久久人妻| 亚洲中文av在线| 3wmmmm亚洲av在线观看| 国产成人a∨麻豆精品| 久久久久久人妻| 亚洲av二区三区四区| 亚洲欧洲国产日韩| 日韩欧美精品免费久久| h视频一区二区三区| 亚洲av日韩在线播放| xxx大片免费视频| 老司机亚洲免费影院| 天天影视国产精品| 婷婷色麻豆天堂久久| 伦理电影免费视频| 久热久热在线精品观看| 久久女婷五月综合色啪小说| 亚洲一级一片aⅴ在线观看| 2021少妇久久久久久久久久久| 亚洲少妇的诱惑av| 久久人妻熟女aⅴ| 啦啦啦在线观看免费高清www| 免费日韩欧美在线观看| 亚洲三级黄色毛片| 在线观看免费日韩欧美大片 | 一级毛片我不卡| 欧美激情 高清一区二区三区| 免费人成在线观看视频色| 久久狼人影院| 大又大粗又爽又黄少妇毛片口| 久久这里有精品视频免费| 欧美日本中文国产一区发布| 日本黄大片高清| 2022亚洲国产成人精品| 最新中文字幕久久久久| 日韩中文字幕视频在线看片| 亚洲情色 制服丝袜| 美女cb高潮喷水在线观看| 99热6这里只有精品| 午夜精品国产一区二区电影| 精品久久久噜噜| 欧美激情极品国产一区二区三区 | 又大又黄又爽视频免费| 亚洲四区av| 日韩电影二区| 国产精品无大码| 啦啦啦中文免费视频观看日本| 视频区图区小说| 能在线免费看毛片的网站| 午夜影院在线不卡| av在线观看视频网站免费| 各种免费的搞黄视频| 七月丁香在线播放| 国产精品久久久久久av不卡| 亚洲国产欧美日韩在线播放| 少妇高潮的动态图| 亚洲精品视频女| 国产精品国产三级国产av玫瑰| 精品少妇久久久久久888优播| 最近2019中文字幕mv第一页| 97在线人人人人妻| 在线精品无人区一区二区三| 美女脱内裤让男人舔精品视频| 国产 一区精品| 亚洲内射少妇av| av在线播放精品| 一边摸一边做爽爽视频免费| 精品人妻一区二区三区麻豆| 青春草视频在线免费观看| 久久鲁丝午夜福利片| 美女大奶头黄色视频| 九九久久精品国产亚洲av麻豆| 国产高清国产精品国产三级| 2021少妇久久久久久久久久久| 18禁动态无遮挡网站| 精品视频人人做人人爽| 国产成人精品在线电影| 日产精品乱码卡一卡2卡三| 久久久久视频综合| 一级毛片我不卡| 伊人久久国产一区二区| 视频在线观看一区二区三区| 亚洲色图综合在线观看| 一级毛片aaaaaa免费看小| 性高湖久久久久久久久免费观看| 精品酒店卫生间| 免费观看性生交大片5| 亚洲精品久久午夜乱码| 国产在视频线精品| 久热久热在线精品观看| 久久精品熟女亚洲av麻豆精品| 国产成人freesex在线| 边亲边吃奶的免费视频| 国产伦精品一区二区三区视频9| 亚洲av电影在线观看一区二区三区| 能在线免费看毛片的网站| 亚洲国产欧美日韩在线播放| 精品一区二区三卡| 精品国产一区二区三区久久久樱花| 欧美bdsm另类| 美女内射精品一级片tv| 国产成人午夜福利电影在线观看| 成人国产av品久久久| 日本黄色日本黄色录像| 中文精品一卡2卡3卡4更新| 在线免费观看不下载黄p国产| 人妻制服诱惑在线中文字幕| 大码成人一级视频| 国产欧美日韩一区二区三区在线 | 一边摸一边做爽爽视频免费| 99视频精品全部免费 在线| 国产精品免费大片| 免费看av在线观看网站| 国产av码专区亚洲av| 夫妻午夜视频| 久久精品国产鲁丝片午夜精品| 国产成人精品无人区| 青春草亚洲视频在线观看| 在线亚洲精品国产二区图片欧美 | 天堂俺去俺来也www色官网| 伦理电影大哥的女人| 亚洲av福利一区| 国产69精品久久久久777片| 亚洲av二区三区四区| 精品久久久久久久久av| 黑人巨大精品欧美一区二区蜜桃 | 精品一区在线观看国产| 一个人免费看片子| 精品久久久久久久久亚洲| 2021少妇久久久久久久久久久| a级片在线免费高清观看视频| 日韩伦理黄色片| 国产老妇伦熟女老妇高清| 人人妻人人爽人人添夜夜欢视频| 国产国拍精品亚洲av在线观看| 日韩一区二区三区影片| 精品熟女少妇av免费看| 亚洲精品av麻豆狂野| 欧美精品高潮呻吟av久久| 日产精品乱码卡一卡2卡三| 国产免费一区二区三区四区乱码| 午夜免费男女啪啪视频观看| 日韩大片免费观看网站| 亚洲国产欧美日韩在线播放| 视频在线观看一区二区三区| 99久国产av精品国产电影| 国产伦理片在线播放av一区| 自拍欧美九色日韩亚洲蝌蚪91| 日本黄色片子视频| 91久久精品电影网| 午夜91福利影院| kizo精华| 国产男人的电影天堂91| 伊人久久国产一区二区| 久久99蜜桃精品久久| 中文字幕精品免费在线观看视频 | 91国产中文字幕| 午夜老司机福利剧场| 少妇的逼水好多| 亚洲欧美成人综合另类久久久| 欧美精品亚洲一区二区| 欧美三级亚洲精品| 99久久综合免费| 欧美变态另类bdsm刘玥| 亚洲精品色激情综合| 91精品国产九色| 在线看a的网站| 久久97久久精品| 日韩欧美一区视频在线观看| 精品国产一区二区久久| 丰满饥渴人妻一区二区三| av.在线天堂| 在线观看美女被高潮喷水网站| 丝袜喷水一区| 日本wwww免费看| 最近的中文字幕免费完整| 久久国产精品大桥未久av| 成人18禁高潮啪啪吃奶动态图 | av福利片在线| 99视频精品全部免费 在线| 久久久久国产网址| 国产黄频视频在线观看| 国产高清国产精品国产三级| 亚洲欧美日韩另类电影网站| 午夜福利,免费看| 欧美人与性动交α欧美精品济南到 | 亚洲精品,欧美精品| 亚洲情色 制服丝袜| 妹子高潮喷水视频| 99国产精品免费福利视频| 久久精品国产自在天天线| 国产成人a∨麻豆精品| 十八禁高潮呻吟视频| 中文字幕人妻丝袜制服| 国产日韩一区二区三区精品不卡 | 亚洲av成人精品一二三区| 午夜激情福利司机影院| 三级国产精品欧美在线观看| 亚洲色图综合在线观看| 伦精品一区二区三区| 18禁裸乳无遮挡动漫免费视频| 晚上一个人看的免费电影| 国产精品 国内视频| 久久午夜福利片| 精品少妇黑人巨大在线播放| 97在线视频观看| 十八禁高潮呻吟视频| 亚洲色图 男人天堂 中文字幕 | 青春草国产在线视频| 亚洲国产精品一区二区三区在线| 亚洲综合精品二区| 黄色一级大片看看| 视频中文字幕在线观看| 蜜桃在线观看..| 777米奇影视久久| 高清毛片免费看| 中文欧美无线码| 国产片特级美女逼逼视频| 免费大片黄手机在线观看| 美女主播在线视频| 2018国产大陆天天弄谢| 免费观看在线日韩| 国产欧美另类精品又又久久亚洲欧美| 人妻少妇偷人精品九色| 老女人水多毛片| 你懂的网址亚洲精品在线观看| √禁漫天堂资源中文www| 亚洲精品456在线播放app| 午夜影院在线不卡| 黄片无遮挡物在线观看| 欧美性感艳星| 久久精品国产鲁丝片午夜精品| 寂寞人妻少妇视频99o| 亚洲欧洲日产国产| 美女国产高潮福利片在线看| h视频一区二区三区| 超色免费av| 在线观看免费高清a一片| 国产亚洲精品第一综合不卡 | 99视频精品全部免费 在线| 精品久久久久久久久av| 国产亚洲一区二区精品| 国产在线免费精品| 黑人高潮一二区| 少妇 在线观看| 国产在线视频一区二区| 日日撸夜夜添| 水蜜桃什么品种好| 不卡视频在线观看欧美| 亚洲av男天堂| 中文字幕亚洲精品专区| 69精品国产乱码久久久| 飞空精品影院首页| av福利片在线| 成年av动漫网址| 日本色播在线视频| 国产精品久久久久久久久免| 黄片无遮挡物在线观看| 亚洲av综合色区一区| 亚洲一区二区三区欧美精品| 秋霞在线观看毛片| 久久久午夜欧美精品| 精品一区二区三区视频在线| 成人毛片a级毛片在线播放| 午夜福利影视在线免费观看| 女人久久www免费人成看片| 免费大片18禁| 亚洲久久久国产精品| 国产精品一区二区在线不卡| 国产精品不卡视频一区二区| 成人18禁高潮啪啪吃奶动态图 | 免费不卡的大黄色大毛片视频在线观看| 精品99又大又爽又粗少妇毛片| 一区在线观看完整版| 麻豆精品久久久久久蜜桃| 久久鲁丝午夜福利片| 久久久国产一区二区| 一本一本综合久久| 欧美日韩国产mv在线观看视频| 亚洲欧美日韩卡通动漫| 另类亚洲欧美激情| 女人精品久久久久毛片| 97精品久久久久久久久久精品| 欧美+日韩+精品| 久久人人爽人人爽人人片va| 丝瓜视频免费看黄片| av播播在线观看一区| 国产片特级美女逼逼视频| www.av在线官网国产| 久久人人爽人人爽人人片va| 91精品国产九色| 亚洲欧美精品自产自拍| 午夜免费鲁丝| 亚洲少妇的诱惑av| 亚洲激情五月婷婷啪啪| 天堂8中文在线网| 中文字幕人妻熟人妻熟丝袜美| 亚洲一区二区三区欧美精品| 亚洲成人手机| 免费人成在线观看视频色| 国产一区二区三区av在线| 国产在视频线精品| 日日摸夜夜添夜夜添av毛片| 亚洲av.av天堂| 久久精品久久久久久久性| 熟女电影av网| 美女xxoo啪啪120秒动态图| 老女人水多毛片| 国产精品无大码| 久久久久精品性色| 女性生殖器流出的白浆| 91精品国产国语对白视频| 女人久久www免费人成看片| 久久国产亚洲av麻豆专区| 五月开心婷婷网| 日日爽夜夜爽网站| 国产高清不卡午夜福利| 免费观看的影片在线观看| 日本-黄色视频高清免费观看| 国产精品久久久久成人av| 精品国产一区二区三区久久久樱花| 丰满少妇做爰视频| 人体艺术视频欧美日本| 热99国产精品久久久久久7| 在线看a的网站| av女优亚洲男人天堂| 美女大奶头黄色视频| 九色亚洲精品在线播放| 制服丝袜香蕉在线| av在线播放精品| 五月玫瑰六月丁香| 日韩欧美精品免费久久| 国产精品不卡视频一区二区| 亚洲欧美成人综合另类久久久| 丰满饥渴人妻一区二区三| 久久99蜜桃精品久久| 婷婷成人精品国产| 亚洲欧洲国产日韩| tube8黄色片| 亚洲一级一片aⅴ在线观看| 午夜影院在线不卡| 国产亚洲午夜精品一区二区久久| 尾随美女入室| 亚洲,一卡二卡三卡| 国产白丝娇喘喷水9色精品| 精品少妇内射三级| 精品国产露脸久久av麻豆| 最近手机中文字幕大全| 成人毛片a级毛片在线播放| 97精品久久久久久久久久精品| 日产精品乱码卡一卡2卡三| 超碰97精品在线观看| 美女主播在线视频| 乱码一卡2卡4卡精品| 国产精品人妻久久久影院| 一级毛片aaaaaa免费看小| 国产亚洲最大av| 伦精品一区二区三区| xxxhd国产人妻xxx| 亚洲经典国产精华液单| 国产精品熟女久久久久浪| 22中文网久久字幕| 国产精品久久久久久精品古装| 最黄视频免费看| 自线自在国产av| 日日撸夜夜添| 一区二区日韩欧美中文字幕 | 国产国拍精品亚洲av在线观看| av视频免费观看在线观看| 久久午夜综合久久蜜桃| 少妇高潮的动态图| 大香蕉97超碰在线| 婷婷色麻豆天堂久久| 成人二区视频| 美女国产视频在线观看| h视频一区二区三区| 国产精品一区二区在线不卡| 岛国毛片在线播放| 国产精品久久久久成人av| 亚洲精品aⅴ在线观看| 在线免费观看不下载黄p国产| 国产精品无大码| 自线自在国产av| 18禁裸乳无遮挡动漫免费视频| 成年女人在线观看亚洲视频| 99九九线精品视频在线观看视频| 色94色欧美一区二区| 少妇猛男粗大的猛烈进出视频| 欧美+日韩+精品| 欧美bdsm另类| 亚洲av福利一区| 亚洲人成网站在线播| 黄片无遮挡物在线观看| 日韩制服骚丝袜av| 曰老女人黄片| 99视频精品全部免费 在线| 校园人妻丝袜中文字幕| 午夜影院在线不卡| 少妇人妻久久综合中文| 亚洲精品456在线播放app| 超色免费av| 少妇 在线观看| 亚洲av欧美aⅴ国产| 夜夜骑夜夜射夜夜干| 女性生殖器流出的白浆| 久久99精品国语久久久| 欧美精品亚洲一区二区| 免费观看a级毛片全部| 99久久人妻综合| 大片电影免费在线观看免费| 亚洲国产精品999| 日日摸夜夜添夜夜爱| 高清午夜精品一区二区三区| 少妇 在线观看| av福利片在线| 黄色欧美视频在线观看| 人人澡人人妻人| 国产精品一区二区在线观看99| 老司机影院成人| 只有这里有精品99| 18禁观看日本| 2018国产大陆天天弄谢| 免费看光身美女| 婷婷色综合www| 亚洲av国产av综合av卡| 国产欧美日韩一区二区三区在线 | 中文字幕人妻熟人妻熟丝袜美| 国产精品三级大全| 黄色一级大片看看| 中文乱码字字幕精品一区二区三区| av在线观看视频网站免费| 高清毛片免费看| 日韩av免费高清视频| 日日摸夜夜添夜夜添av毛片| 国产免费一区二区三区四区乱码| 国产精品 国内视频| 久久久久久久大尺度免费视频| 亚州av有码| 亚洲无线观看免费| 亚洲综合精品二区| 在线观看免费日韩欧美大片 | 男女国产视频网站| 亚洲中文av在线| 久久久久网色| 精品久久久久久久久亚洲| 亚洲av成人精品一区久久| 制服人妻中文乱码| 欧美日韩综合久久久久久| 国产黄色视频一区二区在线观看| 日韩精品有码人妻一区| 日韩一区二区视频免费看| av国产精品久久久久影院| 亚洲欧美精品自产自拍| 国产一级毛片在线| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲,欧美,日韩| 欧美日韩国产mv在线观看视频| 男的添女的下面高潮视频| 在线亚洲精品国产二区图片欧美 | 亚洲在久久综合| 99久久精品国产国产毛片| 亚洲精品,欧美精品| 精品亚洲成a人片在线观看| 欧美少妇被猛烈插入视频| videosex国产| 最近的中文字幕免费完整| 搡女人真爽免费视频火全软件| 99久久精品一区二区三区| 国产免费视频播放在线视频| 精品人妻偷拍中文字幕| 狠狠精品人妻久久久久久综合| 国产国拍精品亚洲av在线观看| 亚洲国产日韩一区二区| 亚洲国产色片| 成人国产麻豆网| 考比视频在线观看| 国产免费一区二区三区四区乱码| 这个男人来自地球电影免费观看 | 国产亚洲午夜精品一区二区久久| 国产精品.久久久| 国产成人aa在线观看| 搡女人真爽免费视频火全软件| av视频免费观看在线观看| 99久久综合免费| 亚洲一区二区三区欧美精品| 极品少妇高潮喷水抽搐| 搡老乐熟女国产| 国产熟女欧美一区二区| 国产视频内射| 97超视频在线观看视频| 久久精品国产亚洲av涩爱| 老司机影院成人| 看非洲黑人一级黄片| 在线观看免费视频网站a站| 久久韩国三级中文字幕| 一级,二级,三级黄色视频| 国产深夜福利视频在线观看| 伦理电影免费视频| 在线天堂最新版资源| 亚洲成人一二三区av| 午夜激情久久久久久久| 看免费成人av毛片| 最近2019中文字幕mv第一页| 美女国产视频在线观看| 亚洲,欧美,日韩| 五月天丁香电影| 国产国语露脸激情在线看| 男女无遮挡免费网站观看| 日日摸夜夜添夜夜爱| 91午夜精品亚洲一区二区三区| 久久精品国产鲁丝片午夜精品| 成人18禁高潮啪啪吃奶动态图 | 美女内射精品一级片tv| 国产极品天堂在线| 老女人水多毛片| 91在线精品国自产拍蜜月| 欧美精品一区二区大全| 久久久久久久久久久丰满| 成人18禁高潮啪啪吃奶动态图 | 国产免费一区二区三区四区乱码| 只有这里有精品99| 成人综合一区亚洲| 少妇猛男粗大的猛烈进出视频| 国产黄色免费在线视频| 18禁观看日本| 另类精品久久| av在线观看视频网站免费| 蜜桃国产av成人99| 国产极品粉嫩免费观看在线 | 伊人久久精品亚洲午夜| a 毛片基地| 99热网站在线观看| 久久精品国产亚洲av涩爱| 午夜激情av网站| 桃花免费在线播放| 日韩,欧美,国产一区二区三区| a 毛片基地| 色吧在线观看| 亚洲精品美女久久av网站| 亚洲欧美一区二区三区国产| 我的女老师完整版在线观看| 国产探花极品一区二区| 狂野欧美白嫩少妇大欣赏| 久久久久久久久大av| 色5月婷婷丁香| 久久久久网色| 美女国产视频在线观看| 日韩,欧美,国产一区二区三区| 少妇人妻精品综合一区二区| 久久久久视频综合| 成人免费观看视频高清| 日本黄色片子视频| 一级a做视频免费观看| 国产av码专区亚洲av| 纯流量卡能插随身wifi吗| 91久久精品国产一区二区三区| 亚洲欧美精品自产自拍| 亚洲图色成人|