Matlab微分方程应用2.pptVIP

油画真假辨别 放射物质衰变原理：物质的放射性正比于现存物质的原子数 衰变史： 鉴别油画的方法： 模型建立： 是否现代膺品的判别 其它鉴定结果： * 药物在体内的分布 药物进入机体形成血药浓度(单位体积血液的药物量) 血药浓度需保持在一定范围内——给药方案设计 药物在体内吸收、分布和排除过程 ——药物动力学 建立房室模型——药物动力学的基本步骤 何为房室系统？ 在用微分方程研究实际问题时，人们常常采用一种叫“房室系统”的观点来考察问题。根据研究对象的特征或研究的不同精度要求，我们把研究对象看成一个整体（单房室系统）或将其剖分成若干个相互存在着某种联系的部分（多房室系统）。 房室具有以下特征：它由考察对象均匀分布而成，房室中考察对象的数量或浓度（密度）的变化率与外部环境有关，这种关系被称为“交换”且交换满足着总量守衡。以下我们将用房室系统的方法来研究药物在体内的分布。 交换 环境 内部 单房室系统 均匀分布 房室——机体的一部分，药物在一个房室内均匀分布(血药浓度为常数)，在房室间按一定规律转移 药物的分解与排泄（输出）速率通常被认为是与药物当前的浓度成正比的，即： 药物分布的单房室模型 单房室模型是最简单的模型，它假设：体内药物在任一时刻都是均匀分布的，设t时刻体内药物的总量为x(t)；系统处于一种动态平衡中，即成立着关系式： 药物的输入规律与给药的方式有关。下面，我们来研究一下在几种 常见的给药方式下体内药体的变化 规律。 机体 环境 药物总量 图3-8 假设药物均匀分布 情况1 快速静脉注射 机体 环境 只输出不输入房室 其解为： 药物的浓度： 与放射性物质类似，医学上将血浆药物浓度衰减一半所需的时间称为药物的血浆半衰期： 负增长率的Malthus模型 在快速静脉注射时，总量为D的药物在瞬间被注入体内。设机体的体积为V，则我们可以近似地将系统看成初始总量为D，浓度为D/V，只输出不输入的房室，即系统可看成近似地满足微分方程： 情况2 恒速静脉点滴 机体 环境 恒定速率输入房室 药物以恒速点滴方式进入体内，即: 则体内药物总量满足： （x(0)=0） 这是一个一阶常系数线性方程，其解为： 或 易见： 称为稳态血药浓度 对于多次点滴，设点滴时间为T1，两次点滴之间的间隔时间设为T2，则在第一次点滴结束时病人体内的药物浓度可由上式得出。其后T2时间内为情况1。故： 0≤t≤T1 T1≤t≤T1 +T2 类似可讨论以后各次点滴时的情况，区别只在初值上的不同。第二次点滴起，患者 体内的初始药物浓度不为零。 情况3 口服药或肌注 y(t) x(t) K1y K1x 环境 机体 外部药物 口服药或肌肉注射时，药物的吸收方式与点滴时不同，药物虽然瞬间进入了体内，但它一般都集中与身体的某一部位，靠其表面与肌体接触而逐步被吸收。设药物被吸收的速率与存量药物的数量成正比，记比例系数为K1，即若记t时刻残留药物量为y(t)，则y满足： D为口服或肌注药物总量 因而： 所以 解得： 从而药物浓度 图3-9给出了上述三种情况下体内血药浓度的变化曲线。容易看出，快速静脉注射能使血药浓度立即达到峰值，常用于急救等紧急情况；口服、肌注与点滴也有一定的差异，主要表现在血药浓度的峰值出现在不同的时刻，血药的有效浓度保持时间也不尽相同。 图3-9 我们已求得三种常见给药方式下的血药浓度C(t)，当然也容易求得血药浓度的峰值及出现峰值的时间，因而，也不难根据不同疾病的治疗要求找出最佳治疗方案。 上述研究是将机体看成一个均匀分布的同质单元，故被称单房室模型，但机体事实上并不是这样。药物进入血液，通过血液循环药物被带到身体的各个部位，又通过交换进入各个器官。因此，要建立更接近实际情况的数学模型就必须正视机体部位之间的差异及相互之间的关联关系，这就需要多房室系统模型。 I II k12 k21 两房室系统 图3-10表示的是一种常见的两房室模型，其间的k12表示由室I渗透到室II的变化率前的系数，而k21则表示由室II返回室I的变化率前的系数，它们刻划了两室间的内在联系，其值应当用实验测定，使之尽可能地接近实际情况。 当差异较大的部分较多时，可以类似建立多房室系统，即N房室系统，包括：中心室(心、肺、肾等)和周边室(四肢、肌肉等) 历史背景： 二战后,荷兰保安机关开始搜捕纳粹分子的合作者,于1945年5月以通敌罪逮捕了一名三流画家H.A.Vanmeegren,此人曾将荷兰著名画家创作的一批名贵油画盗卖给德国.但此人被捕后称自己从未出卖过荷