施肥效果分析大作业二

1、施肥效果分析摘 要 本文对作物生长所需的营养素与作物产量的关系及经济效益做了一定的讨论与分析。鉴于三种营养素混合考虑的复杂性，我们首先进行了单一营养素的分析，基本思想是先把三种营养素中的两种设为定值，把单一营养素的散点图通过MATLAB软件拟合出了施肥量与产量的关系方程及回归曲线图像，然后再由三种营养素的散点图通过SPSS软件拟合出了施肥量与产量的三元非线性方程及回归曲线图像，土豆与生菜的二次方程分别为： 并用LINGO对其应用价值进行了估价，从而，求得了获得最大经济效益所需的各营养素的量及最大经济效益，即最优方案的最优值，如表一所示：表一 因素作物N（kg/ha）P（kg/ha）K（kg/h

2、a）y（t/ha）最大经济效益（元/公顷）土豆292245.854245.252683.6生菜20535934121.111376.5关键词：散点图 非线性回归分析 一、问题重述施肥量(kg/ha)产量(t/ha)0346710113520225933640447115.8121.3625.7232.2934.0339.4543.1543.4640.8330.75施肥量(kg/ha)产量(t/ha)02449739814719624519434215.8121.3625.7232.2934.0339.4543.1543.4640.8330.75施肥量(kg/ha)产量(t/ha)0479314

3、018627937246555865118.9827.3534.8638.5238.4437.7338.4343.8742.7746.22 NP施肥量(kg/ha)产量(t/ha)049981471962943914895876856.399.4812.4614.3817.1021.9422.6421.3422.0724.53施肥量(kg/ha)产量(t/ha)0479314018627937246555865115.7516.7616.8916.2417.5619.2017.9715.8420.1119.40施肥量(kg/ha)产量(t/ha)028568411216822428033639

4、211.0212.7014.5616.2717.7522.5921.6319.3416.1214.11某地区作物生长所需的营养素主要是(N)、钾(K)、磷(P)。某作物研究所在该地区对土豆和生菜做了一定数量的实验，实验数据如下列表格所示，其中ha表示公顷，t表示吨，kg表示公斤。当一个营养素的施肥量变化时，总将另二个营养素的施肥量保持在第七个水平上，如对土豆产量关于N的施肥量分别取为196kg/ha与372kg/ha。试分析施肥量与产量之间关系，并对所得结果从应用价值与如何改进等方面做出估价。二、符号说明 表示氮肥的施肥量； 表示磷肥的施肥量； 表示钾肥的施肥量； 表示氮肥作为变量时土豆或生菜

5、的产量； 表示磷肥作为变量时土豆或生菜的产量； 表示钾肥作为变量时土豆或生菜的产量； 表示氮肥、磷肥、钾肥都为变量时土豆或生菜的产量； 表示氮肥作为变量时每公顷土豆或生生菜的纯收入； 表示磷肥作为变量时每公顷土豆或生菜的纯收入； 表示钾肥作为变量时每公顷土豆或生菜的纯收入； 表示氮肥、磷肥、钾肥都为变量时每公顷土豆或生菜的纯收入；三、模型假设 1气候、温度等外界自然条件适宜作物生长； 2作物水分充足，生长良好； 3其他营养素(除氮、磷、钾)都充足，不对作物生长造成影响； 4忽略所喷洒的农药对作物生长的影响； 5忽略施其他营养素(除氮、磷、钾)所需费用； 6假设土豆和生菜的价格是一定值，不随季节

6、变化；四、模型的建立与求解市场价格1NPK土豆生菜价格(元/千克)101861.41模型一：首先分析施氮肥的量与土豆产量的关系。把磷和钾的施肥量取为定值，分别为196kg/ha和372kg/ha；由题中施氮肥的量和产量的关系画出散点图，再由MATLAB2拟合得回归曲线，曲线图如图1所示：再经最小二乘法3求得拟合曲线方程： .由以氮肥为变量拟合的曲线方程对土豆的应用价值进行估价可知，每公顷土豆的纯属收入：由LINGO4解得：图1模型二：分析施磷肥的量与土豆产量的关系；把氮肥和钾肥的施肥量取为定值，分别为259kg/ha和372kg/ha；由题中施磷肥的量和产量的关系画出散点图，再由MATLAB拟

7、合得回归曲线，曲线图如图2所示：图2再经最小二乘法的拟合曲线方程：由以磷为变量拟合的曲线方程对土豆的应用价值进行估价可知，每公顷土豆的纯属收入：由LINGO解得：模型三：分析施钾肥的量与土豆产量的关系；把氮肥和磷肥的施肥量取为定值，分别为259kg/ha和372kg/ha；由题中施磷肥的量和产量的关系画出散点图，再由MATLAB拟合得回归曲线，曲线图如图3所示： 图3再经最小二乘法的拟合曲线方程：由以钾为变量拟合的曲线方程对土豆的应用价值进行估价可知，每公顷土豆的纯属收入：由LINGO解得：模型四：分析施氮肥、磷肥、钾肥的量与土豆产量的关系（氮肥、磷肥、钾肥都为变量），由题中氮肥、磷肥、钾肥与

