农作物施肥效果分析.docx

上传人:b****8 文档编号:10960530 上传时间:2023-02-24 格式:DOCX 页数:18 大小:242.10KB
下载 相关 举报
农作物施肥效果分析.docx_第1页
第1页 / 共18页
农作物施肥效果分析.docx_第2页
第2页 / 共18页
农作物施肥效果分析.docx_第3页
第3页 / 共18页
农作物施肥效果分析.docx_第4页
第4页 / 共18页
农作物施肥效果分析.docx_第5页
第5页 / 共18页
点击查看更多>>
下载资源
资源描述

农作物施肥效果分析.docx

《农作物施肥效果分析.docx》由会员分享,可在线阅读,更多相关《农作物施肥效果分析.docx(18页珍藏版)》请在冰豆网上搜索。

农作物施肥效果分析.docx

农作物施肥效果分析

农作物施肥效果分析

摘要

本文通过观察实验数据及考虑现实情况,先放弃进行全因素试验,只在部分因素水平下进行试验,观察能否较全面的求出尽可能最大产量对应的N、P、K的施肥量。

即首先进行了忽略相关性,对作物进行单一因素分析,将N肥P肥和K肥三种因素两定一变,通过matlab软件将散点图汇出并通过最小二乘法拟合最相关的函数图像,分析找到产生最佳产量对应的最佳施肥量,建立出模型一。

鉴于三种营养素有相关性,可能存在1+1>2的效果,我们利用多元回归分析模拟出三元二次函数,并利用matlab编程求出最值,从而进行实验数据检验分析。

 

关键词:

最小二乘法回归分析三元二次函数matlab

 

一、问题重述

某研究所为了研究N、P、K三种肥料对于土豆和生菜的作用,分别对每种作物进行了三组实验,实验中将每种肥料的施用量分为10个水平,在考察其中一种肥料的施用量与产量关系时,总是将另二种肥料固定在第7个水平上,实验数据如下列表格所示,其中ha表示公顷,t表示吨,kg表示千克,试建立反映施肥量与产量关系的模型,并从应用价值和如何改进等方面作出评价.

表1-1施肥量对土豆产量的实验数据

 

N

P

K

施肥量

(kg/ha)

产量

(t/ha)

施肥量

(kg/ha)

产量

(t/ha)

施肥量

(kg/ha)

产量

(t/ha)

0

15.18

0

33.46

0

18.98

34

21.36

24

32.47

47

27.35

67

25.72

49

36.06

93

34.86

101

32.29

73

37.96

140

38.52

135

34.03

98

41.04

186

38.44

202

39.45

147

40.09

279

37.73

259

43.15

196

41.26

372

38.43

336

43.46

245

42.17

465

43.87

404

40.83

294

40.36

558

42.77

471

30.75

342

42.73

651

46.22

表1-2施肥量对生菜产量的实验数据

 

N

P

K

施肥量

(kg/ha)

产量

(t/ha)

施肥量

(kg/ha)

产量

(t/ha)

施肥量

(kg/ha)

产量

(t/ha)

0

11.02

0

6.39

0

15.75

28

12.70

49

9.48

47

16.76

56

14.56

98

12.46

93

16.89

84

16.27

147

14.38

140

16.24

112

17.75

196

17.10

186

17.56

168

22.59

294

21.94

279

19.20

224

21.63

391

22.64

372

17.97

280

19.34

489

21.34

465

15.84

336

16.12

587

22.07

558

20.11

392

14.11

685

24.53

651

19.40

二、问题的分析

在实际生活中,化肥对农作物具有二重性。

结合所学的初等函数,可以使用二次函数对数据进行初步的拟合。

做出数据散点图,利用matlab软件使用最小二乘法拟合找到相关函数,并进行初步的分析和判断。

生活中,考虑到N、P、K肥对农作物的影响,具有1+1>2的效果,所以不能仅仅考虑单一因素对其的影响,还需要进一步进行回归分析。

所以我们队决定拟合三项二次函数,将三者因素结合在一起。

进行数据回归分析和观察,并将模型一中的结果进行相关的检验。

三、建立模型

(一)基本假设

1.气候、温度等外界自然条件适宜作物生长,且处于稳定的状态。

2.作物水分充足,生长良好,其他营养素(除氮、磷、钾)都充足,不对作物生长造成影响;且忽略所喷洒的农药对作物生长的影响;

