基于分组序列码本集的大规模终端多址接入方法

1.本发明涉及通信技术领域，具体涉及一种基于分组序列码本集的大规模终端多址接入方法。


背景技术：

2.物联网设备数量的大幅增加对无线系统的设计提出了新的挑战，其中，大规模机器类型通信(massive machine type communication，mmtc)是第五代移动通信系统的基本组成部分之一。基于压缩感知的多用户检测(compressive sensing multiuser detection，csmud)是应对mmtc大规模连接需求的候选者，它促进了免授权非正交码分多址的使用，以容纳大量的物联网设备。csmud对活动用户和数据进行联合检测的能力促进了直接随机访问的可靠检测，以最小的性能损失为代价来提高频谱效率并简化传感器节点处的信息处理。
3.csmud利用mmtc零星传输的特性，从接收信号中引入的稀疏性，进而检测传感器是否处于活动状态(活动检测)以及传感器传输的数据(数据检测)。csmud为用户分配非正交扩展序列作为其签名，其性能往往取决于扩展序列之间的相关性。而随着未来用户数量的增加，扩展序列之间的相关性将越来越高，从而为mmtc网络中的多用户检测问题带来了挑战。有人试过通过增加扩展序列的长度来降低序列之间相关性，但效果并不明显，且会浪费更多的无线资源；如何基于多种传输需求建立多用户检测方案，以提高系统检测的性能，是一个需要解决的技术问题。


技术实现要素：

