频谱分析仪FFT无法正常绘图问题求助

核心问题定位 #

你提到FFT计算结果有真实数值，但绘图显示零值，问题大概率出在FFT结果的后处理或绘图参数设置上，而非FFT计算本身。以下是针对性排查步骤：

1. 检查FFT幅度的处理逻辑 #

FFT输出是复数数组，直接绘图会导致异常，必须先提取幅度并做正确缩放：

• 提取幅度：fft_magnitude = np.abs(fft_result)
• 实信号频谱仅需保留前半部分（避免镜像）：fft_magnitude = fft_magnitude[:N//2]（N为采样点数）
• 幅度归一化：对于实信号，非直流分量需除以N/2，直流分量除以N，确保幅度与原始时域信号匹配：

  fft_magnitude[0] = fft_magnitude[0] / N  # 直流分量
  fft_magnitude[1:] = fft_magnitude[1:] / (N/2)  # 交流分量
  

2. 验证频率轴的正确性 #

频率轴范围不匹配会导致峰值被压缩至不可见：

• 正确计算频率轴：基于采样率Fs，频率范围应为0到Fs/2，步长为Fs/N：

  freq_axis = np.linspace(0, Fs/2, N//2)
  

• 若频率轴设置过大（如远高于信号实际频率），会导致峰值窄到看起来像零，需确认Fs的赋值与实际采样率一致。

3. 排查绘图参数设置 #

• 避免直接绘制复数数组或未处理的FFT结果，必须使用处理后的幅度数组与对应频率轴配对绘图：

  plt.plot(freq_axis, fft_magnitude)
  

• 检查坐标轴刻度：若使用对数纵坐标（plt.yscale('log')），但幅度值过小（如小于1e-6），会显示为零，可切换为线性刻度或调整纵坐标范围：

  plt.ylim(0, np.max(fft_magnitude)*1.2)  # 锁定纵坐标范围到峰值的1.2倍
  

• 确认绘图数据维度一致：fft_magnitude与freq_axis的长度必须相同，否则会导致绘图异常。

4. 快速验证方法 #

打印处理后的fft_magnitude最大值和平均值，若数值远小于坐标轴默认范围（比如最大值为0.1，但坐标轴默认到10），只需手动设置纵坐标范围即可显示出信号。

内容的提问来源于stack exchange，提问作者Skymero