【北京交通大学《数学模型》电子教材 】4.1 颅内压与血流速度的关系问题（精品推荐 免费分享）.doc

第4章 随机模型 4.1 颅内压与血流速度的关系问题 1)?????? 问题的提出 正常人颅内有一定的压力，称为颅内压。颅高压是神经科的急危重症，是神经内外科病人死亡的第一原因。因此颅内压值的测定对临床治疗的决策起至关重要的作用。目前，颅内压的定量测定均使用创伤性的方法。创伤性方法对患者的身体具有损害，且当患者出现急性颅内压增高而需要尽快获得实时的颅压值时，创伤性的方法显然不适宜。 大量动物实验及临床脑血管造影观察表明脑血循环和颅内压的改变关系密切。下表一数据来自北京市人民医院神经外科研究人员对5只兔子的实验数据，其中ICP代表颅内压；V代表脑血流速度。试根据此数据建立兔颅内压增高与脑血流速度的关系的数学模型。 表一 急性ICP增高与对应的V值 0.98 832 1.96 770 3.92 617 4.90 268 ICP(KPa) V(Hz/s) 8.13 0 2）问题的分析与假设 颅内压与很多因素有关，血流速度是其中最重要的一个，简化起见这里只考虑与血流速度的关系。本问题只给出了有限的实验数据点，因此通常可以想到用插值、拟合或回归的方法来求出颅内压与血流速度的关系。考虑到颅内压与血流速度的随机因素，这里选用回归的方法来确定颅内压与血流速度的关系。注意到由于医学上已有的关于颅内压和血流速度的定性关系较为丰富，建立基于微分方程的参数辨识模型也是可能的。因为所给数据由实验得到，难免有误差，这里认为它们是基本反映事实的，而不考虑实验动物或实验仪器本身不正常的特殊情况。 3）模型的建立与求解 将表一中的数据作图，横轴代表颅内压，纵轴代表血流速度，可以得到下面这张图 ?????????? 回归模型 根据回归分析的使用经验可知，回归分析常用于这样的情况：虽然变量Y与变量X之间有一定的关系，但这种关系与通常的函数关系不同，Y的值不能由X的值精确确定。这里的颅内压与血流速度的关系即属于这种情形。事实上影响或反映颅内压变化的因素很多，血流速度只是其中之一，只不过一般医学界认为它是一种主要因素，因此考虑它们二者的关系用回归分析是可行的。 设血流速度为V，颅内压为P，如果用直线回归方法，则有回归式如下： V =+ ( (1) 其中，( 是随机误差，其数学期望为零。确定未知系数、的常用方法是最小二乘法，计算结果如下： = 987．343 ， = -123.163 它的方差分析见下面表二： 表二 线性回归方差分析表 ─────────────────────────────────────- 均差平方和 自由度 均 方 F P ─────────────────────────────────────- 回归变异 472563.8634 1 472563.8634 50.5589 0.005723 误差 28040.33656 3 4515.9830 总变异 500603.2000 4 ───────────────────────────────────── 利用Methematica 4.0 得到如下图形，图中的五个点是已知点: 观察图形可知，直线回归效果并不是很好。根据散点图尝试采用指数回归，可以获得如下回归式： （2） 经计算得到回归方程 其方差分析表为： 表三 指数回归方差分析表 ─────────────────────────────────────- 均差平方和 自由度 均 方 F P ─────────────────────────────────────- 回归变异 487089.38115 1 487089.38115 108.131 0.0018975 误差