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

    考慮J2攝動(dòng)的彈道導(dǎo)彈高精度彈道預(yù)報(bào)和誤差傳播分析

    2016-07-16 03:26:21孟凡坤鄧啟正
    彈道學(xué)報(bào) 2016年2期

    孫 瑜,吳 楠,孟凡坤,鄧啟正

    (1.中國人民解放軍96669部隊(duì),北京 102208;2.信息工程大學(xué) 導(dǎo)航與空天目標(biāo)工程學(xué)院,鄭州 450001)

    ?

    考慮J2攝動(dòng)的彈道導(dǎo)彈高精度彈道預(yù)報(bào)和誤差傳播分析

    孫瑜1,2,吳楠2,孟凡坤2,鄧啟正1

    (1.中國人民解放軍96669部隊(duì),北京 102208;2.信息工程大學(xué) 導(dǎo)航與空天目標(biāo)工程學(xué)院,鄭州 450001)

    摘要:為進(jìn)行高精度的彈道導(dǎo)彈彈道預(yù)報(bào)和誤差傳播分析,分析了平根數(shù)法解攝動(dòng)方程的解形式,提出了一種考慮J2攝動(dòng)并包含短周期項(xiàng)的解析法,并將其用于彈道預(yù)報(bào),結(jié)合無跡變換構(gòu)造了彈道預(yù)報(bào)的誤差傳播分析方法,解決了解析法中開普勒根數(shù)求解過程不可微、偏導(dǎo)數(shù)矩陣難以計(jì)算的問題。仿真結(jié)果表明,將解析法用于彈道預(yù)報(bào)時(shí),J2攝動(dòng)短周期項(xiàng)導(dǎo)致的誤差達(dá)3個(gè)數(shù)量級,不可忽略;該算法預(yù)報(bào)的相對誤差只有不到5%,精度與高精度數(shù)值法相當(dāng),遠(yuǎn)高于傳統(tǒng)解析法,但運(yùn)算速度是后者的7倍,具有更好的效費(fèi)比。

    關(guān)鍵詞:彈道預(yù)報(bào);無跡變換;J2攝動(dòng);解析法;短周期項(xiàng)

    彈道預(yù)報(bào)與誤差傳播分析是戰(zhàn)略預(yù)警系統(tǒng)的重要研究內(nèi)容。傳統(tǒng)的彈道預(yù)報(bào)方法通常有數(shù)值法和解析法兩類。數(shù)值法是在充分考慮導(dǎo)彈飛行過程中各種攝動(dòng)因素的基礎(chǔ)上,建立導(dǎo)彈動(dòng)力學(xué)方程,并采用一定的積分準(zhǔn)則,利用數(shù)值逼近的方法計(jì)算彈道。文獻(xiàn)[1~2]基于導(dǎo)彈動(dòng)力學(xué)方程采用不同的數(shù)值積分方法計(jì)算彈道,驗(yàn)證了數(shù)值法的有效性。然而,這種方法預(yù)報(bào)的精確程度及計(jì)算速度與積分步長的選取有很大關(guān)系。小步長積分可以保證較高的預(yù)報(bào)精度,但是計(jì)算速度較慢;大步長積分計(jì)算速度快,但是預(yù)報(bào)精度不高。解析法是將導(dǎo)彈自由段飛行軌跡近似為橢圓彈道,通過計(jì)算開普勒根數(shù)進(jìn)行彈道預(yù)報(bào)[3]。文獻(xiàn)[4~7]基于橢圓彈道理論,采用解析法實(shí)現(xiàn)彈道預(yù)報(bào),并分別從時(shí)間和預(yù)報(bào)精度等方面提出優(yōu)化算法。盡管采用了優(yōu)化算法,但解析法的精度仍低于數(shù)值法,這是因?yàn)榻馕龇▽⒌厍蛞曌髻|(zhì)點(diǎn),沒有考慮地球是形狀不規(guī)則的扁球體。文獻(xiàn)[8~10]用考慮J2攝動(dòng)的解析法計(jì)算人造衛(wèi)星的軌道,因?yàn)槿嗽煨l(wèi)星的運(yùn)行周期很長,所以在解攝動(dòng)方程的時(shí)候?yàn)榱擞?jì)算簡便,通常忽略短周期項(xiàng)。然而導(dǎo)彈與人造衛(wèi)星截然不同的是導(dǎo)彈的飛行時(shí)間很短,它的短周期項(xiàng)能否忽略值得討論。

    本文提出了一種考慮J2攝動(dòng)并包含短周期項(xiàng)的解析法,并將其用于彈道預(yù)報(bào),與數(shù)值法彈道預(yù)報(bào)進(jìn)行仿真對比預(yù)報(bào)的精度及耗時(shí)。同時(shí),在誤差傳播分析過程中,結(jié)合無跡變換(Unscented Transformation,UT)構(gòu)造了彈道預(yù)報(bào)的誤差傳播分析方法,解決了解析法中開普勒根數(shù)求解過程不可微、偏導(dǎo)數(shù)矩陣難以計(jì)算的問題。

    1考慮J2攝動(dòng)并包含短周期項(xiàng)的彈道預(yù)報(bào)解析法

    彈道導(dǎo)彈關(guān)機(jī)后進(jìn)入自由段依靠慣性飛行,而自由段的高度一般在200 km以上,此時(shí)空氣阻力的影響可以忽略不計(jì),導(dǎo)彈主要外力源為地球中心引力,可以將導(dǎo)彈繞地球的運(yùn)動(dòng)視為二體運(yùn)動(dòng),從理論上可以認(rèn)為其彈道為橢圓彈道,落點(diǎn)為橢圓彈道與地球的交點(diǎn)。解析法通過計(jì)算橢圓彈道的6個(gè)開普勒根數(shù)來預(yù)報(bào)導(dǎo)彈任意時(shí)刻的狀態(tài)。然而地球作為一個(gè)形狀不規(guī)則的扁球體,其重力場對自由段的影響很大,二階帶諧系數(shù)J2對彈道導(dǎo)彈自由段落點(diǎn)偏差的影響最大可達(dá)十幾km,四階帶諧系數(shù)J4對自由段落點(diǎn)偏差的影響可達(dá)數(shù)十m[11-12]。文獻(xiàn)[13]采用數(shù)值積分法,分析了不同積分步長、不同緯度和發(fā)射方位角下J2和J4項(xiàng)以及不同階數(shù)的擾動(dòng)重力對導(dǎo)彈自由段落點(diǎn)的影響。J2項(xiàng)導(dǎo)致的攝動(dòng)加速度約為10-3量級,而J4項(xiàng)攝動(dòng)與大氣阻力、第三體引力、太陽光壓等其他攝動(dòng)均為10-6量級及以下。若考慮J4項(xiàng)攝動(dòng),則大氣阻力、三體引力、太陽光壓等攝動(dòng)均不能忽略,此時(shí)將大幅增加計(jì)算耗時(shí)。因此,綜合考慮預(yù)報(bào)的精度和耗時(shí),本文只考慮J2項(xiàng)攝動(dòng)的影響。

    1.1J2攝動(dòng)運(yùn)動(dòng)方程

    地球是一個(gè)形狀不規(guī)則的扁球體,要精確地積分計(jì)算出其勢函數(shù),必須知道地球的表面形狀及內(nèi)部的密度分布。目前,通常采用球函數(shù)展開式推導(dǎo)出地球勢函數(shù)的標(biāo)準(zhǔn)表達(dá)式[14]:

    (Cnmcosmλ+Snmsinmλ)Pnm(sinφ)

    (1)

    式中:f為萬有引力常數(shù);me為地球質(zhì)量;r為地球外地心距;ae為地球赤道平均半徑;Jn為帶諧系數(shù),Jn=-Cn0;n≠m時(shí),Cnm,Snm稱為田諧系數(shù);n=m時(shí),Cnm,Snm稱為扇諧系數(shù);Pn(sinφ)稱為勒讓德函數(shù);Pnm(sinφ)稱為締合勒讓德函數(shù);φ,λ為地球的地心緯度和經(jīng)度。

    物體繞地球運(yùn)動(dòng)的微分方程可寫成:

    (2)

    (Cnmcosmλ+Snmsinmλ)Pnm(sinφ)

    (3)

    當(dāng)n=2時(shí),即為考慮J2項(xiàng)攝動(dòng),取J2=1.082 63×10-3。

    以6個(gè)開普勒根數(shù)為基本變量的拉格朗日型攝動(dòng)運(yùn)動(dòng)方程[15]為

    (4)

    式中:a為軌道半長軸;e為偏心率;i為軌道傾角;Ω為升交點(diǎn)赤經(jīng);A為近地點(diǎn)幅角;M為平近點(diǎn)角;ω為平均角速度。設(shè)σ=(aeiΩAM)為開普勒根數(shù)構(gòu)成的向量。

    角析法中,平均根數(shù)法解攝動(dòng)運(yùn)動(dòng)方程的特征是將攝動(dòng)變化中不同性質(zhì)的項(xiàng)區(qū)分開,相應(yīng)的冪級數(shù)解的形式為

    (5)

    式中:

    1.2開普勒根數(shù)轉(zhuǎn)移

    對于衛(wèi)星而言,由于其在軌運(yùn)行時(shí)間長,在計(jì)算其攝動(dòng)運(yùn)動(dòng)方程時(shí),為了計(jì)算簡便通常忽略短周期項(xiàng)。對于彈道導(dǎo)彈,其飛行時(shí)間很短,所以在計(jì)算彈道導(dǎo)彈的攝動(dòng)運(yùn)動(dòng)方程的時(shí)候不能忽略短周期項(xiàng)。而對于具體取到多少階的攝動(dòng)解,則要取決于預(yù)報(bào)的精度需求。文獻(xiàn)[13]分析了不同階數(shù)對彈道導(dǎo)彈自由段落點(diǎn)的影響,中低階部分的影響從幾十m到數(shù)百m不等,對于自由段射程約6 000km的彈道導(dǎo)彈,考慮到30階左右可以達(dá)到m級的計(jì)算精度。本文重點(diǎn)在于綜合考慮預(yù)報(bào)的精度和預(yù)警時(shí)間,因此本文只考慮一階攝動(dòng)解:包括所有的一階攝動(dòng)項(xiàng)(長期項(xiàng)、長周期項(xiàng)和短周期項(xiàng))和二階長期項(xiàng)。

    考慮J2項(xiàng)攝動(dòng)的平均根數(shù)法一階解完整結(jié)果如下:

    (6)

    2基于無跡變換的誤差傳播分析

    計(jì)算協(xié)方差矩陣時(shí),純粹的解析法需要計(jì)算笛卡爾坐標(biāo)對開普勒根數(shù)的偏導(dǎo)數(shù)矩陣以及開普勒根數(shù)的轉(zhuǎn)移矩陣,然而在計(jì)算開普勒根數(shù)轉(zhuǎn)移矩陣時(shí)會(huì)出現(xiàn)一個(gè)問題:如果考慮短周期項(xiàng),則開普勒根數(shù)求解過程不可微,偏導(dǎo)數(shù)矩陣難以計(jì)算;如果不考慮短周期項(xiàng),計(jì)算轉(zhuǎn)移矩陣所用的模型就與狀態(tài)方程中所用的模型不一致,得到的解是完全錯(cuò)誤的,這個(gè)結(jié)論在文獻(xiàn)[16]中得到證實(shí)。而UT用采樣點(diǎn)來近似狀態(tài)向量的后驗(yàn)概率密度,結(jié)合UT的解析法在對彈道導(dǎo)彈進(jìn)行誤差傳播分析時(shí),就不需要計(jì)算笛卡爾坐標(biāo)對開普勒根數(shù)的偏導(dǎo)數(shù)矩陣以及開普勒根數(shù)的轉(zhuǎn)移矩陣。

    (7)

    (8)

    通過橢圓彈道理論計(jì)算每個(gè)χ點(diǎn)對應(yīng)的開普勒根數(shù)σ,然后考慮J2攝動(dòng)采用平根數(shù)法解拉格朗日攝動(dòng)運(yùn)動(dòng)方程,計(jì)算出t時(shí)刻對應(yīng)的開普勒根數(shù)σt。通過這些非線性變換得到樣點(diǎn)的輸出:

    Ψj=f(χj)j=0,1,…,12

    (9)

    最后,計(jì)算t時(shí)刻位置速度的均值和協(xié)方差矩陣:

    (10)

    基于UT變換來處理協(xié)方差矩陣的非線性傳遞,可以避開計(jì)算笛卡爾坐標(biāo)對開普勒根數(shù)的偏導(dǎo)數(shù)矩陣以及開普勒根數(shù)的轉(zhuǎn)移矩陣,并且對任意t時(shí)刻的協(xié)方差矩陣可以直接計(jì)算,不需要計(jì)算t時(shí)刻之前的協(xié)方差矩陣,能夠有效提高彈道預(yù)報(bào)效率。

    圖1 系統(tǒng)流程框圖

    3仿真分析

    仿真將通過普通解析法、基于UT的解析法和數(shù)值法對彈道導(dǎo)彈自由段進(jìn)行預(yù)報(bào),忽略空氣阻力等其他攝動(dòng)力,只考慮J2攝動(dòng)的影響。普通解析法是直接解算攝動(dòng)運(yùn)動(dòng)方程,通過計(jì)算開普勒根數(shù)變化率,解算出預(yù)報(bào)終點(diǎn)的位置速度;基于UT的解析法即上文所提的方法;數(shù)值法是通過對目標(biāo)的運(yùn)動(dòng)方程進(jìn)行積分計(jì)算彈道的方法,數(shù)值積分采用四階龍格-庫塔算法,仿真步長設(shè)為1s。首先分析解析法中短周期項(xiàng)對彈道預(yù)報(bào)的影響,然后分析基于UT的解析法計(jì)算的協(xié)方差矩陣和普通解析法進(jìn)行蒙特卡洛仿真得到的均方根誤差的一致性,最后對比分析解析法和數(shù)值法的預(yù)報(bào)精度和運(yùn)算速度。設(shè)在地心慣性坐標(biāo)系下,探測到導(dǎo)彈在自由段500s時(shí)刻的位置為x=2 163 451m,y=-6 342 006m,z=2 191 877m,速度vx=594.31m/s,vy=711.49m/s,vz=5 453.41m/s,協(xié)方差為

    預(yù)報(bào)終點(diǎn)t=1 287s。

    3.1短周期項(xiàng)對彈道預(yù)報(bào)的影響

    有短周期項(xiàng)的解析法與數(shù)值法的均值差變化情況及無短周期項(xiàng)的解析法與數(shù)值法的均值變化情況如圖2所示。

    圖2 均值之差的變化情況

    從圖2中可以看到,無短周期項(xiàng)的解析法與數(shù)值法計(jì)算的均值差在位置的x軸、y軸方向及速度的z軸方向均成發(fā)散趨勢,速度的x軸、y軸方向及位置的z軸方向有成周期變化的趨勢。位置的x軸方向終點(diǎn)處偏差約12km,y軸方向終點(diǎn)處偏差達(dá)到25km,z軸方向偏差最大達(dá)到15km,終點(diǎn)處約14km;速度的x軸方向偏差最大有3m/s,終點(diǎn)處為0.12m/s,y軸方向偏差最大達(dá)到14m/s,終點(diǎn)處為12.38m/s,z軸方向終點(diǎn)處偏差為38.08m/s。而與無短周期項(xiàng)相比,考慮了短周期項(xiàng)的解析法與數(shù)值法之差則小得多,接近于0。t=1 287s時(shí)刻,均值差的具體結(jié)果如表1所示。由此可見,由于彈道導(dǎo)彈飛行時(shí)間短,短周期項(xiàng)對彈道預(yù)報(bào)的影響很大,位置偏差能夠達(dá)到3個(gè)數(shù)量級,所以在進(jìn)行高精度彈道預(yù)報(bào)時(shí)不能忽略短周期項(xiàng)。

    表1 均值之差仿真結(jié)果

    3.2誤差一致性分析

    采用UT變換解析法計(jì)算的標(biāo)準(zhǔn)差和普通解析法1 000次蒙特卡洛仿真得到的均方根誤差的變化情況如圖3所示。

    圖3 普通解析法均方根誤差與基于UT變換解析法的標(biāo)準(zhǔn)差變化情況

    從圖3中可以看到,二者的趨勢基本相同,且預(yù)報(bào)的協(xié)方差位置和速度均略小于均方根誤差。在預(yù)報(bào)的終點(diǎn)處,位置x軸的相對誤差為1%,位置y軸的相對誤差為5.7%,位置z軸的相對誤差為3.03%,速度x軸的相對誤差為1.7%,速度y軸的相對誤差為7.1%,速度z軸的相對誤差為2.6%??梢?二者的相對誤差非常小,而UT變換計(jì)算的協(xié)方差矩陣比蒙特卡洛仿真計(jì)算的均方根誤差略小,是因?yàn)閰f(xié)方差矩陣是截?cái)嗾`差,忽略了高階項(xiàng),屬于正常誤差,可以說二者是一致的。

    3.3彈道預(yù)報(bào)和誤差傳播分析

    圖4給出了考慮J2攝動(dòng)解析法的標(biāo)準(zhǔn)差和數(shù)值法1 000次蒙特卡洛仿真得到的均方根誤差的變化情況。從圖中可以看到,二者的變化趨勢基本一致,解析法的協(xié)方差矩陣在x軸、y軸方向的位置和速度均略大于數(shù)值法的均方根誤差,在z軸方向的位置和速度均略小于數(shù)值法的均方根誤差,說明解析法的精度在x軸和y軸略低于數(shù)值法,而z軸略高于數(shù)值法。在預(yù)報(bào)的終點(diǎn)處,位置x軸的相對誤差為4.6%,位置y軸的相對誤差為4.2%,位置z軸的相對誤差為2.5%,速度x軸的相對誤差為4.8%,速度y軸的相對誤差為3.5%,速度z軸的相對誤差為1.1%。無論位置還是速度,兩種方法的相對誤差都不超過5%,說明解析法和數(shù)值法的預(yù)報(bào)精度相當(dāng),基本在同一數(shù)量級。

    考慮J2攝動(dòng)的解析法的總耗時(shí)為128.075ms,其中每個(gè)χ點(diǎn)的總耗時(shí)為9.851ms。每個(gè)χ點(diǎn)在傳播的過程中,500s初始時(shí)刻的位置速度轉(zhuǎn)換成開普勒根數(shù)的耗時(shí)最大,為4.478ms,占其傳播總耗時(shí)的45%左右;計(jì)算500s時(shí)刻平均根數(shù)的耗時(shí)為2.001ms,占20%左右;計(jì)算長期項(xiàng)耗時(shí)為0.883ms,占9%左右;計(jì)算1 287s時(shí)刻的平均根數(shù)的耗時(shí)為0.616ms,占6%左右;計(jì)算短周期項(xiàng)和長周期項(xiàng)的耗時(shí)分別為0.755ms和0.698ms,分別占8%和7%左右??梢?在解析法中位置速度與開普勒根數(shù)之間的轉(zhuǎn)換耗時(shí)最大。數(shù)值法積分步長為1s,每步長的平均耗時(shí)約為1.048ms,預(yù)報(bào)自由段時(shí)長為788s,數(shù)值法總耗時(shí)為826.602ms,是解析法的7倍左右。

    預(yù)報(bào)飛行時(shí)長800s左右的彈道得到相同級別的精度,數(shù)值法計(jì)算耗時(shí)可以達(dá)到考慮J2攝動(dòng)解析法的7倍左右。原因是數(shù)值法在計(jì)算任意時(shí)刻導(dǎo)彈狀態(tài)的時(shí)候,需要按照一定的步長逐步積分,耗時(shí)與所選取的步長成反比,保證一定的精度就需要選取較小的步長,而以耗時(shí)的增加為代價(jià);而解析法中,在觀測到導(dǎo)彈初始狀態(tài)后,就能夠快速計(jì)算出彈道6個(gè)開普勒根數(shù)的變化情況,只要給出任意時(shí)刻點(diǎn),就可以直接計(jì)算導(dǎo)彈在該時(shí)刻的狀態(tài),所以計(jì)算速度快。

    圖4 數(shù)值法均方根誤差與考慮J2攝動(dòng)解析法標(biāo)準(zhǔn)差的變化情況

    綜上所述,總結(jié)得到如下結(jié)論:

    ①由于彈道導(dǎo)彈飛行時(shí)間短,短周期項(xiàng)對彈道預(yù)報(bào)影響很大,位置偏差能夠達(dá)到3個(gè)數(shù)量級,在進(jìn)行高精度彈道預(yù)報(bào)的時(shí)候不能忽略;

    ②結(jié)合UT的解析法進(jìn)行誤差傳播分析,能夠有效避開計(jì)算笛卡爾坐標(biāo)對開普勒根數(shù)的偏導(dǎo)數(shù)矩陣以及開普勒根數(shù)的轉(zhuǎn)移矩陣;

    ③考慮J2攝動(dòng)的解析法與積分步長為1s的數(shù)值法相比,彈道預(yù)報(bào)的精度在同一數(shù)量級,并且在預(yù)報(bào)800s左右的彈道時(shí),耗時(shí)只有數(shù)值法的1/7左右。

    4結(jié)束語

    本文通過仿真對解析法和數(shù)值法在彈道預(yù)報(bào)中的性能進(jìn)行了對比分析。解析法在進(jìn)行誤差傳播分析時(shí)結(jié)合UT,能夠有效避開計(jì)算笛卡爾坐標(biāo)對開普勒根數(shù)的偏導(dǎo)數(shù)矩陣以及開普勒根數(shù)的轉(zhuǎn)移矩陣。傳統(tǒng)的解析法只是簡單地近似成二體問題,沒有考慮其他攝動(dòng)因素,所以計(jì)算精度不高;而考慮了J2攝動(dòng)之后的解析法,計(jì)算精度和誤差傳播發(fā)散程度都與積分步長為1s的數(shù)值法相當(dāng),在預(yù)報(bào)800s左右的彈道時(shí)運(yùn)算速度是數(shù)值法的7倍左右。但是考慮J2攝動(dòng)的解析法在彈道預(yù)報(bào)中,不能夠像計(jì)算衛(wèi)星軌道一樣簡單地忽略短周期項(xiàng),短周期項(xiàng)對彈道預(yù)報(bào)影響很大,在進(jìn)行高精度的預(yù)報(bào)彈道時(shí)不能忽略。

    參考文獻(xiàn)

    [1]張龍杰,謝曉方.通用變步長彈道仿真方法研究[J].彈道學(xué)報(bào),2011,23(2):42-46.

    ZHANG Long-jie,XIE Xiao-fang.Research on general variable step-size ballistic simulation method[J].Journal of Ballistics,2011,23(2):42-46.(in Chinese)

    [2]程誠,張小兵.某制導(dǎo)炮彈二位兩相流內(nèi)彈道性能分析與數(shù)值模擬研究[J].兵工學(xué)報(bào),2015,36(1):58-63.

    CHENG Cheng,ZHANG Xiao-bing.Two-dimensional numerical simulation on two-phase flow interior ballistic performance of a guided projectile[J].Acta Armamentarii,2015,36(1):58-63.(in Chinese)

    [3]張毅,楊慧耀,李俊莉.彈道導(dǎo)彈彈道學(xué)[M].長沙:國防科技大學(xué)出版社,1999.

    ZHANG Yi,YANG Hui-yao,LI Jun-li.Ballistic missile ballistics[M].Changsha:National University of Defense Technology Press,1999.(in Chinese)

    [4]李森,劉紅軍,程仲.提高彈道導(dǎo)彈落點(diǎn)預(yù)報(bào)精度方法[J].電子科技,2012,25(9):85-87.

    LI Sen,LIU Hong-jun,CHENG Zhong.Improvement of TBM’s impact-point prediction precision[J].Electronic Sci & Tech,2012,25(9):85-87.(in Chinese)

    [5]沈慧娜,徐振來.提高彈道導(dǎo)彈落點(diǎn)預(yù)報(bào)精度的動(dòng)弧平滑平均法[J].現(xiàn)代雷達(dá),2009,31(7):55-58.

    SHEN Hui-na,XU Zhen-lai.Moving arc smoothing averaging algorithm for improving the prediction precision of ballistic missile impact point[J].Modern Radar,2009,31(7):55-58.(in Chinese)

    [6]劉仁,王愛華,郭桂治.基于關(guān)機(jī)點(diǎn)狀態(tài)的戰(zhàn)術(shù)彈道導(dǎo)彈落點(diǎn)估計(jì)[J].空軍工程大學(xué)學(xué)報(bào),2010,11(1):27-30.

    LIU Ren,WANG Ai-hua,GUO Gui-zhi.Impact point estimation of tactical ballistic missile based on the state of burnout point[J].Journal of Air Force Engineering University,2010,11(1):27-30.(in Chinese)

    [7]劉彥君,喬士東,黃金才,等.一種高精度彈道導(dǎo)彈落點(diǎn)預(yù)測方法[J].彈道學(xué)報(bào),2012,24(1):23-26.

    LIU Yan-jun,QIAO Shi-dong,HUANG Jin-cai,et al.A method of impact point prediction of ballistic missile[J].Journal of Ballistics,2012,24(1):23-26.(in Chinese)

    [8]劉林.一種人造地球衛(wèi)星的攝動(dòng)計(jì)算方法[J].天文學(xué)報(bào),1975,16(1):5-80.

    LIU Lin.A calculating method of perturbation of artificial satellite[J].Acta Astronomica Sinica,1975,16(1):5-80.(in Chinese)

    [9]曹文華,李傳榮,李子揚(yáng),等.基于攝動(dòng)模型的衛(wèi)星軌道預(yù)測算法精度評估[J].遙感信息,2012,27(6):21-27.

    CAO Wen-hua,LI Chuan-rong,LI Zi-yang,et al.Accuracy assessment of satellite orbit prediction based onJ2perturbation model[J].Remote Sensing Information,2012,27(6):21-27.(in Chinese)

    [10]張俊.基于開普勒二體運(yùn)動(dòng)修正地球扁率J2攝動(dòng)項(xiàng)算法[J].航天控制,2014,32(6):22-25.

    ZHANG Jun.Correction algorithm ofJ2perturbations of the earth oblateness based on Kepler two-body motion[J].Aerospace Control,2014,32(6):22-25.(in Chinese)

    [11]鄭偉,湯國建.上升段引力常數(shù)變化及其對彈道導(dǎo)彈運(yùn)動(dòng)的影響[J].導(dǎo)彈與航天運(yùn)載技術(shù),2008,295(3):12-14.

    ZHENG Wei,TANG Guo-jian.Gravitation constant variety in ascend phase and its effect on ballistic missile[J].Missile and Space Vehcile,2008,295(3):12-14.(in Chinese)

    [12]鄭偉,湯國建.彈道導(dǎo)彈自由段解算的等高約束解析解[J].宇航學(xué)報(bào),2007,28(2):269-272.

    ZHENG Wei,TANG Guo-jian.Contour restricted analytical solution for free flight trajectory of ballistic missile[J].Journal of Astronautics,2007,28(2):269-272.(in Chinese)

    [13]馮偉,劉根友,郝曉光.重力場對彈道導(dǎo)彈自由段落點(diǎn)影響的仿真分析[J].測繪科學(xué),2009,34(5):118-120.

    FENG Wei,LIU Gen-you,HAO Xiao-guang.Simulated analysis of the effect of gravity field to the fall point in free flight trajectory of ballistic missile[J],Science of Surveying and Mapping,2009,34(5):118-120.(in Chinese)

    [14]賈沛然,陳克俊,何力.遠(yuǎn)程火箭彈道學(xué)[M].長沙:國防科技大學(xué)出版社,1993.

    JIA Pei-ran,CHEN Ke-jun,HE Li.Long-range rocket ballistics[M].Changsha:National University of Defense Technology Press,1993.(in Chinese)

    [15]李濟(jì)生,葉杰,王家松,等.航天器軌道確定[M].北京:國防工業(yè)出版社,2003.

    LI Ji-sheng,YE Jie,WANG Jia-song,et al.Determination of spacecraft[M].Beijing:National Defense Industry Press,2003.(in Chinese)

    [16]MAY J A.A study of the effects of state transition matrix approximations,NASA CP 21123[R].1980.

    Analysis on High Precision Trajectory Prediction and Error Propagation of Ballistic Missile ConsideringJ2Perturbation

    SUN Yu1,2,WU Nan2,MENG Fan-kun2,DENG Qi-zheng1

    (1.Unit 96669 of PLA,Beijing 102208,China;2.College of Navigation and Aerospace Target,Information Engineering University,Zhengzhou 450001,China)

    Abstract:To carry out high-precision trajectory prediction and error propagation analysis of ballistic missile,the solution form of mean element algorithms for perturbation equation was analyzed,and an analytical algorithm consideringJ2perturbation and containing the short period term was proposed for trajectory prediction.By constructing the error propagation analysis method of trajectory prediction based on unscented transform,the problem that the Kepler elements is not differentiable and the partial derivative matrix is difficult to calculate was solved.The simulation results show that the error caused by theJ2perturbation short period term,which can not be ignored,reaches 3 orders of magnitude when analytical algorithm is applied in the trajectory prediction.The relative error predicted by the proposed algorithm is only less than 5%,and the precision is far higher than the precision of traditional analytical algorithm,which is comparable with the precision of high-precision numerical algorithm,but the operation speed is 7 times of the latter,and the algorithm has a better effectiveness-cost ratio.

    Key words:trajectory prediction;unscented transformation;J2perturbation;analytical algorithm;short-period term

    收稿日期:2015-12-30

    作者簡介:孫瑜(1986- ),男,助理工程師,碩士,研究方向?yàn)榭仗炷繕?biāo)預(yù)警信息處理。E-mail:woshizhongsy@163.com。

    中圖分類號:TJ760.1

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1004-499X(2016)02-0018-07

    精品人妻1区二区| www日本黄色视频网| 亚洲国产欧洲综合997久久,| 黄片大片在线免费观看| 制服人妻中文乱码| 熟女电影av网| 午夜福利欧美成人| 欧美日韩乱码在线| 国产日本99.免费观看| 日日夜夜操网爽| 国产一级毛片七仙女欲春2| 国产av一区二区精品久久| 欧美日韩乱码在线| 一卡2卡三卡四卡精品乱码亚洲| 搡老熟女国产l中国老女人| 国产爱豆传媒在线观看 | 午夜老司机福利片| 91在线观看av| 又粗又爽又猛毛片免费看| 哪里可以看免费的av片| bbb黄色大片| 国产探花在线观看一区二区| 国产爱豆传媒在线观看 | 又黄又粗又硬又大视频| 五月伊人婷婷丁香| 香蕉丝袜av| 久9热在线精品视频| 国产在线精品亚洲第一网站| 老司机靠b影院| 中文字幕人妻丝袜一区二区| 国产成人aa在线观看| 一级毛片女人18水好多| 国产亚洲精品久久久久久毛片| 国产主播在线观看一区二区| 日韩大码丰满熟妇| 国产aⅴ精品一区二区三区波| 黄色成人免费大全| 美女高潮喷水抽搐中文字幕| 欧美日韩一级在线毛片| 麻豆一二三区av精品| 两个人的视频大全免费| videosex国产| 最新在线观看一区二区三区| 无限看片的www在线观看| 人人妻人人澡欧美一区二区| 亚洲成人国产一区在线观看| 免费人成视频x8x8入口观看| 可以免费在线观看a视频的电影网站| 可以免费在线观看a视频的电影网站| 国产精品一区二区免费欧美| 国内精品久久久久久久电影| 欧美+亚洲+日韩+国产| 国产av在哪里看| 色综合亚洲欧美另类图片| 日本 av在线| av在线天堂中文字幕| 久久性视频一级片| 国产精品 欧美亚洲| 亚洲自偷自拍图片 自拍| 成人特级黄色片久久久久久久| 999精品在线视频| 日韩欧美国产一区二区入口| av国产免费在线观看| 国产精品一区二区三区四区久久| 久久人妻福利社区极品人妻图片| 欧美成人一区二区免费高清观看 | 美女午夜性视频免费| 亚洲男人天堂网一区| 精品国产乱码久久久久久男人| 一级片免费观看大全| 午夜激情福利司机影院| 天堂影院成人在线观看| 欧美性长视频在线观看| 一进一出好大好爽视频| 麻豆成人午夜福利视频| 欧美最黄视频在线播放免费| 人人妻人人澡欧美一区二区| 亚洲国产欧美一区二区综合| 国产探花在线观看一区二区| 色在线成人网| 久久精品成人免费网站| 可以在线观看毛片的网站| 999久久久国产精品视频| 露出奶头的视频| 三级毛片av免费| 欧美日本视频| 黄色女人牲交| 欧美日韩一级在线毛片| 亚洲欧美日韩东京热| 美女免费视频网站| 欧美成狂野欧美在线观看| 精品国产乱子伦一区二区三区| 黄色 视频免费看| 国产一区二区在线av高清观看| 国产麻豆成人av免费视频| 在线观看舔阴道视频| 欧美日本亚洲视频在线播放| 久久天躁狠狠躁夜夜2o2o| 国产aⅴ精品一区二区三区波| 亚洲精品中文字幕在线视频| 香蕉av资源在线| 18禁美女被吸乳视频| 国产成人啪精品午夜网站| 麻豆成人av在线观看| 亚洲男人天堂网一区| 久99久视频精品免费| 国产亚洲精品av在线| 无遮挡黄片免费观看| 国产高清有码在线观看视频 | 亚洲真实伦在线观看| 黄频高清免费视频| 久久人妻av系列| 老司机靠b影院| 91成年电影在线观看| 国产视频内射| 在线观看66精品国产| 欧美黄色片欧美黄色片| 可以在线观看的亚洲视频| 欧美日本视频| 91麻豆精品激情在线观看国产| 制服诱惑二区| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久大精品| 欧美中文综合在线视频| 精品电影一区二区在线| 亚洲黑人精品在线| 国产av麻豆久久久久久久| 99国产极品粉嫩在线观看| 免费观看人在逋| svipshipincom国产片| 不卡一级毛片| 搡老妇女老女人老熟妇| 黄色a级毛片大全视频| 91字幕亚洲| 国产精品精品国产色婷婷| 国产成人欧美在线观看| 亚洲精品中文字幕在线视频| 日日干狠狠操夜夜爽| 中文在线观看免费www的网站 | 久久草成人影院| 久热爱精品视频在线9| www.999成人在线观看| 99热6这里只有精品| 亚洲欧美日韩东京热| 久久精品国产清高在天天线| 两人在一起打扑克的视频| 成年版毛片免费区| 欧美性猛交黑人性爽| 桃色一区二区三区在线观看| 色综合婷婷激情| 麻豆成人av在线观看| 精品国产超薄肉色丝袜足j| 亚洲自拍偷在线| 久久久久久久午夜电影| 后天国语完整版免费观看| 女人被狂操c到高潮| 免费在线观看日本一区| 午夜福利欧美成人| 国产成人av教育| 成人永久免费在线观看视频| 波多野结衣高清无吗| 99精品欧美一区二区三区四区| 午夜亚洲福利在线播放| 亚洲狠狠婷婷综合久久图片| 国产亚洲欧美98| 又黄又爽又免费观看的视频| 美女免费视频网站| 亚洲av成人av| 久久性视频一级片| 色尼玛亚洲综合影院| 久久欧美精品欧美久久欧美| 免费人成视频x8x8入口观看| 亚洲av片天天在线观看| 宅男免费午夜| 久久久久国产一级毛片高清牌| 亚洲人成网站在线播放欧美日韩| 我要搜黄色片| 看免费av毛片| 久久久久亚洲av毛片大全| 麻豆国产97在线/欧美 | 中文字幕高清在线视频| 欧美精品啪啪一区二区三区| 97碰自拍视频| 国产精品久久久久久精品电影| 叶爱在线成人免费视频播放| 久久这里只有精品19| 啦啦啦韩国在线观看视频| 12—13女人毛片做爰片一| 久久久久久久久久黄片| 这个男人来自地球电影免费观看| 在线播放国产精品三级| 又黄又爽又免费观看的视频| 巨乳人妻的诱惑在线观看| 叶爱在线成人免费视频播放| 成年免费大片在线观看| 亚洲国产日韩欧美精品在线观看 | 男女那种视频在线观看| 黑人巨大精品欧美一区二区mp4| 91麻豆精品激情在线观看国产| 久久精品成人免费网站| 黑人欧美特级aaaaaa片| 欧美成人午夜精品| 亚洲精品在线观看二区| 伦理电影免费视频| 看免费av毛片| 人人妻,人人澡人人爽秒播| 床上黄色一级片| 他把我摸到了高潮在线观看| 这个男人来自地球电影免费观看| 欧美国产日韩亚洲一区| 日韩国内少妇激情av| 亚洲成人免费电影在线观看| 丰满人妻一区二区三区视频av | 黄色毛片三级朝国网站| 亚洲第一电影网av| 午夜视频精品福利| 亚洲av电影在线进入| 国产日本99.免费观看| 无限看片的www在线观看| 久久精品国产亚洲av香蕉五月| 国产成人精品无人区| 国产精品 欧美亚洲| 老司机靠b影院| 国产三级黄色录像| 欧美成狂野欧美在线观看| 可以在线观看的亚洲视频| 亚洲专区字幕在线| 国产私拍福利视频在线观看| 在线观看日韩欧美| 欧美最黄视频在线播放免费| 精品国产亚洲在线| 精品久久久久久久毛片微露脸| 97超级碰碰碰精品色视频在线观看| 亚洲欧美日韩东京热| 色综合亚洲欧美另类图片| 禁无遮挡网站| 亚洲人成伊人成综合网2020| 成人三级做爰电影| 天天一区二区日本电影三级| 精华霜和精华液先用哪个| 免费搜索国产男女视频| 免费在线观看日本一区| 99久久久亚洲精品蜜臀av| 国产成年人精品一区二区| 18禁裸乳无遮挡免费网站照片| 精品国内亚洲2022精品成人| 色在线成人网| 两个人视频免费观看高清| 成年女人毛片免费观看观看9| 搡老熟女国产l中国老女人| 欧美成人午夜精品| 熟女少妇亚洲综合色aaa.| 中文字幕av在线有码专区| 一进一出抽搐gif免费好疼| 久久精品影院6| 高清毛片免费观看视频网站| 久久人妻av系列| 不卡一级毛片| 黄色视频不卡| 免费在线观看亚洲国产| 女人爽到高潮嗷嗷叫在线视频| 色噜噜av男人的天堂激情| 欧美成人一区二区免费高清观看 | 精品久久久久久成人av| 老司机在亚洲福利影院| 国产免费av片在线观看野外av| 国产成人影院久久av| 精品国内亚洲2022精品成人| 婷婷精品国产亚洲av在线| 性色av乱码一区二区三区2| 午夜福利成人在线免费观看| 亚洲精品一区av在线观看| 最近在线观看免费完整版| av福利片在线观看| 人成视频在线观看免费观看| 免费在线观看视频国产中文字幕亚洲| 婷婷精品国产亚洲av| 听说在线观看完整版免费高清| 午夜精品一区二区三区免费看| 国产高清激情床上av| 国产97色在线日韩免费| 又粗又爽又猛毛片免费看| 成人手机av| 国产精品98久久久久久宅男小说| 叶爱在线成人免费视频播放| 久久精品国产综合久久久| 老汉色av国产亚洲站长工具| 黑人欧美特级aaaaaa片| 99热只有精品国产| 久久这里只有精品19| 亚洲一区二区三区不卡视频| 我的老师免费观看完整版| 久久久久久国产a免费观看| 黄色片一级片一级黄色片| 男男h啪啪无遮挡| 亚洲精品中文字幕在线视频| 久久久久久大精品| 99久久综合精品五月天人人| 一级黄色大片毛片| 亚洲精品美女久久久久99蜜臀| 精品无人区乱码1区二区| 免费在线观看视频国产中文字幕亚洲| 国产免费av片在线观看野外av| 国产熟女午夜一区二区三区| 国产熟女xx| 日韩欧美一区二区三区在线观看| 欧美色视频一区免费| 女人爽到高潮嗷嗷叫在线视频| 亚洲片人在线观看| 久久久国产成人免费| 国产熟女xx| 99久久精品热视频| 少妇的丰满在线观看| 男人舔奶头视频| 午夜精品久久久久久毛片777| 精品少妇一区二区三区视频日本电影| 一边摸一边抽搐一进一小说| 久热爱精品视频在线9| 69av精品久久久久久| 国产精品亚洲美女久久久| 亚洲av日韩精品久久久久久密| 欧美一级a爱片免费观看看 | 亚洲av成人精品一区久久| 欧美+亚洲+日韩+国产| 国产亚洲av嫩草精品影院| 国产精品一区二区三区四区免费观看 | 婷婷丁香在线五月| 悠悠久久av| 一个人观看的视频www高清免费观看 | 九色国产91popny在线| 国产精品av视频在线免费观看| 国产v大片淫在线免费观看| 免费看a级黄色片| 国产成人av教育| 欧美一区二区精品小视频在线| 国产av又大| 亚洲激情在线av| 99在线视频只有这里精品首页| 亚洲成人免费电影在线观看| 久久性视频一级片| 色哟哟哟哟哟哟| 男插女下体视频免费在线播放| 三级男女做爰猛烈吃奶摸视频| 久久久久久久精品吃奶| 老熟妇乱子伦视频在线观看| 人成视频在线观看免费观看| 亚洲欧美日韩东京热| 制服丝袜大香蕉在线| 欧美乱码精品一区二区三区| 欧美性猛交╳xxx乱大交人| 9191精品国产免费久久| 欧美日韩亚洲综合一区二区三区_| 亚洲国产欧美网| 久久久国产精品麻豆| 国产99久久九九免费精品| 搡老岳熟女国产| 国产精华一区二区三区| 怎么达到女性高潮| 在线观看66精品国产| 欧美黑人欧美精品刺激| 日本撒尿小便嘘嘘汇集6| 精品高清国产在线一区| 日本 av在线| 国产麻豆成人av免费视频| 曰老女人黄片| 麻豆国产97在线/欧美 | 午夜视频精品福利| 夜夜爽天天搞| 人人妻人人看人人澡| 男女做爰动态图高潮gif福利片| 精品熟女少妇八av免费久了| 亚洲午夜精品一区,二区,三区| 九色成人免费人妻av| 色av中文字幕| 欧美+亚洲+日韩+国产| 一边摸一边抽搐一进一小说| 国产亚洲精品综合一区在线观看 | 国产高清videossex| 欧美黑人精品巨大| 欧美大码av| 宅男免费午夜| 免费在线观看日本一区| 九九热线精品视视频播放| 国产成人精品久久二区二区免费| 最近视频中文字幕2019在线8| 欧美高清成人免费视频www| 国产一区二区三区在线臀色熟女| 国产伦在线观看视频一区| 香蕉丝袜av| 国产伦一二天堂av在线观看| 亚洲七黄色美女视频| 色老头精品视频在线观看| 国产精品久久久久久人妻精品电影| 白带黄色成豆腐渣| 午夜免费激情av| 亚洲精品在线美女| 日本一区二区免费在线视频| 国产亚洲精品第一综合不卡| 人妻久久中文字幕网| 欧美av亚洲av综合av国产av| 蜜桃久久精品国产亚洲av| 99久久精品热视频| 亚洲在线自拍视频| 日韩三级视频一区二区三区| 亚洲中文日韩欧美视频| 亚洲九九香蕉| 免费无遮挡裸体视频| 国产精品久久视频播放| 身体一侧抽搐| x7x7x7水蜜桃| 97人妻精品一区二区三区麻豆| 午夜福利成人在线免费观看| 91麻豆精品激情在线观看国产| 久久热在线av| 国产黄a三级三级三级人| 免费一级毛片在线播放高清视频| 俄罗斯特黄特色一大片| 婷婷亚洲欧美| 午夜福利视频1000在线观看| 老熟妇仑乱视频hdxx| 18禁黄网站禁片午夜丰满| 免费看a级黄色片| 欧美一级毛片孕妇| 亚洲精品久久成人aⅴ小说| www.熟女人妻精品国产| АⅤ资源中文在线天堂| 波多野结衣巨乳人妻| 丝袜美腿诱惑在线| 白带黄色成豆腐渣| 国产精华一区二区三区| 97碰自拍视频| 特大巨黑吊av在线直播| 在线永久观看黄色视频| 中文字幕久久专区| 欧美性猛交╳xxx乱大交人| 亚洲最大成人中文| 18禁裸乳无遮挡免费网站照片| 999精品在线视频| 热99re8久久精品国产| 村上凉子中文字幕在线| 国产成人啪精品午夜网站| 在线免费观看的www视频| 看免费av毛片| 国产精品av久久久久免费| 少妇熟女aⅴ在线视频| 麻豆av在线久日| 亚洲精品中文字幕在线视频| aaaaa片日本免费| www日本黄色视频网| 在线播放国产精品三级| 19禁男女啪啪无遮挡网站| 搡老熟女国产l中国老女人| 美女大奶头视频| 久久久国产成人精品二区| 国产不卡一卡二| 久久人妻av系列| av有码第一页| 老司机深夜福利视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 999精品在线视频| 黄色片一级片一级黄色片| 国产三级在线视频| 色尼玛亚洲综合影院| 天堂动漫精品| 一本综合久久免费| 亚洲精品久久国产高清桃花| ponron亚洲| 性欧美人与动物交配| 啦啦啦免费观看视频1| 欧美日韩亚洲国产一区二区在线观看| 1024香蕉在线观看| 看黄色毛片网站| 欧美成人性av电影在线观看| 全区人妻精品视频| 狂野欧美白嫩少妇大欣赏| 国产成人系列免费观看| 精品久久久久久久久久久久久| 老汉色av国产亚洲站长工具| 特级一级黄色大片| 一本久久中文字幕| 午夜福利在线在线| 精品久久久久久成人av| 长腿黑丝高跟| 亚洲一区二区三区不卡视频| 这个男人来自地球电影免费观看| 好男人电影高清在线观看| 久久精品国产99精品国产亚洲性色| 日韩欧美一区二区三区在线观看| 国产精品野战在线观看| 99国产精品一区二区蜜桃av| 人妻夜夜爽99麻豆av| 亚洲免费av在线视频| 亚洲片人在线观看| 亚洲精品av麻豆狂野| 18禁美女被吸乳视频| 夜夜看夜夜爽夜夜摸| 亚洲人成伊人成综合网2020| 国产精品av视频在线免费观看| 国产探花在线观看一区二区| 18禁美女被吸乳视频| 亚洲国产高清在线一区二区三| 99热6这里只有精品| 色噜噜av男人的天堂激情| 欧美一级a爱片免费观看看 | 夜夜看夜夜爽夜夜摸| 精品一区二区三区四区五区乱码| 国产一区二区在线观看日韩 | 成人手机av| 最好的美女福利视频网| 欧美大码av| 亚洲狠狠婷婷综合久久图片| 在线观看66精品国产| 黄色片一级片一级黄色片| 在线观看一区二区三区| 国产探花在线观看一区二区| 最近最新中文字幕大全免费视频| 999久久久精品免费观看国产| 老司机午夜十八禁免费视频| 免费av毛片视频| 麻豆一二三区av精品| 99久久无色码亚洲精品果冻| 18禁观看日本| 欧美中文综合在线视频| 久久精品综合一区二区三区| 国产精品99久久99久久久不卡| 别揉我奶头~嗯~啊~动态视频| 天天躁狠狠躁夜夜躁狠狠躁| 黑人欧美特级aaaaaa片| 欧美日韩乱码在线| 国产一区二区三区视频了| 欧美中文日本在线观看视频| 深夜精品福利| 国产精品亚洲美女久久久| 黄片大片在线免费观看| 日韩欧美精品v在线| 国产又黄又爽又无遮挡在线| 欧美日韩亚洲综合一区二区三区_| 熟女电影av网| 成人国语在线视频| 国产精品一区二区三区四区久久| 国产精品乱码一区二三区的特点| 欧美成狂野欧美在线观看| 少妇人妻一区二区三区视频| 国产在线精品亚洲第一网站| 性色av乱码一区二区三区2| 亚洲中文字幕一区二区三区有码在线看 | 亚洲av成人一区二区三| 一边摸一边抽搐一进一小说| 国产三级在线视频| 亚洲av成人精品一区久久| 精品国产乱码久久久久久男人| 精品国产亚洲在线| 精品午夜福利视频在线观看一区| 国产黄a三级三级三级人| 曰老女人黄片| 国产成人一区二区三区免费视频网站| 99久久精品热视频| 在线观看舔阴道视频| 男人的好看免费观看在线视频 | 欧美大码av| 香蕉丝袜av| 999精品在线视频| 欧美极品一区二区三区四区| 亚洲国产欧美人成| 久久国产精品影院| 久久久久久亚洲精品国产蜜桃av| 亚洲,欧美精品.| 国产成+人综合+亚洲专区| 亚洲五月婷婷丁香| 巨乳人妻的诱惑在线观看| 丁香六月欧美| 国模一区二区三区四区视频 | 久久中文字幕人妻熟女| 动漫黄色视频在线观看| 色综合欧美亚洲国产小说| 国产精品日韩av在线免费观看| 亚洲专区国产一区二区| 最新在线观看一区二区三区| 国产欧美日韩精品亚洲av| 成年版毛片免费区| 国产日本99.免费观看| 亚洲,欧美精品.| 少妇熟女aⅴ在线视频| 最近在线观看免费完整版| www日本黄色视频网| 丁香六月欧美| 久久久久国内视频| 久久人人精品亚洲av| 法律面前人人平等表现在哪些方面| 美女黄网站色视频| netflix在线观看网站| 香蕉国产在线看| 国产主播在线观看一区二区| 欧美乱码精品一区二区三区| 亚洲欧美日韩高清专用| 欧美黑人精品巨大| 香蕉久久夜色| 久9热在线精品视频| 两个人免费观看高清视频| 国产午夜精品论理片| 真人一进一出gif抽搐免费| 久久精品aⅴ一区二区三区四区| 亚洲午夜理论影院| 国产欧美日韩一区二区精品| 国产一区二区在线观看日韩 | 欧美成狂野欧美在线观看| avwww免费| 欧美乱码精品一区二区三区| 久久久国产成人免费| 丰满的人妻完整版| 99国产精品一区二区三区| 狂野欧美白嫩少妇大欣赏| 亚洲男人天堂网一区| 国语自产精品视频在线第100页| 国产69精品久久久久777片 |