부록 A. 공식 모음

푸리에 해석 · 푸리에 변환

기본 · 약 30분

요약

네 변환의 정의, 변환쌍 표, 성질 표, 주파수(Hz) 관례로 옮기기, 푸리에 급수와 DFT, DFT 실무 공식, 자주 쓰는 계산, 더 읽을 책

자주 쓰는 정의, 변환쌍, 성질을 한곳에 모은 공식 모음입니다. 각 공식 옆에 그것을 다룬 장을 적어 두었으므로, 공식이 왜 성립하는지 궁금하면 그 장으로 가면 됩니다.

모든 공식은 관례 (A) (X(ω)=∫xe−jωt dtX(\omega) = \int xe^{-j\omega t}\, dt) 기준입니다. 괄호 안은 본문의 위치입니다.

A.1 정의

변환 분석 합성
푸리에 변환 (6장) X(ω)=∫−∞∞x(t)e−jωt dtX(\omega) = \int_{-\infty}^{\infty}x(t)e^{-j\omega t}\, dt x(t)=12π∫−∞∞X(ω)ejωt dωx(t) = \frac{1}{2\pi}\int_{-\infty}^{\infty}X(\omega)e^{j\omega t}\, d\omega
푸리에 급수 (5장) ck=1T∫Tx(t)e−jkω0t dtc_k = \frac{1}{T}\int_Tx(t)e^{-jk\omega_0t}\, dt x(t)=∑k=−∞∞ckejkω0tx(t) = \sum_{k=-\infty}^{\infty}c_ke^{jk\omega_0t}, ω0=2πT\omega_0 = \frac{2\pi}{T}
DTFT (11장) Xd(Ω)=∑n=−∞∞x[n]e−jΩnX_d(\Omega) = \sum_{n=-\infty}^{\infty}x[n]e^{-j\Omega n} x[n]=12π∫−ππXd(Ω)ejΩn dΩx[n] = \frac{1}{2\pi}\int_{-\pi}^{\pi}X_d(\Omega)e^{j\Omega n}\, d\Omega
DFT (4장) X[k]=∑n=0N−1x[n]e−j2πkn/NX[k] = \sum_{n=0}^{N-1}x[n]e^{-j2\pi kn/N} x[n]=1N∑k=0N−1X[k]ej2πkn/Nx[n] = \frac{1}{N}\sum_{k=0}^{N-1}X[k]e^{j2\pi kn/N}

A.2 푸리에 변환쌍

