MIT 6.300 — Lecture 3

Complex Fourier Series

§ Lecture: Complex Fourier Series

Annoyingly, we denote j:=1j := \sqrt{-1}. Fourier series look a lot simpler now.

f(t)=k=0(ckcos(kω0t)+dksin(kω0t))=k=akejkω0t.f(t) = \sum_{k = 0}^{\infty} \left( c_k \cos(k \omega_0 t) + d_k \sin(k \omega_0 t) \right) = \sum_{k = -\infty}^{\infty} a_k e^{j k \omega_0t}.

The relationship between the two looks like this:

ck=ak+a(k)  and  dk=j(aka(k))  for k1.      ak={c0 for k=0.12(ckjdk) for k>0.12(c(k)+jd(k)) for k<0.c_k = a_k + a_{(-k)} ~ \text{ and } ~ d_k = j\left(a_k - a_{(-k)}\right) ~ \text{ for } k \geq 1. ~~~ \parallel ~~~ a_k = \begin{cases} c_0 & \text{ for } k = 0. \\ \frac{1}{2}(c_k - jd_k) & \text{ for } k > 0. \\ \frac{1}{2}(c_{(-k)} + jd_{(-k)}) & \text{ for } k < 0. \end{cases}

And more directly, the coefficients {ak}k=\{a_k\}_{k = -\infty}^{\infty} are generated like so.

Theorem. (Complex Fourier Coefficients) Given a signal f(t)f(t) with period TT and fundamental frequency ω0:=2π/T\omega_0 := 2\pi / T, we have:

a=1TTf(t)ejω0tdt.a_{\ell} = \dfrac{1}{T} \int_T f(t) e^{-j \omega_0 \ell t} \, \mathrm{d}t.

Here, T\int_T denotes an integral over any interval of length TT. Note that this holds even for =0\ell = 0.

Proof: It's even simpler than last time.

Tejkω0tejω0tdt=Tejω0(k)tdt={T if k=.0 otherwise.\int_T e^{j k \omega_0 t} \cdot e^{-j\omega_0 \ell t} \, \mathrm{d}t = \int_T e^{j \omega_0 (k - \ell) t} \, \mathrm{d}t = \begin{cases} T & \text{ if } k = \ell. \\ 0 & \text{ otherwise.} \end{cases}

We can do the same “Transforming Fourier Series” analysis as before using complex numbers, too.

Problem. Given the Fourier coefficients {ak}k=\{a_k\}_{k = -\infty}^{\infty} for f(t)f(t), determine the Fourier coefficients for f(2t)f(2t) and f(tT/2)f(t - T/2).

Solution: For f(2t)f(2t), the coefficients remain the same. Meanwhile, for f(tT/2)f(t - T/2), we notice (recalling ω0=2π/T\omega_0 = 2\pi / T) that:

ejkω0(tT/2)=ejkω0T/2ejkω0t=(1)kejkω0t.e^{j k \omega_0 (t - T/2)} = e^{-j k \omega_0 T/2} \cdot e^{j k \omega_0 t} = (-1)^k \cdot e^{j k \omega_0 t}.

So for f(tT/2)f(t - T/2), the odd-indexed coefficients are negated, and all others remain the same.

The usefulness of complex numbers reveals itself when considering more complicated transformations.

Problem. Determine how the Fourier coefficients for f(t)f(t) relate to those of f(tT/8)f(t - T/8).

Solution: Sparing you the computational details, you can verify:

ck={ckk0(mod8).22(ckdk)k1(mod8).dkk2(mod8).22(ckdk)k3(mod8).ckk4(mod8).22(ck+dk)k5(mod8).dkk6(mod8).22(ck+dk)k7(mod8).c_k' = \begin{cases} c_k & k \equiv 0 \pmod{8}. \\ \frac{\sqrt{2}}{2}(c_k - d_k) & k \equiv 1 \pmod{8}. \\ -d_k & k \equiv 2 \pmod{8}. \\ \frac{\sqrt{2}}{2}(-c_k - d_k) & k \equiv 3 \pmod{8}. \\ -c_k & k \equiv 4 \pmod{8}. \\ \frac{\sqrt{2}}{2}(-c_k + d_k) & k \equiv 5 \pmod{8}. \\ d_k & k \equiv 6 \pmod{8}. \\ \frac{\sqrt{2}}{2}(c_k + d_k) & k \equiv 7 \pmod{8}. \end{cases}

The expression for dkd_k' is equally messy. In comparison, ak=ejk(2π/8)aka_k' = e^{-j k(2\pi / 8)} a_k; the kthk^{\text{th}} coefficient rotates clockwise by (45k)(45k)^{\circ}.

§ Recitation: Complex Fourier Series

In some cases, finding the complex Fourier series can be really easy.

Problem. Determine the Fourier series coefficients for f(t)=cos100(t)f(t) = \cos^{100}(t), taking T=2πT = 2\pi.

Solution: This is just (12ejt+12ejt)100\left(\frac{1}{2} e^{jt} + \frac{1}{2}e^{-jt}\right)^{100}, so ak=2100(10050k/2)a_k = 2^{-100} \binom{100}{50 - k / 2} for all k{100,98,,98,100}k \in \{-100, -98, \dots, 98, 100\}, and ak=0a_k = 0 otherwise.

Problem. Determine the Fourier series coefficients for f(t)=1(sin(3π2t)+jcos(3π2t))2f(t) = \dfrac{-1}{\left(\sin\left(\frac{3\pi}{2}t\right) + j\cos\left(\frac{3\pi}{2}t\right)\right)^2}, taking T=4T = 4.

Solution: It's just a6=1a_6 = 1 and everything else is zero.

As we saw in lecture, shifting a signal just multiplies each of its complex Fourier coefficients by a unit-magnitude constant.

Theorem. (Delay Property) If f(t)f(t) has Fourier coefficients {ak}k=\{a_k\}_{k = -\infty}^{\infty}, then f(tτ)f(t - \tau) has Fourier coefficients ak=akejkω0τa_k' = a_k \cdot e^{-jk\omega_0 \tau}.

We can also say things about symmetry and antisymmetry.

Might as well share a fun fact.

Theorem. (Parseval's Theorem) For a signal with Fourier coefficients {ak}k=\{a_k\}_{k = -\infty}^{\infty}, we have 1TTf(t)2dt=k=ak2\frac{1}{T} \int_T |f(t)|^2 \, \mathrm{d}t = \sum_{k = -\infty}^{\infty} |a_k|^2.

Proof: Just use zz=z2z \overline{z} = |z|^2 and the fact that {ejkω0t}k=\{e^{jk \omega_0 t}\}_{k = -\infty}^{\infty} forms an orthonormal basis under g,h:=1TTg(t)h(t)dt\langle g, h \rangle := \frac{1}{T} \int_T g(t) \overline{h(t)} \, \mathrm{d}t.   \blacksquare