- 1、本文档共7页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数学模型综合实验鱼模型
开课学院、实验室:数统学院 实验时间 : 2016年3月13日
课程
名称 数学实验 实验项目
名 称 MATLAB软件入门之绘图与程序设计 实验项目类型 验证 演示 综合 设计 其他 指导
教师 成 绩 实验目的
通过对一个鱼类重量的问题的分析,使学生:
1、了解可以用高等数学的基础知识来解决这种类型的重要问题;
2、体验利用高等数学的思想、方法来分析和解决实际问题;
3、接受写作研究报告的初步训练;
4、激发学生学习高等数学的兴趣。
实验内容
1、应用场景
2、问题分析
3、建立数学模型
一、问题重述
出于保护的目的,垂钓俱乐部鼓励会员钓到鱼之后马上放生。该俱乐部还希望根据钓到的鱼的总重量给予奖励。垂钓者怎样确定所钓到的鱼的重量呢?如果采用便携称,一方面,称起来不方便,另一方面对于小鱼并不准确。请根据某个容易量测的量来估计鱼的重量。
在大奖赛期间收集到得数据为:
长度
14.5
12.5
17.25
14.5
12.625
17.75
14.125
12.625
重量
27
17
41
26
17
48
23
16
如果考虑腰围应该怎么做?可用如下数据来做测试检验。
长度/cm
36.8
31.8
43.8
36.8
32.1
45.1
35.9
32.1
重量/g
765
482
1162
737
482
1389
652
454
腰围/cm
24.8
21.3
27.9
24.8
21.6
31.8
22.9
21.6
二、问题分析
通过测量鱼的长度估计鱼的重量,我们不能陷入对鱼复杂生理结构的研究,否则将复杂化,得不到有使用价值的模型,经过大量的分析研究,我们利用类比方法以及鱼的体形都是相似的,找出小鱼的身长和体重两者关系,或是大鱼的身长,胸围和体重三者的比例关系,进而建立了测量鱼重量的模型。
进行以下模型假设:
1、池塘里的鱼体型都是相似的;
2、每条鱼被钓上的几率是相等的;
3、鱼肉的密度是相等的:
4、不区别鱼的雌雄且鱼的肥瘦均匀;
5、鱼的腰围指鱼身的最大周长;
6、池塘的鱼都是在同一条件下生长;
三、数学模型的建立与求解
1、只有身长时,我们将鱼看作一个矩形,建立身长三次方与重量的关系函数,有身长和腰围时,将鱼看作是上下两个圆锥拼接而成,建立函数。输入数据与程序,描绘出一个散点图
图1 鱼身长与重量一次拟合
图2 鱼身长,腰围与重量一次拟合
2、先建立一个一次拟合查看效果,如图2,图3
3、借助cftool工具箱,反复尝试比对一次二次及以上的拟合函数效果
四、实验结果及分析
经过cftool中的反复比较,主要比较了R-square这个值,发现在只有身长时,三次关系最符合,在有身长和腰围时,二次关系最符合。
所以实验结果如图3,图4所示。
当只有身长时,我们根据(1)式计算:
f(x) = -0.1261*x^3 + 11.59*x^2 + -208*x + 2840 (1)
当有身长和腰围时,我们根据(2)式计算:
f(x) = 0.01265*x^2 + 9.146*x + 7820 (2)
这个数据虽然符合的很好,但估计有一定巧合因素在其中,计算也甚为复杂,所以如果是在现实中使用,也可以考虑一次关系。
图3 鱼身长和重量三次关系函数
图4 鱼身长,腰围和重量二次关系图
五、附录(程序等)
身长和重量的拟合求解程序:
L1=[14.5,12.5,17.25,14.5,12.625,17.75,14.125,12.625];%长度(cm)
W1=[27,17,41,26,17,48,23,16];%重量(g)
X1=L1.^3;%身长的三次方记为X
P1=polyfit(X1,W1,1);
Q1=polyval(P1,X1);
plot(X1,W1,o,X1,Q1,b);
title(Model One);
身长,腰围,重量的拟合求解程序:
L2=[36.8,31.8,43.8,36.8,32.1,45.1,35.9,32.1];% 长度(cm)
W2=[765,482,1162,737,482,1389,652,454];%重量(g)
B2=[24.8,21.3,27.9,24.8,21.6,31.8,22.9,21.6];%腰围(cm)
X2=L2.*(B2.^2);%
文档评论(0)