基于局部最小二乘支持向量機(jī)的潛空間廣義預(yù)測(cè)控制器

張壤文,田學(xué)民，王平

基于局部最小二乘支持向量機(jī)的潛空間廣義預(yù)測(cè)控制器

張壤文,田學(xué)民，王平

(中國(guó)石油大學(xué)(華東) 信息與控制工程學(xué)院，山東 青島 266580)

針對(duì)多變量、非線性、時(shí)變的實(shí)際工業(yè)過程系統(tǒng)，提出了一種基于局部最小二乘支持向量機(jī)的潛空間廣義預(yù)測(cè)控制方法。該方法通過偏最小二乘構(gòu)建潛變量空間，從而將復(fù)雜的多變量系統(tǒng)轉(zhuǎn)變成多個(gè)單變量子系統(tǒng)，然后在每個(gè)采樣點(diǎn)利用即時(shí)學(xué)習(xí)選擇相關(guān)數(shù)據(jù)樣本，在潛空間內(nèi)在線建立每個(gè)單變量子系統(tǒng)(SISO)的局部最小二乘支持向量機(jī)(LSSVM)模型，最后利用廣義預(yù)測(cè)控制器對(duì)這多個(gè)子系統(tǒng)分別實(shí)施預(yù)測(cè)控制。利用即時(shí)學(xué)習(xí)剔除冗余數(shù)據(jù)樣本，提升了LSSVM的魯棒性，并且使其更適用于實(shí)時(shí)建模和控制。利用該控制器對(duì)四容水箱對(duì)象進(jìn)行仿真研究，驗(yàn)證了算法的有效性。

多變量 非線性 偏最小二乘 即時(shí)學(xué)習(xí) 最小二乘支持向量機(jī) 廣義預(yù)測(cè)控制

1 概 述

模型預(yù)測(cè)控制[1]MPC(model predictive control)是目前過程工業(yè)中應(yīng)用最為廣泛也最成功的一種先進(jìn)控制算法。傳統(tǒng)的模型預(yù)測(cè)控制算法包括模型算法控制(MAC)[2]，動(dòng)態(tài)矩陣控制(DMC)[3]和廣義預(yù)測(cè)控制(GPC)[4]等，都是基于線性系統(tǒng)提出的一類控制算法,因而對(duì)于線性系統(tǒng)可以得到令人滿意的控制效果。然而，大量的實(shí)際工業(yè)過程中都存在著非線性，并且其操作點(diǎn)也會(huì)在一定的范圍內(nèi)發(fā)生改變，因而不能簡(jiǎn)單的用線性模型來模擬系統(tǒng)的真實(shí)工況，這也導(dǎo)致了MPC算法無法直接應(yīng)用于非線性系統(tǒng)。

目前大量針對(duì)于非線性系統(tǒng)的預(yù)測(cè)控制方法都是通過引入神經(jīng)網(wǎng)絡(luò)[5-6]、模糊模型[7-8]等經(jīng)驗(yàn)?zāi)Ｐ蛠肀硎鞠到y(tǒng)的非線性特性。但是神經(jīng)網(wǎng)絡(luò)算法是基于經(jīng)驗(yàn)風(fēng)險(xiǎn)最小化原則的,因而存在局部極小點(diǎn)、收斂速度慢以及過學(xué)習(xí)等不足，這些缺陷會(huì)嚴(yán)重影響神經(jīng)網(wǎng)絡(luò)在非線性MPC[9-10]中的應(yīng)用。而模糊模型作為預(yù)測(cè)模型在建模時(shí)往往需要足夠大的樣本容量，且模糊規(guī)則庫(kù)和控制區(qū)域的劃分往往比較復(fù)雜，這也限制了模糊預(yù)測(cè)控制[11-12]在工業(yè)過程中的應(yīng)用。支持向量機(jī)(SVM)[13-14]是一種基于結(jié)構(gòu)風(fēng)險(xiǎn)最小化原則的算法，不僅兼顧了學(xué)習(xí)算法的經(jīng)驗(yàn)風(fēng)險(xiǎn)和泛化能力,并且具有訓(xùn)練速度快等優(yōu)點(diǎn),但其缺點(diǎn)是當(dāng)訓(xùn)練樣本過多時(shí)則會(huì)導(dǎo)致維數(shù)災(zāi)難。最小二乘支持向量機(jī)(LSSVM)[15-16]是一類改進(jìn)的SVM，它采用等式約束代替不等式約束，通過求解一組等式方程得到參數(shù)的解析解，從而能夠克服用二次規(guī)劃方法求解SVM的維數(shù)災(zāi)難問題。但是，LSSVM模型的訓(xùn)練需要進(jìn)行矩陣求逆，當(dāng)累計(jì)數(shù)據(jù)量較大時(shí)，會(huì)導(dǎo)致計(jì)算量和所需內(nèi)存容量大幅增加，從而不適用于在線建模。同時(shí)，由于LSSVM的支持向量是沒有稀疏性的，如果不加區(qū)別的將所有樣本用于建模也會(huì)降低LSSVM的魯棒性。

