主题目录

多项式插值与点值方法

拓展更新于 2026-09-29

本页目录

问题与唯一性

给定 n+1n+1 个互异节点 x0,…,xnx_0,\ldots,x_n 及指定值 y0,…,yny_0,\ldots,y_n,希望找次数不超过 nn 的多项式 PP,满足 P(xi)=yiP(x_i)=y_i。节点和值可以为实数,也可以为复数。

若 P,QP,Q 都符合要求,则 P−QP-Q 的次数至多为 nn,却在这 n+1n+1 个节点为零,故恒为零。这证明至多有一个解;下面还要构造解,才能证明存在。不要把“至多一个”和“存在唯一一个”混为一谈。

拉格朗日插值公式及证明

对 i=0,…,ni=0,\ldots,n,定义基多项式

Li(x)=∏0≤j≤nj≠ix−xjxi−xj.L_i(x)=\prod_{\substack{0\le j\le n\\j\ne i}}\frac{x-x_j}{x_i-x_j}.

所有分母非零,因为节点互异。LiL_i 在自己的节点 xix_i 取值 11,在其他节点取值 00。因此

P(x)=∑i=0nyiLi(x)P(x)=\sum_{i=0}^n y_iL_i(x)

次数不超过 nn,且代入 xkx_k 时只有第 kk 项保留下来,得到 P(xk)=ykP(x_k)=y_k。结合唯一性,这就是所求多项式。n=0n=0 时空乘积取 11,公式给出常数 y0y_0。

插值结果的实际次数可能小于 nn,例如全部 yiy_i 相同。节点重复时公式没有意义:重复的相同条件不增加信息,重复节点给出不同数值则根本无解。

不展开全部系数也能求值

已知次数不超过二的 PP 满足 P(0)=1,P(1)=2,P(2)=5P(0)=1,P(1)=2,P(2)=5。直接插值得

P(x)=(x−1)(x−2)2−2x(x−2)+52x(x−1)=x2+1.P(x)=\frac{(x-1)(x-2)}2-2x(x-2)+\frac52x(x-1)=x^2+1.

若只求 P(3)P(3),可先把 33 代入三个基多项式:系数依次为 1,−3,31,-3,3,故 P(3)=1−6+15=10P(3)=1-6+15=10。计算量来自点值的线性组合,不必先设一般式再解完整方程组。

比较最高次项与构造含已知根的多项式

若 P(x)=anxn+⋯+a0P(x)=a_nx^n+\cdots+a_0,允许 an=0a_n=0,比较插值式中的 xnx^n 系数得到

an=∑i=0nP(xi)∏j≠i(xi−xj).a_n=\sum_{i=0}^n\frac{P(x_i)}{\prod_{j\ne i}(x_i-x_j)}.

特别地,取 P(x)=xkP(x)=x^k,当 0≤k<n0\le k<n 时该和为零,k=nk=n 时为一。k=0k=0 表示常数多项式 11,并不要求对零节点解释 000^0。这些关系可用于点值估计和求系数。

若点值具有其他共同结构,先构造一个含已知根的多项式往往更短。例如已知 P(ti)=1/tiP(t_i)=1/t_i 且 ti≠0t_i\ne0,则 xP(x)−1xP(x)-1 在各 tit_i 为零。先核对其次数,再提出这些一次因式,最后用其他点求常数。完整例子见点值构造。若次数仍有余量,商可能是非常数多项式,不能直接写成常数倍。

插值与余式、差分的联系

设 D(x)=∏i=0n(x−xi)D(x)=\prod_{i=0}^n(x-x_i),节点互异。任意多项式 FF 除以 DD 的余式 RR 次数至多为 nn,且 R(xi)=F(xi)R(x_i)=F(x_i)。所以对 FF 的节点值做插值,得到的正是余式,而不必等于 FF 本身。

节点等距时,可用有限差分与牛顿有限展开组织同一插值问题;重复节点则需要导数值等额外条件,见重根除式的余式。有限个样本总能插值,不代表原函数在未给出的点也服从这个多项式,更不能凭此断言数列的无限规律。

关联阅读

多项式目录 · 全站目录