CN1705242A - 一种分组垂直分层空时检测方法及装置 - Google Patents
一种分组垂直分层空时检测方法及装置 Download PDFInfo
- Publication number
- CN1705242A CN1705242A CN 200410042933 CN200410042933A CN1705242A CN 1705242 A CN1705242 A CN 1705242A CN 200410042933 CN200410042933 CN 200410042933 CN 200410042933 A CN200410042933 A CN 200410042933A CN 1705242 A CN1705242 A CN 1705242A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- symbol
- math
- centerdot
- 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
- 238000000034 method Methods 0.000 title claims abstract description 51
- 239000013598 vector Substances 0.000 claims abstract description 201
- 238000004364 calculation method Methods 0.000 claims abstract description 47
- 238000001514 detection method Methods 0.000 claims description 101
- 239000011159 matrix material Substances 0.000 claims description 63
- 230000008707 rearrangement Effects 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 abstract 6
- 238000010998 test method Methods 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 101000741965 Homo sapiens Inactive tyrosine-protein kinase PRAG1 Proteins 0.000 description 3
- 102100038659 Inactive tyrosine-protein kinase PRAG1 Human genes 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000013139 quantization Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 241001105470 Valenzuela Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012956 testing procedure Methods 0.000 description 1
Images
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种分组垂直分层空时检测方法包括:将接收信号矢量中待检测的符号分为Z组,Z>1,令i=1;A.检测接收信号矢量,计算出第i组内所有符号的估计值,消除该组所有符号的干扰,得到修正的接收信号矢量,再将修正的接收信号矢量作为需检测的接收信号矢量,并令i=i+1,执行步骤A,在执行Z-1次步骤A后,转入B;B.检测接收信号矢量,计算出第Z组内所有符号的估计值。本发明还公开了一种分组垂直分层空时检测装置。本发明将接收信号矢量中待检测的符号分成Z组,组与组间进行串行检测与串行干扰对消,在组内进行并行检测与并行干扰对消,降低系统时延,减少计算量,提高系统容量和速率,并且采用软判决方法得到符号的估计,提高符号估计的性能。
Description
技术领域
本发明涉及无线通信领域多输入多输出(MIMO,Multiple Input MultipleOutput)技术,特别是指一种MIMO系统的分组垂直分层空时检测方法及装置。
背景技术
MIMO技术是无线通信领域的重大技术突破,它能够在不增加带宽的情况下成倍地提高通信系统的容量和频谱利用率。MIMO技术在发送端和接收端采用多天线(天线阵列)同时发送和接收信号。由于各发射天线同时发送的信号占用同一个频带,因而通信带宽并没有增加。每个发送天线和每个接收天线之间存在一个空间信道。如果每个空间信道的信道冲击响应独立,则MIMO系统通过多个发送天线和多个接收天线可以在发送端和接收端之间创造多个并行的独立的空间信道。通过这些并行的空间信道独立地传输信息,MIMO系统的传输数据率必然成倍增加。1998年G.J.Foschini和M.J.Gans充分论证了上述结论,并定量指出:假设MIMO系统有M根发送天线和N根接收天线,在窄带慢衰落信道下,就可以建立N×M阶信道矩阵。该矩阵的元素为独立同分布的复高斯随机变量。MIMO系统可以获得的信道容量将是单输入单输出(SISO,Single Input Single Output)系统的min(m,n)倍,其中min(m,n)表示取m和n中的最小值,且总发射功率保持不变。
1996年G.J.Foshini曾经提出了MIMO系统的D-BLAST贝尔实验室对角分层空时结构(Diagonal BLAST,BLAST:Bell Laboratories LayeredSpace-time Architecture)方法,该方法可以实现MIMO系统的解调,并且可以获得接近理论容量90%的容量。但是,D-BLAST方法复杂度较大,不易实时实现。
1999年,G.D.Golden、G.J.Foshini、R.A.Valenzuela和P.W.Wolniansky提出一种简化的BLAST方法——贝尔实验室垂直分层空时结构(V-BLAST,Vertical BLAST)。该方法已经在实验室实时实现。实验结果表明:发送端和接收端分别采用8根天线和12根天线的MIMO系统在室内环境下,当平均SNR在24dB~34dB的范围内变化时,可以获得20-40bps/Hz的频谱效率。尽管该频谱效率是在室内环境下获得,但这个量级的频谱效率是空前的。
在V-BLAST方法中,假设M为发射天线数目,相应地接收机采用N根天线进行接收,被发送的符号矢量为:a=[a1,a2,…,aM]T,接收信号矢量为:
r1=Ha+v (1)
这里,H为N×M维信道矩阵,v表示N维噪声矢量,其分量是独立同分布的高斯白噪声,σ2为接收信号矢量中每个噪声分量的功率,σs 2是发送符号的功率,r1为N维接收信号矢量。
s≡[k1,k2,…,kM]T是整数1,2,......,M的一个排列,该排列表示符号矢量a中符号分量被检测的顺序。V-BLAST检测算法由r1按照s≡[k1,k2,…,kM]T给定的顺序计算决策统计量yk1,yk2,…,ykM。由决策统计量yk1,yk2,…,ykM可以形成数据符号的估计k1,k2,…,kM。
可以按照某种最优准则计算得到s≡[k1,k2,…,kM]T,这里将按照某种最优准则计算得到的最优排列记为SOPT≡[k1,k2,…,kM]T。下面阐述在给定的s≡[k1,k2,…,kM]T下V-BLAST算法的检测过程。
V-BLAST算法采用线性迫零和符号对消方法,串行地计算决策统计量。步骤如下:
步骤1:采用N维迫零权矢量wk1,通过计算r1分量的线性合并,得到决策统计量yk1:
步骤2:由yk1量化得到k1:
其中,Q(·)是量化函数,该函数与发送端调制方法相对应。a=[a1,a2,…,aM]T中每个符号是信息比特经过调制得到的。
步骤3:假设
从接收信号矢量r1中消除ak1,得到修正的接收信号矢量r2:
这里(H)k1表示矩阵H的第k1列。
按照上述步骤1~3对修正的接收信号矢量r2进行同样的计算,即:由wk2得到yk2,再由yk2得到ak2的估计ak2和r3。按照顺序k3,k4,…,kM计算下去,可以求出其他被发送符号的估计k3,k4,…,kM。
迫零权矢量wk1,wk2,…,wkM的计算是算法的关键。迫零权矢量的具体计算方法与所采用的计算准则有关。通常采用ZF准则和MMSE准则进行迫零权矢量的计算。在采用ZF准则下,矢量wk1,wk2,…,wkM满足如下关系:
因此,在采用ZF准则下,V-BLAST算法可以通过如下递归过程实现:
初始化:
G1=H+ (6a)
i=1 (6b)
进行如下计算:
ki=argmin‖(Gi)j‖2
j=i+1(6i)
递归计算公式(6c)~(6i),直到i=M+1时才停止计算。这时已经通过V-BLAST方法得到了所有被发送符号的估计:k1,k2,…,kM
上述公式中,+表示Moore-Penrose伪逆;令矩阵H中的第k1,k2,…,ki列为0列矢量,其他列矢量保持不变,就得到矩阵H
ki
。(Gi)ki表示矩阵Gi的第ki行。
上述V-BLAST方法的递归实现中,公式(6a)和(6h)是ZF准则下迫零矢量的计算公式。如果采用MMSE准则,公式(6a)和(6h)需要进行相应的更改。无论采用何种准则计算迫零矢量,V-BLAST方法的递归实现过程都如上所述,只是不同准则下公式(6a)和(6h)的表达式要进行相应的变换,与所采用的准则保持一致。
从上述V-BLAST方法的递归实现中可以看到:V-BLAST方法采用串行干扰对消的思想。即:检测一个符号就进行一次干扰对消;通过干扰对消对接收信号矢量进行一次修正,使修正的接收信号矢量不再包含已经检测的符号的信息;其余符号的检测基于修正的接收信号矢量,从而消除了已经检测出来的符号对其他符号检测的影响,提高其余符号的检测性能。
在上述V-BLAST方法的递归实现中,公式(6c)是按照SNR最大准则计算SOPT≡[k1,k2,…,kM]T的公式。采用该公式表明:V-BLAST算法先计算SNR最大的符号的估计后,从接收信号矢量中消除该符号对其他符号检测的影响,再从修正的接收信号矢量中进行SNR次最大的符号的检测,然后再消除该符号对其余符号检测的影响。如此类推,SNR最小的符号最后进行检测。这种检测顺序使SNR最小的符号检测受益最大,从而使整个串行干扰对消方法的性能最优。
V-BLAST方法的复杂度较D-BLAST方法降低了,但是V-BLAST方法也存在不足之处:
(1)V-BLAST方法采用量化函数得到符号的估计,这是对符号的硬判决,参见公式(3)和公式(6f)。当符号硬判决结果不正确时,不仅不能消除干扰,反而会使干扰功率增加。因此,V-BLAST方法的符号估计方式需要改进。
(2)V-BLAST方法采用串行干扰对消,处理时延较大。通常MIMO通信中都需要接收端反馈ACK/NACK信息,由ACK/NACK信息确定是否进行符号重传。因此,接收机端的处理时延越小,ACK/NACK反馈的时延就越小,相应地系统的容量和吞吐量就越大。所以,V-BLAST的串行对消方式需要改进。
发明内容
有鉴于此,本发明的目的在于提供一种分组垂直分层空时检测方法,使其不仅能提高符号估计的性能,而且能减小MIMO系统的处理时延,进而提高系统的容量和速率。
本发明的另一个目的在于提供一种MIMO系统的分组垂直分层空时检测装置,使其能提高符号估计的性能,减小MIMO系统的处理时延,从而提高系统容量和速率。
本发明提供的一种分组垂直分层空时检测方法包括:
将接收信号矢量中待检测的符号分为Z组,Z为大于1的整数,令i=1,该方法还包括:
A.检测接收信号矢量,计算出第i组内所有符号的估计值,消除该组中所有符号的干扰,得到部分修正的接收信号矢量,再将部分经过修正的接收信号矢量作为需要检测的接收信号矢量,令i=i+1,再次执行步骤A,并且直到执行Z-1次步骤A后,执行步骤B;
B.检测接收信号矢量,计算出第Z组内所有符号的估计值,将每组所有符号的估计值的组合作为接收信号矢量的检测结果。
所述将接收信号矢量中待检测的符号分为Z组是利用信道矩阵H按照一定的原则进行分组、重排实现。
所述分组、重排的步骤包括:
A01、将接收信号矢量对应的M个符号a1,a2,…,aM分成Z组,第i(i=1,2,…,Z)组符号数为Li,第i组Li个符号为bi1,bi2,…,biLi,其中,
j=1,2,…,Li;Mi,j(Mi,j∈[1,2,…,M])表示第i组第j个符号在原来符号集合中的序号。
A02、将符号a1,a2,…,aM重排为:611,b12,…,b1L1,b21,…,b2L2,…,bZ1,…,bZLZ;则对应的发送符号矢量转换为:
A03、信道矩阵H=[H1,H2,…,HM],利用与序列a1,a2,…,aM同样的方法对序列H1,H2,…,HM进行重排,得到序列Q1,Q2,…,QM,记Q=[Ql,Q2,L,QM],则r1=Qb+v。
所述第i组的符号的估计值的计算步骤包括:
A11、利用接收信号矢量和矩阵Q的估计计算出迫零权矢量,wi,j,j=1,2,…,Li;
A12、分别利用每个迫零权矢量计算出Li个决策统计量, j=1,2,…,Li;ri是待检测的接收信号矢量。
A13、利用每个决策统计量计算出每个符号的估计值,i,j=D(yi,j),j=1,2,…,Li。
所述迫零权矢量是通过以每组的信道矩阵的估计为参数计算得到的。
所述迫零权矢量采用ZF准则,或MMSE准则计算得到。
在ZF准则下,并且在第一次执行步骤A时,由信道矩阵
按照公式 和公式
计算得到迫零权矢量。
在不重新进行信道估计的情况下,信道矩阵
的前
列为0列矢量,而第(J+1)列至第M列为
的第(J+1)列至第M列;
在MMSE准则下,在第一次执行步骤A时,由信道矩阵
按照公式 和公式
计算得到迫零权矢量,其中,wi,j为迫零权矢量,σs 2是发送符号的功率,σ2为接收信号矢量中噪声分量的功率。
在MMSE准则下,在第i(1<i≤z)组符号检测时,由信道矩阵
按照公式
和公式
计算得到迫零权矢量,其中,σs 2是发送符号的功率,σ2为接收信号矢量中噪声分量的功率。
在重新进行信道估计的情况下,信道矩阵
的前
列为0列矢量,而第(J+1)列至第M列通过信道估计得到,其中,第j∈[J+1,M]列通过对Q矩阵的第j∈[J+1,M]列重新进行信道估计而得到。
步骤A13中每个符号的估计值通过软判决函数确定。
当采用软判决函数确定符号的估计值时,所述符号的估计采用校正因子修正信道估计不理想状态下对符号估计的影响。
所述通过软判决函数确定符号的估计值是通过公式 计算出,
其中,βij称为校正因子,用以校正信道估计不理想造成的符号估计的偏差和干扰对消的偏差;f(yij|Ak)表示发送符号为Ak时接收到yij的概率密度函数。
若发送的符号是复数,则
其中,ARk、AIk分别为Ak的实部和虚部;YRij、YIij分别为yij的实部和虚部;VRij、VIij分别为vij的实部和虚部,分别服从N(0,σij 2)分布,且βij是复数。
若发送符号是实数,则
其中,YRij为yij的实部;VRij为vij的实部,服从N(0,σij 2)分布,且βij是实数。
βij的取值可以是符号的SNR的函数,也可以直接取1。
本发明提供的一种分组垂直分层空时检测装置包括:
符号分组与重排单元、一个以上个组内符号检测与干扰对消单元和组内符号检测单元;其中,
符号分组与重排单元,用于接收信号矢量、将接收信号矢量对应的符号分组、重排为若干组,并将接收信号矢量输出至第一个组内符号检测与干扰对消单元;
每个组内符号检测与干扰对消单元,检测接收信号矢量中对应的组内符号,得到符号的估计值,并将所有检测出的符号进行干扰对消,得到修正的接收信号矢量作为下一个组内符号检测与干扰对消单元/组内符号检测单元的待检测的接收信号矢量。
组内符号检测单元,用于接收最后一个组内符号检测与干扰对消单元发送的修正的接收信号矢量,估计出对应组内符号的估计值。
所述每个组内符号检测与干扰对消单元包括:
权矢量计算单元,用于根据信道矩阵的估计计算出一个以上个迫零权矢量,并将其输出至每个决策统计量计算单元;
一个以上个决策统计量计算单元,用于接收信号矢量,并根据接收到的每个迫零权矢量确定决策统计量,并将其输出至发送符号估计单元;
一个以上个发送符号估计单元,根据接收到的决策统计量估计发送符号,并将发送符号的估计值输出至干扰对消单元;
干扰对消单元,接收信号矢量以及发送符号的估计,并根据发送符号的估计对接收信号矢量进行干扰对消,得到一个修正的接收信号矢量并将其输出至下一个组内符号检测与干扰对消单元/组内符号检测单元。
所述组内符号检测单元包括:
权矢量计算单元,用于根据信道矩阵的估计计算出一个以上个迫零权矢量,并将其输出至每个决策统计量计算单元;
一个以上个决策统计量计算单元,用于接收修正的接收信号矢量,并根据接收到的每个迫零权矢量确定决策统计量,并将其输出至发送符号估计单元;
一个以上个发送符号估计单元,根据接收到的决策统计量估计发送符号,得到一个以上发送符号的估计值。
本发明将所有接收信号矢量中待检测的符号分成Z组。在组与组之间进行串行检测与串行干扰对消,在组内进行并行检测与并行干扰对消。本发明采用软判决方法提高符号估计的性能。本发明可以有效地降低MIMO系统的处理时延,并减少计算量,提高符号估计的性能,使系统的容量和速率得到提高。
附图说明
图1为实现本发明装置的结构示意图;
图2为每个组内符号检测与干扰对消单元的结构示意图;
图3为组内符号检测单元的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合附图对本发明作进一步的详细描述。
本发明的核心思想:接收信号矢量,将所有接收信号矢量中待检测的符号分成Z组。在组与组之间进行串行检测与串行干扰对消,在组内进行并行检测与并行干扰对消。并采用软判决方法提高符号估计的性能。具体的说,在第一组内,并行进行所有符号的检测,得到所有符号的估计,并从接收信号矢量中消除第一组内所有符号的干扰,得到修正的接收信号矢量,该接收信号矢量作为第一组给第二组的输入信号用于第二组符号的检测,使第二组符号的检测不受第一组符号的干扰;将第一组输出的接收信号矢量作为第二组内每个符号检测时采用的接收信号矢量,并行进行所有符号的检测,得到第二组内每个符号的估计,并从第一组输出的接收信号矢量中消除第二组内所有符号的干扰,得到进一步修正的接收信号矢量,该接收信号矢量作为第二组给第三组的输入信号矢量用于第三组符号的检测。依此类推,进行其他后续各组符号的检测。在进行第Z组符号的检测时,把第(Z-1)组输入的接收信号矢量作为第Z组内每个符号检测时采用的接收信号矢量,并行进行所有符号的检测,得到第Z组内每个符号的估计。至此,已经得到Z组中每组符号的估计,完成了所有符号的检测。
实现本发明的装置如图1所示,图中组的总数为Z。从图1可以看到:本发明的装置包括符号分组与重排单元11、(Z-1)个组内符号检测与干扰对消单元12和组内符号检测单元13。其中,符号分组与重排单元11将M个符号按照一定准则分成Z组,并按照分组顺序与组内符号顺序将M个符号进行重新排列,相应地组成信道矩阵的列矢量也进行重排形成新的信道矩阵;(Z-1)个组内符号检测与干扰对消单元分别完成第一组至第(Z-1)组符号的检测与干扰对消,分别得到每组内所有符号的估计和修正的接收信号矢量;组内符号检测单元13完成第Z组符号的检测,得到组内所有符号的估计。
图2为第i个组内符号检测与干扰对消单元12的结构示意图。参见图2,每个组内符号检测与干扰对消单元12包括权矢量计算单元21、决策统计量计算单元22、发送符号估计单元23、干扰对消单元24。
图3为组内符号检测单元13的结构示意图。组内符号检测单元包括权矢量计算单元21、决策统计量计算单元22、发送符号估计单元23。
下面结合图1、图2和图3说明实现本发明方法的具体过程。
本发明实现分组垂直分层空时检测方法的步骤如下:
第一步:符号分组与重排。
步骤101:符号分组与重排单元11由接收信号矢量和信道矩阵H的估计
按照某个准则将M个符号a1,a2,…,aM分成Z组,第i(i=1,2,…,Z)组符号数为Li,第i组Li个符号为bi1,bi2,…,biLi。其中,
表示第i组第j个符号在原来符号集合中序号为Mi,j(Mi,j∈[1,2,…,M]),j=1,2,…,Li,j的取值随组号i的变化而变化。
步骤102:将符号a1,a2,…,aM进行重排,重排结果为:b11,b12,…,b1L1,b21,…,b2L2,…,bZ1,…,bZLZ。因此,对应的发送符号矢量转换为:
公式(1)中的信道矩阵H可以写成列矢量的形式:H=[H1,H2,…,HM]。
接下来,对序列H1,H2,…,HM进行与序列a1,a2,…,aM同样的重排,得到的序列为Q1,Q2,…,QM,记Q=[Q1,Q2,…,QM]。
这样,经过H和a的元素重排,公式(1)可以更新如下:
r1=Qb+v (1a)
第二步:第一组符号的检测
第一个组内符号检测与干扰对消单元12中的权矢量计算单元21根据Q矩阵的估计
计算出一个以上的迫零权矢量,并将每个迫零权矢量依次输入决策统计量计算单元;每个决策统计量计算单元22根据接收信号矢量和迫零权矢量计算出决策统计量,并将计算出的决策统计量发送至与其相连的发送符号估计单元23。发送符号估计单元23分别根据收到的决策统计量,估计出接收信号矢量中对应的第一组内符号的估计值,然后将得到的符号估计值输出至干扰对消单元24进行组内干扰对消,得到修正的接收信号矢量,经过修正的接收信号矢量作为第二组检测的输入信号,输出至第二个组内符号检测与干扰对消单元。
第三步:第二组至第(Z-1)组符号的检测
在第i(1<i≤Z-1)个组内符号检测与干扰对消单元12中的权矢量计算单元21根据Q矩阵的估计
计算出一个以上的迫零权矢量,并将每个迫零权矢量依次输入决策统计量计算单元;每个决策统计量计算单元22根据第(i-1)组输出的经过修正的接收信号矢量ri和迫零权矢量计算出决策统计量,并将计算出的决策统计量发送至与其相连的发送符号估计单元23。发送符号估计单元23根据收到的决策统计量,估计出本组内对应符号的估计值,然后将得到的符号估计值输出至干扰对消单元24进行组内干扰对消,得到修正的接收信号矢量,经过修正的接收信号矢量作为第(i+1)组检测的输入信号,输出至第(i+1)个组内符号检测与干扰对消单元。当(i+1)=Z时,就输出至组内符号检测单元。
这里,第二组至第(Z-1)组符号的检测与第一组符号的检测的区别是:第一组接收信号矢量是未经过修正的,第二组至第(Z-1)组中的接收信号矢量是上一组经过修正的接收信号矢量。
第四步:第Z组符号的检测
在第Z组组内符号检测单元13中的权矢量计算单元21根据Q矩阵的估计
计算出一个以上的迫零权矢量,并将每个迫零权矢量依次输入决策统计量计算单元;每个决策统计量计算单元22根据第(Z-1)组输出的经过修正的接收信号矢量rZ和迫零权矢量计算出决策统计量,并将计算出的决策统计量发送至与其相连的发送符号估计单元23。发送符号估计单元23根据收到的决策统计量,估计出本组内对应符号的估计值,并将该估计值作为第Z组符号估计结果输出。
在上述第一步~第四步中,组内符号检测与干扰对消单元完成功能包括:
第i(i∈[1,2,…,Z-1])个组内符号检测与干扰对消单元12在输入信号矢量为ri、组内符号总数为Li时,完成组内Li个符号bi1,bi2,…,biLi的检测并得到修正的接收信号矢量ri+1。
下面再参见图2详细说明第i个组内符号检测与干扰对消单元的功能。其具体的计算步骤如下:
权矢量计算单元21由修正的接收信号矢量ri和Q矩阵的估计
按照公式计算出各迫零权矢量wi,j,j=1,2,…,Li。
Li个决策统计量计算单元22并行地计算出Li个决策统计量。其中,第j个决策统计量计算单元22由第j个迫零权矢量计算第j个决策统计量:
Li个发送符号估计单元23并行地计算出Li个符号的估计。第j个发送符号估计单元23由第j个决策统计量得到第j个发送符号的估计:
i,j=D(yi,j),j=1,2,…,Li (7b)
在公式(7b)中,以软判决函数D(·)替代了公式(3)和公式(6f)中的量化函数Q(·)。
公式(7a)可以表示如下:
yi,j=bij+vij,j=1,2,…,Li (7a-1)
其中,bij∈{A1,A2,…,AK},{A1,A2,…,AK}是所有可能的发送符号构成的集合;vij是高斯白噪声。
则(7b)具体按照下式计算:
其中,βij称为校正因子,用以校正信道估计不理想造成的符号估计的偏差和干扰对消的偏差;f(yij|Ak)表示发送符号为Ak时接收到yij的概率密度函数。βij和f(yij|Ak)的计算方法如下:
(1)若发送的符号是复数,则
其中,ARk、AIk分别为Ak的实部和虚部;YRij、YIij分别为yij的实部和虚部;VRij、VIij分别为vij的实部和虚部,分别服从N(0,σij 2)分布,且βij是复数。
(2)若发送符号是实数,则
其中,YRij为yij的实部;VRij为vij的实部,服从N(0,σij 2)分布。且βij是实数。
(3)βij的取值与符号的SNR密切相关,是符号的SNR的函数。
为了简化计算,也可以直接令βij=1,即:不考虑信道估计的偏差对符号估计和干扰对消造成的影响。
通常,在接收端符号的信噪比(SNR,Signal to Noise Ratio)要达到一定的数值,才能使解映射和译码后的比特误块率(BLER,Block Error Rate)性能满足业务质量的要求。对于给定的BLER数值,可以通过仿真确定:为使译码以后BLER性能达到要求,符号所需要达到的最低SNR数值SMIN。并取T1=SMIN-δ1,T2=SMIN-δ2。其中,δ1>0,δ2>δ1。
如果第i组符号j的SNR大于阈值T1=SMIM-δ1,就认为该符号的信道估计比较准确,可以近似认为βij=1。
如果第i组符号j的SNR小于或等于阈值T2=SMIN-δ2,就认为该符号的SNR太低,使信道估计非常差,可以近似认为βij=0。
如果第i组符号j的SNR大于阈值T2=SMIN-δ2且小于等于阈值T1=SMIN-δ1,可以通过COSSAP(Communication System Simulation Analysis Package,一种仿真软件包)仿真或MATLAB(矩阵实验室,Matrix Laboratory,一种仿真软件包)仿真优化确定:当SNR在(T2,T1)内变化时βij的具体取值。
δ1、δ2取值直接影响βij对符号估计偏差的校正的精度和SNR优化区间(T2,T1)的大小。可以根据对校正精度和优化计算量的需要确定δ1和δ2的数值。
综上所述,βij=β取值如下:
上式中,αi的取值通过在较小的SNR区间[SNRi-I,SNRi]的COSSAP仿真或MATLAB仿真优化得到。其中,I是对SNR优化区间(T2,T1)量化得到的小区间的个数,根据需要确定。
干扰对消单元24由组内所有符号的估计,得到修正的接收信号矢量,修正的接收信号矢量用于第(i+1)组符号的检测。修正的接收信号矢量为:
其中,
是在第i个组内符号检测时信道矩阵Q的估计。
当i=1时,即在第一个组内符号检测时,
作为Q矩阵的估计,可以由信道估计得到
然后由
通过重排得到
其中,第j∈[J+1,M]列通过对Q矩阵的第j∈[J+1,M]列重新进行信道估计而得到。
组内符号检测单元完成第Z组符号的检测。该单元输入信号矢量为rZ,组内包含LZ个符号:bZ1,bZ2,…,bZLZ。
参见图3,权矢量计算单元21由修正的接收信号矢量rZ和Q矩阵的估计
按照公式计算出各迫零权矢量wZ,j,j=1,2,…,LZ。
LZ个决策统计量计算单元22并行地计算出LZ个决策统计量。
令i=Z,第j个决策统计量计算单元22由第j个迫零权矢量按照公式(7a)计算第j个决策统计量。
LZ个发送符号估计单元23并行地计算出LZ个符号的估计。令i=Z,第j个发送符号估计单元23由第j个决策统计量按照公式(7b)、(7b-1)、(7b-2)、(7b-3)和(7b-4)得到第j个发送符号的估计。
在上述分组垂直分层空时检测方法的实现过程中,迫零权矢量的计算是关键技术,权矢量的计算以信道矩阵为参数,可以采用不同的优化准则计算得到,比如,ZF准则和MMSE准则等。而采用不同准则,迫零权矢量的计算公式不同。下面给出ZF准则和MMSE准则下分组垂直分层空时检测方法中迫零权矢量计算的公式。
在MMSE准则下,在第i(1<i≤Z)组符号检测时,可以由信道矩阵
按照公式(9c)和公式(9d)计算得到迫零权矢量:
在重新进行信道估计的情况下,信道矩阵
的前
列为0列矢量,而第(J+1)列至第M列通过信道估计得到。其中,第j∈[J+1,M]列通过对Q矩阵的第j∈[J+1,M]列重新进行信道估计而得到。
本发明可以预先设置接收信号矢量中被检测符号的组数。在多用户检测中,需要将用户分组,具体的分组方式很多。如:可以根据符号的SNR进行分组。可以根据SNR的均值和方差设置若干个门限,将所有符号分成Z组等。第一组符号的SNR比其他各组符号的SNR都高。第二组符号的SNR比后面各组符号的SNR都高。依此类推,最后一组符号的SNR比前面各组符号的SNR都小。
本发明采用了H和a的重排。在经过分组处理得到发送符号矢量b的估计,可以对矢量b的估计进行反重排,可以得到矢量a=[a1,a2,…,aM]T的估计,该矢量的估计可以用于后续的译码等处理,并可以减少计算量。当然,由于H和a是否进行重排并不影响各组符号的检测结果,因此也可以不将H和a进行重排。
总之,以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。
Claims (21)
1、一种分组垂直分层空时检测方法,其特征在于,将接收信号矢量中待检测的符号分为Z组,Z为大于1的整数,令i=1,该方法还包括以下步骤:
A.检测接收信号矢量,计算出第i组内所有符号的估计值,从接收信号矢量中消除该组中所有符号的干扰,得到部分修正的接收信号矢量,再将部分经过修正的接收信号矢量作为需要检测的接收信号矢量,并令i=i+1,再次执行步骤A,并且直到执行Z-1次步骤A后,执行步骤B;
B.检测接收信号矢量,计算出第Z组内所有符号的估计值,将每组所有符号的估计值的组合作为接收信号矢量的检测结果。
2、根据权利要求1所述的方法,其特征在于,所述将接收信号矢量中待检测的符号分为Z组是利用信道矩阵H按照一定的原则进行分组、重排实现。
3、根据权利要求2所述的方法,其特征在于,所述分组、重排的步骤包括:
A01、将接收信号矢量对应的M个符号a1,a2,…,aM分成Z组,第i(i=1,2,…,Z)组符号数为Li,第i组Li个符号为bi1,bi2,…,biLi,其中,
Mi,j(Mi,j∈[1,2,…,M])表示第i组第j个符号在原来符号集合中的序号;
A02、将符号a1,a2,…,aM重排为:b11,b12,…,b1L1,b21,…,b2L2,…,bZ1,…,bZLZ;则对应的发送符号矢量转换为:
A03、信道矩阵H=[H1,H2,…,HM],利用与序列a1,a2,…,aM同样的方法对序列H1,H2,…,HM进行重排,得到序列Q1,Q2,…,QM,记Q=[Q1,Q2,L,QM],则r1=Qb+v。
4、根据权利要求2所述的方法,其特征在于,所述第i组的符号的估计值的计算步骤包括:
A11、利用接收信号矢量和矩阵Q的估计计算出迫零权矢量,wi,j,j=1,2,…,Li;
A12、分别利用每个迫零权矢量计算出Li个决策统计量, ri是待检测的接收信号矢量;
A13、利用每个决策统计量计算出每个符号的估计值,
5、根据权利要求4所述的方法,其特征在于,所述迫零权矢量是通过以每组的信道矩阵的估计为参数计算得到的。
6、根据权利要求4所述的方法,其特征在于,所述迫零权矢量采用ZF准则,或MMSE准则计算得到。
7、根据权利要求4所述的方法,其特征在于,在ZF准则下,并且在第一次执行步骤A时,由信道矩阵
按照公式
和公式 计算得到迫零权矢量。
9、根据权利要求8所述的方法,其特征在于,
在不重新进行信道估计的情况下,信道矩阵
的前
列为0列矢量,而第(J+1)列至第M列为
的第(J+1)列至第M列;
在重新进行信道估计的情况下,信道矩阵
的前
列为0列矢量,而第(J+1)列至第M列通过信道估计得到,其中,第j∈[J+1,M]列通过对Q矩阵的第j∈[J+1,M]列重新进行信道估计而得到。
10、根据权利要求4所述的方法,其特征在于,在MMSE准则下,在第一次执行步骤A时,由信道矩阵
按照公式 和公式
计算得到迫零权矢量,其中,wi,j为迫零权矢量,σs 2是发送符号的功率,σ2为接收信号矢量中噪声分量的功率。
11、根据权利要求4所述的方法,其特征在于,在MMSE准则下,在第i(1<i≤Z-1)次执行步骤A和执行步骤B时,由信道矩阵
按照公式 和公式
计算得到迫零权矢量,其中,σs 2是发送符号的功率,σ2为接收信号矢量中噪声分量的功率。
13、根据权利要求4所述的方法,其特征在于,步骤A13中每个符号的估计值通过软判决函数确定。
14、根据权利要求13所述的方法,其特征在于,当采用软判决函数确定符号的估计值时,所述符号的估计采用校正因子修正信道估计不理想状态下对符号估计的影响。
15、根据权利要求13所述的方法,其特征在于,所述通过软判决函数确定符号的估计值是通过公式
其中,βij称为校正因子,用以校正信道估计不理想造成的符号估计的偏差和干扰对消的偏差;f(yij|Ak)表示发送符号为Ak时接收到yij的概率密度函数。
16、根据权利要求15所述的方法,其特征在于,若发送的符号是复数,则
17、根据权利要求15所述的方法,其特征在于,若发送符号是实数,则
其中,YRij为yij的实部;VRij为vij的实部,服从N(0,σij 2)分布,且βij是实数。
18、根据权利要求15所述的方法,其特征在于,βij的取值是符号的SNR的函数,或βij=1。
19、一种分组垂直分层空时检测装置,其特征在于,该装置包括:符号分组与重排单元、一个以上个组内符号检测与干扰对消单元和组内符号检测单元;其中,
符号分组与重排单元,用于接收信号矢量、将接收信号矢量对应的符号分组、重排为若干组,并将接收信号矢量输出至第一个组内符号检测与干扰对消单元;
每个组内符号检测与干扰对消单元,检测接收信号矢量中对应的组内符号,得到符号的估计值,并将所有检测出的符号进行干扰对消,得到修正的接收信号矢量作为下一个组内符号检测与干扰对消单元/组内符号检测单元的待检测的接收信号矢量;
组内符号检测单元,用于接收最后一个组内符号检测与干扰对消单元发送的修正的接收信号矢量,估计出对应组内符号的估计值。
20、根据权利要求18所述的装置,其特征在于,所述每个组内符号检测与干扰对消单元包括:
权矢量计算单元,用于根据信道矩阵的估计计算出一个以上个迫零权矢量,并将其输出至每个决策统计量计算单元;
一个以上个决策统计量计算单元,用于接收信号矢量,并根据接收到的每个迫零权矢量确定决策统计量,并将其输出至发送符号估计单元;
一个以上个发送符号估计单元,根据接收到的决策统计量估计发送符号,并将发送符号的估计值输出至干扰对消单元;
干扰对消单元,接收信号矢量以及发送符号的估计,并根据发送符号的估计对接收信号矢量进行干扰对消,得到一个修正的接收信号矢量并将其输出至下一个组内符号检测与干扰对消单元/组内符号检测单元。
21、根据权利要求20所述的装置,其特征在于,所述组内符号检测单元包括:
权矢量计算单元,用于根据信道矩阵的估计计算出一个以上个迫零权矢量,并将其输出至每个决策统计量计算单元;
一个以上个决策统计量计算单元,用于接收修正的接收信号矢量,并根据接收到的每个迫零权矢量确定决策统计量,并将其输出至发送符号估计单元;
一个以上个发送符号估计单元,根据接收到的决策统计量估计发送符号,得到一个以上发送符号的估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100429333A CN100361414C (zh) | 2004-06-01 | 2004-06-01 | 一种分组垂直分层空时检测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100429333A CN100361414C (zh) | 2004-06-01 | 2004-06-01 | 一种分组垂直分层空时检测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1705242A true CN1705242A (zh) | 2005-12-07 |
CN100361414C CN100361414C (zh) | 2008-01-09 |
Family
ID=35577682
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004100429333A Expired - Fee Related CN100361414C (zh) | 2004-06-01 | 2004-06-01 | 一种分组垂直分层空时检测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100361414C (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102231641A (zh) * | 2011-07-21 | 2011-11-02 | 西安电子科技大学 | 多输入多输出mimo逐级并行检测方法 |
CN101572561B (zh) * | 2008-04-30 | 2012-09-05 | 华为技术有限公司 | 多输入多输出系统的信号检测方法和装置 |
CN101383652B (zh) * | 2007-09-04 | 2012-09-26 | 中兴通讯股份有限公司 | 一种多输入多输出系统的信号检测方法及装置 |
CN104869086A (zh) * | 2015-05-27 | 2015-08-26 | 东南大学 | 基于二维压缩感知的mimo-ofdm通信系统下行信道估计方法、装置 |
CN104977579A (zh) * | 2015-07-03 | 2015-10-14 | 中国科学院声学研究所 | 一种基于随机协方差矩阵的多亮点目标空时检测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2001270734A1 (en) * | 2000-07-04 | 2002-01-14 | Thales | Space-time coding digital transmission systems and methods |
CN1155190C (zh) * | 2001-10-09 | 2004-06-23 | 华为技术有限公司 | 垂直的贝尔实验室分层空时编码阵列线性检测方法 |
-
2004
- 2004-06-01 CN CNB2004100429333A patent/CN100361414C/zh not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101383652B (zh) * | 2007-09-04 | 2012-09-26 | 中兴通讯股份有限公司 | 一种多输入多输出系统的信号检测方法及装置 |
CN101572561B (zh) * | 2008-04-30 | 2012-09-05 | 华为技术有限公司 | 多输入多输出系统的信号检测方法和装置 |
US8437418B2 (en) | 2008-04-30 | 2013-05-07 | Huawei Technologies Co., Ltd. | Method and apparatus for detecting signals of multi-input multi-output system |
CN102231641A (zh) * | 2011-07-21 | 2011-11-02 | 西安电子科技大学 | 多输入多输出mimo逐级并行检测方法 |
CN102231641B (zh) * | 2011-07-21 | 2013-08-14 | 西安电子科技大学 | 多输入多输出mimo逐级并行检测方法 |
CN104869086A (zh) * | 2015-05-27 | 2015-08-26 | 东南大学 | 基于二维压缩感知的mimo-ofdm通信系统下行信道估计方法、装置 |
CN104869086B (zh) * | 2015-05-27 | 2017-11-14 | 东南大学 | 基于二维压缩感知的mimo‑ofdm通信系统下行信道估计方法、装置 |
CN104977579A (zh) * | 2015-07-03 | 2015-10-14 | 中国科学院声学研究所 | 一种基于随机协方差矩阵的多亮点目标空时检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100361414C (zh) | 2008-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1287531C (zh) | 自适应均衡装置 | |
CN1943194A (zh) | 校准无线mimo通信系统中的下行链路和上行链路信道响应 | |
CN1256814C (zh) | 无线电接收设备和无线电接收方法 | |
CN101048949A (zh) | 在通用rake接收机中用于参数估计的方法和设备 | |
CN101069375A (zh) | 发送控制帧生成装置、发送控制帧处理装置、发送控制帧生成方法及发送控制帧处理方法 | |
CN1679269A (zh) | 按特征模式应用选择性信道反转的已编码mimo系统 | |
CN1692572A (zh) | 在无线通信中用于多天线的空时发射分集(sttd) | |
CN1808959A (zh) | 一种传输数据的方法及通信系统 | |
CN1203693C (zh) | 自适应均衡方法及自适应均衡器 | |
CN1906864A (zh) | 空分多址的多天线传输 | |
CN1697362A (zh) | 用于在多输入多输出通信系统中接收信号的装置和方法 | |
CN1538700A (zh) | 通过顺序蒙特卡罗的接近最佳的多输入多输出信道检测 | |
CN1889555A (zh) | 空时比特交织调制系统的迭代译码算法及接收系统 | |
CN1155189C (zh) | 双层加权并行干扰对消方法的简化方法 | |
CN2579091Y (zh) | 利用阵列处理进行数据检测的用户装置 | |
CN1649335A (zh) | 相关器和利用其的接收装置 | |
CN1446409A (zh) | 一族线性多用户检测器(muds) | |
CN1943155A (zh) | 用于提供发射序列的估计的设备和方法 | |
CN1905383A (zh) | 一种同频小区信道估计装置及方法 | |
CN1845483A (zh) | Td-scdma系统中接收机的联合检测系统及其数据处理方法 | |
CN1701553A (zh) | 无线通信系统内质量反馈解码 | |
CN1288612A (zh) | 针对直接序列-码分多址信号的周期自适应接收器 | |
CN1140944A (zh) | 数据接收装置 | |
CN1282156A (zh) | Cdma移动通信系统中的cdma接收设备和接收信号功率测量设备 | |
CN1705242A (zh) | 一种分组垂直分层空时检测方法及装置 |
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: 20080109 Termination date: 20200601 |
|
CF01 | Termination of patent right due to non-payment of annual fee |