有关养牛的农场规划问题 数学建模.docx

1、有关养牛的农场规划问题 数学建模【摘要】本题探究的是如何对农场5年的生产计划作出决策，我们将运用目标优化模型进行求解。据题意可知第j年卖0岁母牛的数量Nj（Af=0），第j年种地的亩数Sj（在解题时，我们以买卖饲料的数量以及养牛需要的饲料数量来决定），以及贷款的总额M为三个主要的决策变量。其中第j年卖0岁母牛的数量Nj（Af=0）会影响每一年牛的总数进而影响种地的亩数，同时也会影响农场主的贷款额。所以，本题我们将以第j年牛的总头数Nj，年收入Wj，年成本Cj作为对象，并以5年的总收益的最大值Z作为目标函数进行讨论。为得到五年的净收益总和，我们将用每年的收入Wj减去每年的成本Cj得到每年的净收益

2、并求和来得到。对于成本中牛的数量超过130时的额外投资费用，我们用年初出生的小母牛的头数减去年初卖掉的小母牛的头数乘以2000来得到；在决定贷款额时，我们首先对每年除去还款额后的净收益进行粗略计算得到每一年都是可以盈利的，所以我们将使用第一年的成本Cj作为贷款总额。最后，运用lingo软件对决策变量进行规划得到的结果。我们通过求解得到5年的总收益的最大值Z=686625.6，贷款总额M=416055.8，种粮食亩数S1=80，种甜菜亩数S2=120，各年卖出的0岁母牛头数分别为48,0,15,60,90。其他因素的变化对计划造成的影响，其中银行利率r的变化会对总收益M造成影响，但对其他决策变量

3、影响较小；另外，如果农产品价格和产量以及劳动力价格发生变化，将有可能改变种农产品的亩数和购买农产品的数量，贷款额M和相应的总收益Z的变化，各因素的具体影响我们方式将在模型解答中加以阐述。【关键词】 优化模型 贷款总额M 卖掉的0岁母牛的头数 买、卖的饲料吨数 种植饲料的土地面积 lingo软件 一、问题重述某公司计划承包有200亩土地的农场，建立奶牛场。开始承包时农场有120头母牛，其中20头为不到01岁的幼牛，100头为产奶牛。产奶牛平均每头每年生0.55头公牛，生出后不久即以每头300元卖掉；产奶牛平均每头每年生0.05头的母牛要么出生后以400元卖掉，要么饲养，养至2岁成为产奶牛。幼牛年

4、损失5%；产奶牛年损失2%。产奶牛满12岁就以每头1200元卖掉。这120头牛中从0岁到11岁每个年龄阶段各有10头。应该卖掉的小母牛都已卖掉。只有20头是要饲养成产奶牛的。一头牛产奶年收入为3700元。现在农场最多只能养130头牛。超过此数每多养一头，要投资2000元。每头产奶牛每年消耗0.6吨粮食和0.7吨甜菜。每头小牛每年消耗粮食和甜菜量为奶牛的2/3。粮食和甜菜可以自己种植，每亩产甜菜1.5吨。只有80亩的土地适于种粮食，产量平均0.9吨。从市场购粮食每吨900元，卖出750元。买甜菜每吨700元，卖出500元。养牛和种植所需的劳动量为：每头小牛每年10小时；每头产奶牛每年42小时；种

5、一亩粮食每年需20小时；种一亩甜菜每年需30小时。其它费用：每头幼牛每年500元，产奶牛每头每年1000元；种粮食每亩每年150元，种甜菜每亩每年100元。劳动力成本为每小时费用为10元承包农场需要一笔费用，其中租地费用为每年6万元（每年底付清），以及开始承包时农场已有的120头牛的费用。平均产奶牛每头4000元，小牛每头400元，到承包结束时，农场的牛按此价折价抵卖。 该公司资金全部来自银行贷款，复利计息，每年归还本利和的1/5。五年还清。此外，农场主不希望产奶牛的数目在五年末与现在相比减少超过50%，也不希望增加超过75%。现在要求分析承包人有无盈利的可能性。若有，应如何安排5年的生产，使

