专利名称：一种在多天线数字无线通信系统中检测信号的方法
技术领域：
本发明涉及信号检测技术，特别是指一种在多天线数字无线通信系统中检测信号的方法。
背景技术：
根据信息论，在通信系统的发射端和接收端同时使用多天线阵列可以极大的提高传输比特率。
在发射端和接收端同时使用多天线阵列的具有空-时架构的无线通信系统如图1所示。该系统工作在瑞利散射环境，信道矩阵的各个元素可以近似看作是统计独立的。在图1所示的系统中，一个数据序列分成M个不相关的码元子序列，每个子序列由M个发射天线的一个发射。M个子序列在经过一个信道矩阵为H的信道的影响后，在接收端由N个接收天线接收。发射信号s1，...，sM分别通过M个不同的天线单元a-1，...，a-M发射，相应的接收信号x1，...，xN分别从N个不同的天线单元b-1，...，b-N接收。该系统中，发射天线单元数M最少是2，而接收天线单元数N最少是M。信道矩阵H是一个N×M的矩阵，矩阵中第i行j列的元素表示第i个接收天线和第j个发射天线通过传输信道的耦合。接收信号x1，...，xN在数字信号处理器中被处理以产生恢复的发射信号 此图中也显示了求和成分c-1，c-2，...，c-N，它们代表包含的无法避免的噪声信号w1，w2，...，wN，这些噪声信号分别加入到接收天线单元b-1，b-2，...，b-N接收到的信号中。
由发射天线和接收天线之间的信道系数组成的矩阵为信道矩阵，信道系数是利用接收信号进行信道估计得到。在图1所示系统中的信道矩阵H是一个N×M的矩阵，表示为
信道矩阵H是一个N×M复数矩阵，假定该信道矩阵在K个符号的时期内是常数。向量hn:(n＝1，2，...，N)和h:m(m＝1，2，...，M)的长度分别是M和N。其中，信道矩阵H包含的信道向量h:1至h:M分别表示信道对M个传输信号中每个传输信号的影响。更明确的，信道向量h:m(m＝1，2，...，M)包括信道矩阵项h1m至hNm，分别表示在接收天线单元b-1至b-N中每个接收天线上，信道对发射信号sm的影响。
在图1所示的系统中，发射信号的向量与接收信号的向量之间满足关系式x(k)=Σm=1Mh:msm(k)+w(k)=Hs(k)+w,]]>其中k表示采样时刻，k＝1，2，...，K。用向量形式表示上述关系为x→=Σm=1Mh:msm+w→=Hs→+w→,]]>再把该式写为x→=s1·h:1+s2·h:2+...+sm·h:m+···+sM·h:M+w→]]>的形式，可以清楚的看到各个发射信号对接收信号向量 的影响。
发射信号的最小均方误差(MMSE)估计为s→^=(HH·H+αIM×M)-1HHx→,]]>其中，符号-1表示求矩阵的逆矩阵，α为与发射信号的信噪比相关的常数，α=σw2σs2.]]>本发明所述的估计误差e=s→-s→^]]>的协方差矩阵为把加性高斯白噪声的方差归一化为1的情况下的协方差矩阵，即E{(s→-s→^)(s→-s→^)H}=(HH·H+αIM×M)-1,]]>记为Q，并定义R＝(HH·H+αIM×M)，则有Q＝R-1。从而，对发射信号的估计值可以表示为s→^=QHHx→]]>的形式。
在现有技术中，在以上所述的MIMO系统中实现信号检测的方法是首先利用信道矩阵H得到所有待检测发射信号的估计误差协方差矩阵Q的初始值，然后利用所得到的Q的初始值计算对发射信号的估计值。其中，如何计算Q的初始值将会影响信号检测的计算量和复杂度。下面介绍现有技术中的具体实现方法。
在现有技术中计算Q的初始值是通过以下所述的递推方法计算首先设置一个M×M的对角线上元素为(1/α)的对角矩阵，即Q
＝(1/α)·IM×M，其中的Q
作为递推Q的初始值的第一个递推变量。
然后，按照Q[l]=Q[l-1]-Q[l-1]hl:hl:HQ[l-1]1+hl:HQ[l-1]hl:,]]>l＝1，...，M的公式递推得到所有待检测发射信号对应的估计误差协方差矩阵Q的初始值。
利用Q的初始值计算对发射信号的估计值的方法为假设m为检测信号过程中的迭代变量，则检测第m个发射信号时的M+1-m个待检测发射信号对应的估计误差协方差矩阵记为QM+1-m。
当m＝1时，即当检测第一个发射信号时，此时的QM+1-m为Q的初始值，因此直接利用Q的初始值计算得到第一个被检测发射信号的估计值。然后，利用Q递推得到剩下的M-1个待检测发射信号对应的估计误差协方差矩阵QM-1。
当m＞1时，此时对应M+1-m个待检测发射信号对应的估计误差协方差矩阵QM+1-m已在上一次检测发射信号后计算得到，并利用QM+1-m计算得到第m个被检测发射信号的估计值。然后，利用QM+1-m递推得到剩下的M-m个待检测发射信号对应的估计误差协方差矩阵QM-m。
综上所述，在信号检测时需要的Q的初始值的递推过程中的矩阵从一开始就是M×M的矩阵，由于矩阵的计算复杂度与矩阵维数直接相关，当发射天线个数M大时，递推的计算复杂度会变得更高。而且，在信号检测过程中，每一次都要重新计算当前待检测发射信号对应的估计误差协方差矩阵，这必然导致信号检测计算复杂度的提高。
发明内容
有鉴于此，本发明的主要目的在于提供一种在多天线数字无线通信系统中检测信号的方法，使得检测信号的计算复杂度减少。
为了达到上述目的，本发明提供一种在多天线数字无线通信系统中检测信号的方法，在多入多出MIMO系统中检测至少两个发射信号，所述发射信号由发射端各个不同的发射天线分别发射并经过一个信道到达接收端，该方法包括a.接收端的至少两个接收天线接收所述发射信号，获得至少两个接收信号；b.接收端根据接收信号进行信道估计，得到由发射天线和接收天线之间的信道系数组成的信道矩阵H；c.利用信道矩阵H计算出所有发射信号中的部分发射信号的估计误差协方差矩阵，然后利用信道矩阵H和计算出的部分发射信号的估计误差协方差矩阵，递推求得包括所述部分发射信号且个数多于所述部分发射信号个数的发射信号的估计误差协方差矩阵；d.利用步骤c所得到的估计误差协方差矩阵，检测步骤c中所述的包括部分发射信号且个数多于所述部分发射信号个数的发射信号。
步骤c所述递推求得包括所述部分发射信号且个数多于所述部分发射信号个数的发射信号的估计误差协方差矩阵为递推求得所有发射信号的估计误差协方差矩阵；步骤d所述检测包括部分发射信号且个数多于所述部分发射信号个数的发射信号为检测所有发射信号。
步骤c所述利用部分发射信号的估计误差协方差矩阵递推求得所有发射信号的估计误差协方差矩阵的步骤为以部分发射信号的估计误差协方差矩阵的Sherman-Morrison结果作为子矩阵，递推求得所有发射信号的估计误差协方差矩阵。
所述步骤b和步骤c之间进一步包括设置检测至少两个发射信号的先后顺序；所述步骤c包括c11.利用所设置的检测顺序中最后被检测的第一数目个发射信号对应的信道矩阵，计算所述最后被检测的第一数目个发射信号的估计误差的协方差矩阵；设置用于递推检测发射信号估计误差协方差矩阵大于所述第一数目的第二数目；c12.利用所设置的检测顺序中最后被检测的第二数目个发射信号对应的信道矩阵，并以上一次递推或步骤c11所得到的最后被检测的第一数目个发射信号的估计误差协方差矩阵的Sherman-Morrison结果作为子矩阵，递推计算所述最后被检测的第二数目个发射信号的估计误差协方差矩阵，如果已得到所有被检测发射信号的估计误差协方差矩阵，则结束本流程，否则令第一数目的值等于第二数目的值后，第二数目的值加1或大于1的整数值，返回步骤c12。
所述步骤d包括d1.在待检测发射信号中选择当前被检测的一个发射信号，利用步骤c得到的所有发射信号的估计误差协方差矩阵、信道矩阵H以及接收信号得到对所述当前被检测的一个发射信号的估计值；d2.利用步骤d1得到的当前被检测的一个发射信号的估计值计算对检测后续待检测发射信号的干扰值，消除所述当前被检测的一个发射信号对检测后续待检测发射信号的干扰；d3.重复步骤d1、d2，直到检测到所有待检测发射信号。
所述步骤d1之前进一步包括利用信道矩阵H对接收信号进行预匹配滤波变换；计算信道矩阵H的互相关信道矩阵Φ，Φ＝HH·H；步骤d1所述得到对当前被检测的一个发射信号的估计值的步骤包括利用待检测发射信号的估计误差协方差矩阵和所述接收信号的预匹配滤波结果得到所述当前被检测的一个发射信号的估计值；所述步骤d2包括利用所述当前被检测的一个发射信号的估计值和信道矩阵H的互相关信道矩阵Φ计算已检测的发射信号对检测后续发射信号的干扰值，并从所述接收信号的预匹配滤波结果中消除已检测的发射信号的干扰，得到修正的接收信号的预匹配滤波结果，作为下一次检测信号时的接收信号的预匹配滤波结果。
根据本发明提供的检测信号方法，在计算所有发射信号的估计误差协方差矩阵的初始值时，预先设置在接收端检测发射信号的顺序，并按照所述顺序依次递推从少到多的发射信号对应的估计误差协方差矩阵，最终得到所有发射信号的估计误差协方差矩阵的初始值，所述的递推方法，减少了递推的计算量，从而能够减少检测信号的计算复杂度。在初始值的递推过程中得到的中间结果还可以在检测信号过程中利用，如果预先设置的检测顺序与实际检测顺序相近或相同，则在检测信号过程中减少重新计算待检测发射信号对应的估计误差协方差矩阵的步骤，从而能够减少检测信号的计算复杂度。
图1所示为现有技术中多天线数字无线通信系统框图；图2所示为本发明中计算所有发射信号的估计误差协方差矩阵的初始值的流程图；图3所示为本发明中检测信号的流程图。
具体实施例方式
为使本发明的目的、技术方案和优点更加清楚明白，下面举具体实施例，对本发明作进一步详细的说明。
本发明在检测信号过程中，计算所有发射信号的估计误差协方差矩阵的初始值的主要思想为按照从少到多的发射信号的顺序，先计算得到较少发射信号的估计误差协方差矩阵，然后利用所得到的较少发射信号的估计误差协方差矩阵，递推计算较多发射信号的估计误差协方差矩阵。通过以上递推方法计算估计误差协方差矩阵的初始值，不仅能够减少递推初始值的计算量，而且在检测信号的步骤中还可以充分利用在计算初始值过程中所得到的中间结果，从而进一步减少检测信号的计算量。
下面给出具体实现方法。图2所示为计算所有发射信号的估计误差协方差矩阵的初始值的流程图，包括以下几个步骤步骤201接收端接收到发射端从M个发射天线分别发射的M个发射信号，获得N个接收信号，并根据接收信号进行信道估计，得到由发射天线和接收天线之间的信道系数组成的信道矩阵H。
预先设置所有M个发射信号在接收端被检测的先后顺序，用发射信号序号记为tM，tM-1，...，tm，...，t2，t1。相应的，把信道矩阵H按列重新排序，得到按照所述预先设置的检测发射信号的先后顺序排序的信道矩阵H，记为HM(tM)，HM(tM)=h:t1h:t2...h:tM-1h:tM.]]>利用向量f＝[t1，t2，...，tm，...，tM-1，tM]T记录与信道矩阵HM(tM)对应的发射信号的索引。
步骤202用信道矩阵HM(tM)，先求得HM(tM)的互相关信道矩阵ΦM=(HM(tM))H·HM(tM),]]>再由ΦM求得估计误差的协方差矩阵QM的逆矩阵RM=(HM(tM))H·HM(tM)+αIM×M=ΦM+αIM×M.]]>其中，RM=(HM(tM))H·HM(tM)+αIM×M]]> 其中，*表示对1个复数取共轭。
步骤203计算最后被检测的一个发射信号t1对应的估计误差协方差矩阵，记为Q1(t1)。
对应发射信号t1的信道矩阵为H1(t1)=[h:t1].]]>从步骤202得到的RM中，得到发射信号t1的估计误差协方差矩阵的逆矩阵为R1(t1)=(h:t1)H·h:t1+α=rt1t1,]]>其中，rt1t1就是RM的第1行第1列的元素。
然后利用R1(t1)计算对应最后被检测的发射信号t1的估计误差的协方差矩阵Q1(t1)。由Q1(t1)=(R1(t1))-1]]>得到Q1(t1)=1/rt1t1.]]>下面让递推Q的初始值的变量m设为2，进入步骤204。在下面的递推Q的初始值的步骤中，将最后被检测的m个发射信号tm，...，t2，t1对应的估计误差协方差矩阵记为Qm(tm)。
步骤204判断是否已得到所有被检测发射信号的估计误差协方差矩阵，即判断m是否大于M，如果是，则说明已得到M个被检测发射信号的估计误差协方差矩阵，转到步骤208；否则，递推m个被检测发射信号的估计误差协方差矩阵Qm(tm)的值，执行步骤205、206、207。
步骤205最后被检测的m个发射信号tm，...，t2，t1对应的信道矩阵为Hm(tm)=h:t1h:t2...h:tm,]]>相应的，对应最后被检测的m个发射信号tm，...，t2，t1的估计误差协方差矩阵的逆矩阵为Rm(tm)=(Hm(tm))H·Hm(tm)+αI(m)×(m).]]>Rm(tm)与Rm-1(tm-1)有如下的递推关系Rm(tm)=Rm-1(tm-1)vm-1(tm)(vm-1(tm))Hβ1(tm),]]>其中，Rm-1(tm-1)是上一次递推得到的对应最后被检测的m-1个发射信号tm-1，...，t2，t1的估计误差协方差矩阵的逆矩阵或者是对应最后被检测的1个发射信号t1的估计误差的协方差矩阵的逆矩阵R1(t1)；β1(tm)=h:tmH·h:tm+α=rtmtm;vm-1(tm)=h:t1H·h:tmh:t1H·h:tm...h:tm-1H·h:tm=rt1tmrt2tm...rtm-1tm.]]>容易看到β1(tm)和vm-1(tm)都可以从步骤202中计算的RM中直接得到，更具体的，是RM的第m行第m列的元素，而vm-1(tm)是由RM的第m列的头m-1项组成。从而不需要任何计算，就可以直接得到Rm(tm)。
步骤206求最后被检测的m个发射信号tm，...，t2，t1对应的估计误差协方差矩阵Qm(tm)。
利用对应最后被检测的一个发射信号的估计误差协方差矩阵Q1(t1)或者上一次递推得到的对应最后被检测的m-1个发射信号tm-1，...，t2，t1的估计误差协方差矩阵Qm-1(tm-1)，递推得到Qm(tm)，递推方法如下所述首先计算Qm-1(tm-1)的Sherman-Morrison结果，即使用Sherman-Morrison公式得到(Tm-1(tm1))-1，(Tm-1(tm))-1=Qm-1(tm-1)+Qm-1(tm-1)vm-1(tm)(vm-1(tm))HQm-1(tm-1)β1(tm)-(vm-1(tm))HQm-1(tm-1)vm-1(tm);]]>然后由(Tm-1(tm)-1、vm-1(tm)，以及β1(tm)得到最后被检测的m个发射信号tm，...，t2，t1对应的估计误差协方差矩阵Qm(tm)为，Qm(tm)=(Tm-1(tm))-1-(Tm-1(tm))-1vm-1(tm)β1(tm)-((Tm-1(tm))-1vm-1(tm)β1(tm))H1β1(tm)+(vm-1(tm))H(Tm-1(tm))-1vm-1(tm)(β1(tm))2.]]>步骤207m的值增加1，即m＝m+1，然后返回步骤204。
步骤208得到所有M个发射信号tM，tM-1，...，tm，...，t2，t1对应的估计误差协方差矩阵QM(tM)的值。QM(tM)就是信号检测过程中决定一个最优检测顺序，且依照所述最优检测顺序并使用干扰消除的方法逐次检测各个发射信号时，所使用的估计误差协方差矩阵的初始值，所述初始值记为QM，QM=QM(tM).]]>通过以上步骤，得到了所有被检测发射信号的估计误差协方差矩阵的初始值QM；同时，还得到了检测部分发射信号的估计误差协方差矩阵Qm，其中，m＝1，2，...，M-1，Qm对应发射天线tm，...，t2，t1的估计误差协方差矩阵。
得到所有待检测发射信号的估计误差协方差矩阵的初始值后，进入图3所示的检测信号的流程中，即转到图3的a。
图3为信号检测的流程图，图3所示的信号检测从a开始。
步骤301对接收信号进行预匹配滤波变换，得到接收信号的预匹配滤波结果zM=(HM)H·x→,]]>其中，(HM)H为匹配滤波器，向量 为表示接收信号x1，...，xN的向量。发射信号的索引仍然是向量f＝[t1，t2，...，tm，...，tM-1，tM]T。
下面将检测信号过程中的迭代变量m设为1，进入步骤302。在下面的检测信号的迭代步骤中，在M+1-m个发射信号中检测某一发射信号时，对应M+1-m个待检测发射信号的估计误差协方差矩阵记为QM+1-m。
步骤302判断是否检测最后一个发射信号，即判断m是否等于M，如果是，则执行步骤313；否则执行步骤303。
步骤303在M+1-m个发射信号中确定接收信噪比最好的发射信号，方法是从QM+1-m中查找对角元素最小的行记为第lm行，lm=argminiqM+1-m,ii,]]>其中，qM+1-m，ii是矩阵QM+1-m的第i行第i列的元素即对角线上的元素，i＝1，...，M+1-m。所述的第lm行对应于M+1-m个发射信号中接收信噪比最好的信号，即当前被检测的发射信号。
步骤304在矩阵QM+1-m中交换第lm行和M+1-m行，交换第lm列和M+1-m列；相应的，在矩阵RM+1-m中交换第lm行和M+1-m行，交换第lm列和M+1-m列；相应的，在矩阵ΦM+1-m中交换第lm行和M+1-m行，交换第lm列和M+1-m列；相应的，在对接收信号的预匹配滤波结果向量zM+1-m中交换第lm项和第M+1-m项。在向量f中交换第lm项和第M+1-m项。
步骤305当前被检测的发射信号为向量f中第M+1-m项，记为pm，pm＝f(M+1-m)。计算对当前被检测的发射信号的估计值ypm。发射信号pm的估计值ypm为，ypm=qM+1-mH·zM+1-m.]]>其中，qM+1-m表示QM+1-m的第M+1-m列。
步骤306对所得到的发射信号的估计值ypm进行量化，得到对发射信号的检测结果 步骤307从接收信号的预匹配滤波结果中消除当前检测到的发射信号的影响，得到尚未被检测的M-m个发射信号对应的多个接收信号的预匹配滤波结果，即从对接收信号的预匹配滤波结果向量zM+1-m中删除第M+1-m项得到具有M-m项的(zM+1-m)minus；从(zM+1-m)minus中消除当前被检测到的发射信号的干扰，得到对应所有M-m个未被检测的发射信号的多个接收信号的预匹配滤波结果 其中M-m+1是矩阵ΦM-m+1的第M+1-m列的头M-m行，M-m+1具有M-m项。
步骤308从矩阵ΦM+1-m中删除第M+1-m列和第M+1-m行，得到矩阵ΦM-m；从矩阵RM+1-m中提取头M-m行头M-m列得到矩阵RM-m，从矩阵RM+1-m中提取第M+1-m列的头M-m行得到向量vM-m，从矩阵RM+1-m中提取第M+1-m行第M+1-m列的元素得到一项βpm。
步骤309判断在图2所示求QM的初始值的过程中，是否已得到对应M-m个发射信号的估计误差协方差矩阵QM-m，如果是，则执行步骤310；否则，执行步骤311。
步骤310将递推QM的初始值的过程中得到的QM-m作为下一次迭代所需要的对应M-m个发射信号的估计误差协方差矩阵QM-m。然后执行步骤312。
步骤311利用矩阵QM+1-m计算下一次迭代所需要的对应M-m个发射信号的估计误差协方差矩阵QM-m，计算方法是从矩阵QM+1-m中删除第M+1-m行和第M+1-m列得到矩阵TM-m-1，然后利用TM-m-1、vM-m以及βpm得到QM-m，QM-m=TM-m-1-TM-m-1vM-mvM-mHTM-m-1βpm+vM-mHTM-m-1vM-m.]]>然后执行步骤312。
步骤312m的值增加1，即m＝m+1，返回步骤302。
步骤313当前被检测的发射信号为向量f的最后一项，记为pM。计算最后被检测的发射信号的估计值ypM。发射信号pM的估计值ypM为，ypM=q1H·z1.]]>其中，q1就是Q1。
步骤314对所得到的发射信号的估计值ypM进行量化，得到对发射信号的检测结果 通过以上方法，对发射信号的检测顺序为p1，p2，...，pM，对应的，发射信号的检测结果为 以上所述实施例中可以看到，如果递推发射信号估计误差协方差矩阵Q的初始值时所预先设置的检测发射信号的先后顺序与实际检测发射信号的顺序相同，则不用如步骤311所述的方法计算得到下一次迭代所需的估计误差协方差矩阵，而直接从递推Q的初始值过程中所得到的中间结果得到下一次迭代所需的估计误差协方差矩阵，因此，能够减少很多计算量；或者，如果递推Q的初始值时所预先设置的检测发射信号的先后顺序与实际检测发射信号的顺序相近，也能够省略所述步骤311，因此，也能够减少不少计算量。
从而，在慢衰落信道中，将接收端检测M个发射天线的先后顺序可以设置为最近一次检测的最优顺序，使得在图3所示的检测信号过程中，减少所述步骤311带来的计算量。因为在慢衰落信道中，信道特性变化缓慢，与最近一次的最优检测顺序相比，当前时刻的最优检测顺序变化不大或相同，因此能够很好的利用递推Q的初始值过程中所得到的中间结果减少检测信号的计算量。
在快衰落信道中，也有众多的现有技术，可以通过信道矩阵H，估计一个检测顺序，使得这个检测顺序接近实际使用的最优检测顺序。
在某些应用中，检测所有发射信号的顺序是预先固定的，依照这个预先固定的检测顺序逐个检测发射信号，在此过程中不需要去求最优的检测顺序。在这种情况下，将所述步骤201中预先设置的所有M个发射信号在接收端被检测的先后顺序，就是所述的预先固定的检测顺序，这样，不需要如步骤303所述的确定检测顺序的过程，而是根据所述预先固定的检测顺序确定当前选择哪一个待检测信号加以检测；从而，不需要如步骤304所述的矩阵的行和列的交换；同时，不需要如步骤311所述的计算下一次迭代所需的估计误差协方差矩阵的过程，直接利用递推Q过程中所得到的中间结果检测下一个待检测发射信号。通过以上方法，能够减少很多计算量。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
权利要求
1.一种在多天线数字无线通信系统中检测信号的方法，在多入多出MIMO系统中检测至少两个发射信号，所述发射信号由发射端各个不同的发射天线分别发射并经过一个信道到达接收端，其特征在于，a.接收端的至少两个接收天线接收所述发射信号，获得至少两个接收信号；b.接收端根据接收信号进行信道估计，得到由发射天线和接收天线之间的信道系数组成的信道矩阵H；c.利用信道矩阵H计算出所有发射信号中的部分发射信号的估计误差协方差矩阵，然后利用信道矩阵H和计算出的部分发射信号的估计误差协方差矩阵，递推求得包括所述部分发射信号且个数多于所述部分发射信号个数的发射信号的估计误差协方差矩阵；d.利用步骤c所得到的估计误差协方差矩阵，检测步骤c中所述的包括部分发射信号且个数多于所述部分发射信号个数的发射信号。
2.根据权利要求1所述的方法，其特征在于，步骤c所述递推求得包括所述部分发射信号且个数多于所述部分发射信号个数的发射信号的估计误差协方差矩阵为递推求得所有发射信号的估计误差协方差矩阵；步骤d所述检测包括部分发射信号且个数多于所述部分发射信号个数的发射信号为检测所有发射信号。
3.根据权利要求2所述的方法，其特征在于，步骤c所述利用部分发射信号的估计误差协方差矩阵递推求得所有发射信号的估计误差协方差矩阵的步骤为以部分发射信号的估计误差协方差矩阵的Sherman-Morrison结果作为子矩阵，递推求得所有发射信号的估计误差协方差矩阵。
4.根据权利要求2所述的方法，其特征在于，所述步骤b和步骤c之间进一步包括设置检测至少两个发射信号的先后顺序；所述步骤c包括c11.利用所设置的检测顺序中最后被检测的第一数目个发射信号对应的信道矩阵，计算所述最后被检测的第一数目个发射信号的估计误差的协方差矩阵；设置用于递推检测发射信号估计误差协方差矩阵大于所述第一数目的第二数目；c12.利用所设置的检测顺序中最后被检测的第二数目个发射信号对应的信道矩阵，并以上一次递推或步骤c11所得到的最后被检测的第一数目个发射信号的估计误差协方差矩阵的Sherman-Morrison结果作为子矩阵，递推计算所述最后被检测的第二数目个发射信号的估计误差协方差矩阵，如果已得到所有被检测发射信号的估计误差协方差矩阵，则结束本流程，否则令第一数目的值等于第二数目的值后，第二数目的值加1或大于1的整数值，返回步骤c12。
5.根据权利要求2所述的方法，其特征在于，所述步骤b和步骤c之间进一步包括设置检测至少两个发射信号的先后顺序；所述步骤c包括c21.利用所设置的检测顺序中最后被检测的一个发射信号对应的信道矩阵，计算所述最后被检测的一个发射信号的估计误差的协方差矩阵；c22.利用与所设置的检测顺序中最后被检测的m个发射信号对应的信道矩阵，并以上一次递推或者步骤c21得到的最后被检测的m-1个发射信号的估计误差协方差矩阵的Sherman-Morrison结果作为子矩阵，递推计算所述最后被检测的m个发射信号的估计误差协方差矩阵，如果已得到所有发射信号的估计误差协方差矩阵，则结束本步骤，否则m的值加1，返回步骤c22；其中，m的初始值设为2。
6.根据权利要求5所述的方法，其特征在于，所述发射天线有M个，接收天线有N个，所述M个不同的发射天线分别发射M个发射信号；所述设置检测至少两个发射信号的先后顺序的步骤包括对发射天线发射的M个发射信号重新排序得到所述先后顺序，用发射信号的序列表示为tM，tM-1…，tm，…，t2，t1；所述步骤c21包括利用发射信号t1对应的信道矩阵[h:t1]得到该发射信号的估计误差协方差矩阵的逆矩阵R1(t1)，并根据发射信号t1的估计误差协方差矩阵与R1(t1)矩阵满足的关系得到发射信号t1的估计误差协方差矩阵；所述步骤c22包括利用m个发射信号t1…tm对应的信道矩阵[h:t1h:t2…h:tm]得到所述m个发射信号的估计误差协方差矩阵的逆矩阵Rm(tm)中不包含在Rm-1(tm-1)中的部分项，并根据所述m个发射信号的估计误差协方差矩阵与Rm(tm)满足的关系，以及利用所得到的m-1个发射信号t1…tm-1的估计误差协方差矩阵，递推m个发射信号的估计误差协方差矩阵，如果已得到M个被检测信号的估计误差协方差矩阵，则结束本步骤；否则m的值加1，返回执行步骤c22；其中，h:ti表示信道矩阵H中与发射信号ti对应的列向量，i＝1…M。
7.根据权利要求6所述的方法，其特征在于，步骤c21所述的R1(t1)=(h:t1)H·h:t1+α,]]>其中α为与发射信号的信噪比相关的常数；步骤c22所述Rm(tm)中不包含在Rm-1(tm-1)中的部分项为一个标量β1(tm)和一个向量vm-1(tm)，其中，β1(tm)=h:tmH·h:tm+α,vm-1(tm)=h:t1H·h:tmh:t1H·h:tm···h:tm-1H·htm;]]>步骤c22所述递推得到的m个发射信号的估计误差协方差矩阵为在m-1个发射信号的估计误差协方差矩阵的Sherman-Morrison结果基础上，增加一行和一列得到的矩阵。
8.根据权利要求7所述的方法，其特征在于，步骤c22所述在m-1个发射信号的估计误差协方差矩阵的Sherman-Morrison结果基础上，增加一行和一列得到的矩阵为m-1个发射信号的估计误差协方差矩阵Qm-1(tm-1)的Sherman-Morrison结果为(Tm-1(tm))-1=Qm-1(tm-1)+Qm-1(tm-1)vm-1(tm)(vm-1(tm))HQm-1(tm-1)β1(tm)-(vm-1(tm))HQm-1(tm-1)vm-1(tm),]]>在(Tm-1(tm))-1基础上，增加由矩阵(Tm-1(tm))-1、向量vm-1(tm)以及标量βm-1(tm)组成的一列和一行，得到m个发射信号的估计误差协方差矩阵Qm(tm)，其中，一列和一行相交的项是1β1(tm)+(vm-1(tm))H(Tm-1(tm))-1vm-1(tm)(β1(tm))2,]]>一列的其它项为 一行的其它项为 其中，Qm-1(tm-1)为上一次递推的结果或步骤c21中得到的Q1(t1)。
9.根据权利要求2所述的方法，其特征在于，所述步骤d包括d1.在待检测发射信号中选择当前被检测的一个发射信号，利用步骤c得到的所有发射信号的估计误差协方差矩阵、信道矩阵H以及接收信号得到对所述当前被检测的一个发射信号的估计值；d2.利用步骤d1得到的当前被检测的一个发射信号的估计值计算对检测后续待检测发射信号的干扰值，消除所述当前被检测的一个发射信号对检测后续待检测发射信号的干扰；d3.重复步骤d1、d2，直到检测到所有待检测发射信号。
10.根据权利要求9所述的方法，其特征在于，所述步骤d1之前进一步包括利用信道矩阵H对接收信号进行预匹配滤波变换；计算信道矩阵H的互相关信道矩阵Φ，Φ＝HH·H；步骤d1所述得到对当前被检测的一个发射信号的估计值的步骤包括利用待检测发射信号的估计误差协方差矩阵和所述接收信号的预匹配滤波结果得到所述当前被检测的一个发射信号的估计值；所述步骤d2包括利用所述当前被检测的一个发射信号的估计值和信道矩阵H的互相关信道矩阵Φ，计算已检测的发射信号对检测后续发射信号的干扰值，并从所述接收信号的预匹配滤波结果中消除已检测的发射信号的干扰，得到修正的接收信号的预匹配滤波结果，作为下一次检测信号时的接收信号的预匹配滤波结果。
11.根据权利要求10所述的方法，其特征在于，所述计算信道矩阵H的互相关信道矩阵Φ的步骤包括利用信道矩阵H计算发射信号的估计误差协方差矩阵的逆矩阵R，利用Φ＝HH·H和R＝HH·H+αIM×M的关系，得到Φ。
12.根据权利要求10所述的方法，其特征在于，所述利用信道矩阵H对接收信号进行预匹配滤波变换的步骤包括将信道矩阵H的共轭转置矩阵作为接收信号的预匹配滤波器，对接收信号向量进行预匹配滤波得到接收信号的预匹配滤波结果；步骤d1所述在待检测发射信号中选择当前被检测的一个发射信号为待检测发射信号的估计误差协方差矩阵中对角元素最小的行对应的发射信号为当前被检测的一个发射信号；步骤d1所述利用待检测发射信号的估计误差协方差矩阵和所述接收信号的预匹配滤波结果得到当前被检测的一个发射信号的估计值的步骤包括待检测发射信号的估计误差协方差矩阵中对角元素最小的行，与接收信号的预匹配滤波结果相乘，得到所述当前被检测的一个发射信号的估计值；所述步骤d2包括根据所述当前被检测的发射信号的估计值和所述信道矩阵的互相关信道矩阵Φ中与当前被检测的发射信号对应的元素组成的向量的乘积，得到已检测的发射信号对检测后续发射信号的干扰值，然后从接收信号的预匹配滤波结果中删除已检测发射信号对应的一项，再从所述删除一项后的接收信号的预匹配滤波结果中，消除所述干扰得到修正的接收信号的预匹配滤波结果。
13.根据权利要求12所述的方法，其特征在于，所述步骤d2之后，在检测下一个发射信号之前进一步包括判断是否在步骤c所述递推所有发射信号的估计误差协方差矩阵的过程中，已得到下一次检测信号时的待检测发射信号对应的估计误差协方差矩阵，如果是，则直接执行检测下一个发射信号的步骤；否则，利用本次检测信号时的待检测发射信号对应的估计误差协方差矩阵，递推得到下一次检测信号时的待检测发射信号对应的估计误差协方差矩阵。
14.根据权利要求13所述的方法，其特征在于，所述利用本次检测信号时的待检测发射信号对应的估计误差协方差矩阵递推得到下一次检测信号时的待检测发射信号对应的估计误差协方差矩阵为从本次检测信号时的待检测发射信号对应的估计误差协方差矩阵中，删除对角元素最小的行和列得到的子矩阵的Sherman-Morrison结果，作为下一次检测信号时的待检测发射信号对应的估计误差协方差矩阵。
全文摘要
本发明公开了一种在多天线数字无线通信系统中检测信号的方法，该方法包括接收端的至少两个接收天线接收发射信号，获得至少两个接收信号；接收端根据接收信号进行信道估计，得到由发射天线和接收天线之间的信道系数组成的信道矩阵H；利用信道矩阵H计算出所有发射信号中的部分发射信号的估计误差协方差矩阵，然后利用信道矩阵H和计算出的部分发射信号的估计误差协方差矩阵，递推求得所有发射信号的估计误差协方差矩阵；利用所得到的所有发射信号的估计误差协方差矩阵，检测所有发射信号。根据本发明公开的方法，减少了检测信号的计算复杂度。
文档编号H04L25/02GK1905432SQ20061010404
公开日2007年1月31日 申请日期2006年7月31日 优先权日2006年7月31日
发明者朱胡飞 申请人:华为技术有限公司