matlab编写函数实现傅里叶逆变换

以下是一个 Matlab 函数，用于实现傅里叶逆变换：

function [x, t] = my_ifft(X, fs, n)
% X: 频谱图 STFT(k, l)
% fs: 采样率
% n: FFT 点数
% x: 时域重构信号
% t: 时间向量

% 计算频率向量
df = fs / n;
f = (0:n/2-1) * df;

% 计算时间向量
dt = n / fs;
t = (0:size(X, 2)-1) * dt;

% 进行逆变换
x = zeros(n, size(X, 2));
for i = 1:size(X, 2)
    x(:, i) = ifft(X(:, i), n);
end

% 取实部并归一化
x = real(x(1:n/2, :)) * 2 / n;

end

该函数的输入参数为频谱图 STFT(k, l)、采样率 fs 和 FFT 点数 n，输出参数为时域重构信号 x 和时间向量 t。函数首先计算频率向量和时间向量，然后对每一列进行傅里叶逆变换，并取实部并归一化得到时域重构信号。