《数学实验》8.pptVIP

*/16 一元函数极小值问题 线性规划问题标准形式 线性规划问题求解命令 思考题与练习题 《数学实验》8 最优化方法——从可行方案中寻求系统的最优方案 给定条件下,按某一指标寻找最优方案。可表示为求函数在约束条件下的极大值(或极小值)问题。 无约束极值问题的数学模型 有约束线性规划的数学模型 决策变量: X = [x1, x2, ···, xn ]T 目标函数系数: C=[c1, c2, ···, cn]T 不等式约束矩阵和向量: A, b 温室 例1 梯子问题:在花园靠楼房处有一花棚,花棚伸入花 解:设梯子长度为L, 梯子与地面的 夹角为 数学模型: L=inline(2./cos(alpha)+3./sin(alpha)) ezplot(L,0,pi/2) [x,Lmin]=fminbnd(L,0.8,0.9) x = 0.8528 Lmin = 7.0235 梯子长度基本可行。 园2米,高3米,花棚上方是楼房的窗台,要将梯子从花 园地上放靠在楼房墙上不碰坏花棚,现用7米长的梯子 是否可行？ 无约束优化问题 min f(x) x1 x x2 Matlab求解命令使用格式: x= fminbnd (fun,x1,x2) 或 [x, fval]= fminbnd(fun,x1,x2) 例2 求 f(x) = e-x sin x 在[0, 8]区间中的最大值和最小值 f=inline(exp(-x).*sin(x)); f1=inline(-exp(-x).*sin(x)); ezplot(f,0,8) [xmin,ymin]=fminbnd (f, 0,8) xmax=fminbnd (f1, 0,8) ymax=f(xmax) xmin =3.9270，ymin = -0.0139 xmax =0.7854，ymax = 0.3224 例3.边长3米的正方形铁板,在四个角剪去相等小正方形制成无盖水槽,问如何使水槽容积最大？ 解: 设小正方形边长为 x ,则水槽容积 V=(3 – 2x)2x 最大值问题: max (3 – 2x)2x 最小值问题: min – (3 – 2x)2x f=inline(-x.*(3-2*x).^2); [x,Fmin]=fminbnd(f,0,1.5) xmin=x,Fmax=-Fmin xmin = 0.5000 Fmax = 2.0000 即剪掉的正方形的边长为0.5米时水槽的容积最大,最大容积为2立方米 例4 某建筑公司投标修建办公楼和住宅楼,已知建办公楼利润是50元/平方米，住宅楼利润60元/平方米，要求完成建筑面积不少于5000m2，其中办公楼的面积不能大于5000 m2，住宅楼不能大于3000m2。 假定公司当年建办公楼x1平方米，建住宅楼x2平方米。以所得利润最大为目标，得目标函数 z = 50x1 + 60x2 根据招标单位的要求,得约束条件 x1+x2≥ 5000 x1 ≤ 5000 x2 ≤ 3000 x1≥ 0, x2≥ 0 求解线性规划命令使用格式 (1) x=linprog(C, A, b) [x,fval] = linprog(C, A, b) (2) x=linprog(c,A,b,Aeq,beq) [x,fval]=linprog(c,A,b,Aeq,beq) (3) x=linprog(c,A,b,Aeq,beq,e0,e1) MATLAB中线性规划问题标准形式 C——目标函数系数向量 A——约束条件系数矩阵 b——约束条件常数向量 解线性规划问题标准形式 min -( 50x1+ 60x2) s. t. -x1-x2≤-5000 x1 ≤5000 x2≤3000 x1≥0, x2≥0 ? C = – [ 50 60] C=[50,60]; A=[-1,-1;1,0;0,1;]; b=[-5000;5000;3000]； x=linprog(-C,A,b) z=C*x x = 5000 3000 z= 430000 程序: 结果: ①一组决策变量(x1, x2, ···, xn)表示某套方案; ②一组线性不等式或线性等式为约束条件; ③以决策变量的线性函数作为目标函数. 明确问题的目标 假设一组决策变量 考虑目标函数 考虑约束条件 建立线性规划问题数学模型: 例5 某厂生产两种产