張書花，李艷龍，李 磊,景孟旗

（1,4、河海大學(xué)水利水電學(xué)院，南京，210098；2，3、河海大學(xué)力學(xué)與材料學(xué)院，南京，210098）

0 引言

一般的數(shù)值方法在求解大規(guī)模線性約束優(yōu)化問題時，達(dá)不到計算時間和內(nèi)存的要求。因此本文對大規(guī)模優(yōu)化問題，首先運(yùn)用信賴域方法將原優(yōu)化問題的設(shè)計變量轉(zhuǎn)換為搜索方向，再結(jié)合零空間方法和線搜索技術(shù)對移動漸近函數(shù)進(jìn)行構(gòu)造，運(yùn)用算例驗證了該算法在求解大規(guī)模線性等式約束優(yōu)化問題方面的有效性。通過求解問題得到下降搜索方向,然后計算步長，得到下一迭代點

Svanberg在1987年對優(yōu)化問題首次提出MMA法，用移動漸近線函數(shù)逼近優(yōu)化問題中的目標(biāo)和約束函數(shù)，從而產(chǎn)生子問題。通過求解子問題，最終求得約束優(yōu)化問題的最優(yōu)解。本文在文獻(xiàn)[6,7]的基礎(chǔ)上，對移動漸近線函數(shù)進(jìn)行了深入的研究，構(gòu)造了含參數(shù)的新的移動漸進(jìn)線函數(shù)，并形成新的問題。子問題是嚴(yán)格凸且可分的，不需更新信賴域半徑。在此基礎(chǔ)上，獲得了一個解線性等式約束優(yōu)化問題的新MMA算法。討論了參數(shù)的選取準(zhǔn)則，證明了新算法的全局收斂性，并通過試驗結(jié)果得出算法的有效性。

1 新MA子問題

文獻(xiàn)[6]對MA法提出，并用于線性等式約束優(yōu)化問題(1)。在當(dāng)前迭代點處構(gòu)造的子問題：

2 新MMA算法步驟及收斂性

類似于文獻(xiàn)[6]中定理的證明，得到算法1具有全局收斂性。

3 數(shù)值驗算

在本節(jié)中，進(jìn)行數(shù)值試驗：將算法1與文獻(xiàn)[3]中的MMA算法(以下簡稱為OMMA算法)、投影梯度算法(以下簡稱為PGM算法)進(jìn)行比較。

表1 ： 算法1與PGM算法、OMMA算法的運(yùn)行結(jié)果Table 1： The numerical results of algorithm 1、PGM and OMMA

由計算結(jié)果表1可以看出，對于優(yōu)化問題，算法1的效果優(yōu)于OMMA算法和PGM算法，算法1的運(yùn)行時間也少于OMMA算法和PGM算法，因此，算法1在解大規(guī)模線性約束優(yōu)化問題方面非?？扇?。

4 結(jié)論

為了有效地求解線性等式約束優(yōu)化問題，本文構(gòu)造了一個新的移動漸近線函數(shù)，進(jìn)而建立了新的移動漸近線子問題，并討論了子問題中參數(shù)的選取策略，求出子問題的最優(yōu)解，然后將此最優(yōu)解作為原問題的一個下降方向進(jìn)行線搜索。我們對該算法進(jìn)行數(shù)值試驗，得出算法1是一種有效算法，適合求解大規(guī)模線性等式約束優(yōu)化問題，有可能求得全局最優(yōu)解，新的MMA算法可應(yīng)用于大型工程計算。

[1] Bulteau J P,Vial J P.A restricted trust-region algorithm for unconstrained optimization, J Optim.Theory Appl., 1985, 47：413-435

[2] Shi Z J, Guo J H. A new trust region method for unconstrained optimization. Journal of Computational and Applied Mathematics, 2008, 213：509-520

[3] Fletcher R. Practical methods of optimization.Chichester： John Wiley and Sons,1981.

[4] 袁亞湘,孫文瑜. 最優(yōu)化理論與方法[M]. 北京：科學(xué)出版社，1997

[5] Svanberg K.The method of moving asymptotes—a new method for structural optimization[J].International Journal for Numerical Methods in Engineering, 1987, 24(2)： 359-373

[6] 王海軍. 解非線性最優(yōu)化問題的移動漸近線法及應(yīng)用；[D]南京，南京航空航天大學(xué)；2010

[7] 胡平,賈朝輝,倪勤.一種解無約束優(yōu)化問題的新移動漸近線算法；工程數(shù)學(xué)學(xué)報；2012