matlab 8点fft蝶形图

以下是使用 MATLAB 实现 8 点 FFT 的蝶形图代码：

% 输入信号
x = [1 2 3 4 4 3 2 1];

% 定义旋转因子
W8 = exp(-1j*2*pi/8);

% 定义蝶形运算
butterfly = @(x) [x(1)+x(2) (x(1)-x(2))*W8];

% 第一级蝶形运算
X1 = [butterfly(x(1:2)); butterfly(x(3:4)); butterfly(x(5:6)); butterfly(x(7:8))];

% 第二级蝶形运算
X2 = [butterfly(X1(1,:)); butterfly(X1(2,:)); butterfly(X1(3,:)); butterfly(X1(4,:))];

% 第三级蝶形运算
X3 = [butterfly(X2(1,:)); butterfly(X2(2,:))];

% 输出频谱
X = [X3(1,:) X3(2,:)];

其中，`x` 为输入信号，`W8` 为旋转因子，`butterfly` 定义了蝶形运算，`X1`、`X2` 和 `X3` 分别表示第一级、第二级和第三级蝶形运算的结果，最终输出 `X` 为 8 点 FFT 的频谱。