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

    約束非線性系統(tǒng)穩(wěn)定經(jīng)濟(jì)模型預(yù)測控制

    2016-12-17 08:23:50何德峰
    自動化學(xué)報 2016年11期
    關(guān)鍵詞:閉環(huán)約束穩(wěn)定性

    何德峰

    約束非線性系統(tǒng)穩(wěn)定經(jīng)濟(jì)模型預(yù)測控制

    何德峰1

    考慮約束非線性系統(tǒng),提出一種新的具有穩(wěn)定性保證的經(jīng)濟(jì)模型預(yù)測控制(Economic model predictive control, EMPC)策略.由于經(jīng)濟(jì)性能指標(biāo)的非凸性和非正定性,引入關(guān)于經(jīng)濟(jì)最優(yōu)平衡點的正定輔助函數(shù).利用輔助函數(shù)的最優(yōu)值函數(shù)定義原始EMPC優(yōu)化問題的穩(wěn)定性約束.應(yīng)用終端約束集、終端代價函數(shù)和局部控制器三要素,建立閉環(huán)系統(tǒng)關(guān)于經(jīng)濟(jì)最優(yōu)平衡點的漸近穩(wěn)定性和漸近平均性能.進(jìn)一步,結(jié)合多目標(biāo)理想點概念,將提出的控制策略用于多個經(jīng)濟(jì)性能指標(biāo)的優(yōu)化控制,得到穩(wěn)定多目標(biāo)EMPC策略.最后,以連續(xù)攪拌反應(yīng)器為例,比較仿真結(jié)果驗證本文策略的有效性.

    模型預(yù)測控制,非線性系統(tǒng),經(jīng)濟(jì)優(yōu)化,穩(wěn)定性,多目標(biāo)控制

    近年來,經(jīng)濟(jì)模型預(yù)測控制(Economic model predictive control,EMPC)在學(xué)術(shù)界和工業(yè)界受到了廣泛關(guān)注[1?15].EMPC的一個顯著優(yōu)點是將過程實時控制與經(jīng)濟(jì)性能優(yōu)化集成在一個最優(yōu)控制的框架內(nèi)設(shè)計[1?2].不同于傳統(tǒng)MPC(Model predictive control)需要預(yù)先設(shè)定目標(biāo)值為前提(故也稱為目標(biāo)跟蹤型MPC)[1],EMPC原則上并不需要這樣一個前提[2].由于這類MPC性能函數(shù)通常與過程的經(jīng)濟(jì)(環(huán)保)性能相關(guān),故統(tǒng)稱為經(jīng)濟(jì)MPC.但該“經(jīng)濟(jì)”并不特指某個經(jīng)濟(jì)性能,而泛指一類非正定和(或)非凸的任意性能函數(shù).盡管如此,經(jīng)濟(jì)MPC和傳統(tǒng)MPC的實施機(jī)理一樣,即都采用滾動時域方式實現(xiàn)系統(tǒng)的閉環(huán)狀態(tài)(或輸出)反饋控制[1?3],但直接優(yōu)化經(jīng)濟(jì)性能計算控制量給EMPC帶來更復(fù)雜的穩(wěn)定性問題[1,16].

    常規(guī)EMPC穩(wěn)定性需引入強(qiáng)對偶性或耗散性假設(shè),并附加罰函數(shù)將經(jīng)濟(jì)性能函數(shù)轉(zhuǎn)換為正定性能函數(shù),再結(jié)合終端等式或不等式約束,建立EMPC的遞推可行性和閉環(huán)穩(wěn)定性[4?8,16?17].由于系統(tǒng)與性能函數(shù)不一定滿足正則性(Regulation)條件[6],約束非線性系統(tǒng)通常不滿足對偶性與耗散性假設(shè)[4].進(jìn)一步,由于附加項通常是平衡點偏差的函數(shù),將改變原經(jīng)濟(jì)性能的最優(yōu)路徑,影響原始經(jīng)濟(jì)性能,特別是終端等式約束加重了EMPC在線優(yōu)化的計算量,同時會減小EMPC的吸引域.針對終端約束對經(jīng)濟(jì)性能的影響,文獻(xiàn)[18]應(yīng)用耗散性和能控性假設(shè)建立了無終端約束EMPC穩(wěn)定性,但要求預(yù)測時域足夠大,而這將增加EMPC在線優(yōu)化的計算量.基于能控性假設(shè)和切換控制思想,文獻(xiàn)[19]提出了基于Lyapunov函數(shù)的穩(wěn)定EMPC策略,但為保證閉環(huán)系統(tǒng)的穩(wěn)定性和漸近平均性能,需要在線求解三個非線性最優(yōu)控制問題[20],增加了EMPC運行的復(fù)雜性.現(xiàn)有研究結(jié)果表明, EMPC的經(jīng)濟(jì)性能與穩(wěn)定性存在一定的矛盾[7,16].如何在優(yōu)化經(jīng)濟(jì)性能的同時確保閉環(huán)系統(tǒng)的穩(wěn)定性是EMPC研究的一個關(guān)鍵問題.

    另一方面,在工業(yè)過程經(jīng)濟(jì)目標(biāo)優(yōu)化操作中,通常存在多個相互矛盾的性能指標(biāo)[21?22],且這些性能指標(biāo)具有不同的物理含義.因為加權(quán)函數(shù)法使用簡單,在控制實踐中,通常被用來近似處理多目標(biāo)優(yōu)化控制問題.但各目標(biāo)加權(quán)系數(shù)的選擇與整定需要通過大量的離線實驗完成,特別是在系統(tǒng)約束和目標(biāo)函數(shù)非凸的情況下,權(quán)值選擇與整定會異常困難[22?23].近年來,提出了一些新的多目標(biāo)優(yōu)化MPC設(shè)計,如基于穩(wěn)態(tài)目標(biāo)規(guī)劃計算的雙層多目標(biāo)MPC[23]、參數(shù)規(guī)劃多目標(biāo)線性MPC[24]、混合邏輯整數(shù)規(guī)劃多目標(biāo)MPC[25?27]、考慮性能指標(biāo)優(yōu)先級要求的多目標(biāo)MPC[28?30]以及基于多目標(biāo)理想點跟蹤的多目標(biāo)MPC[31?32]等.由于理想點跟蹤多目標(biāo)MPC策略不需要人為選擇各性能的加權(quán)系數(shù),在不需要系統(tǒng)先驗知識的條件下能自動處理各個性能指標(biāo)的沖突性,成為目前廣受關(guān)注的一種多目標(biāo)MPC策略.但在理想點跟蹤多目標(biāo)MPC中,由于性能指標(biāo)相互沖突,折中性能函數(shù)在穩(wěn)態(tài)平衡點不為零,即折中性能函數(shù)不是穩(wěn)態(tài)平衡點的正定函數(shù),對應(yīng)的控制器通常不穩(wěn)定,因此可看作是一類特殊的EMPC問題,可采用EMPC方法建立理想點跟蹤多目標(biāo)MPC的閉環(huán)穩(wěn)定性.

    本文考慮具有狀態(tài)和控制約束的非線性系統(tǒng),提出一種新的具有遞推可行性和穩(wěn)定性保證的EMPC策略.通過離線計算經(jīng)濟(jì)性能指標(biāo)的最優(yōu)平衡點,引入關(guān)于該點偏差的正定輔助函數(shù).利用輔助函數(shù)的最優(yōu)值函數(shù)定義原始EMPC優(yōu)化問題的一個穩(wěn)定性約束.再應(yīng)用終端約束集、終端代價函數(shù)和局部控制器等三要素[3],建立閉環(huán)系統(tǒng)關(guān)于經(jīng)濟(jì)最優(yōu)平衡點的漸近穩(wěn)定性和漸近平均性能.與現(xiàn)有穩(wěn)定EMPC策略相比,本文的一個主要創(chuàng)新點是EMPC的遞推可行性和漸近穩(wěn)定性與非線性系統(tǒng)的強(qiáng)對偶性、耗散性條件無關(guān),從而簡化EMPC策略的設(shè)計.進(jìn)一步,結(jié)合多目標(biāo)理想點法[31],將上述結(jié)果推廣至多個經(jīng)濟(jì)性能指標(biāo)的優(yōu)化控制問題,提出一種新的具有穩(wěn)定性保證的理想點跟蹤多目標(biāo)EMPC策略.最后以連續(xù)攪拌反應(yīng)器多目標(biāo)控制為例,對比仿真驗證本文結(jié)果的有效性.

    符號說明:集合R≥0和I≥0分別表示非負(fù)實數(shù)集和非負(fù)整數(shù)集.Ia:b={i∈I≥0:a≤i≤b},其中a∈I≥0和b∈I≥0.符號uk表示k∈I時刻的一個信號序列,系統(tǒng)在k∈I≥0時刻的狀態(tài)為系統(tǒng)在零時刻的狀態(tài).給定一個向量符號T表示的轉(zhuǎn)置,表示的p范數(shù),表示在第k時刻對

    未來第k+i時刻的預(yù)測變量.

    1 問題描述

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

    考慮經(jīng)濟(jì)性能函數(shù)Le:X×U→R.該函數(shù)關(guān)于連續(xù)且有界,但可能是非凸的且對于任意不一定是正定的.由文獻(xiàn)[16]可知,極小化這類性能函數(shù)的控制律不一定保證閉環(huán)系統(tǒng)漸近穩(wěn)定,甚至出現(xiàn)震蕩和周期響應(yīng).盡管周期響應(yīng)適用一些間歇工業(yè)過程,但大部分連續(xù)工業(yè)過程要求閉環(huán)系統(tǒng)是漸近穩(wěn)定的.本文目標(biāo)是設(shè)計一個經(jīng)濟(jì)最優(yōu)狀態(tài)反饋控制律使閉環(huán)系統(tǒng)能漸近穩(wěn)定于最優(yōu)經(jīng)濟(jì)平衡點同時滿足系統(tǒng)狀態(tài)和控制約束(2),以及漸近平均性能條件

    2 穩(wěn)定經(jīng)濟(jì)預(yù)測控制

    2.1 預(yù)測控制器

    令N∈I≥1為預(yù)測時域,定義k時刻的N步預(yù)測控制序列和對應(yīng)預(yù)測狀態(tài)序列給定緊集XT?X,如果uk滿足?i∈I0:N?1,以及和則該控制序列稱為系統(tǒng)(1)的一個可行預(yù)測控制序列.

    考慮系統(tǒng)(1)的一個可行預(yù)測控制序列uk及其對應(yīng)的預(yù)測狀態(tài)序列xk,定義如下目標(biāo)函數(shù):

    為計算EMPC控制律,在每個時刻求解如下有限時域最優(yōu)控制問題

    如果優(yōu)化問題(7)在k時刻是可行的,則根據(jù)滾動時域控制原理,定義經(jīng)濟(jì)預(yù)測控制律如下:

    本文將證明值函數(shù)(9)是系統(tǒng)(15)的Lyapunov函數(shù),建立閉環(huán)系統(tǒng)關(guān)于的漸近穩(wěn)定性結(jié)論.

    算法1.(EMPC算法)

    2.2 穩(wěn)定性與平均性能

    顯然,XN滿足XT?XN?X,且XN?XN+1.

    定理1.如果假設(shè)1和2成立,則對于任意給定α≥0,優(yōu)化問題(7)在XN內(nèi)具有遞推可行性,進(jìn)而XN是閉環(huán)系統(tǒng)(15)的一個不變集.

    進(jìn)一步考慮控制序列(16),可得:

    定理2.如果假設(shè)1和2成立,且優(yōu)化問題(7)在初始時刻存在可行解.給定α∈[0,1),則是閉環(huán)系統(tǒng)(15)在XN內(nèi)的漸近穩(wěn)定平衡點,且XN是閉環(huán)系統(tǒng)的一個吸引域.

    證明.任意給定α∈[0,1),由于優(yōu)化問題(7)在初始時刻是可行的.根據(jù)定理1可知,該優(yōu)化問題在任意k時刻存在可行解.

    考慮假設(shè)1和2,將式(19)代入不等式(20)并整理得:

    注3. 由定理2證明過程可知,閉環(huán)系統(tǒng)(15)的穩(wěn)定性是通過不等式(7e)建立的,即利用(Ea,XT,π)滿足假設(shè)1和2,使函數(shù)(9)成為閉環(huán)系統(tǒng)的Lyapunov函數(shù),因此式(7e)是一類基于Lyapunov函數(shù)的穩(wěn)定性約束.在常規(guī)的基于Lyapunov函數(shù)MPC策略[1,38]中,性能函數(shù)是關(guān)于平衡點偏差的正定函數(shù),故Lyapunov穩(wěn)定性約束是通過強(qiáng)制一個已知的正定函數(shù)的導(dǎo)數(shù)小于零構(gòu)造的,而本文采用MPC的三要素原理建立閉環(huán)系統(tǒng)的Lyapunov函數(shù),因此是兩種不同的穩(wěn)定性約束條件.其次,經(jīng)濟(jì)優(yōu)化問題(7)的最優(yōu)解并不能保證是優(yōu)化問題(8)的最優(yōu)解,而通常只是一個可行解;同理,盡管是優(yōu)化問題(8)的最優(yōu)解,但只是優(yōu)化問題(7)的一個可行解.盡管如此,結(jié)合定理1和MPC三要素原理可得,由最優(yōu)解定義的滾動時域控制器同樣能保證的漸近穩(wěn)定性,且XN是對應(yīng)閉環(huán)系統(tǒng)的一個吸引域.最后,控制序列(16)只是優(yōu)化問題(8)的一個可行序列,但不一定是優(yōu)化問題(7)的可行序列,因為不等式(17)使得序列(16)不一定滿足約束(7e).

    定理3.如果假設(shè)1和2成立,且優(yōu)化問題(7)在初始時刻存在可行解.如果給定α∈[0,1),則閉環(huán)系統(tǒng)(15)在吸引域XN內(nèi)滿足漸近平均性能(3).

    即任意充分小ε>0,總存在一個有限時間T?,當(dāng)T>T?時,式(24)成立.從而由極限理論可知,

    結(jié)合式(25),對式(22)兩邊同時取極限,可得:

    即閉環(huán)系統(tǒng)(15)在XN內(nèi)具有漸近平均性能. □

    注 4.閉環(huán)系統(tǒng)的漸近平均性能(26)是建立在無窮時間域上的結(jié)果,因此,在任意一個有限時間區(qū)間的平均性能值不一定小于穩(wěn)態(tài)經(jīng)濟(jì)性能

    3 一類多目標(biāo)經(jīng)濟(jì)預(yù)測控制

    上述穩(wěn)定EMPC策略可以推廣應(yīng)用到多經(jīng)濟(jì)目標(biāo)優(yōu)化控制問題.本節(jié)結(jié)合多目標(biāo)理想點概念[39],給出一種穩(wěn)定多目標(biāo)EMPC策略.

    3.1 多目標(biāo)理想點

    考慮l∈I≥2個連續(xù)且有界的性能函數(shù)Lj: X×U → R,j∈I1:l.令為函數(shù)向量.在多目標(biāo)優(yōu)化問題中,這些性能函數(shù)可能相互沖突,此時不存在唯一的最優(yōu)解同時極小化各個性能函數(shù)j∈I1:l.通常采用Pareto最優(yōu)解定義多目標(biāo)優(yōu)化問題的有效解,但Pareto最優(yōu)解通常不唯一(甚至無窮多個),并在目標(biāo)空間Rl形成一個Pareto面(Pareto front)[39].因此,文獻(xiàn)[28–30]等對具有優(yōu)先級排序的多目標(biāo)優(yōu)化控制問題,采用字典序多目標(biāo)優(yōu)化方法設(shè)計多目標(biāo)MPC控制器;而文獻(xiàn)[31–32]等對一類沖突但無排序要求的多目標(biāo)優(yōu)化控制問題,提出了理想點跟蹤多目標(biāo)MPC控制器設(shè)計等.

    其中,向量p范數(shù)‖·‖p定義了目標(biāo)空間中性能值到理想點L?的距離.則本節(jié)的目的是設(shè)計一個多目標(biāo)狀態(tài)反饋控制器使閉環(huán)系統(tǒng)在約束(2)作用下能漸近穩(wěn)定于Pareto解同時極小化一組性能函數(shù)為此,結(jié)合算法1和多目標(biāo)理想點法,設(shè)計一種理想點跟蹤多目標(biāo)優(yōu)化經(jīng)濟(jì)模型預(yù)測控制器.

    3.2 多目標(biāo)EMPC

    下面給出多目標(biāo)EMPC控制律(30)的實施步驟.

    算法2.(多目標(biāo)EMPC算法)

    1)系統(tǒng)在吸引域XN內(nèi)漸近穩(wěn)定于

    2)系統(tǒng)在XN內(nèi)滿足如下性能:

    證明.考慮性能函數(shù)(29),并利用定理1~3的結(jié)果直接得到定理4. □

    4 實例仿真

    考慮二元物系聚合過程連續(xù)攪拌釜反應(yīng)器[38],其動態(tài)聚合過程可描述為

    其中,CA、T和Tc分別為反應(yīng)物濃度(mol/L)、反應(yīng)器溫度(K)和冷卻劑溫度(K),模型參數(shù)為q=100 L/min,CAf=1 mol/L,UA=5×104J/min·K,ρ=1000 g/L,Cp=0.239 J/g·K, E/R=8750 K,k0=7.2×1010min?1,V=100 L,Tf=350 K,△H=?5×104J/mol.定義狀態(tài)變量和控制變量u=Tc,及狀態(tài)和控制約束

    我們用多目標(biāo)優(yōu)化控制結(jié)果驗證本文結(jié)果的有效性.

    在聚合反應(yīng)過程中,通常希望增強(qiáng)聚合反應(yīng)速率以提高生產(chǎn)效率,即極小化性能函數(shù)

    同時希望降低能源消耗,采用極小化性能函數(shù)

    表示降低能源消耗.因此,一個理想的聚合反應(yīng)過程控制器應(yīng)該同時滿足上述兩個性能要求.

    選擇向量2-范數(shù),定義折中經(jīng)濟(jì)性能函數(shù)

    可以驗證,三要素(40)滿足假設(shè)1和2.

    在仿真中,采用歐拉差分法離散系統(tǒng)(34),取采樣周期為0.1 min,預(yù)測步長N為7,仿真總步長為30.采用MatLab 2007的fmincon函數(shù)優(yōu)化計算最優(yōu)控制問題(7)和(8).下面,先驗證本文結(jié)果的正確性,再與文獻(xiàn)[6]等提出的EMPC策略、常規(guī)目標(biāo)跟蹤型MPC策略做比較以驗證本文結(jié)果的優(yōu)越性.

    4.1 仿真實驗一

    選取系統(tǒng)初始狀態(tài)(0.7 mol/L,320 K),運行算法2,仿真結(jié)果如圖1所示.其中,實線為α=0的仿真結(jié)果;虛線為α=0.3的仿真結(jié)果;點劃線為α=0.7的仿真結(jié)果;點線為α=0.99的仿真結(jié)果.需要指出的是:為了更清楚地顯示值函數(shù)在不同α取值時的差異,圖1右下子圖的時間坐標(biāo)截取為前20步.由圖1可知,對于系數(shù)α∈[0,1),值函數(shù)都是單調(diào)遞減的.因此,由算法2產(chǎn)生的多目標(biāo)EMPC及其閉環(huán)系統(tǒng)漸近穩(wěn)定于Pareto最優(yōu)穩(wěn)態(tài)解且在所有時刻滿足狀態(tài)和控制約束(35).

    盡管對系數(shù)α∈[0,1),閉環(huán)系統(tǒng)漸近穩(wěn)定,但不同值對應(yīng)的閉環(huán)系統(tǒng)的動態(tài)響應(yīng)和經(jīng)濟(jì)性能是不一樣的,圖2給出了上述取值時的性能實時曲線.分析圖1和2可知,α取值越小,值函數(shù)下降越快,閉環(huán)系統(tǒng)趨于最優(yōu)平衡點的時間越短,但經(jīng)濟(jì)性能可能會下降,如圖2左子圖實線所示.反之亦然,如對應(yīng)α=0.99的點線所示,盡管在30步仿真時間內(nèi),閉環(huán)系統(tǒng)還沒有收斂至平衡點,但它的性能指標(biāo)L1可以達(dá)到更小的值.因此,設(shè)計者可以通過調(diào)整α取值對閉環(huán)系統(tǒng)的收斂速度和經(jīng)濟(jì)性能進(jìn)行權(quán)衡,從而提高經(jīng)濟(jì)預(yù)測控制器應(yīng)用的靈活性.

    圖1 狀態(tài)軌跡、控制輸入和值函數(shù)曲線Fig.1 Profiles of states,control and value functions

    4.2 仿真實驗二

    為驗證本文控制策略的優(yōu)越性,對比文獻(xiàn)[6]等提出耗散EMPC策略(簡記為D-EMPC)和目標(biāo)跟蹤MPC策略(簡記為T-MPC).由于T-MPC策略的性能函數(shù)總是設(shè)定值偏差的正定函數(shù)[1?2],因此在該次實驗中,令正定函數(shù)(39)為T-MPC策略的性能函數(shù).考慮初始狀態(tài)(0.7 mol/L,320 K),圖3給出了本文策略和其余兩種策略在相同仿真環(huán)境下的控制結(jié)果.其中,實線表示本文策略(α=0.2)對應(yīng)結(jié)果;點線表示D-MPC策略對應(yīng)結(jié)果;虛線表示T-MPC策略對應(yīng)結(jié)果;右下子圖符號“*”表示性能理想點,“o”表示最優(yōu)平衡點對應(yīng)的經(jīng)濟(jì)性能.

    圖2 性能函數(shù)實時曲線Fig.2 Real-time profiles of cost functions

    圖3 狀態(tài)軌跡、控制輸入和性能函數(shù)相位曲線Fig.3 Profiles of states,control,and cost functions

    由圖3分析可知,本文策略和T-MPC策略對應(yīng)的閉環(huán)系統(tǒng)漸近收斂至而D-MPC策略對應(yīng)的閉環(huán)系統(tǒng)在暫態(tài)響應(yīng)過程中出現(xiàn)震蕩現(xiàn)象(如控制輸入u),同時存在穩(wěn)態(tài)誤差(如反應(yīng)器溫度x2),這主要是因為D-EMPC策略要求被控系統(tǒng)和經(jīng)濟(jì)性能函數(shù)滿足嚴(yán)格耗散性條件.但對于非線性系統(tǒng),嚴(yán)格耗散性條件很難成立[5],如本例中的聚合過程(34)和經(jīng)濟(jì)性能(38)不滿足嚴(yán)格耗散性條件,從而D-EMPC策略不能保證閉環(huán)系統(tǒng)收斂至對應(yīng)的性能相位曲線也不能收斂至最優(yōu)平衡點對應(yīng)的經(jīng)濟(jì)性能“o”,如圖3右下子圖點線所示.另一方面,盡管本文策略和T-MPC策略的閉環(huán)系統(tǒng)都收斂至但收斂路徑不一樣.如果以進(jìn)入穩(wěn)態(tài)解±5%為界,那么本文策略對應(yīng)閉環(huán)系統(tǒng)的收斂過渡時間小于T-MPC策略閉環(huán)系統(tǒng)的收斂過渡時間,從而加快閉環(huán)系統(tǒng)的過渡響應(yīng)過程.

    令A(yù)v{L1}和Av{L2}分別表示性能函數(shù)(36)和(37)的漸近平均性能,則進(jìn)一步比較三種策略取得的漸近平均性能.取足夠大的仿真時間步長100,圖4和表1分別給出了三種策略關(guān)于性能函數(shù)(36)和(37)的瞬時性能和漸近平均性能.圖4中實線為本文策略對應(yīng)的結(jié)果;點線為D-MPC策略對應(yīng)的結(jié)果;虛線表示T-MPC策略對應(yīng)結(jié)果.由圖4分析可知,對于本文策略和T-MPC策略取得的平均性能較優(yōu)于D-MPC策略;對于三種策略取得的平均性能較為接近.進(jìn)一步,以D-MPC策略取得的漸近平均性能值作為參考基準(zhǔn),對于本文策略和T-MPC策略分別增效1.52%和1.69%;對于本文策略和T-MPC策略分別增效0.34%和0.27%.綜合圖3和4和表1結(jié)果,驗證了本文策略在性能優(yōu)化方面的優(yōu)越性.

    圖4 性能函數(shù)的瞬時值曲線Fig.4 Real-time value profiles of cost functions

    表1 漸近平均性能Table 1 Asymptotic average performance

    5 結(jié)語

    本文針對約束非線性系統(tǒng),提出了一種新的穩(wěn)定經(jīng)濟(jì)MPC策略.通過引入關(guān)于經(jīng)濟(jì)最優(yōu)平衡點的輔助正定函數(shù),并利用輔助函數(shù)的最優(yōu)值函數(shù)定義原始EMPC優(yōu)化問題的一類穩(wěn)定性約束.應(yīng)用常規(guī)MPC的終端約束集、終端代價函數(shù)和局部控制器三要素方法,建立了閉環(huán)系統(tǒng)關(guān)于經(jīng)濟(jì)最優(yōu)平衡點漸近穩(wěn)定性和漸近平均性能的充分條件.結(jié)合多目標(biāo)理想點法,將提出EMPC策略推廣至多經(jīng)濟(jì)目標(biāo)函數(shù)優(yōu)化控制問題.建立了多目標(biāo)EMPC策略的漸近穩(wěn)定性充分條件.通過對連續(xù)攪拌反應(yīng)器多目標(biāo)優(yōu)化控制的對比研究,驗證了本文結(jié)果的有效性.

    1 Ellis M,Durand H,Christofides P D.A tutorial review of economic model predictive control methods.Journal of Process Control,2014,24(8):1156?1178

    2 Tran T,Ling K V,Maciejowski J M.Economic model predictive control—a review.In:Proceedings of the 31st International Symposium on Automation and Robotics in Construction.Sydney,Australia:IAARC,2014.35?42

    3 Mayne D Q.Model predictive control:recent developments and future promise.Automatica,2014,50(12):2967?2986

    4 Angeli D,Amrit R,Rawlings J B.On average performance and stability of economic model predictive control. IEEE Transactions on Automatic Control,2012,57(7): 1615?1626

    5 Diehl M,Amrit R,Rawlings J B.A Lyapunov function for economic optimizing model predictive control.IEEE Transactions on Automatic Control,2011,56(3):703?707

    6 Amrit R,Rawlings J B,Angeli D.Economic optimization using model predictive control with a terminal cost.Annual Reviews in Control,2011,35(2):178?186

    7 Fagiano L,Teel A R.Generalized terminal state constraint for model predictive control.Automatica,2013,49(9): 2622?2631

    8 Ferramosca A,Limon D,Camacho E F.Economic MPC fora changing economiccriterion forlinearsystems. IEEE Transactions on Automatic Control,2014,59(10): 2657?2667

    9 Huang R,Harinath E,Biegler L T.Lyapunov stability of economically oriented NMPC for cyclic processes.Journal of Process Control,2011,21(4):501?509

    10 Chu B,Duncan S,Papachristodoulou A,Hepburn C.Using economic model predictive control to design sustainable policies for mitigating climate change.In:Proceedings of the 51st IEEE Conference on Decision and Control.Maui, HI:IEEE,2012.406?411

    11 Hovgaard T G,Larsen L F S,J?rgensen J B.Flexible and cost efficient power consumption using economic MPC a supermarket refrigeration benchmark.In:Proceedings of the 50th IEEE Conference on Decision and Control and European Control Conference.Orlando,FL:IEEE,2011. 848?854

    12 Ma W J,Gupta V.Desynchronization of thermally-coupled first-order systems using economic model predictive control. In:Proceedings of the 51st IEEE Conference on Decision and Control.Maui,HI:IEEE,2012.278?283

    13 Lao L F,Ellis M,Christofides P D.Economic model predictive control of parabolic PDE systems:addressing state estimation and computational efficiency.Journal of Process Control,2014,24(4):448?462

    14 Lao L F,Ellis M,Christofides P D.Smart manufacturing: handling preventive actuator maintenance and economics using model predictive control.AIChE Journal,2014,60(6): 2179?2196

    15 Lao L F,Ellis M,Christofides P D.Economic model predictive control of transport-reaction processes.Industrial& Engineering Chemistry Research,2014,53(18):7382?7396

    17 Rawlings J B,Angeli D,Bates C N.Fundamentals of economic model predictive control.In:Proceedings of the 51st IEEE Conference on Decision and Control.Maui,HI:IEEE, 2012.3851?3861

    19 Heidarinejad M,Liu J F,Christofides P D.Economic model predictive control of nonlinear process systems using Lyapunov techniques.AIChE Journal,2012,58(3):855?870

    20 Ellis M,Christofides P D.On finite-time and infinite-time cost improvement of economic model predictive control for nonlinear systems.Automatica,2014,50(10):2561?2569

    21 Maree J P,Imsland L.On multi-objective economic predictive control of cyclic process operation.Journal of Process Control,2014,24(8):1328?1336

    22 Vallerio M,van Impe J,Logist F.Tuning of NMPC controllers via multi-objective optimisation.Computers and Chemical Engineering,2014,61(11):38?50

    23 Pan Hong-Guang,Gao Hai-Nan,Sun Yao,Zhang Ying, Ding Bao-Cang.The algorithm and software implementation for double-layered model predictive control based on multi-priority rank steady-state optimization.Acta Automatica Sinica,2014,40(3):405?414 (潘紅光,高海南,孫耀,張英,丁寶蒼.基于多優(yōu)先級穩(wěn)態(tài)優(yōu)化的雙層結(jié)構(gòu)預(yù)測控制算法及軟件實現(xiàn).自動化學(xué)報,2014,40(3): 405?414)

    24 Bemporad A,de la Pea D M.Multiobjective model predictive control.Automatica,2009,45(12):2823?2830

    25 Zou Tao,Wang Xi,Li Shao-Yuan.Multi-model predictive control for nonlinear systems based on mixed logic.Acta Automatica Sinica,2007,33(2):188?192 (鄒濤,王昕,李少遠(yuǎn).基于混合邏輯的非線性系統(tǒng)多模型預(yù)測控制.自動化學(xué)報,2007,33(2):188?192)

    26 Wang H Y,Huang D X,Jin Y H.A hybrid model predictive control for handling infeasibility and constraint prioritization.Chinese Journal of Chemical Engineering,2005,13(2): 211?217

    27 Tyler M L,Morari M.Propositional logic in control and monitoring problems.Automatica,1999,35(4):565?582

    28 He De-Feng,Song Xiu-Lan,Huang Hua.Stabilizing constrained model predictive control with prioritized multiobjectives.Control and Decision,2013,28(12):1831?1836 (何德峰,宋秀蘭,黃驊.優(yōu)先級多目標(biāo)穩(wěn)定化約束模型預(yù)測控制.控制與決策,2013,28(12):1831?1836)

    29 He D F,Wang L,Sun J.On stability of multiobjective NMPC with objective prioritization.Automatica,2015,57: 189?198

    30 Ocampo-Martinez C,Ingimundarson A,Puig V,Quevedo J.Objective prioritization using lexicographic minimizers for MPC of sewer networks.IEEE Transactions on Control Systems Technology,2008,16(1):113?121

    31 Zavala V M,Flores-Tlacuahuac A.Stability of multiobjective predictive control:a utopia-tracking approach.Automatica,2012,48(10):2627?2632

    32 He D F,Wang L,Yu L.Multi-objective nonlinear predictive control of process systems:a dual-mode tracking control approach.Journal of Process Control,2015,25(3):142?151

    34 Yu Shu-You,Chen Hong,Zhang Peng,Li Xue-Jun.An LMI optimization approach for enlarging the terminal region of MPC for nonlinear systems.Acta Automatica Sinica,2008, 34(7):798?804 (于樹友,陳虹,張鵬,李學(xué)軍.一種基于LMI的非線性模型預(yù)測控制終端域優(yōu)化方法.自動化學(xué)報,2008,34(7):798?804)

    35 He D F,Yu L,Song X L.Optimized-based stabilization of constrained nonlinear systems:a receding horizon approach. Asian Journal of Control,2014,16(6):1693?1701

    36 Xi Yu-Geng,Li De-Wei.Fundamental philosophy and status of qualitative synthesis of model predictive control.Acta Automatica Sinica,2008,34(10):1225?1234 (席裕庚,李德偉.預(yù)測控制定性綜合理論的基本思路和研究現(xiàn)狀.自動化學(xué)報,2008,34(10):1225?1234)

    37 He De-Feng,Ding Bao-Cang,Yu Shu-You.Review of fundamental properties and topics of model predictive control for nonlinear systems.Control Theory and Applications,2013, 30(3):273?287 (何德峰,丁寶蒼,于樹友.非線性系統(tǒng)模型預(yù)測控制若干基本特點與主題回顧.控制理論與應(yīng)用,2013,30(3):273?287)

    38 Mhaskar P,El-Farra N H,Christofides P D.Stabilization of nonlinear systems with state and control constraints using Lyapunov-based predictive control.System&Control Letters,2006,55(8):650?659

    39 Xu Jiu-Ping,Li Jun.Multiple Objective Decision Making Theory and Methods.Beijing:Tsinghua University Press, 2005. (徐玖平,李軍.多目標(biāo)決策的理論與方法.北京:清華大學(xué)出版社, 2005.)

    何德峰 浙江工業(yè)大學(xué)信息工程學(xué)院教授.2008年獲得中國科學(xué)技術(shù)大學(xué)自動化系博士學(xué)位.主要研究方向為魯棒與預(yù)測控制.

    E-mail:hdfzj@zjut.edu.cn

    (HE De-Feng Professor at the College of Information Engineering,Zhejiang University of Technology.He received his Ph.D.degree from the Department of Automation,University of Science and Technology of China in 2008.His research interest covers robust and predictive control.)

    Stabilizing Economic Model Predictive Control of Constrained Nonlinear Systems

    HE De-Feng1

    This paper proposes a novelty economic model predictive control(EMPC)scheme with guaranteed stability for constrained nonlinear systems.Some auxiliary positive-definite functions at economy optimal equilibrium points are introduced for economic performance functions that are generally non-convex and non-positive-definite.The optimal value function of the auxiliary function is used to define an adjustable stability constraint which is imposed on the original EMPC optimization problem.By the triplet of terminal constrained sets,terminal penalty functions and local controllers, it is proved that the closed-loop system is asymptotically stable at the economically optimal equilibrium point and has an asymptotic average performance.Moreover,the proposed EMPC scheme is extended to optimization control of multiple economic performance functions,together with the concept of multi-objective utopia points.Then a stabilizing multiobjective EMPC scheme is presented.Finally,the effectiveness of the obtained results is illustrated by a simulation example of a continuously stirred tank reactor.

    Model predictive control(MPC),nonlinear systems,economic optimization,stability,multi-objective control Citation He De-Feng.Stabilizing economic model predictive control of constrained nonlinear systems.Acta Automatica Sinica,2016,42(11):1680?1690

    何德峰.約束非線性系統(tǒng)穩(wěn)定經(jīng)濟(jì)模型預(yù)測控制.自動化學(xué)報,2016,42(11):1680?1690

    2015-04-20 錄用日期2016-06-30

    Manuscript received April 20,2015;accepted June 30,2016

    國家自然科學(xué)基金(61374111),浙江省公益技術(shù)應(yīng)用研究計劃項目(2015C31057)資助

    Supported by National Natural Science Foundation of China (61374111),Public Welfare Technology Application Research Project of Zhejiang Province(2015C31057)

    本文責(zé)任編委喬俊飛

    Recommended by Associate Editor QIAO Jun-Fei

    1.浙江工業(yè)大學(xué)信息工程學(xué)院杭州310023

    1.College of Information Engineering,Zhejiang University of Technology,Hangzhou 310023

    DOI 10.16383/j.aas.2016.c150234

    猜你喜歡
    閉環(huán)約束穩(wěn)定性
    “碳中和”約束下的路徑選擇
    約束離散KP方程族的完全Virasoro對稱
    非線性中立型變延遲微分方程的長時間穩(wěn)定性
    單周期控制下雙輸入Buck變換器閉環(huán)系統(tǒng)設(shè)計
    黑龍江電力(2017年1期)2017-05-17 04:25:05
    半動力系統(tǒng)中閉集的穩(wěn)定性和極限集映射的連續(xù)性
    雙閉環(huán)模糊控制在石化廢水處理中的研究
    最優(yōu)價格與回收努力激勵的閉環(huán)供應(yīng)鏈協(xié)調(diào)
    一種基于全閉環(huán)實時數(shù)字物理仿真的次同步振蕩阻尼控制
    適當(dāng)放手能讓孩子更好地自我約束
    人生十六七(2015年6期)2015-02-28 13:08:38
    模糊微分方程的一致穩(wěn)定性
    午夜福利在线观看免费完整高清在 | 国产精品嫩草影院av在线观看 | 宅男免费午夜| 黄色成人免费大全| 色视频www国产| 久久6这里有精品| 日本精品一区二区三区蜜桃| 国内少妇人妻偷人精品xxx网站| 国产精品久久久久久久电影 | 国产免费一级a男人的天堂| 91九色精品人成在线观看| 久久久精品大字幕| 别揉我奶头~嗯~啊~动态视频| 香蕉丝袜av| 国产精品久久久人人做人人爽| 91久久精品电影网| 级片在线观看| 中文字幕精品亚洲无线码一区| 欧美乱妇无乱码| 综合色av麻豆| 91av网一区二区| 国产三级在线视频| 国内精品久久久久精免费| 成人鲁丝片一二三区免费| 欧美zozozo另类| 欧美日韩中文字幕国产精品一区二区三区| 18禁国产床啪视频网站| 人人妻人人澡欧美一区二区| 亚洲av二区三区四区| 午夜免费男女啪啪视频观看 | 色综合亚洲欧美另类图片| 不卡一级毛片| 成人永久免费在线观看视频| 亚洲不卡免费看| 久久午夜亚洲精品久久| 欧美一区二区亚洲| 丰满人妻一区二区三区视频av | 成人国产一区最新在线观看| 亚洲精品一卡2卡三卡4卡5卡| 怎么达到女性高潮| 高清在线国产一区| 欧美日韩亚洲国产一区二区在线观看| bbb黄色大片| 最近最新中文字幕大全免费视频| 十八禁网站免费在线| 国产av在哪里看| 精品久久久久久,| 欧美在线黄色| 亚洲精品色激情综合| 国产精品久久久人人做人人爽| 久久精品人妻少妇| 日本成人三级电影网站| 国产精品野战在线观看| 麻豆成人午夜福利视频| 欧美激情在线99| 免费高清视频大片| eeuss影院久久| 神马国产精品三级电影在线观看| 亚洲成人免费电影在线观看| 91九色精品人成在线观看| 午夜日韩欧美国产| 成人一区二区视频在线观看| 色av中文字幕| 黄色片一级片一级黄色片| 久久久久久久精品吃奶| 变态另类丝袜制服| av国产免费在线观看| 国产高清视频在线播放一区| 婷婷精品国产亚洲av| netflix在线观看网站| 亚洲人成网站高清观看| 人妻久久中文字幕网| 久久久久性生活片| 久久6这里有精品| 淫秽高清视频在线观看| 88av欧美| 男女做爰动态图高潮gif福利片| 免费看a级黄色片| 国产视频一区二区在线看| 特级一级黄色大片| 啦啦啦韩国在线观看视频| 国产精品久久久久久久久免 | 精品一区二区三区av网在线观看| 男插女下体视频免费在线播放| 叶爱在线成人免费视频播放| 国产精品自产拍在线观看55亚洲| 国产精品一区二区三区四区久久| 久久欧美精品欧美久久欧美| 悠悠久久av| 亚洲成人免费电影在线观看| 网址你懂的国产日韩在线| 欧美日本视频| 91在线观看av| 欧美高清成人免费视频www| 中文字幕av在线有码专区| 精品一区二区三区av网在线观看| 老司机福利观看| 亚洲一区高清亚洲精品| 97超级碰碰碰精品色视频在线观看| av女优亚洲男人天堂| 1024手机看黄色片| 亚洲五月婷婷丁香| 啪啪无遮挡十八禁网站| 69人妻影院| 91在线精品国自产拍蜜月 | 国产视频内射| 成人性生交大片免费视频hd| 青草久久国产| 桃红色精品国产亚洲av| 天天躁日日操中文字幕| 久久久久九九精品影院| 久久久久九九精品影院| 桃红色精品国产亚洲av| 天天添夜夜摸| 亚洲精品成人久久久久久| 大型黄色视频在线免费观看| 99久久久亚洲精品蜜臀av| 国产黄片美女视频| 又紧又爽又黄一区二区| 国产一级毛片七仙女欲春2| 婷婷精品国产亚洲av在线| 午夜福利在线观看免费完整高清在 | 99riav亚洲国产免费| 国产高清有码在线观看视频| 听说在线观看完整版免费高清| 国产三级黄色录像| 久久精品综合一区二区三区| a级毛片a级免费在线| 淫秽高清视频在线观看| 90打野战视频偷拍视频| 久久精品影院6| 亚洲无线在线观看| 久久久国产精品麻豆| 久久久精品欧美日韩精品| 九九热线精品视视频播放| 最近视频中文字幕2019在线8| 美女高潮喷水抽搐中文字幕| 国产精品精品国产色婷婷| 精品久久久久久久毛片微露脸| 91久久精品国产一区二区成人 | 欧美成人一区二区免费高清观看| 亚洲aⅴ乱码一区二区在线播放| 国模一区二区三区四区视频| 无人区码免费观看不卡| 日本在线视频免费播放| 丁香六月欧美| 又黄又爽又免费观看的视频| 国产三级中文精品| 一a级毛片在线观看| 九九热线精品视视频播放| 亚洲狠狠婷婷综合久久图片| 日本免费a在线| 日韩欧美精品免费久久 | 一边摸一边抽搐一进一小说| 在线看三级毛片| 成年免费大片在线观看| 免费人成在线观看视频色| 一级黄色大片毛片| 欧美av亚洲av综合av国产av| 男女做爰动态图高潮gif福利片| 69人妻影院| 国产97色在线日韩免费| bbb黄色大片| 成人午夜高清在线视频| xxxwww97欧美| 一本综合久久免费| 啦啦啦免费观看视频1| 亚洲av第一区精品v没综合| 中出人妻视频一区二区| 在线国产一区二区在线| 神马国产精品三级电影在线观看| 欧美性猛交╳xxx乱大交人| 免费在线观看亚洲国产| 人人妻人人看人人澡| 成人国产一区最新在线观看| 国产欧美日韩精品亚洲av| 一个人观看的视频www高清免费观看| 欧美一级毛片孕妇| 麻豆一二三区av精品| 日本黄大片高清| av在线天堂中文字幕| 国产精品嫩草影院av在线观看 | 中文亚洲av片在线观看爽| 黄片小视频在线播放| 午夜免费激情av| 亚洲不卡免费看| 免费av不卡在线播放| 国产真实伦视频高清在线观看 | 亚洲国产精品合色在线| avwww免费| 欧美在线一区亚洲| 亚洲男人的天堂狠狠| 欧美日韩精品网址| www.色视频.com| 亚洲av免费高清在线观看| 中文字幕久久专区| 日韩欧美 国产精品| 有码 亚洲区| 女警被强在线播放| 少妇的逼水好多| 精品人妻1区二区| 成人特级av手机在线观看| 久久久色成人| 精品久久久久久久人妻蜜臀av| 亚洲精品亚洲一区二区| 亚洲内射少妇av| 熟女人妻精品中文字幕| 国产熟女xx| 国产真实伦视频高清在线观看 | 成人国产综合亚洲| 91字幕亚洲| 极品教师在线免费播放| 少妇人妻精品综合一区二区 | 亚洲内射少妇av| x7x7x7水蜜桃| 搞女人的毛片| 国产一区在线观看成人免费| 91久久精品国产一区二区成人 | 精品福利观看| 国产不卡一卡二| 免费人成视频x8x8入口观看| 亚洲久久久久久中文字幕| 国产精品一区二区三区四区免费观看 | 欧美成人免费av一区二区三区| 精品99又大又爽又粗少妇毛片 | 叶爱在线成人免费视频播放| 欧美色欧美亚洲另类二区| 国产成年人精品一区二区| 18禁美女被吸乳视频| 久久草成人影院| 午夜日韩欧美国产| 精品电影一区二区在线| 亚洲精品粉嫩美女一区| 久9热在线精品视频| 国产成人a区在线观看| 国产黄片美女视频| 国产野战对白在线观看| 国产成年人精品一区二区| 天天躁日日操中文字幕| 亚洲人成网站在线播| 啦啦啦韩国在线观看视频| 真人做人爱边吃奶动态| 熟女电影av网| 国产真实伦视频高清在线观看 | 亚洲国产欧美人成| 国产成人aa在线观看| 成年人黄色毛片网站| 成人亚洲精品av一区二区| 少妇的逼水好多| 亚洲中文字幕一区二区三区有码在线看| 一a级毛片在线观看| 岛国在线免费视频观看| 精品一区二区三区视频在线 | 久久精品国产综合久久久| 波野结衣二区三区在线 | 国内揄拍国产精品人妻在线| 亚洲熟妇熟女久久| 亚洲av第一区精品v没综合| 免费观看的影片在线观看| 最新在线观看一区二区三区| 国产色爽女视频免费观看| 一本综合久久免费| 亚洲电影在线观看av| 欧美在线黄色| 成年女人永久免费观看视频| 中亚洲国语对白在线视频| 午夜精品久久久久久毛片777| 亚洲第一欧美日韩一区二区三区| 亚洲欧美精品综合久久99| 亚洲成a人片在线一区二区| 亚洲国产精品久久男人天堂| 久久精品国产综合久久久| 麻豆一二三区av精品| 久久精品国产综合久久久| 搡老熟女国产l中国老女人| 内射极品少妇av片p| 精品国产超薄肉色丝袜足j| 精品久久久久久久人妻蜜臀av| 搞女人的毛片| 国产免费一级a男人的天堂| 亚洲av免费在线观看| 日本免费a在线| av女优亚洲男人天堂| 免费在线观看成人毛片| 九色国产91popny在线| 国产精品一区二区三区四区免费观看 | 少妇人妻精品综合一区二区 | 亚洲va日本ⅴa欧美va伊人久久| 三级国产精品欧美在线观看| 国产精品一及| 最近最新免费中文字幕在线| 亚洲男人的天堂狠狠| 无限看片的www在线观看| 午夜福利高清视频| 亚洲片人在线观看| 国产高清有码在线观看视频| 狂野欧美白嫩少妇大欣赏| 亚洲国产中文字幕在线视频| 成人永久免费在线观看视频| 久久久久久九九精品二区国产| 亚洲av免费在线观看| 美女cb高潮喷水在线观看| 国产成人啪精品午夜网站| 国内精品一区二区在线观看| 亚洲精品一区av在线观看| 亚洲av免费在线观看| 欧美激情在线99| 亚洲欧美精品综合久久99| 久久精品亚洲精品国产色婷小说| 麻豆久久精品国产亚洲av| 丰满的人妻完整版| 亚洲精品影视一区二区三区av| 校园春色视频在线观看| 欧美+日韩+精品| 人妻夜夜爽99麻豆av| 午夜福利欧美成人| 欧美日韩中文字幕国产精品一区二区三区| 最近在线观看免费完整版| 国内精品久久久久精免费| 一级毛片女人18水好多| 中文亚洲av片在线观看爽| avwww免费| 黄色视频,在线免费观看| 国产欧美日韩一区二区精品| 日韩大尺度精品在线看网址| 久久久久久人人人人人| 露出奶头的视频| 91麻豆av在线| 国产黄片美女视频| 又粗又爽又猛毛片免费看| 国产欧美日韩一区二区三| 久久久久国内视频| 黄色丝袜av网址大全| av福利片在线观看| 男女下面进入的视频免费午夜| 亚洲av一区综合| 成年免费大片在线观看| 在线国产一区二区在线| 老汉色∧v一级毛片| 90打野战视频偷拍视频| 日本黄色片子视频| 国产av麻豆久久久久久久| 日韩精品青青久久久久久| 99国产综合亚洲精品| 色视频www国产| 欧美精品啪啪一区二区三区| 国产真人三级小视频在线观看| 叶爱在线成人免费视频播放| 97超级碰碰碰精品色视频在线观看| 欧美乱码精品一区二区三区| 日本与韩国留学比较| 一夜夜www| 精品久久久久久成人av| 变态另类成人亚洲欧美熟女| 久久久久久久久大av| 亚洲最大成人手机在线| 亚洲在线观看片| 91久久精品电影网| 国产精品女同一区二区软件 | 亚洲无线观看免费| 亚洲成人精品中文字幕电影| 欧美zozozo另类| 欧美大码av| 欧美乱色亚洲激情| 亚洲精品成人久久久久久| 精品久久久久久久毛片微露脸| 国产成人av教育| 国产成人a区在线观看| 亚洲国产高清在线一区二区三| 中文字幕人妻丝袜一区二区| 给我免费播放毛片高清在线观看| 日日夜夜操网爽| 99国产极品粉嫩在线观看| a级毛片a级免费在线| 美女 人体艺术 gogo| 真人做人爱边吃奶动态| 757午夜福利合集在线观看| 日韩人妻高清精品专区| 国产精品三级大全| 非洲黑人性xxxx精品又粗又长| 精品久久久久久久人妻蜜臀av| 91久久精品电影网| 亚洲精品国产精品久久久不卡| 精品欧美国产一区二区三| 国产色爽女视频免费观看| 精品国产超薄肉色丝袜足j| 99精品欧美一区二区三区四区| 午夜精品在线福利| 国内精品久久久久精免费| 一区二区三区高清视频在线| 国产国拍精品亚洲av在线观看 | 国产成人欧美在线观看| 少妇高潮的动态图| 国产伦精品一区二区三区视频9 | 亚洲国产日韩欧美精品在线观看 | 18禁裸乳无遮挡免费网站照片| 制服人妻中文乱码| 国产精品亚洲一级av第二区| 亚洲熟妇熟女久久| 色视频www国产| 午夜两性在线视频| 国产一区二区亚洲精品在线观看| aaaaa片日本免费| 久久久久国产精品人妻aⅴ院| 18禁国产床啪视频网站| 真人一进一出gif抽搐免费| 少妇的逼水好多| 国产精品av视频在线免费观看| 18禁裸乳无遮挡免费网站照片| 99久久99久久久精品蜜桃| 亚洲精品日韩av片在线观看 | 两人在一起打扑克的视频| 成人国产综合亚洲| 久久精品影院6| 又爽又黄无遮挡网站| 国产免费一级a男人的天堂| 国产一区二区激情短视频| 人妻久久中文字幕网| 999久久久精品免费观看国产| 乱人视频在线观看| 尤物成人国产欧美一区二区三区| 日韩免费av在线播放| 岛国在线免费视频观看| 午夜精品一区二区三区免费看| 此物有八面人人有两片| АⅤ资源中文在线天堂| 久久久久性生活片| 三级男女做爰猛烈吃奶摸视频| 黄色视频,在线免费观看| 婷婷丁香在线五月| 国产精品综合久久久久久久免费| 69av精品久久久久久| 免费观看精品视频网站| 五月伊人婷婷丁香| 给我免费播放毛片高清在线观看| 欧美日本视频| 国产欧美日韩一区二区三| 国产97色在线日韩免费| 听说在线观看完整版免费高清| 中国美女看黄片| 午夜视频国产福利| 国产av在哪里看| 日日干狠狠操夜夜爽| 级片在线观看| 国内少妇人妻偷人精品xxx网站| 老鸭窝网址在线观看| av片东京热男人的天堂| netflix在线观看网站| 一区福利在线观看| 久久久国产成人免费| 老司机深夜福利视频在线观看| 久久精品国产清高在天天线| 亚洲人成电影免费在线| 国产亚洲精品av在线| 黄色女人牲交| 热99re8久久精品国产| 亚洲成av人片在线播放无| 精品人妻1区二区| 国内毛片毛片毛片毛片毛片| 久久精品国产99精品国产亚洲性色| 天天一区二区日本电影三级| 亚洲av免费在线观看| 成人亚洲精品av一区二区| av在线蜜桃| 性色av乱码一区二区三区2| 美女被艹到高潮喷水动态| 男女做爰动态图高潮gif福利片| 国产精品嫩草影院av在线观看 | 99久久九九国产精品国产免费| 变态另类丝袜制服| 日本三级黄在线观看| av国产免费在线观看| 欧美日韩福利视频一区二区| 午夜两性在线视频| av视频在线观看入口| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产国拍精品亚洲av在线观看 | 搡老熟女国产l中国老女人| 久久精品国产清高在天天线| 校园春色视频在线观看| 午夜福利欧美成人| 免费看日本二区| 亚洲精品一区av在线观看| 成人av在线播放网站| 女同久久另类99精品国产91| 久久天躁狠狠躁夜夜2o2o| 国内精品一区二区在线观看| 免费大片18禁| 欧美乱码精品一区二区三区| 嫩草影院精品99| 亚洲一区二区三区不卡视频| 偷拍熟女少妇极品色| 国产精品99久久久久久久久| 舔av片在线| 欧美乱色亚洲激情| 深爱激情五月婷婷| 91九色精品人成在线观看| 亚洲av一区综合| 国产探花极品一区二区| 亚洲av熟女| 欧美日韩精品网址| 色精品久久人妻99蜜桃| h日本视频在线播放| 久久久久久久精品吃奶| 国产伦在线观看视频一区| 国内毛片毛片毛片毛片毛片| 内射极品少妇av片p| 18禁美女被吸乳视频| 757午夜福利合集在线观看| 两个人视频免费观看高清| 老熟妇乱子伦视频在线观看| 免费看十八禁软件| 亚洲中文字幕日韩| 男人和女人高潮做爰伦理| 我的老师免费观看完整版| 麻豆国产av国片精品| 亚洲中文字幕一区二区三区有码在线看| 亚洲成av人片免费观看| 亚洲中文字幕日韩| 久久性视频一级片| 日韩欧美 国产精品| 非洲黑人性xxxx精品又粗又长| 人人妻,人人澡人人爽秒播| 亚洲人成伊人成综合网2020| 人人妻人人澡欧美一区二区| 两人在一起打扑克的视频| 欧美大码av| 午夜免费成人在线视频| x7x7x7水蜜桃| 亚洲美女黄片视频| 老司机午夜福利在线观看视频| 欧美在线黄色| 亚洲国产精品久久男人天堂| 男人舔奶头视频| 草草在线视频免费看| 久久精品国产99精品国产亚洲性色| 亚洲精品日韩av片在线观看 | 日日干狠狠操夜夜爽| 亚洲国产日韩欧美精品在线观看 | 久久久精品大字幕| av国产免费在线观看| 两人在一起打扑克的视频| 12—13女人毛片做爰片一| 亚洲欧美日韩高清专用| 日韩精品中文字幕看吧| 最近视频中文字幕2019在线8| 国产在视频线在精品| eeuss影院久久| 免费看日本二区| 国产三级在线视频| 午夜福利视频1000在线观看| 少妇熟女aⅴ在线视频| 欧美性猛交黑人性爽| 日日摸夜夜添夜夜添小说| 一进一出好大好爽视频| 亚洲av不卡在线观看| 色在线成人网| 亚洲欧美激情综合另类| 我的老师免费观看完整版| 久久午夜亚洲精品久久| 日韩中文字幕欧美一区二区| 成人无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 欧美成人性av电影在线观看| 人人妻人人看人人澡| 久久久久久国产a免费观看| 日日夜夜操网爽| 国产国拍精品亚洲av在线观看 | 伊人久久精品亚洲午夜| or卡值多少钱| 在线观看舔阴道视频| 免费观看人在逋| 欧美乱码精品一区二区三区| 国产精品精品国产色婷婷| 久久久久久久久中文| 国产亚洲欧美在线一区二区| 可以在线观看毛片的网站| 国产欧美日韩精品亚洲av| 国产 一区 欧美 日韩| 国产精品久久久人人做人人爽| 亚洲久久久久久中文字幕| 亚洲精品日韩av片在线观看 | 欧美另类亚洲清纯唯美| 美女大奶头视频| 国产亚洲欧美98| 日韩欧美国产一区二区入口| 在线观看日韩欧美| 每晚都被弄得嗷嗷叫到高潮| 成年人黄色毛片网站| 99视频精品全部免费 在线| 色精品久久人妻99蜜桃| 婷婷精品国产亚洲av| 亚洲五月天丁香| 成人特级黄色片久久久久久久| 精品一区二区三区视频在线 | 欧美乱妇无乱码| 1000部很黄的大片| 一本一本综合久久| 老汉色∧v一级毛片| 最近在线观看免费完整版| 亚洲一区二区三区色噜噜| 成年人黄色毛片网站| 身体一侧抽搐| 麻豆国产av国片精品| 两个人的视频大全免费| 99精品久久久久人妻精品| 欧美一区二区精品小视频在线| 亚洲成人久久爱视频| 成人国产一区最新在线观看| 国产成人aa在线观看| 91字幕亚洲| 国产中年淑女户外野战色| 日本五十路高清| 亚洲精品456在线播放app | 一进一出抽搐动态|