• <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

    亚洲成人精品中文字幕电影| 又粗又硬又长又爽又黄的视频 | 国产精品蜜桃在线观看 | 国产一区亚洲一区在线观看| 国产精品久久视频播放| 亚洲人成网站在线播| 久久久欧美国产精品| 寂寞人妻少妇视频99o| 国产不卡一卡二| 欧美一区二区亚洲| 99在线视频只有这里精品首页| 午夜久久久久精精品| 三级男女做爰猛烈吃奶摸视频| 亚洲成人精品中文字幕电影| 校园人妻丝袜中文字幕| 国产日本99.免费观看| 99精品在免费线老司机午夜| 国产熟女欧美一区二区| 色综合亚洲欧美另类图片| 日日摸夜夜添夜夜爱| 日本黄色片子视频| 亚洲无线观看免费| 只有这里有精品99| 亚洲最大成人av| 简卡轻食公司| 国产精品电影一区二区三区| 国产乱人视频| 一级二级三级毛片免费看| 欧美xxxx性猛交bbbb| 尤物成人国产欧美一区二区三区| 国产亚洲精品久久久久久毛片| 99热网站在线观看| 久久午夜福利片| 亚洲不卡免费看| www.av在线官网国产| 国产真实伦视频高清在线观看| 男女边吃奶边做爰视频| 国产成人aa在线观看| 青春草视频在线免费观看| 国产午夜福利久久久久久| 亚洲av男天堂| 日韩亚洲欧美综合| 日韩视频在线欧美| 成人毛片60女人毛片免费| 亚洲精品色激情综合| 久久久久九九精品影院| 国产综合懂色| 九九在线视频观看精品| 日日撸夜夜添| 1024手机看黄色片| 联通29元200g的流量卡| 午夜福利高清视频| 色哟哟·www| 能在线免费观看的黄片| 亚洲人成网站在线观看播放| 久久99蜜桃精品久久| 嫩草影院精品99| 99热全是精品| 亚洲欧美日韩东京热| 日韩中字成人| 一区二区三区免费毛片| 国产亚洲av片在线观看秒播厂 | 日韩欧美国产在线观看| 精品熟女少妇av免费看| 蜜臀久久99精品久久宅男| av在线观看视频网站免费| 国产精品福利在线免费观看| 久久久久久大精品| 最近手机中文字幕大全| 中文字幕免费在线视频6| 啦啦啦韩国在线观看视频| 日韩欧美在线乱码| 在线天堂最新版资源| 男女啪啪激烈高潮av片| 亚洲自拍偷在线| 一区福利在线观看| 国产一区二区亚洲精品在线观看| 最后的刺客免费高清国语| 99热这里只有是精品在线观看| 99久久成人亚洲精品观看| 亚洲av中文字字幕乱码综合| 一个人看视频在线观看www免费| 夫妻性生交免费视频一级片| 成年女人看的毛片在线观看| 日产精品乱码卡一卡2卡三| 日韩欧美 国产精品| 一级二级三级毛片免费看| 日本免费a在线| 国产高清有码在线观看视频| 国产亚洲精品久久久com| 免费看av在线观看网站| 亚洲美女搞黄在线观看| 午夜免费男女啪啪视频观看| 亚洲av二区三区四区| ponron亚洲| 99热网站在线观看| a级毛片a级免费在线| 男女边吃奶边做爰视频| 毛片一级片免费看久久久久| 人妻少妇偷人精品九色| 老司机福利观看| 天堂网av新在线| 国产伦在线观看视频一区| 白带黄色成豆腐渣| 91久久精品电影网| av国产免费在线观看| 哪里可以看免费的av片| 欧美在线一区亚洲| 午夜精品在线福利| 午夜福利成人在线免费观看| 国产午夜精品一二区理论片| 一个人看的www免费观看视频| 欧美日韩在线观看h| 最近中文字幕高清免费大全6| 国产精品久久电影中文字幕| av在线蜜桃| 97在线视频观看| 国产亚洲欧美98| 搞女人的毛片| 热99在线观看视频| 亚洲最大成人手机在线| 成人午夜精彩视频在线观看| 欧美日韩精品成人综合77777| 亚洲精品影视一区二区三区av| 久久精品国产亚洲av天美| 久久久精品大字幕| 国产高清不卡午夜福利| 中文精品一卡2卡3卡4更新| 国产精品一区www在线观看| 看片在线看免费视频| 午夜福利在线在线| 卡戴珊不雅视频在线播放| 久久人妻av系列| 精品人妻偷拍中文字幕| 综合色丁香网| 免费黄网站久久成人精品| 黄色日韩在线| 亚洲最大成人av| 欧美成人免费av一区二区三区| 日本成人三级电影网站| 色综合站精品国产| 人妻久久中文字幕网| 亚洲中文字幕日韩| 久久精品综合一区二区三区| 高清毛片免费看| 天堂中文最新版在线下载 | 国产精品一二三区在线看| 少妇高潮的动态图| 亚洲国产欧洲综合997久久,| 女的被弄到高潮叫床怎么办| 麻豆久久精品国产亚洲av| 国产精品爽爽va在线观看网站| 晚上一个人看的免费电影| 三级男女做爰猛烈吃奶摸视频| 欧美精品国产亚洲| 欧美色欧美亚洲另类二区| 成人午夜精彩视频在线观看| 国产成人freesex在线| 国产真实乱freesex| 内地一区二区视频在线| 欧美潮喷喷水| 欧美最新免费一区二区三区| 国产精品久久久久久精品电影| 久久午夜福利片| 午夜久久久久精精品| 亚洲成人av在线免费| 精品一区二区免费观看| 69人妻影院| 久久精品国产99精品国产亚洲性色| 精品少妇黑人巨大在线播放 | 悠悠久久av| 精品久久国产蜜桃| 国产麻豆成人av免费视频| 成人美女网站在线观看视频| 九九在线视频观看精品| 久久久午夜欧美精品| 天堂√8在线中文| 一区二区三区四区激情视频 | 日韩av不卡免费在线播放| 狠狠狠狠99中文字幕| 午夜久久久久精精品| 亚洲精品日韩在线中文字幕 | 99久久久亚洲精品蜜臀av| 免费看美女性在线毛片视频| 成人美女网站在线观看视频| av在线观看视频网站免费| 亚洲精品影视一区二区三区av| 国产精品电影一区二区三区| 色综合亚洲欧美另类图片| 国产麻豆成人av免费视频| 丝袜喷水一区| 日韩国内少妇激情av| 亚洲精品亚洲一区二区| 天天躁日日操中文字幕| 亚洲av中文字字幕乱码综合| 日韩制服骚丝袜av| 在线免费观看不下载黄p国产| 国产精品国产高清国产av| 人人妻人人澡欧美一区二区| 欧美日本亚洲视频在线播放| 国产一级毛片七仙女欲春2| 欧美最黄视频在线播放免费| 少妇的逼水好多| 国产精品一区二区性色av| 在线观看一区二区三区| a级毛片a级免费在线| 99热这里只有是精品在线观看| 一级毛片aaaaaa免费看小| 色哟哟哟哟哟哟| 亚洲国产欧美人成| 99热全是精品| 欧美高清成人免费视频www| 99久久久亚洲精品蜜臀av| 99久久精品热视频| 天美传媒精品一区二区| 国产亚洲欧美98| 在线播放无遮挡| 男人舔奶头视频| 亚洲国产欧美人成| 人妻少妇偷人精品九色| 大又大粗又爽又黄少妇毛片口| 欧美三级亚洲精品| 欧美bdsm另类| 久久这里有精品视频免费| 日日啪夜夜撸| 一区二区三区四区激情视频 | 熟妇人妻久久中文字幕3abv| 欧美又色又爽又黄视频| 狂野欧美激情性xxxx在线观看| 精品人妻熟女av久视频| 欧美变态另类bdsm刘玥| 国产成人91sexporn| 国产中年淑女户外野战色| 亚洲经典国产精华液单| 黄色视频,在线免费观看| av在线老鸭窝| 亚洲美女视频黄频| 桃色一区二区三区在线观看| 久久热精品热| 欧美色视频一区免费| 亚洲最大成人中文| 国产一区二区三区在线臀色熟女| 联通29元200g的流量卡| 久久婷婷人人爽人人干人人爱| 久久精品夜色国产| 亚洲国产高清在线一区二区三| 日韩成人av中文字幕在线观看| 联通29元200g的流量卡| 国产私拍福利视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 91久久精品电影网| 春色校园在线视频观看| 久久精品国产亚洲av涩爱 | 日本熟妇午夜| 欧美一区二区亚洲| 天堂av国产一区二区熟女人妻| 久久久久久伊人网av| 国产一级毛片在线| 欧美+亚洲+日韩+国产| 久久久久九九精品影院| 国产精品.久久久| 激情 狠狠 欧美| 久久久久久久亚洲中文字幕| 久久99热这里只有精品18| 久久99热6这里只有精品| av在线蜜桃| 亚洲精品影视一区二区三区av| 亚洲性久久影院| 国产成人精品婷婷| 亚洲第一电影网av| 婷婷色综合大香蕉| 久久精品国产亚洲网站| 国产精品综合久久久久久久免费| 能在线免费看毛片的网站| 久久精品国产亚洲av涩爱 | 国产精品.久久久| 欧美性猛交╳xxx乱大交人| 精品久久久久久久久久免费视频| 99热网站在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产自在天天线| a级毛色黄片| 国产一级毛片七仙女欲春2| 免费搜索国产男女视频| 中国美白少妇内射xxxbb| 哪里可以看免费的av片| 成人美女网站在线观看视频| 免费电影在线观看免费观看| 久久精品国产鲁丝片午夜精品| 国产一区二区三区av在线 | 搡女人真爽免费视频火全软件| 欧美丝袜亚洲另类| 中文字幕av成人在线电影| 亚洲va在线va天堂va国产| 少妇熟女aⅴ在线视频| 青春草亚洲视频在线观看| 一级毛片aaaaaa免费看小| 色吧在线观看| 高清毛片免费看| 免费观看a级毛片全部| 秋霞在线观看毛片| 久久热精品热| 最近2019中文字幕mv第一页| 午夜福利高清视频| 午夜精品在线福利| 全区人妻精品视频| 国产精品久久久久久av不卡| 国语自产精品视频在线第100页| 午夜亚洲福利在线播放| 欧美人与善性xxx| 亚洲最大成人av| 美女被艹到高潮喷水动态| 丰满人妻一区二区三区视频av| 免费黄网站久久成人精品| 美女大奶头视频| 少妇人妻精品综合一区二区 | 午夜老司机福利剧场| 色综合色国产| 夜夜爽天天搞| av福利片在线观看| av卡一久久| 伊人久久精品亚洲午夜| 偷拍熟女少妇极品色| a级一级毛片免费在线观看| 久久久成人免费电影| 特大巨黑吊av在线直播| 欧美日韩精品成人综合77777| 亚洲一区高清亚洲精品| 国产av在哪里看| 久久精品夜色国产| 神马国产精品三级电影在线观看| 一级黄片播放器| 日韩欧美 国产精品| 久久午夜亚洲精品久久| 亚洲第一区二区三区不卡| 一边摸一边抽搐一进一小说| 美女脱内裤让男人舔精品视频 | 麻豆成人av视频| 波多野结衣高清作品| 国产精品一二三区在线看| 熟女电影av网| 欧美成人一区二区免费高清观看| 人体艺术视频欧美日本| 国产精华一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 一边摸一边抽搐一进一小说| 全区人妻精品视频| 午夜精品在线福利| 悠悠久久av| 国产老妇伦熟女老妇高清| 亚洲国产高清在线一区二区三| 一本久久精品| 色哟哟哟哟哟哟| 在现免费观看毛片| 人妻夜夜爽99麻豆av| 久久久色成人| 在线观看av片永久免费下载| 最近的中文字幕免费完整| 22中文网久久字幕| 久久人人爽人人片av| 午夜福利在线在线| 女人被狂操c到高潮| 国产精品福利在线免费观看| 欧美成人精品欧美一级黄| 国产真实伦视频高清在线观看| 少妇熟女欧美另类| 亚洲欧美日韩无卡精品| 亚洲欧美精品综合久久99| 五月伊人婷婷丁香| 国产精品久久电影中文字幕| 美女国产视频在线观看| 丝袜喷水一区| 变态另类丝袜制服| 欧美变态另类bdsm刘玥| 深爱激情五月婷婷| 日韩国内少妇激情av| 国产高清有码在线观看视频| 国产精品1区2区在线观看.| 夜夜爽天天搞| 你懂的网址亚洲精品在线观看 | 亚洲av成人精品一区久久| 国产精品免费一区二区三区在线| 黄色欧美视频在线观看| 一级毛片电影观看 | 三级经典国产精品| 精品久久久久久久末码| 日本-黄色视频高清免费观看| 一本一本综合久久| 欧美+亚洲+日韩+国产| 免费在线观看成人毛片| 黄片无遮挡物在线观看| 99久国产av精品| 色噜噜av男人的天堂激情| 高清毛片免费看| 亚洲综合色惰| 99热精品在线国产| 噜噜噜噜噜久久久久久91| 美女被艹到高潮喷水动态| 丰满乱子伦码专区| 欧美激情国产日韩精品一区| 国产真实乱freesex| 天天躁夜夜躁狠狠久久av| 精品无人区乱码1区二区| 日本与韩国留学比较| 久久精品久久久久久久性| 成人三级黄色视频| 欧美高清性xxxxhd video| 久久婷婷人人爽人人干人人爱| 天堂影院成人在线观看| 亚洲国产欧美在线一区| 成人一区二区视频在线观看| 狂野欧美激情性xxxx在线观看| 亚洲性久久影院| 亚洲av熟女| 久久久久久久久久久免费av| 毛片女人毛片| 亚洲精品亚洲一区二区| 蜜桃亚洲精品一区二区三区| 热99在线观看视频| 黄色欧美视频在线观看| 最近视频中文字幕2019在线8| 国产伦理片在线播放av一区 | 中文字幕av在线有码专区| 深爱激情五月婷婷| 国语自产精品视频在线第100页| 亚洲va在线va天堂va国产| 麻豆久久精品国产亚洲av| 最近的中文字幕免费完整| 黄片无遮挡物在线观看| 亚洲国产精品久久男人天堂| 国产高清视频在线观看网站| 欧美变态另类bdsm刘玥| 亚洲七黄色美女视频| av天堂中文字幕网| 一本一本综合久久| 日韩av不卡免费在线播放| 美女高潮的动态| 男女下面进入的视频免费午夜| av黄色大香蕉| 成年女人永久免费观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 夜夜爽天天搞| 国产免费一级a男人的天堂| 国产成人aa在线观看| 黄片无遮挡物在线观看| 亚洲国产精品国产精品| 哪里可以看免费的av片| 久久久久久久久久成人| 亚洲婷婷狠狠爱综合网| 三级国产精品欧美在线观看| 久久99热6这里只有精品| 久久精品久久久久久噜噜老黄 | 亚洲国产高清在线一区二区三| 简卡轻食公司| 久久国产乱子免费精品| 老女人水多毛片| 男女边吃奶边做爰视频| 看免费成人av毛片| 午夜免费激情av| 欧美日韩一区二区视频在线观看视频在线 | 非洲黑人性xxxx精品又粗又长| 亚洲激情五月婷婷啪啪| 亚洲av不卡在线观看| 国产成人91sexporn| 国产探花在线观看一区二区| 综合色丁香网| 中文亚洲av片在线观看爽| 亚洲成人久久爱视频| 亚洲国产精品sss在线观看| 九九爱精品视频在线观看| 亚洲自偷自拍三级| 亚洲无线观看免费| 在线观看午夜福利视频| 天天躁夜夜躁狠狠久久av| 亚洲精品456在线播放app| 久久久久久伊人网av| 久久这里有精品视频免费| 日韩成人av中文字幕在线观看| 一级黄色大片毛片| 亚洲国产色片| 精华霜和精华液先用哪个| 老师上课跳d突然被开到最大视频| 国产av不卡久久| 日日摸夜夜添夜夜添av毛片| 欧美xxxx黑人xx丫x性爽| 成人毛片a级毛片在线播放| 色视频www国产| www.色视频.com| 精品国内亚洲2022精品成人| 一个人免费在线观看电影| 亚洲成人av在线免费| 熟女人妻精品中文字幕| 国产一区二区三区av在线 | a级毛片a级免费在线| 岛国毛片在线播放| 国产淫片久久久久久久久| 国产午夜精品久久久久久一区二区三区| 国内久久婷婷六月综合欲色啪| 国产成人精品婷婷| 伊人久久精品亚洲午夜| 日本欧美国产在线视频| 欧美日本亚洲视频在线播放| 国产大屁股一区二区在线视频| 国产成人福利小说| 亚洲最大成人中文| 亚洲欧美日韩高清专用| 国产一区二区三区av在线 | 舔av片在线| av福利片在线观看| 在线免费观看不下载黄p国产| 国内揄拍国产精品人妻在线| 国产探花在线观看一区二区| 两个人视频免费观看高清| 日本黄大片高清| 色视频www国产| 我的女老师完整版在线观看| 国产亚洲精品av在线| 成人无遮挡网站| 成人毛片60女人毛片免费| 蜜桃亚洲精品一区二区三区| 国产高清视频在线观看网站| 国产成人freesex在线| 国产精华一区二区三区| 国产一区二区三区av在线 | 欧美一区二区亚洲| 亚洲国产高清在线一区二区三| 成年免费大片在线观看| 在线免费观看的www视频| 亚洲国产日韩欧美精品在线观看| 成人无遮挡网站| 亚洲精品久久久久久婷婷小说 | 精品久久久久久久人妻蜜臀av| 免费搜索国产男女视频| 久久精品人妻少妇| 久久久久久久久久久丰满| or卡值多少钱| 国产精品久久久久久久电影| 国产视频内射| 亚洲精品色激情综合| 内射极品少妇av片p| 久久午夜亚洲精品久久| 欧美日韩精品成人综合77777| 联通29元200g的流量卡| 卡戴珊不雅视频在线播放| 日本免费一区二区三区高清不卡| 一卡2卡三卡四卡精品乱码亚洲| 精品久久久久久久久av| 不卡视频在线观看欧美| 亚洲国产精品成人久久小说 | 免费人成在线观看视频色| 欧美zozozo另类| 边亲边吃奶的免费视频| 九九爱精品视频在线观看| 亚洲无线在线观看| 国产老妇女一区| 久久精品久久久久久久性| 国产成年人精品一区二区| 欧美激情国产日韩精品一区| 99在线视频只有这里精品首页| 久久精品国产自在天天线| 看片在线看免费视频| 欧美成人一区二区免费高清观看| 亚洲乱码一区二区免费版| 久久久a久久爽久久v久久| 免费大片18禁| 国产女主播在线喷水免费视频网站 | 综合色丁香网| 99久久久亚洲精品蜜臀av| 高清毛片免费看| 中国国产av一级| 亚洲第一区二区三区不卡| 特大巨黑吊av在线直播| 色综合站精品国产| 插逼视频在线观看| 久久久国产成人精品二区| 波多野结衣高清作品| 欧美不卡视频在线免费观看| 美女 人体艺术 gogo| 丰满人妻一区二区三区视频av| 久久精品国产亚洲av涩爱 | eeuss影院久久| 少妇的逼水好多| 午夜免费激情av| 欧美日韩精品成人综合77777| 麻豆一二三区av精品| 插阴视频在线观看视频| 97人妻精品一区二区三区麻豆| 国内精品久久久久精免费| 色播亚洲综合网| 亚洲真实伦在线观看| 日韩视频在线欧美| 亚洲婷婷狠狠爱综合网| 亚洲成人精品中文字幕电影| 啦啦啦韩国在线观看视频| 国产精品久久视频播放| 国产精品无大码| 国产av麻豆久久久久久久| 能在线免费看毛片的网站| 日本色播在线视频| 精品久久久久久久人妻蜜臀av| 午夜免费激情av| 色噜噜av男人的天堂激情| 亚洲欧美成人精品一区二区| 日本欧美国产在线视频| 亚洲欧美成人综合另类久久久 | 亚洲国产精品sss在线观看| 国产精品野战在线观看| 成人无遮挡网站| av国产免费在线观看| 一级毛片aaaaaa免费看小| 性色avwww在线观看| 精华霜和精华液先用哪个| 日韩欧美在线乱码| 亚洲欧美日韩高清专用| 丰满的人妻完整版|