八优化模型总汇.ppt

function [c,ceq,GC,GCeq] = mycon(x) c = ... % Nonlinear inequalities at x ceq = ... % Nonlinear equalities at x if nargout 2 % nonlcon called with 4 outputs GC = ... % Gradients of the inequalities GCeq = ... % Gradients of the equalities end * * 第8章 优化模型 8.1 无约束最优化问题 8.2 线性规划问题 8.3 二次规划问题 8.4 非线性规划问题 8.5 0-1规划问题 8.6 圆形工件检验优化模型 8.7 钢管的订购与运输 8.1 无约束最优化问题 数学描述 优化变量 目标函数 matlab解 一元函数极小 X=fminbnd(‘f’,x1,x2) 无约束极小 X=fminunc(‘f’,x0) (牛顿法) X=fminsearch(‘f’,x1,x2) 注： 得到的只是局部最优解，并非全局最优解 解: f=2*exp(-x).*sin(x); fplot(f,[0,8]); [xmin,ymin]=fminbnd (f, 0,8) f1=-2*exp(-x).*sin(x); [xmax,ymax]=fminbnd (f1, 0,8) 例2 对边长为3米的正方形铁板，在四个 角剪去相等的正方形以制成方形无盖 水槽，问如何剪法使水槽的 容积最大？ 例 1 求 f=2 x e x sin - 在 0 x 8 中的最小值与最大值 解: 设剪去的正方形的边长为x, ，则水槽的容积为： 建立无约束优化模型为： , 0x1.5 先编写M文件如下: function f=myfun(x) f=-(3-2*x).^2*x; 调用fminbnd: [x,fval]=fminbnd(myfun,0,1.5) 运算结果为: x = 0.5000,fyal =2.0000. 即剪掉的正方形的边长为0.5米时水槽的容积最大,最大容积为2立方米. 例3(初值的影响力)设目标函数为 试观察不同的初值得出的最小值。 f=inline(exp(-2*t).*cos(10*t)+exp(-3*(t+2).*sin(2*t)),t) t0=1;[t1 f1]=fminsearch(f,t0) t0=0.1;[t2 f2]=fminsearch(f,t0) syms t; y=exp(-2*t)*cos(10*t)+exp(-3*(t+2))*sin(2*t); ezplot(y,[0,2.5]); axis([0 2.5 -0.6 1]) 数形结合 目标函数 约束条件 如果目标函数和约束条件都是线性函数， 则该模型称为线性规划. 决策变量 8.2 线性规划问题 数学描述 M m eq eq x x x x A B Ax to sub ￡ ￡ = ￡ B : . z = min T f x matlab解 [x,f_opt,flag,c]= linprog(f,A,B,Aeq,Beq,xm,xM,x0,opt) X: 解 f_opt: 最优值 Flag：大于零表示求解成功，否则求解出问题 C：求解信息 X0：搜索点的初值 Opt：最优化控制项 Optimset常用算法选项 约束条件的精度 tolcon 解的精度 tolx 误差函数精度 tolfun 最大允许的迭代次数 maxiter 是否使用大规模问题算法， on,off largescale 目标函数的梯度是否已知，on,off Gradobj 中间结果显示方式off,iter,final display 例4：某车间生产A和B两种产品，为了生产A和B，所需的原料分别为2个和3个单位，所需的工时分别为4个和2个单位。现在可以应用的原料为100个单位，工时为120个单位。每生产一台A和B分别可获得利润6元和4元。应当生产A和B各多少台能获得最大利润？ 总原料100个单位； 总工时120个单位。 分析： 一台A 一台B 2个原料 4个工时 赚6元钱 3个原料 2个工时 赚4元钱 解：设生产A产品 台，生产B产品 台 f=[-6,-4]; A=[2 3;4 2