logo

샘플링 정리 📂푸리에해석

샘플링 정리

빌드업1

어떤 물리적 신호 ff가 시간 t1<t2<t3<t_{1} < t_{2} < t_{3} < \cdots에 따라서 측정되고 있다고 생각해보자. f(t1),f(t2),f(t_{1}), f(t_{2}), \dots를 알고 있다고 해도, 일반적으로는 임의의 f(t)f(t) 값을 알 수는 없다. 하지만 여기서 신호 ff가 특정한 범위 내의 주파수의 신호만 포함한다고 가정해보자. 즉 어떤 상수 Ω\Omega보다 작은 주파수만 포함하는 신호 ff를 생각하자는 것인데, 이러한 ff대역제한 신호band-limited signal, 주파수제한 신호라고 한다.

이를 푸리에 해석의 언어로 말해보면 f^(ω)\hat{f}(\omega)ωΩ|\omega| \ge \Omega인 영역에서는 함숫값이 전부 00이라는 것과 같다. 따라서 ff가 대역제한이라는 조건은 f^(ω)=0 for ωΩ\hat{f} (\omega) = 0\ \text{for } | \omega | \ge \Omega라는 조건과 동치이고, f^L1\hat{f} \in L^{1}라는 것을 의미한다. 이러한 조건에서는 아래와 같은 강력한 정리가 성립한다.

정리

fL2f\in L^{2}이고 f^(ω)=0 for ωΩ\hat{f} (\omega) = 0\ \text{for } | \omega | \ge \Omega라고 가정하자. 그러면 f(t)f(t)nπ/Ω(n=0,±1,±2,)n\pi / \Omega (n=0, \pm 1, \pm 2, \dots)에서의 값들로 결정된다. 즉 다음이 성립한다.

f(t)=n=f(nπΩ)sin(Ωtnπ)Ωtnπ f(t) = \sum \limits_{n = -\infty}^{\infty} f \left( \dfrac{n\pi}{\Omega}\right) \dfrac{\sin (\Omega t -n\pi )}{\Omega t - n\pi}

설명

이를 샘플링 정리sampling theorem라 한다. 샘플링 정리는 가산개의 함숫값 f(nπΩ)f \left( \dfrac{n\pi}{\Omega}\right)만으로도 모든 tt에 대한 f(t)f(t)를 결정할 수 있는 조건을 제시한다.

신호 해석의 관점에서 샘플링 정리를 봤을 때 f^(ω)=0 for ωΩ\hat{f} (\omega) = 0\ \text{for } | \omega | \ge \Omega라는 조건은 신호 ff의 주파수가 제한되어 있다는 말과 같다.

한편 주파수 함수 f^L2\hat{f} \in L^{2}시간제한 신호time-limited signal f(t)f(t)에 대해서도 위와 같은 꼴의 공식이 성립한다.

주파수 샘플링 정리

f^L2\hat{f} \in L^{2}이고 f(t)=0 for tLf (t) = 0\ \text{for } | t | \ge L라고 가정하자. 그러면 f^(ω)\hat{f}(\omega)nπ/L(n=0,±1,±2,)n\pi / L(n=0, \pm 1, \pm 2, \dots)에서의 값들로 결정된다. 즉 다음이 성립한다.

f^(ω)=n=f^(nπL)sin(Lωnπ)Lωnπ \hat{f}(\omega) = \sum \limits_{n = -\infty}^{\infty} \hat{f} \left( \dfrac{n\pi}{L}\right) \dfrac{\sin (L \omega -n\pi )}{L \omega - n\pi}

증명은 시간 샘플링 정리와 같다.

증명

f^L1\hat{f} \in L^{1}이므로 f^\hat{f}푸리에 급수로 표현가능하다. f^(ω)\hat{f} (\omega)[Ω,Ω][-\Omega, \Omega]에서 정의되어 있으므로 복소 푸리에 급수는 다음과 같다.

