a=(0:15)'*(0:15);W=exp(-2*pi*1i*a/16)/sqrt(16);不会太麻烦，两句搞定推广到不同N值情况a=(0:N-1)'*(0:N-1);W=exp(-2*pi*1i*a/N)/sqrt(N);
直接用 fft(eye(n));就OK了，实际上dftmtx（n）就是调用的fft(eye(n));的，你可以看一下dftmtx函数的源代码。
