插值算法之拉格朗日插值总结.pdfVIP

记一下拉格朗日插值公式的推导和一些要点【这里说的都是二维插值，多维上的以此类推】 1、插值问题：在做实验的过程中，往往得到一堆离散的数据，现在想用数学公式模拟这堆 离散数据。 怎么办，数学家们提出了插值问题。 插值问题的提法是这样的给定一堆数据点 (x0, y0), (x1, y1), (x2, y2)...(xn, yn) ，要求一个函数 y = f(x) ，要求该函数经过上面所有的数据 点。 2 、多项式插值及其唯一性：在所有的函数中，多项式函数是最简单的函数，所以只要是人 就会想到用多项式函数来作为插值函数， 好， 以上给定了 n+1 个点， 现在要求一个 n 次多项 式 y = an * x^n + ... a1 * x + a0, 使它们经过这 n+1 个点；通过范德蒙行列式 和 克莱姆法 则，可以判定如果这 n+1 个点的 x 值各不相同，那么这个多项式是唯一的。结果唯一，但是 用直接法很不好求。现在用别的办法来求之。这就是：拉格朗日多项式 3 、拉格朗日多项式的构造，以四个点为例子进行说明 由于函数经过 4 个点 (x0, y0) ，(x1, y1) ，(x2, y2) ，(x3, y3) ，所以可以设函数为： f(x) = b0(x) * y0 + b1(x) * y1 + b2(x) * y2 + b3(x) * y3 注意： b0(x) ， ...，b3(x) 都是 x 的3次多项式，称之为拉格朗日插值基函数。 由于要求当 x 为 x0 时候， f(x) = y0, 所以最简单的做法就是让 b0(x0) = 1, b1(x0) = b2(x0) = b3(x0) = 0; 同理可知，在 x1 ，x2 ，x3 点上，插值基函数的值构造如下： b0(x) b1(x) b2(x) b3(x) x=x0 1 0 0 0 x=x1 0 1 0 0 x=x2 0 0 1 0 x=x3 0 0 0 1 问题 1、根据这些值来确定 b0(x) 的表达式， 由于 b0(x1) = b0(x2) = b0(x3) = 0 ，所以 x1, x2, x3 是 b0(x) 的零点，由于 b0(x) 是三次多项 式，所以设 b0(x) = c0 * (x-x1) * (x-x2) * (x-x3) 由于 b0(x0) = 1, 所以 1 = c0 * (x0-x1) * (x0-x2) * (x0-x3) 得到 c0 = 1/[(x0-x1)(x0-x2)(x0-x3)] 所以： b0(x) = (x-x1)*(x-x2)*(x-x3)/[(x0-x1)*(x0-x2)*(x0-x3)] 同理可求 b1(x) 、b2(x) ，略 问题 2 、根据上面的表格说明插值基函数的一个性质：无论 x 取和值，它们的和都为 1. 【这 个叫做调和函数】 以3 次为例子说明： 将上述表格的每一行分别相加， 得到的事函数： g(x) = b0(x) + b1(x) + b2(x) + b3(x) 在 x0, x1, x2, x3 的值，都为 1. b0(x) + b1(x) + b2(x) + b3(x) x=x0 1+0+0+0 = 1 x=x1 0+1+0+0 = 1 x=x2