logo

중심 B-스플라인 📂푸리에해석

중심 B-스플라인

정의1

mNm\in \mathbb{N}에 대해서, 중심 B-스플라인 BmB_{m}을 다음과 같이 정의한다.

Bm(x):=Tm2Nm(x)=Nm(x+12) B_{m}(x):= T_{-\frac{m}{2}}N_{m}(x)=N_{m}(x+{\textstyle \frac{1}{2}})

이때 TTL2L^{2} 공간트랜슬레이션이다.

설명

다음과 같이 정의할 수도 있다.

B1:=χ[1/2,1/2],Bm+1:=BmB1, mN B_{1}:= \chi_{[-1/2,1/2]},\quad B_{m+1}:=B_{m}*B_{1},\ m\in\mathbb{N}

두 정의는 실제로 같은 함수를 의미한다. 여기서 핵심은 BmB_{m}을 우함수가 되게끔 정의했다는 것이다. B-스플라인에서와 같이 다음의 식이 성립함을 쉽게 알 수 있다.

Bm+1(x)=Bm(xt)B1(t)dt=1212Bm(xt)dt B_{m+1}(x)=\int _{-\infty} ^{\infty} B_{m}(x-t)B_{1}(t)dt=\int_{-{\textstyle \frac{1}{2}}}^{{\textstyle \frac{1}{2}}}B_{m}(x-t)dt

중심 B-스플라인은 B-스플라인을 평행이동한 것에 불과하므로 B-스플라인의 성질을 그대로 갖는다.

성질

(a) suppBm=[m2,m2]\mathrm{supp} B_{m}=[-\frac{m}{2},\frac{m}{2}]

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

(c) m2m\ge 2에 대해서,

kZBm(xk)=1,xR \sum \limits_{k\in \mathbb{Z}}B_{m}(x-k)=1,\quad \forall x\in R

(c’) m=1m=1일 때, 위 식은 xR{±12,±32,}x\in \mathbb{R}\setminus \left\{ \pm\frac{1}{2},\pm \frac{3}{2},\dots \right\}에 대해서 성립한다.

(d) 중심 B-스플라인의 푸리에 변환은 다음과 같다.

Bm^(γ)=(eπiγeπiγ2πiγ)m=(sin(πγ)πγ)m \widehat{B_{m}}(\gamma)=\left( \frac{e^{\pi i \gamma}-e^{-\pi i \gamma}}{2\pi i \gamma} \right)^{m}=\left( \frac{\sin (\pi\gamma)}{\pi \gamma} \right)^{m}

이때 ff의 푸리에 변환 f^\widehat{f}의 정의는 다음과 같다.

f^(γ):=f(x)e2πixγdx \widehat{f}(\gamma):=\int _{-\infty} ^{\infty} f(x)e^{-2\pi i x\gamma}dx

증명

(a)

suppNm=[m2,m2]\mathrm{supp}N_{m}=\left[ -\frac{m}{2}, \frac{m}{2} \right]이고 Bm=Tm2NmB_{m}=T_{- \frac{m}{2}}N_{m}이므로,

suppBm=[0m2,mm2]=[m2,m2] \mathrm{supp} B_{m} = \left[0-\frac{m}{2},m-\frac{m}{2} \right] = \left[ -\frac{m}{2},\frac{m}{2} \right]

(b)

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

이고 Bm=Tm2NmB_{m}=T_{- \frac{m}{2}}N_{m}이므로,

Bm(x)dx=Tm2Bm(x)dx=1 \int _{-\infty} ^{\infty} B_{m}(x)dx=\int _{-\infty} ^{\infty} T_{-\frac{m}{2}}B_{m}(x)dx=1

(c)

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

이므로,

kZBm(xk)=kZTm2Nm(xk)=1,xR \sum \limits_{k\in \mathbb{Z}}B_{m}(x-k)=\sum \limits_{k \in \mathbb{Z}} T_{-\frac{m}{2}}N_{m}(x-k)=1,\quad \forall x\in \mathbb{R}

(c')

m=1m=1일 때, NmN_{m}에 대해서 xRZx\in \mathbb{R}\setminus \mathbb{Z}일 때 성립하고 Bm=Tm2NmB_{m}=T_{- \frac{m}{2}}N_{m}이므로, xR{±12,±32,}x\in \mathbb{R}\setminus \left\{ \pm\frac{1}{2},\pm \frac{3}{2},\dots \right\}에 대해서 성립한다.

(d)

B-스플라인의 푸리에 변환은 다음과 같다.

Nm^(γ)=(1e2πiγ2πiγ)m \widehat{N_{m}}(\gamma)=\left( \frac{1-e^{-2\pi i\gamma}}{2\pi i \gamma} \right)^{m}

그러면 푸리에 변환의 성질에 의해

F[Nm(x+m2)](γ)=e2πim2γNm^(γ) \mathcal{F}\left[ N_{m}(x+\frac{m}{2}) \right] (\gamma)=e^{2\pi i \frac{ m}{2} \gamma}\widehat{N_{m}}(\gamma)

이므로,

Bm^(γ)= F[Bm(x)](γ)=F[Nm(x+m2)](γ)= e2πim2γNm^(γ)= (eπiγ)m(1e2πiγ2πiγ)m= (eπiγeπiγ2πiγ)m= (sin(πγ)πγ)m \begin{align*} \widehat{B_{m}}(\gamma) =&\ \mathcal{F}\left[ B_{m}(x) \right] (\gamma) =\mathcal{F}\left[ N_{m}(x+\textstyle \frac{m}{2}) \right] (\gamma) \\ =&\ e^{2\pi i \frac{m}{2} \gamma}\widehat{N_{m}}(\gamma) \\ =&\ \left(e^{\pi i \gamma}\right)^{m}\left( \frac{1-e^{-2\pi i\gamma}}{2\pi i \gamma} \right)^{m} \\ =&\ \left( \frac{e^{\pi i \gamma}-e^{-\pi i \gamma}}{2\pi i \gamma} \right)^{m} \\ =&\ \left( \frac{\sin (\pi\gamma)}{\pi \gamma} \right)^{m} \end{align*}


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