3.模型一是在单一因素中考虑,认为其他因素对产量的影响达到最大,且N、P、K不相互影响。

4.模型二认为各种肥料施放的效果是互相影响

(二)符号表示

(三)模型

模型一

使用matlab分别将N、P、K对农作物土豆的影响数据描绘出散点图。

如图3-1

图3-1N、P、K对农作物土豆效果的散点图

拟合选择二次函数,得到相对应的图象。

如图3-2

图3-2拟合二次函数

得到的二次函数分别为

产量最大值对应的自变量分别为

经分析可以得到,在仅考虑N、P、K肥单一因素对作物产生影响,即可得到每公顷土豆对其分别使用294.12kg,257.14kg,535.71kg的肥料,可使作物产生最大。

同理得到N、P、K对农作物生菜的拟合二次函数为

产量最大值对应的自变量分别为

即在仅考虑N、P、K肥单一因素对作物产生影响,即可得到每公顷土豆对其分别使用208.33kg,554.55kg,3541.7kg的肥料,可使作物产生最大。

同时通过对拟合函数的分析,含K肥料对作物影响可能会存在单一性,K含量越多,作物产量越大。

模型二

在实际施肥过程中,不同的肥料之间往往相互影响,使得各自的作用并不完全独立。

因此考虑采用三元二次多项式函数来拟合肥料施放实验的数据。

仍然使用matlab自带的语句拟合数据,考虑到施放肥料的量级比起产量太小,所以将肥料施放量级更改为百千克。

对于土豆,得到的函数为

计算得产量最大时,各肥料释放量如下:

对于生菜,得到的函数为

计算得产量最大时,各肥料释放量如下:

使用matlab实现函数求解最值程序如下

(1)土豆产量最值

函数文件

functionf=fun1(x)

f=-3.2578*x

(1)^2+8.2752*x

(1)+15.1518-1.7121*x

(2)^2+2.7090*x

(2)-0.6781*x(3)^2+2.7913*x(3)+2.1748*x

(1)*x

(2)+1.7440*x

(1)*x(3)+0.0198*x

(2)*x(3);

f=-f;

MATLAB求解步骤:

x0=[1;1;1];

Aeq=[];Beq=[];A=[];B=[];

vlb=[0;0;0];vhb=[];

[x,fval]=fmincon('fun1',x0,A,B,Aeq,Beq,vlb,vhb)

Localminimumfoundthatsatisfiestheconstraints.

Optimizationcompletedbecausetheobjectivefunctionisnon-decreasingin

feasibledirections,towithinthedefaultvalueoftheoptimalitytolerance,

andconstraintsaresatisfiedtowithinthedefaultvalueoftheconstrainttolerance.

 

x=

4.7680

3.8671

8.2461

 

fval=

-51.6266

-fval=51.6266

(2)生菜产量最值

函数文件

functionf=fun2(x)

f=-2.1910*x

(1)^2+2.5148*x

(1)+7.1369-0.3423*x

(2)^2+2.0551*x

(2)-0.3018*x(3)^2+1.4469*x(3)+1.1907*x

(1)*x

(2)+0.5801*x

(1)*x(3)-0.0433*x

(2)*x(3);

f=-f;

MATLAB求解步骤:

>>x0=[1;1;1];

Aeq=[];Beq=[];A=[];B=[];

vlb=[0;0;0];vhb=[];

[x,fval]=fmincon('fun2',x0,A,B,Aeq,Beq,vlb,vhb)

Localminimumfoundthatsatisfiestheconstraints.

Optimizationcompletedbecausetheobjectivefunctionisnon-decreasingin

feasibledirections,towithinthedefaultvalueoftheoptimalitytolerance,

andconstraintsaresatisfiedtowithinthedefaultvalueoftheconstrainttolerance.

 

x=

3.8132

9.2929

5.3953

 

fval=

-25.3838

-fval=25.3838

(四)模型扩展

建立模型,还需要考虑肥料的价格及农作物的经济价格,从而可以得到最有利于农民生产农作物的最佳施用量。

使其能够更贴近实际生活,而广泛使用。

四、模型分析及优缺点

模型一无论在拟合还是预测上都较为简单易懂,然而并未考虑各肥料施放量对于结果影响的独立性,仅仅适用于只能施放某一种肥料,且其他两种肥料土壤含量较为理想的情况;

