生物医学数据分析解析课程--医用高等数学实例分析解析资料文档.pptVIP

第3章医用数学基础及其MATLAB实现 知识结构 第3章医用数学基础及其MATLAB实现 第3章医用数学基础及其MATLAB实现 医用高等数学实例分析 函数分析 微分方程 函数分析 函数 函数分析 函数分析 函数文件的建立 函数分析 函数文件的调用 函数分析 函数分析 函数分析 常用重要的函数分析命令 函数分析 函数分析 函数分析 函数分析 多元函数分析 函数分析 函数分析 函数分析 函数拟合 函数分析 最小二乘拟合(Least Squares Fitting) 函数分析 最小二乘拟合(Least Squares Fitting) 函数分析 最小二乘拟合(Least Squares Fitting) 函数分析 函数分析 关于函数拟合 函数分析 关于函数拟合 函数分析 关于函数拟合 微分方程 微分方程数学模型 微分方程 微分方程 微分方程求解 微分方程 微分方程求解 微分方程 微分方程求解 微分方程 例如：Van der Pol Equation 微分方程 微分方程 微分方程 微分方程 微分方程 微分方程 微分方程 微分方程 微分方程 微分方程 微分方程 微分方程 微分方程 微分方程 微分方程 【例3.5】 当k=0.5时 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 【例3.6】细菌繁殖的非理想环境模型：考虑除系统本身的繁殖外有的细菌向系统外迁移或死亡，损耗速率是时间t的线性函数，即At+B，系统内细菌的繁殖率与细菌的数目成正比，若 时细菌的数目为 ，求系统内的细菌繁殖规律。 解：设t时刻的细菌数目为 ，则根据题意有 【例3.6】 解析解： 数值解： function dxdt=func3f4(t,x,k,A,B) dxdt=k*x-(A*t+B); 应掌握带调用参数的函数的微分方程解法 【例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),tspan,x0);%%k=0.5,A=1.8,B=1; figure, plot(t1,x_1n,r-o,t2,x_2n,b-+,t3,x_3n,k*-); legend(k=0.5,A=2.0,B=1,k=0.5,A=2.2,B=1,k=0.5,A=1.8,B=1); title(细菌繁殖模型的解(非理想环境)); xlabel(时间t) ylabel(细菌数目); 【例3.6】 图 不同参数对细菌繁殖系统行为的影响 【例3.7】细菌自然生长模型(Logistic增长模型)：t时刻的细菌数目 由细菌的繁殖率m和死亡率n决定，实践证明： 若 时细菌的数目为 ，求系统内的细菌繁殖规律。 解：由题意可知，单位时间内细菌的净增加数目应为当前的细菌数目乘以(繁殖率－死亡率),即 【例3.7】 解析解： 其各阶导数的求法 syms r k t x0 f=r*x0/(k*x0-(k*x0-r)*exp(-r*t)); f1=diff(f,t)?;%%一阶导数 f2=diff(f1,t)?; %%二阶导数 【例3.7】 tspan=[0,10]; x0=2; r=1; k=0.2; t=linspace(tspan(1),tspan(2),100); x1=r*x0./(k*x0-(k*x0-r).*exp(-r.*t)); r=0.7; x2=r*x0./(k*x0-(k*x0-r).*exp(-r.*t)); r=0.4; x3=r*x0./(k*x0-(k*x0-r).*exp(-r.*t)); r=0.01; x4=r*x0./(k*x0-(k*x0-r).*exp(-r.*t)); jj=find(t5); jj=jj(1);%%找到t=5