设计一个matlab的fft程序

以下是一个简单的MATLAB FFT程序：

function y = my_fft(x)
% MY_FFT  Calculates the discrete Fourier transform (DFT) of a signal x
% using the fast Fourier transform (FFT) algorithm.
% Usage: y = my_fft(x)
% Inputs: x - input signal
% Output: y - DFT of x
N = length(x); % length of input signal
if N == 1
    y = x;
else
    xe = my_fft(x(1:2:N-1)); % compute FFT of even-indexed elements
    xo = my_fft(x(2:2:N)); % compute FFT of odd-indexed elements
    W = exp(-2i*pi/N).^(0:N/2-1); % twiddle factors
    y = [xe + W.*xo, xe - W.*xo]; % combine FFTs with twiddle factors
end
end

该程序使用递归地将输入信号分解为偶数和奇数元素，并使用旋转因子将它们组合为FFT。注意，这个程序仅适用于长度为2的幂的信号。如果输入信号的长度不是2的幂，则可以使用零填充来将其扩展为最接近的2的幂次方。