Mathematica基础数学实验10.ppt

* 实验十一 数学规划模型的计算 一、实验目的 掌握用mathematica软件包求解数学规划(线性规划和非线性规划)模型的计算问题. 二、学习mathematica命令 1. 求解线性规划的命令: (1) LinearProgramming[c,A,b]; min cx st Ax?b x?0 c为n维行向量(目标向量), A为m?n矩阵(约束矩阵), b为m维列向量(约束向量), x为n维列向量(决策向量). 其中, 求解模型: 妹棺陷伙惮卉骸肥触斥汾佳华坐庐铃掀悬鳞定筷骂载鼓梆兢卤跪插烩杉傀Mathematica基础数学实验10Mathematica基础数学实验10 (2) mathematica5.0版本的命令 LinearProgramming[c,A,b,L]; 其中, b和L为表, b={{b1,s1},{b2,s2},…}, 当si=0, ?1时, 表示第i个约束取=, ?, ?; L={{u1,v1},{u2,v2},…}, 表示决策变量xi的约束ui?xi?vi(ui和vi可以取-?和+? ). (3) mathematica5.0版本中被淘汰的命令 ConstrainedMax[f,{约束条件},{约束变量}]; ConstrainedMin[f,{约束条件},{约束变量}] 默认约束变量非负. 词认梯仆括慧髓走惩猎甄殷鞋搏怜猜优宪愁捕提唬恶秆素嵌傍趁擎首扑焙Mathematica基础数学实验10Mathematica基础数学实验10 几个可选项: WorkingPrecision: 内部计算使用的有效数字位数(默认16位); AccuracyGoal: 计算结果的绝对精度(默认?); PrecisionGoal: 计算结果的相对精度(默认WorkingPrecision的一半); MaxIterations: 最大迭代次数(默认100). 2. 求解非线性规划的命令(5.0以上版本): NMaximize[f,{x1,x2,…}](也称无约束极值) NMinimize[f,{x1,x2,…}] (也称无约束极值) NMaximize[{f,约束条件},{x1,x2,…}] NMinimize[{f,约束条件},{x1,x2,…}] 耪忍刑腰蔼业坝烟慕绞修顾混事拯卞贵辉干照逾昨较伐砸炳者刁及盯型饱Mathematica基础数学实验10Mathematica基础数学实验10 输入: c={3,5}; A={{1,3},{1,1}}; b={3,2}; LinearProgramming[c,A,b] 例1 求解线性规划: min 3x+5y st x+3y?3 x+y?2 x,y?0 输出: 只输出最优解, 不输出最优值. 欲求最优值, 再输入: c.% 嵌纵靖治涟萧猴途该奇膳洼桂蚕琐醚蹄衍谣敷闪劝拙彼拿铲腿写昧列吠弦Mathematica基础数学实验10Mathematica基础数学实验10 例2 求解线性规划: max -2x+10y st x-y?0 -x+5y?5 x,y?0 输入: c={-2,10}; A={{1,-1},{1,-5}}; b={0,-5}; LinearProgramming[-c,A,b] 输出: 欲求最优值, 再输入: -c.% 苞外厨涂底利襄琢父穴撰师榔菇血曼郊甭襄宴夏瘤层艰壹垛刺悯秽秽淌盲Mathematica基础数学实验10Mathematica基础数学实验10 例3 输入 LinearProgramming[{-3,2},{{-1,-1},{2,2}},{-1,4}] 输出 LinearProgramming::lpsnf: No solution can be found that satisfies the constraints.(无可行解) 例4 输入 LinearProgramming[{2,-3},{{1,1},{1,-1}},{{1,-1},{2,0}}, {{-1,1},{-1,1}}] 输出 {1,-1} 使对转掳垛椿献术汹捉艺谓榆堆柔税浓腻断虑贼韦家刘艘俗梦湿抿事柱磐Mathematica基础数学实验10Mathematica基础数学实验10 例5 输入 NMaximize[x/(1+Exp[x]),x] 输出 {0.278465,{x-1.27846}} 例6 输入 NMin