LINEAR-ALGEBRA Β· Interactive Practice

Why Inverting the Fourier Matrix Is Almost Free

IKey Formulas

Formula Name Description
w=e2Ο€i/nw = e^{2\pi i / n} Primitive nn-th root of unity Smallest nn with wn=1w^n = 1
(Fn)ij=wij(F_n)_{ij} = w^{ij} Fourier matrix entry Indices i,j=0,…,nβˆ’1i,j = 0, \ldots, n-1
⟨u,v⟩=βˆ‘kuk‾ vk\langle u, v \rangle = \sum_k \overline{u_k}\, v_k Hermitian inner product Conjugate first, then multiply
Fnβˆ’1=1nFnHF_n^{-1} = \frac{1}{n} F_n^{H} Inverse of Fourier matrix Conjugate transpose, scaled by 1/n1/n

IIVisualization 1 β€” Roots of Unity

The powers of w=e2Ο€i/nw = e^{2\pi i/n} land on nn equally spaced points of the unit circle.

πŸ’‘ Multiplying by ww rotates by 2Ο€/n2\pi/n; after nn steps you return to wn=w0=1w^n = w^0 = 1 β€” the cyclic symmetry that fills FnF_n.

IIIVisualization 2 β€” The Fourier Matrix Entries

Every entry (Fn)ij=wij(F_n)_{ij} = w^{ij} is a power of ww; row 0 and column 0 are all ones.

IVVisualization 3 β€” Inner Product Between Columns

Distinct columns of F4F_4 are orthogonal β€” inner product 00; a column with itself gives n=4n = 4.

πŸ’‘ Distinct columns cancel to 00 and equal columns give nn, so FnHFn=nIF_n^H F_n = nI and Fnβˆ’1=1nFnHF_n^{-1} = \tfrac{1}{n} F_n^H.

VQuiz Questions

Question 1

Let w=e2Ο€i/4=iw = e^{2\pi i/4} = i be the primitive 4th root of unity. What is w6w^6?

βœ… Correct! w6=w6β€Šmodβ€Š4=w2=βˆ’1w^6 = w^{6 \bmod 4} = w^2 = -1.

❌ Not quite. Reduce the exponent mod 4 first, since w4=1w^4 = 1.

Show solution

Solution:

Since w4=1w^4 = 1, we can reduce the exponent modulo 4: w6=w6β€Šmodβ€Š4=w2=i2=βˆ’1w^6 = w^{6 \bmod 4} = w^2 = i^2 = -1

This is exactly how we fill in the entries of the Fourier matrix: any power of ww can be reduced to w0,w1,w2,w^0, w^1, w^2, or w3w^3.

Question 2

Consider columns 1 and 3 of F4F_4:

u=(1,i,βˆ’1,βˆ’i)u = (1, i, -1, -i) and v=(1,βˆ’i,βˆ’1,i)v = (1, -i, -1, i).

Using the Hermitian inner product ⟨u,v⟩=βˆ‘kukβ€Ύvk\langle u, v \rangle = \sum_k \overline{u_k} v_k, what is ⟨u,v⟩\langle u, v \rangle?

βœ… Correct! Distinct columns of FnF_n are orthogonal under the Hermitian inner product.

❌ Not quite. You may have skipped the conjugation step. The Hermitian inner product conjugates the first vector.

❌ Not quite. Conjugate uu first, then multiply entrywise and sum.

Show solution

Solution:

First conjugate uu: uΛ‰=(1,βˆ’i,βˆ’1,i)\bar{u} = (1, -i, -1, i).

Now multiply entry-wise with v=(1,βˆ’i,βˆ’1,i)v = (1, -i, -1, i) and sum:

⟨u,v⟩=(1)(1)+(βˆ’i)(βˆ’i)+(βˆ’1)(βˆ’1)+(i)(i)\langle u, v \rangle = (1)(1) + (-i)(-i) + (-1)(-1) + (i)(i) =1+i2+1+i2=1βˆ’1+1βˆ’1=0= 1 + i^2 + 1 + i^2 = 1 - 1 + 1 - 1 = 0

The columns are orthogonal under the Hermitian inner product. This is the key fact behind FnHFn=nIF_n^H F_n = nI.

Question 3

True or False: For the 4Γ—44 \times 4 Fourier matrix F4F_4, the matrix F4HF4F_4^H F_4 equals 4I4I (four times the identity).

βœ… Correct! This is exactly why F4βˆ’1=14F4HF_4^{-1} = \frac{1}{4}F_4^H.

❌ Not quite. Each column has squared length 4, and distinct columns are orthogonal, so F4HF4F_4^H F_4 is diagonal with 4s on the diagonal.

❌ Not quite. Try again β€” the hints above can help.

Show solution

Solution:

The entry (F4HF4)ij(F_4^H F_4)_{ij} is the Hermitian inner product of column ii with column jj of F4F_4.

  • Diagonal (i=ji = j): each column has squared length ∣1∣2+∣i∣2+βˆ£βˆ’1∣2+βˆ£βˆ’i∣2=4|1|^2 + |i|^2 + |-1|^2 + |-i|^2 = 4.
  • Off-diagonal (iβ‰ ji \ne j): distinct columns are orthogonal, so the inner product is 00.

Therefore F4HF4=4IF_4^H F_4 = 4I, and dividing by 4 gives: F4βˆ’1=14F4HF_4^{-1} = \frac{1}{4} F_4^H

Question 4

For a general nΓ—nn \times n Fourier matrix FnF_n, what is Fnβˆ’1F_n^{-1}?

βœ… Correct! Since FnHFn=nIF_n^H F_n = nI, the inverse is 1nFnH\frac{1}{n}F_n^H.

❌ Close, but missing the scale factor. Each column has squared length nn, not 1, so we need the 1/n1/n factor.

❌ Not quite. The key relation is FnHFn=nIF_n^H F_n = nI. Solve for Fnβˆ’1F_n^{-1}.

Show solution

Solution:

The columns of FnF_n are mutually orthogonal under the Hermitian inner product, and each has squared length nn (since ∣wkj∣2=1|w^{kj}|^2 = 1 for each of the nn entries).

Therefore: FnHFn=nIβ€…β€ŠβŸΉβ€…β€ŠFnβˆ’1=1nFnHF_n^H F_n = nI \implies F_n^{-1} = \frac{1}{n} F_n^H

This is why the inverse Discrete Fourier Transform is essentially the same as the forward DFT β€” just conjugate and divide by nn. No Gaussian elimination needed. This same orthogonality powers the Fast Fourier Transform, JPEG compression, and PDE solvers.

Solved: 0 / 4