优化问题的MATLAB求解总汇.ppt

* 1)先建立M文件 fun4.m,定义目标函数: function f=fun4(x); f=exp(x(1)) *(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1); x1+x2=0 s.t. 1.5+x1x2 - x1 - x2 0 -x1x2 –10 0 例 2)再建立M文件mycon.m定义非线性约束： function [g,ceq]=mycon(x) g=[1.5+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10]; ceq=[]; * 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) 4)运算结果为： x = -1.2250 1.2250 fval = 1.8951 The end! 练习 1.（连续投资问题）某部门有现金10万元，在今后五年内考虑给下面项目投资，已知： 项目A，从第一年到第四年每年年初需要投资，并于次年末回收本利115%; 项目B，第三年初需要投资，到第五年末回收本利125%，但最大投资额不超过4万元； 项目C，第二年初需要投资，到第五年末回收本利140%，但最大投资额不超过3万元； 项目D，五年内每年初可购买公债，于当年末归还，并加利息6%； 问：如何投资到第五年末时收益最大？ 练习 2.（供应与选址问题）某公司有6个建筑工地要开工，每个工地的位置（用平面坐标系a，b表示，距离单位：千米 ）及水泥日用量d(吨)由下表给出。目前有两个临时料场位于A(5,1)，B(2,7)，日储量各有20吨。假设从料场到工地之间均有直线道路相连。 （1）试制定每天的供应计划，即从A，B两料场分别向各工地运送多少吨水泥，使总的吨千米数最小。 （2）为了进一步减少吨千米数，打算舍弃两个临时料场，改建两个新的，日储量各为20吨，问应建在何处，节省的吨千米数有多大？ 优化问题的MATLAB求解 一、matlab优化工具箱中常用的功能函数 求解线性规划问题的主要函数是linprog。 求解二次规划问题的主要函数是quadprog。 求解无约束非线性规划问题的主要函数是fminbnd、fminunc和fminsearch。 求解约束非线性规划问题的主要函数是fgoalattain和fminimax。 二、一般步骤 建立目标函数文件 针对实际问题建立优化设计的数学模型 不等式约束条件表示成g(X)≥0的形式 建立调用优化工具函数的命令文件 文件内容：必须的输入参数、描述标函数表达式等 存储：以自定义的目标函数文件名存储在文件夹中 建立约束函数文件 文件内容：必须的输入参数、约束函数表达式等 存储：以自定义的约束函数文件名存储在文件夹中 将优化设计的命令文件复制到MATLAB命令窗口中进行运算求解 分析优化设计的数学模型，选择适用的优化工具函数 文件内容：初始点，设计变量的边界约束条件， 运算结果输出等内容 存储：以自定义的命令文件名存储于文件夹中。 线性规划问题 三、线性规划数学模型 1.主要应用对象： （1）在有限的资源条件下完成最多的任务； （2）如何统筹任务以使用最少资源。 2.数学模型形式： min f TX s.t. AX≤b （线性不等式约束条件） AeqX=beq （线性等式约束条件） lb ≤X ≤ub （边界约束条件） 约束条件 决策变量 目标函数 非负数 线性 3.MATLAB中函数调用格式 [xopt, fopt]=linprog( f, A, b, Aeq, beq, lb, ub, x0, options) 最优解 最优值 目标函数各维变量系数向量 初始点 可选项 四、线性规划例题 生产规划问题：某厂利用a,b,c三种原料生产A,B,C三种产品，已知生产每种产品在消耗原料方面的各项指标和单位产品的利润，以及可利用的数量，试制定适当