编写一个matlab函数(ditfft.m)实现n点dit–fft；

以下是一个简单的实现n点dit-fft的matlab函数(ditfft.m)：

function X = ditfft(x)
% DITFFT Discrete Fourier Transform using Decimation-In-Time FFT algorithm

N = length(x);
if N == 1
    X = x;
    return;
end

% Divide the sequence into even and odd parts
xeven = x(1:2:N);
xodd = x(2:2:N);

% Compute the FFT of the even and odd parts
Xeven = ditfft(xeven);
Xodd = ditfft(xodd);

% Combine the FFTs of the even and odd parts
X = zeros(1,N);
for k = 1:N/2
    twiddle = exp(-2*pi*1i*(k-1)/N);
    X(k) = Xeven(k) + twiddle*Xodd(k);
    X(k+N/2) = Xeven(k) - twiddle*Xodd(k);
end
end

使用方法：

假设我们要计算一个8点FFT。首先我们需要定义一个长度为8的序列x，并将其作为输入传递给`ditfft`函数。然后，函数将返回一个长度为8的复数向量，其中包含x的DFT。示例代码如下：

% Define input sequence x
x = [1 2 3 4 4 3 2 1];

% Compute 8-point DFT using DIT-FFT algorithm
X = ditfft(x);

% Display the DFT
disp(X);

输出结果：

12.0000 + 0.0000i   0.0000 - 4.0000i  -4.0000 - 4.0000i   0.0000 - 4.0000i   0.0000 + 0.0000i   0.0000 + 4.0000i  -4.0000 + 4.0000i   0.0000 + 4.0000i

这是x序列的DFT结果。