- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2.
构造数学模型
我们进行检验,先假设存在其他的决策单元组合的产出不 j0 j0 低于 而且投入尽可能的比 小,构造数学模型如下:
min
n j y rj y rj 0 (r 1,2, , s ) j 1 n s.t. j xij xij 0 (i 1,2, , m) j 1 n j 1, j 0( j 1,2, , n) j 1
上页 下页 返回
修正单纯形法简介
修正单纯形法要点:
寻求初始可行解,方法与单纯形法相同。 其迭代过程如下:
确定换入变量,方法与单纯形法相同。 确定换出变量,方法与单纯形法相同。 确定新的基可行解: 首先导出B-1 然后计算XB= B-1 b 迭代终止原则与单纯形法相同。
第二节
第一节
单纯形法的矩阵描述 及改进单纯形法介绍
单纯形法的矩阵描述
继续
改进单纯形法介绍
返回
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页
单纯形法的矩阵描述
设线性规划问题
不妨设基为
max z CX s.t AX b X 0
则
返回
B P1 P2 Pm A ( P1 P2 Pn ) ( B N ) X (XB XN ) C (CB CN )
令
XN 0
得当前的目标函数值为:
~ 1 z0 CBb CB B b
当前目标值
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
单纯形法的矩阵描述
检验数 ~ N C N CB N
~ ~ (Cm1 Cn ) (C1 Cm )( Pm1 Pn ) ~ n 1 Cm1 C B Pm1 当前检验数 ~ n Cn C B Pn
上页 下页 返回
3、计算步骤
例、解下列线性规划问题:
min 2 x1 3x2 - x3 s.t. x1 x2 x3 10 2 x1 x2 2 x3 8 0 x1 4 0 x2 4 1 x3 2
第三节
可分解的 大规模线性规划
学生讨论报告
返回
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
线性规划应用 ---数据包络分析法
数据包络分析法(Data Envelopment Analysis, 简称 DEA),是著名运筹学家 A.Charnes和W.W.Copper等 学者以“相对效率”概念为基础,以凸分析和线 性规划为工具,根据多指标投入和多指标产出对 相同类型的单位(部门)进行相对有效性或效益 评价的一种新的系统分析方法。它是处理多目标 决策问题的好方法。
矩阵单纯形法计算的描述
当基变量为 X B 时,新的单纯形表
基变量 非基变量
CB
X B B 1b cj zj
当前基解
XB I 0
XN Xs B 1 N B C N CB B 1 N CB B 1
当前检验数
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
变量有界的 大规模线性规划
返回
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
1、基本可行解概念的推广
考虑线性规划问题:
min s.t
z CX AX b l X u
A为m*n,秩为m
基本解X(0) :X(0)为AX=b的一个解,其中m个分量对应A 的列线性无关,其余n-m个分量取上界或下界值。 基本可行解X(0) :基本解X(0) 中m个基变量的值介于上下 界之间。
j
x j ( x1 j , x2 j ,, xmj )T 0, y j ( y1 j , y2 j ,, ysj )T 0,
j 1,2,, n
而且 xij 0, yrj 0, i 1,2,, m; r 1,2,, s 即每个决策单元都有m种类型的输入以及s种类型的输出
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
讨论最优性条件
设x是线性规划(LP)的一个基本可行解,若
对每个取下界值的非基变量,有
zj cj 0
对每个取上界值的非基变量,有
zj cj 0
则x是最优解。
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
非基变量
基变量
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
单纯形法的矩阵描述
约束方程组 XB AX b ( B N ) X N BX B NX N b
1
~ ~ X B B (b NX N ) b NX N ~ ~ 1 其中 b B b, N B1N
令
XN 0
~ 1 XB b B b
得当前的基解为: 当前基解
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
单纯形法的矩阵描述
目标函数 XB z (C B C N ) C B X B C N X N X N 1 1 C B B b (C N C B B N ) X N ~ ~ C B b (C N C B N ) X N
因此,我们假设该决策单元的第i项投入为
j 1 n j 1
n
j
xij (i 1,2,, m)
返回
产出为 j y rj (r 1,2,, s)
且 j
j 1
n
1( j 0)
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
上页 下页 返回
修正单纯形法简介
有关公式:
确定新的换入变量
1
j c j CB B Pj c j Pj
其中 C B B 1 单纯形乘子(行向量) ~ ~ 1 1 Pk B Pk , b B b i
确定新的换出变量
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
xij 为第j个决策单元对第i种类型输入的投入量;
y rj 为第j个决策单元对第r种类型输出的产出量。
这些都是已知的数据。
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页
DEA步骤
1.
假设
现在,我们是要最优化这些决策单元,那么我们假设 一个假想决策单元 j0 满足产出最大,同时投入最小。在此 基础上,我们来判断 j0 是否真的满足该条件。
推广基本可行解集与可行域凸集K的极点集等价
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
2、基本可行解的改进
设X(0)是一个基本可行解
xB ( 0) B 1b B 1 N l B 1 N u 1 N1 2 N2 (0) (0) X x N1 l N1 (0) uN2 xN 2 令c c B , cN1 , c N 2 ) (
特点:
1. 2.
具有一定的输入和输出 在将输入转换成输出的过程中,努力实现自身的决策 目标。
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
重要概念
决策单元的相对有效性
评价的依据是决策单元的“输入”和“ 输出”数据,根据输入和输出数据来评价决 策单元的优劣。 决策单元的相对有效性(即决策单元的 优劣)被称为DEA有效,它用数学规划模型 计算比较决策单元之间的相对效率,为评价 对象作出评价。
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
DEA步骤
1.
假设
设某个DMU的输入向量为 x ( x , x ,x )T , 1 2 m T。 输出向量为 y ( y , y , y ) 1 2 s 则n个 DMU ( 1 j n )对应的输入、输出向量分别为:
上页 下页 返回
修正单纯形法简介
修正单纯形法的优点:
能够从问题的原来参数(A,b,C),
计算出单纯形表中所有的数据,只要导 B 1即可。 出 单纯形表中的任一数字,只要作部分的 矩阵乘法即可获得。
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
上页 下页 返回
换入变量?
max max z j c j , max ( z j c j ) 0 xk jR jR2 1
换出变量?
若k R1 , 令xk lk k, k 非负,且要保证解的可行性, 即 k 必须满足下列条件: (1)保持xB不越下界; (2)保持xB不越上界; (3)保持xk 不越上界。
修正单纯形法简介
原因:
单纯形法的目的是要求问题的最优解, 而在迭代过程中,单纯形表中的某些列与 求最优解关系不大。因此,对单纯形法进 行修正。
思路:
~ ~ 每次迭代关键求出 B , Pk b , Pk , j ,i
1
需要换入的变量对应的列
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述
单 大 纯 规 形 对 模 法 偶 线 矩 问 性 题 阵 规 描 划 述