基于分时复用技术的cholesky矩阵求逆系统技术方案

本发明专利技术公开了一种基于分时复用技术的cholesky矩阵求逆系统，降低了在雷达信号处理中对共轭对称矩阵求逆所使用的硬件资源，同时通过提高矩阵求逆运算的并行度来保证运算速度。本发明专利技术根据公式A = L * D * LH，通过主控状态机将cholesky分解矩阵求逆划为三个步骤进行，分别是cholesky分解、上三角矩阵求逆和矩阵相乘。在主控状态机的控制下，通过复选器，分时复用运算资源和存储资源。该技术方案是基于单端口存储器和单精度浮点复数运算单元的并行流水处理，并支持任意阶数的共轭对称矩阵求逆，有效加快了运算速度，提高了硬件利用率，且具有良好的兼容性。

【技术实现步骤摘要】

本专利技术属于雷达信号处理
，具体是一种基于FPGA的正定复数矩阵求逆的硬件结构。在雷达信号处理中，常对正定复数矩阵求逆，若使用通用矩阵求逆的算法和结构，会使用额外的资源。本专利技术提出的硬件结构专门针对正定复数矩阵求逆，适合雷达信号处理中的求逆操作。
技术介绍
矩阵求逆的方法有很多，包括伴随矩阵法、高斯消元法、初等变换法、矩阵分解法等等。上述各种矩阵求逆的方法中，伴随矩阵法、高斯消元法、初等变换法中需要求大量的行列式，每个行列式都几乎要计算到所有的矩阵元素，计算量大，对存储空间的需求大，在FPGA上实现比较困难。矩阵分解法是将矩阵分解为一些较简单矩阵的乘积，如三角矩阵或酉矩阵，他们都具有某种特性，较容易得到逆矩阵。对分解得到的矩阵求逆后相乘，得到原矩阵的逆矩阵。矩阵分解法分为几种：1)LU分解法，2)QR分解法，3)奇异值分解法，4)cholesky分解法。前三种矩阵分解法，都是对所有可逆复数矩阵进行求解与计算，算法较复杂，资源使用较多；而第四种分解法是针对共轭对称矩阵求逆，算法较简单，资源使用也较少。而在雷达信号处理中，主要是对自相关后得到的共轭对称矩阵进行求逆操作，因此本专利技术选择cholesky分解法。cholesky分解法常用于线性方程的求解、矩阵求逆和行列式计算，在雷达信号处理、科学计算等诸多领域也有广泛应用。根据公式A＝L*D*LH，cholesky分解法是将原矩阵A分解为单位下三角矩阵L(对角线元素为1)、对角矩阵D，上三角矩阵LH(是L的共轭转置矩阵)，原矩阵A对角线上的元素必须是实数，其他元素可以是复数。因三角矩阵几乎一半数据都为0，所以求逆简单，且LH是L的共轭转置矩阵，因此只需求得(LH)-1，取共轭转置即得到了L-1。令Gij＝Lij*Dii,则cholesky分解公式可由此推得：d11＝a11，对i＝2...n有gij=aij-Σk=1j-1gik*ljk*,(j=1,2...i-1);lij=gijdij,(j=1,2...i-1);dii=dii-Σk=1i-1gik*lik*;]]>由于dii出现在分母中，因此只有当dii≠0(i＝1...n)时，分解才能进行到底。此外，在实际的计算过程中，即使不为零，但当它很小的时候，由于计算机的精度限制，会导致下溢，也会使分解运算不稳定。
技术实现思路
专利技术目的：在FPGA芯片上提供一种硬件复杂度低，资源利用率高的正定矩阵浮点求逆结构。在尽可能减少硬件资源消耗的情况下，通过提高矩阵求逆运算的并行度来保证运算速度。技术方案：一种基于分时复用技术的cholesky矩阵求逆系统，包括五个模块：1)主控状态机模块2)资源复用模块3)cholesky分解模块4)上三角矩阵求逆模块5)矩阵相乘模块。所述的主控状态机模块与系统总线相连接，在求逆操作开始前，接收使能拉高，根据待求逆矩阵维度，从外部按行顺序接收原矩阵数据，存放在资源复用模块的存储器中，接收完成后，将接收使能拉低，并通过cholesky分解模块、上三角矩阵求逆模块、矩阵相乘模块分别进行cholesky分解、上三角矩阵求逆、矩阵乘法3个步骤进行cholesky矩阵求逆操作。求逆完成后向总线发出求逆结束信号，通知系统总线从存储模块读取运算结果。所述的资源复用模块内部有3个存储器，1个单精度浮点复数乘法器，1个单精度浮点复数累加器，存储器用于存放求逆过程中各个阶段的矩阵数据，在cholesky分解，上三角求逆，矩阵乘法的运算中，都需进行复数乘法，复数累加操作，且三个运算步骤是顺序执行，因此可复用单精度浮点复数乘法器和单精度浮点复数累加器(见表1)。资源复用模块与主控状态机模块、cholesky分解模块、上三角矩阵求逆模块、矩阵相乘模块相连接。根据主控状态机的状态，资源复用内部的存储资源和运算资源受不同模块控制。所述的cholesky分解模块是通过对原矩阵元素进行乘法，加法，除法运算将原矩阵A分解为单位下三角矩阵L(对角线元素为1)、对角矩阵D，上三角矩阵LH(是L的共轭转置矩阵)。cholesky分解模块内部是一个地址生成状态机和两个单精度浮点除法器，外部与主控状态机和资源复用模块相连，在分解过程中将使用资源复用模块中的2个存储器和浮点乘法器，浮点累加器，2个存储器分别存放原始矩阵A和分解得到的上三角矩阵LH。在接收到主控状态机模块输入的cholesky分解起始信号后，cholesky分解模块的内部地址生成状态机根据矩阵维度开始按cholesky分解公式的计算顺序生成原矩阵读地址和分解后上三角矩阵的写地址。将读出的数据依次送入浮点乘法器，浮点累加器，浮点除法器，除法器结果即为分解后上三角矩阵的元素。在分解完成后，通过复用浮点除法器对分解后得到的对角矩阵的所有对角元素求倒数(即求对角阵的逆)。所述的上三角矩阵求逆模块，因choelsky分解后得到的上三角矩阵的对角线元素全为1，所以可以对上三角矩阵求逆公式进行简化，省去除法操作。上三角矩阵求逆模块内部是一个地址生成状态机，外部与主控状态机和资源复用模块相连，在上三角矩阵求逆过程中将使用资源复用模块中的2个存储器、浮点乘法器和浮点累加器，2个存储器分别存放的是分解得到的上三角矩阵LH和求逆得到的上三角逆矩阵(LH)-1，在接收到主控状态机模块输入的上三角矩阵求逆起始信号后，上三角矩阵求逆模块的内部地址生成状态机根据矩阵维度开始按上三角矩阵求逆公式的计算顺序生成原上三角矩阵读地址和上三角逆矩阵的写地址，将读出的数据依次送入浮点乘法器，浮点累加器，累加器结果即为上三角逆矩阵的元素。所述的矩阵乘法模块内部是一个地址生成状态机，外部与主控状态机和资源复用模块相连，在矩阵乘法过程中将使用资源复用模块中的3个存储器、浮点乘法器和浮点累加器，3个存储器分别存放的是求逆得到的上三角逆矩阵(LH)-1、矩阵(LH)-1和矩阵D-1相乘的结果(LH)-1*D-1和原矩阵的逆A-1。在接收到主控状态机模块输入的矩阵乘法模块起始信号后，通过矩阵乘法模块内部的地址生成状态机控制的存储器读写，先进行上三角矩阵(LH)-1与对角阵D-1相乘，再进行上三角矩阵(LH)-1*D-1与下三角矩阵L-1相乘，得到元矩阵的逆A-1。工作原理：cholesky分解矩阵求逆运算划分为三个步骤进行：(1)cholesky分解，根据公式A＝L*D*LH，将原矩阵A分解为单位下三角矩阵L(对角线本文档来自技高网...

【技术保护点】
一种基于分时复用技术的cholesky矩阵求逆系统，其特征在于，包括五个模块：主控状态机模块 、资源复用模块、cholesky分解模块、上三角矩阵求逆模块和矩阵相乘模块；所述的主控状态机模块与系统总线相连接，在求逆操作开始前，接收使能拉高，根据待求逆矩阵维度，从外部按行顺序接收原矩阵数据，存放在资源复用模块的存储器中，接收完成后，将接收使能拉低，并通过cholesky分解模块、上三角矩阵求逆模块、矩阵相乘模块分别进行cholesky分解、上三角矩阵求逆、矩阵乘法3个步骤进行cholesky矩阵求逆操作；求逆完成后向总线发出求逆结束信号，通知系统总线从存储模块读取运算结果。

【技术特征摘要】
1.一种基于分时复用技术的cholesky矩阵求逆系统，其特征在于，包括五个模块：主控状态机模块、资源复用模块、cholesky分解模块、上三角矩阵求逆模块和矩阵相乘模块；
所述的主控状态机模块与系统总线相连接，在求逆操作开始前，接收使能拉高，根据待求逆矩阵维度，从外部按行顺序接收原矩阵数据，存放在资源复用模块的存储器中，接收完成后，将接收使能拉低，并通过cholesky分解模块、上三角矩阵求逆模块、矩阵相乘模块分别进行cholesky分解、上三角矩阵求逆、矩阵乘法3个步骤进行cholesky矩阵求逆操作；求逆完成后向总线发出求逆结束信号，通知系统总线从存储模块读取运算结果。
2.如权利要求1所述的基于分时复用技术的cholesky矩阵求逆系统，其特征在于，所述的资源复用模块内部有3个存储器，1个单精度浮点复数乘法器，1个单精度浮点复数累加器；资源复用模块与主控状态机模块、cholesky分解模块、上三角矩阵求逆模块、矩阵相乘模块相连接；根据主控状态机的状态，资源复用内部的存储资源和运算资源受不同模块控制。
3.如权利要求1所述的基于分时复用技术的cholesky矩阵求逆系统，其特征在于，所述的cholesky分解模块是通过对原矩阵元素进行乘法，加法，除法运算将原矩阵A分解为单位下三角矩阵L、对角矩阵D，上三角矩阵LH；cholesky分解模块内部是一个地址生成状态机和两个单精度浮点除法器，外部与主控状态机和资源复用模块相连，在分解过程中将使用资源复用模块中的存储资源和浮点乘...

【专利技术属性】
技术研发人员：沙金，陈万雄，胡光辉，沈伟，任道，邵伟，张盟盟，潘红兵，李丽，李伟，
申请(专利权)人：南京大学，
类型：发明
国别省市：江苏;32