多用户mimo-ofdm上行链路资源分配方法

专利名称：：多用户mimo-ofdm上行链路资源分配方法
技术领域：
：本发明涉及通信技术，具体说就是一种多用户MIMO-OF匿上行链路资源分配方法。(二)
背景技术：
：多用户多输入多输出_正交频分复用(MMO-OFDM)系统上行链路资源分配算法，能够有效利用空域、频域、时域资源，具有频谱利用率高、抗多径时延等优点，已成为下一代宽带移动通信技术的热点。由于共道干扰(CCI)的存在，多用户MIMO-OFDM系统上行链路资源分配算法，多采用子载波独占方式。子载波独占方式资源分配，每个子载波仅分配给一个用户，任意子载波不存在共道干扰，属频分复用方式。但是，由于子载波独占方式并没有充分利用空间资源，频谱利用率低。多用户MIMO-OFDM系统上行链路子载波时分复用方式(TDMA)资源分配，各用户分时独占所有子载波。同样，由于没有充分利用空间资源，频谱利用率低。多用户MIM0-0F匿系统上行链路子载波共享方式资源分配，充分利用空间资源，各用户共用时域、频域资源，属空分复用方式，频谱利用率高。但由于共道干扰的存在，系统性能受用户空间相关性的影响严重。典型共道干扰抑制方法是利用波束形成技术的基于多用户信道矩阵求逆的迫零算法，但迫零算法高阶矩阵奇异值分解(SVD)[1]的计算复杂度高，限制了迫零共享方式的发展。文献[2]提出在基站端利用多用户检测的分组共享方法，将空间相关性大(波达角(DOA)<20°)的用户分在一组，每组内用户以独占的方式分配全部子载波，每组间的用户由于空间相关性很小，基站端可以利用匹配滤波器消除共道干扰，从而实现子载波的共享，本质上属于直接共享方式。然而，该方法在用户分组时，同一组用户DOA的夹角要求很大，因此，适用性受限。用户空间相关性对系统性能的影响，以及多用户共道干扰消除的计算复杂度，制约了多用户MIMO-OFmi系统上行链路子载波共享方式的发展。(三)
发明内容本发明的目的在于提供一种适用于多用户MIM0-0F匿系统上行链路子载波、子空间、比特、功率资源分配优化过程的多用户MIMO-OFDM上行链路资源分配方法。本发明的目的是这样实现的所述的一种多用户MIMO-OFDM上行链路资源分配方法，它由以下步骤实现步骤一用户按空间相关性分组；空间相关性高的用户分在一组，DOA夹角小于5度的用户分在同一组；步骤二每组内用户按子载波独占方式分配子载波；依据比例公平原则，对组内用户以独占方式分配子载波，任意子载波每组有且仅有一个用户占用。定义子载波分配因子wmk，为1表示用户k占用子载波m，为0表示用户k不占用子载波m;步骤三各组间共享同一子载波的用户利用波束形成技术，基于零空间交集的连续迫零方法，实现共道干扰抑制的同时完成子空间的分配。对每个子载波，将共享该子载波的用户编号为1K，利用基于零空间交集的连续迫零方法，按如下步骤，以用户编号的正序计算各用户的迫零矩阵Z、发射加权矩阵W、接收合并矩阵R;HZ2，Z3H:AJ,，令R丄H，W!二V。(。)]a2v2h4r2=[U2(1)]hz2，w2v9U3(°)]八,/，令尺3=[U3")]HZ3，W:(1)—对用户k:Zk对用户K:ZK%—pZkHk=[UkAJkH，令Rk(1)(0)—Uk(。)]AkV/，令RkHZk，Wk=VkHzK，wK=vKHZK—i，ZKHK=[UK(1):步骤四、依据等效信道增益，自适应比特功率分配；对每个用户，以单用户方式对所分配的所有子载波进行比特、功率分配，典型算法步骤五、基站端基于连续干扰消除(SIC)方法进行多用户检测；基站端对接收的多用户的信号和进行分离，计算出每个共享用户的信号。对每个对用户K:1对用户k:a:A(y層￡好严,》=[^(1)产4(/^^^+2://》+^"=4^+^KA)-、(r—Z好》对用户1:A2好y)=[171(1)产(好iW义i+")=4zi+尺i"步骤六、各用户依据最大似然法则进行MQAM解码。本发明的核心是针对多用户MIMO-OF匿系统上行链路，在按用户空间相关性对用户分组的基础上，提出一种基于零空间交集的连续迫零算法结合基于连续干扰消除的多用户检测算法的低复杂度的子载波共享方式共道干扰抑制方法，进而提出多用户MIMO-OF匿系统上行链路子载波共享方式的资源分配方法，在提高多用户MIMO-OFDM系统频谱利用率等性能的同时，降低了资源分配算法的计算复杂度。本发明适用于多用户MIM0-0F匿系统上行链路子载波、子空间、比特、功率等资源分配的优化过程。图1为不同资源分配方式系统性能与用户空间相关性曲线图；图2为直接共享资源分配方式系统性能与用户空间相关性曲线图；图3为本发明的多用户MIM0-0FDM上行链路系统结构图；图4为本发明的基于零空间交集的连续迫零算法结合基于连续干扰消除的多用户检测的原理框图；图5为本发明方法与其他方法的用户子载波平均比特数与Eb/No图；图6为本发明的SIC-ZF方法各用户子载波平均比特数与Eb/No图；图7为本发明的SIC-GZF方法各用户子载波平均比特数与Eb/No;图8为本发明方法与其他方法的BER与Eb/No。具体实施例方式下面结合附图举例对本发明作进一步说明。实施例1:本发明一种多用户MIMO-OFDM上行链路资源分配方法，它由以下步骤实现步骤一用户按空间相关性分组；空间相关性高的用户分在一组，D0A夹角小于5度的用户分在同一组；步骤二每组内用户按子载波独占方式分配子载波；依据比例公平原则，对组内用户以独占方式分配子载波，任意子载波每组有且仅有一个用户占用。定义子载波分配因子w二为1表示用户k占用子载波m，为0表示用户k不占用子载波m;步骤三各组间共享同一子载波的用户利用波束形成技术，基于零空间交集的连续迫零方法，实现共道干扰抑制的同时完成子空间的分配。对每个子载波，将共享该子载波的用户编号为1K，利用基于零空间交集的连续迫零方法，按如下步骤，以用户编号的正序计算各用户的迫零矩阵Z、发射加权矩阵W、接收合并矩阵R;对用户1:Z丄=l，H丄=[U/1)，U/0)]AJ,，令R丄=[U/"]h，W丄=V工；对用户2:Z2=[U/H"Z2H2=[U2(1)，U2(°)]A2V/，令R2=[U2(1)]HZ2，W2=V2;对用户3:Z3=[U2(°)]HZ2，Z3H3=[U3(1)，U3(°)]A3V/，令R3=[U3(1)]HZ3，W3=V3;......对用户k:Zk=[Uk—/。)]%—"ZkHk=[Uk("，Uk(。)]AJkH，令Rk=[Uk(1)]HZk，Wk=Vk;......对用户K:ZK=[UK—/。)]%—nZKHK=[UK(1)，UK(°)]AkV/，令RK=[UK(1)]HZK，WK=VK;步骤四、依据等效信道增益，自适应比特功率分配；对每个用户，以单用户方式对所分配的所有子载波进行比特、功率分配，典型算法6为贪婪算法；步骤五、基站端基于连续干扰消除方法进行多用户检测；基站端对接收的多用户的信号和进行分离，计算出每个共享用户的信号。对每个子载波，按如下步骤，以用户编号的逆序计算各共享用户的接收信号；对用户K:<formula>formulaseeoriginaldocumentpage7</formula>对用户k:<formula>formulaseeoriginaldocumentpage7</formula>对用户1:<formula>formulaseeoriginaldocumentpage7</formula>步骤六、各用户依据最大似然法则进行MQAM解码。实施例2:结合图5-图8，针对应用本发明的多用户MIM0-0F匿上行链路系统，满足各用户速率及误比特率要求时，最小化系统发射功率的自适应资源分配问题，进行了MonteCarlo仿真验证，仿真参数设置如表1。表1仿真参数Table1Simulationparameters<table>tableseeoriginaldocumentpage7</column></row><table>8(分成4组)16264160，1，2，4，6，8最小化系统发射功率的数学模型为arg,nZZ尸附subjectto:Z《=&(1)臓^層一对各子载波的共享用户应用所提发明的基于零空间交集的连续迫零算法结合基于连续干扰消除的多用户检测方法，使得各用户的共道干扰得到消除，从而共享同一子载波的各用户发射功率满足4-g(BER一K、(2)即某用户的比特功率分配与其他共享用户的比特功率分配无关，仅与其信道矩阵及小于其序号的共享用户的信道矩阵有关，使复杂的多用户资源分配问题转化为简单的多个并行独立单用户的资源分配问题，数学模型为argmin五《尸二附=1subjectto:Zm^《=&(3)層气S脂一对于MQAM调制方式，发射功率有如下关系4=(rCT2/4)(2《-i)(4)其中，Amk为等效信道增益，由步骤三求得；o2是噪声功率；T为一定BER下的信噪比差(SNRg即)，调整T可改变Pmk，从而满足不同BER要求。图5给出子载波独占方式、基本迫零共享方式、所提发明方法不应用按用户空间相关性分组(SIC-ZF)及应用按用户空间相关性分组(SIC-GZF)，在各用户BER=10—3时，用户平均每子载波的比特数与Eb/No关系图。由图可见，子载波独占方式频谱利用率极低；所提发明方法在不应用按用户空间相关性分组时，在不同频谱利用率下，系统性能相对基本迫零方法均有10dB以上的改进；所提发明方法在应用按用户空间相关性分组时，在频谱利用率范围内，系统性能又有很大的改进。该仿真结果验证了所提发明方法由于充分利用了空间资源，在不同频谱利用率需求时，有效解决了用户空间相关性高对系统性能影响大的问题，改进了系统性能。图6给出所提发明方法不应用按用户空间相关性分组(SIC-ZF)，8用户数接收天线数发射天线数[QQ79]子载波数W保护间隔G/MQAM调制比特数在各用户BER二10—s时，各用户平均每子载波的比特数与Eb/No关系图。由图可见，在不同频谱利用率下，用户1、3、5、7的性能分别优于用户2、4、6、8至少5dB。该仿真结果是由于用户2、4、6、8迫零矩阵的计算分别包含了空间高度相关的用户1、3、5、7，使得连续迫零算法的性能变差所导致。另外，由于用户发射功率的动态范围变大，使得各用户的QoS难以保障，因此，当存在空间相关性高的用户时必须应用按用户空间相关性分组。图7给出所提发明方法应用按用户空间相关性分组(SIC-GZF)，在BER二10—3时，各用户平均每子载波的比特数与Eb/No关系图。由图可见，在不同频谱利用率下，各用户的性能相差最大不超过2dB。这是由于通过按用户空间相关性分组，使得同一子载波不存在空间高度相关的共享用户的结果。与图6比较，应用用户分组，同时兼顾了系统性能与各用户的QoS。图8给出应用基本迫零共享方式、所提发明方法不应用按用户空间相关性分组(SIC-ZF)及应用按用户空间相关性分组(SIC-GZF)，在用户平均每子载波比特数为3时，系统BER与Eb/No关系比较。由图可见在不同BER情况下，所提发明方法不应用按用户空间相关性分组时所需功率相对基本迫零方法至少改进10dB，应用按用户空间相关性分组时所需功率又改进了近5dB。该仿真结果进一步验证了所提发明方法由于充分利用了空间资源，有效解决了用户空间性高对系统性能影响大的问题，从而在不同BER需求时，使系统性能得以改进。实施例3:本发明的计算复杂度分析奇异值分解(SVD)是MM0系统的重要方法，对于一个pXq维、秩r的矩阵，SVD的计算复杂度为min(0(pq2)，0(qp2))，对于低秩稠密矩阵，计算复杂度为0(pqr)[3]。因此，SVD的计算复杂度由矩阵的维度决定。对于多用户MIM0-0FDM系统，高阶矩阵SVD的计算复杂度决定了资源分配算法的计算复杂度。对于基本迫零算法，各用户在各子载波的计算复杂度相同。通过对维度为nKX((K-l)nT)的矩阵5k=Hk—HK]应用SVD方法，计算用户k的迫零矩阵Zk，其计算复杂度为0(X((K-l)riT)2)，求得的迫零矩阵Zk的维度为(nK-(K-l)nT)Xrv通过对维度为(nK-(K-l)nT)XnT的矩阵Z,Hk应用SVD方法，计算用户k的加权矩阵，其计算复杂度为0((nK-(K-l)nT)X(nT)2)。因此，应用SVD方法，计算各用户迫零矩阵及加权矩阵的总计算复杂度为KXNcX(0(nKX((K_l)nT)2)+0((nK-(K_l)nT)X(nT)2))。对于所提发明方法，由于应用了基于零空间交集的连续迫零算法，用户迫零矩阵的计算是通过计算用户加权矩阵时SVD的计算结果迭代相乘求得，其计算复杂度相对于利用SVD方法计算用户加权矩阵的计算复杂度可以忽略。对于用户k，通过对维度为(nK-(k-l)nT)XnT的矩阵ZkHk应用SVD方法，计算用户的加权矩阵，其计算复杂度为0((nK-(k-l)nT)X()2)，因此，所有用户的总计算复杂度为<formula>formulaseeoriginaldocumentpage9</formula>通过比较基本迫零算法与所提发明方法的所有用户SVD的总计算复杂度，可见，当共享子载波的用户数K较大时，基本迫零算法SVD的总计算复杂度远高于所提发明方法的总计算复杂度。实施例4:一种多用户MIM0-0F匿系统上行链路资源分配方法，其技术关键点在于(1)用户按空间相性分组，空间相关性极高的用户分在同一组，保证了各用户组之间不存在空间高度相关的用户；(2)各组内用户由于空间相关性高，采用独占方式分配子载波，各组间用户由于不存在空间相关性极高的用户，采用共享方式分配子载波，解决了共享方式系统性能受用户空间相关性影响大的问题，改进了多用户系统性能；(3)对于各组间共享同一子载波的用户之间产生的共道干扰，利用波束形成技术，通过基于零空间交集的连续迫零算法，实现各用户共道干扰抑制的同时完成各用户子空间的分配。基于零空间交集的连续迫零算法对多用户系统的任意子载波，将按用户空间相关性分组后共享该子载波的用户编号为1K。对任一用户k，假定其后序号的用户信号已知(即对用户k的干扰在基站端可以消除)，为避免其前序号的用户对用户k的干扰，用户k的迫零矩阵应为其前序号所有用户的零空间。基于零空间交集定理，通过计算两个低阶矩阵零空间交集的迫零方法，应用中间计算结果，连续递归计算各用户的迫零矩阵，通过递归迭代使用户迫零矩阵的计算复杂度降低。基于零空间交集的连续迫零算法的用户迫零矩阵的计算为按用户编号的正序，在各用户SVD的计算过程中，一次SVD的计算结果，不仅确定了用户的迫零矩阵，同时确定了每个用户的发射加权矩阵W、接收合并矩阵R、及等效信道增益A，实现了子空间的分配，进一步降低了整个算法的复杂度。(4)各用户跟据连续迫零算法计算的等效信道增益，自适应比特功率分配。对每个用户，以单用户方式对所分配的所有子载波进行比特、功率分配。(5)基于零空间交集的连续迫零算法结合基于连续干扰消除的多用户检测方法利用基于零空间交集的连续迫零算法的计算过程中求得的各用户的接收合并矩阵及发射加权矩阵，按用户编号的逆序计算用户信号。计算任一用户的信号，首先在接收端减去所有已经计算出的大于该用户序号的用户信号经过信道传输后的信号，然后利用该用户的加权矩阵实现线性信号分离；(6)多用户MIMO-OFDM系统上行链路资源分配方法的流程步骤。实施例5:结合图3-图4，一种多用户MIM0-0F匿系统上行链路资源分配方法，它的实现方法为用户按空间相关性分组，相关性极高的用户分在一组；每组内用户依据比例公平原则按独占方式分配子载波，任意子载波每组有且仅有一个用户占用；由于用户分组，任意子载波不存在空间相关性过高的用户，各组间共享子载波的用户利用基于零空间交集的连续迫零算法，实现共道干扰抑制；依据各用户等效信道增益以单用户方式进行自适应比特、功率分配；利用连续干扰消除方法实现多用户检测，进而正确解调。基于该发明的多用户MIM0-0FDM上行链路系统结构图如图3所示。用户k(k二1，2...K)在子载波m(m=1，2...Nc)的nKXnT维上行链路的信道矩阵为<formula>formulaseeoriginaldocumentpage10</formula><formula>formulaseeoriginaldocumentpage10</formula>其中hi,/'"1为第j发射天线到第i接收天线间的信道增益，Amk为由DOA生成的导引向量矩阵。cC、Wmk、Rmk分别为用户k在子载波m的发送数据流、发射加权向量、接收合并向量。为书写方便，以下省略子载波序号m。该发明的核心技术点是基于零空间交集的连续迫零算法结合基于连续干扰消除的多用户检测，其原理框图如图4所示。对多用户系统的任意子载波，将按用户空间相关性分组后，共享该子载波的用户编号为1K。对任一用户k，假定其后序号的用户信号已知(即对用户k的干扰在基站端可以消除)，为避免其前序号的用户对用户k的干扰，用户k的迫零矩阵应为其前序号所有用户的零空间。基于零空间交集的连续迫零算法结合连续干扰消除的多用户检测包括连续迫零及连续干扰消除两个过程，其核心是如何确定每个用户的发射加权矩阵W、及接收合并矩阵R。零空间交集定理令AGCmXn，BGCmXp，Z是B的零空间，X是ZA的零空间；则XZ是[AB]的零空间。依据零空间交集定理，一个高阶矩阵的零空间可以通过将该矩阵按列分成两个行数相同的低阶矩阵，通过递归计算两个低阶矩阵零空间交集的连续迫零的方法得到。对任一用户k，令Hk=[Hk—凡—2.HJ=[Hk—瓦—J，Zk—!为5k—!的零空间，Xk—!为Zk—凡—！的零空间，则Zk二Xk—A—工为!U勺零空间，即Zk为用户l至用户(k-1)的零空间交集，满足ZkHj二0(1《j《k-1)。因此，Zk可通过递归迭代求得，使得用户迫零矩阵的计算复杂度降低。基于零空间交集定理的连续迫零算法，递归计算每个用户的迫零矩阵Zk、发射加权矩阵Wk、接收合并矩阵Rk。对用户1A=1，&=[U/"，U/。)]AJ,，令Ri=[U/"]h，Wi=V丄；对用户2:Z2=[U/。)]h，Z2H2=[U2(1)，U2(°)]A2V/，令R2=[U2(1)]HZ2，W2=V2;......对用户k:Zk=[Uk—/。)]%—"ZkHk=[Uk("，Uk(。)]AJkH，令Rk=[Uk(1)]HZk，Wk=Vk;其中，Ak为等效信道增益Ak(i)组成的对角阵。基于连续干扰消除的多用户检测a:对任一子载波，基站接收信号为「！>严_/^+"，对任一用户k7=1<formula>formulaseeoriginaldocumentpage11</formula><formula>formulaseeoriginaldocumentpage11</formula><formula>formulaseeoriginaldocumentpage11</formula>因此<formula>formulaseeoriginaldocumentpage11</formula>基于以上推导，可以通过先计算编号最大的用户的信号，按编号递减顺序计算其他用户的信号。即计算任一用户的信号，首先在接收端减去所有已经计算出的大于该用户序号的用户信号经过信道传输后的信号，然后利用该用户加权矩阵实现线性信号分离。因此，基于零空间交集的连续迫零算法的用户迫零矩阵的计算为按用户编号的正序，而基于连续干扰消除的多用户检测的计算为按用户编号的逆序。由于整个递归过程中一次SVD的计算结果，即可以用于迫零矩阵的计算，又可以用于多用户的连续干扰消除，因此，该发明方法的计算复杂度明显降低。该发明方法各用户性能分析由于编号最大的用户在发射端乘的迫零矩阵为小于其编号的所有用户的零空间，其等效增益最小，但由于没有引入其他用户的共道干扰，性能得以改进；由于编号最小的用户最后解调，不需乘迫零矩阵，其等效增益最大，但由于噪声的存在，SIC算法引入了其他用户的部分干扰，使得性能变差。因此，整个算法中用户编号顺序对用户的性能影响不大。实施例6:本发明中对任一用户k，假定其后序号的用户信号已知(即对用户k的干扰在基站端可以消除)，为避免其前序号的用户对用户k的干扰，用户k的迫零矩阵应为其前序号所有用户的零空间。基于零空间交集定理，通过计算两个低阶矩阵零空间交集的迫零方法，应用中间计算结果，连续递归计算各用户的迫零矩阵，通过递归迭代使用户迫零矩阵的计算复杂度降低。基于零空间交集的连续迫零算法的用户迫零矩阵的计算为按用户编号的正序，在各用户SVD的计算过程中，一次SVD的计算结果，不仅确定了用户的迫零矩阵，同时确定了每个用户的发射加权矩阵W、接收合并矩阵R、及等效信道增益A，实现了子空间的分配，进一步降低了整个算法的复杂度。某文献针对多用户MIMO-OFDM下行链路，子载波共享方式高阶矩阵迫零算法计算复杂度高的问题，提出了基于干扰用户零空间交集的低复杂度的迫零矩阵计算方法。该方法本质上仍是针对所有干扰用户信道矩阵组成的高阶矩阵进行SVD计算，只是将直接计算高阶矩阵的零空间转化为间接计算多个低阶矩阵的零空间交集，干扰用户零空间交集计算的中间计算结果没有利用，各用户迫零矩阵的计算复杂度相同。1权利要求一种多用户MIMO-OFDM上行链路资源分配方法，其特征在于它由以下步骤实现步骤一用户按空间相关性分组；空间相关性高的用户分在一组，DOA夹角小于5度的用户分在同一组；步骤二每组内用户按子载波独占方式分配子载波；依据比例公平原则，对组内用户以独占方式分配子载波，任意子载波每组有且仅有一个用户占用。定义子载波分配因子wmk，为1表示用户k占用子载波m，为0表示用户k不占用子载波m；步骤三各组间共享同一子载波的用户利用波束形成技术，基于零空间交集的连续迫零方法，实现共道干扰抑制的同时完成子空间的分配。对每个子载波，将共享该子载波的用户编号为1～K，利用基于零空间交集的连续迫零方法，按如下步骤，以用户编号的正序计算各用户的迫零矩阵Z、发射加权矩阵W、接收合并矩阵R；对用户1Z1＝1，H1＝[U1(1)，U1(0)]Λ1V1H，令R1＝[U1(1)]H，W1＝V1；对用户2Z2＝[U1(0)]HZ1，Z2H2＝[U2(1)，U2(0)]Λ2V2H，令R2＝[U2(1)]HZ2，W2＝V2；对用户3Z3＝[U2(0)]HZ2，Z3H3＝[U3(1)，U3(0)]Λ3V3H，令R3＝[U3(1)]HZ3，W3＝V3；......对用户kZk＝[Uk-1(0)]HZk-1，ZkHk＝[Uk(1)，Uk(0)]ΛkVkH，令Rk＝[Uk(1)]HZk，Wk＝Vk；......对用户KZK＝[UK-1(0)]HZK-1，ZKHK＝[UK(1)，UK(0)]ΛKVKH，令RK＝[UK(1)]HZK，WK＝VK；步骤四、依据等效信道增益，自适应比特功率分配；对每个用户，以单用户方式对所分配的所有子载波进行比特、功率分配，典型算法为贪婪算法；步骤五、基站端基于连续干扰消除(SIC)方法进行多用户检测；基站端对接收的多用户的信号和进行分离，计算出每个共享用户的信号。对每个子载波，按如下步骤，以用户编号的逆序计算各共享用户的接收信号；对用户K<mrow><msub><mi>R</mi><mi>K</mi></msub><mi>Y</mi><mo>=</mo><msup><mrow><mo>[</mo><msup><msub><mi>U</mi><mi>K</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>]</mo></mrow><mi>H</mi></msup><msub><mi>Z</mi><mi>K</mi></msub><mrow><mo>(</mo><msub><mi>H</mi><mi>K</mi></msub><msub><mi>W</mi><mi>K</mi></msub><msub><mi>X</mi><mi>K</mi></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>K</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>H</mi><mi>j</mi></msub><msub><mi>W</mi><mi>j</mi></msub><msub><mi>X</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>R</mi><mi>K</mi></msub><mi>n</mi><mo>=</mo><msub><mi>&Lambda;</mi><mi>K</mi></msub><msub><mi>X</mi><mi>K</mi></msub><mo>+</mo><msub><mi>R</mi><mi>K</mi></msub><mi>n</mi></mrow><mrow><msub><mover><mi>X</mi><mo>^</mo></mover><mi>K</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>&Lambda;</mi><mi>K</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>R</mi><mi>K</mi></msub><mi>Y</mi></mrow>......对用户k<mrow><msub><mi>R</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>H</mi><mi>j</mi></msub><msub><mi>W</mi><mi>j</mi></msub><msub><mi>X</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msup><msub><mi>U</mi><mi>k</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>]</mo></mrow><mi>H</mi></msup><msub><mi>Z</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>H</mi><mi>k</mi></msub><msub><mi>W</mi><mi>k</mi></msub><msub><mi>X</mi><mi>k</mi></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>H</mi><mi>j</mi></msub><msub><mi>W</mi><mi>j</mi></msub><msub><mi>X</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>R</mi><mi>k</mi></msub><mi>n</mi><mo>=</mo><msub><mi>&Lambda;</mi><mi>k</mi></msub><msub><mi>X</mi><mi>k</mi></msub><mo>+</mo><msub><mi>R</mi><mi>k</mi></msub><mi>n</mi></mrow><mrow><msub><mover><mi>X</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>&Lambda;</mi><mi>k</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>R</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>H</mi><mi>j</mi></msub><msub><mi>W</mi><mi>j</mi></msub><msub><mi>X</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow>......对用户l<mrow><msub><mi>R</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>2</mn></mrow><mi>K</mi></munderover><msub><mi>H</mi><mi>j</mi></msub><msub><mi>W</mi><mi>j</mi></msub><msub><mi>X</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msup><msub><mi>U</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>]</mo></mrow><mi>H</mi></msup><mrow><mo>(</mo><msub><mi>H</mi><mn>1</mn></msub><msub><mi>W</mi><mn>1</mn></msub><msub><mi>X</mi><mn>1</mn></msub><mo>+</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&Lambda;</mi><mn>1</mn></msub><msub><mi>X</mi><mn>1</mn></msub><mo>+</mo><msub><mi>R</mi><mn>1</mn></msub><mi>n</mi></mrow><mrow><msub><mover><mi>X</mi><mo>^</mo></mover><mn>1</mn></msub><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>&Lambda;</mi><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>R</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>2</mn></mrow><mi>K</mi></munderover><msub><mi>H</mi><mi>j</mi></msub><msub><mi>W</mi><mi>j</mi></msub><msub><mi>X</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow>步骤六、各用户依据最大似然法则进行MQAM解码。全文摘要本发明的目的在于提供一种适用于多用户MIMO-OFDM系统上行链路子载波、子空间、比特、功率资源分配优化过程的多用户MIMO-OFDM上行链路资源分配方法。它由以下步骤实现用户按空间相关性分组；每组内用户按子载波独占方式分配子载波；各组间共享同一子载波的用户利用波束形成技术；依据等效信道增益，自适应比特功率分配；基站端基于连续干扰消除(SIC)方法进行多用户检测；各用户依据最大似然法则进行MQAM解码。本发明在提高多用户MIMO-OFDM系统频谱利用率等性能的同时，降低了资源分配算法的计算复杂度。本发明适用于多用户MIMO-OFDM系统上行链路子载波、子空间、比特、功率等资源分配的优化过程。文档编号H04B1/707GK101765123SQ200910073389公开日2010年6月30日申请日期2009年12月11日优先权日2009年12月11日发明者张成文申请人:哈尔滨工业大学