• <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)引頭的解耦算法
    亚洲精华国产精华精| 亚洲九九香蕉| 国产精品久久久久久精品电影 | 亚洲在线自拍视频| 欧美丝袜亚洲另类 | 午夜精品在线福利| 亚洲国产精品久久男人天堂| 91九色精品人成在线观看| ponron亚洲| 国产精品二区激情视频| 午夜福利在线观看吧| 精品一区二区三区四区五区乱码| 成熟少妇高潮喷水视频| 国产精品久久电影中文字幕| 成熟少妇高潮喷水视频| 免费久久久久久久精品成人欧美视频| 午夜福利成人在线免费观看| 久久国产精品人妻蜜桃| 麻豆成人av在线观看| 国产亚洲精品第一综合不卡| 黄色视频不卡| 黄色片一级片一级黄色片| 看免费av毛片| 国产成人欧美| 国产精品九九99| 久久午夜亚洲精品久久| 黄色丝袜av网址大全| 国产伦一二天堂av在线观看| 最好的美女福利视频网| 国产成人精品在线电影| 日本vs欧美在线观看视频| 国产1区2区3区精品| 亚洲一区二区三区色噜噜| 国产一卡二卡三卡精品| 国产亚洲av嫩草精品影院| 亚洲男人的天堂狠狠| 窝窝影院91人妻| 午夜久久久久精精品| a级毛片在线看网站| 久久久久九九精品影院| 日日干狠狠操夜夜爽| 亚洲国产欧美日韩在线播放| 自线自在国产av| ponron亚洲| 日韩中文字幕欧美一区二区| 中文字幕久久专区| 高清毛片免费观看视频网站| 精品一区二区三区视频在线观看免费| 一二三四在线观看免费中文在| 欧美精品啪啪一区二区三区| 变态另类成人亚洲欧美熟女 | 在线观看66精品国产| 亚洲国产高清在线一区二区三 | 热99re8久久精品国产| 青草久久国产| 免费看美女性在线毛片视频| 在线观看一区二区三区| 午夜福利高清视频| 成人特级黄色片久久久久久久| 久久婷婷成人综合色麻豆| 美女午夜性视频免费| 久久亚洲精品不卡| 岛国视频午夜一区免费看| av有码第一页| 国内毛片毛片毛片毛片毛片| 1024视频免费在线观看| 国内精品久久久久久久电影| 国产免费男女视频| 午夜福利一区二区在线看| 1024视频免费在线观看| 成人欧美大片| 国产av在哪里看| 香蕉丝袜av| 成人国语在线视频| 国产精品免费视频内射| 19禁男女啪啪无遮挡网站| 99香蕉大伊视频| 中文字幕av电影在线播放| 午夜精品国产一区二区电影| 成人免费观看视频高清| 热99re8久久精品国产| 91在线观看av| 精品国产一区二区三区四区第35| 一本久久中文字幕| 丁香六月欧美| 亚洲欧美日韩无卡精品| 一级片免费观看大全| 成熟少妇高潮喷水视频| 动漫黄色视频在线观看| videosex国产| 人妻久久中文字幕网| 国产精品一区二区精品视频观看| 久久精品国产清高在天天线| 国产麻豆成人av免费视频| 亚洲中文日韩欧美视频| 村上凉子中文字幕在线| 亚洲国产欧美一区二区综合| 制服人妻中文乱码| 激情视频va一区二区三区| 国产精品免费一区二区三区在线| 久久九九热精品免费| 精品国产乱子伦一区二区三区| 国产精品综合久久久久久久免费 | 精品第一国产精品| 自线自在国产av| 久久中文字幕一级| 在线播放国产精品三级| 老熟妇仑乱视频hdxx| 亚洲狠狠婷婷综合久久图片| 日韩有码中文字幕| 国产片内射在线| 久久影院123| 国产精华一区二区三区| 久久久久国产一级毛片高清牌| 国内精品久久久久精免费| 搞女人的毛片| 嫩草影视91久久| 他把我摸到了高潮在线观看| 极品教师在线免费播放| 桃红色精品国产亚洲av| 亚洲第一欧美日韩一区二区三区| 我的亚洲天堂| 一本综合久久免费| 日本黄色视频三级网站网址| 欧美在线黄色| 色播亚洲综合网| 亚洲精品国产色婷婷电影| 在线观看www视频免费| 久久久水蜜桃国产精品网| 久久人妻熟女aⅴ| 日本三级黄在线观看| 变态另类成人亚洲欧美熟女 | 精品久久久久久成人av| 久久婷婷成人综合色麻豆| 亚洲人成77777在线视频| 18禁黄网站禁片午夜丰满| svipshipincom国产片| 啦啦啦韩国在线观看视频| 欧美成人午夜精品| 欧美亚洲日本最大视频资源| 村上凉子中文字幕在线| 男女做爰动态图高潮gif福利片 | 日本一区二区免费在线视频| e午夜精品久久久久久久| 亚洲欧美日韩高清在线视频| 国产精品国产高清国产av| av天堂在线播放| 欧美黄色片欧美黄色片| 极品教师在线免费播放| 久久久国产成人免费| 欧美乱码精品一区二区三区| av在线播放免费不卡| 十八禁人妻一区二区| 欧美最黄视频在线播放免费| 国产一区在线观看成人免费| 一级作爱视频免费观看| 老熟妇乱子伦视频在线观看| 国产欧美日韩一区二区三| 天天添夜夜摸| 久久九九热精品免费| 久久精品国产99精品国产亚洲性色 | 国产亚洲av嫩草精品影院| 亚洲成国产人片在线观看| 看片在线看免费视频| 亚洲一区二区三区色噜噜| 欧美老熟妇乱子伦牲交| a级毛片在线看网站| 女人精品久久久久毛片| 好看av亚洲va欧美ⅴa在| 12—13女人毛片做爰片一| 亚洲av美国av| 国产成人影院久久av| 成年女人毛片免费观看观看9| 9191精品国产免费久久| 亚洲精品国产区一区二| 日韩欧美一区二区三区在线观看| 女性生殖器流出的白浆| www.www免费av| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩瑟瑟在线播放| 精品久久久久久,| 欧美日本中文国产一区发布| 淫秽高清视频在线观看| 国产麻豆成人av免费视频| 午夜免费鲁丝| 动漫黄色视频在线观看| 最好的美女福利视频网| 天堂√8在线中文| 一级,二级,三级黄色视频| 欧美成人免费av一区二区三区| 久久国产亚洲av麻豆专区| 91字幕亚洲| 欧美日本中文国产一区发布| 别揉我奶头~嗯~啊~动态视频| 午夜免费鲁丝| 亚洲成人免费电影在线观看| 国产精品香港三级国产av潘金莲| 国产三级黄色录像| 久99久视频精品免费| www国产在线视频色| 国产欧美日韩一区二区三| 99在线视频只有这里精品首页| 免费一级毛片在线播放高清视频 | 最近最新中文字幕大全免费视频| 十分钟在线观看高清视频www| 老熟妇乱子伦视频在线观看| 99热只有精品国产| 免费在线观看黄色视频的| 韩国av一区二区三区四区| 99re在线观看精品视频| 国产精品久久视频播放| 欧美一级毛片孕妇| 欧美性长视频在线观看| av天堂在线播放| 亚洲五月天丁香| 夜夜躁狠狠躁天天躁| 精品不卡国产一区二区三区| 人人妻,人人澡人人爽秒播| 黑丝袜美女国产一区| 久久久精品国产亚洲av高清涩受| 伦理电影免费视频| 日韩免费av在线播放| 国产精品精品国产色婷婷| 亚洲一区高清亚洲精品| www日本在线高清视频| 男人舔女人下体高潮全视频| 如日韩欧美国产精品一区二区三区| 看黄色毛片网站| 97人妻精品一区二区三区麻豆 | 亚洲av熟女| 日韩欧美在线二视频| 50天的宝宝边吃奶边哭怎么回事| 午夜免费成人在线视频| 免费女性裸体啪啪无遮挡网站| 欧美国产日韩亚洲一区| 成人av一区二区三区在线看| 欧美日韩中文字幕国产精品一区二区三区 | 色尼玛亚洲综合影院| 国产一区二区三区视频了| 久久久久久国产a免费观看| 非洲黑人性xxxx精品又粗又长| av电影中文网址| 美女大奶头视频| 国产亚洲欧美在线一区二区| 国产野战对白在线观看| 亚洲精品国产精品久久久不卡| 男人的好看免费观看在线视频 | 亚洲国产欧美日韩在线播放| 91在线观看av| 婷婷精品国产亚洲av在线| 叶爱在线成人免费视频播放| 久久久国产成人精品二区| 久久久久国产精品人妻aⅴ院| 亚洲第一av免费看| 国产精品久久久久久人妻精品电影| 婷婷丁香在线五月| 777久久人妻少妇嫩草av网站| 黄色女人牲交| 亚洲人成电影观看| 啦啦啦观看免费观看视频高清 | 最近最新免费中文字幕在线| 国产成人一区二区三区免费视频网站| 国产精品久久久久久人妻精品电影| 亚洲第一欧美日韩一区二区三区| 亚洲自偷自拍图片 自拍| 日韩欧美国产一区二区入口| videosex国产| 在线av久久热| 9热在线视频观看99| 丰满的人妻完整版| 69av精品久久久久久| 成人三级做爰电影| 激情在线观看视频在线高清| 他把我摸到了高潮在线观看| 欧美日韩福利视频一区二区| 色综合婷婷激情| 色在线成人网| 日韩视频一区二区在线观看| 精品一区二区三区四区五区乱码| 欧美 亚洲 国产 日韩一| 巨乳人妻的诱惑在线观看| 国产亚洲精品久久久久5区| 国产麻豆69| 此物有八面人人有两片| www.自偷自拍.com| 久久精品人人爽人人爽视色| 成人国语在线视频| 欧美日韩精品网址| 亚洲久久久国产精品| 别揉我奶头~嗯~啊~动态视频| 成人精品一区二区免费| 18禁观看日本| 国产一区二区激情短视频| 国产麻豆69| 9191精品国产免费久久| 乱人伦中国视频| 欧美日韩亚洲综合一区二区三区_| 一卡2卡三卡四卡精品乱码亚洲| 黄网站色视频无遮挡免费观看| 国产三级黄色录像| 69精品国产乱码久久久| 国产亚洲精品一区二区www| 久久亚洲精品不卡| 90打野战视频偷拍视频| 人妻久久中文字幕网| 激情视频va一区二区三区| or卡值多少钱| www国产在线视频色| 色婷婷久久久亚洲欧美| 久热这里只有精品99| 午夜福利视频1000在线观看 | 亚洲欧美精品综合一区二区三区| 十八禁网站免费在线| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品合色在线| 国产高清激情床上av| 国产一区二区在线av高清观看| 久9热在线精品视频| 精品国产乱子伦一区二区三区| 黄色 视频免费看| 欧美日韩亚洲综合一区二区三区_| 亚洲男人的天堂狠狠| 999精品在线视频| 嫩草影院精品99| 国产亚洲欧美精品永久| 亚洲成人免费电影在线观看| 日韩三级视频一区二区三区| 99久久精品国产亚洲精品| 一级a爱视频在线免费观看| 99国产精品免费福利视频| 极品人妻少妇av视频| 美女国产高潮福利片在线看| 美女高潮喷水抽搐中文字幕| 精品一区二区三区四区五区乱码| 一边摸一边抽搐一进一出视频| 国产精品一区二区三区四区久久 | 亚洲中文av在线| 十分钟在线观看高清视频www| 精品少妇一区二区三区视频日本电影| videosex国产| 法律面前人人平等表现在哪些方面| 久久久久久人人人人人| 亚洲av第一区精品v没综合| 日日摸夜夜添夜夜添小说| 国产成年人精品一区二区| 级片在线观看| 18禁国产床啪视频网站| 亚洲在线自拍视频| 99在线人妻在线中文字幕| 校园春色视频在线观看| ponron亚洲| 黄频高清免费视频| 9色porny在线观看| 老司机靠b影院| 欧美在线一区亚洲| 看免费av毛片| 久久久久久人人人人人| 免费av毛片视频| 国产精品 国内视频| 法律面前人人平等表现在哪些方面| 欧美日韩瑟瑟在线播放| 99国产精品免费福利视频| 国产97色在线日韩免费| 麻豆成人av在线观看| www.999成人在线观看| 热99re8久久精品国产| 精品福利观看| 1024香蕉在线观看| 国产国语露脸激情在线看| 一二三四在线观看免费中文在| 变态另类丝袜制服| 色av中文字幕| 免费无遮挡裸体视频| 亚洲片人在线观看| 亚洲精华国产精华精| 高清在线国产一区| 国产亚洲精品久久久久久毛片| 国产免费av片在线观看野外av| 欧美日韩一级在线毛片| 真人做人爱边吃奶动态| 在线观看66精品国产| 乱人伦中国视频| 91麻豆精品激情在线观看国产| 欧美日韩中文字幕国产精品一区二区三区 | 12—13女人毛片做爰片一| 丁香欧美五月| 国产精品二区激情视频| 欧美日韩乱码在线| 如日韩欧美国产精品一区二区三区| 首页视频小说图片口味搜索| 国产aⅴ精品一区二区三区波| 色综合婷婷激情| 精品少妇一区二区三区视频日本电影| 嫩草影视91久久| 国产主播在线观看一区二区| АⅤ资源中文在线天堂| 日日干狠狠操夜夜爽| 别揉我奶头~嗯~啊~动态视频| 人妻久久中文字幕网| 精品国内亚洲2022精品成人| 亚洲国产精品合色在线| 日韩免费av在线播放| 久久精品国产亚洲av香蕉五月| 在线天堂中文资源库| 国产成人精品无人区| 国产欧美日韩综合在线一区二区| 国产成人啪精品午夜网站| 99国产精品99久久久久| 国产亚洲av高清不卡| 国产亚洲精品第一综合不卡| 岛国视频午夜一区免费看| 免费高清在线观看日韩| 欧美成人一区二区免费高清观看 | 久久久久久久午夜电影| 午夜福利在线观看吧| 成人国语在线视频| 国产99久久九九免费精品| 老熟妇乱子伦视频在线观看| 亚洲av成人不卡在线观看播放网| 精品人妻在线不人妻| 久久中文字幕人妻熟女| 波多野结衣一区麻豆| 法律面前人人平等表现在哪些方面| 亚洲中文日韩欧美视频| 久久久久久免费高清国产稀缺| 久久中文字幕人妻熟女| √禁漫天堂资源中文www| 动漫黄色视频在线观看| 视频区欧美日本亚洲| 国内精品久久久久久久电影| www日本在线高清视频| 国产91精品成人一区二区三区| 国产区一区二久久| 9191精品国产免费久久| 亚洲专区中文字幕在线| 岛国视频午夜一区免费看| av超薄肉色丝袜交足视频| 最新在线观看一区二区三区| 免费一级毛片在线播放高清视频 | 日韩欧美免费精品| 极品人妻少妇av视频| 亚洲国产高清在线一区二区三 | 久久午夜亚洲精品久久| 色精品久久人妻99蜜桃| 精品欧美一区二区三区在线| 精品久久久精品久久久| 黑人巨大精品欧美一区二区蜜桃| 亚洲激情在线av| av网站免费在线观看视频| 狂野欧美激情性xxxx| 久久精品国产综合久久久| 女人被躁到高潮嗷嗷叫费观| 欧美黑人欧美精品刺激| 十八禁人妻一区二区| www.熟女人妻精品国产| 欧美精品亚洲一区二区| 久久狼人影院| 高清黄色对白视频在线免费看| 最新在线观看一区二区三区| 亚洲av成人一区二区三| 伊人久久大香线蕉亚洲五| 啦啦啦 在线观看视频| 国产熟女午夜一区二区三区| 午夜免费鲁丝| 亚洲精品国产区一区二| 88av欧美| 两性夫妻黄色片| 久久久久久久久免费视频了| 又紧又爽又黄一区二区| 国产精品久久久久久精品电影 | 91精品三级在线观看| 两性夫妻黄色片| 夜夜躁狠狠躁天天躁| 不卡一级毛片| 久久久水蜜桃国产精品网| 男人的好看免费观看在线视频 | ponron亚洲| 欧美激情极品国产一区二区三区| 69精品国产乱码久久久| 日日摸夜夜添夜夜添小说| 久久人人97超碰香蕉20202| 欧美日韩中文字幕国产精品一区二区三区 | 天堂√8在线中文| 99精品久久久久人妻精品| 中文字幕av电影在线播放| 精品电影一区二区在线| 国产精品日韩av在线免费观看 | av免费在线观看网站| 精品人妻在线不人妻| 极品人妻少妇av视频| 亚洲五月色婷婷综合| 久久久国产精品麻豆| 国产高清videossex| 91成人精品电影| 波多野结衣一区麻豆| 国产乱人伦免费视频| 国产一区二区三区在线臀色熟女| 伦理电影免费视频| 不卡一级毛片| 亚洲狠狠婷婷综合久久图片| 狂野欧美激情性xxxx| 亚洲欧美精品综合一区二区三区| 久久中文看片网| 免费在线观看黄色视频的| 亚洲午夜理论影院| 国产真人三级小视频在线观看| 色精品久久人妻99蜜桃| 91老司机精品| 黑人巨大精品欧美一区二区蜜桃| 免费观看人在逋| 国产激情久久老熟女| 日本黄色视频三级网站网址| 在线国产一区二区在线| 天天躁狠狠躁夜夜躁狠狠躁| 电影成人av| 高潮久久久久久久久久久不卡| 999久久久国产精品视频| 十八禁人妻一区二区| 欧美性长视频在线观看| 露出奶头的视频| 天堂影院成人在线观看| 99国产极品粉嫩在线观看| 制服诱惑二区| 曰老女人黄片| 国内毛片毛片毛片毛片毛片| 中文字幕av电影在线播放| 国产私拍福利视频在线观看| 午夜福利成人在线免费观看| 日本五十路高清| 波多野结衣一区麻豆| 老熟妇乱子伦视频在线观看| 老汉色av国产亚洲站长工具| www国产在线视频色| 亚洲中文字幕日韩| 黑丝袜美女国产一区| 亚洲第一电影网av| 黑人巨大精品欧美一区二区蜜桃| 大型av网站在线播放| 亚洲情色 制服丝袜| 婷婷精品国产亚洲av在线| a级毛片在线看网站| 真人做人爱边吃奶动态| 日韩 欧美 亚洲 中文字幕| 国产av又大| 欧美性长视频在线观看| 日本精品一区二区三区蜜桃| 欧美黄色淫秽网站| 日韩中文字幕欧美一区二区| 国产精品自产拍在线观看55亚洲| www.熟女人妻精品国产| 久久欧美精品欧美久久欧美| 最近最新免费中文字幕在线| 一区二区三区国产精品乱码| 精品人妻在线不人妻| 非洲黑人性xxxx精品又粗又长| 久热爱精品视频在线9| 亚洲av电影在线进入| 午夜久久久久精精品| 久久人人精品亚洲av| 成人手机av| 女性被躁到高潮视频| 19禁男女啪啪无遮挡网站| 精品久久久久久成人av| 亚洲 国产 在线| 色精品久久人妻99蜜桃| 亚洲午夜理论影院| 88av欧美| 精品久久久久久成人av| 亚洲av成人不卡在线观看播放网| 久久久久国内视频| 免费在线观看视频国产中文字幕亚洲| 欧美黑人欧美精品刺激| 亚洲av电影不卡..在线观看| 大型黄色视频在线免费观看| 1024视频免费在线观看| 国产伦一二天堂av在线观看| 91精品三级在线观看| 精品久久久久久成人av| 国产av在哪里看| 后天国语完整版免费观看| 精品福利观看| 亚洲九九香蕉| 99精品在免费线老司机午夜| 午夜福利一区二区在线看| 国产区一区二久久| 一进一出抽搐gif免费好疼| a级毛片在线看网站| 搞女人的毛片| 久久久久久大精品| 亚洲在线自拍视频| 久久婷婷人人爽人人干人人爱 | 国产亚洲精品久久久久久毛片| 久久人妻熟女aⅴ| 热re99久久国产66热| 91精品三级在线观看| 久久久久精品国产欧美久久久| 一级a爱片免费观看的视频| 国产极品粉嫩免费观看在线| 动漫黄色视频在线观看| 99riav亚洲国产免费| 国产99白浆流出| 中文字幕人妻熟女乱码| 国产欧美日韩一区二区三区在线| 国产国语露脸激情在线看| 亚洲伊人色综图| www国产在线视频色| 精品国产一区二区三区四区第35| 美女扒开内裤让男人捅视频| 男人舔女人下体高潮全视频| 国产高清videossex| 一区二区日韩欧美中文字幕| 中文字幕色久视频| 亚洲成人精品中文字幕电影| 日韩一卡2卡3卡4卡2021年| 久久人人97超碰香蕉20202|