数值分析复习(四)——多项式插值与函数逼近

本文详细介绍了多项式插值与函数的最佳逼近,包括Lagrange插值多项式、余项表示,以及Hermite插值的原理、应用实例和误差分析。讨论了高阶插值的缺点,如Runge现象,并提出了分段低次插值的解决方案,如线性插值和分段三次Hermite插值。

四、多项式插值与函数最佳逼近

设函数 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 bax0<x1<<xnb 上的值 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=knxkxixxi
  • 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)\ (\xix[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=0n(xxi)=(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)=sin⁡x,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=0x1=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)(02π)(0π)(x2π)(xπ)+f(2π)(2π0)(2ππ)(x0)(xπ)+f(π)(π0)(π2π)(x0)(x2π)=π24x(xπ)

例2:设函数
f(x)=2π∫0xe−t2dt f(x)=\frac{2}{\sqrt{\pi}}\int_0^xe^{-t^2}dt f(x)=π20xet2dt
的函数值已造成函数表。假设在区间 [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(ξ)(xx0)(xx1), ξ[4,6]=π2ex2=π4xex2=π4(2x21)ex2>0, x(4,6)21×f(4)×(54)(56)=0.508×106

差商及 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} 0f[xi]=f(xi)1f[xi,xj]=xjxif(xj)f(xi)2f[xi,xj,xk]=xkxif[xj,xk]f[xi,xj]kf[x0,x1,...,xk1,xk]=xkx0f[x1,...,xk]f[x0,...,xk1]

  • 😄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=0ki=0,i=mk(xmxi)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](xx0)+f[x0,x1,x2](xx0)(xx1)+...+f[x0,x1,...,xn](xx0)(xx1)...(xxn1)

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=0n(xxi)mi+1

  • f∈Cn[a,b]f \in C^n[a, b]fCn[a,b]a≤x0<x1<...<xn≤ba \le x_0 \lt x_1 \lt ... \lt x_n \le bax0<x1<...<xnb,则有

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=1i=1ntiτn={(t1,t2,...,tn)∣ti≥0\tau_n=\{(t_1, t_2, ..., t_n)|t_i \ge 0τn={(t1,t2,...,tn)ti0∑i=1nti≤1}\sum_{i=1}^nt_i \le 1\}i=1nti1}nnn 维单纯形

  • 注意到被积函数是通过一元连续函数 f(n)(x)f^{(n)}(x)f(n)(x)nnn 元线性连续函数 ∑i=0ntixi\sum_{i=0}^nt_ix_ii=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]=lim⁡x1,x2,...,xk→x0f[x0,x1,...,xk]=lim⁡x1,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,...,xkx0limf[x0,x1,...,xk]=x1,x2,...,xkx0limk!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](xx0)+...+f[x0,...,x0](xx0)m0+f[x0,...,x0,x1](xx0)m0+1+...+f[x0,...,x0,x1,x1](xx0)m0+1(xx1)+f[x0,...,x0,x1,...x1,x2](xx0)m0+1(xx1)(m1+1)+...+f[x0,...,x0,...,xn1,...,xn1,xn,...,xn](xx0)m0+1...(xxn1)mn1+1(xxn)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=0n(xxi)mi+1

例题

例1:若 n=0n = 0n=0m0=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)(xx0)+...+k!f(k)(x0)(xx0)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)=3H(0)=4H(1)=5H(1)=6H(1)=7

:列表计算各点差商

kkkxkx_kxkf(xk)f(x_k)f(xk)1阶差商2阶差商3阶差商4阶差商
0034-26−132-\frac{13}{2}213
10324−12-\frac{1}{2}21
215672\frac{7}{2}27
3156
415

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(x0)2(x0)2+6(x0)2(x1)213(x0)2(x1)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)(ξ)(x0)2(x1)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 上的函数值:

xxxx0x_0x0x1x_1x1⋯\cdotsxn−1x_{n-1}xn1xnx_nxn
f(x)f(x)f(x)f(x0)f(x_0)f(x0)f(x1)f(x_1)f(x1)⋯\cdotsf(xn−1)f(x_{n-1})f(xn1)f(xn)f(x_n)f(xn)

