第7章 方程求根与最优化.pptVIP

*/16 */16 第六章 方程求根与最优化 一元方程求根命令 线性规划问题标准形式 线性规划问题求解命令 一元函数极小值问题 ①求多项式零点命令: roots(p) ②求一元函数零点命令: fzero(fun, x0) 多项式 P(x) = a1xn + a2xn-1 + ······ + anx + an+1 系数提取 p =[a1 , a2 , ······ , an+1 ] p=[1 -30 0 2552]; roots(p) ans = 26.3146 11.8615 -8.1761 r x 例1. 球体的吃水深度计算 半径 r =10 cm的球体,密度 ? =0.638.浸入水深度 x = ? 解：重量 体积 x3 –30x2 +2552 =0 fun=inline(x.^3-30*x.^2+2552); d=fzero(fun,10) d = 11.8615 1、一元方程求根命令 x y 例2. 求 x = tan x 正根 f=inline(x-tan(x)) x=fzero(f,3.5) function y=fx(x) y=x-tan(x); x=fzero(fx,3.5) 或 2、最优化方法——从可行方案中寻求系统的最优方案 无约束极值问题的数学模型 线性规划的数学模型 决策变量: X = [x1, x2, ···, xn ]T 目标函数系数: C=[c1, c2, ···, cn]T 不等式约束矩阵和向量: A, b 例2 某建筑公司投标修建办公楼和住宅楼,已知建办公楼利润是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)表示某套方案; ②一组线性不等式或线性等式为约束条件; ③以决策变量的线性函数作为目标函数. 明确问题的目标 假设一组决策变量 考虑目标函数 考虑约束条件 建立线性规划问题数学模型: 例3 某厂生产两种产品,产一吨甲产品用A资源 3吨、 B 甲 乙 上限 资源A 3 2