999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

一類非自治隨機積微分方程的最優控制問題

2025-07-14 00:00:00劉夢杰范虹霞
哈爾濱理工大學學報 2025年2期
關鍵詞:定義系統

中圖分類號:0175.6 文獻標志碼:A 文章編號:1007-2683(2025)02-0150-09

Abstract:Inordertostudytheoptimalcontrolproblemofaclassofnon-autonomousstochasticintegro-diferentialequationin Hilertspace,theexistencetheoremofmildsolutionoftheproblemisestablishedbyusingstochasticanalysis,theoriesofoperator semigroupsandSadovskifixedpointtheorem.Ontisbasis,theexistenceoftheoptialstate-controlpairisdiscussedbyostucting minimization sequences twice.

Keywords:stochasticnalysis;teoryofoperatoremigroups;ntego-diferetialequatin;Sadovskifixedpointtheorem;tial state-control pair

0 引言

最優控制理論是數學控制理論中的一個重要概念,在控制系統中起著關鍵作用,它著重于研究使控制系統的性能指標實現最優化的基本條件和綜合方法,在數學、工程、醫學等領域有著廣泛的應用[-4]。而脈沖微分及積微分方程在很多領域中也起著越來越重要的作用。近年來,脈沖微分及積微分方程受到了廣泛研究[5-8]。另外,由于確定性系統通常在實際生活中會受到白噪聲等隨機擾動的影響,因此,對脈沖隨機微分及積微分方程的研究是非常有必要的[9-15] 。

2020年,Rajesh等[2]利用Krasnoselskii不動點定理研究了如下由混合分數布朗運動驅動的二階脈沖隨機微分方程

mild解及最優狀態-控制對的存在性。

由于具有非局部條件的發展方程比經典的柯西問題更具有一般性,在實際物理問題中有更廣泛的應用[16-19]。2021年,Yang 等[4]利用Krasnoselskii 不動點定理研究了下列具有非局部條件的非自治脈沖積微分方程

x(0)+g(x)=x0

mild解及最優狀態-控制對的存在性。

受上述文獻啟發,本文研究如下一類具有非瞬時脈沖的非自治隨機積微分方程

f(t,x(t),Hx(t))dt+

R(t,x(t))dw(t)

mild解及最優狀態-控制對的存在性問題,其中狀態函數 x(?) 在 H 中取值, H 中的內積與范數分別為(?,?) 和 為一稠定閉線性算子族, .A(t) 在 H 上生成一個發展系統 {U(t, (204s):0?s?t?b} ,其中 D(A) 與 Ψt 無關。 H 為一線性有界算子,控制函數 u∈Uad,Uad 為后文定義的可容許控制集 L20(K,H) 為適當定義的函數 (204號x(s))ds,γk 為一非瞬時脈沖函數,設 0=s0lt;…lt; 為一常數。

1 預備知識

本節將介紹一些數學符號并回顧一些必要的基本概念及引理

設 L(K,H) 為所有 K 到 H 的有界線性算子構成的空間。其中 K,H 為實可分Hilbert空間, K,H L(K,H) 的范數均用 ·表示。設 (Ω,Ft,{Ft}t≥0 P )為一完備賦流概率空間,流 {Ft}t?0 為 F 的一列右連續單調遞增的子 σ- 代數族,并且 {F0} 包含所有概率測度為0的集合。設 Q∈L(K,K) 是由 Qei= λiei 定義的算子,且具有有限跡 oi}i≥1 為一非負實數有界序列, {ei}i≥1 是 K 中一組完備標準正交基。設 βi(t) 是完備賦流概率空間中相互獨立的實值標準布朗運動序列,使得

設 δ∈L(K,H) 并且定義如下:

,則稱 δ 為 Q -Hilbert-Schmidt算子。L20(K,H) 為所有 Q -Hilbert-Schmidt算子 所構成的空間。

