【正文】
fun39。 x0=[1,1]。b=2。,x,options,xl,xu)。 [x,options]=constr(39。inf。2]。 options(13)=0 %等式約束的個(gè)數(shù)為 0 xl=[0。10。 G(1)=(x(1)+x(2))*x(3)120。,x0,…) 具體含義請(qǐng)參見(jiàn) tbp154 求解非線性規(guī)劃問(wèn)題,需要首先編制一個(gè) m文件,描述需要求解的問(wèn)題,也即求解需以“ 被調(diào) ”的形式進(jìn)行 數(shù)學(xué)實(shí)驗(yàn) 非線性規(guī)劃問(wèn)題的范例 圈地 1 某旅游發(fā)展有限公司計(jì)劃開(kāi)發(fā)度假村,公司要求先用一批舊磚建一圈矩形圍墻,以便存放建筑材料.舊磚的數(shù)量是固定的,圍墻的高度不能低于兩米,圍墻圍住的面積越大越好.要你來(lái)設(shè)計(jì)。,x0,OPTIONS,VLB,VUB) 具體含義請(qǐng)參見(jiàn)聯(lián)機(jī) help 數(shù)學(xué)實(shí)驗(yàn) MATLAB非線性規(guī)劃函數(shù) 無(wú)約束優(yōu)化函數(shù)介紹: fminu (fminunc) 第10.6節(jié) 非線性規(guī)劃函數(shù) 08 m in ( )fx調(diào)用語(yǔ)法: fminu [X,OPTIONS] = fminu(39。 相對(duì)而言無(wú)約束優(yōu)化更容易求解一些,一個(gè)基本的思路就是, 化約束優(yōu)化為一系列無(wú)約束 優(yōu)化問(wèn)題,例如:序列無(wú)約束極小化技術(shù)、可變?nèi)莶罘ǎ鹊?~~ 對(duì)于具體算法細(xì)節(jié),我們雖然不會(huì)本課程中深入探究,但能夠從 算法本身的著手改進(jìn) 將會(huì)是很有價(jià)值的工作,比如本章作者開(kāi)發(fā)的逼近精確罰函數(shù)法,感興趣可以仔細(xì)研讀 第9.34節(jié) 非線性規(guī)劃方法概要 06 數(shù)學(xué)實(shí)驗(yàn) MATLAB非線性規(guī)劃函數(shù) 約束優(yōu)化函數(shù)介紹: constr 第10.6節(jié) 非線性規(guī)劃函數(shù) 07 m i n ( ). . 1 ( ) 02 ( ) 0fxs t g xgxvl b x vu b????調(diào)用語(yǔ)法: constr (fmincon) [X,OPTIONS] = constr(39。 MATLAB優(yōu)化工具箱 (Optim)針對(duì)約束和非約束分別由函數(shù) constr和 fminu進(jìn)行求解。 非線性規(guī)劃問(wèn)題同樣是企業(yè)運(yùn)作、科技研發(fā)和工程設(shè)計(jì)的常見(jiàn)問(wèn)題,甚至在某種意義上應(yīng)用面比線性規(guī)劃更廣。x0=zeros(4,1)。 ub=inf*ones(4,1)。beq=[0]。24]。 b=[0。2,3,0,0。200]。1000。 [x,fval,exitflag,output,lambda]=linprog(f,A,b,[],[],lb)。30]。 b=[20。3,2,4。6]。*x 數(shù)學(xué)實(shí)驗(yàn) 線性規(guī)劃的 LP解法 相關(guān)函數(shù)介紹: linprog 第9.2節(jié) 線性規(guī)劃LP解法 10 m i n..(, , , , ,Txfxs t A x bAeq x b e ql b u u bf x b b e q l b u b A Aeq??????沒(méi)有賦空[ ] )其中 和 為向量 和 為矩陣x=linprog(f