CN101739387B - 奇异值分解方法及装置 - Google Patents
奇异值分解方法及装置 Download PDFInfo
- Publication number
- CN101739387B CN101739387B CN2009102120863A CN200910212086A CN101739387B CN 101739387 B CN101739387 B CN 101739387B CN 2009102120863 A CN2009102120863 A CN 2009102120863A CN 200910212086 A CN200910212086 A CN 200910212086A CN 101739387 B CN101739387 B CN 101739387B
- Authority
- CN
- China
- Prior art keywords
- matrix
- treatment circuit
- vector
- coupled
- unitary
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
- H04L25/0248—Eigen-space methods
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/0204—Channel estimation of multiple channels
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Power Engineering (AREA)
- Radio Transmission System (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种奇异值分解方法及装置。所述奇异值分解方法包含:(a)简化矩阵以导出简化矩阵;(b)对简化矩阵执行迭代矩阵乘法以产生迭代矩阵;(c)依据迭代矩阵乘法的迭代次数撷取迭代矩阵的向量;(d)从上述向量中解相关简化矩阵以更新简化矩阵;(e)重复步骤(b)、(c)及(d)直到预设数量的向量被导出;以及(f)存储依据奇异值分解中预设数量的向量确定的第一酉矩阵。以上所述的奇异值分解方法及装置可降低分解信道矩阵的总的计算复杂度且具有快速收敛的特性。
Description
技术领域
本发明有关于奇异值分解(singular value decomposing)方法及相关装置,且特别有关于具有快速收敛(fast-convergence)及较少计算量的奇异值分解方法及装置。
背景技术
如图1所示,在无线通信系统中,多入多出(multiple-input multiple-output,MIMO)收发器利用天线阵列可显著增强传送数据的吞吐量(throughput)。图1是传统MIMO系统的天线阵列的信道的示意图。发送器有mt个天线,而接收器则有mr个天线。接收器每一天线Y1~Ymr可接收自发送器每一天线X1~Xmt传送的信号。换句话说,第一发送器的天线X1可具有mr个发送信道(transmitting channel)h11,h21,...,hmr1,并且第二发送器的天线也具有mr个发送信道h12,h22,...,hmr2等等。因此,MIMO系统中的天线被用来增加多样性以抑制多径衰落(multi-path fading)或空间分离元件(spatially separatedevice)。此外,信号的线性度(linearity of the signal)与发送器与接收器使用的天线的数量成比例。天线间形成的干扰是MIMO系统的主要问题。各种不同方法可被使用,以解决MIMO系统的内干扰(inter-interference)问题。奇异值分解被认为对于消除干扰具有最佳性能。然而,解决奇异值分解的传统算法需要大量迭代运算,其可导致收敛时间(converging time)及硬件复杂度的增加。
发明内容
有鉴于此,特提供以下技术方案:
本发明的实施例提供一种酉矩阵计算装置,用于在对矩阵执行奇异值分解时计算第一酉矩阵。上述酉矩阵计算装置包含第一处理电路、第二处理电路、第三处理电路、第四处理电路、开关电路及存储模块。第一处理电路用于简化矩阵以导出简化矩阵;第二处理电路耦接于第一处理电路,用于对简化矩阵执行迭代矩阵乘法以产生迭代矩阵;第三处理电路耦接于第二矩阵电路,用于依据迭代矩阵乘法的迭代次数撷取迭代矩阵的单位向量;第四处理电路耦接于第三处理电路,用于将单位矩阵与所述单位向量导出的特征向量相减以得到减法矩阵,并将所述减法矩阵与所述简化矩阵相乘以得到解相关矩阵来更新简化矩阵;开关电路耦接于第四处理电路及第二处理电路之间,用于选择性将解相关矩阵耦接至第二处理电路,直到预设数量的单位向量被导出;存储模块耦接于第三处理电路,用于存储依据奇异值分解中的预设数量的单位向量确定的第一酉矩阵。
本发明实施例另提供一种包含所述酉矩阵计算装置的奇异值分解装置,所述奇异值分解装置还包含第六处理电路以及第七处理电路。第六处理电路耦接至所述存储模块,用于将所述矩阵分别与所述多个第一正规化向量相乘以产生多个合成向量;以及第七处理电路耦接至所述第六处理电路,用于导出所述奇异值分解的第二酉矩阵。
本发明的实施例另提供一种包含所述酉矩阵计算装置的奇异值分解装置,所述奇异值分解装置还包含第八处理电路以及第九处理电路。第八处理电路耦接至所述存储模块,用于将所述矩阵分别与所述预设数量的单位向量相乘,以产生多个合成向量;以及第九处理电路耦接至所述第八处理电路,用于导出由分别正规化所述多个合成向量产生的多个正规化向量组成的第二酉矩阵,并存储所述第二酉矩阵。
以上所述的酉矩阵计算装置以及奇异值分解装置可降低分解信道矩阵的总的计算复杂度且具有快速收敛的特性。
附图说明
图1是传统MIMO系统的天线阵列中信道的示意图;
图2是依本发明实施例的信道矩阵的奇异值分解方法的流程图;
图3是依本发明第二实施例的奇异值分解装置的示意图;
图4是依本发明第三实施例的奇异值分解装置的示意图。
具体实施方式
在说明书及权利要求书当中使用了某些词汇来指称特定的组件。所属领域中的技术人员应可理解,制造商可能会用不同的名词来称呼同样的组件。本说明书及权利要求书并不以名称的差异来作为区分组件的方式,而是以组件在功能上的差异来作为区分的基准。在通篇说明书及后续的请求项当中所提及的“包含”是一开放式的用语,故应解释成“包含但不限定于”。另外,“耦接”一词在此包含任何直接及间接的电气连接手段。因此,若文中描述一第一装置耦接于一第二装置,则代表第一装置可直接电气连接于第二装置,或通过其它装置或连接手段间接地电气连接至第二装置。
请参考图2,图2是对信道矩阵H执行奇异值分解方法200的流程图,其中信道矩阵H是通过无线MIMO系统建模的复矩阵(complex matrix)。依据本发明的实施例,信道矩阵H是m*m复矩阵(m≥1)。换句话说,信道矩阵H的建模是通过m个发送器天线及m个接收器天线之间的信道来实现的,但此实施例仅用于说明,并非本发明的限制。换句话说,本领域的技术人员应可理解,本发明也可被用于在任意数量的发送器天线及任意数量的接收器天线之间建模的信道矩阵。此外,本发明的奇异值分解方法200用于将信道矩阵H分解为U∑VT的形式,其中矩阵U、V是酉矩阵(unitary matrix),矩阵∑是对角矩阵(diagonal matrix),并且VT是矩阵V的共轭转置矩阵(transpose conjugatematrix)。只要大体上相同的结果可被获取,图2所示流程图的步骤不必依照显示的精确顺序执行,也不必是连续的;即,中间可插入其他步骤。奇异值分解方法200包含如下步骤:
步骤202:简化信道矩阵H以导出简化矩阵(simplified matrix);
步骤204:对简化矩阵执行迭代矩阵乘法(iterative matrix multiplication)的n次迭代以产生迭代矩阵(iterated matrix);当n次迭代完成后,转至步骤206,其中n≥1,且迭代矩阵乘法的一次迭代包含如下步骤:
步骤2042:将简化矩阵自乘(multiply by itself)以产生迭代矩阵;
步骤2044:确定n次迭代是否已经完成;若是,转至步骤206,若否,转至步骤2046;
步骤2046:再分配一个字长(word-length)用于存储迭代矩阵以产生位移矩阵(shifted matrix);以及
步骤2048:利用位移矩阵更新简化矩阵,转至步骤2042;
步骤206:对迭代矩阵执行终止检测(termination checking)以确定是否迭代矩阵的所有单位向量(unit vectors)均被撷取;若终止检测显示迭代矩阵的所有单位向量均被撷取,转至步骤212,若终止的检测显示并非迭代矩阵的所有单位向量均被撷取,转至步骤208;
步骤208:撷取对应迭代矩阵任意行或列向量(row or column vector)的单位向量,并存储此单位向量,然后继续进行步骤210、212、214以及220;
步骤210:从步骤208中撷取的单位向量中解相关简化矩阵以产生解相关矩阵(de-correlated matrix),并利用解相关矩阵更新简化矩阵,转至步骤204;
步骤212:导出由步骤208中撷取的多个单位向量组成的矩阵的酉矩阵V;
步骤214:利用信道矩阵H与步骤208中撷取的单位向量相乘以产生一个合成向量(resultant vector);
步骤216:撷取并存储对应于合成向量的单位向量,转至步骤218及步骤220;
步骤218:导出由步骤216中产生的多个单位向量组成的矩阵的酉矩阵U;
步骤220:依据步骤208中产生的单位向量及步骤216中产生的单位向量导出并存储矩阵的一个特征值(eigenvalue);
步骤222:利用步骤220中产生的多个特征值产生对角矩阵∑。
首先,由于信道矩阵H中存在三个变量矩阵(variable matrix),当信道矩阵H被接收时,奇异值分解方法200减少信道矩阵H的一个变量矩阵以简化信道矩阵H。因此,在步骤202中,首先对信道矩阵H执行共轭转置操作以产生共轭转置矩阵HT;且随后将共轭转置矩阵HT与信道矩阵H相乘以产生简化矩阵P,如下方程所示:
其中M是信道矩阵H的列数(column number),组成对角矩阵∑的σi是降序排列的第i个奇异值,其意为σ1≥σ2≥...≥σM,且vi是酉矩阵V的第i个列向量(column vector)。请注意,尽管本实施例是通过消除酉矩阵U来简化信道矩阵H的,其并非本发明的限制。本领域的技术人员应可理解,信道矩阵H与共轭转置矩阵HT的相乘可被互换以产生简化矩阵P,从而消除酉矩阵V,也就是,P=HHT=U∑2UT,其也涵盖在本发明的范围内。
当简化矩阵P被获取后,对简化矩阵P执行迭代矩阵乘法以在步骤204中产生迭代矩阵其中在迭代矩阵乘法中执行n次迭代。然而,其并非本发明的限制。在本发明另一实施例中,迭代矩阵乘法对简化矩阵P执行m次迭代以产生迭代矩阵迭代矩阵乘法的说明如下述方程所示:
由于n次迭代后酉矩阵V并不改变,很明显,对简化矩阵P执行n次迭代之后,最大奇异值σ1与其他较小奇异值σ2,...,σM的差距(gap)将会增大。因此,迭代矩阵可被近似为其中仅有与最大奇异值σ1对应的向量(即,v1)被留在迭代矩阵中。迭代矩阵乘法的一次迭代首先将该简化矩阵P自乘以产生一阶简化矩阵(first order simplified matrix)P2(即,步骤2042),然后将一阶简化矩阵P2更新为简化矩阵P,以用于执行下一迭代。然而,在实践中,为避免一阶简化矩阵P2的最高有效位(most significant bit,MSB)的溢出,步骤2046可能要求迭代再分配一个字长,用于存储一阶简化矩阵P2以及用于产生位移矩阵。随后,位移矩阵被更新为简化矩阵P以用于执行下一次迭代直到n次迭代完成。换句话说,在迭代矩阵乘法的迭代中,每一迭代矩阵都被执行动态移位(dynamic shift)操作。更特别的,在步骤2046中,简化矩阵P中的所有元素均应当被二进制移位t位,以使得|2t*{Re(P),Im(P)}|的最大元素可在不丢失其MSB的情形下被表述。尽管本实施例是通过将矩阵P中的所有元素乘以2t来实施,但本发明并非仅限于此。本领域的技术人员应当可以理解,在此步骤中将矩阵P中的所有元素同时乘以其它值也能达成上述目的。此外,步骤2044确定n次迭代是否完成。若n次迭代尚未完成,迭代操作必须执行直到n次迭代完成。当n次迭代完成时,意味着迭代矩阵已被获取,对应于最大奇异值σ1的单位向量可在步骤208中自迭代矩阵中撷取。在单位向量可被撷取之前,在步骤206中对迭代矩阵执行终止检测,以确定对应于所有奇异值的所有单位向量是否均被撷取。若所有单位向量均被撷取,由多个单位向量组成的矩阵的酉矩阵V可在步骤212中导出。否则,迭代矩阵的单位向量在步骤208中被撷取。
对应于最大奇异值σ1的单位向量被获得后,信道H中的相关分量(correlated components)应当被消除,以用于导出下一奇异值及对应的单位向量。因此,在步骤210中,对应于最大奇异值σ1的单位向量自简化矩阵P中被解相关,以产生解相关矩阵Pd并利用解相关矩阵Pd更新简化矩阵P。下列方程阐述本发明的奇异值分解方法200的初始条件:
因此,获取对应于v1~vk的单位向量后,步骤210更进一步将单位矩阵(identity matrix)IM与对应于v1~vk的单位向量所形成的自相关矩阵(auto-correlation matrix)相减,以产生减法矩阵(subtracted matrix)Psub;并且将减法矩阵Psub与简化矩阵P相乘。相应地,解相关矩阵Pd可通过如下推导获得:
方法1
方法2
方法3
方法4
因此,解相关矩阵Pd可在步骤204中应用,以导出对应于vk+1的单位向量。
由于迭代矩阵在步骤204中连续产生,对应于各个奇异值σ1,...,σM的单位向量在步骤208中被逐一获取并存储。当对应于各个奇异值σ1,...,σM的单位向量均被获得后,可对单位向量执行正规化操作(nomalizing operation),以分别产生特征向量(eigenvector)v1,...,vM,其中特征向量v1,...,vM被用于导出矩阵的酉矩阵V(步骤212)。因此,矩阵的酉矩阵V可被获取。请注意,步骤208连同步骤212并非本发明用于获取特征向量v1,...,vM的限制。在本发明的另一实施例中,每一特征向量vi是通过正规化步骤204中得到的迭代矩阵任一列向量或行向量而获得。此外,在本发明的另一实施例中,对应于奇异值σi的特征向量vi可通过简单地消去迭代矩阵的部分而获得。
无论步骤208中对应于最大奇异值σ1的单位向量何时获得,所述单位向量在步骤204中被进一步利用来与信道矩阵H相乘以产生合成向量Pu,如下述方程所示:
类似于步骤208,对应于最大奇异值σ1的单位向量可在步骤216中自合成向量Pu中撷取。由于迭代矩阵在步骤204中连续产生,对应于各个奇异值σ1,...,σM的单位向量在步骤216中被逐一获得并且存储。当对应于各个奇异值σ1,...,σM的单位向量均被获得后,可对单位向量执行正规化操作以分别产生特征向量u1,...,uM,其中特征向量u1,...,uM被用于导出矩阵的酉矩阵U(步骤218)。因此,矩阵的酉矩阵U可被获得。请注意,在本发明的另一实施例中,每一特征向量ui是通过正规化步骤214中得到的合成向量Pu任一列向量或行向量而获得。
此外,无论步骤208及步骤216中对应于最大奇异值σ1的单位向量(分别对应于v1及u1)何时获得,所述单位向量在步骤220中被进一步利用来与信道矩阵H相乘以产生奇异值σi,如下述方程所示:
σi=ui THvi (5)
当所有奇异值σ1~σM均被获得后,对角矩阵∑可在步骤222中被导出。
依据本发明的上述揭露内容,奇异值分解方法200中无除法计算且无平方根计算。换句话说,奇异值分解方法200中无除法操作或者平方根操作,从而可降低分解信道矩阵H的总的计算复杂度(即,只有乘法及加法-multiplication and add,MAC),并且奇异值分解方法200具有快速收敛(fast-convergence)的特性。此外,奇异值分解方法200可适用于各种不同大小的信道矩阵,例如4*4矩阵或者4*3矩阵。此外,因为特征向量vi,ui以及σi可在本发明的实施例中分别导出,与传统方法相比,奇异值分解方法200具有高的平行度(high parallelism)。
请参考图3。图3是依据本发明第二实施例的奇异值分解装置300的示意图。奇异值分解装置300被用于对信道矩阵执行奇异值分解,以产生信道矩阵的酉矩阵。为更清楚地描述本发明,奇异值分解装置300被用于分解信道矩阵H以产生如图2所示的第一实施例的酉矩阵U。换句话说,信道矩阵H是由无线MIMO系统建模的复矩阵,且信道矩阵H是m*m复矩阵(m≥1),其通过m个发送器天线及m个接收器天线之间的信道建模,且两个实施例均用于说明的目的,并非本发明的限制。奇异值分解装置300包含第一处理电路301、第二处理电路302、第三处理电路303、第四处理电路304、开关电路305以及存储模块306。第一处理电路301接收信道矩阵H以简化信道矩阵H,从而导出上述方程(1)中所述的简化矩阵P。第二处理电路302耦接于第一处理电路301,用于对简化矩阵P执行迭代矩阵乘法的n次迭代以产生上述方程(2)中所述的迭代矩阵第三处理电路303耦接于第二处理电路302,用于在迭代矩阵乘法的n次迭代后撷取迭代矩阵的单位向量。第四处理电路304耦接于第三处理电路303,用于从第三处理电路中撷取的单位向量中解相关简化矩阵P,以产生解相关矩阵Pd。开关电路305耦接于第四处理电路304及第二处理电路302之间,用于选择性将解相关矩阵Pd耦接至第二处理电路302,直到预设数量的单位向量被导出;而存储模块306耦接于第三处理电路303,用于存储依据自第三处理电路303中撷取的预设数量的单位向量确定的酉矩阵V。
此外,第一处理电路301包含第一算术单元3011及第二算术单元3012。第一算术单元3011对信道矩阵H执行共轭转置操作以产生共轭转置矩阵HT,且第二算术单元3012耦接于第一算术单元3011,用于将共轭转置矩阵HT与信道矩阵H相乘以产生简化矩阵P。第二处理电路302包含第三算术单元3021及以移位模块(shifting module)3022。第三算术单元3021将解相关矩阵Pd进行n次自乘以产生迭代矩阵移位模块3022耦接于第三算术单元3021,用于再分配一个字长用于存储迭代矩阵以产生位移矩阵,并利用位移矩阵更新解相关矩阵Pd直到迭代次数与预设值相等,也就是n次迭代。第四处理电路304包含第四算术单元3041以及第五算术单元3042。第四算数单元3041将单位矩阵IM与第三处理电路303中撷取的单位向量导出的特征向量(即,)相减,以产生减法矩阵Psub。第五算术单元3042将减法矩阵Psub与简化矩阵P相乘以产生另一解相关矩阵Pd。存储模块306包含第五处理电路3061及存储单元3062。第五处理电路3061分别正规化第三处理电路303中撷取的多个单位向量以导出多个特征向量v1~vm。存储单元3062耦接于第五处理电路3061,用于存储由多个特征向量v1~vm组成的酉矩阵V。
首先,开关电路305切换至第四处理电路304,用于将简化矩阵P耦接至第二处理电路302。随后,开关电路305切换至移位模块3022。从而,第三算术单元3021将该简化矩阵P自乘以产生一阶简化矩阵P2,且一阶简化矩阵P2随后被输入移位模块3022,其中移位模块3022再分配简化矩阵P2的字长,以避免一阶简化矩阵P2的MSB溢出。随后位移矩阵被生成,并且开关电路305将位移矩阵耦接至第三算术单元3021以用于执行下一次迭代,直到n次迭代完成。换句话说,在迭代矩阵乘法的迭代中,移位模块3022对每一迭代矩阵执行动态移位操作。更特别的,简化矩阵P中的所有元素均被二进制移位t位,以使得|2t*{Re(P),Im(P)}|的最大元素可在不丢失其MSB的情形下被表述。尽管本实施例是通过将矩阵P中的所有元素乘以2t来实施,但本发明并非仅限于此。本领域的技术人员应当可以理解,在此步骤中将矩阵P中的所有元素同时乘以其它值也能达成上述目的。当第二处理电路302及开关电路305完成n次迭代时,第三处理电路303进一步对迭代矩阵执行终止检测,以确定对应于所有奇异值σ1,...,σM的所有单位向量是否均已被撷取。若并非对应于所有奇异值σ1,...,σM的所有单位向量均已被撷取,第三处理电路303撷取迭代矩阵的单位向量。此种状况下,当n次迭代完成时,迭代矩阵被获得,并且在第三处理电路303中,对应于最大奇异值σ1的单位向量可自迭代矩阵中撷取。获取对应于v1~vk的单位向量后,第四处理电路304执行下述计算以产生解相关矩阵Pd:
方法1
方法2
方法3
方法4
其中第二算术单元3041将单位矩阵IM与第三处理电路303中撷取的单位向量获取的特征向量(例如,)相减,并且第三算术单元3042将减法矩阵Psub与简化矩阵P相乘以产生另一解相关矩阵Pd。随后开关电路305将解相关矩阵Pd耦接至第二处理电路302,以获取对应于vk+1的单位向量。
当第三处理电路303的终止检测确定所有对应于奇异值的单位向量均已被撷取后,第五处理电路3061分别正规化第三处理电路303撷取的单位向量,以获取多个特征向量v1~vM。存储单元3062耦接至第五处理电路3061,用于存储由多个特征向量v1~vM组成的酉矩阵V。从而,信道矩阵H的酉矩阵V可被获得。请注意,对应于奇异值的特征向量vi也可在未于第五处理电路3061中正规化的情形下获得。换句话说,在本发明的另一实施例中,第五处理电路3061简单消去迭代矩阵的部分,以获得特征向量vi。
请一并参考图3及图4。图4是依据本发明的第三实施例的奇异值分解装置400的示意图。奇异值分解装置400被用于对信道矩阵H执行奇异值分解以产生信道矩阵H的酉矩阵V、酉矩阵U以及对角矩阵∑(如图2所示)。类似地,信道矩阵H是由无线MIMO系统建模的复矩阵,且信道矩阵H是m*m复矩阵(此处m≥1),其通过m个发送器天线及m个接收器天线之间的信道建模,且两个实施例均用于说明的目的,并非本发明的限制。奇异值分解装置400不仅包含奇异值分解装置300中的组件,更包含第六处理电路307、第七处理电路308、第八处理电路309以及第九处理电路3010。第六处理电路307耦接至存储模块306,用于分别将信道矩阵H与特征向量vi相乘以产生合成向量Pu。第七处理电路308耦接至第六处理电路307,用于导出信道矩阵H的酉矩阵U,其中酉矩阵U由分别正规化合成向量Pu产生的多个正规化向量(即,特征向量)ui组成。第八处理电路309耦接至第六处理电路307,用于将信道矩阵H与自存储模块306中获得的特征向量vi相乘,并与第六处理电路307中获得的特征向量ui相乘,以产生对应的奇异值σi。第九处理电路3010耦接至第八处理电路309,用于依据第八处理电路309中获得的奇异值σi导出信道矩阵H的对角矩阵∑。
依据图3所示的第三实施例,一旦特征向量vi产生,第六处理电路307将信道矩阵H与特征向量vi相乘以获取如方程(4)所示的合成向量Pu(即,)。正规化合成向量Pu后,特征向量ui可被获得。请注意,特征向量ui也可在未对合成向量Pu正规化的情形下获得。在本发明的另一实施例中,对应于奇异值的特征向量可通过简单地消去迭代矩阵Pu(也就是,)的部分来获得。
随后,信道矩阵H的酉矩阵U可依据特征向量ui(i=1~M)在第七处理电路中导出。
此外,依据上述揭露内容,一旦特征向量vi在第五处理电路3061中被获得,对应的特征向量ui也在第六处理电路307中被获得,从而对应的奇异值σi可在第八处理电路309中通过如方程(5)中所示,将信道矩阵H与特征向量vi及特征向量ui相乘而被获取。当所有奇异值σ1~σM均被获得后,对角矩阵∑可在第九处理电路3010中被导出。
依据本发明的上述揭露内容,奇异值分解装置400中无除法计算且无平方根计算。换句话说,奇异值分解装置400中无除法操作或者平方根操作,从而可降低分解信道矩阵H的总的计算复杂度(即,只有乘法及加法-MAC),并且具有快速收敛的特性。此外,奇异值分解装置400可适用于各种不同大小的信道矩阵,例如4*4矩阵或者4*3矩阵。此外,因为特征向量vi、ui以及σi可在本发明的实施例中分别导出,与传统装置相比,奇异值分解装置400具有高的平行性。
以上所述仅为本发明的较佳实施例,本领域相关的技术人员依据本发明的精神所做的等效变化与修饰,都应当涵盖在权利要求范围内。
Claims (10)
1.一种酉矩阵计算装置,用于在对矩阵执行奇异值分解时计算第一酉矩阵,其特征在于,所述酉矩阵计算装置包含:
第一处理电路,用于简化所述矩阵以导出简化矩阵;
第二处理电路,耦接于所述第一处理电路,用于对所述简化矩阵执行迭代矩阵乘法以产生迭代矩阵;
第三处理电路,耦接于所述第二处理电路,用于依据所述迭代矩阵乘法的迭代次数撷取所述迭代矩阵的单位向量;
第四处理电路,耦接于所述第三处理电路,用于将单位矩阵与所述单位向量导出的特征向量相减以得到减法矩阵,并将所述减法矩阵与所述简化矩阵相乘以得到解相关矩阵来更新所述简化矩阵;
开关电路,耦接于所述第四处理电路及所述第二处理电路之间,用于选择性将所述解相关矩阵耦接至所述第二处理电路,直到预设数量的单位向量被导出;以及
存储模块,耦接于所述第三处理电路,用于存储依据该奇异值分解中的所述预设数量的单位向量确定的第一酉矩阵。
2.如权利要求1所述的酉矩阵计算装置,其特征在于,所述第一处理电路包含:
第一算术单元,用于对所述矩阵执行共轭转置操作以产生共轭转置矩阵;以及
第二算术单元,耦接于所述第一算术单元,用于将所述矩阵与所述共轭转置矩阵相乘以产生所述简化矩阵。
3.如权利要求1所述的酉矩阵计算装置,其特征在于,所述第二处理电路包含:
第三算术单元,用于将所述简化矩阵自乘以产生该迭代矩阵;
移位模块,耦接于所述第三算术单元,用于再分配一个字长用于存储所述迭代矩阵,以产生位移矩阵,并且利用所述位移矩阵更新所述简化矩阵,直到所述迭代次数与预设值相同。
4.如权利要求1所述的酉矩阵计算装置,其特征在于,所述第四处理电路包含:
第四算术单元,用于将所述单位矩阵与所述第三处理电路中撷取的所述多个单位向量导出的所述特征向量相减,以产生所述减法矩阵;以及
第五算术单元,耦接至所述第四算术单元,用于将所述减法矩阵与所述简化矩阵相乘以得到所述解相关矩阵来更新所述简化矩阵。
5.如权利要求1所述的酉矩阵计算装置,其特征在于,所述矩阵是多入多出信道矩阵。
6.如权利要求1所述的酉矩阵计算装置,其特征在于,所述存储模块包含:
第五处理电路,用于分别正规化所述多个预设数量的单位向量以导出多个第一正规化向量;以及
存储单元,耦接于所述第五处理电路,用于存储由所述多个第一正规化向量组成的所述第一酉矩阵。
7.一种包含权利要求6所述的酉矩阵计算装置的奇异值分解装置,其特征在于,所述奇异值分解装置还包含:
第六处理电路,耦接至所述存储模块,用于将所述矩阵分别与所述多个第一正规化向量相乘以产生多个合成向量;以及
第七处理电路,耦接至所述第六处理电路,用于导出所述奇异值分解的第二酉矩阵,其中所述第二酉矩阵由分别正规化所述多个合成向量产生的多个第二正规化向量组成。
8.如权利要求7所述的奇异值分解装置,其特征在于,在所有所述多个第一正规化向量被第五处理电路导出前,所述第六处理电路将所述矩阵与由所述第五处理电路导出的一个正规化向量相乘。
9.一种包含权利要求1所述的酉矩阵计算装置的奇异值分解装置,其特征在于,所述奇异值分解装置还包含:
第八处理电路,耦接至所述存储模块,用于将所述矩阵分别与所述预设数量的单位向量相乘,以产生多个合成向量;以及
第九处理电路,耦接至所述第八处理电路,用于导出由分别正规化所述多个合成向量产生的多个正规化向量组成的第二酉矩阵,并存储所述第二酉矩阵。
10.如权利要求9所述的奇异值分解装置,其特征在于,在所有预设数量的单位向量被导出前,所述第八处理电路将所述矩阵与所述预设数量的单位向量其中之一相乘,以导出合成向量。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/264,258 | 2008-11-04 | ||
US12/264,258 US8321488B2 (en) | 2008-11-04 | 2008-11-04 | Singular value decomposing method and related singular value decomposing device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101739387A CN101739387A (zh) | 2010-06-16 |
CN101739387B true CN101739387B (zh) | 2012-07-04 |
Family
ID=42131376
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102120863A Expired - Fee Related CN101739387B (zh) | 2008-11-04 | 2009-11-04 | 奇异值分解方法及装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8321488B2 (zh) |
CN (1) | CN101739387B (zh) |
TW (1) | TWI393394B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104618293B (zh) * | 2015-01-27 | 2017-11-28 | 东南大学 | 一种平滑奇异值分解的酉变换矩阵的优化方法 |
CN105763880B (zh) * | 2016-02-19 | 2019-02-22 | 北京大学 | 消除视频编码图像块效应的方法 |
CN107977534A (zh) * | 2017-12-21 | 2018-05-01 | 北京华大九天软件有限公司 | 层次化电路仿真中避免矩阵奇异的方法 |
US11190244B1 (en) * | 2020-07-31 | 2021-11-30 | Samsung Electronics Co., Ltd. | Low complexity algorithms for precoding matrix calculation |
EP4164137A1 (en) * | 2021-10-05 | 2023-04-12 | Nokia Solutions and Networks Oy | Computation of beamforming parameters |
CN114357322B (zh) * | 2022-03-16 | 2022-06-14 | 交通运输部规划研究院 | 一种基于时空分解的船舶轨迹异常检测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1614906A (zh) * | 2003-11-05 | 2005-05-11 | 索尼株式会社 | 无线通信系统、无线通信方法以及无线通信设备 |
CN101156403A (zh) * | 2005-04-01 | 2008-04-02 | 美商内数位科技公司 | 频道矩阵奇异值分解方法及装置 |
CN101262311A (zh) * | 2008-04-18 | 2008-09-10 | 西安电子科技大学 | 多用户mimo系统下行链路发射端信息数据的预处理方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5615288A (en) * | 1993-04-26 | 1997-03-25 | Fuji Xerox Co., Ltd. | Singular value decomposition coding and decoding apparatuses |
US8483761B2 (en) * | 2006-01-18 | 2013-07-09 | Intel Corporation | Singular value decomposition beamforming for a multiple-input-multiple-output communication system |
JPWO2007088579A1 (ja) * | 2006-01-31 | 2009-06-25 | 三菱電機株式会社 | 無線送信装置、無線受信装置および無線通信システム |
JP2008060877A (ja) * | 2006-08-31 | 2008-03-13 | Hitachi Ltd | Mimo無線データ伝送システム |
US7949070B2 (en) * | 2006-11-15 | 2011-05-24 | Sanyo Electric Co., Ltd. | Method and apparatus for executing MIMO eigenmode transmission |
JP4518184B2 (ja) * | 2008-05-02 | 2010-08-04 | ソニー株式会社 | 無線通信装置、無線通信方法、無線通信システム及びコンピュータプログラム |
KR101470501B1 (ko) * | 2008-08-20 | 2014-12-08 | 삼성전자주식회사 | 양자화된 채널 상태 정보에 기반하여 데이터를 전송하는 장치 및 방법 |
-
2008
- 2008-11-04 US US12/264,258 patent/US8321488B2/en active Active
-
2009
- 2009-11-04 CN CN2009102120863A patent/CN101739387B/zh not_active Expired - Fee Related
- 2009-11-04 TW TW098137397A patent/TWI393394B/zh not_active IP Right Cessation
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1614906A (zh) * | 2003-11-05 | 2005-05-11 | 索尼株式会社 | 无线通信系统、无线通信方法以及无线通信设备 |
CN101156403A (zh) * | 2005-04-01 | 2008-04-02 | 美商内数位科技公司 | 频道矩阵奇异值分解方法及装置 |
CN101262311A (zh) * | 2008-04-18 | 2008-09-10 | 西安电子科技大学 | 多用户mimo系统下行链路发射端信息数据的预处理方法 |
Also Published As
Publication number | Publication date |
---|---|
TW201110625A (en) | 2011-03-16 |
TWI393394B (zh) | 2013-04-11 |
US20100111213A1 (en) | 2010-05-06 |
US8321488B2 (en) | 2012-11-27 |
CN101739387A (zh) | 2010-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bellili et al. | Generalized approximate message passing for massive MIMO mmWave channel estimation with Laplacian prior | |
Lapidoth et al. | Capacity bounds via duality with applications to multiple-antenna systems on flat-fading channels | |
CN101739387B (zh) | 奇异值分解方法及装置 | |
CN108964726B (zh) | 一种低复杂度的大规模mimo上行链路传输信道估计方法 | |
CN109150278B (zh) | 基于改进牛顿迭代的massive MIMO信号检测方法 | |
CN109951214B (zh) | 一种适用于大规模mimo系统的信号检测方法 | |
Kim et al. | Practical channel estimation and phase shift design for intelligent reflecting surface empowered MIMO systems | |
Zhang et al. | Leveraging the restricted isometry property: Improved low-rank subspace decomposition for hybrid millimeter-wave systems | |
CN105071843A (zh) | 大规模mimo系统低复杂度多项式展开矩阵求逆方法及应用 | |
CN102271006B (zh) | 无线通信系统中的通信方法和装置 | |
Ghavami et al. | Noncoherent massive MIMO detection by expectation propagation | |
US7561644B2 (en) | Method, system and device for reducing co-channel interference | |
Zheng et al. | Hybrid Driven Learning for Channel Estimation in Intelligent Reflecting Surface Aided Millimeter Wave Communications | |
Moustakas et al. | Random matrix theory of multi-antenna communications: the Ricean channel | |
US11256780B1 (en) | Methods and apparatus for fast eigenvalue decomposition of Hermitian matrices | |
CN109379116A (zh) | 基于切比雪夫加速法与sor算法的大规模mimo线性检测算法 | |
CN106850096A (zh) | 空间相关性估算的方法及装置 | |
Kreyndelin et al. | Decreasing of computational complexity of demodulation algorithms in multi-antenna systems due to application of fast algorithms | |
CN108809389A (zh) | 基于Givens旋转矩阵复用的QR分解块对角化预编码实现方法 | |
Zuo et al. | Distributed channel estimation algorithm for mmWave massive MIMO communication systems | |
KR100965100B1 (ko) | 적응적 빔포밍 장치 및 그 방법 | |
Riegler et al. | Noncoherent SIMO pre-log via resolution of singularities | |
Hong et al. | Channel Estimation-Free Deep Direct Beamforming With Low Complexity in mmWave Massive MIMO | |
Esfandiari et al. | ADMM-based solution for mmWave UL channel estimation with one-bit ADCs via sparsity enforcing and Toeplitz matrix reconstruction | |
Crâşmariu et al. | Optimized block-diagonalization precoding technique using givens rotations QR decomposition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120704 Termination date: 20191104 |
|
CF01 | Termination of patent right due to non-payment of annual fee |