引理 1[2] (20 對于任意的 c?1 和 L20- 可預測過程 v(?) ,有下面不等式成立

特別地,當 c=1 時,有

定義1 若下列兩個條件成立:

1)U(s,s)=I,U(t,τ)U(τ,s)=U(t,s),

0?s?t?b

則稱雙參有界線性算子族 {U(t,s),0?s?t?b} 為發展系統。

引理 2[4] ( 若 {U(t,s),0?s

設 L2(Ω,H) 為所有強可測、平方可積的 H- 值隨機變量 x 所構成的Banach空間,并賦有范數 為數學期望。

設 C(J,L2(Ω,H)) )為所有從 J 到 L2(Ω,H) 的連續映射構成的Banach空間且滿足條件supE 的子集 L02(Ω,H) 定義如下:

L02(Ω,H)={x∈L2(Ω,H)∣x 是 F0 -可測}。

令 PC(J,H) 為所有 Ft← 適定可測 ?H- 值隨機過程 {x(t):t∈J} 構成的空間,使得 x 在 ?t≠?tk 處連續,在 t=tk 處有 x(tk-)=x(tk)?x(tk+) 存在, k= 1,2,…,m 則 為一Banach空間并賦有范數

算子 B∈L(J,L(G,H)),|B| 。為Banach空間 L )上的范數, G 為一自反可分Hilbert空間,控制函數 u∈LF2(J,G) 在 G 中取值。LF2(J,G) 表示所有 G- 值 -適定可測隨機過程構成的空間,其滿足 且賦有范數

設 U 為 G 中一非空有界凸閉集,集

Uad={u∈LF2(J,G)∣u(t)∈U,t∈J} 為可容許控制集,Bu ∈L2(J,H),u∈Uad。

令 rgt;0 為有限常數,并記

記 S(u):={xu∈Ωr:xu 是隨機系統(2)相應于控制函數 u∈Uad 的mild解}且 :u∈Uad,xu∈S(u)} 。

設 xu∈S(u) 是隨機系統(2)相應于 u∈Uad 的mild解,則最優控制問題可以轉化為有限拉格朗日問題:

找到一個可容許狀態-控制對 u∈Uad ,使得

為積分成本函數。

為隨機系統(2)相應于控制函數 u0 Π∈Uad 的mild解。若 滿足式(4),則 為隨機系統(2)一最優狀態-控制對。

引理 3[4] 若 U(t,s) 為一緊發展系統,則算子

相對緊。此外,當 n?∞ 時,若 un∈Uad 弱收斂于 u ,則 。

下面介紹非緊性測度的定義及基本性質。

設 s 是Banach空間 E?B 中的有界子集,則Hausdorff非緊性測度 β 的定義如下:

可以被半徑小于 ε 的球體覆蓋}。

引理 4[20] (20 設 S,T?E?B 是 E?B 中的有界子集。Hausdorff非緊性測度 β 具有如下性質:

1)β(S)=0?S 在 E?B 中相對緊;

2)S?T?β(S)?β(T)

, 為 s 的凸閉包;4)β(S∪T)?max{β(S),β(T)} :

5)對 ?c∈R,β(cS)=∣c∣β(S)

6)β(S+T)?β(S)+β(T) ,其中

S+T={x∣x=y+z,y∈S,z∈T}

7)如果 是以 k 為常數的Lipschitz連續映射,則對任意有界子集 V?EB ,有 β(Z(V)) ?kβ(V) 。

分別為 中有界集的Hausdorff非緊性測度。此外,對于任意有界集 B?PC(J,H) , t∈[0,b] ,有

引理 5[20] 若 B?E?B 有界, E?B 為Banach空間,則存在可數集 B0?B ,使得 β(B)?2β(B0) 。

