四、多项式插值与函数最佳逼近
设函数 y=f(x)y = f(x)y=f(x) 在区间 [a,b][a, b][a,b] 上有定义,且已知在点 a≤x0<x1<⋯<xn≤ba \le x_0 \lt x_1 \lt \cdots \lt x_n \le ba≤x0<x1<⋯<xn≤b 上的值 f(x0), f(x1), ⋯ , f(xn)f(x_0),\ f(x_1),\ \cdots,\ f(x_n)f(x0), f(x1), ⋯, f(xn), 若存在一个次数不超过 nnn 的多项式 pn(x)p_n(x)pn(x), 使
pn(xi)=f(xi) (i=0,1,2,⋯ ,n)(1.1)
p_n(x_i) = f(x_i)\ (i = 0, 1, 2,\cdots,n)\quad (1.1)
pn(xi)=f(xi) (i=0,1,2,⋯,n)(1.1)
成立,则称 pn(x)p_n(x)pn(x) 为 f(x)f(x)f(x) 的 nnn 次插值多项式, 式 (1.1)(1.1)(1.1) 为插值条件,点 x0, x1,⋯ , xnx_0, \ x_1,\cdots ,\ x_nx0, x1,⋯, xn 称为插值节点,称 f(x)f(x)f(x) 为被插值函数
- 满足插值条件 (1.1)(1.1)(1.1) 的 nnn 次多项式 pn(x)p_n(x)pn(x) 是存在唯一的
Lagrange 插值多项式及余项表示
- 📖Lagrange 插值多项式: Ln(x)=∑k=0nf(xk)lk(x)=∑k=0nf(xk)∏i=0,i≠knx−xixk−xiL_n(x)=\sum_{k=0}^nf(x_k)l_k(x)=\sum_{k=0}^nf(x_k)\prod_{i=0,i \ne k}^n\frac{x-x_i}{x_k-x_i}Ln(x)=∑k=0nf(xk)lk(x)=∑k=0nf(xk)∏i=0,i=knxk−xix−xi
- l0(x),l1(x),...,ln(x)l_0(x),l_1(x),...,l_n(x)l0(x),l1(x),...,ln(x) 线性无关,它是 n 次多项式空间 PnP_nPn 的一组基,而 1,x,x2,...,xn\pmb{1,x,x^2,...,x^n}1,x,x2,...,xn1,x,x2,...,xn1,x,x2,...,xn 也是其一组基。l0(x),l1(x),...,ln(x)l_0(x),l_1(x),...,l_n(x)l0(x),l1(x),...,ln(x) 称为 nnn 次 Lagrange 插值基函数
- 插值多项式余项:Rn(x)=f(x)−Ln(x)R_n(x)=f(x)-L_n(x)Rn(x)=f(x)−Ln(x)
- 📖设 f(n)(x)f^{(n)}(x)f(n)(x) 在 [a,b][a,b][a,b] 上连续,f(n+1)(x)f^{(n+1)}(x)f(n+1)(x) 在 (a,b)(a,b)(a,b) 内存在,x0,x1,...,xn∈[a,b]x_0,x_1,...,x_n \in [a,b]x0,x1,...,xn∈[a,b] 为互异节点,Ln(x)L_n(x)Ln(x) 是满足 pn(xi)=f(xi) (i=0,1,2,...,n)p_n(x_i)=f(x_i)\ (i=0,1,2,...,n)pn(xi)=f(xi) (i=0,1,2,...,n) 的插值多项式,则对 ∀x∈[a,b],∃ξ∈(a,b) (ξ\forall x \in [a,b],\exists \xi \in(a,b)\ (\xi∀x∈[a,b],∃ξ∈(a,b) (ξ 依赖于 x)x)x),使得
Rn(x)=f(x)−Ln(x)=f(n+1)(ξ)(n+1)!∏i=0n(x−xi)=f(n+1)(ξ)(n+1)!ωn+1(x)ξ=ξ(x)∈(min{x0, x1, ⋯ , xn}, max{x0, x1, ⋯ , xn}) R_n(x)=f(x)-L_n(x)=\frac{f^{(n+1)}(\xi)}{(n+1)!}\prod_{i=0}^n(x-x_i)=\frac{f^{(n+1)}(\xi)}{(n+1)!}\omega_{n+1}(x) \\ \xi = \xi(x) \in (\min\{ x_0,\ x_1,\ \cdots,\ x_n \},\ \max\{ x_0,\ x_1,\ \cdots,\ x_n \}) Rn(x)=f(x)−Ln(x)=(n+1)!f(n+1)(ξ)i=0∏n(x−xi)=(n+1)!f(n+1)(ξ)ωn+1(x)ξ=ξ(x)∈(min{x0, x1, ⋯, xn}, max{x0, x1, ⋯, xn})
- ξ\xiξ 一般不能求出,因此只能估计误差
- Lagrange 插值的缺点:当节点增加或减少时,插值多项式 Ln(x)L_n(x)Ln(x) 将发生变化,计算不变
例题
例1:已知函数 f(x)=sinx,x∈[0,π]f(x)=\sin{x},x \in [0, \pi]f(x)=sinx,x∈[0,π],以 x0=0,x1=π2,x2=πx_0=0,x_1=\frac{\pi}{2},x_2=\pix0=0,x1=2π,x2=π 为插值节点,求 f(x)f(x)f(x) 的 2 次插值多项式 L2(x)L_2(x)L2(x)。
解:由 Lagrange 插值多项式知
L2(x)=f(0)(x−π2)(x−π)(0−π2)(0−π)+f(π2)(x−0)(x−π)(π2−0)(π2−π)+f(π)(x−0)(x−π2)(π−0)(π−π2)=−4π2x(x−π)
\begin{aligned}
L_2(x)&=f(0)\frac{(x-\frac{\pi}{2})(x-\pi)}{(0-\frac{\pi}{2})(0-\pi)}+f(\frac{\pi}{2})\frac{(x-0)(x-\pi)}{(\frac{\pi}{2}-0)(\frac{\pi}{2}-\pi)}+f(\pi)\frac{(x-0)(x-\frac{\pi}{2})}{(\pi-0)(\pi-\frac{\pi}{2})} \\
&= -\frac{4}{\pi^2}x(x-\pi)
\end{aligned}
L2(x)=f(0)(0−2π)(0−π)(x−2π)(x−π)+f(2π)(2π−0)(2π−π)(x−0)(x−π)+f(π)(π−0)(π−2π)(x−0)(x−2π)=−π24x(x−π)
例2:设函数
f(x)=2π∫0xe−t2dt
f(x)=\frac{2}{\sqrt{\pi}}\int_0^xe^{-t^2}dt
f(x)=π2∫0xe−t2dt
的函数值已造成函数表。假设在区间 [4, 6] 上用线性插值计算 f(x)f(x)f(x) 的近似值,问会有多大的误差?
解:在 [4, 6] 上作 f(x)f(x)f(x) 的线性插值多项式 p1(x)p_1(x)p1(x),则
R1(x)=f(x)−p1(x)=12f′′(ξ)(x−x0)(x−x1), ξ∈[4,6]f(x)′=2πe−x2f(x)′′=−4xπe−x2f(x)′′′=4π(2x2−1)e−x2>0, x∈(4,6),则∣R2(x)∣≤12×∣f′′(4)∣×∣(5−4)(5−6)∣=0.508×10−6
\begin{aligned}
R_1(x) &= f(x)-p_1(x)=\frac{1}{2}f^{''}(\xi)(x-x_0)(x-x_1), \ \xi \in [4, 6] \\
f(x)^{'} &= \frac{2}{\sqrt{\pi}}e^{-x^2} \\
f(x)^{''} &= -\frac{4x}{\sqrt{\pi}}e^{-x^2} \\
f(x)^{'''} &= \frac{4}{\sqrt{\pi}}(2x^2-1)e^{-x^2} > 0, \ x \in (4, 6) ,则\\
|R_2(x)| &\le \frac{1}{2} \times |f^{''}(4)| \times |(5-4)(5-6)|=0.508 \times 10^{-6}
\end{aligned}
R1(x)f(x)′f(x)′′f(x)′′′∣R2(x)∣=f(x)−p1(x)=21f′′(ξ)(x−x0)(x−x1), ξ∈[4,6]=π2e−x2=−π4xe−x2=π4(2x2−1)e−x2>0, x∈(4,6),则≤21×∣f′′(4)∣×∣(5−4)(5−6)∣=0.508×10−6
差商及 Newton 插值多项式
差商
📖设已知函数 f(x)f(x)f(x) 在 n+1n + 1n+1 个互异节点 x0, x1,⋯ , xnx_0,\ x_1, \cdots ,\ x_nx0, x1,⋯, xn 上的函数值为 f(x0), f(x1),⋯ , f(xn)f(x_0),\ f(x_1), \cdots,\ f(x_n)f(x0), f(x1),⋯, f(xn), 称
0阶差商:f[xi]=f(xi)1阶差商:f[xi,xj]=f(xj)−f(xi)xj−xi2阶差商:f[xi,xj,xk]=f[xj,xk]−f[xi,xj]xk−xik阶差商:f[x0,x1,...,xk−1,xk]=f[x1,...,xk]−f[x0,...,xk−1]xk−x0
\begin{aligned}
& 0阶差商:f[x_i] = f(x_i) \\
& 1阶差商:f[x_i, x_j] = \frac{f(x_j)-f(x_i)}{x_j - x_i} \\
& 2阶差商:f[x_i, x_j, x_k] = \frac{f[x_j, x_k] - f[x_i, x_j]}{x_k - x_i} \\
& k阶差商:f[x_0, x_1, ..., x_{k-1}, x_k] = \frac{f[x_1, ..., x_k] - f[x_0, ..., x_{k-1}]}{x_k - x_0}
\end{aligned}
0阶差商:f[xi]=f(xi)1阶差商:f[xi,xj]=xj−xif(xj)−f(xi)2阶差商:f[xi,xj,xk]=xk−xif[xj,xk]−f[xi,xj]k阶差商:f[x0,x1,...,xk−1,xk]=xk−x0f[x1,...,xk]−f[x0,...,xk−1]
- 😄kkk 阶差商可表示成函数值 f(x0), f(x1),⋯ , f(xk)f(x_0),\ f(x_1),\cdots,\ f(x_k)f(x0), f(x1),⋯, f(xk) 的线性组合,即
f[x0,x1,...,xk]=∑m=0kf(xm)∏i=0,i≠mk(xm−xi) f[x_0,x_1,...,x_k]=\sum_{m=0}^k \frac{f(x_m)}{\prod_{i=0,i\ne m}^k(x_m-x_i)} f[x0,x1,...,xk]=m=0∑k∏i=0,i=mk(xm−xi)f(xm)
- 😄kkk 阶差商 f[x0,...,xk]f[x_0, ..., x_k]f[x0,...,xk] 与节点的次序无关
- 📖kkk 阶差商和 kkk 阶导数之间关系
f[x0,x1,...,xk]=f(k)(η)k!, η∈(min{x0,...,xk},max{x0,...,xk}) f[x_0,x_1,..., x_k] = \frac{f^{(k)}(\eta)}{k!},\ \eta \in(min\{x_0,...,x_k\}, max\{x_0, ..., x_k\}) f[x0,x1,...,xk]=k!f(k)(η), η∈(min{x0,...,xk},max{x0,...,xk})
Newton 插值多项式
Ln(x)=f(x0)+f[x0,x1](x−x0)+f[x0,x1,x2](x−x0)(x−x1)+...+f[x0,x1,...,xn](x−x0)(x−x1)...(x−xn−1) \begin{aligned} L_n(x) &= f(x_0)+f[x_0, x_1](x-x_0)+f[x_0, x_1, x_2](x-x_0)(x-x_1) \\ &+...+f[x_0, x_1, ..., x_n](x-x_0)(x-x_1)...(x-x_{n-1}) \end{aligned} Ln(x)=f(x0)+f[x0,x1](x−x0)+f[x0,x1,x2](x−x0)(x−x1)+...+f[x0,x1,...,xn](x−x0)(x−x1)...(x−xn−1)
Hermite 插值
- 📖Hermite 插值多项式:给定 [a,b][a, b][a,b] 中的 n+1n+1n+1 个互异节点 xi(i=0,1,...,n)x_i(i=0,1,...,n)xi(i=0,1,...,n) 上的函数值和直到 mim_imi 阶的导数值 f(xi),f′(xi),...,f(mi)(xi)f(x_i),f^{'}(x_i),...,f^{(m_i)}(x_i)f(xi),f′(xi),...,f(mi)(xi)。令 m=∑i=0n(mi+1)−1m=\sum_{i=0}^n(m_i+1)-1m=∑i=0n(mi+1)−1,若存在一个次数不超过 mmm 的多项式 Hm(x)H_m(x)Hm(x) 使得下述内容成立,则 Hm(x)H_m(x)Hm(x) 为 Hermite 插值多项式
{Hm(x0)=f(x0),Hm′(x0)=f′(x0),...,Hm(m0)(x0)=f(m0)(x0)Hm(x1)=f(x1),Hm′(x1)=f′(x1),...,Hm(m1)(x1)=f(m1)(x1)...(4.1)Hm(xn)=f(xn),Hm′(xn)=f′(xn),...,Hm(mn)(xn)=f(mn)(xn) \left\{ \begin{aligned} &H_m(x_0) = f(x_0),H_m^{'}(x_0)=f^{'}(x_0),...,H_m^{(m_0)}(x_0)=f^{(m_0)}(x_0) \\ &H_m(x_1) = f(x_1),H_m^{'}(x_1)=f^{'}(x_1),...,H_m^{(m_1)}(x_1)=f^{(m_1)}(x_1) \\ &... && (4.1) \\ &H_m(x_n) = f(x_n),H_m^{'}(x_n)=f^{'}(x_n),...,H_m^{(m_n)}(x_n)=f^{(m_n)}(x_n) \\ \end{aligned} \right. ⎩⎪⎪⎪⎪⎨⎪⎪⎪⎪⎧Hm(x0)=f(x0),Hm′(x0)=f′(x0),...,Hm(m0)(x0)=f(m0)(x0)Hm(x1)=f(x1),Hm′(x1)=f′(x1),...,Hm(m1)(x1)=f(m1)(x1)...Hm(xn)=f(xn),Hm′(xn)=f′(xn),...,Hm(mn)(xn)=f(mn)(xn)(4.1)
- 满足上式的 mmm 次多项式 Hm(x)H_m(x)Hm(x) 存在唯一
- 📖设 Hm(x)H_m(x)Hm(x) 为满足 (4.1)(4.1)(4.1) 的 mmm 次 Hermite 插值多项式,f(x)f(x)f(x) 在包含 (n+1)(n + 1)(n+1) 个互异节点 x0, x1,⋯ , xnx_0,\ x_1,\cdots,\ x_nx0, x1,⋯, xn 的区间 [a,b][a, b][a,b] 上具有 mmm 阶连续导数, 在 (a,b)(a, b)(a,b) 内存在 (m+1)(m + 1)(m+1) 阶导数, 则对任意 x∈[a,b]x \in [a, b]x∈[a,b],存在 ξ∈(a,b)\xi \in (a, b)ξ∈(a,b),使得插值余项 Rm(x)R_m(x)Rm(x)
Rm(x)=f(x)−Hm(x)=f(m+1)(ξ)(m+1)!∏i=0n(x−xi)mi+1 R_m(x)=f(x)-H_m(x)=\frac{f^{(m+1)}(\xi)}{(m+1)!}\prod_{i=0}^n(x-x_i)^{m_i+1} Rm(x)=f(x)−Hm(x)=(m+1)!f(m+1)(ξ)i=0∏n(x−xi)mi+1
- 若 f∈Cn[a,b]f \in C^n[a, b]f∈Cn[a,b],a≤x0<x1<...<xn≤ba \le x_0 \lt x_1 \lt ... \lt x_n \le ba≤x0<x1<...<xn≤b,则有
f[x0,x1,...,xn]=∫⋯∫f(n)(t0x0+t1x1+...+tnxn)dt1...dtn f[x_0, x_1,..., x_n] = \int \cdots \int f^{(n)}(t_0x_0+t_1x_1+...+t_nx_n)dt_1...dt_n f[x0,x1,...,xn]=∫⋯∫f(n)(t0x0+t1x1+...+tnxn)dt1...dtn
上述具有 τn\tau_nτn 个积分 ∫\int∫,其中 t0=1−∑i=1ntit_0=1-\sum_{i=1}^nt_it0=1−∑i=1nti,τn={(t1,t2,...,tn)∣ti≥0\tau_n=\{(t_1, t_2, ..., t_n)|t_i \ge 0τn={(t1,t2,...,tn)∣ti≥0,∑i=1nti≤1}\sum_{i=1}^nt_i \le 1\}∑i=1nti≤1} 为 nnn 维单纯形
- 注意到被积函数是通过一元连续函数 f(n)(x)f^{(n)}(x)f(n)(x) 与 nnn 元线性连续函数 ∑i=0ntixi\sum_{i=0}^nt_ix_i∑i=0ntixi 复合而成,所以 f[x0,x1,...,xn]f[x_0, x_1, ..., x_n]f[x0,x1,...,xn] 是 x0,x1,...,xnx_0, x_1, ..., x_nx0,x1,...,xn 的连续函数,因此
f[x0,x0,...,x0]=limx1,x2,...,xk→x0f[x0,x1,...,xk]=limx1,x2,...,xk→x0f(k)(η)k!=f(k)(x0)k! \begin{aligned} f[x_0, x_0, ..., x_0]&=\lim_{x_1, x_2,..., x_k \rightarrow x_0}f[x_0, x_1,...,x_k] \\ &= \lim_{x_1, x_2,..., x_k \rightarrow x_0} \frac{f^{(k)}(\eta)}{k!} = \frac{f^{(k)}(x_0)}{k!} \end{aligned} f[x0,x0,...,x0]=x1,x2,...,xk→x0limf[x0,x1,...,xk]=x1,x2,...,xk→x0limk!f(k)(η)=k!f(k)(x0)
- 📖重节点插值的 mmm 次 Hermite 插值多项式为
Hm(x)=f(x0)+f[x0,x0](x−x0)+...+f[x0,...,x0](x−x0)m0+f[x0,...,x0,x1](x−x0)m0+1+...+f[x0,...,x0,x1,x1](x−x0)m0+1(x−x1)+f[x0,...,x0,x1,...x1,x2](x−x0)m0+1(x−x1)(m1+1)+...+f[x0,...,x0,...,xn−1,...,xn−1,xn,...,xn](x−x0)m0+1...(x−xn−1)mn−1+1(x−xn)mn \begin{aligned} H_m(x) &= f(x_0)+f[x_0, x_0](x-x_0)+...+f[x_0,..., x_0](x-x_0)^{m_0} \\ &+f[x_0,...,x_0,x_1](x-x_0)^{m_0+1} +...+f[x_0,...,x_0,x_1,x_1](x-x_0)^{m_0+1}(x-x_1) \\ &+f[x_0,...,x_0,x_1,...x_1,x_2](x-x_0)^{m_0+1}(x-x_1)^{(m_1+1)} \\ &+...+f[x_0,...,x_0,...,x_{n-1},...,x_{n-1},x_n,...,x_n](x-x_0)^{m_0+1}...(x-x_{n-1})^{m_{n-1}+1}(x-x_n)^{m_n} \end{aligned} Hm(x)=f(x0)+f[x0,x0](x−x0)+...+f[x0,...,x0](x−x0)m0+f[x0,...,x0,x1](x−x0)m0+1+...+f[x0,...,x0,x1,x1](x−x0)m0+1(x−x1)+f[x0,...,x0,x1,...x1,x2](x−x0)m0+1(x−x1)(m1+1)+...+f[x0,...,x0,...,xn−1,...,xn−1,xn,...,xn](x−x0)m0+1...(x−xn−1)mn−1+1(x−xn)mn
- 📖重节点插值余项为
Rm(x)=f(x)−Hm(x)=f(m+1)(ξ)(m+1)!∏i=0n(x−xi)mi+1 R_m(x)=f(x)-H_m(x)=\frac{f^{(m+1)}(\xi)}{(m+1)!}\prod_{i=0}^n(x-x_i)^{m_i+1} Rm(x)=f(x)−Hm(x)=(m+1)!f(m+1)(ξ)i=0∏n(x−xi)mi+1
例题
例1:若 n=0n = 0n=0,m0=km_0 = km0=k,即 1 个 k+1 重插值节点,其插值多项式为
Hk(x)=f(x0)+f′(x0)(x−x0)+...+f(k)(x0)k!(x−x0)k
H_k(x) = f(x_0)+f^{'}(x_0)(x-x_0)+...+\frac{f^{(k)}(x_0)}{k!}(x-x_0)^k
Hk(x)=f(x0)+f′(x0)(x−x0)+...+k!f(k)(x0)(x−x0)k
例2:求 4 次插值多项式 H(x)H(x)H(x),使得 H(0)=3,H′(0)=4,H(1)=5,H′(1)=6,H′′(1)=7H(0)=3,H^{'}(0)=4,H(1)=5,H^{'}(1)=6,H^{''}(1)=7H(0)=3,H′(0)=4,H(1)=5,H′(1)=6,H′′(1)=7
解:列表计算各点差商
| kkk | xkx_kxk | f(xk)f(x_k)f(xk) | 1阶差商 | 2阶差商 | 3阶差商 | 4阶差商 |
|---|---|---|---|---|---|---|
| 0 | 0 | 3 | 4 | -2 | 6 | −132-\frac{13}{2}−213 |
| 1 | 0 | 3 | 2 | 4 | −12-\frac{1}{2}−21 | |
| 2 | 1 | 5 | 6 | 72\frac{7}{2}27 | ||
| 3 | 1 | 5 | 6 | |||
| 4 | 1 | 5 |
故 H(x)=3+4(x−0)−2(x−0)2+6(x−0)2(x−1)−132(x−0)2(x−1)2H(x)=3+4(x-0)-2(x-0)^2+6(x-0)^2(x-1)-\frac{13}{2}(x-0)^2(x-1)^2H(x)=3+4(x−0)−2(x−0)2+6(x−0)2(x−1)−213(x−0)2(x−1)2
其插值余项 R4(x)=f(x)−H4(x)=f(5)(ξ)5!(x−0)2(x−1)3,ξ∈(a,b)R_4(x)=f(x)-H_4(x)=\frac{f^{(5)}(\xi)}{5!}(x-0)^2(x-1)^3, \xi \in (a, b)R4(x)=f(x)−H4(x)=5!f(5)(ξ)(x−0)2(x−1)3,ξ∈(a,b)
高阶插值的缺点及分段低次插值
高次插值的病态性质
当 xxx 在 ±1\plusmn 1±1 附近时,f(x)f(x)f(x) 的值和 L10(x)L_{10}(x)L10(x) 的值相差很大。这种现象称 Runge 现象。其实可以证明,f(x)f(x)f(x) 的 nnn 次插值多项式 Ln(x)L_n(x)Ln(x) 在 [−1,1][−1, 1][−1,1] 上不是一致收敛到 f(x)f(x)f(x)
分段线性插值
给定 f(x)f(x)f(x) 在 n+1n + 1n+1 个节点 a=x0<x1<⋯<xn=ba = x_0 \lt x_1 \lt \cdots \lt x_n = ba=x0<x1<⋯<xn=b 上的函数值:
| xxx | x0x_0x0 | x1x_1x1 | ⋯\cdots⋯ | xn−1x_{n-1}xn−1 | xnx_nxn |
|---|---|---|---|---|---|
| f(x)f(x)f(x) | f(x0)f(x_0)f(x0) | f(x1)f(x_1)f(x1) | ⋯\cdots⋯ | f(xn−1)f(x_{n-1})f(xn−1) | f(xn)f(x_n)f(xn) |
📖记 hi=xi+1−xih_i = x_{i+1} − x_ihi=xi+1−xi,h=max0≤i≤n−1hih = \max_{0 \le i \le n−1}h_ih=max0≤i≤n−1hi,在每个小区间 [xi,xi+1][x_i , x_{i+1}][xi,xi+1] 上作 f(x)f(x)f(x) 的线性插值
L1,i(x)=f(xi)+f[xi,xi+1](x−xi),x∈[xi,xi+1]
L_{1, i}(x) = f(x_i) + f[x_i, x_{i+1}](x-x_i), \quad x \in [x_i, x_{i+1}]
L1,i(x)=f(xi)+f[xi,xi+1](x−xi),x∈[xi,xi+1]
📖其误差为
f(x)−L1,i(x)=12f′′(ξi)(x−xi)(x−xi+1),ξi∈(xi,xi+1)
f(x) - L_{1, i}(x) = \frac{1}{2}f^{''}(\xi_i)(x-x_i)(x-x_{i+1}), \quad \xi_i \in (x_i, x_{i+1})
f(x)−L1,i(x)=21f′′(ξi)(x−xi)(x−xi+1),ξi∈(xi,xi+1)
从而有
maxxi≤x≤xi+1∣f(x)−L1,i(x)∣≤maxxi≤x≤xi+1∣12f′′(ξi)(x−xi)(x−xi+1)∣≤18hi2maxxi≤x≤xi+1∣f′′(x)∣
\max_{x_i \le x \le x_{i+1}}|f(x) - L_{1, i}(x)| \le \max_{x_i \le x \le x_{i+1}}|\frac{1}{2}f^{''}(\xi_i)(x-x_i)(x-x_{i+1})| \le
\frac{1}{8}h_i^2\max_{x_i \le x \le x_{i+1}}|f^{''}(x)|
xi≤x≤xi+1max∣f(x)−L1,i(x)∣≤xi≤x≤xi+1max∣21f′′(ξi)(x−xi)(x−xi+1)∣≤81hi2xi≤x≤xi+1max∣f′′(x)∣
令
L~1(x)={L1,0(x),x∈[x0,x1)L1,1(x),x∈[x1,x2)⋮L1,n−2(x),x∈[xn−2,xn−1)L1,n−1(x),x∈[xn−1,xn]
\tilde{L}_1(x) = \left\{
\begin{aligned}
&L_{1, 0}(x), && x \in [x_0, x_1) \\
&L_{1, 1}(x), && x \in [x_1, x_2) \\
\vdots \\
&L_{1, n-2}(x), && x \in [x_{n-2}, x_{n-1}) \\
&L_{1, n-1}(x), && x \in [x_{n-1}, x_n] \\
\end{aligned}
\right.
L~1(x)=⎩⎪⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎪⎧⋮L1,0(x),L1,1(x),L1,n−2(x),L1,n−1(x),x∈[x0,x1)x∈[x1,x2)x∈[xn−2,xn−1)x∈[xn−1,xn]
📖于是 L~1(xi)=f(xi), (i=0, 1, ..., n)\tilde{L}_1(x_i) = f(x_i),\ (i=0,\ 1,\ ...,\ n)L~1(xi)=f(xi), (i=0, 1, ..., n),即 L~1(x)\tilde{L}_1(x)L~1(x) 是 f(x)f(x)f(x) 的分段线性插值函数,其插值误差为
maxa≤x≤b∣f(x)−L~1(x)∣=max0≤i≤n−1maxxi≤x≤xi+1∣f(x)−L~1(x)∣=max0≤i≤n−1maxxi≤x≤xi+1∣f(x)−L1,i(x)∣≤max0≤i≤n−118hi2maxxi≤x≤xi+1∣f′′(x)∣≤18h2maxa≤x≤b∣f′′(x)∣
\begin{aligned}
\max_{a \le x \le b}|f(x)-\tilde{L}_1(x)| & = \max_{0 \le i \le n-1} \max_{x_i \le x \le x_{i+1}}|f(x)-\tilde{L}_1(x)| \\
&=\max_{0 \le i \le n-1} \max_{x_i \le x \le x_{i+1}}|f(x)-L_{1, i}(x)| \\
&\le \max_{0 \le i \le n-1} \frac{1}{8}h_i^2 \max_{x_i \le x \le x_{i+1}}|f^{''}(x)| \\
&\le \frac{1}{8}h^2 \max_{a \le x \le b}|f^{''}(x)|
\end{aligned}
a≤x≤bmax∣f(x)−L~1(x)∣=0≤i≤n−1maxxi≤x≤xi+1max∣f(x)−L~1(x)∣=0≤i≤n−1maxxi≤x≤xi+1max∣f(x)−L1,i(x)∣≤0≤i≤n−1max81hi2xi≤x≤xi+1max∣f′′(x)∣≤81h2a≤x≤bmax∣f′′(x)∣
- 只要 f(x)f(x)f(x) 在 [a,b][a, b][a,b] 上有 2 阶连续导数,当 h→0h \rightarrow 0h→0 时,余项一致趋于零,即分段线性插值函 数 L~1(x)\tilde{L}_1(x)L~1(x) 一致收敛于 f(x)f(x)f(x)
分段 3 次 Hermite
给定 f(x)f(x)f(x) 在 n+1n + 1n+1 个节点 a=x0<x1<⋯<xn=ba = x_0 \lt x_1 \lt \cdots \lt x_n = ba=x0<x1<⋯<xn=b 上的函数表
| xxx | x0x_0x0 | x1x_1x1 | ⋯\cdots⋯ | xn−1x_{n-1}xn−1 | xnx_nxn |
|---|---|---|---|---|---|
| f(x)f(x)f(x) | f(x0)f(x_0)f(x0) | f(x1)f(x_1)f(x1) | ⋯\cdots⋯ | f(xn−1)f(x_{n-1})f(xn−1) | f(xn)f(x_n)f(xn) |
| f′(x)f^{'}(x)f′(x) | f′(x0)f^{'}(x_0)f′(x0) | f′(x1)f^{'}(x_1)f′(x1) | ⋯\cdots⋯ | f′(xn−1)f^{'}(x_{n-1})f′(xn−1) | f′(xn)f^{'}(x_n)f′(xn) |
📖记 hi=xi+1−xih_i = x_{i+1} − x_ihi=xi+1−xi,h=max0≤i≤n−1hih = \max_{0 \le i \le n−1}h_ih=max0≤i≤n−1hi,在每个小区间 [xi,xi+1][x_i , x_{i+1}][xi,xi+1] 上利用数据
| xxx | xix_ixi | xi+1x_{i+1}xi+1 |
|---|---|---|
| f(x)f(x)f(x) | f(xi)f(x_i)f(xi) | f(xi+1)f(x_{i+1})f(xi+1) |
| f′(x)f^{'}(x)f′(x) | f′(xi)f^{'}(x_i)f′(xi) | f′(xi+1)f^{'}(x_{i+1})f′(xi+1) |
作 3 次 Hermite 插值
H3,i(x)=f(xi)+f′(xi)(x−xi)+f[xi,xi+1]−f′(xi)hi(x−xi)2+f′(xi+1)−2f[xi,xi+1]+f′(xi)hi2(x−xi)2(x−xi+1)
\begin{aligned}
H_{3, i}(x) = &f(x_i) + f^{'}(x_i)(x-x_i)+\frac{f[x_i, x_{i+1}] - f^{'}(x_i)}{h_i}(x-x_i)^2+ \\
&\frac{f^{'}(x_{i+1}) - 2f[x_i, x_{i+1}] + f^{'}(x_i)}{h_i^2}(x-x_i)^2(x-x_{i+1})
\end{aligned}
H3,i(x)=f(xi)+f′(xi)(x−xi)+hif[xi,xi+1]−f′(xi)(x−xi)2+hi2f′(xi+1)−2f[xi,xi+1]+f′(xi)(x−xi)2(x−xi+1)
其插值余项为
f(x)−H3,i(x)=f(4)(ξ)4!(x−xi)2(x−xi+1)2, ξ∈(xi,xi+1)
f(x) - H_{3, i}(x) = \frac{f^{(4)}(\xi)}{4!}(x-x_i)^2(x-x_{i+1})^2, \ \xi \in (x_i, x_{i+1})
f(x)−H3,i(x)=4!f(4)(ξ)(x−xi)2(x−xi+1)2, ξ∈(xi,xi+1)
于是
maxxi≤x≤xi+1∣f(x)−H3,i(x)∣≤14!hi416maxxi≤x≤xi+1∣f(4)(x)∣
\max_{x_i \le x \le x_{i+1}}|f(x) - H_{3,i}(x)| \le \frac{1}{4!}\frac{h_i^4}{16}\max_{x_i \le x \le x_{i+1}} |f^{(4)}(x)|
xi≤x≤xi+1max∣f(x)−H3,i(x)∣≤4!116hi4xi≤x≤xi+1max∣f(4)(x)∣
令
H~3(x)={H3,0(x)x∈[x0,x1)H3,1(x)x∈[x1,x2)⋯H3,n−2(x)x∈[xn−2,xn−1)H3,n−1(x)x∈[xn−1,xn]
\tilde{H}_3(x) =
\left\{
\begin{aligned}
&H_{3,0}(x) && x \in [x_0, x_1) \\
&H_{3,1}(x) && x \in [x_1, x_2) \\
&\cdots \\
&H_{3,n-2}(x) && x \in [x_{n-2}, x_{n-1}) \\
&H_{3,n-1}(x) && x \in [x_{n-1}, x_n] \\
\end{aligned}
\right.
H~3(x)=⎩⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎧H3,0(x)H3,1(x)⋯H3,n−2(x)H3,n−1(x)x∈[x0,x1)x∈[x1,x2)x∈[xn−2,xn−1)x∈[xn−1,xn]
则 H~3(xi)=f(xi)\tilde{H}_3(x_i) = f(x_i)H~3(xi)=f(xi),H~3′(xi)=f′(xi) (i=0,1,...,n)\tilde{H}_3^{'}(x_i) = f^{'}(x_i)\ (i=0,1,...,n)H~3′(xi)=f′(xi) (i=0,1,...,n),即 H~3(x)\tilde{H}_3(x)H~3(x) 满足插值条件,称 f(x)f(x)f(x) 的分段三次插值函数为 H~3(x)\tilde{H}_3(x)H~3(x)
其误差为
maxa≤x≤b∣f(x)−H~3(x)∣=max0≤i≤n−1maxxi≤x≤xi+1∣f(x)−H~3(x)∣=max0≤i≤n−1maxxi≤x≤xi+1∣f(x)−H3,i(x)∣≤max0≤i≤n−114!hi416maxxi≤x≤xi+1∣f(4)(x)∣≤1384h4maxa≤x≤b∣f(4)(x)∣
\begin{aligned}
\max_{a \le x \le b}|f(x)-\tilde{H}_3(x)| & = \max_{0 \le i \le n-1} \max_{x_i \le x \le x_{i+1}}|f(x)-\tilde{H}_3(x)| \\
&=\max_{0 \le i \le n-1} \max_{x_i \le x \le x_{i+1}}|f(x)-H_{3, i}(x)| \\
&\le \max_{0 \le i \le n-1} \frac{1}{4!}\frac{h_i^4}{16} \max_{x_i \le x \le x_{i+1}}|f^{(4)}(x)| \\
&\le \frac{1}{384}h^4 \max_{a \le x \le b}|f^{(4)}(x)|
\end{aligned}
a≤x≤bmax∣f(x)−H~3(x)∣=0≤i≤n−1maxxi≤x≤xi+1max∣f(x)−H~3(x)∣=0≤i≤n−1maxxi≤x≤xi+1max∣f(x)−H3,i(x)∣≤0≤i≤n−1max4!116hi4xi≤x≤xi+1max∣f(4)(x)∣≤3841h4a≤x≤bmax∣f(4)(x)∣
分段三次 Hermite 插值的余项和 f(x)f(x)f(x) 的 4 阶导数有关,当 f(x)f(x)f(x) 在 [a,b][a, b][a,b] 上有 4 阶连续导数,则有 H~3(x)→f(x)\tilde{H}_3(x) \rightarrow f(x)H~3(x)→f(x)
三次样条插值
端点处附加条件(边界条件)
- 已知两端点的一阶导数,即 S′(x0)=f′(x0),S′(xn)=f′(xn)S^{'}(x_0)=f^{'}(x_0),S^{'}(x_n)=f^{'}(x_n)S′(x0)=f′(x0),S′(xn)=f′(xn)
- 已知两端点的二阶导数,即 S′′(x0)=f′′(x0),S′′(xn)=f′′(xn)S^{''}(x_0)=f^{''}(x_0),S^{''}(x_n)=f^{''}(x_n)S′′(x0)=f′′(x0),S′′(xn)=f′′(xn)。当 S′′(x0)=0,S′′(xn)=0S^{''}(x_0)=0,S^{''}(x_n)=0S′′(x0)=0,S′′(xn)=0 时称为自然边界条件
- 周期边界条件,当 f(x0)=f(xn)f(x_0)=f(x_n)f(x0)=f(xn) 时,S′(x0+0)=S′(xn−0),S′′(x0+0)=S′′(xn−0)S^{'}(x_0+0)=S^{'}(x_n-0),S^{''}(x_0+0)=S^{''}(x_n-0)S′(x0+0)=S′(xn−0),S′′(x0+0)=S′′(xn−0)
样条插值的建立
TODO
最佳一致逼近
线性赋范空间
- 设 XXX 是一个集合,如果对 ∀x, y∈X, λ∈R\forall x,\ y \in X, \ \lambda \in R∀x, y∈X, λ∈R, 有 λx∈X, x+y∈X\lambda x \in X,\ x + y \in Xλx∈X, x+y∈X,则称 XXX 是线性空间
- 设 XXX 是 RRR 上的一个线性空间,φ1, φ2,⋯ , φm∈X\varphi_1,\ \varphi_2, \cdots ,\ \varphi_m \in Xφ1, φ2,⋯, φm∈X,如果存在不全为零的数 a1, a2,⋯ , am∈Ra_1,\ a_2, \cdots ,\ a_m \in Ra1, a2,⋯, am∈R 使得
a1φ1+a2φ2+⋯+amφm=0 a_1 \varphi_1 + a_2 \varphi_2 + \cdots + a_m \varphi_m = 0 a1φ1+a2φ2+⋯+amφm=0
则称 φ1, φ2,⋯ , φm\varphi_1,\ \varphi_2, \cdots ,\ \varphi_mφ1, φ2,⋯, φm 线性相关;否则称 φ1, φ2,⋯ , φm\varphi_1,\ \varphi_2, \cdots ,\ \varphi_mφ1, φ2,⋯, φm 线性无关
- 📖设 XXX 是一个线性空间。若对 ∀x∈X\forall x \in X∀x∈X,对应于实数,记为 ∥x∥\parallel x \parallel∥x∥ 且满足下面关系:
- ∀x∈X\forall x \in X∀x∈X,有 ∥x∥≥0\parallel x \parallel \ge 0∥x∥≥0,且 ∥x∥=0⇒x=0\parallel x \parallel = 0 \Rightarrow x=0∥x∥=0⇒x=0
- ∀λ∈R\forall \lambda \in R∀λ∈R, x∈X\ x \in X x∈X,有 ∥λx∥=∣λ∣∥x∥\parallel \lambda x \parallel = |\lambda|\parallel x \parallel∥λx∥=∣λ∣∥x∥
- x,y∈Xx, y \in Xx,y∈X,有 ∥x+y∥≤∥x∥+∥y∥\parallel x+y \parallel \le \parallel x \parallel + \parallel y \parallel∥x+y∥≤∥x∥+∥y∥
则称 ∥⋅∥\parallel \cdot \parallel∥⋅∥ 为 XXX 上的一个范数,对应的空间称线性赋范空间
- 设 XXX 是线性赋范空间,x, y∈Xx,\ y \in Xx, y∈X,称 ∥x−y∥\parallel x - y \parallel∥x−y∥ 为 xxx 和 yyy 之间的距离
- 📖设 XXX 是线性赋范空间,M⊆XM \subseteq XM⊆X 是 XXX 的子空间,f∈Xf \in Xf∈X。若 ∃φ∈M\exists \varphi \in M∃φ∈M 使 ∀ψ∈M\forall \psi \in M∀ψ∈M 有
∥f−φ∥≤∥f−ψ∥ \parallel f - \varphi \parallel \le \parallel f - \psi \parallel ∥f−φ∥≤∥f−ψ∥
则称 φ\varphiφ 是 fff 在 MMM 中的最佳逼近元
最佳一致逼近多项式
记 Mn={pn∣pnM_n = \{p_n | p_nMn={pn∣pn 为次数不超过 nnn 的多项式 }\}},则 Mn⊂C[a,b]M_n \subset C[a, b]Mn⊂C[a,b]
- 设 f∈C[a,b]f \in C[a, b]f∈C[a,b]。若 ∃pn∈Mn\exists p_n \in M_n∃pn∈Mn,使得对 ∀qn∈Mn\forall q_n \in M_n∀qn∈Mn,有 ∥f−pn∥∞≤∥f−qn∥∞\parallel f − p_n \parallel_{\infty} \le \parallel f − q_n\parallel_{\infty}∥f−pn∥∞≤∥f−qn∥∞。则称 pn(x)p_n(x)pn(x) 是 f(x)f(x)f(x) 的 nnn 次最佳一致逼近多项式
- 设 f∈C[a,b]f \in C[a, b]f∈C[a,b],则 fff 在 MnM_nMn 中存在唯一的 nnn 次最佳一致逼近多项式 pn(x)p_n(x)pn(x)
- 设 g∈C[a,b]g \in C[a, b]g∈C[a,b],如果 ∃x0∈[a,b]\exists x_0 \in [a, b]∃x0∈[a,b] 使得 ∣g(x0)∣=∥g∥∞=maxa≤x≤b∣g(x)∣|g(x_0)| = \parallel g\parallel_{\infty} = \max_{a \le x \le b} |g(x)|∣g(x0)∣=∥g∥∞=maxa≤x≤b∣g(x)∣,则称 x0x_0x0 为 g(x)g(x)g(x) 在 [a,b][a, b][a,b] 上的偏差点。当 g(x0)=∥g∥∞g(x_0) = \parallel g \parallel_{\infty}g(x0)=∥g∥∞,x0 称 g(x) 的正偏差点。当 g(x0)=−∥g∥∞g(x_0) = −\parallel g \parallel_{\infty}g(x0)=−∥g∥∞,x0x_0x0 称 g(x)g(x)g(x) 的负偏差点
- 设 f∈C[a,b]f \in C[a, b]f∈C[a,b],pn(x)p_n(x)pn(x) 是 f(x)f(x)f(x) 的 nnn 次最佳一致逼近多项式,则 f−pnf − p_nf−pn 必存在正负偏差点
- 设 f∈C[a,b]f \in C[a, b]f∈C[a,b],pn(x)p_n(x)pn(x) 是 nnn 次多项式,则 pn(x)p_n(x)pn(x) 是 f(x)f(x)f(x) 的 nnn 次最佳一致逼近多项式 ⇔\Leftrightarrow⇔ f(x)−pn(x)f(x) − p_n(x)f(x)−pn(x) 在 [a,b][a, b][a,b] 上至少有 (n+2)(n + 2)(n+2) 个交错偏差点,即存在 (n+2)(n + 2)(n+2) 个点 a≤x0<x1<⋯<xn<xn+1≤ba \le x_0 \lt x_1 \lt \cdots \lt x_n \lt x_{n+1} \le ba≤x0<x1<⋯<xn<xn+1≤b,使得
f(xi)−pn(xi)=(−1)iδ∥f−pn∥∞, i=0,1,...,n+1 f(x_i) - p_n(x_i) = (-1)^i\delta\parallel f - p_n \parallel_{\infty}, \ i = 0,1,..., n+1 f(xi)−pn(xi)=(−1)iδ∥f−pn∥∞, i=0,1,...,n+1
其中 σ=1\sigma = 1σ=1 或 σ=−1\sigma = -1σ=−1
- 📖设 f∈C[a,b]f \in C[a, b]f∈C[a,b],pn(x)p_n(x)pn(x) 是 f(x)f(x)f(x) 的 nnn 次最佳一致逼近多项式。如果 f(n+1)(x)f ^{(n+1)}(x)f(n+1)(x) 在 (a,b)(a, b)(a,b) 内存在且保号,则 f(x)−pn(x)f(x) − p_n(x)f(x)−pn(x) 在 [a,b][a, b][a,b] 内恰有 (n+2)(n + 2)(n+2) 个交错偏差点,且两端点 aaa,bbb 也是偏差点
- 📖设 f∈C[a,b]f \in C[a, b]f∈C[a,b],则 f(x)f(x)f(x) 的 nnn 次最佳一致逼近多项式 pn(x)p_n(x)pn(x) 为 f(x)f(x)f(x) 的某个 nnn 次插值多项式
本文详细介绍了多项式插值与函数的最佳逼近,包括Lagrange插值多项式、余项表示,以及Hermite插值的原理、应用实例和误差分析。讨论了高阶插值的缺点,如Runge现象,并提出了分段低次插值的解决方案,如线性插值和分段三次Hermite插值。
2242

被折叠的 条评论
为什么被折叠?



