用matlab函数形式写一个 DFT 的计算机实现方法

function X = myDFT(x)
% x: 输入的信号序列
% X: 输出的频率域序列

N = length(x); % 信号序列长度
X = zeros(1,N); % 初始化频率域序列

for k = 0:N-1
for n = 0:N-1
X(k+1) = X(k+1) + x(n+1)*exp(-1j*2*pi*k*n/N);
end
end

end