【SoCVista】9.滤波器算法强度缩减.docVIP

第九章、滤波器的算法强度缩减 学过数字信号处理，就不可能不知道滤波器。经典的滤波器分两类：FIR和IIR。FIR是有限长冲击响应滤波器，硬件电路是非递归的；而IIR是无限长冲击响应滤波器，硬件电路存在递归环路，值得注意的是，IIR可以看成是一个FIR和递归环路的级联。本章的目的在于如何并行实现FIR滤波器，在获得高吞吐率的同时尽可能节约面积和功耗。 回想第八章、快速卷积，前面讨论过FIR其实也是一个卷积的过程，只是这个卷积不是短卷积，而是长卷积。根据数字信号处理理论（参见——胡广书《数字信号处理：理论、算法与实现》第二版3.6.2节），长卷积可以用短卷积实现，短卷积正是第八章所解决的问题，如此一来FIR就能高效实现。OK！没错，第八章的知识已经为FIR并行实现提供了一个途径，但是这里要讨论FIR的另一种并行结构，当然了，这种新的FIR并行结构仍然和快速卷积密不可分，稍后便见分晓。 在进入正题之前，要先擦亮眼睛！记住一点，卷积其实是一种多项式乘法，快速卷积给出了“短”多项式乘法的有效实现，只要是多项式乘法，就能用快速卷积知识搞定。废话那么多，无非想告诉大家，接下来讨论的FIR并行实现其实可看成是“短”多项式乘法。 本章的任务：弄清楚以下4个问题， FIR的多相分解。 多相FIR的并行结构，以及如何使用快速卷积构造高效的多相FIR并行结构。 大尺寸并行FIR算法。 转置的并行滤波器。 很多数字信号算法，比如FFT和现在讨论的并行FIR，都让我联想到“分而治之”这个道理。分治就是把一个复杂的问题分成两个或更多的相同或相似的子问题，再把子问题分成更小的子问题……直到最后子问题可以简单的直接求解，原问题的解即子问题的解的合并。 FIR多相分解也是一种“分而治之”处理，其做法就是将“待卷积”的两个序列和按交叉数据分配的原则分为“相同数量的”若干个子序列，比如2相分解如下示， 图 1 2相分解 在频域，可表示为 其中， 同理，有 其中，当为偶数时， 若为奇数，只需修正一下公式(4)的求和上标即可。 如此滤波结果的频域表示为， 如果也做两相分解，则正好有 其中， 将替换为，则公式(7)就是 也就是说的奇偶序列分别对应两个子序列序列的和，即 图 2 2相分解并行结构，没有经过优化，占用资源多 很多朋友看到这可能会很迷惑，或者说有点乱！！其原因是什么呢？这里给出几点关键的提示，大家开动脑筋仔细的推导一遍，相信会有所收获： 为什么序列的奇偶分解，在频域会表示成如公式(1)、公式(3)和公式(6)的形式，其中公式右边的自变量是，而不是“理所当然”的？在公式右边第二项中出现的又是为何？ 以为例，“深入”推导一下， 记的奇偶序列分别为和，则公式(9)可表示如下 因为 显然，有 ，，，注意，的奇偶序列的z变换分别是和，而不是和。我们不能简单的说，是奇偶序列的z变换之和，关系没那么简单，详细思考公式(10)，看看到底意味着什么就清楚了。 为什么会有公式(7)，也就是说，怎么知道属于，而不是？ 如果仔细观察/“深入去想象”和每一项的因子的幂次到底是什么就能明白。其实，每一项因子的幂次必为偶数，而每一项因子的幂次必为奇数。的每一项因子的幂次是偶数，所以属于。大家可以摸清其中的规律，这样在更多相分解（3相或4相等等）时，就不会搞错。 公式(5)的中的意味着2个延时，可是为什么图2中只有一个D？课本上说到“2并行电路的属于快延时，在实际电路中就是一个延时”，此话怎讲？还有公式(5)的中的又表示什么呢，难道是1/2个延时吗？ 实际上，图2是根据公式(8)画出的，而非公式(5)。注意公式(5)(6)(7)(8)，公式(8)也是从公式(5)推导而来。请思考公式(5)和公式(8)的关系，见图3，顶端子图相当于公式(8)所表示的奇偶序列，往下即为公式(7)，再往下的奇序列右移一位就是公式(6)中的，最后是序列相加，结果正好是完整序列。 图 3 序列的合成 题外话：序列的多相分解初看感觉很容易，但细想又没那么简单。能深入理解每个公式的物理意义是最好不过的，如果不能，至少应该记住多相分解的公式！！ 推而广之，FIR的3相分解如下， 简化后，有 三路输出，分别为 写成矩阵形式，有 公式(15)和(16)对应于2相分解中的公式(6)和(7)，注意，公式(15)和(16)的自变量是，只需进行变量代换，将替换为，就能得到FIR 3并行电路的构造式子，如下 根据公式(17)很容易画出最终电路图，见课本图9-2。 值得注意的3点： 公式(16)右边的向量是多项式相乘的结果的各项系数，见公式(14)的第二步。 对这些系数进行“合并同类项”，就能得到输出，也就是公式(14)的最后一步，写成