一讲插值与拟合2课件

1、2013数学建模培训10/10/202212013数学建模培训10/3/20221培训内容与方式10/10/20222培训内容与方式10/3/20222培训内容1. 插值与拟合；2. 灰色系统；3. 层次分析法；(重要)4. 模糊数学方法；(重要)5. 多元统计分析与Spss；(重要)6. 数学规划与Lingo；(重要)7. 图论模型及其Matlab程序；10/10/20223培训内容10/3/202238. 建模论文写作；9. 建模案例讲评。 培训内容覆盖了数学建模中大部分常用的数学方法。通过刻苦钻研，努力学习，掌握了上述方法后，必将大大提高建立数学模型和运用计算机解决实际问题的综合能力，也

2、一定会在全国建模大赛中取得较好成绩。10/10/202248. 建模论文写作；10/3/20224培训方式 由于培训内容众多且有相当难度,而培训时间又较短，所以如果学生不事先认真预习相关内容，那么在课堂上不可能完全听懂老师所讲内容，大部分学生会云山雾罩，一头雾水，从而使得培训效果大打折扣。 综上，强烈建议学生课前认真、10/10/20225培训方式10/3/20225反复研读授课PPT和其它相关资料，然后带着疑问和兴趣再听老师讲解，这样才能保证培训效果。 培训PPT及软件、程序下载邮箱: MM: matlabmaple 培训以讲解为主，学生如有疑问, 可在下午上机辅导时向老师咨询。10/10/

3、20226反复研读授课PPT和其它相关资料，然后带着疑问和兴趣再听老师上机操作与练习10/10/20227上机操作与练习10/3/20227 众所周知，在建模竞赛中，能否熟练使用相关数学软件是能否取得好成绩的关键之一。因此，数学软件的培训应该是建模培训的重要内容。 建模中常用的数学软件有Matlab, Spss, Lingo, Maple等。 由于培训时间有限，在课堂上只能简单介绍Spss和Lingo，而Matlab和10/10/20228 众所周知，在建模竞赛中，能否熟练使用相关数学Maple只能在上机过程中穿插介绍。 几乎所有的数学方法最终都要用数学软件和程序实现，所以上机练习是与课堂讲授

4、同等重要的培训内容。 上机训练的主要目的和内容是通过练习，掌握实现各类数学方法的数学软件和程序。 上机所用数学软件和程序主要由10/10/20229Maple只能在上机过程中穿插介绍。10/3/20229教师提供，学生只要会用即可。 但需要提醒同学们注意的是，现成的软件和程序不可能解决建模中的所有问题。通过上机训练，掌握一些基本编程和计算技能(如用Matlab做数据处理、画图, 用Maple做简单的解析计算)，对于参加建模竞赛是绝对必要的。10/10/202210教师提供，学生只要会用即可。10/3/202210 上机前，培训老师会布置上机练习，提供相关软件或程序，讲解关键步骤和程序语句。上机

5、过程中，辅导老师负责解答学生的疑难问题。题外话 特别要提醒大家的是，国赛完全不同于网络挑战赛，题目的难度和开放性适中，有答题要点、参考方法，10/10/202211 上机前，培训老师会布置上机练习，提供相关软件甚至有参考答案；初评评委大部分由指导教师担任，评阅结果通常比较合理、靠谱。那种连基本建模方法都不懂竟然也能获大奖的现象在国赛中几乎是不可能发生的，全国一等奖的论文必须模型、方法合理，结果基本可信、正确，写作清晰、规范。 一句话，建模不能全靠忽悠。10/10/202212甚至有参考答案；初评评委大部分由指导教师担任，评阅结果通常比第一讲 插值与拟合10/10/202213第一讲 插值与拟合

6、10/3/202213一、引 言10/10/202214一、引 言10/3/202214 插值与拟合属数值分析中函数逼近内容。在数学建模竞赛中，插值与拟合是一种基本的数据分析手段，被公认为建模中的常用算法之一。 本讲第一部分首先介绍了插值问题、插值原理、高次插值的Runge现象，然后讲解了Matlab中的一维和二维插值命令，最后简要介绍了具有10/10/202215 插值与拟合属数值分析中函数逼近内容。在数学建强大插值功能的三维图形(曲面图，等高线等)绘制软件Surfer。 第二部分首先介绍了拟合问题、拟合原理与步骤，然后介绍了Origin中的拟合功能及专业拟合软件TableCurve 2D,

