关于生物心率与体重之间的关系建立初等模型

1、数学模型课程设计业:数学与应用数学姓 名:*学 号:指导教师：*日 期：2013年6月25 0关于生物心率与体重之间的关系建立初等模型一. 问题提出：生物学家认为，对于休息状态的热血动物，消耗的能量主要用于维持体温， 能量与从心脏到全身的血流量成正比，而体温主要通过身体表面散失，建立一个 动物（体重：g）与心率单位：次/min之间关系的模型，并用下面的数据加以检 验.动物体重/g心率/ （次 min"）皿鼠25670家鼠200420兔2 000205小狗5 000120大狗30 00085羊50 00070人70 00072马45000038二. 问题分析1根据生活经验可知，心率大小

2、的影响因素不止一两个，由于很多因素我们 无法控制，如在运动状态下和休息状态下，心率的大小明显不同；而随着年龄的 增长，心率也会改变为了简化问题，我们假设所讨论的结果是在休息的状态 下，显而易见，这样的假设比较粗略。2从问题的提出可以看到，对于热血动物来说，消耗能量与全身血流量成正 比，体温从体表散失，于是有：体表散热率二心跳产热率。而动物消耗的能量主 要用于维持体温，而体内热量通过表面积散失，则有：能量*动物的表面积，能 量*血流量3.建模的目的在于找出体重与心率之间的近似数量关系。显然，由上表可看 出：体重越大，其心率就越小；反z,动物的体重越小，其心率就越大。三. 问题假设木次模型只考虑在

3、外界温度不变动物休息状态下的、身体正常的情况下 的心率问题。假设生物的体积与长度的立方完全成正比，表面积与长度的平方也是完 全成正比。假设体表散热率等于心跳产热率假设每分钟总的血流量等于心率与每分钟心脏到全身的血流量的乘积假设供血能力与体重成正比假设每分钟通过血液提供的能量与每分钟总血流量成正比符号说明符号说明s生物身体的表面积l生物的长度v生物的体积w生物的体重q每分钟为维持体温而流失的能量e每分钟通过血液提供的能量t每分钟总血流量n生物的心率q每次心跳从心脏到全身的血流量k表面积s与长度l的比例系数k2体积v与长度l的比例系数k3体积v与体重w的比例系数能量q与表面积s的比例系数k.q与w

4、, n的比例系数能量e与总血流量t的比例系数k心率q与血流量w的比例系数五. 模型建立根据以往的数学知识可知：表面积soc长度l的平方，体积voc长度l的立 方.即：s = k、c v = k2l（其中仏怠为常数且大于0）。于是，有s = v3生物的体积v与体重w成正比，即v = k3w,将其代入式可得：折=一沪k魏（d 为维持体温而流失的能量q与动物的表面积有关，且是一种正比关系，则有 w = k4s （其中心为常数且大于0） 将式代入式，可得出：q二叭,昕 （哄）1由假设体表散热率等于心跳产热率，即0=，则有：e = /雪叨 （來3 ）3由于假设供血能力与体重成正比，则有： q = kv

5、（其中心为常数且大于0）由于假设每分钟总的血流量等于心率与每分钟心脏到全身的学流量的乘积,则有：t=nq,由式,可得：t = k5nw 由于假设每分钟通过血液提供的能量与每分钟总血流量成正比，则有： e = k6t（其中心为常数且大于0） 将式与式代入式，可得：= k6k5nw «p： （mx k's加 （丛严式就是通过比例关系所推导得出的数学模型。六. 模型求解要利用由w来计算n,则首先得确定其中的参数你込、心、心、%、k6o 由于要确定这六个参数需要进行测量以及查阅相关资料，而迫于时间限制，可以 丄k k用另一-种简便的方法：将式改写成n = kw (其中，k=),这样,

6、(丛严5我们只需要确定参数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。已给数据的散点图)» xlabelc生物的体重')；» ylabel('生物的心率j；已给数据的散点图700600500400300200100000.511.522.533.544.5生物的体重x 1q5由图可看出这些数据组成的散点图是一个曲

