局部脑血流的测定

局部脑血流的测定

摘要

随着医学的快速发展，脑部的研究越来越得到人们的重视。而对脑血流系数的精确测定可以帮助人们快速得到相关指标，所以对医学领域来说很有意义。

问题1，首先我们根据题设中的信息及放射性元素衰减性，列出头部记数率和呼出记数率关于脑部血流系数的微分表达式；再分析试验得到的呼出气记数率的数据，用Matlab 的拟合工具箱拟合得到关于呼出气记数率的函数；最后代入上述微分表达式求解微分方程即可得到关于脑部血流系数的表达式，表达式为)(5

.11000)(5.1Kt t

e e K k t N ----=

。

问题2，先用Matlab 拟合工具箱求出头部记数率的函数，用对比系数法

可得到K 和k 的近似值，即3977.0,5015

.0≈≈k K ，但这种方法是不精确的，只是用于后面方法得到参数的验证；再利用最小二乘法求解拟合后曲线参数的函数以及由问题1得到的关于脑部血流系数的表达式，即可得到脑血流系数：0.5000=K ，4001.0=k ，最后对得到的值进行误差分析，可知脑血流的预测值和实际值很吻合，比较符合题意。

关键词：对比系数法 放射性元素衰减 曲线拟合

一问题重述

用放射性同位素测定大脑局部脑血流量的方法如下：由受试者吸入含有某种放射性同位素的气体，然后将探测器置于受试者头部某固定处，定时测量该处的放射性记数率(简称记数率)，同时测量他呼出气的记数率。

由于动脉血将肺部的放射性同位素传送至大脑，使脑部同位素增加，而脑血流又将同位素带离，使同位素减少。实验证明由脑血流引起局部地区记数率下降的速率与当时该处的记数率成正比，其比例系数反映了该处的脑血流量，被称为脑血流量系数，只要确定该系数即可推算出脑血流量。动脉血从肺输送同位素至大脑引起脑部记数率上升的速度与当时呼出气的记数率成正比。若某受试者的测试数据如附表1所示：

根据以上题目所给的条件及数据，回答以下问题：

1、建立确定脑部血流系数的数学模型；

2、计算上述受试者的脑血流系数。

二模型假设

1.脑部记数率的上升只与从肺部输送的放射性同位素有关；

2.脑部记数率的下降只与当时该处的脑血流量有关；

3.脑血流量在测定期间恒定，心脏博动、被测试者大脑活动、情感波动

等带来的变化忽略不予考虑；

4.每次仪器测量为相互独立事件，各测量值无记忆相关；

5.在吸入气体瞬时，脑中放射物记数率为零；

6.脑血流量与脑血流量系数成单值函数关系，求得后者即可确定前者。

三符号说明

四 问题分析

问题1，首先根据题设可知：一方面，由脑血流引起局部地区记数率下降的速率与当时该处的记数率成正比, 并且其比例系数反映了该处的脑血流量；另一方面，动脉血从肺输送同位素至大脑引起脑部记数率上升的速率与当时呼出气的记数率成正比，由上述两方面可得到头部记数率关于脑部血流系数和呼出气记数率的表达式；再分析试验得到的呼出气记数率的数据，用MATLAB 的拟合工具箱拟合得到关于呼出气记数率的函数；最后代入上述表达式求解微分方程即可得到关于脑部血流系数的表达式。

问题2，利用Matlab 里专门求解拟合后曲线参数的函数，再由问题1得到的关于脑部血流系数的表达式即可得到脑血流系数。最后再对得到的值进行误差分析。

五 模型的建立与求解

5.1求解脑血流系数的数学模型

设某时刻0≥t 时，脑部记数率为)(t N ，在t ∆时刻后记数率为)(t t N ∆+，

由题设及基本假设1和2可知，脑部记数率的增量)()(t N t t N N -∆+=∆只与下面两个因素有关：

(1)动脉血从肺部输送放射性同位素至大脑引起脑部记数率的增量为

1N ∆；

(2)脑血流将放射性同位素带离使得脑部记数率的减量为2N ∆。 (3)根据文献[1]，放射性元素自身有衰减，设其半衰期为τ，由此引起的

记数率下降为∆3N

又由医学实验和假定有：

)()(2

1

t KN dt

dN t kP dt

dN ==， 又

311(t)ln 2(t)2t

N N t τ

οτ∆⎛⎫∆=-⨯⨯⨯∆+∆ ⎪⎝⎭

3ln 2

(t)dN N dt τ

=- 所以考虑△: 时刻内头部放射性元素记数率变化, 有

123(t)N (t)(t)(t)N N N ∆=∆-∆+∆

其变化率为：

3

12(t)(t)dN dN dN dN dt dt dt dt

=-+ 于是得到

3

ln 2

()()dN kP t KN t N dt τ=-- (1)

由于在测试时放射性同位素的半衰期一般很大, 这样会给测量和试验带来严重影响,

因此假定: τ→∞于是( l) 变为 ()()dN

kP t KN t dt =- （2） 分析式(2)，要确定脑血流系数的模型，必须分析)(t P 和)(t N 的实验数据，观察其变化趋势。首先用Matlab 绘出)(t P 和)(t N 的散点图并观察其变化

趋势，)(t P 和)(t N 的散点图如下：

t

P (t )

图5-1 呼出气记数率的散点图

123

45

678910

t

N (t )

图5-2 脑部记数率的散点图

由上图可知，t 和)(t P 有近似于bt

ae 的关系，而t 和)(t N 的关系暂时不能

直接观察出，设bt

ae t P =)(，用MATLAB 的拟合工具箱可得到参数a =1000，

=b -1.5，而其相关系数为1，说明拟合的非常精确。拟合后的图像如下图所示：

1

2345678910

图5-2 拟合函数)(t P 的图像

由基本假设5，即0)0(=N 和式(1)联立可得一带有初值的微分方程：

⎪⎩

⎪⎨⎧=-=-0)0()

(10005.1N t KN ke dt

dN

t (3) 解此微分方程得到

)(5

.11000)(5.1Kt t

e e K k t N ----= (4)

5.2 求解脑血流系数的算法模型

以下为求解)(t N 的算法模型：

算法模型I ：再分析)(t N 的试验数据，同样用Matlab 的拟合工具箱里现有的指数函数可精确拟合出)(t N 的函数，表达式为：

)(t N =t t e e 479.15015.039834032--⨯-⨯，可近似认为5.1479.1,39834032≈≈，则对比系数可得到3977.0,5015.0≈≈k K ，但这种方法是不精确的，只是用于后面