【正文】 , 則計(jì)算結(jié)果為 : x=[3． 0000 5． 0000 4． 0000 7． 0000 1． 0000 0 0 0 0 0 5． 0000 11． 0000 5． 6959 4． 9285 7． 2500 7． 7500]’ fval =89． 8835 exitflag = 1 總的噸千米數(shù) 89． 8835比上面結(jié)果更好． 通過(guò)此例可看出 fmincon函數(shù)在選取初值上的重要性． MATLAB（ gying2） 返回 。 料場(chǎng)位置為(xj， yj)，日儲(chǔ)量為 ej， j=1,2；料場(chǎng) j向工地 i的運(yùn)送量為 Xij． 目標(biāo)函數(shù)為：? ?? ?????216122 )()(mi nj iijijij byaxXf 約束條件為： 2,1 ,6,2,1 ,6121????????jeXidXjiijijij? 當(dāng)用臨時(shí)料場(chǎng)時(shí)決策變量為： Xij， 當(dāng)不用臨時(shí)料場(chǎng)時(shí)決策變量為： Xij， xj， yj． （二）使用臨時(shí)料場(chǎng)的情形 使用兩個(gè)臨時(shí)料場(chǎng) A(5,1)， B(2,7)．求從料場(chǎng) j向工地 i的運(yùn)送量 Xij . 在各工地用量必須滿足和各料場(chǎng)運(yùn)送量不超過(guò)日儲(chǔ)量的條件下，使總的噸千米數(shù)最小，這是線性規(guī)劃問(wèn)題． 線性規(guī)劃模型為： ? ?? ??2161),(m i nj iijXjiaaf2,1 , 6,2,1 , s . t .6121????????jeXidXjiijijij?其中 22 )()(),( ijij byaxjiaa ???? ， i = 1, 2, … , 6, j = 1, 2, 為常數(shù) . 設(shè) X11=X1, X21= X 2, X31= X 3, X41= X 4, X51= X 5, X61= X 6 X12= X 7, X22= X 8, X32= X 9, X42= X 10, X52= X 11, X62= X 12 編寫程序 gying1． m MATLAB（ gying1） 計(jì)算結(jié)果為： x =[ 3． 0000 5． 0000 0． 0000 7． 0000 0． 0000 1． 0000 0． 0000 0． 0000 4． 0000 0． 0000 6． 0000 10． 0000]’ fval = 136． 2275 即由料場(chǎng) A 、 B 向 6 個(gè)工地運(yùn)料方案為： 1 2 3 4 5 6 料場(chǎng) A 3 5 0 7 0 1 料場(chǎng) B 0 0 4 0 6 10 總的噸千米數(shù)為 13 275 . （三）改建兩個(gè)新料場(chǎng)的情形 改建兩個(gè)新料場(chǎng)，要同時(shí)確定料場(chǎng)的位置 (xj,yj)和運(yùn)送量Xij，在同樣條件下使總噸千米數(shù)最?。@是非線性規(guī)劃問(wèn)題．非線性規(guī)劃模型為： ? ?? ?????216122 )()(mi nj iijijij byaxXf 2161. , 1 , 2 , , 6 , 1 , 2ij ijij jiX d iX e j???????? 設(shè) X11=X1, X21= X 2, X31= X 3, X41= X 4, X51= X 5, X61= X 6 X12= X 7, X22= X 8, X32= X 9, X42= X 10, X52= X 11, X62= X 12 x1=X13, y1=X14, x2=X15, y2=X16 （ 1）先編寫 M文件 liaoch． m定義目標(biāo)函數(shù)． MATLAB（ liaoch） (2) 取初值為線性規(guī)劃的計(jì)算結(jié)果及臨時(shí)料場(chǎng)的坐標(biāo) : x0=[3 5 0 7 0 1 0 0 4 0 6 10 5 1 2 7]39。mycon239。fun39。VUB=[5 10]。2． 5]。x(1)^2x(2)^27]。 f=2*x(1)x(2)。mycon39。fun439。vub=[]。beq=[0]。b=[]。1]。 x(1)*x(2)10]。 1 221 2 1 2 2( ) e ( 4 2 4 2 1 )xf x x x x x x? ? ? ? ? x1+x2=0 . +x1x2 x1 x2 0 x1x2 –10 0 ??例 3 2．再建立 M文件 mycon． m定義非線性約束： function [g,ceq]=mycon(x) g=[x(1)+x(2)。,x0,A,b,Aeq,beq,VLB,VUB) 4． 運(yùn)算結(jié)果為： x = 0． 7647 1． 0588 fval = 2． 0294 1． 先建立 M文件 fun4． m定義目標(biāo)函數(shù) : function f=fun4(x)。 [x,fval]=fmincon(39。0]。beq=[]。5]。1 4]。1]。 2．一般非線性規(guī)劃 標(biāo)準(zhǔn)型為： m i n F ( X ) s .t ． AX ? b be qXA e q ?? G ( X ) 0? Ceq ( X )= 0 VL B ? X ? V U B 其中 X為 n維變?cè)蛄浚?G(X)與 Ceq(X)均為非線性函數(shù)組成的向量，其他變量的含義與線性規(guī)劃、二次規(guī)劃中相同．用MATLAB求解上述問(wèn)題，基本步驟分三步： 2. 若約束條件中有非線性約束 : G (X) 0? 或 Ceq(X)= 0 ,則建立 M 文件 nonlcon .m 定義函數(shù) G(X) 與 Ceq (X) : function [G,Ceq]=nonlcon(X) G= … Ceq= … 3． 建立主程序 .求解非線性規(guī)劃的函數(shù)是 fmincon,命令的基本格式如下： (1) x=fmincon(‘fun’,X0,A,b) (2) x=fmincon(‘fun’,X0,A,b,Aeq,beq) (3) x=fmincon(‘fun’,X0,A,b, Aeq,beq,VLB,VUB) (4) x=fmincon(‘fun’,X0,A,b,Aeq,beq,VLB,VUB,’nonlcon’) (5)x=fmincon(‘fun’,X0,A,b,Aeq,beq,VLB,VUB,’nonlcon’,options) (6) [x,fval]= fmincon(… ) (7) [x,fval,exitflag]= fmincon(… ) (8)[x,fval,exitflag,output]= fmincon(… ) 輸出極值點(diǎn) M文件 迭代的初值 參數(shù)說(shuō)明 變量上下限 注意： [1] fmincon函數(shù)提供了大型優(yōu)化算法和中型優(yōu)化算法．默認(rèn)時(shí) : 若在 fun函數(shù)中提供了梯度（ options參數(shù)的 GradObj設(shè)置為’ on’），并且只有上下界存在或只有等式約束， fmincon函數(shù)將選擇大型算法．當(dāng)既有等式約束又有梯度約束時(shí)，使用中型算法． [2] fmincon函數(shù)的中型算法使用的是序列二次規(guī)劃法．在每一步迭代中求解二次規(guī)劃子問(wèn)題，并用 BFGS法更新拉格朗日Hesse矩陣． [3] fmincon函數(shù)可能會(huì)給出局部最優(yōu)解，這與初值 X0的選取有關(guān)． 1． 寫成標(biāo)準(zhǔn)形式 ： . ?????????????????00546322121xxxx?????????????2100xx222121 21212m in xxxxf ?????222121 21212m i n xxxxf ????? 2x1+3x2 6 . x1+4x2 5 x1,x2 0 ??例 2 ?2． 先建立 M文件 fun3． m: function f=fun3(x)。 [x,z]=quadprog(H,c,A,b,Aeq,beq,VLB,VUB) 3． 運(yùn)算結(jié)果 為： x =3 2 z = T111222 2 2 21m in ( , )2 4 62xxz x x ?? ? ? ?? ? ? ???? ? ? ?? ? ? ???? ? ? ?? ? ? ?1212 1 1 21 2 200xxxx??? ? ? ???? ? ? ???? ? ? ?? ? ? ????????????????? ?????. 1． 首先建立 M文件 ,用來(lái)定義目標(biāo)函數(shù) F（ X） : function f=fun(X)。0]。beq=[]。2]。 2 4]。6]。 1 2]。 8． [x,fval,exitflag,output]=quaprog(…)。 6． [x,fval]=quaprog(…)。 4． x=quadprog(H,C,A,b, Aeq,beq ,VLB,VUB,X0)。 2． x=quadprog(H,C,A,b,Aeq,beq)。 y=fun1(x) : x= y = 1221 2 1 2 2( ) ( 4 2 4 2 1 ) e xf x x x x x x? ? ? ? ? ?*非線性規(guī)劃的基本解法 非線性規(guī)劃的基本概念 非線性規(guī)劃 返回 定義 如果目標(biāo)函數(shù)或約束條件中至少有一個(gè)是非線性函數(shù) ,則最優(yōu)化問(wèn)題就叫做 非線性規(guī)劃問(wèn)題 ． 非現(xiàn)性規(guī)劃的基本概念 一般形式 : （ 1） 其中 ， 是定義在 Rn 上的實(shí)值函數(shù)，簡(jiǎn)記 : ? ?Xfm in ji hgf , 其它情況 : 求目標(biāo)函數(shù)的最大值，或約束條件小于等于零兩種情況，都可通過(guò)取其相反數(shù)化為上述一般形式． 1 n j 1 n i 1 n R : h ,R : g ,R : R R R f ? ? ? ? ? n T n R x x x X ? ? , , , 2 1 ? ? ? ? ? ? ? ? ? ? ? ? ? ? . ,..., 2 , 1 0 m。fun139。 M文件 : x0 = [1, 1]。steepdesc39。dfp39。bfgs39。 (默認(rèn)值 ),使用大型算法 LargeScale=‘off39。cubicpoly39。,0,)。 主程序?yàn)?: [x,fval]=fminbnd(39。2*exp(x).*sin (x)39。 fplot(f,[0,8])。,value2,...) 創(chuàng)建名稱為 oldops的參數(shù)的拷貝 ,用指定的參數(shù)值修改 oldops中相應(yīng)的參數(shù) . 返回 用 MATLAB解無(wú)約束優(yōu)化問(wèn)題 1. 一