f^(ω)=cneinπω/Ω(ωΩ) \begin{equation} \hat{f}(\omega)=\sum_{-\infty}^{\infty} c_{-n} e^{-i n \pi \omega / \Omega} \quad(|\omega| \leq \Omega) \label{eq1} \end{equation}

이때 증명 과정에서의 편의를 위해 nnn-n으로 표기하였다. 계수 cnc_{-n}은 다음과 같다.

cn=12ΩΩΩf^(ω)einπω/Ωdω=12Ωf^(ω)einπω/Ωdω=πΩ12πf^(ω)eiω(nπ/Ω)dω \begin{align*} c_{-n} &= \dfrac{1}{2 \Omega} \int_{-\Omega}^{\Omega} \hat{f}(\omega) e^{i n \pi \omega / \Omega} d \omega \\ &= \dfrac{1}{2 \Omega} \int_{-\infty}^{\infty} \hat{f}(\omega) e^{i n \pi \omega / \Omega} d \omega \\ &= \dfrac{\pi}{\Omega} \dfrac{1}{2\pi} \int_{-\infty}^{\infty} \hat{f}(\omega) e^{i \omega (n \pi / \Omega)} d\omega \end{align*}

이때 f^(ω)=0 for ωΩ\hat{f} (\omega) = 0\ for\ | \omega | \ge \Omega이므로 f^(ω)L2\hat{f} (\omega) \in L^{2}이고, 푸리에 역변환 정리에 의해 다음이 성립한다.

f(t)=12πf^(ω)eiωtdω \begin{equation} f(t) =\dfrac{1}{2\pi} \int _{-\infty} ^{\infty}\hat{f}(\omega) e^{i\omega t}d\omega \label{eq2} \end{equation}

따라서 다음의 식을 얻는다.

cn=πΩ12πf^(ω)eiω(nπ/Ω)dω=πΩf(nπΩ) \begin{align*} c_{-n} &= \dfrac{\pi}{\Omega} \dfrac{1}{2\pi} \int_{-\infty}^{\infty} \hat{f}(\omega) e^{i \omega (n \pi / \Omega)} d\omega \\ &= \dfrac{\pi}{\Omega} f \left( \dfrac{n\pi}{\Omega}\right) \end{align*}

이를 (eq1)\eqref{eq1}에 대입하면 다음을 얻는다.

f^(ω)=cneinπω/Ω(ωΩ)=πΩf(nπΩ)einπω/Ω(ωΩ) \hat{f}(\omega) = \sum_{-\infty}^{\infty} c_{-n} e^{-i n \pi \omega / \Omega} \quad(|\omega| \leq \Omega) = \sum_{-\infty}^{\infty} \dfrac{\pi}{\Omega} f \left( \dfrac{n\pi}{\Omega}\right) e^{-i n \pi \omega / \Omega} \quad(|\omega| \leq \Omega)

다시 이를 (eq2)\eqref{eq2}에 대입하면 아래의 식을 얻는다.

f(t)=12πf^(ω)eiωtdω=12πΩΩf^(ω)eiωtdω=12ΩΩΩf(nπΩ)einπω/Ωeiωtdω \begin{align*} f(t) &= \frac{1}{2 \pi} \int_{-\infty}^{\infty} \hat{f}(\omega) e^{i \omega t} d \omega = \frac{1}{2 \pi} \int_{-\Omega}^{\Omega} \hat{f}(\omega) e^{i \omega t} d \omega \\ &=\frac{1}{2 \Omega} \int_{-\Omega}^{\Omega} \sum_{-\infty}^{\infty} f\left(\frac{n \pi}{\Omega}\right) e^{-i n \pi \omega / \Omega} e^{i \omega t} d \omega \end{align*}

함수공간에서는 내적이 정적분으로 정의되고, 내적은 연속 이므로 극한을 밖으로 빼낼 수 있으므로 위 식은 다음과 같다.

