含菁英政策之基因演算法於模式最佳化近似的應用

(1)

工程科技與教育學刊第七卷第三期民國九十九年六月第 437～444 頁

含菁英政策之基因演算法於模式最佳化近似的應用

林家富、劉宏毅、蔡金橤、郭東義 國立高雄應用科技大學化學工程與材料工程系 E-Mail：[email protected]

摘要

本文是應用含菁英政策之基因演算法於複雜系統之模式最佳化近似。此最佳化近似模式以含有時間延遲之有理轉移函數描述，並藉由最小化頻率響應L2_{誤差準則來搜尋最佳化近似模式的參數。以頻率響應做為誤差準} 則之優點是此方法可有效的應用於以有理及無理函數描述的系統。模擬結果顯示，相較於傳統的基因演算法，本文提出的方法可降低基因演算過早收歛入局部最佳解，最佳近似模式有較好的適應函數值。關鍵詞：基因演算法，菁英政策，頻率響應

Abstract

In this paper, the problem of optimally approximating complex systems is solved by a genetic algorithm with elitism. The optimal approximate model described by a rational transfer function with a time delay is sought such that a frequency-domain L2_{-error criterion is minimized. The chosen frequency-domain error criterion makes the proposed}

approach quite efficacious for optimally approximating systems having rational and irrational functions. The simulated results show that the proposed approach is able to avoid falling into the local maximum due to the parameter convergence and therefore the proposed optimal approximants have better fitness than that obtained by a traditional genetic algorithm.

Keyword: genetic algorithm, elitism, frequency response.

1. 前言

隨著控制技術與控制要求的提高，系統的階次也越來越高，所以在模擬與控制器中設計一個動態線性系統，對系統的模式做近似是一件非常重要的工作。舉例來說，一個工廠的回饋系統中，其轉移函數之分母若包含著時延( dead time)[1]項次，這使得在時域分析或是控制器的設計時，會變得更加困難。目前對於解決模式最佳化近似的問題，已經提出許多方法[2]，其可分為兩大類：一個是非性能指標導向的方法，主要是藉由將高階系統截去部分資訊轉變成低階系統，此方式是相當的容易。而將高階系統轉成低階系統的動作，稱之為模式簡化 (model reduction)，其方法有傳統的 Padé 模式簡化法(model reduction method)[3,4]、Routh 近似法[4,5,6]、平衡截斷法(balanced truncation method)[7,8]以及最佳 Hankel 近似法[9,10]等；另一個則是性能指標導向的方法，它藉由最小化某一誤差準則來獲得近似模式，其方法則有H₂-norm[11]、 -norm[12,13]、L₂ H_∞-norm[14]以及 -norm[15] 等，這一類的方法，它們必須藉助一套有效率演算法來計算性能指標。以 L_∞ 2 H -norm、 -norm 為例，常用的方法是以梯度法為基礎(gradient-based)的搜尋法則[16]，但不幸地，梯度法通常都會陷入局部的最佳解。同樣的， 2 L H∞-norm 及 -norm 的困難點在於它們被選擇為性能指標時，就必須求其導數，但這兩者沒有任何一點可微分的性質。所以，仍然缺少計算最小化以 L∞ H_∞-norm 或L_∞-norm 為性能指標之有效率的計算方法。

基因演算法(genetic algorithm, GA)起源於 60 年代對自然和人工自適應系統的研究，最早由 Holland 教授所

(2)

