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

    彈道導(dǎo)彈發(fā)射點(diǎn)位置的彈道重構(gòu)尋優(yōu)估計(jì)

    2018-10-15 03:47:02劉新建胡佳鑫楊樂平
    宇航總體技術(shù) 2018年5期
    關(guān)鍵詞:發(fā)射點(diǎn)彈道導(dǎo)彈彈道

    劉新建,胡佳鑫,楊樂平

    (國防科技大學(xué)空天科學(xué)學(xué)院,長沙410073)

    0 引言

    彈道落點(diǎn)預(yù)報(bào)是彈道導(dǎo)彈預(yù)警和反導(dǎo)中的重要任務(wù),但彈道導(dǎo)彈發(fā)射陣地的位置估計(jì)在彈道導(dǎo)彈防御中也具有重要作用。對于戰(zhàn)略彈道導(dǎo)彈,要準(zhǔn)確判定是哪個(gè)敵對國家進(jìn)行的攻擊,何時(shí)何地發(fā)射的,以便迅速進(jìn)行反擊;對于戰(zhàn)術(shù)彈道導(dǎo)彈,由于發(fā)射也需要調(diào)平和方位瞄準(zhǔn),一般是事先安排了有陣地依托的發(fā)射,也要盡可能準(zhǔn)確知道敵方是何時(shí)何地發(fā)射的,以便告知偵察衛(wèi)星對該區(qū)域掃描拍照,通過衛(wèi)星圖像進(jìn)一步精確鎖定坐標(biāo),以便轟炸發(fā)射陣地或彈道導(dǎo)彈貯存基地。

    Danish[1]較早研究了天基測量戰(zhàn)術(shù)彈道導(dǎo)彈(TBM)主動段彈道的發(fā)射點(diǎn)參數(shù)估計(jì)問題;Isaacson等[2]研究了戰(zhàn)術(shù)彈道導(dǎo)彈的狀態(tài)估計(jì)和軌跡預(yù)報(bào)。國內(nèi)李盾等[3]基于同步軌道預(yù)警衛(wèi)星對助推段彈道導(dǎo)彈的測量信息,利用改進(jìn)的高斯-牛頓法探討了發(fā)射點(diǎn)參數(shù)的最大似然估計(jì),但必須利用主動段目標(biāo)彈道曲線的先驗(yàn)信息;宋偉等[4]、陽曙光等[5]研究了雙星交匯測量戰(zhàn)術(shù)彈道導(dǎo)彈主動段的彈道參數(shù),并利用多項(xiàng)式擬合進(jìn)行彈道外推估計(jì)發(fā)射點(diǎn)參數(shù);王瑞等[6]研究了天基紅外低軌衛(wèi)星系統(tǒng)對自由段空間目標(biāo)定位的最小二乘、卡爾曼跟蹤濾波及自適應(yīng)濾波的方法。但鮮見基于對大氣層外被動段彈道的短弧測量對發(fā)射點(diǎn)位置進(jìn)行估計(jì)的研究,而且這種多項(xiàng)式彈道擬合外推發(fā)射點(diǎn)的方法較為簡單,對TBM主動段可能有效,但對中遠(yuǎn)程、洲際彈道導(dǎo)彈的被動段測量信息反推發(fā)射點(diǎn)未必有效。

    發(fā)射點(diǎn)的位置確定比彈道的落點(diǎn)預(yù)報(bào)更難,前者是逆問題;后者是正問題,相對要簡單。

    本文基于天基或地基設(shè)備的被動段短弧測量數(shù)據(jù),假設(shè)在WGS84系中,得到了來襲彈道導(dǎo)彈大氣層外的部分彈道,提出發(fā)射陣地位置估計(jì)方法。

    1 濾波定軌

    (1)坐標(biāo)參數(shù)轉(zhuǎn)換

    假設(shè)ti∶xi,yi,zi,vxi,vyi,vzi(i=1,2,…,n)為彈道導(dǎo)彈在真空中的一段短弧軌跡測量數(shù)據(jù),分別從WGS84系轉(zhuǎn)換到某一地心慣性系,如J2000 系,得到ti∶Xi,Yi,Zi,Vxi,Vyi,Vzi。

    其中,ti為飛行時(shí)刻,Xi、Yi、Zi為彈道導(dǎo)彈在該時(shí)刻的3個(gè)直角坐標(biāo)位置分量,Vxi、Vyi、Vzi為該時(shí)刻的3個(gè)速度分量。

    (2)均值濾波去噪

    由于每個(gè)測量值含有噪聲,并非是彈道導(dǎo)彈的真實(shí)狀態(tài),因此需要去噪濾波。工程上有多種數(shù)字濾波方法,如卡爾曼濾波、最小二乘濾波、自適應(yīng)濾波、平均值濾波。特別對于地基雷達(dá),由于地球曲率的影響,可見弧段極短,而彈道導(dǎo)彈飛行速度極快,獲得的測量數(shù)據(jù)不會很多,因此可采用單點(diǎn)均值濾波。均值濾波最簡單,沒有延遲,對少量的短弧測量數(shù)據(jù)魯棒可靠。

    所謂均值濾波,就是依次利用測量數(shù)據(jù)ti:Xi,Yi,Zi,Vxi,Vyi,Vzi(i=1,2,…,n),根據(jù)真空段自由飛行動力學(xué)模型式 (2),進(jìn)行數(shù)值積分,獲得測量數(shù)據(jù)段之后任選的定軌時(shí)刻tf的n個(gè)狀態(tài)數(shù)據(jù)。如果測量弧段的數(shù)據(jù)沒有噪聲,那么依次積分,預(yù)報(bào)tf時(shí)刻的六分量狀態(tài)應(yīng)該是完全相同的;但正因?yàn)榇嬖跍y量噪聲,所以依次積分獲得該時(shí)刻的狀態(tài)也存在噪聲,那么求平均值。

    可是tf時(shí)刻的彈道導(dǎo)彈真實(shí)狀態(tài)只有一個(gè),如果測量數(shù)據(jù)的噪聲遵循白噪聲分布,而且采樣點(diǎn)足夠多,則平均值逼近真值,就得到了彈道導(dǎo)彈真實(shí)狀態(tài)tf:Xf,Yf,Zf,Vxf,Vyf,Vzf。

    (3)根據(jù)如下地慣系中空間飛行器的動力學(xué)

    微分方程組[7]

    其中,地球扁率J2=1.08263×10-3, 地球引力常數(shù)μ =3.986005×1014m3/s2, 地球長半軸=6378140m,地心矢徑再由上述定軌點(diǎn)的真實(shí)狀態(tài),進(jìn)行彈道積分預(yù)報(bào)再入落點(diǎn)e的直角坐標(biāo)位置和速度Xe、Ye、Ze、Vxe、Vye、Vze。

    2 發(fā)射點(diǎn)初始解析估計(jì)

    為獲得初始估計(jì),不考慮彈道導(dǎo)彈助推段彈道及再入段大氣阻力的影響,假設(shè)導(dǎo)彈在地面瞬間獲得足夠能量,滿足測量彈道所預(yù)報(bào)落點(diǎn)的位置和速度要求,不妨稱為脈沖發(fā)射點(diǎn)。脈沖發(fā)射點(diǎn)的位置估計(jì)是發(fā)射陣地位置的近似估計(jì),是彈道重構(gòu)尋優(yōu)估計(jì)的基礎(chǔ)。

    圖1 二體脈沖彈道示意圖Fig.1 Impulse ballistic sketch of two-body problem

    1)如圖1所示,根據(jù)再入點(diǎn)e的直角坐標(biāo)和速度Xe、Ye、Ze、Vxe、Vye、Vze, 計(jì)算經(jīng)典軌道根數(shù)a、e、i、Ω、ω、M,并求其真近點(diǎn)角fe[8]。

    2)求二體脈沖發(fā)射點(diǎn)的真近點(diǎn)角fk

    于是不難求得其對應(yīng)的平近點(diǎn)角Mk和偏近點(diǎn)角Ek。

    3)求二體飛行時(shí)間

    其中,a是定軌得到的彈道長半軸,平近點(diǎn)角Me=M。

    4)求脈沖發(fā)射點(diǎn)的經(jīng)緯度

    理論上,這是一個(gè)已知再入終點(diǎn)位置和飛行時(shí)間,反求初始點(diǎn)位置的逆向動力學(xué)問題。這里提出利用二體運(yùn)動軌道的封閉性特點(diǎn),避開逆向積分問題。

    根據(jù)二體運(yùn)動理論,二體脈沖發(fā)射點(diǎn)的直角坐標(biāo)有如下解析公式[9]

    再將其從地慣系轉(zhuǎn)換到WGS84系的xk、yk、zk,由此計(jì)算脈沖發(fā)射點(diǎn)的經(jīng)緯度λk、фk,就是發(fā)射點(diǎn)的初步估計(jì),誤差與主動段的彈道設(shè)計(jì)程序及參數(shù)有關(guān),對于射程超過10000km的洲際彈道導(dǎo)彈,數(shù)值仿真表明本方法求得發(fā)射點(diǎn)最大誤差為200km。

    初始解析估計(jì)的計(jì)算量很小,適合于彈道導(dǎo)彈防御系統(tǒng)實(shí)時(shí)快速預(yù)報(bào),至少能準(zhǔn)確告知彈道導(dǎo)彈來自哪個(gè)敵對國家。

    3 彈道重構(gòu)尋優(yōu)估計(jì)

    由于脈沖點(diǎn)的估計(jì)與發(fā)射陣地的真實(shí)位置存在誤差,需要進(jìn)一步修正。策略是在脈沖發(fā)射點(diǎn)附近進(jìn)行系列彈道重構(gòu),仿真模擬彈道導(dǎo)彈從發(fā)射點(diǎn)垂直起飛、主動段程序轉(zhuǎn)彎、末級瞄準(zhǔn)與閉路控制關(guān)機(jī),優(yōu)選出滿足定軌預(yù)報(bào)再入點(diǎn)約束的重構(gòu)彈道,輸出所對應(yīng)的發(fā)射點(diǎn)經(jīng)緯度。

    3.1 尋優(yōu)方法

    1)以脈沖發(fā)射點(diǎn)及對應(yīng)的自由段飛行時(shí)間為中心,設(shè)定經(jīng)緯度誤差邊界范圍,如±2°(對應(yīng)440km邊長的正方形),自由段時(shí)間偏差,如±100s,劃分三維網(wǎng)格如圖2所示,遍歷搜索可行解與最優(yōu)解,發(fā)射點(diǎn)位置精度與網(wǎng)格大小有關(guān),顯然網(wǎng)格越細(xì),搜索的時(shí)間越長。

    2)根據(jù)脈沖點(diǎn)與再入落點(diǎn)的經(jīng)緯度,由球面三角公式估計(jì)彈道射程。

    3)按照射程大小 (洲際、遠(yuǎn)程、中程、近程和短程),參考對應(yīng)的國內(nèi)外典型彈道導(dǎo)彈的參數(shù),虛擬來襲彈道導(dǎo)彈的動力學(xué)仿真參數(shù),如級數(shù)、各級的推力、質(zhì)量、秒流量、燃燒時(shí)間以及氣動阻力、升力系數(shù)。

    4)按照彈道導(dǎo)彈主動段的彈道設(shè)計(jì)方法,設(shè)計(jì)轉(zhuǎn)彎程序如攻角函數(shù)、各級轉(zhuǎn)彎角速率[8]。

    圖2 脈沖發(fā)射位置點(diǎn)附近的三維網(wǎng)格Fig.2 Three dimensional grid near impulse launch location

    5)為確保再入位置約束,導(dǎo)彈末級采用Lambert制導(dǎo)方法[7],根據(jù)末級當(dāng)前位置矢量,再入落點(diǎn)位置矢量和飛行時(shí)間要求,求取需要速度,利用速度增益制導(dǎo)構(gòu)成閉路關(guān)機(jī)控制,求得導(dǎo)彈末級飛行的俯仰角、偏航角和關(guān)機(jī)時(shí)間。

    3.2 約束條件

    以上給出了主要搜索變量的邊界約束和方法,但要模擬真實(shí)彈道,還需加入彈道重構(gòu)約束條件,即仿真彈道的再入點(diǎn)位置和速度六分量應(yīng)與彈道測量后的定軌值相同或逼近。

    再入點(diǎn)位置誤差約束:Δλ≤ελ,Δφ≤εφ,ΔH ≤εH。

    即經(jīng)度、緯度和高程三分量仿真值與定軌值的誤差小于設(shè)定值,因?yàn)槟┘壊捎昧薒ambert制導(dǎo)方法需要速度閉路制導(dǎo),仿真會自動滿足,簡化了位置約束的求解。

    再入角度誤差約束:ΔΘ≤εΘ,Δσ≤εσ,即速度的傾角和偏角誤差小于設(shè)定值

    3.3 目標(biāo)函數(shù)

    取再入速度模的誤差ΔVe為目標(biāo)函數(shù),應(yīng)取最小。

    4 仿真驗(yàn)證

    這里以某三級彈道導(dǎo)彈A為例,給出數(shù)學(xué)仿真予以說明。

    (1)導(dǎo)彈A的主要參數(shù)[10]

    一級起飛質(zhì)量為34.5t、推力為912kN、工作時(shí)間為61s、推進(jìn)劑秒流量為360kg/s;

    二級起飛質(zhì)量為12t、推力為270kN、工作時(shí)間為65s、推進(jìn)劑秒流量為108kg/s;

    三級起飛質(zhì)量為5t、推力為155kN、工作時(shí)間為59s、推進(jìn)劑秒流量為62kg/s。

    末修級工作時(shí)間由閉路關(guān)機(jī)控制指令決定。但彈道數(shù)值仿真時(shí),為簡化起見,可不設(shè)末修級,第三級工作時(shí)間不固定,由閉路關(guān)機(jī)仿真時(shí)間決定,只要滿足自由段飛行時(shí)間約束和落點(diǎn)位置約束就可以。其他中短程的戰(zhàn)術(shù)彈道導(dǎo)彈類似,也可不設(shè)末修級。

    (2)發(fā)射點(diǎn)和目標(biāo)點(diǎn)參數(shù)假設(shè)

    發(fā)射點(diǎn):經(jīng)度-118°,緯度30°;

    目標(biāo)點(diǎn):經(jīng)度112°,緯度34°;

    兩點(diǎn)之間的射程為11643km。

    (3)仿真過程

    1)在脈沖估計(jì)的基礎(chǔ)上,編寫3自由度仿真優(yōu)化搜索程序。

    2)設(shè)計(jì)一條從發(fā)射點(diǎn)到目標(biāo)點(diǎn)的標(biāo)稱攻擊彈道,獲得中段N個(gè)短弧采樣數(shù)據(jù),例如50個(gè)軌跡點(diǎn),每隔1秒一個(gè)采樣數(shù)據(jù),適當(dāng)加入白噪聲,基于這N個(gè)數(shù)據(jù),完成定軌計(jì)算。

    3)基于以上方法,完成發(fā)射點(diǎn)的二體脈沖估計(jì)和彈道重構(gòu)尋優(yōu)估計(jì)。

    (4)仿真結(jié)果

    二體脈沖估計(jì):

    發(fā)射點(diǎn)位置誤差

    經(jīng)度差:-0.127°(對應(yīng)14km)

    緯度差:1.0015°(對應(yīng)110km)彈道重構(gòu)尋優(yōu)估計(jì):

    發(fā)射點(diǎn)位置誤差

    經(jīng)度差:0.125°(對應(yīng)14km)

    緯度差:0.0015° (對應(yīng)0.2km)

    (5)約束條件逼近度

    速度誤差:0.047m/s

    當(dāng)?shù)厮俣葍A角誤差:0.85°

    速度偏角誤差:0.027°

    經(jīng)度差:0.0045° (對應(yīng)0.5km誤差)

    緯度差:0.0044° (對應(yīng)0.5km誤差)

    圖3給出了上述數(shù)值仿真算例的尋優(yōu)估計(jì)重構(gòu)彈道。

    圖3 彈道重構(gòu)尋優(yōu)估計(jì)數(shù)值仿真軌跡示意圖Fig.3 Numerical simulation sketch of optimal estimation with trajectory regeneration

    可見,對于本算例的二體快速脈沖估計(jì),發(fā)射點(diǎn)的位置誤差為120km,這對于判斷來襲彈道導(dǎo)彈的敵對國家是夠用的。而經(jīng)過彈道重構(gòu)的尋優(yōu)估計(jì),發(fā)射點(diǎn)位置誤差估計(jì)降到了15km左右,提高了可觀的估計(jì)精度,如果發(fā)射點(diǎn)尋優(yōu)網(wǎng)格進(jìn)一步細(xì)分,誤差還可能降低。

    5 結(jié)論

    提出的彈道導(dǎo)彈發(fā)射參數(shù)估計(jì)方法包括脈沖初始估計(jì)和彈道重構(gòu)仿真兩部分,前者實(shí)時(shí)性好,也是后者的基礎(chǔ)。后者以脈沖點(diǎn)為估計(jì)中心,對于戰(zhàn)略彈道導(dǎo)彈,以經(jīng)緯度±2°誤差為位置邊界,以二體飛行時(shí)間tke為中心,±100s為自由段的飛行時(shí)間邊界,劃分三維網(wǎng)格,并以再入點(diǎn)的位置三分量、當(dāng)?shù)厮俣葍A角、速度偏角共5個(gè)彈道重構(gòu)約束條件,以速度誤差為目標(biāo)函數(shù),進(jìn)行三自由度彈道控制仿真與優(yōu)化搜索,求取發(fā)射陣地坐標(biāo),這樣來減小發(fā)射點(diǎn)位置估計(jì)誤差。對于戰(zhàn)術(shù)彈道導(dǎo)彈,網(wǎng)格邊界可通過仿真進(jìn)一步縮小,方法類似。

    這種彈道重構(gòu)方法,不需要完全一致的彈道導(dǎo)彈參數(shù),只要大致接近就可,因?yàn)槭前凑諒椀缹?dǎo)彈攻擊目標(biāo)所采用的典型彈道規(guī)劃方法,基于動力學(xué)與控制仿真,并通過優(yōu)化搜索參數(shù)和Lambert末制導(dǎo)來滿足再入邊界條件的,使得數(shù)學(xué)仿真彈道能夠逼近實(shí)際飛行軌跡。

    只要獲得了來襲彈道測量參數(shù),本方法可以獲得較高置信度的估計(jì),在彈道導(dǎo)彈預(yù)警和防御中有一定參考價(jià)值。

    猜你喜歡
    發(fā)射點(diǎn)彈道導(dǎo)彈彈道
    飛行器虛擬實(shí)驗(yàn)平臺開發(fā)及其在教學(xué)中的應(yīng)用
    讓人膽寒的“遠(yuǎn)程殺手”:彈道導(dǎo)彈
    軍事文摘(2022年20期)2023-01-10 07:18:38
    彈道——打勝仗的奧秘
    快響火箭機(jī)動發(fā)射點(diǎn)位/發(fā)射窗口一體規(guī)劃
    航天控制(2020年5期)2020-03-29 02:10:32
    一維彈道修正彈無線通信系統(tǒng)研制
    電子制作(2019年7期)2019-04-25 13:17:48
    基于分解協(xié)調(diào)法的彈道導(dǎo)彈突防效能控制優(yōu)化
    基于PID控制的二維彈道修正彈仿真
    全固態(tài)PDM 1kW中波廣播發(fā)射機(jī)的維護(hù)
    多場強(qiáng)下調(diào)頻同步廣播的調(diào)試技術(shù)探索與實(shí)踐
    聲屏世界(2015年13期)2015-12-29 05:01:16
    消除彈道跟蹤數(shù)據(jù)中伺服系統(tǒng)的振顫干擾
    亚洲综合精品二区| 亚洲在久久综合| 国产精品欧美亚洲77777| 又粗又硬又长又爽又黄的视频| 各种免费的搞黄视频| 亚洲精品乱久久久久久| 亚洲免费av在线视频| 国产精品三级大全| 国产日韩一区二区三区精品不卡| 久久97久久精品| 精品国产国语对白av| 曰老女人黄片| 欧美日韩精品网址| 国产无遮挡羞羞视频在线观看| 欧美日韩精品网址| 好男人视频免费观看在线| 午夜日韩欧美国产| 日本猛色少妇xxxxx猛交久久| 视频区图区小说| 免费人妻精品一区二区三区视频| 香蕉国产在线看| 超色免费av| 男女边摸边吃奶| 国产有黄有色有爽视频| 叶爱在线成人免费视频播放| 久久99一区二区三区| 国产极品天堂在线| 午夜老司机福利片| 国产亚洲欧美精品永久| 久久国产亚洲av麻豆专区| 丝袜在线中文字幕| 午夜福利免费观看在线| 七月丁香在线播放| 毛片一级片免费看久久久久| www.自偷自拍.com| 亚洲精品美女久久av网站| 在线观看免费日韩欧美大片| 校园人妻丝袜中文字幕| 纯流量卡能插随身wifi吗| 国产日韩欧美在线精品| 亚洲av综合色区一区| 侵犯人妻中文字幕一二三四区| a 毛片基地| 秋霞在线观看毛片| 久久人人97超碰香蕉20202| 久久人人97超碰香蕉20202| 99re6热这里在线精品视频| 欧美精品一区二区大全| 美女大奶头黄色视频| 亚洲色图综合在线观看| 高清av免费在线| 久久97久久精品| 亚洲国产精品成人久久小说| 综合色丁香网| 丝袜脚勾引网站| 亚洲精品国产av成人精品| 丰满饥渴人妻一区二区三| av线在线观看网站| 亚洲成人国产一区在线观看 | 亚洲人成网站在线观看播放| 亚洲第一区二区三区不卡| 婷婷色综合大香蕉| 18禁裸乳无遮挡动漫免费视频| 国产精品一国产av| 如日韩欧美国产精品一区二区三区| 久久这里只有精品19| 中文字幕高清在线视频| 搡老岳熟女国产| 亚洲国产精品一区二区三区在线| 亚洲精品国产一区二区精华液| 久久精品亚洲熟妇少妇任你| 国产午夜精品一二区理论片| 一级片免费观看大全| 欧美另类一区| 丝袜在线中文字幕| 免费观看人在逋| 18禁裸乳无遮挡动漫免费视频| 老司机深夜福利视频在线观看 | 丝袜在线中文字幕| 你懂的网址亚洲精品在线观看| 亚洲免费av在线视频| 亚洲自偷自拍图片 自拍| 天天躁狠狠躁夜夜躁狠狠躁| 欧美在线一区亚洲| 久久午夜综合久久蜜桃| 老汉色∧v一级毛片| 欧美日本中文国产一区发布| 桃花免费在线播放| 国产精品99久久99久久久不卡 | 女人精品久久久久毛片| 亚洲av国产av综合av卡| 最近最新中文字幕大全免费视频 | 亚洲国产av影院在线观看| 国产亚洲欧美精品永久| 我的亚洲天堂| 国产精品嫩草影院av在线观看| 黄网站色视频无遮挡免费观看| 久久久久久人妻| 悠悠久久av| 亚洲精品中文字幕在线视频| 亚洲精品中文字幕在线视频| 999精品在线视频| 毛片一级片免费看久久久久| 欧美日韩视频高清一区二区三区二| 欧美中文综合在线视频| 久热这里只有精品99| 国产毛片在线视频| 日本欧美视频一区| 国产精品一区二区在线不卡| 黄色毛片三级朝国网站| 国产片特级美女逼逼视频| 一级a爱视频在线免费观看| 伊人久久国产一区二区| 欧美黑人欧美精品刺激| 国产有黄有色有爽视频| 亚洲成av片中文字幕在线观看| 欧美日韩一区二区视频在线观看视频在线| 欧美精品人与动牲交sv欧美| 高清黄色对白视频在线免费看| 1024视频免费在线观看| 狂野欧美激情性xxxx| 亚洲欧洲精品一区二区精品久久久 | 国产精品免费大片| 国产黄色视频一区二区在线观看| 国产欧美亚洲国产| 一本—道久久a久久精品蜜桃钙片| 曰老女人黄片| 国产伦理片在线播放av一区| 亚洲欧美成人综合另类久久久| 99久久99久久久精品蜜桃| 国产野战对白在线观看| 国产黄色免费在线视频| 美女福利国产在线| 一区二区日韩欧美中文字幕| 成人手机av| 久久国产亚洲av麻豆专区| videos熟女内射| 满18在线观看网站| 午夜免费鲁丝| 欧美av亚洲av综合av国产av | 国产欧美日韩一区二区三区在线| 国产有黄有色有爽视频| 国产男女内射视频| 久久久久久人人人人人| 永久免费av网站大全| 精品一区在线观看国产| 国产深夜福利视频在线观看| 亚洲精品日韩在线中文字幕| 亚洲成人手机| 国产精品国产av在线观看| 久久99精品国语久久久| 综合色丁香网| 波多野结衣av一区二区av| 交换朋友夫妻互换小说| 91精品伊人久久大香线蕉| 波多野结衣av一区二区av| 久久免费观看电影| 中文字幕高清在线视频| 久久狼人影院| 欧美日韩福利视频一区二区| 成年人午夜在线观看视频| 蜜桃在线观看..| 日韩精品有码人妻一区| 亚洲欧美一区二区三区国产| 女人久久www免费人成看片| 一本—道久久a久久精品蜜桃钙片| 一二三四在线观看免费中文在| 99re6热这里在线精品视频| 国产成人91sexporn| 99久国产av精品国产电影| 亚洲视频免费观看视频| 亚洲久久久国产精品| 亚洲熟女精品中文字幕| 9191精品国产免费久久| 高清视频免费观看一区二区| 久久久久国产一级毛片高清牌| 在线免费观看不下载黄p国产| 看免费成人av毛片| 亚洲国产中文字幕在线视频| 天天躁夜夜躁狠狠躁躁| 久久人妻熟女aⅴ| 日本vs欧美在线观看视频| 亚洲免费av在线视频| 黄片播放在线免费| 成人国产av品久久久| 国产免费又黄又爽又色| 国产成人精品在线电影| 国产精品国产三级专区第一集| 日韩一卡2卡3卡4卡2021年| 高清黄色对白视频在线免费看| av又黄又爽大尺度在线免费看| 久久精品国产综合久久久| 国产精品一二三区在线看| 久久久久久久国产电影| 亚洲精品乱久久久久久| 久久久久精品性色| 久久狼人影院| 国产精品.久久久| a级毛片在线看网站| 视频区图区小说| 黑人猛操日本美女一级片| 国产视频首页在线观看| www.熟女人妻精品国产| 国产精品免费视频内射| 成人午夜精彩视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产伦人伦偷精品视频| 免费黄频网站在线观看国产| 蜜桃在线观看..| 亚洲av欧美aⅴ国产| 中文欧美无线码| 国产成人av激情在线播放| 丰满饥渴人妻一区二区三| 国产精品国产av在线观看| 亚洲欧美精品自产自拍| 少妇的丰满在线观看| av女优亚洲男人天堂| 成人午夜精彩视频在线观看| a级毛片在线看网站| 日韩av免费高清视频| 一二三四在线观看免费中文在| 精品少妇黑人巨大在线播放| 波多野结衣一区麻豆| 国产视频首页在线观看| 晚上一个人看的免费电影| 国产免费现黄频在线看| 亚洲美女视频黄频| av电影中文网址| 日韩制服丝袜自拍偷拍| 成年人免费黄色播放视频| 亚洲国产欧美一区二区综合| 精品午夜福利在线看| 日韩欧美精品免费久久| 久久精品国产a三级三级三级| 成人亚洲欧美一区二区av| 亚洲天堂av无毛| 久久 成人 亚洲| 久久人人97超碰香蕉20202| 久久精品国产亚洲av涩爱| 免费观看人在逋| 久久久久视频综合| 亚洲精品日本国产第一区| 精品亚洲成国产av| 国产熟女欧美一区二区| 亚洲七黄色美女视频| 丁香六月天网| 国产精品久久久av美女十八| 大香蕉久久网| 欧美xxⅹ黑人| 久久久久人妻精品一区果冻| 亚洲熟女毛片儿| 亚洲国产日韩一区二区| 中文字幕色久视频| 一边摸一边做爽爽视频免费| 亚洲国产日韩一区二区| 国产精品亚洲av一区麻豆 | 人妻一区二区av| e午夜精品久久久久久久| 成年人午夜在线观看视频| 精品一品国产午夜福利视频| 久久久精品免费免费高清| 亚洲精品自拍成人| 中文精品一卡2卡3卡4更新| 国产精品久久久久久精品电影小说| 老熟女久久久| 国产欧美日韩一区二区三区在线| 波野结衣二区三区在线| 午夜老司机福利片| 久久精品久久精品一区二区三区| 老司机靠b影院| 叶爱在线成人免费视频播放| 熟女少妇亚洲综合色aaa.| 麻豆精品久久久久久蜜桃| 久久亚洲国产成人精品v| 永久免费av网站大全| 亚洲久久久国产精品| 久久99精品国语久久久| 丰满乱子伦码专区| 国产亚洲一区二区精品| 老司机影院毛片| 欧美另类一区| 嫩草影视91久久| 美女脱内裤让男人舔精品视频| 国产精品一区二区在线不卡| 一级黄片播放器| 成人亚洲欧美一区二区av| 黑人猛操日本美女一级片| 美女视频免费永久观看网站| 波野结衣二区三区在线| 色婷婷av一区二区三区视频| 色婷婷久久久亚洲欧美| 亚洲精华国产精华液的使用体验| 免费观看a级毛片全部| 波多野结衣一区麻豆| 啦啦啦视频在线资源免费观看| 国产精品蜜桃在线观看| 性高湖久久久久久久久免费观看| 性少妇av在线| 99久国产av精品国产电影| 色播在线永久视频| 亚洲成人手机| 国产亚洲av片在线观看秒播厂| 男女边摸边吃奶| 亚洲av电影在线观看一区二区三区| 水蜜桃什么品种好| 中文字幕人妻丝袜一区二区 | 国产黄色免费在线视频| 精品久久蜜臀av无| 亚洲国产精品成人久久小说| 欧美在线黄色| 日韩av不卡免费在线播放| av国产精品久久久久影院| 国产成人系列免费观看| 亚洲天堂av无毛| 一边摸一边抽搐一进一出视频| 日韩欧美一区视频在线观看| 在线观看免费日韩欧美大片| 欧美精品av麻豆av| 日韩伦理黄色片| 日日摸夜夜添夜夜爱| 丝袜脚勾引网站| 国产成人欧美在线观看 | 国产av码专区亚洲av| 亚洲综合色网址| 悠悠久久av| 国产男女超爽视频在线观看| 性少妇av在线| 亚洲欧美激情在线| 老司机影院毛片| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧美中文字幕日韩二区| 国产野战对白在线观看| 五月开心婷婷网| 在线观看三级黄色| 亚洲成人一二三区av| 亚洲中文av在线| 亚洲国产精品一区二区三区在线| 搡老岳熟女国产| 成年美女黄网站色视频大全免费| 久久热在线av| 天天躁狠狠躁夜夜躁狠狠躁| 午夜日本视频在线| 中文精品一卡2卡3卡4更新| 亚洲国产成人一精品久久久| 永久免费av网站大全| 久久久久久久精品精品| 999久久久国产精品视频| 国产有黄有色有爽视频| 亚洲国产精品一区三区| 丝瓜视频免费看黄片| 这个男人来自地球电影免费观看 | 婷婷色av中文字幕| 欧美黑人精品巨大| avwww免费| 天天操日日干夜夜撸| 最近最新中文字幕大全免费视频 | 国产毛片在线视频| 亚洲精品乱久久久久久| 亚洲精品在线美女| 满18在线观看网站| 精品卡一卡二卡四卡免费| 赤兔流量卡办理| 爱豆传媒免费全集在线观看| 三上悠亚av全集在线观看| 国产福利在线免费观看视频| 黑人猛操日本美女一级片| 超碰97精品在线观看| 久久人妻熟女aⅴ| 制服人妻中文乱码| 国产激情久久老熟女| 99精品久久久久人妻精品| 啦啦啦在线免费观看视频4| 女人高潮潮喷娇喘18禁视频| 欧美最新免费一区二区三区| 777米奇影视久久| 亚洲欧美精品综合一区二区三区| 深夜精品福利| 伊人亚洲综合成人网| 视频区图区小说| 大话2 男鬼变身卡| 久久久久久久久免费视频了| 精品一区二区三卡| 一区二区av电影网| 蜜桃国产av成人99| 美女大奶头黄色视频| 久久韩国三级中文字幕| 亚洲,欧美,日韩| 成人三级做爰电影| 国产乱人偷精品视频| 欧美日韩视频精品一区| 下体分泌物呈黄色| 久久精品人人爽人人爽视色| 国产精品欧美亚洲77777| 十八禁高潮呻吟视频| 色精品久久人妻99蜜桃| 好男人视频免费观看在线| 久久毛片免费看一区二区三区| 日韩大片免费观看网站| 久久人人97超碰香蕉20202| 亚洲精品日韩在线中文字幕| 丁香六月欧美| 亚洲人成网站在线观看播放| 久久久久精品久久久久真实原创| 水蜜桃什么品种好| 成人亚洲欧美一区二区av| 欧美精品av麻豆av| 少妇精品久久久久久久| 精品久久久久久电影网| 亚洲av综合色区一区| 中文字幕人妻丝袜制服| www.精华液| 精品亚洲成国产av| 亚洲一码二码三码区别大吗| 十八禁网站网址无遮挡| 一级片'在线观看视频| 女的被弄到高潮叫床怎么办| 最近最新中文字幕大全免费视频 | 人人澡人人妻人| 人人妻,人人澡人人爽秒播 | 亚洲精品日韩在线中文字幕| 最近手机中文字幕大全| www.av在线官网国产| 伦理电影免费视频| 亚洲国产看品久久| 在线观看免费午夜福利视频| 精品国产乱码久久久久久男人| 青春草亚洲视频在线观看| 伊人久久国产一区二区| 一区二区三区精品91| 91精品国产国语对白视频| 久久天堂一区二区三区四区| 十分钟在线观看高清视频www| 亚洲熟女毛片儿| 国产精品久久久久久久久免| 少妇人妻精品综合一区二区| 国产精品二区激情视频| 亚洲国产欧美在线一区| 欧美 日韩 精品 国产| 大陆偷拍与自拍| 最近最新中文字幕免费大全7| 97人妻天天添夜夜摸| 精品国产露脸久久av麻豆| 看免费av毛片| 亚洲视频免费观看视频| 欧美日韩亚洲国产一区二区在线观看 | 9热在线视频观看99| 美女高潮到喷水免费观看| 日韩av免费高清视频| 精品亚洲乱码少妇综合久久| 国产 一区精品| 热99国产精品久久久久久7| 人人妻人人添人人爽欧美一区卜| 侵犯人妻中文字幕一二三四区| 国产成人免费无遮挡视频| 亚洲欧美色中文字幕在线| 久久这里只有精品19| www.自偷自拍.com| 人人妻人人澡人人看| 久久久精品区二区三区| 久久精品国产亚洲av高清一级| 在线观看免费高清a一片| 赤兔流量卡办理| 亚洲伊人久久精品综合| 免费看av在线观看网站| av在线观看视频网站免费| 麻豆精品久久久久久蜜桃| 日韩一本色道免费dvd| 老司机亚洲免费影院| 国产一区二区激情短视频 | 欧美精品人与动牲交sv欧美| 天美传媒精品一区二区| 天天影视国产精品| 男人操女人黄网站| 国产黄色视频一区二区在线观看| 国产av码专区亚洲av| 亚洲成人手机| 中文字幕人妻丝袜一区二区 | 国产日韩欧美亚洲二区| 欧美日韩一级在线毛片| 国产在视频线精品| 五月开心婷婷网| 看免费av毛片| 久久久欧美国产精品| 欧美日本中文国产一区发布| 不卡视频在线观看欧美| 国产一区二区三区av在线| 丰满少妇做爰视频| 午夜福利,免费看| 国产精品一国产av| a级毛片黄视频| 欧美日韩一区二区视频在线观看视频在线| 国产成人a∨麻豆精品| 国产精品二区激情视频| 日韩av免费高清视频| 久久久国产欧美日韩av| 色吧在线观看| 亚洲国产毛片av蜜桃av| 亚洲成国产人片在线观看| 丝瓜视频免费看黄片| 亚洲熟女精品中文字幕| 天堂俺去俺来也www色官网| 亚洲av欧美aⅴ国产| 国产黄色视频一区二区在线观看| 视频在线观看一区二区三区| 水蜜桃什么品种好| 欧美精品人与动牲交sv欧美| 日日爽夜夜爽网站| 国产在线视频一区二区| 午夜免费鲁丝| 国产成人一区二区在线| av天堂久久9| 亚洲成人国产一区在线观看 | av国产久精品久网站免费入址| 日韩一区二区视频免费看| 各种免费的搞黄视频| 国产深夜福利视频在线观看| 91aial.com中文字幕在线观看| 韩国av在线不卡| 母亲3免费完整高清在线观看| 亚洲欧美中文字幕日韩二区| 久久精品人人爽人人爽视色| 久久青草综合色| 免费在线观看完整版高清| 久久热在线av| 国产成人午夜福利电影在线观看| 国产精品一区二区在线不卡| 国产精品欧美亚洲77777| 在现免费观看毛片| 纵有疾风起免费观看全集完整版| 不卡av一区二区三区| 亚洲成人av在线免费| 午夜福利影视在线免费观看| 最近最新中文字幕大全免费视频 | 国产99久久九九免费精品| 十八禁网站网址无遮挡| 日日爽夜夜爽网站| 国产熟女欧美一区二区| 欧美精品亚洲一区二区| 中文欧美无线码| 亚洲专区中文字幕在线 | 亚洲专区中文字幕在线 | 国产一区有黄有色的免费视频| 多毛熟女@视频| 国产一卡二卡三卡精品 | 久久久久国产精品人妻一区二区| 亚洲一码二码三码区别大吗| www日本在线高清视频| 欧美日韩一区二区视频在线观看视频在线| 高清不卡的av网站| 欧美97在线视频| 中国三级夫妇交换| 97精品久久久久久久久久精品| 久久鲁丝午夜福利片| 中文字幕最新亚洲高清| 国产日韩欧美在线精品| 色综合欧美亚洲国产小说| 久久久久久久国产电影| 国产乱来视频区| 男的添女的下面高潮视频| 一区二区三区精品91| 国产乱来视频区| 又大又爽又粗| 精品国产超薄肉色丝袜足j| 日韩大码丰满熟妇| 建设人人有责人人尽责人人享有的| 永久免费av网站大全| 亚洲欧美精品综合一区二区三区| 男女边摸边吃奶| 十分钟在线观看高清视频www| 成年av动漫网址| 久久久久精品久久久久真实原创| 久久久久久人妻| 人人妻人人澡人人看| 免费黄网站久久成人精品| 国产熟女午夜一区二区三区| 侵犯人妻中文字幕一二三四区| 永久免费av网站大全| 国产精品久久久久久精品电影小说| 人成视频在线观看免费观看| 亚洲欧美一区二区三区国产| 亚洲国产欧美一区二区综合| 精品少妇一区二区三区视频日本电影 | 一边摸一边做爽爽视频免费| 国产精品.久久久| 欧美少妇被猛烈插入视频| tube8黄色片| 国产成人啪精品午夜网站| 欧美日韩亚洲高清精品| 视频区图区小说| 国产一区有黄有色的免费视频| 成人国产av品久久久| 啦啦啦在线观看免费高清www| 一个人免费看片子| 性少妇av在线| 无遮挡黄片免费观看| 亚洲国产av新网站| 久久人人97超碰香蕉20202| 亚洲av成人精品一二三区| 韩国精品一区二区三区| 日本91视频免费播放| 久久精品人人爽人人爽视色| 亚洲色图综合在线观看| 女人久久www免费人成看片| 中文字幕最新亚洲高清| 久久天堂一区二区三区四区| av一本久久久久| 热99久久久久精品小说推荐| 久久久久久久久免费视频了| 十八禁网站网址无遮挡| 久久国产亚洲av麻豆专区| 在线天堂中文资源库| 久久 成人 亚洲| 大香蕉久久网| 交换朋友夫妻互换小说| bbb黄色大片|