為了在建模過程中剔除冗余樣本信息，提升LSSVM建模的魯棒性，同時(shí)減小計(jì)算量，利用即時(shí)學(xué)習(xí)(just-in-time learning)[17-18]對(duì)所有樣本進(jìn)行篩選，從中選取相關(guān)樣本集建立當(dāng)前時(shí)刻的局部LSSVM模型，并且隨著系統(tǒng)工作點(diǎn)的變化，在線滾動(dòng)建立系統(tǒng)的若干個(gè)局部模型，以滿足非線性系統(tǒng)實(shí)時(shí)建模和控制的需要。但是目前大量基于即時(shí)學(xué)習(xí)的控制方法[19-20]針對(duì)于多入多出(MIMO)系統(tǒng)的解決方法通常是將其轉(zhuǎn)化為幾個(gè)多入單出(MISO)子系統(tǒng)，然后建立MISO子系統(tǒng)的局部模型并設(shè)計(jì)控制器。

為了避免對(duì)MISO系統(tǒng)的建模同時(shí)也為了減小計(jì)算量，利用偏最小二乘(PLS)[21]建立潛變量空間，從而將復(fù)雜的多變量控制問題轉(zhuǎn)變?yōu)閱巫兞靠刂茊栴}，減小計(jì)算量的同時(shí)也簡(jiǎn)化了控制器。PLS通過構(gòu)造潛變量，可以將原始空間的MIMO 控制系統(tǒng)分解為潛空間內(nèi)多個(gè)的單輸入單輸出(SISO)子系統(tǒng)，進(jìn)而可以實(shí)現(xiàn)自動(dòng)解耦以及回路匹配的特性，簡(jiǎn)化了控制器的設(shè)計(jì)。在潛空間中，由于多個(gè)SISO系統(tǒng)是并行求解控制律，因而不會(huì)由于維數(shù)增加而加大計(jì)算耗時(shí)。然而，傳統(tǒng)的PLS方法僅僅考慮了系統(tǒng)輸入輸出之間的靜態(tài)關(guān)系，這使得它們不適用于系統(tǒng)的動(dòng)態(tài)建模。為了建立系統(tǒng)的動(dòng)態(tài)模型從而實(shí)施預(yù)測(cè)控制， Chi等利用將ARX模型嵌入到PLS內(nèi)模型中來刻畫過程動(dòng)態(tài)的方法設(shè)計(jì)了廣義預(yù)測(cè)控制器[22]，但是這類方法對(duì)于強(qiáng)非線性對(duì)象則不能獲得令人滿意的建模和控制結(jié)果。在此基礎(chǔ)上，Chi等提出了一種在PLS框架下的多模型預(yù)測(cè)控制策略[23]，將非線性過程看成分段線性的系統(tǒng)，通過一組線性模型來建模整個(gè)過程。但多模型方法存在難以確定子模型類型、個(gè)數(shù)和如何選擇合適的模型切換算法的問題。

