自適應(yīng)粒子群優(yōu)化算法在聚丙烯熔融指數(shù)預(yù)報(bào)上的應(yīng)用

1、自適應(yīng)粒子群優(yōu)化算法在聚丙烯熔融指數(shù)預(yù)報(bào)上的應(yīng)用基金項(xiàng)目:  國家自然科學(xué)基金資助項(xiàng)目（編號(hào)：50876093）、浙江省科技廳國際合作項(xiàng)目（編號(hào)：2009C34008）、國家863計(jì)劃項(xiàng)目（編號(hào)：2006AA05Z226） Supported by NSFC of P.R. China (No. 50876093); International Cooperation and Exchange Project of Science and Technology Department of Zhejiang Province (No. 2009C34008); Nationa

2、l HI-TECH Research & Development Program of P.R.China (863 Program, No. 2006AA05Z226)趙成業(yè) 劉興高（工業(yè)控制技術(shù)國家重點(diǎn)實(shí)驗(yàn)室 浙江大學(xué)控制系 浙江 杭州 310027）（E-mail: ）摘 要：針對(duì)丙烯聚合生產(chǎn)控制中聚丙烯熔融指數(shù)在線測(cè)量的控制要求，以及過程變量間相關(guān)性高的特點(diǎn)，提出一種基于自適應(yīng)粒子群優(yōu)化算法和徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)的聚丙烯熔融指數(shù)預(yù)報(bào)新方法。該方法采用變參數(shù)的自適應(yīng)粒子群優(yōu)化算法提高優(yōu)化算法的效率和收斂性，融合了主成分分析、統(tǒng)計(jì)建模以及智能優(yōu)化方法，降低了預(yù)報(bào)模型的復(fù)雜度。提出了一

3、種基于徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)的統(tǒng)計(jì)預(yù)報(bào)模型的參數(shù)優(yōu)化和結(jié)構(gòu)優(yōu)化方法。使用該統(tǒng)計(jì)模型對(duì)工廠實(shí)際生產(chǎn)過程進(jìn)行預(yù)報(bào)，并與國內(nèi)外相關(guān)研究報(bào)道相比較，表明了本文所提出的預(yù)報(bào)方法的有效性和更高的準(zhǔn)確性。關(guān)鍵詞：聚丙烯熔融指數(shù)預(yù)報(bào)，自適應(yīng)粒子群優(yōu)化算法,徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)， Melt index prediction of propylene polymerization based on adaptive particle swarm optimizationZhao Chengye Liu Xinggao (State Key Laboratory of Industrial Control Technol

4、ogy, Control Department, Zhejiang University, Hangzhou 310027)（E-mail: ）Abstract: A high-precision on-line method of predicting melt index of propylene polymerization based on principal component analysis and adaptive particle swarm optimization is proposed to overcome the high correlation character

5、istics and high nonlinear characteristics in the propylene polymerization process. Adaptive particle swarm optimization (APSO) is employed to get better search efficiency and higher precision than classical particle swarm optimization (PSO), and principal component analysis is used to reduce the com

6、plexity of predicting model. A new method of optimizing both structure and parameters of radial basis function (RBF) network is also proposed. The validity of these methods is demonstrated through practical data in real factory, and research result shows higher precision and shorter computing time t

7、han before.Key words: melts index prediction of propylene polymerization, adaptive particle swarm optimization, radial basis function neural network.1 引言粒子群優(yōu)化算法（PSO）是由Kenedy和Eberhart在1995年提出的群智算法，該算法從鳥類遷徙和魚群巡游過程中提煉出一種簡單的群體運(yùn)動(dòng)機(jī)制，用來引導(dǎo)群體中的粒子找到問題的全局最優(yōu)解。PSO算法易于實(shí)現(xiàn)的優(yōu)點(diǎn)使得它在近年來發(fā)展迅速，并被應(yīng)用到許多實(shí)際領(lǐng)域1。作為一種基于群體的迭代進(jìn)化算法

