插值和拟合(讲义)培训教程文件.pptVIP

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)) 程序： 通过此例对最近邻点插值、双线性插值方法和双三次插值方法的插值效果进行比较． figure(1); meshz(x,y,z) xlabel(X) ylabel(Y) zlabel(Z) xi=0:50:5600; yi=0:50:4800; figure(2) z1i=interp2(x,y,z,xi,yi,nearest); surfc(xi,yi,z1i) xlabel(X),ylabel(Y),zlabel(Z) figure(3) z2i=interp2(x,y,z,xi,yi); surfc(xi,yi,z2i) xlabel(X),ylabel(Y),zlabel(Z) figure(4) z3i=interp2(x,y,z,xi,yi,cubic); surfc(xi,yi,z3i) xlabel(X),ylabel(Y),zlabel(Z) figure(5) subplot(1,3,1),contour(xi,yi,z1i,10,r); subplot(1,3,2),contour(xi,yi,z2i,10,r); subplot(1,3,3),contour(xi,yi,z3i,10,r); 插值函数griddata格式为: cz =griddata（x，y，z，cx，cy，‘method’） 用MATLAB作散点数据的插值计算 要求cx取行向量，cy取为列向量． 被插值点 插值方法 插值节点 被插值点的函数值 ‘nearest’最邻近插值 ‘linear’ 双线性插值 ‘cubic’ 双三次插值 v4- MATLAB提供的插值方法 缺省时, 双线性插值 例6. 在某海域测得一些点(x,y)处的水深z由下表给出，船的吃水深度为5英尺，在矩形区域（75，200）×（-50，150）里的哪些地方船要避免进入． x=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5]; y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5]; z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9]; cx=75:0.5:200; cy=-50:0.5:150; cz=griddata(x,y,z,cx,cy,cubic); 4.作出水深小于5的海域范围,即z=5的等高线. . . 3 作海底曲面图 meshz(cx,cy,cz),rotate3d xlabel(X),ylabel(Y),zlabel(Z) figure(2),contour(cx,cy,cz,[-5 -5]);grid hold on plot(x,y,+) xlabel(X),ylabel(Y) 用MATLAB解拟合问题 1.线性最小二乘拟合 2.非线性最小二乘拟合 用MATLAB作线性最小二乘拟合 1. 作多项式f(x)=a1xm+ …+amx+am+1拟合,可利用已有程序: a=polyfit(x,y,m) 2. 对超定方程组 可得最小二乘意义下的解． ，用 3.多项式在x处的值y可用以下命令计算：y=polyval（a，x） 输出拟合多项式系数 a=[a1, …,am , am+1] (数组)） 输入同长度 的数组x，y 拟合多项 式次数 即要求 出二次多项式: 中 的 使得: 例7. 对下面一组数据作二次多项式拟合 1）输入以下命令： x=0:0.1:1; y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; R=[(x.^2)‘ x’ ones(11,1)]; A=R\y 解法1．用解超定方程的方法 2）计算结果： Ａ = -9.8108 20.1293 -0.0317 1）输入以下命令： x=0:0.1:1; y=[-