数学实验第七次教程.ppt

数学实验之七 数据拟合;实验目的 [1] 了解最小二乘拟合的原理，掌握用 MATLAB作线性最小二乘拟合的方法。 [2] 通过实例学习如何用拟合方法解决 实际问题；通过实例理解参数辨识问题的几种方法。 ;主要内容;求电阻R随温度t的变化规律。; 设 R=at+b a,b为待定系数;引例2：血药浓度的变化规律;作图观察;主要内容;曲线拟合的基本原理——问题的提法;问题的数学模型 ;记; 特别，若选定一组函数 r1(x), r2(x), …rm(x), mn, 令 f(x,a1,a2, …am )=a1r1(x)+a2r2(x)+ …+amrm(x) 其中 a1,a2, …am 为待定系数。; 1) 除了最小二乘准则（即各点误差的平方和最小），你认为还可以用怎样的拟合准则？ 2) 比较起来，最小二乘准则有什么优点？;主要内容;最小二乘拟合函数的选取 ;最小二乘拟合函数的选取 ;主要内容;多项式拟合： 作多项式f(x)=a1xm+ …+amx+am+1函数拟合,可利用已有程序polyfit,其调用格式为:;例1. 由数据;2.用命令polyfit(x,y,m)作最小二 乘拟合 3.编写MATLAB程序dianzu1.m，并运行得到：;t=[20.5 32.5 51 73 95.7]; r=[765 826 873 942 1032]; aa=polyfit(t,r,1); a=aa(1) b=aa(2) y=polyval(aa,t); plot(t,r,k+,t,y,r);2. 曲线拟合： 作一般的最小二乘曲线拟合，可利 用已有程序lsqcurvefit,其调用格式为： [a,resnorm,residual]=lsqcurvefit(‘f’, a0, x, y) ;用MATLAB作最小二乘曲线拟合;a0=[1,1,1,0]; xdata=[0:.1:2]; %ydata省略; [a,resnorm,residual,flag,output]=lsqcurvefit(fitf,a0,xdata,ydata) xi=linspace(0,2,200); yi=fitf(a,xi); plot(xdata,ydata,ro,xi,yi) xlabel(x),ylabel(y=f(x)), title(nonlinear curve fitting);主要内容;范例1：静脉注射的给药方案 ;范例1：静脉注射的给药方案 ;若c1=10, c2=25(?g/ml), 给药方案设计归结为根据数据(ti,ci) i=1,…n (d 给定）拟合曲线c(t), 以确定系数k, v.;给药方案 ;主要内容;; 某种医用薄膜在试制时需测定其被物质分子穿透的能力。 测定方法：用面积为S的薄膜将容器分成两部份，在两部分中分别注满该物质的两种不同浓度的溶液。此时该物质分子就会从一侧向另一侧扩散。平均每单位时间通过单位面积薄膜的物质分子量与膜两侧溶液的浓度差成正比，比例系数K称为渗透率。定时测量容器中薄膜某一侧的溶液浓度，以此确定K。;3）薄膜是双向同性的即物质从膜的任何一侧向另一侧渗透的性能是相同的。;VA; 考察时段[t，t+Δt]薄膜的一侧容器中该物质质量的变化。 ;于是有：;分别表示在初始时刻两侧溶;从而：;基于假设和前面的分析，B侧的浓度CB(t)应满足如下微分方程和初始条件：;模型求解方法;;用MATLAB软件进行计算. 1）编写函数M-文件 nongdu.m function f = nongdu(x,tdata) f = x(1)+x(2)*exp(-0.02*x(3)*tdata); 其中 x(1) = a;x(2) = b;x(3) = k; 2) 编写M文件 (baomo.m) tdata = linspace(100,1000,10); cdata =[4.54 4.99 5.35 5.65 5.90 6.10 ... 6.26 6.39 6.50 6.59]; x0 = [0.2,0.05,0.05]; x = lsqcurvefit(‘nongdu’,x0,tdata,cdata);3) 输出结果: x = 0.007 -0.003 0.1012 即 k = 0.1012, a = 0.007, b = -0.003, ;2. 非线性规划法;3. 导函数拟合法;即为求参数K, a使下列误差函数达到最小：;4. 线性化迭代法;4. 线性化迭代法;1）当K的初值取为k=0.3时，出现