fft matlab代码自写

当然可以！以下是一个简单的FFT（快速傅里叶变换）的MATLAB代码示例：

function X = myfft(x)
    N = length(x);
    if N <= 1
        X = x;
    else
        even = myfft(x(1:2:N));
        odd = myfft(x(2:2:N));
        factor = exp(-2i*pi*(0:N/2-1)/N);
        X = [even + factor.*odd, even - factor.*odd];
    end
end

这个代码实现了一个递归的FFT函数 `myfft`，它输入一个向量 `x`，输出它的FFT结果 `X`。代码首先检查向量的长度是否小于等于 1，如果是，则直接返回向量本身。否则，它将向量分成偶数和奇数部分，并对它们进行递归调用 `myfft`。最后，根据傅里叶变换的定义，将两部分结果合并起来。

你可以使用以下代码测试这个函数：

x = [1, 2, 3, 4];
X = myfft(x);
disp(X);

输出结果应该是 `[10, -2+2i, -2, -2-2i]`，与MATLAB内置的 `fft` 函数的结果相同。

希望这可以帮助到你！如果你有任何其他问题，请随时提问。