8、，PSO算法的收斂比較緩慢，而且經(jīng)典的PSO算法在解決復(fù)雜的多維優(yōu)化問題時(shí)容易過早的陷入局部最優(yōu)點(diǎn)。這些弱點(diǎn)限制了PSO算法的應(yīng)用范圍。因此提高PSO算法的收斂速度和防止算法陷入局部最優(yōu)點(diǎn)是PSO算法研究的重要方面，主要措施包括控制算法運(yùn)行時(shí)參數(shù)以及引入附加進(jìn)化操作算子2。RBF神經(jīng)網(wǎng)絡(luò)基于多變量有限點(diǎn)嚴(yán)格插值理論，能夠逼近任意非線性函數(shù)并且具有快速收斂性。文獻(xiàn)3將徑向基函數(shù)網(wǎng)絡(luò)與主元分析分析結(jié)合，建立聚丙烯熔融指數(shù)預(yù)報(bào)模型。文獻(xiàn)4通過構(gòu)造多尺度徑向基函數(shù)(multi-scale RBF)網(wǎng)絡(luò)，提高網(wǎng)絡(luò)的預(yù)報(bào)和泛化能力。文獻(xiàn)5利用主元分析算法和遺傳算法優(yōu)化RBF網(wǎng)絡(luò)。文獻(xiàn)6使用PSO算法優(yōu)化徑

9、向基函數(shù)（RBF）神經(jīng)網(wǎng)絡(luò)，給出了RBF神經(jīng)網(wǎng)絡(luò)編碼的基本原則。文獻(xiàn)7給出了基于PSO等優(yōu)化算法的混合優(yōu)化模型。本文提出以RBF神經(jīng)網(wǎng)絡(luò)為基礎(chǔ)的熔融指數(shù)預(yù)報(bào)建模方法，提出采用主元分析減少數(shù)據(jù)冗余，降低預(yù)報(bào)模型的復(fù)雜度，使用基于自適應(yīng)參數(shù)控制的PSO算法優(yōu)化神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)和參數(shù)，減小人為因素對(duì)建模的影響。通過對(duì)某工廠實(shí)際生產(chǎn)的預(yù)報(bào)以及與國際相關(guān)報(bào)道的比較，表明了該預(yù)報(bào)方法具有較高的精度和較強(qiáng)的泛化能力。2 基于RBF神經(jīng)網(wǎng)絡(luò)的統(tǒng)計(jì)模型基于RBF神經(jīng)網(wǎng)絡(luò)的統(tǒng)計(jì)模型的前端是使用主元分析的數(shù)據(jù)預(yù)處理部分。使用主元分析(PCA)將高維數(shù)據(jù)映射到較低維度，并且保留原數(shù)據(jù)的主要信息。PCA主要包括特征選擇和

10、特征提取兩個(gè)過程。特征選擇的關(guān)鍵是選取特征向量并獲得輸入向量在特征向量上的投影。設(shè)輸入數(shù)據(jù)是均值為零的維隨機(jī)向量，求得其協(xié)方差陣，然后計(jì)算的特征向量矩陣以及相應(yīng)的特征值矩陣。特征提取時(shí)將特征值從大到小排序，保留較大特征值而舍去較小的特征值。最后使用保留的特征值所對(duì)應(yīng)的特征向量重構(gòu)的投影，從而得到維度較低的數(shù)據(jù)。RBF神經(jīng)網(wǎng)絡(luò)具有三層結(jié)構(gòu)，包括輸入層、非線性神經(jīng)元單元層（隱含層）和輸出層。對(duì)于給定輸入向量，RBF網(wǎng)絡(luò)將其映射到輸出向量，的計(jì)算公式為 (1)其中,是徑向基函數(shù)，是隱藏層到輸出層之間的連接權(quán)值，是隱藏層神經(jīng)元數(shù)目，是隱藏層中神經(jīng)元的中心向量。本文使用歐幾里德范數(shù)，并取高斯函數(shù)作為基函

