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

数学模型课程设计专 业： 数学与应用数学 姓 名： * 学 号： * 指导教师： * 日 期： 2013 年 6 月 25日关于生物心率与体重之间的关系建立初等模型一问题提出：生物学家认为，对于休息状态的热血动物，消耗的能量主要用于维持体温，能量与从心脏到全身的血流量成正比，而体温主要通过身体表面散失，建立一个动物（体重：g）与心率单位：次/min之间关系的模型，并用下面的数据加以检验.动物 体重/g心率/（次 ）田鼠25670家鼠200420兔2 000205小狗5 000120大狗30 00085羊50 00070人70 00072马45000038二问题分析1.根据生活经验可知，心率大小的影响因素不止一两个，由于很多因素我们无法控制，如在运动状态下和休息状态下，心率的大小明显不同；而随着年龄的增长，心率也会改变为了简化问题，我们假设所讨论的结果是在休息的状态下，显而易见，这样的假设比较粗略。2.从问题的提出可以看到，对于热血动物来说，消耗能量与全身血流量成正比，体温从体表散失，于是有：体表散热率=心跳产热率。而动物消耗的能量主要用于维持体温，而体热量通过表面积散失，则有：能量动物的表面积，能量血流量3.建模的目的在于找出体重与心率之间的近似数量关系。显然，由上表可看出：体重越大，其心率就越小；反之，动物的体重越小，其心率就越大。三问题假设u 本次模型只考虑在外界温度不变动物休息状态下的、身体正常的情况下 的心率问题。u 假设生物的体积与长度的立方完全成正比，表面积与长度的平方也是完全成正比。 u 假设体表散热率等于心跳产热率u 假设每分钟总的血流量等于心率与每分钟心脏到全身的血流量的乘积u 假设供血能力与体重成正比u 假设每分钟通过血液提供的能量与每分钟总血流量成正比四 符号说明符号说明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*103.*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*103.*W.(-1/3);>> NN = 626.8088 313.4044 145.4694 107.1828 58.9850 49.7498 44.4716 23.9172动物 体重实际心率模型结果田鼠25670626.8088家鼠200420313.4044兔2 000205145.4694小狗5 000120107.1828大狗30 0008558.9850羊50 0007049.7498人70 0007244.4716马4500003823.9172由上表可看出，虽然曲线轨迹比较接近原数据，可是将模型结果与实际数据比较，差别还是很大的。此次模型是先通过比例关系推导出的模型：，然后为了简便，将模型简化为，最后又用最小二乘法求出比例系数K。而且，假设条件比较粗糙，又缺乏足够的数据支持，所以，得到的结果不够满意。八 优缺点与存在问题v 模型的优点：1) 建立模型时假设环境不变，减少影响结果的因素，使得问题简化，便于分析求解。2) 由于本课程是关于心率与体重的关系，所以二者之间的关系不能直接观测出，而建立模型使得问题迎刃而解。3) 利用matlab编程进行求解，所得结果数据比较准确，利用matlab画图利于比较实际数据与模型结果的差异。v 模型的缺点：1) 建立模型时，假设热血动物在休息状态时消耗的能量全部转换为热量，忽略了其他热量的散失方式。2) 假设外界环境不变，但不能确定外界环境是否一直保持不变。3) 假设动物表面积、体积、与长度成正比，是非常粗略的。九 模型推广通过对生物体重与心率这个模型的研究，虽然结果不太令人满意，但是二者之间的大概关系很容易看出来，体重越大，心率越小。在实际医疗问题中，通过控制体重来预防心脏病，预防心脏病的发病率。可以通过进一步优化，利用足够的数据，改进这个模型，使人们能更有效的认识心脏病的发生、发展规律和研究防治措施。十参考文献【1】启源，金星，叶俊.数学模型，第四版。高等教育【2】静，但琦。数学建模与数学实验，高等教育。【3】启源，金星，刑文训，立平。大学数学实验，第2版。清华大学【4】Richard L.Burden，J.Douglas Faires.数值分析，第七版。高等教育