logo

해석학에서 스플라인, B-스플라인 📂푸리에해석

해석학에서 스플라인, B-스플라인

정의1

함수 f:RRf:\mathbb{R} \to \mathbb{R}구간마다 다항식piecewise polynomial이면 ffR\mathbb{R}위에서의 스플라인spline이라고 한다. 다항식이 바뀌는 점을 knot 이라 한다.

설명

정의를 보면 알 수 있듯이 스플라인이 연속함수여야하는 것은 아니다. 아래와 같은 함수 ff는 스플라인의 한 예이다.

f(x)={0x[,0]2x2x(0,1]2xx(1,4]116x3x(4,] f(x) = \begin{cases} 0 & x\in[\infty,0] \\ 2x^{2}&x\in(0,1] \\ 2-x & x\in (1,4] \\ \frac{1}{16}x^{3} & x\in(4,\infty] \end{cases}

untitled2.png

위의 경우에서 x=0x=0, x=1x=1, x=4x=4가 놋이다. B-스플라인 은 스플라인 중에서 좋은 특징을 가진 스플라인이다. BB-스플라인 N1N_{1}을 아래와 같이 구간 [0,1][0,1]에서의 특성함수로 정의한다.

N1(x):=χ[0,1](x),xR N_{1}(x) :=\chi_{[0,1]}(x)\quad , x\in \mathbb{R}

그리고 mNm \in \mathbb{N}에 대해서 B-스플라인 Nm+1N_{m+1}을 아래와 같이 정의한다.

Nm+1(x):=(NmN1)(x) \begin{equation} N_{m+1}(x) := (N_{m} * N_{1})(x)\end{equation}

이때 \ast컨볼루션이다. mm을 B-스플라인 NmN_{m}오더order라고 한다. 정의 (1)(1)에 의해 다음이 성립한다.

Nm=Nm1N1=Nm2N1N1=Nm3N1N1N1=N1N1N1N1m \begin{align*} N_{m} &= N_{m-1} \ast N_{1} \\ &= N_{m-2} \ast N_{1} \ast N_{1} \\ &= N_{m-3} \ast N_{1} \ast N_{1}\ast N_{1} \\ &= \underbrace{N_{1} \ast N_{1} \ast N_{1} \cdots \ast N_{1}}_{m} \end{align*}

또한 N1N_{1}컨볼루션의 정의로 아래의 식이 성립함을 알 수 있다.

Nm+1(x)=Nm(xt)N1(t)dt=01Nm(xt)dt N_{m+1}(x)=\int _{-\infty} ^{\infty}N_{m}(x-t)N_{1}(t)dt=\int_{0}^{1}N_{m}(x-t)dt

아래의 그림은 왼쪽에서부터 N2N_{2}, N3N_{3}의 그래프를 그린 것이다.2020-10-05135620.png

성질

오더가 mNm\in \mathbb{N}B-스플라인은 다음과 같은 성질을 만족한다.

(a) suppNm=[0,m]\mathrm{supp}N_{m}=[0,m] and\text{and} Nm(x)>0 for x(0,m)N_{m}(x)>0 \text{ for } x\in(0,m)

(b) Nm(x)dx=1\displaystyle \int _{-\infty} ^{\infty} N_{m}(x)dx=1

(c) m2m\ge 2에 대해서 아래의 식이 성립한다.

kZNm(xk)=1,xR \begin{equation} \sum \limits_{k \in \mathbb{Z}} N_{m}(x-k)=1,\quad \forall x\in \mathbb{R} \end{equation}

(c’) m=1m=1일 때, 위 식은 xRZx\in \mathbb{R}\setminus \mathbb{Z}에 대해서 성립한다.

같이보기


  1. Ole Christensen, Functions, Spaces, and Expansions: Mathematical Tools in Physics and Engineering (2010), p203-204 ↩︎