用Xilinx Vivado HLS可以快速、高效地实现QRD矩阵分解

用XilinxVivadoHLS可以快速、高效地实现QRD矩阵分解在数字信号处理领域，如自适应滤波、DPD系数计算、MIMODecoder等，常常需要矩阵解方程运算以获得其系数，因此需对矩阵进行求逆运算。然而，由于直接对矩阵求逆会导致庞大的运算量，所以在实际工程中往往需要先将矩阵分解成几个特殊矩阵（正规正交矩阵或上、下三角矩阵以求其逆矩阵需要更小的运算量）的乘积。目前，QRD矩阵分解法是求一般矩阵全部特征值的最有效且广泛应用的方法之一。它是将矩阵分解成一个正规正交矩阵Q与上三角形矩阵R，称为QRD矩阵分解。由于浮点具有更大的数据动态范围，所以在众多多算法中具有只需要一种数据类型的优势，所以很多QRD矩阵分解是基于浮点数据类型的。不过在通信应用中，更多的场景还是复数类型。因此，随着通信技术的日益发展，算法的复杂度越来越高，QRD矩阵的维度也越来越大。如果是用传统的手写RTL，浮点复数超大维度QRD矩阵分解的FPGA实现将变得非常复杂，需要很长的时间来编写RTL代码、仿真和进行验证等工作，使得开发效率不是很高。本文将介绍如何使用XilinxVivadoHLS（Vivado高层次综合）工具实现浮点复数QRD矩阵分解并提升开发效率。使用VivadoHLS可以快速、高效地基于FPGA实现各种矩阵分解算法，降低开发者对算法FPGA的实现难度。其中包括：

使用VivadoHLS开发效率比手写RTL实现快5-10倍，而实现的FPGA资源效率与手写RTL接近

由于C/C++仿真验证比传统FPGARTL要快100倍，VivadoHLS实现可以大大缩短用户的代码开发时间和仿真验证时间，从而大幅提高生产效率。1.VivadoHLS工具介绍VivadoHLS是Xilinx公司2010年收购AutoESL以后重新打造的高层次综合工具，它可以让用户通过添加适当的directives（制导语句）和constrains（约束）,将其C/C++/SystemC代码直接转换成FPGARTL（Verilog,VHDL,SystemC）代码。让用户可以在算法开发环境而非通常的硬件开发环境中只需专注于算法规格和算法的C实现，VivadoHLS工具会自动考虑FPGA微观实现架构，并可生成可综合的FPGARTL代码。如图1所示。

图1VivadoHLS设计介绍VivadoHLSFPGA设计流程：

首先用C/C++/SystemC将算法实现，并编写Ctestbench验证C的功能，确保其功能正确；

然后就可以通过VivadoHLS工具进行C综合,将C转换成RTL；

接着做C/RTL的协同仿真（Co-simulation），以保证生成的RTL代码功能与C的功能完全一致。

最后VivadoHLS生成的RTL代码可直接用于Xilinx设计开发环境下做系统集成、仿真和生成bit文件。如图2所示。

图2VivadoHLS设计流程2.QRD矩阵分解算法

3.QRD矩阵分解VivadoHLS实现C++代码构架QRD矩阵分解C++实现代码的顶层模块是qrd_engine.cpp，它调用cal_core.cpp（核心计算函数）,coef_cal.cpp（系数计算函数），以及浮点加、减、乘、除法等子函数来实现。structcf_t{

floatre;

floatim;

};

voidqrd_engine

(

cf_tin_u[(R_DIM+Y_DIM)/DIV_NUM][DIV_NUM],

cf_tpd_err_in,

floatlamda,

floatlamda_sqrt,

floatdiag[R_DIM],

cf_tr[R_DIM][X_DIM],

cf_tp[R_DIM]

)

{

#pragmaHLSARRAY_PARTITIONvariable=in_ucompletedim=2

#pragmaHLSARRAY_PARTITIONvariable=rcompletedim=2//注：这里对数组加入完全分割directive，目的是提高数据的并行带宽，从而获得并行