医学数据分析及matlab2012-2精选.ppt

讲解人：刘青萍 生物医学数据分析 刘青萍 3.1.2 微分方程在生物医学中的应用 代谢过程 药物动力学工程 细胞动力学工程 生物系统种群 4 1 2 3 微分方程是表示未知函数的导数以及 自变量之间关系的方程。 Click to edit title style 1 2 3 细菌的繁殖 药物动力学模型 流行病数学模型 微分方程数学模型 (1) 寻找改变量，通常是描述方程文字中的变化率 (微商、单位增加量或单位减少量)。 (2) 对问题中的特征进行数学刻画。 (3) 用微元法建立微分方程。 (4) 确定微分方程的定解条件(初、边值条件)。 (5) 求解或讨论方程(数值解或定性理论)。 (6) 模型和结果的讨论与分析。 数学模型的建立流程 例题 y’=y y=c*exp(x) y’’=y y=c1*sinx+c2cosx 微分方程求解 一步法，用于一阶微分方程 取两点斜率 的平均 斜率来计算，精度高 Euler 欧拉法 Runge Kutta 龙格-库塔法 求解 Matlab求解微分方程 dsolve ode45 简单的微分方程或微分方程组，用于线性常系数微分方程求解—符号解 是（4,5）阶的龙格-库塔-芬尔格方法，用于常微分方程的数值求解。 dsolve命令 ode45命令 常用 方法 dsolve命令 dsolve(‘equation’,‘condition’) equation表示微分方程式 Dy表示y’，D2y表示y’’，condition表示初始条件 dsolve(‘equation’) dsolve(‘equation’,‘condition’) 给出微分方程初值问题的解，表示为t的函数 dsolve(‘equation’,‘v’) 给出微分方程的解析解，表示为v的函数 dsolve(‘equation’,‘condition’,‘v’) 例题 计算微分方程dy/dx+3xy=xe-x2的通解。 湖南中医药大学 dsolve(Dy+3*x*y=x*exp(-x^2)) ans = (C2/exp(3*t*x) + 1)/(3*exp(x^2)) dsolve(Dy+3*x*y=x*exp(-x^2),x) ans = C6/exp((3*x^2)/2) + exp(x^2/2)/exp((3*x^2)/2) ode命令 [T,Y]=ode45(‘F’,TSPAN,Y0) F定义此微分方程的形式y’=F(t,y)，函数 F(t,y)应返回一列向量 TSPAN=[T0 TFANAL]表示此微分方程的积分限是从T0到TFANAL，也可以是一些离散的点 Y0表示初始条件 fplot命令 fplot(‘函数名’,lims,’s’) 绘制指定函数的曲线 limits是一个指定x轴范围的向量[xmin xmax]或者是x轴和y轴范围的向量[xmin xmax ymin ymax]。 s可指定线形 例题 例3.5 细菌繁殖模型 在理想环境中，细菌的繁殖率与细菌的数目成正比。若t=0时细菌的数目为x(0)=x0,求细菌的繁殖规律 设t时刻的细菌数目为x(t)，则根据题意有： 设x(t)=u，du=kudt du/u=kdt ∫1/udu= ∫kdt 可求出 ln|u|=kt ，|u|=ekt 又k0，x(0)=x0，即x(t)=x0*ekt dsolve方法求解 dsolve(Dx=k*x,x(0)=x0) ans = x0*exp(k*t) ode方法求解 k=0.5; tspan=[0,10]; x0=10; %%%解析解画图 fplot(@(t) 10*exp(0.5*t),tspan,k-); %%%数值解 [t,x_n]=ode45(@func3f3,tspan,x0); hold on,plot(t,x_n,ro) legend(解析解,数值解); title(细菌繁殖模型的解); xlabel(时间t) ylabel(细菌数目); hold off dsolve(Dx=k*x,x(0)=x0) ans = x0*exp(k*t) 思考题1 按照例3.5的思路自学例3.6 k=0.5; tspan=[0,10]; x0=10; %%数值解 [t1,x_1n]=ode45(@(t,x)func3f4(t,x,0.5,2,1),tspan,x0);??%%k=0.5,A=2,B=1; [t2,x_2n]=ode45(@(t,x)func3f4(t,x,0.5,2.2,1),tspan,x0);%%k=0.5,A=2.2,B=1; [t3,x_3n]=ode45(@(t,x)func3f4(t,x,0.5,1.8,1