• <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)的振顫干擾
    亚洲中文字幕日韩| 人妻少妇偷人精品九色| 久久国产乱子免费精品| 国产精品亚洲一级av第二区| 亚洲国产精品久久男人天堂| 国产精品女同一区二区软件| 日本成人三级电影网站| 国产单亲对白刺激| 久久久精品94久久精品| 老熟妇仑乱视频hdxx| 99精品在免费线老司机午夜| 狂野欧美白嫩少妇大欣赏| 国产精品久久视频播放| 午夜福利在线观看吧| 欧美在线一区亚洲| 中文字幕精品亚洲无线码一区| 日韩欧美国产在线观看| 免费av不卡在线播放| 国产男靠女视频免费网站| 国产一区亚洲一区在线观看| 一本久久中文字幕| 又爽又黄a免费视频| 搞女人的毛片| 日韩欧美精品v在线| 亚洲熟妇熟女久久| 男插女下体视频免费在线播放| 成年女人毛片免费观看观看9| 成人欧美大片| 亚洲成人久久性| 91av网一区二区| 国内揄拍国产精品人妻在线| 欧美一区二区精品小视频在线| 久久亚洲精品不卡| 91久久精品国产一区二区成人| 狂野欧美白嫩少妇大欣赏| 人妻久久中文字幕网| 18+在线观看网站| 寂寞人妻少妇视频99o| 国内揄拍国产精品人妻在线| 亚洲精品影视一区二区三区av| 久久精品夜色国产| 亚洲人成网站在线观看播放| 国产女主播在线喷水免费视频网站 | 国语自产精品视频在线第100页| 亚洲精品在线观看二区| 欧美一区二区国产精品久久精品| 春色校园在线视频观看| 久久人妻av系列| 国产精品一区二区三区四区免费观看 | 日本a在线网址| 秋霞在线观看毛片| 国产精品野战在线观看| 91午夜精品亚洲一区二区三区| 日本精品一区二区三区蜜桃| 老司机影院成人| 久久精品国产亚洲av天美| 免费看光身美女| 婷婷精品国产亚洲av在线| 国产av不卡久久| 18禁裸乳无遮挡免费网站照片| 国产91av在线免费观看| 国产精品不卡视频一区二区| 亚洲自拍偷在线| 欧美一区二区国产精品久久精品| 亚洲人与动物交配视频| 欧美日韩国产亚洲二区| 亚洲精品在线观看二区| 免费大片18禁| 特大巨黑吊av在线直播| 免费无遮挡裸体视频| 深夜a级毛片| 综合色av麻豆| 性欧美人与动物交配| 亚洲精华国产精华液的使用体验 | 国产高清视频在线观看网站| 欧美一区二区亚洲| 日韩高清综合在线| 国产日本99.免费观看| 99热这里只有是精品在线观看| 亚洲成a人片在线一区二区| 国产成人aa在线观看| 国产精品女同一区二区软件| 欧美日韩精品成人综合77777| 大香蕉久久网| 国产精品伦人一区二区| 丝袜美腿在线中文| 成人国产麻豆网| 国产黄片美女视频| 精品一区二区三区av网在线观看| 在线免费观看不下载黄p国产| 一级毛片我不卡| 一夜夜www| 在线国产一区二区在线| 国产高清三级在线| 欧美最黄视频在线播放免费| 日日摸夜夜添夜夜添小说| 精品国内亚洲2022精品成人| 欧美日本视频| 国语自产精品视频在线第100页| 亚洲国产高清在线一区二区三| 91久久精品电影网| 亚洲乱码一区二区免费版| 日本 av在线| 欧美精品国产亚洲| 最近视频中文字幕2019在线8| 一个人观看的视频www高清免费观看| 久久国内精品自在自线图片| 性欧美人与动物交配| 国产色婷婷99| 日本五十路高清| 最近视频中文字幕2019在线8| 麻豆av噜噜一区二区三区| 伦精品一区二区三区| 亚洲成人中文字幕在线播放| 一级毛片我不卡| 极品教师在线视频| 国产色婷婷99| 欧美日韩国产亚洲二区| 一个人看的www免费观看视频| 久久欧美精品欧美久久欧美| 99久久精品热视频| 亚洲国产精品成人综合色| 韩国av在线不卡| 亚洲人成网站高清观看| 国产亚洲精品综合一区在线观看| 人妻制服诱惑在线中文字幕| 精品不卡国产一区二区三区| 日韩一本色道免费dvd| 亚洲美女黄片视频| 国产成人福利小说| 日本黄色片子视频| h日本视频在线播放| 亚洲自拍偷在线| 99久久精品热视频| 女生性感内裤真人,穿戴方法视频| 性色avwww在线观看| 女的被弄到高潮叫床怎么办| 日本熟妇午夜| 久久久久国产网址| 亚洲乱码一区二区免费版| 国产一区二区三区在线臀色熟女| 99在线视频只有这里精品首页| 国产大屁股一区二区在线视频| 免费看光身美女| 欧美+亚洲+日韩+国产| 亚洲av免费在线观看| 欧美成人a在线观看| 韩国av在线不卡| 丝袜喷水一区| 国产探花极品一区二区| 国产免费一级a男人的天堂| 亚州av有码| 少妇的逼好多水| 午夜福利视频1000在线观看| 欧美另类亚洲清纯唯美| 99久久精品一区二区三区| 日本-黄色视频高清免费观看| 久久99热这里只有精品18| 嫩草影院新地址| 久久久国产成人精品二区| 熟妇人妻久久中文字幕3abv| 欧美日韩精品成人综合77777| av视频在线观看入口| 成人亚洲精品av一区二区| 亚洲国产精品合色在线| 精品久久久久久久久亚洲| 国产黄色小视频在线观看| 一本一本综合久久| 日本在线视频免费播放| 网址你懂的国产日韩在线| 男人舔女人下体高潮全视频| 日日摸夜夜添夜夜爱| 内地一区二区视频在线| 欧美人与善性xxx| 国产高潮美女av| 天堂√8在线中文| 久久精品夜色国产| 久久久久久久久中文| 亚洲中文字幕一区二区三区有码在线看| 日韩三级伦理在线观看| 97热精品久久久久久| 直男gayav资源| 插逼视频在线观看| 日韩欧美一区二区三区在线观看| 免费搜索国产男女视频| 99久国产av精品| 国产免费男女视频| 人人妻,人人澡人人爽秒播| 亚洲一级一片aⅴ在线观看| 91在线精品国自产拍蜜月| 国产久久久一区二区三区| 亚洲欧美日韩东京热| 日本色播在线视频| 欧美激情在线99| 久久午夜亚洲精品久久| 女人被狂操c到高潮| 国产成人freesex在线 | 乱人视频在线观看| 日韩高清综合在线| 精品久久久久久成人av| 波多野结衣巨乳人妻| 黄色配什么色好看| 97碰自拍视频| 黑人高潮一二区| 丝袜美腿在线中文| 久久久久性生活片| 亚洲人成网站高清观看| 老熟妇乱子伦视频在线观看| 国产爱豆传媒在线观看| 国产精品一及| 少妇猛男粗大的猛烈进出视频 | 五月玫瑰六月丁香| 国产精品爽爽va在线观看网站| 亚洲精品国产av成人精品 | 免费看a级黄色片| 国产亚洲91精品色在线| 搡老妇女老女人老熟妇| 成人漫画全彩无遮挡| 午夜老司机福利剧场| 国内久久婷婷六月综合欲色啪| 久久久久久久久中文| 日韩欧美 国产精品| 久久久成人免费电影| 日日干狠狠操夜夜爽| 国产精华一区二区三区| 中出人妻视频一区二区| 99热这里只有是精品50| 亚洲av二区三区四区| 久久人人精品亚洲av| 亚洲内射少妇av| 亚洲精品影视一区二区三区av| 女的被弄到高潮叫床怎么办| 欧美日本视频| 亚洲精品在线观看二区| 长腿黑丝高跟| 男人和女人高潮做爰伦理| 亚洲第一电影网av| 看十八女毛片水多多多| 国产不卡一卡二| 俄罗斯特黄特色一大片| 日本一二三区视频观看| 国产aⅴ精品一区二区三区波| 国产精品久久久久久久电影| 国内揄拍国产精品人妻在线| 男人和女人高潮做爰伦理| 尾随美女入室| 九九久久精品国产亚洲av麻豆| 午夜爱爱视频在线播放| 天天一区二区日本电影三级| 久久中文看片网| 国产亚洲精品久久久com| 国产大屁股一区二区在线视频| 人妻制服诱惑在线中文字幕| 婷婷精品国产亚洲av在线| a级毛片免费高清观看在线播放| 最后的刺客免费高清国语| 日韩欧美国产在线观看| 在线观看免费视频日本深夜| 国产私拍福利视频在线观看| 一本久久中文字幕| 亚洲人成网站在线播| 日韩在线高清观看一区二区三区| 婷婷精品国产亚洲av在线| 91在线精品国自产拍蜜月| 午夜福利成人在线免费观看| 国产探花极品一区二区| 久久婷婷人人爽人人干人人爱| 六月丁香七月| 成熟少妇高潮喷水视频| 网址你懂的国产日韩在线| 国产私拍福利视频在线观看| 日本一二三区视频观看| 国产精品野战在线观看| 国产午夜精品久久久久久一区二区三区 | 国产精品久久久久久精品电影| 久久天躁狠狠躁夜夜2o2o| 国产高清三级在线| 亚洲中文字幕一区二区三区有码在线看| 波多野结衣高清无吗| 久久久精品94久久精品| 国产欧美日韩一区二区精品| 精品人妻一区二区三区麻豆 | 国产成人91sexporn| 午夜免费激情av| 午夜久久久久精精品| 国内精品美女久久久久久| 99久久精品热视频| 2021天堂中文幕一二区在线观| 亚洲熟妇中文字幕五十中出| 能在线免费观看的黄片| 国产精品亚洲美女久久久| 3wmmmm亚洲av在线观看| 成人午夜高清在线视频| 欧美3d第一页| 全区人妻精品视频| 3wmmmm亚洲av在线观看| 舔av片在线| 99久国产av精品| 精品久久久噜噜| 老师上课跳d突然被开到最大视频| 免费av不卡在线播放| 亚洲美女视频黄频| 亚洲精品国产av成人精品 | 久久久久国产精品人妻aⅴ院| 亚洲欧美日韩无卡精品| 日本a在线网址| 午夜精品国产一区二区电影 | 国产精品久久视频播放| 男女那种视频在线观看| 免费在线观看影片大全网站| 国产一区二区激情短视频| 99热6这里只有精品| 尾随美女入室| 男女下面进入的视频免费午夜| 一区二区三区免费毛片| 成人美女网站在线观看视频| 亚洲av中文av极速乱| 一个人免费在线观看电影| 如何舔出高潮| 插逼视频在线观看| 综合色av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 国产午夜福利久久久久久| 97在线视频观看| 久久精品国产亚洲网站| 日本成人三级电影网站| 国产一区二区亚洲精品在线观看| 国产在视频线在精品| 欧美日韩精品成人综合77777| 日本撒尿小便嘘嘘汇集6| 国产精品久久久久久久电影| 欧美性猛交╳xxx乱大交人| 欧美成人精品欧美一级黄| 久久鲁丝午夜福利片| 久久精品国产自在天天线| 亚洲欧美清纯卡通| 亚洲天堂国产精品一区在线| 国产高清不卡午夜福利| 久久久久久久久大av| 亚洲四区av| a级一级毛片免费在线观看| 亚洲美女黄片视频| 久久韩国三级中文字幕| 国产精品一区二区免费欧美| 精品人妻偷拍中文字幕| 一级毛片久久久久久久久女| 久久精品国产亚洲av天美| 三级国产精品欧美在线观看| 淫秽高清视频在线观看| 男人狂女人下面高潮的视频| 国产亚洲精品久久久久久毛片| 久久精品影院6| 色视频www国产| 黄色日韩在线| 国产欧美日韩精品亚洲av| 性欧美人与动物交配| 精品熟女少妇av免费看| 国产精品女同一区二区软件| 中文字幕人妻熟人妻熟丝袜美| 色在线成人网| 午夜福利在线观看吧| 亚洲中文日韩欧美视频| 亚洲精品在线观看二区| 免费人成视频x8x8入口观看| 在线播放国产精品三级| 两个人的视频大全免费| 精品久久久久久久久av| 男人狂女人下面高潮的视频| 国产白丝娇喘喷水9色精品| 国产成年人精品一区二区| 99热全是精品| 波多野结衣高清无吗| 久久精品夜色国产| 日产精品乱码卡一卡2卡三| 精品无人区乱码1区二区| 嫩草影视91久久| 久久九九热精品免费| 欧美绝顶高潮抽搐喷水| 3wmmmm亚洲av在线观看| 最近手机中文字幕大全| 国产伦精品一区二区三区四那| 亚洲自拍偷在线| 99久久无色码亚洲精品果冻| 亚洲天堂国产精品一区在线| 国产精品综合久久久久久久免费| 嫩草影院新地址| 人人妻人人看人人澡| 亚洲自偷自拍三级| 晚上一个人看的免费电影| 欧美三级亚洲精品| 亚洲最大成人av| 欧美成人精品欧美一级黄| 你懂的网址亚洲精品在线观看 | 成人高潮视频无遮挡免费网站| 亚洲色图av天堂| 国产精品人妻久久久影院| 极品教师在线视频| 国产高清视频在线播放一区| 黄色视频,在线免费观看| av在线老鸭窝| 日韩一本色道免费dvd| 69人妻影院| 一本精品99久久精品77| 老熟妇仑乱视频hdxx| 亚洲国产精品成人综合色| 国产伦在线观看视频一区| 一级毛片久久久久久久久女| 日韩精品中文字幕看吧| 97热精品久久久久久| 国产熟女欧美一区二区| 午夜精品在线福利| 中文字幕av成人在线电影| 国产麻豆成人av免费视频| 大型黄色视频在线免费观看| 在线看三级毛片| 欧美日韩在线观看h| 99久国产av精品| 一区福利在线观看| 在线播放国产精品三级| 波多野结衣高清作品| 精品国内亚洲2022精品成人| 久久久精品欧美日韩精品| 成人永久免费在线观看视频| 夜夜爽天天搞| 午夜福利在线观看免费完整高清在 | 亚洲欧美清纯卡通| 精品人妻偷拍中文字幕| 成人鲁丝片一二三区免费| 亚洲成人中文字幕在线播放| 久久久久精品国产欧美久久久| 亚洲一级一片aⅴ在线观看| 免费观看精品视频网站| 中文字幕熟女人妻在线| 日韩国内少妇激情av| 成人午夜高清在线视频| 国产熟女欧美一区二区| 国产人妻一区二区三区在| 97碰自拍视频| 亚洲欧美精品综合久久99| 99久久精品热视频| 日韩大尺度精品在线看网址| 淫秽高清视频在线观看| 久久久久国产网址| 午夜老司机福利剧场| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧美中文字幕日韩二区| 人人妻人人澡欧美一区二区| 免费黄网站久久成人精品| 99久久精品热视频| 联通29元200g的流量卡| 久久综合国产亚洲精品| 国产精品1区2区在线观看.| 18禁裸乳无遮挡免费网站照片| 身体一侧抽搐| 在线免费观看不下载黄p国产| 日韩欧美精品v在线| 内地一区二区视频在线| 亚洲国产欧洲综合997久久,| 97超碰精品成人国产| 国产精品精品国产色婷婷| 免费av观看视频| 国产精品女同一区二区软件| 两性午夜刺激爽爽歪歪视频在线观看| 卡戴珊不雅视频在线播放| 在线观看66精品国产| 日韩亚洲欧美综合| 国产高清视频在线播放一区| 99久国产av精品国产电影| 亚洲av二区三区四区| 嫩草影院精品99| 久久草成人影院| 国产高清不卡午夜福利| 国产男人的电影天堂91| 乱系列少妇在线播放| 啦啦啦观看免费观看视频高清| 九九在线视频观看精品| 国产成年人精品一区二区| 久久九九热精品免费| 啦啦啦韩国在线观看视频| 免费在线观看影片大全网站| 亚洲国产欧美人成| 51国产日韩欧美| 亚洲成av人片在线播放无| 国产一区二区三区av在线 | 男女做爰动态图高潮gif福利片| 97热精品久久久久久| 看黄色毛片网站| 丰满的人妻完整版| 天堂影院成人在线观看| 99久久中文字幕三级久久日本| 精品一区二区三区av网在线观看| 美女xxoo啪啪120秒动态图| 欧美最新免费一区二区三区| 一进一出抽搐动态| 最新中文字幕久久久久| 一级毛片aaaaaa免费看小| 日本爱情动作片www.在线观看 | 少妇熟女欧美另类| 亚洲欧美成人综合另类久久久 | 中文在线观看免费www的网站| 97在线视频观看| 婷婷亚洲欧美| 麻豆成人午夜福利视频| 国产一区二区三区av在线 | 国产一级毛片七仙女欲春2| 久久人人爽人人片av| 男女做爰动态图高潮gif福利片| 可以在线观看毛片的网站| av在线播放精品| 亚洲中文字幕一区二区三区有码在线看| 亚洲成a人片在线一区二区| 日韩欧美一区二区三区在线观看| 亚洲在线自拍视频| 精品午夜福利在线看| av专区在线播放| 中文字幕熟女人妻在线| 久久久久久久久久成人| 91在线精品国自产拍蜜月| 大香蕉久久网| 美女内射精品一级片tv| 国内精品一区二区在线观看| 精品免费久久久久久久清纯| 欧美绝顶高潮抽搐喷水| 深爱激情五月婷婷| 亚洲丝袜综合中文字幕| 长腿黑丝高跟| 亚洲第一电影网av| 日日撸夜夜添| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一及| 国产在线精品亚洲第一网站| 亚州av有码| 国产大屁股一区二区在线视频| 亚洲一级一片aⅴ在线观看| 久久久久久久亚洲中文字幕| 哪里可以看免费的av片| 久久久a久久爽久久v久久| 中文字幕av在线有码专区| 国产精品一区二区性色av| 最近的中文字幕免费完整| 亚洲自拍偷在线| 欧美成人a在线观看| 亚洲欧美精品自产自拍| 嫩草影院精品99| 久久鲁丝午夜福利片| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲av涩爱 | 亚洲av成人av| 成人一区二区视频在线观看| 91久久精品国产一区二区三区| 伦精品一区二区三区| 丝袜美腿在线中文| 欧美区成人在线视频| 好男人在线观看高清免费视频| 此物有八面人人有两片| 美女内射精品一级片tv| 两个人的视频大全免费| 国产69精品久久久久777片| 国产伦精品一区二区三区视频9| 国产大屁股一区二区在线视频| 国产精品一及| 欧美一级a爱片免费观看看| 麻豆国产97在线/欧美| 波多野结衣巨乳人妻| 日韩成人av中文字幕在线观看 | 一夜夜www| 精品一区二区免费观看| 丝袜喷水一区| 亚洲精品粉嫩美女一区| 小蜜桃在线观看免费完整版高清| 欧美性猛交╳xxx乱大交人| 国产精品99久久久久久久久| 国产精品久久视频播放| 久久人人精品亚洲av| 午夜精品一区二区三区免费看| 99国产精品一区二区蜜桃av| 亚洲精品日韩在线中文字幕 | 亚洲第一区二区三区不卡| 蜜桃亚洲精品一区二区三区| 国产一区二区激情短视频| 国产高清不卡午夜福利| 欧美中文日本在线观看视频| АⅤ资源中文在线天堂| 婷婷精品国产亚洲av| 夜夜夜夜夜久久久久| 亚洲18禁久久av| 白带黄色成豆腐渣| 在线看三级毛片| 国产精品一区二区三区四区久久| 免费无遮挡裸体视频| 亚洲成a人片在线一区二区| 国产男人的电影天堂91| 网址你懂的国产日韩在线| 97超级碰碰碰精品色视频在线观看| 亚洲18禁久久av| 网址你懂的国产日韩在线| 麻豆精品久久久久久蜜桃| 国产精品国产三级国产av玫瑰| 亚洲欧美日韩东京热| 中文字幕精品亚洲无线码一区| 免费不卡的大黄色大毛片视频在线观看 | 一区二区三区免费毛片| 成人欧美大片| 国产高清视频在线观看网站| 看黄色毛片网站| 欧美成人a在线观看| 悠悠久久av| a级毛色黄片| 九九久久精品国产亚洲av麻豆| 亚洲人成网站高清观看| 麻豆国产97在线/欧美| 久久久久免费精品人妻一区二区| 国产精品久久久久久av不卡| 人妻夜夜爽99麻豆av| 性欧美人与动物交配|