血糖浓度论文

承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是：血液中葡萄糖浓度的稳定性分析所属学院：数学与计算机科学学院参赛队员(打印并签名)：陈志伟P091712725沈良P091712741金华P091712718谢琼辉P091712735田亚玲P091712724日期：2012年4月30日评阅编号(教师评阅时填写)：血液中葡萄糖浓度的稳定性问题摘要本文讨论人体注射葡萄糖溶液时，血液中葡萄糖浓度的稳定情况与人体血液容积的关系。首先，分析在人体注射葡萄糖溶液时血液中葡萄糖浓度的变化率dgdt与注射速率r、人体血液的容积v和血液中葡萄糖浓度g的关系，在人体血液容积不变的前提下建立关于葡萄糖浓度的微分方程模型。其次，考虑到人体的血液容积会因溶液的注射而增加，但又不能无限增加，利用常数变易法可解微分方程模型②的通解，但通解很难实现积分，故考虑实际情况，采用数值逼近法做进一步修改与完善，使其更符合实际。最后，通过对Logistic模型的分析与Mathematica软件的运用，得出人体血液中葡萄糖浓度在不同人体血液容积机理下的稳定情况。关键词微分方程模型常数变易法数值逼近法Logistic模型一、问题重述输液就是让某种液体以稳定的速率进入静脉的过程。当输入葡萄糖溶液时，血液中的葡萄糖浓度的增长率与注射速率成正比，与人体血液溶积成反比，而由于人体组织的吸收作用，血液中葡萄糖浓度的减少率与其本身成正比。根据是否考虑人体血液容积V的变化情况，可分以下3中情况，并在3种假设下建立模型，进而讨论稳定情况。问题1人体血液容积V不变；问题2V随着注入葡萄糖溶液而增加；问题3由于排泄等因素V的增加有极限值。二、问题分析讨论人体血液葡萄糖浓度的稳定情况时，由于血液中葡萄糖浓度的变化率等于浓度增长率与减少率之差，而葡萄糖浓度的增长率又与注射速率r成正比，与人体血液体积成反比；同时，考虑到人体组织的吸收作用，葡萄糖浓度的减少率与本身成正比。针对问题1利用血液中葡萄糖浓度的变化率等于浓度的增长率与减少率之差建立一阶线性微分方程，得到血液中葡萄糖浓度与时间的函数关系，根据自治微分方程的求法得到平衡点；针对问题2由注射速率的微分表达式，结合葡萄糖浓度随时间变化的微分表达式得到葡萄糖浓度与实践的函数关系，利用数学软件画出图像得到平衡点；针对问题3人体血液容积增加有极限值，符合阻滞增长模型，结合葡萄糖浓度随时间变化的微分表达式最终得到平衡点。三、基本假设1、假设注射为静脉注射，注射速率r恒定；2、假设人体所注射的溶液浓度不变；3、假设该人的各项生理机能正常；4、不考虑其他因素对人体内血液中葡萄糖浓度的影响.四、符号表示注射葡萄糖溶液之前，人体血液中的葡萄糖浓度；g人体血液中葡萄糖浓度；r注射速率；V人体血液容积；人体初始血液体积：p人体中葡萄糖浓度的增加率；q人体中葡萄糖浓度的减少率；1k关于g的增加率与r、v之间的关系(1ko)；2k关于g的减少率与g本身的常数(2ko).五、模型建立与求解利用人体血液中葡萄糖浓度的变化率等于浓度增长率与减少率之差，建立以下微分方程模型。㈠人体血液容积V不变设时刻t的葡萄糖浓度g(t)是连续可微函数，因为血液中葡萄糖浓度g(t)的增长率与注射速率r成正比，与人体血液容积V成反比，而由于人体组织的吸收作用，g(t)的减少率与g(t)本身成正比，则血液中葡萄糖浓度g的增长率p=Vrk1，g的减少率2qkg。考虑t到t+的增长率，就有g(t+)-g(t)=(p-q)=(Vrk1-gk2)设t=0时血液中葡萄糖浓度为，即得微分方程①令f(g)=Vrk1-gk2=0,解得g()=Vkrk21,即为平衡点则此方程的解为g(t)=Vkrk21+()tke2则有tlimg(t)=Vkrk21=g()所以即为①的稳定点。㈡人体容积V随着注入溶液而增加当人体容积V随着注入溶液而增加时，V不再是常量，而是一个随时间变化而变化的变量，则可设V=+rt根据模型①的建立过程，同理可以得dtdg=Vrk1-gk2因此，可建立人体血液容积V随着注入溶液而增加情况下的微分方程模型：②对于一阶线性微分方程求解，利用常数变易法，可得通解tktktkegdtertVrketg222001由于微分方程不是自治方程，所以没有平衡状态。㈢由于排泄等因素V的增加有极限值当V的增加有极限值时，符合Logistic模型。因此可建立由于排泄等因素人体血液容积V的增加有极限值情况下的微分方程模型：③根据阻滞增长模型可得解得可以建立微分方程由微分方程稳定性有从而解得平衡点同时，有综上所述，就是稳定点。六、模型结果分析与检验对于模型（1），可以画出图像来分析：对于模型（2），由于右端积分比较难以求出，因此如果像模型（1）一样来求解或对解进行讨论，势必会出现混乱。图2()gt随时间t在0至100分钟内的变化图形和许多微分方程一样，数值逼近看来是唯一可行的求解方法。为了能够用数学软件Mathematica来描绘出的图形，我们分别取单位（1/min），（mg/min），（cm3），1.0（cm3/mg），（mg/cm3），并假设人体含有4600ml的血液，输入的葡萄糖溶液为500ml，并假设输液过程进行了100min.在Mathematica软件中输入程序（程序见附录）得到如下()gt随时间t在0至100分钟内的变化图形：204060801000.250.50.7511.251.51.75)t(gt)t(gtg图1人体中葡萄糖溶液浓度g与时间t的关系12krgkv，即0dgdt时，g为增函数；12krgkv，即0dgdt时，g为减函数。vkrk21t对于模型（3），人体的容量是有限的，静脉中液体的体积不可能无限增长；另一种机理，如排尿，将使体积保持协调，这是人体生理功能的一种实际极限。七、模型的评价与改进模型评价本文研究了人体血液中葡萄糖浓度的稳定性问题，在静脉注射的前提条件下，讨论了人体血液中葡萄糖浓度在不同人体血液容积机理下的稳定情况，建立的三个微分方程模型中的参数可以通过参数估计法得到，能够推广并应用到更一般的实际事例中，不足之处在于模型的求解和参数估计方面很难做到精确。模型改进对于（1）式中两个参数1k，2k，可以用参数估计的方法确定参数的值，并在注入葡萄糖溶液之前测量血糖浓度，确定0g的值，然后在1t，2t时刻对人的血糖浓度进行测量，由方程2102,1,2iktikrggeikv可确定1k，2k的值，进而可以做出图像进行分析。八、模型推广由（1）式表示的微分方程模型，它描述了葡萄糖血液浓度随时间变化的规律，在药物动力学、农业种植和社会经济等方面也有着广泛的应用，因此可以推广到更一般的形式：dxabxdt（4）其中,0ab，初始条件为()oxtx，对于（4）式进行分离变量dxdtbx对上式两端求不定积分得11lnxtcb即btxce再由常数变易法得（4）式的解为btaxceb由初始条件可得00()btbcexa由此得（4）的解为0()()btbtobaxtexeab由（1）和（4）可知()xt具有如下性质：1.稳定点axb；2.当axb时，0dxdt，即表明()xt严格单调递增；当axb时，0dxdt，即表明()xt严格单调递减；3.00lim()bttaxtxeb.因此在药物动力学、农业种植和社会经济等方面，只要符合微分方程模型（1）和关系式（4）的情况，都具有上述3点性质，可以更好的为社会实践、生产服务。参考文献：[1]杨启帆，边馥萍，数学模型，浙江：浙江大学出版社，1990.[2]数学建模，西安：陕西师范大学出版社，2003.[3]姜启源，数学模型，北京：高等教育出版社，2005.[4]东北师范大学微分方程教研室，常微分方程，北京：高等教育出版社，2005.