运筹学概述及线性规划的软件求解.pptVIP

“运筹帷幄之中，决策千里之外” ——运 筹 学 Operations Research (Operational Research) 本次课主要内容 一.运筹学的研究目的 二.运筹学的研究步骤 一.运筹学的研究目的 1. 明确决策目标,紧紧围绕这一目标,强调目标的优化和实现这一目标的方案的优化； 2 .依靠改变资源的应用方式或方法，合理有效运用资源； 3 .决策方案须有定量的依据，且可操作性强。此外，也要尽量反所涉及的某些非定量的依据考虑进去； 4 .表达研究成果的技术是运筹学的重要组成部分，对研究成果的表述要尽量通俗易懂。 二.运筹学的研究步骤 1. 提出问题； 2. 建立模型； 3. 求解模型； 4. 解的检验（灵敏度分析）； 5. 解的控制（改进和修正）； 6. 解的实施。 线性规划模型的一般形式 (标准形式略) max (min) z = c1 x1 + c2 x2 + … + cn xn a11 x1 + a12 x2 + … + a1n xn ≤(≥,=) b1 s.t. a21 x1 + a22 x2 + … + a2n xn ≤ (≥,=) b2 …… …… am1 x1 + am2 x2 + … + amn xn ≤ (≥,=) bm x1 ，x2 ，… ，xn ≥ 0 二. 线性规划的灵敏度分析 在线性规划标准模型 max z = c · x s.t. Ax = b x ≥0 中,我们一般假设 A, b , c 都是常数,但实际中可能多种原因引起它们的变化,现在的问题是:这些系数在什么范围内变化时,使线性规划问题的最优解不变? 这就是灵敏度分析要研究的问题。 线性规划是静态模型 参数发生变化，原问题的最优解还是不是最优？ 哪些参数容易发生变化？ C, b, A 每个参数发生多大的变化不会破坏最优解？ 灵敏度越小，解的稳定性越好 2. 1 价值系数 cj 的灵敏度分析 cj 变动可能由于市场价格的波动，或生产成本的变动 cj 的灵敏度分析是在保证最优解的基变量不变的情况下，分析cj 允许的变动范围?cj cj 的变化会引起检验数的变化，有两种情况 非基变量对应的价值系数变化，不影响其它检验数； 基变量对应的价值系数变化，影响所有非基变量检验数。 1、非基变量对应的价值系数的灵敏度分析 2.2 基变量对应的价值系数的灵敏度分析 由于基变量对应的价值系数在CB中出现，因此它会影响所有非基变量的检验数 只有一个基变量的 cj? 发生变化，变化量为? cj? 令 cj? 在CB中的第k行，研究非基变量xj 机会成本的变化 2.3 右端项 bi 的灵敏度分析 设 XB=B?1b 是最优解，则有XB=B?1b?0 b 的变化不会影响检验数 b 的变化量 ?b 可能导致原最优解变为非可行解 2.4 技术系数 aij 的灵敏度分析 技术系数aij变化的影响比较复杂 对应基变量的 aij ，且资源bi已全部用完 对应基变量的 aij ，但资源bi未用完 对应非基变量的 aij ，且资源bi全用完或未用完 1、对应基变量的 aij ，且资源bi已全部用完 ?aij=0 2、对应基变量的 aij ，但资源bi未用完 ?? ? ?aij? xn+i /xj 上述两个公式不充分，为什么？ B–1发生变化，从而引起非基变量检验数 cj– zj 的变化 3、对应非基变量的 aij 只影响对应非基变量xj的检验数 cj– zj 若? aij 0，不会破坏最优解 若? aij 0，必须保证 cj– zj ? 0 三. 线性规划问题的软件求解 model: sets: row/1..m/: b; !m表示数组维数（约束条件行数） arrange/1..n/: x,c; !n表示数组维数（决策变量个数） link (row, arrange): a; !系数矩阵形状 endsets Data: b=b1,b2,…,bm ! 约束条件右边的常数值; c=c1,c2,…,cn !目标函数系数的实际数值; a=a(1,1),a(1,2),…,a(1,n),…a(m,1),a(m,2),…,a(m,n); !系数矩阵； Enddata [OBJ]max=@sum(a