数学实验施肥效果+代码报告Word下载.docx
- 文档编号:18437500
- 上传时间:2022-12-16
- 格式:DOCX
- 页数:12
- 大小:215.90KB
数学实验施肥效果+代码报告Word下载.docx
《数学实验施肥效果+代码报告Word下载.docx》由会员分享,可在线阅读,更多相关《数学实验施肥效果+代码报告Word下载.docx(12页珍藏版)》请在冰豆网上搜索。
kg/ha
11.02
0
6.39
15.75
28
12.70
49
9.48
47
16.76
56
14.56
98
12.46
93
16.89
84
16.72
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
试分析施肥量与产量之间的关系,并对所得结果从应用价值与如何改进等方面做出估价。
三、模型的建立与问题的分析
问题的分析与模型的建立
由题意知,当一个营养素的施肥量发生改变时,总将另外两个营养素的肥量保持在第7个水平上,因此找出一个营养素的变化与作物产量之间的关系,即利用给出的数据,找出生菜产量与各元素之间的函数关系
。
施用氮肥的平均效率:
(22.59-11.02)/168=0.0689;
施用磷肥的平均效率:
(24.53-6.39)/685=0.0265;
施用钾肥的平均效率:
(20.11-15.75)/558=0.0078;
我们从数据中我们可以发现:
当氮肥增加时产量逐步增加,但当增加到一定程度的时候产量反而减少,这就是农业生产中氮肥过量使用会造成烧苗的原因。
从磷肥与产量的数据可以发现:
磷肥增加时,产量逐步增加,但增加速率随着磷肥的增加而减小。
钾肥与产量的关系与上述两种肥料都不同,从钾肥与产量的数据可以发现:
钾肥的使用对产量的影响不大。
则我们有如下三种方案:
(n,391,372),(224,p,372),
(224,391,k);
由于钾肥与产量的关系与氮肥和磷肥的关系都不同,所以选择用曲线回归的方法建立生菜产量y与施肥量之间函数关系,具体的步骤如下:
(1)作出散点图
生菜产量与施氮肥效应
生菜产量与施磷肥效应
生菜产量与施钾肥效应
(2)建立回归曲线方程
设y与n(氮肥的量)的函数为:
y=an
+bn+c;
利用最小二乘法求下式成立的函数y
=min
则用MTTLAB软件:
我们得到产量y与N(氮肥的量)的函数为:
y(n)=-0.0002n
+0.1013n+10.2294;
同理,我们利用最小二乘法求y与p(磷肥的量)之间的函数关系式为:
y(p)=-0.0001p
+0.0606p+6.8757
而钾肥的使用量与产量的关系来看其比较特殊,故改用非线性曲线拟合的方法对其建立数学函数模型求解;
从图象我们可以看出其符合指数函数;
我们设函数为:
y(k)=a(1-be
),其中a,b,c均为待定系数,我们用数据带入,分别再与原函数图像比较拟合,得到以下函数关系式:
y(k)=25.0467*(1-0.3537*e
)
四模型的应用与改进
在上面我们建立了氮、钾、磷肥与产量之间的关系,利用上述函数关系我们可以进行定量分析计算各种肥料的最佳使用量;
题中已给出了氮、磷、钾与生菜的市场价如图所示
市场价格(元/吨)
商品
氮肥
磷肥
钾肥
生菜
价格
350
320
640
200
由于一种肥料改变时,另外的两种的肥料都保持在第7种水平上,于是我们可以有下面的三种方案:
(n,391,372),(224,p,372),(224,391,k)
为了寻求最佳的施肥效果只需寻找到最大的利润即可:
第一种:
设每公顷加的氮肥n千克时的利润为
Ln(n)=200y(n)-(0.35n+0.32*391+0.64*372)
-=-0.04n
+19.91n+1682.7
由Ln’(n)=-0.08n+19.91=0
n=248.87,而Ln’’(n)=-0.08<
0;
即n=248.87时函数有极大值,而且为问题的最大值,最大的利润为:
Ln_max(248.87)=4160.3(元)
第二种:
设每公顷加的磷肥p千克时的利润为
Lp(p)=200y(p)-(0.35*224+0.32p+0.64*372)
=-0.02p
+11.8p+1058.7
由Lp’(p)=-0.04p+11.8=0
p=295.00
而其二元导数为负值,即p=295.00时即为函数的最大值
其最大利润为:
Lp_max(295.00)=2799.2(元)
第三种:
设每公顷加的钾肥k千克时的利润为
Lk(k)=200y(k)-(0.35*224+0.32*391+0.64k)
=4805.8-1771.8*exp(-6.4114e-004*k))-0.64*k
可求出其最大值为:
Lk(700)=3226.7
综上可知最佳的施肥方案为第一种即(249,391,372)
根据农作物生长的原理,氮、磷、钾3种肥料缺一不可,但又是一个有机的整体,因此,要得到农作物的产量与3种肥料之间的使用量的关系,必须考虑3种肥料之间交互影响的数据,也就是说在设计实验时应采取正交实验,或均匀设计的方法,利用这样实验得到的数据建立农作物产量与3种肥料间的多元函数关系,才能比较准确的找到最佳施肥量。
五、实验结果
数学实验MATLAB施肥效果分析
*********************************************************************
根据离散点图像,用线性拟合的方法得到拟合函数的三个系数:
a=
-2.3815e-004
b=
0.1013
c=
10.2294
生菜的产量与氮肥的施肥量函数关系为:
y_n=
Inlinefunction:
y_n(n)=-0.0002*n^2+0.1013*n+10.2294
使用氮肥时生菜的利润函数关系为:
Ln=
Ln(n)=-0.04*n^2+19.91*n+1682.7
此时氮肥的使用量为:
(吨)时,利润最大
n=
248.8750
使用氮肥的最大利润为:
(元)
maxn=
4.1603e+003
a1=
-5.4528e-005
b1=
0.0606
c1=
6.8757
生菜的产量与磷肥的施肥量函数关系为:
y_p=
y_p(p)=-0.0001*p^2+0.0606*p+6.8757
使用磷肥时生菜的利润函数关系为:
Lp=
Lp(p)=-0.02*p^2+11.8*p+1058.7
此时磷肥的使用量为:
p=
295.0000
使用磷肥的最大利润为:
maxp=
2.7992e+003
根据离散图像,其图像符合一指数函数,用非线性拟合的方法设其拟合函数为:
y_k=
y_k(a,k)=a
(1)*(1-a
(2)*exp(-a(3)*k))
a2=
25.0467
b2=
0.3537
c2=
6.4114e-004
生菜的产量与钾肥的施肥量函数关系为:
y_k(k)=25.0467*(1-0.3537*exp(-6.4114e-004*k))
使用钾肥时生菜的利润函数关系为:
Lk=
Lk(k)=25.0467*(1-0.3537*exp(-6.4114e-004*k))*200-224*0.35-391*0.32
-0.64*k
此时钾肥的使用量为:
k=
700.0000
使用钾肥的最大利润为:
maxk=
3.2267e+003
综上所知,最佳的施肥方案为第一种方案(249,391,372)
离散图和拟合后的图形如下所示:
全部代码:
functionshiyan
clc
disp('
数学实验MATLAB施肥效果分析'
***********************************************************************************************'
n=[0285684112168224280336392];
yn=[11.0212.7014.5616.2717.7522.5921.6319.3416.1214.11];
N=polyfit(n,yn,2);
x1=0:
1:
400;
y1=polyval(N,x1);
subplot(2,2,1);
plot(n,yn,'
g'
x1,y1,'
b'
xlabel('
氮肥的使用量'
ylabel('
生菜产量'
legend('
离散图'
'
拟合图'
4)
'
symsabc
a=N
(1),b=N
(2),c=N(3)
y_n=inline('
-0.0002*n^2+0.1013*n+10.2294'
Ln=inline('
-0.04*n^2+19.91*n+1682.7'
)%为了显示
-(-0.04*n^2+19.91*n+1682.7)'
);
%为了计算最大值,前加‘-’号,同理,以下一样
[x,fval]=fminbnd(Ln,0,400);
(吨)时,利润最大'
n=x
(元)'
maxn=abs(fval)
p=[04998147196294391489587685];
yp=[6.399.4812.4614.3317.1021.9422.6421.3422.0724.53];
P=polyfit(p,yp,2);
x2=0:
700;
y2=polyval(P,x2);
subplot(2,2,2);
plot(p,yp,'
x2,y2,'
磷肥的使用量'
symsa1b1c1
a1=P
(1),b1=P
(2),c1=P(3)
y_p=inline('
-0.0001*p^2+0.0606*p+6.8757'
Lp=inline('
-0.02*p^2+11.8*p+1058.7'
)
-(-0.02*p^2+11.8*p+1058.7)'
[x,fval]=fminbnd(Lp,0,700);
p=x
maxp=abs(fval)
k=[04793140186279372465558651];
yk=[15.7516.7616.8916.2417.5619.2017.9715.8420.1119.40];
symsa2b2c2
y_k=inline('
a
(1)*(1-a
(2)*exp(-a(3)*k))'
a'
k'
A=nlinfit(k,yk,y_k,[14.52,0.06,0.02]);
a2=A
(1)
b2=A
(2)
c2=A(3)
y3=a2*(1-b2*exp(-c2*k));
subplot(2,1,2);
plot(k,yk,'
k,y3,'
)%必须有相同的点,否则运行出现错误
钾肥的使用量'
25.0467*(1-0.3537*exp(-6.4114e-004*k))'
Lk=inline('
25.0467*(1-0.3537*exp(-6.4114e-004*k))*200-224*0.35-391*0.32-0.64*k'
-(25.0467*(1-0.3537*exp(-6.4114e-004*k))*200-224*0.35-391*0.32-0.64*k)'
[x,fval]=fminbnd(Lk,0,700);
k=x
maxk=abs(fval)
综上所知,最佳的施肥方案为第一种方案(249,391,372)'
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 实验 施肥 效果 代码 报告
![提示](https://static.bdocx.com/images/bang_tan.gif)