数学建模培训之四拟合与插值专题ppt培训课件

数学建模培训之四拟合与插值专题ppt培训课件

ID:33491431

大小:1.43 MB

页数:65页

时间:2018-05-23

数学建模培训之四拟合与插值专题ppt培训课件_第1页
数学建模培训之四拟合与插值专题ppt培训课件_第2页
数学建模培训之四拟合与插值专题ppt培训课件_第3页
数学建模培训之四拟合与插值专题ppt培训课件_第4页
数学建模培训之四拟合与插值专题ppt培训课件_第5页
资源描述:

《数学建模培训之四拟合与插值专题ppt培训课件》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、2012数学建模集训班----拟合与插值专题精彩源于坚持,搏过才知其美9/20/2021数学建模在大量的应用领域中,人们经常面临用一个解析函数描述数据(通常是测量值)的任务。对这个问题有两种方法。一种是插值法,数据假定是正确的,要求以某种方法描述数据点之间所发生的情况。另一种方法是曲线拟合或回归。人们设法找出某条光滑曲线,它最佳地拟合数据,但不必要经过任何数据点。本专题的主要目的是:了解插值和拟合的基本内容;掌握用Matlab求解插值与拟合问题的基本命令。9/20/2021数学建模内容提纲1.拟合问题引例及基本理论2.Matlab求解拟合问题3.

2、应用实例4.插值问题引例及基本理论5.Maltab求解插值问题6.应用实例9/20/2021数学建模拟合问题9/20/2021数学建模拟合问题引例1温度t(0C)20.532.751.073.095.7电阻R()7658268739421032已知热敏电阻数据:求600C时的电阻R。设R=at+ba,b为待定系数9/20/2021数学建模拟合问题引例2t(h)0.250.511.523468c(g/ml)19.2118.1515.3614.1012.899.327.455.243.01已知一室模型快速静脉注射下的血药浓度数据(t=0注射300

3、mg)求血药浓度随时间的变化规律c(t).作半对数坐标系(semilogy)下的图形MATLAB(aa1)9/20/2021数学建模曲线拟合问题的提法已知一组(二维)数据,即平面上n个点(xi,yi)i=1,…n,寻求一个函数(曲线)y=f(x),使f(x)在某种准则下与所有数据点最为接近,即曲线拟合得最好。+++++++++xyy=f(x)(xi,yi)ii为点(xi,yi)与曲线y=f(x)的距离9/20/2021数学建模线性最小二乘拟合f(x)=a1r1(x)+…+amrm(x)中函数{r1(x),…rm(x)}的选取1.通过机理分析建

4、立数学模型来确定f(x);++++++++++++++++++++++++++++++f=a1+a2xf=a1+a2x+a3x2f=a1+a2x+a3x2f=a1+a2/xf=aebxf=ae-bx2.将数据(xi,yi)i=1,…n作图,通过直观判断确定f(x):9/20/2021数学建模曲线拟合问题最常用的解法——线性最小二乘法的基本思路第一步:先选定一组函数r1(x),r2(x),…rm(x),m

5、最小二乘准则):使n个点(xi,yi)与曲线y=f(x)的距离i的平方和最小。记问题归结为,求a1,a2,…am使J(a1,a2,…am)最小。9/20/2021数学建模线性最小二乘法的求解:预备知识超定方程组:方程个数大于未知量个数的方程组即Ra=y其中超定方程一般是不存在解的矛盾方程组。如果有向量a使得达到最小,则称a为上述超定方程的最小二乘解。9/20/2021数学建模线性最小二乘法的求解定理:当RTR可逆时,超定方程组(3)存在最小二乘解, 且即为方程组RTRa=RTy------正则(正规)方程组的解:a=(RTR)-1RTy所以,曲

6、线拟合的最小二乘法要解决的问题,实际上就是求以下超定方程组的最小二乘解的问题。其中Ra=y(3)9/20/2021数学建模用MATLAB解拟合问题1、线性最小二乘拟合2、非线性最小二乘拟合9/20/2021数学建模用MATLAB作线性最小二乘拟合1.作多项式f(x)=a1xm+…+amx+am+1拟合,可利用已有命令:a=polyfit(x,y,m)3.对超定方程组可得最小二乘意义下的解。,用2.多项式在x处的值y的计算命令:y=polyval(a,x)输出拟合多项式系数a=[a1,…,am,am+1]’(数组)输入同长度数组X,Y拟合多项式 次

7、数9/20/2021数学建模即要求出二次多项式:中的使得:例对下面一组数据作二次多项式拟合9/20/2021数学建模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=Ry'MATLAB(zxec1)解法1.解超定方程的方法2)计算结果:A=[-9.8108,20.1293,-0.0317]9/20/2021数学建模2)计算结果:A=[-9.8108,20.1293,-0.0317]解

8、法2.用多项式拟合的命令MATLAB(zxec2)1)输入命令:x=0:0.1:1;y=[-0.447,1.978,3.28,6.16,

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。