«上一篇
文章快速检索     高级检索
下一篇»
  哈尔滨工程大学学报  2019, Vol. 40 Issue (4): 827-833  DOI: 10.11990/jheu.201708044
0

引用本文  

禹永植, 张未坤, 郜丽鹏. 多用户MIMO中继系统包含直传链路的联合预编码[J]. 哈尔滨工程大学学报, 2019, 40(4): 827-833. DOI: 10.11990/jheu.201708044.
YU Yongzhi, ZHANG Weikun, GAO Lipeng. Joint precoding for multi-user MIMO relay system with direct links[J]. Journal of Harbin Engineering University, 2019, 40(4): 827-833. DOI: 10.11990/jheu.201708044.

基金项目

国家自然科学基金项目(61571146)

通信作者

禹永植, E-mail:yuyongzhi@hrbeueducn

作者简介

禹永植, 男, 讲师, 博士;
郜丽鹏, 男, 教授, 博士

文章历史

收稿日期:2017-08-18
网络出版日期:2018-09-28
多用户MIMO中继系统包含直传链路的联合预编码
禹永植 , 张未坤 , 郜丽鹏     
哈尔滨工程大学 信息与通信工程学院, 黑龙江 哈尔滨 150001
摘要:针对现有的联合预编码算法忽略多用户之间的直传链路而导致通信性能下降的问题,本文以最小和均方误差为设计准则,提出了一种发射端的所有用户集中优化的联合预编码算法。首先,由于MSMSE优化问题的非凸性,将单流信号输入的多用户等效为一个多流信号输入的基站进行优化处理,其中基站的预编码矩阵由所有发射端的用户的子预编码矩阵构成,原非凸的MSMSE优化问题被转化为发射端用户的子预编码矩阵和中继收发矩阵2个子优化问题分别进行处理。其次,根据最小均方误差接收滤波理论,直接求解出接收滤波矩阵的优化表达式。发射端用户的子预编码矩阵的子优化问题通过求解条件获得,中继收发矩阵的子优化问题被转化成一个标准的半正定问题进行优化处理。最后,发射端用户的子预编码矩阵、中继收发矩阵和接收滤波矩阵联合迭代至收敛,得到优化后的预编码矩阵。实验仿真结果表明:提出的联合预编码算法在系统和均方误差性能和误码率性能上都有明显的提升,并且仿真进一步验证了算法有快速的收敛速度。
关键词发射/接收端多用户    直传链路    SMSE    BER    K.K.T条件    SDP问题    
Joint precoding for multi-user MIMO relay system with direct links
YU Yongzhi , ZHANG Weikun , GAO Lipeng     
College of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China
Abstract: To solve the problem of existing precoding algorithms that fail to improve communication performance because they ignore the direct links between multiple users, we design a joint precoding algorithm based on minimum sum mean-square-error (MSMSE) design criteria.In this algorithm, all users at the transmitter are optimized in a centralized manner.First, as the MSMSE optimization problem is non-convex, all multiple users of a single-stream signal input are equivalent to a base station of multi-stream signal input for optimization, in which the precoding matrix of the base station consists of sub-precoding matrixes of all the transmitting users.The original problem is transformed into two sub-optimization problems to be solved:transmitter user's sub-precoding matrix and relay transmitter/receiver matrix.Second, according to the minimum mean-square-error receiving filter theory, expressions of the receiving filter matrixes are derived directly.The optimal expressions of the transmitter sub-precoding matrices are solved under Karush-Kuhn-Tucker conditions.To optimize the relay transmitter/receiver matrix, we transform the optimization problem into a standard semi-definite programming problem.Finally, the transmitter sub-precoding matrices, relay amplifying matrix, and receiving filter matrixes are optimized iteratively, and the solution is obtained after the optimization of the algorithm.Numerical results show that the proposed algorithm offers good performance in sum mean-square-error and bit-error rate.In addition, the algorithm exhibits fast convergence.
Keywords: multiple users of transmitter/receiver    direct links    minimun sum mean-square-error (SMSE)    bit-error rate (BER)    Karush-Kuhn-Tucker (KKT) conditions    semi-definite programming (SDP) problem    

MIMO技术由于其分集传输和多路复用的特点,在频谱利用率和链接的可靠性方面性能都有显著的提高[1]。而中继协作MIMO技术自从被Sendonaris等首次提出后,因其具有提高系统容量、频谱利用率和增大系统覆盖范围等特点吸引了无线通信领域研究者们的兴趣。另外,放大转发(amplify-and-forward, AF)中继协议因其具有较低的计算复杂度和较小的传输时延而受到较多的关注[2]。因此,结合中继技术和MIMO技术的MIMO中继技术被认为是下一代无线通信领域最具有应用前景的技术之一。特别地,因为考虑直传链路的多用户MIMO AF中继系统可以明显地提高通信质量,所以,越来越多的学者对此产生兴趣。

近年来,多用户MIMO中继系统中的预编码算法的设计问题一直备受关注,预编码算法是解决无线通信中多用户MIMO信道干扰问题的主要方法,也是提高无线通信系统性能的关键技术之一。对于多用户MIMO AF中继系统,Khandaker等[3]针对上行的多用户MIMO中继系统,提出一种基于MMSE设计准则的预编码算法,联合优化了发射端,中继节点和接收端的预编码矩阵;基于最小SMSE设计准则且在没有考虑中继节点的情况下,文献[4]提出一种联合发射预编码矩阵和接收滤波矩阵的迭代算法。文献[3]和[4]都是在发射端为多用户的情况下展开的研究。然而,在满足功率限制的条件下,文献[5]在由一个基站、一个中继节点和多个用户组成的系统模型下,提出了一种联合预编码算法。另外,文献[5]还研究了基站的预编码设计和中继节点的预编码设计这2种特殊情况。为了减小计算复杂度,文献[6]对于多用户对MIMO AF中继系统,提出了2种子优化非迭代设计方案。然而,上述提到的文献均没有考虑发射端和接收端用户之间的直传链路。到目前为止,仅有少量的文献考虑直传链路对系统性能的贡献。在基于MMSE设计准则下,文献[7-8]针对包含直传链路的三节点双向的MIMO AF中继系统,提出了一种预编码算法,但是算法不适合应用到多用户MIMO AF中继系统中。文献[9]基于系统容量和MSE这2个设计准则,设计了一种考虑直传链路的联合预编码算法。不同于文献[9],在包含直传链路的多用户MIMO AF中继系统中,一种加权合速率(weighted sun rate, WSR)最大化的联合预编码算法在文献[10]中被提出,算法考虑将最初的优化问题等效为容易处理的问题模型,且实现简单。