6、得五年的净收益为最大，并以此为基础进一步讨论，若遇到银行利率波动（例如上下波动2个百分点），还贷方式改变（如规定每年还息，改变还本的方式），由于气候等外因变化引起的农产品产量与价格的变化及劳动力市场价格的变动等将会对农场主的五年生产计划及收益产生怎样的影响。二、问题分析对于问题一：本题要求分析农场主5年后能否盈利，并制定5年的生产计划设法使净收益最大，我们认为这是一个研究资源的合理配置和优化问题。由题意可知，农场的收入来源有卖牛奶的收入，卖牛的收入，以及卖出多余的粮食的收入，同时成本来源包括雇佣劳动力产生的费用，地租，养牛成本（养牛总数高出130头时追加的投资费，饲料费，照顾费等），还贷款以及

7、其他一些费用。由于题中所给的资料均以年为单位，因此，可以通过得出年收益的通项并求和就可得到5年的总收益。经分析，我们可以知道，要想规划每年的生产计划就必须要对每年卖牛的数量，种地的亩数（在解题时，我们以买卖饲料的数量以及养牛需要的饲料数量来决定），以及贷款的总额做出决策。由此可见，其中最为总要的是卖牛的数量，因为它会影响一年牛的总数进而影响种地的亩数，同时也会影响农场主的贷款额。本题我们将以牛的头数以及总收益作为主要的建模对象，来对问题进行分析和阐述。对于问题二：银行利率等贷款方式的波动会在一定程度上影响贷款的额度，如果实际利率高，将考虑精简贷款额度，反之，可能适当提高贷款额度以获得更高收益；

8、农产品价格和产量的变化可能导致在考虑种植不同种饲料（粮食和甜菜）的亩数的决策上有所影响；最后，劳动力市场的价格变化可能对不同的牛（12岁母牛，0岁母牛）的养殖数目，种植饲料的亩数，以及间接地对贷款额造成一定影响。此题我们将结合模型中利率，农产品产量和价格，及劳动力价格这几个变量来进行分析。三、模型假设1、幼牛出生发生在年初，损失发生在年末。2、刚出生的小牛无损失，即刚出生就卖掉，12岁的老母牛在年末被卖掉。3、幼牛损失的各年龄和奶牛损失的各年龄是均匀的。4、银行按照复利计息，即农场每年归还1/5(1+r)5M5、相邻两个年龄组的牛在相邻两年之间的变化是连续的。6、在5年末12岁的牛被以4000

9、的价格收购。7、粮食是在每年末才能成熟并收割，即第一年的牛的全部饲料都从市场上购买，以后每年的牛吃上一年种的粮食，不够的再到市场上购买。8、年底才能拿到所有的收入，而支出发生在年初。9、不会卖一岁的母牛四、符号说明1、Wj第j年总收益， Wkj卖奶收益，Wcj卖牛收益，Wqj卖粮食收益2、 Nj 是第j年初母牛的总数量，Nij 是第j年初第i岁母牛的数量，F1j第j年初成年母牛总数，F2j第j年初未成年母牛总数3、Nj(A=12) 卖出12岁母牛总数，Nj(Am=0)卖出公牛总数，Nj(Af=0)卖出幼年母牛总数4、Clj劳动力花费, Csj养牛花费, Rsj地租, Cbj贷款本利和，Caj为

10、第j年牛的头数超过130头时的额外投资，Cj为第j年末总花费5、Pl为劳动力单位时间费用， Pr1为粮食卖出价格，Pr2为粮食买入价格，Ps1为甜菜卖出价格，Ps2为甜菜买入价格6、Wrj是卖粮食重量，Wsj是卖甜菜的重量，Wmj是购买粮食的重量，Wej是购买甜菜的重量7、S1j是第j年种植粮食的面积 S2j是第j年种植甜菜的面积8、O1为每亩粮食产量，O2为每亩甜菜产量9、W3j是第j年需要的粮食吨数，W4j是第j年需要的甜菜吨数10、W5第五年当年净收益，Z为5年的最大净收益11、M为贷款总额，利率r五、模型的建立（一）目标函数： ，M=C1（二）约束条件： 我们将贷款总额M，卖掉的0岁母

11、牛的头数Nj（Af=0），卖掉的粮食吨数Wrj，卖掉的甜菜吨数Wsj，买入的粮食吨数Wmj，买入的甜菜吨数Wej，六个变量作为决策变量，并用(W3j+1+WsjWmj)/O1 表示种植粮食的土地面积S1j，用(Wsj+W4j+1Wej)/O2表示种植甜菜的土地面积S2j，则（1）14年末收益：卖奶收益Wkj=F1j*3700 卖牛收益 Wcj=Nj(A=12)*1200+Nj(Am=0)*300+Nj(Af =0)*400 卖粮食和甜菜的收益Wqj=Wrj*Pr1+Wsj*Ps1总收益Wj=Wkj+Wcj+Wqj= F1j*3700+ Nj(A=12)*1200+Nj(Am=0)*300+Nj

