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

    一種LSTM模型預(yù)測BC值的空間碎片無控再入預(yù)報方法

    2022-02-01 13:29:26蒯家偉趙柯昕孫立剛廖名傳
    宇航學(xué)報 2022年12期
    關(guān)鍵詞:時刻神經(jīng)網(wǎng)絡(luò)誤差

    蒯家偉,趙柯昕,孫立剛,廖名傳

    (1. 南京航空航天大學(xué)航天學(xué)院,南京 211106; 2. 北京電子工程總體研究所,北京 100854;3. 北京遙測技術(shù)研究所,北京 100094; 4. 中國科學(xué)院微小衛(wèi)星創(chuàng)新研究院,上海 201304)

    0 引 言

    低地球軌道(LEO)空間碎片在大氣阻力耗散作用下,軌道能量降低并再入大氣層,進(jìn)而解體隕落??臻g碎片再入有利于降低駐空物體(RSO)碰撞風(fēng)險,維持空間環(huán)境物體數(shù)量穩(wěn)定,由于其中存在大型碎片如運(yùn)載火箭上面級、燃料儲箱等,可能因再入大氣層時燃燒不充分而撞擊地面,對地面人類財產(chǎn)等帶來安全隱患[1],因而對空間碎片的再入時刻預(yù)報對空間環(huán)境的治理與地面安全的保障有著重要意義,為了降低這類事故風(fēng)險需要對空間碎片再入時刻進(jìn)行高精度的預(yù)測?,F(xiàn)階段空間碎片再入時刻預(yù)報誤差主要來源于大氣阻力的估計誤差。大氣阻力攝動fdrag主要與碎片的大氣阻力系數(shù)Cd,面積質(zhì)量比A/m,大氣密度ρ和相對于大氣速度矢量v有關(guān),其表達(dá)式如式(1)所示:

    (1)

    式中:v=|v|,其中大氣密度受太陽和地磁活動影響,存在較強(qiáng)的短期時變特征,較長時間預(yù)報導(dǎo)致軌道誤差較大,因此空間碎片再入時刻精確預(yù)測一直是空間碎片研究的重點和難點[2]。

    國際機(jī)構(gòu)間空間碎片協(xié)調(diào)委員會(IADC)自1998年以來為了提高空間碎片再入預(yù)報水平,組織開展了27次再入隕落預(yù)報聯(lián)測活動,從歷史結(jié)果看總體統(tǒng)計精度水平在30%~10%之間[3]。在預(yù)測方法方面,通過實時衛(wèi)星數(shù)據(jù)矯正密度模型提高大氣密度預(yù)測值[4],由高精度大氣模型得到較精確的相對大氣速度[5]等方法可一定程度上提高短期預(yù)報精度,但因空間碎片外型、材質(zhì)、運(yùn)動復(fù)雜多樣,這類方法難以適用??臻g碎片大氣阻力系數(shù)和面質(zhì)比參數(shù)主要受大氣密度、材質(zhì)、碎片姿態(tài)的影響,直接測量阻力系數(shù)和面質(zhì)比參數(shù)十分困難。因此在難以確定面積質(zhì)量的情形下,通常軌道預(yù)報將阻力系數(shù)和面質(zhì)比耦合成一個系數(shù),稱為彈道系數(shù)(Ballistic coefficient,BC),本文在公式中以β表示。

    (2)

    BC的值可以在具有雷達(dá)、光學(xué)等觀測條件下通過精密定軌進(jìn)行估算,從而提高再入預(yù)測精度。而當(dāng)信息不充分時,要提升空間碎片再入彈道系數(shù)估算精度是很困難的。對于無觀測設(shè)備、或者觀測能力不足的情況下,BC的估算還可以通過機(jī)構(gòu)(space-track.org)發(fā)布的兩行根數(shù)進(jìn)行擬合預(yù)測,這是當(dāng)前IADC國際聯(lián)測中各機(jī)構(gòu)基于space-track.org數(shù)據(jù)的主要方法。以Pardini和Anselmo提出基于BC序貫修正預(yù)測BC值的再入時刻預(yù)測方法為例,其BC值的不確定度為±20%[6],再入時刻預(yù)測偏差接近20%。

    近年來隨著人工智能(AI)技術(shù)的發(fā)展,深度學(xué)習(xí)(DL)如深度神經(jīng)網(wǎng)絡(luò)(DNN)、卷積神經(jīng)網(wǎng)絡(luò)(CNNs)、長短周期記憶神經(jīng)網(wǎng)絡(luò)(Long short term memory neural networks,LSTM)等方法被廣泛應(yīng)用于參數(shù)預(yù)測、建模。航天任務(wù)與空間科學(xué)中人工智能算法研究近年來也是熱門課題,涉及地磁活動性指數(shù)、熱層大氣密度、電離層電子含量等在內(nèi)的多種參量都使用深度學(xué)習(xí)法進(jìn)行了較高精度地預(yù)測或建模[7],包括與本文密切相關(guān)的對大氣密度存在重要影響的參數(shù)地磁指數(shù)的預(yù)測,Kp指數(shù)[8]和Dst指數(shù)[9]都可以通過神經(jīng)網(wǎng)絡(luò)算法進(jìn)行預(yù)報;Chen等[10]開發(fā)出R-DCGAN模型實現(xiàn)了TEC Map缺失數(shù)據(jù)的補(bǔ)全,以及全球預(yù)測和未來預(yù)測,受此啟發(fā),Pan等[11]設(shè)計開發(fā)了DCGAN-PB;隨后,Liu等[12]和Tang等[13]也分別使用長短周期記憶神經(jīng)網(wǎng)絡(luò)預(yù)測了常態(tài)及磁暴期間電離層TEC。張君彪等[14]提出了一種基于集合經(jīng)驗?zāi)B(tài)分解和注意力長短時記憶神經(jīng)網(wǎng)絡(luò)的高超聲速滑行飛行器軌跡智能預(yù)測方法,對氣動加速度進(jìn)行實時估計,進(jìn)而智能預(yù)測飛行軌跡。此外在空間物體有控再入以及路徑規(guī)劃方面,深度學(xué)習(xí)的應(yīng)用也取得了較多的研究成果,周大鵬等[15]設(shè)計了一種智能神經(jīng)網(wǎng)絡(luò)自適應(yīng)復(fù)合控制方法,基于誤差反饋學(xué)習(xí)準(zhǔn)則在線更新神經(jīng)網(wǎng)絡(luò)權(quán)重以自適應(yīng)提升空天飛機(jī)再入精度。相比于傳統(tǒng)方法,深度學(xué)習(xí)可以從大量觀測樣本中提取有價值的特征,并利用多神經(jīng)網(wǎng)絡(luò)層重現(xiàn)相關(guān)特征,從而實現(xiàn)狀態(tài)預(yù)測。

    本工作提出了一種利用兩行軌道根數(shù)(Two line element, TLE)數(shù)據(jù)獲得精確彈道系數(shù)BC的深度學(xué)習(xí)方法,采用長短期記憶神經(jīng)網(wǎng)絡(luò)LSTM模型對其進(jìn)行預(yù)測,繼而利用高精度動力學(xué)模型計算再入時刻。本文首先介紹了利用簡化通用攝動模型(SGP4)從TLE數(shù)據(jù)中獲取精度較高BC值的算法,并實現(xiàn)其計算過程與結(jié)果,建立BC值訓(xùn)練集,然后詳細(xì)闡述了LSTM預(yù)測模型結(jié)構(gòu)的再入時刻預(yù)報模型的建立過程,最后對訓(xùn)練好的模型性能進(jìn)行驗證。使用LSTM模型剔除了BC值中存在的誤差因素,可以在空間碎片再入前中長期獲得比較準(zhǔn)確的再入時刻。不僅如此,這種方法在一次學(xué)習(xí)之后可以快速地對大批量空間碎片再入時刻進(jìn)行動態(tài)跟蹤預(yù)測,提高預(yù)測精度,并可大幅減少積分計算量,顯著提升運(yùn)算效率。

    1 彈道系數(shù)的計算

    1.1 TLE數(shù)據(jù)的預(yù)處理

    TLE數(shù)據(jù)中包含了空間碎片的軌道根數(shù)、平運(yùn)動以及和大氣阻力相關(guān)的B*數(shù)據(jù)。在大批量使用TLE數(shù)據(jù)時,大概率會遇到空間碎片TLE數(shù)據(jù)異常情況,包括重復(fù)歷元、數(shù)據(jù)格式錯亂、編號異常等。因此需要對整個TLE數(shù)據(jù)進(jìn)行預(yù)處理,以去除異常和錯誤的數(shù)據(jù)[16]。預(yù)處理過程主要包括去除重復(fù)TLE數(shù)據(jù)、平均運(yùn)動項異常數(shù)據(jù)、傾角異常、B*為負(fù)值的數(shù)據(jù)等,由于TLE格式可讀性問題,預(yù)處理過后的TLE數(shù)據(jù)依然會存在某些異常情況,會明顯反應(yīng)在后續(xù)計算結(jié)果中。

    1.2 彈道系數(shù)BC估計方法

    以Saunders等[17]和Gondelach等[18]提出利用TLE數(shù)據(jù)修正BC值的方法作為本文研究的比對基準(zhǔn)。其基本原理是通過利用TLE配套的SGP4模型從TLE數(shù)據(jù)中獲得碎片的軌道狀態(tài)、物理參數(shù)(特別是B*數(shù)據(jù))等信息,建立連續(xù)TLE數(shù)據(jù)間的平運(yùn)動參數(shù)變率迭代估計任意時刻的BC值。

    由于TLE數(shù)據(jù)中不包含短周期項的變化,數(shù)據(jù)中半長軸的變化可以假設(shè)為僅受大氣阻力影響下的長期項變化??紤]到太陽光壓導(dǎo)致的半長軸長周期變化可能包括在TLE數(shù)據(jù)中,但由于太陽光壓與低軌衛(wèi)星再入過程中大氣阻力攝動相比,其影響要小得多,因此本文忽略光壓的影響。由此可假設(shè)基于TLE數(shù)據(jù)平運(yùn)動參數(shù)計算得出的結(jié)果與由軌道外推計算結(jié)果之間的任何差異都是因為大氣阻力彈道系數(shù)BC值的不正確猜測引起的。利用連續(xù)TLE數(shù)據(jù)計算BC值的過程如下:

    1)利用TLE數(shù)據(jù)中的平運(yùn)動參數(shù)計算兩個TLE數(shù)據(jù)的半長軸的變化量Δa。

    (3)

    Δa=a2-a1

    (4)

    式中:a為軌道半長徑;μ為地球引力常數(shù);n0為空間物體軌道平均軌道運(yùn)動速度,也就是TLE的平運(yùn)動項。a1與a2分布表示前后兩次TLE數(shù)據(jù)的半長徑。

    2)使用高精度的軌道動力學(xué)模型進(jìn)行外推,計算兩個TLE時刻(tTLE1,tTLE2)內(nèi)僅由大氣阻力引起的軌道半長軸的變化。其中p為軌道半通徑,θ為真近點角,e為偏心率,p=a(1-e2),fRdrag,fTdrag分別是由阻力在徑向和橫向分量產(chǎn)生的加速度。

    (5)

    (6)

    3)利用TLE數(shù)據(jù)中的B*參數(shù)確定BC的初值。B*參數(shù)是SGP4模型中的近似反應(yīng)阻力的系數(shù),可以通過(7)式計算得到BC的初值β1。第二個TLE根數(shù)的BC值β2可以通過第一個BC的值β1和Δaprop之間的線性關(guān)系計算得到。當(dāng)Δaprop小于10-4千米時,BC值收斂。如果在后續(xù)計算過程中出現(xiàn)了過早再入的情況,可以將β1的值設(shè)置為原來值的90%,但是需要更多次的迭代方能收斂。

    (7)

    4)對BC值進(jìn)行迭代計算,直到迭代收斂。βn是迭代過程中第n個BC值。

    (8)

    式中:Δadiff=Δa-Δaprop。

    在默認(rèn)情況下,使用正向傳播來進(jìn)行BC值的估算,取最早TLE時刻的狀態(tài),將它外推到最新TLE數(shù)據(jù)時刻的狀態(tài)。也可以使用反向傳播,從最新的TLE時刻的狀態(tài)向后傳播到前一個TLE時刻的狀態(tài),反向傳播可以防止預(yù)測時空間碎片過早的再入。在臨近再入時刻估計時十分有用,因為不準(zhǔn)確的BC估計值可能導(dǎo)致再入預(yù)報異常。以上估計過程其本質(zhì)原理類似Newton-Raphson方法。如果TLE數(shù)據(jù)充分,且TLE數(shù)據(jù)在時序上連續(xù)且間隔不大,不失為一種較優(yōu)的方法,然后對于中短期預(yù)報來說,由于其數(shù)據(jù)使用率較高,如果其中出現(xiàn)異常點,其預(yù)測值會突然拉偏,導(dǎo)致預(yù)報誤差突然增大??臻g物體的軌道是攝動情況下連續(xù)可微的擬二次曲線,因此有必要在預(yù)測過程中考慮短期的數(shù)據(jù)情況的變化,同時又能兼顧長期的動力學(xué)演化趨勢,因此本文參考了長短期記憶網(wǎng)絡(luò)預(yù)測模型LSTM在空間任務(wù)中的研究結(jié)論,將其應(yīng)用于彈道系數(shù)BC值的預(yù)測。

    2 長短期記憶神經(jīng)網(wǎng)絡(luò)LSTM模型原理

    為了解決原有的循環(huán)神經(jīng)網(wǎng)絡(luò)(RNN)存在梯度爆炸或梯度消失的問題[19],Hochreiter等[20]于1997年提出了長短期記憶神經(jīng)網(wǎng)絡(luò)LSTM,并被Graves等[21]進(jìn)行了改進(jìn)和推廣。LSTM神經(jīng)網(wǎng)絡(luò)經(jīng)過了特殊的設(shè)計,使得其適合在重要事件之間存在很長時間滯后時從經(jīng)驗中學(xué)習(xí),以對時間序列進(jìn)行分類、處理和預(yù)測,目前在多個領(lǐng)域得到了廣泛的應(yīng)用,其單元內(nèi)部包含遺忘門、輸入門和輸出門,具體結(jié)構(gòu)如圖1所示。

    圖1 LSTM單元內(nèi)部結(jié)構(gòu)Fig.1 Internal structure of LSTM unit

    遺忘門用來確定狀態(tài)傳遞的過程中上一時刻的單元狀態(tài)ct-1有多少保留到當(dāng)前時刻的狀態(tài)ct。讀取上一個時刻輸出值ht-1和當(dāng)前時刻輸入的BC值數(shù)據(jù)xt,并利用設(shè)定的連接權(quán)重矩陣Wf和sigmoid函數(shù)σ計算得到經(jīng)過處理的需要保留下來的歷史信息狀態(tài),如式(9)所示,其中bf為遺忘門的偏移值。將會輸出一個0至1之間的數(shù)值給歷史狀態(tài)ct-1,1表示完全保留,0表示完全遺棄。

    ft=σ(Wf[ht-1,xt]+bf)

    (9)

    輸入門利用連接權(quán)重參數(shù)Wi、偏移值bi等參數(shù)來確定當(dāng)前時刻輸入的BC值數(shù)據(jù)xt有多少可以作為輸入門的輸出it,如式(10)所示。再利用計算新狀態(tài)的連接權(quán)重矩陣WC、偏移值bC確定臨時的當(dāng)前時刻新狀態(tài),如式(11)。再利用遺忘門和輸入門的數(shù)據(jù)更新狀態(tài),將上一時刻狀態(tài)ct-1更新到當(dāng)前時刻狀態(tài)ct,如式(12)。

    it=σ(Wi·[ht-1,xt]+bi)

    (10)

    (11)

    (12)

    輸出門利用連接權(quán)重矩陣Wo、偏移值bo確定當(dāng)前時刻的輸出值ht,如式(13)、式(14)所示。

    ot=σ(Wo[ht-1,xt]+bo)

    (13)

    ht=ot·tanh(Ct)

    (14)

    3 再入預(yù)報流程與性能評估

    3.1 再入時刻預(yù)測流程

    首先建立預(yù)測模型的訓(xùn)練數(shù)據(jù)集,別選取近三年再入大氣層的80個和20個空間碎片的TLE數(shù)據(jù)作為預(yù)測模型的訓(xùn)練集和測試集。TLE數(shù)據(jù)和再入時刻都是從北美防空司令部空間目標(biāo)監(jiān)視編目軌道根數(shù)發(fā)布網(wǎng)站(https://www.space-track.org)上獲取的,該網(wǎng)站發(fā)布了空間物體歷史TLE數(shù)據(jù)以及再入日期,確切的再入時刻未給出,因此最終再入時刻僅能精確到天。

    LSTM模型預(yù)測彈道系數(shù)BC值的具體流程如圖2所示。首先對所有TLE數(shù)據(jù)進(jìn)行預(yù)處理,剔除其中誤差較大及錯誤的數(shù)據(jù),利用1.2節(jié)的方法計算出每一TLE時刻的BC值,隨后將得到的BC值隨時間變化數(shù)據(jù)代入LSTM模型,確定模型中的特征參數(shù),由此得到了BC值的精確預(yù)測模型,再將測試集代入預(yù)測模型中,得到對未來BC值的預(yù)測,即可計算出精確的再入時刻,最后對再入時刻進(jìn)行評估。

    圖2 空間碎片再入時刻預(yù)測流程Fig.2 The reentry time prediction scheme of space debris

    3.2 彈道系數(shù)BC的估計

    圖3為編號為32712的空間碎片再入大氣層前200天分別利用1.2節(jié)所述的迭代過程和由TLE數(shù)據(jù)中B*參數(shù)兩種方法計算得到的BC值隨時間的變化曲線。圖中橫軸為再入前的時間,縱軸為BC值的大小。圖示TLE數(shù)據(jù)是經(jīng)過預(yù)處理的,可以看出兩種方法計算的BC值具有相同的變化趨勢,但是有一個偏移量,同時BC值中仍然存在一些異常點,如果進(jìn)一步利用近地點半徑進(jìn)行濾波,這些異常值會被剔除。

    圖3 編號32712空間碎片BC值與再入時刻關(guān)系Fig.3 The relationship between BC value and reentry time of the space debris No. 32712

    利用1.2節(jié)的迭代方法,選取了100個已再入隕落的空間碎片作為目標(biāo)集,計算了目標(biāo)集對應(yīng)的從再入前200天到再入時刻彈道系數(shù)BC的最小值、平均值和最大值。計算結(jié)果表明基于迭代修正的BC值預(yù)測再入時刻精度在10%的置信度為60%~70%,預(yù)測精度在20%的置信度為95%。根據(jù)歐洲航天局的DISCOS數(shù)據(jù)庫提供的信息,彈道系數(shù)BC估計值精度20%的置信度為85%??紤]到大氣密度10~20%誤差因素,表明BC可以利用1.2節(jié)的流程建立訓(xùn)練集。因此以此為參考對象建立訓(xùn)練集樣本,利用LSTM模型預(yù)測BC值,進(jìn)而預(yù)報再入時刻,并評估預(yù)報精度與模型性能。

    3.3 建立預(yù)測BC值的LSTM模型

    LSTM模型的建立主要分3個步驟:首先需要對BC值數(shù)據(jù)集進(jìn)行整理;其次是完成對模型的構(gòu)建,包括確定用于訓(xùn)練和預(yù)測所需要的特征參數(shù)和確定LSTM神經(jīng)網(wǎng)絡(luò)模型的具體參數(shù);最后是利用測試集數(shù)據(jù)代入得到預(yù)測結(jié)果,并對預(yù)測模型準(zhǔn)確性進(jìn)行驗證,如圖4所示。

    圖4 LSTM預(yù)測模型框架Fig.4 The framework of LSTM prediction model

    采用多時刻輸入單一時刻輸出的模式對彈道系數(shù)BC值進(jìn)行預(yù)測,利用40天的TLE數(shù)據(jù)作為預(yù)測模型的輸入,預(yù)測結(jié)果作為輸出。設(shè)置兩層LSTM預(yù)測模型,并使用全連接層對LSTM層進(jìn)行線性化。訓(xùn)練過程中,反向傳播的誤差項是損失函數(shù)對輸出值的倒數(shù),損失函數(shù)是均方根誤差函數(shù),給定各層神經(jīng)元個數(shù)、學(xué)習(xí)率以及訓(xùn)練迭代步數(shù),利用梯度下降算法更新網(wǎng)絡(luò)權(quán)重并獲得最終的隱藏層。隨后將測試集進(jìn)行迭代預(yù)測,并輸出預(yù)測結(jié)果。

    建立彈道系數(shù)預(yù)測模型時需要考慮5個參數(shù),即訓(xùn)練步長、輸出層維數(shù)、隱藏層層數(shù)、隱藏層維數(shù)和輸出變量維數(shù)。訓(xùn)練步長的設(shè)置要考慮該步長包含的信息量是否可以滿足模型訓(xùn)練,若設(shè)置的過小,則不足以使模型在訓(xùn)練中學(xué)習(xí)BC數(shù)據(jù)的變化規(guī)律,同時考慮到模型訓(xùn)練難易程度,設(shè)置過大可能使得數(shù)據(jù)序列中包含過多無用數(shù)據(jù),導(dǎo)致預(yù)測性能下降并增加訓(xùn)練時間。本工作測模型中訓(xùn)練步長設(shè)置為20步。

    針對于本文的隕落時間預(yù)測場景,輸入層維數(shù)的設(shè)置根據(jù)雙行軌道根數(shù)中軌道和動力學(xué)特征數(shù)據(jù)來確定,在本模型中使用TLE數(shù)據(jù),所以輸入層維數(shù)設(shè)置為1。多隱藏層的網(wǎng)絡(luò)從理論上來說有足夠的表達(dá)能力可以擬合任何復(fù)雜的函數(shù),但網(wǎng)絡(luò)結(jié)構(gòu)越復(fù)雜越難以訓(xùn)練,并且可能出現(xiàn)過擬合的情況,即在訓(xùn)練數(shù)據(jù)上表現(xiàn)很好,但泛化能力差。隱藏層維數(shù)即為隱藏層神經(jīng)元個數(shù),對于較復(fù)雜、非線性程度大的序列,設(shè)置的神經(jīng)元個數(shù)越多擬合效果越好,但花費(fèi)的時間也越多,同時可能出現(xiàn)過擬合的現(xiàn)象。本模型中隱藏層由兩層LSTM層和全連接層構(gòu)成,第二層的LSTM層神經(jīng)元個數(shù)是第一層個數(shù)的兩倍,第三層的全連接層神經(jīng)元個數(shù)和第一層LSTM層個數(shù)相等,由于本文主要是利用算法,算法具體細(xì)節(jié)不作贅述。

    3.4 LSTM預(yù)測模型的性能評估

    使用平均絕對百分比誤差(MAPE)作為評估參數(shù)與預(yù)測模型輸出值之間的預(yù)測準(zhǔn)確性度量,用以量化訓(xùn)練后預(yù)測模型的性能。MAPE定義為測試數(shù)據(jù)集中每個數(shù)據(jù)點上的BC預(yù)測值βpre與計算值βcal之間的真實誤差之和與計算值βcal之和的比值,下式中以ξ標(biāo)志,其數(shù)學(xué)表達(dá)式為:

    (15)

    式中:n是參與訓(xùn)練的有效數(shù)據(jù)數(shù)量。利用預(yù)測模型預(yù)測出結(jié)果越精確,式(15)將越逼近其下邊界0,MAPE的數(shù)值越低說明預(yù)測模型的性能越好。這里選擇的度量標(biāo)準(zhǔn)與預(yù)測模型的目標(biāo)一致,該目標(biāo)是通過最大程度減少訓(xùn)練數(shù)據(jù)中的總誤差,使所有數(shù)據(jù)點盡可能平滑。

    3.5 LSTM模型預(yù)測BC值效能分析

    采用編號20230的空間碎片,利用其再入前200天至160天之間40天的TLE數(shù)據(jù)作為預(yù)測模型的輸入,預(yù)測再入前160天的BC值隨時間的變化關(guān)系。

    圖5橫軸表示再入預(yù)報所用TLE根數(shù)時間與www.space-track.org公布的再入時間的間隔,即當(dāng)前歷元距離再入日期的天數(shù),從再入前160天至再入時刻為止,縱軸表示了BC值的大小。其中橙色折線的中心點表示了10天內(nèi)的BC計算值的平均值,藍(lán)色折線的中心點表示了10天內(nèi)的BC預(yù)測值的平均值,誤差線表示了10天內(nèi)的最大值和最小值。總數(shù)據(jù)統(tǒng)計預(yù)測模型的性能指標(biāo)MAPE達(dá)到13.9%,在某些時間段,MAPE更低,其中有一些未能預(yù)處理掉的異常數(shù)據(jù)點也清晰地反映在圖上,如果剔除一些異常點,預(yù)測結(jié)果的MAPE值會更低。數(shù)據(jù)結(jié)果表明訓(xùn)練好的預(yù)測模型可以較好地學(xué)習(xí)BC值變化趨勢,表明本文建立的LSTM預(yù)測模型是有效的。

    圖5 編號20230的空間碎片再入前160天的BC預(yù)測值和計算值與預(yù)測時間的關(guān)系Fig.5 The relationship between BC predicted/calculated values and predicted times for 160 days prior to reentry of the space debris No.20230

    3.6 再入預(yù)報能分析

    BC值預(yù)測的目標(biāo)是為了更好地預(yù)報空間碎片再入時刻,或者是在軌壽命。在獲得BC預(yù)測值之后,基于最新的TLE根數(shù),考慮地球非球形攝動、大氣阻力攝動、太陽光壓攝動和日月三體引力攝動,利用RKF7/8階數(shù)值積分進(jìn)行軌道外推,外推截止軌道高度設(shè)定80 km即認(rèn)為目標(biāo)已經(jīng)發(fā)生了再入,采集此時的歷元作為再入時刻。

    為了評估再入時刻的預(yù)報精度,使用時間的相對誤差估計的精度,再入預(yù)報精度ε定義為實際預(yù)測時間誤差與最后所用軌道根數(shù)歷元距離隕落日期的時間之比值。

    (16)

    式中:tpre是預(yù)測的再入日期,tact是space-track.org網(wǎng)站上公布的再入日期,tTLElast是最后用于預(yù)測再入時刻的TLE軌道根數(shù)歷元。

    以編號42727目標(biāo)為例,取再入前160天內(nèi)所有有效的TLE數(shù)據(jù),利用當(dāng)天的TLE數(shù)據(jù)計算出衛(wèi)星的狀態(tài),并預(yù)測出從TLE數(shù)據(jù)時刻到再入時刻的BC變化趨勢。利用衛(wèi)星軌道參數(shù)和BC變化計算出衛(wèi)星再入時刻并與space-track.org網(wǎng)站給出的真實時刻計算再入時刻誤差,如圖6所示??梢园l(fā)現(xiàn),利用本文方法得到平均預(yù)報誤差為7.3%,再入前160天時預(yù)測誤差在4%以下,再入前30天時預(yù)測誤差在10%以下。但是再入前20天預(yù)測的效果沒有長期預(yù)測的理想,最大預(yù)測誤差可達(dá)到51%,這是因為space-track.org給出的目標(biāo)再入時刻只精確到天,依據(jù)式(16)的標(biāo)準(zhǔn),越到末端誤差會急劇增大,導(dǎo)致即使預(yù)測是準(zhǔn)確的,再入前10天預(yù)測時每差24小時就帶有近10%的誤差,再入最后2天每一至兩個小時就接近10%誤差。同時,短期預(yù)測誤差率大也可能是因為接近再入大氣層時氣動力快速變化導(dǎo)致無法準(zhǔn)確預(yù)測與計算,因此后續(xù)評估主要對大于30天臨近再入時間的再入預(yù)報有較好的指導(dǎo)意義,但本文同時也給出了小于30天的預(yù)測結(jié)果。

    圖6 編號42727空間碎片再入時刻預(yù)報結(jié)果Fig.6 Re-entry time forecast results of the space debris No.42727

    選取5個已再入物體作為測試目標(biāo),利用本文的方法預(yù)測其再入前160天至再入時刻的BC變化,基于再入前160天,90天和30天的TLE數(shù)據(jù)計算再入時刻,并計算再入時刻誤差。計算結(jié)果如表1所示。

    表1 測試目標(biāo)再入時刻預(yù)報結(jié)果Table 1 Re-entry time forecast results of the test targets

    對選取的90個目標(biāo)于再入前10天、20天、30天、60天、90天和160天時預(yù)測再入時刻并得到累計誤差。所得到的再入預(yù)報時間累積誤差和95%置信區(qū)間結(jié)果如圖7所示。圖7表明,再入前160天時的預(yù)測最準(zhǔn)確,總誤差低于5%,再入前10、20天預(yù)測的結(jié)果不如整體結(jié)果準(zhǔn)確是由于公布的再入時刻精確到天導(dǎo)致??傮w來看,再入前160天,90天,30天分別得到的預(yù)測誤差在5%以下,8%以下和12%以下,算例表明利用長短期記憶神經(jīng)網(wǎng)絡(luò)先預(yù)測出物體BC變化規(guī)律再進(jìn)行再入時刻預(yù)報可以得到較好的空間碎片無控再入預(yù)報時刻。

    圖7 90個目標(biāo)的再入前10、20、30、60、90和160天的再入時刻預(yù)報誤差的累積分布Fig.7 Cumulative distributions of reentry time prediction errors of 90 objects for 10, 20, 30, 60, 90, and 180 days prior to reentry

    4 結(jié) 論

    空間碎片無控再入時刻預(yù)報是一項國際上非常關(guān)注的重要工作,其受到大氣密度、軌道數(shù)據(jù)精度、彈道系數(shù)估計等諸多因素的制約,以至于準(zhǔn)確預(yù)測空間碎片的再入時刻是十分困難的。在已有大氣阻力彈道系數(shù)BC值預(yù)測方法的基礎(chǔ)上,建立了基于長短期記憶神經(jīng)網(wǎng)絡(luò)深度學(xué)習(xí)LSTM模型預(yù)測BC值,進(jìn)而預(yù)測空間碎片再入時刻的方法,獲得了較好的研究結(jié)果。從選取的樣本訓(xùn)練模型以及測試結(jié)果看,在95%的置信度內(nèi),90天預(yù)測精度達(dá)到8%,30天預(yù)測精度達(dá)到12%,考慮到大氣密度模型帶來的不確定影響,結(jié)果表明LSTM模型預(yù)測方法在空間碎片再入預(yù)報具有很好應(yīng)用價值。后續(xù)我們將進(jìn)一步精細(xì)化考慮空間碎片類型,將相近類型的空間碎片進(jìn)行聚類,建立更精確數(shù)據(jù)樣本,提高再入或者在軌壽命預(yù)測精度。

    猜你喜歡
    時刻神經(jīng)網(wǎng)絡(luò)誤差
    冬“傲”時刻
    捕獵時刻
    角接觸球軸承接觸角誤差控制
    哈爾濱軸承(2020年2期)2020-11-06 09:22:26
    Beidou, le système de navigation par satellite compatible et interopérable
    神經(jīng)網(wǎng)絡(luò)抑制無線通信干擾探究
    電子制作(2019年19期)2019-11-23 08:42:00
    壓力容器制造誤差探究
    基于神經(jīng)網(wǎng)絡(luò)的拉矯機(jī)控制模型建立
    九十億分之一的“生死”誤差
    山東青年(2016年2期)2016-02-28 14:25:41
    復(fù)數(shù)神經(jīng)網(wǎng)絡(luò)在基于WiFi的室內(nèi)LBS應(yīng)用
    街拍的歡樂時刻到來了
    很黄的视频免费| 国产av不卡久久| 午夜久久久在线观看| 国产精品永久免费网站| 夜夜夜夜夜久久久久| 日韩高清综合在线| 久久国产亚洲av麻豆专区| 伊人久久大香线蕉亚洲五| 老司机午夜福利在线观看视频| 亚洲激情在线av| 欧美乱色亚洲激情| 亚洲va日本ⅴa欧美va伊人久久| 麻豆成人av在线观看| 日韩欧美 国产精品| 国产午夜福利久久久久久| 国产欧美日韩一区二区三| 91在线观看av| 国产精品自产拍在线观看55亚洲| 国产私拍福利视频在线观看| 国产黄色小视频在线观看| 母亲3免费完整高清在线观看| 欧美午夜高清在线| 人妻丰满熟妇av一区二区三区| 中文字幕人妻熟女乱码| 色精品久久人妻99蜜桃| 亚洲欧美日韩无卡精品| 51午夜福利影视在线观看| 三级毛片av免费| 法律面前人人平等表现在哪些方面| 欧洲精品卡2卡3卡4卡5卡区| 成人手机av| 精品一区二区三区视频在线观看免费| 淫妇啪啪啪对白视频| 婷婷丁香在线五月| 1024手机看黄色片| 亚洲av电影不卡..在线观看| 婷婷亚洲欧美| 国产爱豆传媒在线观看 | 久久久久久久久久黄片| 国产精品野战在线观看| 精品久久久久久久久久免费视频| 亚洲国产欧美一区二区综合| 久久人妻福利社区极品人妻图片| 一个人观看的视频www高清免费观看 | 国产成人精品久久二区二区91| 亚洲精品色激情综合| 搡老岳熟女国产| 成人av一区二区三区在线看| 国产一区在线观看成人免费| 国产一区在线观看成人免费| 成人av一区二区三区在线看| 两性夫妻黄色片| 免费电影在线观看免费观看| 亚洲国产高清在线一区二区三 | 精品国产一区二区三区四区第35| 午夜福利在线观看吧| 熟女少妇亚洲综合色aaa.| 看免费av毛片| 怎么达到女性高潮| 久久这里只有精品19| 午夜激情av网站| 中文字幕av电影在线播放| 欧美日韩福利视频一区二区| 久久久国产成人免费| 欧美一级a爱片免费观看看 | 国产精品免费一区二区三区在线| 国产亚洲av嫩草精品影院| 国产片内射在线| 国产精品国产高清国产av| 一进一出抽搐gif免费好疼| 校园春色视频在线观看| 日韩国内少妇激情av| 国产区一区二久久| 99久久无色码亚洲精品果冻| 国产精华一区二区三区| √禁漫天堂资源中文www| 成人18禁在线播放| 欧美乱妇无乱码| 男人舔女人的私密视频| 麻豆一二三区av精品| 99久久国产精品久久久| 真人一进一出gif抽搐免费| 久久国产精品人妻蜜桃| 亚洲国产高清在线一区二区三 | 夜夜夜夜夜久久久久| 9191精品国产免费久久| 久热爱精品视频在线9| 亚洲一区二区三区不卡视频| 99精品久久久久人妻精品| 香蕉丝袜av| 国产欧美日韩一区二区三| av中文乱码字幕在线| 人妻丰满熟妇av一区二区三区| 国产成+人综合+亚洲专区| 国产高清视频在线播放一区| 亚洲片人在线观看| av欧美777| 99re在线观看精品视频| 亚洲av成人一区二区三| 桃色一区二区三区在线观看| 久久久国产成人免费| 免费看日本二区| 青草久久国产| 国产精品美女特级片免费视频播放器 | 免费观看人在逋| 亚洲成a人片在线一区二区| √禁漫天堂资源中文www| 亚洲avbb在线观看| 中文字幕人妻丝袜一区二区| 老熟妇仑乱视频hdxx| 黑人操中国人逼视频| 国产精品野战在线观看| 手机成人av网站| 国产爱豆传媒在线观看 | √禁漫天堂资源中文www| 久久久久久大精品| 少妇被粗大的猛进出69影院| 久久精品国产亚洲av香蕉五月| 国产私拍福利视频在线观看| 国产伦人伦偷精品视频| 美女高潮喷水抽搐中文字幕| 最近在线观看免费完整版| 日日夜夜操网爽| 欧美又色又爽又黄视频| 51午夜福利影视在线观看| 国产日本99.免费观看| 国产精品久久电影中文字幕| 亚洲五月婷婷丁香| 长腿黑丝高跟| 亚洲欧美日韩无卡精品| www日本黄色视频网| 国产亚洲av高清不卡| 午夜福利18| 欧美日韩亚洲综合一区二区三区_| 欧美黄色淫秽网站| 日日摸夜夜添夜夜添小说| 好看av亚洲va欧美ⅴa在| 精品高清国产在线一区| 精品久久久久久久末码| 亚洲av日韩精品久久久久久密| 欧美zozozo另类| 亚洲精品久久成人aⅴ小说| 久久国产精品男人的天堂亚洲| 久久久久久久精品吃奶| 欧美久久黑人一区二区| 美国免费a级毛片| 欧美激情久久久久久爽电影| 亚洲天堂国产精品一区在线| 亚洲成av人片免费观看| tocl精华| 中文字幕精品亚洲无线码一区 | 很黄的视频免费| 亚洲最大成人中文| 啦啦啦 在线观看视频| 久久这里只有精品19| 色婷婷久久久亚洲欧美| 日韩免费av在线播放| 亚洲激情在线av| 成人三级黄色视频| 国产精品1区2区在线观看.| 亚洲在线自拍视频| 国产亚洲精品一区二区www| 最新美女视频免费是黄的| 欧美最黄视频在线播放免费| 黄色成人免费大全| 听说在线观看完整版免费高清| 可以免费在线观看a视频的电影网站| 精品久久久久久久人妻蜜臀av| 国产精品乱码一区二三区的特点| 亚洲一区二区三区不卡视频| 性欧美人与动物交配| 久9热在线精品视频| 99热这里只有精品一区 | 人人妻,人人澡人人爽秒播| 亚洲欧美一区二区三区黑人| www国产在线视频色| 韩国av一区二区三区四区| 又黄又爽又免费观看的视频| 曰老女人黄片| 国产人伦9x9x在线观看| 国产精品亚洲av一区麻豆| 狂野欧美激情性xxxx| 国产成人一区二区三区免费视频网站| 国产精品野战在线观看| 在线永久观看黄色视频| 日韩一卡2卡3卡4卡2021年| 亚洲av片天天在线观看| netflix在线观看网站| 黄色成人免费大全| 在线视频色国产色| 欧美一区二区精品小视频在线| 亚洲国产高清在线一区二区三 | 亚洲一区二区三区色噜噜| 精品国产亚洲在线| 国产精品乱码一区二三区的特点| 国产高清视频在线播放一区| 母亲3免费完整高清在线观看| 成人一区二区视频在线观看| 亚洲国产欧美日韩在线播放| 国产免费av片在线观看野外av| 久久精品aⅴ一区二区三区四区| 熟女电影av网| 手机成人av网站| 亚洲全国av大片| 99精品欧美一区二区三区四区| 久久精品国产亚洲av香蕉五月| 久久久久久久精品吃奶| 91在线观看av| 国产黄a三级三级三级人| 国产成人一区二区三区免费视频网站| 淫妇啪啪啪对白视频| 天天一区二区日本电影三级| 亚洲一区二区三区不卡视频| 啦啦啦韩国在线观看视频| 看黄色毛片网站| 国产精品98久久久久久宅男小说| 国产99白浆流出| 国产精品,欧美在线| 国产精品美女特级片免费视频播放器 | www.www免费av| 午夜成年电影在线免费观看| 中文字幕久久专区| 少妇熟女aⅴ在线视频| 美国免费a级毛片| 欧美日韩精品网址| 天天躁夜夜躁狠狠躁躁| 中文在线观看免费www的网站 | 国产亚洲精品av在线| 久久香蕉精品热| 免费看a级黄色片| 日韩免费av在线播放| 中文字幕人成人乱码亚洲影| 他把我摸到了高潮在线观看| 哪里可以看免费的av片| 午夜福利免费观看在线| 久久午夜亚洲精品久久| 亚洲最大成人中文| 俄罗斯特黄特色一大片| 丁香六月欧美| 国产欧美日韩一区二区三| 午夜福利视频1000在线观看| √禁漫天堂资源中文www| 久久久久国产精品人妻aⅴ院| 日本熟妇午夜| 亚洲色图av天堂| 午夜视频精品福利| 香蕉av资源在线| 亚洲精品色激情综合| 久久人妻福利社区极品人妻图片| 丁香六月欧美| 国产私拍福利视频在线观看| 日韩欧美在线二视频| av天堂在线播放| 亚洲成av人片免费观看| 又黄又爽又免费观看的视频| 黄片小视频在线播放| 麻豆av在线久日| 狂野欧美激情性xxxx| e午夜精品久久久久久久| 国产熟女午夜一区二区三区| 国产精品 欧美亚洲| 国产亚洲精品av在线| 好男人电影高清在线观看| 丝袜美腿诱惑在线| 亚洲在线自拍视频| 午夜福利高清视频| 欧美激情极品国产一区二区三区| 国产精品久久视频播放| 久久精品91无色码中文字幕| 亚洲精华国产精华精| 后天国语完整版免费观看| 观看免费一级毛片| 女性被躁到高潮视频| 黄色女人牲交| 真人一进一出gif抽搐免费| 亚洲精品av麻豆狂野| 国产男靠女视频免费网站| 香蕉丝袜av| 免费在线观看成人毛片| 很黄的视频免费| 色综合欧美亚洲国产小说| 亚洲av成人一区二区三| 一级a爱视频在线免费观看| 18禁观看日本| 日本a在线网址| 天堂影院成人在线观看| 老司机在亚洲福利影院| 97碰自拍视频| 中文字幕人成人乱码亚洲影| 免费电影在线观看免费观看| 免费看美女性在线毛片视频| 欧美丝袜亚洲另类 | 午夜福利在线在线| 欧美性猛交黑人性爽| 亚洲人成电影免费在线| 欧美国产日韩亚洲一区| 欧美中文综合在线视频| 可以在线观看毛片的网站| 在线观看免费视频日本深夜| 丁香欧美五月| 黄频高清免费视频| 成人18禁在线播放| 999精品在线视频| 国产黄a三级三级三级人| 免费在线观看亚洲国产| 欧美激情久久久久久爽电影| 日韩欧美一区二区三区在线观看| 母亲3免费完整高清在线观看| 窝窝影院91人妻| 国产不卡一卡二| av超薄肉色丝袜交足视频| 成人欧美大片| 日韩免费av在线播放| 久久久精品欧美日韩精品| 欧美成人性av电影在线观看| 给我免费播放毛片高清在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 精品久久久久久久久久免费视频| 国产精品一区二区精品视频观看| 久久精品夜夜夜夜夜久久蜜豆 | 久久精品人妻少妇| 精品人妻1区二区| 欧美午夜高清在线| 国产高清有码在线观看视频 | 亚洲成av片中文字幕在线观看| 中文字幕人妻熟女乱码| 国产一区二区三区视频了| 一级片免费观看大全| 国产成人av激情在线播放| 午夜福利在线在线| 日日干狠狠操夜夜爽| 成人18禁在线播放| 国产99久久九九免费精品| 女性被躁到高潮视频| 人人妻,人人澡人人爽秒播| 国内毛片毛片毛片毛片毛片| bbb黄色大片| 亚洲国产欧美网| 日韩 欧美 亚洲 中文字幕| 非洲黑人性xxxx精品又粗又长| 在线观看www视频免费| 又大又爽又粗| 久热爱精品视频在线9| 欧美在线一区亚洲| 日日夜夜操网爽| 欧美在线黄色| 国产99久久九九免费精品| 国产在线精品亚洲第一网站| 观看免费一级毛片| 无限看片的www在线观看| 国产精品影院久久| 国产av一区二区精品久久| 久久亚洲真实| 黄片播放在线免费| 国产91精品成人一区二区三区| www.熟女人妻精品国产| 男人舔女人的私密视频| cao死你这个sao货| 中文在线观看免费www的网站 | bbb黄色大片| 国产精品久久电影中文字幕| 日韩精品中文字幕看吧| avwww免费| 午夜精品在线福利| 91在线观看av| 色综合亚洲欧美另类图片| 在线观看免费日韩欧美大片| 午夜久久久在线观看| 日日夜夜操网爽| 色尼玛亚洲综合影院| 成年版毛片免费区| 久久精品国产亚洲av香蕉五月| 一本一本综合久久| 亚洲一区中文字幕在线| 亚洲av第一区精品v没综合| 精品久久久久久久人妻蜜臀av| 丰满的人妻完整版| 国产成人一区二区三区免费视频网站| 欧美日本视频| 一边摸一边做爽爽视频免费| 午夜福利18| 极品教师在线免费播放| 国产片内射在线| 999精品在线视频| 少妇裸体淫交视频免费看高清 | 国产成人精品久久二区二区免费| 不卡一级毛片| 香蕉丝袜av| 亚洲欧美日韩高清在线视频| 怎么达到女性高潮| 欧美日韩乱码在线| 欧美黑人巨大hd| 欧美色视频一区免费| 51午夜福利影视在线观看| www国产在线视频色| 国内少妇人妻偷人精品xxx网站 | 国产熟女xx| 日韩欧美三级三区| 国产久久久一区二区三区| 精品人妻1区二区| 在线观看免费午夜福利视频| av免费在线观看网站| 少妇裸体淫交视频免费看高清 | 国产日本99.免费观看| 亚洲 欧美 日韩 在线 免费| 国产高清有码在线观看视频 | 国产熟女午夜一区二区三区| 亚洲精品在线美女| 大香蕉久久成人网| 不卡av一区二区三区| 老汉色av国产亚洲站长工具| 国产片内射在线| av在线播放免费不卡| 午夜激情福利司机影院| 少妇 在线观看| 中文字幕人妻丝袜一区二区| 国产97色在线日韩免费| 成人av一区二区三区在线看| 麻豆av在线久日| 99国产精品一区二区三区| 老熟妇乱子伦视频在线观看| 在线视频色国产色| 91国产中文字幕| 又黄又爽又免费观看的视频| 久久中文看片网| 午夜免费激情av| 精品久久久久久久久久免费视频| 日韩免费av在线播放| 欧美成人性av电影在线观看| 久久久久久久久免费视频了| 最好的美女福利视频网| 一本久久中文字幕| 看黄色毛片网站| 国产精品免费一区二区三区在线| 亚洲成国产人片在线观看| 午夜福利高清视频| 午夜老司机福利片| 精品电影一区二区在线| 青草久久国产| 久久国产精品影院| 日日干狠狠操夜夜爽| 真人做人爱边吃奶动态| 免费av毛片视频| 国产成人欧美在线观看| 黄色 视频免费看| 91成年电影在线观看| 日本撒尿小便嘘嘘汇集6| 女性被躁到高潮视频| 国产私拍福利视频在线观看| 亚洲人成伊人成综合网2020| 久久香蕉国产精品| 久久中文字幕一级| 在线播放国产精品三级| 国产精品99久久99久久久不卡| 亚洲精品美女久久久久99蜜臀| 久久久精品国产亚洲av高清涩受| 欧美性长视频在线观看| 哪里可以看免费的av片| 99热这里只有精品一区 | 51午夜福利影视在线观看| 色综合婷婷激情| 亚洲精品美女久久久久99蜜臀| 在线播放国产精品三级| 国产野战对白在线观看| 国产一区二区三区在线臀色熟女| 在线观看免费午夜福利视频| 亚洲 欧美一区二区三区| 精品国产一区二区三区四区第35| 麻豆国产av国片精品| 欧美乱色亚洲激情| 国产精品二区激情视频| 天堂影院成人在线观看| 女生性感内裤真人,穿戴方法视频| 日日夜夜操网爽| 久久精品影院6| av欧美777| 日韩高清综合在线| 免费一级毛片在线播放高清视频| 老熟妇仑乱视频hdxx| 亚洲av日韩精品久久久久久密| e午夜精品久久久久久久| 亚洲国产高清在线一区二区三 | 69av精品久久久久久| 91九色精品人成在线观看| 日韩欧美免费精品| 午夜日韩欧美国产| av有码第一页| 正在播放国产对白刺激| 一区二区三区高清视频在线| a级毛片a级免费在线| 精品欧美国产一区二区三| 日日摸夜夜添夜夜添小说| 亚洲激情在线av| 久久久久国产一级毛片高清牌| 亚洲 欧美 日韩 在线 免费| 国产成人精品无人区| 91在线观看av| 母亲3免费完整高清在线观看| 精品少妇一区二区三区视频日本电影| 狠狠狠狠99中文字幕| 欧美国产精品va在线观看不卡| 淫秽高清视频在线观看| 观看免费一级毛片| 免费在线观看影片大全网站| 亚洲av日韩精品久久久久久密| 99久久国产精品久久久| 国产亚洲精品久久久久久毛片| 国产伦人伦偷精品视频| 亚洲av成人一区二区三| 国产一卡二卡三卡精品| 欧美日韩亚洲国产一区二区在线观看| 久久人人精品亚洲av| 两人在一起打扑克的视频| 一进一出好大好爽视频| 精品久久久久久久人妻蜜臀av| 欧美在线一区亚洲| 亚洲三区欧美一区| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美精品综合一区二区三区| videosex国产| xxx96com| 久久久久久九九精品二区国产 | 国产成年人精品一区二区| 国产成人欧美在线观看| 女人高潮潮喷娇喘18禁视频| 中文字幕精品亚洲无线码一区 | 露出奶头的视频| 国产精品电影一区二区三区| 国产精品二区激情视频| 精品久久蜜臀av无| 亚洲五月色婷婷综合| 久久天堂一区二区三区四区| 在线观看免费午夜福利视频| 午夜两性在线视频| 夜夜爽天天搞| 搞女人的毛片| 国产黄色小视频在线观看| 免费在线观看视频国产中文字幕亚洲| 少妇裸体淫交视频免费看高清 | 国产麻豆成人av免费视频| 99国产精品一区二区三区| av有码第一页| 日韩欧美国产在线观看| 亚洲成人精品中文字幕电影| 国产成人欧美在线观看| 亚洲欧美日韩无卡精品| 日本撒尿小便嘘嘘汇集6| 日韩精品青青久久久久久| x7x7x7水蜜桃| 人人妻人人澡欧美一区二区| 黑人欧美特级aaaaaa片| 久久久精品欧美日韩精品| 亚洲一区二区三区不卡视频| 亚洲av第一区精品v没综合| 在线十欧美十亚洲十日本专区| 俄罗斯特黄特色一大片| 免费一级毛片在线播放高清视频| 国产精品久久久久久精品电影 | 亚洲一区二区三区色噜噜| 亚洲精华国产精华精| 少妇熟女aⅴ在线视频| 欧美在线一区亚洲| 757午夜福利合集在线观看| а√天堂www在线а√下载| 日韩精品中文字幕看吧| 色精品久久人妻99蜜桃| 波多野结衣高清无吗| 中文字幕最新亚洲高清| 久久中文字幕人妻熟女| 身体一侧抽搐| 亚洲avbb在线观看| 亚洲熟妇中文字幕五十中出| 波多野结衣巨乳人妻| 成人亚洲精品一区在线观看| 亚洲午夜理论影院| 黄色视频,在线免费观看| 波多野结衣巨乳人妻| 久久久水蜜桃国产精品网| 成人18禁高潮啪啪吃奶动态图| 一区二区日韩欧美中文字幕| 国产成人啪精品午夜网站| 国产又色又爽无遮挡免费看| 熟女少妇亚洲综合色aaa.| 他把我摸到了高潮在线观看| 亚洲欧洲精品一区二区精品久久久| 岛国视频午夜一区免费看| 夜夜躁狠狠躁天天躁| 深夜精品福利| 中文字幕av电影在线播放| 成在线人永久免费视频| 岛国在线观看网站| 午夜福利在线观看吧| 一区二区三区激情视频| 91国产中文字幕| 悠悠久久av| 99精品久久久久人妻精品| 亚洲精品在线美女| 一级作爱视频免费观看| 一区二区三区高清视频在线| 真人一进一出gif抽搐免费| 麻豆久久精品国产亚洲av| 欧美色欧美亚洲另类二区| 村上凉子中文字幕在线| 老司机福利观看| 国产精品久久久久久精品电影 | 国产主播在线观看一区二区| 国产成人av激情在线播放| 午夜福利免费观看在线| 成人国产综合亚洲| 嫁个100分男人电影在线观看| 老汉色∧v一级毛片| 国产精品美女特级片免费视频播放器 | 久久欧美精品欧美久久欧美| aaaaa片日本免费| av在线播放免费不卡|