本文考虑了所有通信节点均配备多天线且包含直传链路的发射/接收端多用户MIMO AF中继系统模型,在所有发射端的用户和中继节点满足功率限制条件下,研究了一种发射端的用户集中优化的联合预编码问题。目前,没有文献对包含直传链路的发射/接收端多用户MIMO AF中继系统的预编码算法进行研究,而且多用户之间的通信更加具有实际的应用意义。基于这些,本文首先将系统所有发射端的用户看成一个整体,即将发射单流信号的多用户等效为一个发射多流信号的基站,以一种集中优化的方式处理所有发射端用户的子预编码矩阵。其次,提出一种有效的联合预编码算法,结合标准的SDP凸优化模型和迭代思想,简化了算法设计的数学推导过程,算法交替迭代至所有预编码矩阵的优化值。最后,本文提出的算法与文献[6]的2个非迭代设计和文献[10]的迭代算法进行性能比较。

1 多用户MIMO AF中继系统模型

本文考虑了一个包含直传链路的多用户MIMO AF中继系统模型,模型由K个发射端的用户、K个接收端的用户和一个中继节点组成。如图 1所示,发射端的用户和接收端的用户配备相同的天线数为Nk(k=1, 2, …, K),中继节点配备Nr根天线。在该通信系统中,第k个发射端的用户通过中继节点或者直接将信息发送给第k个接收端的用户。对于任意一个接收端的用户来说,最终在接收端接收到的信息不仅会有相应的发射端的用户发送的有用信息,还会有来自发射端其他用户的干扰信息。因此,本文将所有发射端的用户等效为一个传输多信息流的基站,基站将所有发射端的用户发送的单信息流进行同时传输。为简化分析,本文假设系统工作在半双工模式下,信息由发射端的用户通过2个传输时隙发送到接收端的用户,并假设每个时隙内的信道状态信息(channel state information, CSI)保持不变。

Download:
图 1 包含直传链路的多用户MIMO AF中继系统模型 Fig. 1 Multiuser MIMO AF relay system model with direct links

在第1个传输时隙内,第k个发射端的用户将原始的单流信号skCNk×1(k=1, 2, …, K)经过发射预编码矩阵${\mathit{\boldsymbol{B}}_k} \in {\rm{ }}{{\bf{C}}^{{N_b} \times {N_k}}}({N_b} = \sum\limits_{k = 1}^K {{N_k}} )$处理后得到发送信号xk=BkskCNb×1,信号xk被同时发送给中继节点和第k个接收端的用户。其中,原始信号sk满足的功率限制条件为Rsk=Ε[skskH]=INk。令${\mathit{\boldsymbol{s}}^{\rm{T}}} \buildrel \Delta \over = {\rm{ }}[\mathit{\boldsymbol{s}}_1^{\rm{T}} \cdots \mathit{\boldsymbol{s}}_k^{\rm{T}} \cdots {\rm{ }}\mathit{\boldsymbol{s}}_K^{\rm{T}}] \in {{\bf{C}}^{1 \times {N_b}}}$为所有发射端的用户的信号,即等效基站的原始多流信号,且信号s满足的功率限制条件为Ε[ssH]=INb。第1个时隙内,中继节点和第k个接收端的用户接收到的信号分别为yrCNr×1yk1CNk×1,具体表达式为:

$ {\mathit{\boldsymbol{y}}_r} = \mathit{\boldsymbol{HBs}} + {\mathit{\boldsymbol{n}}_r} = \mathit{\boldsymbol{H}}\sum\limits_{i = 1,i \ne k}^K {{\mathit{\boldsymbol{B}}_k}{\mathit{\boldsymbol{s}}_k}} + {\mathit{\boldsymbol{n}}_r} $ (1)
$ \mathit{\boldsymbol{y}}_k^1 = \underbrace {{\mathit{\boldsymbol{T}}_k}{\mathit{\boldsymbol{B}}_k}{\mathit{\boldsymbol{s}}_k}}_{有用信号} + \underbrace {{\mathit{\boldsymbol{T}}_k}\sum\limits_{i = 1,i \ne k}^K {{\mathit{\boldsymbol{B}}_i}{\mathit{\boldsymbol{s}}_i}} }_{干扰信号} + \mathit{\boldsymbol{n}}_k^1 $ (2)

式中:HCNr×Nb是等效后的基站到中继节点的MIMO信道矩阵;$\mathit{\boldsymbol{B}} \buildrel \Delta \over = [{\mathit{\boldsymbol{B}}_1}, \cdots, {\mathit{\boldsymbol{B}}_k}, \cdots, {\mathit{\boldsymbol{B}}_K}] \in {{\bf{C}}^{{N_b} \times {N_b}}}$是等效后基站的预编码矩阵,矩阵B由所有发射端用户的子预编码矩阵构成。等效后的基站到第k个接收端用户之间的MIMO信道矩阵为TkCNk×Nb(k=1, 2, …, K)。另外,nrCNr×1nk1CNk×1分别为中继节点和第k个接收端用户处的复加性高斯白噪声(AWGN),且满足均值为零、方差矩阵分别为Ε[nrnrH]=σr2INr和Ε[nk1(nk1)H]=σk2INkσr2σk2分别为中继节点和第k个接收端用户处的噪声功率。

在第2个传输时隙内,中继节点通过中继预处理矩阵FCNr×Nr对接收信号yr进行处理,并转发给所有的接收端用户,第k个接收端的用户接收到的信号yk2可表示为:

$ \mathit{\boldsymbol{y}}_k^2 = \underbrace {{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{FH}}{\mathit{\boldsymbol{B}}_k}{\mathit{\boldsymbol{s}}_k}}_{有用信号} + \underbrace {\sum\limits_{i = 1,i \ne k}^K {{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{FH}}{\mathit{\boldsymbol{B}}_i}{s_i}} }_{干扰信号} + {\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{F}}{\mathit{\boldsymbol{n}}_r} + \mathit{\boldsymbol{n}}_k^2 $ (3)

式中:GkCNk×Nr(k=1, 2, …, K)为中继节点到第k个接收端用户的MIMO信道矩阵;nk2CNk×1为第k个接收端用户在第2个传输时隙内的复AWGN,且满足均值为零、方差矩阵为Ε[nk2(nk2)H]=σk2INkPr为中继节点处最大的发射功率。那么,中继节点处的功率限制条件可具体表达为:

