基于GNU Radio的解调后信号自动频率测量方案问询

GNU Radio原生支持解调后信号的自动频率测量功能，无需额外依赖第三方工具，常见实现方案如下：

• 方案1：直接调用原生Frequency Detector模块
  该模块为GNU Radio内置信号处理模块，支持输入解调后的实/复信号，直接输出归一化频率值，后续可对接QT GUI Number Sink模块实时展示数值结果。需按当前链路采样率做单位换算得到实际频率：

  实际频率值 = 模块输出归一化频率 × 采样率 / 2

• 方案2：FFT+峰值检测组合方案（精度更高）
  适用于对测量精度要求更高，或需要同步观测频谱分布的场景，链路搭建逻辑如下：

  解调后信号 → 流式FFT模块 → 复数转幅度模块 → 峰值检测模块 → 频率换算模块 → 数值显示模块
  峰值检测得到峰值对应的FFT bin索引后，可通过自定义Python Block完成频率换算，示例代码如下：

  import numpy as np
  from gnuradio import gr
  
  class freq_calculator(gr.sync_block):
      def __init__(self, samp_rate=1e6, fft_size=1024):
          gr.sync_block.__init__(
              self,
              name="freq_calculator",
              in_sig=[np.float32],
              out_sig=[np.float32]
          )
          self.samp_rate = samp_rate
          self.fft_size = fft_size
  
      def work(self, input_items, output_items):
          peak_bin = np.argmax(input_items[0])
          actual_freq = (peak_bin - self.fft_size/2) * self.samp_rate / self.fft_size
          output_items[0][:] = actual_freq
          return len(output_items[0])
  

  换算公式为：实际频率 =（峰值bin索引 - FFT长度/2）× 采样率 / FFT长度

• 方案3：过零检测方案（延迟最低）
  若解调后输出为单音正弦信号，可选用该方案实现最低延迟的频率测量，调用内置Zero Crossing Counter模块统计单位时间内的信号过零次数，频率值直接为过零次数/2，无需复杂换算。

• 注意事项：所有测量方案执行前，建议先对解调后信号做低通滤波，滤除带外噪声，避免噪声干扰检测结果，提升测量准确率。

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