一种低复杂度的快速并行矩阵求逆方法技术

本发明专利技术公开了一种低复杂度的快速并行矩阵求逆方法，该发明专利技术包括以下步骤：第一步，对于给定矩阵A，令矩阵E为与矩阵A同阶的单位矩阵，将矩阵A与矩阵E组成扩展矩阵B，对矩阵B使用改进的Givens旋转(MSGR，Modified Squared Givens Rotations)，求得上三角矩阵U和其中矩阵U和的定义如下：使用SGR(Squared Givens Rotations)方法对矩阵A的QR分解进行变形与原始QR分解的关系为QA＝QDR，DR＝diag(R)，U＝DRR，DU＝diag(DRR)。MSGR方法可以去掉Givens旋转过程中的平方根运算同时减少除法运算，显著降低算法复杂度；第二步，回代法求上三角矩阵U的逆矩阵U-1；第三步，矩阵求逆本发明专利技术去掉了大量的除法以及平方根运算，降低了算法复杂度，可以应用于无线通信、信号处理以及数值计算等领域的矩阵求逆问题。

【技术实现步骤摘要】
一种低复杂度的快速并行矩阵求逆方法
本专利技术涉及一种去平方根运算、减少除法运算的低复杂度的快速并行矩阵求逆方法，主要涉及无线通信、信号处理以及数值计算等领域。
技术介绍
随着无线通信技术的不断发展，用户对通信系统的可靠性和有效性要求越来越高。为了提高系统的传输速率、增加系统的频谱效率，多输入多输出(MultipleInputMultipleOutput，MIMO)技术作为一种关键技术得到了广泛研究。IEEE802.11n、IEEE802.11ac、IEEE802.11ad以及3GPP-LTE等无线通信标准，都采用了MIMO技术。MIMO通信系统，在接收端的设计相当复杂。接收端信道估计会涉及到大量的矩阵求逆运算，尤其随着发射天线以及接收天线数目的增加，矩阵的维数也随着增加，矩阵求逆的运算复杂度也会大量增加，导致矩阵求逆的算法设计以及硬件实现难度大大提高，这也成为当前MIMO系统亟待解决的一个技术难题。因此，提出一个有效的矩阵求逆算法，可以提高通信系统的效率、降低硬件实现的成本，对于通信系统的发展具有很重要的现实意义。Givens旋转在数值线性代数中的主要用途是在向量或矩阵中介入零，常用于矩阵的QR分解，Givens旋转具有稳定性高、易于映射到脉动阵列、便于并行处理等优点，改进的Givens旋转，去掉Givens旋转过程中的平方根以及除法运算有利于降低算法复杂度。
技术实现思路
专利技术目的：本专利技术目的在于提供一种基于改进Givens旋转的矩阵求逆方法，去掉Givens旋转过程中的平方根、除法运算，大大降低算法复杂度。为更好的理解本专利技术的技术方案，首先对技术理论基础说明如下：对于矩阵A进行QR分解，可以分解为一个酉矩阵和一个上三角矩阵的乘积，即A＝QR，其中Q是一个酉矩阵，R是一个上三角矩阵，使用SGR(SquaredGivensRotations)方法，矩阵A的QR分解可以变形为：(公式1)其中QA＝QDR，DR＝diag(R)，U＝DRR为上三角矩阵，则有：(公式2)本专利技术的主要目标是提供一种省略平方根运算的同时减少除法运算的低复杂度并行矩阵求逆方法，主要是获得上三角矩阵U及相应的逆矩阵和技术方案：本专利技术提供一种低复杂度的快速并行矩阵求逆方法，应用于MIMO通信系统接收端的信道估计和接收端的信号均衡处理，对于一个发送天线数为M，接收天线数为N的MIMO通信系统，其接收机信号可以表示为r＝Hs+n，其中，r表示接收信号，是维数为M的列向量；s表示发送信号，是维数为N的列向量；H表示信道矩阵，是维数为M×N的矩阵；n表示加性高斯白噪声，是维数为M的列向量。接收端的均衡器从接收到的信号r估计出发送信号s，常见的均衡算法包括迫零算法和最小均方误差算法。基于迫零算法的均衡表达式为其中表示基于迫零算法的均衡器对发送信号的估计量；基于最小均方误差算法的均衡表达式为其中表示基于最小均方误差算法的均衡器对发送信号的估计量，上标H表示矩阵的共轭转置，IM表示M维单位矩阵，表示噪声的平均功率。将待求逆矩阵(HHH)或矩阵记为A，E为与A同阶的单位矩阵，分别表示如下：(公式3)本专利技术所提供的方法主要包括以下步骤：步骤1：采用MSGR(ModifiedSquaredGivensRotations)方法求解矩阵A按(公式1)分解时的上三角矩阵U以及步骤2：采用回代法求上三角矩阵U的逆矩阵U-1；步骤3：将步骤2中获得的U-1和步骤1中获得的相乘，得到矩阵A的逆矩阵。所述步骤1的具体步骤如下：步骤1.1：定义矩阵B为矩阵A与E组成的扩展矩阵，即B表示为：(公式4)步骤1.2：记矩阵B的维数为M1×M2，如果M1＜2，则执行步骤1.6；如果M1≥2，对矩阵B调用矩阵消零模块，矩阵消零模块输出矩阵(公式5)其中符号×表示非零元素值，符号表示用向量标量对表示矩阵的各行元素，即用<m,n>表示矩阵的第1行元素，用<pi,qi>表示矩阵的第i行元素，i＝2,3,…,M1，各向量标量对的取值由矩阵消零模块计算得到。步骤1.3：定义向量用来保存向量m，定义标量用来保存n，即令步骤1.4：将矩阵B重新赋值且维数减1，即M1＝M1-1，M2＝M2-1且令表示矩阵B由矩阵的第2行到第M1行以及第2列到第M2列元素组成。步骤1.5：对矩阵B重复步骤1.2到步骤1.4。步骤1.6：执行完上述步骤，矩阵B的维数为1×(N+1)，定义向量hN，令hN＝B*(1,1)·B(1,1:N+1)，其中B(1,1)表示矩阵B的第一行第一列位置的元素值，B(1,1:N)表示矩阵B的第一行元素，定义标量gN＝1。步骤1.7：令hx(j)表示向量hx的第j个元素，x＝1,2,…,N，j＝1,2,…,2N-x+1，则上三角矩阵U和可以分别表示成：(公式6)所述步骤2采用回代法求上三角矩阵U的逆矩阵U-1的具体方法为：上三角矩阵U的逆矩阵也是上三角矩阵，U-1的形式可以表示为：(公式7)采用下面的回代公式，按照w11,w12,…w1n,w22,…,w2n,w33,…,wnn的顺序求出W中各个元素的值：(公式8)上述步骤1.2中的矩阵消零模块对复数矩阵B进行消零运算的具体方法为：对任意M1×M2维复数矩阵B，记为：其中bk表示矩阵B的第k行元素，pk表示维数与bk相同的行向量，qk表示标量，符号表示分别用向量标量对<pk,qk>表示矩阵B的第k行；pk初始化为矩阵B的第k行元素组成的行向量，即pk＝bk，qk初始化为1，k＝1,2,…,M1。定义向量s，t为M1维列向量且元素初始值为s(l)＝1和t(l)＝1，其中s(l)和t(l)分别表示向量s和t的第l个元素，l＝1,2,…,M1。对复数矩阵B进行消零运算的具体步骤为：步骤1.2.1：旋转行向量b1和b2将b21消为0。1)将b1用向量m和标量n表示为对m，n，进行初始化：p1(r)表示向量p1的第r个元素，如果p1(1)＝0，则m＝s(1)p1，n＝t(1)q1，如果p1(1)≠0，则m＝s(1)p1(1)*p1，2)使用MSGR旋转，对各变量进行更新：约定表示x更新后的值，则向量b1和b2更新为：下面分情况给出以及和的表达式：当m(1)≠0时当m(1)＝0，且p2(1)＝0时当m(1)＝0，且p2(1)≠0时(公式9)3)溢出处理：exp(c)表示浮点数c的指数值，令x，y为两个复数，约定对于复数num，real(num)表示该复数的实部，imag(num)表示该复数的虚部，则x和y的实部和虚部可以分别用x1，x2，y1，y2表示如下：x1＝real(x)，x2＝imag(x)，y1＝real(y)，y2＝imag(y)，令约定语句IFTHEN如果条件A成立，则执行B操作。对复数对<x,y>进行溢出处理可以用以下语句表示：IFTHENIFTHEN根据上述方法分别对复数对(或复数向量对)进行溢出处理。4)将更新后的以及和分别赋值给m，n，p2，q2，以及s(2)和t(2)，即令步骤1.2.2：将b31，b41，……，位置元素消为0。将b2，p2，q2，s(2)，t(2)替换为br，pr，qr，s(r)，t(r)，r＝3,4,…,M1，重复步骤1.2.1中2)到4)的操作，可本文档来自技高网...

