lesson9插值与数据拟合建模.pptx

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