数学建模-数值分析实验-给药方案的设计(共10页)

1、给药方案的设计 摘 要一种新药用于临床之前，必须设计给药方案。要设计给药方案,必须知道给药后血药浓度随时间变化的规律。众所周知，快速静脉注射后，血药浓度立即上升；然后迅速下降。当血药浓度太低时，达不到预期的治疗效果；当血药浓度太高，又可能导致药物中毒或副作用太强。临床上，每种药物有一个最小有效浓度和一个最大有效浓度。设计给药方案时，要使血药浓度保持在最小有效浓度和最大有效浓度之间。在本文中先根据时刻与血药浓度的关系，建立时刻与血药浓度之间的数学模型，然后分别用线性最小二乘拟合和非线性最小二乘拟合两种方法通过MATLB软件对给出的时刻及血药浓度相应数值进行曲线拟合，得到模型中各参数的值，从而确定

2、函数方程，然后设计给药方案，根据已知的最小有效浓度和最大有效浓度计算初次剂量和每次注射剂量的大小，以及间隔时间的长短，从而确定给药方案。另外在非线性最小二乘进行曲线拟合时，分别使用isqcurvefit函数和isqnonlin函数进行拟合，发现拟合效果一致。关键词：线性最小二乘拟合、非线性最小二乘拟合、isqcurvefit函数、isqnonlin函数问题重述临床上，每种药物有一个最小有效浓度c1和一个最大有效浓度c2。设计给药方案时，要使血药浓度 保持在c1c2之间。本题设c1=10，c2=25(ug/ml).要设计给药方案,必须知道给药后血药浓度随时间变化的规律。从实验和理论两方面着手：在

3、实验方面,对某人用快速静脉注射方式一次注入该药物300mg后,在一定时刻t(小时)采集血药,测得血药浓度c(ug/ml)如下表:t(h)0.250.511.523468c(ug/ml)19.2118.1515.3614.1012.899.327.455.243.01问题分析要想设计出给药方案，就要得知时刻与血药浓度的关系，根据题目中给出的时刻及血药浓度相应数值数据，可以用两种方法，一是线性最小二乘拟合，二是非线性最小二乘拟合，通过利用MATLB软件进行曲线拟合，找到时刻与血药浓度的函数方程，从而确定给药方案。模型假设1、 假设机体是一个房室，室内血药浓度均匀（即一室模型）。2、 假设排除速率和

4、血药浓度成正比，比例系数为k,并且k>0.3、 假设血液容积为v，t=0时的注射剂量为d，则t=0时的血药浓度为d/v。符号说明c 血药浓度t 时刻d 初始时的注射剂量v 血液容积k 排除速率和血药浓度的比例系数D 每次注射剂量 初次剂量T 间隔时间 最小有效浓度 最大有效浓度模型建立与求解由假设2可得 (1)由假设3可得 (2)由（1）（2）式可得： (3)方法一：用线性最小二乘拟合c(t)对式(3)两边同时取对数得: (4)用MATLAB做线性最小二乘拟合得：k=0.2347，v=15.0219。程序如下：clear all;close all;clct=0.25 0.5 1 1.5

5、 2 3 4 6 8;c=19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01;y=log(c);p=polyfit(t,y,1); C=polyval(p,t); C=exp(C);v=exp(log(300)-p(2); k=-p(1); plot(t,c,'or',t,C,'b')从而拟合曲线： （5）见下图：图1 线性最小二乘拟合曲线由，求得：即：第一次注射375.5475mg,以后每隔3.9041h注射225.3285mg。方法二：用非线性最小二乘拟合1、利用MATLAB的isqcurvefit函数做非线性

6、最小二乘拟合可得：k=0.2420，v=14.8212。误差平方和为norm=1.0659。MATLAB程序如下：function f=fun(x,xdata)f=(300/x(1)*exp(-xdata*x(2);endclear all;close all;clc;clear all;close all;clc;t=0.25 0.5 1 1.5 2 3 4 6 8;c=19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01;c0=1 1;x,norm,res=lsqcurvefit('fun',c0,t,c)f=fun(x,t);

7、plot(t,c,'or',t,f,'b')从而拟合曲线： （6）拟合曲线见下图：图 2非线性最小二乘拟合曲线求得：即：第一次注射370.5300mg,以后每隔3.7863h注射222.3180mg。2、利用MATLAB的isqnonlin函数做非线性最小二乘拟合知：k=0.2420，v=14.8212。误差平方和为norm=1.0659。MATLAB程序如下：function f=fun(x,xdata,c)f=(300/x(1)*exp(-xdata*x(2)-c;endclear all;close all;clc;t=0.25 0.5 1 1.5 2 3 4 6 8;c=19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01;c0=1 1;x,norm,res=ls