12、(Af =0)*400+ Wrj*750+Wsj*50014年末花费：劳动力花费 Clj=(F1j*42+F2j*10+S1j*20+S2j*30)*P1 养牛花费Csj= F1j *1000+ F2j *500+Wmj*Pr2+Wej*Ps2 地租 Rsj=60000+S1j*150+S2j*100 贷款本利和Cbj = 1/5M*(1+r)5 r=12%(年利率) 超过130头年的额外投资 2=j=5第一年超过130头牛的额外投资 总花费Cj= Clj+Csj+ Rsj+ Cbj +Caj注：W3j = (F1j + F2j *2/3)*0.6，W4j =（F1j + F2j *2/3）*

13、0.7) （2）第五年净收益 （211岁成年母牛总数）， （12岁母牛总数）Nj(A=12)=N11j-1*98%Nj(Am=0)= Nj(Am=0) 是卖的0岁公牛总数 F11=100,F21=20+,Wr1=0,Ws1=0,Wm1=68,We1=80,r=12%,Pl=10,O1=0.9,O2=1.5,Pr1=900,Pr2=750,Ps1=700,Ps2=500, N1(Am=0)=55,N1(Af=12)=9.8六、模型求解（一）问题一的求解为简化模型，我们增设以下假设：假设一：前四年每一年都种80亩粮食，120亩甜菜，第五年种200亩甜菜。理由：我们从运行结果上可以看到第二年的牛的总

14、头数已经超过了130，经计算，当农场里恰好有100头成年母牛，30头未成年母牛时，农场中80亩地种植的粮食恰好够吃，且种粮食比买粮食更便宜，并且种植甜菜和粮食都能够获得收益。所以，在用程序进行计算时，我们对种植粮食做了规划，200亩地全部种完，由于粮食最多只能种80亩，所以，前4年80亩地种粮食，120亩地种甜菜；由于第五年吃的粮食来源于第四年的收成，而第5年的粮食到年底才有收成，由于甜菜每亩的收益比粮食大，所以第五年200亩的土地全部中甜菜。假设二：不购买甜菜理由：经计算，甜菜足够支撑起在题目范围内的母牛的增长率假设三：第一年的总成本C1等于贷款总额M理由：因为经过估算，几乎每一年都可以得到

15、正的净收益，并且，由于第一年的资金比较充裕至少会大于第一年的支出，因为年底才会有收成，所以，第一年的净收益就等于第一年的收入W1，会比较大足矣垫付之后每一年的亏损，所以，为达到支出最小化，我们只借相当于第一年的总成本额度的贷款。假设四：小母牛在年初产下理由：由于年初生牛和年末生牛的计算结果相差不大，为方便模型计算，设置在年初生牛。将已知数代入lingo，运行结果如下： Global optimal solution found. Objective value: 751029.9 Objective bound: 751029.9 Infeasibilities: 0.7331380E-10

16、Extended solver steps: 0 Total solver iterations: 10 Variable Value Reduced Cost Z1 2139.106 0.000000 Z2 126925.6 0.000000 Z3 348.9236 0.000000 Z4 397370.5 0.000000 Z5 330750.0 0.000000 X5 178.0000 -2134.520 S1 519798.7 0.000000 X1 100.0000 0.000000 R1 14.51357 0.000000 Y1 60.48643 0.000000 X2 98.00

17、000 0.000000 R2 48.30230 0.000000 Y2 53.55980 0.000000 X3 132.0000 588.9184 R3 0.8840093 0.000000 Y3 77.03380 0.000000 X4 125.0000 358.4026 R4 68.13019 0.000000 Y4 0.000000 1593.333 R5 97.90000 0.000000 Y5 0.000000 0.000000 S2 325892.8 0.000000 S3 547304.4 0.000000 S4 316200.0 0.000000 S5 424080.0 0

18、.000000 X6 174.9994 0.000000 Row Slack or Surplus Dual Price 1 751029.9 1.000000 2 0.000000 1.000000 3 0.000000 1.000000 4 0.000000 1.000000 5 0.000000 1.000000 6 0.000000 1.000000 7 0.000000 -0.7623417 8 0.000000 0.000000 9 0.000000 0.000000 10 0.000000 0.000000 11 0.000000 0.000000 12 0.000000 -81