為了處理多變量非線性系統(tǒng)的實(shí)時(shí)建模和控制問題，筆者提出了一種基于局部LSSVM的動(dòng)態(tài)PLS建模方法并在此基礎(chǔ)上設(shè)計(jì)了預(yù)測(cè)控制器。該方法利用PLS降維、解耦的同時(shí)建立潛變量空間，從而將MIMO非線性系統(tǒng)轉(zhuǎn)化為多個(gè)SISO子系統(tǒng)，然后利用潛空間內(nèi)過程的輸入輸出數(shù)據(jù)建立原始數(shù)據(jù)庫(kù)，在每個(gè)采樣時(shí)刻，根據(jù)當(dāng)前狀態(tài)，在數(shù)據(jù)庫(kù)中找到相關(guān)樣本集從而建立每個(gè)子系統(tǒng)的局部LSSVM模型。為了避免在系統(tǒng)運(yùn)行的每個(gè)采樣周期求解非線性規(guī)劃問題的復(fù)雜計(jì)算,在每個(gè)采樣時(shí)刻將非線性模型線性化,以解決非線性預(yù)測(cè)模型在線優(yōu)化的困難，最后利用廣義預(yù)測(cè)控制算法實(shí)現(xiàn)多步預(yù)測(cè)，求解控制策略?？紤]到實(shí)時(shí)性問題，為了加快在線建模速度，降低建模復(fù)雜度，將系統(tǒng)當(dāng)前時(shí)刻狀態(tài)和上一時(shí)刻狀態(tài)的相似度進(jìn)行比較，如果大于某一閾值，則繼續(xù)沿用上一時(shí)刻的模型。同時(shí)，提出了一種選取相似度閾值的自適應(yīng)方法，避免了閾值選取的盲目性，仿真實(shí)例驗(yàn)證了算法的有效性。

2 基于LSSVM的動(dòng)態(tài)PLS建模方法

2.1 偏最小二乘方法

在標(biāo)準(zhǔn)的PLS中，輸入數(shù)據(jù)為X=(xij)l×m，輸出數(shù)據(jù)為Y=(yij)l×n，l,m,n分別定義了樣本數(shù)、輸入變量和輸出變量的維數(shù)，外部模型可通過對(duì)輸入數(shù)據(jù)X和輸出數(shù)據(jù)Y進(jìn)行迭代提取主元來獲得：

(1)

(2)

式中：a——主元個(gè)數(shù)；T，U——得分矩陣；P，Q——相對(duì)應(yīng)的負(fù)荷矩陣；ti，ui——得分矩陣T和U的第i(0

則PLS回歸模型最終可以被表達(dá)為

(3)

式中：B=diag(b1,b2, …,ba)，bi為回歸系數(shù)。

實(shí)際上，PLS將多變量回歸問題轉(zhuǎn)變?yōu)橐幌盗袉巫兞炕貧w問題，在此基礎(chǔ)上，非方和病態(tài)問題就可以很容易地被解決。但是在傳統(tǒng)的PLS建模方法中僅僅考慮到了輸入輸出之間的穩(wěn)態(tài)關(guān)系，忽略了其動(dòng)態(tài)關(guān)系，所以在傳統(tǒng)PLS的基礎(chǔ)上進(jìn)行改進(jìn)使其體現(xiàn)系統(tǒng)的動(dòng)態(tài)特性具有重要意義。

2．2 最小二乘支持向量機(jī)

f(x)=wTφ(x)+b

(4)

來估計(jì)未知的非線性函數(shù)。則LSSVM的目標(biāo)函數(shù)[19]被定義為

(5)

s.t.yi=wTφ(xi)+b+ei

i=1, …,N

式中： 目標(biāo)函數(shù)的第一項(xiàng)對(duì)應(yīng)模型泛化能力；第二項(xiàng)代表模型的精確性；γ——是模型泛化能力和精度之間的1個(gè)折中參數(shù)，正常數(shù)；ei——第i個(gè)數(shù)據(jù)的實(shí)際輸出和預(yù)測(cè)輸出間的誤差。

為了解決該優(yōu)化問題，對(duì)LSSVM的目標(biāo)函數(shù)構(gòu)建Lagrange函數(shù)并分別對(duì)其求微分，最終得到LSSVM模型估計(jì)為

(6)

式中：αi——Lagrange乘子；K(,)——核函數(shù)，本文選用徑向基函數(shù)作為核函數(shù)：

(7)

2．3 即時(shí)學(xué)習(xí)方法

即時(shí)學(xué)習(xí)算法(JITL)將過程的輸入輸出數(shù)據(jù)儲(chǔ)存在數(shù)據(jù)庫(kù)中，當(dāng)新的過程數(shù)據(jù)到來時(shí)，在數(shù)據(jù)庫(kù)中找到與之相關(guān)的數(shù)據(jù)，并基于這些相關(guān)樣本數(shù)據(jù)建立局部模型。

