三次样条插值

最新推荐文章于 2020-12-17 08:58:34 发布

怡乐长吟

最新推荐文章于 2020-12-17 08:58:34 发布

阅读量997

点赞数

分类专栏：逼近论文章标签：逼近论数值逼近样条插值

本文链接：https://blog.csdn.net/RJ_MC_2016/article/details/78088891

版权

逼近论专栏收录该内容

5 篇文章 0 订阅

订阅专栏

$\newcommand{\nj}[1]{{\left\langle#1\right\rangle}}\newcommand{\kh}[1]{{\left(#1\right)}}\newcommand{\zkh}[1]{{\left[#1\right]}}\newcommand{\dkh}[1]{{\left\{{#1}\right \}}}\newcommand{\norm}[1]{\left|#1\right|}\newcommand{\normm}[1]{\left\Vert#1\right\Vert}\newcommand{\calA}{{\mathcal{A}}}\newcommand{\calC}{{\mathcal{C}}}\newcommand{\calH}{{\mathcal{H}}}\newcommand{\calN}{{\mathcal{N}}}\newcommand{\calQ}{{\mathcal{Q}}}\newcommand{\calR}{{\mathcal{R}}}\newcommand{\calW}{{\mathcal{W}}}\newcommand{\calX}{{\mathcal{X}}}\newcommand{\calZ}{{\mathcal{Z}}}\newcommand{\bbC}{\mathbb{C}}\newcommand{\bbN}{\mathbb{N}}\newcommand{\bbR}{\mathbb{R}}\newcommand{\bbZ}{\mathbb{Z}}\newcommand{\sin}{{\rm sin}}\newcommand{\cos}{{\rm cos}}\newcommand{\Id}{{\rm Id}}$

局部插值方法

　　指插值函数满足 $s\kh{x}=\sum_{i=1}^nl_i\kh{x}f\kh{x_i}$ ，其中 $l_i\kh{x_j}=\delta_{ij}$ ，且 $l_i$ 仅在 $x_i$ 某邻域内非零。例如分片线性、局部分片三次多项式插值。

三次样条插值

　　给定插值点均匀分布，间隔为 $h$ 。要求 $s$ 为有界分片三次函数，且插值函数二次导数连续。则有 $s$ 每个分片 $s_i$

⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ s j (x) = s j (x j + 1) = s'' j (x j + 1) = f (x j) + s' (x j) (x - x j) + c 2 j (x - x j) 2 + c 3 j (x - x j) 3 f (x j + 1) s'' j + 1 (x j + 1)

$\left\{ \begin{align*} s_j\kh{x}={}&f\kh{x_j}+s'\kh{x_j}\kh{x-x_j}+c^2_j\kh{x-x_j}^2+c^3_j\kh{x-x_j}^3\\ s_j\kh{x_{j+1}}={}&f\kh{x_{j+1}}\\ s''_j\kh{x_{j+1}}={}&s''_{j+1}\kh{x_{j+1}} \end{align*} \right.$

得

⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ c 2 j = c 3 j = 3 f ( x j + 1 ) - f ( x j ) h 2 - 2 s ' ( x j ) - s ' ( x j + 1 ) h 2 f ( x j ) - f ( x j + 1 ) h 3 + s ' ( x j ) + s ' ( x j + 1 ) h 2

$\left\{ \begin{align*} c^2_j={}&3\frac{f\kh{x_{j+1}}-f\kh{x_j}}{h^2}-\frac{2s'\kh{x_{j}}-s'\kh{x_{j+1}}}{h}\\ c^3_j={}&2\frac{f\kh{x_{j}}-f\kh{x_{j+1}}}{h^3}+\frac{s'\kh{x_{j}}+s'\kh{x_{j+1}}}{h^2} \end{align*} \right.$

即

c 2 j + 1 = c 2 j + 3 c 3 j h ⟹ s' (x j - 1) + 4 s' (x j) + s' (x j + 1) = 3 h [f (x j + 1) - f (x j - 1)] .

$c_{j+1}^2=c_j^2+3c_j^3h\Longrightarrow s'\kh{x_{j-1}}+4s'\kh{x_{j}}+s'\kh{x_{j+1}}=\frac{3}{h}\zkh{f\kh{x_{j+1}}-f\kh{x_{j-1}}}.$

　　代入 $s\kh{x}=\sum_{i=1}^nl_i\kh{x}f\kh{x_i}$ ，其中 $l_i\kh{x_j}=\delta_{ij}$ 。则有

l' i (x j - 1) + 4 l' i (x j) + l' i (x j + 1) = 3 h [δ i, j + 1 - δ i, j - 1] .

$l'_i\kh{x_{j-1}}+4l'_i\kh{x_{j}}+l'_i\kh{x_{j+1}}=\frac{3}{h}\zkh{\delta_{i,j+1}-\delta_{i,j-1}}.$

　　由上式右端当 $j>i+1$ 时为0，易得对 $j>i$ ， $l'_i\kh{x_j}=\alpha\kh{-2+\sqrt{3}}^{j-i}+\beta\kh{-2-\sqrt{3}}^{j-i}$ 。 $l_i$ 有界，故 $\beta=0$ 。 $j<i$ 同理。从而代入得

l' i (x j) = ⎧ ⎩ ⎨ ⎪ ⎪ ⎪ ⎪ - 3 h (- 2 + 3 \sqrt) j - i 0 3 h (- 2 + 3 \sqrt) j - i j < i j = i j > i

$l'_i\kh{x_j}=\begin{cases} -\frac{3}{h}\kh{-2+\sqrt{3}}^{j-i}&j<i\\ 0&j=i\\ \frac{3}{h}\kh{-2+\sqrt{3}}^{j-i}&j>i \end{cases}$

结合 $l_i\kh{x_j}=\delta_{ij}$ 可解出 $l_i$ 方程。

= = : = = max x \in [0, h] ∥ f ∥ \infty \leq 1 ∣ ∣ ∣ ∣ \sum j \in N l j (x) f (x j) ∣ ∣ ∣ ∣ max x \in [0, h] \sum j \in N ∣ ∣ l j (x) ∣ ∣ max x \in [0, h] \sum j \in Z (- 1) j [l - j (x) + l j + 1 (x)] max x \in [0, h] p (x) 1 + 3 3 \sqrt 4 .

$\begin{align*} &\max_{\substack{x\in\zkh{0,h}\\\normm{f}_\infty\leq1}}\norm{\sum_{j\in\bbN}l_j\kh{x}f\kh{x_j}}\\ ={}&\max_{x\in\zkh{0,h}}\sum_{j\in\bbN}\norm{l_j\kh{x}}\\ ={}&\max_{x\in\zkh{0,h}}\sum_{j\in\bbZ}\kh{-1}^j\zkh{l_{-j}\kh{x}+l_{j+1}\kh{x}}\\ :={}&\max_{x\in\zkh{0,h}}p\kh{x}\\ ={}&\frac{1+3\sqrt{3}}{4}. \end{align*}$