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

    一種基于神經(jīng)網(wǎng)絡(luò)的中制導(dǎo)改進算法

    2016-12-23 01:53:07魏倩蔡遠(yuǎn)利
    西安交通大學(xué)學(xué)報 2016年7期
    關(guān)鍵詞:訓(xùn)練樣本算例制導(dǎo)

    魏倩,蔡遠(yuǎn)利

    (西安交通大學(xué)電子與信息工程學(xué)院,710049,西安)

    ?

    一種基于神經(jīng)網(wǎng)絡(luò)的中制導(dǎo)改進算法

    魏倩,蔡遠(yuǎn)利

    (西安交通大學(xué)電子與信息工程學(xué)院,710049,西安)

    針對地球扁率影響下的大氣層外導(dǎo)彈中段制導(dǎo)問題,提出了基于BP神經(jīng)網(wǎng)絡(luò)的模型預(yù)測改進算法,并且創(chuàng)新使用軌道偏差解析解來構(gòu)造訓(xùn)練樣本集。首先,利用極點變換方法把彈體受到的J2項攝動引力優(yōu)化分解為與運動軌跡相關(guān)的擾動函數(shù);然后,采用偏差狀態(tài)空間的轉(zhuǎn)移矩陣,建立起導(dǎo)彈在J2項攝動作用下的軌道偏差公式;最后,利用偏差公式構(gòu)造取值廣泛的訓(xùn)練樣本集并訓(xùn)練BP神經(jīng)網(wǎng)絡(luò),從而建立起關(guān)于虛擬目標(biāo)信息的預(yù)測模型,計算出中段制導(dǎo)控制所需的增益速度矢量。該模型的優(yōu)點是利用極點變換和狀態(tài)轉(zhuǎn)移矩陣直接求解J2項攝動偏差,避免了進行大規(guī)模的數(shù)值積分運算;神經(jīng)網(wǎng)絡(luò)擁有強大的學(xué)習(xí)能力,保證了預(yù)測模型的全面性及精確性;BP神經(jīng)網(wǎng)絡(luò)可以預(yù)先離線訓(xùn)練、學(xué)習(xí),大大縮短了計算時間。與傳統(tǒng)Lambert迭代補償修正方法相比,改進型BP神經(jīng)網(wǎng)絡(luò)補償算法可以同時滿足實時計算速度及計算精度的雙重要求,具有較強的實際工程意義。

    極點變換;攝動偏差;J2項攝動;BP神經(jīng)網(wǎng)絡(luò)

    在大氣層外的遠(yuǎn)距離導(dǎo)彈攻防作戰(zhàn)中,彈體的飛行過程分為主動助推段、中段制導(dǎo)段和末制導(dǎo)段。中段制導(dǎo)段是指飛行器在主動段到達(dá)關(guān)機點后,末段制導(dǎo)開始之前,飛行器進行較長時間的無控滑行過程。在整個中段制導(dǎo)段,飛行器僅采用脈沖推力等方式進行少數(shù)幾次飛行狀態(tài)修正。因此,整個中段制導(dǎo)階段可以抽象為固定時間的Lambert制導(dǎo)問題,經(jīng)典方法求解法有傳統(tǒng)高斯解法、普適變量法、Battin-Vaughan算法等[1-5],若考慮地球扁率所產(chǎn)生的J2項引力攝動(引力場函數(shù)僅包含2階帶諧項),這些方法會產(chǎn)生較大的軌道偏差。

    現(xiàn)階段對于J2攝動下的Lambert制導(dǎo)問題的研究,主要為引力差修正法和補償Lambert制導(dǎo)法[6-11]。引力差修正法是在中段制導(dǎo)過程中,分階段對于引力差所產(chǎn)生的軌道偏差進行彈道修正,其優(yōu)點是制導(dǎo)精度對于飛行狀態(tài)參數(shù)的魯棒性較好,缺點是引力差的計算進行了大量的近似處理,且修正過程持續(xù)整個中制導(dǎo)階段。補償Lambert制導(dǎo)法是在經(jīng)典Lambert制導(dǎo)基礎(chǔ)上,補償修正J2引力攝動對飛行彈道的影響,是一種簡單、實用的制導(dǎo)方法。但是,現(xiàn)有的補償制導(dǎo)法在計算精度和運行時間兩方面難以取得平衡,限制了上述方法在實際工程中的應(yīng)用。

    本文提出了基于BP神經(jīng)網(wǎng)絡(luò)的模型預(yù)測制導(dǎo)算法,并使用偏差解析解來構(gòu)造取值廣泛的訓(xùn)練樣本集,這可大大縮短計算時間,滿足實時在線計算的要求。

    1 改進Lambert制導(dǎo)問題描述

    在研究Lambert制導(dǎo)時,通常假設(shè)飛行器可以瞬間獲得轉(zhuǎn)移軌道所需要的速度增益矢量,不涉及脈沖推力控制問題。

    圖1 經(jīng)典Lambert制導(dǎo)下目標(biāo)點P2與虛擬目標(biāo)點

    2 軌道偏差解析解

    2.1 J2項引力攝動的分解

    彈體受到的J2項引力函數(shù)可表示為[1]

    (1)

    式中:r為地心距;φ為地心緯度;ae=6 378 137;J=1.082 6×10-3。

    圖2 極點變換的幾何關(guān)系

    關(guān)于極點Op(λp,φp)可由初始點P1位置信息得到,即

    (2)

    因此,標(biāo)準(zhǔn)橢圓軌道上任意點Q(λ,φ),變換出的新坐標(biāo)(η,σ)關(guān)系式為

    (3)

    J2項引力攝動的加速度,關(guān)于導(dǎo)彈軌跡運動的分量表示為

    (4)

    其中

    (5)

    2.2 軌道偏差解析解

    為避免直接求解微分方程帶來的困難,本文采用偏差狀態(tài)的轉(zhuǎn)移矩陣直接解算出軌道偏差的解析解。

    (6)

    當(dāng)考慮到實際軌道與標(biāo)準(zhǔn)橢圓軌道同為地心距r時,軌道狀體參數(shù)的偏差稱為等地心距偏差。設(shè)等地心距偏差構(gòu)成的狀態(tài)向量為

    (7)

    (8)

    其中

    (9)

    圖3 等地心距偏差與等地心角偏差的關(guān)系

    因此,Yr的攝動偏差方程表示為

    (10)

    由于J2項引力擾動,導(dǎo)彈的軌道偏差可分解為東向和北向兩個方向

    (11)

    式中:Ap2為目標(biāo)點P2的方位角。

    (12)

    3 改進型Lambert制導(dǎo)算法

    3.1 虛擬目標(biāo)點預(yù)測模型

    考慮到Lambert制導(dǎo)問題對于預(yù)測精度的要求,引入多隱含層來提高網(wǎng)絡(luò)預(yù)測精度。通過對比分析,采用網(wǎng)絡(luò)結(jié)構(gòu)為3-9-9-3的4層BP神經(jīng)網(wǎng)絡(luò)滿足模型的預(yù)測精度。

    使用BP神經(jīng)網(wǎng)絡(luò)建模時,訓(xùn)練樣本的選取十分重要。訓(xùn)練樣本滿足2個條件:①樣本數(shù)量要足夠多,能夠反映出關(guān)系的復(fù)雜程度;②樣本需要具有代表性,盡可能包含各種特征模式。

    在大氣層外導(dǎo)彈的取值范圍內(nèi)以任務(wù)目標(biāo)點P2為中心選取采樣數(shù)據(jù),構(gòu)成采樣矩陣Ti,j,k(λTi,φTj,tTk),滿足條件

    (13)

    式中:Δλ、Δφ、Δt為小量的采樣間隔;i=0,±1,±2,…,±nλ,j=0,±1,±2,…,±nφ,k=0,±1,±2,…,±nt,nλ、nφ和nt分別決定在λ、φ和tF上采樣規(guī)模的大小。

    圖4 神經(jīng)網(wǎng)絡(luò)訓(xùn)練原理圖

    因為Δλ、Δφ取值較小,采樣樣本集中分布在目標(biāo)點周圍一定的范圍內(nèi),此樣本集減小了J2項引力攝動對位置、飛行時間等的復(fù)雜影響,使得樣本的輸出信息也在小范圍內(nèi)變化,小范圍變化的各種因素使得映射關(guān)系接近于線性關(guān)系。同時,此采樣集兼顧考慮了時間、空間等因素不同數(shù)量級變化對于飛行軌道偏差的實際影響。

    3.2 制導(dǎo)方法流程

    4 算法驗證與分析

    假設(shè)遠(yuǎn)距離飛行的導(dǎo)彈只受到J2項引力攝動作用的影響,且不進行導(dǎo)彈末段制導(dǎo)修正,通過仿真實驗對以上改進型Lambert制導(dǎo)算法的正確性、可靠性和精度進行驗證。

    4.1 仿真算例分析

    取導(dǎo)彈關(guān)機點P1在地心東經(jīng)緯度為0°、北緯經(jīng)度為10°、高度為1 200 km處,用2個具體算例來對比分析。在算例1中,假設(shè)導(dǎo)彈的任務(wù)目標(biāo)信息為(46.2°,7°,1 150 s)。在考慮到J2項引力攝動下,導(dǎo)彈的無修正落點為(46.260 4°,7.006 5°,1 148 s),落點誤差6 754 m。通過BP神經(jīng)網(wǎng)絡(luò)預(yù)測出虛擬目標(biāo)點(46.139 7°,6.993 6°,1 151.6 s),BP改進算法的修正效果良好,修正落點誤差僅為7.6 m。在算例2中,假設(shè)導(dǎo)彈的任務(wù)目標(biāo)為(35°,17°,800 s)。在考慮到J2項引力攝動下,導(dǎo)彈的無修正落點為(34.959°,16.985°,796.6 s),落點誤差為4 844 m。通過BP神經(jīng)網(wǎng)絡(luò)預(yù)測出虛擬目標(biāo)點(35.040 8°,17.015 3°,803.3 s),修正落點誤差僅為2.972 m。

    通過對2個算例數(shù)據(jù)對比可看出:在J2項引力攝動的影響下,導(dǎo)彈的位置偏差為104m,并且偏差隨著飛行時間的延長、飛行距離的增長而不斷變大;改進型算法利用BP神經(jīng)網(wǎng)絡(luò)預(yù)測出虛擬目標(biāo)點,可很好修正J2項引力攝動引起的軌道位置偏差。

    根據(jù)BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練樣本的選取規(guī)則,算例1中,將任務(wù)目標(biāo)點的信息(46.2°,7°,1 150 s)作為BP神經(jīng)網(wǎng)絡(luò)的訓(xùn)練樣本的采樣中心點,nλ=10、nφ=10、nt=5、Δλ=0.57°、Δφ=0.57°、Δt=20 s,采樣樣本個數(shù)為4 851,隨機抽取4 751個樣本作為訓(xùn)練樣本,剩余的100個樣本作為測試數(shù)據(jù)。

    圖5 算例1的測試誤差率

    利用算例1中的訓(xùn)練數(shù)據(jù)對預(yù)先訓(xùn)練完成的BP神經(jīng)網(wǎng)絡(luò)(BPnet)進行第2次更新訓(xùn)練,迭代次數(shù)為100,均方誤差為10-6。由于BPnet網(wǎng)絡(luò)經(jīng)過預(yù)先訓(xùn)練過程,二次訓(xùn)練的迭代次數(shù)為3,訓(xùn)練誤差為9.37×10-7,運行時間為1.3 s,測試樣本數(shù)據(jù)的誤差為9.6×10-6,測試誤差率如圖5所示。對于算例2進行相同的神經(jīng)網(wǎng)絡(luò)訓(xùn)練,得到二次訓(xùn)練的迭代次數(shù)為1,訓(xùn)練誤差為8.6×10-7,運行時間為1.64s,測試樣本誤差為8.5×10-6,測試數(shù)據(jù)的測試誤差率如圖6所示。對于BPnet網(wǎng)絡(luò)的二次更新訓(xùn)練,使得靜態(tài)BP神經(jīng)網(wǎng)絡(luò)能夠滿足不同任務(wù)目標(biāo)點的相關(guān)要求,并且在同一任務(wù)目標(biāo)點周圍一定區(qū)域中,BP神經(jīng)網(wǎng)絡(luò)具有較好的預(yù)測能力。

    通過2個算例的仿真計算可看出,本文討論的基于BP神經(jīng)網(wǎng)絡(luò)的改進型Lambert制導(dǎo)算法,可很好修正J2項攝動產(chǎn)生的軌道偏差,實現(xiàn)精確的飛行制導(dǎo)。

    圖6 算例2的測試誤差率

    4.2 算法對比與分析

    忽略關(guān)機點測量誤差,且不進行末制導(dǎo)修正,在相同的運算環(huán)境下分別用BP神經(jīng)網(wǎng)絡(luò)改進法及經(jīng)典虛擬目標(biāo)補償修正法進行對比仿真分析。不同算法的仿真精度、計算結(jié)果以及所耗費的運行時間如表1所示。

    表1 不同算法的仿真效果對比

    經(jīng)典的虛擬目標(biāo)點補償法[1]采用數(shù)值積分,整個飛行仿真周期為50 s,計算出位置偏差并進行迭代修正,次數(shù)一般不小于3,數(shù)值積分所耗費的時間相當(dāng)可觀,并且軌道積分的路徑差別也影響計算精度。改進型算法使用BP神經(jīng)網(wǎng)絡(luò)建立目標(biāo)點與實際位置之間的映射關(guān)系,從而可精準(zhǔn)預(yù)測出虛擬目標(biāo)點信息。

    仿真對比實驗結(jié)果表明,BP神經(jīng)網(wǎng)絡(luò)改進型Lambert制導(dǎo)策略在計算精度相當(dāng)?shù)那闆r下,具有運算量小、方法簡便、易于工程實現(xiàn)等優(yōu)點。

    5 結(jié) 論

    考慮J2項攝動影響的Lambert問題是一個典型的多變量非線性問題,常規(guī)方法不能滿足計算時間和運算精度的雙重要求。本文利用偏差解析解構(gòu)造出基于BP神經(jīng)網(wǎng)絡(luò)的預(yù)測模型,并計算得到補償制導(dǎo)的控制參數(shù),其最主要的優(yōu)點是:利用狀態(tài)轉(zhuǎn)移矩陣直接求解J2項攝動引起的軌道偏差,避免了進行大規(guī)模的數(shù)值積分運算;神經(jīng)網(wǎng)絡(luò)的強大學(xué)習(xí)能力,保證了預(yù)測模型的全面性以及精確性;BP神經(jīng)網(wǎng)絡(luò)的離線訓(xùn)練、學(xué)習(xí),大大縮短了計算時間。此改進型制導(dǎo)算法可同時滿足實時計算效率及計算精度的雙重要求,具有較大的實際工程意義。

    [1] BATTIN R H. Introduction to the mathematics and methods of astrodynamics [M]. Reston, VA, USA: AIAA, 1999: 325-417.

    [2] AHN J, lEE S. Lambert algorithm using analytic gradients [J]. Journal of Guidance, Control, and Dynamics, 2013, 36(6): 1751-1761.

    [3] AHN J, BANG J, LEE S. Acceleration of zero-revolution Lambert’s algorithms using table-based initialization [J]. Journal of Guidance, Control, and Dynamics, 2015, 38(2): 335-342.

    [4] ZHANG Gang, CAO Xibin, ZHOU Di. Two-impulse cotangent rendezvous between coplanar elliptic and hyperbolic orbits [J]. Journal of Guidance, Control, and Dynamics, 2014, 37(3): 965-969.

    [5] WAILLIEZ S E. On Lambert’s problem and the elliptic time of flight equation: a simple semi-analytical inversion method [J]. Advances in Space Research, 2014, 53(5): 890-898.

    [6] 高海燕, 蔡遠(yuǎn)利. 高超聲速飛行器的滑模預(yù)測控制方法 [J]. 西安交通大學(xué)學(xué)報, 2014, 48(1): 67-72. GAO Haiyan, CAI Yuanli. Sliding mode predictive control for hypersonic vehicle [J]. Journal of Xi’an Jiaotong University, 2014, 48(1): 67-72.

    [7] 徐明, 譚田, 李志武, 等. Lambert轉(zhuǎn)移中途修正的全局概率最優(yōu)策略 [J]. 北京航空航天大學(xué)學(xué)報, 2012, 38(5): 574-578. XU Ming, TAN Tian, LI Zhiwu, et al. Optimal correction strategy during Lambert transfer from view of probability [J]. Journal of Beijing University of Aeronautics and Astronautics, 2012, 38(5): 574-578.

    [8] 胡正東, 郭才發(fā), 曹淵, 等. 軌道轟炸飛行器過渡段軌道設(shè)計與制導(dǎo) [J]. 固體火箭技術(shù), 2009, 32(5): 473-479. HU Zhengdong, GUO Caifa, CAO Yuan, et al. Transition trajectory planning and guidance for orbital bombing vehicle [J]. Journal of Solid Rocket Technology, 2009, 32(5): 473-479.

    [9] 柴華, 仲明, 梁彥剛. 采用狀態(tài)轉(zhuǎn)移矩陣的攔截中段制導(dǎo)方法 [J]. 國防科技大學(xué)學(xué)報, 2015, 37(4): 137-142. CHAI Hua, ZHONG Ming, LIANG Yangang. Midcourse guidance of interception using state transition matrix [J]. Journal of National University of Defense Technology, 2015, 37(4): 137-142.

    [10]鄭偉, 湯國建. 擾動引力場中彈道導(dǎo)彈飛行力學(xué) [M]. 北京: 國防工業(yè)出版社, 2009: 32-45.

    [11]YAMADA K, KIMURA M. New state transition matrix for formation flying inJ2-perturbed elliptic orbits [J]. Journal of Guidance, Control, and Dynamics, 2012, 35(2): 536-547.

    [12]鄧一兵, 胡偉, 高峰, 等. 遺傳神經(jīng)網(wǎng)絡(luò)在載人飛船環(huán)控決策系統(tǒng)中的應(yīng)用研究 [J]. 西安交通大學(xué)學(xué)報, 2010, 44(7): 64-69. DENG Yibing, HU Wei, GAO Feng, et al. Application of genetic neural network to decision support system for environmental control and life support system [J]. Journal of Xi’an Jiaotong University, 2010, 44(7): 64-69.

    [本刊相關(guān)文獻(xiàn)鏈接]

    張小棟,郭晉,李睿,等.表情驅(qū)動下腦電信號的建模仿真及分類識別.2016,50(6):1-8.[doi:10.7652/xjtuxb201606001]

    姜濤,黃偉,王安麟.多路閥閥芯節(jié)流槽拓?fù)浣Y(jié)構(gòu)組合的神經(jīng)網(wǎng)絡(luò)模型.2016,50(6):36-41.[doi:10.7652/xjtuxb201606 006]

    宋青松,田正鑫,孫文磊,等.用于孤立數(shù)字語音識別的一種組合降維方法.2016,50(6):42-46.[doi:10.7652/xjtuxb2016 06007]

    陳江城,張小棟.人體下肢行走關(guān)節(jié)連續(xù)運動表面肌電解碼方法.2016,50(6):61-67.[doi:10.7652/xjtuxb201606010]

    陳斌,胡平舸,屈丹.子空間域相關(guān)特征變換與融合的語音識別方法.2016,50(4):60-67.[doi:10.7652/xjtuxb201604010]

    王安麟,孟慶華,曹巖,等.液力變矩器的葉片數(shù)神經(jīng)網(wǎng)絡(luò)模型.2015,49(7):11-16.[doi:10.7652/xjtuxb201507003]

    馬忠麗,文杰,梁秀梅,等.無人艇視覺系統(tǒng)多類水面目標(biāo)特征提取與識別.2014,48(8):60-66.[doi:10.7652/xjtuxb2014 08011]

    侯雨伸,王秀麗.氣象過程信息挖掘與輸電線路覆冰預(yù)測.2014,48(6):43-49.[doi:10.7652/xjtuxb201406008]

    王萬召,王杰.采用限定記憶極限學(xué)習(xí)機的過熱汽溫逆建模研究.2014,48(2):32-37.[doi:10.7652/xjtuxb201402006]

    馮斌,梅雪松,楊軍,等.數(shù)控機床摩擦誤差自適應(yīng)補償方法研究.2013,47(11):65-69.[doi:10.7652/xjtuxb201311012]

    (編輯 趙煒)

    A Modified Algorithm on the Midcourse Guidance Based on BP Neural Network

    WEI Qian,CAI Yuanli

    (School of Electronic and Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China)

    Aiming at solving the midcourse guidance problem of exo-atmosphere missile under the influence of earth’s oblateness perturbation, a prediction model for missile guidance was proposed based on BP neural network. This method provides a new training sample set constructed by the analytical formula of trajectory deviation. First, by using the pole transform method, the missile’sJ2perturbed gravity is decomposed into the disturbing function related to its flight trajectory. Then, with the state space matrix method, the analytic solution of trajectory deviation with theJ2perturbation is calculated. Finally, using trajectory deviation function to construct a wide range of training sample set, the BP neural network of prediction model is established. The neural network can forecast the virtual target point information, so as to calculate the vector of gained velocity for midcourse guidance control. Using the modified algorithm, the trajectory deviation ofJ2perturbation can be directly solved by pole transform and state transition matrix, avoiding large-scale numerical calculation. The BP neural network has powerful learning and training ability, ensuring the comprehensiveness and accuracy of the prediction model and saving calculation time by the off-line training and learning before the simulation tests. In comparison with traditional correction method of Lambert guidance, this modified algorithm can satisfy the requirements on both efficiency and accuracy of real-time computation, being of practical engineering significance.

    pole transform; perturbation deviation;J2perturbations; BP neural network

    2015-12-04。 作者簡介:魏倩(1984—),女,博士生;蔡遠(yuǎn)利(通信作者),男,教授,博士生導(dǎo)師。 基金項目:國家自然科學(xué)基金資助項目(61308120,61463029)。

    10.7652/xjtuxb201607019

    V448.2

    A

    0253-987X(2016)07-0125-06

    猜你喜歡
    訓(xùn)練樣本算例制導(dǎo)
    人工智能
    寬帶光譜成像系統(tǒng)最優(yōu)訓(xùn)練樣本選擇方法研究
    基于MPSC和CPN制導(dǎo)方法的協(xié)同制導(dǎo)律
    基于在線軌跡迭代的自適應(yīng)再入制導(dǎo)
    融合原始樣本和虛擬樣本的人臉識別算法
    基于稀疏重構(gòu)的機載雷達(dá)訓(xùn)練樣本挑選方法
    帶有攻擊角約束的無抖振滑模制導(dǎo)律設(shè)計
    基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
    互補問題算例分析
    基于CYMDIST的配電網(wǎng)運行優(yōu)化技術(shù)及算例分析
    精品久久国产蜜桃| 99久久精品一区二区三区| 99久久人妻综合| 天堂√8在线中文| 免费看光身美女| 秋霞伦理黄片| 国产大屁股一区二区在线视频| 天堂俺去俺来也www色官网 | 国产免费视频播放在线视频 | 中文字幕av在线有码专区| 真实男女啪啪啪动态图| 欧美潮喷喷水| 国产免费又黄又爽又色| 黄色欧美视频在线观看| 啦啦啦韩国在线观看视频| 日韩欧美一区视频在线观看 | 日韩三级伦理在线观看| 国产亚洲91精品色在线| 啦啦啦中文免费视频观看日本| 美女cb高潮喷水在线观看| 又爽又黄无遮挡网站| 黄色欧美视频在线观看| 欧美+日韩+精品| 亚洲真实伦在线观看| 国产亚洲最大av| 中文乱码字字幕精品一区二区三区 | 日韩国内少妇激情av| 汤姆久久久久久久影院中文字幕 | 亚洲天堂国产精品一区在线| 国产精品久久久久久久电影| 七月丁香在线播放| 男插女下体视频免费在线播放| 亚洲精品一区蜜桃| 国产老妇女一区| 婷婷六月久久综合丁香| 欧美成人一区二区免费高清观看| 十八禁国产超污无遮挡网站| 成人毛片60女人毛片免费| 亚洲自偷自拍三级| 久99久视频精品免费| 啦啦啦韩国在线观看视频| 两个人的视频大全免费| 精品一区在线观看国产| 女人被狂操c到高潮| 中国美白少妇内射xxxbb| 婷婷色麻豆天堂久久| 一级毛片电影观看| 啦啦啦啦在线视频资源| 日韩,欧美,国产一区二区三区| 国产精品无大码| 18禁在线播放成人免费| 欧美高清性xxxxhd video| 亚洲色图av天堂| 免费播放大片免费观看视频在线观看| 有码 亚洲区| 免费观看性生交大片5| 国产精品久久久久久精品电影小说 | 亚洲最大成人av| 非洲黑人性xxxx精品又粗又长| 亚洲三级黄色毛片| 亚洲精品日韩在线中文字幕| 日本av手机在线免费观看| 一级av片app| 特级一级黄色大片| 九九久久精品国产亚洲av麻豆| 欧美成人a在线观看| 蜜桃久久精品国产亚洲av| 免费大片黄手机在线观看| 日韩中字成人| 亚洲va在线va天堂va国产| 中文在线观看免费www的网站| 色网站视频免费| 99久久精品热视频| 精品久久久精品久久久| 成年女人在线观看亚洲视频 | 精品一区在线观看国产| 久久草成人影院| 日本熟妇午夜| 国产成人91sexporn| 能在线免费看毛片的网站| 免费黄频网站在线观看国产| 日韩av不卡免费在线播放| 亚洲精品国产成人久久av| 非洲黑人性xxxx精品又粗又长| 欧美日本视频| 男女下面进入的视频免费午夜| 成人午夜高清在线视频| 老师上课跳d突然被开到最大视频| av.在线天堂| 亚洲国产欧美在线一区| 亚洲精品成人av观看孕妇| 中文资源天堂在线| 久久久久精品久久久久真实原创| 日本一二三区视频观看| 人人妻人人看人人澡| 自拍偷自拍亚洲精品老妇| 视频中文字幕在线观看| 熟妇人妻不卡中文字幕| 亚洲av日韩在线播放| 国产亚洲午夜精品一区二区久久 | 亚洲av中文字字幕乱码综合| 2021天堂中文幕一二区在线观| 亚洲精品久久午夜乱码| 老师上课跳d突然被开到最大视频| 国产一级毛片七仙女欲春2| 精品久久久久久久久av| 欧美3d第一页| 人体艺术视频欧美日本| 淫秽高清视频在线观看| 欧美极品一区二区三区四区| 最近中文字幕2019免费版| 日韩三级伦理在线观看| 超碰97精品在线观看| 黄色一级大片看看| 在线观看三级黄色| 久久精品夜色国产| 女性生殖器流出的白浆| 亚洲综合色惰| 日韩欧美精品免费久久| 亚洲精品在线美女| 18禁观看日本| 精品视频人人做人人爽| 黄频高清免费视频| 久久久精品国产亚洲av高清涩受| 亚洲成国产人片在线观看| 精品国产国语对白av| 久久精品国产鲁丝片午夜精品| 99精国产麻豆久久婷婷| 亚洲国产av新网站| 久久精品国产综合久久久| 九九爱精品视频在线观看| 女性生殖器流出的白浆| 免费播放大片免费观看视频在线观看| 天美传媒精品一区二区| 黄色 视频免费看| 亚洲精品美女久久av网站| 如日韩欧美国产精品一区二区三区| 国产一区二区 视频在线| 国产精品99久久99久久久不卡 | 日本欧美视频一区| 久久国产精品大桥未久av| 日本wwww免费看| 亚洲 欧美一区二区三区| 国产日韩欧美视频二区| 国产成人欧美| 欧美av亚洲av综合av国产av | 亚洲人成77777在线视频| 麻豆av在线久日| 久久久久久伊人网av| 国产亚洲午夜精品一区二区久久| 一本久久精品| 妹子高潮喷水视频| 青春草视频在线免费观看| 亚洲经典国产精华液单| 黄片无遮挡物在线观看| 久久久精品94久久精品| 一本久久精品| 男人添女人高潮全过程视频| 精品国产乱码久久久久久小说| 久久国内精品自在自线图片| 国产精品国产av在线观看| 成人黄色视频免费在线看| 亚洲国产最新在线播放| 尾随美女入室| 桃花免费在线播放| 99热全是精品| 久久人妻熟女aⅴ| av一本久久久久| 丝袜在线中文字幕| 久久久久久久久久久久大奶| 亚洲欧美一区二区三区久久| a级片在线免费高清观看视频| 亚洲美女视频黄频| 少妇熟女欧美另类| 90打野战视频偷拍视频| 啦啦啦在线观看免费高清www| 人妻系列 视频| 成人免费观看视频高清| 欧美国产精品一级二级三级| 伊人久久国产一区二区| 在线观看免费日韩欧美大片| www日本在线高清视频| 日韩在线高清观看一区二区三区| 成人国语在线视频| 中文字幕亚洲精品专区| 男的添女的下面高潮视频| 久久青草综合色| 国产日韩欧美在线精品| 精品国产乱码久久久久久小说| 老司机影院成人| 中文字幕人妻丝袜一区二区 | 国产精品久久久久成人av| 人妻 亚洲 视频| 免费播放大片免费观看视频在线观看| 国产乱来视频区| 国产av精品麻豆| 中文字幕最新亚洲高清| 欧美在线黄色| 免费播放大片免费观看视频在线观看| 黄网站色视频无遮挡免费观看| 黑丝袜美女国产一区| 老司机亚洲免费影院| 久久久久精品人妻al黑| 国产视频首页在线观看| 极品少妇高潮喷水抽搐| 美女主播在线视频| 麻豆乱淫一区二区| 成年美女黄网站色视频大全免费| 国产1区2区3区精品| 国产色婷婷99| 国产麻豆69| 国产 精品1| 国产精品.久久久| 2018国产大陆天天弄谢| 中文字幕制服av| 交换朋友夫妻互换小说| 美国免费a级毛片| 1024香蕉在线观看| 大码成人一级视频| 乱人伦中国视频| 午夜日本视频在线| av不卡在线播放| 欧美日本中文国产一区发布| 激情五月婷婷亚洲| 免费少妇av软件| 欧美日韩精品成人综合77777| 亚洲国产色片| 人体艺术视频欧美日本| 国产免费视频播放在线视频| 深夜精品福利| 国产视频首页在线观看| 成年女人在线观看亚洲视频| 最近最新中文字幕免费大全7| 激情视频va一区二区三区| 欧美激情高清一区二区三区 | 狠狠精品人妻久久久久久综合| 高清黄色对白视频在线免费看| av一本久久久久| av电影中文网址| 欧美国产精品一级二级三级| 国产 精品1| 1024香蕉在线观看| 久久精品国产自在天天线| 午夜福利视频精品| 久久精品国产鲁丝片午夜精品| 亚洲视频免费观看视频| 久久狼人影院| 久久久久久久久久人人人人人人| 国产精品免费大片| 热re99久久精品国产66热6| 亚洲人成77777在线视频| 欧美国产精品一级二级三级| 丝袜在线中文字幕| 精品国产一区二区三区久久久樱花| 777米奇影视久久| 久久久精品区二区三区| 国产一区二区三区综合在线观看| 色吧在线观看| 亚洲欧美精品综合一区二区三区 | 不卡视频在线观看欧美| 国产极品粉嫩免费观看在线| 天天躁日日躁夜夜躁夜夜| 日本欧美视频一区| 1024视频免费在线观看| 校园人妻丝袜中文字幕| 99久久人妻综合| 满18在线观看网站| 多毛熟女@视频| 国产日韩欧美视频二区| 91精品三级在线观看| 观看av在线不卡| 少妇猛男粗大的猛烈进出视频| 天天躁夜夜躁狠狠躁躁| 欧美黄色片欧美黄色片| 制服诱惑二区| 国产精品三级大全| 国产黄色视频一区二区在线观看| 午夜免费观看性视频| 国产高清不卡午夜福利| 看非洲黑人一级黄片| 久久ye,这里只有精品| 激情视频va一区二区三区| 国产一区二区三区综合在线观看| 日韩精品有码人妻一区| 国产成人免费无遮挡视频| 香蕉精品网在线| 九草在线视频观看| 久久久久国产一级毛片高清牌| 九九爱精品视频在线观看| 国产精品成人在线| 一二三四在线观看免费中文在| 黄频高清免费视频| 精品卡一卡二卡四卡免费| 人人妻人人爽人人添夜夜欢视频| 国产人伦9x9x在线观看 | 久久影院123| 日本欧美视频一区| xxx大片免费视频| 国产在线一区二区三区精| 观看av在线不卡| 建设人人有责人人尽责人人享有的| 亚洲欧美精品自产自拍| 久久国内精品自在自线图片| 免费观看av网站的网址| 美女xxoo啪啪120秒动态图| 久久人人爽人人片av| 久热这里只有精品99| 视频在线观看一区二区三区| 视频区图区小说| 美女福利国产在线| 久久这里有精品视频免费| a级毛片在线看网站| 国产精品一区二区在线不卡| 久久久精品区二区三区| 国产成人精品福利久久| 一本色道久久久久久精品综合| av天堂久久9| 欧美 日韩 精品 国产| 国产探花极品一区二区| 亚洲激情五月婷婷啪啪| 国产av国产精品国产| 蜜桃国产av成人99| 女性生殖器流出的白浆| 国产97色在线日韩免费| 最近手机中文字幕大全| 两个人免费观看高清视频| 日本午夜av视频| 春色校园在线视频观看| 最近2019中文字幕mv第一页| 亚洲第一青青草原| 国产一区二区三区av在线| 亚洲av男天堂| 少妇被粗大猛烈的视频| av女优亚洲男人天堂| 午夜福利,免费看| 亚洲精品视频女| 亚洲精品久久午夜乱码| 国产一区有黄有色的免费视频| 中文欧美无线码| 一区二区三区乱码不卡18| 国产精品麻豆人妻色哟哟久久| 十八禁网站网址无遮挡| 国产精品一国产av| 精品久久久精品久久久| 国产成人精品一,二区| 婷婷色麻豆天堂久久| 国产老妇伦熟女老妇高清| 欧美日韩一区二区视频在线观看视频在线| 一本久久精品| 午夜福利在线观看免费完整高清在| av电影中文网址| 亚洲成人av在线免费| 亚洲精品久久午夜乱码| 十八禁网站网址无遮挡| 国产成人免费无遮挡视频| 欧美激情极品国产一区二区三区| 男男h啪啪无遮挡| 欧美激情极品国产一区二区三区| 毛片一级片免费看久久久久| 午夜福利在线免费观看网站| 性色avwww在线观看| 亚洲av电影在线进入| 亚洲成人手机| 日本免费在线观看一区| 欧美国产精品va在线观看不卡| 制服丝袜香蕉在线| 曰老女人黄片| 汤姆久久久久久久影院中文字幕| 妹子高潮喷水视频| 日韩伦理黄色片| 纯流量卡能插随身wifi吗| 最黄视频免费看| 超碰97精品在线观看| 高清黄色对白视频在线免费看| 亚洲av男天堂| 一级毛片 在线播放| 欧美精品亚洲一区二区| 亚洲精品美女久久久久99蜜臀 | 久久久久久久大尺度免费视频| 欧美最新免费一区二区三区| 成人18禁高潮啪啪吃奶动态图| 精品亚洲乱码少妇综合久久| 制服诱惑二区| 一二三四在线观看免费中文在| 欧美日韩综合久久久久久| 国产av国产精品国产| 赤兔流量卡办理| 国产国语露脸激情在线看| 老司机影院成人| 99re6热这里在线精品视频| 各种免费的搞黄视频| 亚洲欧美日韩另类电影网站| 精品少妇黑人巨大在线播放| 国产 一区精品| 一本大道久久a久久精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久这里只有精品19| 中文字幕精品免费在线观看视频| 欧美变态另类bdsm刘玥| 亚洲欧洲国产日韩| 精品少妇一区二区三区视频日本电影 | 国产成人一区二区在线| 国产亚洲午夜精品一区二区久久| 国产爽快片一区二区三区| 大香蕉久久成人网| 欧美日韩精品网址| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久久av不卡| 免费在线观看完整版高清| 有码 亚洲区| 精品亚洲成a人片在线观看| 中文字幕人妻熟女乱码| 少妇 在线观看| 边亲边吃奶的免费视频| 日韩免费高清中文字幕av| 久久精品人人爽人人爽视色| 久久青草综合色| 久久av网站| 男女边吃奶边做爰视频| 久久久a久久爽久久v久久| 日韩精品有码人妻一区| 国产成人欧美| 日韩人妻精品一区2区三区| 国产精品久久久久成人av| 免费少妇av软件| 国产精品无大码| 久久精品国产综合久久久| 欧美+日韩+精品| 最近最新中文字幕免费大全7| av女优亚洲男人天堂| 美女脱内裤让男人舔精品视频| 免费av中文字幕在线| 天天躁夜夜躁狠狠躁躁| 成人免费观看视频高清| 国产亚洲欧美精品永久| 亚洲欧美色中文字幕在线| 日韩制服骚丝袜av| 婷婷色av中文字幕| 天堂中文最新版在线下载| av国产精品久久久久影院| 一区福利在线观看| 亚洲熟女精品中文字幕| 欧美 亚洲 国产 日韩一| 日韩欧美一区视频在线观看| 亚洲伊人久久精品综合| 国产亚洲一区二区精品| 香蕉丝袜av| 免费观看性生交大片5| 久久这里有精品视频免费| 亚洲av.av天堂| 香蕉精品网在线| 两性夫妻黄色片| 久久精品夜色国产| 十八禁高潮呻吟视频| 捣出白浆h1v1| 丰满乱子伦码专区| 波多野结衣av一区二区av| 观看av在线不卡| 性色av一级| 亚洲国产精品一区三区| 亚洲欧美色中文字幕在线| 两性夫妻黄色片| 免费在线观看视频国产中文字幕亚洲 | 午夜精品国产一区二区电影| 国产高清不卡午夜福利| 久久久精品国产亚洲av高清涩受| 天天躁夜夜躁狠狠躁躁| 自拍欧美九色日韩亚洲蝌蚪91| 久久久精品免费免费高清| 人人妻人人添人人爽欧美一区卜| 美女脱内裤让男人舔精品视频| 人妻少妇偷人精品九色| 国产熟女午夜一区二区三区| 国产野战对白在线观看| 午夜久久久在线观看| 91精品三级在线观看| 午夜91福利影院| 国产成人精品久久二区二区91 | 亚洲av在线观看美女高潮| 久久久久国产一级毛片高清牌| 亚洲精华国产精华液的使用体验| av卡一久久| 男女免费视频国产| 热re99久久精品国产66热6| 日韩制服骚丝袜av| 久久久精品免费免费高清| av在线播放精品| 黄频高清免费视频| 国产亚洲一区二区精品| 七月丁香在线播放| kizo精华| 午夜日韩欧美国产| 交换朋友夫妻互换小说| 飞空精品影院首页| 国产男人的电影天堂91| 国产探花极品一区二区| 亚洲av免费高清在线观看| 亚洲精品,欧美精品| 国产在线一区二区三区精| 色哟哟·www| 久久鲁丝午夜福利片| 亚洲四区av| 91国产中文字幕| 新久久久久国产一级毛片| 91久久精品国产一区二区三区| 久久鲁丝午夜福利片| 99久国产av精品国产电影| 国产精品偷伦视频观看了| 精品人妻一区二区三区麻豆| 伦理电影免费视频| 亚洲人成77777在线视频| 国产黄色视频一区二区在线观看| 成人国产麻豆网| 精品视频人人做人人爽| 日日爽夜夜爽网站| 男女下面插进去视频免费观看| 婷婷成人精品国产| 亚洲av福利一区| 十八禁网站网址无遮挡| 在线 av 中文字幕| 69精品国产乱码久久久| 国产成人精品一,二区| 国产综合精华液| 在线天堂最新版资源| 麻豆乱淫一区二区| 国产成人免费观看mmmm| 久久av网站| 夜夜骑夜夜射夜夜干| 大香蕉久久成人网| 亚洲国产毛片av蜜桃av| 黄片播放在线免费| 欧美xxⅹ黑人| 欧美 日韩 精品 国产| 自拍欧美九色日韩亚洲蝌蚪91| 老汉色av国产亚洲站长工具| 中文乱码字字幕精品一区二区三区| av免费观看日本| 在线观看一区二区三区激情| 欧美日韩视频高清一区二区三区二| 肉色欧美久久久久久久蜜桃| 热re99久久精品国产66热6| 亚洲成国产人片在线观看| 亚洲成av片中文字幕在线观看 | 成人亚洲精品一区在线观看| 免费黄色在线免费观看| 久久久久久伊人网av| 欧美激情 高清一区二区三区| 亚洲精品国产av成人精品| 两性夫妻黄色片| 国产又爽黄色视频| 最近中文字幕高清免费大全6| 国产成人精品无人区| 久久久久久人人人人人| √禁漫天堂资源中文www| 女人久久www免费人成看片| 九色亚洲精品在线播放| 一区二区av电影网| 国产不卡av网站在线观看| 9191精品国产免费久久| 久久99蜜桃精品久久| 国产日韩欧美在线精品| 国精品久久久久久国模美| 精品人妻偷拍中文字幕| 国产乱人偷精品视频| 黑人欧美特级aaaaaa片| 在线看a的网站| 女人被躁到高潮嗷嗷叫费观| √禁漫天堂资源中文www| 久久国产精品大桥未久av| 一区二区三区激情视频| 自拍欧美九色日韩亚洲蝌蚪91| 岛国毛片在线播放| 777久久人妻少妇嫩草av网站| 日本wwww免费看| 久久亚洲国产成人精品v| 久久国内精品自在自线图片| 少妇被粗大猛烈的视频| 午夜av观看不卡| 国产精品香港三级国产av潘金莲 | 色视频在线一区二区三区| 国产成人精品久久二区二区91 | 欧美97在线视频| 99国产综合亚洲精品| 日本欧美国产在线视频| 国产黄色免费在线视频| 婷婷成人精品国产| 日韩精品有码人妻一区| 在线观看人妻少妇| 高清av免费在线| 黄片小视频在线播放| 一二三四中文在线观看免费高清| 国产有黄有色有爽视频| 午夜激情av网站| 午夜av观看不卡| 97在线人人人人妻| 国产成人精品婷婷| 国产毛片在线视频| 午夜日韩欧美国产| 天天操日日干夜夜撸| 啦啦啦啦在线视频资源| 国产麻豆69| 午夜福利乱码中文字幕| 久久久久久久精品精品| 欧美日韩综合久久久久久| 国产黄色视频一区二区在线观看| 日本色播在线视频| 国产亚洲av片在线观看秒播厂| 国产精品免费大片| 精品人妻熟女毛片av久久网站| 午夜老司机福利剧场| 国产成人精品久久二区二区91 | 久久影院123| 美女脱内裤让男人舔精品视频| 久热久热在线精品观看| 国产视频首页在线观看|