假設(shè)即時(shí)學(xué)習(xí)的當(dāng)前數(shù)據(jù)庫(kù)包含了N個(gè)過程數(shù)據(jù)(y(i),x(i))i = 1 - N，令當(dāng)前數(shù)據(jù)為xq。為了建立非線性過程在當(dāng)前操作條件下的局部LSSVM模型，首先需通過當(dāng)前數(shù)據(jù)信息在數(shù)據(jù)庫(kù)中選擇相似的回歸向量，則同時(shí)包含距離和角度測(cè)量的相似性指標(biāo)si定義為

(8)

if cos(θi)≥0

式中：θi——Δxq與Δxi之間的角度，Δxq=xq-xq-1， Δxi=xi-xi-1；si——界定于0～1的值，當(dāng)si越接近于1，則表明xi與xq越相似。

當(dāng)計(jì)算完所有的si后，將其按降序排列，然后選取相似度較大的樣本數(shù)據(jù)建立當(dāng)前時(shí)刻的局部模型。

2．4 相似度閾值的自適應(yīng)更新

考慮到對(duì)每個(gè)xq作預(yù)測(cè)時(shí)都需要從數(shù)據(jù)庫(kù)中選擇相關(guān)數(shù)據(jù)來建模，當(dāng)數(shù)據(jù)庫(kù)龐大時(shí)，這種建模很耗時(shí)；而且當(dāng)系統(tǒng)的動(dòng)態(tài)特性變化緩慢或者是在某一區(qū)間穩(wěn)定工作時(shí)，模型變化不大，無需再?gòu)臄?shù)據(jù)庫(kù)中選擇相關(guān)數(shù)據(jù)來建模，可以沿用上一時(shí)刻的局部模型?；诖怂枷耄P者利用當(dāng)前時(shí)刻數(shù)據(jù)信息與上一時(shí)刻數(shù)據(jù)信息相似度比較的方法來降低在線建模的復(fù)雜度。利用下式計(jì)算當(dāng)前時(shí)刻的狀態(tài)信息和上一時(shí)刻狀態(tài)信息的相似度snew, new - 1：

(9)

如果snew, new - 1≥slim，則認(rèn)為當(dāng)前時(shí)刻的模型不變，繼續(xù)沿用上一時(shí)刻的模型，否則重新從數(shù)據(jù)庫(kù)中選擇相關(guān)數(shù)據(jù)來建模。這里slim為相似度閾值,這樣在保證模型精度的前提下，減小了建模過程的在線計(jì)算量，加快了系統(tǒng)的響應(yīng)速度。

為了獲得相似度閾值slim，通過計(jì)算當(dāng)前時(shí)刻的數(shù)據(jù)信息和所有建模樣本數(shù)據(jù)之間的相似度的最小值作為當(dāng)前時(shí)刻的slim。隨著時(shí)間的推移，相似度閾值slim也隨之實(shí)時(shí)更新，避免了人為選定閾值的盲目性。

2．5 基于局部LSSVM的動(dòng)態(tài)PLS建模方法

為了構(gòu)造動(dòng)態(tài)PLS模型，將局部LSSVM模型嵌入到PLS內(nèi)模型的動(dòng)態(tài)描述中，其動(dòng)態(tài)內(nèi)部模型為

ui=Hi(φi)

(10)

則整個(gè)動(dòng)態(tài)PLS模型為

(11)

(12)

式中：Hi(·)——在潛變量空間中第i個(gè)子系統(tǒng)具有LSSVM形式的動(dòng)態(tài)模型。

通過上述LSSVM建模方法，建立當(dāng)前時(shí)刻潛空間中輸入和輸出的動(dòng)態(tài)局部LSSVM模型為

(13)

式中：φi——回歸數(shù)據(jù)向量，φi(k)=[ui(k-1), …,ui(k-na),ti(k-1), …,ti(k-nb)]，其包含了第i個(gè)潛變量空間內(nèi)輸入輸出變量當(dāng)前和歷史的信息，且na和nb為潛變量子系統(tǒng)的階數(shù)?？梢姡瑸榱说玫接?xùn)練數(shù)據(jù)，首先需構(gòu)造回歸向量φi。