7、线，其数学模型是n = kw 若对 两边取对数，则有：lg(/t) = lg(a：)-|lg(w),令 y=lg(n),a=lg(k),& = -|,由此可 得模型为：y = ax,是一个直线，可用最小二乘法拟合出a的值。下面用mat lab 3画出改造数据后的散点图。» w=25,200,2000,5000,30000,50000,70000,450000;» n二670,420,205,120,85,70,72,38;» x=loglo(w);» y=logl0(n);» plot(x,y,'r*');»

8、titled改造数据后的散点图)改造数据后的散点图32.51 511.522.533.544.55由此可以看出，与前面的分析符合，改造后的数据大致与一个直线拟合，由之前 学过的数值分析，可直接利用最小二乘法的结果求出a的值。若要用最小二乘法 拟合出数据为(xpx)i的逼近直线歹=筑+ax，贝山(m /=!丿 /=1 丿mm加加*=i=lz=l/=1所以，此次模型所要求得的:加mmmz/=!/=1» w=25,200,2000,5000,30000,50000,70000,450000;» n=670,420,205,120,85,70,72,381;» x=log

9、lo(w);» y=logl0(n);» a=(sum(x.a2)*sum(y)-sum(x.*y)*sum(x)/(8*(sum(x.a2)-sum(x)a2)3.2631» plot(x,y,r*',x,a-l/3*x,'g,)32.82.62.42.221.81.61.411.522.533.544.555.56所以:k = 10" =1.8328x10'，h =1.8328xlo3vv_5 下面就根据所求出的模型作图，与题目所给的数据比较。» w=25,200,2000,5000,30000,50000,7000

10、0,450000;» n=670,420,205,120,85,70,72,38;» n=1.8328*10a3.*w.a(-1/3)；» plot(w,n,ro,w,n,b,)» xlabelc生物的体重')；» ylabef生物的心率*);» title。拟合后的曲线与原数据的比较)拟合后的曲线与原数据的比较700600500400300200100000.511.522.533.544.5生物的体重x 105七. 结论解释由上面的图可看出，第5, 6, 7这儿个点拟合的不够好，下面是将生物的体重代_1入模型n = 1.8

11、328x1o3w'3与实际数据作比较得出的表格» w=25,200,2000,5000,30000,50000,70000,450000;» n=670,420,205,120,85,70,72,38;» n= 1.8328* 10a3.*w.a(-l/3);»nn =626.8088313.4044145.4694107.182858.985049.749844.471623.9172动物体重实际心率模型结果田鼠25670626.8088家鼠200420313.4044兔2 000205145.4694小狗5 000120107.1828人狗3

12、0 0008558.9850羊50 0007049.7498人70 0007244.4716马4500003823.9172由上表可看岀，虽然曲线轨迹比较接近原数据，可是将模型结果与实际数据比较，差别还是很大的。此次模型是先通过比例关系推导出的模型：2 坐 加,然后为了简 （丛严x便，将模型简化为n = kw'最后又用最小二乘法求出比例系数ko而且，假 设条件比较粗糙，又缺乏足够的数据支持，所以，得到的结果不够满意。八. 优缺点与存在问题模型的优点：1）建立模型时假设环境不变，减少影响结果的因素，使得问题简化，便于分析 求解。2）由于本课程是关于心率与体重的关系，所以二者之间的关系不能直接观测 出，而建立模型使得问题迎刃而解。3）利用matlab编程进行求解，所得结果数据比较准确，利用matlab画图利于 比较实际数据与模型结果的差异。模型的缺点：1）建立模型时，假设热血动物在休息状态时消耗的能量全部转换为热量，忽略 了其他热量的散失方式。2）假设外界环境不变，但不能确定外界环境是否一直保持不变。3）假设动物表面积、体积、与长度成正比，是非常粗略的。九. 模型推广通过对牛物体重与心率这个模型的研究，虽然结果不太令人满意，但是二者 之间的大概关系很容易看岀来，体重越人，心率越小。在实际医疗问题中，通过 控制体重來预防心脏病，预防心脏病的发病率。可以通过