matlab第九章 数学规划求解.pptVIP

第九章 数学规划求解 线性规划、二次规划、无约束非线 性规划、有约束非线性规划、极大极小化问题 一、线性规划 基本模型 % min f(x) % s.t .(约束条件)： Ax=b % (等式约束条件)： Aeqx=beq % lb=x=ub 调用格式 x=linprog(f,A,b,Aeq,beq)x=linprog(f,A,b,Aeq,beq,lb,ub)x=linprog(f,A,b,Aeq,beq,lb,ub,x0)x=linprog(f,A,b,Aeq,beq,lb,ub,x0)[x,fval]=linprog(f,A,b,Aeq,beq,lb,ub,x0)[x,fval,exitflag]=linprog(f,A,b,Aeq,beq,lb,ub,x0,)[x,fval,exitflag,output]=linprog(f,A,b,Aeq,beq,lb,ub,x0)[x,fval,exitflag,output,lambda]=linprog(f,A,b,Aeq,beq,lb,ub,x0) 说明：f:是优化参数x的系数矩阵：lb,ub:设置优化参数x的上下界：fval:返回目标函数在最优解x点的函数值：exitflag:返回算法的终止标志：output:返回优化算法的信息的一个数据结构。 例：求使函数f(x)=-5x1-4x2-6x3取最小值的x值，且满足约束条件：x1-x2+x3=203x1+2x2+4x3=423x1+2x2=30x1=0,x2=0,x3=0求解的MATLAB代码：f=[-5;-4;-6];A=[1,-1,1;3,2,4;3,2,0];b=[20;42;30];lb=zeros(3,1);[x,fval]=linprog(f,A,b,[],[],lb) 例1： c=[-1,-3]; A=[-1,1;1,1;2,-1]; b=[3;7;8]; lb=[0;0]; ub=1e+10*[1;1]; x=linprog(c,A,b,[],[],lb,ub) x = 2.0000 5.0000 f=c*x f = -17.0000 例2： c=[-0.5,-1,-0.5,-1]; A=[1,1,1,1;0.2,0.5,1,2;2,1,0.5,0.2]; b=[10;10;10];Aeq=[1,1,1,-2];beq=6; lb=[0;0;0;0]; ub=1e+10*[1;1;1;1]; x0=[0;0;0;0]; x=linprog(c,A,b,Aeq,beq,lb,ub,x0) x = 0.0000 8.6667 0.0000 1.3333 f=c*x f = -10.0000 二、二次规划 二次规划的标准形式为：min??f=(1/2)XTHX+fTX约束条件：A.x=b，Aeq.x=beq，lb?x?ub其中：c、b、beq、lb、ub、x是矢量，H、A、Aeq为矩阵。在MATLAB中可以使用quadprog函数来求最小值。 例3：计算下面二次规划问题 minf（x）= (1/2)x1^2+x2^2- x1x2-2x1-6x2约束条件： x1+x2 ?2；-x1+2x2 ?2 ；2x1+x2 ?3；0? x1，0? x2 解：把二次规划问题写成标准形式：(1/2)XTHX+fTX这里：H=? 1? -1???? f=?? -2?????? X=? x1??????????????? -1?? 2???????? ? -6????????????? x2 H=[1 -1;-1 2];f=[-2;-6];A=[1 1;-1 2;2 1];b=[2;2;3];lb=zeros(2,1);[x,fval,exitflag,output,lambda]=quadprog(H,f,A,b,[],[],lb) 程序与编程 指定次数重复循环语句 For vi=s1:s3:s2 Sta End 不定次数重复循环语句 While va Sta End 条件转移语句 简单if语句 If exp Sta End 双目if语句 If exp Sta1 Else Sta2 end 嵌套if语句 If exp1 Sta1 Elseif exp2 Sta2 Elseif expn Stan Else Stan+1 End 开关切换语句 Switch kv Case kv1 Sta1 Case kv2 Sta2 Case kvn Stan Otherwise Break end 其他语句 Break Return pause * * x