8、土豆产量对应关系，再由SPSS拟合得回归曲线，曲线图如图4所示：施肥量与产量之间的关系NPK图4再经SPSS4软件得拟合曲线方程： 由以氮、磷、钾为变量拟合的曲线方程对土豆的应用价值进行估价可知，由LINGO解得每公顷土豆的纯属收入：，此时的 N（kg/ha）P（kg/ha）K（kg/ha）Y（t/ha）最大经济效益（元/公顷）N（PK变）279.819637243.2452094P（NK变）259163.737241.65450547K（NP变）52919621244.7847061NPK均变292245.854245.252683.6表二综上：模型一、模型二、模型三、模型四的最优解，如表二

9、所示：模型五：首先分析施氮肥的量与生菜产量的关系；把磷和钾的施肥量取为定值，分别为196kg/ha和372kg/ha；由题中施氮肥的量和产量的关系画出散点图，再由MATLAB拟合得回归曲线，曲线图如图5所示：图5再经最小二乘法求得拟合曲线方程： 由以氮肥为变量拟合的曲线方程对生菜的应用价值进行估价可知，每公顷生菜的纯属收入：由LINGO解得：模型六：分析施磷肥的量与生菜产量的关系；把氮肥和钾肥的施肥量取为定值，分别为259kg/ha和372kg/ha；由题中施磷肥的量和产量的关系画出散点图，再由MATLAB拟合得回归曲线，曲线图如图6所示：图6再经最小二乘法的拟合曲线方程：由以磷为变量拟合的曲

10、线方程对生菜的应用价值进行估价可知，每公顷生菜的纯属收入：由LINGO解得：模型七：分析施钾肥的量与生菜产量的关系；把氮肥和磷肥的施肥量取为定值，分别为259kg/ha和372kg/ha；由题中施磷肥的量和产量的关系画出散点图，再由MATLAB拟合得回归曲线，曲线图如图7所示： 图7再经最小二乘法的拟合曲线方程：由以钾为变量拟合的曲线方程对生菜的应用价值进行估价可知，每公顷生菜的纯属收入：由LINGO解得：模型八：分析施氮肥、磷肥、钾肥的量与生菜产量的关系（氮肥、磷肥、钾肥都为变量），由题中氮肥、磷肥、钾肥与生菜产量对应关系，再由MATLAB拟合得回归曲线，曲线图如图8所示：图8再经SPSS软

11、件得拟合曲线方程： 由以氮、磷、钾为变量拟合的曲线方程对生菜的应用价值进行估价可知，由LINGO解得每公顷生菜的纯属收入：此时综上：模型五、模型六、模型七、模型八的最优解，如下表三所示：表三N（kg/ha）P（kg/ha）K（kg/ha）y（t/ha）最大经济效益（元/公顷）N（PK变）187.539137220.549397.5P（NK变）224185.737216.416881.4K（NP变）224391221.220.498548.456NPK均变20535934121.111376.5五、模型的优缺点及改进优点：对于每种作物，我们先把单一营养素通过MATLAB绘出了模型图像，可以使我们

12、一目了然的看出氮、磷、钾三种营养素分别对于产量以及应用价值的关系，更加直观，更加易于理解；然后，研究讨论了三种营养素均变化时对产量的影响，得到了多元纯二次模型，而且，从上述模型和公式中可以看出产量及应用价值是受三种元素共同作用时，有最佳匹配关系，从而，求得了获得最大经济效益所需的各营养素的量及最大经济效益，通过一定的假设我们可以将一些问题简单化，从而对于实际生产有一定的指导作用。缺点：市场价格是受地点、季节和供求关系影响的，没有考虑人工、农药和灌溉的费用以及一些自然因素和人为因素对产量的影响。改进：根据不同地点、不同月份的市场价格，把每个月的单价进行修改，考虑一些人为和自然因素对产量的影响，考虑人工费、水费以及农药的费用六、参考文献【1】高全福，作物及肥料市场价格,网址:,200.8.11.【2】 刘卫国编, MATLAB程序设计教程,北京:中国水利水电出版社,2010.2.【3】 盛骤、谢式千、潘承毅编,概率论与数理统计,北京:高等教育出版社,2008.6.【4】 谢金星、薛毅编,优化建模与LINDO/LINGO软件,北京:清华大学出版,2008.7.【5】 林杰斌、林川雄、刘明德编,SPSS 12统计建模与应用实务,北京:中