$ {\rm{tr}}\left( {\mathit{\boldsymbol{F}}\left( {\mathit{\boldsymbol{HB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}{\mathit{\boldsymbol{H}}^{\rm{H}}} + \sigma _r^2{\mathit{\boldsymbol{I}}_{{N_r}}}} \right){\mathit{\boldsymbol{F}}^{\rm{H}}}} \right) \le {P_r} $ (4)

结合表达式(2)和(3),第k个接收端的用户在2个传输时隙内接收的信号yk1yk2可以表达为:

$ \begin{array}{*{20}{c}} {{\mathit{\boldsymbol{y}}_k} = \left[ {\begin{array}{*{20}{c}} {\mathit{\boldsymbol{y}}_k^1}\\ {\mathit{\boldsymbol{y}}_k^2} \end{array}} \right] = \underbrace {\left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{T}}_k}}\\ {{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{FH}}} \end{array}} \right]{\mathit{\boldsymbol{B}}_k}{\mathit{\boldsymbol{s}}_k}}_{有用信号} + \underbrace {\sum\limits_{i = 1,i \ne k}^K {\left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{T}}_k}}\\ {{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{FH}}} \end{array}} \right]{\mathit{\boldsymbol{B}}_i}{\mathit{\boldsymbol{s}}_i}} }_{干扰信号} + }\\ {\underbrace {\left[ {\begin{array}{*{20}{c}} \mathit{\boldsymbol{I}}&{\bf{0}}&{\bf{0}}\\ {\bf{0}}&{{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{F}}}&\mathit{\boldsymbol{I}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\mathit{\boldsymbol{n}}_k^1}\\ {{\mathit{\boldsymbol{n}}_r}}\\ {\mathit{\boldsymbol{n}}_k^2} \end{array}} \right]}_{噪声}} \end{array} $ (5)

式中:本文定义${\mathit{\boldsymbol{\overline H}} _k} = {[{\mathit{\boldsymbol{T}}_k}\;\;{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{FH}}]^{\rm{T}}}$ (k=1, 2, …, K)为等效后的基站到第k个接收端用户之间的等效信道矩阵;${\mathit{\boldsymbol{\overline n}} _k} = {[\mathit{\boldsymbol{n}}_k^1\;\;{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{F}}{\mathit{\boldsymbol{n}}_r} + {\boldsymbol{n}}_k^2]^{\rm{T}}}$ (k=1, 2, …, K)为等效后的基站到第k个接收端用户之间的等效噪声。那么,表达式(5)可进一步表达为:

$ {\mathit{\boldsymbol{y}}_k} = {{\mathit{\boldsymbol{\bar H}}}_k}{\mathit{\boldsymbol{B}}_k}{\mathit{\boldsymbol{s}}_k} + \sum\limits_{i = 1,i \ne k}^K {{{\mathit{\boldsymbol{\bar H}}}_k}{\mathit{\boldsymbol{B}}_i}{\mathit{\boldsymbol{s}}_i}} + {{\mathit{\boldsymbol{\bar n}}}_k} $ (6)

在第k个接收端用户处用滤波矩阵WkCNk×2Nk(k=1, 2, …, K)对信号sk进行估计:

$ {{\mathit{\boldsymbol{\bar s}}}_k} = {\mathit{\boldsymbol{W}}_k}{\mathit{\boldsymbol{y}}_k} $ (7)

式中第k个接收端用户处的滤波矩阵Wk由第1个传输时隙的滤波矩阵W1kCNk×Nk和第2个传输时隙的滤波矩阵W2kCNk×Nk组成。根据文献[11]的MMSE接收滤波理论,第k个接收端用户处的滤波矩阵Wk的表达式可表达为:

$ \begin{array}{*{20}{c}} {{\mathit{\boldsymbol{W}}_k} = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{W}}_{1k}}}&{{\mathit{\boldsymbol{W}}_{2k}}} \end{array}} \right] = \arg \mathop {\min }\limits_{{W_k}} {\rm{E}}\left[ {\left\| {{\mathit{\boldsymbol{W}}_k}{\mathit{\boldsymbol{y}}_k} - {\mathit{\boldsymbol{s}}_k}} \right\|_2^2} \right] = }\\ {\mathit{\boldsymbol{B}}_k^{\rm{H}}\mathit{\boldsymbol{\bar H}}_k^{\rm{H}}{{\left( {{{\mathit{\boldsymbol{\bar H}}}_k}\mathit{\boldsymbol{B}}{\mathit{\boldsymbol{B}}^{\rm{H}}}\mathit{\boldsymbol{\bar H}}_k^{\rm{H}} + {\mathit{\boldsymbol{Z}}_k}\mathit{\boldsymbol{Z}}_k^{\rm{H}}} \right)}^{ - 1}}} \end{array} $ (8)

式中:ZkZkH=diag(σk2INk, σk2INk+σr2GkFFHGkH)。根据表达式(5)和表达式(8),第k个接收端用户处的均方误差矩阵MSEk的具体表达式为:

