Skip to main content

Derive the discrete-time Fourier series

Short answer

For period N, only N distinct exponentials exist (e^(jk2πn/N) repeats every N). Synthesis: x[n]=∑_{k=<N>}a_k e^(jk2πn/N). Analysis: a_k=(1/N)∑ x[n]e^(-jk2πn/N), obtained using geometric series orthogonality: ∑ e^(j(k-r)2πn/N) equals N if k≡r mod N, else 0.

Derive the exponential form of the Fourier series (synthesis and analysis equations) for a discrete time periodic signal x[n] with period N.

A discrete-time signal x[n]x[n] is periodic with period NN if x[n+N]=x[n]x[n+N] = x[n]; its fundamental frequency is ω0=2π/N\omega_0 = 2\pi/N. The DTFS writes it as a sum of harmonically related complex exponentials ϕk[n]=ejkω0n\phi_k[n] = e^{jk\omega_0 n}.

Only NN distinct exponentials

ϕk+N[n]=ej(k+N)(2π/N)n=ejk(2π/N)n ej2πn=ϕk[n]\phi_{k+N}[n] = e^{j(k+N)(2\pi/N)n} = e^{jk(2\pi/N)n}\, e^{j2\pi n} = \phi_k[n]

since ej2πn=1e^{j2\pi n} = 1 for integer nn. So only NN different exponentials exist, e.g. k=0,1,…,N−1k = 0, 1, \dots, N-1. The sum therefore runs over any NN consecutive values of kk, written k=⟨N⟩k = \langle N\rangle.

Synthesis equation

x[n]=∑k=⟨N⟩ak ejkω0n=∑k=⟨N⟩ak ejk(2π/N)nx[n] = \sum_{k=\langle N\rangle} a_k\, e^{jk\omega_0 n} = \sum_{k=\langle N\rangle} a_k\, e^{jk(2\pi/N)n}

Orthogonality result

For integer rr, using the finite geometric series:

∑n=0N−1ejr(2π/N)n={N,r=0,±N,±2N,…1−ej2πr1−ej2πr/N=0,otherwise\sum_{n=0}^{N-1} e^{jr(2\pi/N)n} = \begin{cases} N, & r = 0, \pm N, \pm 2N, \dots \\ \dfrac{1 - e^{j2\pi r}}{1 - e^{j2\pi r/N}} = 0, & \text{otherwise} \end{cases}

Analysis equation

Multiply the synthesis equation by e−jr(2π/N)ne^{-jr(2\pi/N)n} and sum over one period:

∑n=⟨N⟩x[n]e−jr(2π/N)n=∑n=⟨N⟩∑k=⟨N⟩akej(k−r)(2π/N)n=∑k=⟨N⟩ak∑n=⟨N⟩ej(k−r)(2π/N)n\begin{aligned} \sum_{n=\langle N\rangle}x[n]e^{-jr(2\pi/N)n} &= \sum_{n=\langle N\rangle}\sum_{k=\langle N\rangle}a_k e^{j(k-r)(2\pi/N)n} \\ &= \sum_{k=\langle N\rangle}a_k \sum_{n=\langle N\rangle}e^{j(k-r)(2\pi/N)n} \end{aligned}

The inner sum is NN when k=rk = r (within one period of kk) and 0 otherwise, so

∑n=⟨N⟩x[n]e−jr(2π/N)n=Nar\sum_{n=\langle N\rangle}x[n]e^{-jr(2\pi/N)n} = N a_r ak=1N∑n=⟨N⟩x[n] e−jk(2π/N)na_k = \frac{1}{N}\sum_{n=\langle N\rangle}x[n]\,e^{-jk(2\pi/N)n}

The DTFS pair

EquationFormula
Synthesisx[n]=∑k=⟨N⟩akejk(2π/N)nx[n] = \sum_{k=\langle N\rangle}a_k e^{jk(2\pi/N)n}
Analysisak=1N∑n=⟨N⟩x[n]e−jk(2π/N)na_k = \frac{1}{N}\sum_{n=\langle N\rangle}x[n]e^{-jk(2\pi/N)n}

Key points:

  • Both sums are finite (NN terms), so there is no convergence problem, unlike the CT Fourier series.
  • The coefficients are periodic: ak+N=aka_{k+N} = a_k.
  • a0=1N∑x[n]a_0 = \frac{1}{N}\sum x[n] is the average (dc) value.

Example: x[n]=sin⁡(2πn/5)x[n] = \sin(2\pi n/5), N=5N = 5: x[n]=12jej(2π/5)n−12je−j(2π/5)nx[n] = \frac{1}{2j}e^{j(2\pi/5)n} - \frac{1}{2j}e^{-j(2\pi/5)n}, so a1=12ja_1 = \frac{1}{2j}, a−1=−12ja_{-1} = -\frac{1}{2j}, and these repeat every 5 values of kk.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Signals and Systems chapter 2 questions →