林家富、劉宏毅、蔡金橤、郭東義 438 提出[17]，主要是依據生物學家查理．達爾文適者生存理論所發展出來的一套演算法則，其基本精神在於模仿生物界物競天擇、適者生存的自然演化過程，所謂『適者生存，不適者淘汰』即是自然界為繁衍出更優良的下一代，會隨著環境的變化而自我調適，使子代(offspring)的基因排列優於母代(parents)的基因排列，而產生比母代更適應環境的子代。因此，基因演算法可以說是一種仿傚生物演化來求解最佳化問題的工具。雖然基因演算的搜尋技術是以隨機搜尋為架構，但其採用多點搜尋的方式在參數空間中同時尋求問題最佳解且運用演化過程所提供的資訊，因此是一個強有力的數值最佳化之演算法[18]。本文是應用基因演算法於複雜系統之最佳化近似，我們於基因演算法加入菁英政策(elitist strategy)，使其在選擇複製、交配、突變等運算子中，都能一直保有原本最佳的基因，讓其參數搜尋時，可加速基因演算收斂速度和降低基因演算過早收歛落入局部最佳解，找到近似於原系統最佳化的參數解。

2. 模式最佳化近似之問題描述

對於控制系統的族群中，最能充分表示一個複雜動態系統( complex dynamic system )的就是包含有時間延遲 ( time delay )的低階轉移函數。考慮一個單輸入單輸出的高階系統或非有理函數 G(s)，主要希望能將其表示為 如下m 階具時延的近似模式H s_m( )： 1 0 1 1 1 0 1 1 ( ) d m s m m m m m b b s b s H s e a a s a s s τ − − − − − + + = ⋅ + + + L L (1)

而H s_m( )包含了G(s)的特徵。在本文中，將依據如下的L₂誤差性能指標(L₂-error performance index ):

2 0 | ( ) ( ) | N m i m i J G jω H jω = =

∑

