施肥效果分析 摘要
本文研究了营养素对作物的产量的影响，分析了不同营养素对不同作物生长产量的差异，建立了施肥效果模型。并采用控制变量法和计算机数据拟合法建立了营养素对作物生长影响的模型。
根据研究所所得的营养素与作物产量的数据，运用MATLAB得到营养素与作物产量关系的散点图。进一步运用拟合工具进行拟合数据，得到多项式的二次，三次函数和正弦函数一项，两项和三项函数。利用方差比较，得到N在三次多项式时拟合度最好，而P和K在二次多项式时拟合度最好。
本文最后总结了模型的优点和不足之处，并对施肥效果改进意见。
关键词：散点图，方差比较，拟合方程，控制变量
一．问题重述
作物生长所需的营养素主要是氮（N）、钾（K）、磷（P）。为研究三种营养素对作物生长的影响，某作物研究所在该地区选取土豆与生菜做了一定数量的实验，实验过程中当一个营养素的施肥量变化时，总将另二个营养素的施肥量保持在最适宜植物生长状态。分析数据得出施肥量与产量之间关系，并对所得结果从应用价值与如何改进等方面作出估价。
二．问题分析
氮元素可促进植株茎叶的生长，更好的进行光合作用。磷元素具有一部分促根发育的作用还具有促进开花的作用。钾元素主要是促进果实的干物质积累，用来膨大果实。增加产量。由施肥量与产量的关系表格可得营养素对土豆生菜的产量有明显的促进作用。
根据农业期刊《Biology and fertility of soils》，一般来说，产量W可以用营养素施肥量的多项正弦函数表示，故做拟合曲线并代入试验数据求得关系表达式；同时联想到Logistic函数的导函数曲线为二次多项式（也是随着自变量先增后减），因此作一次二次以及三次多项式拟合，并进行比较。
三．基本假设
①每次试验独立且试验条件（如环境条件，种植密度，土壤条件）相同； ②由于数据由研究所提供，所以假设试验数据不会出现较大误差；
③三种元素的使施用量同作物产量有一定的函数关系，同一种元素对不同作物的作用表现为同一类的函数关系；
④忽略土壤中原有的N、P、K对作物生长的影响； ⑤三种元素对作物增长的作用是相互独立的；
四．名词解释和符号说明
名词解释:种植密度：单位面积作物种植量
符号说明：①pi（i=1,2,3.....)多项式系数
②ai，bi，ci正弦函数各项系数和常数项 五.模型建立和求解
采用MATLAB2012b中配置的curve fitting tool（曲线拟合工具），直接输入数据，进行曲线拟合。
二、N对作物生长的影响
1.N对土豆生长的影响 数据散点图
根据图像可进行大致估计，其函数关系应为二次（或更高次）多项式函数，此函数关系亦有可能为三角函数。根据实际，现假设其函数关系可能为二次多项式函数、三次多项式函数、一项正弦函数、两项正弦函数。根据已有数据建立回归模型，进行曲线拟合。（以下关于其他几组数据的讨论皆为这几种模型，不再赘述） （1） 二次多项式函数回归模型
结果如下Linear model Poly2: f(x) = p1*x^2 + p2*x + p3
Coefficients (with 95% confidence bounds):（95%致信区间内的拟合常数） p1 = -0.0003395 (-0.0003886, -0.0002905) p2 = 0.1971 (0.1736, 0.2207) p3 = 14.74 (12.63, 16.85)
Goodness of fit: SSE: 11.33 R-square: 0.9863 Adjusted R-square: 0.9824 RMSE: 1.272 （2）三次多项式函数回归模型
Linear model Poly3:
f(x) = p1*x^3 + p2*x^2 + p3*x + p4 Coefficients (with 95% confidence bounds):
p1 = -2.666e-07 (-5.998e-07, 6.655e-08) p2 = -0.0001532 (-0.0003899, 8.347e-05) p3 = 0.1645 (0.1189, 0.2102) p4 = 15.71 (13.5, 17.91)
Goodness of fit: SSE: 6.914 R-square: 0.9916 Adjusted R-square: 0.9875 RMSE: 1.073 (3)单项正弦函数回归模型
General model Sin1:
f(x) = a1*sin(b1*x+c1)
Coefficients (with 95% confidence bounds): a1 = 43.76 (42.27, 45.25)
b1 = 0.004162 (0.003909, 0.004416) c1 = 0.3586 (0.3058, 0.4115)
Goodness of fit: SSE: 10.32 R-square: 0.9875 Adjusted R-square: 0.984 RMSE: 1.214 (4)两项正弦函数回归模型
General model Sin2:
f(x) = a1*sin(b1*x+c1) + a2*sin(b2*x+c2) Coefficients (with 95% confidence bounds):
a1 = 92.31 (-1.453e+04, 1.472e+04) b1 = 0.006413 (-0.1468, 0.1596) c1 = -0.02898 (-27.93, 27.88)
a2 = 50.29 (-1.459e+04, 1.469e+04) b2 = 0.008166 (-0.2279, 0.2442) c2 = 2.783 (-42.74, 48.31)
Goodness of fit:
SSE: 3.436 R-square: 0.9958 Adjusted R-square: 0.9906 RMSE: 0.9268 2．N对生菜生长的影响
（1）二次多项式函数回归模型