normaldft.m

标准dft

function [result]=normalDFT(x);
k=0:15;n=0:15; 
result=x*(exp(-j*pi/8)).^(n'*k);