插值与拟合-MATLAB学习之四剖析.pptVIP

流量曲线见图 n=(3,4) n=(5,6) 用给定的多项式，如y=x3-6x2+5x-3，产生一组数据(xi,yi，i=1,2,…,n),再在yi上添加随机干扰(可用rand产生(0,1)均匀分布随机数,或用rands产生N(0,1)分布随机数)，然后用xi和添加了随机干扰的yi作的3次多项式拟合，与原系数比较。 如果作2或4次多项式拟合，结果如何？ 练习1 练习2 山区地貌：在某山区测得一些地点的高程如下表：(平面区域1200=x=4000,1200=y=3600)，试作出该山区的地貌图和等高线图，并对几种插值方法进行比较。 　将估计水塔的流量问题整理成一个完善的建模论文（word 文档） 练习3 谢谢大家 其它问题(自学) 电阻问题 温度t(0C) 20.5 32.7 51.0 73.0 95.7 电阻R(?) 765 826 873 942 1032 例. 由数据 拟合R=a1t+a2 方法1.用命令 polyfit(x,y,m) 得到 a1=3.3940, a2=702.4918 方法2.直接用 结果相同。 一室模型：将整个机体看作一个房室，称中心室，室内血药浓度是均匀的。快速静脉注射后，浓度立即上升；然后迅速下降。当浓度太低时，达不到预期的治疗效果；当浓度太高，又可能导致药物中毒或副作用太强。临床上，每种药物有一个最小有效浓度c1和一个最大有效浓度c2。设计给药方案时，要使血药浓度 保持在c1~c2之间。本题设c1=10，c2=25(ug/ml). 拟 合 问 题 实 例 2 给药方案 —— 一种新药用于临床之前，必须设计给药方案. 药物进入机体后血液输送到全身，在这个过程中不断地被吸收、分布、代谢，最终排出体外，药物在血液中的浓度，即单位体积血液中的药物含量，称为血药浓度。 在实验方面,对某人用快速静脉注射方式一次注入该药物300mg后,在一定时刻t(小时)采集血药,测得血药浓度c(ug/ml)如下表: t (h) 0.25 0.5 1 1.5 2 3 4 6 8 c (?g/ml) 19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01 要设计给药方案,必须知道给药后血药浓度随时间变化的规律。从实验和理论两方面着手： 给药方案 1. 在快速静脉注射的给药方式下，研究血药浓度（单位体积血液中的药物含量）的变化规律。 t c2 c c1 0 ? 问题 2. 给定药物的最小有效浓度和最大治疗浓度，设计给药方案：每次注射剂量多大；间隔时间多长。 分析 理论：用一室模型研究血药浓度变化规律 实验：对血药浓度数据作拟合，符合负指数变化规律 3.血液容积v, t=0注射剂量d, 血药浓度立即为d/v. 2.药物排除速率与血药浓度成正比，比例系数 k(0) 模型假设 1. 机体看作一个房室，室内血药浓度均匀——一室模型 模型建立 在此，d=300mg，t及c（t）在某些点处的值见前表，需经拟合求出参数k、v 用线性最小二乘拟合c(t) 计算结果： d=300; t=[0.25 0.5 1 1.5 2 3 4 6 8]; c=[19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01]; y=log(c); a=polyfit(t,y,1) k=-a(1) v=d/exp(a(2)) 程序： 用非线性最小二乘拟合c(t)-用lsqcurvefit 2、主程序lihe2.m如下 clear tdata=[0.25 0.5 1 1.5 2 3 4 6 8]; cdata=[19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01]; x0=[10,0.5]; x=lsqcurvefit(curvefun3,x0,tdata,cdata); f=curvefun3(x,tdata) x 1、用M-文件curvefun3.m定义函数 function f=curvefun3(x,tdata) d=300 f=(x(1)\d)*exp(-x(2)*tdata) % x(1)=v; x(2)=k 给药方案 设计 c c2 c1 0 ? t 设每次注射剂量D, 间隔时间? 血药浓度c(t) 应c1? c(t) ? c2 初次剂量D0 应加大 给药方案记为： 2、 1、 计算结果： 给药方案： c