探究鲈鱼的质量估计的方法及应用Word格式.docx

1、符号说明单位x1鱼的身长cmx2鱼的胸围y鱼的质量gV鱼的体积cm3R鱼的最大周长所对应半径同时将题中所给的统计数据由左向右依次标号，鱼的身长依次为：x11,x12x18，鱼的质量依次为：y1,y2y8，鱼的胸围依次为：x21,x22x28。3 模型建立模型一 机理分析法分别探究身长和胸围对鱼质量的影响4.1.1 一次函数模型假设鱼质量和身长以及鱼质量和胸围的关系都为一次函数关系，具体记为：y=a1x1+b1y=a2x2+b2其中：a1，b1，a2，b2都是参数。参数可以用题中所给的统计数据应用最小二乘法拟合得到，具体操作如下：求a1,b1,a2,b2使目标函数满足条件,目标函数为：miny（

2、xk）-yk2 （n=1,28）其中y（xk）是鱼质量的估计值。通过matlab对模型的参数进行参数估计，求解模型得：a1=65.3b1=-1637.3a2=92b2=-1497.5程序详见附表1。 4.1.2 二次函数模型假设鱼质量和身长以及鱼质量和胸围的关系都为二次函数关系，具体记为：y=a1x12+b1x2+c1y=a2x22+b2x2+c2a1,b1,c1,a2,b2,c2都是参数。参数估计方法同一次函数模型。通过matlab对模型的参数进行参数估计，可得：a1=1.6247b1=-59.3124c1=709.7392a2=1.3b2=157.9c2=-2344.8程序详见附表2。4.

3、1.3 三次函数模型假设鱼质量和身长以及鱼质量和胸围的关系都为三次函数关系，具体记为：y=a1x13+b1x12+c1x1+d1y=a2x23+b2x22+c2x2+d2a1,b1,c1,d1,a2,b2,c2,d2都是参数。参数估计方法同上，可得：a1=1b1=-80c1=3008d1=-37262a2=-1b2=90c2=-2228d2=18113程序详见附表3。 误差分析利用方差对上述结果进行误差分析，其实质上就是比较三种函数模型的目标函数值的大小，目标函数数值小的误差小。由此可得:方法身长误差胸围误差一次函数关系1.6543e+0041.6258e+004二次函数关系9.8461e+0

4、032.1530e+007三次函数关系2.0881e+0095.0456e+007程序详见附表4。显然，二次函数关系模型对鱼质量和身长的关系估计的更准确，一次函数关系对鱼质量和胸围的关系估计的更准确。模型二 机理分析法探究身长和胸围共同对鱼质量的影响由模型一知胸围和身长对质量都有影响，因此在模型二中考虑二者共同对质量的影响，采用两种方法。4.2.1 加权平均法对模型一的进一步处理对模型一中胸围和身长分别对质量的影响加权平均。加权平均法是一种根据各种因素对结果影响大小对因素加不同权重的方法。公式：y=am+bn（其中a,b分别为权重且满足a+b=1，m,n分别为影响因素，若影响因素不止两个，可以

5、按此格式继续扩展）对于此研究课题影响因素分别为身长和胸围，在模型一中已得出身长对质量的影响函数为：Y1=y=1.6247x12-59.3124x1+709.7392胸围对质量的影响函数为：Y2=y=92x2-1497.5对两种因素加权重，此课题不妨取权重分别为0.5（根据需要也可取权重分别为0.4,0.6或0.6,0.4等）得出加权模型：Y=0.5Y1+0.5Y2在概率论与数理统计中我们已经学过，方差具有一定的计算法则，例如：D（ax+b）=a2D（x）。我们很容易得到该模型的误差，误差为：6.5260e+0034.2.2 三维空间基本曲线模型假设鱼质量和鱼身长、鱼胸围的关系为三维空间基本曲线

6、函数关系，具体记为：y=ax12+bx22+cx1x2+dx1+ex2+f其中a,b,c,d,e,f均为参数。参数估计方法同模型一，模型求解得：a=-0.0001b=0.7686c=1.1436d=0.0432e=0.0289f=0.0011程序详见附表5。利用方差对上述结果进行误差分析，误差结果为：8.1896e+007模型三 机理分析法探究体积对鱼质量的影响关于鱼的最大周长到底处于鱼体的何位置，在我们所建立的鱼的近似几何体（如图）中，x2=2R;R=V1=R2l1; V2=l2;V =l1 +l2=（l1+l2）R2l;所以鱼的体积和鱼的胸围处于鱼体何处无关。则：V = = ;我们已假设鱼

7、的质量均匀，则设：m=aV+ba，b为正参数。参数估计方法同上，模型求解得：a=1.1587b=42.4279程序详见附表6。利用方差对上述结果进行误差分析，误差分析结果为：7.0276e+0044 模型评价5.1模型比较模型误差模型一9.8461e+003 1.6258e+004加权平均模型三维空间基本曲线模型体积质量模型由上表可知，加权平均模型是最优之选。5.2优点评述1） 模型不仅可以应用于本题的背景，其实，在鱼苗鱼种的生产中，需要对鱼体的成长情况观察了解，随着望向培育鱼种工艺的应用及有关的实验项目的展开，对鱼体长度和重量的测定更成为一种经常性的工作内容。在工作量较大数量较多时，适用常规

