FPGA 的 DSP：Verilog 中的简单 FIR 滤波器

本项目介绍如何用 Verilog 实现一个带有预生成系数的简单 FIR 滤波器。

Things used in this project 、

Story

简陋的 FIR 滤波器是 FPGA 数字信号处理中最基本的构建模块之一，因此了解如何利用给定的抽头数和相应的系数值组装一个基本模块非常重要。因此，在这个关于在 FPGA 上入门 DSP 基础知识的实用方法迷你系列中，我将从一个简单的 15 抽头低通滤波器 FIR 开始，先在 Matlab 中生成初始系数值，然后将这些数值转换为 Verilog 模块中的使用值。

有限脉冲响应或 FIR 滤波器的定义是，滤波器的脉冲响应在一定时间内趋于零值，因此它是有限的。脉冲响应归零所需的时间与滤波器的阶（抽头数）直接相关，而阶（抽头数）就是 FIR 底部传递函数多项式的阶数。FIR 的传递函数不包含反馈，因此如果输入一个值为 1 的脉冲，然后是一堆零值，那么输出将只是滤波器的系数值。

任何滤波器的作用都是对信号进行调节，主要侧重于选择滤除或允许哪些频率通过。最简单的例子之一就是低通滤波器，它允许低于某一阈值（截止频率）的频率通过，同时大大衰减高于该阈值的频率，如下图所示。

本项目的主要重点是在 HDL（特别是 Verilog，但其概念可以很容易地转换为 VHDL）中实现 FIR，可将其分解为三个主要逻辑组件：将每个采样时钟送入的循环缓冲器，该缓冲器可适当考虑串行输入的延迟；每个抽头系数值的乘法器；以及每个抽头输出相加结果的累加器寄存器。

由于我的重点是 FPGA 逻辑中的 FIR 机制，因此我只是使用 Simulink 和 Matlab 中的 FDA 工具为低通滤波器插入一些简单的参数，然后使用生成的系数值为 Verilog 模块计算适当的寄存器值（稍后完成）。

我选择实施一个简单的 15 抽头低通滤波器 FIR，采样频率为 1Ms/s，通带频率为 200 kHz，截止频率为 355kHz，从而得到以下系数：

Create Design File for FIR Module

在新的 Vivado 项目中从头开始，使用流程导航器窗口中的添加源选项为 FIR 模块创建新的设计源。

在确定了 FIR 的阶数（抽头数）和系数值之后，下一组必须定义的参数是输入采样、输出采样和系数本身的位宽。

对于这个 FIR，我选择将输入采样和系数寄存器设置为 16 位宽，输出采样寄存器设置为 32 位宽，因为两个 16 位值的乘积是一个 32 位值（两个值相乘的宽度相加得出乘积的宽度，因此如果我选择 16 位输入采样和 8 位抽头，那么输出采样的宽度将是 24 位）。

这些值也都是带符号的，因此 MSB 被用作符号位，其余较低的位数则是值必须包含的位数（在选择输入采样寄存器的初始宽度时务必牢记这一点）。要在 Verilog 中将这些值设置为有符号数据类型，需要使用关键字 signed：

reg signed [15:0] register_name;

接下来要解决的问题是如何在 Verilog 中处理系数值，需要将十进制点值转换为定点值。由于所有的系数值都小于 1，因此寄存器的全部 15 位（总共 16 位中的 MSB 为带符号位）都可以用于小数位。通常，你必须决定寄存器中的整数部分和小数部分各占多少位。因此，转换小数抽头的数学方法是：（小数系数值）*(2^(15))，如果系数值为负数，该乘积的任何十进制值都将被四舍五入，并计算该值的两位数：

现在，我们终于可以专注于 FIR 模块的逻辑了，首先是循环缓冲器，它引入串行输入采样流，并为滤波器的 15 个抽头创建一个包含 15 个输入采样的数组。

    always @ (posedge clk)
        begin
            if(enable_buff == 1'b1)
                begin
                    buff0 <= in_sample;
                    buff1 <= buff0;        
                    buff2 <= buff1;         
                    buff3 <= buff2;      
                    buff4 <= buff3;      
                    buff5 <= buff4;       
                    buff6 <= buff5;    
                    buff7 <= buff6;       
                    buff8 <= buff7;       
                    buff9 <= buff8;       
                    buff10 <= buff9;        
                    buff11 <= buff10;       
                    buff12 <= buff11;       
                    buff13 <= buff12;       
                    buff14 <= buff13;    
                end
        end

接下来，乘法阶段将每个样本乘以每个系数值：

always @ (posedge clk)
        begin
            if (enable_fir == 1'b1)
                begin
                    acc0 <= tap0 * buff0;
                    acc1 <= tap1 * buff1;
                    acc2 <= tap2 * buff2;
                    acc3 <= tap3 * buff3;
                    acc4 <= tap4 * buff4;
                    acc5 <= tap5 * buff5;
                    acc6 <= tap6 * buff6;
                    acc7 <= tap7 * buff7;
                    acc8 <= tap8 * buff8;
                    acc9 <= tap9 * buff9;
                    acc10 <= tap10 * buff10;
                    acc11 <= tap11 * buff11;
                    acc12 <= tap12 * buff12;
                    acc13 <= tap13 * buff13;
                    acc14 <= tap14 * buff14;
                end
        end 

乘法阶段的结果值通过加法累积到寄存器中，最终成为滤波器的输出数据流。

 /* Accumulate stage of FIR */   
    always @ (posedge clk) 
        begin
            if (enable_fir == 1'b1)
                begin
                    m_axis_fir_tdata <= acc0 + acc1 + acc2 + acc3 + acc4 + acc5 + acc6 + acc7 + acc8 + acc9 + acc10 + acc11 + acc12 + acc13 + acc14;
                end
        end 

最后，逻辑的最后一部分是 FIR 模块的数据流接口。AXI 流接口是最常见的接口之一，因此我选择了它来实现。其关键在于有效信号和就绪信号，这两个信号可以控制上下游设备之间的数据流。这意味着 FIR 模块需要向其下游设备提供一个有效信号，以表明其输出为有效数据，同时在下游设备取消就绪信号时，能够暂停（但仍保留）其输出。FIR 模块还必须能够在其主站接口上以同样的方式与上游设备进行通信。

以下是 FIR 模块的逻辑设计概览：

请注意有效信号和就绪信号是如何设置 FIR 输入循环缓冲器和乘法级的使能值的，以及数据或系数通过的每个寄存器是如何声明为带符号的。

FIR 模块 Verilog 代码：

`timescale 1ns / 1ps

module FIR(
    input clk,
    input reset,
    input signed [15:0] s_axis_fir_tdata, 
    input [3:0] s_axis_fir_tkeep,
    input s_axis_fir_tlast,
    input s_axis_fir_tvalid,
    input m_axis_fir_tready,
    output reg m_axis_fir_tvalid,
    output reg s_axis_fir_t