本题探究的是如何对农场5年的生产计划作出决策,我们将运用目标优化模型进
行求解据题意可知第j年卖0岁母牛的数量1(A=0),第j年种地的亩数S(在
解题时,我们以买卖饲料的数量以及养牛需要的饲料数量来决定),以及贷款的
总额M为三个主要的决策变量。其中第j年卖0岁母牛的数量I,(A=0)会影
响每一年牛的总数进而影响种地的亩数,同时也会影响农场主的贷款额。所以,
本题我们将以第j年牛的总头数,年收入W,年成本C作为对象,并以5年
的总收益的最大值Z作为目标函数进行讨论。为得到五年的净收益总和,我们
将用每年的收入W减去每年的成本C得到每年的净收益并求和来得到。对于成
本中牛的数量超过130时的额外投资费用,我们用年初出生的小母牛的头数减去
年初卖掉的小母牛的头数乘以2000来得到;在决定贷款额时,我们首先对每年
除去还款额后的净收益进行粗略计算得到每一年都是可以盈利的,所以我们将使
用第一年的成本C作为贷款总额。最后,运用ingo软件对决策变量进行规划得
到的结果。我们通过求解得到5年的总收益的最大值2=6866256,贷款总额
M=4160558,种粮食亩数s=80,种甜菜亩数S=120,各年卖出的0岁母牛头
数分别为48,0,15,60,90。其他因素的变化对计划造成的影响,其中银行利率r的
变化会对总收益M造成影响,但对其他决策变量影响较小;另外,如果农产品
价格和产量以及劳动力价格发生变化,将有可能改变种农产品的亩数和购买农产
品的数量,贷款额M和相应的总收益Z的变化,各因素的具体影响我们方式将
在模型解答中加以阐述。
【关键词】优化模型贷款总额M卖掉的0岁母牛的头数买、卖的饲料吨
数种植饲料的土地面积ingo软件
一、问题重述
某公司计划承包有200亩土地的农场，建立奶牛场。开始承包时农场有120头母牛，其中20头为不到0~1岁的幼牛，100头为产奶牛。产奶牛平均每头每年生0.55头公牛，生出后不久即以每头300元卖掉产奶牛平均每头每年生0.05头的母牛要么出生后以400元卖掉，要么饲养，养至2岁成为产奶牛。幼牛年损失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小时；种一亩粮食每年需20小时；种一亩甜菜每年需30小时。
其它费用：每头幼牛每年500元，产奶牛每头每年1000元；种粮食每亩每年150元，种甜菜每亩每年100元。劳动力成本为每小时费用为10元承包农场需要一笔费用，其中租地费用为每年6万元（每年底付清），以及开始承包时农场已有的120头牛的费用。平均产奶牛每头4000元，小牛每头400元，到承包结束时，农场的牛按此价折价抵卖。
该公司资金全部来自银行贷款，复利计息，每年归还本利和的1/5。五年还清。此外，农场主不希望产奶牛的数目在五年末与现在相比减少超过50%，也不希望增加超过75%。
现在要求分析承包人有无盈利的可能性。若有，应如何安排5年的生产，使得五年的净收益为最大，并以此为基础进一步讨论，若遇到银行利率波动（例如上下波动2个百分点），还贷方式改变（如规定每年还息，改变还本的方式），由于气候等外因变化引起的农产品产量与价格的变化及劳动力市场价格的变动等将会对农场主的五年生产计划及收益产生怎样的影响。
二、问题分析
对于问题一：本题要求分析农场主5年后能否盈利，并制定5年的生产计划设法使净收益最大，我们认为这是一个研究资源的合理配置和优化问题。由题意可知，农场的收入来源有卖牛奶的收入，卖牛的收入，以及卖出多余的粮食的收入，同时成本来源包括雇佣劳动力产生的费用，地租，养牛成本（养牛总数高出130头时追加的投资费，饲料费，照顾费等），还贷款以及其他一些费用。由于题中所给的资料均以年为单位，因此，可以通过得出年收益的通项并求和就可得到5年的总收益。经分析，我们可以知道，要想规划每年的生产计划就必须要对每年卖牛的数量，种地的亩数（在解题时，我们以买卖饲料的数量以及养牛需要的饲料数量来决定），以及贷款的总额做出决策。由此可见，其中最为总要的是卖牛的数量，因为它会影响一年牛的总数进而影响种地的亩数，同时也会影响农场主的贷款额。本题我们将以牛的头数以及总收益作为主要的建模对象，来对问题进行分析和阐述。
对于问题二：银行利率等贷款方式的波动会在一定程度上影响贷款的额度，如果实际利率高，将考虑精简贷款额度，反之，可能适当提高贷款额度以获得更高收益农产品价格和产量的变化可能导致在考虑种植不同种饲料（粮食和甜菜）的亩数的决策上有所影响；最后，劳动力市场的价格变化可能对不同的牛（12岁母牛，0岁母牛）的养殖数目，种植饲料的亩数，以及间接地对贷款额造成一定影响。此题我们将结合模型中利率，农产品产量和价格，及劳动力价格这几个变量来进行分析。
三、模型假设
1、幼牛出生发生在年初，损失发生在年末。
2、刚出生的小牛无损失，即刚出生就卖掉，12岁的老母牛在年末被卖掉。
3、幼牛损失的各年龄和奶牛损失的各年龄是均匀的。
4、银行按照复利计息，即农场每年归还1/5（1+r）~5M
5、相邻两个年龄组的牛在相邻两年之间的变化是连续的。
6、在5年末12岁的牛被以4000的价格收购。
7、粮食是在每年末才能成熟并收割，即第一年的牛的全部饲料都从市场上购买，以后每年的牛吃上一年种的粮食，不够的再到市场上购买。
8、年底才能拿到所有的收入，而支出发生在年初。
9、不会卖一岁的母牛
四、符号说明
1、w第；年总收益，W，卖奶收益，Wj卖牛收益，w卖粮食收益2、N；是第j年初母牛的总数量，N；是第j年初第i岁母牛的数量，F1；第j年初成年母牛总数，F，，第；年初未成年母牛总数3、N（A=12）卖出12岁母牛总数，N（Am-0）卖出公牛总数，N（A=0）卖出幼年母牛总数
4、C劳动力花费，C，养牛花费，R地租，C，贷款本利和，C，为第；年牛的头数超过130头时的额外投资，C.为第；年末总花费5、P1为劳动力单位时间费用，P1为粮食卖出价格，P2为粮食买入价格，P1为甜菜卖出价格，P。为甜菜买入价格
6、W，是卖粮食重量，W，是卖甜菜的重量，Wai是购买粮食的重量，W，是购买甜菜的重量
7、S1是第；年种植粮食的面积S；是第；年种植甜菜的面积8、O1为每亩粮食产量，o，为每亩甜菜产量
9、w；是第；年需要的粮食吨数，W；是第；年需要的甜菜吨数10、w；第五年当年净收益，z为5年的最大净收益11、M为贷款总额，利率r
五、模型的建立
4
（一）目标函数；Z=max2（W，-C，）+W3，i=1
M=C1
（二）约束条件：50≤>N*98%+N6*95%≤175
Nay≥0N2≥950≤S；≤800≤S&；0≤82，+S；≤2002%-C）+Al，>=0
H我们将贷款总额M，卖掉的0岁母牛的头数N（A=0），卖掉的粮食吨数Wr，卖掉的甜菜吨数Ws，买入的粮食吨数W，买入的甜菜吨数w。，六个变量作为决策变量，并用（W2u+W.-Wm）/o，表示种植粮食的土地面积s1，用（W，+W4j1-W，）/O，表示种植甜菜的土地面积S2，则
（1）1-4年末收益：卖奶收益W一F1*3700卖牛收益Wcj=N；（A=12）*1200+N；（Am=0）*300+N（A=0）*400卖粮食和甜菜的收益 W。=W，*P，+W，*P1总收益W=W计W，+W于F*3700+N（A=12）*1200+N，（Am-0）*300+N，（Ar
=0）*400+W，*750+W，*5001-4年未花费：劳动力花费C=（F1*42+F2*10+S1*20+S2*30）*P1养牛花费C、=F1；*1000+F2，*500+WmP2+W，p。
地租R=60000+S1；*150+S，*100贷款本利和C；=1/5[M*（1+r）5]r=12%（年利率）超过130头年的额外投资C=Q2N，*055-N，（4=00）*20002=<js=5第一年超过130头牛的额外投资
Ca=CNa*055-N（4=00+120-130*2000豆
总花费C=C1+C+R+Cu+Cj
注：W3=（F1j+F2；*2/3）*0.6，Wj=（F1；+F2j*2/3）*0.7）
（2）第五年净收益AW，=200*1O，*500-000+30*）]
+（2N5*98%+M5*95%）*400022N4*98%*11/2*300+4000+2N6*3700-2N
*1000-N5*500-60000一。*P2-W。*E2-1/TM*I+r）5]
w2.“%+N“95%
NMo=20：*11/2*98%-N，（4=00
NM=NMy:*95%
（2-11岁成年母牛总数）凡，=2N，=Ny+N
（12岁母牛总数）N（A=12）=N1j1*98%
N（Am-0）=2N1*11/2*98%N（Am-0）是卖的0岁公牛总数F1u=100，F1=20+，W1=0，W1=0，W-68，W1=80，r=12%，P=10，01=0.9，0，=1.5，P1=900，P2=750，P1=700，P2=500，N1（Am-0）=55，N1（A=12）=9.8
六、模型求解
（一）问题一的求解
为简化模型，我们增设以下假设：
假设一：前四年每一年都种80亩粮食，120亩甜菜，第五年种200亩甜菜。理由我们从运行结果上可以看到第二年的牛的总头数已经超过了130，经计算，当农场里恰好有100头成年母牛，30头未成年母牛时，农场中80亩地种植的粮食恰好够吃，且种粮食比买粮食更便宜，并且种植甜菜和粮食都能够获得收益。所以，在用程序进行计算时，我们对种植粮食做了规划，200亩地全部种完，由于粮食最多只能种80亩，所以，前4年80亩地种粮食，120亩地种甜菜；由于第五年吃的粮食来源于第四年的收成，而第5年的粮食到年底才有收成，由于甜菜每亩的收益比粮食大，所以第五年200亩的土地全部中甜菜。
假设二：不购买甜菜
理由：经计算，甜菜足够支撑起在题目范围内的母牛的增长率假设三：第一年的总成本C1等于贷款总额M理由：因为经过估算，几乎每一年都可以得到正的净收益，并且，由于第一年的资金比较充裕至少会大于第一年的支出，因为年底才会有收成，所以，第一年的净收益就等于第一年的收入W1，会比较大足矣垫付之后每一年的亏损，所以，为达到支出最小化，我们只借相当于第一年的总成本额度的贷款。
假设四：小母牛在年初产下
理由：由于年初生牛和年末生牛的计算结果相差不大，为方便模型计算，设置在年初生牛。
将已知数代人ingo，运行结果如下：
Global optimal solution found.
objective bound：751029.9
Infeasibilities：0.7331380E-10
Extended solver steps：0
Total solver iterations：10
绘制成表格如下：
|年份|年初产|年初卖|年初卖出|年初|卖出的|年初卖牛后|年末支出的|当年一年
奶牛头|出的幼|后剩余的|牛总|12岁母|超过130头|费用|的年末净
数|年母牛|幼母牛数|数|牛头数|的牛的数目||利润
|第一|100|1560160|9.8|30519798.7|2139.106
|第二|98 |48|54|152 |9.6|22|325892.8 |126925.6
|第三|132|1|77 |209 |9.4|79 |547304.4 |348.9236
|第四|125 680125|9.2|5|316200 397370.5
|第五|178 |98|0|178 |0|48|424080|330750
|年
即可知：
第五年末盈利Z=2139.106+126925.6+348.9236+397370.5+330750=751029.9元
其中5年的生产计划如下
贷款总额M=519798.7，前四年种粮食亩数S1=80，种甜菜亩数S2=120，第五年
种植200亩甜菜，各年卖出的0岁母牛头数分别为Nj（A-0）=15，Nj（A-0）
=48，Nj（A=0）=1，Nj（Af=0）=68，Nj（Af=0）=98.
由上式可以看到，对于银行利率r上升会导致年贷款本利和C1，的上升，进而导致年成本C的上升，从而引起总收益Z的降低，但由于贷款总额M是由第一年的成本C决定，与银行利率没有关系，所以贷款总额M将不会变化，并且银行利率与其他决策变量没有直接或间接地关系，所以不会对其造成影响；如果其他因素不变，当甜菜的亩产o，低于0.56时应不种值，当0.56=<0，<=0.8时，应当种植，当o，>0.8时，才能从种植甜菜上盈利；同理，当粮食的亩产O，<=0.39时不宜种植，当0.39=<O，<=0.46时，前四年应种满80亩，当O1>0.46时，才能从种植粮食上盈利。另外，当亩产降低时将会增加成本，从而增加贷款额度M，减少总收益Z。
其他因素不变，当粮食的卖出价格P>388.9时，且此时若粮食买入价格Pos=P1，前四年应该种满80亩，并使牛吃的所有都从市场上购买，当粮食的买入价格，P1<P。时，牛先吃自己种的粮食，不够的从市场上买；当粮食的卖出价格P1<388.9时，且此时若粮食买入价格Po<=P1，应不种，并使牛吃的所有都从市场上购买，当粮食的买入价格，Pi<P。时，牛先吃自己种的粮食，不够的从市场上买。同理，266.7为甜菜的买卖价格的临界点。此外，如果价格变动带来粮食收入的升高，但会将减少贷款额度M，增加总收益z，卖奶牛的计划不会发生改变；反之，则会增加贷款额度M，减少总收益z，卖奶牛的计划不会有改变。
如果其他因素不变，劳动力单位时间费用增加则成本增加，总收益z减少，贷款总额M增加；当P1<31.7时，前四年种满80亩粮食，种120亩甜菜，并先吃自己种的，不够再买；当31.7<P<33时，不种甜菜只种粮食，种80亩粮食，先吃自己种的，不够再买；当p>33时，不种农作物，所有饲料都从市场买。同时卖奶牛的头数不发生变化。
七、模型评价
（一）优点：本模型适用广泛，对其他一些生物资源的开发具有同等可靠性，特别是对于线性的优化问题，ingo软件配合目标优化模型的准确性较高。此外，本模型在解决每一个小问题时，根据题设对模型做了额外的假设，简化了解题步骤，并对程序和模型本身进行了二次简化，使结果更加简单易懂，清晰明了。模型中我们建立了多个决策变量，我们尽量使用了比较直观的符号来表示变量名称，并简化了变量的表达式，使得每个变量之间的关系一目了然。而且，本题主要以讨论牛的头数为出发点，层层深入，逻辑性更显严密。
（二）缺点：在模型的优化前期，变量个数相对较多，需要耗费比较大的时间来分析变量之间的关系并对模型进行简化；由于结构非常的严密，逻辑性很强，变量很多，所以，读者需要密切注意每一个变量。
八、模型推广
以本模型主要解决如何达到制定的长远的目标优化问题，以此为理论基础，适当改变某些参数就能更多地推广应用到多种经济问题中、工作计划，科研，教育等多个领域及其它饲养场中的一些简单优化问题；如果对模型进一步改进，引入非线性的变量可能能够解决一些复杂的优化问题。例如，如何在边际效益递减的情况下，获得效益的最大化，又或者是诸如国家经济的发展计划的一些变量较复杂，变量间关系较为复杂，呈现非正规线性关系的优化问题。
参考文献
[1]姜启源谢金星叶俊，《数学模型》（第四版），北京高等教育出版社，2011年；
[2]韩中庚，《数学建模方法及其应用》北京市：高等教育出版社，2009年；
[3]赵可培，《运筹学》，上海市：上海财经大学出版社，2000年；
[4]吴祈宗，《运筹学与最优化方法》北京市：机械工业出版社，2003年；
[5]谢金星薛毅，《优化建模与LINDO/LINGO软件》，北京：清华大学出版社，
2005年。
附录
Lingo程序如下：
mode1：
max=z1+z2+z3+z4+z5-100*4000-20*400+x5*4000*0.98-（（1+0.12）
入5-1）*s1；
z1=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-（×1*0.55-r1-10）*2000-
60000；
z2=3700*x2+×2*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-×2*42*10-x2*1000-（×2*0.55-
r2）*2000-60000；
z3=3700*×3+x3*0.55*300+r3*400+10*0.98*0.98*0.98*1200-（0.6
*×3+0.4*y3-72）*900-20*80*10-150*80-120*30*10-120*100+（180
-0.7*×3-0.7*2/3*y3）*500-y3*100-y3*500-×3*42*10-x3*1000-（x
3*0.55-r3）*2000-60000；
z4=3700*x4+x1*0.55*300+r4*400+10*0.98*0.98*0.98*0.98*1200
-（0.6*×4+0.4*y4-72）*900-20*80*10-150*80-120*30*10-120*100
+（180-0.7*×4-0.7*2/3*y4）*500-y4*100-y4*500-x4*42*10-x4*10
00-60000+200*1.5*500；
z5=3700*X5+x5*0.55*300+r5*400-（0.6*X5+0.4*y5-72）*900-200*
30*10-200*100+（180-0.7*x5-0.7*2/3*y5）*500-y5*100-y5*500-x
5*42*10-85*1000-60000；
s1=（0.6*x1+0.4*y1）*900+20*80*10+150*80+120*30*10+120*100+
y1*100+y1*500+x1*42*10+x1*1000+（x1*0.55-r1-10）*2000+60000
+（0.7*x1+0.7*2/3*y1）*700；
s2=（0.6*×2+0.4*y2-72）*900+20*80*10+150*80+120*30*10+120*1
00+y2*100+y2*500+×2*42*10+×2*1000+（×2*0.55-r2）*2000+60000
s3=（0.6*×3+0.4*y3-72）*900+20*80*10+150*80+120*30*10+120*1
00+y3*100+y3*500+×3*42*10+×3*1000+（×3*0.55-r3）*2000+60000；
s4=(0.6*×4+0.4*y4-72)*900+20*80*10+150*80+120*30*10+120*1
00+y4*100+y4*500+x4*42*10+x4*1000+60000;
s5=(0.6*×5+0.4*y5-72)*900+200*30*10+200*100+y5*100+y5*500
+x5*42*10+x5*1000+60000;
x1=100;
×2=x1*0.98;
X3=×2*0.98-10*0.98*0.98+10*0.95*0.95+(x1*0.55-r1)*0.95*0.
95;
x4=×3*0.98-10*0.98*0.98*0.98+(X2*0.55-r2)*0.95*0.95;
x5=×4*0.98-10*0.98*0.98*0.98*0.98+(×3*0.55-r3)*0.95*0.95;
X6=X5*0.98+(×4*0.55-r4)*0.95*0.95;
r1>=0;
r2>=0;
r3>=0;
r1<=x1*0.55;
r2<=X2*0.55;
r3<=×3*0.55;
r4<=×4*0.55;
r5=X5*0.55;
y1=x1*0.55-r1+20;
y2=(x1*0.55-r1)*0.95+10*0.95+×2*0.55-r2;
y3=(×2*0.55-r2)*0.95+×3*0.55-r3;
y4=（x3*0.55-r3）*0.95-r4；
y5=0；
X6>=50；
X6<=175；
@gin（×1）；@gin（x2）；0gin（x3）；@gin（×4）；@gin（x5）；
end
程序中的变量名一览：
X1第一年初产奶牛的总数y1第一年初剩余小母牛的总数
X2第二年初产奶牛的总数y2第二年初剩余小母牛的总数
X3第三年初产奶牛的总数
y3第三年初剩余小母牛的总数
X4第四年初产奶牛的总数y4第四年初剩余小母牛的总数
X5第五年初产奶牛的总数y5第五年初剩余小母牛的总数
X6第五年未产奶牛的总数
Z1第一年当年的净收益s1第一年当年的成本
Z2第二年当年的净收益
s2第二年当年的成本
Z3第三年当年的净收益
s3第三年当年的成本
24第四年当年的净收益s4第四年当年的成本
Z5第五年当年的净收益s5第年年当年的成本
r1第一年初卖掉的小母牛的总数
R2第二年初卖掉的小母牛的总数
R3第三年初卖掉的小母牛的总数
R4第四年初卖掉的小母牛的总数
R5第五年初卖掉的小母牛的总数
网友评论