数学建模之减肥问题的数学模型

东北大学秦皇岛分校数学模型课程设计报告减肥问题的数学建模学院数学与统计学院专业信息与计算科学学号5133117姓名楚文玉指导教师张尚国刘超成绩教师评语：指导教师签字：2016年01月09日数学与统计学院课程设计（实习）报告第1页摘要肥胖已成为公众日益关注的卫生健康问题.肥胖是与目前严重危害人类健康疾病，如糖尿病、高血压、冠心病、血脂异常、胆囊疾病、痛风、骨关节病、阻塞性睡眠呼吸暂停、某些癌症等的发病有明显相关的危险因素之一.但是实际情况却是人们不会理性的对待自己的身体状况，经常使用一些不健康的方式减肥，到最后适得其反，给自己的身体造成很大的伤害.本文特别的从数学模型的角度来考虑和认识问题，通过该模型的建立，科学的解释了肥胖的机理，引导群众合理科学的减肥.本文建立了减肥的数学模型，从数学的角度对有关身体肥胖的规律做进一步的探讨和分析.在研究此问题时，体重的实时变化数据是我们研究的核心数据，这就会使我们联系到变量的变化率或导数，这样所得到变量之间的关系式就是微分方模型.微分方程模型反映的是变量之间的间接关系，因此，在研究体重，能量与运动之间的关系时，得到直接关系就得求解微分方程.本文利用了微分方程模型求解减肥的实际问题，根据基本规律写出了平衡关系式[()()][()()]tttDABRtt再利用一定的转换条件进行转化为简单明了的式子，求解出模型关系式()(1)dtdtateed然后根据建立的模型表达式来解决一些实际的减肥问题，给出数学模型所能解答的一些实际建议.关键字：微分方程模型能量守恒能量转换系数1问题重述1.1课题的背景随着社会的进步和发展，人们的生活水平在不断提高，饮食营养摄入量的改善和变化、生活方式的改变，使得肥胖成了社会关注的一个问题.为此，联合国世界卫生组织曾颁布人体体重指数（简记BMI）：体重（单位：kg）除以身高（单位：m）的平方，规定BMI在18.5至25为正常，大于25为超重，超过30则为肥胖.据悉我国有关机构针对东方人的特点，拟将上述规定中的25改为24，30改为29.无论从健康的角度，是数学与统计学院课程设计（实习）报告第2页从审美的角度，人们越来越重视减肥，大量的减肥机构和商品出现，不少自感肥胖的人加入了减肥的行列，盲目的减肥，使得人们感到不理想，如何对待减肥问题，不妨通过组建模型，从数学的角度，对有关的规律作一些探讨和分析.根据背景知识，我们知道任何人通过饮食摄取的能量不能低于用于维持人体正常生理功能所需要的能量，因此作为人体体重极限值的减肥效果指标一定存在一个下限1，当1*时表明能量的摄入过低并致使维持他本人正常的生理功能的所需不足，这种减肥所得到的结果不能认为是有效的，它将危机人的身体健康，是危险的，称1为减肥的临界指标.另外，人们认为减肥所采取的各种体力运动对能量的消耗也有一个所能承受的范围，记为10RR，当能量的摄取量高于体重0时，这是体重不会从0减少，所以可以看到单一的措施达不到减肥效果.1.2具体的问题和相关数据现有五个人，身高、体重和BMI指数分别如下表1.1所示，体重长期不变，试为他们按照以下方式制定减肥计划，使其体重减至自己的理想目标，并维持下去：表1.1身高，体重和BMI指数表人数编号12345身高1.71.681.641.721.71体重100112113114124BMI34.633.535.234.835.6理想目标7580808590每天摄入能量28572543273426892776题目具体要求如下：（1）在基本不运动的情况下安排计划，每天吸收的热量保持下限，减肥达到目标；（2）若是加快进程，增加运动，重新安排计划，经过查找资料得到以下各项运动每小时每kg体重的消耗的热量如下表1.2所示：表1.2每小时每kg体重的热量消耗运动跑步跳舞乒乓自行车（中速）游泳（50m/min）热量消耗7.03.04.42.57.9（3）给出达到目标后维持体重的方案.数学与统计学院课程设计（实习）报告第3页2模型假设与符号说明2.1问题分析本问题要建立减肥的数学模型，减肥是一个比较长期和不定的过程，因此要用数学的方法对减肥这一问题建模，就需要选定一个测量肥胖的标准量.因为人体的脂肪是能量的主要贮存和提供的方式，而且也是减肥的主要目标.因此，我们以人体脂肪的重量作为体重的标志.已知脂肪的能量转换率为100﹪，每千克脂肪可以转换为8000kcal，称D为脂肪的能量转换系数.肥胖主要是体现在人的身体上，减肥其实就是将人的体重降下来，所以归根到底，研究减肥就是要研究体重的变化，因此在减肥过程中我们要对人的体重进行持续的检测，忽略个体间的差异（年龄、性别、健康状况等）对减肥的影响，可以将人体的体重看成是时间t的函数()t.在减肥的过程中，无论是由于进食摄取能量导致体重的增加，还是由于体力活动消耗能量致使体重的减少，异或还有其他一些不可预知的因素，这都是一个渐变的过程,所以认定()t是连续光滑的.所以我们认为能量的摄取和消耗都是随时发生的，而不同的活动对能量的消耗是不同的.所以我们在建模的过程中需要设定一个参数用来表示某种活动消耗的人体能量.记r为某一种活动每小时所消耗的能量，记b为1kg体重每小时所消耗的能量.2.2模型假设1．假设以人体脂肪的重量作为体重的标志.2．假设体重随时间的变化()t是连续而且充分光滑的.3．假设在单位时间人体的能量消耗与其体重成正比.4．假设人体每天摄入的能量是一定的.记为A.5．正常代谢引起的减少正比于体重，每人每千克体重消耗热量一般为28.75~45.71kcal，且因人而异.6．假设在研究减肥的过程中，我们忽略个体间的差异对减肥的影响.7．人体每天摄入量是一定的，为了安全和健康，每天吸收热量不要小于1429kcal.8．假设单位时间内人体由于基础代谢和食物特殊动力作用所消耗的能量正比于人的体重.数学与统计学院课程设计（实习）报告第4页2.3符号说明D：脂肪的能量转化系数.()t：人体的体重关于时间t的的函数..r：每千克体重每小时运动所消耗的能量(/)/kcalkgh.b：每千克体重每小时所消耗的能量(/)/kcalkgh.0A：每天摄入的能量.1W：五个人理想的体重目标向量.A：五个人每天分别摄入的能量..W：五个人减肥前的体重.B：每人每天每千克体重基础代谢的能量消耗.3模型建立与求解3.1一般模型建立如果以1天为时间的计量单位，于是每天基础代谢的能量消耗量应=24(/)Bbkcald，由于人的某种运动一般不会是全天候的，不妨假设每天运动h小时，则每天由于运动所消耗的能量应为=(/)Rrhkcald.按照假设2，体重随时间的变化()t是连续而且充分光滑的，我们可以在任何一个时间段内考虑由于能量的摄入与消耗引起人的体重的变化.按照能量的平衡原理，任何时间段内由于体重的改变所引起的人体内能量的变化应该等于这段时间内摄入的能量与消耗的能量之差.我们选取某一段时间(,)ttt，在时间段(,)ttt内考虑能量的改变:设体重改变的能量变化为W,则有=[(+)()]WtttD（3.1）设摄入与消耗的能量之差为M，则有[()()]MABRtt（3.2）根据能量平衡原理有MW（3.3）得：[()()][()()]tttDABRtt（3.4）取0t，可得数学与统计学院课程设计（实习）报告第5页0dd(0)adt=（3.5）其中/aAD，()/dBRD，0t（模型开始考察时刻），即减肥问题的数学模型模型求解得()(1)dtdtateed（3.6）/aAD表示由于能量的摄入而增加的体重，而()/dBRD表示由于能量的消耗而失掉的百分数（每单位体重中由于基础代谢和活动而消耗掉的那部分）.3.2针对实际问题的模型建立1.由一般模型的建立已经知道减肥问题的数学模型为微分方程模型（3.6），利用此方法可求解出每个人要达到自己的理想体重的天数.首先确定此人每天每千克体重基础代谢的能量消耗B，因为没有运动，所以有0R，根据式（3.6）式，得ABW（3.7）从而得到每人每天每千克体重基础代谢的能量消耗从假设5可知，这些人普遍属于代谢消耗相当弱的人，加上吃得比较多，有没有运动，所以会长胖，进一步，由()t(五人的理想体重)，W（五人减肥前的体重），D=8000kcal/kg（脂肪的能量转换系数），根据式（3.6）式有001/lnln/adDBAtdadBBA（3.8）将A（五个人每天分别摄入的能量）的值代入上式时，就会得出五个人要达到自己的理想体重时的天数，如下表3.1所示表3.1达到理想体重所需天数表人12345天数194372313266298Matlab源程序：R=0;D=8000;%能量转换系数数学与统计学院课程设计（实习）报告第6页W1=[7580808590];%理想的体重目标A=[28572543273426892776];%每人每天摄入的能量W=[100112113114124];%每人的体重n=length(W);B=A./W%每人每天每千克体重基础代谢的能量消耗a=A./Dd=(B+R)./Dfori=1:nt(i)=-(D/B(i))*log((W1(i)*B(i)-A0)/(W(i)*B(i)-A0));%减肥所需要的时间end2.为加快进程，增加运动，结合查找资料得到各项运动每小时每kg体重消耗的热量表2，再结合假设3，取1hh，Rrhr，根据式（4.6）有001/()lnln/()adDBRAtdadBRBRA（3.9）将A（五个人每天分别摄入的能量）的值代入上式时，取不同的r，得到一组数据，在运动的情况下，我们选取的是一个小时，得到了每个人在不同运动强度下，要达到自己的理想目标所需的天数，如下表3.2所示：表3.2不同运动强度下达到理想体重所需天数运动跑步跳舞乒乓自行车游泳时间/天122155141160116187261229274176173232207243164148198177206140163220196230154Matlab源程序：h=1;r=[7.03.04.42.57.9];R=h.*r;n1=length(R);D=8000;%能量转换系数W1=[7580808590];%理想的体重目标数学与统计学院课程设计（实习）报告第7页A=[28572543273426892776];%每人每天摄入的能量W=[100112113114124];%每人的体重n=length(W);B=A./W;%每人每天每千克体重基础代谢的能量消耗forj=1:n1fori=1:nt=(i,j)=-(D./(B(i)+R(j))*log((W1(i).*(B(i)+R(j))-A0)./(W(i).*(B(i)+R(j))-A0)));%减肥所需要的时间endend3.要使体重稳定在一个定值，则有*ABR（3.10）根据自己的不同理想目标和B（每人每天每千克体重基础代谢的能量消耗），在不同小时下的能量消耗表：（1）在1h的情况下运动所消耗的能量，如下表3.3表3.31h的情况下运动所消耗的能量运动跑步跳舞乒乓自行车游泳消耗能量(kcal)2667.002367.8002472.8002330.2002735.3002376.4002056.4002168.4002016.4002448.4002495.6002175.6002287.6002135.6002567.6002600.0002260.0002379.0002217.5002676.5002644.8002284.8002410.8002239.8002725.800（2）在2h的情况下运动所消耗的能量，如下表3.4表3.42h的情况下运动所消耗的能量运动跑步跳舞乒乓自行车游泳消耗能量(kcal)3198.002592.8002802.8002517.