引理 6[20] (204號 若 B?C([a1,a2],EB) 有界且等度連續,則 β(B(t)) 在 [a1,a2] 上連續,且有 βc(B) β=maxt∈[a1,a2]β(B(t)) 。

引理 7[20] (20號 設 D={xn}?C([a1,a2],EB) 為 可數集,若存在函數 m∈L1([a1,a2],R+) ,使得對 任意的 n∈N+ ,幾乎處處 ,有 ,則 β(D(t)) 在 [a1,a2] 上是 Lebesgue可積的,并且

定義2 若 Ft -適定可測隨機過程 x∈PC(J, H )滿足以下條件:

1)對幾乎處處的 t∈[0,b],x(t)∈H 有 路徑(左極限存在右連續);

2)對任意 t∈[0,b] ,有

t∈?k=1m(sk,tk+1]

則稱 x(?) 為隨機系統(2)的mild解。

主要結果的證明基于下面的不動點定理

引理 8[4] 設 s 是Banach空間 X 中的一個非空子集,若對任意有界集 D?S ,滿足 α(Q(D) )lt;α(D) ,則連續映射 稱為凝聚的。

引理9[20] (Sadovskii不動點定理)設 X 為一Banach空間, W?X 為一有界凸閉集。若 是凝聚映射,則 Q 在 W 中至少存在一個不動點。

2 存在性結果

為了得到隨機系統(2)mild解的存在性結果,還需如下假設:

(H1) {U(t,s),0?s 滿足

1)函數 h 連續;

2)存在常數 Lh,Dhgt;0 ,使得對于任意的 x,y∈ H ,有

E E

(H3)函數 滿足

1)f(t,?,?) 對于 t∈J 幾乎處處連續且對于 x ,y∈H×H,t 到 f(t,x,y) 是強可測的;

2)存在連續函數 ?:J?R+ ,使得對于幾乎處處 t∈J ,任意 x,y∈H ,有

3)存在函數 lf∈L1(J,R+) ,使得對于任意有界集 A,B?H 且對幾乎處處 t∈J ,有

β(f(t,A,B))?lf(t)(β(A)+β(B))

(H4)函數 滿足

1)g(t,s,?):H?H 對 Φ(t,s)Φ∈Λ 幾乎處處連續且 對 x∈H 強可測;

2)存在連續函數 ,有

E 其中

3)存在函數 η1,η2∈L1(J,R+) ,使得對于任意有界集 C?H 且對幾乎處處 (t,s)∈Λ ,有

β(g(t,s,C))?η1(t)η2(s)β(C) 其中 Λ={(t,s),0?s?t?b} 。

(H5)脈沖函數 滿足

1)函數 γ?k,k=1,2,…,m 連續;

2)存在常數 Lγk,Dγkgt;0 ,使得對于任意 x,y∈ H ,有

其中

(H6)函數 滿足

1)R(t,?) 對于幾乎處處 t∈J 連續且對于所有x∈H,t 到 R(t,x) 是強可測的;

2)存在連續函數 和一個連續非減函數 Dq∈([0,+∞),(0,+∞)) ,使得對于 (Γt,x) ∈J×H ,有

3)存在正常數 L?R ,使得對于任意有界集 D? H ,有

β(R(t,D))?L?Rβ(D)

(H7) 成立。

(H8)對于任意有界集 B1r ,存在一個可數集 {xn}n?1?B1 ,使得

其中 t∈J

定理1 若假設(H1) ~ (H8)成立,且

則隨機系統(2)在 J 上至少有一個mild解。

證明:定義 且(ψx)(t)=(ψ1x)(t)+(ψ2x)(t) ,這里

顯然 x(t) 為隨機系統(2)的mild解當且僅當 x(t) 為算子 ψ 的不動點,下證算子 ψ 存在一個不動點。

第一步 證明存在 rgt;0 ,使得 ,如若不然,對任意 rgt;0 ,存在 ,使得

當 t∈[0,t1] 時,由假設 (H2)~(H4) 及(H6)可知

