下面给程序包含在您的工作目录当你执行计划
函数X = dft (xn N)
%计算序列的DFT x (n)
L =长度(xn);%的序列长度
%检查DFT的长度
如果(N < L)
错误(必须> ' N = L ')
结束
N-L x1 = [xn 0 (1)];%附加零
%计算抚弄因素
为k = 0:1: n - 1;
为n = 0:1: n - 1;
p = exp(我* 2 *π* n * k / n);
x2 (k + 1, n + 1) = p;
结束
结束
X = x1 * x2;