[发明专利]一种基于FPGA的Jacobi变换的高效实现方法

说明书

技术领域

本发明属于信号处理技术领域，具体的说是涉及一种基于FPGA的Jacobi变换的高效实现方法。

背景技术

在信号处理中，Jacobi变换是一个应用广泛的矩阵变换，可以用求解矩阵的奇异值分解SVD、求逆等，而这些基本的矩阵分解算法在科学计算、信号处理领域具有广泛的应用，如数据压缩、噪声去除、数值分析，包括近几年兴起的机器学习、深度学习其基本核心操作也包括矩阵奇异值分解与求逆等变换。实现这些矩阵分解算法常用方法有Gauss变换、Householder变换、Jacobi变换等，其中Jacobi变换是精度较高的方法并且很适合在FPGA中实现。因此基于FPGA的Jacobi变换的高效实现技术在实际工程中具有很高的应用价值。

目前工程上多数采用如附图1所示的并行结构来计算Jacobi变换矩阵中的这四个元素值，其中G为Jacobi旋转矩阵。虽然这样的结构可以使用流水线方式计算但在实际计算过程中，每次变换都会影响矩阵A中第i行、第j行、第i列、第j列的元素，所以需要等待这些元素计算完成才能进行下一次Jacobi变换。并且传统的并行计算方案需要多个CORDIC模块，从输入数据有效到输出数据有效需要进行两次CORDIC计算。因为需要等待其他元素计算完成并不能有效使用这些CORDIC模块的计算能力，所以这样的方案消耗较多的硬件资源，反而却不能很好的利用起并行结构可以进行流水线方式高速计算的优点。同时因为FPGA资源消耗较多，在一些低端的FPGA芯片中，不能实现基于传统并行计算方案实现Jacobi旋转的信号处理算法。

发明内容

本发明所要解决的技术问题是，针对上述问题，提供一种FPGA资源消耗比较小的Jacobi变换技术，可以在资源量比较小的FPGA芯片中实现基于Jacobi变换的信号处理算法，如单基线相位干涉仪等。与传统计算方案相比在完成计算需要时钟数基本不增加的前提下减少了CORDIC模块，减少了FPGA逻辑资源的消耗。

本发明的技术方案是：

一种基于FPGA的Jacobi变换的高效实现方法，其特征在于，包括以下步骤：

a.确定Jacobi变换前后矩阵A中元素之间的关系：

设矩阵A∈R^n×n，单次Jacobi旋转的计算公式如下公式1所示：

其中，G∈R^n×n为Jacobi旋转矩阵，G的对角线元素除了g_ii＝cosθ、g_jj＝cosθ其他全为1，其非对角线元素除了g_ij＝-sinθ、g_ji＝sinθ外其他全为0：

则完成一次Jacobi变换之后的与A的a_ii、a_jj、a_ij、a_ji的关系确定为如下公式2所示：

b.构建用于Jacobi变换的Cordic模块：

构建一个流水线的旋转Cordic模块，输入为(xi,yi,θ)，输出为(xo,yo)，输入与输出关系如下公式3所示：

c.采用步骤b中构建的Cordic模块计算Jacobi变换后矩阵A的元素

c1.将Jacobi变换前矩阵A的元素(a_ji,a_ii,θ)、(a_jj,a_ij,θ)在两个时钟内分别依次输入到Cordic模块中，即依次计算公式4和公式5：

获得中间变量：(x₁,y₁,x₂,y₂)；

c2.将步骤c1中获得的中间变量组合为(y₂,y₁,θ)、(x₂,x₁,θ)，并在两个时钟内分别依次输入到Cordic模块中，即依次计算公式6和公式7

得到Jacobi变换后矩阵A的元素完成Jacobi变换。

本发明的有益效果为：与传统方法相比，本发明的方法使用一种高效的串行控制方案算法，只使用一个Cordic模块实现Jacobi变换的计算，并且在减少硬件资源消耗的同时，完成计算所需要的时间与并行算法一样，即这样的计算方案算法具有计算速度快、节省硬件资源消耗的优点，在实际工程中具有很高的应用价值。

附图说明