建立動(dòng)態(tài)PLS模型的框圖如圖1所示，其中，Hi(·)為在潛變量空間中第i個(gè)具有LSSVM形式的動(dòng)態(tài)局部模型。

圖1 動(dòng)態(tài)PLS模型示意

可見，原來的MIMO過程建模問題實(shí)際上就分解成了潛空間中多個(gè)SISO子模型的辨識(shí)，這使得過程的辨識(shí)和控制變得簡(jiǎn)單。

綜上所述，基于局部LSSVM的潛空間建模一般流程如圖2所示。

圖2 基于局部LSSVM潛空間建模設(shè)計(jì)流程

3 基于局部LSSVM的潛空間預(yù)測(cè)控制方法

3.1 模型實(shí)時(shí)線性化

為了避免在系統(tǒng)的每個(gè)采樣周期求解非線性規(guī)劃問題，可在每個(gè)采樣時(shí)刻將在線得到的被控系統(tǒng)局部LSSVM模型線性化。

采用u表示潛空間輸出，t表示潛空間輸入，φ表示潛空間回歸向量：

φk=[u(k-1), …,u(k-na),
t(k-1), …,t(k-nb)]

設(shè)當(dāng)前采樣時(shí)刻為k，把LSSVM模型在φk處泰勒展開，得到線性化模型：

P-α1φ(1)-…-αnaφ(na)+
β1φ(na+1)+…+βnbφ(na+nb)

(14)

式中：P——預(yù)測(cè)時(shí)域，且只與當(dāng)前采樣點(diǎn)有關(guān)，且

則被控系統(tǒng)當(dāng)前時(shí)刻的差分方程可以表示為

u(k)=P+β1t(k-1)+…+
βnbt(k-nb)-α1u(k-1)-…-

αnau(k-na)

(15)

即：

A(z-1)u(k)=B(z-1)t(k-1)+P

式中：A(z-1)=1+α1z-1+…+αnaz-na;B(z-1)=β1+β2z-1+…+βnbz-(nb-1)。

顯而易見，式(15)中的參數(shù)是通過把輸入映射到更高維的特征空間而得到的，因而具有更大的自由度和更高的推廣性能來描述訓(xùn)練數(shù)據(jù)。

3.2 廣義預(yù)測(cè)控制

在上述動(dòng)態(tài)建模的基礎(chǔ)上，將原始空間的MIMO系統(tǒng)轉(zhuǎn)化為潛空間里多個(gè)SISO子系統(tǒng)，則可對(duì)每一個(gè)SISO回路進(jìn)行預(yù)測(cè)控制[4]。

根據(jù)當(dāng)前時(shí)刻的差分方程，可將潛空間中每個(gè)回路的動(dòng)態(tài)內(nèi)模型表示為如下形式：

(16)

則預(yù)測(cè)模型[4]可以表示為

(17)

式中：Ej，F(xiàn)j——滿足Diophantin方程，并記Gj=EjB，可根據(jù)已知的輸入輸出信息及未來的輸入值，預(yù)測(cè)未來的輸出。

則相應(yīng)目標(biāo)函數(shù)可以被寫為

(18)

式中：u(k+j)——潛空間內(nèi)預(yù)測(cè)模型的輸出；usp(k+j)——潛空間內(nèi)參考軌跡；M——控制時(shí)域；λ——控制作用加權(quán)系數(shù)，表明了控制作用變化的限制；H——輸入變量的約束矩陣?？梢酝ㄟ^在每一時(shí)刻求解該二次規(guī)劃問題得到未來的輸入序列。

對(duì)于無約束問題，則通過最小化目標(biāo)函數(shù)，得到潛空間最優(yōu)控制序列為

Δt=(GTG+λI)-1GT(ysp-f)

(19)

即潛空間最優(yōu)控制量為

t(k)=t(k-1)+gT(usp-f)

(20)

式中：gT——矩陣(GTG+λI)-1GT的第1行。

基于局部LSSVM的潛空間預(yù)測(cè)控制器設(shè)計(jì)框圖如圖3所示。

圖3 局部LSSVM-LVMPC示意

綜上所述，基于局部LSSVM的潛空間預(yù)測(cè)控制器設(shè)計(jì)流程如下：