11、數(shù)。RBF神經(jīng)網(wǎng)絡(luò)的主要參數(shù)包括：隱含層節(jié)點(diǎn)個(gè)數(shù)，中心向量，高斯函數(shù)寬度向量，權(quán)值向量，因此網(wǎng)絡(luò)可以由長度為的編碼序列表示。當(dāng)是可變量時(shí)，網(wǎng)絡(luò)的編碼長度將不確定。使用粒子群算法優(yōu)化網(wǎng)絡(luò)時(shí)，可以統(tǒng)一定義編碼長度為，其中是的最大值，并將長度小于該值的編碼將擴(kuò)展到該長度。定義RBF網(wǎng)絡(luò)的適應(yīng)度函數(shù) (2) 其中，是隱含層節(jié)點(diǎn)個(gè)數(shù)，是網(wǎng)絡(luò)的輸出誤差，是可調(diào)參數(shù)。3 經(jīng)典粒子群算法粒子群算法通過模擬粒子群體在解空間上的運(yùn)動(dòng)來尋找優(yōu)化解?？紤]一個(gè)維空間上的優(yōu)化問題，初始化該問題的一個(gè)解集群體，其中個(gè)體是該問題的一個(gè)有效解，是解集的大小。定義個(gè)體的速度為，按照(3)(4)式迭代更新解集和個(gè)體的速度，直到整個(gè)

12、群體的取值滿足給定優(yōu)化條件。 (3) (4)其中，是學(xué)習(xí)速率參數(shù)，和是加速度參數(shù)，和是0與1之間的隨機(jī)數(shù)， 和分別是粒子到達(dá)過的最優(yōu)解和全局最優(yōu)解，是迭代次數(shù)。算法詳細(xì)步驟如下。1) 確定問題的優(yōu)化要求、編碼方式以及算法常量；2) 算法初始化：隨機(jī)初始化個(gè)體的位置、速度和適應(yīng)度；3) 按照(3)式更新個(gè)體的速度，按照(4)式更新個(gè)體的位置；4) 更新每個(gè)個(gè)體到達(dá)過的局部最優(yōu)點(diǎn)；5) 如果需要更新全局最優(yōu)點(diǎn)，則在更新后轉(zhuǎn)到3)，否則轉(zhuǎn)到6)；6) 如果迭代結(jié)束，則輸出全局最優(yōu)點(diǎn)作為問題的最優(yōu)解，否則返回3)；4 自適應(yīng)粒子群算法PSO算法的收斂性能和穩(wěn)定性分析由Clerc和Kennedy給出8。

13、目前更多的研究集中在如何改進(jìn)算法的性能，主要工作包括：算法參數(shù)的選擇，引入附加操作算子，改進(jìn)算法拓?fù)浣Y(jié)構(gòu)。一種改進(jìn)方法是引入慣性學(xué)習(xí)速率，使得學(xué)習(xí)速率的值隨迭代過程線性減小9。 (5)其中是當(dāng)前迭代次數(shù)，是預(yù)定義的最大迭代次數(shù)。最大學(xué)習(xí)速率和最小學(xué)習(xí)速率的取值分別可以是0.9和0.49。加速度參數(shù)和是影響算法的重要參數(shù)，實(shí)現(xiàn)變加速度參數(shù)的PSO算法對(duì)于改進(jìn)算法的性能有重要意義。研究表明在算法初期采用較大的和較小的，并在算法運(yùn)行過程中逐漸減小并增大，能夠提高算法性能10。本文采用一種基于進(jìn)化狀態(tài)估計(jì)（ESE）的自適應(yīng)變參數(shù)PSO算法(APSO）。ESE將PSO算法的進(jìn)化過程分為四個(gè)狀態(tài)：探索期、

14、開拓期、聚合期以及跳出期，并在算法運(yùn)行過程中動(dòng)態(tài)估計(jì)進(jìn)化狀態(tài)，調(diào)整算法參數(shù)。1) 定義粒子與群體中其他粒子的平均距離為 (6)其中，是種群規(guī)模。2) 定義反映群體分布特性的進(jìn)化因子， (7)其中，為當(dāng)前全局最優(yōu)解與其他粒子的平均距離。和分別是中的最小值和最大值。3) 分別定義隸屬于四個(gè)狀態(tài)的隸屬度函數(shù)，如(8)式至(11)式。直觀來說，當(dāng)值取值適中時(shí)算法處于探索期，值取值較小時(shí)算法處于開拓期，值取值接近于0時(shí)算法處于聚合期，值取值接近于1時(shí)算法處于跳出期。 (8)，(9) (10)，(11)在自適應(yīng)PSO算法中，學(xué)習(xí)速率和加速度參數(shù)取決于隨著進(jìn)化因子，而不是算法的迭代次數(shù)。初始化加速度參數(shù)為0

