数据插值与拟合讲课文档.pptVIP

数据插值与拟合 第一页，共二十五页。 （优选）第九讲数据插值与拟合 第二页，共二十五页。 插值函数一般是已知函数的线性组合或者称为加权平均．插值在工程实践和科学实验中有着非常广泛而又十分重要的应用，例如，信息技术中的图像重建、图像放大中为避免图像的扭曲失真的插值补点、建筑工程的外观设计。化学工程实验数据与模型的分析、天文观测数据、地理信息数据的处理如(天气预报）以及社会经济现象的统计分析等等． 插值则要求函数在每个观测点处一定要满足 第三页，共二十五页。 1、船在该海域会搁浅吗 在某海域测得一些点(x，y)处的水深z（单位：英尺）由下表给出，水深数据是在低潮时测得的．船的吃水深度为5英尺，问在矩形区域(75,200)*(-50，150）里的哪些地方船要避免进人． 一、实例及其模型 第四页，共二十五页。 分析 由于测量点是散乱分布的，先在平面上作出测量点的分布图，再利用二维插值方法补充一些点的水深，然后作出海底曲面图和等高线图，并求出水深小于5的海域范围． 第五页，共二十五页。 在化学反应中，为研究某化合物的浓度随时间的变化规律，测得一组数据如表 2、浓度的变化规律 表中的数据反映了浓度随时间变化的函数关系，它是一种离散关系若需要推断20，40分钟时的浓度值，能否用一个显函数y=f(t)来拟合表中的离散数据，然后再计算浓度值f(20), f（40）? 第六页，共二十五页。 问题分析 （1）首先将这些离散数据分布在直角坐标系下，由此可发现浓度与时间之间呈现什么规律．这种数据分布在直角坐标系下的图形被称为散点图； （2）根据散点图，判段它接近于哪类函数曲线， 即确定函数形式 （3）函数形式确定以后，关键是要确定函数中含有的 待定参数。 最常用的确定待定系数的方法是，曲线拟合的最小二乘法 第七页，共二十五页。 二、 插值与拟合 1、插值方法 （1）分段线性插值 分段线性插值的提法如下： 第八页，共二十五页。 （2）分段三次埃尔米特插值 在插值问题中，如果除了插值节点的函数值给定外，还要求在节点的导数值为给定值，即插值问题变为 相当于在每一小段上应满足四个条件（方程），可以确定四个待定参数．三次多项式正好有四个系数，所以可以考虑用三次多项式函数作为插值函数，这就是分段三次埃尔米特插值，它与分段线性插值一起都称为分段多项式插值 第九页，共二十五页。 （3）三次样条插值 第十页，共二十五页。 2、曲线拟合的最小二乘法 给定平面上的点 进行曲线拟合有多种方法，最小二乘法是解决曲线拟合最常用的一种方法 最小二乘法的原理是求f(x),使 达到最小 简单地说，最小二乘法准则就是使所有散点到曲线的距离平方和最小 第十一页，共二十五页。 线性最小二乘法 拟合函数可由一些简单的“基函数”（例如幂函数，三角 函数等等） 来线性表示 现在要确定系数 使 达到极小为此 第十二页，共二十五页。 三、插值的matlab实现 1、一维插值 MATLAB中的插值函数为interp1，其调用格式为 其中x，y为插值点，yi为在被插值点xi处的插值结果，x, y为向量。 注意：所有的插值方法都要求x是单调的，并且xi不能够超过x的范围。 第十三页，共二十五页。 MATLAB提供的插值方法有几种 表示采用的插值方法 ：分段线性插值 ：三次Hermite插值（立方插值） ：三次分段样条插值 ：最近点等值方式 缺省时表示线性插值 第十四页，共二十五页。 例1 在一 天24小时内，从零点开始每间隔2小时测得的环境温度数据分别为 12，9，9，1，0，18 ，24，28，27，25，20，18，15，13， 推测中午（即13点）时的温度． x=0：2：24； y=[12 9 9 10 18 24 28 27 25 20 18 15 13]； x1＝13 ； y1＝interp1（x，y，x1，‘spline’） 若要得到一天24小时的温度曲线 x=0：2：24； y=[12 9 9 10 18 24 28 27 25 20 18 15 13] xi＝0：1／3600：24； yi=interp1(x,y,xi,’spline’ ); plot(x, y, ’o’, xi, yi) 第十五页，共二十五页。 2、高维插值 N维插值函数interpN（） 其中N可以为2，3，……，如N＝2为二维插值，调用格式为 其中 x，y，z为插值节点，zi为被插值点(xi，yi)处的插值结果 且， xi， yi为被插值节点构成的新的网格数据 ‘me