x(t)x(t) X(ω)X(\omega) 비고
δ(t)\delta(t) 11
δ(t−t0)\delta(t - t_0) e−jωt0e^{-j\omega t_0}
11 2πδ(ω)2\pi\delta(\omega)
ejω0te^{j\omega_0t} 2πδ(ω−ω0)2\pi\delta(\omega - \omega_0)
cos⁡ω0t\cos\omega_0t π[δ(ω−ω0)+δ(ω+ω0)]\pi[\delta(\omega - \omega_0) + \delta(\omega + \omega_0)]
sin⁡ω0t\sin\omega_0t πj[δ(ω−ω0)−δ(ω+ω0)]\frac{\pi}{j}[\delta(\omega - \omega_0) - \delta(\omega + \omega_0)]
u(t)u(t) πδ(ω)+1jω\pi\delta(\omega) + \frac{1}{j\omega}
sgn⁡(t)\operatorname{sgn}(t) 2jω\frac{2}{j\omega}
1πt\frac{1}{\pi t} −jsgn⁡(ω)-j\operatorname{sgn}(\omega) 힐베르트 변환
rect⁡(t/τ)\operatorname{rect}(t/\tau) 2sin⁡(ωτ/2)ω=τsinc⁡ωτ2π\frac{2\sin(\omega\tau/2)}{\omega} = \tau\operatorname{sinc}\frac{\omega\tau}{2\pi} 폭 τ\tau
tri⁡(t/τ)\operatorname{tri}(t/\tau) τsinc⁡2ωτ2π\tau\operatorname{sinc}^2\frac{\omega\tau}{2\pi} 밑변 2τ2\tau
sin⁡Wtπt\frac{\sin Wt}{\pi t} rect⁡ω2W\operatorname{rect}\frac{\omega}{2W} 이상적 저역통과
e−atu(t)e^{-at}u(t) 1a+jω\frac{1}{a + j\omega} Re⁡a>0\operatorname{Re}a > 0
tne−atu(t)t^ne^{-at}u(t) n!(a+jω)n+1\frac{n!}{(a + j\omega)^{n+1}} Re⁡a>0\operatorname{Re}a > 0
e−a∣t∣e^{-a\lvert t\rvert} 2aa2+ω2\frac{2a}{a^2 + \omega^2} a>0a > 0
1a2+t2\frac{1}{a^2 + t^2} πae−a∣ω∣\frac{\pi}{a}e^{-a\lvert\omega\rvert} a>0a > 0
e−t2/(2σ2)e^{-t^2/(2\sigma^2)} σ2π e−σ2ω2/2\sigma\sqrt{2\pi}\, e^{-\sigma^2\omega^2/2}
e−atcos⁡(ω0t)u(t)e^{-at}\cos(\omega_0t)u(t) a+jω(a+jω)2+ω02\frac{a + j\omega}{(a + j\omega)^2 + \omega_0^2} a>0a > 0
e−atsin⁡(ω0t)u(t)e^{-at}\sin(\omega_0t)u(t) ω0(a+jω)2+ω02\frac{\omega_0}{(a + j\omega)^2 + \omega_0^2} a>0a > 0
∑nδ(t−nT)\sum_n\delta(t - nT) 2πT∑kδ(ω−2πkT)\frac{2\pi}{T}\sum_k\delta\left(\omega - \frac{2\pi k}{T}\right) 임펄스열
∑kckejkω0t\sum_kc_ke^{jk\omega_0t} 2π∑kckδ(ω−kω0)2\pi\sum_kc_k\delta(\omega - k\omega_0) 주기신호

A.3 성질

성질 시간 영역 주파수 영역
선형성 αx+βy\alpha x + \beta y αX+βY\alpha X + \beta Y
시간 이동 x(t−t0)x(t - t_0) e−jωt0X(ω)e^{-j\omega t_0}X(\omega)
주파수 이동 ejω0tx(t)e^{j\omega_0t}x(t) X(ω−ω0)X(\omega - \omega_0)
변조 x(t)cos⁡ω0tx(t)\cos\omega_0t 12X(ω−ω0)+12X(ω+ω0)\frac{1}{2}X(\omega - \omega_0) + \frac{1}{2}X(\omega + \omega_0)
스케일링 x(at)x(at) 1∣a∣X(ω/a)\frac{1}{\lvert a\rvert}X(\omega/a)
켤레 x(t)∗x(t)^* X(−ω)∗X(-\omega)^*
쌍대성 X(t)X(t) 2πx(−ω)2\pi x(-\omega)
미분 x(n)(t)x^{(n)}(t) (jω)nX(ω)(j\omega)^nX(\omega)
적분 ∫−∞tx(s) ds\int_{-\infty}^tx(s)\, ds X(ω)jω+πX(0)δ(ω)\frac{X(\omega)}{j\omega} + \pi X(0)\delta(\omega)
tt 곱 tnx(t)t^nx(t) jnX(n)(ω)j^nX^{(n)}(\omega)
합성곱 (x∗y)(t)(x * y)(t) X(ω)Y(ω)X(\omega)Y(\omega)
곱 x(t)y(t)x(t)y(t) 12π(X∗Y)(ω)\frac{1}{2\pi}(X * Y)(\omega)
상관 ∫x(t+τ)y(t)∗ dt\int x(t + \tau)y(t)^*\, dt X(ω)Y(ω)∗X(\omega)Y(\omega)^*
표본화 x(t)∑nδ(t−nTs)x(t)\sum_n\delta(t - nT_s) 1Ts∑kX(ω−kωs)\frac{1}{T_s}\sum_kX(\omega - k\omega_s)
주기화 ∑nx(t+nT)\sum_nx(t + nT) 2πT∑kX(kω0)δ(ω−kω0)\frac{2\pi}{T}\sum_kX(k\omega_0)\delta(\omega - k\omega_0)
  • 파세발: ∫xy∗ dt=12π∫XY∗ dω\int xy^*\, dt = \frac{1}{2\pi}\int XY^*\, d\omega, ∫∣x∣2 dt=12π∫∣X∣2 dω\int|x|^2\, dt = \frac{1}{2\pi}\int|X|^2\, d\omega.
  • 넓이: X(0)=∫x(t) dtX(0) = \int x(t)\, dt, x(0)=12π∫X(ω) dωx(0) = \frac{1}{2\pi}\int X(\omega)\, d\omega.
  • 실수 신호: X(−ω)=X(ω)∗X(-\omega) = X(\omega)^*. 실수 우함수는 실수 우함수로, 실수 기함수는 순허수 기함수로 갑니다.
  • 불확정성: Δt Δω≥12\Delta t\, \Delta\omega \ge \frac{1}{2}.

