CN1901388B - 一种低复杂度的线性迭代多用户检测方法 - Google Patents
一种低复杂度的线性迭代多用户检测方法 Download PDFInfo
- Publication number
- CN1901388B CN1901388B CN200510035893A CN200510035893A CN1901388B CN 1901388 B CN1901388 B CN 1901388B CN 200510035893 A CN200510035893 A CN 200510035893A CN 200510035893 A CN200510035893 A CN 200510035893A CN 1901388 B CN1901388 B CN 1901388B
- Authority
- CN
- China
- Prior art keywords
- matrix
- centerdot
- detection method
- linear
- sgn
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 46
- 238000000034 method Methods 0.000 claims abstract description 7
- 230000006978 adaptation Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 6
- 238000004422 calculation algorithm Methods 0.000 description 29
- 238000004891 communication Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 230000005764 inhibitory process Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Landscapes
- Mobile Radio Communication Systems (AREA)
- Noise Elimination (AREA)
Abstract
本发明提供了一种低复杂度的线性迭代多用户检测方法,该方法通过一系列循环迭代简单地避免了对矩阵求逆的计算。由于该线性循环叠代方法减少了计算中的舍入误差,从而在大大降低数据计算量的同时,提高了计算的精度,降低了由多址干扰(MAI)造成的噪声,最终很好地恢复了多用户信息。
Description
技术领域:
本发明涉及一种低复杂度的线性迭代多用户检测方法,属于使用电磁波技术的数字通信领域,尤其涉及直接序列码分多址(DS-CDMA)***中多用户检测的装置和方法,这些装置和方法能实时清除存在于DS-CDMA***中的多址干扰(MAI)。
背景技术:
在CDMA通信***中,制约***性能的主要因素包括同信道干扰、码间干扰等。码间干扰约束了***通信的速率,主要采用均衡或分集技术来抑制;而主要由同小区的移动台、相邻小区的同频基站或移动台以及其他通信***等多址干扰(MAI)产生的同信道干扰则限制了***通信的容量。目前,消除多址干扰的技术是多用户检测技术,该技术在传统检测技术基础上,充分利用了产生多址干扰的所有用户信号信息,从而具有优良的抗干扰性能,降低了***对功率控制精度的要求,显著地提高了***容量。
多用户检测算法主要分为四类:最优检测算法、线性检测算法、干扰对消检测算法、盲检测算法。最优多用户检测算法基于最大似然估计准则,虽然性能最佳,但是运算复杂度太高;相应的改进算法确实也在一定程度上降低了计算量,但仍难于实用。干扰对消检测算法主要分为串行干扰抵消算法和并行干扰抵消算法,其基本思想是通过再造部分多址干扰(PMAI)并抵消从而达到抑制MAI的目的,干扰抵消器的性能依赖于再造PMAI的准确度和干扰抵消权值。在无法改善信道估计准确度时,部分干扰抵消接收机的性能在很大程度上受到干扰抵消权值的约束。盲检测器主要包括基于最小能量输出(MOE)准则的最小均方/递归二乘(RLS)/QR-RLS盲检测器、基于信道子空间搜索的盲检测器、基于窄带***和阵列处理的盲检测器、混合型半盲检测器四类。盲检测器的主要思想是通过子空间跟踪技术获得信号子空间并利用它来消除未知用户造成的干扰。一般而言,盲检测器在AWGN信道下具有优良的检测性能。性能次优的线性检测算法主要分为解相关检测算法和线性最小均方误差(LMMSE)检测算法两大类。解相关算法能完全消除MAI,但对噪声有所放大。LMMSE算法降低了线性映射对噪声的放大,但MAI消除不完全,同时必须对接收机信号的幅度进行估计。
影响这两种线性多用户检测算法成为实用的关键原因是矩阵求逆的问题。一方面,在多变的移动环境中,求解相关矩阵R的逆阵不易。尤其是在CDMA***中,为了更好地减少干扰,一般都要采用语音激活和变速率编码技术,这使得用户信号之间的相关随时变化,因而R的逆阵也是时变的,这使得实时地求解R的逆阵变得非常困难。另一方面,求逆矩阵的最快速算法的乘法运算量也是K3阶,当用户数K较多时,它的计算复杂度更是难以接受的。
发明内容:
本发明的目的是:针对现有技术的不足,提供一种计算复杂度更低的、新的适用于码分多址***的快速线性迭代多用户检测方法。
为实现上述目的,本发明的技术方案是:采用线性迭代的方法取代现有技术中的对相关矩阵所做的求逆运算,从而将计算复杂度大大降低。该算法通过简单的循环叠代快速地实现了对多址干扰的抑制。
下面介绍现有线性多用户检测方法与本发明的原理:
1、现有线性多用户检测方法
相关器输出离散信号的矩阵表示为:
y=RAb+N (1)
其中,R=S·SH,相关矩阵R的大小为K×K,它代表两个码字之间的相关值;S为一矩阵,大小为K×K,由R计算可得;b为传输的比特流矩阵,A为用户的信号幅度矩阵,N为噪声矩阵,大小均为K×K。
常规线性多用户检测算法包括常规多用户检测算法、解相关算法和最小均方估计多用户检测算法三种,均可同一表示为:
其中,L是K×K矩阵。
在常规多用户检测算法中,假设I为K×K大小的单位阵,有L=I,此时:
在解相关检测器中,采用逆矩阵L=R-1相乘以解调出用户信号,有:
而常规的最小均方估计检测器采用逆矩阵L=(R+σ2(A2)-1)-1,有:
式中,σ2为噪声的功率。
2、快速线性叠代多用户检测算法
常规线性多用户检测器须对相关阵R求逆,这造成了算法计算复杂度高,没法实现实时性。针对这一缺点,提出了快速线性叠代多用户检测算法,该算法不仅计算复杂度低,而且性能更优。
假设α为自适应调整因子,可以为1.5、2.0、4.0中任意值,该算法采用逆矩阵L=(R+ασ2(A2)-1)-1相乘以解调出用户信号,有:
对于相关阵R∈CK×K且非奇异,则(R+ασ2(A2)-1)可写成为:
R+ασ2(A2)-1=E+A (7)
其中E为单位阵。因此,有:
其中,(E+A)*表示矩阵(E+A)的伴随矩阵。
假设
(E+A)*=P(K)+P(K-1)+...+P(1)+P(1) (9)
式中,P(K)、P(K-1)、...、P(2)、P(1)为K阶常数矩阵。
同时,假设
|E+A|=1+c[K]+c[K-1]+...+c[2]+c[1] (10)
式中,c[K]、c[K-1]、...、c[2]、c[1]为常数。
定义n阶方阵A主对角线上元素的和为矩阵A的迹,记为tr(A)。
根据矩阵之间的关系,有更新循环叠代公式:
式中,tr(·)表示矩阵的迹;其中初值满足:
P(K)=E (13)
AP(1)=c[1]E
于是,有:
此时,
式中P(i)、c[i]、由更新叠代公式(11)、(12)、(13)获得。显然,快速线性叠代多用户检测算法不须计算相关矩阵R的逆阵。
本发明的具体操作步骤如下:
步骤1:在接收端,可获得k个匹配滤波器组输出的向量:
y=RBb+N (16)
其中,R=S·SH,相关矩阵R的大小为K×K,它代表两个码字之间的相关值。
步骤2:对于最小均方估计多用户检测器,选定自适应调整因子α,实测噪声功率σ2和信号幅度B。计算出:
R+ασ2(B2)-1=E+A (17)
获得矩阵A。其中,E为K×K的单位阵。
步骤3:设定初值:
P(K)=E (18)
AP(1)=c[1]E (19)
步骤4:根据矩阵之间关系,利用更新循环叠代公式:
获得P(i)和c[i](i=1,2,...,K)的值。
步骤5:利用最小均方估计检测器可获得信号:
本发明在CDMA通信***的接收端通过一系列线性迭代运算取代对相关矩阵的求逆运算,尤其是在用户较多时,大大降低了计算复杂度,达到了快速实现的目的。
同时,经过CDMA通信***链路仿真也证明,与其他常规方法相比,该算法采用循环叠代求逆将运算复杂度大大降低,而且循环迭代减少了计算中的舍入误差,从而使误比特率和抗远近效应的性能均更优。
附图说明:
图1是说明根据本发明的CDMA***中多用户检测装置的框图。
图2、图3均为用户1的误比特率情况。如图所示,改进的快速线性叠代检测器的性能比解相关检测器和最小均方估计检测器误比特率性能更佳,具有更强的抵消多址干扰的能力。通过图3可以看出,当α取0.5和0.9时,其误码率均比最小均方估计E检测器和解相关检测器性能更佳;而且当α取2.0和4.0时,其误码率均比α取0.5和0.9时性能更佳,显然,α取1.5、2.0、4.0中的一个值就可。
图4为用户1的误比特率随“远近”效应的变化情况。如图所示,当α取1.5、2.0、4.0时,本专利方法的抗“远近”效应性能均比0.9时更佳。所以,α取1.5、20、4.0中的一个值就可,改进的快速线性叠代检测器性能比解相关检测器和最小均方估计检测器具有更佳的抗“远近”效应性能。
具体实施方式:
下面通过具体的实施对本发明的技术方案作进一步的描述。
采用Matlab仿真软件仿真,仿真***为同步/异步DS-CDMA***,AWGN信道,BPSK调制,载波相位为0,用长度为10的伪随机序列扩频,对修改前后的解相关多用户检测器分别从误比特率、抗“远近”效应和计算复杂度三个方面进行了分析。假设用户1为期望用户,用户数为10。
步骤1:在接收端,可获得10个匹配滤波器组输出的向量:
y=RBb+N
其中,R=S·SH,相关矩阵R的大小为10×10,它代表两个码字之间的相关值。
步骤2:对于最小均方估计多用户检测器,选定自适应调整因子α,实测噪声功率σ2和信号幅度B。计算出:
R+ασ2(B2)-1=E+A
获得矩阵A。其中,E为K×K的单位阵。
步骤3:设定初值:
P(K)=E
AP(1)=c[1]E
步骤4:根据矩阵之间关系,利用更新循环叠代公式:
获得P(i)和c[i](i=1,2,...,K)的值。
步骤5:线性叠代多用户检测器可获得信号:
Claims (3)
1.一种低复杂度的线性迭代多用户检测方法,其特征在于包括如下步骤:
步骤1在接收端,可获得i个匹配滤波器组输出的向量:
y=RBb+N (1)
其中,R=S·SH,相关矩阵R的大小为K×K,它代表两个码字之间的相关值;S为一矩阵,大小为K×K,由R计算可得;b为传输的比特流矩阵,A为用户的信号幅度矩阵,N为噪声矩阵,大小均为K×K;
步骤2:对于线性迭代多用户检测器,选定自适应调整因子α为1.5、2.0、4.0中任意值,实测噪声功率σ2和信号幅度B,计算出:
R+ασ2(B2)-1=E+A (2)
获得矩阵A,其中,E为K×K的单位阵;
步骤3:设定初值:
P(K)=E (3)
AP(1)=c[1]E (4)
步骤4:根据矩阵之间关系,利用更新循环叠代公式:
获得P(i)和c[i](i=1,2,...,K)的值;式中,tr(·)表示矩阵的迹;
步骤5:进行信号检测:
由此,可获得信号。
2.如权利要求1所说,一种低复杂度的线性迭代多用户检测方法,其特征在于:该方法采用了简单的线性迭代取代了现有技术中对相关矩阵的求逆运算。
3.如权利要求1所说,一种低复杂度的线性迭代多用户检测方法,其特征在于:可用于DS-CDMA***或CDMA***中多用户检测去除多址干扰。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200510035893A CN1901388B (zh) | 2005-07-20 | 2005-07-20 | 一种低复杂度的线性迭代多用户检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200510035893A CN1901388B (zh) | 2005-07-20 | 2005-07-20 | 一种低复杂度的线性迭代多用户检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1901388A CN1901388A (zh) | 2007-01-24 |
CN1901388B true CN1901388B (zh) | 2010-05-05 |
Family
ID=37657159
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200510035893A Expired - Fee Related CN1901388B (zh) | 2005-07-20 | 2005-07-20 | 一种低复杂度的线性迭代多用户检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1901388B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101807939B (zh) * | 2009-02-12 | 2013-10-09 | 中国电信股份有限公司 | 多用户检测方法和*** |
CN106357318B (zh) * | 2016-10-31 | 2019-10-11 | 东南大学 | 收敛速率可调的大规模mimo迭代检测方法 |
CN108512580B (zh) * | 2018-03-02 | 2021-05-25 | 东南大学 | 适用于低精度量化的大规模多用户mimo迭代检测方法 |
-
2005
- 2005-07-20 CN CN200510035893A patent/CN1901388B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN1901388A (zh) | 2007-01-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101128992B (zh) | 使用并行剩余补偿的自适应加权干扰消除的***、装置和方法 | |
Strom et al. | Propagation delay estimation in asynchronous direct-sequence code-division multiple access systems | |
CN1902834B (zh) | 用于码分多址通信的方法、装置和*** | |
CN105791184B (zh) | 基于相似度的非周期长码直扩信号扩频序列盲估计方法 | |
CN101316115A (zh) | 导频序列信号的检测方法、设备及*** | |
CN101800568A (zh) | 自适应并行干扰抵消多用户检测器的算法及实现 | |
CN1615598B (zh) | 低复杂度多用户检测器及cdma接收机***中为用户产生解扩序列的方法 | |
CN1901388B (zh) | 一种低复杂度的线性迭代多用户检测方法 | |
CN101420248A (zh) | 一种td-scdma终端频偏估计的方法及装置 | |
Orsak et al. | Efficient importance sampling techniques for simulation of multiuser communication systems | |
CN101001092A (zh) | 低复杂度的解相关多用户检测方法与装置 | |
KR100909519B1 (ko) | 타임-슬롯 cdma 시스템에서의 다중-코드-셋 채널 추정 방법 | |
CN101222240A (zh) | 异步ds-cdma盲多用户检测的粒子滤波方法 | |
CN110995364B (zh) | 一种提升双差分扩频水声通信***通信速率的通信方法 | |
CN1992558B (zh) | 一种无线通信***以及无线通信方法 | |
Ertin et al. | Iterative techniques for DS/CDMA multipath channel estimation | |
CN106603201B (zh) | 一种基于采样点处理的多用户联合检测方法 | |
Laurenson et al. | A computationally efficient multipath channel simulator for the COST 207 models | |
CN101257324A (zh) | Td-scdma***中的线性联合信道估计方法 | |
CN1953343A (zh) | 线性联合检测输出生成信道译码器软比特的方法 | |
Sengupta et al. | Subspace‐based tracking of multipath channel parameters for CDMA systems | |
Jeney et al. | Blind adaptive stochastic neural network for multiuser detection | |
Linlin et al. | Performance analysis of Volterra-based nonlinear adaptive blind multiuser detectors for DS-CDMA systems | |
CN100539456C (zh) | 多用户联合检测方法 | |
Sengupta et al. | Fixed point error analysis of multiuser detection and synchronization algorithms for CDMA communication systems |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100505 Termination date: 20100720 |