A 题施肥效果分析第三组摘要本文就施肥量与农作物产量的关系进行研究分析，运用逐步回归的思想建立数学模型，最初方案从比较简单的模型入手，逐步优化最终得到各营养素施肥量的最优配比，并使其应用价值得到推广。
思路：考虑到多种肥料对产量的影响复杂且不易得到其关系，因此我们采用控制单一变量的方法，做出散点图并进行逐步回归，得到产量与各个单一变量的关系式，在此基础上将多个变量做多元回归最终得到施肥量与产量的函数关系。
综合考虑各个方案选择最佳方案作为最终模型，并加以推广。
最终的最优模型为：当氮的施肥量为290.2542、磷的施肥量为303、钾的施肥量为536.0742时土豆产量达到最优解；当氮的施肥量为290.2542、磷的施肥量为290.2542、钾的施肥量为290.2542时生菜产量达到最优解。
我们通过运用excel、MATLAB、SPSS等软件做出散点图并进行曲线的拟合，用LINGO 软件规划求得给定目标函数在限制条件下的最优解,用SPSS进行一元回归。
综合以上所给出的最终模型，各营养素施肥量的最优配比对农作物产量的提高有很好的应用价值。
模型改进：因为获得的实际数据较少，使模型的精确度受到影响，采用的数学模型因此不够精准，改进建议是收集更多实际数据，统计分析，改进模型。
关键词：散点图逐步回归目标函数目录一、问题重述................................................................................. 错误！未定义书签。
二、符号说明................................................................................. 错误！未定义书签。
三、模型假设................................................................................. 错误！未定义书签。
四、问题分析................................................................................. 错误！未定义书签。
五、模型的建立与求解................................................................. 错误！未定义书签。
六、模型的改进与评价................................................................. 错误！未定义书签。
七、参考文献................................................................................. 错误！未定义书签。
八、附录......................................................................................... 错误！未定义书签。
一、问题重述某地区作物生长所需的营养素主要是氮N、钾K、磷P。
某研究所在该地区对土豆与生菜做了一定数量的实验，实验数据如下表所示，其中ha表示公顷，t表示吨，kg表示公斤。
当一个营养素的施肥量变化时，P与K的施肥量分别取为ha196与kg/ kg/372.ha试分析施肥量与产量之间的关系，并对所得结果从应用价值与如何改进等方面做出估价。
土豆：生菜：二、模型假设1.农作物实验是在相同的实验条件下进行，产量的变化仅由施肥量引起；2.土壤中有一定量的天然肥，满足其在不施肥下也能生长；3.当一个营养素的施肥量变化时，另两个营养素的施肥量总保持在第七水平上不变；4.每次实验独立进行，互不影响；5.剔除所给数据中偏差较大的点后进行曲线拟合。
三、符号说明x：每公顷土地中氮元素的施肥量1x：每公顷土地中磷元素的施肥量2x：每公顷土地中钾元素的施肥量3y：每公顷土地中施加氮后作物的产量1y：每公顷土地中施加磷后作物的产量2y：每公顷土地中施加钾后作物的产量3四、问题分析这是一个由多个变量控制的数学模型，以各营养素施肥量配比使产量达到较好状态为目标，因此需首先求得各个变量的最优解从而使产量达到最佳状态。
对题目中给出的数据进行处理，找出产量与单一营养素在其他营养素施肥量不变条件下的数据关系，绘出土豆和生菜的产量与施肥量的散点图观察其分布规律，选用适当的拟合曲线，并根据曲线类型得出函数关系式，运用SPSS对所得函数进行回归分析和假设检验最终得出各营养素的最终配比。
五、模型的建立与求解方案一：用excel表格做出施肥量与产量的拟合曲线如下图所示：1.施肥量对土豆产量的影响：氮元素：磷元素：钾元素：2.施肥量对生菜产量的影响：氮元素：磷元素：钾元素：综上第一套方案对曲线的拟合程度来看，虽然拟合程度很高，但是运用excel做出的拟合曲线并不能反映出施肥量对农作物产量的影响总体规律，只能是在给定的数据范围内达到拟合的要求，并不能准确的把握施肥量与产量之间的关系，故不具有应用价值，不可采用。
方案二：用SPSS做施肥量与产量的拟合曲线：1.施肥量对土豆产量的影响：氮元素：氮元素施肥量与土豆产量的关系式是：742.14197.02+=xy用SPSS对模型做回归分析所得数据：2R ：0.986，估计值的标准误差：1.272 所得模型符合总体要求。
磷元素：磷元素施肥量与土豆产量的关系式：005.3210493.6114.073+⨯+=-x x y 用SPSS 对模型做回归分析所得数据：2R ：0.898，估计值的标准误差：1.424 所得模型符合总体要求。
钾元素：钾元素施肥量与土豆产量的关系式：549.2010966.3168.072+⨯+=-x x y 用SPSS 对模型做回归分析所得数据：2R ：0.933，估计值的标准误差：2.578 所得模型符合总体要求。
2.施肥量对生菜产量的影响： 氮元素：氮元素施肥量与生菜产量的函数关系式：229.10101.0+=x y 用SPSS 对模型做回归分析所得数据：2R ：0.925，估计值的标准误差：1.167 所得模型符合总体要求。
磷元素：磷元素施肥量肥与生菜产量的函数关系式：685.510056.188.073+⨯+=-x x y 用SPSS 对模型做回归分析所得数据：2R ：0.978，估计值的标准误差：1.149 所得模型符合总体要求。
钾元素：钾与生菜产量的函数关系式：69.1510573.51047.5018.0853+⨯+⨯-=--x x x y 用SPSS 对模型做回归分析所得数据：2R ：0.515，估计值的标准误差：1.328 所得模型有误差。
综上对方案二的研究分析可以看出虽然SPSS 对曲线的拟合度较高，但是其所拟合的曲线并不符合实际规律，因为根据实际经验，肥料在一定量下并不是施的越多，作物产量越大，而是到超过一定量时，会导致土壤中溶液浓度过大而发生“烧苗”现象，经过综合实际情况与所得结果将施肥量与产量的关系定性为二次函数关系式，由此得出改进方案。
改进方案：SPSS 拟合的二次曲线对自变量系数偏差较大，在回归分析时误差较大，因此我们采用软件MATLAB 求解得到二次函数关系方程式，所得方程式用LINGO 软件进行分析并求解其最优解。
土豆产量与各个营养素施肥量的关系图：氮与土豆产量函数关系式7416.141971.0103953.324++⨯-=-x x y磷与土豆产量函数关系式7375.320659.0100825.124++⨯-=-x x y 钾与土豆产量函数关系式4144.240750.0109953.625++⨯-=-x x y用LINGO 软件求所得函数关系式的最优解：当氮的施肥量为290.2542时使得土豆产量达到最优解为43.34615； 当磷的施肥量为303时使得土豆产量达到最优解为42.7423；当钾的施肥量为536.0742时使得土豆产量达到最优解为44.51718.生菜产量与各个营养素施肥量的关系图：氮与生菜产量函数关系式2294.101013.0103815.224++⨯-=-x x y磷与生菜产量函数关系式8820.60606.0104577.525++⨯-=-x x y钾与生菜产量函数关系式8362.150106.0100768.726++⨯-=-x x y用LINGO软件求所得函数关系式的最优解：当氮的施肥量为290.2542时使得土豆产量达到最优解为43.34615当磷的施肥量为290.2542时使得土豆产量达到最优解为43.34615当钾的施肥量为290.2542时使得产量达到最优解为43.34615六、模型结果的改进与评价优点：1）本模型运用回归分析的方法求解，理论可得最优解。
2）问题的考虑是逐层深入，使得问题的各个因素考虑比较全面。
3）模型是分为两个独立的模型进行逐步回归。
4）我们运用的 MATLAB、SPSS等软件使得结果更加精确，误差主要取决于Lingo 的软件的精度。
缺点：1）在实际工作中，三种肥料之间除了与产量有直接的关系外，还有彼此之间的交互作用。
因此，本模型只是初步探讨，要得到三种营养素之间与产量的准确关系应在实验之初采用正交实验与均匀设计的方法，得到更有价值的实验数据，建立更有可信度的自变量与因变量之间的函数关系式。
2）整个方案中对于问题考虑的不全面，导致最终结果存在误差。
七、参考文献[1]崔群法、张瑜、杨光霞、李乃文.《数据建模与应用》.清华大学出版社[2]张举刚、张彩霞、李荣平、张小平.《统计学》.河北人民出版社[3]马昌风.《最优化方法及其MATLAB程序设计》.科学出版社[4]蔡建琼、于惠芳.《SPSS统计分析实例精选》.清华大学出版社八、附录x=[0 34 67 101 135 202 259 336 404 471];y=[15.18 21.36 25.72 32.29 34.03 39.45 43.15 43.46 40.83 30.75];subplot(2,3,1);plot(x,y,'o'),xlabel('N'),ylabel('CL'),title('N施肥量与土豆产量的关系');p=polyfit(x,y,2);m=p(1);b=p(2);c=p(3);x1=[0:5:500];y1=m*x1.^2+b*x1+c;subplot(2,3,4);plot(x1,y1,x,y,'o'),xlabel('N'),ylabel('CL'),x=[0 24 49 73 147 196 245 342];y=[33.46 32.47 36.06 37.96 40.09 41.26 42.17 42.73];subplot(2,3,2);plot(x,y,'o'),xlabel('P'),ylabel('CL'),title('P施肥量与土豆产量的关系');p=polyfit(x,y,2);m=p(1);b=p(2);a=p(3);x1=[0:5:500];y1=m*x1.^2+b*x1+a;subplot(2,3,5);plot(x1,y1,x,y,'o'),xlabel('N'),ylabel('CL'),x=[0 47 93 140 186 279 372 465 558 651];y=[18.98 27.35 34.86 38.52 38.44 37.73 38.43 43.87 42.77 46.22];subplot(2,3,3);plot(x,y,'o'),xlabel('K'),ylabel('CL'),title('K施肥量与土豆产量的关系') p=polyfit(x,y,2);m=p(1);b=p(2);a=p(3);x1=[0:5:700];y1=m*x1.^2+b*x1+a;subplot(2,3,6);plot(x1,y1,x,y,'o'),xlabel('K'),ylabel('CL'),x=[0 28 56 84 112 168 224 280 336 392];y=[11.02 12.70 14.56 16.27 17.75 22.59 21.63 19.34 16.12 14.11];subplot(2,3,1);plot(x,y,'o'),xlabel('N'),ylabel('CL'),title('N施肥量与生菜产量的关系');p=polyfit(x,y,2);m=p(1);b=p(2);a=p(3);x1=[0:5:500];y1=m*x1.^2+b*x1+a;subplot(2,3,4);plot(x1,y1,x,y,'o'),xlabel('N'),ylabel('CL'),x=[0 49 98 147 196 294 391 489 587 685];y=[6.39 9.48 12.46 14.38 17.10 21.94 22.64 21.34 22.07 24.53];subplot(2,3,2);plot(x,y,'o'),xlabel('P'),ylabel('CL'),title('P施肥量与生菜产量的关系');p=polyfit(x,y,2);m=p(1);b=p(2);a=p(3);x1=[0:5:700];y1=m*x1.^2+b*x1+a;subplot(2,3,5);plot(x1,y1,x,y,'o'),xlabel('P'),ylabel('CL'),x=[0 47 93 140 186 279 372 558 651];y=[15.75 16.76 16.89 16.24 17.56 19.20 17.97 20.11 19.40];subplot(2,3,3);plot(x,y,'o'),xlabel('K'),ylabel('CL'),title('K施肥量与生菜产量的关系');p=polyfit(x,y,2);m=p(1);b=p(2);a=p(3);x1=[0:5:700];y1=m*x1.^2+b*x1+a;subplot(2,3,6);plot(x1,y1,x,y,'o'),xlabel('K'),ylabel('CL'),x=[0 34 67 101 135 202 259 336 404 471];y=[15.18 21.36 25.72 32.29 34.03 39.45 43.15 43.46 40.83 30.75];subplot(2,3,1);plot(x,y,'o'),xlabel('N'),ylabel('CL'),title('N施肥量与土豆产量的关系');p=polyfit(x,y,2);m=p(1);b=p(2);c=p(3);x1=[0:5:500];y1=m*x1.^2+b*x1+c;subplot(2,3,4);plot(x1,y1,x,y,'o'),xlabel('N'),ylabel('CL'),x=[0 24 49 73 147 196 245 342];y=[33.46 32.47 36.06 37.96 40.09 41.26 42.17 42.73];subplot(2,3,2);plot(x,y,'o'),xlabel('P'),ylabel('CL'),title('P施肥量与土豆产量的关系');p=polyfit(x,y,2);m=p(1);b=p(2);a=p(3);x1=[0:5:500];y1=m*x1.^2+b*x1+a;subplot(2,3,5);plot(x1,y1,x,y,'o'),xlabel('N'),ylabel('CL'),x=[0 47 93 140 186 279 372 465 558 651];y=[18.98 27.35 34.86 38.52 38.44 37.73 38.43 43.87 42.77 46.22];subplot(2,3,3);plot(x,y,'o'),xlabel('K'),ylabel('CL'),title('K施肥量与土豆产量的关系') p=polyfit(x,y,2);m=p(1);b=p(2);a=p(3);x1=[0:5:700];y1=m*x1.^2+b*x1+a;subplot(2,3,6);plot(x1,y1,x,y,'o'),xlabel('K'),ylabel('CL'),x=[0 28 56 84 112 168 224 280 336 392];y=[11.02 12.70 14.56 16.27 17.75 22.59 21.63 19.34 16.12 14.11];subplot(2,3,1);plot(x,y,'o'),xlabel('N'),ylabel('CL'),title('N施肥量与生菜产量的关系');p=polyfit(x,y,2);m=p(1);b=p(2);a=p(3);x1=[0:5:500];y1=m*x1.^2+b*x1+a;subplot(2,3,4);plot(x1,y1,x,y,'o'),xlabel('N'),ylabel('CL'),x=[0 49 98 147 196 294 391 489 587 685];y=[6.39 9.48 12.46 14.38 17.10 21.94 22.64 21.34 22.07 24.53];subplot(2,3,2);plot(x,y,'o'),xlabel('P'),ylabel('CL'),title('P施肥量与生菜产量的关系');p=polyfit(x,y,2);m=p(1);b=p(2);a=p(3);x1=[0:5:700];y1=m*x1.^2+b*x1+a;subplot(2,3,5);plot(x1,y1,x,y,'o'),xlabel('P'),ylabel('CL'),x=[0 47 93 140 186 279 372 558 651];y=[15.75 16.76 16.89 16.24 17.56 19.20 17.97 20.11 19.40];subplot(2,3,3);plot(x,y,'o'),xlabel('K'),ylabel('CL'),title('K施肥量与生菜产量的关系');p=polyfit(x,y,2);m=p(1);b=p(2);a=p(3);x1=[0:5:700];y1=m*x1.^2+b*x1+a;subplot(2,3,6);plot(x1,y1,x,y,'o'),xlabel('K'),ylabel('CL'),。