1、用MATLAB優(yōu)化工具箱解線性規(guī)劃min z=cX 1、模型：命令：x=linprog（c，A，b）2、模型：命令：x=linprog（c，A，b，Aeq,beq）注意：若沒有不等式：存在，則令A(yù)= ，b= .若沒有等式約束, 則令A(yù)eq= , beq= .3、模型：命令：1 x=linprog（c，A，b，Aeq,beq, VLB，VUB） 2 x=linprog（c，A，b，Aeq,beq, VLB，VUB, X0）注意：1 若沒有等式約束, 則令A(yù)eq= , beq= . 2其中X0表示初始點4、命令：x,fval=linprog()返回最優(yōu)解及處的目標(biāo)函數(shù)值fval.例1 max 解編

2、寫M文件小xxgh1.m如下：c=-0.4 -0.28 -0.32 -0.72 -0.64 -0.6; A=0.01 0.01 0.01 0.03 0.03 0.03;0.02 0 0 0.05 0 0;0 0.02 0 0 0.05 0;0 0 0.03 0 0 0.08; b=850;700;100;900; Aeq=; beq=; vlb=0;0;0;0;0;0; vub=;x,fval=linprog(c,A,b,Aeq,beq,vlb,vub)例2 解: 編寫M文件xxgh2.m如下： c=6 3 4; A=0 1 0; b=50; Aeq=1 1 1; beq=120; vlb=3

