pub fn dft_matrix(n: usize) -> Vec<Vec<Complex>>
Compute the N×N DFT matrix W where W[j][k] = exp(-2πi·j·k/N).
The matrix form of the DFT: X = W·x.