实验五医用数学在MATLAB中的实现1.docVIP

实验名称:医用数学在MATLAB中的实现1 时间：2013-4-15 地点：三机房 教 师 评 语 成绩 教师签名 一、实验目的 Matlab提供了强大的数据分析功能，本次实验旨在熟悉和掌握应用Matlab实现函数分析、微分方程求解的功能。 二、实验原理 1.函数分析 （1）自变量为单个变量时，称为一元函数，当自变量为多个变量时，称为多元函数，表示为：y=f(x1,…,xn),x1 D1,…,xn Dn 在进行函数分析的时候，将每个函数编写成M函数文件即可。 （2）函数拟合 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,就会弹出一个曲线拟合工具界面界面窗口。 2.微分方程在生物医学中的应用 （1）对于简单的微分方程或微分方程组，可以用symbolic math toolbox 里德dsolve命令来获得解析解。 Dx(t)/dt=kx(t) 使用dsolve命令如下： f=dsolve(‘Dx=k*x’,’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为对应时间点的待求变量值。 3.常用统计量函数 算数均值函数mean，几何均值函数geomean,中位数函数median,极差函数range,四分位数间距的函数iqr，方差函数var，标准差函数std，斜坡函数skewness,峭度函数kurtosis,正态分布概率函数normcdf。 三、实验仪器和设备 1、计算机一台 以及MATLAB R2009a软件 四、实验内容和步骤 （习题第三章） 五．程序流程图及程序清单 fplot((1/sqrt(2*pi))*exp(-x.^2/2),[-6.5,6.5]) x=-6.5:0.01:6.5; y=(1/sqrt(2*pi))*exp(-x.^2/2); j=find(y==max(y)) x(j) y(j) j = 651 ans = 0 ans = 0.3989 x=randn(1000,1) xmean=mean(x) xvar=var(x) xmean=median(x) xmin=min(x) xmax=max(x) 运行结果： xmean = -0.0431 xvar = 0.8902 xmean = -0.0131 xmin = -2.6442 xmax = 2.7316 normcdf(16.1,19.5,2.3) 1-normcdf(22.9,19.5,2.3) normcdf(14.6,19.5)-normcdf(23.9,19.5,2.3) normcdf(16.1,19.5,2.3) 1-normcdf(22.9,19.5,2.3) normcdf(14.6,19.5)-normcdf(23.9,19.5,2.3) ns = 0.0697 ans = 0.0697 ans = -0.9721 六、实验总结 Matlab是一种使用方便，帮助全面的工具箱，它的编程语言简单易学，科学计算能力强大，数据可视化功能强大，使用m文件，可以让有些复杂的数据结构运算变得简单明了，因此学好matlab软件变得极其的重要，而它也得到了广泛的应用。 Mean： 算数均值函数， geomean: 几何均值函数 median: 中位数函数 range: 极差函数 iqr: 四分位数间距的函数 var: 方差函数 std: 标准差函数 skewness: 斜坡函数 kurtosi: 峭度函数 normcdf: