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

    近地軌道航天器編隊(duì)僅測(cè)距相對(duì)導(dǎo)航方法

    2022-09-23 06:18:34龔柏春張偉夫
    宇航學(xué)報(bào) 2022年8期
    關(guān)鍵詞:鏡像測(cè)距編隊(duì)

    龔柏春,王 沙,張偉夫,周 亮,李 爽

    (1. 南京航空航天大學(xué)航天學(xué)院,南京 210016;2. 上海衛(wèi)星工程研究所,上海 201109;3. 航天科工空間工程發(fā)展有限公司,北京 100039)

    0 引 言

    編隊(duì)航天器因?yàn)榫哂邢噍^于單個(gè)航天器的許多優(yōu)勢(shì),成為學(xué)術(shù)界的研究熱點(diǎn)。例如,在太空任務(wù)中,單個(gè)航天器發(fā)生故障將直接導(dǎo)致任務(wù)失敗,而編隊(duì)航天器的冗余度設(shè)計(jì)可以在編隊(duì)某個(gè)成員發(fā)生故障時(shí)由備用航天器代替完成任務(wù)。除此之外,編隊(duì)航天器可以攜帶不同的載荷,完成單個(gè)航天器無(wú)法完成的多角度成像、通信中繼等任務(wù)。編隊(duì)航天器的飛行任務(wù)比較復(fù)雜,需要在協(xié)同飛行的過(guò)程中控制編隊(duì)成員之間的相對(duì)距離和方位,使整個(gè)編隊(duì)保持確定的幾何形態(tài),而完成編隊(duì)隊(duì)形控制的前提是編隊(duì)成員之間的精確相對(duì)導(dǎo)航。而要實(shí)現(xiàn)相對(duì)導(dǎo)航,集群成員之間首先需要進(jìn)行相對(duì)測(cè)量,其次需要建立可靠通信連接實(shí)現(xiàn)信息的實(shí)時(shí)交換。受航天器載荷體積和功耗的限制,光學(xué)相機(jī)僅測(cè)角和無(wú)線電僅測(cè)距成為航天器自主相對(duì)導(dǎo)航的主要發(fā)展方向。

    杜榮華等在僅有目標(biāo)視線角測(cè)量的條件下,提出了一種空間非合作目標(biāo)交會(huì)接近的自主檢測(cè)和跟蹤算法。周文雅等研究了一種改進(jìn)的基于圖像特征點(diǎn)的快速定位和簡(jiǎn)易旋轉(zhuǎn)算法,能夠提高雙目視覺無(wú)人機(jī)之間的相對(duì)定位精度。金澤明等聯(lián)合擴(kuò)展卡爾曼濾波算法和擴(kuò)展卡爾曼粒子濾波算法提高了空間非合作目標(biāo)單目視覺位姿估計(jì)的精度和計(jì)算效率。然而無(wú)源相機(jī)僅測(cè)角相對(duì)導(dǎo)航受光照條件制約,且相機(jī)有視場(chǎng)角限制,在測(cè)量過(guò)程中需要不斷調(diào)整來(lái)跟蹤目標(biāo)的方位,因此常用于非合作航天器的相對(duì)導(dǎo)航。

    合作航天器利用數(shù)據(jù)鏈測(cè)距實(shí)現(xiàn)相對(duì)導(dǎo)航時(shí),不需要考慮光學(xué)測(cè)量或無(wú)線電測(cè)量時(shí)的傳感器指向問(wèn)題,相對(duì)導(dǎo)航設(shè)備質(zhì)量輕、功耗低。Burchett等研究了飛行器在三維空間的僅測(cè)距相對(duì)導(dǎo)航算法,得出系統(tǒng)是局部弱可觀的結(jié)論。Wang等和Christian分析了Clohessy-Wiltshire(CW)動(dòng)力學(xué)模型下的僅測(cè)距模糊軌道,得到僅測(cè)距相對(duì)導(dǎo)航存在三種鏡像模糊軌道和四種變形模糊軌道的結(jié)論。Wang等也分析了TH動(dòng)力學(xué)模型下的僅測(cè)距模糊軌道,指出僅測(cè)距相對(duì)導(dǎo)航理論上存在三種鏡像模糊軌道,仿真中存在額外的四種變形模糊軌道。Gong等研究了近程圓或近圓軌道航天器編隊(duì)利用傳感器偏于質(zhì)心安裝提供的桿臂效應(yīng)來(lái)避免僅測(cè)距收斂于模糊軌道的方法,但是Wang團(tuán)隊(duì)和Gong等所采用的CW和TH方程均沒有考慮航天器在軌的各種攝動(dòng)因素。Shalaby等研究了利用在同一個(gè)飛行器安裝多個(gè)用于測(cè)距的通信信標(biāo),提升了僅測(cè)距相對(duì)導(dǎo)航系統(tǒng)的可觀測(cè)度。但是安裝在同一架飛行器上的信標(biāo)基線受飛行器實(shí)際尺寸的限制,因此該方案更適合于近距離或航天器交會(huì)對(duì)接場(chǎng)景下的相對(duì)導(dǎo)航。Xiong等通過(guò)在編隊(duì)中引入一架無(wú)人機(jī)作為參考信標(biāo),文獻(xiàn)[12-15]利用多個(gè)固定點(diǎn)作為參考信標(biāo),研究了超寬帶測(cè)距相對(duì)定位方法。Martel等利用超寬帶至少六次的距離測(cè)量信息和V-SLAM測(cè)量的各飛行器從各原點(diǎn)出發(fā)的位移信息,提出了一種獨(dú)特的四自由度相對(duì)位姿定位方案。

    本文針對(duì)近地軌道航天器編隊(duì)存在的自主相對(duì)導(dǎo)航需求,研究利用低成本數(shù)據(jù)鏈測(cè)距實(shí)現(xiàn)航天器相對(duì)導(dǎo)航的方法。論文推導(dǎo)了在地球非球形引力攝動(dòng)下的航天器高精度相對(duì)運(yùn)動(dòng)模型,分析了基于該動(dòng)力學(xué)模型設(shè)計(jì)的僅測(cè)距相對(duì)導(dǎo)航系統(tǒng)的可觀測(cè)性。此外,論文設(shè)計(jì)了以四個(gè)航天器為例的編隊(duì)分布式估計(jì)策略并研究了相對(duì)導(dǎo)航的誤差傳播規(guī)律,為大型航天器編隊(duì)的相對(duì)導(dǎo)航算法設(shè)計(jì)奠定了基礎(chǔ)。

    1 坐標(biāo)系定義

    本節(jié)定義描述航天器絕對(duì)運(yùn)動(dòng)的地心赤道慣性坐標(biāo)系和相對(duì)運(yùn)動(dòng)的航天器軌道坐標(biāo)系,第2節(jié)將在此基礎(chǔ)上推導(dǎo)航天器在J2攝動(dòng)下的相對(duì)運(yùn)動(dòng)方程。

    地心赤道慣性坐標(biāo)系:原點(diǎn)位于地心,軸指向春分點(diǎn),軸垂直于赤道平面并且與地球自轉(zhuǎn)角速度方向一致,軸位于赤道平面內(nèi)與軸、軸構(gòu)成右手直角坐標(biāo)系。

    航天器軌道坐標(biāo)系:原點(diǎn)位于航天器質(zhì)心,軸沿徑向方向由地心指向航天器質(zhì)心,軸垂直于軌道平面指向軌道角速度方向,位于軌道平面內(nèi)與軸、軸組成右手正交坐標(biāo)系。

    2 相對(duì)運(yùn)動(dòng)動(dòng)力學(xué)模型

    相對(duì)運(yùn)動(dòng)模型的精度直接影響著航天器相對(duì)導(dǎo)航的精度,已有的航天器相對(duì)導(dǎo)航研究所使用的動(dòng)力學(xué)模型都是基于無(wú)攝動(dòng)二體動(dòng)力學(xué)推導(dǎo)的CW方程或者TH方程,具有一定的局限性。對(duì)于近地軌道航天器,特別是中低軌道航天器,地球非球形引力攝動(dòng)是主要攝動(dòng)源之一,而J2攝動(dòng)又是地球非球形引力攝動(dòng)中的最主要部分。因此為了實(shí)現(xiàn)近地軌道航天器的高精度相對(duì)導(dǎo)航,本文首先要建立航天器在J2攝動(dòng)下的相對(duì)運(yùn)動(dòng)模型。

    2.1 J2攝動(dòng)相對(duì)運(yùn)動(dòng)模型

    記航天器在系下的位置為=[,,],則J2攝動(dòng)下航天器的動(dòng)力學(xué)方程如下

    (1)

    (2)

    (3)

    (4)

    (5)

    (6)

    (7)

    (8)

    (9)

    (10)

    (11)

    (12)

    (13)

    2.2 相對(duì)運(yùn)動(dòng)模型精度分析

    理論上,對(duì)于近地軌道航天器,J2攝動(dòng)相對(duì)運(yùn)動(dòng)模型因?yàn)榭紤]了地球非球形引力攝動(dòng)源中最主要的J2攝動(dòng)項(xiàng),因此比CW和TH相對(duì)運(yùn)動(dòng)模型具有更高的精度。為了對(duì)比三種相對(duì)運(yùn)動(dòng)模型的精度,進(jìn)行如下仿真,仿真中航天器絕對(duì)運(yùn)動(dòng)學(xué)只考慮地心引力和J2攝動(dòng)項(xiàng)。根據(jù)第1節(jié)航天器軌道坐標(biāo)系的定義,可以利用航天器在系下的位置換算得到系下的相對(duì)位置,并以此作為三種不同動(dòng)力學(xué)模型計(jì)算結(jié)果的參考,追蹤器和目標(biāo)器軌道六根數(shù)設(shè)置分別如表1中航天器1和航天器2所示。表1中,,,,,,分別代表軌道的半長(zhǎng)軸、偏心率、升交點(diǎn)赤經(jīng)、軌道傾角、近地點(diǎn)輻角和真近點(diǎn)角。航天器1和2相對(duì)距離為10 km到12 km,仿真步長(zhǎng)為10 s,總仿真時(shí)長(zhǎng)為2個(gè)軌道周期,約3.4 h。

    表1 航天器編隊(duì)軌道六根數(shù)設(shè)置Table 1 Classical orbital elements setting for spacecraft formation

    仿真結(jié)果中,隨著時(shí)間的推移,根據(jù)相對(duì)運(yùn)動(dòng)模型遞推的狀態(tài)誤差逐漸增加,這是因?yàn)檫@三種相對(duì)運(yùn)動(dòng)模型本質(zhì)上都是只保留了泰勒展開后一階項(xiàng)的線性相對(duì)運(yùn)動(dòng)模型。其中,CW相對(duì)運(yùn)動(dòng)模型狀態(tài)遞推誤差在0.5和1.5個(gè)軌道周期時(shí)達(dá)到最大值0.23%,TH相對(duì)運(yùn)動(dòng)模型狀態(tài)遞推誤差在2個(gè)軌道周期時(shí)達(dá)到最大值0.16%。J2攝動(dòng)相對(duì)運(yùn)動(dòng)模型精度最高,狀態(tài)遞推誤差在2個(gè)軌道周期時(shí)達(dá)到最大值0.06%。

    3 數(shù)據(jù)鏈測(cè)距模型

    合作航天器有兩種模式可以實(shí)現(xiàn)數(shù)據(jù)鏈測(cè)相對(duì)距離,其中一種是飛行時(shí)間差(Time difference of arrival, TDOA)測(cè)距模式。記Δ是航天器2信號(hào)接收時(shí)間與航天器1信號(hào)發(fā)送時(shí)間之差,由此計(jì)算航天器之間距離

    (14)

    式中:代表光速。TDOA測(cè)距是最簡(jiǎn)單的測(cè)距方法,航天器之間的通信可以是單向的,但是它要求航天器時(shí)鐘必須是同步的。另一種是飛行時(shí)間(Time-of-flight, TOF)測(cè)距模式,Δ是航天器2接收信號(hào)時(shí)間與航天器1發(fā)送信號(hào)時(shí)間之差,Δ是航天器1接收信號(hào)時(shí)間與航天器2發(fā)送信號(hào)時(shí)間之差,航天器相對(duì)距離為

    (15)

    利用TOF測(cè)距時(shí),航天器之間通信是雙向的,可以避免航天器鐘差帶來(lái)的測(cè)距誤差,因此本文將采用TOF測(cè)距模式。若時(shí)間測(cè)量誤差滿足高斯白噪聲分布,由于和Δ具有線性關(guān)系,則數(shù)據(jù)鏈測(cè)距可以建模為

    (16)

    4 可觀測(cè)性分析

    可觀測(cè)性分析是指在不考慮系統(tǒng)噪聲和觀測(cè)噪聲的條件下,根據(jù)歷史測(cè)量信息能否唯一確定系統(tǒng)狀態(tài)初值,它決定了濾波算法能否有效地估計(jì)出系統(tǒng)的狀態(tài)。本節(jié)將根據(jù)第2節(jié)建立的動(dòng)力學(xué)模型和第3節(jié)建立的觀測(cè)模型,分析相對(duì)導(dǎo)航系統(tǒng)可能存在的模糊軌道。因?yàn)镴2攝動(dòng)相對(duì)運(yùn)動(dòng)模型比較復(fù)雜,非線性系統(tǒng)可觀測(cè)性分析的李導(dǎo)數(shù)方法不再適用,所以本節(jié)將結(jié)合理論和數(shù)值仿真的方法分析系統(tǒng)的可觀測(cè)性。

    4.1 僅測(cè)距模糊軌道

    僅測(cè)距相對(duì)導(dǎo)航因?yàn)槿鄙傧鄬?duì)方位信息,系統(tǒng)可觀測(cè)性不足,往往存在模糊軌道。模糊軌道與真實(shí)軌道滿足相同動(dòng)力學(xué)模型并且具有相同測(cè)距輸出,它與真實(shí)軌道的區(qū)別在于系統(tǒng)狀態(tài)的初值不同。Wang等的研究結(jié)果表明,僅測(cè)距相對(duì)導(dǎo)航可能存在的模糊軌道一共7種,見表2。其中,(a)、(b)和(c)類模糊軌道因?yàn)楹驼鎸?shí)軌道形狀大小相同,被稱為鏡像模糊軌道;(d)、(e)、(f)和(g)類模糊軌道和真實(shí)軌道形狀和大小不同,被稱為變形模糊軌道。

    表2 僅測(cè)距模糊軌道初始狀態(tài)Table 2 Initial state of range-only ambiguous relative orbits

    4.2 可觀測(cè)性理論分析

    對(duì)于如下標(biāo)準(zhǔn)線性時(shí)變系統(tǒng)

    (17)

    根據(jù)式(17)設(shè)計(jì)的僅測(cè)距相對(duì)導(dǎo)航系統(tǒng),記模糊軌道初始狀態(tài)為′=-??梢宰C明,在動(dòng)力學(xué)模型的約束下,該模糊軌道在任意時(shí)刻的狀態(tài)量′()滿足

    ′()=-()

    (18)

    它與真實(shí)軌道狀態(tài)()是相反數(shù)關(guān)系,因此與真實(shí)軌道具有同樣測(cè)距輸出,是如表2所示的(c)類鏡像模糊軌道。CW方程、TH方程和本文所建立的J2攝動(dòng)相對(duì)運(yùn)動(dòng)方程都可以表示為式(17),因此基于這三種相對(duì)運(yùn)動(dòng)模型的僅測(cè)距相對(duì)導(dǎo)航系統(tǒng)都存在該類鏡像模糊軌道。但是對(duì)于式(17)所示的線性時(shí)變系統(tǒng),因?yàn)榭刂屏康妮斎?初始狀態(tài)為′=-的模糊軌道將不一定滿足式(18),因此可能避免模糊軌道。這種思想與Woffinden等通過(guò)軌道機(jī)動(dòng)來(lái)解決僅測(cè)角相對(duì)導(dǎo)航的不可觀測(cè)問(wèn)題很相似。

    (19)

    4.3 可觀測(cè)性數(shù)值分析

    Wang等在論文中通過(guò)數(shù)值仿真發(fā)現(xiàn),保留CW方程或者TH方程推導(dǎo)過(guò)程中泰勒展開的二階項(xiàng)、三階項(xiàng),可以提高僅測(cè)距相對(duì)導(dǎo)航系統(tǒng)的可觀測(cè)性,避免濾波收斂于模糊軌道。下面通過(guò)數(shù)值仿真分析式(9)中的J2攝動(dòng)項(xiàng),是否同樣有抑制濾波發(fā)散的作用。對(duì)于表1中的航天器1和2,使用無(wú)味卡爾曼濾波(Unscented Kalman filter, UKF)進(jìn)行相對(duì)運(yùn)動(dòng)狀態(tài)的估計(jì),仿真中假設(shè)追蹤器在系下的位置和速度精確已知。濾波關(guān)鍵參數(shù)設(shè)置見表3,總仿真時(shí)長(zhǎng)為10000 s,約2個(gè)軌道周期。

    表3 濾波算法參數(shù)設(shè)置Table 3 Simulation parameters setting for Kalman filter

    定義僅測(cè)距相對(duì)導(dǎo)航收斂域?yàn)槟苁篂V波收斂于真實(shí)軌道的初始估計(jì)誤差集合,收斂域越大,表示系統(tǒng)的可觀測(cè)性越好。但是在比較考慮J2攝動(dòng)和不考慮J2攝動(dòng)(采用TH模型)系統(tǒng)的收斂域時(shí),因?yàn)椴荒芨F舉所有的初始估計(jì)誤差,因此選取如表4 所示的8個(gè)方向進(jìn)行比較。

    表4 初始估計(jì)誤差的八個(gè)方向Table 4 Initial error setting for the eight directions

    (20)

    圖1是兩種相對(duì)運(yùn)動(dòng)模型的距離估計(jì)誤差圖,由圖1(a)可知,僅測(cè)距相對(duì)導(dǎo)航不同方向的收斂性是不一樣的。通過(guò)對(duì)比圖1(a)和圖1(b)可知,J2攝動(dòng)模型和TH模型在選定的8個(gè)方向上具有相同的收斂性。另外,J2攝動(dòng)模型具有更高的相對(duì)導(dǎo)航精度。例如在方向1上,當(dāng)=03時(shí),TH模型對(duì)應(yīng)的RMSE為245 m,而J2攝動(dòng)模型對(duì)應(yīng)的RMSE為19 m。因此,可以預(yù)測(cè)J2攝動(dòng)模型與TH模型存在相同的收斂域,都至少存在如表2所示的(a)、(b)和(c)三種鏡像模糊軌道,下面將對(duì)該結(jié)論進(jìn)行驗(yàn)證。

    圖1 航天器距離估計(jì)誤差Fig.1 Graphs of relative distance estimation error

    4.4 可觀測(cè)性仿真校驗(yàn)

    無(wú)味卡爾曼濾波4組狀態(tài)初值設(shè)置見表2中真實(shí)軌道和鏡像模糊軌道狀態(tài)初值,其余濾波參數(shù)設(shè)置見表3。雙星編隊(duì)的軌道參數(shù)設(shè)置見表1中航天器1和2,通過(guò)設(shè)置不同濾波初值可以分別得到真實(shí)軌道和鏡像模糊軌道。圖2是J2攝動(dòng)下僅測(cè)距相對(duì)導(dǎo)航系統(tǒng)的鏡像模糊軌道與真實(shí)軌道。

    圖2 僅測(cè)距相對(duì)導(dǎo)航鏡像模糊軌道Fig.2 Mirror ambiguous orbits for range-only relative navigation

    從圖2可以看出,因?yàn)镴2攝動(dòng)項(xiàng)的存在,藍(lán)色實(shí)線表示的真實(shí)軌道不再閉合。但是三種鏡像模糊軌道與真實(shí)軌道在三維空間內(nèi)仍然具有相同的形狀和大小,呈現(xiàn)“鏡像”特性。驗(yàn)證了之前的可觀測(cè)性分析結(jié)論:

    1) 對(duì)于式(17)所示的僅測(cè)距相對(duì)導(dǎo)航系統(tǒng),始終存在鏡像模糊軌道;

    2) 根據(jù)式(9)與TH方程設(shè)計(jì)的僅測(cè)距相對(duì)導(dǎo)航系統(tǒng)可觀測(cè)性一樣,都至少存在三種鏡像模糊軌道。

    5 編隊(duì)航天器分布式相對(duì)導(dǎo)航

    Wang團(tuán)隊(duì)研究了由三個(gè)航天器組成編隊(duì)的分布式相對(duì)導(dǎo)航算法,并利用一致性濾波來(lái)提高系統(tǒng)的可觀測(cè)性。三個(gè)航天器是相對(duì)簡(jiǎn)單的航天器編隊(duì),本節(jié)以四個(gè)航天器為例,研究不同分布式估計(jì)策略的相對(duì)導(dǎo)航精度與導(dǎo)航誤差傳播特性。

    5.1 分布式估計(jì)策略

    設(shè)無(wú)向圖由頂點(diǎn)集合和邊集合組成,頂點(diǎn)表示編隊(duì)中的航天器,邊表示兩個(gè)航天器相鄰,它們之間相互通信并估計(jì)彼此的相對(duì)運(yùn)動(dòng)狀態(tài)。為了實(shí)現(xiàn)編隊(duì)中任意兩個(gè)成員的相對(duì)導(dǎo)航,無(wú)向圖必須是連通的。這樣即使兩個(gè)航天器在無(wú)向圖中不相鄰,它們之間的相對(duì)運(yùn)動(dòng)狀態(tài)也可以通過(guò)如下坐標(biāo)轉(zhuǎn)換矩陣與矢量運(yùn)算得到

    (21)

    (22)

    (23)

    在某次更新相對(duì)運(yùn)動(dòng)狀態(tài)時(shí),無(wú)向圖的邊表示執(zhí)行一次濾波算法,記計(jì)算量為。若無(wú)向圖G的邊個(gè)數(shù)為,則編隊(duì)總的計(jì)算量為。圖3是四個(gè)航天器組成編隊(duì)的分布式估計(jì)策略,圖3(a)是最簡(jiǎn)單的情形,其計(jì)算量也最少(為3),圖3(b)、(c)和(d)估計(jì)策略的計(jì)算量分別為4,5和6。以圖3(a) 為例,航天器2相對(duì)于航天器1的運(yùn)動(dòng)信息可以通過(guò)濾波算法獲得,而航天器3和航天器4相對(duì)于航天器1的運(yùn)動(dòng)信息需要通過(guò)式(21)所示的關(guān)系計(jì)算得到。

    圖3 四個(gè)航天器組成編隊(duì)的分布式估計(jì)策略Fig.3 Decentralized estimation strategy for spacecraft formation

    5.2 相對(duì)導(dǎo)航誤差傳播規(guī)律

    為了研究編隊(duì)相對(duì)導(dǎo)航誤差傳播規(guī)律,本節(jié)采用蒙特卡洛打靶法。軌道和濾波參數(shù)設(shè)置分別見表1和表3,濾波初值誤差δ為相對(duì)運(yùn)動(dòng)狀態(tài)初值的0.1倍。另外,與4.3節(jié)可觀測(cè)性數(shù)值分析的仿真中假設(shè)追蹤器在系下的位置和速度精確已知不同,本節(jié)與第6節(jié)的仿真中均為表1所示四個(gè)航天器在系下的初始位置和速度添加均值為0、方差分別為10 m和0.1 m·s的高斯誤差。

    以圖3(a)最簡(jiǎn)單三條邊估計(jì)策略為例,圖4是航天器2、3、4相對(duì)于航天器1的200次蒙特卡洛打靶仿真的距離估計(jì)誤差統(tǒng)計(jì)結(jié)果。圖4中藍(lán)色雙劃線、橙色實(shí)線和紅色點(diǎn)劃線分別代表估計(jì)誤差的+3σ邊界、均值和-3σ邊界,邊界值表示導(dǎo)航精度,邊界范圍表示導(dǎo)航穩(wěn)定性。8000 s時(shí),航天器1和2相對(duì)距離估計(jì)誤差在13 m到18 m之間。相比之下,8000 s時(shí)航天器1和3相對(duì)距離估計(jì)誤差較大,在42 m到74 m之間。8000 s時(shí)航天器1和4相對(duì)距離估計(jì)誤差最大,在77 m到150 m之間。

    圖4 相鄰和非相鄰成員相對(duì)距離估計(jì)誤差Fig.4 Relative distance estimation error for nearby and no-nearby members

    (24)

    (25)

    式中:表示相關(guān)的狀態(tài)量估計(jì)誤差。另外,不同航天器進(jìn)行相對(duì)導(dǎo)航時(shí)的濾波算法和數(shù)據(jù)鏈測(cè)距結(jié)果是相互獨(dú)立的,相對(duì)導(dǎo)航誤差也是相互獨(dú)立的。式(24)和式(25)的導(dǎo)航誤差是近似線性疊加的,這揭示了編隊(duì)航天器相對(duì)導(dǎo)航的誤差傳播規(guī)律。因此,在實(shí)現(xiàn)編隊(duì)中兩個(gè)航天器的相對(duì)導(dǎo)航時(shí),應(yīng)盡可能利用這兩個(gè)航天器的相對(duì)測(cè)距信息并進(jìn)行濾波。

    6 編隊(duì)航天器一致性算法

    式(21)是航天器相對(duì)運(yùn)動(dòng)狀態(tài)需要滿足的物理約束,通過(guò)在濾波狀態(tài)更新算法中引入該等式確定的先驗(yàn)估計(jì),可以設(shè)計(jì)一致性濾波算法。一致性濾波算法與UKF算法結(jié)合,便得到一致性無(wú)味卡爾曼濾波算法(CUKF)。因?yàn)闀r(shí)間更新方程和UKF相同,下面僅給出一致性濾波算法的狀態(tài)更新方程

    (26)

    以圖3(b)左邊四邊形環(huán)估計(jì)策略為例,研究CUKF算法下航天器相對(duì)導(dǎo)航系統(tǒng)可觀測(cè)性。仿真的軌道和濾波參數(shù)設(shè)置見表1和表3,令航天器1對(duì)航天器2的初始相對(duì)導(dǎo)航誤差為δ=2δ,δ見表4。圖5對(duì)比了UKF和CUKF濾波結(jié)果的相對(duì)距離估計(jì)誤差,由圖可知,在CUKF算法下僅測(cè)距相對(duì)導(dǎo)航系統(tǒng)是收斂的,驗(yàn)證了一致性算法對(duì)于提高系統(tǒng)可觀測(cè)性的作用。

    圖5 UKF與CUKF相對(duì)距離估計(jì)誤差Fig.5 Relative distance estimation error under UKF and CUKF

    7 結(jié) 論

    本文研究了近地軌道航天器在不完備測(cè)量即僅有相對(duì)距離情形下的相對(duì)導(dǎo)航問(wèn)題,建立了航天器在地球非球形引力J2攝動(dòng)下的僅測(cè)距相對(duì)導(dǎo)航系統(tǒng),并對(duì)系統(tǒng)進(jìn)行了可觀測(cè)性分析和分布式估計(jì)策略設(shè)計(jì)。論文的主要研究成果包括四個(gè)方面:

    1) 推導(dǎo)了航天器在地球非球形引力攝動(dòng)下線性的J2攝動(dòng)相對(duì)運(yùn)動(dòng)模型,該模型適用于航天器運(yùn)行在圓或橢圓軌道的情形;

    2) 可觀測(cè)性分析結(jié)果表明,利用J2攝動(dòng)相對(duì)運(yùn)動(dòng)模型設(shè)計(jì)的僅測(cè)距相對(duì)導(dǎo)航系統(tǒng)不可觀測(cè),存在三種鏡像模糊軌道;

    3) 對(duì)于四個(gè)航天器組成的編隊(duì),設(shè)計(jì)了六種不同的分布式估計(jì)策略,理論上分析并通過(guò)數(shù)值仿真校驗(yàn)了航天器相對(duì)導(dǎo)航的誤差具有近似線性傳播的特性;

    4) 對(duì)比了一致性算法和無(wú)味卡爾曼濾波算法的蒙特卡洛打靶結(jié)果,驗(yàn)證了一致性算法對(duì)于提高系統(tǒng)可觀測(cè)性的作用。

    下一步的研究將從利用航天器軌道機(jī)動(dòng)提高僅測(cè)距相對(duì)導(dǎo)航系統(tǒng)可觀測(cè)性的角度出發(fā),研究?jī)H測(cè)距編隊(duì)控制問(wèn)題。

    猜你喜歡
    鏡像測(cè)距編隊(duì)
    2023年1月25日,美軍一次演習(xí)期間,空軍正在進(jìn)行編隊(duì)飛行
    軍事文摘(2023年5期)2023-03-27 08:56:26
    鏡像
    類星體的精準(zhǔn)測(cè)距
    科學(xué)(2020年3期)2020-01-06 04:02:51
    鏡像
    小康(2018年23期)2018-08-23 06:18:52
    基于事件驅(qū)動(dòng)的多飛行器編隊(duì)協(xié)同控制
    淺談超聲波測(cè)距
    電子制作(2017年7期)2017-06-05 09:36:13
    鏡像
    小康(2015年4期)2015-03-31 14:57:40
    鏡像
    小康(2015年6期)2015-03-26 14:44:27
    基于PSOC超聲測(cè)距系統(tǒng)設(shè)計(jì)
    基于預(yù)測(cè)控制的無(wú)人機(jī)編隊(duì)內(nèi)部避碰
    搡老熟女国产l中国老女人| 国产精品亚洲av一区麻豆| 国产日本99.免费观看| 久久久久精品国产欧美久久久| 国产亚洲精品av在线| 国内精品美女久久久久久| 久久99热这里只有精品18| 露出奶头的视频| 香蕉国产在线看| 高清在线国产一区| 精品不卡国产一区二区三区| 国产精品久久久人人做人人爽| 最近最新中文字幕大全免费视频| 很黄的视频免费| 免费在线观看日本一区| 黄色 视频免费看| 午夜福利欧美成人| 我要搜黄色片| 久久久久国内视频| 欧美极品一区二区三区四区| 日韩欧美在线乱码| 久久久久久人人人人人| 亚洲av电影在线进入| 亚洲av第一区精品v没综合| 日韩有码中文字幕| 法律面前人人平等表现在哪些方面| 国产激情欧美一区二区| 欧美不卡视频在线免费观看| 日韩欧美精品v在线| 天堂√8在线中文| 亚洲国产欧美一区二区综合| 日本三级黄在线观看| 国产亚洲精品综合一区在线观看| 在线观看舔阴道视频| 日本三级黄在线观看| 国产精品精品国产色婷婷| 亚洲精品久久国产高清桃花| 免费看十八禁软件| 国产亚洲av嫩草精品影院| 午夜福利在线观看吧| 禁无遮挡网站| 中文字幕人妻丝袜一区二区| 蜜桃久久精品国产亚洲av| 我要搜黄色片| 嫩草影院精品99| 亚洲中文av在线| 色综合站精品国产| 90打野战视频偷拍视频| 国产三级黄色录像| 亚洲av美国av| 国产主播在线观看一区二区| 美女高潮喷水抽搐中文字幕| 欧美日韩中文字幕国产精品一区二区三区| 最近视频中文字幕2019在线8| 精品国产乱码久久久久久男人| 亚洲av第一区精品v没综合| 免费看十八禁软件| 精品无人区乱码1区二区| 欧美中文综合在线视频| 国产aⅴ精品一区二区三区波| 窝窝影院91人妻| netflix在线观看网站| 少妇的逼水好多| 小蜜桃在线观看免费完整版高清| 亚洲中文日韩欧美视频| 一本综合久久免费| 白带黄色成豆腐渣| www.www免费av| 91av网站免费观看| 亚洲av第一区精品v没综合| 国产成人精品无人区| 91在线精品国自产拍蜜月 | 久久天堂一区二区三区四区| 黄片小视频在线播放| 日韩欧美三级三区| 亚洲熟妇熟女久久| 国产三级黄色录像| 免费电影在线观看免费观看| 91av网站免费观看| 老司机午夜十八禁免费视频| 国产aⅴ精品一区二区三区波| 在线十欧美十亚洲十日本专区| 色哟哟哟哟哟哟| 97碰自拍视频| 午夜福利欧美成人| 中文字幕熟女人妻在线| 亚洲七黄色美女视频| 香蕉国产在线看| www.精华液| 美女被艹到高潮喷水动态| 国产精品98久久久久久宅男小说| 成人高潮视频无遮挡免费网站| 两个人的视频大全免费| 精品国产三级普通话版| 女生性感内裤真人,穿戴方法视频| 色播亚洲综合网| 国产欧美日韩一区二区精品| 最新在线观看一区二区三区| 级片在线观看| 精品久久久久久久人妻蜜臀av| 欧美高清成人免费视频www| 亚洲欧美日韩东京热| 欧美三级亚洲精品| 少妇人妻一区二区三区视频| 琪琪午夜伦伦电影理论片6080| 成人午夜高清在线视频| 毛片女人毛片| 两人在一起打扑克的视频| 欧美大码av| 一本综合久久免费| 国内精品美女久久久久久| 国产精品久久久av美女十八| 男女床上黄色一级片免费看| 日本免费a在线| 一本一本综合久久| 国产成人啪精品午夜网站| 免费看日本二区| 午夜久久久久精精品| 欧美国产日韩亚洲一区| 免费在线观看日本一区| 欧美xxxx黑人xx丫x性爽| 这个男人来自地球电影免费观看| 国产单亲对白刺激| 欧美中文日本在线观看视频| 亚洲人成伊人成综合网2020| 欧美日韩福利视频一区二区| 免费看十八禁软件| 18美女黄网站色大片免费观看| 日韩av在线大香蕉| 亚洲欧美一区二区三区黑人| 欧美3d第一页| 啦啦啦韩国在线观看视频| 久久久久久久久中文| 婷婷精品国产亚洲av| 午夜福利欧美成人| 露出奶头的视频| 搡老妇女老女人老熟妇| 性色av乱码一区二区三区2| 免费看美女性在线毛片视频| 国产精品自产拍在线观看55亚洲| 久久精品综合一区二区三区| 亚洲av日韩精品久久久久久密| 免费高清视频大片| 亚洲国产欧洲综合997久久,| 丰满人妻一区二区三区视频av | 不卡av一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| av福利片在线观看| 男人舔奶头视频| 无遮挡黄片免费观看| 精品国产三级普通话版| 琪琪午夜伦伦电影理论片6080| 亚洲黑人精品在线| 97超级碰碰碰精品色视频在线观看| 久久国产精品人妻蜜桃| 精品国产亚洲在线| 最好的美女福利视频网| 91麻豆精品激情在线观看国产| 欧美一区二区国产精品久久精品| 成人特级av手机在线观看| 一夜夜www| 日本a在线网址| 成年人黄色毛片网站| 亚洲第一欧美日韩一区二区三区| 美女扒开内裤让男人捅视频| 日韩欧美在线二视频| 可以在线观看毛片的网站| 成年版毛片免费区| 国产极品精品免费视频能看的| 国产亚洲av高清不卡| 欧美日韩黄片免| av中文乱码字幕在线| 国产又色又爽无遮挡免费看| 波多野结衣高清作品| 搡老妇女老女人老熟妇| 这个男人来自地球电影免费观看| 亚洲欧美日韩高清专用| 亚洲欧美激情综合另类| 国产精品亚洲av一区麻豆| 精品午夜福利视频在线观看一区| 亚洲天堂国产精品一区在线| 免费观看精品视频网站| 欧美绝顶高潮抽搐喷水| 人妻夜夜爽99麻豆av| 中文字幕人成人乱码亚洲影| 亚洲国产精品久久男人天堂| 19禁男女啪啪无遮挡网站| 午夜亚洲福利在线播放| 国产成人aa在线观看| 老司机午夜福利在线观看视频| 久久人妻av系列| 国产1区2区3区精品| 亚洲专区字幕在线| 亚洲成人免费电影在线观看| 久久久久亚洲av毛片大全| 神马国产精品三级电影在线观看| 热99在线观看视频| 天堂动漫精品| 天堂√8在线中文| 国产一区二区三区在线臀色熟女| www.www免费av| 免费在线观看成人毛片| 欧美大码av| 久久久久国产精品人妻aⅴ院| 91字幕亚洲| 亚洲精品色激情综合| 久久中文字幕人妻熟女| 午夜福利视频1000在线观看| 性欧美人与动物交配| 精品无人区乱码1区二区| 亚洲欧美精品综合久久99| 欧洲精品卡2卡3卡4卡5卡区| 久久热在线av| 夜夜躁狠狠躁天天躁| 99久久99久久久精品蜜桃| 国产精品1区2区在线观看.| 身体一侧抽搐| 亚洲国产精品sss在线观看| 久久精品国产亚洲av香蕉五月| 亚洲熟妇中文字幕五十中出| 美女大奶头视频| 一二三四社区在线视频社区8| 日本a在线网址| 一级作爱视频免费观看| 欧美极品一区二区三区四区| 91在线观看av| 一区二区三区激情视频| 一级a爱片免费观看的视频| 哪里可以看免费的av片| 午夜亚洲福利在线播放| 性色avwww在线观看| 18禁黄网站禁片午夜丰满| 男女那种视频在线观看| 天天一区二区日本电影三级| 亚洲av成人不卡在线观看播放网| 亚洲精品美女久久av网站| 后天国语完整版免费观看| www.熟女人妻精品国产| 欧美日本视频| 久久久久久九九精品二区国产| 久久久久久久久久黄片| 久久亚洲真实| 中出人妻视频一区二区| 中文字幕高清在线视频| 亚洲 国产 在线| 五月玫瑰六月丁香| 51午夜福利影视在线观看| 亚洲成av人片免费观看| 看黄色毛片网站| 午夜福利高清视频| 丰满的人妻完整版| 国产精品爽爽va在线观看网站| 丝袜人妻中文字幕| 免费在线观看成人毛片| 国产精品一区二区免费欧美| 亚洲av五月六月丁香网| 少妇丰满av| 精品免费久久久久久久清纯| 好男人在线观看高清免费视频| 日韩欧美在线二视频| 天天一区二区日本电影三级| 午夜影院日韩av| 午夜福利在线观看吧| 色av中文字幕| 在线国产一区二区在线| 午夜福利高清视频| 最近视频中文字幕2019在线8| 1024手机看黄色片| 亚洲天堂国产精品一区在线| 久久香蕉精品热| 1024香蕉在线观看| 性色avwww在线观看| 国产 一区 欧美 日韩| 一级黄色大片毛片| 色吧在线观看| 久久久久国内视频| 最新美女视频免费是黄的| 999久久久精品免费观看国产| 国产精品一区二区三区四区久久| 1024手机看黄色片| 观看免费一级毛片| 欧美日韩综合久久久久久 | 国产成人啪精品午夜网站| 午夜久久久久精精品| 999久久久精品免费观看国产| 19禁男女啪啪无遮挡网站| 久久午夜亚洲精品久久| 啦啦啦观看免费观看视频高清| 欧美日韩黄片免| 国产精品一及| 色精品久久人妻99蜜桃| 精品午夜福利视频在线观看一区| 看片在线看免费视频| 悠悠久久av| 亚洲av美国av| 国模一区二区三区四区视频 | 久久欧美精品欧美久久欧美| 偷拍熟女少妇极品色| 国产成人aa在线观看| 天堂网av新在线| 午夜两性在线视频| 亚洲精品粉嫩美女一区| 免费看十八禁软件| 97超级碰碰碰精品色视频在线观看| 国内精品久久久久精免费| 久久国产精品影院| 黄色女人牲交| 久久亚洲真实| 亚洲av成人精品一区久久| www日本黄色视频网| 久久天躁狠狠躁夜夜2o2o| 国产精品 欧美亚洲| 12—13女人毛片做爰片一| 国产成人福利小说| 亚洲国产色片| 看黄色毛片网站| 97碰自拍视频| 噜噜噜噜噜久久久久久91| 午夜福利在线观看吧| 久久这里只有精品19| 国产aⅴ精品一区二区三区波| 国产黄片美女视频| 日韩成人在线观看一区二区三区| 香蕉丝袜av| 女人被狂操c到高潮| 亚洲午夜理论影院| 亚洲av成人不卡在线观看播放网| 亚洲av电影在线进入| 狂野欧美白嫩少妇大欣赏| 天堂av国产一区二区熟女人妻| 国产精品久久久久久久电影 | 熟女电影av网| avwww免费| 中国美女看黄片| 伦理电影免费视频| 精品一区二区三区视频在线观看免费| 亚洲熟妇中文字幕五十中出| 欧美日韩一级在线毛片| 国产伦精品一区二区三区四那| 欧美一区二区精品小视频在线| 午夜福利免费观看在线| 国内精品久久久久久久电影| 欧美日韩国产亚洲二区| 亚洲国产欧美一区二区综合| 极品教师在线免费播放| 日本撒尿小便嘘嘘汇集6| 久久久久免费精品人妻一区二区| 18禁黄网站禁片免费观看直播| 嫩草影院入口| 午夜久久久久精精品| 麻豆成人av在线观看| 香蕉av资源在线| 久久久国产成人免费| 久久久久久久久中文| 91字幕亚洲| av视频在线观看入口| 国产综合懂色| 亚洲av成人不卡在线观看播放网| 亚洲九九香蕉| 日韩有码中文字幕| 中文在线观看免费www的网站| 亚洲精品在线美女| 国产精品综合久久久久久久免费| 亚洲精品粉嫩美女一区| 色视频www国产| 欧美日韩精品网址| 91麻豆精品激情在线观看国产| 两性夫妻黄色片| 国产精品电影一区二区三区| 丁香欧美五月| 国模一区二区三区四区视频 | 日韩欧美国产一区二区入口| 亚洲av电影在线进入| 91麻豆av在线| 麻豆成人av在线观看| 欧美性猛交╳xxx乱大交人| 熟女少妇亚洲综合色aaa.| 日本黄大片高清| 亚洲av电影在线进入| 国产麻豆成人av免费视频| 99久久99久久久精品蜜桃| 午夜福利18| 黄频高清免费视频| 熟女少妇亚洲综合色aaa.| 婷婷精品国产亚洲av| 1024香蕉在线观看| 午夜精品在线福利| 真人做人爱边吃奶动态| 婷婷精品国产亚洲av在线| or卡值多少钱| 桃色一区二区三区在线观看| 日本熟妇午夜| 在线永久观看黄色视频| 男人的好看免费观看在线视频| 亚洲欧美日韩高清在线视频| 免费在线观看影片大全网站| 欧美一区二区国产精品久久精品| 在线看三级毛片| 国产精品乱码一区二三区的特点| 午夜a级毛片| 国模一区二区三区四区视频 | 亚洲美女视频黄频| 日本一二三区视频观看| 美女被艹到高潮喷水动态| 日韩成人在线观看一区二区三区| 禁无遮挡网站| 亚洲av片天天在线观看| 精华霜和精华液先用哪个| 欧美成人一区二区免费高清观看 | 国产一区在线观看成人免费| 欧美黑人巨大hd| 欧美一级a爱片免费观看看| 最好的美女福利视频网| 夜夜爽天天搞| 久久精品亚洲精品国产色婷小说| 亚洲激情在线av| 精品无人区乱码1区二区| 国产精品1区2区在线观看.| or卡值多少钱| 伦理电影免费视频| 最好的美女福利视频网| 丝袜人妻中文字幕| 一级毛片高清免费大全| 国产又色又爽无遮挡免费看| 国产成人aa在线观看| 一二三四在线观看免费中文在| 国产av不卡久久| 亚洲欧美激情综合另类| 亚洲国产精品合色在线| 国产欧美日韩精品亚洲av| 亚洲av电影不卡..在线观看| 人妻夜夜爽99麻豆av| 香蕉丝袜av| 午夜久久久久精精品| 国产精品影院久久| 免费看十八禁软件| 久久精品国产综合久久久| 免费av毛片视频| 操出白浆在线播放| 国产精品国产高清国产av| 久久精品影院6| 免费看日本二区| 日本黄大片高清| 国产伦人伦偷精品视频| 国产高清三级在线| 村上凉子中文字幕在线| 久久久久国内视频| 亚洲国产欧美网| 久久婷婷人人爽人人干人人爱| 黄色片一级片一级黄色片| 一二三四在线观看免费中文在| 国产精品九九99| 蜜桃久久精品国产亚洲av| 国内少妇人妻偷人精品xxx网站 | 国产激情偷乱视频一区二区| 亚洲国产看品久久| 国产成人av教育| 在线a可以看的网站| 国产精品自产拍在线观看55亚洲| 国产高清videossex| 男人舔女人的私密视频| 嫩草影院入口| 黑人操中国人逼视频| 国产私拍福利视频在线观看| 亚洲第一欧美日韩一区二区三区| 日韩免费av在线播放| 又大又爽又粗| 少妇熟女aⅴ在线视频| 亚洲乱码一区二区免费版| 在线观看日韩欧美| 国产午夜精品论理片| 国产欧美日韩精品亚洲av| 99久久综合精品五月天人人| 97碰自拍视频| 一本一本综合久久| 欧美日韩精品网址| or卡值多少钱| 国产探花在线观看一区二区| av欧美777| www国产在线视频色| 97碰自拍视频| 美女cb高潮喷水在线观看 | 亚洲无线在线观看| 亚洲人成网站高清观看| 黄色女人牲交| 校园春色视频在线观看| 亚洲乱码一区二区免费版| 久99久视频精品免费| 亚洲国产高清在线一区二区三| 欧美在线一区亚洲| 久久精品国产综合久久久| 女生性感内裤真人,穿戴方法视频| 夜夜躁狠狠躁天天躁| 久久99热这里只有精品18| 国产v大片淫在线免费观看| 午夜福利在线在线| 麻豆成人午夜福利视频| 国产亚洲精品久久久久久毛片| 亚洲欧美日韩卡通动漫| 久久精品亚洲精品国产色婷小说| 亚洲无线在线观看| 一级毛片精品| 一二三四社区在线视频社区8| 制服人妻中文乱码| 在线永久观看黄色视频| 精品99又大又爽又粗少妇毛片 | 亚洲国产精品sss在线观看| 国产成人aa在线观看| 国产99白浆流出| 麻豆av在线久日| 在线观看66精品国产| 亚洲欧美一区二区三区黑人| 91在线观看av| 欧美性猛交黑人性爽| 亚洲欧美日韩高清在线视频| 久久精品国产清高在天天线| 成人三级做爰电影| 欧美日本亚洲视频在线播放| bbb黄色大片| 此物有八面人人有两片| 欧美日韩一级在线毛片| 悠悠久久av| 亚洲五月婷婷丁香| 在线观看66精品国产| 夜夜爽天天搞| 一二三四在线观看免费中文在| 亚洲aⅴ乱码一区二区在线播放| 国产精品亚洲美女久久久| 国产精品九九99| 99热精品在线国产| АⅤ资源中文在线天堂| 亚洲人成网站高清观看| 国产免费av片在线观看野外av| 婷婷丁香在线五月| 成人三级黄色视频| 欧美色欧美亚洲另类二区| 大型黄色视频在线免费观看| 午夜精品久久久久久毛片777| 亚洲一区二区三区色噜噜| 午夜福利视频1000在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲五月天丁香| 国产精品99久久99久久久不卡| 亚洲七黄色美女视频| 亚洲最大成人中文| 18禁美女被吸乳视频| 韩国av一区二区三区四区| 国产一区二区三区在线臀色熟女| 国产欧美日韩精品一区二区| 久久精品综合一区二区三区| 真人做人爱边吃奶动态| 床上黄色一级片| 中文字幕久久专区| 免费电影在线观看免费观看| 亚洲欧洲精品一区二区精品久久久| 国产精品av视频在线免费观看| 色av中文字幕| 一区福利在线观看| 在线播放国产精品三级| 免费av毛片视频| 午夜久久久久精精品| 1024香蕉在线观看| 精品久久久久久久久久免费视频| 天堂动漫精品| 欧美最黄视频在线播放免费| 国产av不卡久久| 精品国内亚洲2022精品成人| 19禁男女啪啪无遮挡网站| 成人一区二区视频在线观看| 国产黄a三级三级三级人| 国产欧美日韩精品一区二区| 综合色av麻豆| 熟妇人妻久久中文字幕3abv| 亚洲精品中文字幕一二三四区| 国产精品乱码一区二三区的特点| а√天堂www在线а√下载| 搞女人的毛片| 狠狠狠狠99中文字幕| 成人高潮视频无遮挡免费网站| 久久久久久久精品吃奶| 日韩精品中文字幕看吧| 免费看十八禁软件| 九九在线视频观看精品| 国产精华一区二区三区| 一进一出抽搐gif免费好疼| 国产成人aa在线观看| 精品电影一区二区在线| 亚洲欧美激情综合另类| 久久久国产欧美日韩av| 色综合欧美亚洲国产小说| 国产高清激情床上av| 国产人伦9x9x在线观看| 天堂动漫精品| 国产激情欧美一区二区| 亚洲中文av在线| 国产熟女xx| cao死你这个sao货| 免费无遮挡裸体视频| 母亲3免费完整高清在线观看| 麻豆成人av在线观看| 俺也久久电影网| 国产极品精品免费视频能看的| 午夜精品久久久久久毛片777| www.自偷自拍.com| 国产视频内射| 男女下面进入的视频免费午夜| 天堂影院成人在线观看| 国产91精品成人一区二区三区| 精品一区二区三区av网在线观看| 中文字幕人妻丝袜一区二区| 欧美一级a爱片免费观看看| 日韩三级视频一区二区三区| 99久久99久久久精品蜜桃| 国产精品av久久久久免费| 激情在线观看视频在线高清|