19、67.200 13 0.000000 9121.901 14 0.000000 7496.838 15 0.000000 5237.673 16 0.000000 3560.480 17 0.000000 0.000000 18 14.51357 0.000000 19 48.30230 0.000000 20 0.8840093 0.000000 21 40.48643 0.000000 22 5.597695 0.000000 23 71.71599 0.000000 24 0.6198088 0.000000 25 0.000000 400.0000 26 0.000000 -1707.

20、546 27 0.000000 -1193.333 28 0.000000 -1193.333 29 0.000000 400.0000 30 0.000000 -1193.333 31 124.9994 0.000000 32 0.6225200E-03 0.00000 截图如下：绘制成表格如下：年份年初产奶牛头数年初卖出的幼年母牛年初卖出后剩余的幼母牛数年初牛总数卖出的12岁母牛头数年初卖牛后超过130头的牛的数目年末支出的费用当年一年的年末净利润第一年10015601609.830519798.72139.106第二年9848541529.622325892.8126925.6第三年13

21、21772099.479547304.4348.9236第四年1256801259.25316200397370.5第五年178980178048424080330750即可知：第五年末盈利Z=2139.106+126925.6+348.9236+397370.5+330750=751029.9元其中5年的生产计划如下贷款总额M=519798.7，前四年种粮食亩数S1=80，种甜菜亩数S2=120，第五年种植200亩甜菜，各年卖出的0岁母牛头数分别为Nj（Af=0）=15, Nj（Af=0）=48, Nj（Af=0）=1, Nj（Af=0）=68, Nj（Af=0）=98.（二）问题二的求解：

22、为简化模型，现增设以下假设：假设一：农场主养奶牛能够盈利假设二：在其他因素不变的情况下，即指处讨论的那个因素变量外，其他变量的取值都来源于第一题所给的数值。贷款总额M，卖掉的0岁母牛的头数Nj（Af=0），卖掉的粮食吨数Wrj，卖掉的甜菜吨数Wsj，买入的粮食吨数Wmj，买入的甜菜吨数Wej，六个最主要的决策变量，并用(W3j+1+WsjWmj)/O1 表示种植粮食的土地面积S1j，用(Wsj+W4j+1Wej)/O2表示种植甜菜的土地面积S2j。对他们造成影响的因素有以下几个Pl为劳动力单位时间费用， Pr1为粮食卖出价格，Pr2为粮食买入价格，Ps1为甜菜卖出价格，Ps2为甜菜买入价格，利

23、率r,O1为每亩粮食产量，O2为每亩甜菜产量 ，M=C1 Nj（Af=0）， S1j=(W3j+1+WsjWmj)/O1 ， S2j=(Wsj+W4j+1Wej)/O2W3j = (F1j + F2j *2/3)*0.6，W4j =（F1j + F2j *2/3）*0.7)， 卖牛收益 Wcj=Nj(A=12)*1200+Nj(Am=0)*300+Nj(Af =0)*400年末收益：卖奶收益Wkj=F1j*3700 卖牛收益 Wcj=Nj(A=12)*1200+Nj(Am=0)*300+Nj(Af =0)*400 卖粮食和甜菜的收益Wqj=Wrj*Pr1+Wsj*Ps1总收益Wj=Wkj+Wc

24、j+Wqj= F1j*3700+ Nj(A=12)*1200+Nj(Am=0)*300+Nj(Af =0)*400+ Wrj*750+Wsj*500年末花费：劳动力花费 Clj=(F1j*42+F2j*10+S1j*20+S2j*30)*P1 养牛花费Csj= F1j *1000+ F2j *500+Wmj*Pr2+Wej*Ps2 地租 Rsj=60000+S1j*150+S2j*100 贷款本利和Cbj = 1/5M*(1+r)5 r=12%(年利率) 超过130头年的额外投资 2=j=5第一年超过130头牛的额外投资 总花费Cj= Clj+Csj+ Rsj+ Cbj +Caj 由上式可以看

25、到，对于银行利率r上升会导致年贷款本利和Cbj的上升，进而导致年成本Cj的上升，从而引起总收益Z的降低，但由于贷款总额M是由第一年的成本C1决定，与银行利率没有关系，所以贷款总额M将不会变化，并且银行利率与其他决策变量没有直接或间接地关系，所以不会对其造成影响；如果其他因素不变，当甜菜的亩产O2低于0.56时应不种值，当0.56=O20.8时，才能从种植甜菜上盈利；同理，当粮食的亩产O1=0.39时不宜种植，当0.39=O10.46时，才能从种植粮食上盈利。另外，当亩产降低时将会增加成本，从而增加贷款额度M，减少总收益Z。 其他因素不变，当粮食的卖出价格Pr1388.9时，且此时若粮食买入价格

