CN109639607A - 一种低复杂度的非正交多址接入系统信号检测方法 - Google Patents
一种低复杂度的非正交多址接入系统信号检测方法 Download PDFInfo
- Publication number
- CN109639607A CN109639607A CN201910055664.0A CN201910055664A CN109639607A CN 109639607 A CN109639607 A CN 109639607A CN 201910055664 A CN201910055664 A CN 201910055664A CN 109639607 A CN109639607 A CN 109639607A
- Authority
- CN
- China
- Prior art keywords
- matrix
- column
- row
- vector
- indicates
- 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.)
- Granted
Links
Classifications
-
- 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/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03891—Spatial equalizers
- H04L25/03898—Spatial equalizers codebook-based design
- H04L25/0391—Spatial equalizers codebook-based design construction details of matrices
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
- H04L1/0045—Arrangements at the receiver end
- H04L1/0047—Decoding adapted to other signal detection operation
- H04L1/0048—Decoding adapted to other signal detection operation in conjunction with detection of multiuser or interfering signals, e.g. iteration between CDMA or MIMO detector and FEC decoder
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
- H04L1/0045—Arrangements at the receiver end
- H04L1/0047—Decoding adapted to other signal detection operation
- H04L1/005—Iterative decoding, including iteration between signal detection and decoding operation
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)
Abstract
本发明公开了一种低复杂度的非正交多址接入系统信号检测方法,属于无线通信技术领域。本发明依据信号的SNR值排序,通过对矩阵进行排序QR检测,得到一个正交矩阵和一个上三角矩阵,利用列向量之间的对称性,将矩阵中的一对相邻列与另一对相邻列交换,减少了排序QR分解算法中的正交计算,再根据迭代公式对接收端输出的接收信号矩阵进行检测,以得到发射信号的估计值。在系统性能没有降低的前提下,本发明降低了现有算法的复杂度,本发明的计算复杂度几乎是传统方法的50%。
Description
技术领域
本发明涉及到无线通信技术,特别是涉及一种低复杂度的非正交多址接入系统信号检测方法。
背景技术
在2014年,IMT-2020推进组提出了5G白皮书。在第五代通信系统中,将提出更高的需求,5G相比与4G将会满足1000倍甚至更多的容量需求、10到100倍的用户设备连接,5到15倍频谱效率提升。面对5G海量连接场景,现有的4GOFDMA技术有一定的局限性。为了进一步提升系统容量和效率,基于非正交多址接入(Non-Orthogonal Multiple Access,简称NOMA)的无线传输技术成为未来移动通信系统的热门技术。
虽然NOMA技术提高了无线通信的容量和效率,但是它也意味着更复杂的配置,这对于数据检测来说是棘手的。随着用户和基站端天线数目的增加,依赖于最大似然(ML)准则或最大后验(MAP)准则的最优检测方法复杂度降成指数增长,限制了他们的实际应用。为了降低相关计算复杂度,许多线性检测方案被相机提出。其中,迫零算法和最小均方误差算法都是性能接近最优的选择,但是他们包含直接矩阵求逆运算,复杂度较高。
为了避免线性检测器的主要复杂度的来源,即矩阵求逆运算,研究人员开始求助于近似方法或迭代方法。QR分解信号检测算法属于线性检测算法,利用接收信号组合中的串行干扰消除技术,它有效的避免了多次矩阵求伪逆的运算,从而减小了运算复杂度;排序的QR分解信号检测算法则在减小运算复杂度的同时改善了系统的性能。通过分析己有的基于串行干扰消除原则的排序QR分解检测算法,虽然它们能保障信号输出端检测到的信号尽可能的大,但是算法复杂度较高,需要进行多次正交运算。
发明内容
针对现有技术的不足,本发明的目的在于提供一种低复杂度的非正交多址接入系统信号检测方法,用于解决现有技术中存在的算法复杂度较高,需要进行多次正交运算的问题。
为解决现有技术问题,本发明采取的技术方案为:
一种低复杂度的非正交多址接入系统信号检测方法,所述方法包括以下步骤:
获取加性高斯白噪声信道的复信道系数矩阵、接收信号向量、发送信号向量、发送天线数和接收天线数,对所述复信道系数矩阵进行实数化处理,得到实数矩阵,所述实数矩阵的行数和列数分别是复信道系数矩阵的两倍;
根据发送天线数、接收天线数对所述实数矩阵进行对称的排序QR分解,得到一个正交矩阵Q和一个上三角矩阵R;
根据所述发送信号向量对所述接收信号向量进行实数化处理、检测,得到发射信号估计值。
进一步的,所述复信道系数矩阵采用公式(1)进行实数化处理:
式中:G为复信道系数矩阵;H为实数矩阵。
进一步的,所述排序QR的分解方法包括如下步骤:
令Q=H,R=0,排列向量p=[1,(1+C),2,(2+C),...,C,2C],其中C表示实数矩阵H的列数,C的大小与发送天线数相等;
将矩阵Q中的矢量按照向量p重新排列,此时矩阵Q中相邻两列向量对称,即Q2i=[-QM+1:2M,2i-1 T,Q1:M,2i-1 T]T,且Q2i-1与Q2i正交,
其中M表示矩阵Q的列数,M的大小与接收天线数相等,(·)T表示矩阵转置,Q2i表示矩阵Q的第2i列,QM+1:2M,2i-1表示取矩阵Q第2i-1列中第M+1行到2M行的元素,Q1:M,2i-1表示取矩阵Q第2i-1列中第1行到M行的元素,Q2i-1表示矩阵Q的第2i-1列;
计算Q中奇数列的二范数模值的平方,记为 其中表示对矩阵Q的第2i-1列做转置;i表示第一迭代计数器,N表示迭代总次数,与发送天线数相等;
计算矩阵Q和矩阵R中的元素。
进一步的,计算矩阵Q和矩阵R中的元素过程包括:
步骤2.41:设置迭代总次数N,令第二迭代计数器j=0;
步骤2.42:k=argminl=j,...,Nnorml,其中k表示二范数模值平方最小向量的坐标,l表示第三迭代计数器,norml表示向量norm的第l列;
步骤2.43:交换R,Q,p,norm的2j-1与2k-1项,交换R,Q,p的2j与2k项;
步骤2.44:按照式(2)计算;
其中,R2j-1,2j-1表示矩阵R的第2j-1行第2j-1列的元素、normj表示向量norm的第j列、Q2j-1表示矩阵Q的第2j-1列、R2j,2j表示矩阵R的第2j行第2j列的元素、R2j-1,2j表示矩阵R的第2j-1行第2j列的元素、Q2j表示矩阵Q的第2j列、QM+1:2M,2j-1 T表示对矩阵Q第2j-1列中第M+1行到2M行的元素所组成的向量做转置、Q1:M,2j-1 T表示对矩阵Q第2j-1列中第1行到M行的元素所组成的向量做转置;
R2j,2j与R2j-1,2j计算过程如下:
因为Q2j-1与Q2j正交,所以||Q2j-1||与||Q2j||相等,R2j-1,2j-1=||Q2j-1||,将Q2j-1的标准化值映射到Q2j上得到R2j-1,2j,R2j-1,2j表示为:
步骤2.45:设置迭代次数N-j,令第四迭代计数器n=j+1;
步骤2.46:按照式(3)计算
其中R2j,2n,R2j,2n-1与Q2n计算过程如下:
R2j-1,2n-1的计算过程同R2j,2n,Q2n=Q2n-R2j-1,2nQ2j-1-R2j,2nQ2j,由以上结果得出:Q2n与Q2n-1对称,Q2n=[-QM+1:2M,2n-1 T,Q1:M,2n-1 T]T;
其中,R2j-1,2n-1表示矩阵R的第2j-1行第2n-1列的元素、表示对矩阵Q的第2j-1列做转置、Q2n-1表示矩阵Q的第2n-1列、R2j-1,2n表示矩阵R的第2j-1行第2n列的元素、Q2n表示矩阵Q的第2n列、R2j,2n-1表示矩阵R的第2j行第2n-1列的元素、R2j,2n表示矩阵R的第2j行第2n列的元素、normn表示向量norm的第n列、QM+1:2M,2n-1 T表示对矩阵Q第2n-1列中第M+1行到2M行的元素所组成的向量做转置、Q1:M,2n-1 T表示对矩阵Q第2n-1列中第1行到M行的元素所组成的向量做转置、Q1:M,2j T表示对矩阵Q第2j列中第1行到M行的元素所组成的向量做转置、QM+1:2M,2j T表示对矩阵Q第2j列中第1行到M行的元素所组成的向量做转置、Q1:M,2n T表示对矩阵Q第2n列中第1行到M行的元素所组成的向量做转置、QM+1:2M,2n T表示对矩阵Q第2n列中第1行到M行的元素所组成的向量做转置;
步骤2.47:令n=n+1,并返回至步骤2.45,直至迭代达到预设次数N-j;
步骤2.48:令j=j+1,并返回至步骤2.41,直至迭代达到预设次数N为止,则得到正交矩阵Q和上三角矩阵R。
进一步的,所述接收信号向量实数化处理过程包括:
采用公式(5)表示接收信号向量:
r=Gs+v (5)
其中r为接收信号向量,s为发送信号向量,v为加性高斯白噪声向量;
对公式(5)等式两边的参数分别进行实数化,得到:
y=Hx+w (6)
y表示实数化接受信号向量,x表示实数化发送信号向量,w表示实数化加性高斯白噪声向量,即:
进一步的,所述接收信号向量检测过程包括如下步骤:
步骤3.1:所述实数化接受信号向量表示为:
y=Hx+w=QRx+w (8)
在等式两边分别乘上QH,获得检测向量
其中z=QHw;
步骤3.2:计算信号估计值。
进一步的,所述信号估计值的计算包括如下步骤:
步骤3.21:设置迭代总次数为N,第五迭代计数器a=N;
步骤3.22:设置迭代总次数为N-(a-1),第六迭代计数器b=a+1;
步骤3.23:按照下式计算干扰总和:
其中Ra,a表示矩阵R的第a行第a列的元素;表示第b个发送信号的估计值;
步骤3.24:令b=b+1,返回步骤3.22,至迭代达到预设次数N-(a-1)为止;
步骤3.25:按照下式计算信号估计值:
其中SGN[·]表示对信号进行硬判决;表示被检测向量;表示第a个发送信号的估计值;
步骤3.26:令a=a+1,并返回至步骤3.21,直至迭代达到预设次数N为止。
与现有技术相比,本发明具有如下的有益效果:
本发明改进的Gram-Schmidt(MGS)过程的每一步都有两个列向量正交,对其中一个列向量进行对称处理,可以有效地消除其中一个列向量正交所需要的计算量;利用列向量之间的对称性,将矩阵Q中的一对相邻列与另一对相邻列交换,减少了现有排序QR分解算法中的正交计算,在系统性能没有明显降低的前提下,降低了现有算法的复杂度,通过该算法本发明系统信号检测方法的计算复杂度几乎是传统的50%。
附图说明
图1是用户数为2且基站端天线数和用户端天线数均为16时,本发明方法的检测结果与理论最优检测方法最大合并比(MRC)、未排序的QR检测和已排序的QR检测的误码率表现比较图。
具体实施方式
下面对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
一种低复杂度的非正交多址接入系统信号检测方法,所述方法包括以下步骤:
步骤1:获取加性高斯白噪声信道的复信道系数矩阵、接收信号向量、发送信号向量、发送天线数和接收天线数,对所述复信道系数矩阵进行实数化处理,得到实数矩阵,所述实数矩阵的行数和列数分别是复信道系数矩阵的两倍;
所述复信道系数矩阵采用公式(1)进行实数化处理:
式中:G为复信道系数矩阵;H为实数矩阵。
步骤2:根据发送天线数、接收天线数对所述实数矩阵进行对称的排序QR分解,得到一个正交矩阵Q和一个上三角矩阵R;
步骤2.1:令Q=H,R=0,排列向量p=[1,(1+C),2,(2+C),...,C,2C],其中C表示实数矩阵H的列数,C的大小与发送天线数相等;
步骤2.2:将矩阵Q中的矢量按照向量p重新排列,此时矩阵Q中相邻两列向量对称,即Q2i=[-QM+1:2M,2i-1 T,Q1:M,2i-1 T]T,且Q2i-1与Q2i正交,
其中M表示矩阵Q的列数,M的大小与接收天线数相等,i表示第一迭代计数器,(·)T表示矩阵转置;Q2i表示矩阵Q的第2i列;QM+1:2M,2i-1表示取矩阵Q第2i-1列中第M+1行到2M行的元素;Q1:M,2i-1表示取矩阵Q第2i-1列中第1行到M行的元素;Q2i-1表示矩阵Q的第2i-1列;
步骤2.3:计算Q中奇数列的二范数模值的平方,记为 其中(·)T表示矩阵转置,表示对矩阵Q的第2i-1列做转置;N表示总次数,与发送天线数相等;
步骤2.4:计算矩阵Q和矩阵R中的元素。
步骤2.41:设置迭代总次数N,令第二迭代计数器j=0;
步骤2.42:k=argminl=j,...,Nnorml,其中k表示二范数模值平方最小向量的坐标,l表示第三迭代计数器,argmin()表示得到最小值的坐标;norml表示向量norm的第l列;
步骤2.43:交换R,Q,p,norm的2j-1与2k-1项,交换R,Q,p的2j与2k项;
步骤2.44:按照式(2)计算;
其中,R2j-1,2j-1表示矩阵R的第2j-1行第2j-1列的元素;normj表示表示向量norm的第j列;(normj)1/2表示对normj开方;Q2j-1表示矩阵Q的第2j-1列;R2j,2j表示矩阵R的第2j行第2j列的元素;R2j-1,2j表示矩阵R的第2j-1行第2j列的元素;Q2j表示矩阵Q的第2j列;QM+1:2M,2j-1 T表示对矩阵Q第2j-1列中第M+1行到2M行的元素所组成的向量做转置;Q1:M,2j-1 T表示对矩阵Q第2j-1列中第1行到M行的元素所组成的向量做转置。
R2j,2j与R2j-1,2j计算过程如下:
因为Q2j-1与Q2j正交,所以||Q2j-1||与||Q2j||相等,R2j-1,2j-1=||Q2j-1||,将Q2j-1的标准化值映射到Q2j上得到R2j-1,2j,所以R2j-1,2j可以表示为:
步骤2.45:设置迭代次数N-j,令第四迭代计数器n=j+1;
步骤2.46:按照式(3)计算
其中R2j,2n,R2j,2n-1与Q2n计算过程如下:
R2j-1,2n-1的计算过程同R2j,2n,Q2n=Q2n-R2j-1,2nQ2j-1-R2j,2nQ2j,由以上结果可得Q2n与Q2n-1对称,所以Q2n=[-QM+1:2M,2n-1 T,Q1:M,2n-1 T]T;其中,R2j-1,2n-1表示矩阵R的第2j-1行第2n-1列的元素;表示对矩阵Q的第2j-1列做转置;Q2n-1表示矩阵Q的第2n-1列;R2j-1,2n表示矩阵R的第2j-1行第2n列的元素;Q2n表示矩阵Q的第2n列;R2j,2n-1表示矩阵R的第2j行第2n-1列的元素;R2j,2n表示矩阵R的第2j行第2n列的元素;normn表示向量norm的第n列;QM+1:2M,2n-1 T表示对矩阵Q第2n-1列中第M+1行到2M行的元素所组成的向量做转置;Q1:M,2n-1 T表示对矩阵Q第2n-1列中第1行到M行的元素所组成的向量做转置;Q1:M,2j T表示对矩阵Q第2j列中第1行到M行的元素所组成的向量做转置;QM+1:2M,2j T表示对矩阵Q第2j列中第1行到M行的元素所组成的向量做转置;Q1:M,2n T表示对矩阵Q第2n列中第1行到M行的元素所组成的向量做转置;QM+1:2M,2n T表示对矩阵Q第2n列中第1行到M行的元素所组成的向量做转置。
步骤2.47:令n=n+1,并返回至步骤2.45,直至迭代达到预设次数N-j;
步骤2.48:令j=j+1,并返回至步骤2.41,直至迭代达到预设次数N为止,则得到正交矩阵Q和上三角矩阵R。
步骤3:根据所述发送信号向量对所述接收信号向量进行实数化处理、检测,得到发射信号估计值。
所述接收信号向量实数化处理过程包括:
采用公式(5)表示接收信号向量:
r=Gs+v (5)
其中r为接收信号向量,s为发送信号向量,v为加性高斯白噪声向量;
其中r为维数为M×1的接收信号向量,G为维数为M×C的矩阵,s为维数为C×1的发送信号向量,v为维数为M×1的加性高斯白噪声向量,其每一项均服从均值为0,方差为σ2,其协方差矩阵为E{nHn}=N0I,N0为噪声的单边带功谱密度,I为单位矩阵,且M与接收天线数相等,C与发送天线数相等。
对公式(5)等式两边的参数分别进行实数化,得到:
y=Hx+w (6)
y表示实数化接受信号向量,x表示实数化发送信号向量,w表示实数化加性高斯白噪声向量,即:
所述接收信号向量检测过程包括如下步骤:
步骤3.1:所述实数化接受信号向量表示为:
y=Hx+w=QRx+w (8)
在等式两边分别乘上QH,检测向量可以表示为
其中(·)H表示对矩阵做共轭转置,其中z=QHw;
步骤3.2:计算信号估计值;
步骤3.21:设置迭代总次数为N,第五迭代计数器为a=N;
步骤3.22:设置迭代总次数为N-(a-1),第六迭代计数器为b=a+1;
步骤3.23:按照下式计算干扰总和:
其中Ra,b表示矩阵R的第a行第b列的元素;表示第b个发送信号的估计值;
步骤3.24:令b=b+1,返回步骤3.22,至迭代达到预设次数N-(a-1)为止;
步骤3.25:按照下式计算信号估计值:
其中SGN[·]表示对信号进行硬判决,表示第a个发送信号的估计值,Ra,a表示矩阵R的第a行第a列的元素;
步骤3.26:令a=a+1,并返回至步骤3.21,直至迭代达到预设次数N为止。
最后,对本发明的方法与经典排序QR算法进行对比:
表1
通过计算算法所需的数值运算的数量来分析传统算法和本发明提出的方法中算法计算复杂度,结果如表1所示。两种算法都基于MGS正交化过程进行排序QR检测;然而,由于每一步都有两个列向量正交化,因此在该算法中MGS过程的步数减半。此外,利用两列向量之间的对称性,可以完全消除其中一个列向量正交所需要的计算。因此,每个步骤所涉及的计算复杂度完全没有增加。结果表明,与传统算法相比,该方法中的算法综合计算复杂度降低了近50%,通过该算法本发明系统信号检测方法的计算复杂度几乎是传统的50%。
通过Matlab仿真来分析传统算法与本发明方法中的算法的性能,结果如图1所示。图1表明对称SQRD算法在相同信噪比情况下性能优于理论最优算法MRC和不排序的QRD算法,即误比特率较低。与SQRD算法相比,在低信噪比情况下误比特率低于SQRD算法,在高信噪比情况下,误比特率没有明显升高。因此,本发明在性能降低较少的前提下,通过该算法本发明系统信号检测方法的计算复杂度降低了50%。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (7)
1.一种低复杂度的非正交多址接入系统信号检测方法,其特征在于,所述方法包括以下步骤:
获取加性高斯白噪声信道的复信道系数矩阵、接收信号向量、发送信号向量、发送天线数和接收天线数,对所述复信道系数矩阵进行实数化处理,得到实数矩阵,所述实数矩阵的行数和列数分别是复信道系数矩阵的两倍;
根据发送天线数、接收天线数对所述实数矩阵进行对称的排序QR分解,得到一个正交矩阵Q和一个上三角矩阵R;
根据所述发送信号向量对所述接收信号向量进行实数化处理、检测,得到发射信号估计值。
2.根据权利要求1所述的一种低复杂度的非正交多址接入系统信号检测方法,其特征在于,所述复信道系数矩阵采用公式(1)进行实数化处理:
式中:G为复信道系数矩阵;H为实数矩阵。
3.根据权利要求2所述的一种低复杂度的非正交多址接入系统信号检测方法,其特征在于,所述排序QR的分解方法包括如下步骤:
令Q=H,R=0,排列向量p=[1,(1+C),2,(2+C),...,C,2C],其中C表示实数矩阵H的列数,C的大小与发送天线数相等;
将矩阵Q中的矢量按照向量p重新排列,此时矩阵Q中相邻两列向量对称,即Q2i=[-QM+1:2M,2i-1 T,Q1:M,2i-1 T]T,且Q2i-1与Q2i正交,
其中M表示矩阵Q的列数,M的大小与接收天线数相等,(·)T表示矩阵转置,Q2i表示矩阵Q的第2i列,QM+1:2M,2i-1表示取矩阵Q第2i-1列中第M+1行到2M行的元素,Q1:M,2i-1表示取矩阵Q第2i-1列中第1行到M行的元素,Q2i-1表示矩阵Q的第2i-1列;
计算Q中奇数列的二范数模值的平方,记为 其中表示对矩阵Q的第2i-1列做转置;i表示第一迭代计数器,N表示迭代总次数,与发送天线数相等;
计算矩阵Q和矩阵R中的元素。
4.根据权利要求3所述的一种低复杂度的非正交多址接入系统信号检测方法,其特征在于,计算矩阵Q和矩阵R中的元素过程包括:
步骤2.41:设置迭代总次数N,令第二迭代计数器j=0;
步骤2.42:k=argminl=j,...,Nnorml,其中k表示二范数模值平方最小向量的坐标,l表示第三迭代计数器,norml表示向量norm的第l列;
步骤2.43:交换R,Q,p,norm的2j-1与2k-1项,交换R,Q,p的2j与2k项;
步骤2.44:按照式(2)计算;
其中,R2j-1,2j-1表示矩阵R的第2j-1行第2j-1列的元素、normj表示向量norm的第j列、Q2j-1表示矩阵Q的第2j-1列、R2j,2j表示矩阵R的第2j行第2j列的元素、R2j-1,2j表示矩阵R的第2j-1行第2j列的元素、Q2j表示矩阵Q的第2j列、QM+1:2M,2j-1 T表示对矩阵Q第2j-1列中第M+1行到2M行的元素所组成的向量做转置、Q1:M,2j-1 T表示对矩阵Q第2j-1列中第1行到M行的元素所组成的向量做转置;
R2j,2j与R2j-1,2j计算过程如下:
因为Q2j-1与Q2j正交,所以||Q2j-1||与||Q2j||相等,R2j-1,2j-1=||Q2j-1||,将Q2j-1的标准化值映射到Q2j上得到R2j-1,2j,R2j-1,2j表示为:
步骤2.45:设置迭代次数N-j,令第四迭代计数器n=j+1;
步骤2.46:按照式(3)计算
其中R2j,2n,R2j,2n-1与Q2n计算过程如下:
R2j-1,2n-1的计算过程同R2j,2n,Q2n=Q2n-R2j-1,2nQ2j-1-R2j,2nQ2j,由以上结果得出:Q2n与Q2n-1对称,Q2n=[-QM+1:2M,2n-1 T,Q1:M,2n-1 T]T;
其中,R2j-1,2n-1表示矩阵R的第2j-1行第2n-1列的元素、表示对矩阵Q的第2j-1列做转置、Q2n-1表示矩阵Q的第2n-1列、R2j-1,2n表示矩阵R的第2j-1行第2n列的元素、Q2n表示矩阵Q的第2n列、R2j,2n-1表示矩阵R的第2j行第2n-1列的元素、R2j,2n表示矩阵R的第2j行第2n列的元素、normn表示向量norm的第n列、QM+1:2M,2n-1 T表示对矩阵Q第2n-1列中第M+1行到2M行的元素所组成的向量做转置、Q1:M,2n-1 T表示对矩阵Q第2n-1列中第1行到M行的元素所组成的向量做转置、Q1:M,2j T表示对矩阵Q第2j列中第1行到M行的元素所组成的向量做转置、QM+1:2M,2j T表示对矩阵Q第2j列中第1行到M行的元素所组成的向量做转置、Q1:M,2n T表示对矩阵Q第2n列中第1行到M行的元素所组成的向量做转置、QM+1:2M,2n T表示对矩阵Q第2n列中第1行到M行的元素所组成的向量做转置;
步骤2.47:令n=n+1,并返回至步骤2.45,直至迭代达到预设次数N-j;
步骤2.48:令j=j+1,并返回至步骤2.41,直至迭代达到预设次数N为止,则得到正交矩阵Q和上三角矩阵R。
5.根据权利要求4所述的一种低复杂度的非正交多址接入系统信号检测方法,其特征在于,所述接收信号向量实数化处理过程包括:
采用公式(5)表示接收信号向量:
r=Gs+v (5)
其中r为接收信号向量,s为发送信号向量,v为加性高斯白噪声向量;
对公式(5)等式两边的参数分别进行实数化,得到:
y=Hx+w (6)
y表示实数化接受信号向量,x表示实数化发送信号向量,w表示实数化加性高斯白噪声向量,即:
6.根据权利要求5所述的一种低复杂度的非正交多址接入系统信号检测方法,其特征在于,所述接收信号向量检测过程包括如下步骤:
步骤3.1:所述实数化接受信号向量表示为:
y=Hx+w=QRx+w (8)
在等式两边分别乘上QH,获得检测向量
其中z=QHw,(·)H表示对矩阵做共轭转置;
步骤3.2:计算信号估计值。
7.根据权利要求6所述的一种低复杂度的非正交多址接入系统信号检测方法,其特征在于,所述信号估计值的计算包括如下步骤:
步骤3.21:设置迭代总次数为N,第五迭代计数器a=N;
步骤3.22:设置迭代总次数为N-(a-1),第六迭代计数器b=a+1;
步骤3.23:按照下式计算干扰总和:
其中Ra,a表示矩阵R的第a行第a列的元素;表示第b个发送信号的估计值;
步骤3.24:令b=b+1,返回步骤3.22,至迭代达到预设次数N-(a-1)为止;
步骤3.25:按照下式计算信号估计值:
其中SGN[·]表示对信号进行硬判决;表示被检测向量;表示第a个发送信号的估计值;
步骤3.26:令a=a+1,并返回至步骤3.21,直至迭代达到预设次数N为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910055664.0A CN109639607B (zh) | 2019-01-21 | 2019-01-21 | 一种低复杂度的非正交多址接入系统信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910055664.0A CN109639607B (zh) | 2019-01-21 | 2019-01-21 | 一种低复杂度的非正交多址接入系统信号检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109639607A true CN109639607A (zh) | 2019-04-16 |
CN109639607B CN109639607B (zh) | 2021-09-03 |
Family
ID=66062264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910055664.0A Active CN109639607B (zh) | 2019-01-21 | 2019-01-21 | 一种低复杂度的非正交多址接入系统信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109639607B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113141233A (zh) * | 2021-03-11 | 2021-07-20 | 西安电子科技大学 | 一种信道矩阵排序qr分解处理方法、系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105337699A (zh) * | 2015-11-19 | 2016-02-17 | 北京邮电大学 | 一种应用于非正交多址接入系统的信号检测方法及装置 |
CN105978835A (zh) * | 2016-05-06 | 2016-09-28 | 华东交通大学 | 一种低复杂度的分布式干扰对齐方法 |
CN106411796A (zh) * | 2015-07-31 | 2017-02-15 | 电信科学技术研究院 | 一种非正交多址接入中多终端信号检测方法及基站 |
CN107592183A (zh) * | 2017-09-19 | 2018-01-16 | 北京邮电大学 | 一种应用于非正交多址系统的多用户信号检测方法及装置 |
CN108768477A (zh) * | 2018-05-24 | 2018-11-06 | 天津工业大学 | 一种基于模型扰动与格基约减的mimo检测方法 |
WO2018203627A1 (ko) * | 2017-05-02 | 2018-11-08 | 엘지전자(주) | 무선 통신 시스템에서 신호를 송수신하는 방법 및 이를 위한 장치 |
-
2019
- 2019-01-21 CN CN201910055664.0A patent/CN109639607B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106411796A (zh) * | 2015-07-31 | 2017-02-15 | 电信科学技术研究院 | 一种非正交多址接入中多终端信号检测方法及基站 |
CN105337699A (zh) * | 2015-11-19 | 2016-02-17 | 北京邮电大学 | 一种应用于非正交多址接入系统的信号检测方法及装置 |
CN105978835A (zh) * | 2016-05-06 | 2016-09-28 | 华东交通大学 | 一种低复杂度的分布式干扰对齐方法 |
WO2018203627A1 (ko) * | 2017-05-02 | 2018-11-08 | 엘지전자(주) | 무선 통신 시스템에서 신호를 송수신하는 방법 및 이를 위한 장치 |
CN107592183A (zh) * | 2017-09-19 | 2018-01-16 | 北京邮电大学 | 一种应用于非正交多址系统的多用户信号检测方法及装置 |
CN108768477A (zh) * | 2018-05-24 | 2018-11-06 | 天津工业大学 | 一种基于模型扰动与格基约减的mimo检测方法 |
Non-Patent Citations (2)
Title |
---|
FAN JIANG, YAN ZHANG, AND CHENG LI: ""A New SQRD-Based Soft Interference Cancelation Scheme in Multi-User MIMO SC-FDMA System"", 《IEEE COMMUNICATIONS LETTERS》 * |
张晓然等: ""一种用于检测中的快速分解算法"", 《信号处理》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113141233A (zh) * | 2021-03-11 | 2021-07-20 | 西安电子科技大学 | 一种信道矩阵排序qr分解处理方法、系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109639607B (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105915477B (zh) | 基于gs方法的大规模mimo检测方法及硬件架构 | |
CN110086555B (zh) | 一种大规模mimo系统中的分组导频分配方法及其分配装置 | |
CN107210783A (zh) | 采用协作中继进行联合mimo传输和压缩来降低干扰的系统和方法 | |
CN113517941B (zh) | 一种大规模mimo系统信道估计与迭代检测仿真方法和系统 | |
CN106028451A (zh) | 一种应用于noma中的用户分组机制 | |
CN109743086A (zh) | 一种大规模mimo系统的信道估计方法 | |
CN109347532A (zh) | 广义空间调制系统中改进的gomp检测算法 | |
Dong et al. | Improved joint antenna selection and user scheduling for massive MIMO systems | |
CN109995403A (zh) | 大规模mimo系统中基于模拟退火思想改进的las检测算法 | |
Tu et al. | An efficient massive MIMO detector based on second-order Richardson iteration: From algorithm to flexible architecture | |
CN105634568A (zh) | 一种基于大规模mimo系统信号检测的llr计算方法 | |
CN109818887B (zh) | 基于evd-ilsp的半盲信道估计方法 | |
CN109639607A (zh) | 一种低复杂度的非正交多址接入系统信号检测方法 | |
CN104994041B (zh) | 三维大规模天线系统基于导频的波达波离角联合估计方法 | |
CN109981151A (zh) | 大规模mimo系统中改进的高斯树近似消息传递检测算法 | |
CN109617577A (zh) | 一种基于压缩感知信号检测的无线光空间调制方法 | |
CN107017923A (zh) | 一种基于无线输能的非正交多址接入通信系统 | |
CN108736934A (zh) | 一种高效的大规模mimo系统信号检测方法 | |
CN107181705A (zh) | 一种半盲信道估计方法及系统 | |
CN114389756B (zh) | 基于分组ml检测和并行迭代干扰抵消的上行mimo检测方法 | |
CN106789781A (zh) | 基于Givens变换块对角化预编码的干扰消除方法 | |
CN106027434A (zh) | 一种基于csm的预编码方法 | |
CN107346985B (zh) | 一种结合发射天线选择技术的干扰对齐方法 | |
CN105429687A (zh) | 一种最小化干扰功率与维度的干扰对齐方法 | |
CN100373811C (zh) | 一种分布式多入多出系统的信号检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |