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

    有界擾動下約束非線性系統(tǒng)魯棒經(jīng)濟(jì)模型預(yù)測控制

    2022-03-10 11:04:38何德峰王青松
    自動化學(xué)報 2022年2期
    關(guān)鍵詞:魯棒閉環(huán)擾動

    何德峰 韓 平 王青松

    近年來,經(jīng)濟(jì)模型預(yù)測控制(Economic model predictive control,EMPC)在工業(yè)界和學(xué)術(shù)界引起了廣泛關(guān)注[1-2].作為一種新近發(fā)展的先進(jìn)控制技術(shù),EMPC 有望成為解決復(fù)雜系統(tǒng)節(jié)能、降耗和增效優(yōu)化控制問題的重要手段,目前已應(yīng)用于能源、造紙、車輛等系統(tǒng)的能效優(yōu)化控制[3-9].除具有常規(guī)模型預(yù)測控制 (Model predictive control,MPC)的顯式處理約束和多變量控制的優(yōu)點外,EMPC 還能優(yōu)化 “經(jīng)濟(jì)”類目標(biāo)函數(shù),通常這類函數(shù)不是設(shè)定值跟蹤偏差的正定函數(shù),而是系統(tǒng)狀態(tài)和控制變量的非凸或非正定函數(shù)[1-9].因此,把以設(shè)定值跟蹤偏差的正定函數(shù)為優(yōu)化目標(biāo)的常規(guī)MPC 稱為目標(biāo)跟蹤MPC,而不以跟蹤偏差的正定函數(shù)為優(yōu)化目標(biāo)的MPC 統(tǒng)稱為經(jīng)濟(jì)MPC[2].現(xiàn)有研究表明:經(jīng)濟(jì)最優(yōu)性目標(biāo)與閉環(huán)系統(tǒng)的穩(wěn)定性目標(biāo)具有一定的沖突性[1-2],因此近年來EMPC 的穩(wěn)定性綜合策略得到了廣泛研究.

    為建立EMPC 關(guān)于經(jīng)濟(jì)平衡點的穩(wěn)定性,一種主要方法是構(gòu)造基于經(jīng)濟(jì)優(yōu)化目標(biāo)函數(shù)的Lyapunov 函數(shù)[5,10-17].例如,使用終端等式約束和強(qiáng)對偶性假設(shè),定義經(jīng)濟(jì)目標(biāo)函數(shù)的旋轉(zhuǎn)代價函數(shù)并將其作為閉環(huán)系統(tǒng)的一個Lyapunov 函數(shù)[10],而引入廣義終端等式約束[15],建立了經(jīng)濟(jì)性能變化下的遞推可行性與閉環(huán)系統(tǒng)的有界穩(wěn)定性[11].進(jìn)一步,采用嚴(yán)格耗散性條件、終端不等式約束和終端代價函數(shù)代替,降低了EMPC 穩(wěn)定性綜合策略的保守性[12-13].在無終端約束EMPC 策略中,閉環(huán)軌跡在足夠長的預(yù)測時域情況下收斂到平衡點的鄰域[16-17].雖然無終端約束增大了閉環(huán)系統(tǒng)的吸引域,但長時域預(yù)測將大大增加了在線優(yōu)化的計算負(fù)擔(dān).進(jìn)一步,EMPC穩(wěn)定性和經(jīng)濟(jì)性是一對存在沖突的控制目標(biāo)[18-19],且穩(wěn)定性和經(jīng)濟(jì)性目標(biāo)無法統(tǒng)一度量,難以通過權(quán)重標(biāo)定.對此,從多目標(biāo)優(yōu)化控制角度,考慮非線性系統(tǒng)強(qiáng)對偶性或耗散性條件難以滿足情況,文獻(xiàn)[20-21]通過構(gòu)造穩(wěn)定性收縮約束,建立閉環(huán)系統(tǒng)關(guān)于最優(yōu)經(jīng)濟(jì)平衡點的漸近穩(wěn)定性.

    實際系統(tǒng)總是存在不確定擾動,現(xiàn)有EMPC 策略通常難以保證受擾系統(tǒng)的可行性和穩(wěn)定性.對于目標(biāo)跟蹤MPC,目前已有較多魯棒穩(wěn)定性結(jié)果[22-32],主要包括本質(zhì)魯棒MPC[22]、Tube 魯棒MPC[24-25]以及min-max MPC[26-32]等,其中min-max MPC 采用微分對策原理,在使最壞擾動輸入情況下系統(tǒng)的性能指標(biāo)上界達(dá)到最小.相比于本質(zhì)魯棒MPC和Tube 魯棒MPC,min-max MPC 能大大降低魯棒MPC 的保守性,但會增加優(yōu)化問題的在線計算量[1].為降低min-max MPC 的在線計算量,文獻(xiàn)[30]采用仿射輸入結(jié)構(gòu),使MPC 含有抑制擾動的閉環(huán)成分和易于求解的開環(huán)優(yōu)化.另一方面,輸入到狀態(tài)穩(wěn)定性(Input-to-state stability,ISS)成為分析不確定系統(tǒng)魯棒穩(wěn)定性的一個有效工具[23,27-32],并應(yīng)用到了EMPC 魯棒性研究,如文獻(xiàn)[33-34]采用強(qiáng)對偶性假設(shè)和約束緊縮方法,證明了周期性擾動下線性系統(tǒng)EMPC 閉環(huán)收斂性,文獻(xiàn)[35]獲得了非線性系統(tǒng)EMPC 的有界穩(wěn)定性結(jié)果,提高了經(jīng)濟(jì)性能優(yōu)化的靈活性,文獻(xiàn)[36]將穩(wěn)定性目標(biāo)和經(jīng)濟(jì)性目標(biāo)相加,證明EMPC 線性系統(tǒng)關(guān)于經(jīng)濟(jì)目標(biāo)的最大值是ISS 的,文獻(xiàn)[37]施加保證魯棒穩(wěn)定性的顯式收縮約束,提出兩種非線性系統(tǒng)魯棒EMPC算法,文獻(xiàn)[38]提出Lipschitz連續(xù)非線性系統(tǒng)的隱式收縮魯棒EMPC 策略,提高了系統(tǒng)的平均經(jīng)濟(jì)性能.

    本文針對含有未知有界擾動的不確定非線性系統(tǒng),提出一種新的具有遞推可行性以及ISS 保證的魯棒EMPC 策略.該策略明確考慮經(jīng)濟(jì)最優(yōu)性和魯棒穩(wěn)定性控制目標(biāo)的矛盾特點,采用微分對策原理在線滾動優(yōu)化計算這對沖突目標(biāo)的min-max 問題.離線計算最優(yōu)經(jīng)濟(jì)平衡點,并利用狀態(tài)與該平衡點的偏差定義魯棒穩(wěn)定性目標(biāo)函數(shù),而經(jīng)濟(jì)目標(biāo)函數(shù)則由系統(tǒng)的經(jīng)濟(jì)性能給定.通過特殊設(shè)計EMPC 優(yōu)化問題的隱式收縮約束,并在魯棒穩(wěn)定性目標(biāo)優(yōu)化問題中引入一個新約束,保證EMPC 優(yōu)化的遞推可行性和閉環(huán)系統(tǒng)關(guān)于不確定擾動輸入的ISS.相比現(xiàn)有魯棒EMPC 策略,本文首先建立了約束非線性系統(tǒng)具有ISS 的魯棒EMPC 策略;其次,EMPC 遞推可行性和魯棒穩(wěn)定性無需強(qiáng)對偶性或耗散性假設(shè)條件,從而擴(kuò)大了魯棒 EMPC 的應(yīng)用范圍;最后,采用微分對策原理得到了保守性更低的容許擾動上界.采用一個受擾非線性連續(xù)攪拌釜反應(yīng)器(Continuous stirred tank reactor,CSTR)的仿真實例,驗證本文提出策略的有效性與優(yōu)越性.

    符號說明:Z+表示非負(fù)整數(shù)集,Ia:b表示集合{i ∈Z+:a ≤i ≤b,a ∈Z+,b ∈Z+},I≥j表示集合{i ∈Z+:i ≥j,j ∈Z+}.x+表示x的下一時刻狀態(tài),|x|表示x的歐幾里得范數(shù),‖x‖=sup{|x(k)|,k ∈Z+},x(i|k) 表示在第k時刻對未來第k+i時刻的預(yù)測變量.連續(xù)函數(shù)h1:R+→R+稱為K類函數(shù)系指該函數(shù)單調(diào)遞增,且h1(0)=0 ;函數(shù)h2:R+→R+稱為K∞類 函數(shù)系指該函數(shù)是K類函數(shù),且當(dāng)s →∞時,有h(s)→∞;函數(shù)h3:R+×R+→R+稱為KL函 數(shù)系指對于任意固定的t≥0,h3(·,t) 是K類函數(shù),而對于任意固定的s>0,h3(s,·) 是單調(diào)遞減,且當(dāng)t→∞時,h3(s,t)→0 .

    1 問題描述

    考慮不確定離散時間非線性系統(tǒng)

    其中,x(k)∈Rn,u(k)∈Rm,w(k)∈Rq和z(k)∈Rs分別表示系統(tǒng)在k時刻的狀態(tài)、控制輸入、擾動輸入和輔助輸出.假設(shè)f1(·),f2(·),f3(·)和g(·) 分別為定義在 Rn上的光滑函數(shù),滿足f1(0)=0,f2(0)=0,f3(0)=0 以及g(0)=0 .系統(tǒng)(1) 在k時刻的狀態(tài)x(k)=φ(k;x0,u,w) 由初始狀態(tài)x0,控制序列u={u(0),u(1),···}和擾動序列w={w(0),w(1),···}表示.假設(shè)系統(tǒng)(1)的狀態(tài)完全可測,且系統(tǒng)狀態(tài)和控制輸入滿足約束

    其中,集合X?Rn和U?Rm均為凸的緊集,且它們內(nèi)部包含某些平衡點 (xe,ue) .

    假設(shè) 1.擾動輸入w滿足如下約束

    其中,W為包含平衡點 (xe,ue) 的緊集.滿足約束(3)的擾動稱為容許擾動.注意w表示參數(shù)不確定性、模型失配以及持續(xù)外部擾動等多種有界不確定性[29],其上界表示為‖w‖=sup{|wt|:wt∈W,t ∈Z+}.

    考慮系統(tǒng)(1)的經(jīng)濟(jì)性能函數(shù)Le:X×U →R,基于系統(tǒng)(1)的名義模型x+=f1(x)+f2(x)u離線計算如下最優(yōu)經(jīng)濟(jì)平衡點

    不失一般性,后文假設(shè)系統(tǒng)(1)的最優(yōu)經(jīng)濟(jì)平衡點為原點.

    下面介紹輸入到狀態(tài)穩(wěn)定(ISS)相關(guān)理論,首先回顧魯棒不變集的概念.考慮一般離散不確定非線性系統(tǒng)x+=f(x,w) .

    定義 1[27].考慮系統(tǒng)x+=f(x,w)和集 Θ?Rn,如果對于任意x∈Θ和w∈W,該系統(tǒng)滿足x+=f(x,w)∈Θ,則稱 Θ 為該系統(tǒng)的一個魯棒不變集.

    定義 2[27].考慮系統(tǒng)x+=f(x,w) 及其內(nèi)含原點的魯棒不變集 Θ?X,如果對于任意初始狀態(tài)x0∈Θ 及w∈W,存在K∞類函數(shù)α和KL類函數(shù)β使系統(tǒng)滿足

    則該系統(tǒng)在 Θ 內(nèi)是ISS 的.

    引理 1[29].考慮系統(tǒng)x+=f(x,w) 及其內(nèi)含原點的魯棒不變集 Θ?X,如果對于任意x∈Θ和w ∈W,存在連續(xù)函數(shù)V:Rn→R+滿足

    其中,ξ1,ξ2和ξ3為K∞類函數(shù),ρ1和ρ2為K類函數(shù),則該系統(tǒng)在 Θ 內(nèi)ISS,稱V為該系統(tǒng)的ISSLyapunov 函數(shù).

    注 1.由定義2 可知,當(dāng)系統(tǒng)不受擾動或僅受衰減擾動作用時,系統(tǒng)最終在原點處漸近穩(wěn)定;當(dāng)受持續(xù)有界擾動作用時,系統(tǒng)有界穩(wěn)定,且狀態(tài)軌跡最終收斂的范圍與持續(xù)擾動的上界有關(guān).

    本文目標(biāo)是針對不確定非線性系統(tǒng)(1),通過極小化經(jīng)濟(jì)性能函數(shù)在線計算魯棒EMPC 控制器,要求相應(yīng)閉環(huán)系統(tǒng)滿足約束條件(2),且閉環(huán)系統(tǒng)的最優(yōu)經(jīng)濟(jì)平衡點相對于容許擾動(3)具有ISS.

    2 魯棒經(jīng)濟(jì)模型預(yù)測控制

    考慮有限預(yù)測步長N∈I≥1,定義k∈Z+時刻的N步控制序列u(k)={u(0|k),u(1|k),···,u(N-1|k)}、擾動序列w(k)={w(0|k),w(1|k),···,w(N-1|k)}以及對應(yīng)的預(yù)測狀態(tài)序列x(k)={x(1|k),x(2|k),···,x(N|k)},并且x(i|k)=φ(i;x(k),u(k),w(k)),其中x(k) 為當(dāng)前時刻系統(tǒng)狀態(tài).考慮系統(tǒng)(1)的經(jīng)濟(jì)性能函數(shù)Le,定義N步經(jīng)濟(jì)目標(biāo)函數(shù)

    其中,x(0|k)=x(k)

    為計算魯棒EMPC 控制器,在每個時刻優(yōu)化經(jīng)濟(jì)目標(biāo)函數(shù)(7).由于擾動的存在,在魯棒EMPC 中求解如下min-max 經(jīng)濟(jì)最優(yōu)控制問題:

    其中,(u*(k),w*(k)) 為k時刻的經(jīng)濟(jì)最優(yōu)解,對應(yīng)最優(yōu)預(yù)測狀態(tài)序列x*(k) ;x(0|k)=x(k) 為初始條件;x(N|k)∈Ω 為終端狀態(tài)約束,終端約束集 Ω?X且內(nèi)含原點;(8e)為待設(shè)計的魯棒穩(wěn)定性收縮約束,λ∈[0,1) 為收縮因子.進(jìn)一步定義如下關(guān)于最優(yōu)經(jīng)濟(jì)平衡點的魯棒穩(wěn)定性目標(biāo)函數(shù):

    其中,函數(shù)L:X×U →R+,Lw:W →R+和E:X →R+為連續(xù)有界函數(shù).為書寫簡潔,令VR(x)=VR(x,u,w) .

    假設(shè) 2.存在K∞類函數(shù)αl,αw和βw,使魯棒穩(wěn)定性目標(biāo)函數(shù)VR(x) 滿足:L(x,u)≥αl(|x|),?x ∈X,?u∈U;αw(|w|)≤Lw(w)≤βw(|w|),?w ∈W.

    假設(shè) 3.在終端約束集 Ω 內(nèi)存在局部反饋控制律u=π(x) 使 得π(x)∈U,?x ∈X和如下不等式成立:

    其中,αf(|x|)≤E(x)≤βf(|x|),αf和βf為K∞類函數(shù).

    注意,局部控制律u=π(x) 可采用如H∞控制方法求解[26,30].為此,定義對稱矩陣

    其中,P為對稱正定矩陣,I為單位矩陣.令R=R(0) .不失一般性,令擾動抑制水平γ=1 .

    引理 2[30].假設(shè)矩陣P滿足Riccati 不等式方程

    注 2.如果假設(shè)3 成立,則對于任意w∈W,Ω為閉環(huán)系統(tǒng)(13)的一個魯棒不變集,同時E(x) 為該系統(tǒng)在 Ω 內(nèi)的一個ISS-Lyapunov 函數(shù).文獻(xiàn)[30]給出了保證 Ω 魯棒不變性的一個充分條件,即容許擾動w的上界

    其中,lx和lw分別為閉環(huán)系統(tǒng)(13)關(guān)于x和w的局部Lipschitz 常數(shù),λmax和λmin分別為矩陣P的最大和最小特征值.如擾動w滿足式(15),則 Ω 為閉環(huán)系統(tǒng)(13)的魯棒不變集,且π(x)∈U,?x ∈Ω .

    下面構(gòu)造收縮約束函數(shù)η.令x(k)為當(dāng)前時刻系統(tǒng)狀態(tài),求解如下有限時域魯棒性最優(yōu)控制問題

    其中,(u*(i|k-1),w*(i|k-1)) 是 優(yōu)化問題(8)在k-1 時刻最優(yōu)解 (u*(k-1),w*(k-1)) 的 分量.將 (u*(k-1),w*(k-1))和(u0(k),w0(k)) 分別代入VR(x),據(jù)此定義函數(shù)

    如果優(yōu)化問題(8)在k時刻可行,則根據(jù)滾動時域控制原理,將u*(k) 的第1 個分量定義為魯棒EMPC 控制律,即u(k)=umpc(x(k))=u*(0|k),對應(yīng)閉環(huán)系統(tǒng)為

    計算魯棒EMPC 控制器的算法總結(jié)如下:

    算法 1.魯棒EMPC 算法

    步驟 1.初始化.設(shè)定N∈I≥1,λ ∈[0,1),λ1≥1,Le(x,u),L(x,u),Lw(w)和E(x) ;離線計算π(x)以及 Ω ;令k=0,考慮初始狀態(tài)x0,令η充分大;求解優(yōu)化問題(8),得到經(jīng)濟(jì)最優(yōu)解 (u*(0),w*(0)),并將u*(0) 的首個分量u*(0|0) 作用于系統(tǒng)(1).

    步驟 2.在k時刻,利用k-1 時刻經(jīng)濟(jì)最優(yōu)解(u*(k-1),w*(k-1)) 構(gòu) 造(u1(k),w1(k)),更 新(16e);求解優(yōu)化問題(16),得到當(dāng)前時刻魯棒最優(yōu)解 (u0(k),w0(k)) .

    步驟 3.計算和,更新η;求解優(yōu)化問題(8),得到當(dāng)前時刻經(jīng)濟(jì)最優(yōu)解(u*(k),w*(tk)) .

    步驟 4.將u*(k) 的首個分量u*(0|k) 作用于系統(tǒng)(1).

    步驟 5.令k=k+1,測量系統(tǒng)(1)的狀態(tài);返回步驟2.

    為更清晰地描述算法1 的運(yùn)行過程,圖1 給出了每個時刻兩個優(yōu)化問題的求解順序.

    由算法1和圖1 可知,在初始時刻k=0,由于約束(8e)不起作用,故在初始時刻無需求解優(yōu)化問題(16),并可令η(x0,λ)→∞.

    圖1 魯棒EMPC 算法運(yùn)行過程示意圖Fig.1 A schematic diagram of the robust EMPC algorithm

    3 輸入到狀態(tài)穩(wěn)定性

    閉環(huán)系統(tǒng)(19)具有ISS 性質(zhì)的前提是需要保證算法1 具有遞推可行性,即兩個優(yōu)化問題在每個時刻都至少存在一組可行解(不一定最優(yōu)),使得對于任意容許擾動,經(jīng)濟(jì)優(yōu)化問題(8)和魯棒性優(yōu)化問題(16)的所有約束均滿足.注意,(8b)~(8d)與(16b)~(16d)具有相同約束.

    考慮約束不確定系統(tǒng)(1)~(3)和狀態(tài)τ=x(0|k)∈X.對于優(yōu)化問題(8),定義N步可行控制序列集為

    其中,N步擾動約束集WN=W ×W ×···×W.

    定義 3.考慮系統(tǒng)(1)和狀態(tài)τ=x(0|k)∈X,如果可行控制序列集UN(τ) 非空,則稱τ為該系統(tǒng)的初始可行狀態(tài).所有初始可行狀態(tài)τ的集合稱為系統(tǒng)的初始可行狀態(tài)集XN.

    定理 1.如果假設(shè)1~3 成立,且容許擾動滿足式(15),則優(yōu)化問題(16)具有遞推可行性.

    證明.考慮序列對 (u1(k),w1(k)) 并代入系統(tǒng)(1)得到狀態(tài)序列x1(k)={x*(1|k-1),···,x*(N|k-1),x1(N|k)},其中x*(N|k-1)∈Ω 是對應(yīng)于u*(k-1)的終端預(yù)測狀態(tài).由假設(shè)3和注2 可知,Ω 是閉環(huán)系統(tǒng)(13)的魯棒不變集,故對于任意w1(N -1|k)∈W,u1(N -1|k)=π(x*(N|k-1))∈U和x1(N|k)=φ(N;x*(N|k-1),u1(N -1|k),w1(N -1|k))∈Ω 成立.利用序列對 (u1(k),w1(k)) 構(gòu)造優(yōu)化問題(16)在k時刻的備選解

    其對應(yīng)狀態(tài)序列x2(k)={x1(1|k),···,x1(N|k),x2(N|k)},其中x1(N|k)∈Ω.又 Ω 是魯棒不變的,則對于任意w2(N -1|k)∈W,u2(N -1|k)=π(x1(N|k))∈U和x2(N|k)=φ(N;x1(N|k),u2(N-1|k),w2(N -1|k))∈Ω 成立,因此約束(16b)~(16d)滿足.

    再分別 將 (u1(k),w1(k))和(u2(k),w2(k)) 代入函數(shù)VR(x) 中,計算得

    故約束(16e)滿足,從而 (u2(k),w2(k)) 為問題(16)在k時刻的可行解,即優(yōu)化問題(16)具有遞推可行性.□

    定理 2.如果假設(shè)1~3 成立,且容許擾動滿足式(15),則優(yōu)化問題(8)具有遞推可行性.

    故約束(8e) 成立.由定理1 可知,(u0(k),w0(k))作為優(yōu)化問題(16)在k時刻的最優(yōu)解滿足約束(8b)~(8d).因此,在k時刻總存在可行控制序列u0(k) 滿足優(yōu)化問題(8)的所有約束,從而優(yōu)化問題(8)具有遞推可行性.□

    根據(jù)定理1和定理2,本文魯棒EMPC 策略中的雙目標(biāo)優(yōu)化問題皆具有遞推可行性.下面將基于定理1和定理2 給出閉環(huán)系統(tǒng)(19)的ISS 結(jié)果.

    定理 3.如果假設(shè)1~3 成立,且容許擾動滿足式(15),則當(dāng)優(yōu)化問題(8)在初始時刻存在可行解時,閉環(huán)系統(tǒng)(19)在魯棒不變集XN內(nèi)相對于擾動具有ISS.

    證明.由定理1和定理2 可知,當(dāng)優(yōu)化問題(8)在初始時刻存在可行解時,該優(yōu)化問題在任意k時刻都是可行的,則對任意x(k)∈XN和w(k)∈W,閉環(huán)系統(tǒng)(19)滿足x(k+1)∈XN.由定義3 可知,XN為閉環(huán)系統(tǒng)(19)的一個魯棒不變集.

    成立.則聯(lián)立不等式(26),(32)和(36),由引理1 可知,值函數(shù)VN(x(k))為閉環(huán)系統(tǒng)(19)的一個ISSLyapunov 函數(shù).因此,閉環(huán)系統(tǒng)(19)在XN內(nèi)相對于擾動具有輸入到狀態(tài)穩(wěn)定性.□

    4 實例仿真

    考慮不確定非線性連續(xù)攪拌釜反應(yīng)器(CSTR)

    其中,cA和cB分別為組分A和B的濃度,Q為反應(yīng)器進(jìn)料流量,cAf和cBf分別為進(jìn)料中組分A和B的濃度,Δc為進(jìn)料中組分A 濃度的不確定波動,體積V和反應(yīng)動力學(xué)參數(shù)k0.當(dāng)忽略進(jìn)料中組分濃度的波動,該模型廣泛用于名義穩(wěn)定EMPC 綜合策略的驗證[10,13,21].這里假設(shè)進(jìn)料中組分A的濃度波動是有界的,用于驗證本文魯棒EMPC 的有效性.取模型參數(shù)[10]:cAf=1.0 mol/l,cBf=0,V=10 l和k0=1.2 l/(mol·min).

    令 [cA,cB]T為系統(tǒng)狀態(tài)[x1,x2]T,Q為控制輸入u,w為擾動 Δc.進(jìn)一步,定義狀態(tài)約束和控制約束

    選擇采樣周期Ts=0.1 min,利用歐拉差分法離散化連續(xù)時間模型(37),得CSTR 非線性離散時間模型

    其中,k為采樣時刻.定義CSTR 經(jīng)濟(jì)性能指標(biāo)[10]

    進(jìn)一步,為保證終端狀態(tài)集 Ω 的魯棒不變性,由式(15)計算容許擾動上界為0.143 6.

    令預(yù)測步長N=5,仿真總步長Tsim=100,系數(shù)λ1=1,并考慮不同收縮因子λ∈[0,1) 的控制效果.選擇系統(tǒng)初始狀態(tài)x0=(0.15,0.7)T,仿真結(jié)果如圖2和圖3 所示,圖2 對應(yīng)持續(xù)擾動w(k)=0.1436 sin(k/2)的仿真結(jié)果,圖3 對應(yīng)衰減擾動w(k)=0.1436 exp(-k/10) 的仿真結(jié)果.分析圖2 可知,在持續(xù)擾動下,閉環(huán)系統(tǒng)最終在最優(yōu)經(jīng)濟(jì)平衡點附近有界穩(wěn)定;分析圖3 可知,在衰減擾動下,閉環(huán)狀態(tài)軌跡漸近收斂于.根據(jù)圖2和圖3 可知,對于不同收縮因子λ∈[0,1),由算法1 得到的魯棒EMPC 控制器使得閉環(huán)系統(tǒng)在相對于容許擾動總是ISS 的,但不同λ對應(yīng)的閉環(huán)狀態(tài)軌跡的動態(tài)響應(yīng)不同.以衰減擾動w(k)=0.1436 exp(-k/10)為例,由圖3 可知,因子λ越小,閉環(huán)狀態(tài)軌跡的收斂過渡時間Ttr越短,具體數(shù)據(jù)如表1 最右列所示.

    圖2 持續(xù)擾動 w (k)=0.1436 sin(k/2) 的仿真結(jié)果Fig.2 Simulation results under continuous disturbance w(k)=0.1436 sin(k/2)

    圖3 持續(xù)擾動 w (k)=0.1436 exp(-k/10) 的仿真結(jié)果Fig.3 Simulation results under continuous disturbance w(k)=0.1436 exp(-k/10)

    定義閉環(huán)系統(tǒng)T步平均經(jīng)濟(jì)性能指標(biāo)

    取T=100,表1 給出了不同λ值對應(yīng)的閉環(huán)系統(tǒng)平均經(jīng)濟(jì)性能.分析表1 可以看出,因子λ越大,閉環(huán)系統(tǒng)平均經(jīng)濟(jì)性能越好.這表明,無論是名義系統(tǒng)還是不確定系統(tǒng),閉環(huán)系統(tǒng)的經(jīng)濟(jì)最優(yōu)性和穩(wěn)定性是互相沖突的雙控制目標(biāo)[10,21].在本文魯棒EMPC 策略中,可以通過調(diào)節(jié)收縮因子λ對經(jīng)濟(jì)性控制目標(biāo)和魯棒穩(wěn)定性控制目標(biāo)進(jìn)行權(quán)衡,從而實現(xiàn)經(jīng)濟(jì)性和魯棒穩(wěn)定性綜合控制效果.

    表1 平均經(jīng)濟(jì)性能和收斂過渡時間Table 1 Average economic performance and transient time

    進(jìn)一步,為了驗證本文魯棒EMPC 策略(簡記為MM-EMPC)的優(yōu)越性,對比文獻(xiàn)[38]中魯棒收縮EMPC 策略(簡記為RC-EMPC),研究不同容許擾動上界對系統(tǒng)穩(wěn)定性的影響.對約束CSTR 系統(tǒng)(37)和(38),在系統(tǒng)參數(shù)相同情況下,RC-EMPC 策略容許擾動上界d1=0.0136,MM-EMPC 容許擾動上界d2=0.1436.選擇持續(xù)擾動w(k)=0.04sin(k/5),令初始狀態(tài)x0=(0.5,0.5)T,仿真步長Tsim=500.在相同仿真環(huán)境下,兩種魯棒EMPC 策略的控制結(jié)果如圖4 所示,其中,實線表示MM-EMPC 策略仿真結(jié)果,虛線表示RC-EMPC 策略仿真結(jié)果.

    圖4 閉環(huán)狀態(tài)軌跡和控制輸入曲線Fig.4 Closed-loop state trajectorand control input profiles

    從圖4 可以看出,MM-EMPC 策略對應(yīng)的閉環(huán)狀態(tài)軌跡和控制輸入曲線均收斂于附近,而RC-EMPC 策略對應(yīng)的閉環(huán)軌跡雖然收斂,但存在較大的穩(wěn)態(tài)誤差,且MM-EMPC 閉環(huán)系統(tǒng)的動態(tài)響應(yīng)更加快速.進(jìn)一步,考慮衰減擾動w(k)=0.04 exp(-k/5),選擇4 個不同的初始狀態(tài)對兩種控制策略進(jìn)行仿真對比分析,控制結(jié)果如圖5 所示.綜合圖4和圖5 分析可知,當(dāng)容許擾動上界d ∈(d1,d2) 時,MM-EMPC 策略下的閉環(huán)系統(tǒng)相對于擾動是ISS 的,而RC-EMPC 卻無法保證閉環(huán)系統(tǒng)相對于擾動具有ISS,其中一個重要原因是RCEMPC 對大擾動會丟失遞推可行性.因此,相較于RC-EMPC 策略,本文策略在保證閉環(huán)系統(tǒng)ISS 同時,能夠獲得更大的容許擾動上界,從而降低魯棒EMPC 控制器的保守性.

    圖5 不同初始狀態(tài)的閉環(huán)系統(tǒng)相軌跡Fig.5 Phase trajectories of the closed-loop system from different initial states

    5 結(jié)束語

    本文針對有界擾動下的約束不確定仿射輸入非線性系統(tǒng),提出了一種新的魯棒EMPC 策略.基于微分對策原理分別對經(jīng)濟(jì)目標(biāo)函數(shù)和關(guān)于最優(yōu)經(jīng)濟(jì)平衡點的魯棒穩(wěn)定性目標(biāo)函數(shù)進(jìn)行優(yōu)化,利用得到的魯棒穩(wěn)定性目標(biāo)最優(yōu)值函數(shù)構(gòu)造隱式收縮約束,保證了雙控制目標(biāo)優(yōu)化問題的遞推可行性,并建立了閉環(huán)系統(tǒng)在最優(yōu)經(jīng)濟(jì)平衡點處相對于擾動的輸入到狀態(tài)穩(wěn)定性結(jié)果.通過對不確定CSTR 經(jīng)濟(jì)優(yōu)化控制的對比仿真實驗,驗證了本文策略的有效性和優(yōu)越性.

    盡管基于微分對策的min-max 魯棒EMPC 在理論上能有效提高魯棒EMPC 的性能,但minmax 優(yōu)化在線計算復(fù)雜,將阻礙魯棒EMPC 在快速響應(yīng)系統(tǒng)中的應(yīng)用.因此,降低min-max 魯棒EMPC 的在線優(yōu)化計算量和設(shè)計更高效的魯棒EMPC 策略(如Tube 魯棒EMPC、本質(zhì)魯棒EMPC 等)將是后續(xù)研究重點.

    猜你喜歡
    魯棒閉環(huán)擾動
    Bernoulli泛函上典則酉對合的擾動
    (h)性質(zhì)及其擾動
    基于學(xué)習(xí)的魯棒自適應(yīng)評判控制研究進(jìn)展
    目標(biāo)魯棒識別的抗旋轉(zhuǎn)HDO 局部特征描述
    單周期控制下雙輸入Buck變換器閉環(huán)系統(tǒng)設(shè)計
    黑龍江電力(2017年1期)2017-05-17 04:25:05
    小噪聲擾動的二維擴(kuò)散的極大似然估計
    雙閉環(huán)模糊控制在石化廢水處理中的研究
    基于Cauchy魯棒函數(shù)的UKF改進(jìn)算法
    用于光伏MPPT中的模糊控制占空比擾動法
    目標(biāo)軌跡更新的點到點魯棒迭代學(xué)習(xí)控制
    久久久久精品国产欧美久久久 | 十八禁高潮呻吟视频| 国产精品二区激情视频| 久久精品成人免费网站| www.熟女人妻精品国产| 亚洲一区中文字幕在线| 在线观看一区二区三区激情| 成人亚洲欧美一区二区av| 如日韩欧美国产精品一区二区三区| 欧美国产精品一级二级三级| 在线 av 中文字幕| 久久亚洲精品不卡| 伊人亚洲综合成人网| 亚洲欧美成人综合另类久久久| 我的亚洲天堂| 黄色 视频免费看| 午夜福利,免费看| 日韩av在线免费看完整版不卡| av又黄又爽大尺度在线免费看| 精品久久久精品久久久| 欧美av亚洲av综合av国产av| 国产福利在线免费观看视频| 多毛熟女@视频| 成年人午夜在线观看视频| 亚洲美女黄色视频免费看| 国产精品人妻久久久影院| 午夜91福利影院| 亚洲 国产 在线| 捣出白浆h1v1| 女警被强在线播放| 久久久久久久久免费视频了| 悠悠久久av| 在线精品无人区一区二区三| 在线观看免费午夜福利视频| 国精品久久久久久国模美| 欧美亚洲日本最大视频资源| 国产亚洲欧美在线一区二区| 99热网站在线观看| 我要看黄色一级片免费的| 精品人妻1区二区| 日韩熟女老妇一区二区性免费视频| 日本午夜av视频| 久久久久久久久免费视频了| 嫩草影视91久久| www.av在线官网国产| 如日韩欧美国产精品一区二区三区| 亚洲情色 制服丝袜| 亚洲熟女毛片儿| 热99久久久久精品小说推荐| 日本av免费视频播放| 久久国产精品人妻蜜桃| 黄网站色视频无遮挡免费观看| 少妇被粗大的猛进出69影院| 午夜av观看不卡| 色视频在线一区二区三区| 啦啦啦啦在线视频资源| 女性生殖器流出的白浆| 国产国语露脸激情在线看| 久久国产亚洲av麻豆专区| 中文字幕亚洲精品专区| 久久国产精品大桥未久av| 菩萨蛮人人尽说江南好唐韦庄| 热re99久久国产66热| 黑人巨大精品欧美一区二区蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 日韩av免费高清视频| 亚洲欧美中文字幕日韩二区| 亚洲国产精品999| 亚洲人成电影观看| 日本av手机在线免费观看| 欧美亚洲日本最大视频资源| 老司机影院毛片| 99久久精品国产亚洲精品| 99九九在线精品视频| 999久久久国产精品视频| 美女高潮到喷水免费观看| 18禁黄网站禁片午夜丰满| 99国产精品99久久久久| 乱人伦中国视频| 国产av精品麻豆| 欧美黄色淫秽网站| 久久国产精品影院| 黄色视频不卡| 国产精品亚洲av一区麻豆| 一本大道久久a久久精品| 少妇精品久久久久久久| 亚洲精品自拍成人| 国产99久久九九免费精品| 久久久久久久大尺度免费视频| 国产一区二区激情短视频 | 日日摸夜夜添夜夜爱| 伊人久久大香线蕉亚洲五| 日韩熟女老妇一区二区性免费视频| 日韩制服骚丝袜av| 无限看片的www在线观看| 搡老岳熟女国产| 亚洲av日韩在线播放| 国产亚洲欧美精品永久| 国产精品一二三区在线看| 脱女人内裤的视频| 91字幕亚洲| 久久女婷五月综合色啪小说| 一边亲一边摸免费视频| 老司机靠b影院| 国产亚洲欧美在线一区二区| 亚洲精品乱久久久久久| 一级黄片播放器| 99re6热这里在线精品视频| 中文字幕人妻熟女乱码| 亚洲欧美一区二区三区久久| 国产真人三级小视频在线观看| 老鸭窝网址在线观看| 搡老乐熟女国产| 国产无遮挡羞羞视频在线观看| 欧美日本中文国产一区发布| 一本—道久久a久久精品蜜桃钙片| 国产老妇伦熟女老妇高清| 亚洲精品久久成人aⅴ小说| 亚洲一区二区三区欧美精品| 真人做人爱边吃奶动态| 精品国产一区二区三区四区第35| 最新的欧美精品一区二区| 亚洲精品久久午夜乱码| 欧美精品高潮呻吟av久久| 亚洲国产精品一区二区三区在线| 久久ye,这里只有精品| 又紧又爽又黄一区二区| 成人亚洲精品一区在线观看| av网站在线播放免费| 97精品久久久久久久久久精品| 欧美黄色片欧美黄色片| bbb黄色大片| 中文字幕精品免费在线观看视频| 久久亚洲国产成人精品v| 中文字幕人妻熟女乱码| www日本在线高清视频| 考比视频在线观看| 亚洲国产av新网站| 国产精品99久久99久久久不卡| 欧美老熟妇乱子伦牲交| 亚洲国产av新网站| 国产精品一区二区精品视频观看| 黄色一级大片看看| 一边亲一边摸免费视频| 久久久精品区二区三区| 肉色欧美久久久久久久蜜桃| 欧美久久黑人一区二区| 成人国产一区最新在线观看 | 亚洲精品第二区| 国产欧美日韩精品亚洲av| 久久这里只有精品19| 99国产精品一区二区三区| 亚洲人成网站在线观看播放| 日韩制服骚丝袜av| 韩国高清视频一区二区三区| 男女之事视频高清在线观看 | 国产免费又黄又爽又色| 在线观看免费日韩欧美大片| 欧美国产精品va在线观看不卡| 国产精品国产av在线观看| av福利片在线| 亚洲精品一卡2卡三卡4卡5卡 | 国产免费视频播放在线视频| 男女午夜视频在线观看| 久久久亚洲精品成人影院| 久久影院123| 在线观看www视频免费| 亚洲三区欧美一区| 国产精品国产三级专区第一集| 精品国产乱码久久久久久男人| 电影成人av| 久久久久国产一级毛片高清牌| 女人爽到高潮嗷嗷叫在线视频| 欧美+亚洲+日韩+国产| 国产精品九九99| 大香蕉久久成人网| 十分钟在线观看高清视频www| 在线观看免费高清a一片| 每晚都被弄得嗷嗷叫到高潮| 新久久久久国产一级毛片| 久久亚洲国产成人精品v| 九草在线视频观看| 亚洲七黄色美女视频| 国产精品人妻久久久影院| 精品国产超薄肉色丝袜足j| 老司机影院毛片| 欧美少妇被猛烈插入视频| 亚洲专区中文字幕在线| 亚洲一区二区三区欧美精品| 国产精品久久久久成人av| 熟女av电影| 亚洲av男天堂| 久久久久久免费高清国产稀缺| 久久人妻福利社区极品人妻图片 | 亚洲精品久久成人aⅴ小说| 最黄视频免费看| 免费日韩欧美在线观看| 亚洲av日韩精品久久久久久密 | 精品人妻一区二区三区麻豆| 欧美日韩视频精品一区| 老司机深夜福利视频在线观看 | 午夜av观看不卡| 欧美 亚洲 国产 日韩一| 欧美老熟妇乱子伦牲交| 大片电影免费在线观看免费| 日本五十路高清| 日韩av在线免费看完整版不卡| 热re99久久国产66热| 少妇被粗大的猛进出69影院| 成人国语在线视频| 天天躁夜夜躁狠狠躁躁| 中文欧美无线码| 一本大道久久a久久精品| 深夜精品福利| 免费av中文字幕在线| 欧美+亚洲+日韩+国产| 男女边摸边吃奶| 亚洲中文日韩欧美视频| www.999成人在线观看| 国产精品人妻久久久影院| 美女视频免费永久观看网站| 亚洲av日韩在线播放| 啦啦啦在线免费观看视频4| 国产淫语在线视频| 国产成人免费观看mmmm| 超碰97精品在线观看| 国产亚洲精品久久久久5区| 欧美久久黑人一区二区| 久久久精品免费免费高清| 国产一区二区在线观看av| 狂野欧美激情性bbbbbb| 欧美xxⅹ黑人| 欧美乱码精品一区二区三区| 亚洲中文av在线| av一本久久久久| 丁香六月天网| 国产成人一区二区三区免费视频网站 | 又紧又爽又黄一区二区| 国产福利在线免费观看视频| 美国免费a级毛片| 亚洲av美国av| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久人人人人人| 国产精品成人在线| 久久国产精品影院| 欧美精品一区二区大全| 国产精品熟女久久久久浪| 九草在线视频观看| 久久影院123| 亚洲第一av免费看| 少妇人妻久久综合中文| 国产欧美日韩一区二区三区在线| 日韩av不卡免费在线播放| 好男人视频免费观看在线| 亚洲精品一二三| 操美女的视频在线观看| 欧美性长视频在线观看| 丁香六月天网| 大香蕉久久网| 又大又爽又粗| 久久99热这里只频精品6学生| 亚洲av成人不卡在线观看播放网 | 91国产中文字幕| 王馨瑶露胸无遮挡在线观看| 精品视频人人做人人爽| 电影成人av| 久久热在线av| 自线自在国产av| 国产免费福利视频在线观看| 亚洲精品国产一区二区精华液| 韩国高清视频一区二区三区| av有码第一页| 少妇的丰满在线观看| 老司机亚洲免费影院| 午夜免费观看性视频| 99国产综合亚洲精品| 电影成人av| 精品欧美一区二区三区在线| 一级毛片我不卡| 老熟女久久久| 人人妻人人澡人人爽人人夜夜| 欧美日韩一级在线毛片| 我的亚洲天堂| 国产爽快片一区二区三区| 亚洲国产精品999| 一级,二级,三级黄色视频| 亚洲情色 制服丝袜| 久久国产精品影院| 久久精品aⅴ一区二区三区四区| 免费在线观看完整版高清| 99香蕉大伊视频| 精品久久久精品久久久| 夫妻性生交免费视频一级片| 久久国产亚洲av麻豆专区| 亚洲国产精品999| 国产成人啪精品午夜网站| 欧美av亚洲av综合av国产av| 极品人妻少妇av视频| 午夜视频精品福利| 麻豆国产av国片精品| 麻豆乱淫一区二区| 亚洲精品日韩在线中文字幕| 欧美精品啪啪一区二区三区 | 伊人亚洲综合成人网| 最近手机中文字幕大全| 国产精品一区二区在线观看99| 校园人妻丝袜中文字幕| 精品人妻熟女毛片av久久网站| 一级黄色大片毛片| 天堂中文最新版在线下载| 亚洲国产欧美网| 老汉色av国产亚洲站长工具| 制服人妻中文乱码| 国产伦理片在线播放av一区| 欧美精品av麻豆av| 久久久欧美国产精品| 中文精品一卡2卡3卡4更新| 新久久久久国产一级毛片| 久久精品久久久久久噜噜老黄| 亚洲精品中文字幕在线视频| 成年人黄色毛片网站| 欧美成人午夜精品| 母亲3免费完整高清在线观看| 亚洲欧美清纯卡通| 国产亚洲精品久久久久5区| 99久久99久久久精品蜜桃| 美女高潮到喷水免费观看| 黄色a级毛片大全视频| 18禁国产床啪视频网站| 国产精品一区二区在线不卡| 久久久精品免费免费高清| 一级片'在线观看视频| 国产精品麻豆人妻色哟哟久久| 国产成人91sexporn| 天天躁夜夜躁狠狠躁躁| 国产精品秋霞免费鲁丝片| 午夜精品国产一区二区电影| 18禁黄网站禁片午夜丰满| 悠悠久久av| 成人黄色视频免费在线看| 日韩一区二区三区影片| 久久九九热精品免费| 国产亚洲av高清不卡| 久久久亚洲精品成人影院| 丝袜脚勾引网站| 满18在线观看网站| 一边摸一边抽搐一进一出视频| 久久人妻熟女aⅴ| 一级毛片我不卡| 视频区欧美日本亚洲| 啦啦啦视频在线资源免费观看| 日本黄色日本黄色录像| 真人做人爱边吃奶动态| 老司机影院毛片| 丝袜美足系列| 精品国产乱码久久久久久男人| 一区二区三区乱码不卡18| 人人妻人人澡人人爽人人夜夜| 真人做人爱边吃奶动态| 国产av精品麻豆| 99九九在线精品视频| 国产高清不卡午夜福利| 高清欧美精品videossex| 叶爱在线成人免费视频播放| 欧美性长视频在线观看| av网站免费在线观看视频| 丝袜脚勾引网站| 电影成人av| 一级a爱视频在线免费观看| 男女下面插进去视频免费观看| 老司机亚洲免费影院| 在线看a的网站| 不卡av一区二区三区| 免费观看av网站的网址| 男人爽女人下面视频在线观看| 永久免费av网站大全| 下体分泌物呈黄色| 国产精品国产三级专区第一集| 国产日韩欧美视频二区| 亚洲天堂av无毛| 国产高清视频在线播放一区 | 大片免费播放器 马上看| 国产三级黄色录像| 精品久久久精品久久久| 久久99一区二区三区| 亚洲国产精品成人久久小说| 这个男人来自地球电影免费观看| 亚洲成人免费电影在线观看 | 久久久国产欧美日韩av| 色综合欧美亚洲国产小说| avwww免费| 一级毛片电影观看| 在线精品无人区一区二区三| 啦啦啦视频在线资源免费观看| videosex国产| 色播在线永久视频| 你懂的网址亚洲精品在线观看| 不卡av一区二区三区| cao死你这个sao货| 亚洲三区欧美一区| 日韩伦理黄色片| 自线自在国产av| 亚洲av男天堂| 成年动漫av网址| √禁漫天堂资源中文www| 久久久欧美国产精品| 欧美久久黑人一区二区| 女性生殖器流出的白浆| 亚洲欧美清纯卡通| 日韩人妻精品一区2区三区| 欧美成人午夜精品| 亚洲人成电影免费在线| 日本91视频免费播放| 777米奇影视久久| 久久久精品免费免费高清| 国产日韩一区二区三区精品不卡| 国产在线视频一区二区| 精品久久久精品久久久| 伊人亚洲综合成人网| 亚洲成人免费电影在线观看 | 欧美精品高潮呻吟av久久| 大片电影免费在线观看免费| 亚洲第一青青草原| 国产成人免费无遮挡视频| 免费高清在线观看日韩| 成年人午夜在线观看视频| 亚洲精品自拍成人| 国产高清视频在线播放一区 | 婷婷色综合www| 天天操日日干夜夜撸| 亚洲五月色婷婷综合| 国产一区二区三区av在线| 成年av动漫网址| 一级毛片 在线播放| 精品国产一区二区三区四区第35| 中国国产av一级| 19禁男女啪啪无遮挡网站| 欧美日韩国产mv在线观看视频| 男女免费视频国产| 国产精品一二三区在线看| 精品熟女少妇八av免费久了| 在线观看人妻少妇| 久久久久精品国产欧美久久久 | 国产伦人伦偷精品视频| 久久久久国产精品人妻一区二区| 久9热在线精品视频| av天堂在线播放| 久久久国产一区二区| kizo精华| 亚洲中文日韩欧美视频| 永久免费av网站大全| 建设人人有责人人尽责人人享有的| 国产一级毛片在线| 亚洲欧美成人综合另类久久久| 欧美av亚洲av综合av国产av| 国产高清国产精品国产三级| videosex国产| 久久久久网色| 制服诱惑二区| 国产亚洲av片在线观看秒播厂| 亚洲精品国产一区二区精华液| 欧美黄色片欧美黄色片| 精品国产国语对白av| 亚洲七黄色美女视频| 日韩一区二区三区影片| 国产成人av教育| 巨乳人妻的诱惑在线观看| 国产高清videossex| 精品少妇黑人巨大在线播放| 日韩免费高清中文字幕av| 日韩一本色道免费dvd| 悠悠久久av| 在线av久久热| 亚洲成人国产一区在线观看 | 欧美在线一区亚洲| 亚洲av综合色区一区| 国产麻豆69| 国产真人三级小视频在线观看| 精品第一国产精品| 最近中文字幕2019免费版| 亚洲午夜精品一区,二区,三区| 香蕉国产在线看| 男女午夜视频在线观看| 亚洲精品国产色婷婷电影| 一二三四社区在线视频社区8| 99九九在线精品视频| 一级片免费观看大全| 晚上一个人看的免费电影| 国产成人欧美在线观看 | 男男h啪啪无遮挡| 亚洲一码二码三码区别大吗| 亚洲精品国产区一区二| 少妇的丰满在线观看| 欧美老熟妇乱子伦牲交| 欧美日韩国产mv在线观看视频| 亚洲国产欧美一区二区综合| 中文字幕制服av| 婷婷丁香在线五月| 国产精品二区激情视频| 国产精品久久久久成人av| 超色免费av| 丝瓜视频免费看黄片| 美女扒开内裤让男人捅视频| 九草在线视频观看| 欧美激情 高清一区二区三区| 午夜福利视频精品| 国产精品 国内视频| 三上悠亚av全集在线观看| 午夜激情久久久久久久| 日韩,欧美,国产一区二区三区| 我要看黄色一级片免费的| 亚洲五月婷婷丁香| 在线观看人妻少妇| 久久人人爽av亚洲精品天堂| 肉色欧美久久久久久久蜜桃| 亚洲精品在线美女| 日韩中文字幕视频在线看片| 深夜精品福利| 99九九在线精品视频| 中国美女看黄片| 亚洲图色成人| 成人影院久久| 叶爱在线成人免费视频播放| 日韩中文字幕欧美一区二区 | 午夜福利一区二区在线看| 一区二区三区四区激情视频| 天天操日日干夜夜撸| 亚洲国产成人一精品久久久| 国产精品免费大片| 视频区欧美日本亚洲| 国产黄色视频一区二区在线观看| 精品国产一区二区久久| 赤兔流量卡办理| 王馨瑶露胸无遮挡在线观看| 日韩大码丰满熟妇| 国产野战对白在线观看| 黄片播放在线免费| 免费观看av网站的网址| 婷婷成人精品国产| 大片电影免费在线观看免费| 国产一区二区激情短视频 | 久久久欧美国产精品| netflix在线观看网站| 国产精品国产三级专区第一集| 国产日韩欧美亚洲二区| 亚洲情色 制服丝袜| 黄频高清免费视频| 丝袜在线中文字幕| 熟女少妇亚洲综合色aaa.| 国产亚洲欧美在线一区二区| 九色亚洲精品在线播放| 一区二区av电影网| 蜜桃在线观看..| 午夜福利影视在线免费观看| 大香蕉久久成人网| 天天影视国产精品| 脱女人内裤的视频| 一级片免费观看大全| 精品亚洲乱码少妇综合久久| 搡老乐熟女国产| 免费观看人在逋| 丝袜在线中文字幕| www.精华液| 老司机影院毛片| 天堂俺去俺来也www色官网| 最新在线观看一区二区三区 | av一本久久久久| 国产av精品麻豆| 亚洲av综合色区一区| 大片电影免费在线观看免费| 久久精品人人爽人人爽视色| 麻豆乱淫一区二区| 久久久久国产一级毛片高清牌| 欧美激情极品国产一区二区三区| 老司机在亚洲福利影院| 少妇人妻久久综合中文| 国产精品二区激情视频| 9热在线视频观看99| 亚洲伊人久久精品综合| 亚洲成av片中文字幕在线观看| 亚洲一区中文字幕在线| 国产麻豆69| 深夜精品福利| 下体分泌物呈黄色| 操美女的视频在线观看| 99国产精品免费福利视频| a 毛片基地| 亚洲av男天堂| 日本黄色日本黄色录像| 久久久久国产精品人妻一区二区| 成人国语在线视频| 黄色一级大片看看| 亚洲欧美激情在线| 大码成人一级视频| 亚洲成色77777| 另类精品久久| 亚洲国产中文字幕在线视频| 国产av精品麻豆| 飞空精品影院首页| 国产视频首页在线观看| 日本a在线网址| 999久久久国产精品视频| 久久影院123| 纵有疾风起免费观看全集完整版| av线在线观看网站| 国产精品一区二区免费欧美 | 亚洲精品国产色婷婷电影| 男女之事视频高清在线观看 | 国产精品一国产av| 80岁老熟妇乱子伦牲交| 中文字幕精品免费在线观看视频| 日韩伦理黄色片| 亚洲国产精品999| 老司机在亚洲福利影院| 日韩大片免费观看网站|