A.4 관례 (B): 주파수 f (Hz)로 쓸 때

X^(f)=∫x(t)e−j2πft dt=X(2πf)\hat{X}(f) = \int x(t)e^{-j2\pi ft}\, dt = X(2\pi f)입니다. 델타함수는 δ(ω)=δ(2πf)=12πδ(f)\delta(\omega) = \delta(2\pi f) = \frac{1}{2\pi}\delta(f)로 바뀐다는 점에 주의합니다.

관례 (A): X(ω)X(\omega) 관례 (B): X^(f)\hat{X}(f)
역변환 12π∫X(ω)ejωt dω\frac{1}{2\pi}\int X(\omega)e^{j\omega t}\, d\omega ∫X^(f)ej2πft df\int\hat{X}(f)e^{j2\pi ft}\, df
11 2πδ(ω)2\pi\delta(\omega) δ(f)\delta(f)
cos⁡2πf0t\cos 2\pi f_0t π[δ(ω−ω0)+δ(ω+ω0)]\pi[\delta(\omega - \omega_0) + \delta(\omega + \omega_0)] 12[δ(f−f0)+δ(f+f0)]\frac{1}{2}[\delta(f - f_0) + \delta(f + f_0)]
rect⁡(t/τ)\operatorname{rect}(t/\tau) τsinc⁡(ωτ/2π)\tau\operatorname{sinc}(\omega\tau/2\pi) τsinc⁡(fτ)\tau\operatorname{sinc}(f\tau)
e−πt2e^{-\pi t^2} e−ω2/4πe^{-\omega^2/4\pi} e−πf2e^{-\pi f^2}
미분 jωX(ω)j\omega X(\omega) j2πfX^(f)j2\pi f\hat{X}(f)
쌍대성 X(t)↔2πx(−ω)X(t) \leftrightarrow 2\pi x(-\omega) X^(t)↔x(−f)\hat{X}(t) \leftrightarrow x(-f)
곱 12πX∗Y\frac{1}{2\pi}X * Y X^∗Y^\hat{X} * \hat{Y}
파세발 ∫∣x∣2=12π∫∣X∣2 dω\int\lvert x\rvert^2 = \frac{1}{2\pi}\int\lvert X\rvert^2\, d\omega ∫∣x∣2=∫∣X^∣2 df\int\lvert x\rvert^2 = \int\lvert\hat{X}\rvert^2\, df
임펄스열 (간격 TT) 2πT∑kδ(ω−2πk/T)\frac{2\pi}{T}\sum_k\delta(\omega - 2\pi k/T) 1T∑kδ(f−k/T)\frac{1}{T}\sum_k\delta(f - k/T)
불확정성 Δt Δω≥12\Delta t\, \Delta\omega \ge \frac{1}{2} Δt Δf≥14π\Delta t\, \Delta f \ge \frac{1}{4\pi}

시간 이동, 스케일링, 합성곱 정리는 ω\omega를 2πf2\pi f로 바꾸는 것 말고는 달라지지 않습니다.

A.5 푸리에 급수와 DFT

