施肥效果分析.docx
- 文档编号:11638155
- 上传时间:2023-06-01
- 格式:DOCX
- 页数:19
- 大小:228.71KB
施肥效果分析.docx
《施肥效果分析.docx》由会员分享,可在线阅读,更多相关《施肥效果分析.docx(19页珍藏版)》请在冰点文库上搜索。
施肥效果分析
《数学建模》课程设计
报告
课题名称:
施肥效果分析
系(院):
理学院
数学与应用数学
2010.2011学年二学期
专业:
班级:
学生姓名:
学号:
指导教师:
开课时间:
*****
无名
091*********
许建强
摘要对土豆和生菜分别绘制出他们的产量与三种营养元素之间关系的散点图,拟合两变量之间的关系式。
首先分别确定产量与施肥量之间的函数曲线类型,然后根据曲线类型对所求函数的对应关系进行假设,并利用已知数据计算出所需参数,最终确定变量之间的函数关系,得到最佳施肥量和最优产量。
关键词:
施肥方案散点图曲线拟合matlab
一、问题重述:
某地区作物生长所需的营养素主要是氮(N)、钾(K)、磷(P)。
某作物研究所在某地区对土豆与生菜做了一定数量的实验,实验数据如下列表所示,其中ha表示
公顷,t表示吨,kg表示公斤。
当一个营养素的施肥量变化时,总将另两个营养素
的施肥量保持在第七个水平上,如对土豆产量关于N的施肥量做实验时,P与K的
施肥量分别取为196kg/ha与372kg/ha。
若氮(N)、钾(K)、磷(P)和土豆、生菜的市场价格如表1所示:
表1市场价格(元/吨)
商品
N
P
K
十豆
生菜二
价格
350
320
640
800
200
试分析施肥量与产量之间关系,并对所得结果从应用价值与如何改进等方面做出估计
表2土豆产量与施肥量的关系
施肥量(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
39.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.43J
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
表3生菜产量与施肥量的关系
施肥量(N)
「产量J
施肥量(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.33
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
【设计任务】
(1)根据题目要求建立模型并求解:
(2)模型的应用与改进
由于当一种肥料施肥量改变时,另外的两种肥料都保持在第7个水平上,于是有如
下3个方案:
(n,245,465),(259,p,465),(259,245,k)。
对上述方案分别求出最大利润,然后进行比较就可得到最佳施肥方案。
二、问题分析:
利用散点图对所拟合问题的曲线类型做出判断。
当需要拟合的两变量之间的函数关系式,首先要确定所求函数对应曲线的类型,然后根据曲线类型对所求函数的对应关系进行假设,并利用已知数据计算出所需参数,最终确定变量之间的函数关系。
我们可以分别绘制出土豆和生菜的产量与施肥量的散点图,从图像的角度判断函数关系,再根据题目所给数据确定最终的函数。
三、模型的建立与求解:
散点图:
40
20
50
00
50
42
086
4UJUJ
WKW+1
50
++
+
+
200250300350400450500
施肥量N
00
150200250300350
施肥量P
300400500600700
施肥量K
00
20
00
■■■・■■一■o
45403530252015
MKW41
+
+
++.
++-
+
150200250300350400
施肥量N
00
50
42086420
22211111
*<櫟州
24
22
20
18-
16-
14-
12
10
8
6o100200
300400500600700
施肥量P
20.6
20
19.5
19
18.5
175
17
1B5
16
155
0100200300400
施肥量K
所用matlab程序为:
k1=xlsread('《数学建模课程设计》实验报告\shuju'
y31=xlsread('《数学建模课程设计》实验报告\shuju'
plot(k1,y31,'+')
500600700
'sheetl',‘$L$3:
$L$12'
'sheetl','$M$3:
$M$12'
);
);
土豆产量与施肥量的关系由散点图猜测生菜产量y与施肥量N的关系式为:
y=务n2b,n•c1
2
y与磷肥的量P的函数为:
y=a2pb2pC2
y与钾肥的量K的函数为:
y=a31「b3eqk
由matlab解出:
a仁-0.0003
b1=0.1971
c1=14.7416
a2=-0.0001
b2=0.0719
c2=32.9161
a3=42.7
b3=0.56
c3=0.01
土豆产量与施肥量的关系图:
数据拟合图
数据槪合图
《数学建模课程设计》实验报告\shuju','sheet1','$B$3:
$B$12');
p1=xlsread('
p2=p”2;
《数学建模课程设计》
实验报告\shuju'
'sheet1'
y21=xlsread(
'《数学建模课程设计》
实验报告
\shuju'
'sheet1'
k1=xlsread('
《数学建模课程设计》
实验报告
\shuju'
'sheet1'
y31=xlsread(
'《数学建模课程设计》
实验报告
\shuju'
'sheet1'
'$C$3:
$C$12');
'$D$3:
$D$12');
'$E$3:
$E$12');
'$F$3:
$F$12');
c=ones(10,1);
d1(:
1)=n2;
d1(:
2)=n1;
d1(:
3)=c;
x1=inv(d1'*d1)*d1'*y11
d2(:
1)=p2;
d2(:
2)=p1;
d2(:
3)=c;
x2=inv(d2'*d2)*d2'*y21
x0=[420.550.05];
x3=lsqnonlin('shujunihe',x0)
n=0:
0.001:
393;
p=0:
0.001:
686;
k=0:
0.001:
652;
y1=x1
(1)*n.*n+x1
(2)*n+x1(3);
y2=x2
(1)*p.*p+x2
(2)*p+x2(3);
y3=x3
(1)*(1-x3
(2)*exp(x3(3)*k));
plot(k1,y31,'+',k,y3)
上述文件保存为qimobaogao.m
functionf=shujunihe(x)
'$E$3:
$LE$12');
'$F$3:
$F$12');
c1=xlsread('《数学建模课程设计》实验报告\shuju','sheet1'c2=xlsread('《数学建模课程设计》实验报告\shuju','sheet1'f=c2-x
(1)*(1-x
(2)*exp(x(3)*c1));
上述文件保存为shujunihe.m
用matlab解出最大利润为:
y=37693
最佳施肥方案为第一个方案(328.44,245,465)所用程序为:
clear
clc
a1=-0.0003;b1=0.1971;c1=14.742;
a2=-0.0001;b2=0.0719;c2=32.916;
a3=42.7;b3=0.56;c3=0.01;
n=0:
0.01:
393;
p=0:
0.01:
686;
k=0:
0.01:
652;
y1=(a1*n.*n+b1*n+c1)*800;
y1仁max(y1)
fori=1:
length(n)
ifabs(y1(i)-y11)<=0.001q1=n(i)
break
end
end
y2=(a2*n.*n+b2*n+c2)*800;
y22=max(y2)
fori=1:
length(p)
ifabs(y2(i)-y22)<=0.001q2=p(i)
break
end
endy3=a3*(1-b3*exp(-c3*k));y33=max(y3)
fori=1:
length(k)
ifabs(y3(i)-y33)<=0.001
q3=k(i)breakend
end
运行后的结果如图:
生菜产量与施肥量关系:
yll=
3.7693b4OD4
323.4430y22=
了.S672e+004
q?
=
y33=
42.S648
q3=
649.2100
a3=15.8878
关系图为:
24
22
20
咄16
14
12
W
050100150200250300350400
將肥量N
数据拟合
26
24
22
20
18
16
14
12
10
8
0100200300400500600700
施肥量P
数据拟合图
clear
轴忙Mw
clc
n1=xlsread(
n2=n1.A2;
y11=xlsread(p1=xlsread(p2=p1.A2;
y21=xlsread(
k1=xlsread(
y31=xlsread(c=ones(10,1);d1(:
1)=n2;d1(:
2)=n1;
《数学建模课程设计》
实验报告
\shuju'
'sheetl'
'$H$3:
$H$12'
《数学建模课程设计》
实验报告
\shuju'
'sheetl'
'$l$3:
$l$12'
《数学建模课程设计》
实验报告
\shuju'
'sheetl'
'$J$3:
$J$12'
《数学建模课程设计》
实验报告
\shuju'
'sheetl'
'$K$3:
$K$12'
《数学建模课程设计》
实验报告
\shuju'
'sheetl'
'$L$3:
$L$12'
《数学建模课程设计》
实验报告
\shuju'
'sheetl'
'$M$3:
$M$12'
);
);
);
);
);
);
d1(:
3)=c;
x1=inv(d1'*d1)*d1'*y11
d2(:
1)=p2;
d2(:
2)=p1;
d2(:
3)=c;
x2=inv(d2'*d2)*d2'*y21x0=[420.550.05];
x3=lsqnonlin(
n=0:
0.001:
393;
'shujunihe'
xO)
p=0:
0.001:
686;
k=0:
0.001:
652;
y1=x1
(1)*n.*n+x1
(2)*n+x1(3);
y2=x2
(1)*p.*p+x2
(2)*p+x2(3);
y3=x3
(1)*(1-x3
(2)*exp(x3(3)*k));
plot(k1,y31,'+',k,y3)
上述文件保存为qimobaogao.m
functionf=shujunihe(x)
'$L$3:
$L$12');
'$M$3:
$M$12');
c1=xlsread('《数学建模课程设计》实验报告\shuju','sheet1'
c2=xlsread('《数学建模课程设计》实验报告\shuju','sheet1'
f=c2-x
(1)*(1-x
(2)*exp(x(3)*c1));
上述文件保存为shujunihe.m
用matlab解出最大利润为:
y=18445
最佳施肥方案为第一个方案(253.18,245,465)所用程序为:
clear
clc
a1=-0.0002;b1=0.1013;c1=10.2294;
a2=-0.0001;b2=0.0606;c2=6.8757;
a3=15.8878;b3=-0.0440;c3=0.0026;
n=0:
0.01:
393;
p=0:
0.01:
686;
k=0:
0.01:
652;
y1=(a1*n.*n+b1*n+c1)*800;
y11=max(y1)
fori=1:
length(n)
ifabs(y1(i)-y11)<=0.001
q1=n(i)
break
end
end
y2=(a2*n.*n+b2*n+c2)*800;
y22=max(y2)
fori=1:
length(p)
ifabs(y2(i)-y22)<=0.001
q2=p(i)
break
end
end
y3=a3*(1-b3*exp(c3*k));
y33=max(y3)
fori=1:
length(k)
ifabs(y3(i)-y33)<=0.001
q3=k(i)breakend
end
yl1=
1.8445e-F004
Ql=
253.1800
y22=
1.2845e-F004
Q2=
302・8900
二
19-6961q3=
运行结果如图:
651.9000
四、模型的评价与推广:
4.模型优缺点
4.1模型优点
本模型利用Matlab编程,曲线估计较成功地解决了施肥最佳方案问题,方法简练,道理清晰,结果可信。
曲线估计得到较合适的曲线,最终得到拟合曲线函数表达式。
4.1模型缺点
在实际工作中,三种肥料之间除了与产量有直接的数量关系外,还有彼此之间的交互作用。
因此,本模型只是一个初步的探讨,要得到三种营养素与产量之间的准确关系,应该在实验之初就采取正交实验或均匀设计的方法,得到更有价
值的实验数据,从而更好的把握变量间的数量关系,以达到直到农业生产实践的目的。
五、参考文献:
1熊卫国•数学实验教程[M].广东:
中山大学出版社.2006.
2李玉莉.MATLAB函数速查手册[M].北京:
化学工业出版社.2010
3姜启源谢金星叶俊.数学模型[M].北京:
高等教育出版社.2010
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 施肥 效果 分析
![提示](https://static.bingdoc.com/images/bang_tan.gif)