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

    抗擾動(dòng)的捷聯(lián)慣導(dǎo)系統(tǒng)回溯參數(shù)辨識(shí)對(duì)準(zhǔn)法

    2015-06-15 19:19:45張朝飛羅建軍侯永利
    關(guān)鍵詞:捷聯(lián)慣導(dǎo)對(duì)準(zhǔn)

    張朝飛,羅建軍,侯永利

    (1. 西北工業(yè)大學(xué) 航天學(xué)院,西安 710072;2. 航天飛行動(dòng)力學(xué)技術(shù)重點(diǎn)實(shí)驗(yàn)室,西安 710072;3. 中國(guó)航天第十六研究所,西安 710100)

    抗擾動(dòng)的捷聯(lián)慣導(dǎo)系統(tǒng)回溯參數(shù)辨識(shí)對(duì)準(zhǔn)法

    張朝飛1,2,羅建軍1,2,侯永利3

    (1. 西北工業(yè)大學(xué) 航天學(xué)院,西安 710072;2. 航天飛行動(dòng)力學(xué)技術(shù)重點(diǎn)實(shí)驗(yàn)室,西安 710072;3. 中國(guó)航天第十六研究所,西安 710100)

    針對(duì)捷聯(lián)慣導(dǎo)系統(tǒng)參數(shù)辨識(shí)對(duì)準(zhǔn)法航向角估計(jì)時(shí)間長(zhǎng)且收斂速度受到計(jì)算初值影響大的問題,提出了將對(duì)準(zhǔn)用的陀螺和加速度計(jì)采樣數(shù)據(jù)存儲(chǔ)后按時(shí)間進(jìn)行正向和逆向排列,形成正向序列和逆向序列。由正向?qū)Ш降慕馑愎酵茖?dǎo)出逆向?qū)Ш降慕馑愎?,利用公式的一致性提出了回溯誤差參數(shù)辨識(shí)對(duì)準(zhǔn)法。利用抗擾動(dòng)的慣性系對(duì)準(zhǔn)算法為回溯參數(shù)辨識(shí)法提供對(duì)準(zhǔn)初值,應(yīng)用回溯參數(shù)辨識(shí)法對(duì)誤差參數(shù)反復(fù)進(jìn)行逆向和正向辨識(shí)后對(duì)對(duì)準(zhǔn)結(jié)果進(jìn)行修正,能夠有效提高對(duì)準(zhǔn)精度。實(shí)驗(yàn)結(jié)果表明,利用相同時(shí)間的采樣數(shù)據(jù),該算法航向角的對(duì)準(zhǔn)誤差僅為改進(jìn)前算法的25%。

    捷聯(lián)慣導(dǎo)系統(tǒng);抗擾動(dòng);回溯;參數(shù)辨識(shí)對(duì)準(zhǔn)

    初始對(duì)準(zhǔn)技術(shù)是捷聯(lián)慣性導(dǎo)航系統(tǒng)中一項(xiàng)重要技術(shù),對(duì)準(zhǔn)的快速性和高精度是其主要的指標(biāo)。目前常用的對(duì)準(zhǔn)方法有多位置對(duì)準(zhǔn)法和單位置對(duì)準(zhǔn)法。其中,多位置對(duì)準(zhǔn)算法能夠提高對(duì)準(zhǔn)參數(shù)的可觀測(cè)性,但是需要增加高精度的轉(zhuǎn)位機(jī)構(gòu)才能實(shí)現(xiàn)[1-4]。單位置對(duì)準(zhǔn)算法有:① 經(jīng)典三階調(diào)平回路及方位羅經(jīng)對(duì)準(zhǔn)法[5-6],該算法收斂緩慢,適用于艦船等長(zhǎng)時(shí)間對(duì)準(zhǔn),不利于快速對(duì)準(zhǔn);② 卡爾曼濾波對(duì)準(zhǔn)算法[7-8],該算法狀態(tài)方程維數(shù)高,計(jì)算量大且收斂性受模型參數(shù)的影響較大;③ 參數(shù)辨識(shí)對(duì)準(zhǔn)法[9],該算法原理簡(jiǎn)單,易于實(shí)現(xiàn),因此廣泛應(yīng)用于捷聯(lián)慣導(dǎo)系統(tǒng)的初始對(duì)準(zhǔn)中,但該算法誤差參數(shù)收斂速度慢,航向角對(duì)準(zhǔn)時(shí)間受速度誤差的影響,辨識(shí)精度差。

    因此,本文在文獻(xiàn)[9]的基礎(chǔ)上,改進(jìn)了將陀螺和加速度計(jì)的采樣數(shù)據(jù)按時(shí)間順序先進(jìn)行粗對(duì)準(zhǔn)再進(jìn)行精對(duì)準(zhǔn)的單向?qū)?zhǔn)方式,提出了一種抗擾動(dòng)的捷聯(lián)慣導(dǎo)系統(tǒng)回溯參數(shù)辨識(shí)對(duì)準(zhǔn)法。

    通常提高對(duì)準(zhǔn)精度的方法是減少對(duì)準(zhǔn)時(shí)的干擾和增加對(duì)準(zhǔn)時(shí)間,使得初始對(duì)準(zhǔn)誤差在導(dǎo)航誤差中的所占的比重隨之增加,從而易于被辨識(shí)。但對(duì)準(zhǔn)時(shí)的干擾通常難以避免,而增加對(duì)準(zhǔn)時(shí)間也與對(duì)準(zhǔn)的快速性要求相矛盾。本文通過充分利用對(duì)準(zhǔn)時(shí)間內(nèi)的采樣數(shù)據(jù)所包含的導(dǎo)航誤差信息,而非單純地增加對(duì)準(zhǔn)時(shí)間,將對(duì)準(zhǔn)用的采樣數(shù)據(jù)存儲(chǔ)起來,按時(shí)間順序進(jìn)行正向和逆向排序后形成正向序列和逆向序列;然后利用抗擾動(dòng)的慣性系對(duì)準(zhǔn)算法對(duì)正向序列進(jìn)行粗對(duì)準(zhǔn),得到誤差較小的粗對(duì)準(zhǔn)結(jié)果;其次根據(jù)回溯參數(shù)辨識(shí)法對(duì)逆向序列和正向序列迭代進(jìn)行對(duì)準(zhǔn),最終達(dá)到加快誤差參數(shù)收斂速度和提高航向角對(duì)準(zhǔn)精度的目的。

    1 慣性系粗對(duì)準(zhǔn)

    載體坐標(biāo)系(b系)選擇前右下坐標(biāo)系。導(dǎo)航坐標(biāo)系(n系)為當(dāng)?shù)氐乩碜鴺?biāo)系,本文中為東北天坐標(biāo)系。慣性坐標(biāo)系(ib0系)在粗對(duì)準(zhǔn)開始時(shí)刻t0將b系在慣性空間凝固為ib0系,即。初始對(duì)準(zhǔn)姿態(tài)矩陣為其中:

    式中:L和λ分別為當(dāng)?shù)氐木暥群徒?jīng)度,ωie為地球自轉(zhuǎn)角速率,t為對(duì)準(zhǔn)時(shí)間。C利用陀螺的輸出利用等效旋轉(zhuǎn)矢量法進(jìn)行實(shí)時(shí)更新,因此求解 C的核心問題轉(zhuǎn)換為 C的計(jì)算問題。取 t<t ≤ t <t ,其中t為012m0對(duì)準(zhǔn)起始時(shí)刻,tm為對(duì)準(zhǔn)結(jié)束時(shí)刻,計(jì)算得:

    其中:

    2 回溯參數(shù)辨識(shí)法精對(duì)準(zhǔn)

    2.1 正向?qū)Ш剿惴ń榻B

    等效旋轉(zhuǎn)矢量微分方程的近似形式[10-11]為

    式中:Φ為載體系從tk-1時(shí)刻到tk時(shí)刻角位置變化所對(duì)應(yīng)的等效旋轉(zhuǎn)矢量,ω為該時(shí)間段內(nèi)的旋轉(zhuǎn)角速率。

    在靜基座或者晃動(dòng)基座對(duì)準(zhǔn)時(shí),忽略圓錐誤差補(bǔ)償項(xiàng),采用單子樣算法進(jìn)行姿態(tài)解算,得:

    其中:

    Ts= tk- tk-1為姿態(tài)更新周期。

    速度微分方程為

    求解得:

    位置微分方程為

    求解得:

    2.2 逆向?qū)Ш剿惴ㄍ茖?dǎo)

    若正向?qū)Ш介L(zhǎng)度為m,如圖1所示。以m時(shí)刻的姿態(tài)、速度取反和位置作為逆向?qū)Ш降某跏甲藨B(tài)、初始速度和初始位置,將陀螺采樣和地球自轉(zhuǎn)角速率取反,并且對(duì)陀螺和加速度計(jì)的采樣數(shù)據(jù)進(jìn)行逆序處理,形成逆向數(shù)據(jù)序列。

    圖1 逆向?qū)Ш剿惴ㄔ韴DFig.1 Principle of negative-direction navigation algorithm

    由公式(3)可以推導(dǎo)出,在逆向?qū)Ш竭^程中:

    若不考慮計(jì)算誤差,逆向?qū)Ш綍r(shí)各時(shí)刻的姿態(tài)信息、位置信息與正向?qū)Ш皆趯?duì)應(yīng)時(shí)刻上的姿態(tài)信息、位置信息相同,速度信息大小相同,方向相反。

    由上述推導(dǎo)可知,正向?qū)Ш胶湍嫦驅(qū)Ш降慕馑愎较嗤?,因此其系統(tǒng)誤差方程也相同,逆向參數(shù)辨識(shí)法精對(duì)準(zhǔn)的方程與正向參數(shù)辨識(shí)精對(duì)準(zhǔn)的方程完全相同,僅對(duì)準(zhǔn)過程中的導(dǎo)航信息不同。

    2.3 回溯參數(shù)辨識(shí)法精對(duì)準(zhǔn)

    完成慣性系粗對(duì)準(zhǔn)后,捷聯(lián)數(shù)學(xué)平臺(tái)的失準(zhǔn)角φE、 φN、 φU均為小角度。考慮陀螺常值漂移ε和加速度計(jì)零偏▽,捷聯(lián)慣導(dǎo)系統(tǒng)的失準(zhǔn)角誤差方程可以寫為[9]

    式中,uE、uN、uU為中間控制量:

    進(jìn)一步,可求得水平速度誤差方程為

    式中,VDE、VDN分別為對(duì)做簡(jiǎn)諧波動(dòng)的等效干擾加速度的積分,沒有隨時(shí)間增長(zhǎng)的趨勢(shì)??蓪⑹?8)改寫成:

    式(11)中,Ts為采樣周期, a1E=▽E- gφN0、

    速度誤差中辨識(shí)出,得:

    由以上分析可以得出抗擾動(dòng)的回溯參數(shù)辨識(shí)對(duì)準(zhǔn)法的具體步驟如下:

    Step 1. 先保存慣組從t0到tm時(shí)刻陀螺和加速度計(jì)的輸出,進(jìn)行慣性系粗對(duì)準(zhǔn),得到逆向?qū)?zhǔn)的初值C,開始進(jìn)行回溯參數(shù)辨識(shí)法精對(duì)準(zhǔn);

    Step 2. 以結(jié)束時(shí)刻tm為起始時(shí)刻,對(duì)陀螺和加速度計(jì)輸出數(shù)據(jù)進(jìn)行逆序處理,以慣性系粗對(duì)準(zhǔn)得到的姿態(tài)矩陣 C作為姿態(tài)矩陣的初值,根據(jù)逆向參數(shù)辨識(shí)法進(jìn)行精對(duì)準(zhǔn),得到精對(duì)準(zhǔn)后的姿態(tài)矩陣;

    Step 4. 執(zhí)行Step 2和Step 3,即完成一次回溯參數(shù)辨識(shí)法精對(duì)準(zhǔn),可以反復(fù)進(jìn)行多次,直到姿態(tài)矩陣解算得到的航向角逐漸收斂。

    整個(gè)過程如圖2所示。

    圖2 抗擾動(dòng)回溯參數(shù)辨識(shí)對(duì)準(zhǔn)法步驟Fig.2 Alignment approach of anti-interference backtracking parameter identification

    3 實(shí)驗(yàn)驗(yàn)證

    實(shí)驗(yàn)選用航天十六所生產(chǎn)的某型激光捷聯(lián)慣組,激光陀螺儀等效漂移(3σ)為 0.003 (o)/h,加速度計(jì)等效零偏(3σ)為5×10-5g0,慣組固定于實(shí)驗(yàn)車上,在不同的時(shí)間進(jìn)行6次實(shí)驗(yàn)。實(shí)驗(yàn)車靜止于平坦路面上,發(fā)動(dòng)機(jī)處于啟動(dòng)狀態(tài),上電30 min后進(jìn)行數(shù)據(jù)采集,期間人為進(jìn)行開關(guān)車門及上下車的動(dòng)作,采集 300 s數(shù)據(jù)后停止施加干擾并關(guān)閉發(fā)動(dòng)機(jī),再靜態(tài)采集300 s數(shù)據(jù)。實(shí)驗(yàn)結(jié)束后用G2000型陀螺經(jīng)緯儀對(duì)激光捷聯(lián)慣組的航向角進(jìn)行測(cè)量,作為航向角的真值。

    分別采用抗擾動(dòng)的回溯參數(shù)辨識(shí)法和文獻(xiàn)[9]中的參數(shù)辨識(shí)法對(duì)采樣數(shù)據(jù)進(jìn)行對(duì)準(zhǔn),表1為利用不同方法進(jìn)行對(duì)準(zhǔn)的實(shí)驗(yàn)結(jié)果,其中抗擾動(dòng)回溯參數(shù)辨識(shí)法為對(duì)采樣數(shù)據(jù)進(jìn)行一次回溯參數(shù)辨識(shí)精對(duì)準(zhǔn)的結(jié)果。

    由表1中可以看出,利用300 s的采樣數(shù)據(jù)進(jìn)行對(duì)準(zhǔn),抗干擾回溯參數(shù)辨識(shí)法的對(duì)準(zhǔn)誤差為參數(shù)辨識(shí)法誤差的25%,亦小于參數(shù)辨識(shí)法600 s的對(duì)準(zhǔn)誤差。圖3~圖5為利抗擾動(dòng)的回溯參數(shù)辨識(shí)法和用參數(shù)辨識(shí)法進(jìn)行對(duì)準(zhǔn)時(shí)的對(duì)準(zhǔn)誤差,由圖中可以看出,抗擾動(dòng)的回溯參數(shù)辨識(shí)法的對(duì)準(zhǔn)誤差明顯小于參數(shù)辨識(shí)法。

    表1 三種對(duì)準(zhǔn)方法的航向角對(duì)準(zhǔn)誤差Tab.1 Azimuth alignment errors of three alignment algorithms

    圖3 俯仰角誤差對(duì)比Fig.3 Comparison on pitch errors

    圖4 橫滾角誤差對(duì)比Fig.4 Comparison on roll errors

    圖5 航向角誤差對(duì)比Fig.5 Comparison on azimuth errors

    4 結(jié) 論

    慣性系的粗對(duì)準(zhǔn)能夠有效地隔離外界擾動(dòng)對(duì)捷聯(lián)慣導(dǎo)系統(tǒng)初始對(duì)準(zhǔn)的影響。因此,在有擾動(dòng)的環(huán)境下,利用慣性系粗對(duì)準(zhǔn)能夠?yàn)榛厮輩?shù)辨識(shí)精對(duì)準(zhǔn)提供較好的初始對(duì)準(zhǔn)值。另外,粗對(duì)準(zhǔn)的精度與對(duì)準(zhǔn)時(shí)間有很大關(guān)系,在一定程度上,對(duì)準(zhǔn)的時(shí)間越長(zhǎng),粗對(duì)準(zhǔn)的精度越高??箶_動(dòng)的回溯參數(shù)辨識(shí)法將對(duì)準(zhǔn)用的采樣數(shù)據(jù)按時(shí)間進(jìn)行正向和逆向處理,可以將整個(gè)正向的時(shí)間序列用于粗對(duì)準(zhǔn),單向的參數(shù)辨識(shí)對(duì)準(zhǔn)法僅能提供較短時(shí)間的采樣數(shù)據(jù)用于粗對(duì)準(zhǔn),因此,相對(duì)于單向的參數(shù)辨識(shí)對(duì)準(zhǔn)法而言,回溯參數(shù)辨識(shí)法可以有效地提高粗對(duì)準(zhǔn)的精度。

    回溯參數(shù)辨識(shí)法的本質(zhì)是提高對(duì)準(zhǔn)用的采樣數(shù)據(jù)的利用率,反復(fù)對(duì)采樣數(shù)據(jù)進(jìn)行逆向和正向?qū)Ш?,然后利用?dǎo)航過程中產(chǎn)生的速度誤差去估計(jì)對(duì)準(zhǔn)失準(zhǔn)角。理論分析可得,經(jīng)過多次的回溯,初始對(duì)準(zhǔn)可以達(dá)到極限精度。但是,根據(jù)實(shí)驗(yàn)驗(yàn)證,初始對(duì)準(zhǔn)的精度并不總隨著回溯次數(shù)增加而提高,當(dāng)回溯超過一定次數(shù)后,導(dǎo)航解算的速度誤差不再減小,而對(duì)準(zhǔn)精度也就不再提高。因此,在實(shí)際的應(yīng)用過程中,可以根據(jù)捷聯(lián)慣導(dǎo)系統(tǒng)的性能進(jìn)行有限次的回溯,以達(dá)到對(duì)準(zhǔn)精度和計(jì)算復(fù)雜度的平衡。

    本文中的方法不僅適用于靜基座對(duì)準(zhǔn),對(duì)于動(dòng)基座或行進(jìn)間對(duì)準(zhǔn),也有一定的參考價(jià)值,可以將引入的外界參考信息作為導(dǎo)航解算的真值,從導(dǎo)航誤差中辨識(shí)出初始對(duì)準(zhǔn)的誤差。

    (References):

    [1] 付強(qiáng)文, 秦永元, 張金紅, 等. 捷聯(lián)慣導(dǎo)系統(tǒng)快速最小二乘精對(duì)準(zhǔn)方法[J]. 中國(guó)慣性技術(shù)學(xué)報(bào), 2012, 20(3): 278-282. Fu Qiang-wen, Qin Yong-yuan, Zhang Jin-hong, et al. Rapid recursive least-square fine alignment method for SINS[J]. Journal of Chinese Inertial Technology, 2012, 20(3): 278-282.

    [2] Ma T, Gao Y. Application and comparison of two methods for alignment of FOG SINS[C]//2010 8th World Congress on Intelligent Control and Automation. IEEE, 2010: 3581-3584.

    [3] Shang S T, Gao W S. Application of adaptive Kalman filter technique in initial alignment of single-axial rotation strap-down inertial navigation system[C]//Advanced Materials Research. 2012: 617-621.

    [4] 劉永紅, 劉明雍, 謝波. 捷聯(lián)慣導(dǎo)系統(tǒng)雙位置快速抗干擾對(duì)準(zhǔn)方法[J]. 中國(guó)慣性技術(shù)學(xué)報(bào), 2014, 22(3): 296-300. Liu Yong-hong, Liu Ming-yong, Xie Bo. Rapid antiinterference method of two-position alignment of strapdown inertial navigation system[J]. Journal of Chinese Inertial Technology, 2014, 22(3): 296-300.

    [5] Silson P M G. Coarse alignment of a ship's strapdown inertial attitude reference system using velocity loci[J]. IEEE Transactions on Instrumentation and Measurement, 2011, 60(6): 1930-1941.

    [6] 嚴(yán)恭敏. 捷聯(lián)慣導(dǎo)系統(tǒng)動(dòng)基座初始對(duì)準(zhǔn)及其它相關(guān)問題研究[D]. 西安: 西北工業(yè)大學(xué), 2008. Yan Gong-min. On SINS in-movement initial alignment and some other problems[D]. Xi’an: Northwestern Polytechnical University, 2008.

    [7] Song L, Zhang C, Lu J. Self-alignment of full skewed RSINS: observability analysis and full-observable Kalman filter[J]. Journal of Systems Engineering and Electronics, 2014, 25(1): 104-114.

    [8] Ma L H, Wang K L, Li H. Gyrocompass alignment method of sins based on Kalman filtering pretreatment and dynamic gain adjustment on a rocking base[J]. Information Technology Journal, 2013, 12(4): 777-783.

    [9] 秦永元. 慣性導(dǎo)航[M]. 第二版. 北京: 科學(xué)出版社, 2014: 325-328.

    [10] Savage P G. Explicit frequency-shaped coning algorithms for pseudoconing environments[J]. Journal of Guidance, Control, and Dynamics, 2011, 34(3): 774-782.

    [11] Savage P G. Coning algorithm design by explicit frequency shaping[J]. Journal of Guidance, Control, and Dynamics, 2010, 33(4): 1123-1132.

    Anti-interference backtracking parameter identification alignment approach for strapdown inertial navigation system

    ZHANG Zhao-fei1,2, LUO Jian-jun1,2, HOU Yong-li3
    (1. School of Astronautics, Northwestern Polytechnical University, Xi’an 710072, China; 2. Science and Technology on Aerospace Flight Dynamics Laboratory, Xi’an 710072, China; 3. The 16th Institute, China Aerospace Science and Technology Corporation, Xi’an 710100, China)

    Since the parameter identification alignment approach of strapdown inertial navigation system has such problems as the convergent time is long and the azimuth error is affected by the initial value of recursive calculation, an identification alignment algorithm for backtracking parameters is proposed by the consistency of the formulae, in which the positive and negative direction sequences are obtained by storing and processing the alignment data of gyros and accelerometers. Based on the anti-interference coarse alignment in inertial reference frame, the negative and positive direction identification alignment is used repeatedly in the fine alignment algorithm of backtracking parameter identification alignment. Experiment results confirm that, by the new alignment algorithm, the azimuth error can be reduced by 75% without increasing alignment time.

    strapdown inertial navigation system; anti-interference; backtracking; parameter identification alignment

    V249.3

    :A

    2015-07-10;

    :2015-11-11

    國(guó)家自然科學(xué)基金(61004124)

    張朝飛(1987—),男,博士研究生,從事捷聯(lián)慣導(dǎo)研究。E-mail: fly060fly@qq.com

    聯(lián) 系 人:羅建軍(1965—),男,教授,博士生導(dǎo)師。E-mail: jjluo@nwpu.edu.cn

    1005-6734(2015)06-0733-06

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

    猜你喜歡
    捷聯(lián)慣導(dǎo)對(duì)準(zhǔn)
    自適應(yīng)模糊多環(huán)控制在慣導(dǎo)平臺(tái)穩(wěn)定回路中的應(yīng)用
    對(duì)準(zhǔn)提升組織力的聚焦點(diǎn)——陜西以組織振興引領(lǐng)鄉(xiāng)村振興
    無人機(jī)室內(nèi)視覺/慣導(dǎo)組合導(dǎo)航方法
    彈道導(dǎo)彈的捷聯(lián)慣性/天文組合導(dǎo)航方法
    一種改進(jìn)的速度加姿態(tài)匹配快速傳遞對(duì)準(zhǔn)算法
    基于Bagging模型的慣導(dǎo)系統(tǒng)誤差抑制方法
    捷聯(lián)慣性/天文/雷達(dá)高度表組合導(dǎo)航
    半捷聯(lián)雷達(dá)導(dǎo)引頭視線角速度提取
    INS/GPS組合系統(tǒng)初始滾轉(zhuǎn)角空中粗對(duì)準(zhǔn)方法
    基于多線程的慣導(dǎo)邏輯仿真器設(shè)計(jì)
    久久 成人 亚洲| av天堂在线播放| 在线观看免费午夜福利视频| 久久九九热精品免费| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品一区二区精品视频观看| 一级黄片播放器| 亚洲情色 制服丝袜| 亚洲情色 制服丝袜| 欧美人与性动交α欧美精品济南到| 老鸭窝网址在线观看| 国产三级黄色录像| 亚洲三区欧美一区| 美女福利国产在线| 亚洲成人手机| 美女主播在线视频| 丝袜美腿诱惑在线| 亚洲图色成人| 纯流量卡能插随身wifi吗| 国产深夜福利视频在线观看| 在现免费观看毛片| 一级毛片我不卡| 国产成人a∨麻豆精品| 一区在线观看完整版| 亚洲欧美成人综合另类久久久| 国产精品二区激情视频| 丁香六月天网| 一级黄色大片毛片| 一本—道久久a久久精品蜜桃钙片| 精品国产乱码久久久久久男人| 国产欧美日韩精品亚洲av| 一级a爱视频在线免费观看| 丝袜脚勾引网站| 亚洲国产最新在线播放| av片东京热男人的天堂| 亚洲天堂av无毛| 亚洲色图 男人天堂 中文字幕| 性色av乱码一区二区三区2| 久久天躁狠狠躁夜夜2o2o | 久久国产精品大桥未久av| 欧美激情高清一区二区三区| 黄片小视频在线播放| 亚洲人成网站在线观看播放| 国产亚洲av高清不卡| 亚洲色图综合在线观看| 99国产精品99久久久久| 久久九九热精品免费| 国产一区二区在线观看av| 亚洲精品日韩在线中文字幕| 大香蕉久久网| 91精品三级在线观看| av在线老鸭窝| 在线观看一区二区三区激情| 久久亚洲精品不卡| 久久久国产欧美日韩av| 人人妻,人人澡人人爽秒播 | 中文字幕亚洲精品专区| 精品人妻一区二区三区麻豆| 我要看黄色一级片免费的| 婷婷色综合www| 国产成人欧美在线观看 | 在线天堂中文资源库| 看免费成人av毛片| 狠狠精品人妻久久久久久综合| 精品福利观看| 丝袜在线中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品国产a三级三级三级| 国产又爽黄色视频| 一本—道久久a久久精品蜜桃钙片| 99香蕉大伊视频| 女警被强在线播放| 久久久久国产一级毛片高清牌| 国产一区二区在线观看av| 成年女人毛片免费观看观看9 | 美女脱内裤让男人舔精品视频| 性高湖久久久久久久久免费观看| 黑丝袜美女国产一区| 午夜免费男女啪啪视频观看| 精品少妇黑人巨大在线播放| 欧美人与善性xxx| 亚洲欧洲日产国产| 久久国产精品人妻蜜桃| 亚洲国产欧美在线一区| 精品国产乱码久久久久久小说| 久久久国产精品麻豆| 国产精品国产三级国产专区5o| 国产深夜福利视频在线观看| 欧美黑人精品巨大| 久久精品熟女亚洲av麻豆精品| 国产亚洲av片在线观看秒播厂| 成年动漫av网址| 中文字幕制服av| 欧美成人午夜精品| 999精品在线视频| 亚洲,一卡二卡三卡| 中文字幕高清在线视频| 成人国语在线视频| 成人国产av品久久久| 高潮久久久久久久久久久不卡| 成人亚洲欧美一区二区av| 少妇裸体淫交视频免费看高清 | 狂野欧美激情性bbbbbb| av一本久久久久| 精品熟女少妇八av免费久了| 午夜老司机福利片| 国产成人免费观看mmmm| 少妇的丰满在线观看| 99香蕉大伊视频| 亚洲伊人久久精品综合| 国产有黄有色有爽视频| 日韩av不卡免费在线播放| 美女大奶头黄色视频| 色网站视频免费| 国产精品九九99| 免费观看人在逋| 在线观看免费日韩欧美大片| 亚洲欧美精品自产自拍| 日韩人妻精品一区2区三区| 日韩人妻精品一区2区三区| 性少妇av在线| 亚洲国产精品一区三区| 美女福利国产在线| 下体分泌物呈黄色| 亚洲熟女毛片儿| 亚洲国产欧美日韩在线播放| 999久久久国产精品视频| 国产精品久久久人人做人人爽| 黄色视频不卡| 女人久久www免费人成看片| 欧美日韩综合久久久久久| 精品国产一区二区三区四区第35| 国产一区二区三区综合在线观看| 国产精品香港三级国产av潘金莲 | 色精品久久人妻99蜜桃| 亚洲色图 男人天堂 中文字幕| 久久性视频一级片| 亚洲一区二区三区欧美精品| 国产视频首页在线观看| 亚洲精品av麻豆狂野| 九色亚洲精品在线播放| 美女国产高潮福利片在线看| 91九色精品人成在线观看| 日韩一区二区三区影片| a级毛片黄视频| 免费在线观看影片大全网站 | 成年人免费黄色播放视频| 99精国产麻豆久久婷婷| 国产免费现黄频在线看| 色精品久久人妻99蜜桃| 伊人亚洲综合成人网| 超色免费av| 国产精品av久久久久免费| 国产真人三级小视频在线观看| 成年人免费黄色播放视频| 大香蕉久久网| 日韩大码丰满熟妇| 免费在线观看黄色视频的| 母亲3免费完整高清在线观看| 99精品久久久久人妻精品| 制服诱惑二区| 中国国产av一级| 在线av久久热| 人人妻人人澡人人看| 久久久久视频综合| 国产激情久久老熟女| 亚洲精品一二三| 日韩精品免费视频一区二区三区| 性高湖久久久久久久久免费观看| 国产精品一国产av| 日韩中文字幕视频在线看片| 亚洲,一卡二卡三卡| e午夜精品久久久久久久| 久久影院123| 侵犯人妻中文字幕一二三四区| 赤兔流量卡办理| 母亲3免费完整高清在线观看| 制服诱惑二区| 国产精品一区二区免费欧美 | 日韩中文字幕欧美一区二区 | 中文字幕最新亚洲高清| 水蜜桃什么品种好| 男女无遮挡免费网站观看| 免费女性裸体啪啪无遮挡网站| 日韩制服骚丝袜av| 国产免费视频播放在线视频| 亚洲精品一卡2卡三卡4卡5卡 | 99热网站在线观看| 校园人妻丝袜中文字幕| 青青草视频在线视频观看| 青春草亚洲视频在线观看| 多毛熟女@视频| netflix在线观看网站| 高清欧美精品videossex| 欧美日韩黄片免| 赤兔流量卡办理| 国产免费福利视频在线观看| 99久久精品国产亚洲精品| 国产成人欧美在线观看 | 亚洲精品成人av观看孕妇| 一二三四在线观看免费中文在| 国产高清视频在线播放一区 | 嫩草影视91久久| 91字幕亚洲| 精品一区二区三卡| 一个人免费看片子| 91字幕亚洲| 女性生殖器流出的白浆| 国产亚洲欧美在线一区二区| 人体艺术视频欧美日本| 少妇人妻久久综合中文| 观看av在线不卡| 黄色怎么调成土黄色| 国产老妇伦熟女老妇高清| 我的亚洲天堂| 手机成人av网站| 国产成人欧美| 深夜精品福利| 国产精品一国产av| 最黄视频免费看| 啦啦啦中文免费视频观看日本| 99热全是精品| 亚洲av综合色区一区| 美女午夜性视频免费| 久久久久国产精品人妻一区二区| 午夜激情久久久久久久| 亚洲综合色网址| 曰老女人黄片| a 毛片基地| 亚洲精品久久成人aⅴ小说| 99国产综合亚洲精品| 久久久亚洲精品成人影院| 性少妇av在线| 少妇人妻 视频| 9191精品国产免费久久| 高清欧美精品videossex| 午夜视频精品福利| 国产淫语在线视频| 亚洲一码二码三码区别大吗| 亚洲av男天堂| 大片电影免费在线观看免费| www日本在线高清视频| 极品少妇高潮喷水抽搐| 欧美精品亚洲一区二区| 老司机在亚洲福利影院| 青春草视频在线免费观看| 久久人妻熟女aⅴ| 欧美激情高清一区二区三区| 亚洲成国产人片在线观看| av在线播放精品| 2018国产大陆天天弄谢| 97人妻天天添夜夜摸| 中文字幕高清在线视频| 十八禁高潮呻吟视频| 成年人免费黄色播放视频| 欧美在线一区亚洲| 成人三级做爰电影| 久久青草综合色| 一级黄色大片毛片| 亚洲第一av免费看| 大型av网站在线播放| 桃花免费在线播放| 亚洲免费av在线视频| 亚洲精品日本国产第一区| 丝袜脚勾引网站| 我要看黄色一级片免费的| 久久久国产欧美日韩av| 欧美日本中文国产一区发布| 成人亚洲精品一区在线观看| 亚洲成国产人片在线观看| 丰满迷人的少妇在线观看| 侵犯人妻中文字幕一二三四区| 精品国产一区二区三区四区第35| 天天躁夜夜躁狠狠久久av| 制服人妻中文乱码| 啦啦啦啦在线视频资源| 一二三四社区在线视频社区8| 亚洲激情五月婷婷啪啪| 久久久精品94久久精品| 女性被躁到高潮视频| 天天躁夜夜躁狠狠躁躁| 好男人视频免费观看在线| 亚洲国产欧美日韩在线播放| 捣出白浆h1v1| 美女主播在线视频| 人人澡人人妻人| 国产欧美日韩综合在线一区二区| 狂野欧美激情性xxxx| 啦啦啦在线免费观看视频4| 国产精品免费视频内射| avwww免费| 成人手机av| 欧美日韩视频高清一区二区三区二| svipshipincom国产片| 在现免费观看毛片| 国产又爽黄色视频| 大香蕉久久网| 国产97色在线日韩免费| 久久青草综合色| 麻豆乱淫一区二区| 老熟女久久久| 欧美xxⅹ黑人| 久久这里只有精品19| 老司机影院成人| 青春草视频在线免费观看| 国产精品av久久久久免费| 自拍欧美九色日韩亚洲蝌蚪91| 麻豆乱淫一区二区| 国产日韩欧美亚洲二区| 在线看a的网站| 亚洲人成电影免费在线| 99九九在线精品视频| 王馨瑶露胸无遮挡在线观看| 久久精品熟女亚洲av麻豆精品| 只有这里有精品99| 首页视频小说图片口味搜索 | 精品亚洲成国产av| 国产男女超爽视频在线观看| 欧美中文综合在线视频| 久久精品国产综合久久久| 新久久久久国产一级毛片| 久久人妻福利社区极品人妻图片 | 午夜精品国产一区二区电影| 日本黄色日本黄色录像| 国产熟女午夜一区二区三区| 99久久精品国产亚洲精品| 亚洲国产欧美网| 三上悠亚av全集在线观看| 黄网站色视频无遮挡免费观看| 精品久久久久久久毛片微露脸 | 无遮挡黄片免费观看| 日韩视频在线欧美| 国产精品国产三级专区第一集| 亚洲人成电影免费在线| 欧美性长视频在线观看| 在线观看人妻少妇| 国产精品偷伦视频观看了| 国产成人av激情在线播放| 大片电影免费在线观看免费| 免费在线观看影片大全网站 | 美女大奶头黄色视频| 高清不卡的av网站| 久久久久网色| 巨乳人妻的诱惑在线观看| 日韩精品免费视频一区二区三区| 精品熟女少妇八av免费久了| 久久 成人 亚洲| 国产一区二区 视频在线| 亚洲 国产 在线| 九草在线视频观看| 国产精品一国产av| 超碰成人久久| 九色亚洲精品在线播放| 一区二区日韩欧美中文字幕| 丁香六月天网| 久久99精品国语久久久| www.熟女人妻精品国产| 亚洲少妇的诱惑av| 日韩av在线免费看完整版不卡| 2018国产大陆天天弄谢| 国产欧美日韩一区二区三区在线| 女性生殖器流出的白浆| 一区二区av电影网| 亚洲情色 制服丝袜| 在线亚洲精品国产二区图片欧美| 久久精品国产a三级三级三级| 天天操日日干夜夜撸| 黄色片一级片一级黄色片| 国产一区二区 视频在线| 男人爽女人下面视频在线观看| 婷婷色av中文字幕| 波野结衣二区三区在线| 欧美国产精品一级二级三级| 亚洲av综合色区一区| 欧美国产精品va在线观看不卡| av在线老鸭窝| 亚洲国产精品国产精品| 少妇裸体淫交视频免费看高清 | 日韩熟女老妇一区二区性免费视频| 麻豆av在线久日| 午夜日韩欧美国产| 午夜免费男女啪啪视频观看| 超色免费av| 中文字幕精品免费在线观看视频| 国产在线观看jvid| 国产成人精品久久二区二区91| 国产成人免费无遮挡视频| 国产在线一区二区三区精| 国产成人精品久久久久久| 亚洲精品国产色婷婷电影| 狂野欧美激情性bbbbbb| 久久中文字幕一级| 纵有疾风起免费观看全集完整版| 久久 成人 亚洲| 建设人人有责人人尽责人人享有的| 最近中文字幕2019免费版| 国产成人免费观看mmmm| 国产精品久久久久久精品电影小说| 亚洲av美国av| 国产免费福利视频在线观看| 九草在线视频观看| 亚洲色图 男人天堂 中文字幕| xxxhd国产人妻xxx| 亚洲人成网站在线观看播放| av网站免费在线观看视频| 欧美国产精品一级二级三级| 成年av动漫网址| 欧美精品一区二区大全| 91精品国产国语对白视频| 一个人免费看片子| 免费在线观看黄色视频的| 美女扒开内裤让男人捅视频| 女人高潮潮喷娇喘18禁视频| 男女免费视频国产| 久久热在线av| 久久久久久久大尺度免费视频| 日韩一本色道免费dvd| 亚洲国产欧美在线一区| 国产爽快片一区二区三区| 国产免费现黄频在线看| 成人影院久久| 久久精品熟女亚洲av麻豆精品| 亚洲国产日韩一区二区| av网站在线播放免费| 妹子高潮喷水视频| 精品一区二区三卡| 黄色片一级片一级黄色片| 欧美国产精品一级二级三级| 女人高潮潮喷娇喘18禁视频| 汤姆久久久久久久影院中文字幕| 高清欧美精品videossex| 成人免费观看视频高清| 亚洲人成电影免费在线| 一级黄色大片毛片| 99精国产麻豆久久婷婷| 亚洲国产中文字幕在线视频| av在线app专区| 电影成人av| 2018国产大陆天天弄谢| 无遮挡黄片免费观看| 亚洲av国产av综合av卡| 国产亚洲午夜精品一区二区久久| 一边摸一边做爽爽视频免费| 欧美精品一区二区大全| 大片免费播放器 马上看| 亚洲国产毛片av蜜桃av| 久久久久精品国产欧美久久久 | 精品欧美一区二区三区在线| 丝袜在线中文字幕| 亚洲精品日韩在线中文字幕| 国产成人a∨麻豆精品| 成年人免费黄色播放视频| 久久久久视频综合| 三上悠亚av全集在线观看| 国产福利在线免费观看视频| 一本—道久久a久久精品蜜桃钙片| 18禁观看日本| 纵有疾风起免费观看全集完整版| 无限看片的www在线观看| 另类亚洲欧美激情| bbb黄色大片| 婷婷成人精品国产| 人成视频在线观看免费观看| 一区二区三区乱码不卡18| 国产日韩欧美在线精品| 日韩制服骚丝袜av| 亚洲国产成人一精品久久久| 天堂中文最新版在线下载| 黄片播放在线免费| 成人午夜精彩视频在线观看| 少妇人妻久久综合中文| 欧美人与性动交α欧美软件| 午夜福利免费观看在线| 亚洲精品久久午夜乱码| 99re6热这里在线精品视频| 中文字幕最新亚洲高清| 欧美97在线视频| 中文字幕色久视频| av线在线观看网站| 男人操女人黄网站| 亚洲五月色婷婷综合| 男的添女的下面高潮视频| 老司机深夜福利视频在线观看 | 叶爱在线成人免费视频播放| 欧美中文综合在线视频| 搡老岳熟女国产| 人人妻人人爽人人添夜夜欢视频| 成人影院久久| 国产精品久久久人人做人人爽| 午夜久久久在线观看| 嫩草影视91久久| 视频区图区小说| 亚洲精品中文字幕在线视频| 欧美精品一区二区免费开放| 精品福利观看| 80岁老熟妇乱子伦牲交| 亚洲国产av影院在线观看| 欧美日韩亚洲高清精品| 永久免费av网站大全| 精品福利观看| 日韩精品免费视频一区二区三区| 国产男人的电影天堂91| 在线观看一区二区三区激情| 少妇的丰满在线观看| 亚洲国产精品国产精品| 最新的欧美精品一区二区| 在线 av 中文字幕| 久久人妻熟女aⅴ| 亚洲av成人精品一二三区| 一级黄色大片毛片| 在线观看免费视频网站a站| 亚洲欧洲日产国产| 狠狠精品人妻久久久久久综合| 亚洲黑人精品在线| 97人妻天天添夜夜摸| 少妇裸体淫交视频免费看高清 | av国产久精品久网站免费入址| 晚上一个人看的免费电影| 国产精品一国产av| 侵犯人妻中文字幕一二三四区| 最新的欧美精品一区二区| 国产黄色视频一区二区在线观看| 老汉色av国产亚洲站长工具| 亚洲,一卡二卡三卡| 性高湖久久久久久久久免费观看| 最新在线观看一区二区三区 | 亚洲伊人色综图| 一级毛片 在线播放| 91成人精品电影| 女警被强在线播放| av天堂在线播放| 午夜av观看不卡| 欧美精品啪啪一区二区三区 | 一本大道久久a久久精品| 美女视频免费永久观看网站| 日本av免费视频播放| 色播在线永久视频| 国产亚洲欧美在线一区二区| av天堂久久9| 18禁裸乳无遮挡动漫免费视频| 欧美日韩视频高清一区二区三区二| 日本黄色日本黄色录像| 国产av国产精品国产| 国产亚洲av片在线观看秒播厂| 又大又黄又爽视频免费| 五月天丁香电影| 久久精品国产亚洲av高清一级| 好男人视频免费观看在线| 国产av精品麻豆| 国产男人的电影天堂91| 男人舔女人的私密视频| 这个男人来自地球电影免费观看| 性色av一级| 一区在线观看完整版| 99国产精品99久久久久| 亚洲av日韩精品久久久久久密 | 精品人妻1区二区| 午夜老司机福利片| 波多野结衣av一区二区av| 中文字幕人妻熟女乱码| 国产视频一区二区在线看| 嫁个100分男人电影在线观看 | 婷婷色综合大香蕉| 热re99久久精品国产66热6| 99热国产这里只有精品6| 亚洲欧美精品综合一区二区三区| 国产黄频视频在线观看| 两个人看的免费小视频| 男女床上黄色一级片免费看| xxxhd国产人妻xxx| 丁香六月天网| 日本猛色少妇xxxxx猛交久久| 国产伦理片在线播放av一区| 精品久久蜜臀av无| 久久天堂一区二区三区四区| 成人影院久久| 美女中出高潮动态图| 免费久久久久久久精品成人欧美视频| 18禁国产床啪视频网站| 又大又黄又爽视频免费| tube8黄色片| 2018国产大陆天天弄谢| 国产高清视频在线播放一区 | 国产一区二区在线观看av| 成年人午夜在线观看视频| 18在线观看网站| 男人舔女人的私密视频| 色网站视频免费| 97人妻天天添夜夜摸| 如日韩欧美国产精品一区二区三区| 国产欧美日韩一区二区三区在线| 亚洲天堂av无毛| 丰满饥渴人妻一区二区三| 亚洲欧美成人综合另类久久久| 久久久久久久大尺度免费视频| 亚洲av男天堂| 亚洲 国产 在线| 亚洲精品乱久久久久久| 人妻 亚洲 视频| 亚洲伊人久久精品综合| 久久精品aⅴ一区二区三区四区| 国产成人啪精品午夜网站| 久久精品国产亚洲av涩爱| 成人国语在线视频| 久久狼人影院| 大码成人一级视频| 久久av网站| 久久久精品国产亚洲av高清涩受| 午夜影院在线不卡| 久久久久久久久免费视频了| 国产高清不卡午夜福利| 一级片'在线观看视频| 99re6热这里在线精品视频| 精品少妇久久久久久888优播| 少妇裸体淫交视频免费看高清 |