當 t∈(tk,sk] 時, k=1,2,…,m ,由假設(H5)可知 Dγk(1+r)

當 t∈(sk,tk+1] 時, k=1,2,…,m ,同理,由假設(H2) ~ (H7)可知

綜上,對 t∈[0,b] 有

其中

G* 與 r 無關,對上式兩邊同時除以 r 且當 r?∞ 時求極限,可得

式(6)與式(5)矛盾,因此 O第二步 證明算子 Lipschitz連續。當 t∈[0,t1] 時,對任意 x,y∈Ωr ,由假設(H2)

及 PC(J,H) 空間的定義可知,

E

當 t∈(tk,sk] 時, k=1,2,…,m ,由假設(H5)可知

P

(204號

當 t∈Γ(sk,tk+1] 時 k=1,2,…,m ,由假設(H5)

可知

E |ψ1x(t)-ψ1y(t)|2?MU2LγkE|x(t)-y(t)|2?

綜上,由(H5)可知,對于任意的 x,y∈Ωr

第三步 證明算子 ψ2 在 Ωr 上連續。

設 {xn}n=1∈Ωr 是一個序列且 Ωr ,由 f,g,R 的連續性可知,對任意 s∈J 有

由引理1,假設(H3)及(H7),對于任意 s∈J 有P 4?(s)r(1+bm) E 當 t∈[sk,tk+1] 時, k=0,1,…,m 有

由Lebesgue 控制收斂定理, 0,n∞ ,故 ψ2 在 Ωr 上連續。

第四步 證明算子 等度連續。

對于 x∈Ω,sk''?tk+1,k=0,1,…,m 及 ξgt;0 足夠小有

由假設(H2)~(H4)、(H6)及引理1可知,

當 t'-t'?0,ξ?0 時及由假設(H1)知, I1,I2,I3 →0。

綜上, 0 ,故 等度連續。

第五步 證明 是凝聚算子。

對任意有界集 D∈Ωr ,由引理5可知,存在一個可數集 D0={xn}?D ,使得

βPC(ψ2(D))?2βPC(ψ2(D0))

由于 ψ2(D0)?ψ2(Ωr) 有界且等度連續,由引理6可知

由引理1及引理4可知,有

因此,由引理4及假設(H3)~(H4)、(H6)及(H8)可知,當 t∈[sk,tk+1] 時 k=0,1,…,m ,

對 t∈J ,有

由引理4及第三步可知,對任意有界集 D?Ωr , ,有

為凝聚算子。由Sadovskii不動點定理可得算子 ψ 在 Ωr 上至少有一個不動點,即隨機系統(2)在區間 [0,b] 上至少有一個mild解。

3 最優控制

本節通過構造兩次極小化序列研究相應于隨機系統(2)的有限拉格朗日問題最優控制-狀態對的存在性問題。為了得到隨機系統(2)最優狀態-控制對,有如下假設:

(H9)函數 滿足:

1)函數 1可測;

2)函數 Y(t,?,?),t∈J 在 H×G 上依序列下半連續;

3)函數 Y(t,x,?),t∈J,x∈H 在 G 上為凸的;

4)存在常數 及非負函數 μ∈ L1(J,R+) ,使得

x∈H,u∈G

(H10) {U(t,s),0?s

定理2 若假設(H2)\~(H7)、(H9)及(H10)成立,則相應于隨機系統(2)的有限拉格朗日問題至少存在一個最優狀態-控制對,即存在一個可容許狀態控制對 ,使得

證明:固定 u∈Uad ,定義

?(u):=infxu∈S(u)?(xu,u)

第一步 證明存在 ,使得

不失一般性,假設inf 。由假設xu∈S(u)

(H9)4)可知 為一常數。另一

方面,由下確界定義可知,存在一序列 {xnu}n=1?

S(u) 使得當 n∞ 時, Φ(xnu,u)?Φ(u) 且滿足