푸리에 급수 (주기 TT) DFT (길이 NN)
직교성 ∫Tej(k−m)ω0t dt=Tδkm\int_Te^{j(k - m)\omega_0t}\, dt = T\delta_{km} ∑nej2π(k−m)n/N=Nδkm\sum_ne^{j2\pi(k - m)n/N} = N\delta_{km}
이동 x(t−t0)→e−jkω0t0ckx(t - t_0) \to e^{-jk\omega_0t_0}c_k x[n−m]→e−j2πkm/NX[k]x[n - m] \to e^{-j2\pi km/N}X[k]
미분 x′(t)→jkω0ckx'(t) \to jk\omega_0c_k
합성곱 x⊛h→X[k]H[k]x \circledast h \to X[k]H[k]
실수 신호 c−k=ck∗c_{-k} = c_k^* X[N−k]=X[k]∗X[N - k] = X[k]^*
파세발 1T∫T∣x∣2 dt=∑k∣ck∣2\frac{1}{T}\int_T\lvert x\rvert^2\, dt = \sum_k\lvert c_k\rvert^2 ∑n∣x[n]∣2=1N∑k∣X[k]∣2\sum_n\lvert x[n]\rvert^2 = \frac{1}{N}\sum_k\lvert X[k]\rvert^2

주기 2π2\pi인 기본 파형의 급수(5장 5.4절)는 다음과 같습니다.

사각파(±1): 4π∑k=1,3,5,…sin⁡ktk,톱니파(x=t, ∣t∣<π): 2∑k=1∞(−1)k+1sin⁡ktk,삼각파(x=1−2∣t∣/π): 8π2∑k=1,3,5,…cos⁡ktk2\text{사각파}(\pm 1):\ \frac{4}{\pi}\sum_{k = 1, 3, 5, \ldots}\frac{\sin kt}{k}, \qquad \text{톱니파}(x = t,\ |t| < \pi):\ 2\sum_{k=1}^{\infty}(-1)^{k+1}\frac{\sin kt}{k}, \qquad \text{삼각파}(x = 1 - 2|t|/\pi):\ \frac{8}{\pi^2}\sum_{k = 1, 3, 5, \ldots}\frac{\cos kt}{k^2}

A.6 DFT 실무 공식

fk=kNfs (k<N/2),Δf=fsN=1NTs,fmax⁡=fs2,fa=∣f−fsround⁡(f/fs)∣f_k = \frac{k}{N}f_s\ (k < N/2), \qquad \Delta f = \frac{f_s}{N} = \frac{1}{NT_s}, \qquad f_{\max} = \frac{f_s}{2}, \qquad f_a = \left|f - f_s\operatorname{round}(f/f_s)\right| X(ωk)≈Tse−jωkt0X[k],사인파 진폭≈2∣X[k]∣∑nw[n],FFT 연산량∼Nlog⁡2NX(\omega_k) \approx T_se^{-j\omega_kt_0}X[k], \qquad \text{사인파 진폭} \approx \frac{2|X[k]|}{\sum_nw[n]}, \qquad \text{FFT 연산량} \sim N\log_2N

A.7 자주 쓰는 계산

