1、第 1 页 共 30 页第一部分 线性规划问题的求解一、两个变量的线性规划问题的图解法：概念准备：定义：满足所有约束条件的解为可行解；可行解的全体称为可行（解）域。定义：达到目标的可行解为最优解。图解法：图解法采用直角坐标求解：x 1横轴；x 2竖轴。1、将约束条件（取等号）用直线绘出；2、确定可行解域；3、绘出目标函数的图形（等值线） ，确定它向最优解的移动方向；注：求极大值沿价值系数向量的正向移动；求极小值沿价值系数向量的反向移动。4、确定最优解及目标函数值。参考例题：（只要求下面这些有唯一最优解的类型）例 1：某厂生产甲、乙两种产品，这两种产品均需在 A、B 、C 三种不同的设备上加工，
2、每种产品在不同设备上加工所需的工时不同，这些产品销售后所能获得利润以及这三种加工设备因各种条件限制所能使用的有效加工总时数如下表所示：A B C 利润（万元）甲乙3 5 99 5 37030有效总工时 540 450 720 问：该厂应如何组织生产，即生产多少甲、乙产品使得该厂的总利润为最大？设备消耗产品第 2 页 共 30 页（此题也可用“单纯 形法”或化“对 偶问题”用大 M 法求解）第 3 页 共 30 页解：设 x1、x 2 为生产甲、乙产品的数量。max z = 70x1+30x2s.t. 0723945502112xx，可行解域为 oabcd0，最优解为 b 点。由方程组解出 x1
3、=75，x 2=15720394551xX *= =（75，15） T2xmax z =Z *= 7075+3015=5700、第 4 页 共 30 页例 2：用图解法求解max z = 6x1+4x2s.t. 078022112xx，解：可行解域为 oabcd0，最优解为 b 点。由方程组解出 x1=2，x 2=6810212xX *= =（2，6） T1xmax z = 62+46=36、第 5 页 共 30 页例 3：用图解法求解min z =3x 1+x2s.t. 08215341121xx，解：可行解域为 bcdefb，最优解为 b 点。由方程组 解出 x1=4，x 2=125241
4、x 54X *= =（4， ） T2xmin z =34+ =1151、第 6 页 共 30 页二、标准型线性规划问题的单纯形解法：一般思路：1、用简单易行的方法获得初始基本可行解；2、对上述解进行检验，检验其是否为最优解，若是，停止迭代，否则转入 3；3、根据 L 规则确定改进解的方向；4、根据可能改进的方向进行迭代得到新的解；5、根据检验规则对新解进行检验，若是最优解，则停止迭代，否则转入 3，直至最优解。具体做法（可化归标准型的情况）：设已知max z = c1x1+ c2x2+ cnxns.t. njxbxaaxbxaaj mmn， .21021 22212 11对第 i 个方程加入松
5、弛变量 xn+i，i =1，2，m ，得到 njx bxxaabxaaj mnmmnn， .21021 22212 11列表计算，格式、算法如下：第 7 页 共 30 页c1 c2 cn+mCB XB bx1 x2 xn+mLcn+1 xn+1 b1 a11 a12 a1 n+mc n+2 xn+2 b2 a21 a22 a2 n+m. cn+m xn+m bn am1 am2 am n+mz1 z2 zn+m1 2 n+m注： zj =cn+1 a1j+ cn+2 a2j + cn+m amj= ， （j=1，2，n+m）miijnac1j =cjzj ,当 j 0 时，当前解最优。注：由
6、max j确定所对应的行的变量为“入基变量”；由 L= 确定所对应的行的变量为“出基变量”，行、列0minikab交叉处为主元素，迭代时要求将主元素变为 1，此列其余元素变为 0。例 1：用单纯形法求解 （本题即是本资料 P2“图 解 法 ”例 1 的单纯形解法 ； 也可化 “对偶 问 题 ”求解）max z =70x1+30x2s.t. 0723945502112xx，解：加入松弛变量 x3，x 4，x 5，得到等效的标准模型：max z =70x1+30x2+0 x3+0 x4+0 x5第 8 页 共 30 页s.t. 5,.21,0720394551231jxxxj列表计算如下：70 3
7、0 0 0 0CB XB b x1 x2 x3 x4 x5 L0 x3 540 3 9 1 0 0 540/3 =1800 x4 450 5 5 0 1 0 450/5 =900 x5 720 （9） 3 0 0 1 720/9 =800 0 0 0 070 30 0 0 00 x3 300 0 8 1 0 - 1/3 300/8 =37.50 x4 50 0 （10/3） 0 1 - 5/9 50/10/3 =1570 x1 80 1 1/3 0 0 1/9 80/1/3 =24070 70/3 0 0 70/90 20/3 0 0 70/90 x3 180 0 0 1 12/5 130 x
8、2 15 0 1 0 3/10 - 1/670 x1 75 1 0 0 - 1/10 1/670 30 0 2 20/357000 0 0 -2 20/3X *=（75，15，180 ，0，0） Tmax z =7075+3015=5700第 9 页 共 30 页例 2：用单纯形法求解max z =7x1+12x2s.t. 031325460921xx，解：加入松弛变量 x3，x 4，x 5，得到等效的标准模型：max z =7x1+12x2+0 x3+0 x4+0 x5s.t. 5,.21,0303546914231jxxxj列表计算如下：第 10 页 共 30 页7 12 0 0 0CB
9、XB b x1 x2 x3 x4 x5 L0 x3 360 9 4 1 0 0 360/4 =900 x4 200 4 5 0 1 0 200/5 =400 x5 300 3 （10） 0 0 1 300/10 =300 0 0 0 07 12 0 0 00 x3 240 78/10 0 1 0 - 2/5 240/78/10 =2400/780 x4 50 （5/2） 0 0 1 - 1/2 50/5/2 =2012 x2 30 3/10 1 0 0 1/10 30/3/10 =10018/5 12 0 0 6/517/5 0 0 0 6/50 x3 84 0 0 1 78/25 29/25
10、7 x1 20 1 0 0 2/5 - 1/512 x2 24 0 1 0 3/25 4/287 12 0 34/25 11/354280 0 0 34/25 11/35X *=（20，24，84 ，0，0） Tmax z =720+1224=428三、非标准型线性规划问题的解法：1、一般地，对于约束条件组：若为“” ，则加松弛变量，使方程成为“” ；若为“” ，则减松弛变量，使方程成为“” 。我们在前面标准型中是规定目标函数求极大值。如果在实际问题中遇到的是求极小值，则为非标准型。可作如下处理：第 11 页 共 30 页由目标函数 min z= 变成等价的目标函数 max（z ）=njjxc
11、1 njjxc1)(令z=z /，min z=max z /2、等式约束大 M 法：通过加人工变量的方法，构造人造基，从而产生初始可行基。人工变量的价值系数为M，M 是很大的正数，从原理上理解又称为“惩罚系数” 。 （课本P29）类型一 ：目标函数仍为 max z，约束条件组与。例 1：max z =3x1+5x2s.t. 01823411xx，解：加入松弛变量 x3，x 4，得到等效的标准模型：max z =3x1+5x2s.t. 4,321,08324131jxxj其中第三个约束条件虽然是等式，但因无初始解，所以增加一个人工变量x5，得到： max z =3x1+5x2Mx 5第 12 页
12、 共 30 页s.t. 5,.21,08324131jxxxj单纯形表求解过程如下：3 5 0 0 MCB XB bx1 x2 x3 x4 x5L0 x3 4 （1） 0 1 0 0 4/1 =40 x4 12 0 2 0 1 0 M x5 18 3 2 0 0 1 18/3 =6 3M 2M 0 0 M3M3 5 2M 0 0 03 x1 4 1 0 1 0 0 0 x4 12 0 2 0 1 0 12/2 =6 M x5 6 0 （2） 3 0 1 6/2 =33 2M 3 3M 0 M0 5 3 3M 0 03 x1 4 1 0 1 0 0 4/1 =40 x4 6 0 0 （3） 1
13、1 6/3 =25 x2 3 0 1 3/2 0 1/2 3/(2/3) =9/23 5 9/2 0 5/20 0 9/2 0 M5/2x1 2 1 0 0 1/3 1/330 x3 2 0 0 1 1/3 1/3第 13 页 共 30 页5 x2 6 0 1 0 1/2 03 5 0 3/2 1360 0 0 3/2 M1X *=（2，6，2，0 ） Tmax z =32+56=36类型二 ：目标函数 min z，约束条件组与。例 2：用单纯形法求解min z =4x1+3x2s.t. 0123642211x，解：减去松弛变量 x3，x 4，并化为等效的标准模型：max z/ =4x 13x
14、 2s.t. 4,321,0362132jxxj增加人工变量 x5、x 6，得到：max z/ =4x 13x 2Mx 5Mx 6s.t 6,.21,012342641 532jxxxj第 14 页 共 30 页单纯形表求解过程如下：第 15 页 共 30 页4 0 0 M MCB XB bx1 x2 x3 x4 x5 x6 LM x5 16 2 （4） 1 0 1 0 16/4=4M x6 12 3 2 0 1 0 1 12/2=65M 6 M M M M M5M 4 6M3 M M 0 03 x2 4 1/2 1 1/4 0 1/4 0 4/1/2=8M x6 4 （2） 0 1/2 1
15、1/2 1 4/2=22M 3/2 3 3/4M/2 M M/2 3/4 M2M5 /2 0 M/23/4 M 3/43M/2 03 x2 3 0 1 3/8 1/4 3/8 1/44 x1 2 1 0 1/4 1/2 1/4 1/24 3 1/8 5/4 1/8 5/4170 0 1/8 5/4 M 1/8M5/4X *=（2，3，0，0 ） Tmin z =max z / =（17）=17第 16 页 共 30 页四、对偶问题的解法：什么是对偶问题？1、在资源一定的条件下，作出最大的贡献；2、完成给定的工作，所消耗的资源最少。引例 （与本资料 P2 例 1 “图 解 法 ”、 P7 例 1
16、 “单纯 形 法 ”同） ：某工厂生产甲、乙两种产品，这些产品均需在 A、B、C 三种不同的设备上加工，每种产品在不同设备上加工时需要不同的工时，这些产品售后所能获得的利润值以及这三种加工设备因各种条件下所能使用的有效总工时数如下表：A B C 利润（万元）甲乙3 5 99 5 37030有效总工时 540 450 720 问：该厂应如何组织生产，即生产多少甲、乙产品使得该厂的总利润为最大？解：原问题设 x1、x 2 为生产甲、乙产品的数量。max z = 70x1+30x2s.t. 0723945502112xx，设备消耗产品第 17 页 共 30 页将这个原问题化为它的对偶问题设 y1、y
17、 2、y 2 分别为设备 A、B、C 单位工时数的加工费。min w = 540y1+450y2+720y3s.t. 3210059793132， iyyi用大 M 法，先化为等效的标准模型：max w/ =540y 1450y 2720y 3s.t. 5,.21,030359793142jyyj增加人工变量 y6、y 7，得到：max z/ =540y 1450y 2720y 3My 6My 7s.t 5,.21,0303597931 642jy yyj大 M 法单纯形表求解过程如下：第 18 页 共 30 页540 450 720 0 0 M MCB XB by1 y2 y3 y4 y5
18、y6 y7 LM y6 70 3 5 9 1 0 1 0 70/3M y7 30 （9） 5 3 0 1 0 1 30/9=10/3 12M 10M 12M M M M M12M540 10M450 12M720 M M 0 0M y6 60 0 10/3 （8） 1 1/3 1 1/3 60/8=2.5540 y1 10/3 1 5/9 1/3 0 1/9 0 1/9 10/3/1/3=10-300+10/3M -8M180 M M/3+60 M M/3600 -150+10/3M 8M-540 M M/3 60 0 M/3+60720 y3 15/2 0 5/12 1 1/8 1/24 1
19、/8 1/24 15/2/5/12=18540 y1 5/6 1 （5/12） 0 1/24 1/8 1/24 1/8 5/6/5/12=2540 572 720 135/2 475/12 135/2 75/20 125 0 135/2 475/12 135/2M 75/2My3 20/3 1 0 1 1/6 1/6 1/6 1/6720450 y2 2 12/5 1 0 1/10 3/10 1/10 3/10360 450 720 75 15 75 155700180 0 0 75 15 75M 15M该对偶问题的最优解是 y*=（0，2， ，0，0） T3最优目标函数值 min w =（5
20、700）=5700第 19 页 共 30 页五、运输规划问题：运输规划问题的特殊解法“表上作业法”解题步骤：1、找出初始调运方案。即在(mn)产销平衡表上给出 m+n-1 个数字格。 （最小元素法）2、 （对空格）求检验数。判别是否达到最优解。如已是最优解，则停止计算，否则转到下一步。 （闭回路法）3、对方案进行改善，找出新的调运方案。 （根据检验结果选择入基变量，用表上闭回路法调整即迭代计算，得新的基本可行解）4、重复 2、3，再检验、再迭代，直到求得最优调运方案。类型一 ：供求平衡的运输规划问题（又称“供需平衡” 、 “产销平衡” ）引例：某钢铁公司有三个铁矿和四个炼铁厂，铁矿的年产矿石量
21、分别为100 万吨、80 万吨和 50 万吨，炼铁厂年需矿石量分别为 50 万吨、70 万吨、80万吨和 30 万吨，这三个铁矿与四个炼铁厂的距离如下：B1 B2 B3 B4A1A2A315 20 3 3070 8 14 2012 3 20 15问：该公司应如何组织运输，既满足各炼铁厂需要，又使总的运输费用为最小（按吨 .公里计）？解：用“表上作业法”求解。先用最低费用法（最小元素法）求此问题的初始基础可行解：炼铁距离铁矿厂第 20 页 共 30 页B1 B2 B3 B4产量Si15 20 67 3 30 65A120 80 10070 8 14 44 20 A230 20 308012 53
22、 3 20 33 25 10A3 50 50销量dj 50 70 80 30230230初始方案：Z=1520+380+7030+820+2030+350=3550（吨 .公里）对的初始可行解进行迭代（表上闭回路法） ，求最优解：销地费用产 地2080B1B3A1 B2203030B1B3A2 B250A3第 21 页 共 30 页B1 B2 B3 B4产量Si15 20 14 3 30 12A120 80 10070 53 8 14 9 20 A2 50 308012 3 20 20 25 10A330 20 50销量dj 50 70 80 30230230用表上闭回路法调整后，从上表可看出
23、，所有检验数 0，已得最优解。最优方案：Z=1520+380+850+2030+1230+320=1960（吨 .公里）解法分析： 如何求检验数并由此确定入基变量？有数字的空格称为 “基格 ”、打 的空格称为 “空格 ”，标号为偶数的顶点称为偶点、标号为奇数的顶点称为奇点，出发点算 0 故为偶点。找出所有空格的闭回路后计算它们的检验数 ，必须 0，才偶 点奇 点 ijijij cij得到最优解。否则，应选所有 中正的最大者对应的变量 xj 为入基变量进行迭代（调整） 。 检验后调整运输方案的办法是：在空格的闭回路中所有的偶点均加上奇点中的最小运量，所有的奇点均减去奇点中的最小运量。 重复以上两
24、步，再检验、再调整，直到求得最优运输方案。类型二 ：供求不平衡的运输规划问题销地费用产 地2080B1B3A15030B2B4A23020B1B2A3第 22 页 共 30 页若 ，则是供大于求（供过于求）问题，可设一虚销地 Bn+1，njjmids11令 ci,n+1=0，d n+1= ，转化为产销平衡问题。若 ，则是供njjmids11 njjmids11小于求（供不应求）问题，可设一虚产地 Am+1，令 cm+1,j=0，s m+1=，转化为产销平衡问题。minjjsd11（，2，m；，2，n）六、工作指派问题：工作指派问题的数学模型 假定有 n 项工作需要完成，恰好有 n 个人每人可去
25、完成其中一项工作，效果要好。工作指派问题的特殊解法 “匈牙利法” （考！）解题步骤：1、使系数矩阵（效率矩阵）各行、各列出现零元素。作法：行约简系数矩阵各行元素减去所在行的最小元素，列约简再从所得矩阵的各列减去所在列最小元素。2、试求最优解。如能找出 n 个位于不同行不同列的零元素，令对应的 xij= 1，其余 xij = 0，得最优解，结束；否则下一步。作法：由独立 0 元素的行（列）开始，独立 0 元素处画（ ）标记 ，在有（ ）的行列中划去（也可打*）其它 0 元素；再在剩余的 0 元素中重复此做法，直至不能标记（ ）为止。3、作能覆盖所有 0 元素的最少数直线集合。作法： 对没有（ ）
26、的行打号；对已打号的行中所有 0 元素的所在列打号；再对打有号的列中 0 元素的所在行打号；重复、直到得不出新的打号的行(列)为止； 对没有打号的行画一横线，对打号的列画一纵线，这就得到覆盖所有 0 元素的最少直线数。未被直线覆盖的最小元素为 cij,在未被直线覆盖处减去 cij,在直线交叉处加上 cij。4、重复 2、3，直到求得最优解。第 23 页 共 30 页类型一 ：求极小值的匈牙利法：（重点掌握这种基本问题）例 1：有甲、乙、丙、丁四个人，要派去完成 A、B、C、D 四项工作，他们完成的工时如下表：A B C D甲乙丙丁6 12 13 410 3 12 147 14 13 168 8
27、 12 10试问：应如何分配任务可使总工时为最少？解：用“匈牙利法”求解。已知条件可用系数矩阵（效率矩阵）表示为：（c ij）= 10286347106 2409671822097152 2)0(97)(15)(82*使总工时为最少的分配任务方案为：甲D，乙B，丙A，丁C此时总工时数 W=4+3+7+12=26例 2：求效率矩阵列约简任务工时人行约简标号甲乙丙丁A B C D第 24 页 共 30 页的最优解。54326578910解： 54326578910 32103所画（）0 元素少于 n，未21003 210)(13)(*得到最优解，需要继续变换矩阵（求能覆盖所有 0 元素的最少数直线
28、集合）：2100)(13)(2*未被直线覆盖的最小元素为 cij=1,在未被直线覆盖处减去 1,在直线交叉处加上 1。1024 10022)()(4*）（ ）（得最优解： 010类型二 ：求极大值的匈牙利法：min z=max （z ）标号列约简行约简标号第 25 页 共 30 页（c ij）（M c ij）=（b ij） ， （c ij）中最大的元素为 Mmax z= =jijixj iji xM)(= j iji c)(jijic第一部分到 此结束第 26 页 共 30 页第二部分 动态规划只要求掌握动态规划的最短路问题用“图上标号法”解决：具体解题步骤请参看教材 P103（这是本套资料少
29、见的与教材完全相同的算法类型之一，务必看书掌握）学员们只有完全理解了这种作法（思路：逆向追踪）才有可能做题，考试时数字无论如何变化都能作出正确求解！第二 部分到此 结束第 27 页 共 30 页第三部分 网络分析一、求最小生成树（最小支撑树、最小树）问题：破圈法任取一个圈，从圈中去掉一条权最大的边（如果有两条或两条以上的边都是权最大的边，则任意去掉其中一条） 。在余下的图中，重复这个步骤，直到得到一个不含圈的图为止，这时的图便是最小树。参考例题：例：求下图的最小生成树：67941 5 10v2v1v3 v5v4v6328第 28 页 共 30 页解：用“破圈法”求得最小生成树为：已得最小树，此
30、时权 w=1+2+4+5+9=21 为最小。二、最短路问题：（有向图）TP 标号法（狄克斯托算法）具体解题步骤请参看教材 P125（这是本套资料少见的与教材完全相同的算法类型之一，务必看书掌握）学员们只有完全理解了这种作法（思路：顺向追踪）才有可能做题，考试时数字无论如何变化都能作出正确求解！参考例题：例：教材 P124 图 48 的例子（略）三、网络最大流问题：寻求网络最大流的标号法（福特富克尔逊算法）具体解题步骤请参看教材 P130。 （教材用有序数对（ fij，c ij）表示（流量，容量） ，但我们解题算法格式按南邮要求，用大多数运筹学书籍中的标准格式，即用有序数对（c ij，f ij）
31、表示（容量，流量） 。解法本质是相同的。）学员们只有完全理解了这种作法（思路：标号检查、迭代调整）才有可能做题，考试时数字无论如何变化都能作出正确求解！941 5v2v1v3 v5v4v62第 29 页 共 30 页参考例题：运筹学参考综合习题 （我站搜集信息自编，非南邮综合练习题，仅供参考。另挂于网上）第 三部分到 此结束第 30 页 共 30 页第四部分 存储论（简介）随机存储模型参考例题：例：一食品店要决定每天牛奶的进货量，该店根据过去的销售经验，知道需求量的概率分布如下：需求 x（箱） 25 26 27 28p（x） 0.1 0.3 0.5 0.1若进货每箱 80 元，售价 100 元
32、，又若当天不能售出因牛奶变质而全部损失，试确定每天进货量。解法一：已知 C=80， p=100，g=0，需求 x（箱） 25 26 27 28概率 p（x） 0.1 0.3 0.5 0.1累计概率 0.1 0.4 0.9 1根据单周期随机型存储模型 （“报童模型” ）之离散型随机存储模型公式，可得0.2018)(*sgpCQxp即可以确定进货 26 箱，获利的期望值最大。解法二：k=100 80=20，h=80.= =0.2kh102Px=25=0.10.2Px=25+Px=26=0.4可以确定进货 26 箱，获利的期望值最大。 。注：该问题即 “报童问题 ”（关于报童进报纸问题的数学模型） ，即是相当于将上题的进 “牛奶 ”改为进“报纸 ”等等，解法思路是完全一致的，请注意！ 第四部 分到 此结束