f(t)=12ΩΩΩf(nπΩ)einπω/Ωeiωtdω=12Ωf(nπΩ)einπω/Ω,eiωt=12Ωf(nπΩ)einπω/Ω,eiωt=12ΩΩΩf(nπΩ)einπω/Ωeiωtdω=12Ωf(nπΩ)ΩΩeinπω/Ωeiωtdω \begin{align*} f(t) &= \frac{1}{2 \Omega} \int_{-\Omega}^{\Omega} \sum_{-\infty}^{\infty} f\left(\frac{n \pi}{\Omega}\right) e^{-i n \pi \omega / \Omega} e^{i \omega t} d \omega \\ &= \frac{1}{2 \Omega} \left\langle \sum_{-\infty}^{\infty} f\left(\frac{n \pi}{\Omega}\right) e^{-i n \pi \omega / \Omega}, e^{-i \omega t} \right\rangle \\ &= \frac{1}{2 \Omega} \sum_{-\infty}^{\infty} \left\langle f\left(\frac{n \pi}{\Omega}\right) e^{-i n \pi \omega / \Omega}, e^{-i \omega t} \right\rangle \\ &= \frac{1}{2 \Omega} \sum_{-\infty}^{\infty} \int_{-\Omega}^{\Omega} f\left(\frac{n \pi}{\Omega}\right) e^{-i n \pi \omega / \Omega} e^{i \omega t} d \omega \\ &= \frac{1}{2 \Omega} \sum_{-\infty}^{\infty} f\left(\frac{n \pi}{\Omega}\right) \int_{-\Omega}^{\Omega} e^{-i n \pi \omega / \Omega} e^{i \omega t} d \omega \end{align*}

위의 적분을 계산해보면 다음과 같다.

ΩΩeinπω/Ωeiωtdω=ΩΩei(Ωtnπ)ω/Ωdω=ei(Ωtnπ)ω/Ωi(Ωtnπ)/ΩΩΩ=1i(Ωtnπ)/Ω[ei(Ωtnπ)ei(Ωtnπ)]=2(Ωtnπ)/Ωei(Ωtnπ)ei(Ωtnπ)2i=2sin(Ωtnπ)(Ωtnπ)/Ω \begin{align*} \int_{-\Omega}^{\Omega} e^{-i n \pi \omega / \Omega} e^{i \omega t} d \omega &= \int_{-\Omega}^{\Omega} e^{i(\Omega t- n \pi )\omega/\Omega} d \omega \\ &= \left. \frac{e^{i(\Omega t-n \pi) \omega / \Omega}}{i(\Omega t-n \pi) / \Omega}\right|_{-\Omega} ^{\Omega} \\ &= \frac{1}{i(\Omega t-n \pi) / \Omega}\left[ e^{i(\Omega t-n \pi)} - e^{-i(\Omega t-n \pi)} \right] \\ &= \frac{2}{(\Omega t-n \pi) / \Omega}\dfrac{ e^{i(\Omega t-n \pi)} - e^{-i(\Omega t-n \pi)}}{2i} \\ &= \frac{2 \sin (\Omega t-n \pi)}{(\Omega t-n \pi) / \Omega} \end{align*}

대입해서 정리하면 아래와 같다.

f(t)=12Ωf(nπΩ)ΩΩeinπω/Ωeiωtdω=12Ωf(nπΩ)2sin(Ωtnπ)(Ωtnπ)/Ω=f(nπΩ)2sin(Ωtnπ)Ωtnπ \begin{align*} f(t) &= \frac{1}{2 \Omega} \sum_{-\infty}^{\infty} f\left(\frac{n \pi}{\Omega}\right) \int_{-\Omega}^{\Omega} e^{-i n \pi \omega / \Omega} e^{i \omega t} d \omega \\ &= \frac{1}{2 \Omega} \sum_{-\infty}^{\infty} f\left(\frac{n \pi}{\Omega}\right) \frac{2 \sin (\Omega t-n \pi)}{(\Omega t-n \pi) / \Omega} \\ &= \sum_{-\infty}^{\infty} f\left(\frac{n \pi}{\Omega}\right) \frac{2 \sin (\Omega t-n \pi)}{\Omega t-n \pi} \end{align*}


  1. Gerald B. Folland, Fourier Analysis and Its Applications (1992), p230-231 ↩︎