Matlab求解优化问题.docxVIP

预备知识：M文件简介 在MATLAB中，用户能够利用 Edtior（编写器）成立M文件，此后在命令窗口中的“ ”提示符 下键入M文件的主文件名，回车履行. MATLAB的M文件有两类：命令文件和函数文件。将原来要在MATLAB环境下直接输入的语句，放 在一个以.m 为后缀的文件中，这一文件就称为命令文件；函数文件由五部分构成：函数定义行、 H1 行、函数帮助文本、函数体、说明， MATLAB的内部函数都是由函数文件定义的。 1.11优化(最值、数学规划) 在数学上，优化问题包括最值问题和数学规划问题等，后者又包括线性规划、整数规划（含 0-1 规划）、二次规划等. 在MATLAB中，求解最值问题的命令主要有： fminbnd (f,x1,x2) 求一元函数f在区间[x1,x2] 上的最小值点 [x,fval]=fminbnd(f,x1,x2) 求一元函数f 在区间[x1,x2] 上的最小值点和最小值 fminsearch(’f’,x0) 求多元函数 f在点x0周边的最小值点 [x,fval]=fminsearch( ’f’,x0) 求多元函数f 在点x0周边的最小值点和最小值 例 求函数f(x) x2 3x 2在区间[5,5]上的最小值点和最小值. [x,fval]=fminbnd(x^2+3*x+2,-5,5) x=-1.5000 fval=-0.2500 例 2 2 求函数 f(x1,x2)x1x2 在点(1,1)周边的最小值点和最小值. x1 x2 [x,fval]=fminsearch(x(1)*x(2)+2/x(1)+2/x(2),[11]) x= 1.25991.2599fval= 4.7622 在MATLAB中，求解数学规划问题的命令主要有： min z cTx s.t. Ax b （1）线性规划 Aeqx beq lb x ub 命令：[x,fval]=linprog (c,A,b,Aeq,beq,lb,ub) 在上述命令中，当某些参数空缺时，可用 []代替或省略，下同。 min z 2x1 3x2 6x3 5x4 s.t. x1 x2 2x3 4x4 0 例 求解线性规划问题 x2 x3 x4 0. x1 x2 x3 x4 1 x1,x2,x3,x40 c=[2;3;6;5]; A=[1-1-2-4;0-1-11]; b=[0;0]; Aeq=[1111]; beq=[1]; lb=[0;0;0;0]; [x,fval]=linprog(c,A,b,Aeq,beq,lb,[])Optimizationterminated. x=0.5000 0.5000 0.0000 0.0000 fval= 2.5000 minzcTx 2）0-1规划s.t.Axb Aeqxbeq x 0,1 命令：[x,fval]= bintprog (c,A,b,Aeq,beq) min zx1 1.2x2 0.8x3 s.t. 2.1x1 2x2 1.3x3 5 例 求解0-1 规划问题 0.8x1 x2 5 x1 2.5x2 2x3 . 8 2x2 8 x1,x2,x3 0,1 c=[-1;-1.2;-0.8]; A=[2.1,2,1.3;0.8,1,0;1,2.5,2;0,2,0]; b=[5;5;8;8]; [x,fval]=bintprog(c,A,b,[],[])Optimizationterminated. x= 1 1 0 fval= -2.2000 minz1xTHxcTx 2 3）二次规划s.t.Axb Aeqx beq lb x ub 命令：[x,fval]=quadprog (H,c,A,b,Aeq,beq,lb,ub) min z 1x12 x22 x1x2 2x16x2 2 例 求解二次规划问题 s.t. x1 x2 2 . x1 2x2 2 2x1x2 3 x1,x2 0 z 1x12 x22 x1x22x16x2 1 x1 2 2 H=[1-1;-12]; c=[-2;-6]; A=[11;-12]; b=[22]; Aeq=[21]; beq=[3]; lb=[0;0]; [x,fval]=quadprog(H,c,A,b,Aeq,beq,lb,[]) Warning:Large-scalemethodformulation, switchingtomedium-scalemethod. Inquadprogat242Optimizationterminated.x= 1.0000 1.0000fval= -7.5000  x2 1 1 x1 x1 1 2 x2 26 x2 doesnotcurrentlysolvethisproblem min z f(x) s.t. Ax b beq （4）非