数学建模五步法与灵敏度分析

灵敏度分析 简介 研究与分析一个系统 或模型 的状态或输出变化对系统参数或周围条件 变化的敏感程度的方法 在最优化方法中经常利用灵敏度分析来研究原始数据 不准确或发生变化时最优解的稳定性 通过灵敏度分析还可以决定哪些参数对 系统或模型有较大的影响 因此 灵敏度分析几乎在所有的运筹学方法中以及 在对各种方案进行评价时都是很重要的 用途 主要用于模型检验和推广 简单来说就是改变模型原有的假设条件之后 所得到的结果会发生多大的变化 举例 建模五步法 一头猪重 200 磅 每天增重 5 磅 饲养每天需花费 45 美分 猪的市场价格为 每磅 65 美分 但每天下降 1 美分 求出售猪的最佳时间 建立数学模型的五个步骤 1 提出问题 2 选择建模方法 3 推到模型的数学表达式 4 求解模型 5 回答问题 第一步 提出问题 将问题用数学语言表达 例子中包含以下变量 猪的重量 w 磅 从现在 到出售猪期间经历的时间 t 天 t 天内饲养猪的花费 C 美元 猪的市场价 格 p 美元 磅 出售生猪所获得的收益 R 美元 我们最终要获得的净收益 P 美元 还有一些其他量 如猪的初始重量 200 磅 建议先写显而易见的部分 猪从 200 磅按每天 5 磅增加 w 磅 200 磅 5 磅 天 t 天 饲养每天花费 45 美分 C 美元 0 45 美元 天 t 天 价格 65 美分按每天 1 美分下降 p 美元 磅 0 65 美元 磅 0 01 美元 磅 t 天 生猪收益 R 美元 p 美元 磅 w 磅 净利润 P 美元 R 美元 C 美元 用数学语言总结和表达如下 参数设定 t 时间 天 w 猪的重量 磅 p 猪的价格 美元 磅 C 饲养 t 天的花费 美元 R 出售猪的收益 美元 P 净收益 美元 假设 w 200 5t C 0 45t p 0 65 0 01t R p w P R C t 0 目标 求 P 的最大值 第二步 选择建模方法 本例采用单变量最优化问题或极大 极小化问题 第三步 推导模型的数学表达式子 P R C 1 R p w 2 C 0 45t 3 得到 R p w 0 45t p 0 65 0 01t 4 w 200 5t 5 得到 P 0 65 0 01t 200 5t 0 45t 令 y P 是需最大化的目标变量 x t 是自变量 现在我们将问题转化为集合 S x x 0 上求函数的最大值 y f x 0 65 0 01x 200 5x 0 45x 1 1 第四步 求解模型 用第二步中确定的数学方法解出步骤三 例子中 要求 1 1 式中定义的 y f x 在区间 x 0 上求最大值 下图给出了 1 1 的图像和导数 应用几 何画板绘制 在 x 8 为全局极大值点 此时 f 8 133 20 因此 8 133 20 为 f 在整个实轴上的全局极大值点 同时也是区间 x 0 上的最 大值点 第五步 回答问题 根据第四步 8 天后出售生猪的净收益最大 可以获得净收益 133 20 美元 只 要第一步中的假设成立 这一结果正确 数学建模五步方法总结 第一步 提出问题 1 列出问题中涉及的变量 包括适当的单位 2 注意不要混淆变量和常量 3 列出你对变量所做的全部假设 包括等式和不等式 4 检查单位从而保证你的假设有意义 5 用准确的数学术语给出问题的目标 第二步 选择建模方法 1 选择解决问题的一个一般的求解方法 2 一般地 这一步的成功需要经验 技巧和熟悉相关文献 第三步 推导模型的数学表达式 1 将第一步中得到的问题重新表达成第二步选定的建模方法所需要的形式 2 将第一步中的一些变量名改成与第二步所用的记号一致 3 记下任何补充假设 这些假设是为了使第一步中描述的问题与第二步中选 定的数学结构相适应而做出的 第四步 求解模型 1 将第二步中所选用的一般求解过程应用于第三步得到表达式的特定问题 2 注意你的数学推导 检查是否有错误 你的答案是否有意义 3 采用适当的技术 计算机代数系统 图形工具 数值计算的软件等 都能 扩大你能解决问题的范围 并能减少计算错误 第五步 回答问题 1 用非技术性的语言将第四步的结果重新表述 2 避免数学符号和术语 3 能理解出处提出的问题的人就应该能理解你给出的答案 灵敏度分析 数据是由测量 观察有时甚至完全猜测得到的 因此 我们要考虑数据不 准确的可能性 上例中 生猪现在的重量 现在的价格 每天饲养花费都很容易测量 而 且有相当大的确定性 但是猪的生长率则不那么确定 而价格的下降率则确定 性更低 记 r 为价格的下降率 现在假设 r 的实际值不同 对几个不同的 r 值 重复前面的求解过程 我们会对问题的解关于 r 的敏感程度有所了解 下表给 出了几个不同 r 值求出的计算结果 根据表格绘制图形 我们可以看到售猪的 最优时间对参数 r 很敏感 r 美元 天 x 天 0 00815 0 0 00911 1 0 0108 0 0 0115 5 0 0123 3 对灵敏度的更系统的分析是将 r 视为未知参数 按前面的步骤求解 写出 p 0 65 rt 得到 y f x 0 65 rx 200 5x 0 45x 使得导数为 0 得 到 x 7 500r 25r 当 x 0 时 只要 0 r 0 时 得 到 g 3 769 我们将灵敏度数据用相对改变量表示 例如 r 下降 10 导致了 x 增加了 39 而 g 下降了 10 导致了 x 下降了 34 如果 x 的改变量 x 则 x x 表示相对改变量 如果 r 改变了 r 导致 了 x 有 x 的改变量 则相对改变量的比值为 x x r r 令 r 0 我们有 x x r r dx dr r x 我们称这个极限值 为 x 对 r 的灵敏度 即为 S x r 在售猪问题中 r 0 01 和 x 8 得到 dx dr 7 25r2 2800 因此 S x r dx dr r x 2800 0 01 8 7 2 即若 r 增加 2 则 x 下降 7 由于 dx dg 245 2g2 4 9 我们有 S x g dx dg g x 4 9 5 8 3 0625 于是猪的生长率增加 1 会导致大约等待 3 的时间再将猪售出 灵敏度分析的成功应用要有较好的判断力 通常即不可能对模型中的每个 参数都计算灵敏度分析 也没有特别的要求 我们需要选择那些有较大不确定 性的参数进行灵敏度分析 对灵敏度系数的解释还要依赖与参数的不确定程度 主要问题是数据的不确定程度影响答案的置信度 在这个问题中 我们通常认 为猪的生长率 g 比价格下降率 r 更可靠 如果我们观察了猪或者其他类似动物 在过去的生长情况 则 g 有 25 的误差会是很不寻常的 但对 r 的估计有 25 的 误差则不足为奇 数学模型的稳健性 一个数学模型称为稳健的 是指即使这个模型不完全精确 由其导出的结 果也是正确的 在实际问题中 我们不会有绝对准确的信息 即使能够建立一 个完美的精确模型 我们也可能采取较为简单和易于处理的方法 出于数学处 理的方便和简化的目的 常常要做一些假设 建模者有责任要考察这些假设是 否太特殊 以致使模型的结果无效 上例中我们主要是假设猪的重量和每磅的价格都是时间线性函数 假设一 年后 猪的重量为 200 5 365 2025 磅 卖出收益为 0 65 0 01 365 3 美元 磅 一个更为实际的模型应该考虑到这些函数的非线性性 又考虑到随着时间的推 移不确定性的增加 考察售猪问题中的线性假设 基本方程为 P pw 0 45t 如果模型初始数据 和假设没有与实际相差太远 则售猪的最佳时间应该有令 P 求导为 0 得到 计 算后有 p w pw 0 45 得到只要猪价比饲养的费用增长快 就应暂时不卖出 其中 p w 为价格下降带来的损失 pw 为猪增重而增加的价值 考虑更一般的 模型的情况 猪的未来增长和价格的未来变化并