FPGA：卷积编码及维特比译码仿真

FPGA：卷积编码及维特比译码仿真

本篇记录一下在FPGA中完成卷积编码和维特比译码的过程，通过代码解释编码的过程和译码的过程，便于理解，同时也方便移植到其他工程中。

1. 准备工作

• 卷积编译码IP核—convolutionIP核和viterbiIP核
• 卷积编码原理知识—网上有很多关于原理性的解释
• 卷积编码对照仿真过程—可以对照之前的MATLAB仿真程序便于更好的理解编译码的过程，但是有一点不同，在FPGA仿真代码中利用自然数编码，最后译码回自然数，方便对照正确性。

2. 编码过程

首先添加convolutionIP核，可以看到首先有如下的需要配置的参数。

按照上图的设置方式为2，1，7形式的卷积编码器，这个没有什么需要注意的，按照常规设置就可以了。

然后添加viterbiIP核，打开之后可以看到如下的设置参数。

在第一个Viterbi Type栏中可以可以选择Standard，然后约束长度和编码过程的设置一致，回溯深度可以利用公式计算5*(约束长度-1),这个是最小的设置的值，回溯深度至少要大于这个值。

在第二个配置页中，按照如下的参数进行设置，这个可以参考MATLAB仿真。

这个页面配置的传入的bit数据，因为编码设置的是2,1,7所以这个设置为2，

3. 仿真过程

整个仿真过程利用0-15的自然数循环，对其进行卷积编码，然后通过viterbi译码还原出原来的自然数。

程序的设置流程思路：

• 为了方便移植到其他的工程中，同时为了转换数据速率，在自然数信源和卷积编码中间添加了一个fifo
• 同样的操作，在译码结束之后添加了一个fifo,与前一个fifo对称。
• 在利用译码的过程中，有一点注意事项：编码之后的2bit数据输入到viterbi译码的IP核中时，需要在第0位和第8位填充，构成两bit。这个技术手册中有说明。

下面给出两个fifo的参数设置，自然数位宽为[0:3],首先是信源和编码之间的fifo。

然后是译码之后的fifo。

4. 完整代码

