数学建模—农作物施肥的优化设计.doc

大学生数学建模题目：施肥效果分析学院 电气工程学院班级 组号 姓名 姓名 姓名 姓名 姓名 农作物施肥得优化设计摘要本文在合理得假设之下，通过对实验数据得分析，建立了能够反映施肥量与农作物产量得关系模型，据此求得在保证一定产量得同时,施用肥料最少。首先就是对实验数据进行了较为直观得分析，可知肥、肥、K肥施加不同量均对土豆、生菜得产量造成一定影响，且施N肥过多会烧苗，会使土豆与生菜减产。其次,模型一，我们对实验数据运用el进行拟合,得到各肥料得施肥量与产量得拟合曲线,从而获得对应函数表达式。但由于无法对模型进行误差分析，我们再次运用一元多项式回归方法建立模型进行求解，此时得到不同肥料得施肥量与产量得关系。然后,模型二，利用atlb软件建立模型，求出肥、P肥、肥得施肥量关于土豆及生菜得最优解:当氮得施肥量为9、252时使得土豆产量达到最优解为3、34615；当磷得施肥量为30时使得土豆产量达到最优解为4、742；当钾得施肥量为6、0742时使得土豆产量达到最优解为4、51718.当氮得施肥量为20、242时使得生菜产量达到最优解为、345；当磷得施肥量为20、252时使得生菜产量达到最优解为43、4615；当钾得施肥量为29、2时使得生菜产量达到最优解为4、3465。最后我们就应用价值方面对模型做出改进。由于实验数据中各个自变量与因变量之间并不就是一一对应得关系，所以没有得出各肥料得施肥量与产量得交叉关系,仅得到单一变量得对应关系。关键字:一元多项式回归 Excl拟合 Matlab一、问题得提出某地区作物生长所需得营养素主要就是氮(N）、钾(K)、磷(P)。某作物研究所在某地区对土豆与生菜做了一定数量得实验,实验数据如下列表所示，其中ha表示公顷,表示吨，kg表示公斤.当一个营养素得施肥量变化时，总将另两个营养素得施肥量保持在第七个水平上,如对土豆产量关于N得施肥量做实验时， P与 K得施肥量分别取为 16kha与37kha。试分析施肥量与产量之间关系，并对所得结果从应用价值与如何改进等方面做出估计。土豆: N P K施肥量(gha)产量(/a）施肥量（kg/a)产量（tha)施肥量(kg/ha）产量(ha)15、80、4018、3421、36243、47、356725、7293、069334、610132、2937、961409、21353、39841、04183、4420239、514740、9797、725943、15164、37238、43363、462454、746543、87404、32944、65582、774730、73242、735146、22生菜： N K施肥量(k/a)产量(t/ha）施肥量（kg/ha)产量（th)施肥量(kga）产量(t/a)01、00、39015、5281、99、4846、76561、56981、469316、898416、27174、33416、2111、751617、1018617、56162、92942、219、202241、63912、637217、728、448921、3455、8366、125822、07520、1139214、11652、53519、40二、问题得分析利用散点图对所拟合问题得曲线类型做出判断。当需要拟合得两变量之间得函数关系式，首先要确定所求函数对应曲线得类型,然后根据曲线类型对所求函数得对应关系进行假设,并利用已知数据计算出所需参数，最终确定变量之间得函数关系.我们可以分别绘制出土豆与生菜得产量与施肥量得散点图,从图像得角度判断函数关系，再根据题目所给数据确定最终得函数.三、问题得假设与符号说明2、1模型得合理假设（1） 土壤本身已含有一定数量得氮、磷、钾肥,即具有一定得天然肥力。（2） 每次实验就是独立进行得，互不影响。(）研究所得实验就是在相同得实验条件(实验结果不受温度，水,光照等因素影响)下进行得,产量得变化就是由施肥量得变化引起得。(）当一个营养素得施肥量变化时，另两个营养素得施肥量总保持在第七水平上不变。()所给数据中无较大偏差点,无需剔除。2、2模型得符号说明:土豆产量 :对于土豆氮得施肥量:对于土豆磷得施肥量:对于土豆钾得施肥量:生菜产量：对于生菜氮得施肥量:对于生菜磷得施肥量：对于生菜钾得施肥量四、模型得建立与求解土豆得产量与施肥量得散点图如下：土豆产量-氮施肥量散点图土豆产量-磷施肥量散点图土豆产量钾施肥量散点图生菜得产量与施肥量得散点图如下：生菜产量氮施肥量散点图生菜产量磷施肥量散点图生菜产量钾施肥量散点图所用malab程序为:1=lsred( ：数学建模课程设计实验报告shuu，sheet1，：$L12）；y31=xlsread( E：数学建模课程设计实验报告huju,shet，M3:12）;plt(k1,y31，+)由散点图猜测土豆产量与氮施肥量得关系式为:;与磷施肥量得关系式为:；与钾施肥量得关系式为：由matb解出： a0、00 b1=0、7 c1=14、746 =-、001 b2=0、7 2=32、61 a=42、 30、56 c3=0、0土豆产量与施肥量得关系图：最佳施肥方案为第一个方案(328、，245，65）所用程序为：cearclc1=0、003; b10、97； c1=14、2；a2=-0、001; b2=、071； c2=32、916;a3=42、7； b30、6； c、01；n=0：0、01：393;p=0：0、01：686；k=0：、01：652;y1（a1n、n+b*n+1）*80;y1=ma(y1)oi=1:legh（n） if abs(y（i)y11)=0、01 q1n(） brea edendy2=（a2、*n+b2nc)*80;y22mx（y2）or 1:egh（p) f as(y2(i)22）=0、00 q2=p（i） e endedy=3*（1-b3ep（-c3*k)）；y3=max(y)fo i=1:ngh(k） if as（3(）y33)=0、001 q=k(i） brk ede运算后得结果如下: 生菜产量与施肥量关系：由散点图猜测生菜产量与施肥量得关系式为：与磷肥得量得函数为:与钾肥得量得函数为： 由atab解出：a1=、0002 1=0、101 c=0、229a=0、0001 =0、066 =、877a3=15、87 b-0、40 3=0、0026关系图为: 最佳施肥方案为第一个方案（23、1，245,465）所用程序为：clerlc1=0、0002； b=0、101； c110、2294;a=0、000; b=、060； c2=6、8757;a3=1、8878; b3=-、0440； c=0、026;n=0：0、0：93；=0:0、1：66;=：、01:62；y=（a1*n、*n+b1*n1)*0;y1=m（y)for i1：egth（n） if ab（y1（）11)=0、0 q1=（i) break enddy2=（a2、*nb2n+c2）00;y22=max（y2）for =:lengh（) fas（y（i)-y2）、001 2=p（) break edend3=a(-bxp（c*k);3=max(3)or 1：ength(） if bs(y3（i）y33)=0、001 q3k（) brak ende运行结果如图：五、模型得优缺点与改进方向5、1模型得优点（1)本模型运用回归分析得方法求解，理论可得最优解.（)模型就是独立得模型进行逐步回归.（3）利用Matab编程,曲线估计较成功地解决了施肥最佳方案问题，方法简练,道理清晰，结果可信。曲线估计得到较合适得曲线,最终得到拟合曲线函数表达式。5、模型得缺点在实际工作中，三种肥料之间除了与产量有直接得数量关系外，还有彼此之间得交互作用.因此，本模型只就是一个初步得探讨，要得到三种营养素与产量之间得准确关系，应该在实验之初就采取正交实验或均匀设计得方法，得到更有价值得实验数据，从而更好得把握变量间得数量关系,以达到直到农业生产实践得目得。5、3模型得改进该模型只对N、P、得施肥量与产量进行了分析，还应考虑N、P、K得肥料售价与土豆、生菜每吨得售价，从而获得更高得收益。根据实际情况,当施肥料带来得收入比用于购买肥料得费用多时，应该多施肥，否则少施肥。参考文献1 熊卫国、数学实验教程、广东：中山大学出版社、2002 李玉莉、MATL函数速查手册M、北京：化学工业出版社、203 姜启源谢金星、数学模型M、北京、高等教育出版社、2010 马莉、Malab数学建模与实验、清华出版社5冯杰、数学建模原理与案例、科学出版6 董臻圃、数学建模方法与实践、国防工业出版社