lesson9（插值与数据拟合建模）.ppt
插值与数据拟合建模,插值与数据拟合就是通过一些已知数据去确定某类函数的参数或寻找某个近似函数，使所得的函数与已知数据具有较高的精度，并且能够使用数学分析的工具分析数据所反映的对象的性质 几种常用的方法： 1、一般插值法 2、样条插值法 3、最小二乘曲线 4、曲面的拟合,一、基本概念,1、插值问题： 不知道某一函数f(x)在待定范围a,b上 的具体表达式，而只能通过实验测量得到该 函数在一系列点ax1, x2 , ., xn b上的值 y0, y1, y2, ., yn，需要找一个简单的函数P(x) 来近似地代替f(x)，要求满足： P(xi)=yi (i=1,2,.,n) 2、插值函数：P(x) ， 3、插值法：求插值函数P(x)的方法,二、常用插值函数,1、多项式函数 2、样条函数,1、多项式插值方法,（1）n次代数插值 （2）拉格朗日插值 几点说明： （1）拉格朗日插值基函数仅与节点有关，而 与被插值函数f(x)无关。 （2）拉格朗日插值多项式仅由数对(xi,yi)(i 1,2,n)确定，而与数对排列次序无关。 （3）多项式插值除了上述插值法外还有其它 插值法，如newton插值法、hermite插 值法等。,二、常用插值函数,1、多项式函数 2、样条函数,2、样条插值方法,（1）样条函数m次半截幂函数,（2）k次B样条或k次基本样条函数的定义,（一）广泛使用的样条函数,（1）广泛采用：二次样条、三次样条及B样条。 （2）力学意义： A:二次样条在力学上解释为集中力偶作用 下的弹性细梁挠度曲线。 B:弹性细梁受集中载荷作用形成的挠度曲 线，在小挠度的情况下，恰好表示为三 次样条函数，集中载荷的作用点，恰好 就是三次样条函数的节点。,（1）二次样条的定义,设a,b 的一个划分：a=x0x1, x2 , ., xn= b，函数f ( x )各节点的值分别为： f ( xi )=yi (i=1,2,.,n) 如果二次样条函数：,满足： S ( xi )=yi (i=1,2,.,n),（2）三次样条函数的定义,设a,b 的一个划分：a=x0x1, x2 , ., xn= b， 函数f ( x )各节点的值分别为： f ( xi )=yi (i=1,2,.,n) 如果三次样条函数：,3,满足： S ( xi )=yi (i=1,2,.,n),三、曲线拟合最小二乘法,最小二乘法的一般提法是： 对给定的一组数据(x i , y i )(i0，1，2，m)，要求在函数类=1， 2， n中找一个函数yS*(x)，使误差平方和最小：,其中：,例1,在一个SCS系统中，根据实验所得输出信号与时间关系如下表所示，求输出信号 y与时间t的拟台曲线yf ( t )。,案例1 估计水箱的水流量模型,长度单位：E(3024cm) 容积单位：G(=3.785L(升) 某些镇的用水管理机构需估计公众的用水速度(单位是G/h)和每天总用水量的数据许多地方没有测量流入或流出水箱流量的设备，而只能测量水箱中的水位(误差不超过5%)当水箱水位低于某最低水位L时，水泵抽水，灌入水箱，直至水位达到最高水位H为止。但这也无法测量水泵的流量，因此在水泵启动时不易建立水箱中水位和水泵工作时用水量之间关系。水泵一天灌水12次，每次约2h。试估计在任意时刻(包括水泵灌水期间)t流出水箱的流量f(t)，并估计一天的总用水量。,案例1数据,表中给出了某镇中某一天的真实用水数据，表中测量时间以s为单位，水位以10-2 E为单位例如3316s以后，水箱中的水深降至3110 E时，水泵自动启动把水输入水箱；而当水位回升至35.5E时水泵停止工作。,1、模型假设（1）,1）除了问题中特别说明的数据以外，其它给 定的数据其测量误差不超过5 2）影响水箱流量的唯一因素是该镇公众对水 的普通需要所给的数据反映该镇在通常 情况下一天的用水量，不包括任何非常情 况如水箱中水的短缺、水管破裂、自然 灾害等 一个物理定理,1、模型假设（2）,3）水泵的灌水速度为常数，不随时间变化也 不是巳灌水量的函数，因此假设水泵大约 在水位27E时开始灌水，在水位355E时 停止灌水同时假设水泵不会损坏或不需 要维护 4）从水箱中流出的最大流速小于水泵的灌水 速度为了满足公众的用水需求不让水棺 中的水流尽,1、模型假设（3）,5）每天的用水量分布都是相似的因为公众 对水的消耗量是以全天的活动(诸如洗澡、 做饭、洗衣服等)为基础的，所以每日用水 类型是相似的 6）水箱的流水速度可用光滑曲线来近似每 个用户的用水需求量与整个镇的用水需求 量相比是微不足道的，而且它与整个镇需 求量的增减情况是不相似的,2、问题分析,3、水流量与时间关系,中心差分公式：,用来计算每一组的中间数据点的水流量。对每一组数据的前两个和最后两个数据点，采用向前和向后差分公式,4、数值决定的水流量-时间图,5、三次样条拟合,6、对水泵两段充水时间的处理,（1）第一次平均水流量：,（2）第二次平均水流量：,（3）平均水流量：,7、一天总用水量,8、检验,以不同的时间为起点得到的一天总用水量相差多少,一个物理定理,由托里查里(Torricelli)定律知，从水箱中流出水的最大速度与水位高的平方成正比对于所给的数据，其水位的最大高度为355E，最小高度为27E，因此对两个高度的最大流速比为 ，这个数字已很接近l，所以可以假定水位对流速没有影响类似地，还假设大气情况、温度变化等对水流速均无直接影响。,
