5장. 푸리에 급수

푸리에 해석 · 푸리에 변환

심화 · 약 50분 · 연습문제 6 · 그림 5

요약

주기함수의 복소지수 전개, 계수 공식, 실수 신호의 세 가지 표기, 계수 계산 예제, 대칭성, 수렴과 깁스 현상, 파세발 등식

주기적인 신호(예: 일정한 음높이의 소리)는 기본 주파수와 그 정수배 주파수의 사인파를 더해 만들 수 있습니다. 이것이 푸리에 급수입니다. 이 장에서는 계수를 내적으로 구하는 방법, 불연속점 근처의 튀어나옴(깁스 현상), 급수가 언제 원래 함수로 수렴하는지를 다룹니다.

주기함수를 복소지수함수의 무한합으로 씁니다. 4장의 CN\mathbb{C}^N을 함수공간 L2[0,T]L^2[0, T]로, 합 ∑n\sum_n을 적분 ∫0Tdt\int_0^T dt로 바꾼 것입니다.

5.1 설정

주기 TT인 함수 x(t)x(t)를 생각합니다. ω0=2π/T\omega_0 = 2\pi/T를 기본 각주파수라 합니다. 주기함수는 한 주기만 알면 되므로 L2[0,T]L^2[0, T]의 벡터로 보고, 내적은 ⟨x,y⟩=∫Tx(t)y(t)∗ dt\langle x, y \rangle = \int_T x(t)y(t)^*\, dt입니다.

주기 TT를 갖는 복소지수함수는 ejωte^{j\omega t} 중 ejωT=1e^{j\omega T} = 1인 것, 즉 ω=kω0\omega = k\omega_0 (kk 정수)인 것들뿐입니다.

ek(t)=ejkω0t,k=0,±1,±2,…e_k(t) = e^{jk\omega_0 t}, \qquad k = 0, \pm 1, \pm 2, \ldots

k=±1k = \pm 1을 기본파, ∣k∣≥2|k| \ge 2를 고조파(harmonic)라 합니다. 2장 정리 2.9에 의해

⟨ek,em⟩=Tδkm\langle e_k, e_m \rangle = T\delta_{km}

입니다. {ek}\{e_k\}는 L2[0,T]L^2[0, T]의 직교계입니다. DFT와 다른 점은 kk가 유한 개가 아니라 모든 정수라는 것뿐입니다(연속 함수에서는 ek+N=eke_{k+N} = e_k 같은 중복이 없습니다).

5.2 계수 공식

정의 5.1 (푸리에 급수). 주기 TT인 함수 xx의 푸리에 계수(Fourier coefficient)와 푸리에 급수(Fourier series)는

ck=⟨x,ek⟩∥ek∥2=1T∫Tx(t) e−jkω0t dt,x(t)=∑k=−∞∞ckejkω0tc_k = \frac{\langle x, e_k \rangle}{\lVert e_k\rVert^2} = \frac{1}{T}\int_T x(t)\, e^{-jk\omega_0 t}\, dt, \qquad x(t) = \sum_{k=-\infty}^{\infty} c_k e^{jk\omega_0 t}

입니다.

계수 공식은 직교기저의 좌표 공식(3장 따름정리 3.18)입니다. 유도를 다시 써 보면, x(t)=∑kckejkω0tx(t) = \sum_k c_ke^{jk\omega_0 t}라 가정하고 양변에 e−jmω0te^{-jm\omega_0 t}를 곱해 한 주기 적분합니다.

∫Tx(t)e−jmω0t dt=∑kck∫Tej(k−m)ω0t dt=∑kckTδkm=Tcm\int_T x(t)e^{-jm\omega_0 t}\, dt = \sum_k c_k\int_T e^{j(k - m)\omega_0 t}\, dt = \sum_k c_kT\delta_{km} = Tc_m

mm번째 항만 살아남습니다. 무한합과 적분의 순서를 바꾼 것은 5.6절에서 정당화합니다.

계수의 뜻. c0=1T∫Tx dtc_0 = \frac{1}{T}\int_T x\, dt는 한 주기 평균(직류 성분)입니다. ckc_k는 xx에 e−jkω0te^{-jk\omega_0 t}를 곱해 평균 낸 것입니다. xx 안에 ejkω0te^{jk\omega_0 t} 성분이 있으면 e−jkω0te^{-jk\omega_0 t}를 곱했을 때 그 성분만 회전이 멈춰 상수가 되고, 나머지 성분은 여전히 돌고 있어 평균이 0이 됩니다.

