施肥效果分析

数学模型课程实验实验07施肥效果分析姓名:孟蕾学号:201510802008所在院系:数学与计算机学院专业:信息与计算科学指导老师：林宗兵摘要：对土豆和生菜分别绘制出他们的产量与三种营养元素之间关系的散图，拟合两变量之间的关系式。首先分别确定产量与施肥量之间的函数曲线类型，然后根据曲线类型对所求函数的对应关系进行假设，并利用已知数据计算出所需参数，最终确定变量之间的函数关系，得到最佳施肥量和最优产量。关键词：施肥方案散点图曲线拟合matlab一、问题重述：某地区作物生长所需的营养素主要是氮（N）、钾（K）、磷（P）。某作物研究所在某地区对土豆与生菜做了一定数量的实验，实验数据如下列表所示，其中ha表示公顷，t表示吨，kg表示公斤。当一个营养素的施肥量变化时，总将另两个营养素的施肥量保持在第七个水平上，如对土豆产量关于N的施肥量做实验时，P与K的施肥量分别取为196kg/ha与372kg/ha。若氮（N）、钾（K）、磷（P）和土豆、生菜的市场价格如表1所示：表1市场价格（元/吨）商品NPK土豆生菜价格350320640800200试分析施肥量与产量之间关系，并对所得结果从应用价值与如何改进等方面做出估计。表2土豆产量与施肥量的关系施肥量(N)(kg/ha)产量(t/ha)施肥量(P)(kg/ha)产量(t/ha)施肥量(K)(kg/ha)产量(t/ha)015.18033.46018.983421.362432.474727.356725.724936.069334.8610132.297337.9614039.5213534.039841.0418638.4420239.4514740.0927937.7325943.1519641.2637238.4333643.4624542.1746543.8740440.8329440.3655842.7747130.7534242.7365146.22表3生菜产量与施肥量的关系施肥量(N)(kg/ha)产量(t/ha)施肥量(P)(kg/ha)产量(t/ha)施肥量(K)(kg/ha)产量(t/ha)011.0206.39015.752812.70499.484716.765614.569812.469316.898416.2714714.3314016.2411217.7519617.1018617.5616822.5929421.9427919.2022421.6339122.6437217.9728019.3448921.3446515.8433616.1258722.0755820.1139214.1168524.5365119.40设计任务根据题目要求建立模型并求解：模型的应用与改进由于当一种肥料施肥量改变时，另外的两种肥料都保持在第7个水平上，于是有如下3个方案：(n,245,465)，(259,p,465)，(259,245,k)。对上述方案分别求出最大利润，然后进行比较就可得到最佳施肥方案。二、问题分析：利用散点图对所拟合问题的曲线类型做出判断。当需要拟合的两变量之间的函数关系式，首先要确定所求函数对应曲线的类型，然后根据曲线类型对所求函数的对应关系进行假设，并利用已知数据计算出所需参数，最终确定变量之间的函数关系。我们可以分别绘制出土豆和生菜的产量与施肥量的散点图，从图像的角度判断函数关系，再根据题目所给数据确定最终的函数。三、模型的建立与求解:散点图：土豆产量12233445ChOChOUlOUIOOUJ土豆产量涝密合KJ.§2吕3OOs5OOmoo苜

莆留融K++++□—aO百|<|O2CJ1O杪|O++++CJ1o■'土豆产量挣：舀醇所用matlab程序为：k1=xlsread（'E:\《数学建模课程设计》实验报告\shuju','sheet1','$L$3:$L$12'）;y31=xlsread（'E:\《数学建模课程设计》实验报告\shuju','sheet1','$M$3:$M$12'）;plot（k1,y31,'+'）土豆产量与施肥量的关系由散点图猜测生菜产量y与施肥量n的关系式为：＞二叩2+bn+c1y与磷肥的量p的函数为：＞=ap2+b2p+c2

