医用高等数学在MATLAB中的实现.docxVIP

评分 实验目的 Matlab提供了强大的数据分析功能，本次实验旨在熟悉和掌握应用 Matlab实现 函数分析、微分方程求解的功能。 二、 实验环境 1、 硬件配置：计算机 2、 软件环境：MATLAB R2010a 三、 实验原理 函数分析 自变量为单个变量时，称为一元函数，当自变量为多个变量时，称为多元 函数，表示为：y=f(x1，…;xn),x1 D1,…;xn Dn 在进行函数分析的时候，将每个函数编写成 M函数文件即可。 函数拟合 Matlab提供了一个使用函数polyfit,用于最小二乘拟合如下形式的线性多项 式： y(x)=p1*x^ n+p2*x^( n-1)+ +pn *x+p( n+1) 用法：p=polyfit(x,y,n). 其中x,y为输入的自变量和因变量向量，n为待拟合多项式的最高次数，返 回值p为包含n+1个以降幕顺序排列的各阶洗漱的行向量。 另外MATLAB中提供了专门用于函数拟合的工具 cftool，可以拟合各种函数 形式，只要在MATLAB命令窗口输入cftool,就会弹出一个曲线拟合工具界 面界面窗口。 微分方程在生物医学中的应用 Z. 《 《MATLAB 基础》实验报告 第 第 PAGE #页 共6页 (1)对于简单的微分方程或微分方程组，可以用 symbolic math toolbox里德 dsolve命令来获得解析解。 Dx(t)/dt=kx(t) 使用dsolve命令如下： f=dsolve(Dx=k*x,(0)=x0)(始终的 D 表示 d/dt) (2)大部分长微分方程，可以使用 ode45命令来进行述职的求解。 Dx(t)/dt=kx(t) 使用 ode45 命令如下：［T Y］=ode45(odefun.tepan,y0) 其中，odefun表示待求解的微分方程。Tspan表示述职求解的时间范围， 如［0,10］表示0到10秒。Y0表示待求标量的小护士值。返回值 T为数值求 解时间范围内的一系列采样点，y为对应时间点的待求变量值。 常用统计量函数 算数均值函数 mean，几何均值函数geomean,中位数函数median,极差函数 range,四分位数间距的函数iqr，方差函数var，标准差函数std，斜坡函数 skewness,肖度函数kurtosis正态分布概率函数normcdf。 四、实验结果与分析 例题3.7 x0=2; % 当 t=0 时x0=2 r=1; k=0.2; t=li nspace(tspa n( 1),tspa n( 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=fin d(t5); jj=jj(1);%% 找到t=5附近的时刻对应的位置 figure, plot(t,x1,-,LineWidth,2); % 画岀 r=1 时的曲线 text(t(jj),x1(jj)+0.2,r=1,k=0.2); % 文字说明图形表达式 1 内容 hold on, plot(t,x2,-,L in eWidth,2); % 画岀r=0.5 时的曲线； text(t(jj),x2(jj)+0.2,r=0.5,k=0.2); % 文字说明图形表达式 2内容 plot(t,x3,-,L in eWidth,2); % 画岀r=0.4 时的曲线 text(t(jj),x3(jj)+0.2,r=0.4,k=0.2); % 文字说明图形表达式 3内容 plot(t,x4,-,L in eWidth,2); % 画岀r=0.01 时的曲线 text(t(jj),x4(jj)+0.2,r=0.01,k=0.2); % 文字说明图形表达式 4内容 axis（［0,10,0,6］） % 取横坐标范围0到10,纵坐标范围0到6 title（ 不同参数条件下的logistic 模型的解）;% 添加标题 xlabel（ 时间t） % 给横坐标添加标题 ylabel（ 细菌数目）;% 给纵坐标添加标题 hold off; %%%e nd exam37.m 不同参数条件下的logistic模型的解 r=1,k=0.2 r=0.5,k=0.2 目数菌细 目数菌细 r=0.4,k=0.2 r=0.01,k=0.2 8 9 10 习题3.1 1.画出医学上常用的正态分布曲线，并确定其极值点 f(x)(x u)22 2 f(x) (x u)2 2 2 ,其中 u,