정리 5.2 (유한 푸리에 합은 최선의 근사). SK(t)=∑k=−KKckejkω0tS_K(t) = \sum_{k=-K}^K c_ke^{jk\omega_0 t}는, ∑k=−KKakejkω0t\sum_{k=-K}^K a_ke^{jk\omega_0 t} 꼴의 함수 중에서 평균제곱 오차 ∫T∣x(t)−∑akejkω0t∣2 dt\int_T |x(t) - \sum a_ke^{jk\omega_0 t}|^2\, dt를 최소로 하는 유일한 것입니다.

증명. SKS_K는 xx의 span⁡{e−K,…,eK}\operatorname{span}\{e_{-K}, \ldots, e_K\} 위로의 정사영입니다. 3장의 최선근사 정리 3.23입니다. ■\blacksquare

5.3 실수 신호: 세 가지 표기

x(t)x(t)가 실수이면 c−k=1T∫Txejkω0t dt=ck∗c_{-k} = \frac{1}{T}\int_T xe^{jk\omega_0 t}\, dt = c_k^*입니다. ±k\pm k항을 묶으면 ckejkω0t+ck∗e−jkω0t=2Re⁡{ckejkω0t}c_ke^{jk\omega_0 t} + c_k^*e^{-jk\omega_0 t} = 2\operatorname{Re}\{c_ke^{jk\omega_0 t}\}이므로 다음 세 꼴을 얻습니다.

x(t)=∑k=−∞∞ckejkω0t(복소지수형)=c0+∑k=1∞2∣ck∣cos⁡(kω0t+∠ck)(진폭-위상형)=a02+∑k=1∞(akcos⁡kω0t+bksin⁡kω0t)(삼각함수형)\begin{aligned} x(t) &= \sum_{k=-\infty}^{\infty} c_ke^{jk\omega_0 t} && \text{(복소지수형)} \\ &= c_0 + \sum_{k=1}^{\infty} 2|c_k|\cos(k\omega_0 t + \angle c_k) && \text{(진폭-위상형)} \\ &= \frac{a_0}{2} + \sum_{k=1}^{\infty}\left(a_k\cos k\omega_0 t + b_k\sin k\omega_0 t\right) && \text{(삼각함수형)} \end{aligned} ak=2T∫Tx(t)cos⁡kω0t dt,bk=2T∫Tx(t)sin⁡kω0t dt,ck=ak−jbk2 (k≥1),c0=a02a_k = \frac{2}{T}\int_T x(t)\cos k\omega_0 t\, dt, \qquad b_k = \frac{2}{T}\int_T x(t)\sin k\omega_0 t\, dt, \qquad c_k = \frac{a_k - jb_k}{2}\ (k \ge 1), \qquad c_0 = \frac{a_0}{2}

ck=ak−jbk2c_k = \frac{a_k - jb_k}{2}는 e−jkω0t=cos⁡kω0t−jsin⁡kω0te^{-jk\omega_0 t} = \cos k\omega_0 t - j\sin k\omega_0 t를 계수 공식에 넣으면 바로 나옵니다. 세 꼴은 같은 내용입니다. 계산과 이론은 복소지수형이 가장 쉽고, 물리적 해석(“kk번째 고조파의 진폭은 2∣ck∣2|c_k|”)은 진폭-위상형이 가장 쉽습니다.

5.4 예제: 계수를 실제로 구하기

예제 5.3 (사각파). 주기 2π2\pi (ω0=1\omega_0 = 1), x(t)=1x(t) = 1 (0<t<π0 < t < \pi), x(t)=−1x(t) = -1 (−π<t<0-\pi < t < 0)이라 합시다.

