非线性规划11.pptVIP

min z=5x1+6x2+-12x3+2x1^2+4x2^2+6x3^2 - 2x1x2-6x1x3+8x1x3 s.t. x1+2x2+x3≥6; x1+x2+x3≤16; -x1+2x2≤4; x1≥0, x2≥0, x3≥0 . 写成矩阵形式: min 0.5 x’Hx+c’x s.t. Ax≤b; x≥0 . 第31页/共54页 H=[4 -2 -6; -2 8 8; -6 8 12]; f=[5 6 -12]’; A=[-1 -2 -1; 1 1 1; -1 2 0]; b=[-6 16 4]’; lb=[0 0 0]; [x,z]=quadprog(H, f, A, b, [], [], lb); 或options=optimset(largescale,off); [x,z]=quadprog(H, f, A, b, [], [], lb,[],[],options) x = 3.2500 0 2.7500 z = -3.8750 输出结果： 可能会提示你改用中等规模的算法，更正如下 第32页/共54页 3、约束非线性规划（多变量）情形 标准模型 Min f(X) s.t. c(X) ≤0, ceq(X)=0 AX ≤b, Aeq.X=beq, lb ≤X ≤ub 求解格式 第33页/共54页 ① 建立m文件函数fun.m function f=fun(x) f=f(x); ③ 为函数fmincon的其余输入变量赋值，然后调用该函数求出约束规划问题的解。 解释 ② 建立m文件函数nonlcon.m function [c,ceq]= nonlcon(x) c=c(x);ceq=ceq(x) 第34页/共54页 例：求解以下约束非线性规划： Max f(x) = x1 x2 s.t. 2(x1 + x2) x3 ≤500 x3≥2 xj≥0，j = 1,2 第35页/共54页 MATLAB程序 ① function f=fun2(x) f=-x(1)*x(2); ② function [c,ceq]=nlcon(x) c=(x(1)+x(2))*x(3)-250; ceq=[]; ③ x0=[10 10 2]; L=[0 0 2]; [x,fval]=fmincon(fun2,x0,[],[],[],[],L,[], nlcon) 第36页/共54页 计算结果 x = 62.5000 62.5000 2.0000 fval = -3.9063e+003 第37页/共54页 投资决策问题 一般模型与算法概述 Matlab软件求解简介 实 验 内 容 主 要 内 容 范例：选址问题 第38页/共54页 供应与选址 某公司有6个建筑工地要开工，每个工地的位置（用平面坐标a，b表示，距离单位：千米）及水泥日用量d 吨由下表给出。目前有两个临时料场位于A(5,1), B(2,7)，日储量各有20吨。假设从料场到工地均有直线道路相连， 问题1：试制定每天A、B 两料场向各工地供应水泥的供应计划，使总的吨千米数最小。 a 1.25 8.75 0.5 5.75 3 7.25 b 1.25 0.75 4.75 5 6.5 7.75 d 3 5 4 7 6 11 第39页/共54页 问题2： 为进一步减少吨千米数，打算舍弃两个临时料场，改建两个新的，日储量仍各为20吨，问应建在何处，节省的吨千米数有多大？ 第40页/共54页 建立规划模型 记工地的位置为(ai, bi)，水泥日用量为di, i = 1,…,6, 料场位置为(xj, yj), 日储量为ej, j=1,2；从料场j向工地i的运送量为cij。 （工地日用量） （料场日储量） 第41页/共54页 问题1：使用临时料场，即料场位置 (xj, yj)为已知,决策变量为cij，上述模型为线性规划模型。 记决策变量X=[c11,c21,…,c61,c12,…c62] a0=[1.25 8.