【技术保护点】
一种低复杂度的快速并行矩阵求逆方法，应用于MIMO通信系统接收端的信道估计和接收端的信号均衡，对于一个发送天线数为M，接收天线数为N的MIMO通信系统，其接收机信号可以表示为r＝Hs+n，其中，r表示接收信号，是维数为M的列向量；s表示发送信号，是维数为N的列向量；H表示信道矩阵，是维数为M×N的矩阵；n表示加性高斯白噪声，是维数为M的列向量；接收端的均衡器从接收到的信号r估计出发送信号s，均衡算法包括迫零算法和最小均方误差算法，基于迫零算法的均衡表达式为其中表示基于迫零算法的均衡器对发送信号的估计量；基于最小均方误差算法的均衡表达式为其中表示基于最小均方误差算法的均衡器对发送信号的估计量，上标H表示矩阵的共轭转置，IM表示M维单位矩阵，表示噪声的平均功率，将待求逆矩阵(HHH)或矩阵记为A，其特征在于，所述矩阵求逆方法将矩阵A的QR分解变形为其中QA＝QDR，DR＝diag(R)，U＝DRR，Q是一个酉矩阵，R是一个上三角矩阵，所述矩阵求逆包括如下步骤：(1)将矩阵A与同阶单位矩阵E组成的扩展矩阵的各行元素分别用向量标量对表示，基于向量标量对进行迭代消零运算来求解矩阵A按分解的上三角矩阵U以及(2)计算上三角矩阵U的逆矩阵U‑1；(3)将步骤(2)中获得的U‑1和步骤(1)中获得的相乘，得到矩阵A的逆矩阵。...