ck=12π[∫0πe−jkt dt−∫−π0e−jkt dt]=12π[1−e−jkπjk−ejkπ−1jk]=2−2(−1)k2πjk={2jπk(k 홀수)0(k 짝수)c_k = \frac{1}{2\pi}\left[\int_0^\pi e^{-jkt}\, dt - \int_{-\pi}^0 e^{-jkt}\, dt\right] = \frac{1}{2\pi}\left[\frac{1 - e^{-jk\pi}}{jk} - \frac{e^{jk\pi} - 1}{jk}\right] = \frac{2 - 2(-1)^k}{2\pi jk} = \begin{cases} \dfrac{2}{j\pi k} & (k\ \text{홀수}) \\ 0 & (k\ \text{짝수}) \end{cases}

입니다(e±jkπ=(−1)ke^{\pm jk\pi} = (-1)^k이고, 평균이 0이므로 c0=0c_0 = 0). ±k\pm k를 묶으면 2jπk(ejkt−e−jkt)=4πksin⁡kt\frac{2}{j\pi k}(e^{jkt} - e^{-jkt}) = \frac{4}{\pi k}\sin kt이므로

x(t)=4π(sin⁡t+sin⁡3t3+sin⁡5t5+⋯ )x(t) = \frac{4}{\pi}\left(\sin t + \frac{\sin 3t}{3} + \frac{\sin 5t}{5} + \cdots\right)

입니다.

예제 5.4 (톱니파). 주기 2π2\pi, x(t)=tx(t) = t (−π<t<π-\pi < t < \pi)라 합시다. k≠0k \ne 0일 때 부분적분(u=tu = t, v′=e−jktv' = e^{-jkt})으로

∫−ππte−jkt dt=[te−jkt−jk]−ππ+1jk∫−ππe−jkt dt⏟0=π(−1)k+π(−1)k−jk=j2π(−1)kk\int_{-\pi}^{\pi} te^{-jkt}\, dt = \left[\frac{te^{-jkt}}{-jk}\right]_{-\pi}^{\pi} + \frac{1}{jk}\underbrace{\int_{-\pi}^{\pi} e^{-jkt}\, dt}_{0} = \frac{\pi(-1)^k + \pi(-1)^k}{-jk} = \frac{j2\pi(-1)^k}{k}

입니다. 따라서 ck=j(−1)kkc_k = \frac{j(-1)^k}{k} (k≠0k \ne 0), c0=0c_0 = 0이고, c−k=−ckc_{-k} = -c_k이므로 ckejkt+c−ke−jkt=2jcksin⁡ktc_ke^{jkt} + c_{-k}e^{-jkt} = 2jc_k\sin kt입니다.

x(t)=2(sin⁡t−sin⁡2t2+sin⁡3t3−⋯ )x(t) = 2\left(\sin t - \frac{\sin 2t}{2} + \frac{\sin 3t}{3} - \cdots\right)

사각파의 푸리에 부분합을 K = 1, 3, 9, 49인 네 그래프로 그렸습니다. 회색 선이 원래 사각파입니다. K가 클수록 평평한 부분은 잘 맞지만 불연속점 옆에 튀어나온 봉우리가 남습니다.

그림 5.1. 사각파의 푸리에 부분합. KK는 포함한 최고 고조파 차수. 항이 늘수록 평평한 부분은 잘 맞지만 불연속점 옆의 튀어나옴은 사라지지 않습니다(깁스 현상, 5.6절).

예제 5.5 (삼각파). 주기 2π2\pi, x(t)=1−2∣t∣πx(t) = 1 - \frac{2|t|}{\pi} (−π≤t≤π-\pi \le t \le \pi)라 합시다. 우함수이므로 ck=12π∫−ππxcos⁡kt dt=1π∫0π(1−2tπ)cos⁡kt dtc_k = \frac{1}{2\pi}\int_{-\pi}^{\pi} x\cos kt\, dt = \frac{1}{\pi}\int_0^\pi (1 - \frac{2t}{\pi})\cos kt\, dt입니다. k≠0k \ne 0이면 ∫0πcos⁡kt dt=0\int_0^\pi \cos kt\, dt = 0이고, 부분적분으로 ∫0πtcos⁡kt dt=[tsin⁡ktk]0π−1k∫0πsin⁡kt dt=(−1)k−1k2\int_0^\pi t\cos kt\, dt = \left[\frac{t\sin kt}{k}\right]_0^\pi - \frac{1}{k}\int_0^\pi \sin kt\, dt = \frac{(-1)^k - 1}{k^2}입니다. 따라서

ck=−2π2⋅(−1)k−1k2={4π2k2(k 홀수)0(k 짝수),x(t)=8π2(cos⁡t+cos⁡3t9+cos⁡5t25+⋯ )c_k = -\frac{2}{\pi^2} \cdot \frac{(-1)^k - 1}{k^2} = \begin{cases} \dfrac{4}{\pi^2k^2} & (k\ \text{홀수}) \\ 0 & (k\ \text{짝수}) \end{cases}, \qquad x(t) = \frac{8}{\pi^2}\left(\cos t + \frac{\cos 3t}{9} + \frac{\cos 5t}{25} + \cdots\right)

입니다.

예제 5.6 (펄스열). 주기 TT, 한 주기에서 x(t)=1x(t) = 1 (∣t∣<τ/2|t| < \tau/2), 나머지는 0 (τ<T\tau < T)이라 합시다. 2장 예제 2.7의 적분으로

ck=1T∫−τ/2τ/2e−jkω0t dt=1T⋅2sin⁡(kω0τ/2)kω0(k≠0),c0=τTc_k = \frac{1}{T}\int_{-\tau/2}^{\tau/2} e^{-jk\omega_0 t}\, dt = \frac{1}{T} \cdot \frac{2\sin(k\omega_0\tau/2)}{k\omega_0} \quad (k \ne 0), \qquad c_0 = \frac{\tau}{T}

입니다. TckTc_k는 함수 2sin⁡(ωτ/2)ω\frac{2\sin(\omega\tau/2)}{\omega}를 ω=kω0\omega = k\omega_0에서 읽은 값입니다. 다음 장의 출발점이 됩니다.

최고 고조파 차수 K를 바꾸면 부분합의 모양과 최대 넘침이 바뀝니다

이 위젯은 JavaScript가 필요합니다

5.5 대칭성과 기본 성질

계산을 줄여 주는 사실들입니다. xx는 실수, 주기 TT라 합시다.

  • 우함수(even function)이면 bk=0b_k = 0이고 ckc_k는 실수입니다(cos 항만). 기함수(odd function)이면 ak=0a_k = 0이고 ckc_k는 순허수입니다(sin 항만). 사각파와 톱니파는 기함수, 삼각파는 우함수였습니다.
  • 반파 대칭 x(t+T/2)=−x(t)x(t + T/2) = -x(t)이면 짝수 차 고조파가 없습니다. ckc_k의 적분을 반 주기씩 나누면 뒤쪽 절반이 앞쪽의 −e−jkπ=−(−1)k-e^{-jk\pi} = -(-1)^k배이므로 kk가 짝수일 때 상쇄됩니다. 사각파와 삼각파가 그렇습니다.

정리 5.7. x(t)x(t)의 푸리에 계수가 ckc_k일 때 다음이 성립합니다.

  1. (시간 이동) x(t−t0)x(t - t_0)의 계수는 e−jkω0t0cke^{-jk\omega_0 t_0}c_k입니다.
  2. (미분) xx가 연속이고 구간별로 매끄러우면 x′(t)x'(t)의 계수는 jkω0ckjk\omega_0c_k입니다.

증명. (1) 1T∫Tx(t−t0)e−jkω0t dt\frac{1}{T}\int_T x(t - t_0)e^{-jk\omega_0 t}\, dt에서 s=t−t0s = t - t_0으로 치환하면 e−jkω0t0⋅1T∫Tx(s)e−jkω0s dse^{-jk\omega_0 t_0} \cdot \frac{1}{T}\int_T x(s)e^{-jk\omega_0 s}\, ds입니다.

(2) 부분적분하면 1T∫0Tx′(t)e−jkω0t dt=1T[x(t)e−jkω0t]0T+jkω0T∫0Tx(t)e−jkω0t dt\frac{1}{T}\int_0^T x'(t)e^{-jk\omega_0 t}\, dt = \frac{1}{T}\Big[x(t)e^{-jk\omega_0 t}\Big]_0^T + \frac{jk\omega_0}{T}\int_0^T x(t)e^{-jk\omega_0 t}\, dt입니다. xx와 e−jkω0te^{-jk\omega_0 t} 모두 주기 TT이므로 경계항은 0입니다. ■\blacksquare

4장 DFT의 이동 정리(따름정리 4.11)와 같은 모양입니다. (2)는 급수를 항별로 미분한 것과 같습니다: ddtejkω0t=jkω0ejkω0t\frac{d}{dt}e^{jk\omega_0 t} = jk\omega_0e^{jk\omega_0 t}.

예제 5.8. 삼각파(예제 5.5)를 미분하면 x′(t)=∓2πx'(t) = \mp\frac{2}{\pi}로, 사각파의 −2π-\frac{2}{\pi}배입니다. 사각파의 계수 2jπk\frac{2}{j\pi k}에서 삼각파의 계수는 ck=1jk⋅(−2π)⋅2jπk=4π2k2c_k = \frac{1}{jk} \cdot \left(-\frac{2}{\pi}\right) \cdot \frac{2}{j\pi k} = \frac{4}{\pi^2k^2}입니다. 예제 5.5의 직접 계산과 일치합니다.

5.6 수렴

급수 ∑kckejkω0t\sum_k c_ke^{jk\omega_0 t}가 정말 x(t)x(t)와 같을까요? 답은 “어떤 뜻에서 같은가”에 달려 있습니다.

5.6.1 평균제곱 수렴

정리 5.9 (완비성). x∈L2[0,T]x \in L^2[0, T]이면 lim⁡K→∞∫T∣x(t)−SK(t)∣2 dt=0\displaystyle\lim_{K \to \infty}\int_T |x(t) - S_K(t)|^2\, dt = 0입니다.

유한 에너지를 가진 모든 주기신호에 대해 부분합의 평균제곱 오차는 0으로 갑니다. 즉 {ek}\{e_k\}는 완비이고 L2[0,T]L^2[0, T]의 기저입니다. 증명은 “연속인 주기함수는 삼각다항식으로 고르게 근사할 수 있습니다”(페예르의 정리)는 사실과 “L2L^2 함수는 연속함수로 근사할 수 있습니다”는 사실을 합쳐서 합니다. 여기서 다루는 범위를 넘으므로 결과만 받아들입니다.

5.6.2 점별 수렴

정리 5.10 (디리클레). xx가 구간별로 매끄러우면(유한 개의 점을 빼고 미분가능하고, 그 점들에서 좌우 극한과 좌우 미분계수가 존재하면) 모든 tt에서

lim⁡K→∞SK(t)=x(t+)+x(t−)2\lim_{K \to \infty} S_K(t) = \frac{x(t^+) + x(t^-)}{2}

입니다. 연속인 점에서는 x(t)x(t)로, 불연속점에서는 좌극한과 우극한의 평균으로 수렴합니다.

사각파의 급수에 t=0t = 0을 넣으면 모든 항이 0이고, 이는 1+(−1)2=0\frac{1 + (-1)}{2} = 0과 일치합니다.

증명의 뼈대 (xx가 tt에서 미분가능한 경우). 계수 공식을 부분합에 대입하고 합과 적분의 순서를 바꿉니다(유한합이므로 문제없습니다).

SK(t)=∑k=−KK[1T∫Tx(s)e−jkω0s ds]ejkω0t=1T∫Tx(s) DK(t−s) ds,DK(u)=∑k=−KKejkω0uS_K(t) = \sum_{k=-K}^K\left[\frac{1}{T}\int_T x(s)e^{-jk\omega_0 s}\, ds\right]e^{jk\omega_0 t} = \frac{1}{T}\int_T x(s)\, D_K(t - s)\, ds, \qquad D_K(u) = \sum_{k=-K}^K e^{jk\omega_0 u}

DKD_K를 디리클레 핵(Dirichlet kernel)이라 합니다. 등비수열의 합으로 닫힌 꼴을 구하면 (r=ejω0ur = e^{j\omega_0 u})

DK(u)=r−K1−r2K+11−r=r−(K+12)−rK+12r−12−r12=sin⁡((K+12)ω0u)sin⁡(ω0u/2)D_K(u) = r^{-K}\frac{1 - r^{2K+1}}{1 - r} = \frac{r^{-(K + \frac{1}{2})} - r^{K + \frac{1}{2}}}{r^{-\frac{1}{2}} - r^{\frac{1}{2}}} = \frac{\sin\left((K + \frac{1}{2})\omega_0u\right)}{\sin(\omega_0u/2)}

입니다. 1T∫TDK du=1\frac{1}{T}\int_T D_K\, du = 1이므로(k=0k = 0항만 남습니다) x(t)=1T∫Tx(t)DK(u) dux(t) = \frac{1}{T}\int_T x(t)D_K(u)\, du로 쓸 수 있고, u=t−su = t - s로 치환하면

SK(t)−x(t)=1T∫Tx(t−u)−x(t)sin⁡(ω0u/2)⏟g(u)sin⁡((K+12)ω0u)duS_K(t) - x(t) = \frac{1}{T}\int_T \underbrace{\frac{x(t - u) - x(t)}{\sin(\omega_0u/2)}}_{g(u)}\sin\left((K + \tfrac{1}{2})\omega_0u\right)du

입니다. xx가 tt에서 미분가능하면 u→0u \to 0일 때 g(u)→−2x′(t)/ω0g(u) \to -2x'(t)/\omega_0로 유한하므로 gg는 얌전한 함수입니다. 얌전한 함수에 아주 빨리 진동하는 사인파를 곱해 적분하면 양과 음이 상쇄되어 0으로 갑니다(리만–르베그 보조정리, 6장 정리 6.3). 따라서 SK(t)→x(t)S_K(t) \to x(t)입니다.

T = 2π일 때 디리클레 핵 D₅와 D₂₀의 그래프. 원점에 높이 2K + 1인 봉우리가 서고, K가 클수록 봉우리가 높고 좁아집니다.

그림 5.2. 디리클레 핵(T=2πT = 2\pi). KK가 커질수록 원점에 높고 좁은 봉우리가 서고 넓이는 TT로 일정합니다. 부분합은 xx를 이 핵으로 “평균 낸” 것입니다.

5.6.3 깁스 현상

불연속점 근처에서는 수렴이 고르지 않습니다. 그림 5.1에서 불연속점 옆의 튀어나온 봉우리는 KK를 키워도 높이가 줄지 않고 폭만 좁아집니다.

정리 5.11 (깁스 현상). 사각파(높이 ±1\pm 1, 도약 2)의 부분합의 최댓값은 K→∞K \to \infty일 때

2π∫0πsin⁡uu du=1.17898…\frac{2}{\pi}\int_0^\pi \frac{\sin u}{u}\, du = 1.17898\ldots

로 수렴합니다. 일반적으로 도약의 약 8.95%만큼 넘칩니다.

증명. 홀수 고조파 mm개의 합 S(t)=4π∑i=1msin⁡(2i−1)t2i−1S(t) = \frac{4}{\pi}\sum_{i=1}^m \frac{\sin(2i - 1)t}{2i - 1}을 미분하면 S′(t)=4π∑i=1mcos⁡(2i−1)t=2π⋅sin⁡2mtsin⁡tS'(t) = \frac{4}{\pi}\sum_{i=1}^m \cos(2i - 1)t = \frac{2}{\pi} \cdot \frac{\sin 2mt}{\sin t}입니다(등비수열의 합으로 확인할 수 있습니다). t>0t > 0에서 처음으로 S′=0S' = 0이 되는 점은 t=π/2mt = \pi/2m이고 거기가 첫 봉우리입니다. 그 높이는

S(π2m)=2π∫0π/2msin⁡2mtsin⁡t dt=2π∫0πsin⁡u2msin⁡(u/2m) du→m→∞2π∫0πsin⁡uu duS\left(\frac{\pi}{2m}\right) = \frac{2}{\pi}\int_0^{\pi/2m}\frac{\sin 2mt}{\sin t}\, dt = \frac{2}{\pi}\int_0^\pi \frac{\sin u}{2m\sin(u/2m)}\, du \xrightarrow{m \to \infty} \frac{2}{\pi}\int_0^\pi \frac{\sin u}{u}\, du

입니다(u=2mtu = 2mt로 치환했고, m→∞m \to \infty일 때 2msin⁡(u/2m)→u2m\sin(u/2m) \to u). 적분값은 수치적으로 1.85194다. ■\blacksquare

모순은 없습니다. 각 고정된 t≠0t \ne 0에서는 봉우리가 지나가 버린 뒤 SK(t)→x(t)S_K(t) \to x(t)이고, 봉우리의 폭이 0으로 가므로 평균제곱 오차도 0으로 갑니다. 실용적 교훈은 이렇습니다. 불연속이 있는 신호를 유한한 대역폭으로 자르면 약 9%의 넘침과 울림(ringing)이 반드시 생깁니다. 필터 설계와 영상 압축에서 실제로 보이는 현상입니다.

불연속점 t = 0 부근을 확대한 그래프. K = 9, 49, 199인 부분합의 첫 봉우리가 모두 높이 1.179에 있고, K가 클수록 불연속점 쪽으로 다가갑니다.

그림 5.3. 불연속점 t=0t = 0 부근의 부분합. 봉우리의 높이는 1.179로 일정하고 위치만 t=π/(K+1)t = \pi/(K + 1)로 불연속점에 다가갑니다.

5.6.4 매끄러움과 계수의 감소 속도

세 예제를 비교하면 사각파와 톱니파(불연속)는 ∣ck∣∼1/k|c_k| \sim 1/k, 삼각파(연속, 꺾임)는 ∣ck∣∼1/k2|c_k| \sim 1/k^2입니다. 미분 성질(정리 5.7)이 이유를 설명합니다. xx가 pp번 미분가능하고 x(p)x^{(p)}가 유한 에너지이면 x(p)x^{(p)}의 계수 (jkω0)pck(jk\omega_0)^pc_k는 베셀 부등식에 의해 제곱합이 유한하므로, ckc_k는 1/kp1/k^p보다 빨리 줄어야 합니다.

함수가 매끄러울수록 고주파 성분이 적습니다. 불연속이면 ∣ck∣∼1/k|c_k| \sim 1/k, 꺾임이면 1/k21/k^2, 무한히 매끄러우면 어떤 거듭제곱보다 빠르게 감소합니다.

왼쪽은 삼각파와 K = 1, 3까지의 부분합으로, 두 항만으로 거의 맞습니다. 오른쪽은 사각파와 삼각파의 계수 크기를 로그-로그 눈금으로 그린 것으로, 기울기가 각각 −1과 −2인 직선 위에 점이 놓입니다.

그림 5.4. 왼쪽: 삼각파는 두 항만으로 거의 맞습니다. 오른쪽: 사각파와 삼각파의 계수 크기(로그-로그 눈금). 기울기가 −1과 −2다.

이 원리는 푸리에 변환에서도 그대로 성립하며(시간 영역의 매끄러움 ↔ 주파수 영역의 빠른 감쇠), 창함수 설계(11장)의 근거가 됩니다.

5.7 파세발 등식

정리 5.12 (푸리에 급수의 파세발 등식). x,y∈L2[0,T]x, y \in L^2[0, T]의 푸리에 계수를 ckc_k, dkd_k라 하면

1T∫Tx(t) y(t)∗ dt=∑k=−∞∞ckdk∗,1T∫T∣x(t)∣2 dt=∑k=−∞∞∣ck∣2\frac{1}{T}\int_T x(t)\, y(t)^*\, dt = \sum_{k=-\infty}^{\infty} c_kd_k^*, \qquad \frac{1}{T}\int_T |x(t)|^2\, dt = \sum_{k=-\infty}^{\infty}|c_k|^2

입니다.

증명. 정규직교계 ek/Te_k/\sqrt{T}에 대한 xx의 좌표는 ⟨x,ek⟩/T=Tck\langle x, e_k \rangle/\sqrt{T} = \sqrt{T}c_k입니다. 3장 정리 3.23(3)에서 ∥x∥2=∥SK∥2+∥x−SK∥2=T∑∣k∣≤K∣ck∣2+∥x−SK∥2\lVert x\rVert^2 = \lVert S_K\rVert^2 + \lVert x - S_K\rVert^2 = T\sum_{|k| \le K}|c_k|^2 + \lVert x - S_K\rVert^2이고, 완비성에 의해 마지막 항이 0으로 갑니다. 내적에 대한 식도 같은 방법으로 나옵니다. ■\blacksquare

좌변은 신호의 평균 전력입니다. 평균 전력은 각 고조파의 전력 ∣ck∣2|c_k|^2의 합입니다. 주파수가 다른 성분끼리는 직교하므로 전력에 교차항이 없습니다.

예제 5.13 (∑1/n2\sum 1/n^2). 사각파는 12π∫∣x∣2=1\frac{1}{2\pi}\int |x|^2 = 1이고 ∑k∣ck∣2=2∑k 홀수≥14π2k2\sum_k |c_k|^2 = 2\sum_{k\ \text{홀수} \ge 1}\frac{4}{\pi^2k^2}입니다. 따라서 ∑k 홀수1k2=π28\sum_{k\ \text{홀수}}\frac{1}{k^2} = \frac{\pi^2}{8}입니다. 전체 합을 Z=∑n≥11/n2Z = \sum_{n \ge 1} 1/n^2이라 하면 짝수 항의 합은 Z/4Z/4이므로 Z=π28+Z4Z = \frac{\pi^2}{8} + \frac{Z}{4}, 즉

∑n=1∞1n2=π26\sum_{n=1}^\infty \frac{1}{n^2} = \frac{\pi^2}{6}

입니다.

예제 5.14 (전력 분포). 사각파의 전력 중 기본파(k=±1k = \pm 1)가 차지하는 비율은 2⋅4π2=8π2≈81%2 \cdot \frac{4}{\pi^2} = \frac{8}{\pi^2} \approx 81\%, 3차까지 ≈90%\approx 90\%, 5차까지 ≈93%\approx 93\%입니다.

연습문제

문제 1 · 확인 문제

5.1. x(t)=cos⁡2tx(t) = \cos^2 t의 푸리에 계수를 적분 없이 구합니다. 주기를 π\pi (ω0=2\omega_0 = 2)로 봅니다.

계수 값
c0c_0 [?]
c1c_1 [?]
c−1c_{-1} [?]
c2c_2 [?]
힌트 1

cos⁡t=12(ejt+e−jt)\cos t = \frac{1}{2}(e^{jt} + e^{-jt})를 제곱합니다.

정답
1/2
1/4
1/4
0
풀이

cos⁡2t=12+14ej2t+14e−j2t\cos^2 t = \frac{1}{2} + \frac{1}{4}e^{j2t} + \frac{1}{4}e^{-j2t}입니다. 주기 π\pi, ω0=2\omega_0 = 2로 보면 c0=12c_0 = \frac{1}{2}, c±1=14c_{\pm 1} = \frac{1}{4}, 나머지는 0입니다. 이미 복소지수의 합이면 계수는 읽기만 하면 됩니다.

문제 2 · 확인 문제

5.2. 전파 정류 사인파 x(t)=∣sin⁡t∣x(t) = |\sin t| (주기 π\pi)의 푸리에 계수 ckc_k를 구하고 실수형 급수로 씁니다. c0c_0와 c1c_1을 소수 넷째 자리까지 적습니다.

계수 값
c0c_0 [?]
c1c_1 [?]
힌트 1

ω0=2\omega_0 = 2입니다. sin⁡t=12j(ejt−e−jt)\sin t = \frac{1}{2j}(e^{jt} - e^{-jt})를 넣고 e±jπ=−1e^{\pm j\pi} = -1을 씁니다.

정답
0.6366
-0.2122
풀이

ω0=2\omega_0 = 2이고 ck=1π∫0πsin⁡t e−j2kt dtc_k = \frac{1}{\pi}\int_0^\pi \sin t\, e^{-j2kt}\, dt입니다. sin⁡t=12j(ejt−e−jt)\sin t = \frac{1}{2j}(e^{jt} - e^{-jt})를 넣고 e±jπ=−1e^{\pm j\pi} = -1을 쓰면

∫0πsin⁡t e−j2kt dt=12j[−2j(1−2k)−−2−j(1+2k)]=11−2k+11+2k=21−4k2\int_0^\pi \sin t\, e^{-j2kt}\, dt = \frac{1}{2j}\left[\frac{-2}{j(1 - 2k)} - \frac{-2}{-j(1 + 2k)}\right] = \frac{1}{1 - 2k} + \frac{1}{1 + 2k} = \frac{2}{1 - 4k^2}

입니다. 따라서 ck=2π(1−4k2)c_k = \frac{2}{\pi(1 - 4k^2)}이고 ∣sin⁡t∣=2π−4π∑k=1∞cos⁡2kt4k2−1|\sin t| = \frac{2}{\pi} - \frac{4}{\pi}\sum_{k=1}^\infty \frac{\cos 2kt}{4k^2 - 1}입니다. c0=2π≈0.6366c_0 = \frac{2}{\pi} \approx 0.6366, c1=−23π≈−0.2122c_1 = -\frac{2}{3\pi} \approx -0.2122입니다. 정류된 파형의 평균(직류)은 2/π≈0.6372/\pi \approx 0.637이고, 남은 리플의 가장 낮은 주파수는 원래 사인파의 2배입니다.

문제 3 · 확인 문제

5.3. 삼각파(예제 5.5)에 파세발 등식을 적용하여 ∑k=1,3,5,…1k4\displaystyle\sum_{k = 1, 3, 5, \ldots}\frac{1}{k^4}을 구합니다. 소수 넷째 자리까지 적습니다.

힌트 1

평균 전력은 12π∫−ππ(1−2∣t∣π)2 dt=13\frac{1}{2\pi}\int_{-\pi}^{\pi}(1 - \frac{2|t|}{\pi})^2\, dt = \frac{1}{3}입니다.

정답
1.0147
풀이

평균 전력은 12π∫−ππ(1−2∣t∣π)2 dt=13\frac{1}{2\pi}\int_{-\pi}^{\pi}(1 - \frac{2|t|}{\pi})^2\, dt = \frac{1}{3}입니다. 계수 쪽은 2∑k 홀수16π4k42\sum_{k\ \text{홀수}}\frac{16}{\pi^4k^4}입니다. 따라서 ∑k 홀수1k4=π496≈1.0147\sum_{k\ \text{홀수}}\frac{1}{k^4} = \frac{\pi^4}{96} \approx 1.0147입니다.

문제 4 · 확인 문제

5.4 (가). 주기 TT, 폭 τ=T/4\tau = T/4인 펄스열(예제 5.6)에서 0이 되는 고조파를 모두 고릅니다.

  1. k=2k = 2
  2. k=4k = 4
  3. k=6k = 6
  4. k=8k = 8
힌트 1

ck∝sin⁡(kω0τ/2)=sin⁡(kπ/4)c_k \propto \sin(k\omega_0\tau/2) = \sin(k\pi/4)입니다.

정답
2,4
풀이

ck∝sin⁡(kω0τ/2)=sin⁡(kπ/4)c_k \propto \sin(k\omega_0\tau/2) = \sin(k\pi/4)이므로 k=4,8,12,…k = 4, 8, 12, \ldots가 0입니다.

문제 5

5.4 (나). 같은 펄스열에서 전체 전력 중 직류 성분의 비율은 몇 %일까요?

힌트 1

전체 전력은 1T∫T∣x∣2=τT\frac{1}{T}\int_T |x|^2 = \frac{\tau}{T}, 직류 전력은 c02c_0^2입니다.

정답
25
풀이

전체 전력은 1T∫∣x∣2=τT=14\frac{1}{T}\int |x|^2 = \frac{\tau}{T} = \frac{1}{4}이고 직류 전력은 c02=116c_0^2 = \frac{1}{16}입니다. 비율은 25%다.

문제 6

예제 5.14에서 사각파의 전력 중 기본파(k=±1k = \pm 1)가 차지하는 비율은 몇 %일까요? 정수로 반올림합니다.

정답
81
풀이

8π2≈0.8106\frac{8}{\pi^2} \approx 0.8106이므로 약 81%다.