【课件】优化建模.ppt

You can evaluate the constraints at the solution by entering [c,ceq] = confun(x) This returns very small numbers close to zero, such as c = 1.0e-007 * -0.9032 0.9032 ceq = [] Note that both constraint values are, to within a small tolerance, less than or equal to zero; that is, x satisfies 微分方程组数值解：食饵-捕食者模型 xdot=diag([r-a*x(2),-d+b*x(1)])*x 1. 求解与绘图 ts=0:0.1:15; x0=[25,2]; [t,x]=ode45('shier',ts,x0); [t,x] plot(t,x),grid,gtext('x1(t)'),gtext('x2(t)'), {plot(x(:,1),x(:,2)),grid,xlabel('x1'),ylabel('x2')} 2. M-文件 shier.m function xdot=shier(t,x) r=1;d=0.5;a=0.1;b=0.02; xdot=diag([r-a*x(2),-d+b*x(1)])*x; X(t)-图 x1-x2 图 例1 加工奶制品的生产计划 1桶牛奶 3公斤A1 12小时 8小时 4公斤A2 或 获利24元/公斤 获利16元/公斤 50桶牛奶 时间480小时 至多加工100公斤A1 制订生产计划，使每天获利最大 35元可买到1桶牛奶，买吗？若买，每天最多买多少? 可聘用临时工人，付出的工资最多是每小时几元? A1的获利增加到 30元/公斤，应否改变生产计划？ 每天： 1桶牛奶 3公斤A1 12小时 8小时 4公斤A2 或 获利24元/公斤 获利16元/公斤 x1桶牛奶生产A1 x2桶牛奶生产A2 获利 24×3x1 获利 16×4 x2 原料供应 劳动时间 加工能力 决策变量 目标函数 每天获利 约束条件 非负约束 线性规划模型(LP) 时间480小时 至多加工100公斤A1 50桶牛奶 每天 Max 72x1+ 64x2 s. t. x1+ x2≤50 12x1+ 8x2≤480 3x1 ≤100 x1, x2 ≥0 MATLAB程序: c=-[72,64]; a=[1,1;12,8;3,0]; b=[50,480,100]; v1=[0,0]; x=lp(c,a,b,v1) z=-c*x 解：x =[20.0000,30.0000]；z =3360 20桶牛奶生产A1, 30桶生产A2，利润3360元。 LINDO 公司软件产品简要介绍 美国芝加哥(Chicago)大学的Linus Schrage教授于1980年前后开发, 后来成立 LINDO系统公司（LINDO Systems Inc.）， 网址： LINDO: Linear INteractive and Discrete Optimizer (V6.1) LINGO: Linear INteractive General Optimizer (V8.0) LINDO API: LINDO Application Programming Interface (V2.0) What’s Best!: (SpreadSheet e.g. EXCEL) (V7.0) 演示(试用)版、学生版、高级版、超级版、工业版、扩展版… （求解问题规模和选件不同） 优化模型 实际问题中 的优化模型 x~决策变量 f(x)~目标函数 gi(x)?0~约束条件 数学规划 线性规划(LP) 二次规划(QP) 非线性规划(NLP) 纯整数规划(PIP) 混合整数规划(MIP) 整数规划(IP) 0-1整数规划 一般整数规划 连续规划 LINDO和LINGO软件能求解的优化模型