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函数的源代码。