牛顿插值多项式.ppt

分段线性插值 Runge现象 给定函数 取等距插值节点 建立10次插值多项式 三次Hermite插值 问题: 已知函数f(x)在两个节点 x0, x1上 的函数值分别为 y0, y1 及一阶导数值分别为 m0, m1 构造一个插值函数H3(x) ，使满足条件 1°H3(x)是次数?3的多项式 2°H3(x0)=y0 , H3(x1)=y1 ,H3 (x0)=m0 H3(x1)=m1 称这类插值问题为三次Hermite插值问题. 首先求做三次多项式h0(x),h1(x),h0(x),h1(x),使其满足 h0(x0)=1, h0(x1)=0, h0(x0)=0, h0(x1)=0 h1(x0)=0, h1(x1)=1, h1(x0)=0, h1(x1)=0 h0(x0)=0, h0(x1)=0, h0(x0)=1, h0(x1)=0 h1(x0)=0, h1(x1)=0, h1(x0)=0, h1(x1)=1 设 由h0(x0)=1 ，得a=1, 再由h0(x0)=0 ，得 ，于是 同理有 设 由h0(x0)=1 ，得a=1, 于是 同理有 显然H3(x)= y0h0(x)+y1h1(x)+m0h0(x)+m1h1(x) 其中 三次Hermite插值多项式的余项 定理3 设H3(x) 是以x0, x1为插值节点的三次Hermite插值多项式，若f(x)?C3[a, b] ， f(4)(x)在(a, b)上存在，其中 [a, b]是包含(x0, x1)的任一区间，则对任意给定的x?[a,b] ，总存在一点??(a, b)（?依赖于x）使 （2）多项式拟合 有时所给数据点的分布并不一定近似地呈一条直线,这时仍用直线拟合显然是不合适的,可用多项式拟合。对于给定的一组数据 寻求次数不超过m (mN ) 的多项式， 来拟合所给定的数据，与线性拟合类似，使偏差的 平方和 为最小 由于Q可以看作是关于 ( j=0,1,2,…, m)的多元函数, 故上述拟合多项式的构造问题可归结为多元函数的极值问题。令 得 即有 这是关于系数 的线性方程组，通常称为正规方程组（法方程组）。可以证明，正规方程组有惟一解。 例设某实验数据如下： 1 2 3 4 5 6 0 1 2 3 4 5 5 2 1 1 2 3 用最小二乘法求一个多项式拟合这组数据 解：将已给数据点描在坐标系中，可以看出这些点 接近一条抛物线，因此设所求的多项式为 由法方程组， 经计算得 N=6, 其法方程组为 解之得 所求的多项式为 （3）可化为线性拟合的非线性拟合 有些非线性拟合曲线可以通过适当的变量替换转化为线性曲线，从而用线性拟合进行处理，对于一个实际的曲线拟合问题，一般先按观测值在直角坐标平面上描出散点图，看一看散点的分布同哪类曲线图形接近，然后选用相接近的曲线拟合方程。再通过适当的变量替换转化为线性拟合问题，按线性拟合解出后再还原为原变量所表示的曲线拟合方程。 表5-4列举了几类经适当变换后化为线性拟合求解的曲线拟合方程及变换关系 表5-4 曲线拟合方程 变换关系 变换后线性拟合方程 几种常见的数据拟合情况。图 ( a ) 表示数据接近于 直线，故宜采用线性函数 拟合；图(b)数 据分布接近于抛物线。可采拟合；二次多项式 拟合； (a) (b) 图 ( c ) 的数据分布特点是开始曲线上升较快随后逐 渐变慢,宜采用双曲线型函数 或指数型函 数 图 ( d ) 的数据分布特点是开始曲线下降快,随后逐渐变慢,宜采用 或 或 等数据拟合。 ( c ) ( d ) 例 设某实验数据如下: 1 2 3