15、.9，在算法運(yùn)行時(shí)使用(12)式計(jì)算。 (12)設(shè)置加速度參數(shù)和的默認(rèn)值為1，并且定義和在算法不同階段的變化規(guī)律，見表1。和的最大取值必須被限定在一定范圍內(nèi)。表1 和的取值在算法不同階段的變化規(guī)律階段的變化的變化探索期增加減小開拓期緩慢增加緩慢減小聚合期緩慢減小緩慢增加跳出期減小增加5 熔融指數(shù)預(yù)報(bào)實(shí)例熔融指數(shù)（MI）是聚丙烯生產(chǎn)控制的主要質(zhì)量指標(biāo)，它決定了產(chǎn)品的牌號(hào)。目前采用離線化驗(yàn)的方法來獲得MI需要2個(gè)小時(shí)甚至更長，無法滿足在線控制的要求。本文采用統(tǒng)計(jì)建模方法對(duì)熔融指數(shù)建立預(yù)報(bào)模型。根據(jù)反應(yīng)機(jī)理和流程工藝，選擇9個(gè)過程變量作為模型的輸入變量，分別是溫度、壓力、液位、氫氣氣相百分?jǐn)?shù)、3股丙

16、烯進(jìn)料流速和2股催化劑進(jìn)料流速?？偣灿?0組輸入輸出數(shù)據(jù)。使用主元分析方法將輸入數(shù)據(jù)維數(shù)減低至6維，經(jīng)過歸一化處理后輸入RBF網(wǎng)絡(luò)。數(shù)據(jù)集被劃分為訓(xùn)練數(shù)據(jù)集和測(cè)試數(shù)據(jù)集，為了驗(yàn)證建模方法的穩(wěn)定性，設(shè)計(jì)了6組不同的訓(xùn)練集和測(cè)試集的組合，分別進(jìn)行建模和預(yù)報(bào)。統(tǒng)計(jì)模型中的結(jié)構(gòu)和參數(shù)完全由優(yōu)化算法計(jì)算得到。本文使用了變參數(shù)的APSO優(yōu)化算法，并且對(duì)比使用標(biāo)準(zhǔn)PSO算法。APSO算法的參數(shù)設(shè)定為：初始化學(xué)習(xí)速度，初始化加速度參數(shù)，迭代次數(shù)為，群體規(guī)模為120。經(jīng)典PSO算法的參數(shù)設(shè)定為：學(xué)習(xí)速度，加速度參數(shù)，迭代次數(shù)為，群體規(guī)模為120。另外，RBF神經(jīng)網(wǎng)絡(luò)統(tǒng)計(jì)模型的適應(yīng)度函數(shù)參數(shù)=10，隱含層中心節(jié)點(diǎn)

17、個(gè)數(shù)在之間取值。實(shí)驗(yàn)結(jié)果表明：第一，APSO算法相比經(jīng)典PSO算法的收斂速度更快，迭代計(jì)算次數(shù)顯著減少，見圖2；第二，APSO算法的多項(xiàng)參數(shù)能夠隨著具體優(yōu)化問題的不同自動(dòng)調(diào)節(jié)，減少了人為因素對(duì)建模的影響；第三，使用APSO算法建立的熔融指數(shù)預(yù)報(bào)模型具有相對(duì)較高的精度。表2給出了本文的研究結(jié)果與國際上相關(guān)研究報(bào)道結(jié)果的比較，數(shù)據(jù)表明本文提出的預(yù)報(bào)模型比國際上目前最好的報(bào)道結(jié)果5在平均相對(duì)誤差和均方根誤差上都有了顯著改進(jìn)。表2 本文建模方法與參考文獻(xiàn)中的結(jié)果比較文獻(xiàn)建模方法 MRERMSE本文APSO-RBF0.66%0.0237Han, 2004 3Weighted LS-SVM0.15Cao,

