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

    基于偽譜凸優(yōu)化和L1罰函數(shù)的彈道規(guī)劃方法研究

    2022-03-28 07:43:06王慶海王中原尹秋霖
    彈道學報 2022年1期
    關(guān)鍵詞:偽譜信賴彈道

    王慶海,陳 琦,王中原,尹秋霖

    (南京理工大學 能源與動力工程學院,江蘇 南京 210094)

    軌跡優(yōu)化方法按照是否得到解析解被分為2大類:解析法和數(shù)值解法。解析法基于最優(yōu)控制原理,針對實際問題模型,根據(jù)極值條件和邊界條件推導出解析解。史金光等利用龐特里亞金極小值原理,設(shè)計了滑翔制導炮彈以最大射程為指標的最優(yōu)控制參數(shù)。解析法對于簡單系統(tǒng)較為有效,但難以求解復雜的多約束非線性優(yōu)化問題。數(shù)值解法通過某種離散方法,對時間、狀態(tài)和控制量進行離散,將非線性最優(yōu)控制問題轉(zhuǎn)化為有限參數(shù)規(guī)劃問題,然后采用合適的參數(shù)規(guī)劃算法求出最優(yōu)解。直接打靶法是數(shù)值解法中最常用的一種。涂良輝等采用直接打靶法規(guī)劃再入飛行器無動力跳躍軌跡。李瑜采用直接打靶法,對助推-滑翔導彈的最大射程和可達區(qū)域進行了分析。直接打靶法只對控制量進行離散,其優(yōu)化參數(shù)少,便于理解,但是直接打靶法需要使用高階積分算法求解狀態(tài)變量,計算量大,效率不高。根據(jù)FRABIEN的研究,直接打靶法對初始猜測敏感度很高。偽譜法是另一種最常用的數(shù)值解法。陳琦等針對滑翔制導炮彈,采用Gauss偽譜法離散彈道,以最短飛行時間為性能指標,利用序列二次規(guī)劃算法(SQP)規(guī)劃了滑翔段彈道。黃詰等通過Radau偽譜法離散導彈軌跡,分別取最大終端速度和最大落角為性能指標,采用SQP方法規(guī)劃出最優(yōu)攻擊軌跡。張科南等利用SQP方法,為高超聲速飛行器規(guī)劃了多約束突防彈道。偽譜離散與序列二次規(guī)劃法結(jié)合的偽譜法在求解非線性最優(yōu)化問題時十分有效,但其求解速度較慢。造成這一問題的主要原因是序列二次規(guī)劃方法在每一次迭代中不僅要更新Jacobian矩陣,還需要更新Hessian矩陣。

    凸問題在數(shù)學上已經(jīng)被證明,使用內(nèi)點法可以保證在多項式時間內(nèi)獲得全局最優(yōu)解。由于這一固有的優(yōu)點,近年來,利用凸優(yōu)化求解軌跡優(yōu)化問題成為研究熱點。陳嘉澍等采用Chebyshev偽譜法離散高超聲速飛行器動態(tài)方程,結(jié)合無損凸化方法,將原問題轉(zhuǎn)化為標準二階錐規(guī)劃問題求解。ACIKMESE等以最小燃料消耗為指標,利用序列凸優(yōu)化方法規(guī)劃最優(yōu)火星著陸軌跡。同樣以火星著陸為研究背景,文獻[13-14]以最小著陸誤差為性能指標,利用凸優(yōu)化算法規(guī)劃著陸軌跡。王金波等利用翻轉(zhuǎn)Radau偽譜法離散一級火箭回收軌跡,將空心非凸推力約束轉(zhuǎn)換為二階錐約束,以燃料消耗量為性能指標,通過序列凸優(yōu)化方法得到最優(yōu)回收軌跡。文獻[16]采用梯形離散方法離散連續(xù)變量,利用序列凸優(yōu)化方法規(guī)劃了空地導彈的最大存速軌跡,但精度較差。ROH等采用梯形離散方法,離散導彈簡化三自由度模型,將動態(tài)方程線性凸化,采用帶L1懲罰的序列凸優(yōu)化算法(LPSCP)實現(xiàn)實時彈道規(guī)劃。

    文獻[9]提出的(LPSCP)算法信賴域?qū)挾容^大,且恒定;算法中L1懲罰系數(shù)初值較大,且呈指數(shù)增長。該算法在解決簡單軌跡規(guī)劃問題時表現(xiàn)良好,但面對較復雜的彈道規(guī)劃問題時表現(xiàn)較差。較大的信賴域?qū)挾仍斐删€性近似精度低;指數(shù)增長的懲罰系數(shù)會在迭代過程中引起浮點數(shù)計算誤差,尤其制導炮彈軌跡優(yōu)化問題中包含大量小數(shù)量級決策變量,嚴重時會導致無法求解。

    基于上述考慮,本文提出一種改進的L1懲罰序列凸規(guī)劃算法(ILPSCP),算法中采用指數(shù)衰減策略對相對信賴域?qū)挾冗M行更新,并采用指數(shù)增長帶上界的L1懲罰系數(shù)。仿真結(jié)果顯示:與傳統(tǒng)L1懲罰序列凸規(guī)劃算法(LPSCP)相比,本文提出的ILPSCP算法收斂速度更快,穩(wěn)定性更強。

    1 建模

    本文以控制能量最優(yōu)為性能指標,規(guī)劃制導炮彈打擊固定目標的滑翔彈道。為實現(xiàn)一定的打擊效能,要求落角不大于給定上界值,末速度應不小于給定下界值。

    1.1 制導炮彈縱向平面內(nèi)滑翔段彈道模型

    本文對制導炮彈縱向二維平面內(nèi)方案彈道進行規(guī)劃,對應的動態(tài)方程為

    (1)

    式中:為速度,為彈道傾角,為阻力,為升力。阻力和升力的計算公式為

    式中:為特征面積,0為零升阻力系數(shù),為誘導阻力系數(shù),為升力系數(shù)導數(shù),為攻角。

    1.2 模型約束

    假設(shè)制導炮彈在滑翔起控點速度為,起控點彈道傾角為,起控點坐標為(,),目標點坐標為(,)。為實現(xiàn)較好的打擊效能,末速度約束下界為,落角約束上界為。

    于是建立初始狀態(tài)約束為

    (2)

    末狀態(tài)約束:

    (3)

    考慮彈體穩(wěn)定性,攻角不能過大,建立控制量約束:

    ≤≤

    (4)

    1.3 最優(yōu)化模型

    參考文獻[5],以滑翔能量消耗最低為性能指標,對應二維方案彈道的性能指標函數(shù)為

    需要注意的是,初始時間和終止時間是可變的,也是問題的決策變量。

    綜上,制導炮彈滑翔軌跡最優(yōu)化問題模型為

    2 偽譜離散和凸化

    2.1 模型一般化

    為了方便后續(xù)公式推導,本文先建立制導炮彈滑翔段控制能量最優(yōu)軌跡優(yōu)化模型對應的一般化模型。一般化模型中狀態(tài)變量、控制變量以及狀態(tài)函數(shù)都以向量形式表示。

    狀態(tài)向量×1,控制向量=×1,本文中=4,=1。狀態(tài)向量:

    =()

    控制向量為

    =

    于是動態(tài)方程可以寫為

    (5)

    式中:×1為狀態(tài)函數(shù)向量。

    目標函數(shù)可以寫為

    (6)

    式中:()為控制向量的Euclidean范數(shù)。狀態(tài)向量、控制向量和時間變量的上下界約束可以寫為

    (7)

    (8)

    ≤≤

    (9)

    式中:L表示下邊界約束,U表示上邊界約束。

    綜上,制導炮彈滑翔段控制能量最優(yōu)軌跡規(guī)劃模型對應的一般化最優(yōu)化模型為

    2.2 Radau偽譜離散

    偽譜法采用全局多項式插值,相對于傳統(tǒng)的等距梯形離散方法,偽譜離散可以獲得更高的精度。Radau偽譜法使用階Radau勒讓德多項式的根為節(jié)點,其節(jié)點是非對稱、非等間距的,采用它的根作為插值節(jié)點可以避免等距高階多項式插值經(jīng)常出現(xiàn)的龍格現(xiàn)象。相比Gauss偽譜離散,Radau偽譜離散沒有附加的積分約束,更為簡潔。Radau多項式是階勒讓德多項式和-1階勒讓德多項式-1的差:

    ()=()--1(),∈[-1,1]

    階Radau多項式有個根,這個根包含=-1,<1??梢钥闯?配置點包含初始點,但不包含終止點。

    定義∈[-1,1] 為偽譜時間,∈[,]為物理時間。偽譜時間和物理時間之間存在如下映射關(guān)系:

    (10)

    Radau偽譜法的節(jié)點包含階Radau多項式的個根和一個附加終止點,共計+1個節(jié)點。

    定義為第個節(jié)點對應的拉格朗日多項式函數(shù)的基函數(shù),寫作:

    于是,狀態(tài)向量的第個分量()可以用經(jīng)過這+1個節(jié)點的拉格朗日多項式函數(shù)來近似。

    (11)

    對式(11)兩邊求導,可得:

    于是第個狀態(tài)變量在第個配置點處微分的值為

    (12)

    式中:(+1)×為狀態(tài)矩陣,為第個狀態(tài)變量在第個離散節(jié)點的狀態(tài)。在這里定義控制矩陣為×,為第個控制變量在第個離散節(jié)點的控制量。×(+1)被稱為微分矩陣,其每一行對應著相應配置點,每一列對應著相應拉格朗日基函數(shù)。

    對式(10)兩邊求導,并代入式(5)可得:

    定義偽譜域函數(shù):

    即:

    離散形式為

    (13)

    式中:×

    =(,,,)

    式中矩陣下標對應第個離散節(jié)點。結(jié)合式(12)和式(13)可得:

    =(,,,)

    (14)

    由式(14)可知,本文研究的最優(yōu)化問題的決策變量包含:狀態(tài)變量、控制量、起始時間和結(jié)束時間。為了方便數(shù)值求解,需要根據(jù)標準數(shù)值規(guī)劃模型,將所有決策變量轉(zhuǎn)換成列向量。

    定義決策變量列向量為

    =((vec())(vec()))

    式中:vec表示矩陣列變換操作,×1,為決策變量總數(shù)。

    =(+1)++2

    式中:(+1)對應狀態(tài)決策變量總數(shù),對應控制決策變量總數(shù),常數(shù)2對應初末時間,為控制變量數(shù)目。

    定義狀態(tài)函數(shù)列向量為×1:

    ()=vec(())

    于是式(14)可以改寫為

    =()

    (15)

    =(diag(,)()×(+2))

    Radau偽譜法的離散求積公式為

    式中:為被積函數(shù)在第個配置點的值;為對應配置點的積分權(quán)系數(shù)。于是式(6)的離散形式為

    式中:表示第個控制變量。引入輔助變量建立二階錐約束:

    (16)

    (17)

    2.3 線性凸化

    (18)

    =

    (19)

    式中:(·)×,(·)×1。

    彈道規(guī)劃問題中,決策變量分為狀態(tài)變量、控制變量和時間變量3個大類。每一個大類之間數(shù)量級相差較大,而大類之中決策變量的數(shù)量級也相差懸殊。本文采用比例縮放的思想,指定不同決策變量的縮放比例系數(shù)。比例縮放系數(shù)定義為決策變量上邊界值與下邊界值的差,即:

    設(shè)定統(tǒng)一的相對信賴域?qū)挾?從而實現(xiàn)簡單地調(diào)整一個相對信賴域?qū)挾?完成對所有決策變量信賴域?qū)挾鹊恼{(diào)整。

    相對信賴域約束定義為

    (20)

    式中:相對信賴域?qū)挾葹橐粋€足夠小的數(shù)。

    同理,定義相對誤差為

    式(18)、式(20)的物理意義為:新軌跡點在原軌跡點附近,由于軌跡由軌跡點連接而成,即新軌跡在原軌跡附近,從變分學角度看,新軌跡在原軌跡的鄰域內(nèi)。

    式(17)為線性目標函數(shù),式(19)為線性等式約束方程,式(16)和式(20)為錐約束,于是式(16)~式(17)、式(19)~式(20)組成了標準凸優(yōu)化模型,稱該凸優(yōu)化模型為原子問題:

    3 改進L1懲罰序列凸規(guī)劃算法

    由文獻[9]可知,采用普通序列凸規(guī)劃算法(SCP)求解彈道規(guī)劃問題時,其對初始預測十分敏感。如果不能提供合適的初始預測,子問題的解空間將為空,無法求解。合適的初始預測依靠豐富的經(jīng)驗,當問題較為復雜時,經(jīng)驗預測也難以適用。

    L1懲罰序列凸規(guī)劃算法(LPSCP)成功解決了序列凸規(guī)劃算法(SCP)對初始預測敏感這個問題。LPSCP算法通過引入松弛輔助變量,,,,在原問題的基礎(chǔ)上建立L1懲罰子問題。子問題定義為

    式中:(·)×1(·)×1為一組用于松弛線性等式約束方程的輔助變量;×1×1為一組用于松弛信賴域不等式約束的輔助變量;為違反等式約束的懲罰系數(shù),為違反信賴域不等式約束的懲罰系數(shù);=1,2,…,。

    2個懲罰系數(shù)的獲得方式為

    式中:為違反原子問題約束時,隨著迭代次數(shù)呈指數(shù)增長的總懲罰系數(shù),增長速率由系數(shù)確定?!蔥0,1]為一個權(quán)系數(shù),代表對2種約束懲罰權(quán)重。越大,對違反動態(tài)方程等式約束的懲罰越大,對違反信賴域不等式約束懲罰越少,反之同理。當=1時,代表只對動態(tài)約束采取L1懲罰;當=0時,代表只對信賴域約束采取L1懲罰。的具體取值根據(jù)實際情況而定,本文研究的制導炮彈軌跡優(yōu)化問題動態(tài)約束直接影響初始預測的敏感性,相比之下更為重要,所以懲罰權(quán)系數(shù)取值最好不小于0.5。傳統(tǒng)L1懲罰序列凸規(guī)劃(LPSCP)算法流程如圖1所示。

    圖1 傳統(tǒng)LPSCP算法流程圖

    從式(18)、式(20)可以看出,信賴域?qū)挾炔粌H影響線性近似的精度,還將決定每一次迭代優(yōu)化的最大收斂步長。相對信賴域?qū)挾仍叫?則最大收斂步長越小,但線性近似精度越高;相對信賴域?qū)挾仍酱?則迭代最大收斂步長越大,但線性近似精度越低。傳統(tǒng)序列凸規(guī)劃算法使用固定信賴域?qū)挾?如果相對信賴域?qū)挾忍?則無法滿足線性近似條件,會導致算法失效。為了兼具收斂速度和近似精度,改進的L1懲罰序列凸規(guī)劃(ILPSCP)算法采用指數(shù)衰減策略更新相對信賴域?qū)挾?

    =max(×,)

    式中:為初始相對信賴域?qū)挾?由于采用了指數(shù)衰減策略,可以取得適當大一些,有利于前期迭代快速收斂;為衰減系數(shù),越大相對信賴域?qū)挾人p越慢,的取值一般在08~098之間;為迭代次數(shù);為相對信賴域?qū)挾认陆?。算法的終止條件為前后兩次迭代決策變量之間的相對差值不大于相對誤差容錯系數(shù)。信賴域?qū)挾仍O(shè)置下界的原因在于相對信賴域?qū)挾缺仨毚笥诘扔谙鄬φ`差容錯系數(shù),即:

    否則當<,會導致解的精度并未達到預期,迭代過程便提前結(jié)束。

    理論上,當總懲罰系數(shù)足夠大時,子問題和子問題求得的解相同。然而在研究一個數(shù)值算法時,有必要考慮計算機浮點數(shù)計算精度。計算機運算中,數(shù)量級差距太大的數(shù)進行運算時可能會出現(xiàn)誤差,這對于決策變量本身數(shù)量級差距較大的彈道規(guī)劃問題是不利的,懲罰系數(shù)太大時,將無法求解。如圖2所示,當懲罰系數(shù)達到10時,求得的所有的決策變量都在0附近,顯然不是期望的彈道,是無效解。

    圖2 懲罰系數(shù)cex=107 時的仿真結(jié)果

    綜合考慮,為懲罰系數(shù)設(shè)定一個上界,使懲罰系數(shù)滿足:

    =min(×,)

    改進的L1懲罰序列凸規(guī)劃(ILPSCP)算法流程如圖3所示。

    圖3 ILPSCP算法流程圖

    4 仿真和分析

    本節(jié)以制導炮彈縱向平面內(nèi)滑翔段彈道規(guī)劃問題為例,分別采用傳統(tǒng)L1懲罰序列凸規(guī)劃算法(LPSCP)和本文提出的改進L1懲罰序列凸規(guī)劃算法(ILPSCP)進行仿真對比。為驗證本文提出的ILPSCP算法的準確性,使用非線性問題最優(yōu)化通用工具箱GPOPS2進行仿真,作為參照組。

    4.1 仿真參數(shù)設(shè)置

    制導炮彈參數(shù)如表1所示。表2為各決策變量的上下邊界值。

    表1 制導炮彈參數(shù)

    表2 決策變量的上下邊界值

    直接打靶法對初始預測敏感,需要遴選較好的初始預測才能求解,然而本文提出的ILPSCP算法具備強大的冷啟動能力。正常情況下,彈道的始末端點坐標預測值會選擇制導炮彈滑翔起點坐標和目標的坐標,這種初始預測經(jīng)過拉格朗日插值得到的初始預測彈道軌跡比較接近期望的最優(yōu)彈道軌跡,有利于算法的求解。表3為本文仿真使用的始末變量預測值,下標“0”對應變量起始值,下標“f”對應變量末值,下標“g”表示預測值。為了驗證ILPSCP算法的冷啟動能力,表3提供的始末端點坐標預測值極其惡劣,是2個重合的點,并且這個重合點即不是制導炮彈滑翔起點坐標,也不是目標坐標。

    表3 始末變量初始預測值

    表4中為彈道規(guī)劃實例參數(shù),主要包含各狀態(tài)變量的初始點約束、終點約束以及規(guī)劃算法中需要的各參數(shù)。為保證制導炮彈在飛行末段具備一定的機動能力,設(shè)置飛行速度下界為200 m/s;考慮到制導炮彈飛行穩(wěn)定性;攻角絕對值最大值為π/12;其他參數(shù)見表4。

    表4 實例仿真參數(shù)表

    4.2 仿真結(jié)果與分析

    仿真結(jié)果如圖4~圖9所示。圖4為飛行軌跡對比圖;圖5為飛行速度對比圖;圖6為飛行過程中彈道傾角對比圖;圖7為攻角曲線對比圖,也是該模型的控制曲線;圖8為本文提出的ILPSCP算法和傳統(tǒng)LPSCP算法的目標函數(shù)收斂曲線對比圖;圖9為GPOPS2和ILPSCP算法解的節(jié)點分布圖。

    圖4 飛行軌跡對比圖

    圖5 速度曲線對比圖

    圖6 彈道傾角曲線對比圖

    圖7 攻角曲線對比圖

    圖8 目標函數(shù)收斂曲線圖

    圖9 離散節(jié)點分布對比圖

    圖4~圖7中,給出了傳統(tǒng)LPSCP算法第199次和第200次迭代的結(jié)果,可以看出,這兩組曲線都沒能和參照組曲線重合,即沒有收斂到最優(yōu)解。結(jié)合圖8可以看出,LPSCP算法對應的目標函數(shù),僅需20次迭代便趨于一個值附近,并在該值附近震蕩。但在經(jīng)過50次、100次甚至150次迭代之后,該目標函數(shù)曲線仍然持續(xù)震蕩,難以收斂,最后在0.607和0.690之間持續(xù)跳變??梢灶A計,即使進行更多迭代仍然難以收斂。形成鮮明對比的是,圖8中ILPSCP算法對應的目標函數(shù)曲線在經(jīng)過30次迭代后便基本穩(wěn)定,之后變化量極小。

    由圖4~圖7和圖9可以看出,ILPSCP算法的節(jié)點和GPOPS2的節(jié)點雖然不同,但ILPSCP算法規(guī)劃出的各變量曲線與參照組GPOPS2高度重合,這證明本文提出的ILPSCP算法是有效的。

    本文的彈道規(guī)劃模型的性能指標為攻角平方積分,可以發(fā)現(xiàn),圖7中LPSCP算法的2條攻角曲線的平方積分比其他2種算法的更小,從圖8中也可以看出,LPSCP算法的性能指標小于ILPSCP算法,這與本文描述的ILPSCP算法優(yōu)于LPSCP算法這一觀點不符。出現(xiàn)這種現(xiàn)象的原因是:LPSCP算法使用固定相對信賴域=03,較大的信賴域?qū)挾葘τ诤唵螁栴}影響不大,但對于本文較復雜的彈道規(guī)劃問題,會造成線性近似假設(shè)條件不成立,近似誤差大。過大的線性近似誤差導致近似模型無法反映制導炮彈物理模型,所以,對于制導炮彈模型,圖7中LPSCP算法的攻角曲線和圖8中LPSCP算法的目標函數(shù)曲線都是無效的。過大的線性近似誤差,也是導致LPSCP算法無法收斂的原因。采用=03的LPSCP算法在2個無效解之間跳變,而=03的ILPSCP算法收斂到最優(yōu)解,顯然ILPSCP算法比LPSCP算法更適合解復雜的彈道規(guī)劃問題。

    LPSCP算法在=03時便因為信賴域?qū)挾冗^大無法求解,而ILPSCP算法在更大的初始相對信賴域?qū)挾惹闆r下,仍然能穩(wěn)定求解。增大后,ILPSCP算法收斂到指定精度(最優(yōu)解與=03時相同)所需迭代次數(shù),如表5所示。隨著初始相對信賴域的增大,所需迭代次數(shù)增加,但相比信賴域?qū)挾鹊脑鲩L倍數(shù),其增長幅度較小。相比于傳統(tǒng)的LPSCP算法,即使設(shè)置了更大的初始相對信賴域?qū)挾?ILPSCP算法仍然能快速穩(wěn)定收斂,這進一步證明了ILPSCP算法的優(yōu)越性。

    表5 不同δ0下ILPSCP算法迭代次數(shù)

    5 結(jié)束語

    本文推導了時間自由且以控制能量最優(yōu)為指標的彈道規(guī)劃一般化模型,針對該模型采用偽譜離散和線性凸化完成凸優(yōu)化建模。針對傳統(tǒng)L1序列凸規(guī)劃算法(LPSCP)在解決復雜的初末時間自由、非線性強的彈道規(guī)劃問題時存在穩(wěn)定性較差等問題,本文提出了一種改進的L1偽譜凸規(guī)劃算法(ILPSCP)。ILPSCP算法中引入指數(shù)衰減的相對信賴域?qū)挾群蛶辖绲腖1懲罰系數(shù)。仿真結(jié)果證明,ILPSCP算法收斂速度快,精度高,冷啟動能力強,很好地解決了傳統(tǒng)L1序列凸規(guī)劃算法(LPSCP)存在的問題。

    猜你喜歡
    偽譜信賴彈道
    彈道——打勝仗的奧秘
    矩陣偽譜的新定位集及其在土壤生態(tài)系統(tǒng)的應用
    一維彈道修正彈無線通信系統(tǒng)研制
    電子制作(2019年7期)2019-04-25 13:17:48
    淺談行政法的信賴利益保護原則
    信賴利益保護原則的中國化
    行政法論叢(2018年1期)2018-05-21 00:41:50
    紊流環(huán)境下四維軌跡優(yōu)化的偽譜方法研究
    飛行力學(2017年1期)2017-02-15 03:00:58
    一種改進的自適應信賴域算法
    基于PID控制的二維彈道修正彈仿真
    制導與引信(2016年3期)2016-03-20 16:02:02
    偽譜法及其在飛行器軌跡優(yōu)化設(shè)計領(lǐng)域的應用綜述*
    消除彈道跟蹤數(shù)據(jù)中伺服系統(tǒng)的振顫干擾
    久久ye,这里只有精品| 视频区欧美日本亚洲| 日韩熟女老妇一区二区性免费视频| 天天操日日干夜夜撸| 久久中文字幕一级| 另类亚洲欧美激情| 欧美精品高潮呻吟av久久| 久久国产精品人妻蜜桃| 亚洲免费av在线视频| 国产精华一区二区三区| 国产精品一区二区在线不卡| 激情视频va一区二区三区| av视频免费观看在线观看| 亚洲av第一区精品v没综合| 很黄的视频免费| 日日爽夜夜爽网站| 12—13女人毛片做爰片一| 亚洲,欧美精品.| 日韩人妻精品一区2区三区| 久久精品国产99精品国产亚洲性色 | 村上凉子中文字幕在线| 亚洲伊人色综图| 村上凉子中文字幕在线| 三上悠亚av全集在线观看| 国产一区在线观看成人免费| 狠狠狠狠99中文字幕| 十分钟在线观看高清视频www| 一区二区三区国产精品乱码| 国产野战对白在线观看| 亚洲七黄色美女视频| 最新在线观看一区二区三区| www日本在线高清视频| 精品久久久久久,| 99精国产麻豆久久婷婷| e午夜精品久久久久久久| 免费观看精品视频网站| 国产一区二区三区视频了| 午夜影院日韩av| 丁香六月欧美| 天天躁日日躁夜夜躁夜夜| 亚洲在线自拍视频| 十分钟在线观看高清视频www| 亚洲成a人片在线一区二区| 女人被狂操c到高潮| 久久久国产成人精品二区 | 男女高潮啪啪啪动态图| 色播在线永久视频| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲精品第一综合不卡| 大型av网站在线播放| 精品欧美一区二区三区在线| 国产精品免费视频内射| 中文字幕人妻丝袜一区二区| cao死你这个sao货| 国产极品粉嫩免费观看在线| 亚洲aⅴ乱码一区二区在线播放 | 少妇粗大呻吟视频| 变态另类成人亚洲欧美熟女 | 精品福利永久在线观看| 欧美日韩亚洲综合一区二区三区_| 天天添夜夜摸| 久久中文字幕人妻熟女| 99riav亚洲国产免费| 91国产中文字幕| 国产97色在线日韩免费| 国产成人欧美在线观看 | 久热爱精品视频在线9| 三上悠亚av全集在线观看| 极品少妇高潮喷水抽搐| 婷婷丁香在线五月| 国产精品偷伦视频观看了| 成年人黄色毛片网站| 久99久视频精品免费| 欧美一级毛片孕妇| 首页视频小说图片口味搜索| 视频区图区小说| 大片电影免费在线观看免费| 最近最新中文字幕大全电影3 | 久久人妻熟女aⅴ| 久久天堂一区二区三区四区| 午夜福利乱码中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 亚洲中文字幕日韩| 欧美午夜高清在线| 亚洲成av片中文字幕在线观看| 精品国产超薄肉色丝袜足j| 纯流量卡能插随身wifi吗| 在线视频色国产色| 欧美午夜高清在线| 亚洲精品美女久久久久99蜜臀| 午夜福利一区二区在线看| 99精品欧美一区二区三区四区| 在线天堂中文资源库| 亚洲欧美一区二区三区黑人| 性色av乱码一区二区三区2| 午夜福利在线观看吧| 伦理电影免费视频| 纯流量卡能插随身wifi吗| 一个人免费在线观看的高清视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美性长视频在线观看| 精品乱码久久久久久99久播| 性少妇av在线| 又紧又爽又黄一区二区| 欧美激情极品国产一区二区三区| 成人国语在线视频| 在线观看免费视频日本深夜| 国产精品九九99| av片东京热男人的天堂| 久热爱精品视频在线9| 亚洲av美国av| 国产精品亚洲av一区麻豆| 女性被躁到高潮视频| 欧美黑人精品巨大| 欧美日韩一级在线毛片| 女性被躁到高潮视频| 欧美在线黄色| 亚洲精品美女久久av网站| 国产精品免费一区二区三区在线 | 久久香蕉国产精品| 99国产精品免费福利视频| av有码第一页| 黑丝袜美女国产一区| 露出奶头的视频| 操美女的视频在线观看| 一级a爱视频在线免费观看| 欧美黄色片欧美黄色片| 99热网站在线观看| 建设人人有责人人尽责人人享有的| 午夜福利,免费看| 欧美色视频一区免费| 色综合婷婷激情| 波多野结衣av一区二区av| 精品福利永久在线观看| 免费观看精品视频网站| 丁香欧美五月| 99国产精品一区二区蜜桃av | 黄色a级毛片大全视频| 亚洲精品av麻豆狂野| 韩国av一区二区三区四区| 免费日韩欧美在线观看| 一二三四在线观看免费中文在| 午夜福利一区二区在线看| 后天国语完整版免费观看| 亚洲国产欧美网| 搡老岳熟女国产| 亚洲五月天丁香| 亚洲在线自拍视频| 国产在线观看jvid| 69精品国产乱码久久久| 成年人午夜在线观看视频| 国产欧美日韩一区二区三| 一级作爱视频免费观看| 亚洲一区中文字幕在线| 99在线人妻在线中文字幕 | 精品人妻在线不人妻| 波多野结衣一区麻豆| 国产激情欧美一区二区| 国产又色又爽无遮挡免费看| 国产欧美日韩综合在线一区二区| 老熟妇乱子伦视频在线观看| 五月开心婷婷网| 久久精品国产99精品国产亚洲性色 | 国产人伦9x9x在线观看| 欧美人与性动交α欧美软件| 久久香蕉精品热| 啦啦啦免费观看视频1| 宅男免费午夜| 日本精品一区二区三区蜜桃| 日本黄色日本黄色录像| 久9热在线精品视频| 丝袜美足系列| 亚洲精品久久成人aⅴ小说| 国产av又大| 日韩免费高清中文字幕av| 久久精品国产亚洲av香蕉五月 | www.自偷自拍.com| 99riav亚洲国产免费| 热re99久久国产66热| 精品视频人人做人人爽| 亚洲人成伊人成综合网2020| 动漫黄色视频在线观看| 亚洲国产精品一区二区三区在线| 免费在线观看日本一区| 日韩欧美一区视频在线观看| 精品无人区乱码1区二区| a级毛片黄视频| 亚洲国产欧美一区二区综合| 热99re8久久精品国产| 国产激情久久老熟女| 欧美日韩国产mv在线观看视频| 一边摸一边做爽爽视频免费| tube8黄色片| 99在线人妻在线中文字幕 | 中文亚洲av片在线观看爽 | 黄色毛片三级朝国网站| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区三区综合在线观看| 很黄的视频免费| 国产亚洲精品一区二区www | bbb黄色大片| 91大片在线观看| 欧美人与性动交α欧美软件| 成人国产一区最新在线观看| 人妻一区二区av| 亚洲精品自拍成人| 久久热在线av| 在线观看免费视频网站a站| 一区在线观看完整版| 黑人操中国人逼视频| 每晚都被弄得嗷嗷叫到高潮| 啦啦啦免费观看视频1| 日韩熟女老妇一区二区性免费视频| 欧美日韩成人在线一区二区| www.熟女人妻精品国产| 好看av亚洲va欧美ⅴa在| 不卡一级毛片| 久久精品国产清高在天天线| 亚洲aⅴ乱码一区二区在线播放 | 动漫黄色视频在线观看| 国产精品 国内视频| 日日摸夜夜添夜夜添小说| 午夜福利乱码中文字幕| 女人被狂操c到高潮| 不卡一级毛片| 丝袜人妻中文字幕| 99国产精品一区二区三区| 99re6热这里在线精品视频| 亚洲免费av在线视频| 久久久精品区二区三区| 老鸭窝网址在线观看| 成人精品一区二区免费| 亚洲久久久国产精品| 老司机福利观看| 每晚都被弄得嗷嗷叫到高潮| 在线观看免费日韩欧美大片| av免费在线观看网站| av有码第一页| 精品久久久久久久久久免费视频 | 久久亚洲精品不卡| 久久久国产成人精品二区 | 日本黄色日本黄色录像| 久久ye,这里只有精品| 午夜影院日韩av| 又大又爽又粗| 国产成+人综合+亚洲专区| 亚洲熟妇中文字幕五十中出 | 成人永久免费在线观看视频| 美女高潮到喷水免费观看| 久久久国产一区二区| 亚洲全国av大片| 麻豆av在线久日| 91九色精品人成在线观看| 人妻一区二区av| 国产亚洲欧美98| 欧美成人午夜精品| 国产成人精品久久二区二区91| 在线观看免费高清a一片| 一进一出好大好爽视频| 欧美成人免费av一区二区三区 | 一边摸一边抽搐一进一小说 | 午夜福利在线观看吧| 午夜福利影视在线免费观看| 欧美日本中文国产一区发布| 欧美久久黑人一区二区| 一区二区日韩欧美中文字幕| 色综合欧美亚洲国产小说| 久久久久久久国产电影| 1024视频免费在线观看| 国产精品一区二区精品视频观看| 欧美日韩乱码在线| 免费一级毛片在线播放高清视频 | 村上凉子中文字幕在线| 国产国语露脸激情在线看| 麻豆乱淫一区二区| 久久精品91无色码中文字幕| 亚洲一区中文字幕在线| 精品免费久久久久久久清纯 | 久久99一区二区三区| 亚洲成人手机| 久久久久久久午夜电影 | 欧美日韩黄片免| 免费高清在线观看日韩| 国产区一区二久久| 中文字幕av电影在线播放| 在线观看66精品国产| 久热这里只有精品99| 满18在线观看网站| 男女高潮啪啪啪动态图| 日日摸夜夜添夜夜添小说| 亚洲国产欧美日韩在线播放| 啦啦啦免费观看视频1| 国产一卡二卡三卡精品| 久久精品91无色码中文字幕| 黄色丝袜av网址大全| 国产视频一区二区在线看| 国产主播在线观看一区二区| 午夜福利在线观看吧| 人人澡人人妻人| 亚洲国产看品久久| 亚洲第一欧美日韩一区二区三区| 午夜福利视频在线观看免费| 最新在线观看一区二区三区| 男女午夜视频在线观看| 亚洲免费av在线视频| 大型黄色视频在线免费观看| 在线观看日韩欧美| 久久婷婷成人综合色麻豆| 超色免费av| 老熟妇乱子伦视频在线观看| av网站免费在线观看视频| 久久精品亚洲精品国产色婷小说| 欧美激情高清一区二区三区| 咕卡用的链子| 国产男女超爽视频在线观看| 亚洲欧美一区二区三区久久| 成在线人永久免费视频| 精品久久久精品久久久| 久久中文字幕人妻熟女| 女人爽到高潮嗷嗷叫在线视频| 成人免费观看视频高清| 热re99久久国产66热| 欧美激情极品国产一区二区三区| 中文字幕av电影在线播放| 91精品三级在线观看| 91大片在线观看| 国产蜜桃级精品一区二区三区 | 亚洲片人在线观看| 欧美成人免费av一区二区三区 | 老司机靠b影院| 欧美在线一区亚洲| 美女国产高潮福利片在线看| 亚洲国产精品一区二区三区在线| 一边摸一边抽搐一进一出视频| 亚洲性夜色夜夜综合| 国内久久婷婷六月综合欲色啪| 777久久人妻少妇嫩草av网站| 91字幕亚洲| 免费观看人在逋| 亚洲情色 制服丝袜| 九色亚洲精品在线播放| 国产不卡av网站在线观看| 桃红色精品国产亚洲av| 国产成人精品久久二区二区免费| 日本欧美视频一区| 久久草成人影院| 男女午夜视频在线观看| 久久精品国产99精品国产亚洲性色 | 夜夜夜夜夜久久久久| 黑人巨大精品欧美一区二区蜜桃| 老司机午夜十八禁免费视频| 色在线成人网| 少妇的丰满在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 99久久国产精品久久久| 欧美日韩精品网址| 天天影视国产精品| 天天添夜夜摸| 久久人人爽av亚洲精品天堂| 久久狼人影院| 亚洲av熟女| 亚洲精品av麻豆狂野| 久久久久久久国产电影| 天堂中文最新版在线下载| 热re99久久国产66热| 色94色欧美一区二区| 亚洲成人手机| 免费观看人在逋| 高清欧美精品videossex| 欧美日韩乱码在线| 一边摸一边抽搐一进一出视频| 满18在线观看网站| 一二三四在线观看免费中文在| 色精品久久人妻99蜜桃| 国产精品亚洲一级av第二区| 在线观看午夜福利视频| 亚洲三区欧美一区| 精品福利永久在线观看| 高潮久久久久久久久久久不卡| 9热在线视频观看99| 日韩 欧美 亚洲 中文字幕| 久久人妻熟女aⅴ| 久久人妻av系列| 一级黄色大片毛片| 男女床上黄色一级片免费看| 如日韩欧美国产精品一区二区三区| 亚洲欧美一区二区三区久久| 80岁老熟妇乱子伦牲交| 久久久国产成人免费| 超碰成人久久| 亚洲精品久久成人aⅴ小说| 女警被强在线播放| 女人被躁到高潮嗷嗷叫费观| 无人区码免费观看不卡| 波多野结衣av一区二区av| 久久久久国产一级毛片高清牌| 国产真人三级小视频在线观看| 成人三级做爰电影| 少妇被粗大的猛进出69影院| 精品久久久久久久久久免费视频 | 黄色视频不卡| 黑人操中国人逼视频| 深夜精品福利| 十分钟在线观看高清视频www| 欧美 亚洲 国产 日韩一| 国产av一区二区精品久久| 免费在线观看黄色视频的| 国产精品99久久99久久久不卡| 五月开心婷婷网| 久久影院123| 久久国产精品男人的天堂亚洲| 99香蕉大伊视频| www日本在线高清视频| 国产精品自产拍在线观看55亚洲 | 交换朋友夫妻互换小说| 女同久久另类99精品国产91| 波多野结衣一区麻豆| 亚洲专区中文字幕在线| 亚洲精品乱久久久久久| av视频免费观看在线观看| 国产片内射在线| 亚洲成人国产一区在线观看| 日日爽夜夜爽网站| 中文字幕另类日韩欧美亚洲嫩草| 美女国产高潮福利片在线看| 久久精品亚洲av国产电影网| 大型黄色视频在线免费观看| 中文字幕色久视频| 国产av一区二区精品久久| 黑人欧美特级aaaaaa片| 久久久久久久午夜电影 | 亚洲欧美色中文字幕在线| 美女扒开内裤让男人捅视频| 国产一卡二卡三卡精品| 搡老熟女国产l中国老女人| 成年人免费黄色播放视频| 精品高清国产在线一区| 狠狠狠狠99中文字幕| 一级a爱视频在线免费观看| 中文欧美无线码| 18禁观看日本| 欧美一级毛片孕妇| 宅男免费午夜| 热re99久久国产66热| 午夜亚洲福利在线播放| 亚洲欧美日韩另类电影网站| 麻豆国产av国片精品| 国产在线一区二区三区精| 久久久久国产一级毛片高清牌| 久久精品国产亚洲av香蕉五月 | 女性生殖器流出的白浆| 最近最新中文字幕大全电影3 | 久久人人97超碰香蕉20202| 99riav亚洲国产免费| 成人永久免费在线观看视频| 久久久水蜜桃国产精品网| 高清在线国产一区| 亚洲欧美色中文字幕在线| 日本撒尿小便嘘嘘汇集6| 18禁观看日本| 久久久久国产一级毛片高清牌| 黄色女人牲交| 丁香六月欧美| 在线看a的网站| 一级黄色大片毛片| 国产精品美女特级片免费视频播放器 | 91麻豆精品激情在线观看国产 | 国产成人免费观看mmmm| www.自偷自拍.com| 别揉我奶头~嗯~啊~动态视频| 91在线观看av| 亚洲精品美女久久av网站| 亚洲 欧美一区二区三区| 五月开心婷婷网| 国产99久久九九免费精品| 国产91精品成人一区二区三区| 两人在一起打扑克的视频| 变态另类成人亚洲欧美熟女 | 中文亚洲av片在线观看爽 | 一级毛片高清免费大全| 国产一区二区激情短视频| 久久精品亚洲精品国产色婷小说| 丝袜人妻中文字幕| 999久久久国产精品视频| x7x7x7水蜜桃| 极品少妇高潮喷水抽搐| 老司机午夜十八禁免费视频| 国产精品亚洲一级av第二区| ponron亚洲| tocl精华| 久久久精品区二区三区| 91国产中文字幕| 久久久久久免费高清国产稀缺| 91国产中文字幕| 建设人人有责人人尽责人人享有的| 视频在线观看一区二区三区| 亚洲欧美日韩另类电影网站| 国产精品久久久久久人妻精品电影| 国产1区2区3区精品| 国产亚洲欧美在线一区二区| 丰满迷人的少妇在线观看| 久久狼人影院| 久久久国产成人精品二区 | 欧美亚洲日本最大视频资源| 亚洲人成电影观看| 亚洲欧美日韩另类电影网站| 午夜福利在线观看吧| a级毛片黄视频| 国产精品久久久久久人妻精品电影| 免费日韩欧美在线观看| 狠狠婷婷综合久久久久久88av| 99久久99久久久精品蜜桃| 老司机影院毛片| 在线观看免费高清a一片| 大香蕉久久网| 亚洲精品久久午夜乱码| 国产亚洲欧美在线一区二区| 亚洲色图综合在线观看| 高清视频免费观看一区二区| 亚洲精品久久午夜乱码| 亚洲av成人av| 我的亚洲天堂| 亚洲熟妇中文字幕五十中出 | 日日摸夜夜添夜夜添小说| 久久久国产精品麻豆| 国产1区2区3区精品| 久久久久久久午夜电影 | 曰老女人黄片| 欧美在线一区亚洲| 国产精品国产高清国产av | 亚洲午夜精品一区,二区,三区| 亚洲黑人精品在线| 九色亚洲精品在线播放| 免费观看人在逋| 咕卡用的链子| 香蕉丝袜av| 国产精品久久久久成人av| 天堂中文最新版在线下载| 国产免费男女视频| 亚洲av第一区精品v没综合| 成年人午夜在线观看视频| 久久精品国产亚洲av高清一级| 欧美日韩av久久| 1024视频免费在线观看| 一二三四在线观看免费中文在| 最近最新中文字幕大全免费视频| 午夜激情av网站| 丝袜人妻中文字幕| 悠悠久久av| 精品久久久久久电影网| 久久国产精品男人的天堂亚洲| 亚洲精品国产一区二区精华液| 久久天堂一区二区三区四区| 久久人妻av系列| 成年人免费黄色播放视频| 看片在线看免费视频| 亚洲av日韩精品久久久久久密| 无限看片的www在线观看| 久久久久久久久久久久大奶| 国产亚洲一区二区精品| 亚洲熟女精品中文字幕| 欧美+亚洲+日韩+国产| 精品久久久久久电影网| 在线视频色国产色| 亚洲国产中文字幕在线视频| 一本一本久久a久久精品综合妖精| 高潮久久久久久久久久久不卡| 男女床上黄色一级片免费看| 黄片播放在线免费| 亚洲第一青青草原| xxx96com| 久久久水蜜桃国产精品网| 99精品在免费线老司机午夜| 亚洲精品美女久久av网站| 欧美一级毛片孕妇| 巨乳人妻的诱惑在线观看| 可以免费在线观看a视频的电影网站| 99久久人妻综合| 成人国语在线视频| 巨乳人妻的诱惑在线观看| www.精华液| a级片在线免费高清观看视频| 日日爽夜夜爽网站| 欧美激情高清一区二区三区| 欧美精品亚洲一区二区| 国产精华一区二区三区| www.自偷自拍.com| 精品欧美一区二区三区在线| 热99久久久久精品小说推荐| 人妻丰满熟妇av一区二区三区 | 精品国产一区二区久久| 免费看a级黄色片| 久久久国产成人精品二区 | 在线观看免费视频网站a站| 精品第一国产精品| 女性生殖器流出的白浆| 热99国产精品久久久久久7| 国产黄色免费在线视频| 每晚都被弄得嗷嗷叫到高潮| 精品无人区乱码1区二区| 国产一区有黄有色的免费视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲成av片中文字幕在线观看| 亚洲精品av麻豆狂野| 亚洲性夜色夜夜综合| 欧美日韩av久久| 如日韩欧美国产精品一区二区三区| 丝瓜视频免费看黄片| 亚洲第一欧美日韩一区二区三区| 他把我摸到了高潮在线观看| 亚洲欧美日韩高清在线视频| 99国产综合亚洲精品| 97人妻天天添夜夜摸| 在线观看免费日韩欧美大片| 午夜精品国产一区二区电影|