3、0,0,20; vub=; x,fval=linprog(c,A,b,Aeq,beq,vlb,vub例3 （任務(wù)分配問題）某車間有甲、乙兩臺機床，可用于加工三種工件。假定這兩臺車床的可用臺時數(shù)分別為800和900，三種工件的數(shù)量分別為400、600和500，且已知用三種不同車床加工單位數(shù)量不同工件所需的臺時數(shù)和加工費用如下表。問怎樣分配車床的加工任務(wù)，才能既滿足加工工件的要求，又使加工費用最低？解設(shè)在甲車床上加工工件1、2、3的數(shù)量分別為x1、x2、x3，在乙車床上加工工件1、2、3的數(shù)量分別為x4、x5、x6。可建立以下線性規(guī)劃模型：編寫M文件xxgh3.m如下:f = 13 9 10 11

4、 12 8;A = 0.4 1.1 1 0 0 0 0 0 0 0.5 1.2 1.3;b = 800; 900;Aeq=1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1;beq=400 600 500;vlb = zeros(6,1);vub=;x,fval = linprog(f,A,b,Aeq,beq,vlb,vub)例4某廠每日8小時的產(chǎn)量不低于1800件。為了進行質(zhì)量控制，計劃聘請兩種不同水平的檢驗員。一級檢驗員的標(biāo)準(zhǔn)為：速度25件/小時，正確率98%，計時工資4元/小時；二級檢驗員的標(biāo)準(zhǔn)為：速度15小時/件，正確率95%，計時工資3元/小時。檢驗員每錯檢一次，

5、工廠要損失2元。為使總檢驗費用最省，該工廠應(yīng)聘一級、二級檢驗員各幾名？解設(shè)需要一級和二級檢驗員的人數(shù)分別為x1、x2人,則應(yīng)付檢驗員的工資為：因檢驗員錯檢而造成的損失為：故目標(biāo)函數(shù)為：約束條件為：線性規(guī)劃模型：編寫M文件xxgh4.m如下：c = 40;36;A=-5 -3;b=-45;Aeq=;beq=;vlb = zeros(2,1);vub=9;15; %調(diào)用linprog函數(shù)：x,fval = linprog(c,A,b,Aeq,beq,vlb,vub)結(jié)果為：x = 9.0000 0.0000fval =360即只需聘用9個一級檢驗員。Matlab優(yōu)化工具箱簡介1.MATLAB求解優(yōu)

6、化問題的主要函數(shù)2.優(yōu)化函數(shù)的輸入變量使用優(yōu)化函數(shù)或優(yōu)化工具箱中其它優(yōu)化函數(shù)時, 輸入變量見下表:3. 優(yōu)化函數(shù)的輸出變量下表:4控制參數(shù)options的設(shè)置Options中常用的幾個參數(shù)的名稱、含義、取值如下:(1)Display: 顯示水平.取值為off時,不顯示輸出; 取值為iter時,顯示每次迭代的信息;取值為final時,顯示最終結(jié)果.默認(rèn)值為final.(2)MaxFunEvals: 允許進行函數(shù)評價的最大次數(shù),取值為正整數(shù).(3) MaxIter: 允許進行迭代的最大次數(shù),取值為正整數(shù)控制參數(shù)options可以通過函數(shù)optimset創(chuàng)建或修改。命令的格式如下：(1) optio

7、ns=optimset(optimfun) 創(chuàng)建一個含有所有參數(shù)名,并與優(yōu)化函數(shù)optimfun相關(guān)的默認(rèn)值的選項結(jié)構(gòu)options.（2）options=optimset(param1,value1,param2,value2,.) 創(chuàng)建一個名稱為options的優(yōu)化選項參數(shù),其中指定的參數(shù)具有指定值,所有未指定的參數(shù)取默認(rèn)值.(3)options=optimset(oldops,param1,value1,param2, value2,.) 創(chuàng)建名稱為oldops的參數(shù)的拷貝,用指定的參數(shù)值修改oldops中相應(yīng)的參數(shù).例：opts=optimset(Display,iter,TolFun

8、,1e-8) 該語句創(chuàng)建一個稱為opts的優(yōu)化選項結(jié)構(gòu),其中顯示參數(shù)設(shè)為iter, TolFun參數(shù)設(shè)為1e-8.用Matlab解無約束優(yōu)化問題 一元函數(shù)無約束優(yōu)化問題常用格式如下：（1）x= fminbnd (fun,x1,x2)（2）x= fminbnd (fun,x1,x2 ，options)（3）x，fval= fminbnd（.）（4）x，fval，exitflag= fminbnd（.）（5）x，fval，exitflag，output= fminbnd（.）其中（3）、（4）、（5）的等式右邊可選用（1）或（2）的等式右邊。 函數(shù)fminbnd的算法基于黃金分割法和二次插值法，它

9、要求目標(biāo)函數(shù)必須是連續(xù)函數(shù)，并可能只給出局部最優(yōu)解。例1 求在0x 0，且a11 a12；同理， p2=b2-a21x1-a22x2，b2，a21，a22 02成本與產(chǎn)量成負(fù)指數(shù)關(guān)系甲的成本隨其產(chǎn)量的增長而降低,且有一個漸進值,可以假設(shè)為負(fù)指數(shù)關(guān)系,即:同理， 模型建立總利潤為： z(x1,x2)=(p1-q1)x1+(p2-q2)x2若根據(jù)大量的統(tǒng)計數(shù)據(jù),求出系數(shù)b1=100,a11=1,a12=0.1,b2=280,a21=0.2,a22=2,r1=30,1=0.015,c1=20, r2=100,2=0.02,c2=30,則問題轉(zhuǎn)化為無約束優(yōu)化問題：求甲,乙兩個牌號的產(chǎn)量x1，x2，使總

10、利潤z最大.為簡化模型,先忽略成本,并令a12=0,a21=0,問題轉(zhuǎn)化為求: z1 = ( b1 - a11x1 ) x1 + ( b2 - a22x2 ) x2 的極值. 顯然其解為x1 = b1/2a11 = 50, x2 = b2/2a22 = 70,我們把它作為原問題的初始值.模型求解1.建立M-文件fun.m: function f = fun(x) y1=(100-x(1)- 0.1*x(2)-(30*exp(-0.015*x(1)+20)*x(1); y2=(280-0.2*x(1)- 2*x(2)-(100*exp(-0.02*x(2)+30)*x(2); f=-y1-y2;

11、2.輸入命令: x0=50,70; x=fminunc(fun,x0), z=fun(x)3.計算結(jié)果: x=23.9025, 62.4977, z=6.4135e+003 即甲的產(chǎn)量為23.9025,乙的產(chǎn)量為62.4977,最大利潤為6413.5.非線性規(guī)劃1、 二次規(guī)劃用MATLAB軟件求解,其輸入格式如下: 1.x=quadprog(H,C,A,b); 2.x=quadprog(H,C,A,b,Aeq,beq); 3.x=quadprog(H,C,A,b,Aeq,beq,VLB,VUB); 4.x=quadprog(H,C,A,b, Aeq,beq ,VLB,VUB,X0); 5.x=

12、quadprog(H,C,A,b, Aeq,beq ,VLB,VUB,X0,options); 6.x,fval=quaprog(.); 7.x,fval,exitflag=quaprog(.); 8.x,fval,exitflag,output=quaprog(.);例1 min f(x1,x2)=-2x1-6x2+x12-2x1x2+2x22 s.t. x1+x22 -x1+2x22 x10, x20 1、寫成標(biāo)準(zhǔn)形式：s.t.2、 輸入命令： H=1 -1; -1 2; c=-2 ;-6;A=1 1; -1 2;b=2;2; Aeq=;beq=; VLB=0;0;VUB=; x,z=qu

13、adprog(H,c,A,b,Aeq,beq,VLB,VUB)3、運算結(jié)果為： x =0.6667 1.3333 z = -8.2222一般非線性規(guī)劃標(biāo)準(zhǔn)型為：min F(X) s.t AX=bG(X)Ceq(X)=0 VLBXVUB其中X為n維變元向量，G(X)與Ceq(X)均為非線性函數(shù)組成的向量，其它變量的含義與線性規(guī)劃、二次規(guī)劃中相同.用Matlab求解上述問題，基本步驟分三步：1. 首先建立M文件fun.m,定義目標(biāo)函數(shù)F（X）:function f=fun(X);f=F(X);2. 若約束條件中有非線性約束:G(X)或Ceq(X)=0,則建立M文件nonlcon.m定義函數(shù)G(X)

14、與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= fm

15、incon(.)(7) x,fval,exitflag= fmincon(.) (8)x,fval,exitflag,output= fmincon(.)注意：1 fmincon函數(shù)提供了大型優(yōu)化算法和中型優(yōu)化算法。默認(rèn)時，若在fun函數(shù)中提供了梯度（options參數(shù)的GradObj設(shè)置為on），并且只有上下界存在或只有等式約束，fmincon函數(shù)將選擇大型算法。當(dāng)既有等式約束又有梯度約束時，使用中型算法。2 fmincon函數(shù)的中型算法使用的是序列二次規(guī)劃法。在每一步迭代中求解二次規(guī)劃子問題，并用BFGS法更新拉格朗日Hessian矩陣。3 fmincon函數(shù)可能會給出局部最優(yōu)解，這與初值

16、X0的選取有關(guān)。例2s.t.1、寫成標(biāo)準(zhǔn)形式： s.t. 2、先建立M-文件 fun3.m: function f=fun3(x); f=-x(1)-2*x(2)+(1/2)*x(1)2+(1/2)*x(2)23、再建立主程序youh2.m： x0=1;1; A=2 3 ;1 4; b=6;5; Aeq=;beq=; VLB=0;0; VUB=; x,fval=fmincon(fun3,x0,A,b,Aeq,beq,VLB,VUB)4、運算結(jié)果為： x = 0.7647 1.0588 fval = -2.0294例31先建立M文件 fun4.m,定義目標(biāo)函數(shù): function f=fun4(

17、x); f=exp(x(1) *(4*x(1)2+2*x(2)2+4*x(1)*x(2)+2*x(2)+1);2再建立M文件mycon.m定義非線性約束： function g,ceq=mycon(x) g=x(1)+x(2);1.5+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10;3主程序youh3.m為:x0=-1;1;A=;b=;Aeq=1 1;beq=0;vlb=;vub=;x,fval=fmincon(fun4,x0,A,b,Aeq,beq,vlb,vub,mycon)3. 運算結(jié)果為：x = -1.2250 1.2250 fval = 1.8951例4資金使用問題設(shè)有400萬元資金, 要求4年內(nèi)使用完, 若在一年內(nèi)使用資金x萬元, 則可得效益萬元(效益不能再使用),當(dāng)年不用的資金可存入銀行, 年利率為10%. 試制定出資金的使用計劃, 以使4年效益之和為最大.設(shè)變量表示第i年所使用的資金數(shù),則有1先建立M文件 fun44.m,定義目標(biāo)函數(shù):function f=fun44(x)f=-(sqrt(x(1)+sqrt(x(2)+sqrt(x(3)+sqrt(x(4);