局部脑血流地测定数学建模.pdf

局部脑血流的测定 一 . 问题简介 脑血流量是诊断和治疗脑梗塞 ,脑出血 ,动脉瘤和先天性动 ,静脉血管畸形等脑血管疾病 的主要依据。测定脑血流量可为研究人脑在不同的病理和生理条件下的功能提供客观指标， 它对研究脑循环药物的药理作用也很有帮助。 所以人们长期致力于寻找有效地测定脑血流量 的方法。 近年来出现了以放射性同位素作示踪剂测定人脑局部血流量的方法。这种方法大致可描 述如下 :由受试者吸入某种放射性同位素的气体，然后将探测器置于受试者头部某固定处， 定时测量该处放射性同位素的计数率 (简称计数率 ) ，同时测量他呼出气的计数率。 由于动脉血将肺部的放射性同位素输送至大脑， 使脑部同位素增加， 而脑血流又将同位 素带离， 使同位素减少， 实验证明由脑血流引起局部地区计数率下降的速率与当时该处的记 数率成正比， 其比例系数反映了该处的脑血流量， 被称为脑血流量系数， 只要确定该系数即 可推算出脑血流量。 动脉血从肺输送同位素至大脑引起脑部计数率上升的速率与当时呼出气 的计数率成正比。 若某受试者的测试数据如下 : 时 间(分) 1.00 1.25 1.50 1.75 2.00 2.25 2.50 2.75 3.00 3.25 3.50 3.75 头部记数率 1534 1528 1468 1378 1272 1162 1052 947 848 757 674 599 呼出气记数率 2231 1534 1054 724 498 342 235 162 111 76 52 36 时 间(分) 4.00 4.25 4.50 4.75 5.00 5.25 5.50 5.75 6.00 6.25 6.50 6.75 头部记数率 531 471 417 369 326 288 255 225 199 175 155 137 呼出气记数率 25 17 12 8 6 4 3 2 1 1 1 1 时 间(分) 7.00 7.25 7.50 7.75 8.00 8.25 8.50 8.75 9.00 9.25 9.50 9.75 10.0 头部记数率 121 107 94 83 73 65 57 50 44 39 35 31 27 呼出气记数率 0 0 0 0 0 0 0 0 0 0 0 0 0 试建立确定脑血流系数的数学模型并计算上述受试者的脑血流系数。 备注 :该题目是上海市 (1990 年 ) 大学生数学建模竞赛 A 题。 二 . 模型的假定 1. 脑部计数率 (记为 ht()) 的上升只与肺部的放射性同位素有关，上升速度与呼出气的记 数率 (记为 p()t ) 成正比，比例系数记为 k ； 2. 脑部记数率 ht() 的下降只与该处脑血流量有关，其下降速度正比于 ht() ，比例系数为 脑血流系数，记为 K ，这里忽略了放射性元素的衰变和其它因素； 3. 脑血流量在测定期间恒定，心脏博动，被测试者大脑活动，情感波动等带来的变化 忽略不予考虑； 4. 每次仪器测量为相互独立事件，各测量值无记忆相关； 5. 放射性同位素在人体内传递是从吸入气体 (含有放射物 )开始的， 并假定一次吸入， 因 此认为同位素在肺中瞬时达到最大浓度； 6. 在吸入气体瞬时，脑中放射物记数率为零； 7. 脑血流量与脑血流系数 K 成单值函数关系，求得后者即可确定前者。 三 . 模型的建立与分析 由于已知脑局部同位素的增减与已测定的头部记数率 ht() 和呼出气记数率 p()t 成正比 关系 ,于是很自然地确定以脑部同位素量 ,即脑部记数率作为讨论对象 . 。 1. 原始模型的建立 设某时刻 t ≥ 0 时