18、 1999 11魯棒自適應(yīng)RBF0.62Shi, 2006 12ICA-MS-RBF2.98%0.0794Lou, 2007 5PCA - GA - RBF0.84%圖1給出了本文提出的模型在不同批次數(shù)據(jù)上的預(yù)報(bào)結(jié)果，其中實(shí)線為熔融指數(shù)的實(shí)測(cè)值，虛線為統(tǒng)計(jì)模型的預(yù)報(bào)值。圖中顯示了預(yù)報(bào)模型能夠相當(dāng)精確地逼近測(cè)試數(shù)據(jù)，表明了本文建立的預(yù)報(bào)模型具有較高的預(yù)報(bào)精度和較強(qiáng)的泛化能力。圖1(a) APSO算法得到的統(tǒng)計(jì)模型在訓(xùn)練數(shù)據(jù)上的預(yù)報(bào)結(jié)果圖1(b) APSO算法得到的統(tǒng)計(jì)模型在測(cè)試數(shù)據(jù)上的預(yù)報(bào)結(jié)果圖2 對(duì)比了APSO算法和PSO算法在相同訓(xùn)練數(shù)據(jù)上的收斂速度，個(gè)體適應(yīng)度的計(jì)算公式參見(2)式。APS

19、O算法的收斂速度明顯快于PSO算法。圖3顯示了APSO算法運(yùn)行過程中學(xué)習(xí)速度值的變化曲線，可以看到學(xué)習(xí)速率的變化取決于進(jìn)化狀態(tài)而不是迭代次數(shù)。圖2 APSO算法和PSO算法收斂速度比較圖3 APSO算法中學(xué)習(xí)速率的變化曲線6 結(jié)論本文提出一種基于自適應(yīng)粒子群優(yōu)化算法和徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)的熔融指數(shù)預(yù)報(bào)新方法。該方法融合自適應(yīng)PSO算法的高效全局尋優(yōu)能力和徑向基函數(shù)網(wǎng)絡(luò)的精確擬合能力，相比以往的預(yù)報(bào)方法，它簡化了建模的復(fù)雜度，又提高了模型的預(yù)報(bào)精度和泛化能力；相比以往使用經(jīng)典PSO算法，它提高了優(yōu)化過程的效率和精度。通過對(duì)實(shí)際丙烯聚合過程數(shù)據(jù)進(jìn)行建模和預(yù)報(bào)，證明了使用該方法建立的聚丙烯熔融指數(shù)預(yù)報(bào)

20、模型具有準(zhǔn)確性和可靠性，對(duì)于改進(jìn)聚丙烯生產(chǎn)的過程控制具有重要意義。參考文獻(xiàn)1 J. Kennedy, R. C. Eberhart, and Y. H. Shi, Swarm Intelligence. San Mateo, CA: Morgan Kaufmann, 2001.2 R. C. Eberhart and Y. Shi, “Guest editorial,” IEEE Trans. Evol. Comput.Special Issue Particle Swarm Optimization, vol. 8, no. 3, pp. 201203, Jun. 2004.3 In-Su

21、Han, Chonghun Han, Chang-Bock Chung. Melt Index Modeling with Support Vector Machines, Partial Least Squares, and Artificial Neural Networks J. Journal of Applied Polymer, 2005, Vol. 95: 967-974.4 Stephen A. Billings, Hua-Liang Wei, Michael A. Balikhin. Generalized multiscale radial basis function n

22、etworks J. Neural Networks, 2007, 20: 10811094.5 樓巍, 劉興高. 基于PCA - GA - RBF網(wǎng)絡(luò)的聚丙烯熔融指數(shù)預(yù)報(bào)模型 J. 石油化工高等學(xué)校學(xué)報(bào), 2007, 20(3):82-85.6 Hsuan-Ming Feng. Self-generation RBFNs using evolutional PSO learning J. Neurocomputing, 2006, 70: 2412517 Cheng-Jian Lin, Shang-Jin Hong. The design of neuro-fuzzy networks using particle swarm optimization and recursive singular value decomposition J. Neurocomputing, 2007