插值与拟合分析.pptVIP

求得：a = 93.718387 b1 = -0.089032 b2 = -1.166774 所以该种材料的气孔率可近似为 Vp = 93.718387 - 0.089032 T - 1.166774 t §2.2.2 多项式拟合 y(x) = b0 + b1x + b2x2 +…+ bnxn 二次多项式拟合：y = b0 + b1x + b2x2 §2.2.3 指数拟合 有些指数拟合可以转化成线性拟合 两边取对数，得： 若令 则 这是个线性模型。 如： 数据拟合中，如何选择数学模型很重要，通常需要根据物理意义，或所给数据的分布情况来选取合适的数学模型。 y = 4.14896 + 1.14359 t －0.04832 t2 t 1 2 3 4 6 8 10 12 14 16 y 4 6.41 8.01 8.79 9.53 9.86 10.33 10.42 10.53 10.61 《插值与拟合》部分的参考书目： 《数值分析》 “插值与拟合” 所用软件：Matlab “拟合”所用软件： Origin 第二章 数据处理技术 ——插值与拟合 §2.1 插值方法 引例 这就是本章要讨论的“插值问题” 深度（m） 466 741 950 1422 1634 水温（oC） 7.04 4.28 3.40 2.54 2.13 【例2-1】已经测得在某处海洋不同深度处的水温如下： 根据这些数据，希望合理地估计出其它深度（如500米， 600米，1000米…）处的水温。 函数插值的定义 粗略地说，函数插值是对函数的离散数据建立简单的数学模型。 设 y = f(x) 是区间 [a,b] 上的连续函数，记作 。已知 f 在 [a,b] 上n+1个互异点 a≤x0,x1,…,xn-1,xn≤b, xi ≠ xj (i≠j) 处的值 yi = f(xi), i=0,1,2, …,n 若有不超过n次的多项式 满足 (2.1) 则称Ln(x) 为函数 f(x) 在区间 [a,b] 上通过点列 的插值多项式。 其中，[a,b] 称为插值区间， 称为插值节点， 求函数值f(x) 的点x (x≠xi) 称为插值点， f(x) 称为被插函数， Ln(x) 称为插值函数， 式（2.1)称为插值条件。 简单地说，插值法就是根据一组数据点（x1, y1）,（x2, y2）,…, （xn, yn）建立一个便于计算的初等函数或 曲线 y = f(x)，使它通过这些给定的数据点： f(x1) = y1， f(x2) = y2，…， f(xn) = yn 用这种方法所得的近似公式叫插值公式，已知的数据点叫节点。 ? Lagrange插值 ? Newton插值 ? Hermite插值 ? 样条插值 插值方法 拉格朗日插值法 怎样构造插值函数 Ln(x) ？ 从 n = 1，n = 2 推广到一般情况。 n = 1 线性插值 (2.2) (2.3) n = 2 抛物线插值 n = 2 时，构造通过三个点 (x0,y0)， (x1,y1) 和 (x2,y2) 的多项式如下： 这样的L2(x) 满足插值条件L2(x0)= y0， L2(x1)= y1， L2(x2)= y2。 它的几何意义是通过三个插值点的抛物线。 (2.4) 式 2-4 也可写成： 其中 (2.5) 推广到一般情形——拉格朗日插值公式 其中 k = 0, 1, 2, …, n 拉格朗日 插值基函数 拉格朗日 插值公式 (2-6) (2-7) 【例 2-2】已知函数 f(x)的三个点 (0,1), (-1,5) 和 (2,-1), 写出 拉格朗日插值基函数，并用公式(2-5)求2次插值多项式L2(x)。 解： 这里 n = 2，由式 (2-5) 得知三个拉格朗日插值基函数为 x0 = 0 x1 = -1 x2 = 2 y0 = 1 y1 = 5 y2 = -1 代入式 (2-5)得 【例 2-3】已知CO2在1673K、1773K、1873K时的热容分别为2.3136、2.3354、2.3563 kJ/(Nm3?K)。求CO2在1708Ｋ时的热容。 解：令 抛物线插值法 分别将值带入l0(x)、l1(x) 和 l2(x)。 x0 = 1673 x1 = 1773 x2 = 1873 y0 = 2.3136 y1 = 2.3354 y2 = 2.3563 【例 2-4】已知某液体的粘度－温度关系，E330＝60 Pa·S， E350＝30 Pa·