4.为解决上述问题，本发明提出了一种基于分组序列码本集的大规模终端多址接入方法，该方法将所有用户分组形成多个子系统，在每个子系统中分配特定的序列码本集，在提高系统容量的基础上进一步降低扩展序列之间的相关性，从而提高活动用户与数据联合检测的准确性；同时将csmud与多载波系统结合，利用多用户信号的稀疏结构，通过增加总体负载进一步提高频谱效率。请注意，本发明中终端、传感器节点和用户可互换使用，接收端与基站可互换使用。
5.本发明所述的基于分组序列码本集的大规模终端多址接入方法，包括以下步骤：
6.步骤1、对于拥有大规模终端的mmtc系统，将全部用户进行分组形成多个用户子系统，定义每个子系统的用户集合；
7.步骤2、为每个子系统设计唯一的基本扩展序列，根据该序列生成每个子系统的序列码本集，为所述子系统中的每个用户分配其序列码本集中的一个唯一的扩展序列；
8.步骤3、将每个子系统的用户数据经过调制、扩展序列扩频后相加，形成该子系统的和信号；每个子系统的和信号通过正交频分复用(ofdm)技术映射到多个子载波上形成时频数据块；不同子系统的时频数据块重叠复用部分子载波，通过上行信道传输至基站；
9.步骤4、基站接收多个子系统的上行传输信号并恢复出各个子系统的和信号；对于
每个子系统的和信号，采用基于序列码本块的组正交匹配追踪(codebooksequence block based group orthogonal matching pursuit，cb-gomp)算法对其中的活动用户和用户数据进行联合检测，完成大规模用户的随机接入与数据传输。
10.进一步的，步骤1中，所述mmtc系统中共有k个用户，k个用户位于同一个基站的覆盖范围内，每个用户向无线网络发起接入请求和数据传输；将全部k个用户划分为nb组，形成nb个子系统；子系统标号表示为b，它满足：
11.b∈b＝{1,...,nb}
ꢀꢀꢀ
(1)
12.其中，b表示子系统标号集合；
13.每个子系统包含kb个用户，令子系统b中的用户标号为kb，它满足：
14.kb∈{1,...,kb}
ꢀꢀꢀ
(2)
15.其中，表示取上整；
16.用户的数据传输处于偶发状态，整个系统用户总数为k＝nbkb。
17.进一步的，步骤1中，假设每个用户的活动概率为pa，pa＝1，即用户以概率pa进行传输，以1-pa的方式保持静默；因此，在任一给定的时刻，k个用户中只有小部分用户处于活动状态；当一个用户处于活动状态时，它在传输时隙中连续发送短包数据。
18.进一步的，步骤2中，所述每个子系统的基本扩展序列设计方法如下：
19.用sb表示子系统b的基本扩展序列，从单位圆上随机抽取nb个样本，其中第i个样本表示为
[0020][0021]
其中，u服从区间[0,1]上的均匀分布，表示样本对应的相位；的实部和虚部数值分别转换为长度为的二进制码字，最后将它们拼接得到长度为nc的二进制序列，即为子系统b的基本扩展序列sb；
[0022]
nb个子系统的基本扩展序列的组合矩阵表示为：
[0023][0024]
进一步的，步骤2中，所述子系统的序列码本集由一组扩展序列组成，每个扩展序列均由该子系统的基本扩展序列根据特定移位模式进行循环移位后生成，每个子系统所属序列码本集的生成方如下；
[0025]
步骤2-1、生成子系统的特定移位模式；
[0026]
令子系统b，b∈b＝{1,...,nb}的特定移位模式pb为：
[0027][0028]
其中，p
b,1
＝0，kb表示子系统b中的用户数目；d
b,k
表示子系统b中用户k在pb中的索引，即向量pb的第d
b,k
个元素；表示pb的第d
b,k
个元素的数值；
[0029]
子系统b的移位模式pb的第m个索引值p
b,m
,m≠1由下式计算得到：
[0030][0031]
其中，nc表示基本扩展序列sb的长度；表示基本扩展序列sb循环右移j位得到的新扩展序列；表示基本扩展序列sb循环右移p
b,i
位得到的新扩展序列；表示和的相关运算值。
[0032]
步骤2-2、生成子系统的序列码本集；
[0033]
根据子系统b的特定移位模式pb中的kb个元素对基本扩展序列sb分别进行kb次循环右移，得到子系统b中全部kb个用户的对应扩展序列：
[0034][0035]
其中，sb表示子系统b对应的基本扩展序列，pb表示子系统b对应的移位模式，它是长度为kb的向量；d
b,k
表示子系统b中用户k在pb中的索引，即向量pb的第d
b,k
个元素；表示pb的第d
b,k
个元素的数值；是序列sb循环右移位后生成的新序列；
[0036]
这些扩展序列的集合即为子系统b的序列码本集；
[0037]
子系统b中所有用户的扩展序列构成集合即为系统b的序列码本集，定义序列码本矩阵
[0038]
进一步的，步骤3中，所述用户数据的调制、扩频、复用过程如下：
[0039]
步骤3-1、单个子系统的用户数据经调制与扩频后形成和信号；
[0040]
每个子系统中每个用户对其待发送数据进行数字调制后形成l个调制符号，若某一用户无数据发送，则认为其传输l个0；不失一般性，令子系统b，b∈b＝{1,...,nb}中所有用户的发送符号矩阵为
[0041][0042]
其中，表示大小为kb×
l的复数集合；为db的第kb行向量，它表示子系统b中第kb个用户发送的l个连续的时间符号；如果某一用户ku处于静默状态，则db中的第ku行填充l个0；db的每个列向量包含来自kb个用户的符号；
[0043]
用户利用自己的特定扩展序列对已调符号进行扩频，每个已调符号形成长度为nc的扩频序列；每个子系统中所有用户的数据经过调制、扩频后得到的序列相加，得到该子系统的和信号；令子系统b，b∈b＝{1,...,nb}的和信号矩阵为即
[0044]
xb＝dbsbꢀꢀꢀ
(9)
[0045]
其中，为子系统b的码本矩阵；表示大小为nc×
l的复数域空间；
[0046]
每个子系统的和信号通过ofdm技术映射到ns个子载波中形成一个时频资源块；令子系统b，b∈b＝{1,...,nb}的和信号经ofdm技术变换后的频域符号为其中表示大小为ns×
l的复数域空间；
[0047]
步骤3-2、多个子系统的和信号通过载波部分重叠方式复用传输；
[0048]
系统中共有nb个子系统，因此共有nb个时频资源块，每个时频资源块包含ns个连续子载波；令相邻子系统复用no个子载波，0≤no≤ns，nb个子系统的和信号按照载波复用图样映射叠加为系统的上行传输信号表示为：
[0049][0050]
其中，表示大小为n
max
×
l的复数集合；通过重叠复用子载波，系统占据的子载波总数n
max
计算为：
[0051][0052]
定义系统承载的用户总数与子载波总数的比值为系统载波负载，通过子系统之间的子载波复用，系统载波负载计算为：
[0053][0054]
其中，k＝nbkb；
[0055]
该信号经天线馈送至上行信道中传输给基站。
[0056]
进一步的，步骤4中，基站接收多个子系统的上行传输信号并恢复出各个子系统的和信号，所述上行信号接收处理步骤如下：
[0057]
步骤4-1、基站接收多个子系统的上行传输信号；
[0058]
考虑平坦衰落信道，基站接收的频域信号表示为：
[0059]
y＝hx+n
ꢀꢀꢀ
(12)
[0060]
其中，表示大小为n
max
×
l的复数域空间；矩阵表示信道系数矩阵，它是对角矩阵，即
[0061][0062]
其中，表示大小为n
max
×nmax
的复数域空间；hs表示子载波s的平坦衰落系数；n表示信道加性白噪声矩阵，它服从均值为0方差为的复高斯分布，即
[0063]
通过ofdm技术解调并消除信道干扰后，为了简化系统模型，使ns＝nc，将nb个子系统在频域上的和信号转换为时域和信号
[0064]
步骤4-2、基站从上行传输信号中恢复每个子系统的和发送信号；
[0065]
基站接收到多个子系统中活跃用户的上行信号后，通过ofdm技术反变换恢复出单个子系统对应的接收信；令基站接收恢复出的子系统b的和信号为yb表示为：
[0066]
yb＝hbxb+nbꢀꢀꢀ
(14)
[0067]
其中，表示大小为nc×
l的复数域空间；矩阵表示子系统b，b∈b＝{1,...,nb}和信号矩阵；表示信道系数矩阵，它是一个对角矩阵，表示为：
[0068]
hb＝diag{hb}
ꢀꢀꢀ
(15)
[0069]
其中，向量表示子系统b中每个终端用户的平坦衰落系数的集合，表示子系统b中终端用户k的平坦衰落系数；nb表示均值为0方差为的加性高斯白噪声矩阵，即
[0070]
步骤4-3、基站采用cb-gomp算法从每个子系统和信号中联合检测活跃用户及用户数据。
[0071]
进一步的，利用cb-gomp算法迭代检测每个子系统中的活跃用户并恢复其发送信号，以子系统b的序列码本接收信号y为输入，子系统中用户的估计矩阵作为输出；所述cb-gomp算法的具体步骤如下：
[0072]
步骤4-3-1、初始化算法参数；
[0073]
令子系统b的支持集的估计，即子系统b中活跃用户序号集合为γb，其初始值为空集，即令检测子系统码本索引的迭代次数为g，检测子系统b中活跃用户数据的
迭代次数q的初值为0；基站接收到全部子系统的信号y为：
[0074][0075][0076]
其中，y
l
表示当前帧的第l个符号；nc为扩展序列的长度，l为一帧内发送的符号数目；将y
l
分割为nb个子信号，如式(17)所示；表示y
l
分组得到的第g个子信号；残差r表示接收信号与算法恢复信号的差，残差r
l
表示接收信号y
l
与算法恢复信号的差，1≤l≤l；表示迭代次数为q-1时的残差；表示残差r
l
的初值，其中初始化为g初始化为1；执行步骤4-3-2；
[0077]
步骤4-3-2、检测活动用户的子系统码本索引；
[0078]
对接收到的一帧l个符号，即利用式(18)分别求出序列码本集sb各列即与l组初始残差的内积值之和，比较nb组求和结果，得到求和结果最大的一组所对应的子系统的码本索引ig后，执行步骤4-3-3；
[0079][0080]
步骤4-3-3、检测活动用户索引；利用ig对应的序列码本集与残差进行组正交匹配追踪gomp算法求解；
[0081]
活动用户的检测基于从接收信号中获取与码本扩展序列的最大相关值对应的用户索引，通过公式(19)迭代求出序列码本集各列与l组残差值的内积值之和，比较kb组求和结果，得到求和结果数值最大的一组所对应的活动用户的索引，利用公式(20)更新子系统ig中活动用户的支持集；执行步骤4-3-4；
[0082][0083][0084]
步骤4-3-4、检测子系统活动用户数据；
[0085]
利用公式(21)，使用最小二乘估计来估计当前迭代时子系统ig中用户数据的集合其中表示伪逆计算；遍历循环l组残差后，执行步骤4-3-5；
[0086][0087]
步骤4-3-5、更新残差；利用当前检测用户的扩展序列，根据公式(22)，从中减
去估计数据的贡献来更新残差根据公式(23)，残差能量小于阈值时停止更新，得到当前检测子系统ig中用户数据的估计矩阵并保存到系统用户数据的估计矩阵中；遍历更新l组残差后，当前子系统的活动用户数据检测完毕后，执行步骤4-3-6；
[0088][0089]
||r||＜γ
ꢀꢀꢀ
(23)
[0090]
步骤4-3-6、对下一个子系统的活动用户和数据进行迭代检测；将g更新为g+1,则分组信号更新为利用已检测的子系统的支持集和序列码本集，根据公式(24)，从而消除已检测的子系统的活动用户数据对待测子系统的活动用户数据的干扰；消除干扰后的信号为其中可以表示为：
[0091][0092]
迭代更新初始残差返回步骤4-3-2，对待测子系统的活动用户以及数据进行检测，更新子系统i
g+1
用户数据的估计矩阵已检测的子系统不再被重复检测；直到迭代次数达到g＞nb时，活动用户和数据的联合检测结束，得到系统用户数据的估计矩阵
[0093]
本发明的有益效果为：(1)本发明所述的方法改善了用户之间扩展序列的相关性问题，基于分组序列码本集方法，基本扩展序列形成了特定的循环移位版本，从而减小信令开销，在提高系统容量的基础上进一步降低扩展序列之间的相关性，从而进一步提高活动用户与数据联合检测的准确性；
[0094]
(2)本发明所述的方法将csmud与多载波系统结合，利用多用户信号的稀疏结构，通过增加总体负载进一步提高了频谱效率，并通过cb-gomp算法提高csmud的性能，本发明所述的方法能适应于mmtc场景下的大规模连接，相比于传统的方法，具有更高的实用价值。
附图说明
[0095]
图1为本发明实施例中所述的基于分组序列码本集的大规模终端多址接入方法的流程图；
[0096]
图2为本发明实施例中系统场景图；
[0097]
图3为本发明实施例中扩展序列的分配图；
[0098]
图4为本发明实施例中传感器节点的处理过程图；
[0099]
图5为本发明实施例中子系统csmud的资源块图；
[0100]
图6为本发明实施例中cb-gomp算法流程图。。
具体实施方式
[0101]
为了使本发明的内容更容易被清楚地理解，下面根据具体实施例并结合附图，对本发明作进一步详细的说明。
[0102]
本发明所述的基于分组序列码本集的大规模终端多址接入方法的流程图如图1所示，所述方法包括以下步骤：
[0103]
步骤1中，考虑一个典型的上行mmtc系统场景，如图2所示；所述mmtc系统中共有k个用户，k个用户位于同一个基站的覆盖范围内，每个用户可以向无线网络发起接入请求和数据传输；将全部k个用户划分为nb组，形成nb个子系统；子系统标号表示为b，它满足：
[0104]
b∈b＝{1,...,nb}
ꢀꢀꢀ
(1)
[0105]
其中，b表示子系统标号集合；
[0106]
每个子系统包含kb个用户，令子系统b中的用户标号为kb，它满足：
[0107]
kb∈{1,...,kb}
ꢀꢀꢀ
(2)
[0108]
其中，表示取上整；
[0109]
用户的数据传输处于偶发状态，整个系统用户总数为k＝nbkb。
[0110]
假设每个用户的活动概率为pa，pa＝1，即用户以概率pa进行传输，以1-pa的方式保持静默；因此，在任一给定的时刻，k个用户中只有小部分用户处于活动状态；当一个用户处于活动状态时，它在传输时隙中连续发送短包数据。
[0111]
步骤2中，所述每个子系统的基本扩展序列设计方法如下：
[0112]
用sb表示子系统b的基本扩展序列，从单位圆上随机抽取nb个样本，其中第i个样本表示为
[0113][0114]
其中，u服从区间[0,1]上的均匀分布，表示样本对应的相位；的实部和虚部数值分别转换为长度为的二进制码字，最后将它们拼接得到长度为nc的二进制序列，即为子系统b的基本扩展序列sb；
[0115]
nb个子系统的基本扩展序列的组合矩阵表示为：
[0116][0117]
步骤2中，所述子系统的序列码本集由一组扩展序列组成，每个扩展序列均由该子系统的基本扩展序列根据特定移位模式进行循环移位后生成，每个子系统所属序列码本集的生成方如下；
[0118]
步骤2-1、生成子系统的特定移位模式；
[0119]
令子系统b，b∈b＝{1,...,nb}的特定移位模式pb为：
[0120][0121]
其中，p
b,1
＝0，kb表示子系统b中的用户数目；d
b,k
表示子系统b中用户k在pb中的索
引，即向量pb的第d
b,k
个元素；表示pb的第d
b,k
个元素的数值；
[0122]
子系统b的移位模式pb的第m个索引值p
b,m
,m≠1由下式计算得到：
[0123][0124]
其中，nc表示基本扩展序列sb的长度；表示基本扩展序列sb循环右移j位得到的新扩展序列；表示基本扩展序列sb循环右移p
b,i
位得到的新扩展序列；表示和的相关运算值；
[0125]
步骤2-2、生成子系统的序列码本集；
[0126]
根据子系统b的特定移位模式pb中的kb个元素对基本扩展序列sb分别进行kb次循环右移，得到子系统b中全部kb个用户的对应扩展序列：
[0127][0128]
其中，sb表示子系统b对应的基本扩展序列，pb表示子系统b对应的移位模式，它是长度为kb的向量；d
b,k
表示子系统b中用户k在pb中的索引，即向量pb的第d
b,k
个元素；表示pb的第d
b,k
个元素的数值；是序列sb循环右移位后生成的新序列；
[0129]
这些扩展序列的集合即为子系统b的序列码本集；
[0130]
如图3所示，子系统b中所有用户的扩展序列构成集合即为系统b的序列码本集，定义序列码本矩阵
[0131]
步骤3中，所述用户数据的调制、扩频、复用过程如图4所示，具体步骤如下：
[0132]
步骤3-1、单个子系统的用户数据经调制与扩频后形成和信号；
[0133]
每个子系统中每个用户对其待发送数据进行数字调制后形成l个调制符号，若某一用户无数据发送，则认为其传输l个0；不失一般性，令子系统b，b∈b＝{1,...,nb}中所有用户的发送符号矩阵为
[0134][0135]
其中，表示大小为kb×
l的复数集合；为db的第kb行向量，它表示子系统b中第kb个用户发送的l个连续的时间符号；如果某一用户ku处于静默状态，则db中的第ku行填充l个0；db的每个列向量包含来自kb个用户的符号；
[0136]
用户利用自己的特定扩展序列对已调符号进行扩频，每个已调符号形成长度为nc的扩频序列；每个子系统中所有用户的数据经过调制、扩频后得到的序列相加，得到该子系统的信号；令子系统b，b∈b＝{1,...,nb}的信号矩阵为即
[0137]
xb＝dbsbꢀꢀꢀ
(9)
[0138]
其中，为子系统b的码本矩阵；表示大小为nc×
l的复数域空间；
[0139]
每个子系统的和信号通过ofdm技术映射到ns个子载波中形成一个时频资源块，单个子系统csmud的时频资源块如图5所示；令子系统b，b∈b＝{1,...,nb}的和信号经ofdm技术变换后的频域符号为其中表示大小为ns×
l的复数域空间；
[0140]
步骤3-2、多个子系统的和信号通过载波部分重叠方式复用传输；
[0141]
系统中共有nb个子系统，因此共有nb个时频资源块，每个时频资源块包含ns个连续子载波；令相邻子系统复用no个子载波，0≤no≤ns，nb个子系统的和信号按照载波复用图样映射叠加为系统的上行传输信号表示为：
[0142][0143]
其中，表示大小为n
max
×
l的复数集合；通过重叠复用子载波，系统占据的子载波总数n
max
计算为：
[0144][0145]
定义系统承载的用户总数与子载波总数的比值为系统载波负载，通过子系统之间的子载波复用，系统载波负载计算为：
[0146][0147]
其中，k＝nbkb；
[0148]
该信号经天线馈送至上行信道中传输给基站。
[0149]
步骤4中，基站接收多个子系统的上行传输信号并恢复出各个子系统的和信号，所述上行信号接收处理步骤如下：
[0150]
步骤4-1、基站接收多个子系统的上行传输信号；
[0151]
考虑平坦衰落信道，基站接收的频域信号表示为：
[0152]
y＝hx+n
ꢀꢀꢀ
(12)
[0153]
其中，表示大小为n
max
×
l的复数域空间；矩阵表示信道系数矩阵，它是对角矩阵，即
[0154][0155]
其中，表示大小为n
max
×nmax
的复数域空间；hs表示子载波s的平坦衰落系数；n表示信道加性白噪声矩阵，它服从均值为0方差为的复高斯分布，即
[0156]
通过ofdm技术解调并消除信道干扰后，为了简化系统模型，使ns＝nc，将nb个子系统在频域上的和信号转换为时域和信号
[0157]
步骤4-2、基站从上行传输信号中恢复每个子系统的和发送信号；
[0158]
基站接收到多个子系统中活跃用户的上行信号后，通过ofdm技术反变换恢复出单个子系统对应的接收信；令基站接收恢复出的子系统b的信号为yb表示为：
[0159]
yb＝hbxb+nbꢀꢀꢀ
(14)
[0160]
其中，表示大小为nc×
l的复数域空间；矩阵表示子系统b，b∈b＝{1,...,nb}和信号矩阵；表示信道系数矩阵，它是一个对角矩阵，表示为：
[0161]
hb＝diag{hb}
ꢀꢀꢀ
(15)
[0162]
其中，向量表示子系统b中每个终端用户的平坦衰落系数的集合，表示子系统b中终端用户k的平坦衰落系数；nb表示均值为0方差为的加性高斯白噪声矩阵，即
[0163]
步骤4-3、基站采用cb-gomp算法从每个子系统和信号中联合检测活跃用户及用户数据。
[0164]
如图6所示，利用cb-gomp算法迭代检测每个子系统中的活跃用户并恢复其发送信号，以子系统b的序列码本接收信号y为输入，子系统中用户的估计矩阵作为输出；所述cb-gomp算法的具体步骤如下：
[0165]
步骤4-3-1、初始化算法参数；
[0166]
令子系统b的支持集的估计，即子系统b中活跃用户序号集合为γb，其初始值为
空集，即令检测子系统码本索引的迭代次数为g，检测子系统b中活跃用户数据的迭代次数q的初值为0；基站接收到全部子系统的信号y为：
[0167][0168][0169]
其中，y
l
表示当前帧的第l个符号；nc为扩展序列的长度，l为一帧内发送的符号数目；将y
l
分割为nb个子信号，如式(17)所示；表示y
l
分组得到的第g个子信号；残差r表示接收信号与算法恢复信号的差，残差r
l
表示接收信号y
l
与算法恢复信号的差，1≤l≤l；表示迭代次数为q-1时的残差；表示残差r
l
的初值，其中初始化为g初始化为1；执行步骤4-3-2；
[0170]
步骤4-3-2、检测活动用户的子系统码本索引；
[0171]
对接收到的一帧l个符号(即利用式(18)分别求出序列码本集sb各列(即)与l组初始残差的内积值之和，比较nb组求和结果，得到求和结果最大的一组所对应的子系统的码本索引ig后，执行步骤4-3-3；
[0172][0173]
步骤4-3-3、检测活动用户索引；利用ig对应的序列码本集与残差进行组正交匹配追踪gomp算法求解；
[0174]
活动用户的检测基于从接收信号中获取与码本扩展序列的最大相关值对应的用户索引，通过公式(19)迭代求出序列码本集各列与l组残差值的内积值之和，比较kb组求和结果，得到求和结果数值最大的一组所对应的活动用户的索引，利用公式(20)更新子系统ig中活动用户的支持集；执行步骤4-3-4；
[0175][0176][0177]
步骤4-3-4、检测子系统活动用户数据；
[0178]
利用公式(21)，使用最小二乘估计来估计当前迭代时子系统ig中用户数据的集合其中表示伪逆计算；遍历循环l组残差后，执行步骤4-3-5；
[0179]
[0180]
步骤4-3-5、更新残差；利用当前检测用户的扩展序列，根据公式(22)，从中减去估计数据的贡献来更新残差根据公式(23)，残差能量小于阈值时停止更新，得到当前检测子系统ig中用户数据的估计矩阵并保存到系统用户数据的估计矩阵中；遍历更新l组残差后，当前子系统的活动用户数据检测完毕后，执行步骤4-3-6；
[0181][0182]
||r||＜γ
ꢀꢀꢀ
(23)
[0183]
步骤4-3-6、对下一个子系统的活动用户和数据进行迭代检测；将g更新为g+1,则分组信号更新为利用已检测的子系统的支持集和序列码本集，根据公式(24)，从而消除已检测的子系统的活动用户数据对待测子系统的活动用户数据的干扰；消除干扰后的信号为其中可以表示为：
[0184][0185]
迭代更新初始残差返回步骤4-3-2，对待测子系统的活动用户以及数据进行检测，更新子系统i
g+1
用户数据的估计矩阵已检测的子系统不再被重复检测；直到迭代次数达到g＞nb时，活动用户和数据的联合检测结束，得到系统用户数据的估计矩阵
[0186]
以上所述仅为本发明的较佳实施方式，本发明的保护范围并不以上述实施方式为限，但凡本领域普通技术人员根据本发明所揭示内容所作的等效修饰或变化，皆应纳入权利要求书中记载的保护范围内。