在GNU无线电中实现自定义信号源

Question

我正在尝试在GNU无线电中实现自定义信号源，该电台必须发出带有频带的白色高斯噪声，围绕频率f，具有双面带宽bw ，如图所示，图片

为此，我应该实现带通滤波器并将其应用于高斯随机变量。滤波器的输出读取

龙头是滤波器系数的位置，并且 是过滤器输入。在我的情况下，TAPS是通过使用band-pass滤波器的TAPS gnu无线电块获得的，并且将输入设置为随机高斯变量。

为了获取过滤器的输出，我应该编写工作函数看起来像

def work(self, input_items, output_items):
    samp_rate = self.get_samp_rate()
    f0 = self.get_f0()
    BW = self.get_BW()
    my_BP_filter = filter.fir_filter_fff(
                   1,
                   firdes.band_pass(
                       1,
                       samp_rate,
                       f0 - 0.5*BW,
                       f0 + 0.5*BW,
                       BW*1e-2,
                       firdes.WIN_HAMMING,
                       6.76
                       )
                   )
    taps = my_BP_filter.taps()
    ntaps = len(taps)

    min_array = min(output_items[0], taps, key=len)
    max_array = max(output_items[0], taps, key=len)
    min_len, max_len = len(min_array), len(max_array)

    noise = np.random.normal(loc=0.0, scale=1.0, size=min_len)

    t = 0.0
    for idx_out in range(len(output_items[0])):
        for idx_taps, val_taps in enumerate(taps):
            t += val_taps * noise[len(noise)-1-idx_taps]
        output_items[0][idx_out] = t

    return len(output_items[0])

，但代码非常慢。我如何具有与显示图片相似的输出（通过以下结构获得）？

Answer 1

您肯定（Capital D）不想在Work中反复重新恢复A（新的）过滤器。

请记住，在连续的大量数据上，工作被称为数千次。

因此，那是一个根本的建筑错误：您需要一个过滤器，它在工作呼叫之间保持其状态。

另外，您绝不能编写取决于数据的“视图”的信号处理，即无尽的样本流如何恰好到达您的块。想象一下这个很长的数字流：

… … -2 -1 0 1 2 3 4 5 … …

您的工作被列出了，例如，您的第一个工作将值-4000到-2000，第二个，第二个至-1990，第三个-1989至2000年，2001年至2466年，您会明白这一想法。

或者，它完全不同，并且在第一个呼叫工作时，您会在第二-2000至-1500的第二个-1499至-500上看到-4000至-2001第四-499至0，在第五1至2466中。

然后，您的min_aray，max_len取决于信号的打包方式。 这很糟糕，因为它改变了您的数学。

我会建议实施什么，但是整个块对我来说确实没有意义。这是不必要的 - 在噪声源上使用现有的带通滤波器也是如此 - 如果我推荐如何做到这一点，我建议您自己实施相同的东西，没有什么不同：两个块，一个产生噪声，使用现有过滤器的源代码为例，使用正确实现的带通滤波器将其构成第二个块。