2) 在潛空間內(nèi)，利用輸入輸出潛變量建立a個(gè)子系統(tǒng)的初始數(shù)據(jù)庫(kù)。

3) 根據(jù)當(dāng)前狀態(tài)信息，利用即時(shí)學(xué)習(xí)選取數(shù)據(jù)庫(kù)中相似度較大的樣本數(shù)據(jù)，建立當(dāng)前時(shí)刻每個(gè)子系統(tǒng)的局部LSSVM模型，同時(shí)計(jì)算相似度閾值slim。

4) 在潛空間內(nèi)，對(duì)每個(gè)子模型分別實(shí)施GPC從而得到最優(yōu)輸入潛變量。

5) 將每個(gè)子系統(tǒng)的最優(yōu)輸入潛變量投影回原始空間，從而得到實(shí)際的最優(yōu)控制輸入并應(yīng)用于實(shí)際過程。

6) 得到實(shí)際輸出，投影到潛空間，得到潛空間輸入輸出，更新數(shù)據(jù)庫(kù)，同時(shí)計(jì)算當(dāng)前時(shí)刻與上一時(shí)刻的snew, new -1，若snew, new -1≥slim，轉(zhuǎn)至步驟4)，反之，轉(zhuǎn)至步驟3)。

4 仿真驗(yàn)證

四容水箱裝置示意如圖4所示。

由圖4可以看出，該裝置由2個(gè)水泵、2個(gè)分配閥和4個(gè)互連的儲(chǔ)水罐組成。p1，p2分別為加在2個(gè)水泵上的電壓；v1，v2分別為2個(gè)分配閥的流量分配系數(shù)，且v1，v2均在[0, 1]內(nèi)變化；h1，h2，h3和h4分別為4個(gè)水箱的液位高度。系統(tǒng)的輸入即控制變量為p1和p2；系統(tǒng)的輸出即被控變量為2個(gè)下水箱的液位h1和h2。當(dāng)水泵的輸入p1變化時(shí)，不僅會(huì)使得下水箱液位h1變化，還會(huì)通過分配閥使得上水箱h4變化從而影響下水箱液位h2，所以2個(gè)下水箱液位h1和h2呈相互耦合的關(guān)系。

圖4 四容水箱裝置示意

該系統(tǒng)的動(dòng)態(tài)平衡方程[24]為

(21)

(22)

(23)

(24)

p1 min≤p1≤p1 max,p2 min≤p2≤p2 max

式中：A1，A2，A3，A4——4個(gè)水箱的橫截面積；a1，a2，a3，a4——4個(gè)水箱流出管道的橫截面積。

利用隨機(jī)白噪聲作為輸入對(duì)系統(tǒng)進(jìn)行開環(huán)測(cè)試，選擇N=500建立數(shù)據(jù)庫(kù)樣本。LSSVM 的核函數(shù)選用徑向基函數(shù)，模型參數(shù)γ=100,σ2=65。通過多次仿真實(shí)驗(yàn)獲得局部LSSVM學(xué)習(xí)集的規(guī)模l=150。

為了驗(yàn)證本文算法的有效性，將本文方法與基于全局LSSVM模型的預(yù)測(cè)控制算法進(jìn)行比較，所得階躍跟蹤仿真曲線如圖5所示。

由圖5可以看出，對(duì)于該兩入兩出非線性對(duì)象，基于局部LSSVM建模的潛空間GPC控制器通過將其解耦為2個(gè)SISO子系統(tǒng)從而進(jìn)行實(shí)時(shí)建模與預(yù)測(cè)控制，該控制器使得每個(gè)輸出都能較快地跟蹤給定值的變化而變化，快速性好，超調(diào)量較小，且由于耦合作用引起的波動(dòng)較小，控制效果要優(yōu)于基于全局LSSVM 的預(yù)測(cè)控制算法。

圖5 階躍跟蹤比較曲線示意

在階躍跟蹤過程中，兩回路過程的參數(shù)變化曲線如圖6及圖7所示。

為了測(cè)試算法的魯棒性，減小流量分配系數(shù)v2，從而使得流入下水箱2的流量減小，流入上水箱3的流量增大，從而影響下水箱1的流量，因此模型發(fā)生改變，控制效果如圖8所示。

圖6 p1-t1回路參數(shù)變化曲線示意