− _i (2) 其中ω 為頻率點，i=1,2,… ,N，搜尋 G(s)的最佳化近似式_i H s_m( )。亦即，本文目標希望在ω 點時，G(s)_i 與H sm( )的L2誤差性能指標越小越好。倘若G(s)為漸近穩定(asymptotically stable )，則原系統 G(s)和近式模式 ( ) m H s 滿足 ) 0 ( ) 0 ( G Hm = (3) 藉由此限制條件，主要是可確保輸入訊號為單位步階響應(unit-step response)時，G(s)與H s_m( )具有相同的穩態響應值。而參數的最佳化選擇，則是以最小化式(2)之性能指標為準則並藉由傳統的基因演算法與具菁英政策之基因演算法去尋找出最佳之參數，a_i、、b_i τ ，i=0,1,… ,m-1。 _d

(3)

含菁英政策之基因演算法於模式最佳化近似的應用 443 比較圖9 和上述的誤差積分平方值，本文提出的方法與原系統時域響應之近似情形明顯優於傳統基因演算法。 * () ()_m Gj H j ωω − Frequecy( )ω 10-2 10-1 100 101 102 103 10-7 10-6 10-5 10-4 10-3 10-2 二階近似式(GA) 二階近似式(GAE) 三階近似式(GA) 三階近似式(GAE) 圖8 原系統與最佳近似模式的頻域誤差大小響應圖 Time (second) * () () m yt y t − 0 20 40 60 80 -2 -1 0 1 2 3 4 5 6x 10 -3 二階近似式(GA) 二階近似式(GAE) 三階近似式(GA) 三階近似式(GAE) 圖9 原系統式與最佳近似模的時域單位步階響應誤差圖

5. 結論

本文所提出的含菁英政策之基因演算法尋找最佳近似模式，於性能指標、頻率響應和時域響應，均較傳統的基因演算法為佳。但是，和傳統的基因演算法相同，在不同參數搜尋範圍就會得到不同之近似效果，當參數搜尋範圍設定的太小或太大，就都會造成基因演算法在運算時，而引起了參數的局部最佳化，導致運算提早收歛，而無法得到最佳的近似解。整體而言，基因演算法的交配率及突變率、參數搜尋範圍和菁英政策保留最佳染色體個數等參數之設定，均會影響到基因演算法應用於模式是最佳近似之結果。

參考文獻

[1] Hugo A. J., Taylor P. A., and Wright J. D., “Approximate dynamic models for recycle systems,” Ind. End. Chem. Res., vol. 35, pp. 485-487, 1996.

[2] Fortuna L., Nunnari G., and Gallo A., Model Order Reduction Techniques With Applications in Electrical Engineering. London, U.K.: Springer-Verlag, 1992.

[3] Bultheel A. and Van Barel M., “Padé techniques for model reduction in linear system theory: A survey,” J. Comput. Appl. Math., vol. 14, pp. 401-438, 1986.

[4] Mittal S. K., Chandra D., and Dwivedi B., “Improved Routh-Padé approximants Using Vector genetic algorithm to unstable systems,” Int. J. of Engr. and Appl. Sci., vol. 1, issue 2, pp. 1-14, 2009.

[5] Hutton H. F. and Friedland B. F., “Routh approximations for reducing order of linear time-invariant systems,” IEEE Trans. Auto. Control, vol. 20, pp. 329-337, 1975.

[6] Panda S., Tomar S. K., Prasad R., and Ardil C., “Reduction of Linear Time-Invariant Systems Using Routh-Approximation and PSO,” Int. J. of Appl. Math. and Comp. Sci., vol. 5, pp. 82-89, 2009.

[7] Moore B. C., “Principal component analysis in linear systems: Controllability, observability and model reduction,” IEEE Trans. Automat. Contr., vol. 26, pp. 17–31, 1981.

[8] Zilouchian A., “Balanced structures and model reduction of unstable systems,” in Proc. IEEE Southeastcon, vol. 2, pp. 1198-1201, 1991.

[9] Kung S. Y. and Lin D. W., “Optimal Hankel-norm model reductions: Multivariable systems,” IEEE Trans. Automat. Contr., vol. AC–26, pp. 832–852, 1981.

(4)

林家富、劉宏毅、蔡金橤、郭東義 444

[10] Chui C. K., Li X., and Ward J. D., “System reduction via truncated Hankel matrices,” Math. Contr., Signals, Syst., vol. 4, no. 2, pp. 161-175, 1991.

[11] Collins E. G., Ying S. S., Haddad W. M., and Richter S., “An efficient, numerically robust homotopy algorithm for H2 model

reduction using the optimal projection equations,” Math. Modeling Syst.: Meth., Tools Applicat. Eng. Rel. Sci., vol. 2, pp. 101-133, 1996.

[12] Spanos J. T., Milman M. H., and Mingori D. L., “A new algorithm for L2_{optimal model reduction,” Automatica, vol. 28, pp.}

897-909, 1992.

[13] Krajewski W., Lepschy A., Mian G. A., and Viaro U., “Optimality conditions in multivariable L2_{model reduction,” J. Franklin}

Inst., vol. 330, pp. 431-439, 1993.

[14] Grigoriadis K. M., “Optimal H∞_{model reduction via linear matrix in equalities: Continuous- and discrete-time cases,” Syst.} Contr. Lett., vol. 26, pp. 321-333, 1995.

[15] Kavranoglu D. and Bettayeb M., “Characterization and computation of the solution to the optimal L∞_{approximation problem,”} IEEE Trans. Automat. Contr., vol. 39, pp. 1899-1904, Sept. 1994.

[16] Aplevich J. D., “Gradient methods for optimal linear system reduction,” Int. J. Contr., vol. 18, pp. 767-772, 1973. [17] Holland J. H., Adaptation in Nature and Artificial Systems, University of Michigan Press, 1975.

[18] Goldberg D. E., Genetic Algorithms in Search, Optimization and Machine Learning, Addison Wesley, New York, 1989. [19] Cheng S. L. and Hwang C., “Optimal approximation of linear systems by differential evolution algorithm,” IEEE Trans. on

Systems, Man, and Cybernetics-Part A: Systems and Humans, vol. 31, pp. 698-707, 2001.

[20] Hwang C. and Chuang Y. H., “Computation of optimal reduced-order models with time delay,” Chem. Eng. Sci., vol. 49, pp. 3219-3296, 1994.

含菁英政策之基因演算法於模式最佳化近似的應用