农作物施肥效果分析报告

《农作物施肥效果分析报告》由会员分享，可在线阅读，更多相关《农作物施肥效果分析报告（19页珍藏版）》请在金锄头文库上搜索。

1、农作物施肥效果分析摘要我们通过研究氮、磷、钾三种肥料对土豆和生菜的作用，来建立施肥量与产量关系 的模型。通过回归分析的方法，将所给的数据进行MATLAB工具箱拟合，并利用残差分 析的方法，建立反映施肥量与产量关系的模型并检验分析,找到产量的最优解以及氮、 磷、钾三种肥料的最优配合比，在耕地面积一定的情况下研究土豆或生菜可以达到得最 大收益值。由此我们建立的土豆产量模型为生菜产量模型为求解得到土豆产量的最值，当，时，得出，氮磷钾肥料的最优配合比为1.5:1:1.74，土豆是喜钾作物。我们可以得出生 菜的最值，当，时，得，可以看出生菜是喜磷作物。在应用方面，为了直观的展示最大的利润以及最优配合比，

2、设计了一个GUI人机交 互界面，这样可以清晰明了表示获得的最大收益值。关键词：回归分析MATLAB拟合残差分析最优配合比GUI人机交互界面一问题重述俗话说“民以食为天”，我们的生活与农作物的供应息息相关。近年来，随 着人口增多，耕地减少，所以化肥对农作物的生长、提高农作物的产量具有重要 的意义。农作物除了吸收水分和空气中二氧化碳以获得碳、氢、氧等元素外，还 必须从土壤再吸收氮、磷、钾和其他矿质养分，并在太阳能的帮助下合成有机物 质，以建造自己的有机机体，但土壤中的常量营养元素氮、磷、钾和其他矿质养 分一般不能满足作物生长的需求，需要施用含氮、磷、钾的化肥来补充。在本问 题中，某研究所通过研究氮

3、、磷、钾三种肥料对土豆和生菜的作用，来建立施肥 量与产量关系的模型。实验中将每种肥料的施用量分为 10个水平，在考察其中一种肥料的施用量 与产量的关系时，把另两种肥料固定在第 7 个水平上，通过回归分析的方法，将 所给的数据进行MATLAB拟合，从而建立反映施肥量与产量关系的模型，找到产量 的最优解以及氮、磷、钾三种肥料的最优配合比，在耕地面积一定的情况下研究 土豆和生菜可以达到的最大收益值，并从实际情况出发，评价该模型的优缺点。 二问题假设及符号假设2.1 假设在不同的实验条件下，保持水分、温度、湿度、光照、土壤状况等外界 条件一致。2.2 假设施加的化肥完全进入土壤，没有挥发作用。2.3

4、假设在模型一、模型二中氮、磷、钾三种肥料对土豆和生菜的影响是相互独立的。2.4 符号假设：氮、磷、钾施用量分别为 、 、 ；土豆的产量为，生菜的产量为；氮肥 每公斤价格为x,磷肥每公斤价格为y,钾肥每公斤价格为z； 土豆每吨利润为a， 生菜每吨利润为b；其他的固定成本为m；总收益为s。三模型的建立和求解3.1 模型一的建立与求解为了大致分析土豆的产量、生菜的产量与 、 、 的关系,首先利用表中 的数据分别做出与 、与 、与 的散点图并进行拟合、建立模型,如下图所示。图一 对 的散点图从图一可以发现，随着 的增加，先增加在减小，由 MATLAB 拟合为二次函 数。其数学模型为5求参数可得；求解最

5、值得，当 时， ；与实际最值得对比发现，模型给出 的最值小于实际给出的最值，这是因为，给出的模型是为了探究一般情况，故取 该曲线的最值点，而实验给出的最值点不在该曲线上，因而舍弃了实验给出的最 值点。图二 对 的散点图同理给出对 的散点图，在拟合过程中我们发现指数函数拟合效果较二次函数拟合效果好，所以给出指数函数数学模型求解函数最值得，当 时 述原因。，最值的误差来源同样是因为上图三 对 的散点图给出对 的散点图，并得出数学模型由图像可得最值， 时， 。下面我们研究生菜的产量与 、 、 的关系，并做出与 、与 和与 的散 点图并进行拟合、建立模型。图四 对 的散点图给出对 的散点图，并得出数学