由定理1第二步至第五步可知,集 {xnu}n=1 在PC(J,H) 中相對緊。由相對緊的定義可知,存在一子序列,仍記為 且存在 ,使得

此外,由假設(H2)~(H6)及Lebesgue控制收斂定理可知,

因此, 。由于 PC(J,H) 連續嵌入 L1(J,H) 中,由假設(H9)及Balder定理可得

。對于 u∈Uad,? 在 處可以取到最小值。

第二步 證明存在 u0∈Uad ,使得

不失一般性,假設 infu∈Uad?(u)lt;∞ ,由假設(H9)4)可知 。由下確界定義可知,存在一序列 {un}n=1?Uad 使得

由于 {un}n=1 在 LF2(J,G) 中有界,存在一子序列仍記為 {un}n=1 ,有

un?u0∈LF2(J,G),n?∞c

由于 Uad 為一凸閉包,由Marzur引理,可得 u0∈ Uad 。對于任意的 n∈N+ ,由定理1可知,存在 S(un) 滿足 ,則 且滿足

由(H10)及引理3可知, {Bun(?)}n=1 在 PC(J H) )中相對緊。由定理1中第二步及第五步可得 在 PC(J,H) 中相對緊。因此存在一子序列,仍記為 且存在 使得

此外,由假設(H2)~(H6)及Lebesgue控制收斂定理可知

t∈∪k=1m(sk,tk+1]

由上式可得 為可容許-狀態控制對,由于 PC(J,H) 連續嵌入 L1(J,H) 中,由假設

(H9)及Balder定可得

因此, 在 u0∈Uad 處可取得最小值。因此,有

因此,受隨機系統(2)控制的有限拉格朗日問題至少有一個最優狀態-控制對 。證明完畢。

4舉例

考慮如下具有非局部條件的非自治隨機積微分方程

其中 是一致Holder連續, w(t) 是 H 上的一維標準Wiener過程,并定義于一個概率空間 (Ω,F,P) , ψ∈PC(J,R+) 。假設 A(t):H?H 定義如下:

絕對連續, x'∈H,x(0)= (2號 x(π)=0} 。由文[12]可知 {A(t):0?t?b} 生成 一個等度連續族 {U(t,s):0?s?t?b} 且滿足定 義1。設

對于任意 t∈[0,b] ,定義

因此,問題(7)的成本函數可以改寫為系統(2)中成本函數的抽象形式

定理3 若

成立,則問題(7)在區間 [0,b] 上至少一個mild解。

證明:根據 h,f,g,R,γk 的定義可知

易知假設(H1)~(H8)成立,因此定理1中所有條件都成立,其中

則問題(7)在 [0,b] 上至少存在一個mild解且相應的有限拉格朗日問題至少存在一個最優狀態-控制對。

5結語

文[4]在發展系統 U(t,s) 緊的條件下,研究了問題 (1)mild 解的存在性,本文在文[4]的基礎上考慮了隨機效應及非瞬時脈沖的影響,在發展系統U(t,s) 非緊的情形下,利用非緊性測度理論和Sadovskii不動點定理,建立了隨機系統(2)mild解的存在性結果。所獲結論在一定程度上推廣和發展了文[4]的主要結果。

參考文獻:

[1]YAN Zuomao,LU Fangxia.The Optimal Control of a New Class of Impulsive Stochastic Neutral Evolution Integro-differential Equations with Infinite Delay[J]. International Journal of Control,2016,89(8) :1592.

[2]DHAYAL R,MALIK M,ABBAS S,et al. Optimal Controls for Second-order Stochastic Differential Equations Driven by Mixed-fractional Brownian Motion Withimpulses [J].Mathematical Methods in the Applied Sciences, 2020,43(7) :4107.

[3]HUANG Hai,FU Xianlong. Optimal Control Problems for a Semi-linear Integro-differential Evolution System with Infinite Delay[J].Optimal Control Applications and Methods,2021,43(2):459.

