FFT函数如何实现时域信号到频域谱图的转换及频谱搬移操作？

本文共计884个文字，预计阅读时间需要4分钟。

如果您需要对时域信号进行频域分析，但无法通过FFT函数准确获得清晰的频谱结构，可能是因为以下原因：

一、正确构造输入信号与采样参数

fft函数依赖于信号的离散采样点数和采样频率，若采样率过低或信号长度非2的整数次幂，会导致频谱分辨率下降或泄漏。需确保采样满足奈奎斯特准则，并合理设定N点长度。

1、定义采样频率 fs = 1000（单位：Hz），保证大于信号最高频率的两倍。

2、设定信号持续时间 T = 1（单位：秒），生成时间向量 t = linspace(0, T, fs*T)。

3、构建合成信号，例如包含50 Hz与120 Hz正弦分量的叠加：x = sin(2*pi*50*t) + 0.5*sin(2*pi*120*t)。

4、确定FFT点数 N = 1024，若原始长度不足则补零，超过则截断。

二、调用fft函数并计算幅值谱

直接调用fft函数返回复数序列，其模长反映各频率分量强度，但默认输出为零频在首、正负频连续排列的顺序，不便于直观观察。

1、执行变换：X = fft(x, N)。

2、计算单边幅值谱：P2 = abs(X/N)。

3、提取单边谱：P1 = P2(1:N/2+1)，并将除直流与奈奎斯特点外的所有点乘以2：P1(2:end-1) = 2*P1(2:end-1)。

三、生成对应频率轴并绘制原始频谱

频率轴需与FFT输出点一一对应，否则横坐标失真。未搬移时，零频位于左侧，正频居中，负频在右，不符合人眼习惯。

1、构造频率向量：f = fs*(0:(N/2))/N。

2、使用plot(f, P1)绘制单边幅值谱图。

3、添加坐标轴标签：xlabel('Frequency (Hz)')，ylabel('|P1(f)|')。

四、执行fftshift实现频谱搬移

fftshift将零频分量平移到频谱中心，使负频左、零频中、正频右，符合标准频谱表示方式，尤其适用于含载波信号或带通信号分析。

1、对原始FFT结果应用搬移：X_shifted = fftshift(X)。

2、计算搬移后幅值：P_shifted = abs(X_shifted)/N。

3、构造搬移后频率轴：f_shifted = fs*(-N/2:N/2-1)/N。

4、绘制搬移谱：plot(f_shifted, P_shifted)。

五、手动实现频谱搬移（绕过fftshift）

当环境不支持fftshift函数（如部分嵌入式MATLAB精简版或自定义FFT实现）时，可通过索引重排模拟等效效果，增强代码可移植性。

1、分离原始FFT输出的前后两半：X_first = X(1:N/2)，X_second = X(N/2+1:end)。

2、拼接为新序列：X_manual = [X_second, X_first]。

3、计算幅值：P_manual = abs(X_manual)/N。

4、生成对应频率轴：f_manual = [fs/2:fs/N:fs-1, -fs/2:fs/N:-fs/N]（注意浮点精度处理）。