圖7 p2-t2回路參數(shù)變化曲線示意

圖8 模型失配下局部LSSVM與全局LSSVM比較

在模型失配情況下，筆者比較了兩種算法的RMSE，見表1所列。

表1 局部LSSVM與全局LSSVM模型預(yù)測(cè)比較

圖8比較了兩種方法在t=50， 350，750 s時(shí)不斷減小v2所得到的控制曲線。由圖8可以看出，在模型失配情況下，基于全局LSSVM的預(yù)測(cè)控制器由于對(duì)模型參數(shù)是離線辨識(shí)的，無法適應(yīng)改變的過程條件，導(dǎo)致輸出不能較快地跟蹤給定值的變化而變化。而本文所提方法能夠通過在線實(shí)時(shí)更新模型參數(shù)，使得其適用于不斷變化的過程條件。從表1也可以清晰地看出本文所提方法處理模型失配的效果要優(yōu)于基于全局LSSVM的預(yù)測(cè)控制器。

5 結(jié)束語

筆者提出了一種基于局部LSSVM模型的潛空間預(yù)測(cè)控制方法，該方法首先利用PLS建立潛空間，將多變量控制問題轉(zhuǎn)化為單變量控制問題，然后根據(jù)當(dāng)前時(shí)刻的狀態(tài)信息利用JITL篩選出相關(guān)樣本數(shù)據(jù)，建立潛空間內(nèi)的局部LSSVM模型，最后利用GPC對(duì)每個(gè)子系統(tǒng)分別進(jìn)行多步預(yù)測(cè)，求解最優(yōu)控制策略。為了加快在線建模的速度，比較當(dāng)前時(shí)刻狀態(tài)與上一時(shí)刻狀態(tài)的相似度，當(dāng)模型變化不大時(shí)，繼續(xù)沿用上一時(shí)刻的局部模型。與基于全局LSSVM模型的潛空間預(yù)測(cè)控制方法的仿真比較驗(yàn)證了所提方法針對(duì)于非線性多變量系統(tǒng)具有更強(qiáng)的魯棒性，能夠更穩(wěn)定地實(shí)現(xiàn)非線性預(yù)測(cè)控制。仿真實(shí)例驗(yàn)證了算法的有效性。

[1] Camacho D E F, Bordons D C. Model Predictive Control[J]. Preprint, 1999, 127(02)： 6137-6142.

[2] Rouhani R, Mehra R K. Model Algorithmic Control (MAC); Basic Theoretical Properties[J]. Automatica, 1982, 18(04)： 401-414.

[3] Cutler C R, Ramakar B L. Dynamic Matrix Control — A Computer Control Algorithm[C]//Joint Automatic Control Conference. San Francisco： CA, 1980.

[4] Clarke D W, Mohtadi C, Tuffs P S. Generalized Predictive Control — Part I. The Basic Algorithm[J]. Automatica, 1987, 23(87)： 137-148.

[5] Ahmad Z, Mat Noor R, Zhang J. Multiple Neural Networks Modeling Techniques in Process Control： a Review[J]. Asia‐Pacific Journal of Chemical Engineering, 2009, 4(04)： 403-419.

[6] Zhang L, Suganthan P N. A Survey of Randomized Algorithms for Training Neural Networks[J]. Information Sciences, 2016, 364： 146-155.

[7] Lemos A, Caminhas W, Gomide F. Multivariable Gaussian Evolving Fuzzy Modeling System[J]. IEEE Transactions on Fuzzy Systems, 2011, 19(19)： 91-104.

[8] Andrade L P D, Espíndola R P, Pereira G C, et al. Fuzzy Modeling of Plankton Networks[J]. Ecological Modelling, 2016, 337： 149-155.

[9] Li S, Li Y. Nonlinear Model Predictive Control of an Intensified Continuous Reactor Using Neural Networks[C]//Control Conference (CCC). 2015 34th Chinese. IEEE, 2015.

[10] Vasicˇkaninov A, Bako?ov M. Control of a Heat Exchanger Using Neural Network Predictive Controller Combined With Auxiliary Fuzzy Controller[J]. Applied Thermal Engineering, 2015, 89： 1046-1053.

[11] Liu M, Tian X. A Fuzzy Predictive Control Based on T-S Model [J].Microcomputer Information, 2006, 16： 13-15.

