数模培训数据拟合方法.pptx

数据拟合与线性最小二乘法;教学内容; 设 R=at+b a,b为待定系数; 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;数据拟合问题的提法;数据拟合问题的求解思路;线性最小二乘法原理;线性最小二乘法原理;线性最小二乘法原理;线性最小二乘法原理;算例;算例;求解法方程组 ;算例;算例;算例;算例;算例;算例;算例;用Matlab进行数据拟合;例1 已知观测数据点如表所示;t=0:0.1:1.2 s=polyval(p3,t) s1=polyval(p6,t) hold on plot(t,s,r-,linewidth,2) plot(t,s,b--,linewidth,2) grid;例2 用切削机床进行金属品加工时, 为了适当地调整机床, 需要测定刀具的磨损速度. 在一定的时间测量刀具的厚度, 得数据如表所示:;解: 描出散点图, 在命令窗口输入:;解: 描出散点图, 在命令窗口输入:;例2 用切削机床进行金属品加工时, 为了适当地调整机床, 需要测定刀具的磨损速度. 在一定的时间测量刀具的厚度, 得数据如表所示:;例3 一个15.4cm×30.48cm的混凝土柱在加压实验中的应力-应变关系测试点的数据如表所示;已知应力-应变关系可以用一条指数曲线来描述, 即假设;在命令窗口输入:;已知应力-应变关系可以用一条指数曲线来描述, 即假设;在实际应用中常见的拟合曲线有:;2. 非线性曲线拟合: lsqcurvefit.;例4 已知观测数据点如表所示;例4 已知观测数据点如表所示;编写下面的程序调用拟合函数.;例4 已知观测数据点如表所示;f(x)= 3ex+ 4.03x2 + 0.94 x3.;作业：用函数f(x)=a1*exp(-a2*x)+a3*exp(-a4*x)拟合 下列数据点： xdata=[0:.1:2] ydata=[5.8955 3.5639 2.5173 1.9790 1.8990 1.3938 1.1359 1.0096 1.0343 0.8435 0.6856 0.6100 0.5392 0.3946 0.3903 0.5474 0.3459 0.1370 0.2211 0.1704 0.2636];作业：考察某种纤维的强度与其拉伸倍数的关系,下表是实际测定的24个纤维样品的强度与相应的拉伸倍数的记录:;;x1=[2:16]; y1=[6.42,8.2,9.58,9.5,9.7,10,9.93,9.99,10.49,10.59,10.6,10.8,10.6,10.9,10.76]; b01=[0.1435,0.084]; %初始参数值 fun1=inline(x./(b(1)+b(2)*x),b,x); % 定义函数 [b1,r1,j1]=nlinfit(x1,y1,fun1,b01); y=x1./(0.1152+0.0845*x1); %根据b1写出具体函数 plot(x1,y1,*,x1,y,-or);;;作业 给药方案;例5 给药方案;例5 给药方案;2. 给药方案设计;2. 给药方案设计;2. 由实验数据作曲线拟合以确定参数;问题化为由数据 ti , yi ( i=1,…,8 ) 拟合直线;拟合曲线为:;3. 结论;范例：薄膜渗透率的测定;二、问题分析 考察时段[t，t+Δt]薄膜两侧容器中该物质质量的变化。 ; 以容器A侧为例，在时段[t，t+Δt]物质质量的增量为：;于是有：;从而：;三、数学模型 假设： 1）薄膜两侧的溶液始终是均匀的； 2）平均每单位时间通过单位面积薄膜的物质分 子量与膜两侧溶液的浓度差成正比。 3）薄膜是双向同性的即物质从膜的任何一侧向 另一侧渗透的性能是相同的。 基于假设和前面的分析，B侧的浓度CB(t)应满足如下微分方程和初始条件：;四、求解方法：;;用MATLAB软件进行计算. 1）编写函数M-文件 nongdu.m function f = nongdu(x,tdata) f = x(1)+x(2)*exp(-0.02