我试图使代码非常简单:
首先创建这一维DFT函数
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
函数(Xk) = dft (xn)
N =长度(xn);
n = 0:1: n - 1;% n行向量
k = 0:1: n - 1;% k行向量
WN = exp (1 j * 2 *π/ N);%抚弄因子(w)
nk = n ' * k;%创建一个N×N矩阵nk的价值观
WNnk = WN。^ nk;% DFT矩阵
Xk = (WNnk * xn);
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
现在,使用这个函数的第二个维度:
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
函数= dft2 (x)
y = 0(大小(x));
日元= y;
C =大小(x, 2);%的列数
c = 1: c
y (:, c) = dft (x (:, c));
结束
R =大小(x, 1);%的行数
为r = 1: r
日元(r:) = dft (y (r:)。');
结束
= y₁;
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

注意:不要忘记保存dft的函数的名称。m和dft2。米都在同一个文件夹中。现在让你当前的工作目录的文件夹,现在,你可以叫dft一维或二维dft

享受
Montadar