蒋杰的线性规划问题报告.pdfVIP

数学实验报告 题目 线性规划问题 班级 理科数学91 学号 姓名 蒋 杰 日期 2010/10/23 实验题目1 1. 森林管理问题：森林中的树木每年都要有一批被砍伐出售.为了使这片森林不被耗尽且 每年都有所收获,每当砍伐一棵树时,应该就地补种一棵幼苗,使森林树木的总数保持不 变.被出售的树木,其价值取决于树木的高度.开始时森林中的树木有着不同的高度.我 们希望能找到一个方案,在维持收获的前提下,如何砍伐树木,才能使被砍伐的树木获得 最大的经济价值? 思考题1.已知森林具有6 年的生长期g 1 0.28 ，g 2 0.32 ，g 3 0.25 ，g 4 0.23 ， g 5 0.37 ，p 2 50 元，p 3 100 元，p 4 150 元，p 5 200元，p 6 250 元， 求出对其最优的采伐策略。 实验目的 通过对森林采伐问题的操作运算，达到对如何解线性规划问题的掌握，并对实际问题中的简 单的规划问题能够予以解答。 问题分析 如老师的 ppt 所示的分析，即构建出线性方程及其约束条件，由此可以用 matlab 中的 linprog 命令求得该线性方程的最优解。 算法设计 建立如下的线性方程组： max 14* x116* x2 12.5* x3 11.5* x4 18.5* x5 0* x6 0.28* x10.32* x 2 0; 0.32*x 2 0.25* x3 0; 0.25* x3 0.23* x 4 0; 0.23* x 4 0.37* x5 0; x1x 2x 3x 4x 5x 6 s xj 0, j 1,2,3,4,5,6 其中s 是森林中的总的数目数。 程序 format long c=[-14,-16,-12.5,-11.5,-18.5,0]; A=[-0.28,0.32,0,0,0,0; 0,-0.32,0.25,0,0,0; 0,0,-0.25,0.23,0,0; 0,0,0,-0.23,0.37,0]; b=[0;0;0;0]; aeq=[1,1,1,1,1,1]; beq=[10000]; vlb=[0;0;0;0;0;0]; vub=[]; [x,fval]=linprog(c,A,b,aeq,beq,vlb,vub) max=-fval; fprintf(%d是最大收益\n,max) 计算结果 x = 1.0e+003 * 5.333333333369788 4.666666666325552 0.000000000209711 0.000000000072755 0.000000000018550 0.000000000003646 fval = -1.493333333321872e+005 1.493333e+005 是最大收益 结果分析 如上的结果是可以让人可信的，可以看到，当想 x1=5333.3，x2=4666.7，x3=x4=x5=x6=0 时，该方程组取得最大值14933 元的收益。 结论 可以获得的最大收益为14933 元。 改进： 将假设没有死亡改进为存在死亡，即每一种都有一定的死亡率，为此做出的模型应该更加符 合实际中的情况。 实验题目 加上死亡的比率参数即可：设 s1=0.05，s2=0.03，s3=0.02,s4=0.01，s5=0.01，s6=0.01 分别表示六个等级上的数目的死亡率 问题分析 同样利用上述的方法得到平衡的解，然后就可以建立一个线性规划问题，最后予以解决。 算法设计 max 14* x1 13* x2 11.5* x3 9.5* x4 16* x5 0* x6 0.28* x10.05* x1 0; 0.28* x10.32* x 2 0.03* x 2 0; 0.32* x 2 0.25* x3 0.02* x3 0; 0.25* x3