专利名称：用于码分多址系统的检测器的制作方法
技术领域：
本发明一般涉及经受多径衰落和多址联接干扰(MAI)的直接序列(DS)码分多址(CDMA)系统，尤其涉及在这种系统中使用的信号检测技术。
在如建议用于下一代移动无线标准的无线CDMA系统中，在物理层的主要障碍是由不正交于所需用户的同信道用户引起的多径衰落和多址联接干扰。在R.Price和P.E.Green的1958年3月的IRE会刊，第46卷第555到第570页中“用于多径信道的通讯技术”描述的瑞克接收机，通过所需信号的相干组合可分辨的多径复制试图消除多径衰落。在S.Verdu的剑桥大学出版社，New York，1998年的“多用户检测”中描述的多用户检测，通过当检测所需的用户时主动考虑它的存而提出MAI问题。
最近，人们更加感兴趣于使用阵列处理来进一步改进接收机性能。这些技术的焦点在于在基站接收机使用多个天线以提供天线增益和/或分集增益并且允许空间处理的可能性。通过将这些空间域技术与时间域技术如瑞克检测和多用户检测结合起来，得到的时空检测器与传统的只是时间域的检测器相比显示出能够改进CDMA系统能力。例如，参见A.Paulraj和C.Papadias的IEEE信号处理杂志，第14卷第六部分，第49-83页，1997年10月的“用于无线通讯的时空处理”。第一代时空CDMA检测器使用阵列处理和瑞克检测或者多用户检测。可见A.Naguib和A.Paulraj的在IEEE杂志通讯的选择区域，第14卷第9部分，第1770-1783页，1996年12月的“具有M进制正交调制和小区基站天线阵列的无线CDMA的性能”，以及参见S.Miller和S.Schwartz的IEEE通讯学报，第43卷第2/3/4部分，第396-411页，1995年2月、3月、4月的“用于异步高斯多址联接信道的集成的空间时间检测器”。新近的时空CDMA检测器结合了所有三种处理技术。参见H.Huang、S.Schwartz，S.Verdu的在IEEE国际专题讨论会刊的信息理论，第380页，1995年的“用于多用户检测的组合的多径和空间分辨率潜力和问题”；或者M.Nagatsuka和R.Kohno的IEICE通讯学报，第E78-B卷第11部分，第1489-1497页，1995年11月的“对于DS/CDMA使用阵列天线的空间和时间优化的多用户接收机”。
尽管前述的系统可以满意地工作，也可以特别是对于结合所有三种处理技术的时空检测器进行改进。尽管Nagatsuka和Kohno论文中的检测器在最大似然意义上是优化的，但它的计算复杂性对于用户数量成指数关系。因此对于实际系统来说实现起来太复杂。在Huang、Schwartz和Verdu的论文中，实现了性能和复杂性之间的折衷，但是这种检测器不能自适应地实现，因为它使用迫零准则(zero-forcing criteria)。自适应实现允许接收机考虑未知的干扰源因此改进检测器性能并且提高系统性能。例如，基站接收机可以考虑来自相邻小区或来自嵌入微小区的干扰，而手机接收机可以考虑来自不能明确地解调的信号的干扰。
本发明利用一个瑞克(rake)接收机、阵列处理和多用户检测来检测DS-CDMA信号。当结合阵列处理时，瑞克接收机常常被称为时空瑞克接收机。本发明在多用户检测器中使用最小均方差(MMSE)准则。这个准则允许以线性检测器的形式相对简单地实现并且也允许自适应地实现。自适应实现在实际情况下是有用的，在这些情况下对各种接收信号的了解是有限的。如上所述，使用考虑未知干扰的自适应检测器可以改进上行链路和下行链路的能力。
对于脉冲幅度调制(PAM)数据信号，本发明的两个实施例提供了在性能和自适应实现复杂性之间折衷的选择。如果完全了解用户的信号参数，则实现第一个实施例较好。另一方面，对于自适应实现要求更明确的信道信息。假设接收机包括P≥1个天线并且接收信号包括K个DS-CDMA数据信号，每个具有L个延迟的/加权的多径复制。在每个天线，滤波器组匹配于具有它们的多径定时延时的KL个扩展码。滤波器输出根据估计信道(多径和阵列)参数的复数共轭被加权并且被组合以形成一个K向量，这里每个分量相应于K个代码中的一个。采用每个分量的实数部分。例如，可以从序列信号或导频信号获得多径延迟和信道参数的估计。常规的时空瑞克接收机将每个向量分量传递到一个判决装置用于估计相应的用户PAM数据码元。然而，除非用户的信号在时空域正交，否则向量分量被来自其他用户的多址联接干扰损害。本发明建议在判决装置之前使用一个线性组合器以抑制多址联接干扰。这个线性组合器通过一个乘以实数K向量的实数K×K矩阵WA表示。该矩阵将它具有实数K向量的乘积和数据码元K向量之间的均方差减至最小。因为矩阵使用最小均方差(MMSE)准则，所以可以使用公知的自适应算法自适应地获得它。最后的K向量输出的每个分量相应于K个代码中的一个并且被传送到一个判决限制器或者一个译码器用于进一步处理。
如果准确地知道信道参数，第二个实施例在一定程度上具有较差的性能。然而，在信道不匹配的实际情况下，实现第二个实施例常常较好。在自适应实现方面，这个实施例要求较少的信息(它不要求精确的信道估计)但可能适应得较慢。正如第一个实施例的情况，第二个实施例的前端同样包括KL个匹配的滤波器组用于后面有加权和组合的P个天线中的一个。然而，不像第一个实施例，实数运算符的次序和线性组合器(K×K矩阵相乘)交换。这时，加权、组合和线性组合器可以由单个复数K×KLP矩阵WB表示。在产生的K向量中分量的实数部分传送到一个判决装置或一个译码器用于进一步处理。在自适应实现中，自适应算法可以用于获得矩阵WB，因此不明确地要求信道估计。
对于正交PAM(QAM)数据信号，在检测器中不要求实数运算符，因为信号构象是两维的。在这种情况下，可能使用第二个实施例的变型，它不使用实数运算符但直接跟随着具有适当判决装置或译码器的线性组合器。
尽管MMSE检测技术本身是强有力的，但通过使用它们以及其他多用户检测技术如干扰对消可以进一步增强它们的性能，这里干扰不是如线性多用户检测器所做的投射出去，而是明确地从接收信号中提取出来。干扰对消可以在MMSE线性组合器之前或之后发生。
MMSE检测器的两个实施例可以推广到实际的系统，这些系统除了利用数据承载信道以外，还利用可以用作导频信号的辅助信道。同样，MMSE检测器可以在信号以不同扩展因子发射的系统中工作。
相对于不使用多用户检测或使用多用户检测但不自适应地减少未知干扰的常规信号检测器，本发明的CDMA信号检测器可以提供明显改进的性能。
下面结合附图来描述，本发明的前面和其他特性可以变得更加明显，其中
图1示出一个现有技术单用户时空瑞克接收机，它使用一个多元天线阵列、在每个天线的一个N芯片滤波器组以及一个信道加权器和组合器；图2示出本发明的通用实施例；图3示出一个多用户时空MMSE检测器的第一个实施例(检测器A)；图4示出一个多用户时空MMSE检测器的第二个实施例(检测器B)；图5示出第二个实施例的一个自适应实现；图6示出在第一个实施例的线性组合器之前的干扰对消；图7示出在第二个实施例的线性组合器之前的干扰对消；图8示出在本发明的多用户时空MMSE检测器的线性组合器之后的干扰对消；及图9示出根据本发明的后置组合器干扰对消接收机的方块图。
考虑一个K用户系统，这里第k个用户(k＝1…K)用一个N芯片扩展序列sk(t)调制它的数据序列bk(t)。对于检测器A，数据序列是脉冲幅度调制(PAM)的；对于检测器B，数据序列可以是PAM或正交PAM(QAM)的。发射的信号受到信道中的频率选择性衰落并且作为具有复数衰落信道系数ck，l(t)…ck，L(t)的L时间可分辨多径分量到达接收机。假设该接收机是一个P元线性阵列。如果每个可分辨多径分量作为具有一个相应于线性阵列的角度θk，l(t)的平面型波前到达并且如果阵列空间足够的靠近(例如λ/2)使得对于一个给定波前在阵列元之间存在完全的相关，则相对于第一个元的第p个元的相位偏移是hk，l，p＝exp(π(p-1)sinθk，l(t))。对于给定码元周期(并且忽略码元间干扰)在第p个天线接收的信号是rp(t)=Σk=1KΣl=1Lhk,l,p(t)Ck,l(t)Sk(t)(t-τk,l)Akb(t-τk,l)+np(t)---(1)]]>这里Ak是用户k的幅度，τk，l是对于第k个用户的第l个多径的延迟，而np(t)是考虑小区外干扰和背景噪声的加性高斯噪声处理。现在我们做出下列假设来简化分析。
(a)位同步地接收信号。
(b)与码元周期相比较时间扩展是小的使得可以忽略码元间干扰；以及(c)在一个码元周期相位偏移和信道系数是恒定的。
第一个假设以后将省略。在这些假设下，对于方程式(1)的接收信号，芯片匹配滤波器输出是一个复数N向量rp=Σk=1KΣl=1Lhk,l,pCk,lSk,lAkb+np]]>这里sk，l是相应于sk(t-τk，l)的芯片匹配滤波器N向量，而np是相应于高斯噪声的复数N向量。扩展码被归一化为具有单位能量｜｜sk，l｜｜＝1。假设扩展码是随机的。然而，这些MMSE检测器的自适应实现要求使用在几个码元周期之后重复的短扩展码。下列符号将在以后使用hk，l＝[hk，l，l…hk，l，p]T阵列系数的复数P向量H＝[hl，l…hl，L…hK，l…hK，L]复数P×KL阵列矩阵HD＝diag[hl，l…hl，L…hK，l…hK，L]复数KLP×KL阵列矩阵ck＝[ck，l…ck，L]T信道系数的复数L向量
C＝diag(c1…ck)复数KL×K信道矩阵R复数KL×KL相关矩阵由R((k1-1)K+l1,(k2-1)K+l2)=Sk1,l1HSk2,l2]]>定义A＝diag(A1…Ak) 实数K×K幅度矩阵b＝[bl…bk]T数据位的实数K向量Iuu×u单位矩阵1u1的u向量0u×u0的u×u矩阵0u0的u向量对于阵列和信道系数的相应估计值将用在该值的码元上的符号“^”表示。噪声向量是一个零均值复数高斯向量，它的分布可以根据它的(按照分量)实数和虚数分量写出Re(n)Im(n)~η(0N0N,σ2IN0NxN0NxNIN|)]]>这里我们定义对于矩阵和向量的实数和虚数操作为Re(X)＝(X+X*)/2和Im(X)＝(X-X*)/2，这里*表示复数共轭。因此Re(n)和Im(n)是零均值高斯随机向量，它的分量具有方差σ2并且相互间不相关。
常规的(时空瑞克)接收机，(图1)在这个阵列多径信道的上下文中，现有技术单用户检测器是一个匹配于所需用户的复合阵列多径扩展码信号的相关器。这个检测器不考虑干扰的存在；然而，如果不存在任何干扰或者如果它们正交于阵码空间中所需的用户，则它是最大似然检测器。如图1所示，检测器包括在每个天线的一组相关器10，匹配于KL个多径扩展码sl，l…sl，L…sK，l…sK，L。在图1最左边方框中的符号<sk，l，·>表示采用sk，l(相应于第l个多径延迟的第k个用户的扩展码)和到该方框的输入向量之间的内积。在符号<sk，l，·>中的点表示输入。使用单独的定时估计算法获得定时估计。上面左边方框的输出(它相应于向量Z1的第一个分量)是<sl，l，rl>＝sl，lHrl，这里H上标表示一个复数向量的厄米特转置(进行每个元的复数共轭，然后进行生成向量的转置)。用于至少一个多信号子集、这些信号的多个多径分量以及多个接收机天线的匹配滤波器输出在图中用(z1…zp)表示。
在参考标号11，相关器输出由它们的相应信道估计的复数共轭加权。更准确地说，在用于用户k的第l个多径的第p个天线的相关器输出由相应信道(阵列/多径)系数估计的复数共轭
加权。使用单独的信道估计算法获得这些估计。在每个天线，加上用于第k个用户的L个分量(12)，然后用于第k个用户的生成的P个分量加到13k上(k＝1…K)。每个分量随后传送到一个判决装置或一个译码器。对于PAM数据信号，用于第k个用户(14k，k＝1…K)的判决装置给出作为它输出的码元，该码元对于第k个向量分量实数部分在欧几里德距离意义上最近。对于QAM数据信号，不要求实数运算符。
图1中，在匹配滤波器输出端的KLP-向量z＝[z1T…zpT]T可以写为z=R~HDCAb+ns---(2)]]>这里我们定义KLP×KLP矩阵
和KLP×KL矩阵，表示两个矩阵之间的直积(克罗内克积)运算，而复数噪声向量是具有下面分布的复数高斯随机向量Re(ns)Im(ns)~η(0N0N,σ2Re(R~)Im(R~)-Im(R~)Re(R~)|)]]>使用
(这里·表示两个相同尺寸矩阵之间按分量逐个作出的乘积)，在信道组合器输出端的K向量和到判决装置141-14K的输入可以简明地写成Re{yconv}=Re{(H^HC^)Hz}=MAb+nconv---(3)]]>
这里H表示复数共轭转置，·表示两个相同尺寸矩阵之间按分量逐个作出的乘积，并且
。噪声向量完全是实数并且具有分布nconv~η(0N,σ2M^)]]>这里
。如果该数据是BPSK调制的，则对于用户k的位判决是简单的Re{yconv}的第k个分量的硬限制k＝sgn(Re{yconv，k})。常规瑞克接收机的相应的误码率Pkconv是Pkconv(σ)=Q[(MAb)kσM(k,k)^]---(4)]]>这里X(k，k)是矩阵X的第(k，k)个元素。根据BPSK调制数据给出在各种检测器之间的性能比较，但是一般也可以对于任何QAM调制数据进行这些性能比较。
通用时空线性多用户瑞克接收机(图2)上述现有技术时空瑞克接收机在给定用户意义上是单用户接收机，解调仅仅使用来自该用户的信息。因为它不考虑来自其他用户的干扰的存在，它的性能受到损失。图2示出一个通用的多用户时空瑞克接收机。它使用一个信息组合器考虑多址联接干扰。给出这种检测器的两种版本。第一种称为检测器A，在线性组合器之前使用实数运算符。在这种情况下，图2中参考号20包括跟随着一个线性组合器的一组实数运算符，分量141…14k中的每一个是一个限制器，它确定对于它输出的最接近估计码元。第二种称为检测器B，使用跟随着实数运算符由复数矩阵表示的线性组合器。在这种情况下，参考号20仅仅包括线性组合器，而判决装置包括实数运算符和限制器。尽管类似，但用于检测器A和B的线性组合器导致了不同的性能和自适应实现。现在我们更充分地描述两个检测器实施例。
时空线性多用户检测器A(图3)在方程式(3)给出用于K个用户的一组充分的统计量作为K向量Re{yconv}。线性MMSE检测器的目的是将K×K线性组合器应用于这个向量使得在它生成向量和数据向量b之间的均方差减至最小。换句话说，找到实数K×K矩阵WA使得WA=argminW{E[||WRe{yconv}-b||2]}----(5)]]>这个检测器实质上在消除用nconv表示的残留高斯噪声和用
的非对角线项表示的多址联接干扰中取得适当的平衡。可以示出该解决方案为WA＝E[b(Re{yconv})H]{E[Re{yconv}Re{yconv}]}-1(6)=AM(MA2MH+σ2M^)-1]]>图3示出一个根据本发明的检测器A的方块图。它包含P个天线和一组匹配的滤波器10。方框31表示一个由方程式(6)给出的K×K矩阵WA定义的线性组合器。该矩阵可以如下解释它是3项的乘积(a)对角幅度矩阵(A)；(b)等效的阵列/信道/代码相关矩阵(M)；以及矩阵(
)的逆矩阵，它包括阵列/信道/代码相关矩阵(M)、幅度矩阵的平方(A2)和阵列/信道/代码相关矩阵(MH)的厄米特共轭的乘积与由背景噪声方差(σ2)加权的估计阵列/信道/代码相关矩阵
的和。这个矩阵使在方程式(5)给出的均方差准则减至最小，并且因为在(5)中均方差表示的凸度，MMSE解决方案(6)表示一个总的最小值。一组判决限制器(321-32k)跟随着线性组合器。
MMSE检测器的一个具有吸引力的特性是它能够使用公知的自适应算法(33)如最小均方或循环最小二乘来自适应地实现。使用一个序列信号(例如，用户数据位)，这些技术可以用于自适应地获得WA。如果认为WA的直接计算太复杂，则自适应实现是一个选择。注意自适应实现要求了解K个用户的扩展码、多径延迟和信道参数。
假设阵列和信道估计是准确的
，我们可以重新写出(6)中WA表达式WA＝A-1[Re{CH[(HHH)·R]C}+σ2A-2]-1另一个吸引人的特性是随着残留高斯噪声趋近于零，MMSE检测器趋近于解相关(迫零)检测器(高到一个标量因子)limσ→0(WA)=A-1(Re{CH[(HHH)·R]C})---(7)]]>换句话说，从(3)我们可以得到limσ→0(WAy)=A-1(Re{CH[(HHH)·R]C})-1Re{CH[(HHH)·R]C}Ab]]>＝b并且多址联接干扰被迫使为零。在解相关检测器A(DD-A)的实现中，我们以WA＝A-1Re-1{M}代替方程式(6)给出的WA。我们将假设CH[(HHH)·R]C是正定的并且因此是可逆的；所以Re{CH[(HHH)·R]C}也是可逆的。
现在让我们计算这第一个MMSE检测器(MMSE A)的误码率。在限制器输入端的K向量是yA＝WARe{yconv}＝WAMAb+nA这里nA~η(0K,σ2WAM^WAH)]]>因此使用MMSE-A检测器用于第k个用户的误码率是PkMMSE-A(σ)=Q[(WAMAb)kWAM^WA(k,k)Hσ]---(8)]]>设WA＝A-1(Re{CH[(HHH)·R]C})-1并且假设是理想的信道估计，对于DD-A的误码率是PkDD-A(σ)=Q[Ak[(Re{CH[(HHH)·R]C})-1](k,k)σ]---(9)]]>这种MMSE检测器的一个缺点是信道和阵列系数的估计不能结合到用于获得WA的自适应算法中。必须使用一些其他方式(如序列信号或导频信号)明确地获得估计。这导致我们询问“可能设计一个不要求预先了解阵列和信道参数的自适应MMSE吗？”回答是“可以”，正如我们将在下面部分看到的。
在参考标号11，相关器输出由它们的相应信道估计的复数共轭加权。更准确地说，在用于用户k的第l个多径的第p个天线的相关器输出由相应信道(阵列/多径)系数估计的复数共轭
加权。使用单独的信道估计算法获得这些估计。在每个天线，加上用于第k个用户的L个分量(12)，然后用于第k个用户的生成的P个分量加到13k上(k＝1…K)。每个分量随后传送到一个判决装置或一个译码器。对于PAM数据信号，用于第k个用户(14k，k＝1…K)的判决装置给出作为它输出的码元，该码元对于第k个向量分量实数部分在欧几里德距离意义上最近。对于QAM数据信号，不要求实数运算符。
图1中，在匹配滤波器输出端的KLP-向量z＝[z1T…zpT]T可以写为z=R~HDCAb+ns---(2)]]>这里我们定义KLP×KLP矩阵
和KLP×KL矩阵，表示两个矩阵之间的直积(克罗内克积)运算，而复数噪声向量是具有下面分布的复数高斯随机向量Re(ns)Im(ns)~η(0N0N,σ2Re(R~)Im(R~)-Im(R~)Re(R~)|)]]>使用
(这里·表示两个相同尺寸矩阵之间按分量逐个作出的乘积)，在信道组合器输出端的K向量和到判决装置141-14K的输入可以简明地写成Re{yconv}=Re{(H^HC^)Hz}=MAb+nconv---(3)]]>
这里H表示复数共轭转置，·表示两个相同尺寸矩阵之间按分量逐个作出的乘积，并且
。噪声向量完全是实数并且具有分布nconv~η(0N,σ2M^)]]>这里
。如果该数据是BPSK调制的，则对于用户k的位判决是简单的Re{yconv}的第k个分量的硬限制k＝sgn(Re{yconv，k})。常规瑞克接收机的相应的误码率Pkconv是Pkconv(σ)=Q[(MAb)kσM(k,k)^]---(4)]]>这里X(k，k)是矩阵X的第(k，k)个元素。根据BPSK调制数据给出在各种检测器之间的性能比较，但是一般也可以对于任何QAM调制数据进行这些性能比较。
通用时空线性多用户瑞克接收机(图2)上述现有技术时空瑞克接收机在给定用户意义上是单用户接收机，解调仅仅使用来自该用户的信息。因为它不考虑来自其他用户的干扰的存在，它的性能受到损失。图2示出一个通用的多用户时空瑞克接收机。它使用一个信息组合器考虑多址联接干扰。给出这种检测器的两种版本。第一种称为检测器A，在线性组合器之前使用实数运算符。在这种情况下，图2中参考号20包括跟随着一个线性组合器的一组实数运算符，分量141…14k中的每一个是一个限制器，它确定对于它输出的最接近估计码元。第二种称为检测器B，使用跟随着实数运算符由复数矩阵表示的线性组合器。在这种情况下，参考号20仅仅包括线性组合器，而判决装置包括实数运算符和限制器。尽管类似，但用于检测器A和B的线性组合器导致了不同的性能和自适应实现。现在我们更充分地描述两个检测器实施例。
在参考标号11，相关器输出由它们的相应信道估计的复数共轭加权。更准确地说，在用于用户k的第l个多径的第p个天线的相关器输出由相应信道(阵列/多径)系数估计的复数共轭
加权。使用单独的信道估计算法获得这些估计。在每个天线，加上用于第k个用户的L个分量(12)，然后用于第k个用户的生成的P个分量加到13k上(k＝1…K)。每个分量随后传送到一个判决装置或一个译码器。对于PAM数据信号，用于第k个用户(14k，k＝1…K)的判决装置给出作为它输出的码元，该码元对于第k个向量分量实数部分在欧几里德距离意义上最近。对于QAM数据信号，不要求实数运算符。
图1中，在匹配滤波器输出端的KLP-向量z＝[z1T…zpT]T可以写为z=R~HDCAb+ns---(2)]]>这里我们定义KLP×KLP矩阵
和KLP×KL矩阵，表示两个矩阵之间的直积(克罗内克积)运算，而复数噪声向量是具有下面分布的复数高斯随机向量Re(ns)Im(ns)~η(0N0N,σ2Re(R~)Im(R~)-Im(R~)Re(R~)|)]]>使用
(这里·表示两个相同尺寸矩阵之间按分量逐个作出的乘积)，在信道组合器输出端的K向量和到判决装置141-14K的输入可以简明地写成Re{yconv}=Re{(H^HC^)Hz}=MAb+nconv---(3)]]>
对于异步多速率情况的MMSE检测器扩展使用迄今为止，我们已经假设这些线性多用户检测器是位同步操作。现在我们示出它们的操作是如何可以扩展到位异步情况的。首先，我们考虑在同步单速率环境中解相关检测器的操作。同样的原因应用于MMSE检测器；然而，为了教学目的，我们聚焦于解相关检测器。假定有两个用户具有各自的扩展码s1(t)和s2(t)。下面的表A示出这两个代码在时间上的关系，细的垂直线表示码元边界。接收的信号是r(t)＝A1s1(t)b1+A2s2(t)b2+n(t)，这里Ak是第k个用户的幅度，bk是第k个用户的数据位，而n(t)是加性高斯噪声。
表-A2个用户同步单速率系统
表-B2个用户异步单速率系统
表-C2个用户异步多速率系统如果用户1是所需的用户，它的对于一个给定码元间隔的解相关检测器是匹配的滤波器用于投射到s2(t)零空间的代码s1(t)。对于表-B所示的异步单速率情况，用于用户1的解相关检测器投射到由用户2代码的线性组合跨越的零空间，该代码在那个码元间隔期间重叠。如果我们设s2R(t)是以s1(t)重叠s2(t)部分的并且设s2L(t-T)是以s1(t)重叠s2(t-T)的部分，然后假设BPSK数据调制，零空间由s2R(t)+s2L(t-T)和s2R(t)-s2L(t-T)跨越。这个想法可以扩展到表-C所示的异步多速率情况。使用类似的函数定义，用于用户1的解相关器位于由s2R(t)+s2(t-T)+s2L(t-2T)、s2R(t)+s2(t-T)-s2L(t-2T)、s2R(t)-s2(t-T)+s2L(t-2T)、以及s2R(t)-s2(t-T)-s2L(t-2T)跨越的零空间。由于数据率之间的不同增加，线性多用户检测器的效率降低了，因为它将限定于一个更加特定的子空间中。上面的原因可以通过类似的组合干扰波形的适当子空间应用于MMSE检测器。
MMSE检测器和干扰对消这里描述的MMSE技术可以与称为干扰对消的非线性多用户检测技术一起使用。而线性多用户检测器依靠子空间投影减轻干扰，干扰对消器直接去掉干扰。可以使用两类干扰对消前置组合器，这里对消发生在线性组合器之前，以及后置组合器，这里对消发生在线性组合器之后。前置组合器对消用于去除没有由线性组合器考虑的干扰。例如，我们可以选择去除小区内一个高功率、高速率信号而不是以线性组合器考虑它。另一方面，后置组合器对消用于改进由线性组合器做出的码元估计。使用来自线性多用户检测器的初步码元估计，信号被重构并且从接收的信号提取出来以形成用于第二阶段码元估计的增强信号。在两种干扰对消技术中，通过抑制干扰可以潜在地增加性能。
图7示出根据图3中所示的检测器A的前置组合器干扰对消器的方块图。在到线性组合器WA的第k个输入端，存在着由于用户j＝1…K，j≠k的干扰，它将以WA考虑。然而，还可能有其他没有考虑的干扰，例如来自小区内的高功率干扰。如果已知在所有天线的扩展码、多径延迟、数据位和信道参数则这种干扰可以从到WA的输入中提取出来。在这种情况下，可能重构对于干扰的基带信号并且可以估计它对于WA的第k个输入的贡献(将重构信号通过用于用户k的处理链)以及在指定的点70k(k＝1…K)的贡献可以从对于这个的其他干扰中减去或者可以类似地计算和减去其他输入。
图8示出根据图4的检测器B的前置组合器干扰对消器的方块图。使用如上所述对于检测器A前置组合器干扰对消器的相同方法，没有由WB考虑的干扰可以在指定的点80从它的每个输入减去。
图9示出根据本发明的后置组合器干扰对消接收机的方块图。这种接收机包括P个天线并且具有三级，表示为90、91和92。第一级90接收基带信号rp并且以图1的MMSE-A检测器或图2的MMSE-B检测器对于所有K个用户进行初步码元估计。第二级91使用来自第一级90的初步码元估计并且了解K个用户的扩展码、延迟和信道参数以重构对于每个用户的基带接收信号。第二级91从接收的信号r1…rp减去相对于所需用户k的多址联接干扰以形成下面对于用户k的增强的接收信号rk,p(IC)≡rp-Σj=1,j≠kKΣl=1LAjh^j,l.pC^j,lSj,lb^jp=1,···,P]]>注意也可以去除对于用户k的多径干扰。我们将不详细地描述这种选择，因为它的相对增益一般被忽略，除非用户K的数量非常小。理想地，如果初步的码元估计和信道估计是理想的，则在增强的接收信号中没有任何多址联接干扰。然而，一般来说这不是实际的情况。在第三级92，用图1所示的常规时空接收机处理增强的接收信号以产生对于K个用户的最后码元估计k(IC)。
通过无限的重复级91和92可以继续用于多次迭代上述的后置组合器干扰对消过程。由于在第一次两个或三个迭代之后重复的迭代减弱了，可以增强性能。为了减少复杂性，这两类干扰对消可以在一个用户子集上完成。
上面已经描述的可以认为是本发明原理的说明。特别应该注意到本发明的装置和方法可以通过各种技术实现，例如通过使用大规模集成电路、专用集成电路和/或存储程序的通用计算机或专用计算机或者使用任何一种计算机可读介质的微处理器。其他的实施例对于本领域的技术人员也是很明显的，但它们都不会背离本发明的精神和范围。
权利要求
1.一种用于检测分别由多个直接序列扩谱信号调制的数据信号的方法，包括步骤a)产生匹配的滤波器输出，用于至少一个多信号的子集、多个这些信号的多径分量以及多个接收机天线；b)得到相干信道估计；c)使用相干信道估计，加权和组合匹配的滤波器输出；以及d)当抑制多址联接干扰时，对于至少一个扩谱信号子集估计对应的数据码元。
2.如权利要求1所述的方法，其中所述相干信道估计是从所述匹配的滤波器输出中得到的。
3.如权利要求1所述的方法，其中所述相干信道估计是从一个辅助信号中得到的。
4.如权利要求1所述的方法，其中通过一个线性组合器抑制多址联接干扰并且其中判定装置最优地估计给定线性组合器输出对应的数据码元。
5.如权利要求4所述的方法，其中产生匹配的滤波器输出、组合匹配的滤波器输出以及利用线性组合器的步骤使用修改的抽头加权滤波器实现。
6.如权利要求4所述的方法，其中线性组合器使用最小均方差准则。
7.如权利要求6所述的方法，其中最小均方差准则对于第k个用户指定一个线性组合器，它使组合器输出和第k个用户数据码元之间的期望平方差减至最小。
8.如权利要求7所述的方法，其中线性组合器由复数块-托普利兹相关矩阵
、对角幅度矩阵A、复数阵列/信道矩阵G以及背景噪声方差σ2形成矩阵WB。
9.如权利要求7所述的方法，其中线性组合器由复数块-托普利“三圆环轮式磁能动力体”演化而来。定、转子轮均是中间板式磁体结构，对于转子轮从板的
剖面来看，无论长、短轴永磁体，在板两侧均成“八”字形对称布设；而定子轮板两侧的永磁块则均成倒“八”字形对称布设，两端盖所带半个定子轮永磁块是半个倒“八”字形，即保证气隙6两侧永磁块N磁极平行相对应，正是由于气隙6两侧定转子永磁块N磁极相对于各自板斜向约45°左右角度平行相对，使磁斥力作用线沿长轴或短轴轴线作用于转子轮板中心面内，“八”字形作用线夹角约为90°左右，据力的合成定理，合力作用线正是“八”字形夹角的角分线，该角分线与转子轮板中心面重合，与合力作用交点和圆心的连线垂直，换句话说，合力作用线的方向就是合力作用点的切线方向，转子轮中各对永磁块受两侧定子轮永磁块磁斥力的作用线交点的切线方向相同，各转子轮中磁斥力作用线交点的切线方向也相同，从而推动各转子轮同步向一个方向转动。</p><p>周样由图12看到，定转子间气隙21相等且在保证动平衡条件下尽量近。转子轮5与机壳1内周的间隙也相等。三个转子轮5是一样的，两个定子轮6也是相同的，但三个转子轮中，起始磁块端点17与起始纵轴线18的距离是不相等的，这主要是因为磁互相作用力是短程力，为了保证转子组受到的磁推动力是连续的，就需要三个转子按“接力”式排列成递推组合结构，为此用A、B、C三个剖面即图13、14、15、16、17、18、19分别展示三个转子轮中的磁块排列情况。其中图13展示了左转子轮它是由与轴连接的花键圆盘12、幅条13，轮径向中心圆环板14，内周圆环形磁轭圈15、外周圆环磁轭圈16组成轮基体，在“C型”磁轭槽中，圆环板14两侧布设对称的永磁块10和11，其中永磁块11是转子轮右侧、定子轮左侧、右端盖轮通用的，永磁块10是转子轮左侧、定子轮右侧、左端盖轮通用的，并且永磁块10和11在板14两侧成“八”字形布设，两磁块长轴延长线交于板中心面内一点，延长线所夹锐角为90°左右，而磁块长轴与板14成45°夹角。本例中左转轮起始磁块端点17与起始纵轴线18的左偏距是3.5毫米，主视图中10所示为轮式磁能动力机转子轮永磁块在轮主视图中的标准画法，20是其简易画法，后同，不累述。</p><p>图15和16及图17展示了中转子轮，主视图中起始磁块端点17与起始纵轴“三圆环轮式磁能动力体”演化而来。定、转子轮均是中间板式磁体结构，对于转子轮从板的
剖面来看，无论长、短轴永磁体，在板两侧均成“八”字形对称布设；而定子轮板两侧的永磁块则均成倒“八”字形对称布设，两端盖所带半个定子轮永磁块是半个倒“八”字形，即保证气隙6两侧永磁块N磁极平行相对应，正是由于气隙6两侧定转子永磁块N磁极相对于各自板斜向约45°左右角度平行相对，使磁斥力作用线沿长轴或短轴轴线作用于转子轮板中心面内，“八”字形作用线夹角约为90°左右，据力的合成定理，合力作用线正是“八”字形夹角的角分线，该角分线与转子轮板中心面重合，与合力作用交点和圆心的连线垂直，换句话说，合力作用线的方向就是合力作用点的切线方向，转子轮中各对永磁块受两侧定子轮永磁块磁斥力的作用线交点的切线方向相同，各转子轮中磁斥力作用线交点的切线方向也相同，从而推动各转子轮同步向一个方向转动。</p><p>周样由图12看到，定转子间气隙21相等且在保证动平衡条件下尽量近。转子轮5与机壳1内周的间隙也相等。三个转子轮5是一样的，两个定子轮6也是相同的，但三个转子轮中，起始磁块端点17与起始纵轴线18的距离是不相等的，这主要是因为磁互相作用力是短程力，为了保证转子组受到的磁推动力是连续的，就需要三个转子按“接力”式排列成递推组合结构，为此用A、B、C三个剖面即图13、14、15、16、17、18、19分别展示三个转子轮中的磁块排列情况。其中图13展示了左转子轮它是由与轴连接的花键圆盘12、幅条13，轮径向中心圆环板14，内周圆环形磁轭圈15、外周圆环磁轭圈16组成轮基体，在“C型”磁轭槽中，圆环板14两侧布设对称的永磁块10和11，其中永磁块11是转子轮右侧、定子轮左侧、右端盖轮通用的，永磁块10是转子轮左侧、定子轮右侧、左端盖轮通用的，并且永磁块10和11在板14两侧成“八”字形布设，两磁块长轴延长线交于板中心面内一点，延长线所夹锐角为90°左右，而磁块长轴与板14成45°夹角。本例中左转轮起始磁块端点17与起始纵轴线18的左偏距是3.5毫米，主视图中10所示为轮式磁能动力机转子轮永磁块在轮主视图中的标准画法，20是其简易画法，后同，不累述。</p><p>图15和16及图17展示了中转子轮，主视图中起始磁块端点17与起始纵轴iv)用于形成K个生成和的K向量的装置；b)应用一个线性组合器到K向量以抑制多址联接干扰的装置；以及c)用于处理生成信号以产生各个数据码元估计的装置。
26.如权利要求25所述的接收机，其中用于线性组合器的所述装置包括形成下面各项乘积的装置对角幅度矩阵(A)；等效的阵列/信道/代码相关矩阵(M)；以及矩阵(
)的逆矩阵，这里MH是阵列/信道/代码相关矩阵的厄米特共轭，(
)是估计的阵列/信道/代码相关矩阵，而(σ2)是背景噪声方差。
27.如权利要求25所述的接收机，其中在线性组合器之前提取K个生成和的实数部分。
28.接收机包括一组相关器，用于去扩展承载各自数据b1，b2，…，bK的K个扩展代码数据信号s1，s2，…，sK，每个具有最多L个可分辨多径分量以形成一个K向量；应用一个线性组合器到K向量以抑制多址联接干扰的装置；以及用于处理线性组合器的输出以产生各个数据码元估计的装置。
29.如权利要求28所述的接收机，其中用于线性组合器的所述装置包括形成下面各项的乘积对角幅度矩阵(A)；阵列/信道矩阵的厄米特转置；以及矩阵(
)的逆矩阵，这里
是复数块-托普利兹相关矩阵，σ2是背景噪声方差，而IKLP是维数为KLP×KLP的单位矩阵。
30.一种解调在一个P个天线阵列接收的、经受L路径的频率选择性衰落的K个直接序列扩展代码数据信号的方法包括以匹配于在P个天线的每一个接收的K个信号的每一个的L个多径延迟复制的滤波器组去扩展；以一个阵列/信道系数估计的复数共轭加权去扩展信号；求和加权的、去扩展的信号以形成一个K向量；应用一个线性转换到K向量以抑制多址联接干扰；以及处理生成的信号以产生各个数据码元的估计。
31.如权利要求30所述的方法，其中所述线性转换是下面各项的乘积对角幅度矩阵(A)；等效的阵列/信道/代码相关矩阵(M)；以及矩阵(
)的逆矩阵，这里MH是阵列/信道/代码相关矩阵的厄米特共轭，(
)是估计的阵列/信道/代码相关矩阵，而(σ2)是背景噪声方差。
32.一种解调在一个P个天线阵列接收的、经受L路径的频率选择性衰落的K个直接序列扩展代码数据信号的方法包括a)以匹配于在P个天线的每一个接收的K个信号的每一个的L个多径延迟复制的滤波器组去扩展以产生一个K向量；b)应用一个线性转换到K向量以组合匹配的滤波器输出并且抑制多址联接干扰；以及c)处理生成的信号以产生各个数据码元的估计。
33.如权利要求32所述的方法，其中所述用于线性组合器的装置包括下面各项的乘积d)对角幅度矩阵(A)；e)阵列/信道矩阵的厄米特转置(GH)；以及f)矩阵(
)的逆矩阵，这里
是复数块-托普利兹相关矩阵，σ2是背景噪声方差，而IKLP是维数为KLP×KLP的单位矩阵。
34.一种使用直接序列扩谱信号的接收机包括一个具有阵列处理的瑞克接收机；以及一个接收所述瑞克接收机输出的多用户检测器。
35.如权利要求34所述的接收机，其中具有阵列处理的所述瑞克接收机包括一组匹配于多个扩谱信号的多径定时延迟的滤波器；以及一个用相应阵列/信道估计加权匹配滤波器输出的电路。
36.一个用于检测通过多个直接序列扩谱信号分别调制的数据信号的系统，包括多个空间分集天线；与所述天线有关的多个匹配滤波器，用于对至少一个多信号子集、这些信号的多个多径分量以及多个接收机天线提供多个输出(z1…zp)；一个用于从匹配的滤波器输出获得相干信道估计(
)的电路；一个第二电路，使用相干信道估计用于加权和组合匹配的滤波器输出；以及一个第三电路，当抑制多址联接干扰时用于估计至少一个扩谱信号子集的各自数据码元。
37.一个包括执行软件的处理器的接收机，所述软件是在计算机可读介质中并且用于检测通过多个直接序列扩谱信号分别调制的数据信号，所述软件使得所述处理器完成至少下面的步骤a)产生匹配的滤波器输出，用于至少一个多信号子集、这些信号的多个多径分量以及多个接收机天线；b)从匹配的滤波器输出得到相干信道估计；c)使用相干信道估计加权和组合匹配的滤波器输出；以及d)当抑制多址联接干扰时估计用于至少一个扩谱信号子集的各个数据码元。
全文摘要
本发明的线性时空多用户检测器结合了阵列处理、瑞克检测和多用户检测并且满足最小均方差准则。两个实施例都可以自适应地实现以便考虑未知的干扰源。线性时空检测器的第一个实施例的自适应实现要求信道系数的精确估计。另一方面,第二个实施例的自适应实现不要求这些估计,因为它们可以自适应地获得。这两种检测器都提供明显的性能优点,它们超过单用户时空检测器和超过不考虑来自未知源的CDMA干扰的多用户检测器。
文档编号H04B1/707GK1247417SQ9911187
公开日2000年3月15日 申请日期1999年8月3日 优先权日1998年8月4日
发明者霍华德·C·黄, 劳伦斯·欧格内·梅伦德尔, 康斯坦迪诺斯·巴西尔·巴巴迪亚斯 申请人:朗迅科技公司