• <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ù)及算例分析
    一本久久精品| 日韩高清综合在线| 日本一二三区视频观看| 国产美女午夜福利| 好男人在线观看高清免费视频| 淫秽高清视频在线观看| 男人舔女人下体高潮全视频| 日韩一本色道免费dvd| 亚洲欧美日韩无卡精品| 成人三级黄色视频| 日韩三级伦理在线观看| 青青草视频在线视频观看| 女人十人毛片免费观看3o分钟| 国产黄色小视频在线观看| 九九热线精品视视频播放| 天美传媒精品一区二区| 亚洲va在线va天堂va国产| 欧美色视频一区免费| 久久人人爽人人片av| 日本黄色片子视频| 久久精品国产鲁丝片午夜精品| 国产黄色小视频在线观看| av卡一久久| 麻豆成人午夜福利视频| 青春草国产在线视频| 久久人人爽人人片av| 蜜桃亚洲精品一区二区三区| 夜夜爽夜夜爽视频| 国内精品一区二区在线观看| 精品熟女少妇av免费看| 国产极品精品免费视频能看的| 亚洲成人av在线免费| 大香蕉久久网| 国产精品1区2区在线观看.| 国产激情偷乱视频一区二区| 亚洲经典国产精华液单| 中文天堂在线官网| 久久亚洲精品不卡| 爱豆传媒免费全集在线观看| 国产成年人精品一区二区| 日日摸夜夜添夜夜添av毛片| 嫩草影院新地址| 两个人的视频大全免费| 伦理电影大哥的女人| 欧美变态另类bdsm刘玥| 中文在线观看免费www的网站| 中文字幕精品亚洲无线码一区| 听说在线观看完整版免费高清| 激情 狠狠 欧美| 亚洲欧洲日产国产| 国产精品爽爽va在线观看网站| 久久久国产成人免费| 国产日韩欧美在线精品| 国产午夜福利久久久久久| 国内揄拍国产精品人妻在线| 六月丁香七月| 国产在视频线在精品| 亚洲国产成人一精品久久久| a级毛片免费高清观看在线播放| 日本免费在线观看一区| 亚洲av中文字字幕乱码综合| 观看免费一级毛片| 国产亚洲精品av在线| 亚洲经典国产精华液单| 白带黄色成豆腐渣| 亚洲中文字幕一区二区三区有码在线看| 国产午夜福利久久久久久| 久久热精品热| 我的老师免费观看完整版| 看十八女毛片水多多多| 可以在线观看毛片的网站| 国产v大片淫在线免费观看| 久久精品国产亚洲网站| 日韩 亚洲 欧美在线| 国产亚洲av嫩草精品影院| 一二三四中文在线观看免费高清| 校园人妻丝袜中文字幕| 最近视频中文字幕2019在线8| 久久精品影院6| 青春草亚洲视频在线观看| 秋霞在线观看毛片| 看免费成人av毛片| 亚洲国产精品专区欧美| 午夜精品在线福利| 97超碰精品成人国产| 亚洲精品日韩在线中文字幕| 日本av手机在线免费观看| 妹子高潮喷水视频| 久久国产精品大桥未久av| 菩萨蛮人人尽说江南好唐韦庄| 嫩草影院入口| 国产xxxxx性猛交| 亚洲一码二码三码区别大吗| 丰满迷人的少妇在线观看| 亚洲精品第二区| 尾随美女入室| 免费播放大片免费观看视频在线观看| 18禁在线无遮挡免费观看视频| 亚洲图色成人| 在线天堂最新版资源| 成人国语在线视频| 韩国av在线不卡| 国产成人一区二区在线| 少妇人妻 视频| 亚洲国产最新在线播放| 成年女人在线观看亚洲视频| 日本av免费视频播放| 国产av国产精品国产| 中文字幕亚洲精品专区| 一级毛片黄色毛片免费观看视频| 国产亚洲午夜精品一区二区久久| 国产无遮挡羞羞视频在线观看| 亚洲成av片中文字幕在线观看 | 国产黄色视频一区二区在线观看| 欧美人与善性xxx| 中文字幕人妻熟女乱码| 免费久久久久久久精品成人欧美视频 | 亚洲,欧美,日韩| av网站免费在线观看视频| 亚洲国产av新网站| 男女高潮啪啪啪动态图| 久久久精品94久久精品| 久久ye,这里只有精品| www日本在线高清视频| 精品一品国产午夜福利视频| 日韩中字成人| 国产精品熟女久久久久浪| 在线观看一区二区三区激情| 国产黄频视频在线观看| 91精品伊人久久大香线蕉| 国产亚洲精品第一综合不卡 | 国产乱来视频区| 亚洲性久久影院| 久久精品国产自在天天线| 国产成人a∨麻豆精品| videos熟女内射| 满18在线观看网站| 中国国产av一级| 久热久热在线精品观看| 国产av精品麻豆| 日日撸夜夜添| 欧美成人午夜精品| 久久久久久久久久久免费av| 一边亲一边摸免费视频| 黄色 视频免费看| 伦理电影大哥的女人| 一本—道久久a久久精品蜜桃钙片| 丝袜脚勾引网站| 少妇被粗大的猛进出69影院 | 亚洲国产精品成人久久小说| 欧美+日韩+精品| 人人妻人人澡人人看| 一本大道久久a久久精品| 80岁老熟妇乱子伦牲交| 国产欧美日韩综合在线一区二区| 亚洲国产色片| 黄色配什么色好看| 大片电影免费在线观看免费| 亚洲精品国产av蜜桃| 国产熟女欧美一区二区| 国产精品女同一区二区软件| 国产色爽女视频免费观看| 日韩av不卡免费在线播放| 久久99热这里只频精品6学生| 亚洲图色成人| 蜜臀久久99精品久久宅男| av国产精品久久久久影院| 涩涩av久久男人的天堂| 日本与韩国留学比较| 99久久综合免费| 免费人成在线观看视频色| 在线观看免费日韩欧美大片| 精品久久久精品久久久| 国产毛片在线视频| 777米奇影视久久| 国产乱来视频区| 国产成人精品在线电影| av天堂久久9| 色5月婷婷丁香| 男女啪啪激烈高潮av片| 18禁在线无遮挡免费观看视频| 久久婷婷青草| av片东京热男人的天堂| 国产无遮挡羞羞视频在线观看| 黑丝袜美女国产一区| 黑人猛操日本美女一级片| 飞空精品影院首页| 日韩一本色道免费dvd| 看免费av毛片| 一区二区日韩欧美中文字幕 | 国产1区2区3区精品| 丝袜喷水一区| 看非洲黑人一级黄片| 国产精品久久久久成人av| 人妻一区二区av| av片东京热男人的天堂| 久热久热在线精品观看| 日韩av在线免费看完整版不卡| 一区二区三区四区激情视频| 午夜免费男女啪啪视频观看| 免费播放大片免费观看视频在线观看| 97精品久久久久久久久久精品| 少妇高潮的动态图| 中文字幕最新亚洲高清| 欧美97在线视频| 日本av免费视频播放| av卡一久久| 亚洲精品久久成人aⅴ小说| 免费在线观看黄色视频的| 免费日韩欧美在线观看| 国产精品.久久久| 国产激情久久老熟女| 久热久热在线精品观看| 这个男人来自地球电影免费观看 | 国产日韩一区二区三区精品不卡| 国产高清三级在线| 亚洲av成人精品一二三区| 亚洲欧美一区二区三区国产| 免费观看在线日韩| 黄色毛片三级朝国网站| 18禁国产床啪视频网站| 国产精品国产三级国产av玫瑰| 99热全是精品| 成人毛片60女人毛片免费| 亚洲综合色网址| 国产精品嫩草影院av在线观看| 国产精品不卡视频一区二区| 午夜av观看不卡| 卡戴珊不雅视频在线播放| 精品视频人人做人人爽| 成人国产av品久久久| 男人添女人高潮全过程视频| 日韩av不卡免费在线播放| 国产视频首页在线观看| 性色av一级| 成人亚洲精品一区在线观看| 美女主播在线视频| 亚洲精品久久午夜乱码| 天天躁夜夜躁狠狠躁躁| 精品少妇黑人巨大在线播放| 少妇人妻久久综合中文| 国产精品一区www在线观看| 亚洲丝袜综合中文字幕| 国产一区二区三区综合在线观看 | 免费看av在线观看网站| 亚洲一区二区三区欧美精品| 国产片特级美女逼逼视频| 欧美日本中文国产一区发布| 两个人看的免费小视频| 人妻 亚洲 视频| 日本vs欧美在线观看视频| 国产69精品久久久久777片| 亚洲国产精品一区二区三区在线| av卡一久久| 尾随美女入室| 国产视频首页在线观看| 又黄又爽又刺激的免费视频.| 五月伊人婷婷丁香| 一级黄片播放器| 免费观看在线日韩| 一本久久精品| 老司机影院毛片| 人人妻人人添人人爽欧美一区卜| kizo精华| √禁漫天堂资源中文www| av线在线观看网站| 五月伊人婷婷丁香| 亚洲精品中文字幕在线视频| 亚洲精品日韩在线中文字幕| 97人妻天天添夜夜摸| 成人免费观看视频高清| 国产高清国产精品国产三级| 赤兔流量卡办理| 各种免费的搞黄视频| 飞空精品影院首页| 视频在线观看一区二区三区| 曰老女人黄片| 国产成人精品无人区| videossex国产| 免费看光身美女| 免费av中文字幕在线| 在线观看人妻少妇| 99热国产这里只有精品6| 日韩 亚洲 欧美在线| 色网站视频免费| 看免费av毛片| 少妇的逼水好多| 国产欧美另类精品又又久久亚洲欧美| av在线app专区| 看非洲黑人一级黄片| 韩国av在线不卡| 寂寞人妻少妇视频99o| 亚洲 欧美一区二区三区| 18禁国产床啪视频网站| av天堂久久9| 国产淫语在线视频| 亚洲精品乱久久久久久| 晚上一个人看的免费电影| 视频在线观看一区二区三区| 最近的中文字幕免费完整| 午夜福利在线观看免费完整高清在| 黄色毛片三级朝国网站| 成年人免费黄色播放视频| 97在线视频观看| 亚洲欧美色中文字幕在线| 精品一区二区三区四区五区乱码 | 视频在线观看一区二区三区| 久久精品夜色国产| 91午夜精品亚洲一区二区三区| 欧美 亚洲 国产 日韩一| 色94色欧美一区二区| 亚洲精品久久午夜乱码| 亚洲精品美女久久av网站| 久久影院123| 免费人成在线观看视频色| 各种免费的搞黄视频| 国产乱人偷精品视频| 欧美人与性动交α欧美精品济南到 | 国产免费一区二区三区四区乱码| 亚洲av电影在线进入| av国产精品久久久久影院| 亚洲性久久影院| 久久ye,这里只有精品| 午夜福利乱码中文字幕| 高清av免费在线| 乱码一卡2卡4卡精品| 美女主播在线视频| 婷婷成人精品国产| 亚洲美女黄色视频免费看| 欧美xxxx性猛交bbbb| 男女下面插进去视频免费观看 | 两性夫妻黄色片 | 纵有疾风起免费观看全集完整版| a级毛片黄视频| 激情五月婷婷亚洲| 国产午夜精品一二区理论片| 国产精品免费大片| 日韩一区二区视频免费看| 美女国产视频在线观看| 精品亚洲成a人片在线观看| 99热国产这里只有精品6| 如何舔出高潮| 99热国产这里只有精品6| 欧美 日韩 精品 国产| 国产一区二区三区av在线| 建设人人有责人人尽责人人享有的| 22中文网久久字幕| 最近的中文字幕免费完整| 国产成人精品无人区| 在现免费观看毛片| 天天操日日干夜夜撸| 久久久久网色| 亚洲一码二码三码区别大吗| 国产精品99久久99久久久不卡 | 免费观看av网站的网址| 又黄又粗又硬又大视频| 国产免费福利视频在线观看| 国产精品嫩草影院av在线观看| 永久网站在线| 在线 av 中文字幕| 人成视频在线观看免费观看| 欧美日韩av久久| 欧美精品人与动牲交sv欧美| 免费黄色在线免费观看| 亚洲精品久久久久久婷婷小说| 国产黄色视频一区二区在线观看| 免费高清在线观看视频在线观看| 丰满饥渴人妻一区二区三| 免费高清在线观看视频在线观看| 成人二区视频| 国产一区二区在线观看日韩| 免费黄网站久久成人精品| 亚洲欧美成人综合另类久久久| 亚洲美女搞黄在线观看| 十分钟在线观看高清视频www| 日本91视频免费播放| 成人国产av品久久久| 丰满饥渴人妻一区二区三| 免费人妻精品一区二区三区视频| 国产一级毛片在线| 又大又黄又爽视频免费| 伊人久久国产一区二区| 各种免费的搞黄视频| 亚洲第一av免费看| 亚洲综合色网址| 成年动漫av网址| 免费av中文字幕在线| 18+在线观看网站| 青春草亚洲视频在线观看| 人妻一区二区av| 日本爱情动作片www.在线观看| 成年人免费黄色播放视频| 久久 成人 亚洲| 人人妻人人澡人人爽人人夜夜| 大片免费播放器 马上看| 欧美最新免费一区二区三区| 夜夜骑夜夜射夜夜干| 日本-黄色视频高清免费观看| 老司机影院成人| 一级片免费观看大全| 51国产日韩欧美| 日韩av在线免费看完整版不卡| 国产日韩欧美视频二区| 中文字幕人妻熟女乱码| 精品人妻一区二区三区麻豆| 亚洲av福利一区| 精品一品国产午夜福利视频| 成人漫画全彩无遮挡| 亚洲性久久影院| 国产av国产精品国产| 90打野战视频偷拍视频| 国产精品久久久久成人av| 亚洲精品aⅴ在线观看| 国产黄色视频一区二区在线观看| 亚洲美女黄色视频免费看| 女的被弄到高潮叫床怎么办| 亚洲国产精品专区欧美| 男女无遮挡免费网站观看| 一级a做视频免费观看| 少妇被粗大的猛进出69影院 | 成人午夜精彩视频在线观看| 两个人看的免费小视频| 国产白丝娇喘喷水9色精品| 精品卡一卡二卡四卡免费| 在线天堂最新版资源| 五月天丁香电影| 国产永久视频网站| 免费看av在线观看网站| 亚洲国产av影院在线观看| 国产精品麻豆人妻色哟哟久久| 欧美精品国产亚洲| kizo精华| 国产免费视频播放在线视频| 肉色欧美久久久久久久蜜桃| 日韩 亚洲 欧美在线| 中文欧美无线码| 久久久久久久国产电影| 黑人高潮一二区| 久久这里有精品视频免费| 中国国产av一级| 九草在线视频观看| 国产在视频线精品| 日韩人妻精品一区2区三区| 国产69精品久久久久777片| 国产极品粉嫩免费观看在线| 国产高清三级在线| 成人影院久久| 国产一区二区激情短视频 | 国产日韩欧美亚洲二区| 亚洲色图 男人天堂 中文字幕 | 精品久久久久久电影网| 国产黄色视频一区二区在线观看| 全区人妻精品视频| 好男人视频免费观看在线| 精品少妇黑人巨大在线播放| 美女xxoo啪啪120秒动态图| 日韩中字成人| 亚洲欧洲日产国产| av在线app专区| 视频区图区小说| 黄网站色视频无遮挡免费观看| 成人毛片60女人毛片免费| 亚洲欧美清纯卡通| 久久女婷五月综合色啪小说| av国产精品久久久久影院| 久久久久久久久久成人| 国产亚洲最大av| 精品人妻在线不人妻| 国产男人的电影天堂91| 少妇精品久久久久久久| 精品久久久久久电影网| 国产爽快片一区二区三区| 婷婷色综合大香蕉| 青春草国产在线视频| 22中文网久久字幕| 曰老女人黄片| 边亲边吃奶的免费视频| 国产老妇伦熟女老妇高清| 国产成人精品久久久久久| 亚洲精品国产av成人精品| 国产av精品麻豆| 精品久久久精品久久久| 日韩一本色道免费dvd| 亚洲国产欧美日韩在线播放| 黑人高潮一二区| 久久亚洲国产成人精品v| 毛片一级片免费看久久久久| 男的添女的下面高潮视频| 丝袜在线中文字幕| 精品久久蜜臀av无| 久久久久久久久久成人| 一本大道久久a久久精品| 婷婷色麻豆天堂久久| 蜜桃国产av成人99| 国产精品人妻久久久影院| 久久亚洲国产成人精品v| 美女内射精品一级片tv| 一级片免费观看大全| 欧美日韩亚洲高清精品| 五月玫瑰六月丁香| 国产成人一区二区在线| 中文字幕精品免费在线观看视频 | 亚洲国产精品成人久久小说| 国产熟女午夜一区二区三区| 在线观看一区二区三区激情| 久久久久人妻精品一区果冻| 免费高清在线观看日韩| 亚洲精品成人av观看孕妇| 少妇人妻 视频| 久久这里有精品视频免费| 国精品久久久久久国模美| 久久99一区二区三区| 999精品在线视频| 人人澡人人妻人| 伦理电影大哥的女人| 久久久久精品性色| av福利片在线| 亚洲精品久久成人aⅴ小说| 啦啦啦中文免费视频观看日本| 免费观看在线日韩| 亚洲精品中文字幕在线视频| 亚洲av福利一区| 国产一区二区三区综合在线观看 | 两性夫妻黄色片 | 精品国产一区二区久久| 80岁老熟妇乱子伦牲交| 美女内射精品一级片tv| 丝袜人妻中文字幕| 熟妇人妻不卡中文字幕| 国产精品国产三级国产av玫瑰| 熟女av电影| videossex国产| 最新的欧美精品一区二区| 伊人久久国产一区二区| 成人18禁高潮啪啪吃奶动态图| 最近中文字幕2019免费版| 人妻系列 视频| av免费在线看不卡| 人人妻人人爽人人添夜夜欢视频| 大香蕉97超碰在线| 69精品国产乱码久久久| 国产色婷婷99| 久久国产亚洲av麻豆专区| 一级片免费观看大全| 黑人巨大精品欧美一区二区蜜桃 | 久久精品人人爽人人爽视色| 最近的中文字幕免费完整| 菩萨蛮人人尽说江南好唐韦庄| 免费看不卡的av| 精品国产国语对白av| 视频在线观看一区二区三区| 新久久久久国产一级毛片| 久久热在线av| 久久人妻熟女aⅴ| 啦啦啦啦在线视频资源| 日韩三级伦理在线观看| 在线天堂中文资源库| 97人妻天天添夜夜摸| 国产精品秋霞免费鲁丝片| 一个人免费看片子| 色哟哟·www| 日韩 亚洲 欧美在线| 国产伦理片在线播放av一区| 国产高清国产精品国产三级| 久久亚洲国产成人精品v| 校园人妻丝袜中文字幕| 国产在线免费精品| 久久精品国产综合久久久 | 观看美女的网站| 色94色欧美一区二区| 国产亚洲精品久久久com| 人成视频在线观看免费观看| 国产日韩欧美视频二区| 69精品国产乱码久久久| 男女午夜视频在线观看 | 亚洲精品久久成人aⅴ小说| 成人漫画全彩无遮挡| 亚洲中文av在线| 亚洲四区av| 日韩av不卡免费在线播放| 另类精品久久| 国产麻豆69| 日韩视频在线欧美| 国产日韩一区二区三区精品不卡| 国产成人精品福利久久| 精品国产国语对白av| 日韩人妻精品一区2区三区| 男女午夜视频在线观看 | av.在线天堂| 午夜精品国产一区二区电影| 激情五月婷婷亚洲| 欧美日韩一区二区视频在线观看视频在线| 欧美亚洲日本最大视频资源| 男女边吃奶边做爰视频| 国产男女超爽视频在线观看| 亚洲精品中文字幕在线视频| 多毛熟女@视频| 在线天堂最新版资源| 欧美日韩亚洲高清精品| 国产视频首页在线观看| 永久网站在线| 亚洲图色成人| 黄色怎么调成土黄色| 中文天堂在线官网| videos熟女内射| 日日爽夜夜爽网站| 久久99精品国语久久久| 亚洲三级黄色毛片| 美女内射精品一级片tv| 亚洲色图 男人天堂 中文字幕 | 国产欧美另类精品又又久久亚洲欧美| 黄色一级大片看看| 一边摸一边做爽爽视频免费| 日韩三级伦理在线观看| 国产成人精品在线电影|