DFT

gistfile1.txt

function dft{T}(x::Vector{T})
    N = length(x)
    X = Array(Complex{T}, N)
    for k = 1:N
        s = 0.0 + 0.0im
        for n = 1:N
            s += x[n] * exp(2π * (-im) * k * n / N)
        end
        X[k] = s
    end
    X ./= N
    X
end