MATLAB快速傅里叶变换（FFT）示例：不同采样点数下的频谱分析

% 定义输入序列x1和采样点数N1、N2、N3 x1 = [1 1 1 1]; N1 = 8; N2 = 16; N3 = 32;

% 对输入序列进行FFT计算，得到频域表示 y1 = (fft(x1, N1)); y2 = (fft(x1, N2)); y3 = (fft(x1, N3));

% 绘制三幅子图，分别展示不同采样点数下的频域表示 subplot(3, 1, 1); % 第一幅子图在3x1的网格中的位置 stem(0:N1-1, abs(y1)); % 绘制幅度谱，横坐标为0到N1-1 xlabel('x'); % 设置横坐标标签 ylabel('y1'); % 设置纵坐标标签 title('x1的8点DFT的结果'); % 设置子图标题

subplot(3, 1, 2); % 第二幅子图在3x1的网格中的位置 stem(0:N2-1, abs(y2)); % 绘制幅度谱，横坐标为0到N2-1 xlabel('x'); % 设置横坐标标签 ylabel('y2'); % 设置纵坐标标签 title('x1的16点DFT的结果'); % 设置子图标题

subplot(3, 1, 3); % 第三幅子图在3x1的网格中的位置 stem(0:N3-1, abs(y3)); % 绘制幅度谱，横坐标为0到N3-1 xlabel('x'); % 设置横坐标标签 ylabel('y3'); % 设置纵坐标标签 title('x1的32点DFT的结果'); % 设置子图标题