3y与钾肥的量K的函数为:由matlab解出：a1=-0.0003b1=0.1971c1=14.7416a2=-0.0001b2=0.0719c2=32.9161a3=42.7b3=0.56c3=0.01*化回田*比回七土豆产量与施肥量的关系图:*化回田*比回七n1=xlsread（'E:\《数学建模课程设计》实验报告\shuju','sheet1','$A$3:$A$12'）;n2=n1.A2;y11=xlsread('E:\《数学建模课程设计》实验报告\shuju','sheet1','$B$3:$B$12');p1=xlsread('E:\《数学建模课程设计》实验报告\shuju','sheet1','$C$3:$C$12');p2=p1.A2;y21=xlsread('E:\《数学建模课程设计》实验报告\shuju','sheet1','$D$3:$D$12');k1=xlsread('E:\《数学建模课程设计》实验报告\shuju','sheet1','$E$3:$E$12');y31=xlsread('E:\《数学建模课程设计》实验报告\shuju','sheet1','$F$3:$F$12');c=ones(10,1);d1(:,1)=n2;d1(:,2)=n1;d1(:,3)=c;x1=inv(d1'*d1)*d1'*y11d2(:,1)=p2;d2(:,2)=p1;d2(:,3)=c;x2=inv(d2'*d2)*d2'*y21x0=[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.mfunctionf=shujunihe(x)c1=xlsread('E:\《数学建模课程设计》实验报告\shuju','sheet1','$E$3:$LE$12');c2=xlsread('E:\《数学建模课程设计》实验报告\shuju','sheet1','$F$3:$F$12');f=c2-x(1)*(1-x(2)*exp(x(3)*c1));上述文件保存为shujunihe.m用matlab解出最大利润为：y=37693最佳施肥方案为第一个方案(328.44,245,465)所用程序为：clearclca1=-0.0003;b1=0.1971;c1=14.742;a2=-0.0001;b2=0.0719;c2=32.916;a3=42.7;n=0:0.01:393;p=0:0.01:686;k=0:0.01:652;b3=0.56;c3=0.01;y1=(a1*n.*n+b1*n+c1)*800;y11=max(y1)fori=1:length(n)ifabs(y1(i)-y11)<=0.001q1=n(i)breakendendy2=(a2*n.*n+b2*n+c2)*800;y22=max(y2)fori=1:length(p)ifabs(y2(i)-y22)<=0.001q2=p(i)breakendendy3=a3*(1-b3*exp(-c3*k));y33=max(y3)fori=1:length(k)ifabs(y3(i)-y33)<=0.001q3=k(i)breakendendyll=3.7693e+004qi=328.4400y22=3.6672e+004q2=359.3900y33=42.6648运行后的结果如图：必二生菜产量与施肥量关系：649.2100由散点图猜测生菜产量y与施肥量N的关系式为：>二叩2+叩+C1y与磷肥的量P的函数为：＞ap2+bp+c：a\-bf3)y与钾肥的量K的函数为：'由matlab解出：a1=-0.0002b1=0.1013c1=10.2294a2=-0.0001b2=0.0606c2=6.8757a3=15.8878关系图为：b3-0.0440c3=0.0026

所用matlab程序为：clearclcn1=xlsread（'E:\《数学建模课程设计》实验报告\shuju','sheet1','$H$3:$H$12'）;n2=n1.A2;y11=xlsread（'E:\《数学建模课程设计》实验报告\shuju','sheet1','$I$3:$I$12'）;p1=xlsread（'E:\《数学建模课程设计》实验报告\shuju','sheet1','$J$3:$J$12'）;p2=p1.A2;y21=xlsread（'E:\《数学建模课程设计》实验报告\shuju','sheet1','$K$3:$K$12'）;k1=xlsread（'E:\《数学建模课程设计》实验报告\shuju','sheet1','$L$3:$L$12'）;y31=xlsread（'E:\《数学建模课程设计》实验报告\shuju','sheet1','$M$3:$M$12'）;c=ones（10,1）;d1(:,1)=n2;d1(:,2)=n1;d1(:,3)=c;x1=inv(d1'*d1)*d1'*y11d2(:,1)=p2;d2(:,2)=p1;d2(:,3)=c;x2=inv(d2'*d2)*d2'*y21x0=[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.mfunctionf=shujunihe(x)c1=xlsread('E:\《数学建模课程设计》实验报告\shuju','sheet1','$L$3:$L$12');c2=xlsread('E:\《数学建模课程设计》实验报告\shuju','sheet1','$M$3:$M$12');f=c2-x(1)*(1-x(2)*exp(x(3)*c1));上述文件保存为shujunihe.m用matlab解出最大利润为：y=18445最佳施肥方案为第一个方案(253.18,245,465)所用程序为：clearclca1=-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.001q1=n(i)breakendendy2=(a2*n.*n+b2*n+c2)*800;y22=max(y2)fori=1:length(p)ifabs(y2(i)-y22)<=0.001q2=p(i)breakendendy3=a3*(1-b3*exp(c3*k));y33=max(y3)fori=1:length(k)ifabs(y3(i)-y33)<=0.001q3=k(i)breakendendyll=1.8445e+004qi=253.1800y22二1.2845e+004q2二运行结果如图：302.8900y33二四、模型的评价与推卮61本模型利用Matlab编程，曲线估计较成功地解决了施肥最佳方案问题，方法简练,道理清晰，结果可信。■曲线估计得到较合适的