8、测量和称重法，不过那难度大，而且容易导致被测由于受伤，甚至死亡。我们一估计出鱼体长度胸围质量之间存在着某种程度的统计关系，大大方便了研究工作，减少不必要的损失。2） 本文提出的三种估计模型具有普遍适用性。3） 本文通过对估计模型一精确率的探索提出新的加权平均估计模型，具有创新性。4） 对于用长度估计鱼的质量的要求，我们提出三种不同的模型并作了对比和评述，使公平性得到了更好的体现。5.3模型改进1） 理论上讲，模型三的截距应为零，因为我们已经假设鱼的质量是均匀的，密度是一定的，质量等于体积和密度乘积。可以从三方面改进：a） 将该截距看成是残差，通过数据代入我们可以发现在鱼体积不太小的情况下估计还

9、是可靠的。b） 模型三之所以误差稍大，可能因为我们对鱼的几何体模型构造的不准确，可以通过对生物学相关书籍的查阅来完善，进而构造出更精确的立体几何模型。c） 如果条件允许，可以用排水法测定体积，是用一般具有刻度的玻璃量筒，先置水于量筒容量三分之一或二分之一的某一刻度，然后将鱼逐尾投入，并分别从页面上升程度几下增加的毫升读数，此读数记为该鱼的体积值。再用经验获得的鱼的密度值根据公式计算鱼的质量。该种改进方法，鱼不离水，操作便捷，可行性实用性强。2） x1/x2是鱼身长与最大胸围之比，如果x1/x2太大，流体力学角度来讲鱼在水中遇到的阻力增大；如果x1/x2太小，其自身的生长不能达到一种自然吻合，无

10、疑是不利于生存，因此在查阅生物学的有关书籍以及从达尔文进化论角度来讲,可以假定，经过长期进化，对于每一种动物而言x1/x2已经达到其最适合的数值，换句话说，x1/x2应视为与这种鱼的尺寸无关的常数，于是可得到：x1x2由于池塘里的鲈鱼体形都是相似的，对于两条鱼而言，由数学中相似原理得：sx22其中s为鱼的横截面积，又由于：ss1 其中s1为鱼的平均横截面积，则有：sx12由体积公式：Vs1*x1可得：Vs*x1综上比例关系可得:Vx13即得模型：mx13为了检验我们设：m=alb，其中a,b为待定参数，又：log10m=a+blog10l利用最小二乘法根据所给数据拟合上式得到m=0.01322

11、 l3.0265可以看出模型与这个结果吻合的相当好。详细图表：参考文献1 姜启源 谢金星 叶俊 编，数学模型（第四版），北京：高等教育出版社，2011年1月2 季之源 戴俊杰 顾嘉宾 周玉溪，计算幼鱼长度和重量的体积测量法及其应用，江苏省兴化县水产科学研究所，1989年附 录附录1 一次函数参数估计程序x=31.8,32.1,32.1,35.9,36.8,36.8,43.8,45.1;y=482,482,454,652,765,737,1162,1389;p=polyfit（x,y,1）x1=31.8:0.1:45.1;y1=polyval（p,x1）;plot（x,y,o,x1,y1）;ti

12、tle（）x=21.3,21.6,21.6,22.9,24.8,24.8,27.9,31.8;x1=21.3:31.8;程序的运行结果：p =1.0e+003 * 0.0653 -1.6373 p =1.0e+003 * 0.0920 -1.4975 附录2二次函数参数估计程序p=polyfit（x,y,2）p =1.6247 -59.3124 709.7392 p =1.0e+003 * -0.0013 0.1579 -2.3448附录3三次函数参数估计程序p=polyfit（x,y,3）p =1.0e+004 * 0.0001 -0.0080 0.3008 -3.7262 p =1.0e+

13、004 * -0.0001 0.0090 -0.2228 1.8113 附录4模型一误差分析程序x1=31.8,32.1,32.1,35.9,36.8,36.8,43.8,45.1;y1=65.3*x1-1637.3;k=sum（y1-y）.2）y1=1.6247*（x1.2）-59.3124*x1+709.7392;y1=x1.3-80*（x1.2）+3008*x1-37262;x2=21.3,21.6,21.6,22.9,24.8,24.8,27.9,31.8;y1=92*x2-1497.5;y1=1.3*（x2.2）+157.9*x2-2344.8;y1=-1*（x2.3）+90*（x2

14、.2）-2228*x2+18113;附录5三维空间曲线基本函数参数估计程序plot3（x1,x2,y,o-hold ony1=1.7044*（x1.2）+0.7686*（x2.2）+1.1436*（x1.*x2）+0.0432*x1+0.0289*x2+0.0011plot3（x1,x2,y1,rp-）;a=1,1,1,1,1,1,1,1;wei=x1.2;x2.2;x1.*x2;x1;x2;a;wei/yans = 1.7044 0.7686 1.1436 0.0432 0.0289 0.0011附录6三维空间基本曲线误差分析程序附录6体积与质量关系参数估计程序v=（x1.*（x2.2）/（

15、12*pi）;p=polyfit（v,y,1）y1=polyval（p,v）;plot（v,y1,p =1.1587 42.4279附录6体积与质量关系误差分析程序v=3.826971324962187e+02,3.972660168319077e+02,3.972660168319077e+02,4.993836215124780e+02,.6.003714913553345e+02,6.003714913553345e+02,9.043809345407842e+02,1.209761232302733e+03;y1=1.1587*v-42.4279;附录6模型改进相关程序1.利用最小二乘法在matlab软件中求参数的程序如下： l=36.8,31.8,43.8,36.8,32.1,45.1,35.9,32.1; m=765,482,1162,737,482,1389,652,454; plot（l,m, c=log（l）; b=log（m）; p=polyfit（c,b,l）p=3.0265 -4.32572.用matlab画出拟合图形程序如下： p=polyfit（c,b,1） l1=31.8:0.5: m1=polyval（p,l1）; fplot（l（3.0265）,31.8 45.1）