$ \begin{array}{*{20}{c}} {{\rm{MS}}{{\rm{E}}_k} = {\rm{E}}\left[ {\left( {{\mathit{\boldsymbol{W}}_k}{\mathit{\boldsymbol{y}}_k} - {\mathit{\boldsymbol{s}}_k}} \right){{\left( {{\mathit{\boldsymbol{W}}_k}{\mathit{\boldsymbol{y}}_k} - {\mathit{\boldsymbol{s}}_k}} \right)}^{\rm{H}}}} \right] = }\\ {{\mathit{\boldsymbol{I}}_{{N_k}}} - {\mathit{\boldsymbol{W}}_k}{{\mathit{\boldsymbol{\bar H}}}_k}{\mathit{\boldsymbol{B}}_k} - \mathit{\boldsymbol{B}}_k^{\rm{H}}\mathit{\boldsymbol{\bar H}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_k^{\rm{H}} + }\\ {{\mathit{\boldsymbol{W}}_k}{{\mathit{\boldsymbol{\bar H}}}_k}\mathit{\boldsymbol{B}}{\mathit{\boldsymbol{B}}^{\rm{H}}}\mathit{\boldsymbol{\bar H}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_k^{\rm{H}} + {\mathit{\boldsymbol{W}}_k}{\mathit{\boldsymbol{Z}}_k}\mathit{\boldsymbol{Z}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_k^{\rm{H}}} \end{array} $ (9)
2 基于SMSE最小化的联合预编码设计

这部分主要设计了一种发射端的用户集中优化的联合预编码算法,在保证系统SMSE最小化的前提下,将所有传输单流信号的发射端用户等效为一个传输多流信号的基站,通过预编码算法获得所有发射端用户的子预编码矩阵Bk,接收端用户的滤波矩阵Wk和中继收发矩阵F的优化值。具体优化问题为:

$ \mathop {\min }\limits_{\scriptstyle{B_k},F,{W_k}\atop {\scriptstyle {k = 1,2, \cdots ,K}}} \sum\limits_{k = 1}^K {{\rm{tr}}\left( {{\rm{MS}}{{\rm{E}}_k}} \right)} $ (10)
$ {\rm{s}}.{\rm{t}}.\;\;\;{\rm{tr}}\left( {{\mathit{\boldsymbol{B}}_k}\mathit{\boldsymbol{B}}_k^{\rm{H}}} \right) \le {P_s} $ (11)
$ {\rm{tr}}\left( {\mathit{\boldsymbol{F}}\left( {\mathit{\boldsymbol{HB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}{\mathit{\boldsymbol{H}}^{\rm{H}}} + \sigma _r^2{\mathit{\boldsymbol{I}}_{{N_r}}}} \right){\mathit{\boldsymbol{F}}^{\rm{H}}}} \right) \le {P_r} $ (12)

显然,对于矩阵变量BkFWk,优化问题(10)~(12)是一个非凸的问题,所以不能直接获得每个优化变量的全局优化值。对此,我们设计了一个有效的发射端的用户集中优化的迭代算法,通过求解下面几个子凸优化问题获得优化变量BkFWk的局部优化表达式。最后,通过算法交替迭代至收敛,获得每个优化变量矩阵的局部优化值。

2.1 发射端用户的子预编码矩阵的优化设计

给定矩阵变量FWk后,将表达式(9)代入到目标函数表达式(10)中,原多优化变量的优化问题,就可转化为只有发射端用户的子预编码矩阵Bk(k=1, 2, …, K)一个优化变量的优化问题,仅以Bk为优化变量的优化方程可以表示为:

$ \begin{array}{*{20}{c}} {\mathop {\min }\limits_{\begin{array}{*{20}{c}} {{B_k}}\\ {k = 1,2, \cdots ,K} \end{array}} \sum\limits_{k = 1}^K {{\rm{tr}}\left( {{\mathit{\boldsymbol{I}}_{{N_k}}} - {\mathit{\boldsymbol{W}}_k}{{\mathit{\boldsymbol{\bar H}}}_k}{\mathit{\boldsymbol{B}}_k} - \mathit{\boldsymbol{B}}_k^{\rm{H}}\mathit{\boldsymbol{\bar H}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_k^{\rm{H}} + } \right.} }\\ {\left. {{\mathit{\boldsymbol{W}}_k}{{\mathit{\boldsymbol{\bar H}}}_k}\mathit{\boldsymbol{B}}{\mathit{\boldsymbol{B}}^{\rm{H}}}\mathit{\boldsymbol{\bar H}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_k^{\rm{H}} + {\mathit{\boldsymbol{W}}_k}{\mathit{\boldsymbol{Z}}_k}\mathit{\boldsymbol{Z}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_k^{\rm{H}}} \right)} \end{array} $ (13)
$ {\rm{s}}.{\rm{t}}.{\rm{tr}}\left( {{\mathit{\boldsymbol{B}}_k}\mathit{\boldsymbol{B}}_k^{\rm{H}}} \right) \le {P_s} $ (14)

对于上述在给定矩阵变量FWk下的关于矩阵变量Bk的优化问题,暂时忽略中继节点处的功率限制条件,因为该条件可以通过后续的迭代过程得到满足。显而易见,问题(13)~(14)是关于矩阵变量Bk的凸优化问题。

为了求解上述关于矩阵变量Bk的凸优化问题,本文采用传统的凸优化求解方法K.K.T条件求解其优化值[12]。首先构造有关Bk的拉格朗日函数为:

$ \begin{array}{*{20}{c}} {{\mathit{\boldsymbol{L}}_{{B_k}}} = \sum\limits_{k = 1}^K {{\rm{tr}}\left( {{\mathit{\boldsymbol{I}}_{{N_k}}} - {\mathit{\boldsymbol{W}}_k}{{\mathit{\boldsymbol{\bar H}}}_k}{\mathit{\boldsymbol{B}}_k} - \mathit{\boldsymbol{B}}_k^{\rm{H}}\mathit{\boldsymbol{\bar H}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_k^{\rm{H}} + } \right.} }\\ {\left. {{\mathit{\boldsymbol{W}}_k}{{\mathit{\boldsymbol{\bar H}}}_k}\mathit{\boldsymbol{B}}{\mathit{\boldsymbol{B}}^{\rm{H}}}\mathit{\boldsymbol{\bar H}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_k^{\rm{H}}} \right) + {\lambda _{{B_k}}}\left( {{\rm{tr}}\left( {{\mathit{\boldsymbol{B}}_k}\mathit{\boldsymbol{B}}_k^{\rm{H}}} \right) - {P_s}} \right)} \end{array} $ (15)

式中:λBk≥0(k=1, 2, …, K)是与功率限制条件(14)有关的拉格朗日系数。由于表达式(9)中的WkZkZkHWkH与矩阵变量Bk无关,因此在关于矩阵变量Bk的拉格朗日函数中忽略不考虑。通过拉格朗日乘数法求得的矩阵变量Bk的优化表达式为:

$ {\mathit{\boldsymbol{B}}_k} = \sum\limits_{j = 1}^K {{{\left( {\mathit{\boldsymbol{\bar H}}_j^{\rm{H}}\mathit{\boldsymbol{W}}_j^{\rm{H}}{\mathit{\boldsymbol{W}}_j}{{\mathit{\boldsymbol{\bar H}}}_j} + {\lambda _{{B_k}}}{\mathit{\boldsymbol{I}}_{{N_b}}}} \right)}^{ - 1}}\mathit{\boldsymbol{\bar H}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_k^{\rm{H}}} $ (16)

在求解表达式(16)的过程中,应用了矩阵运算规则:$\frac{\partial \text{tr}({{Z}^{\dagger }}\mathit{\boldsymbol{A}})}{\partial {{Z}^{*}}}=\mathit{\boldsymbol{A}}$$\frac{\partial \text{tr}(\mathit{\boldsymbol{AZ}}{{\mathit{\boldsymbol{Z}}}^{\dagger }}\mathit{\boldsymbol{B}})}{\partial {{Z}^{*}}}=\mathit{\boldsymbol{BAZ}}$。其中,未知的拉格朗日系数λBk通过二分法算法获得。

2.2 中继收发矩阵的优化设计

给定矩阵变量BkWk(k=1, 2, …, K),将表达式(9)代入目标函数表达式(10)中,原优化问题就可转化为只有中继收发矩阵F一个优化变量的问题,仅以F为优化变量的优化方程可进一步转化为:

$ \begin{array}{*{20}{c}} {\mathop {\min }\limits_F \sum\limits_{k = 1}^K {{\rm{tr}}\left( {{\mathit{\boldsymbol{I}}_{{N_k}}} - {\mathit{\boldsymbol{W}}_{1k}}{\mathit{\boldsymbol{T}}_k}{\mathit{\boldsymbol{B}}_k} - \mathit{\boldsymbol{B}}_k^{\rm{H}}\mathit{\boldsymbol{T}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_{1k}^{\rm{H}} - {\mathit{\boldsymbol{W}}_{2k}}{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{FH}}{\mathit{\boldsymbol{B}}_k} - } \right.} }\\ {\mathit{\boldsymbol{B}}_k^{\rm{H}}{\mathit{\boldsymbol{H}}^{\rm{H}}}{\mathit{\boldsymbol{F}}^{\rm{H}}}\mathit{\boldsymbol{G}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_{2k}^{\rm{H}} + {\mathit{\boldsymbol{W}}_{2k}}{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{FHB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}\mathit{\boldsymbol{T}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_{1k}^{\rm{H}} + {\mathit{\boldsymbol{W}}_{1k}}{\mathit{\boldsymbol{T}}_k}}\\ {\mathit{\boldsymbol{B}}{\mathit{\boldsymbol{B}}^{\rm{H}}}{\mathit{\boldsymbol{H}}^{\rm{H}}}{\mathit{\boldsymbol{F}}^{\rm{H}}}\mathit{\boldsymbol{G}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_{2k}^{\rm{H}} + {\mathit{\boldsymbol{W}}_{2k}}{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{FHB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}{\mathit{\boldsymbol{H}}^{\rm{H}}}{\mathit{\boldsymbol{F}}^{\rm{H}}}\mathit{\boldsymbol{G}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_{2k}^{\rm{H}} + }\\ {{\mathit{\boldsymbol{W}}_{1k}}{\mathit{\boldsymbol{T}}_k}\mathit{\boldsymbol{B}}{\mathit{\boldsymbol{B}}^{\rm{H}}}\mathit{\boldsymbol{T}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_{1k}^{\rm{H}} + \sigma _r^2{\mathit{\boldsymbol{W}}_{2k}}{\mathit{\boldsymbol{G}}_k}\mathit{\boldsymbol{F}}{\mathit{\boldsymbol{F}}^{\rm{H}}}\mathit{\boldsymbol{G}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_{2k}^{\rm{H}} + \sigma _k^2}\\ {\left. {{\mathit{\boldsymbol{W}}_{1k}}\mathit{\boldsymbol{W}}_{1k}^{\rm{H}} + \sigma _k^2{\mathit{\boldsymbol{W}}_{2k}}\mathit{\boldsymbol{W}}_{2k}^{\rm{H}}} \right)} \end{array} $ (17)
$ {\rm{s}}.{\rm{t}}.{\rm{tr}}\left( {\mathit{\boldsymbol{F}}\left( {\mathit{\boldsymbol{HB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}{\mathit{\boldsymbol{H}}^{\rm{H}}} + \sigma _r^2{\mathit{\boldsymbol{I}}_{{N_r}}}} \right){\mathit{\boldsymbol{F}}^{\rm{H}}}} \right) \le {P_r} $ (18)

与优化问题(13)~(14)一样,优化问题(17)~(18)是关于矩阵变量F的凸优化问题。Matlab仿真软件下的CVX优化工具箱可以高效地求解复杂的凸优化问题,K.K.T条件在求解凸优化问题时需要做大量的数学推导,算法复杂度高。因为CVX工具箱只可以优化特定的凸优化问题,所以需要将优化问题(17)~(18)转化成一个SDP问题后,再去优化矩阵变量F。首先,根据下面的运算规则[13]

$ {\rm{tr}}\left( {{\mathit{\boldsymbol{C}}^{\rm{T}}}\mathit{\boldsymbol{D}}} \right) = {\left( {{\rm{vec}}\left( \mathit{\boldsymbol{C}} \right)} \right)^{\rm{T}}}{\rm{vec}}\left( \mathit{\boldsymbol{D}} \right) $ (19)
$ {\rm{tr}}\left( {{\mathit{\boldsymbol{C}}^{\rm{H}}}\mathit{\boldsymbol{DCA}}} \right) = {\rm{vec}}{\left( \mathit{\boldsymbol{C}} \right)^{\rm{H}}}\left( {{\mathit{\boldsymbol{A}}^{\rm{T}}} \otimes \mathit{\boldsymbol{D}}} \right){\rm{vec}}\left( \mathit{\boldsymbol{C}} \right) $ (20)

${\mathit{\boldsymbol{M}}_k} \buildrel \Delta \over = {\mathit{\boldsymbol{W}}_{2k}}{\mathit{\boldsymbol{G}}_k}$Vk=W1kTkf=vec(F),关于中继收发矩阵F的目标函数表达式(17)可进行分解为:

$ {\rm{tr}}\left( {{\mathit{\boldsymbol{M}}_k}\mathit{\boldsymbol{F}}{\mathit{\boldsymbol{F}}^{\rm{H}}}\mathit{\boldsymbol{M}}_k^{\rm{H}}} \right) = {\mathit{\boldsymbol{f}}^{\rm{H}}}\left( {{\mathit{\boldsymbol{I}}_{{N_r}}} \otimes \left( {\mathit{\boldsymbol{M}}_k^{\rm{H}}{\mathit{\boldsymbol{M}}_k}} \right)} \right)\mathit{\boldsymbol{f}} $ (21)
$ \begin{array}{*{20}{c}} {{\rm{tr}}\left( {{\mathit{\boldsymbol{M}}_k}\mathit{\boldsymbol{FH}}{\mathit{\boldsymbol{B}}_k}} \right) = {\rm{tr}}\left( {\mathit{\boldsymbol{H}}{\mathit{\boldsymbol{B}}_k}{\mathit{\boldsymbol{M}}_k}\mathit{\boldsymbol{F}}} \right) = }\\ {{{\left( {{\rm{vec}}\left( {{{\left( {\mathit{\boldsymbol{H}}{\mathit{\boldsymbol{B}}_k}{\mathit{\boldsymbol{M}}_k}} \right)}^{\rm{T}}}} \right)} \right)}^{\rm{T}}}\mathit{\boldsymbol{f}}} \end{array} $ (22)
$ \begin{array}{*{20}{c}} {{\rm{tr}}\left( {{\mathit{\boldsymbol{M}}_k}\mathit{\boldsymbol{FHB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}\mathit{\boldsymbol{V}}_k^{\rm{H}}} \right) = {{\left( {{\rm{vec}}\left( {{{\left( {\mathit{\boldsymbol{HB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}\mathit{\boldsymbol{V}}_k^{\rm{H}}{\mathit{\boldsymbol{M}}_k}} \right)}^{\rm{T}}}} \right)} \right)}^{\rm{T}}}\mathit{\boldsymbol{f}}}\\ {{\rm{tr}}\left( {{\mathit{\boldsymbol{M}}_k}\mathit{\boldsymbol{FHB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}{\mathit{\boldsymbol{H}}^{\rm{H}}}{\mathit{\boldsymbol{F}}^{\rm{H}}}\mathit{\boldsymbol{M}}_k^{\rm{H}}} \right) = {\mathit{\boldsymbol{f}}^{\rm{H}}}\left( {{{\left( {\mathit{\boldsymbol{HB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}{\mathit{\boldsymbol{H}}^{\rm{H}}}} \right)}^{\rm{T}}}} \right.}\\ {\left. { \otimes \left( {\mathit{\boldsymbol{M}}_k^{\rm{H}}{\mathit{\boldsymbol{M}}_k}} \right)} \right)\mathit{\boldsymbol{f}}} \end{array} $ (24)

$\mathit{\boldsymbol{R}} \buildrel \Delta \over = \mathit{\boldsymbol{HB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}{\mathit{\boldsymbol{H}}^{\rm{H}}} + \sigma _r^2{\mathit{\boldsymbol{I}}_{{N_r}}}$,中继功率限制条件的表达式(18)左边一项可以通过矩阵运算法则(19)~(20)进一步转化为:

$ {\rm{tr}}\left( {\mathit{\boldsymbol{FR}}{\mathit{\boldsymbol{F}}^{\rm{H}}}} \right) = {\mathit{\boldsymbol{f}}^{\rm{H}}}\left( {{\mathit{\boldsymbol{R}}^{\rm{T}}} \otimes {\mathit{\boldsymbol{I}}_{{N_r}}}} \right)\mathit{\boldsymbol{f}} $ (25)

另外,为简化最终的优化问题,定义下面几个变量替换:

$ \begin{array}{*{20}{c}} {{\mathit{\boldsymbol{U}}_k} \buildrel \Delta \over = \left( {{{\left( {\mathit{\boldsymbol{HB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}{\mathit{\boldsymbol{H}}^{\rm{H}}}} \right)}^{\rm{T}}} + \sigma _r^2{\mathit{\boldsymbol{I}}_{{N_r}}}} \right) \otimes \left( {\mathit{\boldsymbol{M}}_k^{\rm{H}}{\mathit{\boldsymbol{M}}_k}} \right)}\\ {{\mathit{\boldsymbol{a}}_k} \buildrel \Delta \over = {\rm{vec}}\left( {{{\left( {\mathit{\boldsymbol{H}}{\mathit{\boldsymbol{B}}_k}{\mathit{\boldsymbol{M}}_k}} \right)}^{\rm{T}}}} \right)}\\ {{\mathit{\boldsymbol{b}}_k} \buildrel \Delta \over = {\rm{vec}}\left( {{{\left( {\mathit{\boldsymbol{HB}}{\mathit{\boldsymbol{B}}^{\rm{H}}}\mathit{\boldsymbol{V}}_k^{\rm{H}}{\mathit{\boldsymbol{M}}_k}} \right)}^{\rm{T}}}} \right)}\\ {{t_k} \buildrel \Delta \over = {\rm{tr}}\left( {{\mathit{\boldsymbol{I}}_{{N_r}}} - {\mathit{\boldsymbol{W}}_{1k}}{\mathit{\boldsymbol{T}}_k}{\mathit{\boldsymbol{B}}_k} - \mathit{\boldsymbol{B}}_k^{\rm{H}}\mathit{\boldsymbol{T}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_{1k}^{\rm{H}} + {\mathit{\boldsymbol{W}}_{1k}}{\mathit{\boldsymbol{T}}_k}\mathit{\boldsymbol{B}}} \right.}\\ {\left. {{\mathit{\boldsymbol{B}}^{\rm{H}}}\mathit{\boldsymbol{T}}_k^{\rm{H}}\mathit{\boldsymbol{W}}_{1k}^{\rm{H}} + \sigma _k^2{\mathit{\boldsymbol{W}}_{1k}}\mathit{\boldsymbol{W}}_{1k}^{\rm{H}} + \sigma _k^2{\mathit{\boldsymbol{W}}_{2k}}\mathit{\boldsymbol{W}}_{2k}^{\rm{H}}} \right)} \end{array} $ (26)

将表达式(21)~(26)代入到优化问题(17)~(18)中,则与中继收发矩阵F相关的优化问题可进一步转化为:

$ \mathop {\min }\limits_f \sum\limits_{k = 1}^K {\left( {{\mathit{\boldsymbol{f}}^{\rm{H}}}{\mathit{\boldsymbol{U}}_k}\mathit{\boldsymbol{f}} + \mathit{\boldsymbol{a}}_k^{\rm{H}}\mathit{\boldsymbol{f}} + {\mathit{\boldsymbol{f}}^{\rm{H}}}{\mathit{\boldsymbol{a}}_k} + \mathit{\boldsymbol{b}}_k^{\rm{H}}\mathit{\boldsymbol{f}} + {\mathit{\boldsymbol{f}}^{\rm{H}}}{\mathit{\boldsymbol{b}}_k} + {\mathit{\boldsymbol{t}}_k}} \right)} $ (27)
$ {\rm{s}}.{\rm{t}}.{\mathit{\boldsymbol{f}}^{\rm{H}}}\left( {{\mathit{\boldsymbol{R}}^{\rm{T}}} \otimes {\mathit{\boldsymbol{I}}_{{N_r}}}} \right)\mathit{\boldsymbol{f}} \le {P_r} $ (28)

为进一步将非凸的优化问题转化为标准的SDP问题,根据舒尔补定理,引入辅助变量pk(k=1, 2, …, K)且满足pk≥tr(MSEk),基于MSMSE准则的关于矩阵变量F的优化问题可进一步转化为如下等价的SDP问题:

$ \mathop {\min }\limits_{\mathop {f,{p_k},}\limits_{k = 1,2, \cdots ,K} } \sum\limits_{k = 1}^K {{p_k}} $ (29)
$ {\rm{s}}.{\rm{t}}.\left[ {\begin{array}{*{20}{c}} {{P_r}}&{{\mathit{\boldsymbol{f}}^{\rm{H}}}}\\ \mathit{\boldsymbol{f}}&{{{\left( {{\mathit{\boldsymbol{R}}^{\rm{T}}} \otimes {\mathit{\boldsymbol{I}}_{{N_r}}}} \right)}^{ - 1}}} \end{array}} \right] \ge 0 $ (30)
$ \left[ {\begin{array}{*{20}{c}} {{p_k} - \mathit{\boldsymbol{a}}_k^{\rm{H}}\mathit{\boldsymbol{f}} - {\mathit{\boldsymbol{f}}^{\rm{H}}}{\mathit{\boldsymbol{a}}_k} - \mathit{\boldsymbol{b}}_k^{\rm{H}}\mathit{\boldsymbol{f}} - {\mathit{\boldsymbol{f}}^{\rm{H}}}{\mathit{\boldsymbol{b}}_k} - {\mathit{\boldsymbol{t}}_k}}&{{\mathit{\boldsymbol{f}}^{\rm{H}}}\mathit{\boldsymbol{U}}_k^{{\rm{1}}/2}}\\ {\mathit{\boldsymbol{U}}_k^{{\rm{H}}/2}\mathit{\boldsymbol{f}}}&{{\mathit{\boldsymbol{I}}_{{N_r}}}} \end{array}} \right] \ge 0 $ (31)

上述优化问题(29)~(31)是CVX可以解决的标准SDP问题,CVX优化工具箱[14]会识别标准的SDP问题且自动调用内点法对该问题进行求解,简化了以往求解此类优化问题时所需要的大量数学运算。

2.3 联合迭代算法

综上所述,当接收端的用户采用MMSE接收滤波矩阵Wk(k=1, 2, …, K)时,最初的基于MSMSE准则的联合优化问题被转化为2个子优化问题,然后单独优化发射端用户的子预编码矩阵Bk(k=1, 2, …, K)和中继收发矩阵F的优化形式。前面已经具体介绍了每个子凸优化问题的求解思路和设计过程,在包含直传链路的发射/接收端多用户的MIMO AF中继系统中,对于预编码矩阵BkFWk的优化问题,一种基于MSMSE准则的联合迭代预编码算法如下。

基于MSMSE准则的联合迭代算法

1) 初始化:BkF的初始值满足功率限制条件,设最大迭代次数为Imax,迭代终止门限为ε,迭代次数初始化n=1;

2) 给定Bk(n)F(n),根据表达式(8)计算所有接收端用户的接收滤波矩阵Wk(n)

3) 给定F(n)Wk(n),通过求解表达式(16)更新所有发射端用户的子预编码矩阵Bk(n+1)

4) 根据更新后的Bk(n+1)Wk(n),通过求解表达式(29)~(31)的SDP问题获得更新后的中继收发矩阵F(n+1)

5) 判断如果‖Bk(n+1)-Bk(n)2ε且‖F(n+1)-F(n)‖≤ε或者n>Imax,结束迭代;否则,跳转2)继续迭代直到满足收敛条件。

因为上述算法在每次更新WkBkF时,对应的表达式(8)、式(16)和式(29)~(31)都是非负的且单调递减,所以上述联合迭代预编码算法满足收敛条件,且算法总是收敛到优化值点。另外,算法在每一次迭代过程中,主要的算法复杂度式求解表达式(29)~(31)的SDP问题来更新F,然而求解BkWk的算法复杂度与求解SDP问题的复杂度相比可以忽略不计,因此每次迭代过程的算法复杂度可以近似为O((Nr2+1)4.5)[15]

3 仿真结果

在数值仿真过程中,假设所有的MIMO信道(包括HGkTk(k=1, 2, …, K))均服从准静态独立同分布的瑞利衰落模型,所有的噪声(包括nrnk(k=1, 2, …, K))均为零均值的加性白噪声,且所有通信节点的噪声功率归一化,即σr2=σk2=1。为简化分析,在下面的论述中用SNR表示所有通信链路的信噪比,即SNR=SNRsr=SNRrd=SNRsd。根据系统模型中对所有节点配备天线数的情况说明,下面给出具体的仿真参数设置:发射端的用户数和接收端的用户数均为K=4,且所有的用户都配备2根天线,即Nk=2,中继节点配备的天线数为Nr=8。

为进一步验证所提算法性能上的优越性,首先将提出的算法与下面的4种优化算法在系统SMSE性能和系统BER性能上进行了比较分析。其中,4种比较算法如下:

1) MSMSE w/o DL算法:一种基于MSMSE设计准则的联合迭代算法,该算法没有考虑发射端与接收端之间的直传链路,系统模型是多用户通过中继节点的协作通信;

2) WSR算法[10]:一种基于WSR准则的联合迭代算法,该算法考虑发射端与接收端之间的直传链路,系统模型是一个基站和多个用户之间通过中继节点的协作通信;

3) ZF-CI算法[6]:迫零信道反转算法,一种多用户对间点对点中继系统中的非迭代算法。其中,${\mathit{\boldsymbol{F}}^{({\rm{ZF}})}} = \sum\limits_{k = 1}^K {\mathit{\boldsymbol{B}}_k^{\left( {{\rm{ZF}}} \right)}{\mathit{\boldsymbol{W}}_k}}, \mathit{\boldsymbol{B}}_k^{\left( {{\rm{ZF}}} \right)} = \mathit{\boldsymbol{Q}}_k^{\left( {{\rm{ZF}}} \right)}\mathit{\boldsymbol{U}}_k^{\left( {{\rm{ZF}}} \right)}\mathit{\boldsymbol{ \boldsymbol{\varGamma} }}_k^{\left( {{\rm{ZF}}} \right)1/2}$

4) MMSE-CI算法[6]:最小均方误差信道反转算法,一种多用户对间点对点中继系统中的非迭代算法。其中,${\mathit{\boldsymbol{F}}^{({\rm{MMSE}})}} = \sum\limits_{k = 1}^K {\mathit{\boldsymbol{B}}_k^{\left( {{\rm{MMSE}}} \right)}{\mathit{\boldsymbol{W}}_k}}, \mathit{\boldsymbol{B}}_k^{\left( {{\rm{MMSE}}} \right)} = \mathit{\boldsymbol{Q}}_k^{\left( {{\rm{MMSE}}} \right)}\mathit{\boldsymbol{U}}_k^{\left( {{\rm{MMSE}}} \right)}\mathit{\boldsymbol{ \boldsymbol{\varGamma} }}_k^{\left( {{\rm{MMSE}}} \right)1/2}$

图 2为包含直传链路的发射/接收端多用户MIMO AF中继系统中的SMSE性能随SNR变化的算法比较曲线图。从仿真图中可以明显地看出,本文提出的迭代算法在整个SNR范围内,SMSE性能最优,并且随着SNR的增加,算法的SMSE值接近于0。在所有的算法仿真曲线中,WSR算法与提出算法的SMSE性能相近,文献[6]中的非迭代ZF-CI算法的SMSE性能最差。在低SNR情况下,例如SNR=0时,提出算法的SMSE值低于ZF-CI算法的SMSE值约为2.5。另外,从仿真图中可以得出2个结论:1)提出算法与没有考虑直传链路的MSMSE w/o DL算法相比,算法可以明显地降低系统的SMSE值;2)尽管文献[6]中提出的2个非迭代算法,ZF-CI算法和MMSE-CI算法减小了迭代算法的计算复杂度,但它们的系统SMSE性能曲线却明显高于提出的算法,尤其在SNR取值较低的情况下。

Download:
图 2 本文算法与其他算法的SMSE性能比较 Fig. 2 SMSE comparisons between the proposed algorithm and the other algorithms

图 3为提出算法与其他4种算法在低SNR情况下的BER性能曲线图,信号调制方式为QPSK。由仿真结果图可知,提出算法的系统BER性能明显优于其他4种算法的BER性能,其他4种算法的BER性能相近,并且提出算法的BER性能随着SNR的增加显著提高。

Download:
图 3 本文算法与其他算法的BER性能比较(QPSK) Fig. 3 BER comparisons between the proposed algorithm and the other algorithms (QPSK modulation)

随着SNR值的增加系统的SMSE值在减小,并且SNR取值的不同将直接影响着算法的收敛速度。例如,在低SNR情况下,提出算法在6次迭代以后收敛;在高SNR的情况下,提出算法在4次迭代以内就可以达到收敛。因此,提出算法仅仅需要较少的迭代次数就可以实现较好的SMSE性能,同时也表明了提出算法有较低的计算复杂度和较小的处理延迟。

在不同的SNR情况下,提出算法的SMSE性能的收敛情况如图 4所示。仿真图给出了提出的算法分别在SNR=10 dB、SNR=20 dB和SNR=30 dB 3种情况下,系统SMSE性能随算法迭代次数变化的仿真曲线。从仿真图中可以看出,随着SNR值的增加系统的SMSE值在减小,即系统的SMSE值与SNR值成反比。并且SNR取值的不同将直接影响着算法的收敛速度,在低SNR情况下,提出算法在6次迭代以后收敛;在高SNR的情况下,提出算法在4次迭代以后能够快速地收敛。因此,提出算法仅仅需要较少的迭代次数就可以实现较好的SMSE性能。同时也表明了提出算法有较低的计算复杂度和较小的延迟。

Download:
图 4 不同SNR条件下,本文算法的SMSE性能随迭代次数变化 Fig. 4 SMSE of the proposed algorithm versus number of iterations under different SNR

为进一步验证提出算法的优越性,在不同用户数的情况下,即K的取值分别为K=2和K=4。图 5图 6分别给出了提出算法的系统SMSE性能随着SNR变化的曲线图和提出算法的系统SMSE性能随着迭代次数变化的曲线图。从图 5中可以得出,当SNR增加时,系统的SMSE性能曲线有明显的下降趋势,并且不同用户数下的SMSE取值最终都趋近于0;从图 6中可以得出,用户数多的算法其收敛速度也是最快的。

Download:
图 5 不同用户对数情况下,本文算法的SMSE性能随SNR变化 Fig. 5 SMSE performance comparisons of the proposed algorithm versus SNR under different user pairs
Download:
图 6 不同用户对数情况下,本文算法的SMSE性能随迭代次数变化 Fig. 6 SMSE performance comparisons of the proposed algorithm versus number of iterations under different user pairs
4 结论

对于包含直传链路的发射/接收端多用户MIMO AF中继系统,提出一种发射端用户集中优化的联合预编码算法。

1) 算法首先将发射单流信号的多用户等效为一个发射多流信号的基站,然后根据等效的通信模型,基于MSMSE准则的算法将原来三变量非凸的优化问题转化为单变量的凸优化问题进行求解。

2) 对于发射端的用户子预编码矩阵Bk(k=1, 2, …, K)的优化问题,通过传统的K.K.T条件求解所有的Bk的优化表达式。

3) 中继收发矩阵F的子优化问题被转化成为一个标准的SDP问题,对于中继收发矩阵F的SDP问题,通过CVX优化工具箱自动识别并求解出局部最优解,最后所有优化矩阵联合迭代至收敛获得最优解。

4) 考虑到现有的文献忽略直传链路对系统性能的贡献,本文将直传链路的作用考虑到发射/接收端多用户MIMO AF中继系统中,可以明显地提高了系统的性能。仿真结果表明,提出算法的SMSE性能和BER性能均优于现有的联合预编码算法,而且算法具有快速的收敛速度。

参考文献
[1]
SHEKHAR, SINGHAL A, SHARMA R, et al. Study of analysis of multiple input and multiple outputs (MIMO) technology in wireless communication[C]//Proceedings of 2016 International Conference on Micro-Electronics and Telecommunication Engineering (ICMETE). Ghaziabad, India: IEEE, 2016: 658-662. (0)
[2]
SETIAWAN D P, ZHAO Huaan. Performance analysis of hybrid AF and DF protocol for relay Networks[C]//Proceedings of 2017 International Conference on Control, Electronics, Renewable Energy and Communications (ICCEREC). Yogyakarta, Indonesia: IEEE, 2017: 207-211. (0)
[3]
KHANDAKER M R A, RONG Yue. Joint transceiver optimization for multiuser MIMO relay communication systems[J]. IEEE transactions on signal processing, 2012, 60(11): 5977-5986. DOI:10.1109/TSP.2012.2212013 (0)
[4]
SHEN Hui, LI Bin, TAO Meixia, et al. MSE-based transceiver designs for the MIMO interference channel[J]. IEEE transactions on wireless communications, 2010, 9(11): 3480-3489. DOI:10.1109/TWC.2010.091510.091836 (0)
[5]
LIU Jun, LIU Zhi, QIU Zhengding. Joint MMSE transceiver design for multiuser non-regenerative MIMO relay downlink systems[C]//Proceedings of the 7th International Wireless Communications and Mobile Computing Conference. Istanbul, Turkey: IEEE, 2011: 877-882. (0)
[6]
SHIN J, MOON J, AHN J. MMSE-based filter design for multi-user peer-to-peer MIMO amplify-and-forward relay systems[C]//Proceedings of 2013 IEEE International Conference on Communications. Budapest, Hungary: IEEE, 2013: 5505-5510. (0)
[7]
HE Zhiqiang, JIANG Weipeng, RONG Yue. Robust design for amplify-and-forward MIMO relay systems with direct link and imperfect channel information[J]. IEEE transactions on wireless communications, 2015, 14(1): 353-363. DOI:10.1109/TWC.2014.2347274 (0)
[8]
HE Zhiqiang, ZHANG Jinnian, LIU Wanning, et al. New results on transceiver design for two-hop amplify-and-forward MIMO relay systems with direct link[J]. IEEE transactions on signal processing, 2016, 64(20): 5232-5241. DOI:10.1109/TSP.2016.2582465 (0)
[9]
WAN Haibin, CHEN Wen. Joint source and relay design for multiuser MIMO nonregenerative relay networks with direct links[J]. IEEE transactions on vehicular technology, 2012, 61(6): 2871-2876. DOI:10.1109/TVT.2012.2196807 (0)
[10]
SUN Qi, LI Lihua. Weighted sum rate maximization for downlink multiuser relay network with direct link[J]. Wireless personal communications, 2014, 75(1): 369-384. DOI:10.1007/s11277-013-1367-6 (0)
[11]
MEHANA A H, NOSRATINIA A. MMSE receive filtering for precoded MIMO systems[C]//Proceedings of 2013 Asilomar Conference on Signals, Systems and Computers. Pacific Grove, CA, USA: IEEE, 2013: 2057-2061. (0)
[12]
BOYD S, VANDENBERGHE L. Convex optimization[M]. Cambridge: Cambridge University Press, 2004. (0)
[13]
ZHANG Xianda. Matrix analysis and applications[M]. Cambridge: Cambridge University Press, 2017: 71-100. (0)
[14]
GRANT M, BOYD S. CVX: Matlab software for disciplined convex programming[EB/OL]. http://stanford.edu/boyd/cvx. (0)
[15]
LOBO M S, VANDENBERGHE L, BOYD S, et al. Applications of second-order cone programming[J]. Linear algebra and its applications, 1998, 284(1/2/3): 193-228. (0)