pub fn naive_dft(xs: &[Complex]) -> Vec<Complex>
Apply the DFT matrix to a vector (O(N²) naive DFT).
Returns X[k] = Σ x[n] · exp(-2πi·k·n/N).