第八章地理系统的线性规划.pptVIP

如果记 以及 则 (二) 单纯形法的计算步骤 第1步，找出初始可行基，建立初始单纯形表。 第2步，判别检验所有的检验系数 （1）如果所有的检验系数 , 则由最优性判定定理知，已获最优解，即此时的基本可行解就是最优解。 （2）若检验系数中，有些为正数，但其中某一正的检验系数所对应的列向量的各分量均非正，则线性规划问题无解。 （3）若检验系数中，有些为正数，且它们所对应的列向量中有正的分量，则需要换基、进行迭代运算。 第3步，选主元。 在所有大于零的检验数中选取最大的一个b0s，对应的非基变量为xs，对应的列向量为 若 则确定brs为主元项。 第4步，在基B中调进Ps，换出Pjr，得到一个新的基 第5步，在单纯形表上进行初等行变换，使第s列向量变为单位向量，又得一张新的单纯形表。 第6步，转入上述第2步。 例1：用单纯形方法求解线性规划问题 解：首先引入松弛变量 ，把原问题化为 标准形式 具体步骤如下： 第1步，确定初始单纯形表。 x1 x2 x3 x4 -z 0 2 3 0 0 x3 12 1 3 1 0 x4 9 2 1 0 1 第2步，判别。在初始单纯形表中b01=2, b02=3，所以B1不是最优基，进行换基迭代。 第3步，选主元。 根据选主元法则，确定主元项 。 第4步 ，换基，得一新基 。 第5步，进行初等行变换， 得B2下的新单纯形表 x1 x2 x3 x4 -z -12 1 0 -1 0 x2 4 1/3 1 1/3 0 x4 5 5/3 0 -1/3 1 第6步，因检验系数有正数b01=1，重复以上步骤可得对应于 B3=[p2,p3]的单纯形表，检验各检验数可知得最优解X1=3,X2=3, X3=0, X4=0:目标函数最大值为 Z=15。 x1 x2 x3 x4 -z -15 0 0 -4/5 -3/5 x2 3 0 1 2/5 -1/5 x1 3 1 0 -1/5 3/5 六、应用实例: 农场种植计划模型 某农场I、II、III等耕地的面积分别为100 hm2、300 hm2和200 hm2，计划种植水稻、大豆和玉米，要求3种作物的最低收获量分别为190 000 kg、130 000 kg和350 000 kg。I、II、III等耕地种植3种作物的单产如表5.1.4所示。若3种作物的售价分别为水稻1.20元/kg，大豆1.50元/ kg，玉米0.80元/kg。那么，（1）如何制订种植计划，才能使总产量最大？（2）如何制订种植计划，才能使总产值最大？ 表 不同等级耕地种植不同作物的单产(单位:kg / hm2) ? I等耕地 II等耕地 III等耕地 水稻 11 000 9 500 9 000 大豆 8 000 6 800 6 000 玉米 14 000 12 000 10 000 表5.1.4 对于上面的农场种植计划问题，我们可以用线性规划方法建立模型。 根据题意，决策变量设置如表所示, 表中Xij表示在第j等级的耕地上种植第i种作物的面积。 3种作物的产量可以用表表示。 作物种类 总产量 水稻 大豆 玉米 表 作物计划种植面积（单位:hm2) 表 3种作物的总产量（单位:kg） ? I等耕地 II等耕地 III等耕地 水稻 大豆 玉米 根据题意，约束方程如下,耕地面积约束 最低收获量约束 非负约束 （1）追求最大总产量的目标函数为 调用Matlab软件系统优化工具箱中的linprog函数，进行求解运算，可以得到一个最优解（如表5.1.7所示）。在该方案下，最优值，即最大总产量为6 892 200 kg。从表中可以看出，如果以追求总产量最大为种植计划目标，那么，玉米的种植面积在I、II、III等耕地上都占绝对优势。 表 追求总产量最大的计划方案（单位：hm2） ? I等耕地 II等耕地 III等耕地 水稻 0 0 21.111 1 大豆 0 0 21.666 7 玉米 100 300 157.222 2 (2) 追求最大总产值的目标函数为 进行求解运算，可以得到一个最优解（如下所示）。在该方案下，最优值，即最大总产值为6 830 500元。从表中可以看出，如果以追求总产值最大为种植计划目标，那么，水稻的种植面积在I、II、III等耕地上都占绝对优势。 表 追求总产值最大的计划方案（单位：