第二次作业wyc.docVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第二次作业wyc

数学实验报告 题目线性规划问题 班级 理科数学91 学号 姓名 王一辰 日期 2010-10-24 实验题目 (1).森林管理问题: 森林管理问题:森林中的树木每年都要有一批被砍伐出售.为了使这片森林不被耗尽且每年都有所收获,每当砍伐一棵树时,应该就地补种一棵幼苗,使森林树木的总数保持不变.被出售的树木,其价值取决于树木的高度.开始时森林中的树木有着不同的高度.我们希望能找到一个方案,在维持收获的前提下,如何砍伐树木,才能使被砍伐的树木获得最大的经济价值? 已知森林具有6年的生长周期=0.28, =0.32, =0.25, =0.23, =0.37, =50元,=100元,=150元,=200元,=250元,求出对其进行最优采伐的策略。 (2).利润最大化问题: 某厂生产一种产品有甲、乙两个牌号,讨论在产销平衡的情况下如何确定各自的产量,使总利润最大. 所谓产销平衡指工厂的产量等于市场上的销量。 二.问题分析 问题一 1.首先,对于简单的模型,假设如下 (1)每年对森林中树木砍伐一次,留下的树木和补种的幼苗,经过一年的生长期后,与上一次砍伐前的高度状态相同(根据保持持续收获的要求)。 (2)在一年的生长期内树木最多只能生长一个高度级,即第k类的树木可能进入k+1类,也可能停留在k类中。 (3)树木没有死亡。 2.建立恰当的数学模型 步骤1:把森林中的树木按照高度分为n类,第k(1kn)类树木的高度为[,],每棵的经济价值为,而第1类树木即幼苗的高度为[0,],每棵的经济价值=0;第n类树木高度为[,)。 步骤2:记(k= 1,2,…,n)是第t年第k类树木的数量, (k= 1,2,…,n)是第k类树木被砍伐的数量。 步骤3:设是经过一年从第k类中长高到第k+1类中的树木比例,其中01。而1-是经过一年仍然留在第k类中的树木比例。 步骤4:设森林中树木的总数是s。 步骤5:令为数量向量,为收获变量, 为种植矩阵,生长矩阵。则 = 由于树木的生长必须维持平衡关系:生长期末的状态减去收获采伐的量后再加上补种的幼苗数应等于生长期开始的量。即 因为,,所以 步骤6:利用收获向量和价值向量得所收获树木的价值为 步骤7:综上所述,此模型为: , 符号说明如下; c=[,,,,,0]; [ 1 1 1 1 1 1 ] 3.模型的改进 在本题中,我们一开始假设树木没有死亡,而树木在生长环境,气候,天敌等诸多因素的影响下会发生死亡,为了更接近实际情况,我们对于每一类树木引进死亡率。假设在每个单位时间内这六类树木的死亡率分别为,,,,,则按照上述思想,可以得到改进的线性规划模型为: 按照相同的思路求解,MATLAB语句中只需修改相应的参量。 问题二 1.问题的假设 (1)价格与销量成线性关系。 (2)成本与产量成负指数关系。 2.模型的建立 (1)价格与销量关系式的建立 利润既取决于销量和价格,也依赖于产量和成本。按照市场规律,甲的价格会随其销量的增长而降低,同时乙的销量的增长也会使甲的价格有稍微的下降,可以简单地假设价格与销量成线性关系: (2)成本与产量关系式的建立 甲,乙的成本随其产量的增长而降低,且有一个渐进值,可以假设为负指数关系: (3) 总利润的的计算式 总利润,对于本题,不妨取: , ,, ,. 三.程序 问题一 (1)简单模型 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; max (2)改进的模型 format long c=[-14,-15.5,-11.5,-10.0,-16.5,0] A=[-0.28,0.33 ,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.38,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; max 问题

文档评论(0)

2017ll + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档