模型二较为复杂,但是拟合与预测较精准,效果较好,适用情况广泛。

不过对于各肥料施放效果独立性的体现,较为初级。

模型一二都是基于肥料施放的效果与肥料施放量的函数曲线变化率不断减少这一特征,选择了二次多项式函数来拟合数据,实际上对肥料施放效果的产生原理生疏,无法利用合乎原理的函数,是本模型的最大缺憾。

五、参考文献

赵东方.数学模型与竞赛.湖北:

华东师范大学出版社.2014

中国牧草产业网,施肥量确定方法之四—肥料效应函数法

姜启源谢金星叶俊.数学模型[M].北京:

高等教育出版社.2011

戴朝寿.数理统计简明教程.北京:

高等教育出版社.2010

周品赵新芳.MATLAB数理统计分析.北京:

国防工业出版社.2009

 

六、附录

1.土豆对作物残差图

2.对生菜的产量影响拟合二次函数图像。

3.模型二拟合函数代码

>>[beta,r,j]=nlinfit(yy,pp,@zuoye1,a)

警告:

TheJacobianatthesolutionisill-conditioned,andsomemodelparametersmay

notbeestimatedwell(theyarenotidentifiable).Usecautioninmaking

predictions.

>Innlinfit(line376)

beta=

1至8列

-3.25788.275215.1518-1.71212.7090-0.67812.79132.1748

9至10列

1.74400.0198

r=

0.1512

0.2391

-0.5935

1.3687

-0.7460

-0.7173

0.6986

1.2845

2.1103

-1.6180

0.9256

-1.9855

-0.1868

0.1948

1.9028

-1.1154

-1.1914

-0.7053

-2.1171

1.4402

-5.5235

-0.4569

4.1101

5.0595

2.6166

-1.9939

-4.0214

-0.1360

-1.6176

2.6238

j=

1至8列

001.00003.84161.960013.83843.72000

0.11560.34001.00003.84161.960013.83843.72000.6664

0.44890.67001.00003.84161.960013.83843.72001.3132

1.02011.01001.00003.84161.960013.83843.72001.9796

1.82251.35001.00003.84161.960013.83843.72002.6460

4.08042.02001.00003.84161.960013.83843.72003.9592

6.70812.59001.00003.84161.960013.83843.72005.0764

11.28963.36001.00003.84161.960013.83843.72006.5856

16.32164.04001.00003.84161.960013.83843.72007.9184

22.18414.71001.00003.84161.960013.83843.72009.2316

6.70812.59001.00000013.83843.72000

6.70812.59001.00000.05760.240013.83843.72000.6216

6.70812.59001.00000.24010.490013.83843.72001.2691

6.70812.59001.00000.53290.730013.83843.72001.8907

6.70812.59001.00000.96040.980013.83843.72002.5382

6.70812.59001.00002.16091.470013.83843.72003.8073

6.70812.59001.00003.84161.960013.83843.72005.0764

6.70812.59001.00006.00252.450013.83843.72006.3455

6.70812.59001.00008.64362.940013.83843.72007.6146

6.70812.59001.000011.69643.420013.83843.72008.8578

6.70812.59001.00003.84161.9600005.0764

6.70812.59001.00003.84161.96000.22090.47005.0764

6.70812.59001.00003.84161.96000.86490.93005.0764

6.70812.59001.00003.84161.96001.96001.40005.0764

6.70812.59001.00003.84161.96003.45961.86005.0764

6.70812.59001.00003.84161.96007.78412.79005.0764

6.70812.59001.00003.84161.960013.83843.72005.0764

6.70812.59001.00003.84161.960021.62254.65005.0764

6.70812.59001.00003.84161.960031.13645.58005.0764

6.70812.59001.00003.84161.960042.38016.51005.0764

9至10列

07.2912

1.26487.2912

2.49247.2912

3.75727.2912

5.02207.2912

7.51447.2912

9.63487.2912

12.49927.2912

15.02887.2912

17.52127.2912

9.63480

9.63480.8928

9.63481.8228

9.63482.7156

9.63483.6456

9.63485.4684

9.63487.2912

9.63489.1140

9.634810.9368

9.634812.7224

00

1.21730.9212

2.40871.8228

3.62602.7440

4.81743.6456

7.22615.4684

9.63487.2912

12.04359.1140

14.452210.9368

16.860912.7596

%以上为拟合函数代码

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 工作范文 > 制度规范

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1