6、模型求解最值得，当 时， ；最值的误差是因为拟合方程与 实验数据不完全吻合。图五 对 的散点图给出对 的散点图，并得出数学模型求解最值得，当 时， ，最值误差来源仍是上述原因图六 对 的散点图给出对 的散点图，并得出数学模型显然由图六所示，该模型为线性模型，故无最值的存在。3.2 模型二的建立与求解 模型一只是给出了产量与单一变量（氮、磷、钾的三种施用量）之间的关系，我们在此基础上进行改进，建立产量与多变量之间的数学模型。首先先建立土豆产量与三种肥料施用量关系的模型。综合上面的分析，土豆 的产量与氮肥施用量的关系是二次函数模型；土豆的产量与磷肥施用量的关系是 四次函数模型；土豆的产量与钾肥施用

7、量的关系是指数函数模型，为建立三者统 一的回归模型，必须简化计算，将土豆的产量与磷肥施用量的关系、土豆的产量 与钾肥施用量的关系均转化为二次函数模型。结合土豆产量的模型一建立如下回 归模型其中 是当给出氮肥施用量、磷肥施用量、钾肥施用量时， 土豆产量的平均值，其中参数， ， ， ， ， ，称为回归系数，应大致服从 均值为 0 的正态分布。我们利用MATLAB统计工具箱中的命令regress求解，使用格式为b,bint,r,rint,stats=regress（w1,X，0.05）其中输入为上述回归模型中的数据（ 维向量，），为对应于回归系数的（，）的数据矩阵（ 矩阵，其中第一列为全1向量），0

8、.05为置信水平；输出b为的估计值，bint为b的置信区间，r为残差 向量， rint 为 r 的置信区间， stats 为回归模型的检验统计量，有4个值，第1 个是回归方程的决定系数（R是相关系数），第2个是F统计量的值，第3个是 与F统计量对应的概率值p，第4个是剩余方差。得到上述回归模型的回归系 数估计值及其置信区间（置信水平0.05）、检验统计量、的结果见表参数参数估计值参数置信区间-12.8361，0.1903，0.0842，0.0735，-0.0003，-0.0002，由此可以得出该回归模型为同理，可得生菜与三种肥料施用量关系的模型二模型求解可得各参数的值参数参数估计值参数置信区间

9、-0.4938,0.0756,0.0234,0.0067,-0.0002,-0.0000,表2由此可以得出该回归模型为四模型的检验和分析在模型一的建立与求解过程中已对该模型进行了检验与分析，现在我们着重 检验分析模型二。图七土豆产量残差图首先我们分析土豆产量的残差图，从图七中发现在土豆产量的 30 个数据中， 仅有2 个数据（图中红线所示）的置信区间不包含零点，对于较多的点来说，残 差的置信区间均包含零点，这说明回归模型能够良好的符合原始数据。图八生菜产量残差图我们从图八中发现在生菜产量的 30个数据中，仅有 3个数据（图中红线所 示）的置信区间不包含零点，对于较多的点来说，残差的置信区间均包

10、含零点， 这说明回归模型能够良好的符合原始数据。现在我们对模型二的结果进行分析，由表1显示， 指因变量的 91.90%可由模型确定且，因而土豆产量的模型从整体来看是可用的。由表2显示指因变量的78.13%可由模型确定且，因而生菜产量的模型从整体来看也是可用的。表1的回归系数给出了土豆产量模型中的， ， ， ， ， ，的估计值，观察他们的置信区间发现，均不含零点，所以说明各项均对模型来说是显著的。表1的回归系数给出了生菜产量模型中的， ， ， ， ， ，的估计值，观察他们的置信区间发现，除外其他置信区间内均不含零点，所以说明 0 次项对 模型不显著。用以下的MATLAB程序求出土豆产量的最值ba

