数学建模_血液中葡萄糖浓度的稳定性问题

血液中葡萄糖浓度的稳定性问题

摘要本文讨论人体注射葡萄糖溶液时，血液中葡萄糖浓度的稳定情况与人体血液容积的关系。

首先，分析在人体注射葡萄糖溶液时血液中葡萄糖浓度的变化率dg

dt

与注射速率r、人

体血液的容积v和血液中葡萄糖浓度g的关系，在人体血液容积不变的前提下建立关于葡萄糖浓度的微分方程模型。

其次，考虑到人体的血液容积会因溶液的注射而增加，但又不能无限增加，利用常数变易法可解微分方程模型（2）的通解，但通解很难实现积分，故考虑实际情况，采用数值逼近法做进一步修改与完善，使其更符合实际。

最后，通过对Logistic模型的分析与Mathematica软件的运用，得出人体血液中葡萄糖浓度在不同人体血液容积机理下的稳定性。

一、问题重述

输液就是让某种液体以稳定的速率进入静脉的过程。当输入葡萄糖溶液时，血液中的葡萄糖浓度的增长率与注射速率成正比，与人体血液溶积成反比，而由于人体组织的吸收作用，血液中葡萄糖浓度的减少率与其本身成正比。根据是否考虑人体血液容积V的变化情况，可分以下3中情况，并在3种假设下建立模型，进而讨论稳定情况。

问题1 人体血液容积V不变；

问题2 V随着注入葡萄糖溶液而增加；

问题3 由于排泄等因素V的增加有极限值。

二、问题分析

讨论人体血液葡萄糖浓度的稳定情况时，由于血液中葡萄糖浓度的变化率等于浓度增长率与减少率之差，而葡萄糖浓度的增长率又与注射速率r成正比，与人体血液体积成反比；同时，考虑到人体组织的吸收作用，葡萄糖浓度的减少率与本身成正比。

针对问题1 利用血液中葡萄糖浓度的变化率等于浓度的增长率与减少率之差建立一阶线性微分方程，得到血液中葡萄糖浓度与时间的函数关系，根据自治微分方程的求法得到平衡点；

针对问题2 由注射速率的微分表达式，结合葡萄糖浓度随时间变化的微分表达式得到葡萄糖浓度与实践的函数关系，利用数学软件画出图像得到平衡点；

针对问题3 人体血液容积增加有极限值，符合阻滞增长模型，结合葡萄糖浓度随时间变化的微分表达式最终得到平衡点。

三、基本假设

1、假设注射为静脉注射，注射速率r恒定；

2、假设人体所注射的溶液浓度不变；

3、假设该人的各项生理机能正常；

4、不考虑其他因素对人体内血液中葡萄糖浓度的影响.

四、符号表示

0g 注射葡萄糖溶液之前，人体血液中的葡萄糖浓度； g 人体血液中葡萄糖浓度； r 注射速率； v 人体血液容积；

p 人体中葡萄糖浓度的增加率； q 人体中葡萄糖浓度的减少率；

1k 关于g 的增加率与r 、v 之间的关系(1k o >)； 2k 关于g 的减少率与g 本身的常数 (2k o >).

五、模型建立与求解

利用人体血液中葡萄糖浓度的变化率等于浓度增长率与减少率之差，建立以下微分方程模型。

5.1 人体血液容积V 不变

分析人体血液容积不变的情况，利用一阶线性非齐次方程的求解方法得到人体血液中

葡萄糖浓度随时间变化的函数。先表示血液中葡萄糖浓度g 的增长率1k r

P v

= 由于人体组织的吸收作用，则g 的减少率 2q k g =

注射葡萄糖溶液之前，人体血液中的葡萄糖浓度0g 为微分方程的初值条件，即

0(0)g g =

因此，可建立人体血液容积V 不变情况下的微分方程模型：

(1)

令错误！未找到引用源。 得平衡点 错误！未找到引用源。 而满足微分方程的初值解为

错误！未找到引用源。

又有 错误！未找到引用源。 所以 错误！未找到引用源。 即为(1)的稳定点。 5.2 人体容积V 随着注入溶液而增加

当人体血液容积V 变化不能忽略时，不再是常量，则可设

0v v rt =+

根据5.1中微分方程模型（1）的建立过程，同理可以得

12dg k r k g dt v

=-

因此，可建立人体血液容积V随着注入溶液而增加情况下的微分方程模型：

(2) 对于一阶线性微分方程求解，利用常数变易法，可得通解

由于微分方程不是自治方程，所以没有平衡状态。

5.3 由于排泄等因素V的增加有极限值

当V的增加有极限值时，符合Logistic模型。因此可建立由于排泄等因素人体血液容积V的增加有极限值情况下的微分方程模型：

错误！未找到引用源。(3)

根据阻滞增长模型可得

解得错误！未找到引用源。

可以建立微分方程

由微分方程稳定性有

错误！未找到引用源。

从而解得平衡点错误！未找到引用源。

同时，有

综上所述，错误！未找到引用源。就是稳定点。

六、模型结果分析与检验

对于模型（1），可以画出图像来分析：

对于模型（2），由于右端积分比较难以求出，因此如果像模型（1）一样来求解或

对解进行讨论，势必会出现混乱。

图2 ()g t 随时间t 在0至100分钟内的变化图形

和许多微分方程一样，数值逼近看来是唯一可行的求解方法。 对于模型（3），人体的容量是有限的，静脉中液体的体积不可能无限增长；另一种机理，如排尿，将使体积保持协调，这是人体生理功能的一种实际极限。

七、模型的评价与改进

g

图1 人体中葡萄糖溶液浓度g 与时间t 的关系

12k r g k v <，即0dg dt >时，g 为增函数； 12k r g k v >

，即0dg dt

<时，g 为减函数。

v

k r k 21

t

为了能够用数学软件Mathematica 来描绘出的图形，我们分别取单位（1/min ），（mg/min ）， （cm3 ），1.0（cm3 /mg ），（mg/ cm3），并假设人体含有4600ml 的血液，输入的葡萄糖溶液为500ml ，并假设输液过程进行了100min.在Mathematica 软件中输入程序（程序见附录）得到如下()g t 随时间t 在0至

100分钟内的变化图形：

20

40

60

80

0.25

0.50.7511.25

1.51.75)

t (g t

)t (g t