← 返回目录
第 1 章:插值法
1.1 Lagrange 插值
构造辅助函数 \(\varphi(t) = f(t) - L_n(t) - K\cdot\omega_{n+1}(t)\),其中 \(\omega_{n+1}(t)=\prod_{j=0}^n(t-x_j)\),\(K\) 由 \(\varphi(x)=0\) 确定,即 \(K = \frac{f(x)-L_n(x)}{\omega_{n+1}(x)}\)。
\(\varphi\) 在 \(x_0,\dots,x_n,x\) 共 \(n+2\) 个点处为零,由 Rolle 定理反复应用,\(\exists\,\xi\) 使 \(\varphi^{(n+1)}(\xi)=0\)。
\(\varphi^{(n+1)}(\xi) = f^{(n+1)}(\xi) - 0 - K\cdot(n+1)! = 0\),故 \(K = \frac{f^{(n+1)}(\xi)}{(n+1)!}\)。
取 \(f\equiv1\),其插值多项式为常数 1,由唯一性 \(\sum\ell_k\equiv1\)。
类似地,取 \(f(x)=x^m\)(\(m\le n\)),其 Lagrange 插值多项式就是 \(x^m\) 本身(次数 \(\le n\),在节点处精确),故 \(\sum_{k=0}^n x_k^m\ell_k(x)=x^m\)。
1.2 Newton 差商插值
\(f[x_0,x_1] =\)
\(f[x_0,\dots,x_k] =\)
\[f[x_0,x_1] = \frac{f(x_1)-f(x_0)}{x_1-x_0}\]
\[f[x_0,\dots,x_k] = \frac{f[x_1,\dots,x_k]-f[x_0,\dots,x_{k-1}]}{x_k - x_0}\]
用归纳法或直接计算:\(f[x_0,\dots,x_n] = \frac{(-1)^n}{\prod_{k=0}^n(x_k-a)}\)。
验证:\(n=1\):\(f[x_0,x_1]=\frac{\frac{1}{x_1-a}-\frac{1}{x_0-a}}{x_1-x_0}=\frac{-(x_1-x_0)}{(x_1-a)(x_0-a)(x_1-x_0)}=\frac{-1}{(x_0-a)(x_1-a)}\),与公式一致,且关于 \(x_0,x_1\) 对称。
对称性:\(\frac{(-1)^n}{\prod(x_k-a)}\) 显然关于 \(x_0,\dots,x_n\) 的任意置换不变。
1.3 Hermite 插值
在节点 \(x_i\) 处要求 \(H^{(j)}(x_i)=f^{(j)}(x_i)\),\(j=0,\dots,m_i-1\)(\(m_i\) 为重数)。总条件数 \(N=\sum m_i\),\(H\) 为次数 \(\le N-1\) 的多项式。
1.4 三次样条插值
令 \(h=g-S\),则 \(h(x_i)=0\)。展开:\(\int(g'')^2=\int(S'')^2+2\int S''h''+\int(h'')^2\)。
对交叉项分部积分两次:\(\int_a^b S''h''=[S''h']_a^b-[S'''h]_a^b+\int_a^b S^{(4)}h\,dx\)。
三项均为零:\(S''(a)=S''(b)=0\)(自然边界);\(h(x_i)=0\) 且 \(S'''\) 分段常数故 \([S'''h]_a^b=0\);\(S^{(4)}\equiv0\) 逐段成立。
故 \(\int(g'')^2=\int(S'')^2+\int(h'')^2\ge\int(S'')^2\)。
1.5 综合练习
由差商递推定义:\(f[x_0,\dots,x_{n-1},t,t+h] = \frac{f[x_0,\dots,x_{n-1},t+h]-f[x_0,\dots,x_{n-1},t]}{h} = \frac{g(t+h)-g(t)}{h}\)。
令 \(h\to0\),左边趋于 \(f[x_0,\dots,x_{n-1},t,t]\)(差商关于节点连续),右边趋于 \(g'(t)\),故 \(g'(t)=f[x_0,\dots,x_{n-1},t,t]\)。
令 \(r=p-q\),则 \(\deg r\le 2n+1\),且 \(r(x_i)=0\),\(r'(x_i)=0\),\(i=0,\dots,n\)。
故每个 \(x_i\) 是 \(r\) 的至少二重零点,\(r\) 有至少 \(2(n+1)=2n+2\) 个零点(计重数)。
但 \(\deg r\le 2n+1 < 2n+2\),故 \(r\equiv0\),即 \(p\equiv q\)。