关于生物心率与体重之间的关系建立初等模型.docx
《关于生物心率与体重之间的关系建立初等模型.docx》由会员分享，可在线阅读，更多相关《关于生物心率与体重之间的关系建立初等模型.docx（11页珍藏版）》请在冰豆网上搜索。
关于生物心率与体重之间的关系建立初等模型
《数学模型》课程设计
专业：
数学与应用数学
姓名：
****
学号：
******
指导教师：
****
日期：
2013年6月25日
关于生物心率与体重之间的关系建立初等模型
一．问题提出：
生物学家认为，对于休息状态的热血动物，消耗的能量主要用于维持体温，能量与从心脏到全身的血流量成正比，而体温主要通过身体表面散失，建立一个动物（体重：
g）与心率单位：
次/min之间关系的模型，并用下面的数据加以检验.
动物
体重/g
心率/（次•
）
田鼠
25
670
家鼠
200
420
兔
2000
205
小狗
5000
120
大狗
30000
85
羊
50000
70
人
70000
72
马
450000
38
二．问题分析
1.根据生活经验可知，心率大小的影响因素不止一两个，由于很多因素我们无法控制，如在运动状态下和休息状态下，心率的大小明显不同；而随着年龄的增长，心率也会改变……为了简化问题，我们假设所讨论的结果是在休息的状态下，显而易见，这样的假设比较粗略。
2.从问题的提出可以看到，对于热血动物来说，消耗能量与全身血流量成正比，体温从体表散失，于是有：
体表散热率=心跳产热率。
而动物消耗的能量主要用于维持体温，而体内热量通过表面积散失，则有：
能量∝动物的表面积，能量∝血流量
3.建模的目的在于找出体重与心率之间的近似数量关系。
显然，由上表可看出：
体重越大，其心率就越小；反之，动物的体重越小，其心率就越大。
三．问题假设
◆本次模型只考虑在外界温度不变动物休息状态下的、身体正常的情况下的心率问题。
◆假设生物的体积与长度的立方完全成正比，表面积与长度的平方也是完全成正比。
◆假设体表散热率等于心跳产热率
◆假设每分钟总的血流量等于心率与每分钟心脏到全身的血流量的乘积
◆假设供血能力与体重成正比
◆假设每分钟通过血液提供的能量与每分钟总血流量成正比
四．符号说明
符号
说明
S
生物身体的表面积
L
生物的长度
V
生物的体积
W
生物的体重
Q
每分钟为维持体温而流失的能量
E
每分钟通过血液提供的能量
T
每分钟总血流量
n
生物的心率
q
每次心跳从心脏到全身的血流量
表面积S与长度L的比例系数
体积V与长度L的比例系数
体积V与体重W的比例系数
能量Q与表面积S的比例系数
q与W，n的比例系数
能量E与总血流量T的比例系数
K
心率q与血流量W的比例系数
五．模型建立
✧根据以往的数学知识可知：
表面积S∝长度L的平方，体积V∝长度L的立方.即：
。
于是，有
………………………………………………①
✧生物的体积V与体重W成正比，即
，将其代入①式可得：
……………………………………②
✧为维持体温而流失的能量Q与动物的表面积有关，且是一种正比关系，则有
………………………………③
将③式代入②式，可得出：
…………………………………………………④
✧由假设体表散热率等于心跳产热率，即Q=E，则有：
…………………………………………………⑤
✧由于假设供血能力与体重成正比，则有：
……………………………⑥
✧由于假设每分钟总的血流量等于心率与每分钟心脏到全身的学流量的乘积，则有：
T=nq,由⑥式，可得：
…………………………………………………⑦
✧由于假设每分钟通过血液提供的能量与每分钟总血流量成正比，则有：
……………………………⑧
将⑤式与⑦式代入⑧式，可得：
即：
………………………………………………⑨
⑨式就是通过比例关系所推导得出的数学模型。
六．模型求解
要利用⑨由W来计算n，则首先得确定其中的参数
。
由于要确定这六个参数需要进行测量以及查阅相关资料，而迫于时间限制，可以用另一种简便的方法：
将⑨式改写成
（其中，
），这样，我们只需要确定参数K就可以得出n与W之间的比较粗糙的关系。
>>W=[25,200,2000,5000,30000,50000,70000,450000];
>>n=[670,420,205,120,85,70,72,38];
>>plot（W,n,'ro'）;
>>title（'已给数据的散点图'）
>>xlabel（'生物的体重'）;
>>ylabel（'生物的心率'）;
由图可看出这些数据组成的散点图是一个曲线，其数学模型是
，若对两边取对数，则有：
，令y=lg（n）,a=lg（K）,
由此可得模型为：
，是一个直线，可用最小二乘法拟合出a的值。
下面用matlab画出改造数据后的散点图。
>>W=[25,200,2000,5000,30000,50000,70000,450000];
>>n=[670,420,205,120,85,70,72,38];
>>x=log10（W）;
>>y=log10（n）;
>>plot（x,y,'r*'）;
>>title（'改造数据后的散点图'）
由此可以看出，与前面的分析符合，改造后的数据大致与一个直线拟合，由之前学过的数值分析，可直接利用最小二乘法的结果求出a的值。
若要用最小二乘法拟合出数据
的逼近直
，则：
所以，此次模型所要求得的：
>>W=[25,200,2000,5000,30000,50000,70000,450000];
>>n=[670,420,205,120,85,70,72,38];
>>x=log10（W）;
>>y=log10（n）;
>>a=（sum（x.^2）*sum（y）-sum（x.*y）*sum（x））/（8*（sum（x.^2））-sum（x）^2）
a=
3.2631
>>plot（x,y,'r*',x,a-1/3*x,'g'）
所以：
，
下面就根据所求出的模型作图，与题目所给的数据比较。
>>W=[25,200,2000,5000,30000,50000,70000,450000];
>>n=[670,420,205,120,85,70,72,38];
>>N=1.8328*10^3.*W.^（-1/3）;
>>plot（W,n,'ro',W,N,'b'）
>>xlabel（'生物的体重'）;
>>ylabel（'生物的心率'）;
>>title（'拟合后的曲线与原数据的比较'）
七．结论解释
由上面的图可看出，第5，6，7这几个点拟合的不够好，下面是将生物的体重代入模型
与实际数据作比较得出的表格
>>W=[25,200,2000,5000,30000,50000,70000,450000];
>>n=[670,420,205,120,85,70,72,38];
>>N=1.8328*10^3.*W.^（-1/3）;
>>N
N=
626.8088313.4044145.4694107.182858.985049.749844.471623.9172
动物
体重
实际心率
模型结果
田鼠
25
670
626.8088
家鼠
200
420
313.4044
兔
2000
205
145.4694
小狗
5000
120
107.1828
大狗
30000
85
58.9850
羊
50000
70
49.7498
人
70000
72
44.4716
马
450000
38
23.9172
由上表可看出，虽然曲线轨迹比较接近原数据，可是将模型结果与实际数据比较，差别还是很大的。
此次模型是先通过比例关系推导出的模型：
，然后为了简便，将模型简化为
，最后又用最小二乘法求出比例系数K。
而且，假设条件比较粗糙，又缺乏足够的数据支持，所以，得到的结果不够满意。
八．优缺点与存在问题
v模型的优点：
1）建立模型时假设环境不变，减少影响结果的因素，使得问题简化，便于分析求解。
2）由于本课程是关于心率与体重的关系，所以二者之间的关系不能直接观测出，而建立模型使得问题迎刃而解。
3）利用matlab编程进行求解，所得结果数据比较准确，利用matlab画图利于比较实际数据与模型结果的差异。
v模型的缺点：
1）建立模型时，假设热血动物在休息状态时消耗的能量全部转换为热量，忽略了其他热量的散失方式。
2）假设外界环境不变，但不能确定外界环境是否一直保持不变。
3）假设动物表面积、体积、与长度成正比，是非常粗略的。
九．模型推广
通过对生物体重与心率这个模型的研究，虽然结果不太令人满意，但是二者之间的大概关系很容易看出来，体重越大，心率越小。
在实际医疗问题中，通过控制体重来预防心脏病，预防心脏病的发病率。
可以通过进一步优化，利用足够的数据，改进这个模型，使人们能更有效的认识心脏病的发生、发展规律和研究防治措施。
十．参考文献
【1】姜启源，谢金星，叶俊.数学模型，第四版。
高等教育出版社
【2】赵静，但琦。
《数学建模与数学实验》，高等教育出版社。
【3】姜启源，谢金星，刑文训，张立平。
《大学数学实验》，第2版。
清华大学出版社
【4】RichardL.Burden，J.DouglasFaires.数值分析，第七版。
高等教育出版社
升级会员 