ejθ=cos⁡θ+jsin⁡θ,cos⁡θ=12(ejθ+e−jθ),sin⁡θ=12j(ejθ−e−jθ)e^{j\theta} = \cos\theta + j\sin\theta, \qquad \cos\theta = \tfrac{1}{2}(e^{j\theta} + e^{-j\theta}), \qquad \sin\theta = \tfrac{1}{2j}(e^{j\theta} - e^{-j\theta}) ejπk=(−1)k,ej2πk=1,e±jπ/2=±j,∣ejθ∣=1e^{j\pi k} = (-1)^k, \qquad e^{j2\pi k} = 1, \qquad e^{\pm j\pi/2} = \pm j, \qquad |e^{j\theta}| = 1 ∑n=0N−1rn=1−rN1−r,∑k=−KKejkθ=sin⁡((K+12)θ)sin⁡(θ/2)\sum_{n=0}^{N-1}r^n = \frac{1 - r^N}{1 - r}, \qquad \sum_{k=-K}^{K}e^{jk\theta} = \frac{\sin((K + \frac{1}{2})\theta)}{\sin(\theta/2)} ∫0∞e−(a+jω)t dt=1a+jω,∫−aae−jωt dt=2sin⁡ωaω,∫0∞tne−at dt=n!an+1\int_0^\infty e^{-(a + j\omega)t}\, dt = \frac{1}{a + j\omega}, \qquad \int_{-a}^{a}e^{-j\omega t}\, dt = \frac{2\sin\omega a}{\omega}, \qquad \int_0^\infty t^ne^{-at}\, dt = \frac{n!}{a^{n+1}} ∫−∞∞e−at2 dt=πa,∫−∞∞sin⁡tt dt=π,∫−∞∞dta2+t2=πa\int_{-\infty}^{\infty}e^{-at^2}\, dt = \sqrt{\frac{\pi}{a}}, \qquad \int_{-\infty}^{\infty}\frac{\sin t}{t}\, dt = \pi, \qquad \int_{-\infty}^{\infty}\frac{dt}{a^2 + t^2} = \frac{\pi}{a} ∫−∞∞ejωt dt=2πδ(ω),∫−∞∞δ(t−t0)ϕ(t) dt=ϕ(t0),δ(at)=δ(t)∣a∣\int_{-\infty}^{\infty}e^{j\omega t}\, dt = 2\pi\delta(\omega), \qquad \int_{-\infty}^{\infty}\delta(t - t_0)\phi(t)\, dt = \phi(t_0), \qquad \delta(at) = \frac{\delta(t)}{|a|}

더 읽을 책

  • 미적분. J. Stewart, Calculus. 1장의 내용을 처음부터 연습문제와 함께 익히려면 이 책의 극한, 미분, 적분, 급수 장으로 충분합니다.
  • 선형대수. G. Strang, Introduction to Linear Algebra. 행렬을 열의 선형결합으로 보는 3장의 관점이 이 책의 방식입니다. 증명 중심으로 내적공간과 고유벡터를 다시 보려면 S. Axler, Linear Algebra Done Right.
  • 신호와 시스템. A. V. Oppenheim, A. S. Willsky, Signals and Systems, 2nd ed. 전기공학의 표준 교재이고 여기의 표기와 관례가 같습니다. 6–11장의 내용을 훨씬 많은 예제와 문제로 다룹니다.
  • 푸리에 변환 자체. B. Osgood, Lectures on the Fourier Transform and Its Applications. 델타함수와 표본화를 분포의 관점에서 깔끔하게 설명합니다. 관례 (B)를 씁니다. R. N. Bracewell, The Fourier Transform and Its Applications, 3rd ed.은 그림으로 된 변환쌍 사전이 유명합니다.
  • 수학적 엄밀성. E. M. Stein, R. Shakarchi, Fourier Analysis: An Introduction. 여기서 생략한 증명(완비성, 수렴, 역변환 정리)이 모두 있습니다.
  • 이산 신호처리. A. V. Oppenheim, R. W. Schafer, Discrete-Time Signal Processing, 3rd ed. DTFT, 창함수, 필터 설계, FFT를 다룹니다. 실용적인 입문서로는 S. W. Smith, The Scientist and Engineer's Guide to Digital Signal Processing.

쓰이는 곳. 논문을 읽다 변환쌍이나 성질이 기억나지 않을 때 여기서 찾으면 됩니다. 정규화 상수(2π2\pi를 어디에 두는지)가 논문마다 다르므로, 6장의 관례 절과 함께 보면 좋습니다.

참고 문헌

  1. Calculus, J. Stewart
  2. Introduction to Linear Algebra, G. Strang
  3. Linear Algebra Done Right, S. Axler
  4. Signals and Systems, 2nd ed., A. V. Oppenheim, A. S. Willsky
  5. Lectures on the Fourier Transform and Its Applications, B. Osgood
  6. The Fourier Transform and Its Applications, 3rd ed., R. N. Bracewell
  7. Fourier Analysis: An Introduction, E. M. Stein, R. Shakarchi
  8. Discrete-Time Signal Processing, 3rd ed., A. V. Oppenheim, R. W. Schafer
  9. The Scientist and Engineer's Guide to Digital Signal Processing, S. W. Smith