26、Pr2=Pr1，前四年应该种满80亩，并使牛吃的所有都从市场上购买，当粮食的买入价格，Pr1Pr2时，牛先吃自己种的粮食，不够的从市场上买；当粮食的卖出价格Pr1388.9时，且此时若粮食买入价格Pr2=Pr1，应不种，并使牛吃的所有都从市场上购买，当粮食的买入价格，Pr1Pr2时，牛先吃自己种的粮食，不够的从市场上买。同理，266.7为甜菜的买卖价格的临界点。此外，如果价格变动带来粮食收入的升高，但会将减少贷款额度M，增加总收益Z，卖奶牛的计划不会发生改变；反之，则会增加贷款额度M，减少总收益Z，卖奶牛的计划不会有改变。 如果其他因素不变，劳动力单位时间费用增加则成本增加，总收益Z减少，贷款

27、总额M增加；当P131.7时，前四年种满80亩粮食，种120亩甜菜，并先吃自己种的，不够再买；当31.7P33时，不种农作物，所有饲料都从市场买。同时卖奶牛的头数不发生变化。七、模型评价（一）优点：本模型适用广泛，对其他一些生物资源的开发具有同等可靠性，特别是对于线性的优化问题，lingo软件配合目标优化模型的准确性较高。此外，本模型在解决每一个小问题时，根据题设对模型做了额外的假设，简化了解题步骤，并对程序和模型本身进行了二次简化，使结果更加简单易懂，清晰明了。模型中我们建立了多个决策变量，我们尽量使用了比较直观的符号来表示变量名称，并简化了变量的表达式，使得每个变量之间的关系一目了然。而且

28、，本题主要以讨论牛的头数为出发点，层层深入，逻辑性更显严密。（二）缺点：在模型的优化前期，变量个数相对较多，需要耗费比较大的时间来分析变量之间的关系并对模型进行简化；由于结构非常的严密，逻辑性很强，变量很多，所以，读者需要密切注意每一个变量。八、模型推广以本模型主要解决如何达到制定的长远的目标优化问题，以此为理论基础，适当改变某些参数就能更多地推广应用到多种经济问题中、工作计划，科研，教育等多个领域及其它饲养场中的一些简单优化问题；如果对模型进一步改进，引入非线性的变量可能能够解决一些复杂的优化问题。例如，如何在边际效益递减的情况下，获得效益的最大化，又或者是诸如国家经济的发展计划的一些变量较

29、复杂，变量间关系较为复杂，呈现非正规线性关系的优化问题。参考文献1姜启源 谢金星 叶俊，数学模型（第四版），北京：高等教育出版社，2011年；2韩中庚，数学建模方法及其应用，北京市：高等教育出版社 ，2009年；3赵可培，运筹学，上海市：上海财经大学出版社 ，2000年；4吴祈宗，运筹学与最优化方法，北京市：机械工业出版社 ，2003年；5谢金星 薛毅，优化建模与LINDO/LINGO软件，北京：清华大学出版社，2005年。 附录：Lingo程序如下：model:max=z1+z2+z3+z4+z5-100*4000-20*400+x5*4000*0.98-(1+0.12)5-1)*s1;z1

30、=3700*x1+x1*0.55*300+r1*400+10*0.98*1200-(0.6*x1+0.4*y1)*900-20*80*10-150*80-120*30*10-120*100+(0.7*x1+0.7*2/3*y1)*500-y1*100-y1*500-x1*42*10-x1*1000-(x1*0.55-r1-10)*2000-60000;z2=3700*x2+x2*0.55*300+r2*400+10*0.98*0.98*1200-(0.6*x2+0.4*y2-72)*900-20*80*10-150*80-120*30*10-120*100+(180-0.7*x2-0.7*2/3*y2)*500-y2*100-y2*500-x2*42*10-x2*1000-(x2*0.55-r2)*2000-60000;z3=3700*x3+x3*0.55*300+r3*400+10*0.98*0.98*0.98*1200-(0.6*x3+0.4*y3-72)*900-20*80*10-150*80-120*30*10-120*100+(180-0.7*x3-0.7*2/3*y3)*500-y3*100-y3*500-x3*42*10-x3*1000-(x3*0.55-r3)*2000-60000;z4=3700*x4+x1*0.55*300+r4*400+10*0.