[4]YANG He, ZHAO Yanxia. Existence and Optimal Controls of Non-autonomous Impulsive integro-differential Evolution Equation with Nonlocal Conditions[J].Chaos, Solitons and Fractal Chaos,Solitons and Fractals: the Interdisciplinary Journal of Nonlinear Science,and Nonequilibrium and Complex Phenomena,2021,148:1.

[5]REN Yong,HU Lanying, SAKTHIVEL R. Controllability of Impulsive Neutral Stochastic Functional Differential Inclusions with Infinite Delay[J]. Journal of Computational and Applied Mathematics,2010, 235(8):2603.

[6]DENG Sufang,SHU Xiaobao,MAO Jianzhong. Existence and Exponential Stability for Impulsive Neutral Stochastic Functional Diferential Equations Driven by fBm with Noncompact Semigroup Via Monch Fixed Point[J]. Journal of Mathematical Analysisand Applications, 2018,467(1) :398.

[7]ARTHI G,PARK JH, JUNG H Y. Existence and Exponential Stability for Neutral Stochastic Integro-differential Equations with Impulses Driven by a Fractional Brownian Motion[J]. Communications in Nonlinear Science and Numerical Simulation, 2016,32:145.

[8]ABADA N,CHAHDANE H,HAMMOUCHE H. Existence Results for Impulsive Partial Functional Fractional Diffrential Equation with State Dependent Delay[J]. Nonlinear Analysis : Problems,Applications and Computational Methods,2021,168:1.

[9]DHAYAL R,MALIK M,ABBA S. Existence,Stability and Controllability Results of Stochastic Differential Equations with Non-instantaneous Impulses[J].International Journal of Control. 2022.95(7) :1719.

[10] ZHOU Xia, ZHOU Dongpeng, ZHONG Shouming. Existence and Exponential Stability in the Pth Moment for Impulsive Neutral Stochastic Integro-differential Equations Driven by Mixed Fractional Brownian Motion[J]. Journal of Inequalities and Applications,2019,2019(1):1.

[11]DIEYE M,ABDOUL M,EZZINBI K. On Exponential Stability of Mild Solutions for Some Stochastic Partial Integro-differential Equations[J]. Statistics and Probability Letters,2017,123:61.

[12] SINGH V, CHAUDHARY R, PANDEY D N. Approximate Controllability of Second-order Nonautonomous Stochastic Impulsive Differential Systems[J]. Stochastic Analysis and Applications,2020,39(2): 339.

[13]DIOP A,DIOP M A,EZZINBI K,et al.Optimal Controls Problems for Some Impulsive Stochastic Integro-differential Equations with State-dependent Delay[J]. Stochastics,2022,94(8) : 1186.

[14]HUANG Hao,WU Zheng,HU Ling,et al.Existence and Controllability of Second-order Neutral Impulsive Stochastic Evolution Integro-differential Equations with Statedependent Delay[J]. Journal of Fixed Point Theory and Applications,2018,

[15]DUAN Pengju,REN Yong. Solvability and Stability for Neutral Stochastic Integro-diferential Equations Driven by Fractional Brownian Motion with Impulses[J]. Mediterranean Journal of Mathematics,2018,15(6): 1.

[16]WANG Rongnian,XIAO Tijun,LIANG Jin. A Note on the Fractional Cauchy Problems with Nonlocal Conditions [J].Applied Mathematics Letters,2011,24(8) :1435.

[17] SAKTHIVEL R,REN Yong,Debbouche A,et al. Approximate Controllability of Fractional Stochastic Differential Inclusions with Nonlocal Conditions[J].Applicable Analysis,2016,95(11) :2361.

[18]YAN Zuomao,LU Fangxia. Exponential Stability for Nonautonomous Impulsive Neutral Partial Stochastic Evolution Equations with Delay[J]. International Journal of Control,2019,92(9) : 2037.

