MATLAB优化工具箱__线性规划,非线性规划.pptVIP

  • 4
  • 0
  • 约2.74千字
  • 约 19页
  • 2018-12-28 发布于浙江
  • 举报

MATLAB优化工具箱__线性规划,非线性规划.ppt

例题建模 典型的优化问题 三大要素: 目标:种植A,B两种作物获得最大利润 决策:安排种植A,B两种作物各所少亩 约束:肥料限制;种植土地面积限制 建立模型 max 6x1+4x2 s.t. 2x1+5x2 ≤100 4x1+2x2 ≤120 Matlab优化工具箱 线性规划:linprog 非线性规划:fminbnd,fminsearch,fmincon Matlab求解线性规划模型 函数linprog linprog 语法: 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,options) [x,fval] = linprog(...) [x,fval,exitflag] = linprog(...) [x,fval,exitflag,output] = linprog(...) [x,fval,exitflag,output,lambda] = linprog(...) linprog 例题的求解程序 模型: max 6x1+4x2 s.t. 2x1+5x2 ≤100 4x1+2x2 ≤120 Matlab求解程序: A=[2 5;4 2]; b=[100 120]; f=-[6 4]; [optx ,funvalue,exitflag]=linprog(f,A,b,[],[],[0 0],[inf,inf]) 程序运行结果 输出: Optimization terminated successfully. optx = 25.0000 10.0000 funvalue = -190.0000 exitflag = 1 解释得出实际问题的解 当分别种植A、B两种作物为25亩、10亩时,预计共获得利润190(百元)。 fmincon函数求解形如下面的有约束非线性规划模型 一般形式: 输入参数语法: x = fmincon(fun,x0,A,b) x = fmincon(fun,x0,A,b,Aeq,beq) x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub) x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon) x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options) x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options,P1,P2, ...) [x,fval,exitflag,output,lambda]=fmincon(fun,x0,...) 输入参数的几点说明 模型中如果没有A,b,Aeq,beq,lb,ub的限制,则以空矩阵[ ]作为 参数传入; nonlcon:如果包含非线性等式或不等式约束,则将这些函数 fmincon示例 求解步骤: (1)编写目标函数文件(这里文件名用myobjfun2.m) (2)编写约束条件函数文件(mymodelcons.m) (3)编写调用fmincon主程序(mymain2.m) fmincon示例程序 学习小结 最优化问题建模的关键是先要确定三要素,再转化为数学表达式(数学模型)。 学习中既要初步掌握最优化问题的建模步骤,也要善于运用Matlab的优化工具箱求解优化模型。 有些模型可以采用多个Matlab函数求解,可以比较结果,加深认识。 思考题 一幢楼房的后面是一个很大的花园。在花园中紧靠着楼房建有一个温室,温室高10英尺,延伸进花园7英尺。清洁工要打扫温室上方的楼房的窗户。他只有借助于梯子,一头放在花园中,一头靠在楼房的墙上,攀援上去进行工作。他只有一架20米长的梯子,你认为他能否成功?能满足要求的梯子的最小长度是多少?请就以上问题建立数学模型,并编程求解。 提示: * 数学实验 例子:某农场种植两种作物A、B,需要甲、乙两种化肥。种植每亩作物A和作物B分别需用的化肥数,可得利润及农场现有化肥数量如下表所示: 问在现有条件下,如何安排种植,才能使利润最大? 求解下列形式的线性规划模型: linprog输入参数说明: f, A, b, Aeq, beq lb,ub 边界设置 说明: 如果x(i)无边界,则 lb(i) = -inf, ub(i) = inf 输出参数说明: x 决策变量取值 fval 目标函数最优值 exitflag 0 成功找到最优解 0 达到最大迭代次数

文档评论(0)

1亿VIP精品文档

相关文档