《插值拟合课件.pptx》由会员分享，可在线阅读，更多相关《插值拟合课件.pptx（48页珍藏版）》请在三一办公上搜索。
1、例：在1-12的11小时内，每隔1小时测量一次温度，测得的温度依次为：5，8，9，15，25，29，31，30，22，25，27，24。试估计每隔1/10小时的温度值。,第1页/共48页,拉格朗日插值,分段线性插值,三次样条插值,一、插值的定义,二、插值的方法,三、用Matlab解插值问题,一维插值,第2页/共48页,一维插值的定义,第3页/共48页,第4页/共48页,第5页/共48页,称为拉格朗日插值基函数。,已知函数f(x)在n+1个点x0,x1,xn处的函数值为 y0,y1,yn 。求一n次多项式函数Pn(x)，使其满足： Pn(xi)=yi,i=0,1,n.,解决此问题的拉格朗日插值多
2、项式公式如下:,其中Li(x) 为n次多项式：,拉格朗日(Lagrange)插值,第6页/共48页,拉格朗日(Lagrange)插值,特别地:,两点一次(线性)插值多项式:,三点二次(抛物)插值多项式:,第7页/共48页,For example,取最接近x=5的点， x0=1， x1=4，x2=9为插值节点，运用插值公式L(5)=2.27.,第8页/共48页,For example,取最接近x=5的点， x0=1， x1=4，x2=9为插值节点，运用插值公式L(5)=2.27.,第9页/共48页,采用拉格朗日多项式插值：选取不同插值节点个数n+1，其中n为插值多项式的次数，当n分别取2,4,6
3、,8,10时，绘出插值结果图形。,例,第10页/共48页,第11页/共48页,分段线性插值,计算量与n无关;n越大，误差越小.,第12页/共48页,例 已知飞机下轮廓线上数据如下，求x每改变0.1时的y值。,第13页/共48页,比分段线性插值更光滑。,在数学上，光滑程度的定量描述是：函数(曲线)的k阶导数存在且连续，则称该曲线具有k阶光滑性。 光滑性的阶次越高，则越光滑。是否存在较低次的分段多项式达到较高阶光滑性的方法？三次样条插值就是一个很好的例子。,三次样条插值,第14页/共48页,三次样条插值,第15页/共48页,第16页/共48页,第17页/共48页,用MATLAB作插值计算,一维插值
4、函数：,yi=interp1(x，y，xi，method),nearest ：最邻近插值linear ： 线性插值；spline ： 三次样条插值；cubic ： 立方插值。缺省时： 分段线性插值。,注意：所有的插值方法都要求x是单调的，并且xi不能够超过x的范围。,第18页/共48页,2.拟合的基本原理,1. 拟合问题引例,拟 合,第19页/共48页,拟 合 问 题 引 例 1,求600C时的电阻R。,设 R=at+ba,b为待定系数,第20页/共48页,函数插值与曲线拟合都是要根据一组数据构造一个函数作为近似，由于近似的要求不同，二者的数学方法上是完全不同的。,实例：下面数据是某次实验所得
5、，希望得到X和 f之间的关系？,问题：给定一批数据点，需确定满足特定要求的曲线或曲面,解决方案：,若不要求曲线（面）通过所有数据点，而是要求它反映对象整体的变化趋势，这就是数据拟合，又称曲线拟合或曲面拟合。,若要求所求曲线（面）通过所给所有数据点，就是插值问题；,拟合与插值的关系,第21页/共48页,线性插值、样条插值与曲线拟合结果：,第22页/共48页,曲线拟合问题最常用的解法线性最小二乘法的基本思路,第一步:先选定一组函数 r1(x), r2(x), rm(x), mn, 令 f(x)=a1r1(x)+a2r2(x)+ +amrm(x) （1） 其中 a1,a2, am 为待定系数。,第二
6、步: 确定a1,a2, am 的准则（最小二乘准则）： 使n个点（xi,yi) 与曲线 y=f(x) 的距离i 的平方和最小 。,记,问题归结为，求 a1,a2, am 使 J(a1,a2, am) 最小。,第23页/共48页,线性最小二乘法的求解：预备知识,超定方程组：方程个数大于未知量个数的方程组,超定方程一般是不存在解的矛盾方程组。,如果有向量a使得 达到最小，则称a为上述超定方程的最小二乘解。,第24页/共48页,定理：当RTR可逆时，超定方程组（3）存在最小二乘解，且即为方程组 RTRa=RTy的解：a=(RTR)-1RTy,所以，曲线拟合的最小二乘法要解决的问题，实际上就是求以下超
7、定方程组的最小二乘解的问题。,线性最小二乘法的求解,第25页/共48页,线性最小二乘拟合 中函数r1(x), rm(x)的选取,1. 通过机理分析建立数学模型来确定 f(x)；,2. 将数据 (xi,yi) i=1, n 作图，通过直观判断确定 f(x)：,第26页/共48页,1.线性最小二乘拟合,2.非线性最小二乘拟合,二、用MATLAB解拟合问题,第27页/共48页,1) 作多项式f(x)=a1xm+ +amx+am+1拟合,可利用已有程序:,a=polyfit(x,y,m),2) 对超定方程组,3)多项式在x处的值y可用以下命令计算： y=polyval（a，x）,1. 用MATLAB作
8、线性最小二乘拟合,第28页/共48页,例 对下面一组数据作二次多项式拟合,第29页/共48页,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,解法1用解超定方程的方法,2）计算结果： = -9.8108 20.1293 -0.0317,第30页/共48页,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; A=poly
9、fit(x,y,2) z=polyval(A,x); plot(x,y,k+,x,z,r) %作出数据点和拟合曲线的图形,2）计算结果： = -9.8108 20.1293 -0.0317,解法2用多项式拟合的命令,第31页/共48页,建模案例：黄河小浪底调水调沙问题,2004年6月至7月黄河进行了第三次调水调沙试验，特别是首次由小浪底、三门峡和万家寨三大水库联合调度，采用接力式防洪预泄放水，形成人造洪峰进行调沙试验获得成功。整个试验期为20多天，小浪底从6月19日开始预泄放水，至到7月13日恢复正常供水结束。小浪底水利工程按设计拦沙量为75.5亿立方米，在这之前，小浪底共积泥沙达14.15亿
10、吨。,第32页/共48页,建模案例：黄河小浪底调水调沙问题,这次调水调沙试验一个重要的目的就是由小浪底上游的三门峡和万家寨水库泄洪，在小浪底形成人造洪峰，冲刷小浪底库区沉积的泥沙，在小浪底水库开闸泄洪以后，从6月27日开始三门峡水库和万家寨水库陆续开闸放水，人造洪峰于29日先后到达小浪底，7月3日达到最大流量2700 ，使小浪底水库的排沙量也不断地增加。表7-1是由小浪底观测站从6月29日到7月10日检测到的试验数据。,第33页/共48页,第34页/共48页,1. 模型假设,1) 水流量和排沙量都是连续的，不考虑上游泄洪所带来的含沙量和外界带来的含沙量。,2) 时间是连续变化的，所取时间点依次
11、为1,2,3，,24,单位时间为12h。,第35页/共48页,2. 模型建立与求解,问题一：给出估算任意时刻的排沙量及总排沙量的方法,根据试验数据，要计算任意时刻的排沙量，就要确定出排沙量随时间变化的规律，可以通过插值来实现，考虑到实际中排沙量应该是时间的连续函数，为了提高精度，我们采用三次样条函数来进行插值。,第36页/共48页,1）建立模型,第37页/共48页,第38页/共48页,模型求解,第39页/共48页,计算结果,第40页/共48页,问题二：确定排沙量与水流量的变化关系。,第41页/共48页,散点图,第42页/共48页,误差计算方法,第43页/共48页,第一阶段多项式拟合,第44页/共48页,第45页/共48页,第二阶段多项式拟合,第46页/共48页,第47页/共48页,谢谢您的观看！,第48页/共48页,
