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

    基于偽地球坐標(biāo)系的慣導(dǎo)系統(tǒng)全球?qū)Ш剿惴?/h1>
    2017-11-17 02:13:39李光春高延濱
    中國慣性技術(shù)學(xué)報 2017年4期
    關(guān)鍵詞:系統(tǒng)

    劉 猛,李光春,高延濱,李 瑩,王 帝

    (哈爾濱工程大學(xué) 自動化學(xué)院,哈爾濱 150001)

    基于偽地球坐標(biāo)系的慣導(dǎo)系統(tǒng)全球?qū)Ш剿惴?/p>

    劉 猛,李光春,高延濱,李 瑩,王 帝

    (哈爾濱工程大學(xué) 自動化學(xué)院,哈爾濱 150001)

    由于極區(qū)特殊的地理、電磁條件,慣性導(dǎo)航因其自主性和信息完備性使之成為極區(qū)導(dǎo)航的首選。然而在考慮全球執(zhí)行能力時,現(xiàn)有常用的任何一種力學(xué)編排方案都不能單獨的實現(xiàn)全球?qū)Ш健Mǔ2捎媒M合編排的方式,這樣則不利于慣導(dǎo)算法的全球統(tǒng)一?;诖颂岢隽嘶趥蔚厍蜃鴺?biāo)系的全球?qū)Ш剿惴?。該方案在全球?qū)Ш綍r可以實現(xiàn)慣導(dǎo)算法的內(nèi)在統(tǒng)一,并可保證物理平臺的平穩(wěn)切換,從而實現(xiàn)平臺慣導(dǎo)與捷聯(lián)慣導(dǎo)系統(tǒng)編排方案的統(tǒng)一。另外,該方案也更方便同其它局部慣導(dǎo)系統(tǒng)進行交互通信,僅不同的參數(shù)轉(zhuǎn)換單元是必需的。同時,簡單的切換邏輯也可以減小程序設(shè)計的復(fù)雜度和降低計算機負擔(dān)。最后通過仿真證明了該算法的可行性。

    偽地球坐標(biāo)系;全球?qū)Ш?;極區(qū)導(dǎo)航;捷聯(lián)慣導(dǎo)系統(tǒng);平臺慣導(dǎo)系統(tǒng)

    由于極區(qū)特殊的地理、電磁條件,慣性導(dǎo)航因其自主性和信息完備性使之成為極區(qū)導(dǎo)航的首選[1]。然而由于極區(qū)經(jīng)線快速收斂,且最終匯聚于一點,使建立相對于經(jīng)線的航向方向越來越困難,從而造成慣導(dǎo)系統(tǒng)方位陀螺施矩困難、計算溢出和誤差放大等問題,故傳統(tǒng)指北方位系統(tǒng)的慣導(dǎo)算法在極區(qū)是失效的[2]。

    為了克服上述困難,在極地地區(qū)通常采用游移方位慣導(dǎo)力學(xué)編排的方式。然而由于游移角的定義在地理極點是無意義的,游移方位系統(tǒng)在地理極點附近也無法正常工作[3]。因此,格網(wǎng)坐標(biāo)系和橫向地球坐標(biāo)系被提出用于解決極區(qū)導(dǎo)航問題。針對格網(wǎng)導(dǎo)航,文獻[4]對其進行了詳細研究,將導(dǎo)航參數(shù)中的姿態(tài)和速度信息在格網(wǎng)坐標(biāo)系下表示,位置信息在ECEF坐標(biāo)系中給出,而橫向地球坐標(biāo)系統(tǒng)的定義更相似于傳統(tǒng)地球坐標(biāo)系,且其慣導(dǎo)編排也更直觀和方便,從而可能使橫向慣導(dǎo)系統(tǒng)更受青睞。文獻[5]對橫向慣導(dǎo)力學(xué)編排進行了詳細的分析;文獻[6]和[7]對橫向?qū)Ш剿惴ǖ男U妥枘徇M行設(shè)計來提高其長期導(dǎo)航精度和抑制其周期性振蕩;文獻[8]和[9]在橢球地球模型下對橫向慣導(dǎo)算法進行了優(yōu)化設(shè)計,結(jié)果表明橫向慣導(dǎo)系統(tǒng)可以解決極區(qū)導(dǎo)航問題。然而由于格網(wǎng)慣導(dǎo)和橫向慣導(dǎo)力學(xué)編排未考慮不同導(dǎo)航坐標(biāo)系間的物理切換[10-11],格網(wǎng)坐標(biāo)系編排和橫向地球坐標(biāo)系編排通常僅適用于捷聯(lián)慣導(dǎo)系統(tǒng)。具體來說,當(dāng)導(dǎo)航坐標(biāo)系由當(dāng)?shù)氐乩碜鴺?biāo)系切換到橫向地理坐標(biāo)系或格網(wǎng)坐標(biāo)系時,由于導(dǎo)航坐標(biāo)系不是連續(xù)性變化的,則需以突變的形式進行切換。大階躍響應(yīng)式切換可能會造成物理平臺追蹤切換的施矩困難,從而影響導(dǎo)航精度或直接造成導(dǎo)航系統(tǒng)的發(fā)散。

    另一方面,在考慮全球執(zhí)行能力時,現(xiàn)有常用的任何一種力學(xué)編排都不能單獨實現(xiàn)全球?qū)Ш?,因此通常采用指?格網(wǎng)編排或橫向編排的組合編排方式來實現(xiàn)其全球?qū)Ш降哪芰ΑA硗?,對于目前常用的編排方案,每一種編排方式都是一套自身完備的導(dǎo)航方案,新的力學(xué)編排的引入會不利于慣導(dǎo)算法的全球統(tǒng)一,而且由于其復(fù)雜的切換邏輯,也會增加程序設(shè)計的復(fù)雜度和計算機負擔(dān),因此文獻[10]和[11]提出了以游移編排為核心的間接格網(wǎng)慣導(dǎo)算法和間接橫向慣導(dǎo)算法來實現(xiàn)其全球?qū)Ш降哪芰?。然而游移方位角的引入增加了慣導(dǎo)編排的復(fù)雜度,使其物理意義不明確,而且需要計算游移方位角。

    針對上述問題,本文采用偽地球坐標(biāo)系慣導(dǎo)力學(xué)編排,提出了基于偽地球坐標(biāo)系的慣導(dǎo)系統(tǒng)全球?qū)Ш剿惴?,將整個慣導(dǎo)系統(tǒng)的實現(xiàn)分成“核心部分”和“轉(zhuǎn)換單元”兩部分完成?!昂诵牟糠帧币砸环N簡單統(tǒng)一的導(dǎo)航編排方案來實現(xiàn)慣導(dǎo)系統(tǒng)全球范圍內(nèi)的導(dǎo)航計算,并可以保證物理平臺追蹤切換的平穩(wěn)性,使其不僅適用于捷聯(lián)慣導(dǎo)系統(tǒng),而且也適用于平臺慣導(dǎo)系統(tǒng),同時實現(xiàn)了慣導(dǎo)算法的全球范圍內(nèi)在統(tǒng)一?!稗D(zhuǎn)換單元”則用于修正和轉(zhuǎn)換導(dǎo)航參數(shù),從而得到期望的導(dǎo)航參數(shù),用于人機交互,以及實現(xiàn)與其它不同的局部慣導(dǎo)系統(tǒng)的交互通信。最后通過仿真驗證了該算法的可行性。

    1 偽地球坐標(biāo)系定義

    偽地球坐標(biāo)系是根據(jù)導(dǎo)航坐標(biāo)系切換位置來進行構(gòu)建,其本質(zhì)是一種廣義的地球坐標(biāo)系。假設(shè)導(dǎo)航坐標(biāo)系切換的當(dāng)?shù)亟?jīng)度和地心緯度分別為則偽地球坐標(biāo)系可以通過原地球坐標(biāo)系旋轉(zhuǎn)得到。如圖1所示,在偽地球坐標(biāo)系中,其坐標(biāo)原點仍在地球中心,Np是偽北極點,Oxp軸與原Oze重合,軸是指向載體初始位置在赤道平面的投影點,Ozp在赤道平面內(nèi)垂直于Oyp軸,且偽赤道與當(dāng)?shù)匚恢盟诘淖游缇€圈是重合的,其具體定義詳見文獻[12]。則原地球坐標(biāo)系和偽地球坐標(biāo)系之間的旋轉(zhuǎn)矩陣為

    圖1 偽地球坐標(biāo)系定義Fig.1 Definition of pseudo-Earth frame

    在偽地球坐標(biāo)系統(tǒng)中,載體位置用偽經(jīng)緯度和偽高度來表示,即λp、Lp、hp。偽地理坐標(biāo)系的定義與傳統(tǒng)的地心地理坐標(biāo)系( 地心地理坐標(biāo)系是根據(jù)當(dāng)?shù)氐匦慕?jīng)緯度構(gòu)建的)是相似的[12],如圖1所示。則在導(dǎo)航坐標(biāo)系切換位置,載體的偽位置與當(dāng)?shù)氐匦奈恢靡约皞蔚乩碜鴺?biāo)系與地心地理坐標(biāo)系存在如下關(guān)系[12]:

    式中,h0為在切換位置的載體的高度。

    由公式(2)可知,在導(dǎo)航坐標(biāo)系切換后,初始偽緯度總為 0。則在偽地球坐標(biāo)系中,載體將總是從偽赤道附近開始航行,這樣導(dǎo)航坐標(biāo)系切換后在偽地球坐標(biāo)系進行導(dǎo)航時,可以長時間保持在中低偽緯度區(qū)域,從而避免導(dǎo)航坐標(biāo)系的頻繁切換。另外,由公式(3)可知,在切換位置,偽地理坐標(biāo)系(即偽東-北-天)與傳統(tǒng)南-東-天地心地理坐標(biāo)系是一致的,并且其切換矩陣是一個已知固定的常值矩陣,因此在偽地球坐標(biāo)系下進行力學(xué)編排,偽地理坐標(biāo)系選為導(dǎo)航系,則當(dāng)導(dǎo)航坐標(biāo)系從地心地理坐標(biāo)系切換到偽地理坐標(biāo)系或從之前一個偽地理坐標(biāo)系切換到新建的偽地理坐標(biāo)系時,導(dǎo)航坐標(biāo)系是連續(xù)變化的,從而可以保證物理平臺的平穩(wěn)切換。另外,已知固定的常值切換矩陣,進行簡單的切換邏輯對于捷聯(lián)慣導(dǎo)系統(tǒng)也可以減小程序設(shè)計的復(fù)雜度和降低計算機負擔(dān)。

    2 基于偽地球坐標(biāo)系的全球?qū)Ш剿惴?/h2>

    2.1 基于偽地球坐標(biāo)系的全球?qū)Ш椒桨?/h3>

    為了有效解決慣導(dǎo)系統(tǒng)全球?qū)Ш降膯栴},應(yīng)探索出一種全球?qū)Ш椒桨福瑢θ蚍秶鷥?nèi)導(dǎo)航算法進行統(tǒng)一,使它不僅適用于捷聯(lián)慣導(dǎo)系統(tǒng),而且也適用于平臺慣導(dǎo)系統(tǒng),同時也方便與外部其它局部導(dǎo)航系統(tǒng)進行交互通信?;诖耍疚奶岢隽艘环N基于偽地球坐標(biāo)系的全球?qū)Ш椒桨?,以“躲進小樓成一統(tǒng),管他冬夏與春秋?!钡闹髦妓枷雭韺崿F(xiàn)慣導(dǎo)系統(tǒng)的全球?qū)Ш?,將整個慣導(dǎo)系統(tǒng)的實現(xiàn)分成“核心部分”和“轉(zhuǎn)換單元”兩部分完成,從而保證慣導(dǎo)算法全球范圍內(nèi)的內(nèi)在統(tǒng)一。

    圖2 基于偽地球坐標(biāo)系慣導(dǎo)系統(tǒng)的全球?qū)Ш椒桨傅姆娇驁DFig.2 Scheme of global navigation for inertial navigation system based on pseudo-Earth frame

    如圖2所示,根據(jù)載體活動范圍,選擇合適的切換位置構(gòu)建偽地球坐標(biāo)系,在合適的時刻完成導(dǎo)航坐標(biāo)系的切換,從而可以直接避免在導(dǎo)航過程中出現(xiàn)極區(qū)和偽極區(qū)導(dǎo)航問題。慣導(dǎo)系統(tǒng)核心部分是在球形地球模型下,根據(jù)偽地球坐標(biāo)系力學(xué)編排來完成導(dǎo)航解算和平臺追蹤的。根據(jù)公式(3),對于平臺慣導(dǎo)系統(tǒng),通過切換追蹤軸來實現(xiàn)平臺追蹤的平穩(wěn)切換,對于捷聯(lián)慣導(dǎo)系統(tǒng),以固定的切換邏輯和簡單的變換方式來實現(xiàn)不同導(dǎo)航坐標(biāo)系下導(dǎo)航參數(shù)的切換。然后根據(jù)計算出的偽地球坐標(biāo)系下的導(dǎo)航參數(shù),利用參數(shù)轉(zhuǎn)換單元對導(dǎo)航參數(shù)進行修正和轉(zhuǎn)換,從而得到期望的導(dǎo)航參數(shù)并保證導(dǎo)航精度,進而用于人機交互以及實現(xiàn)與其它局部慣導(dǎo)系統(tǒng)的交互通信,例如,局部指北向力學(xué)編排慣導(dǎo)系統(tǒng)、局部格網(wǎng)慣導(dǎo)系統(tǒng)、局部橫向慣導(dǎo)系統(tǒng)等。因此,該導(dǎo)航方案在全球?qū)Ш綍r可以實現(xiàn)慣導(dǎo)算法的內(nèi)在統(tǒng)一,并適用于平臺慣導(dǎo)與捷聯(lián)慣導(dǎo)系統(tǒng),從而實現(xiàn)慣導(dǎo)編排的統(tǒng)一,同時也是非常方便同其它局部慣導(dǎo)系統(tǒng)進行交互通信。

    對于切換位置的選擇可以分為以下兩種情況:

    1)載體運動范圍固定且運動區(qū)域不太大(即以載體初始位置為中心,活動范圍大小與中低緯度區(qū)域相等)。可以根據(jù)載體的初始位置來直接構(gòu)建偽地球坐標(biāo)系,然后通過跟蹤偽地理坐標(biāo)系來實現(xiàn)全球?qū)Ш浇馑悖瑥亩梢员苊庠趯?dǎo)航過程中出現(xiàn)導(dǎo)航坐標(biāo)系的切換。

    2)載體運動區(qū)域較大(例如,極限區(qū)域,即在整個全球范圍內(nèi))。當(dāng)僅采用根據(jù)初始位置構(gòu)建的偽地球坐標(biāo)系時,同樣會面臨著偽極區(qū)導(dǎo)航的問題,因此在載體運動過程中應(yīng)選擇合適的位置構(gòu)建偽地球坐標(biāo)系進行導(dǎo)航坐標(biāo)系的切換,從而完成慣導(dǎo)系統(tǒng)的全球?qū)Ш?,并可避免出現(xiàn)極區(qū)導(dǎo)航問題。

    綜上所述,本文給出的基于偽地球坐標(biāo)系的全球?qū)Ш椒桨冈O(shè)計具有以下優(yōu)點:

    1)基于載體活動范圍的偽地球坐標(biāo)系構(gòu)建規(guī)則,可以最大程度避免導(dǎo)航坐標(biāo)系的頻繁切換,進而可避免由導(dǎo)航坐標(biāo)系切換帶來的可能潛在問題。

    2)以球形地球模型下偽慣導(dǎo)力學(xué)編排作為慣導(dǎo)系統(tǒng)的核心,用一種簡單的編排方案來實現(xiàn)慣導(dǎo)系統(tǒng)全球?qū)Ш剿惴ǖ膬?nèi)在統(tǒng)一,并可保證物理平臺的平穩(wěn)切換,使其不僅適用于捷聯(lián)慣導(dǎo)系統(tǒng),而且也適用于平臺慣導(dǎo)系統(tǒng)。球形地球模型的導(dǎo)航計算還可以避免復(fù)雜的曲率半徑計算[8],并有利于不同導(dǎo)航坐標(biāo)系之間的導(dǎo)航參數(shù)轉(zhuǎn)換。另外,偽地球坐標(biāo)系統(tǒng)的構(gòu)建與傳統(tǒng)地球坐標(biāo)系以及橫向坐標(biāo)系的構(gòu)建非常相似,從而使定義簡單易懂,也更有利于該全球?qū)Ш椒桨傅膶嵤?/p>

    3)轉(zhuǎn)換單元可以保證慣導(dǎo)系統(tǒng)的導(dǎo)航精度并獲得期望的導(dǎo)航參數(shù),從而使其能更方便地與外部其它局部導(dǎo)航系統(tǒng)進行交互通信。當(dāng)需要不同的導(dǎo)航參數(shù)時,僅僅不同的參數(shù)轉(zhuǎn)換單元是必要的,慣導(dǎo)系統(tǒng)的核心部分則無需改變,從而以“躲進小樓成一統(tǒng),管他冬夏與春秋?!钡姆绞綄崿F(xiàn)了捷聯(lián)慣導(dǎo)與平臺慣導(dǎo)算法的全球范圍內(nèi)在統(tǒng)一,并方便進行人機交互和通信。

    2.2 偽地理坐標(biāo)系和當(dāng)?shù)氐乩碜鴺?biāo)系之間的關(guān)系

    設(shè)S為地球上一點,S點的地心緯度和地理緯度分別為Lc和L,則地理緯度與地心緯度的偏差可以表示為[14]

    其中,f為地球旋轉(zhuǎn)橢球體模型的扁率,f=1/298.3。

    則以地心緯度構(gòu)建的地心地理坐標(biāo)系(T系)與以地理緯度構(gòu)建的當(dāng)?shù)氐乩碜鴺?biāo)系(t系)在x軸上存在一個旋轉(zhuǎn)角σ,并且σ=ΔL,從而T系和t系之間的關(guān)系可以表示為[14]

    進而根據(jù)坐標(biāo)系間的變換關(guān)系,偽地理坐標(biāo)系和地理坐標(biāo)系之間的變換矩陣可以表示為

    式中,

    式(8)可以根據(jù)載體偽位置坐標(biāo)的偽經(jīng)緯度求得。假設(shè)載體的位置為S,且其離地高度為h,則根據(jù)投影關(guān)系容易得到載體相對于偽地球坐標(biāo)系的位置矢量pr,它與相對于地球坐標(biāo)系的位置矢量er存在如下關(guān)系:

    其中,

    R為地球半徑。

    由式(1)(9)則容易得到載體當(dāng)?shù)氐匦慕?jīng)緯度為

    其中,

    根據(jù)式(10)可知,式(8)是容易求得的,進而也可以得到偽地理坐標(biāo)系和地理坐標(biāo)系之間的變換矩陣。

    2.3 球形地球模型下的偽慣導(dǎo)力學(xué)編排

    偽地球坐標(biāo)系慣導(dǎo)力學(xué)編排與傳統(tǒng)指北方位和橫向慣導(dǎo)力學(xué)編排是相似的。在偽地球坐標(biāo)系中,偽地理坐標(biāo)系被選為導(dǎo)航系,因此偽地理坐標(biāo)系的旋轉(zhuǎn)角速度將發(fā)生變化。偽地理坐標(biāo)系的旋轉(zhuǎn)角速度是由地球自轉(zhuǎn)引起的旋轉(zhuǎn)角速度和載體運動引起的相對地球的旋轉(zhuǎn)角速度兩部分組成。根據(jù)定義,由地球自轉(zhuǎn)而引起的偽地理坐標(biāo)系旋轉(zhuǎn)角速度可以表示為

    根據(jù)公式(1)和(7),由地球自轉(zhuǎn)而引起的偽地理坐標(biāo)系的旋轉(zhuǎn)角速度可以很容易得到。另外,由載體運動引起的偽地理坐標(biāo)系相對于偽地球坐標(biāo)系的角速度與傳統(tǒng)指北慣導(dǎo)系統(tǒng)在表達形式上是一樣的,可以表示為

    根據(jù)式(11)(12),容易得到慣導(dǎo)系統(tǒng)控制平臺的旋轉(zhuǎn)角速度,可以表示為

    則對于平臺慣導(dǎo)系統(tǒng),其平臺指令角速度為

    而對于捷聯(lián)慣導(dǎo)系統(tǒng),其姿態(tài)更新矩陣為

    另外,我們還可以推導(dǎo)出如下慣導(dǎo)系統(tǒng)的比力方程和位置方程:

    2.4 轉(zhuǎn)換單元

    根據(jù)上述的慣導(dǎo)力學(xué)編排方案,偽慣導(dǎo)系統(tǒng)導(dǎo)航計算采用的是理想的球形地球模型。因此為了保證導(dǎo)航精度并獲得期望的導(dǎo)航參數(shù),必須有一個導(dǎo)航參數(shù)轉(zhuǎn)換單元對偽地球坐標(biāo)系下導(dǎo)航參數(shù)進行修正和轉(zhuǎn)換。由于受到篇幅限制,并且常規(guī)地球坐標(biāo)系的導(dǎo)航信息是人所熟知和默認的,在此我們僅給出用于得到傳統(tǒng)橢球模型下當(dāng)?shù)氐乩碜鴺?biāo)系導(dǎo)航參數(shù)的修正與轉(zhuǎn)換。

    另外,根據(jù)公式(6)還可以求得當(dāng)?shù)氐乩碜鴺?biāo)系下的姿態(tài)矩陣及載體速度:

    根據(jù)上述導(dǎo)航參數(shù)的轉(zhuǎn)換關(guān)系,可以通過轉(zhuǎn)換單元對偽地球坐標(biāo)系的導(dǎo)航參數(shù)進行轉(zhuǎn)換,得到傳統(tǒng)橢球模型下當(dāng)?shù)氐乩碜鴺?biāo)系導(dǎo)航參數(shù)且無精度犧牲,從而用于人機交互和通信等。

    綜上所述,基于偽地球坐標(biāo)系的慣導(dǎo)系統(tǒng)全球?qū)Ш竭^程可以表示為圖3所示的流程圖。另外,值得注意的是,T系也可以認為是一個特殊的偽地球地理坐標(biāo)系,此時矩陣是一個單位矩陣。

    圖3 基于偽地球坐標(biāo)系的慣導(dǎo)系統(tǒng)全球?qū)Ш搅鞒虉DFig.3 Flow chart of global navigation for INS based on pseudo-Earth frame

    3 誤差分析

    根據(jù)2.3節(jié)偽地球坐標(biāo)系慣導(dǎo)系統(tǒng)的導(dǎo)航方程,偽慣導(dǎo)編排下的靜態(tài)誤差方程很容易得到。忽略重力誤差并略去二階小量,則小失準(zhǔn)角偽導(dǎo)航靜態(tài)誤差方程為:

    由式(23)~(31)和文獻[5]可知,偽慣導(dǎo)靜態(tài)誤差方程和橫向慣導(dǎo)的靜態(tài)誤差方程在數(shù)學(xué)形式上是一致的,因此他們具有相同的傳播特性。根據(jù)文獻[5],偽地球坐標(biāo)系慣導(dǎo)系統(tǒng)同樣存在著舒勒、傅科以及地球三種周期振蕩,并且其偽緯度誤差和偽經(jīng)度誤差是隨時間發(fā)散的。另外由式(2)可知,在構(gòu)建偽地球坐標(biāo)系的切換位置時,偽緯度Lp總是為0,因此偽方位誤差方程式(22)可以重寫為:

    由式(32)偽方位誤差方程可知,在切換位置,線性運動誤差與方位誤差之間實現(xiàn)了解耦。因此如果依據(jù)載體初始位置來構(gòu)建偽地球坐標(biāo)系進行初始對準(zhǔn),則可以消除導(dǎo)航坐標(biāo)系旋轉(zhuǎn)角速度誤差(即或在傳統(tǒng)編排下為對方位對準(zhǔn)的影響。因此基于偽地球坐標(biāo)系的初始對準(zhǔn)也具有優(yōu)良性質(zhì)[12]。另外,由公式(3)可知,偽地理坐標(biāo)系(即偽東-北-天)與傳統(tǒng)南-東-天地心地理坐標(biāo)系是一致的,因此,即使由傳統(tǒng)的初始對準(zhǔn)算法完成慣導(dǎo)系統(tǒng)的初始對準(zhǔn)也無需進行重新設(shè)計。其次,由于偽地球坐標(biāo)系是根據(jù)載體位置構(gòu)建的,因此它不受載體初始位置的影響,也可用于極區(qū)初始對準(zhǔn),此結(jié)論已在文獻[12]和[13]給出了證明。因此基于偽地球坐標(biāo)的全球?qū)Ш椒桨笇τ趯崿F(xiàn)慣導(dǎo)系統(tǒng)的全球?qū)Ш绞强尚星液侠淼摹?/p>

    4 仿真與分析

    為了驗證基于偽地球坐標(biāo)系全球?qū)Ш剿惴ǖ目尚行院秃侠硇裕竟?jié)通過低緯度導(dǎo)航和極區(qū)導(dǎo)航仿真來對其進行驗證。仿真參數(shù)設(shè)置如下:陀螺常值漂移為0.01 (°)/h,隨機游走系數(shù)為 0.001 (°)/h1/2;加速度計常值零偏為100 μg,隨機游走系數(shù)為10 μg/Hz1/2;載體的真實姿態(tài)角θ=γ=ψ=0°;導(dǎo)航系統(tǒng)的初始姿態(tài)誤差參數(shù)為初始速度誤差為為了防止慣導(dǎo)系統(tǒng)高度通道發(fā)散,在導(dǎo)航解算過程中假設(shè)載體天向速度為零。

    4.1 低緯度導(dǎo)航仿真

    載體初始位置L0=15.6778°,λ0=126.6778°;載體以20 m/s的航速沿著載體艏向方向航行3.7 h。根據(jù)本文基于偽地球坐標(biāo)系的全球?qū)Ш剿惴ǎ瑧T導(dǎo)系統(tǒng)的核心部分首先在南-東-天地心地理坐標(biāo)系(T系)進行導(dǎo)航計算,在0.5 h后通過簡單的變換切換到偽地理坐標(biāo)系(P系)進行導(dǎo)航解算。在此過程中利用傳統(tǒng)地理坐標(biāo)系導(dǎo)航參數(shù)轉(zhuǎn)換單元來修正和轉(zhuǎn)換導(dǎo)航參數(shù)。圖4給出了傳統(tǒng)指北編排導(dǎo)航算法和本文全球?qū)Ш剿惴ㄏ碌恼`差曲線。圖5給出了傳統(tǒng)慣導(dǎo)系統(tǒng)和本文慣導(dǎo)系統(tǒng)核心部分計算出的載體姿態(tài)。圖4和圖5中,虛線代表傳統(tǒng)指北編排算法,實線代表基于偽地球坐標(biāo)系的全球?qū)Ш剿惴ā?/p>

    分析圖4所示的導(dǎo)航誤差曲線可知:傳統(tǒng)慣導(dǎo)算法與本文的全球慣導(dǎo)算法具有幾乎相同的導(dǎo)航精度。在3.7 h時刻,利用兩種導(dǎo)航算法計算的姿態(tài)和速度之間的差別分別是 0.034′、0.065′、0.006′和 0.105 m/s、0.003 m/s。對于位置誤差,基于偽地球坐標(biāo)系全球?qū)Ш剿惴ǖ慕?jīng)度誤差甚至更小。另外即使在0.5 h時刻進行了一次導(dǎo)航坐標(biāo)系的切換,本文全球?qū)Ш剿惴ǖ膶?dǎo)航誤差仍是平穩(wěn)周期性振蕩,沒有任何跳躍和突變。因此本文基于偽地球坐標(biāo)系的全球?qū)Ш剿惴ㄊ强尚泻秃侠淼?,并實現(xiàn)了全球?qū)Ш剿惴ǖ膬?nèi)在統(tǒng)一。

    圖4 傳統(tǒng)指北編排和偽地球坐標(biāo)系編排下導(dǎo)航誤差曲線Fig.4 Error curves of navigation based on traditional north-oriented and pseudo-Earth mechanizations

    圖5 傳統(tǒng)指北編排和偽地球坐標(biāo)系慣導(dǎo)核心的姿態(tài)解算Fig.5 Vehicle attitude with traditional north-oriented mechanization and the core of INS based on pseudo-Earth frame

    其次,分析圖5可知,在導(dǎo)航坐標(biāo)系切換的時刻,偽地理坐標(biāo)系與傳統(tǒng)的地心地理坐標(biāo)系(南-東-天)是一致的,而切換后慣導(dǎo)系統(tǒng)計算的載體姿態(tài)也是連續(xù)變化的。因此基于偽地球坐標(biāo)系導(dǎo)航編排可以實現(xiàn)導(dǎo)航坐標(biāo)系的連續(xù)變化,從而保證物理平臺導(dǎo)航追蹤的平穩(wěn)切換。因此基于偽地球坐標(biāo)系的全球?qū)Ш剿惴梢詫崿F(xiàn)捷聯(lián)慣導(dǎo)與平臺慣導(dǎo)系統(tǒng)的編排方案統(tǒng)一。

    4.2 極區(qū)導(dǎo)航仿真

    設(shè)載體初始位置為L0=87°、λ0=126°,載體以30 m/s 航速沿著載體艏向方向航行4.7 h,則根據(jù)載體的初始位置構(gòu)建偽地球坐標(biāo)系進行導(dǎo)航計算,得出的導(dǎo)航誤差曲線如圖6所示。

    由圖6可知,在極區(qū)使用偽地球坐標(biāo)系力學(xué)編排來進行導(dǎo)航解算,導(dǎo)航誤差也是平滑、連續(xù)地進行振蕩,無跳變現(xiàn)象發(fā)生,因此可以有效解決極區(qū)導(dǎo)航的問題,并滿足極區(qū)導(dǎo)航精度的要求。同時根據(jù)基于偽地球坐標(biāo)系的全球?qū)Ш椒桨福搨蔚厍蜃鴺?biāo)系是根據(jù)載體的初始位置來構(gòu)建的,從而可以最大程度地減小導(dǎo)航坐標(biāo)系的頻繁切換,也可以避免在切換中出現(xiàn)可能的潛在問題。因此本文基于偽地球坐標(biāo)系的全球?qū)Ш椒桨甘强尚星矣幸饬x的,并具有優(yōu)良的性能。

    圖6 偽地球坐標(biāo)力學(xué)編排下的導(dǎo)航誤差Fig.6 Errors of navigation with pseudo-Earth mechanization

    5 結(jié) 論

    為了有效地解決慣導(dǎo)系統(tǒng)全球?qū)Ш剑òO區(qū)導(dǎo)航)的問題,采用了偽地球坐標(biāo)系的慣導(dǎo)力學(xué)編排,從全球?qū)Ш浇嵌瘸霭l(fā)提出了基于偽地球坐標(biāo)系的全球?qū)Ш剿惴?。以“躲進小樓成一統(tǒng),管他冬夏與春秋?!钡闹髦妓枷?,將整個慣導(dǎo)系統(tǒng)的實現(xiàn)分成“核心部分”和“轉(zhuǎn)換單元”兩部分完成?!昂诵牟糠帧币郧蛐蔚厍蚰P拖聜螒T導(dǎo)力學(xué)編排來完成慣導(dǎo)系統(tǒng)在全球范圍內(nèi)的導(dǎo)航計算,從而可以保證物理平臺追蹤切換的平穩(wěn)性,使其不僅適用于捷聯(lián)慣導(dǎo)系統(tǒng)而且也適用于平臺慣導(dǎo)系統(tǒng)。同時“轉(zhuǎn)換單元”用于修正和轉(zhuǎn)換導(dǎo)航參數(shù),從而得到期望的導(dǎo)航參數(shù),而且也是更方便同其它局部慣導(dǎo)系統(tǒng)進行交互通信。當(dāng)需要不同導(dǎo)航系統(tǒng)的導(dǎo)航參數(shù)時,僅不同的參數(shù)轉(zhuǎn)換單元是必要的,慣導(dǎo)系統(tǒng)的核心單元將無需改變。

    最后通過仿真對其進行了驗證,結(jié)果表明,該導(dǎo)航算法是可行和合理的,同時可以解決極區(qū)導(dǎo)航的問題,并對全球?qū)Ш剿惴▽崿F(xiàn)了全球范圍內(nèi)在統(tǒng)一。

    (References):

    [1]Dyer G C. Polar navigation - A new transverse mercator technique[J]. Journal of Navigation, 1971, 24(4): 484-495.

    [2]Salychev O S. Applied inertial navigation: Problems and solutions. Moscow: BMSTU Press, 2004.

    [3]Cheng J, Wang T, Guan D, et al. Polar transfer alignment of shipborne SINS with a large misalignment angle[J]. Measurement Science & Technology, 2016, 27(3): 035101.

    [4]周琪, 秦永元, 付強文, 等. 極區(qū)飛行格網(wǎng)慣性導(dǎo)航算法原理[J]. 西北工業(yè)大學(xué)學(xué)報, 2013, 31(2): 210-217.Zhou Q, Qin Y Y, Fu Q W, et al. Grid mechanization in inertial navigation systems for transpolar aircraft[J].Journal of Northwestern Polytechnical University, 2013,31(2): 210-217.

    [5]李倩, 孫楓, 奔粵陽, 等. 基于橫坐標(biāo)系的捷聯(lián)慣導(dǎo)系統(tǒng)極區(qū)導(dǎo)航方法[J]. 中國慣性技術(shù)學(xué)報, 2014, 22(3):288-295.Li Q, Sun F, Ben Y Y, et al. Polar navigation of strapdown inertial navigation system based on transversal frame in polar region[J]. Journal of Chinese Inertial Technology,2014, 22(3): 288-295.

    [6]Li Q, Ben Y Y, Yu F, et al. System reset of transversal strapdown INS for ship in polar region[J]. Measurement,2015, 60: 247-257.

    [7]Li Q, Ben Y Y, Sun F, et al. Transversal strapdown INS and damping technology for marine in polar region[C]//IEEE/ION Position, Location & Navigation Symposium.2014: 1365-1370.

    [8]Li Q, Ben Y Y, Yu F, et al. Transversal strapdown INS based on reference ellipsoid for vehicle in polar region[J].IEEE Transactions on Vehicular Technology, 2016, 65(9):7791-7795.

    [9]Yao Y. Transverse navigation under the ellipsoidal earth model and its performance in both polar and non-polar areas[J]. Journal of Navigation, 2016, 69(2): 335-352.

    [10]周琪, 岳亞洲, 張曉冬, 等. 極區(qū)飛行間接格網(wǎng)慣性導(dǎo)航算法[J]. 中國慣性技術(shù)學(xué)報, 2014, 22(1): 18-17.Zhou Q, Yue Y Z, Zhang X D, et al. Indirect grid inertial navigation mechanization for transpolar aircraft[J]. Journal of Chinese Inertial Technology, 2014, 22(1): 18-17.

    [11]姚逸卿, 徐曉蘇, 童金武. 極區(qū)間接橫向慣性導(dǎo)航方法[J]. 中國慣性技術(shù)學(xué)報, 2015, 23(1): 29-34.Yao Y Q, Xu X S, Tong J W. Indirect transverse inertial navigation algorithm in polar region[J]. Journal of Chinese Inertial Technology, 2015, 23(1): 29-34.

    [12]劉猛, 高延濱, 李光春, 等. 基于重構(gòu)偽地球坐標(biāo)系的捷聯(lián)慣導(dǎo)初始對準(zhǔn)算法[J]. 中國慣性技術(shù)學(xué)報, 2016,24(6): 710-715.Liu M, Gao Y B, Li G C, et al. Initial alignment algorithm for SINS based on reconstructed pseudo-Earth frame[J].Journal of Chinese Inertial Technology, 2016, 24(6): 710-715.

    [13]Gao Y, Meng L, Li G, et al. Initial alignment for SINS based on pseudo-earth frame in polar regions[J]. Sensors,2017, 17(6): 1416.

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

    Algorithm of global navigation for inertial navigation system based on pseudo-Earth frame

    LIU Meng, LI Guang-chun, GAO Yan-bin, Li Ying, WANG Di
    (College of Automation, Harbin Engineering University, Harbin 150001, China)

    Inertial navigation system (INS), as an autonomous and self-contained system, is an optimal choice for the navigation in polar regions due to the existed special geographical and electromagnetic factors.Nevertheless, any of the existing conventional mechanizations could not achieve the global navigation by itself when the INS is considered to navigate in global regions. Therefore a composite mechanization is usually applied to achieve the global navigation, but it is adverse to the unification of the navigation algorithms in global regions. To solve this problem, a global navigation algorithm based on pseudo-Earth frame is proposed. The global navigation algorithm can not only be applied into strapdown INS but also can be applied into the platform INS, and can achieve the intrinsic unity of the navigation algorithms. Moreover,it is convenient to communicate with other local INS and has low complexity of switching logic, thereby decreasing the complexity of programming and reducing the burden of computer. Simulation results verify the feasibility and effectiveness of the proposed algorithm.

    pseudo-Earth frame; global navigation; polar navigation; strapdown inertial navigation system;platform inertial navigation system

    1005-6734(2017)04-0442-07

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

    U666.1

    A

    2017-04-10;

    2017-07-28

    中國國家科學(xué)技術(shù)部國際科技合作項目(2014DFR10010)

    劉猛(1988—),男,博士研究生,從事慣性技術(shù)研究。E-mail: liumeng_0304@163.com

    聯(lián) 系 人:高延濱(1963—),男,教授,博士生導(dǎo)師。E-mail: gaoyanbin@hrbeu.edu.cn

    猜你喜歡
    系統(tǒng)
    Smartflower POP 一體式光伏系統(tǒng)
    WJ-700無人機系統(tǒng)
    ZC系列無人機遙感系統(tǒng)
    北京測繪(2020年12期)2020-12-29 01:33:58
    基于PowerPC+FPGA顯示系統(tǒng)
    基于UG的發(fā)射箱自動化虛擬裝配系統(tǒng)開發(fā)
    半沸制皂系統(tǒng)(下)
    FAO系統(tǒng)特有功能分析及互聯(lián)互通探討
    連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
    一德系統(tǒng) 德行天下
    PLC在多段調(diào)速系統(tǒng)中的應(yīng)用

    一区二区三区免费毛片| 国内精品宾馆在线| 日日摸夜夜添夜夜爱| 亚洲第一区二区三区不卡| 日本91视频免费播放| 一区二区三区免费毛片| 女性被躁到高潮视频| 黑人欧美特级aaaaaa片| 国产熟女欧美一区二区| 日日啪夜夜爽| 日本爱情动作片www.在线观看| 婷婷色麻豆天堂久久| 欧美97在线视频| 久久国内精品自在自线图片| 18禁裸乳无遮挡动漫免费视频| 少妇高潮的动态图| 久久久久精品性色| 亚洲av免费高清在线观看| 又粗又硬又长又爽又黄的视频| av播播在线观看一区| 成人毛片a级毛片在线播放| 亚洲在久久综合| 黄色欧美视频在线观看| 免费黄网站久久成人精品| 秋霞在线观看毛片| 99久国产av精品国产电影| 综合色丁香网| 多毛熟女@视频| 中文欧美无线码| 伦理电影大哥的女人| 女的被弄到高潮叫床怎么办| 人妻少妇偷人精品九色| 午夜福利视频精品| 国产综合精华液| 国产在线一区二区三区精| 男人爽女人下面视频在线观看| 精品卡一卡二卡四卡免费| 丰满迷人的少妇在线观看| 精品午夜福利在线看| 国产国语露脸激情在线看| 亚洲精品亚洲一区二区| 制服丝袜香蕉在线| 伦理电影大哥的女人| 欧美+日韩+精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美亚洲日本最大视频资源| 妹子高潮喷水视频| 简卡轻食公司| 能在线免费看毛片的网站| 国精品久久久久久国模美| 一级毛片电影观看| 亚洲欧洲日产国产| 日韩成人av中文字幕在线观看| av国产精品久久久久影院| 一级毛片电影观看| 亚洲伊人久久精品综合| 日韩av在线免费看完整版不卡| 欧美日韩综合久久久久久| 亚洲无线观看免费| 精品一品国产午夜福利视频| 99久久人妻综合| 91精品国产九色| 最近的中文字幕免费完整| 亚洲综合精品二区| 亚洲色图综合在线观看| 91久久精品国产一区二区三区| 精品一区二区免费观看| 母亲3免费完整高清在线观看 | av网站免费在线观看视频| 亚洲经典国产精华液单| 国产成人一区二区在线| 久久狼人影院| 99久久人妻综合| 纯流量卡能插随身wifi吗| 国产黄色视频一区二区在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲欧洲精品一区二区精品久久久 | 人妻人人澡人人爽人人| 日本爱情动作片www.在线观看| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品第二区| 亚洲av免费高清在线观看| av国产久精品久网站免费入址| 国产高清国产精品国产三级| 午夜激情久久久久久久| 99热国产这里只有精品6| 天天躁夜夜躁狠狠久久av| av网站免费在线观看视频| av播播在线观看一区| 亚洲第一区二区三区不卡| 国产白丝娇喘喷水9色精品| 人妻 亚洲 视频| 国产成人精品一,二区| 亚洲美女搞黄在线观看| 亚洲国产成人一精品久久久| 亚洲精品中文字幕在线视频| 亚洲欧美日韩另类电影网站| 精品一区二区三区视频在线| 精品国产国语对白av| 亚洲欧美日韩卡通动漫| 国产精品一区二区在线不卡| 国产男人的电影天堂91| 日韩三级伦理在线观看| 一本—道久久a久久精品蜜桃钙片| 又大又黄又爽视频免费| 欧美日韩视频精品一区| 18在线观看网站| 欧美精品一区二区大全| 亚洲国产精品一区二区三区在线| 狠狠婷婷综合久久久久久88av| 成人亚洲精品一区在线观看| 18+在线观看网站| 人妻 亚洲 视频| 国产av一区二区精品久久| 赤兔流量卡办理| 男女无遮挡免费网站观看| 在线观看一区二区三区激情| 国产乱人偷精品视频| 中文字幕亚洲精品专区| 全区人妻精品视频| 亚洲在久久综合| 91精品一卡2卡3卡4卡| 亚洲精品久久成人aⅴ小说 | 国产av一区二区精品久久| 国产又色又爽无遮挡免| 国产一区有黄有色的免费视频| 黑人巨大精品欧美一区二区蜜桃 | 亚洲成人一二三区av| 国产日韩欧美视频二区| 久久99热6这里只有精品| 国产老妇伦熟女老妇高清| 日本欧美国产在线视频| 精品久久久精品久久久| 能在线免费看毛片的网站| 丰满少妇做爰视频| 国产精品一区二区在线不卡| 丰满迷人的少妇在线观看| 在线观看美女被高潮喷水网站| 色视频在线一区二区三区| 观看美女的网站| 高清不卡的av网站| 精品国产国语对白av| 啦啦啦视频在线资源免费观看| 国产亚洲午夜精品一区二区久久| 最新中文字幕久久久久| 亚洲国产精品专区欧美| 免费不卡的大黄色大毛片视频在线观看| 欧美日韩成人在线一区二区| 亚洲精品自拍成人| 久久久亚洲精品成人影院| 只有这里有精品99| 91国产中文字幕| 成人黄色视频免费在线看| 久久久亚洲精品成人影院| 亚洲精品乱码久久久久久按摩| 久久99蜜桃精品久久| 22中文网久久字幕| av在线观看视频网站免费| 这个男人来自地球电影免费观看 | 精品人妻偷拍中文字幕| 国产黄频视频在线观看| 3wmmmm亚洲av在线观看| 日韩伦理黄色片| 国产黄色免费在线视频| 国产 一区精品| 69精品国产乱码久久久| 一级片'在线观看视频| 青春草视频在线免费观看| 内地一区二区视频在线| 人妻系列 视频| 亚洲欧美清纯卡通| 欧美少妇被猛烈插入视频| 这个男人来自地球电影免费观看 | 亚洲av综合色区一区| 久久精品夜色国产| 亚洲人与动物交配视频| 一本一本综合久久| 成人无遮挡网站| 卡戴珊不雅视频在线播放| 国产精品久久久久久精品古装| 亚洲欧美色中文字幕在线| 丝瓜视频免费看黄片| 999精品在线视频| 亚洲精品456在线播放app| 午夜激情福利司机影院| 精品少妇内射三级| 好男人视频免费观看在线| 91久久精品国产一区二区成人| 人人澡人人妻人| 欧美日韩视频高清一区二区三区二| 亚洲欧美成人综合另类久久久| av福利片在线| 制服诱惑二区| 亚洲精品乱码久久久久久按摩| 亚洲av日韩在线播放| 三级国产精品欧美在线观看| 欧美日韩一区二区视频在线观看视频在线| 熟女电影av网| 在线看a的网站| 亚洲精品aⅴ在线观看| 亚洲国产精品一区二区三区在线| 三级国产精品欧美在线观看| 纵有疾风起免费观看全集完整版| 午夜福利在线观看免费完整高清在| av.在线天堂| 午夜福利网站1000一区二区三区| 亚洲精品,欧美精品| 最新中文字幕久久久久| 九九久久精品国产亚洲av麻豆| 日本爱情动作片www.在线观看| 91精品伊人久久大香线蕉| 国产老妇伦熟女老妇高清| 99热6这里只有精品| 国产深夜福利视频在线观看| 成人午夜精彩视频在线观看| 成人免费观看视频高清| 日韩一本色道免费dvd| 国产伦精品一区二区三区视频9| a级毛片免费高清观看在线播放| 黄色怎么调成土黄色| 热99国产精品久久久久久7| 另类精品久久| 日韩视频在线欧美| 国产在线视频一区二区| 91久久精品电影网| 九九久久精品国产亚洲av麻豆| 男女国产视频网站| 91精品伊人久久大香线蕉| 亚洲五月色婷婷综合| 新久久久久国产一级毛片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩熟女老妇一区二区性免费视频| 精品人妻偷拍中文字幕| 日本wwww免费看| 亚洲av福利一区| 日日爽夜夜爽网站| 午夜激情av网站| 亚洲少妇的诱惑av| 成年人免费黄色播放视频| 高清午夜精品一区二区三区| 国产熟女午夜一区二区三区 | 国产成人aa在线观看| 熟妇人妻不卡中文字幕| 亚洲第一av免费看| 久久久久国产网址| 大话2 男鬼变身卡| 欧美bdsm另类| 又粗又硬又长又爽又黄的视频| 亚洲,一卡二卡三卡| 在现免费观看毛片| 午夜精品国产一区二区电影| 久久这里有精品视频免费| 国产精品一区二区在线观看99| 2022亚洲国产成人精品| 亚洲综合色惰| 欧美性感艳星| 热99国产精品久久久久久7| 亚洲欧美日韩另类电影网站| 丁香六月天网| 亚洲无线观看免费| 久久人妻熟女aⅴ| 少妇丰满av| 欧美国产精品一级二级三级| 亚洲婷婷狠狠爱综合网| 韩国高清视频一区二区三区| 一个人免费看片子| 欧美激情极品国产一区二区三区 | 在线看a的网站| 国产伦精品一区二区三区视频9| 欧美少妇被猛烈插入视频| 一本久久精品| 女的被弄到高潮叫床怎么办| 久久久精品免费免费高清| 国产成人精品在线电影| 欧美日韩一区二区视频在线观看视频在线| 男女边摸边吃奶| 日韩伦理黄色片| 一级片'在线观看视频| 精品少妇久久久久久888优播| 一本大道久久a久久精品| 99热全是精品| 极品少妇高潮喷水抽搐| 蜜桃在线观看..| 亚洲欧美一区二区三区黑人 | 中文乱码字字幕精品一区二区三区| 国产欧美日韩一区二区三区在线 | av免费在线看不卡| 日本av手机在线免费观看| 国产精品99久久99久久久不卡 | 久久精品国产亚洲av天美| 精品人妻熟女毛片av久久网站| 亚洲av中文av极速乱| 少妇的逼好多水| 欧美 亚洲 国产 日韩一| 内地一区二区视频在线| 免费看光身美女| 国产精品成人在线| 久久久久人妻精品一区果冻| 曰老女人黄片| 国产日韩欧美在线精品| 全区人妻精品视频| 99久久中文字幕三级久久日本| 毛片一级片免费看久久久久| 丝袜在线中文字幕| 久久亚洲国产成人精品v| 伊人亚洲综合成人网| 一级片'在线观看视频| 欧美 亚洲 国产 日韩一| 免费黄频网站在线观看国产| 欧美最新免费一区二区三区| 哪个播放器可以免费观看大片| 制服人妻中文乱码| 亚洲五月色婷婷综合| 国产午夜精品一二区理论片| 男女边摸边吃奶| 啦啦啦在线观看免费高清www| 国产亚洲精品久久久com| 在线观看一区二区三区激情| 制服丝袜香蕉在线| 国产片内射在线| 亚洲精品视频女| 五月伊人婷婷丁香| 欧美日韩视频高清一区二区三区二| 99热6这里只有精品| 久热这里只有精品99| 午夜福利网站1000一区二区三区| 亚洲精品乱久久久久久| 美女内射精品一级片tv| 如何舔出高潮| 少妇被粗大的猛进出69影院 | 亚洲图色成人| 国产在线免费精品| 国产精品无大码| 人人澡人人妻人| 在线观看www视频免费| 99久久精品一区二区三区| 伦理电影大哥的女人| 免费观看a级毛片全部| 18禁在线播放成人免费| 亚洲五月色婷婷综合| 国产午夜精品一二区理论片| 高清黄色对白视频在线免费看| xxx大片免费视频| 老司机影院成人| tube8黄色片| 亚洲欧美一区二区三区黑人 | 纯流量卡能插随身wifi吗| 80岁老熟妇乱子伦牲交| 免费av不卡在线播放| 一区二区三区精品91| 少妇的逼好多水| 插阴视频在线观看视频| 免费黄网站久久成人精品| 91精品一卡2卡3卡4卡| 在线看a的网站| 亚洲国产精品一区三区| 特大巨黑吊av在线直播| 大陆偷拍与自拍| 久久人人爽人人爽人人片va| 欧美日韩av久久| 热99久久久久精品小说推荐| 国产黄色免费在线视频| 狂野欧美激情性xxxx在线观看| 久久久久国产精品人妻一区二区| 国产又色又爽无遮挡免| 91国产中文字幕| 国产高清有码在线观看视频| 高清毛片免费看| 美女脱内裤让男人舔精品视频| 97超视频在线观看视频| av女优亚洲男人天堂| 欧美变态另类bdsm刘玥| av在线观看视频网站免费| 日本爱情动作片www.在线观看| 亚洲人成网站在线播| 九草在线视频观看| 亚洲综合色惰| 日韩亚洲欧美综合| 在线观看免费视频网站a站| 最新的欧美精品一区二区| 日韩成人av中文字幕在线观看| 亚洲成人一二三区av| 亚洲中文av在线| 成人国产av品久久久| 精品视频人人做人人爽| 一级,二级,三级黄色视频| 亚洲情色 制服丝袜| 蜜桃久久精品国产亚洲av| 女性生殖器流出的白浆| 成人国语在线视频| 国产黄色视频一区二区在线观看| 国产伦理片在线播放av一区| 精品国产一区二区久久| 51国产日韩欧美| 国产熟女欧美一区二区| 青青草视频在线视频观看| 亚洲精品中文字幕在线视频| 成年人免费黄色播放视频| 国产极品粉嫩免费观看在线 | 久久久精品区二区三区| 免费大片黄手机在线观看| 精品久久国产蜜桃| 国产乱人偷精品视频| 国产成人精品无人区| 日韩强制内射视频| 欧美 日韩 精品 国产| 国产av精品麻豆| 又粗又硬又长又爽又黄的视频| 99热6这里只有精品| 久久鲁丝午夜福利片| 色吧在线观看| 久久久久国产精品人妻一区二区| 中国美白少妇内射xxxbb| 又大又黄又爽视频免费| 欧美 亚洲 国产 日韩一| h视频一区二区三区| 如日韩欧美国产精品一区二区三区 | 汤姆久久久久久久影院中文字幕| 18禁在线播放成人免费| 久久精品国产亚洲av天美| 两个人的视频大全免费| 有码 亚洲区| 超碰97精品在线观看| 国产av精品麻豆| √禁漫天堂资源中文www| 超碰97精品在线观看| 人体艺术视频欧美日本| av在线观看视频网站免费| 国产男女内射视频| 美女国产高潮福利片在线看| 久久毛片免费看一区二区三区| 在线观看免费日韩欧美大片 | 国产淫语在线视频| 高清毛片免费看| 人人妻人人澡人人爽人人夜夜| 中国国产av一级| 丝袜喷水一区| 欧美亚洲 丝袜 人妻 在线| 一级毛片我不卡| 午夜免费男女啪啪视频观看| 3wmmmm亚洲av在线观看| 免费日韩欧美在线观看| 熟女av电影| 亚洲国产av影院在线观看| 熟女人妻精品中文字幕| 亚洲国产日韩一区二区| 亚洲国产欧美在线一区| 欧美3d第一页| 精品熟女少妇av免费看| 黑人欧美特级aaaaaa片| 免费黄频网站在线观看国产| 视频在线观看一区二区三区| 精品国产国语对白av| 另类亚洲欧美激情| 精品久久久久久电影网| 久久精品夜色国产| 久久女婷五月综合色啪小说| √禁漫天堂资源中文www| 午夜精品国产一区二区电影| 永久免费av网站大全| 啦啦啦中文免费视频观看日本| 亚洲av成人精品一区久久| 如日韩欧美国产精品一区二区三区 | a级毛片免费高清观看在线播放| 欧美亚洲日本最大视频资源| 久久人妻熟女aⅴ| 最近中文字幕高清免费大全6| 成人综合一区亚洲| 久久精品国产自在天天线| 国产精品嫩草影院av在线观看| 99国产精品免费福利视频| 啦啦啦在线观看免费高清www| 女的被弄到高潮叫床怎么办| 人妻夜夜爽99麻豆av| 啦啦啦视频在线资源免费观看| 成年女人在线观看亚洲视频| www.色视频.com| a级毛色黄片| 亚洲欧美一区二区三区黑人 | 99精国产麻豆久久婷婷| 久久97久久精品| 国产精品人妻久久久久久| 免费观看a级毛片全部| 最近的中文字幕免费完整| 肉色欧美久久久久久久蜜桃| 国产av国产精品国产| 97超碰精品成人国产| 中国三级夫妇交换| 超碰97精品在线观看| 天天操日日干夜夜撸| 亚洲高清免费不卡视频| 另类亚洲欧美激情| 亚洲伊人久久精品综合| 亚洲精品第二区| 久久人妻熟女aⅴ| 一二三四中文在线观看免费高清| 国产av一区二区精品久久| 一级毛片 在线播放| 超碰97精品在线观看| 91久久精品电影网| 精品国产国语对白av| 人人妻人人爽人人添夜夜欢视频| 中文精品一卡2卡3卡4更新| 久久久久久久大尺度免费视频| 日韩av在线免费看完整版不卡| 亚洲综合精品二区| 岛国毛片在线播放| 久久午夜综合久久蜜桃| 欧美一级a爱片免费观看看| 最黄视频免费看| 97精品久久久久久久久久精品| 高清午夜精品一区二区三区| 久久97久久精品| 国产精品国产av在线观看| 久久久久视频综合| 久久久久久久国产电影| 在线观看免费日韩欧美大片 | 看免费成人av毛片| 国产精品一区二区在线不卡| 成人综合一区亚洲| 亚洲av男天堂| 久久久久久人妻| 国产 一区精品| 秋霞在线观看毛片| 国产女主播在线喷水免费视频网站| 精品卡一卡二卡四卡免费| 久久精品国产亚洲av涩爱| 狂野欧美激情性xxxx在线观看| 精品一区二区三卡| av免费在线看不卡| 九九爱精品视频在线观看| 制服诱惑二区| 熟女人妻精品中文字幕| 性高湖久久久久久久久免费观看| 搡老乐熟女国产| 两个人免费观看高清视频| 欧美精品国产亚洲| 免费黄频网站在线观看国产| 国产国语露脸激情在线看| 欧美亚洲日本最大视频资源| 国产熟女午夜一区二区三区 | 欧美精品一区二区大全| 精品国产国语对白av| av有码第一页| 久久99热这里只频精品6学生| 久久亚洲国产成人精品v| 多毛熟女@视频| 亚洲精品久久成人aⅴ小说 | 日韩欧美精品免费久久| 精品酒店卫生间| 久久青草综合色| 亚洲在久久综合| 18禁在线播放成人免费| 在线播放无遮挡| 97在线视频观看| 久久这里有精品视频免费| 高清av免费在线| 一级爰片在线观看| 久久久久久伊人网av| 乱人伦中国视频| 国产极品粉嫩免费观看在线 | 麻豆成人av视频| 18禁在线播放成人免费| 97超碰精品成人国产| 国产高清国产精品国产三级| 日韩一区二区视频免费看| 久久久国产一区二区| 亚洲三级黄色毛片| 午夜激情久久久久久久| 成人综合一区亚洲| 777米奇影视久久| 久久狼人影院| 一边亲一边摸免费视频| 日韩电影二区| 有码 亚洲区| 午夜视频国产福利| av一本久久久久| 亚洲精华国产精华液的使用体验| 欧美少妇被猛烈插入视频| 赤兔流量卡办理| 国产亚洲最大av| 日韩在线高清观看一区二区三区| 搡老乐熟女国产| 亚洲无线观看免费| 国产综合精华液| 不卡视频在线观看欧美| 免费观看性生交大片5| 人妻制服诱惑在线中文字幕| 一区二区三区四区激情视频| 精品酒店卫生间| av国产久精品久网站免费入址| 国产精品一区二区在线观看99| 久久国产亚洲av麻豆专区| 日韩熟女老妇一区二区性免费视频| 午夜福利视频精品| 波野结衣二区三区在线| 性色avwww在线观看| 国产成人精品无人区| 国产精品嫩草影院av在线观看| 免费观看av网站的网址| 免费观看的影片在线观看| 亚洲丝袜综合中文字幕| 18禁裸乳无遮挡动漫免费视频| 在线观看人妻少妇| 亚洲av男天堂| 午夜福利视频在线观看免费| 日韩大片免费观看网站| 国产色爽女视频免费观看| 最近最新中文字幕免费大全7| 亚洲国产精品国产精品| 色婷婷av一区二区三区视频| 新久久久久国产一级毛片| 七月丁香在线播放| 亚洲欧美色中文字幕在线| 自拍欧美九色日韩亚洲蝌蚪91|