7、 3D。 第三部分给出了一个应用插值和拟合的建模实例。10/10/202216强大插值功能的三维图形(曲面图，等高线等)绘制软件Surfe 本讲首先要理解插值问题和拟合问题的特点，其次要了解插值和拟合的原理与方法，还要熟悉Matlab插值与绘图的相关命令，掌握Matlab编程的基本知识与技能，特别要熟练掌握利用Surfer, Origin, TableCurve等软件进行插值和拟合的技能。 下面给出本讲学习大纲，以方便10/10/202217 本讲首先要理解插值问题和拟合问题的特点，其次大家学习。 1. 插值问题和拟合问题的特点及区别；在实际中如何正确地判断、选择插值或拟合方法？ 2. 高次插

8、值的Runge现象及避免方法； 3. Matlab一维插值(interp1), 二维插值(interp2), 散乱点插值(griddat10/10/202218大家学习。10/3/202218a)及相关命令(surf, mesh, meshgrid, contour)；4. Surfer绘图的步骤；5. 了解Surfer中几种插值的意义；6. Origin和TableCurve拟合的步骤；7. Matlab基础知识, 如数组及运算、调用，循环与控制语句，绘图相关命令，函数(m文件)的定义和调用等。10/10/202219a)及相关命令(surf, mesh, meshgrid, 二、插 值10

9、/10/202220二、插 值10/3/2022201. 插值问题 例1 在一天24小时内，从零点开始每间隔2小时测得的环境温度数据分别为12，9，9，10，18 ，24，28， 27，25，20，18，15，13，推测中午1点温度，并做出24小时温度变化曲线图。10/10/2022211. 插值问题10/3/202221 例2 已知飞机下轮廓线上数据如下，画出飞机下轮廓线。 机翼下轮廓线xy10/10/202222 例2 已知飞机下轮廓线上数据如下，画出飞机下 例3 测得平板表面3*5网格点处的温度分别为： 82 81 80 82 84 79 63 61 65 81 84 84 82 85

10、86 做出平板表面的温度分布曲面z=f(x,y)的图形及等温线，并求出温度最高和最低点。 10/10/202223 例3 测得平板表面3*5网格点处的温度分别为 上述问题可归结为“已知函数在某区间(域)内若干点处的值, 求函数在该区间(域)内其它点处的值”，这种问题适宜用插值方法解决。 一维插值问题可描述为：已知函数在x0, x1, , xn处的值y0, y1, yn，求简单函数 p(x)，使 p(xi) = yi。 通常取 p(x)为多项式。10/10/202224 上述问题可归结为“已知函数在某区间(域)内若 可以用范德蒙行列式和克莱姆法则证明(习题集第一章最后一题)： 在x0, x1,

11、, xn处取值y0, y1, , yn的多项式存在且唯一，即插值问题的解唯一存在。 常用的插值方法有Lagrange插值法和Newton插值法。10/10/202225 可以用范德蒙行列式和克莱姆法则证明(习题集第2. 高次插值的Runge现象 在研究插值问题的初期，所有人都想当然地认为插值多项式的次数越高，插值精度越高。 Runge 通过对一个例子的研究发现，上述结论仅仅在插值多项式的次数不超过七时成立；插值多项式的次数超过七时，插值多项式会出现严重10/10/2022262. 高次插值的Runge现象10/3/202226的振荡现象，称之为Runge现象。10/10/202227的振荡现象

12、，称之为Runge现象。10/3/202227 因此，在实际中不应使用七次以上的插值。 避免Runge现象的常用方法是：将插值区间分成若干小区间，在小区间内用低次(二次，三次)插值，即分段低次插值，如样条函数插值。10/10/202228 因此，在实际中不应使用七次以上的插值。10/样条插值结果10/10/202229样条插值结果10/3/202229三、Matlab插值10/10/202230三、Matlab插值10/3/202230 Maple和Matlab都可以进行插值计算，Maple的一维插值计算较为便捷，而Matlab的二维插值功能较强,还能进行散乱点插值。 本节主要介绍Matlab

13、的一维和二维插值命令，大家务必要通过上机操作熟悉这些命令，同时还要初步掌握Matlab的基础知识与技能。10/10/202231 Maple和Matlab都可以进行插值计算，1. 一维插值 一维插值命令是interp1, 其基本格式为yi=interp1(x,y,xi, method)。 x,y为插值点，xi,yi为被插值点和插值结果，x,y和xi,yi通常为向量；method表示插值方法：nearest最邻近插值，linear线性插值， spline三次样条插值，cubic立10/10/2022321. 一维插值10/3/202232方插值，缺省为线性插值。例1Matlab程序 x=0:2:

14、24; y=12 991018 2428272520 18 15 13; x1=13; y1=interp1(x,y,x1,spline) xi=0:1/3600:24; yi=interp1(x,y,xi, spline); plot(x,y, *,xi,yi)10/10/202233方插值，缺省为线性插值。10/3/202233 请理解掌握程序中的每个语句，并改变插值方法，观察图形变化。10/10/20223410/3/202234例2Matlab程序function planex0=0 3 5 7 9 11 12 13 14 15 ;y0=0 1.2 1.7 2.0 2.1 2.0 1.

15、8 1.2 1.0 1.6 ;x=0:0.1:15;y1=lagrange(x0,y0,x);y2=interp1(x0,y0,x);y3=interp1(x0,y0,x,spline);subplot(3,1,1);10/10/202235例2Matlab程序10/3/202235plot(x0,y0,k+,x,y1,r);grid;title(lagrange);subplot(3,1,2);plot(x0,y0,k+,x,y2,r);grid;title(piecewise linear);subplot(3,1,3);plot(x0,y0,k+,x,y3,r);10/10/202236

16、plot(x0,y0,k+,x,y1,r);10/3grid;title(spline);function y=lagrange(x0,y0,x)n=length(x0);m=length(x);for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; 10/10/202237grid;10/3/202237 for j=1:n if j=k p=p*(z-x0(j)/(x0(k)-x0(j); end end s=p*y0(k)+s; end y(i)=s;end10/10/202238 for j=1:n10/3/20223810/10/20223910/3/20

17、2239 例2的程序较复杂，说明如下: (1) 程序依次用Lagrange、分段线性和三次样条三种插值方法进行了计算，其中Lagrange高次插值明显出现了Runge现象； (2) 因为Matlab没有Lagrange高次插值功能，所以程序中单独编写了高次插值函数lagrange，然后调用；10/10/202240 例2的程序较复杂，说明如下:10/3/202 (3) 程序涉及到了数组、循环和条件语句、子函数的定义与调用以及一些绘图命令(subplot,grid,title)等。 请大家通过上机理解、掌握上述命令，特别是函数的定义及调用。10/10/202241 (3) 程序涉及到了数组、循环

18、和条件语句、子2. 二维插值 二维插值命令是interp2, 基本格式为zi=interp2(x,y,z,xi,yi,method)。 二维插值命令的使用较复杂。 x,y,z为插值点，z可以理解为被插值函数在(x,y)处的值；xi,yi为被插值点, zi为输出的插值结果，可理解为插值函数在(xi,yi)处的值；x,y为向10/10/2022422. 二维插值10/3/202242量，xi,yi为向量或矩阵，而z和zi则为矩阵。 method表示插值方法：neares t最邻近插值, linear双线性插值, spline双三次样条插值，cubi c双立方插值，黙认双线性插值。10/10/202

19、243量，xi,yi为向量或矩阵，而z和zi则为矩阵。10/3/2例3Matlab程序x=1:5;y=1:3;temps=82 81 80 82 84;79 63 61 65 81;84 84 82 85 86;figure(1);mesh(x,y,temps);xi=1:0.2:5;yi=1:0.2:3;10/10/202244例3Matlab程序10/3/202244zi=interp2(x,y,temps,xi,yi,cubic);figure(2);mesh(xi,yi,zi);figure(3);contour(xi,yi,zi,20,r);i,j=find(zi=min(min(z

20、i);x=xi(j),y=yi(i),zmin=zi(i,j)i,j=find(zi=max(max(zi);x=xi(j),y=yi(i),zmax=zi(i,j)10/10/202245zi=interp2(x,y,temps,xi,yi,c 上述程序较复杂，说明如下： (1) interp2中的xi为行向量, 而yi为列向量, 其实xi和yi行列不同即可。 (2) plot3(空间曲线), mesh(空间曲面), surf (空间曲面), contour(等高线)是三维作图中的常用命令。mesh和surf的区别是：mesh画的是曲面网格图，而surf画的是曲面表面图。10/10/2022

21、46 上述程序较复杂，说明如下：10/3/2022 contour(x,y,z,n)的功能是作出由点(x,y,z)插值而成曲面的n条等高线。 用meshc和surfc可在曲面下方画等高线。meshz和surfz是画垂帘图。 (3) 程序的最后部分为求最高(低)点，请各位通过上机揣摩min, max特别是find的功能。 (4) 将程序中 figure 语句去除，通10/10/202247 contour(x,y,z,n)的功能是作出过观察结果，体味figure的作用。 例4 在某山区测得一些地点的高程如下表。平面区域为 0 x5600, 0y12时的数据出现了巨大误差。 样条插值的第一个图表明

22、样条插值可以很好地解决第一个问题，但第二个图显示，用样条插值外推t12时的数据也会产生较大误差。 综上，第二问不宜用插值方法。10/10/202295 高次插值的第一个图中出现的明显的Runge现例6多项式(9次)拟合图形R=0.9442810/10/202296例6多项式(9次)拟合图形10/3/202296例7Origin拟合结果Equation y = A1*exp(-x/t1) + A2*exp(-x/t2) + y0Adj. R-Square0.99864ValueStandard ErrorBy011.324160.71521BA1-2.766120.23031Bt111.2397

23、18.0397BA2-9.351890.51014Bt21.511270.1676510/10/202297例7Origin拟合结果Equation y = A1*e例7Origin拟合图形10/10/202298例7Origin拟合图形10/3/2022983. TableCurve简介 TableCurve2D,3D是美国SPSS公司开发的曲线和曲线拟合软件。 TableCurve最独特的功能是能自动找出与数据最匹配的曲线(面)。 TableCurve2D可以从3600多个方程中找出与数据最匹配的曲线。 TableCurve3D可以从4亿多个方10/10/2022993. TableCur

24、ve简介 10/3/20程中找出与数据最匹配的曲面。 与Origin一样，TableCurve也可以自创曲线和曲面方程。 TableCurve使用简单, 只需经过读入数据(通常为Excel文件)，选取数据列，即可拟合出所有曲线或曲面。 这里要提醒各位的是，有些拟合方程虽然效果很好，但方程过于复杂,10/10/2022100程中找出与数据最匹配的曲面。10/3/2022100不便于进一步分析与处理，因此需要综合考虑，选取合适的拟合方程。 下面给出例6的TableCurve2D拟合结果以及2013网络挑战赛A题数据的TableCurve3D拟合结果。10/10/2022101不便于进一步分析与处理

25、，因此需要综合考虑，选取合适的拟合方程例6 TableCurve2D拟合结果10/10/2022102例6 TableCurve2D拟合结果10/3/2022102013A TableCurve3D拟合结果10/10/20221032013A TableCurve3D拟合结果10/3/202七、MCM/ICM1991A题水塔流量的估计10/10/2022104七、MCM/ICM1991A题水塔流量的估计10/3/2 某社区自来水由一个高12.2米, 直径17.4米圆柱形水塔提供。当水位降至8.2米时, 水泵自动启动加水；水位升高到10.8 米时，水泵停止工作；一般水泵每天工作两次。下表给出了某

26、天不同时间水位数据，其中有三次观察时水泵正在供水，无水位记录。 试建立适当的数学模型，计算任10/10/2022105 某社区自来水由一个高12.2米, 直径17.意时刻的水流速度，估计一天的用水量和水泵的工作功率。时刻(h)水位(cm) 0 0.92 1.84 2.95 3.87 4.98 5.90 7.01 7.93 8.97968 948 931 913 898 881 869 852 839 822时刻(h)水位(cm)9.98 10.92 10.95 12.03 12.95 13.88 14.98 15.90 16.83 17.93 / / 1082 1050 1021 994 965 941 918 892时刻(h)水位(cm)19.04 19.96 20.84 22.01 22.96 23.88 24.99 25.91 866 843 822 / 10.82 1059 1035 101810/10/2022106意时刻的水流速度，估计一天的用水量和水泵的工作功率。时刻(h1. 问题分析 题目中