1、第七章第七章 线性规划模型的建立与应用线性规划模型的建立与应用学习目的与要求 线性规划是经济领域广泛应用的一种经济分析方法。讲授本章目的是使同学掌握线性规划分析法的基本原理，掌握图解法和单纯形解法的程序及运算，并借助电化教学，能够初步应用线性规划法解决最低成本的农业生产资源最优配合方式和最大收益的生产结构问题。第七章第七章 线性规划模型的建立与应用线性规划模型的建立与应用一、线性规划的概念二、线性规划三要素三、技术经济研究中运用线性规划方法的特点及局限性四、线性规划模型的基本结构五、线性规划模型的一般形式六、线性规划模型的基本假设 第一节 线性规划模型的基本原理 线性规划是指如何最有效或最佳地
2、谋划经济活动。它所研究的问题有两类：一类是指一定资源的条件下，达到最高产量、最高产值、最大利润；一类是，任务量一定，如何统筹安排，以最小的消耗取完成这项任务。如最低成本问题、最小投资、最短时间、最短距离等问题。前者是求极大值问题，后者是求极小值问题。总之，线性规划是一定限制条件下，求目标函数极值的问题。第一节 线性规划模型的基本原理 一、线性规划的概念 经济大词典定义线性规划：一种具有确定目标，而实现目标的手段又有一定限制，且目标和手段之间的函数关系是线性的条件下，从所有可供选择的方案中求解出最优方案的数学方法。第一节 线性规划模型的基本原理 一、线性规划的概念二、线性规划三要素1.目标函数最
3、优化单一目标 多重目标问题如何处理？2.实现目标的多种方法 若实现目标只有一种方法不存在规划问题。3.生产条件的约束资源是有限的 资源无限不存在规划问题。第一节 线性规划模型的基本原理 三、技术经济研究中运用线性规划方法的特点及局限性 第一节 线性规划模型的基本原理 特点：1.可以使研究对象具体化、数量化。可以对所研究的技术经济问题做出明确的结论；2.线性3.允许出现生产要素的剩余量4.有一套完整的运算程序三、技术经济研究中运用线性规划方法的特点及局限性 第一节 线性规划模型的基本原理 局限性：1.线性规划它是以价格不变和技术不变为前提条件的，不能处理涉及到时间因素的问题。因此，线性规划只能以
4、短期计划为基础。2.在生产活动中，投入产出的关系不完全是线性关系，由于在一定的技术条件下，报酬递减规律起作用，所以要满足线性假定是不可能的。在线性规划解题中，常常把投入产出的非线性关系转化为线性关系来处理，以满足线性的假定性，客观上产生误差。3.线性规划本身只是一组方程式，并不提供经济概念，它不能代替人们对现实经济问题的判断。四、线性规划模型的基本结构1.决策变量 未知数。它是通过模型计算来确定的决策因素。又分为实际变量求解的变量和计算变量，计算变量又分松弛变量（上限）和人工变量（下限）。2.目标函数经济目标的数学表达式。目标函数是求变量的线性函数的极大值和极小值这样一个极值问题。3.约束条件
5、实现经济目标的制约因素。它包括：生产资源的限制（客观约束条件）、生产数量、质量要求的限制（主观约束条件）、特定技术要求和非负限制。第一节 线性规划模型的基本原理 四、线性规划模型的基本结构 Min Z=10 x1+20 x2 s.t.x1+x210 3x1+x215 x1+6x215 x10,x20 约束条件目标函数第一节 线性规划模型的基本原理 五、线性规划模型的一般形式Max Z=c1x1+c2x2+c3x3+cnxn a11x1+a12x2+a1nxn b1 (1)a21x1+a22x2+a2nxn b2 (2)am1x1+am2x2+amnxn bm (m)x1，x2，xn0第一节 线
6、性规划模型的基本原理 极大值模型njxbxaxcxcxcZjinjjijnn,3,2,1,0max12211其简缩形式为其简缩形式为 第一节 线性规划模型的基本原理 极大值模型五、线性规划模型的一般形式Min Z=c1x1+c2x2+c3x3+cnxn a11x1+a12x2+a1nxn b1 (1)a21x1+a22x2+a2nxn b2 (2)am1x1+am2x2+amnxn bm (m)x1，x2，xn0第一节 线性规划模型的基本原理 极小值模型njxbxaxcxcxcZjinjjijnn,3,2,1,0min12211其简缩形式为其简缩形式为 第一节 线性规划模型的基本原理 极小值模
7、型其简缩形式为其简缩形式为 第一节 线性规划模型的基本原理 极大值模型可用向量表示：01jnjjjxbxPCXzMaxnxxxX21mjjjjaaaP21bmbbb21 C=(c1,c2,cn)六、线性规划模型的基本假设1.线性线性 目标函数和约束条件目标函数和约束条件2.可分性可分性 活动对资源的可分性活动对资源的可分性3.可加性可加性 活动所耗资源的可加性，资源总需要活动所耗资源的可加性，资源总需要量为多种活动所需资源数量的总和。量为多种活动所需资源数量的总和。4.明确性明确性 目标的明确性目标的明确性5.单一性单一性 期望值的单一性期望值的单一性6.独立性独立性 变量是独立的表示各种作业
8、对资源都变量是独立的表示各种作业对资源都是互竟关系，没有互助关系是互竟关系，没有互助关系7.非负性非负性第二节 线性规划模型的建立与图解法求解一、建模二、线性规划的求解图解法一、建模例例1某饲料公司用甲、乙两种原料配制饲料，甲乙两种原料的营养成份及配合饲料中所含各营养成份最低量由表1给出。已知单位甲、乙原料的价格分别为10元和20元，求满足营养需要的饲料最小成本配方。一、建模 设配合饲料中，用甲x1单位，用乙x2单位，则配合饲料的原料成本函数，即决策的目标函数为Z=10 x1+20 x2。考虑三种营养含量限制条件后，可得这一问题的线性规划模型如下：Min Z=10 x1+20 x2 x1+x2
9、10 3x1+x215 x1+6x215 x10,x20一、建模例例2某农户计划用12公顷耕地生产玉米，大豆和地瓜，可投入48个劳动日，资金360元。生产玉米1公顷，需6个劳动日，资金36元，可获净收入200元；生产1公顷大豆，需6个劳动日，资金24元，可获净收入150元；生产1公顷地瓜需2个劳动日，资金18元，可获净收入1200元，问怎样安排才能使总的净收入最高。设种玉米，大豆和地瓜的数量分别为x1、x2和x3公顷，根据问题建立线性规划问题模型如下：一、建模Max Z=200 x1+150 x2+100 x3 x1+x2+x312(1)6x1+6x2+2x348(2)36x1+24x2+18
10、x3360(3)x10，x20，x30 一、建模 例例33某农户有耕地20公顷，可采用甲乙两种种植方式。甲种植方式每公顷需投资280元，每公顷投工6个，可获收入1000元，乙方式每公顷需投资150元，劳动15个工日，可获收入1200元，该户共有可用资金4200元、240个劳动工日。问如何安排甲乙两种方式的生产，可使总收入最大?解：设甲方式种x1公顷，乙方式种x2公顷，总收入为Z，则有：一、建模Max Z=1000 x1+1200 x2 280 x1+150 x24200 6x1+15x2240 x1+x220 x10,x20 二、线性规划的求解图解法（一）可行解（二）可行域（三）最优解（四）最
11、优性定理（五）最大化问题的图解法（六）最小化问题的图解法 二、线性规划的求解图解法（一）可行解 线性规划问题的可行解是指，满足规划中所有约束条件及非负约束的决策变量的一组取值，其仅与约束条件有关而与目标函数值的大小无关。（二）可行域 可行域是由所有可行解构成的集合。根据线性规划的基本理论，任一个线性规划问题的可行域，都是一个有限或无限的凸多边形，凸多边形的每个角，称为可行域的极点。（三）最优解 线性规划的最优解是指，使目标函数值达到最优(最大或最小)的可行解。一个线性规划问题可以是有解的，也可能是无解的，最优解的个数可能是惟一的，也可能是有无穷多个，即决策变量有许多组不同的取值，都使目标函数达
12、到同一个最优值。二、线性规划的求解图解法 （四）最优性定理 若一个线性规划问题有最优解，则最优解一定可以在可行域的某个极点上找到一个最优解。同时仍有可能有其他最优解存在，但它们也只可能存在于可行域的其他极点或是边界上。如果我们的目的是找出一个最优解而不是全部最优解，这一定理实际上是把寻找的范围，从可行域中的无穷多个可行点，缩小到可行域的有限几个极点上。二、线性规划的求解图解法 （五）最大化问题的图解法第一步，找出问题的可行域第二步，在可行域中寻求最优解，方法有两种：A.查点法 B.图解法二、线性规划的求解图解法 O 20 40 x120ABCD280 x1+150 x2=42006x1+15x
13、2=240 x1+x2=20 x2Z=1000 x1+1200 x2A(0,16)B(6.7,13.3)C(9.2,10.8)D(15,0)ZA=19200ZB=22660ZC=22160ZD=15000二、线性规划的求解图解法 （五）最小化问题的图解法w 例：Min Z=10 x1+20 x2w s.t.x1+x210w 3x1+x215w x1+6x215w x10,x201515105105OABCDx2x1x1+6x2=15可行域3x1+x2=15x1+x2=1010 x1+20 x20A(0,15)B(2.5,7.5)C(9,1)D(15,0)ZA=300ZB=175ZC=110ZD
14、=150第三节 单纯形法 单纯形方法是一种较为完善的、步骤化的线性规划问题求解方法。它的原理涉及到较多的数学理论上的推导和证明，我们在此仅介绍这种方法的具体操作步骤及每一步的经济上的含义。为更好地说明问题，我们仍结合实例介绍这种方法 第三节 单纯形法一、线性规划的标准型二、线性规划问题的解三、单纯形法 四、单纯型表第三节 单纯形法一线性规划的标准型LP目标函数有的要求实现最大化，有的要求实现最小化，约束条件可以是“=”、“”，这种多样性给讨论问题带来不便。为了便于讨论，我们规定线性规划问题的标准形式为：Max Z=c1x1+c2x2+c3x3+cnxn a11x1+a12x2+a1nxn=b1
15、 (1)a21x1+a22x2+a2nxn=b2 (2)am1x1+am2x2+amnxn=bm (m)x1，x2，xn0 第三节 单纯形法其简缩形式为 一线性规划的标准型njxbxaxcxcxcZjinjjijnn,3,2,1,0max12211用向量表示 01jnjjjxbxPnxxxX21mjjjjaaaP21bmbbb21 其中 C=(c1,c2,cn)向量Pj是其对应变量xj 的系数向量。第三节 单纯形法一线性规划的标准型用矩阵描述CXzMax0XbAX mnmnmmnbbbbbPPPaaaaaaA321212111211;第三节 单纯形法二线性规划问题的解 可行解可行解最优解最优解
16、 基基 设A为约束方程组的mn阶系数矩阵，其秩为m。B是矩阵A中mm阶非奇异子矩阵（）,则称B是线性规划问题的一个基。不失一般性可设 0BmmmmmmPPPaaaaaaB212111211称Pj为基向量，与基变量Pj相对应的变量为基变量。否则为非基变量。w 为了进一步讨论线性规划问题的解，我们来研究约束方程组求解的问题。假设方程组系数矩阵Z的秩为m，因m小于n故它有无穷多个解。假设前m个变量的系数列向量是线性独立的，这时线性规划模型可写成:二线性规划问题的解 nmnnnmmmmmmmmmmmmmxaaaxaaabbbxaaaxaaaxaaa211112112121122212112111nmj
17、jjmjjjxPbxP11或 021nmmxxxTmxxxX)0,0,(21设非基变量用高斯消去法，可求出一个解称X为基本解基本解基本可行解基本可行解 满足非负条件的基本解二线性规划问题的解 w 例例3某工厂在计划期内安排生产x1 x2两种产品，这些产品分别需要在A、B、C、D四种不同的设备上加工。按工艺规定，产品x1和产品x2在各设备上加工的台时数见下表。已知各设备在计划期内有效台时数分别是12、8、16和12。（一台设备工作一小时称为一台时）该工厂每生产一件产品x1可得利润2元，每生产一件产品x2可得利润3元，问如何安排生产计划，才能得到利润最多？三单纯形法 设备产品ABCDx12140
18、x22204三单纯形法w（一）求解过程 w（二）求解过程小结三单纯形法Max Z=2x1+3 x2 2x1+2x212 x1+2x28 4x1 16 4x2 12 x10，x20 引入松弛变量x3 A设备闲置台时数x4 B设备闲置台时数x5 C设备闲置台时数x6D设备闲置台时数将线性规划化为标准型.（8.1）三单纯形法 求解过程Max Z=2x1+3 x2+x3+x4+x5+x6 2x1+2x2+x3 =12 x1+2x2 +x4 =8 4x1 +x5 =16 4x2 +x6=12 x10，x20,x30，x40,x50，x60（8.2）三单纯形法 求解过程 x3,x4,x5,x6的系数列向量
19、p3,p4,p5,p6是线性独立的，这些列向量构成一个基 系数矩阵 100040010004001021000122654321PPPPPPA10000100001000016543PPPPB三单纯形法 求解过程x3 =122x12x2 x4 =8x12x2 x5 =164x1 x6 =124x2 把上式带入目标函数得到 Z=0+2x1+3 x2 （8.4）当非基变量x1=x2=0,便得z=0，这时得到一个基本可行解X(0)对应于B的变量x3,x4,x5,x6为基变量，从标准型我们可以得到：（8.3）三单纯形法 求解过程TxxxxxxX121681200121681200)0(6)0(5)0(
20、4)0(3)0(2)0(1)0(这个基本可行解表示：工厂没有安排生产产品；设备的有效台时数没有被利用，所以构成的利润为0。从分析目标函数的表达式可以看到，非基变量x1,x2系数都是正数，若将非基变量换成基变量，目标函数就会增加。所以，只要在目标函数的表达式中还存在正系数的非基变量，这表示目标函数还有增加的可能，就需要将非基变量换成基变量。一般选择正系数最大的那个非基变量。可按以下方法来确定换出变量。三单纯形法 求解过程 现分析（8.4），将x2定为换入变量后，必须从x3,x4,x5,x6中换出一个，并保证其余的都是非负，即x3,x4,x5,x60 当x10，由（8.3）式得到 x3 =122x
21、2 0 x4 =82x20 （8.5）x5 =160 x6 =124x2 0 从（8.5）式中可以看出，只有选择 Z=0+2x1+3 x2 （8.4）3412,28,212min2x时，才能使（8.5）式成立。因当x2=3时，基变量x6=0这就决定用x2去替换x6。三单纯形法 求解过程为了求得以x3,x4,x5,x2为基变量的一个基本可行解和进一步分析问题，需将（8.5）中的x2位置与x6的位置兑换。得到x3 2x2 =122x1 x4 2x2 =8x1 （8.6）x5 =164x1 4x2 =12 x6 用高斯消去法，将（8.6）式中的x2的系数列向量变为单位列向量。x3 =62x1+1/2
22、x6 x4 =2x1+1/2x6 （8.7）x5 =164x1 x2 =31/4x6三单纯形法 求解过程w 再将（8.7）代入（8.1）目标函数得到：w Z=9+2x1-3/4 x6 （8.8）w 当非基变量x1=x6=0，得到Z=9，并得到另一个基本可行解TTxxxxxxX)0,16,2,6,3,0(),()1(6)1(5)1(4)1(3)1(2)1(1)1(三单纯形法 求解过程 从目标函数的表达式（8.8）中可看到，非基变量x1的系数是正的，说明目标函数值还可以增大，X(1)不一定是最优解。于是用上述方法，确定换入换出变量，继续迭代，再得到另一个基本可行解X(2)再经过一次迭代，又得到一个
23、基本可行解 这时得到的目标函数的表达式是：Z=14-1.5x4-0.125 x5 目标函数值达到最大，X(3)是线性规划的最优解。TX)0,8,0,2,3,2()2(TX)4,0,0,0,2,4()3(三单纯形法 求解过程w 1.人造基、初始基本可行解 w 2.最优性检验 三单纯形法 求解过程小结w 1.人造基、初始基本可行解 w 1.1若从线性规划问题的 Pj中能直接观察到存在m个线性独立的单位向量，经过重新安排次序便得到一个可行基jnjjxczMax101jnjjjxbxP10001000121MPPPB三单纯形法 求解过程小结w 1.人造基、初始基本可行解 w 1.2“”标准化的方法，引
24、入非负的松弛变量重新对xj及aij编号，经整理则可得到下列方程 Max Z=c1x1+c2x2+c3x3+cnxn x1 +a1m+1xm+1+a1m+2xm+2+a1nxn =b1 x2 +a2m+1x m+1+a2m+2x m+2+a2nxn=b2 （8.9）xm +amm+1x m+1+amm+2x m+2+amnxn=bm x1，x2，xn0显然得到一个单位阵 10001000121MPPPB三单纯形法 求解过程小结我们就将B作为可行基。将（8.9）每个等式进行移项得 x1=b1 -a1m+1xm+1-a1m+2xm+2-a1nxn x2=b2-a2m+1x m+1-a2m+2x m+
25、2-a2nxn （8.10）xm=bm-amm+1x m+1-amm+2x m+2-amnxn x1，x2，xn0令x m+1=x m+2=x n=0,由（8.10）可得xi=bi(I=1,2,m)得到一个初始基本可行解TmnmTmnmbbbxxxX)0,0,()0,0,(2121个个，三单纯形法 求解过程小结2.最优性检验w 得到初始可行解后，要检验一下是否是最优解，如果是，则停止迭代，如果不是，则继续迭代。但每次迭代后都要检验一下是否是最优解，为此需要建立一个判别准则。w 一般情况下，经过迭代后式变成w （i=1,2,3,m）w 将上式代入目标函数，整理后得nmjjijiixabx1三单纯
26、形法 求解过程小结nmjjijmiijimiiixaccbcz111)(ijmiijimiiaczbcz110,j=m+1,n nmjjjjxzczz10)(),1(nmjzcjjjnmjjjxzz10三单纯形法 求解过程小结w 2.1最优解判别定理：w 若 为对应于B的基本可行解，且对于一切j=m+1,n有 ，则X（0）为最优解。w 无有限最优解判别定理：w 若 为对应于B的基本可行解，有一个 并且对于一切i=1,2,3,m有，那么该线性规划没有有限最优解。w 2.2换入变量的确定w 2.3换出变量的确定 ，w 为换入变量。TmbbbX)0,0(,2,1)0(0j0kmTmbbbX)0,0(
27、,2,1)0(0,kmia为换入变量则对应的kkjx)0max(ckcikikabaabRi0mincx三单纯形法 求解过程小结三单纯形表例1例1例1例1例2例2例2w 目标函数系数灵敏度分析w 右边值灵敏度分析第四节 灵敏度分析目标函数系数灵敏度分析w 最优解不变的条件下，允许C的变化范围，最优解不变的前提是j 0w 假设玉米价值系数C1发生了变化，其变化量为1 x x1 1x x2 2x x3 3x x4 4x x5 5x x6 61 10 00 0 x x3 36 60 00 01 13 3/2 2-0 0.2 25 50 02 20 00 0+1 1x x1 16 61 11 10 0
28、-0 0.5 51 1/4 40 00 0 x x6 63 36 60 0-1 12 20 0-9 9-4 4.5 51 12 20 00 0+1 11 15 50 01 10 00 00 00 00 01 18 80 00 02 20 00 0+1 12 20 00 0+1 11 10 00 05 50 0-0 0.5 51 12 25 5+0 0.2 25 51 10 01 18 80 00 00 0-5 50 0-1 10 0-5 50 0+0 0.5 51 1-2 25 5-0 0.2 25 51 10 0实实际际活活动动松松弛弛活活动动目目标标系系数数行行c cj j机机会会成成本本
29、行行Z Zj j第第三三单单纯纯形形表表C CB BX XB Bb b检检验验数数行行j j-50 1 1001 -501 1001 100 -50-1 0 -50+0.51 0 -25-0.251 0目标函数系数灵敏度分析w 假设大豆价值系数C2发生了变化，其变化量为2 2 -50-502 0目标函数系数灵敏度分析w 假设地瓜价值系数C3发生了变化，其变化量为3-100/3 3 1003 -100/33 100 -50-1.53 0-250.253 0表8-3 目标系数的允许变动范围活动目标系数当前值可减上限 可增上限 可变范围玉米种植x120050100150300大豆种植x2150无穷大
30、50-200地瓜种植x3100100/3100200/3200 当仅有一种目标系数在允许范围内变动时，最优方案不会变动，但最优目标值会随之变化。右边值敏感性分析w 由线性规划的原理可知，影子价格不变的条件是最优解的松弛变量矩阵与右边值矩阵的乘积大于和等于0，即：0360481212/9904/12/104/12/301PBX当右边值发生变化时，如耕地变化，此时，影子价格不变的条件是 右边值敏感性分析耕地0360481212/9904/12/104/12/3101PBX得到 6+3/210 6-1/210 -414 36-910因此耕地影子价格不变的耕地数量范围为：8，160360481212/
31、9904/12/104/12/3201PBX右边值敏感性分析劳动力得到 61/420 6+1/420 -2428 36-9/220因此劳动力影子价格不变的劳动力数量范围为：24，56右边值敏感性分析资金0360481212/9904/12/104/12/3301PBX得到 60 60 3 36 3630因此资金影子价格不变的资金数量范围为：324，表8-4 常数项的允许变动范围资源现有数量可减上限可增上限可变范围影子价格耕地124481650劳动48248245625资金36036+324+0常数项的允许变动范围这一结果也还有另外一种意义，即它给出了资源影子价格（边际产出率）的有效范围。对耕地而言，当投入使用的数量在8-16公顷之间变化时，其边际产出率都是50元，即每增加或减少1公顷耕地，农户将增加或减少50元净收入。