11、nan a二(x)12.8361-0.1903*x(1)-0.0842*x(2)-0.0735*x(3)+0.0003*(x(1厂2)+0.0002*(x(2厂2)+0.0001*(x(3厂2);x,fval=fminsearch(banana,0,0,0)求解得到最值，当，时，得出，由此结果我们可以得到氮磷钾肥料的最优配合比为1.5:1:1.74。我们可以发现，的值略小于实验数据中给出的最大值，下面我们分析产生 误差的几点原因。第一，给出的实验数据较少，不能很好的给出拟合模型。 第二，在模型二中，只考虑了氮、磷、钾三种肥料独立影响的情况，而忽略 了三者之间的交互影响。第三，在模型二中，缺乏正

12、交试验数据，无法给出较准确地模型。 综上所述，在模型二中的值应当较实际值偏小，但是我们可以看出土豆对钾肥 的喜好氮肥的喜好磷肥的喜好，土豆是喜钾作物。在生菜的模型二中，产量与磷肥钾肥成线性关系，因而无法求出最值，但我 们可以看到产量与氮肥成二次关系，因而氮肥影响更加显著，所以随着氮肥的增 加生菜产量先增加后减少，但磷肥钾肥却成线性关系，虽然符合题中所给出的实 验数据，但不符合实际情况。因此我们只能由实际数据得到最值，当 ， 时，得，可以看出生菜是喜磷作物。下面我们探究生菜产量的模型二与实际情况误差较大的原因 第一，在给出磷肥钾肥的实验数据中，我们可以发现，产量基本是随着施用 量的增加而增加，因

13、而拟合图像只能是线性关系，对于这种情况，应 该增添磷肥钾肥的实验数据来完善模型二。第二，在模型二中，只考虑了氮、磷、钾三种肥料独立影响的情况，而忽略 了三者之间的交互影响。我们通过上述分析发现，土豆与生菜的模型二中，均未对交互项的影响做出 探究，而氮磷钾三者的交互作用不可忽略，因此我们针对这一点对模型二进行改 进。五 模型的改进在土豆产量的模型二中增加四项后得我们把上述模型命名为模型三，在模型三中氮磷钾三种肥料不是相互独立的，其 中，为氮肥、磷肥之间的交互影响项，为氮肥、钾肥之间的交互影响项，为磷肥、 钾肥之间的交互影响项，为氮肥、磷肥、钾肥之间的交互影响项。模型三中的 、 、 称为交互项影响

14、系数，四个系数中的正值代表肥料的效果是相互促进的， 下面直接给出生菜产量的模型三由于所给数据都是固定两个水平，研究另一个的影响因素，题中没有给出正 交实验数据，因而模型三中的交互项系数无法求出。但是模型三的优点比较明显，可以较准确地拟合出实际情况，并可以给出较 准确的最值，以及可以比较几个交互影响的权重，若想得出模型，建议完善实验 数据并给出多组正交实验数据。六 模型的应用价值现在我们主要分析模型二的应用价值以及合理化建议，首先我们研究该模型 的经济价值。土豆的经济收益模型是生菜的经济收益模型是我们可以根据产量和三种肥料之间的关系，某时间内土豆生菜每吨的利润， 每公斤化肥的成本以及固定成本等，

15、利用上述两个经济收益模型，计算出在一定 耕地面积下，我们可以获得的最大的土豆生菜种植的经济收益。为此我们设计了一个 GUI 人机交互界面来预测最大收入201.6111127.1666 ZOO.B33322445 5.3O6利筒通过该模型可以指导我们合理使用化肥，降低对土壤酸度的影响，避免土壤板结、肥力下降，减少对土壤的污染，从而保护生态环境。七 参考文献1 姜启源,谢金星，叶俊.数学模型.4 版.高等教育出版社，2011.2 陈桂明，戚红雨，潘伟.MARLAB数理统计（6.x）.科学出版社，2002.3 辛磊.如何做残差分析.http:/www.doci 附录8.1 氮肥磷肥钾肥作用下产量与施肥量散点图的程序x1=0 34 67 101 135 202 259 336 404 471;y1=15.18 21.36 25.72