📖记 hi=xi+1−xih_i = x_{i+1} − x_ihi=xi+1xih=max⁡0≤i≤n−1hih = \max_{0 \le i \le n−1}h_ih=max0in1hi,在每个小区间 [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](xxi),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)(xxi)(xxi+1),ξi(xi,xi+1)
从而有
max⁡xi≤x≤xi+1∣f(x)−L1,i(x)∣≤max⁡xi≤x≤xi+1∣12f′′(ξi)(x−xi)(x−xi+1)∣≤18hi2max⁡xi≤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)| xixxi+1maxf(x)L1,i(x)xixxi+1max21f(ξi)(xxi)(xxi+1)81hi2xixxi+1maxf(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,n2(x),L1,n1(x),x[x0,x1)x[x1,x2)x[xn2,xn1)x[xn1,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)分段线性插值函数,其插值误差为
max⁡a≤x≤b∣f(x)−L~1(x)∣=max⁡0≤i≤n−1max⁡xi≤x≤xi+1∣f(x)−L~1(x)∣=max⁡0≤i≤n−1max⁡xi≤x≤xi+1∣f(x)−L1,i(x)∣≤max⁡0≤i≤n−118hi2max⁡xi≤x≤xi+1∣f′′(x)∣≤18h2max⁡a≤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} axbmaxf(x)L~1(x)=0in1maxxixxi+1maxf(x)L~1(x)=0in1maxxixxi+1maxf(x)L1,i(x)0in1max81hi2xixxi+1maxf(x)81h2axbmaxf(x)

  • 只要 f(x)f(x)f(x)[a,b][a, b][a,b] 上有 2 阶连续导数,当 h→0h \rightarrow 0h0 时,余项一致趋于零,即分段线性插值函 数 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 上的函数表

xxxx0x_0x0x1x_1x1⋯\cdotsxn−1x_{n-1}xn1xnx_nxn
f(x)f(x)f(x)f(x0)f(x_0)f(x0)f(x1)f(x_1)f(x1)⋯\cdotsf(xn−1)f(x_{n-1})f(xn1)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)⋯\cdotsf′(xn−1)f^{'}(x_{n-1})f(xn1)f′(xn)f^{'}(x_n)f(xn)

📖记 hi=xi+1−xih_i = x_{i+1} − x_ihi=xi+1xih=max⁡0≤i≤n−1hih = \max_{0 \le i \le n−1}h_ih=max0in1hi,在每个小区间 [xi,xi+1][x_i , x_{i+1}][xi,xi+1] 上利用数据