`timescale 1ns / 1ps

module conv_encoder(
    input clk,          //时钟
    input rst_n         // 复位  高电平复位
    );

parameter K = 1;   //  对应MATLAB仿真中的k和n的值，这个在IP核设置中已经有体现
parameter N = 2;   //
parameter L = 7;  // 编码之后的数据长度

reg [5:0] datain_num; // 每一组编码的原始数据个数
reg [3:0] datain;
//reg [5:0] dataout_num;  //输出编码数据的个数

// 定义viterbi IP核需要用到的信号
wire vit_datain_valid;   // 当vit ip的ready信号有效同时 conv的输出有效的时候，这个信号有效
wire vit_datain_ready;   //
wire vit_dataout_valid;
wire vit_dataout_ready;
wire vit_data_out;  // 译码结果输出

// 信源处fifo的相关信号线，该fifo是用来把多位数转为比特流传入到卷积编码中。
wire fifo_encode_empty;
wire fifo_encode_full;
wire fifo_encode_out;
wire fifo_decode_empty;
wire fifo_decode_full;
wire [3:0] fifo_decode_out;

wire rd_en; //第一个fifo的读写控制信号
wire wr_en;
reg wren;

always@(posedge clk)begin
    if(~rst_n)begin
        wren <= 1'b0;
    end
    else begin
        if(fifo_encode_full==1'b1)begin
            wren <= 1'b0;
        end
        else begin
            wren <= 1'b1;
        end
    end
end
assign wr_en = wren & (!fifo_encode_full);  // fifo没有满就往fifo中写数据 
wire rd_en2; // 第二个fifo的读写控制信号
wire wr_en2;


assign rd_en2 = !fifo_decode_empty;  // 最后一个fifo 非空就可以读数据
assign wr_en2 = (!fifo_decode_full) & vit_dataout_valid;
assign vit_dataout_ready = 1'b1;  // 最后直接进fifo了，这里不做特殊控制了，直接常为1，就可以仅看vit_dataout_valid信号了

// 编码信号的控制型号 valid 和 ready
wire conv_datain_ready;
reg conv_datain_valid;
wire [1:0] conv_dataout;
wire conv_dataout_valid;
wire conv_dataout_ready;


// 设计输入数据 这里面需要对输入的数据转化弄成位的形式，
// 这个和卷积编码的参数设置有关，2,1,7,接受一个输入bit生成两个bit

always@(posedge clk)begin
    if(~rst_n)begin
        datain <= 4'b0;
        datain_num <= 6'b0;
    end
    else begin
        if(wr_en == 1'b1)begin // 数据只在wr_en有效的情况下才逐渐累加，这个是为了通过连续的数值检验译码正确性
            datain <= datain + 4'b1;
            datain_num <= datain_num + 6'b1;
        end
        else begin
            datain <= datain;
            datain_num <= datain_num;
        end
    end
end

// 在这里添加转换bit 可以用fifo实现，同时能够控制速率 

assign rd_en = conv_datain_ready;
fifo_encode u1 (
  .clk(clk),      // input wire clk
  .srst(~rst_n),    // input wire srst
  .din(datain),      // input wire [3 : 0] din
  .wr_en(wr_en),  // input wire wr_en
  .rd_en(rd_en),  // input wire rd_en
  .dout(fifo_encode_out),    // output wire [0 : 0] dout
  .full(fifo_encode_full),    // output wire full
  .empty(fifo_encode_empty)
);

always@(posedge clk)begin
    if(~rst_n)begin
        conv_datain_valid <= 1'b0;
    end
    else begin
        conv_datain_valid <= rd_en;
    end
end

assign conv_dataout_ready = 1'b1;

convolution_0 conv (
  .aclk(clk),                              // input wire aclk
  .aresetn(rst_n),                        // input wire aresetn
  .s_axis_data_tdata(fifo_encode_out),    // input wire [7 : 0] s_axis_data_tdata  注意这个只能接受1bit有效数据
  .s_axis_data_tvalid(conv_datain_valid),  // input wire s_axis_data_tvalid
  .s_axis_data_tready(conv_datain_ready),  // output wire s_axis_data_tready
  .m_axis_data_tdata(conv_dataout),    // output wire [7 : 0] m_axis_data_tdata
  .m_axis_data_tvalid(conv_dataout_valid),  // output wire m_axis_data_tvalid
  .m_axis_data_tready(conv_dataout_ready)  // input wire m_axis_data_tready
);

assign vit_datain_valid = conv_dataout_valid & conv_dataout_ready;
    
// 编码完成后进行译码过程
viterbi_0 viterbi (
  .aclk(clk),                                // input wire aclk
  .aresetn(rst_n),                          // input wire aresetn
  .s_axis_data_tdata({7'b0,conv_dataout[1],7'b0,conv_dataout[0]}),      // input wire [15 : 0] s_axis_data_tdata 这个经过编码之后的数据需要按照此种方式输入。
  .s_axis_data_tvalid(vit_datain_valid),    // input wire s_axis_data_tvalid
  .s_axis_data_tready(vit_datain_ready),    // output wire s_axis_data_tready
  .m_axis_data_tdata(vit_data_out),      // output wire [7 : 0] m_axis_data_tdata
  .m_axis_data_tvalid(vit_dataout_valid),    // output wire m_axis_data_tvalid
  .m_axis_data_tready(vit_dataout_ready)   // input wire m_axis_data_tready
);


// 添加一个把数据恢复成自然数的fifo
fifo_decode u2 (
  .clk(clk),      // input wire clk
  .srst(~rst_n),    // input wire srst
  .din(vit_data_out),      // input wire [3 : 0] din
  .wr_en(wr_en2),  // input wire wr_en
  .rd_en(rd_en2),  // input wire rd_en
  .dout(fifo_decode_out),    // output wire [0 : 0] dout
  .full(fifo_decode_full),    // output wire full
  .empty(fifo_decode_empty)  // output wire empty
);

endmodule

下面为tb文件:

`timescale 1ns / 1ps

module conv_tb();

reg l_clk;
reg rst_n;

conv_encoder conv_test_ins(
    .clk(l_clk),          //时钟
    .rst_n(rst_n)         // 复位  高电平复位
//    input [7:0] data_in,  // 输入的待编码数据
//    output [7:0] dataout      // 输出的解码数据
    );

initial l_clk = 1;
always #5 l_clk= !l_clk;  //15.625   


initial begin
    rst_n <= 0;
    #40;
    rst_n <= 1;
    #320;
    //#50000000;
    #320;
//    $stop;
end
endmodule

5. 结果分析

以下为仿真结果图：

最下面的红色线是译码之后经过fifo速率转换之后的结果，可以看到是从0依次递增的自然数，一直到15，然后循环下去。这个有效是和rd_en2这个信号保持一致的，但是看这组红线，译码的结果是两个0，这是因为fifo读数据的时候是有一个clk的时钟延时的，这也是为什么在给出fifo参数设置的时候把latency=1用红线框起来。所以把rd_en2延迟一个系统时钟对照fifo_decode_out看数据就正确了。在使用的时候看接下来数据处理的过程需求，采用合适的操作，这里只做简单的仿真验证。