[19] CHEN Pengyu,LI Yongxiang, ZHANG Xuping. Cauchy Problem for Stochastic Non-autonomous Evolution Equations Governed by Noncompact Evolution Families[J]. Discrete and Continuous Dynamical Systems Series B, 2021,26(3): 1531.

[20] CHEN Pengyu, ZHANG Xuping,LI Yongxiang. Non-autonomous Parabolic Evolution Equations with Non-instantaneous Impulses Governed by Noncompact Evolution Families[J]. Journal of Fixed Point Theory and Applications,2019,21(3) :1.

(編輯:溫澤宇)

猜你喜歡
定義系統
Smartflower POP 一體式光伏系統
工業設計(2022年8期)2022-09-09 07:43:20
WJ-700無人機系統
永遠不要用“起點”定義自己
海峽姐妹(2020年9期)2021-01-04 01:35:44
ZC系列無人機遙感系統
北京測繪(2020年12期)2020-12-29 01:33:58
定義“風格”
基于PowerPC+FPGA顯示系統
半沸制皂系統(下)
連通與提升系統的最后一塊拼圖 Audiolab 傲立 M-DAC mini
成功的定義
山東青年(2016年1期)2016-02-28 14:25:25
修辭學的重大定義
當代修辭學(2014年3期)2014-01-21 02:30:44
主站蜘蛛池模板: 91日本在线观看亚洲精品| 91精品国产自产在线老师啪l| 亚洲视频免费播放| 亚洲日韩每日更新| 国产成人高清在线精品| 99无码中文字幕视频| 91丝袜美腿高跟国产极品老师| 日韩大片免费观看视频播放| 久久永久视频| 丁香六月综合网| 成人免费午夜视频| 日韩亚洲综合在线| 日本中文字幕久久网站| 免费jjzz在在线播放国产| 精品国产99久久| 国产精品视频公开费视频| 久久激情影院| 国产91视频免费| 国产成人精品一区二区三区| 国产在线观看91精品| 国产成人精品高清在线| 国产精品亚洲综合久久小说| 国产剧情一区二区| 国内精品手机在线观看视频| 国产不卡一级毛片视频| 中文字幕乱码二三区免费| 亚洲欧美成人综合| 国产亚洲精品资源在线26u| 青青草国产免费国产| 国产剧情国内精品原创| 无码专区国产精品一区| 一级香蕉视频在线观看| 老司机精品一区在线视频| 国产打屁股免费区网站| 婷婷五月在线| 亚洲精品黄| 综合久久五月天| 久久综合婷婷| 香蕉久久永久视频| 尤物亚洲最大AV无码网站| 99在线观看国产| 国产在线精品网址你懂的| 老司机精品99在线播放| 亚洲国产高清精品线久久| 亚洲天堂免费在线视频| 国产永久在线视频| 亚洲成A人V欧美综合| 少妇被粗大的猛烈进出免费视频| 成年人福利视频| 久久先锋资源| 免费观看三级毛片| 久久九九热视频| 国模粉嫩小泬视频在线观看| 欧美一级高清片久久99| 国产午夜在线观看视频| 视频二区国产精品职场同事| 欧美三级视频在线播放| 久久超级碰| 手机精品福利在线观看| 精品视频一区在线观看| 亚洲视频在线观看免费视频| 亚洲欧美另类色图| 日本一区二区三区精品国产| 久久国产精品影院| 欧美激情视频二区| 亚洲视频影院| 综合网天天| 亚洲v日韩v欧美在线观看| 国产在线欧美| 欧美精品啪啪一区二区三区| 久久综合亚洲鲁鲁九月天| 婷婷丁香在线观看| 91网在线| 福利小视频在线播放| 操国产美女| 97se亚洲综合不卡| 尤物视频一区| 免费无码AV片在线观看中文| 最新国产高清在线| 国产91在线免费视频| 久久综合久久鬼| 国产精品女主播|