【技术特征摘要】
1.一种低复杂度的快速并行矩阵求逆方法，应用于MIMO通信系统接收端的信道估计和接收端的信号均衡，对于一个发送天线数为M，接收天线数为N的MIMO通信系统，其接收机信号可以表示为r＝Hs+n，其中，r表示接收信号，是维数为M的列向量；s表示发送信号，是维数为N的列向量；H表示信道矩阵，是维数为M×N的矩阵；n表示加性高斯白噪声，是维数为M的列向量；接收端的均衡器从接收到的信号r估计出发送信号s，均衡算法包括迫零算法和最小均方误差算法，基于迫零算法的均衡表达式为其中表示基于迫零算法的均衡器对发送信号的估计量；基于最小均方误差算法的均衡表达式为其中表示基于最小均方误差算法的均衡器对发送信号的估计量，上标H表示矩阵的共轭转置，IM表示M维单位矩阵，表示噪声的平均功率，将待求逆矩阵(HHH)或矩阵记为A，其特征在于，所述矩阵求逆方法将矩阵A的QR分解变形为其中QA＝QDR，DR＝diag(R)，U＝DRR，Q是一个酉矩阵，R是一个上三角矩阵，所述矩阵求逆包括如下步骤：(1)将矩阵A与同阶单位矩阵E组成的扩展矩阵的各行元素分别用向量标量对表示，所述向量标量对表示形式为<a,b>，其中a表示向量，b表示标量，基于向量标量对进行迭代消零运算来求解矩阵A按分解的上三角矩阵U以及(2)计算上三角矩阵U的逆矩阵U-1；(3)将步骤(2)中获得的U-1和步骤(1)中获得的相乘，得到矩阵A的逆矩阵。2.根据权利要求1所述的低复杂度的快速并行矩阵求逆方法，其特征在于：所述步骤(1)求解矩阵A按分解的上三角矩阵U以及的具体步骤包括：(1.1)定义矩阵B为矩阵A与E组成的扩展矩阵，E为与A同阶的单位矩阵，矩阵B表示为：(1.2)记矩阵B的维数为M1×M2，如果M1＜2，则执行步骤(1.6)；如果M1≥2，对矩阵B进行消零运算，输出矩阵其中符号×表示非零元素值，符号表示用向量标量对表示矩阵的各行元素，即用<m,n>表示矩阵的第1行元素，用<pi,qi>表示矩阵的第i行元素，i＝2,3,…,M1；(1.3)定义向量令定义标量令(1.4)将矩阵B重新赋值，令并将维数减1，即M1＝M1-1，M2＝M2-1；(1.5)对矩阵B重复步骤(1.2)到步骤(1.4)；(1.6)执行完上述步骤，矩阵B的维数为1×(N+1)，定义向量hN，令hN＝B*(1,1)·B(1,1:N+1)，其中B(1,1)表示矩阵B的第一行第一列位置的元素值，B(1,1:N)表示矩阵B的第一行元素，定义标量gN＝1；(1.7)令hx(j)表示向量hx的第j个元素，x＝1,2,…,N，j＝1,2,…,2N-x+1，则得到上三角矩阵和3.根据权利要求2所述的低复杂度的快速并行矩阵求逆方法，其特征在于，所述步骤(1.2)中对矩阵B进行消零运算的具体步骤包括：(1.2.1)将b1用向量m和标量n表示为对m，n，进行初始化：p1(r)表示向量p1的第r个元素，如果p1(1)＝0，则m＝s(1)p1，n＝t(1)...

【专利技术属性】
技术研发人员：何世文，余登高，黄永明，王海明，杨绿溪，张军，
申请(专利权)人：江苏中兴微通信息科技有限公司，
类型：发明
国别省市：江苏;32