1、价格350 320640 800200试分析施肥量与产量之间关系,并对所得结果从应用价值与如何改进等方面做出估计。 表2 土豆产量与施肥量的关系施肥量(N)(kg/ha) 产量 (t/ha) 施肥量(P) 施肥量(K) 0 15.18 33.46 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.4
2、3 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 生菜产量与施肥量的关系11.02 6.39 15.75 28 12.70 9.48 16.76 56 14.56 12.46 16.89 84 16.27 14.33 16.24 112 17.75 17.10 17.56 168 22.59 21.94 19.20 224 21.63 391 22.64 17.97 280 19.34 489 21.34 15.84 16.12 587 22.07 20.
3、11 392 14.11 685 24.53 19.40 设计任务 (1)根据题目要求建立模型并求解:(2)模型的应用与改进由于当一种肥料施肥量改变时,另外的两种肥料都保持在第7个水平上,于是有如下3个方案:(n,245,465),(259,p,465),(259,245,k)。对上述方案分别求出最大利润,然后进行比较就可得到最佳施肥方案。二、 问题分析:利用散点图对所拟合问题的曲线类型做出判断。当需要拟合的两变量之间的函数关系式,首先要确定所求函数对应曲线的类型,然后根据曲线类型对所求函数的对应关系进行假设,并利用已知数据计算出所需参数,最终确定变量之间的函数关系。我们可以分别绘制出土豆和生
4、菜的产量与施肥量的散点图,从图像的角度判断函数关系,再根据题目所给数据确定最终的函数。三、 模型的建立与求解:散点图:所用matlab程序为:k1=xlsread( E:数学建模课程设计实验报告shuju,sheet1$L$3:$L$12);y31=xlsread($M$3:$M$12plot(k1,y31,+)土豆产量与施肥量的关系由散点图猜测生菜产量y与施肥量N的关系式为:y与磷肥的量P的函数为:y与钾肥的量K的函数为: 由matlab解出: a1=-0.0003 b1=0.1971 c1=14.7416 a2=-0.0001 b2=0.0719 c2=32.9161 a3=42.7 b3
5、=0.56 c3=0.01土豆产量与施肥量的关系图:所用的matlab程序为:clearclcn1=xlsread($A$3:$A$12n2=n1.2;y11=xlsread($B$3:$B$12p1=xlsread($C$3:$C$12p2=p1.2;y21=xlsread($D$3:$D$12$E$3:$E$12$F$3:$F$12c=ones(10,1);d1(:,1)=n2;,2)=n1;,3)=c;x1=inv(d1*d1)*d1*y11d2(:,1)=p2;,2)=p1;x2=inv(d2*d2)*d2*y21x0=42 0.55 0.05;x3=lsqnonlin (shujun
6、ihe,x0)n=0:0.001:393;p=0:686;k=0: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);,k,y3)上述文件保存为qimobaogao.mfunction f=shujunihe(x)c1=xlsread($LE$12c2=xlsread(f=c2-x(1)*(1-x(2)*exp(x(3)*c1);上述文件保存为shujunihe.m用matlab解出最大利润为:y=37693最佳施肥方案为第一个方案(328.44,245,465)所
7、用程序为: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;0.01:y1=(a1*n.*n+b1*n+c1)*800;y11=max(y1)for i=1:length(n) if abs(y1(i)-y11)=0.001 q1=n(i) break endendy2=(a2*n.*n+b2*n+c2)*800;y22=max(y2)length(p) if abs(y2(i)-y22) q2=p(i)y3=a3*(1-b3*exp(-c3*k);y33=ma
8、x(y3)length(k) if abs(y3(i)-y33) q3=k(i)运行后的结果如图:生菜产量与施肥量关系: 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关系图为:$H$3:$H$12$I$3:$I$12$J$3:$J$12$K$3:$K$12y=18445最佳施肥方案为第一个方案(253.18,245,465)a1=-0.0002; b1=0.1013; c1=10.2294; b2=0.0606; c2=6.8757;a3=15.8878
9、; b3=-0.0440; c3=0.0026;y3=a3*(1-b3*exp(c3*k);运行结果如图:四、模型的评价与推广:本模型利用Matlab编程,曲线估计较成功地解决了施肥最佳方案问题, 方法简练, 道理清晰, 结果可信。曲线估计得到较合适的曲线,最终得到拟合曲线函数表达式。在实际工作中, 三种肥料之间除了与产量有直接的数量关系外,还有彼此之间的交互作用。因此, 本模型只是一个初步的探讨, 要得到三种营养素与产量之间的准确关系, 应该在实验之初就采取正交实验或均匀设计的方法, 得到更有价值的实验数据, 从而更好的把握变量间的数量关系, 以达到直到农业生产实践的目的。五、参考文献:熊卫国 . 数学实验教程M. 广东: 中山大学出版社 . 2006.李玉莉 . MATLAB函数速查手册M. 北京:化学工业出版社 . 2010姜启源 谢金星 叶俊 .数学模型M . 北京:高等教育出版社 . 2010
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1