傅里叶变换FFT时域转频域详解

本文深入解析了傅里叶变换（FFT）如何将时域信号精准转换为频域表示，并重点攻克频谱显示不清晰这一常见痛点——从采样参数设置、零填充策略到频谱搬移的完整闭环，手把手演示如何通过奈奎斯特准则设计信号、调用fft计算幅值谱、构建正确频率轴、使用fftshift居中零频，甚至提供不依赖内置函数的手动索引重排方案，兼顾原理透彻性与工程实用性，助你真正看懂、画准、用好频谱图。

如果您对时域信号进行频域分析，但无法通过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]（注意浮点精度处理）。

理论要掌握，实操不能落！以上关于《傅里叶变换FFT时域转频域详解》的详细介绍，大家都掌握了吧！如果想要继续提升自己的能力，那么就来关注golang学习网公众号吧！