xxxxix_ixixi+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)(xxi)+hif[xi,xi+1]f(xi)(xxi)2+hi2f(xi+1)2f[xi,xi+1]+f(xi)(xxi)2(xxi+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)(ξ)(xxi)2(xxi+1)2, ξ(xi,xi+1)
于是
max⁡xi≤x≤xi+1∣f(x)−H3,i(x)∣≤14!hi416max⁡xi≤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)| xixxi+1maxf(x)H3,i(x)4!116hi4xixxi+1maxf(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,n2(x)H3,n1(x)x[x0,x1)x[x1,x2)x[xn2,xn1)x[xn1,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)

其误差为
max⁡a≤x≤b∣f(x)−H~3(x)∣=max⁡0≤i≤n−1max⁡xi≤x≤xi+1∣f(x)−H~3(x)∣=max⁡0≤i≤n−1max⁡xi≤x≤xi+1∣f(x)−H3,i(x)∣≤max⁡0≤i≤n−114!hi416max⁡xi≤x≤xi+1∣f(4)(x)∣≤1384h4max⁡a≤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} axbmaxf(x)H~3(x)=0in1maxxixxi+1maxf(x)H~3(x)=0in1maxxixxi+1maxf(x)H3,i(x)0in1max4!116hi4xixxi+1maxf(4)(x)3841h4axbmaxf(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)=0S(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(xn0)S(x0+0)=S(xn0)
样条插值的建立

TODO

最佳一致逼近

线性赋范空间
  • XXX 是一个集合,如果对 ∀x, y∈X, λ∈R\forall x,\ y \in X, \ \lambda \in Rx, yX, λR, 有 λx∈X, x+y∈X\lambda x \in X,\ x + y \in XλxX, x+yX,则称 XXX 是线性空间
  • XXXRRR 上的一个线性空间,φ1, φ2,⋯ , φm∈X\varphi_1,\ \varphi_2, \cdots ,\ \varphi_m \in Xφ1, φ2,, φmX,如果存在不全为零的数 a1, a2,⋯ , am∈Ra_1,\ a_2, \cdots ,\ a_m \in Ra1, a2,, amR 使得

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 XxX,对应于实数,记为 ∥x∥\parallel x \parallelx 且满足下面关系:
    • ∀x∈X\forall x \in XxX,有 ∥x∥≥0\parallel x \parallel \ge 0x0,且 ∥x∥=0⇒x=0\parallel x \parallel = 0 \Rightarrow x=0x=0x=0
    • ∀λ∈R\forall \lambda \in RλR x∈X\ x \in X xX,有 ∥λx∥=∣λ∣∥x∥\parallel \lambda x \parallel = |\lambda|\parallel x \parallelλx=λx
    • x,y∈Xx, y \in Xx,yX,有 ∥x+y∥≤∥x∥+∥y∥\parallel x+y \parallel \le \parallel x \parallel + \parallel y \parallelx+yx+y

则称 ∥⋅∥\parallel \cdot \parallelXXX 上的一个范数,对应的空间称线性赋范空间

  • XXX 是线性赋范空间,x, y∈Xx,\ y \in Xx, yX,称 ∥x−y∥\parallel x - y \parallelxyxxxyyy 之间的距离
  • 📖设 XXX 是线性赋范空间,M⊆XM \subseteq XMXXXX子空间f∈Xf \in XfX。若 ∃φ∈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φfffMMM 中的最佳逼近元

最佳一致逼近多项式

Mn={pn∣pnM_n = \{p_n | p_nMn={pnpn 为次数不超过 nnn 的多项式 }\}},则 Mn⊂C[a,b]M_n \subset C[a, b]MnC[a,b]

  • f∈C[a,b]f \in C[a, b]fC[a,b]。若 ∃pn∈Mn\exists p_n \in M_npnMn,使得对 ∀qn∈Mn\forall q_n \in M_nqnMn,有 ∥f−pn∥∞≤∥f−qn∥∞\parallel f − p_n \parallel_{\infty} \le \parallel f − q_n\parallel_{\infty}fpnfqn。则称 pn(x)p_n(x)pn(x)f(x)f(x)f(x)nnn最佳一致逼近多项式
  • f∈C[a,b]f \in C[a, b]fC[a,b],则 fffMnM_nMn 中存在唯一nnn 次最佳一致逼近多项式 pn(x)p_n(x)pn(x)
  • g∈C[a,b]g \in C[a, b]gC[a,b],如果 ∃x0∈[a,b]\exists x_0 \in [a, b]x0[a,b] 使得 ∣g(x0)∣=∥g∥∞=max⁡a≤x≤b∣g(x)∣|g(x_0)| = \parallel g\parallel_{\infty} = \max_{a \le x \le b} |g(x)|g(x0)=g=maxaxbg(x),则称 x0x_0x0g(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)=gx0x_0x0g(x)g(x)g(x)负偏差点
  • f∈C[a,b]f \in C[a, b]fC[a,b]pn(x)p_n(x)pn(x)f(x)f(x)f(x)nnn 次最佳一致逼近多项式,则 f−pnf − p_nfpn 必存在正负偏差点
  • f∈C[a,b]f \in C[a, b]fC[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 bax0<x1<<xn<xn+1b,使得

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δfpn, i=0,1,...,n+1

其中 σ=1\sigma = 1σ=1σ=−1\sigma = -1σ=1

  • 📖设 f∈C[a,b]f \in C[a, b]fC[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) 个交错偏差点,且两端点 aaabbb 也是偏差点
  • 📖设 f∈C[a,b]f \in C[a, b]fC[a,b],则 f(x)f(x)f(x)nnn 次最佳一致逼近多项式 pn(x)p_n(x)pn(x)f(x)f(x)f(x) 的某个 nnn 次插值多项式
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值