利用matlab编写离散傅里叶变换程序

以下是一个简单的 Matlab 程序，用于计算离散傅里叶变换：

function y = my_dft(x)
% 离散傅里叶变换函数
% 输入参数 x：原始信号序列
% 输出参数 y：信号在频域的幅度

% 获取信号长度
N = length(x);

% 定义离散傅里叶变换系数矩阵
W_N = zeros(N,N);
for k = 0:N-1
    for n = 0:N-1
        W_N(k+1,n+1) = exp(-1i*2*pi*k*n/N);
    end
end

% 计算离散傅里叶变换
y = W_N * x';

% 取模计算信号在频域的幅度
y = abs(y);
end

在这个函数中，我们首先获取输入信号的长度，然后定义一个 N x N 的矩阵 W_N，用于存储离散傅里叶变换的系数。接下来，我们使用两个循环来计算 W_N 矩阵中的每个元素，其中 k 和 n 分别表示频域的索引和时域的索引。最后，我们将信号序列与 W_N 矩阵相乘，得到信号在频域的表示，并使用 abs 函数计算信号在频域的幅度。

你可以调用这个函数并传入你自己的信号序列来计算其离散傅里叶变换。例如：

x = [1, 2, 3, 4, 5, 6, 7, 8];
y = my_dft(x);
disp(y);

这个程序将输出信号 x 的离散傅里叶变换的幅度。