ZYNQ实验--CIC插值滤波器实验

一、CIC滤波器介绍

  CIC (Cascaded Integrator-Comb) 滤波器是一种常用的数字信号处理滤波器，主要用于降采样（decimation）和升采样（interpolation）操作。它具有简单的硬件实现、高效的运算速度以及适用于需要快速处理的应用场景等优点。半带滤波器具有较好的通带平坦特性,但考虑到面积和功耗,通常只能实现偶数倍的采样率变化。CIC滤波器在信号带宽为8分之一采样率时的滤波效果较好，因此在实际使用中为保证整体滤波的带内平坦和较好的滤波效果，CIC滤波会结合半带滤波器实现高倍采样率变化。

  在理论上，多级CIC（Cascaded Integrator-Comb）滤波器可以通过级联多个单级CIC滤波器来实现。然而，根据Noble恒等式的观点，即“先进行抽取或插值，再进行线性滤波”与“先进行线性滤波，再进行抽取或插值”是等价的，我们可以将多个单级CIC滤波器的结构进行变换，可以节省资源并提高计算速度，更为实用。

1.1 CIC滤波器基本结构

  CIC滤波器分为3部分：积分器，抽取/插值器，梳状器。本实验在Vivado 2018.3 环境中完成并进行仿真。根据插值滤波器的结构，使用Verilog HDL分别完成积分模块、插值模块、梳状模块的设计。下图分别为CIC抽取和CIC插值滤波器的基本结构图，图中参数：R为抽取或插值系数，M为差分延迟，在工程中一般取1或者2，N为滤波器级数。

1.2 CIC滤波器位宽确定

  上文图中R抽取或插值系数，M差分延迟(M为CIC阶数)，N滤波器级数和Bin输入数据的位宽将决定积分器Bout不发生溢出的位宽值。 计算公式如下：

1. 抽取滤波器
   $$B_{\mathrm{out}}=\lceil N{\log }_{2}RM+B_{in}\rceil$$
2. 插值滤波器
   $$B_{\mathrm{out}}=\lceil {\log }_2\frac{(RM{)}^N}{R}+B_{in}\rceil$$
   与CIC抽取滤波器不同，插值并不要求所有滤波器有相同的位宽，这里使用的最大位宽公式。

说明： M这个参数有的文章认为是CIC滤波器的阶数，但有的文章计算公式中有M却并为说明M是什么参数并且设计中也没有体现，结合找到的参考资料在本设计中认为M为差分延迟参数，设置为1。CIC抽取滤波可以在中间级进行舍入操作，但是CIC插值滤波不能进行舍入操作，因为积分器在梳状器之后，梳状器引入的量化误差会在积分器中累加引起滤波器的不稳定。

二、CIC插值滤波器实现

2.1位宽确定

  实验设计的插值滤波器插值系数 $R=96，M=1，N=4，Bin=16bit$，根据计算公式得
$$B_{\mathrm{out}}=\lceil {\log }_2\frac{(RM{)}^N}{R}+B_{in}\rceil =\lceil {\log }_2\frac{(96\ast 1{)}^4}{96}+16\rceil =36bit$$
MATLAB的FDATOOL工具观察CIC滤波器的频谱。

//模块定义
module CIC#(
    parameter STAGES = 4, // 滤波器阶数
    parameter DATA_WIDTH = 16, // 数据宽度
    parameter INDATA_WIDTH = 36, // 中间数据宽度
    parameter Ntimer = 96 // 插值倍数
)
(
    input clk_in,  // 输入数据时钟
    input clk_out, // 输出数据时钟（Ntimer倍于输入数据时钟）
    input reset,   // 复位信号
    input signed [DATA_WIDTH-1:0] data_in, // 输入数据
    output reg signed [DATA_WIDTH-1:0] data_out // 输出数据
);
// 积分器的寄存器
reg signed [INDATA_WIDTH-1:0] integrator [0:STAGES-1]; 
// 梳状器的寄存器
reg signed [INDATA_WIDTH-1:0] comb [0:STAGES-1]; 
reg signed [INDATA_WIDTH-1:0] combd [0:STAGES-1]; 
// 插值的寄存器
reg signed [INDATA_WIDTH-1:0] interpolation = 0;
reg [7:0] cont;
// 输出缓冲
reg signed [INDATA_WIDTH-1:0] output_buffer = 0;