[12] Wu X, Shen J, Li Y, et al. Fuzzy Modeling and Predictive Control of Superheater Steam Temperature for Power Plant[J]. ISA Transactions, 2015, 56： 241-251.

[13] Iplikci S. A Support Vector Machine Based Control Application to the Experimental Three-tank System[J]. Isa Transactions, 2010, 49(03)： 376-386.

[14] Feng K, Lu J, Chen J. Nonlinear Model Predictive Control Based on Support Vector Machine and Genetic Algorithm[J]. Chinese Journal of Chemical Engineering, 2015, 23(12)： 2048-2052.

[15] Samui P, Kothari D P. Utilization of a Least Square Support Vector Machine (LSSVM) for Slope Stability Analysis[J]. Scientia Iranica, 2011, 18(01)： 53-58.

[16] Huang L, Wang Z, Ji X. LS-SVM Generalized Predictive Control Based on PSO and Its Application of Fermentation Control[C]//Proceedings of the 2015 Chinese Intelligent Systems Conference. Springer Berlin Heidelberg, 2016： 605-613.

[17] Cheng C, Chiu M S. A New Data-based Methodology for Nonlinear Process Modeling[J]. Chemical Engineering Science, 2004, 59(13)： 2801-2810.

[18] Su Q L, Kano M, Chiu M S, et al. A New Strategy of Locality Enhancement for Justin-Time Learning Method[J]. Computer Aided Chemical Engineering, 2012, 31(31)： 1662-1666.

[19] Kansha Y, Chiu M S. Adaptive Generalized Predictive Control Based on JITL Technique[J]. Journal of Process Control, 2009, 19(07)： 1067-1072.

[20] Yang X, Jia L, Chiu M S. Adaptive Decentralized PID Controllers Design Using JITL Modeling Methodology[J]. Journal of Process Control, 2012, 22(09)： 1531-1542.

[21] Wold S, Sj?str?m M, Eriksson L. PLS-regression： a Basic Tool of Chemometrics[J]. Chemometrics and Intelligent Laboratory Systems, 2001, 58(02)： 109-130.

[22] Chi Q, Fei Z, Zhao Z, et al. A Model Predictive Control Approach With Relevant Identification in Dynamic PLS Framework[J]. Control Engineering Practice, 2014, 22： 181-193.

[23] Chi Q, Liang J. A Multiple Model Predictive Control Strategy in the PLS Framework[J]. Journal of Process Control, 2015, 25(25)： 129-141.

[24] Kumar E G, Mithunchakravarthi B, Dhivya N. Enhancement of PID Controller Performance for a Quadruple Tank Process with Minimum and Non-Minimum Phase Behaviors[J]. Procedia Technology, 2014, 14： 480-489.

A Generalized Predictive Controller Based on Local Least Squares Support Vector Machines in Latent Space

Zhang Rangwen，Tian Xuemin，Wang Ping

(College of Information and Control Engineering, China University of Petroleum(East China),Qingdao, 266580,China)

s: A generalized predictive control approach based on local least squares support vector machines(LSSVM) model in latent space is proposed for industrial process systems associating with multivariate, nonlinearity and time-varying characteristics. The method constructs latent variable space with partial least squares algorithm. Complicated multi-variable control system is converted into several single input single output subsystems(SISO). Relevant data samples at present are selected by Just-in-time learning(JITL) at every sampling instant. Local LSSVM model for each SISO is online constructed in latent space. Predictive control is implemented to these subsystems separately with generalized predictive controller. Redundant data samples are removed by JITL. Robustness of LSSVM is improved. It is much more applicable for real time modeling and controlling. The proposed predictive control is applied to a quadruple tank process forsimulation. The algorithm efficiency is verified.

multivariate; nonlinearity; partial least squares; just-in-time learning; least squares support vector machines; generalized predictive control

國(guó)家自然科學(xué)基金資助項(xiàng)目(61273160，61403418)。

張壤文(1992—),女，在讀碩士研究生，主要研究方向?yàn)轭A(yù)測(cè)控制及自適應(yīng)控制。

TP273

A

1007-7324(2017)02-0020-07

稿件收到日期： 2016-09-09，修改稿收到日期： 2016-12-19。