生物医学数据分析课件医用高等数学实例分析.pptVIP

微分方程 【例3.5】细菌繁殖模型：在理想环境中，细菌的繁殖率与细菌的数目成正比。若 时细菌的数目为 ，求细菌的繁殖规律。 function dxdt=func3f3(t,x) dxdt=0.5*x [T，Y] =ode45(@func3f3,[0,10],10) 式中的D表示 第三十一页，共四十六页。 微分方程 【例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章医用数学基础及其MATLAB实现 知识结构 第一页，共四十六页。 第3章医用数学基础及其MATLAB实现 The Tower of Babel 沟通的重要性 第二页，共四十六页。 第3章医用数学基础及其MATLAB实现 医用高等数学实例分析 函数分析 微分方程 第三页，共四十六页。 定义：设在某个变化过程中存在两个变量x和y，若对于变量x的每一个允许取的值，变量y按照某一对应关系，都有唯一确定的值与之对应，则称变量y为变量x的函数，记为 其中：x为自变量，y为因变量，f表示二者间的对应函数关系，D为x的定义域，即允许取值的范围。 函数分析 函数 第四页，共四十六页。 【例3.1】 静脉注射G钠盐100000单位后，血清中的药物浓度C为时间t的函数C(t)： 其中时间t的单位为小时(h)，C的单位为单位每毫升(单位/mL)。试绘制C随时间t变化的曲线。 函数分析 第五页，共四十六页。 函数分析 函数文件的建立 %% func3f1.m function y=func3f1(x) if (x=0)(x0.25) y=14.5*x; end if (x=0.25)(x1) y=4.66-4.26*x; end if (x=1)(x=3) y=0.6-0.2*x; end if (x3)||(x0) error(自变量不在该函数定义域内取值) end %% end func3f1.m 第六页，共四十六页。 函数分析 函数文件的调用 直接调用 输出参数＝函数名(输入参数) y=func3f1(2) y=0.2 间接调用 函数句柄是用于间接调用一个函数的Matlab值或数据类型。在调用其它函数时可以传递函数句柄，也可在数据结构中保存函数句柄备用。 通过命令形式 fhandle = @functionname 可以创建函数句柄 trigFun=@sin，或匿名函数sqr = @(x) x.^2;。 trigFun(1) y=feval(@func3f1,2) 第七页，共四十六页。 函数分析 %%%exam31.m %%方法1 t=0:0.05:3; for k=1:length(t) C(k)=feval(@fun1,t(k));nd figure,plot(t,C,o-), title(血清中药物浓度随时间变化图); xlabel(t/h); ylabel(C/unit/ml); 第八页，共四十六页。 函数分析 %%%exam31.m %%方法2 figure,fplot(@func3f1,[0 3]); title(血清中药物浓度随时间变化图) xlabel(t/h);ylabel(C/unit/ml); 第九页，共四十六页。 函数分析 常用重要的函数分析命令 result = func(x0) %求函数值 fplot(func,[xl xd]) %在区间[xl xd]作图 xsolve = fzero(func,x0) %求x0附近的零点 Xmin = fminbnd(func,xl,xd) %求区间[xl xd]内极小值 ezplot %函数作图 第十页，共四十六页。 【例3.2】已知口服一定剂量的药物后，血药浓度C与时间t的关