// 将输出缓冲的值映射到输出端口
always @(posedge clk_out) begin
    data_out <= output_buffer[INDATA_WIDTH-1:INDATA_WIDTH-15]; // 可能需要调整以适应实际的位宽和动态范围
end

endmodule

2.2 梳状器模块(fs)

// 梳状器（由输入时钟驱动）
always @(posedge clk_in or posedge reset) begin
    if (reset) begin
        for (i = 0; i <STAGES; i = i + 1) begin
            comb[i] <= 0;
            combd[i]<= 0;
        end
    end 
    else begin
    	//梳状器操作
        comb[0] <= {{(INDATA_WIDTH-16){data_in[15]}},data_in};
        for (i = 0; i <STAGES; i = i + 1) begin
            combd[i] <= comb[i];
        end
        for (i = 1; i <STAGES; i = i + 1) begin
            comb[i] <= comb[i-1]-combd[i-1];
        end
    end
end

一级comb的RTL图，实现输入与输入延时一个单位后相减的操作。

2.3 插值器模块(N*fs)

对梳状器结果进行插值操作，即间隔N点补零

// 插值器(输出时钟驱动)
always @(posedge clk_out or posedge reset) begin
    if (reset) begin
          interpolation <= 0;
          cont<=0;
    end 
    else begin
         cont<=cont+1;
         if(cont==Ntimer-1) begin //N倍插值
            interpolation <=comb[STAGES-1];
            cont<=0;
         end
         else
            interpolation <=0;
    end
end

2.4 积分器模块(N*fs)

// 积分器逻辑（由输出时钟驱动）
always @(posedge clk_out or posedge reset) begin
    if (reset) begin
        for (i = 0; i <STAGES; i = i + 1) begin
            integrator[i] <= 0;
        end
        output_buffer <= 0;
    end 
    else begin
        integrator[0] <= interpolation;
        for (i = 1; i <STAGES; i = i + 1) begin //积分器
            integrator[i] <= integrator[i] + integrator[i-1];
        end
        output_buffer<=integrator[STAGES-1];
    end
end

一级integartor的RTL图

三、CIC插值滤波器实现

实验仅进行行为仿真验证

3.1 testbench

module testCIC(
    );
reg clk,reset;
reg sclk;
wire [7:0]RA;
wire [15:0]RD;
wire [15:0]RO;

initial begin
    clk  = 1'b0;
    sclk  = 1'b0;
    #100 reset=1'b1;
    #500 reset=1'b0;
    end
    
always begin
   #1 clk =~clk;
end

always begin
   #96 sclk =~sclk;
end

//正弦信号输出
count count1(   
    .CLK(sclk),
    .RST(reset),
    .RA(RA));
DDS DDS1(
    .CLK(sclk),           
    .RA(RA),       
    .RDQ(RD));
//CIC插值滤波
CIC cic1(
     .clk_in(sclk),  
     .clk_out(clk), 
     .reset(reset),  
     .data_in(RD), 
     .data_out(RO) 
);
endmodule

3.2 仿真分析

上变频时钟周期为2ns，输入时钟为192ns，图中data_in和data_out周期比例为96，符合预期设计

CIC插值滤波后的波形未出现波形失真，但存在幅度减小的情况(这里暂未弄清楚是滤波增益造成还是输入数据的问题)

中间梳妆器和积分器的过程就不再展示，可自行进行实验观测。

参考文章

Vivado CIC IP核手册
数字下变频和抽取滤波器—陈老湿·通信MATLAB仿真
《Xilinx FPGA 数字信号处理系统设计指南》-- 何宾 张艳辉 编著

本文目的主要为学习记录，如果您在本文中发现了任何错误或疑问，请随时反馈会及时做出相应的修改。