CN106982085B - 多小区mimo系统中块对角化辅助的鲁棒收发机设计方法 - Google Patents

多小区mimo系统中块对角化辅助的鲁棒收发机设计方法 Download PDF

Info

Publication number
CN106982085B
CN106982085B CN201610028295.2A CN201610028295A CN106982085B CN 106982085 B CN106982085 B CN 106982085B CN 201610028295 A CN201610028295 A CN 201610028295A CN 106982085 B CN106982085 B CN 106982085B
Authority
CN
China
Prior art keywords
receiver
transmitter
kth
interference
channel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610028295.2A
Other languages
English (en)
Other versions
CN106982085A (zh
Inventor
李从改
何晨
蒋铃鸽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201610028295.2A priority Critical patent/CN106982085B/zh
Publication of CN106982085A publication Critical patent/CN106982085A/zh
Application granted granted Critical
Publication of CN106982085B publication Critical patent/CN106982085B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • H04B7/0456Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/06Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0613Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
    • H04B7/0615Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
    • H04B7/0617Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radio Transmission System (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明提供一种多小区MIMO系统中块对角化辅助的鲁棒收发机设计方法,包括步骤1:设置多小区MIMO系统的运行模式;步骤2:定义单个用户或接收机的信干噪比,计算第k个用户或接收机的信干噪比;步骤3:各发射机根据本地的估计信道进行块对角化预编码;步骤4:计算信道存在有界误差时信干噪比的下界值;步骤5:以信干噪比的下界值作为最坏情况的信干噪比,引入松弛变量,以凸优化法转化为求解发射机或接收机的波束成形向量问题;步骤6:通过交替优化方法得到鲁棒的波束成形向量。本发明对有界误差具有较强鲁棒性,每个发射机仅使用本地的CSI进行块对角化预编码,降低了反馈开销,并保证了最小的SINR性能。

Description

多小区MIMO系统中块对角化辅助的鲁棒收发机设计方法
技术领域
本发明涉及的是一种无线通信技术领域的方法,具体是一种多小区MIMO系统中块对角化辅助的鲁棒收发机设计方法。
背景技术
近年来,随着无线通信技术的飞速发展,无线通信大业务量、高速率和高频谱效率的要求日益迫切。在未来无线通信系统中,发射机均采用全频率复用,使得多小区多输入多输出(Multiple-Input Multiple-Output,MIMO)系统中的小区间干扰更为复杂。作为抑制多小区MIMO系统中小区间干扰的关键技术,多个基站协作的收发机联合设计方法成为研究热点之一。根据小区内用户数量,从信息论角度,多小区MIMO系统分为多小区MIMO干扰信道和多小区MIMO干扰广播信道。在多小区MIMO干扰信道中,每个小区内有一个多天线用户;在多小区MIMO干扰广播信道中,每个小区内有多个多天线用户。这里考虑多小区MIMO干扰信道。
通常,联合收发机设计的目标是优化某个特定的系统效用函数,包括和速率最大化、用户公平策略等。
经过对现有技术的检索,L.Li,L.Qiu and G.Wei等人2012年在ElectronicsLetters发表的名称为“Decentralised precoding scheme for MIMO interferencechannel(MIMO干扰信道中分布式的预编码策略)”一文中把和速率最大化作为优化目标。和速率最大化方法把更多的资源用于信道情况好的用户能获得较高的系统容量,但是牺牲了用户的公平性。从用户公平性考虑,一类重要的联合收发机设计是基于SINR的优化方案。Y.-F.Liu,Y.-H.Dai,Z.-Q.Luo等人2013年在IEEE Trans.Signal Process上发表的“Max-min fairness linear tranceiver design problem for a multi-user MIMOinterference channel(多用户MIMO干扰信道中最大-最小公平性线性收发机设计问题)”一文中把基于用户公平性的最大-最小信干噪比优化问题作为优化目标。Y.-F.Liu,M.Hongand Y.-H.Dai等人2013年在IEEE Signal Process.Lett上发表了“Max-min fairnesslinear tranceiver design problem for a multi-user SIMO interference channelis polynomial time solvable(多用户SIMO干扰信道中的最大-最小公平性线性收发机设计问题是多项式时间可解的)”一文中证明在接收波束成形向量固定的单输入多输出干扰信道中,Max-Min SINR预编码设计问题能够在多项式时间(polynomial time)求解。M.Bengtsson,B.Ottersten等2002年IEEE Trans.Wireless Comm上发表了“Jointtransmitter receiver diversity for efficient space division multi-access(高效空分多址中的联合收发机分集)”一文指出联合发送接收波束成形问题可通过交替优化算法求解。R.Mochaourab,P.Cao and E.Jorswieck等人2014年在IEEE SignalProcess.Letters上发表了“Alternating rate profile optimization in singlestream MIMO interference channels(单流MIMO干扰信道中的交替速率优化)”一文中提出在单流MIMO信道中采用交替优化的方式实现max-min SINR优化。
但是,实际系统中,发射机和接收机都很难获得理想的信道状态信息(ChannelState Information,CSI),鲁棒的波束成形算法设计是实际场景中获得有效增益的一个关键因素。在MIMO干扰信道中,SINR信号与干扰加噪声比(Signal to Interference plusNoise Ratio)依赖于发射机和接收机两方面的策略,与MISO信道相比,其耦合性更强,表达式更复杂,因此MISO情况下已有的方案不能直接扩展到MIMO场景中。J.Wang,M.Bengtssonand B.Ottersten等人发表的“Robust MIMO precoding for several classes ofchannel uncertainty(针对几类非理想信道的鲁棒MIMO预编码)”一文中研究了单小区MIMO信道中鲁棒的预编码设计方法。在K对准静态的MIMO干扰信道中,E.Chiu,V.K.N.Lau,H.Huang等人2010年在IEEE Trans.Wireless Commun上发表的“Robust transceiverdesign for K-pairs quasi-static MIMO interference channels via semi-definiterelaxation(K对准静态的MIMO干扰信道中由半正定松弛实现的鲁棒收发机设计)”一文中考虑存在有界误差情况下SINR公平性优化问题,通过半正定松弛提出了鲁棒的波束成形算法;该算法要求CSI误差完全限定在有界的误差的误差界上。R.S.Guiazon,K.K.Wong andD.Wisely等人2014年在IEEE Wireless Commun.Lett发表的“Capacity analysis ofinterference alignment with bounded CSI uncertainty(基于有界误差条件下干扰对齐的容量分析)”一文给出了MIMO干扰信道存在有界误差情况下,采用干扰对齐波束成形算法得到的容量下界。M.J.Rahman,M.Noori and L.Lampe等人2015年在IEEE WCNC上发表的“A low-complexity design for robust SINR fairness in MIMO interferencenetworks”一文中提出了一种基于干扰对齐的鲁棒波束成形算法。众所周知,基于干扰对齐的算法存在一定的限制,如对发送天线的维数要求较高,且仅在高信噪比情况下最优。Q.Zhang,C.He and L.Jiang等人2013年在IEEE GLOBLECOM上发表的“Robust per-streamMSE based transceiver design for MIMO interference channel(MIMO干扰信道中基于每个流MSE的鲁棒收发机设计)”一文中考虑了基于MSE优化的鲁棒波束成形方法。
发明内容
本发明在现有技术基础上,提出一种多小区MIMO系统中块对角化辅助的鲁棒收发机设计方法。
根据本发明提供的多小区MIMO系统中块对角化辅助的鲁棒收发机设计方法,包括以下步骤:
步骤1:设置多小区MIMO系统的运行模式;
步骤2:定义单个用户或接收机的信干噪比,计算第k个用户或接收机的信干噪比,k为多小区内任一用户或接收机的序号;
步骤3:各发射机根据本地的估计信道进行块对角化预编码,消除部分对非同一小区接收机的干扰;
步骤4:计算信道存在有界误差时信干噪比的下界值;
步骤5:以信干噪比的下界值作为最坏情况的信干噪比,引入松弛变量,以凸优化法转化为求解发射机或接收机的波束成形向量问题;
步骤6:通过交替优化方法得到鲁棒的波束成形向量,所述交替优化方法是指:首先固定发送波束成形向量,通过求解SOCP问题获得接收波束成形向量;然后固定接收波束成形向量,求解SOCP问题得到新的接收波束成形向量。
优选地,所述步骤1包括:在多小区MIMO系统中,每个小区有一个发射机并服务一个接收机,每个发射机的天线数M,K个发射机采用全频率复用进行协作,每个接收机有N个接收天线,第k个发射机的功率约束Pk,第k个发射机发送单流数据给第k个接收机,第k个接收机接收到的加性高斯白噪声向量服从零均值、协方差为
Figure GDA0002522585600000032
I的复高斯分布,其中:k=1,...,K。
优选地,所述步骤2包括:
第k个接收机的信干噪比定义为
Figure GDA0002522585600000031
式中,Hki为第i个发射机到第k个接收机的实际信道状态信息,其中:k,i=1,...,K,k=1,...,K,fk为第k个发射机的发送波束成形向量,fi表示第i个发射机的发送波束成形向量,uk为第k个接收机的接收波束成形向量,
Figure GDA00025225856000000415
表示uk的共轭转置,σk表示第k个接收机接收到的加性高斯白噪声的协方差,其中:k=1,...,K;
优选地,所述步骤3包括:
令||uk||2=1,定义信道存在有界误差时第k个接收机的信干噪比为:
Figure GDA0002522585600000041
式中,
Figure GDA0002522585600000042
表示信道存在有界误差时第k个用户的信干噪比,
Figure GDA0002522585600000043
表示第i个发送机到第k个接收机的估计信道矩阵,Δki表示第i个发送机到第k个接收机有界误差信道矩阵,
Figure GDA0002522585600000044
表示向量或矩阵的集合;
各发射机根据本地的估计信道对发送数据进行块对角化预编码,简化SINR效用函数;具体地,
所述的块对角化预编码是指:一种迫零算法,其中,迫零限制为:
Figure GDA0002522585600000045
定义
Figure GDA0002522585600000046
其中
Figure GDA0002522585600000047
L=N(K-1);
式中:Gk表示第k个接收机对其他小区接收机的干扰信道,
Figure GDA0002522585600000048
表示第k个发射机到第K个接收机的估计信道,
Figure GDA0002522585600000049
表示L行M列的复数矩阵,L表示干扰信道的行数;
令fk=Qkvk,其中Qk是Gk零空间的正交基,vk表示待求的预编码向量;对Gk进行奇异值分解获得零奇异值对应的右奇异向量Qk
Figure GDA00025225856000000410
其中,
Figure GDA00025225856000000411
Figure GDA00025225856000000412
Figure GDA00025225856000000413
式中:Uk表示对Gk进行SVD分解获得的左特征向量,
Figure GDA00025225856000000414
表示对Gk进行SVD分解获得的右特征向量,
Figure GDA0002522585600000051
表示Uk的共轭转置,I表示单位矩阵,
Figure GDA0002522585600000052
表示M行M列的复数矩阵,Vk表示
Figure GDA0002522585600000053
的共轭转置,Vk(:,L+1:M)表示矩阵Vk第L+1到第M列的全部元素;
Figure GDA0002522585600000054
表示M行M-L列的复数矩阵;
简化的SINR效用函数是指:
Figure GDA0002522585600000055
式中:
Figure GDA0002522585600000056
表示第k个发射机到第k个接收机的估计信道,Δkk表示第k个发射机到第k个接收机的误差信道,vk表示待求的预编码向量。
优选地,所述步骤4包括:
定义信道存在有界误差时信干噪比的下界为
Figure GDA0002522585600000057
根据内积的性质,利用柯西不等式分析代价函数ζk,把误差参量形式的SINR下界表达式转化成基于误差界的显式表达式;
Figure GDA0002522585600000058
其中:
Figure GDA0002522585600000059
式中:
Figure GDA00025225856000000510
表示
Figure GDA00025225856000000511
的第j个行向量,
Figure GDA00025225856000000512
表示归一化的
Figure GDA00025225856000000513
(·)+表示非负数运算,εki表示误差界,
Figure GDA00025225856000000514
表示Frobenius范数。
优选地,所述步骤5包括:多小区MIMO系统中最坏情况的Max-Min SINR优化问题是指:
Figure GDA00025225856000000515
subject to||uk||2=1,||fk||2≤Pk,||Δki||F≤εki
转化后的优化问题是指:
Figure GDA0002522585600000061
式中:a、tk、bki为引入的松弛变量。
优选地,所述步骤6中的SOCP问题是指:
Figure GDA0002522585600000062
Figure GDA0002522585600000063
其中:a、tk、bki为引入的松弛变量。
与现有技术相比,本发明具有如下的有益效果:
本发明提供的多小区MIMO系统中块对角化辅助的鲁棒收发机设计方法考虑每个发射机的功率约束,优化信道信息存在有界误差情况下多小区MIMO系统中最小的SINR,各发射机根据本地的估计信道进行块对角化预编码,消除部分对其他小区接收机的干扰,简化SINR效用函数;根据内积的性质和柯西-许瓦尔兹不等式,求出最坏情况下SINR的下界;引入松弛变量,利用凸优化理论和交替优化方法,把非凸的优化问题转化成SOCP问题计算收发端的波束成形向量,对有界误差具有较强鲁棒性,每个发射机仅使用本地的CSI进行块对角化预编码,降低了反馈开销,并保证了最小的SINR性能。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为多小区MIMO系统的示意图;
图2为场景K=3,M=6,N=2时(K表示发射机数目,M表示每个发射机的天线数,N表示每个接收机的天线数),不同误差界条件下分别采用本发明提供的方法和现有技术中基于干扰对齐的鲁棒算法(Lower bound)的收敛性比较图。
图3为场景K=3,M=6,N=2时,第一种误差界条件下分别采用本发明提供的方法和现有技术中的基于干扰对齐的鲁棒算法(Lower bound)以及Min-Max MSE的鲁棒算法(MSE-based)的最小SINR性能比较图。
图4为场景K=3,M=6,N=2时,第二种误差界条件下分别采用本发明提供的方法和现有技术中的基于干扰对齐的鲁棒算法(Lower bound)以及Min-Max MSE的鲁棒算法(MSE-based)的最小SINR性能比较图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
根据本发明提供的多小区MIMO系统中块对角化辅助的鲁棒收发机设计方法,在有界误差情况下考虑每个发射机的功率约束,以多小区MIMO系统最小的SINR最大化为优化目标,利用块对角化预编码简化SINR效用函数,引入松弛变量,把简化后的Max-Min优化问题转化二阶锥规划(Second-Order Cone Program,SOCP)形式,利用交替优化方法分别获得鲁棒的发送和接收波束成形方法。
具体地,包括以下步骤:
步骤S1:设置系统参数:多小区MIMO系统,每个小区有一个发射机,每个发射机的天线数M,K个发射机采用全频率服用进行协作,每个接收机有N个接收天线,下行通信表示为
Figure GDA0002522585600000071
第k个发射机的功率约束Pk,第k个发射机发送单流数据给第k个接收机,第k个接收机接收到的向量服从零均值、协方差为
Figure GDA0002522585600000072
的复高斯分布,其中:k=1,...,K,第i个发送机到第k个接收机的估计信道矩阵
Figure GDA0002522585600000073
有界误差信道矩阵Δki,其中:k,i=1,...,K;
步骤S2:定义SINRk为第k个接收机的SINR,
Figure GDA0002522585600000081
其中:Hki为第i个发射机到第k个接收机的实际信道状态信息,其中:k,i=1,...,K,k=1,...,K,fk为第k个发射机的发送波束成形向量,uk为第k个接收机的接收波束成形向量,其中:k=1,...,K;
步骤S3:令||uk||2=1,定义信道存在有界误差时第k个接收机的信干噪比为
Figure GDA0002522585600000082
各发射机根据本地的估计信道矩阵
Figure GDA0002522585600000083
进行块对角化预编码,消除部分对其他小区接收机的干扰,简化SINR效用函数;
所述的块对角化预编码是指:
块对角化预编码是一种迫零算法,其迫零限制为
Figure GDA0002522585600000084
定义
Figure GDA0002522585600000085
其中
Figure GDA0002522585600000086
L=N(K-1)。令fk=Qkvk,其中Qk是Gk零空间的正交基。对Gk进行奇异值分解(Singular Value Decomposition,SVD)可以获得Qk,即
Figure GDA0002522585600000087
其中,
Figure GDA0002522585600000088
Figure GDA0002522585600000089
Figure GDA00025225856000000810
是零奇异值对应的右奇异向量。
简化的SINR效用函数是指:
Figure GDA00025225856000000811
其中:发送波束成形向量fk=Qkvk,Qk是块对角化预编码矩阵,vk是所需的发送波束成形向量,Σk是对Gk进行SVD分解获得的半正定L*M阶对角矩阵。
步骤S4:定义信道存在有界误差时信干噪比的下界为
Figure GDA00025225856000000812
根据内积的性质,利用柯西-许瓦尔兹不等式分析代价函数ζk,把误差参量形式的SINR下界表达式转化成基于误差界的显式表达式;
所述的内积的性质是指:
两个向量
Figure GDA0002522585600000091
Figure GDA0002522585600000092
的内积与迹存在如下关系:aHb=tr(baH)。
柯西-许瓦尔兹不等式(Cauchy-Schwarz inequality)是指:
对于
Figure GDA0002522585600000093
Figure GDA0002522585600000094
其内积满足如下性质
Figure GDA0002522585600000095
等式成立的条件是当且仅当x和y线性无关。
基于误差界的显式表达式是指:
Figure GDA0002522585600000096
其中:
Figure GDA0002522585600000097
Figure GDA0002522585600000098
Figure GDA0002522585600000099
的第j个行向量,j=1,...,N。
步骤S5:信道信息存在有界误差情况下,用代价函数ζk作为最坏情况的SINR,引入松弛变量,根据凸优化知识,对多小区MIMO系统中最坏情况的Max-Min SINR优化问题进行转化;
所述的最坏情况的Max-Min SINR优化问题是指:
Figure GDA00025225856000000910
subject to||uk||2=1,||fk||2≤Pk,||Δki||F≤εki
转化后的优化问题是指:
Figure GDA00025225856000000911
步骤S6:如果固定发送或接收的波束成形向量,转化后的优化问题变成二阶锥规划(second-order cone program,SOCP)问题。采用交替优化方法,首先固定发送波束成形向量,通过求解SOCP问题获得接收波束成形向量;然后固定接收波束成形向量,求解SOCP问题得到新的接收波束成形向量。由此,通过迭代得到了一种鲁棒的收发机设计方法。
所述的固定发送波束成形向量获得的SOCP问题是指:
Figure GDA0002522585600000101
固定接收波束成形向量获得的SOCP问题是指:
Figure GDA0002522585600000102
其中:a、tk、bki为引入的松弛变量。
更进一步地,根据内积的性质和柯西-许瓦尔兹不等式,求出最坏情况下SINR的下界;引入松弛变量,利用凸优化理论和交替优化方法,把非凸的优化问题转化成SOCP问题计算收发端的波束成形向量,对有界误差具有较强鲁棒性,所述信道信息存在有界误差的多小区MIMO系统的Max-Min SINR优化问题是:
Figure GDA0002522585600000103
subject to||uk||2=1,||fk||2≤Pk,||Δki||F≤εki
其中:
Figure GDA0002522585600000104
Figure GDA0002522585600000105
为第i个发射机到第k个接收机的估计信道矩阵,Δki是信道误差矩阵且||Δki||F≤εki,fk为第k个发射机的发送波束成形向量,uk为第k个接收机的接收波束成形向量,Pk是第k个发射机的功率约束,
Figure GDA0002522585600000106
是第k个接收机接收到的复高斯白噪声的方差,k=1,...,K。
具体地,本发明提供的实施例包括以下步骤:
第一步:根据步骤S1设置系统参数;本实施例中,所用的仿真场景为K=3,M=6,N=2;
Figure GDA0002522585600000111
是服从独立同分布的复高斯分布的随机矩阵,||Δki||F≤εki
Figure GDA0002522585600000112
其中:k,i=1,...,K;
且信噪比
Figure GDA0002522585600000113
P=1,
Figure GDA0002522585600000114
第二步:定义SINRk为第k个接收机的SINR,
Figure GDA0002522585600000115
式中:Hki为第i个发射机到第k个接收机的实际信道状态信息,其中:k,i=1,...,K,k=1,...,K,fk为第k个发射机的发送波束成形向量,uk为第k个接收机的接收波束成形向量,其中:k=1,...,K;
第三步:令||uk||2=1,定义信道存在有界误差时第k个接收机的信干噪比为
Figure GDA0002522585600000116
各发射机根据本地的估计信道矩阵
Figure GDA0002522585600000117
进行块对角化预编码,消除部分对其他小区接收机的干扰,简化SINR效用函数,
Figure GDA0002522585600000118
其中:发送波束成形向量fk=Qkvk,Qk是块对角化预编码矩阵,vk是所需的发送波束成形向量;定义
Figure GDA0002522585600000119
其中
Figure GDA00025225856000001110
L=N(K-1)。令fk=Qkvk,其中Qk是Gk零空间的正交基。对Gk进行奇异值分解可以获得Qk,即
Figure GDA00025225856000001111
其中,
Figure GDA00025225856000001112
Figure GDA00025225856000001113
Figure GDA00025225856000001114
Figure GDA00025225856000001115
Figure GDA00025225856000001116
是零奇异值对应的右奇异向量。
第四步:定义信道存在有界误差时信干噪比的下界为
Figure GDA00025225856000001117
根据内积的性质,利用柯西-许瓦尔兹不等式分析代价函数ζk,把误差参量形式的SINR下界表达式转化成基于误差界的显式表达式,
Figure GDA0002522585600000121
其中:(x)+=max{0,x},
Figure GDA0002522585600000122
Figure GDA0002522585600000123
Figure GDA0002522585600000124
的第j个行向量,j=1,...,N。
第五步:信道信息存在有界误差情况下,用代价函数ζk作为最坏情况的SINR,引入松弛变量,根据凸优化知识,对多小区MIMO系统中最坏情况的Max-Min SINR优化问题进行转化,
Figure GDA0002522585600000125
第六步:采用交替优化方法,首先固定发送波束成形向量,通过求解SOCP问题获得接收波束成形向量;然后固定接收波束成形向量,求解SOCP问题得到新的接收波束成形向量。由此,通过迭代得到了一种鲁棒的收发机设计方法。
所述的固定发送波束成形向量获得的SOCP问题是指:
Figure GDA0002522585600000126
固定接收波束成形向量获得的SOCP问题是指:
Figure GDA0002522585600000127
其中:a、tk、bki为引入的松弛变量。在优化过程中,a>0采用二分法进行更新,选定最小值amin和最大值amax,则最优的SINR表示为aopt∈[amin,amax]。令
Figure GDA0002522585600000131
Figure GDA0002522585600000132
Figure GDA0002522585600000133
时算法停止迭代,其中δ是给定的终止阈值。amin和amax的初始值选择:
amin=0,
Figure GDA0002522585600000134
具体地,如图2所示为场景K=3,M=6,N=2时,不同误差界条件下分别采用本发明提供的方法和现有技术中基于干扰对齐的鲁棒算法(Lower bound)的收敛性比较图。图2展示了SNR=15dB时提出的BDA算法的收敛性。从图中可以看出,最小的SINR值如预期单调增加,并且大部分性能增长是在前几次迭代中实现的。在非理想信道条件下,BDA算法的最小SINR性能优于下界。因此,仿真结果证实了提出算法的鲁棒性。
具体地,如图3和图4所示为场景K=3,M=6,N=2时,两种误差界条件下分别采用本发明提供的方法和现有技术中的基于干扰对齐的鲁棒算法(Lower bound)以及Min-MaxMSE的鲁棒算法(MSE-based)的最小SINR性能比较图。由图可知,BDA算法在中高信噪比范围内有更好的最小SINR性能,而MSE-based算法更适用于低信噪比的情况。在多小区MIMO系统中,最小SINR等价于最小速率,因此BDA算法具有可行性和鲁棒性。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。

Claims (1)

1.一种多小区MIMO系统中块对角化辅助的鲁棒收发机设计方法,其特征在于,包括以下步骤:
步骤1:设置多小区MIMO系统的运行模式;
步骤2:定义单个用户或接收机的信干噪比,计算第k个用户或接收机的信干噪比,k为多小区内任一用户或接收机的序号;
步骤3:各发射机根据本地的估计信道进行块对角化预编码,消除部分对非同一小区接收机的干扰;
步骤4:计算信道存在有界误差时信干噪比的下界值;
步骤5:以信干噪比的下界值作为最坏情况的信干噪比,引入松弛变量,以凸优化法转化为求解发射机或接收机的波束成形向量问题;
步骤6:通过交替优化方法得到鲁棒的波束成形向量,所述交替优化方法是指:首先固定发送波束成形向量,通过求解SOCP问题获得接收波束成形向量;然后固定接收波束成形向量,求解SOCP问题得到新的接收波束成形向量;
所述步骤1包括:在多小区MIMO系统中,每个小区有一个发射机并服务一个接收机,每个发射机的天线数M,K个发射机采用全频率复用进行协作,每个接收机有N个接收天线,第k个发射机的功率约束Pk,第k个发射机发送单流数据给第k个接收机,第k个接收机接收到的加性高斯白噪声向量服从零均值、协方差为
Figure FDA0002522585590000011
I的复高斯分布,其中:k=1,...,K;
所述步骤2包括:
第k个接收机的信干噪比定义为
Figure FDA0002522585590000012
式中,Hki为第i个发射机到第k个接收机的实际信道状态信息,其中:k,i=1,...,K,k=1,...,K,fk为第k个发射机的发送波束成形向量,fi表示第i个发射机的发送波束成形向量,uk为第k个接收机的接收波束成形向量,
Figure FDA0002522585590000013
表示uk的共轭转置,σk表示第k个接收机接收到的加性高斯白噪声的协方差,其中:k=1,...,K;
Figure FDA0002522585590000021
表示向量或矩阵的集合;
所述步骤3包括:
令||uk||2=1,定义信道存在有界误差时第k个接收机的信干噪比为:
Figure FDA0002522585590000022
式中,
Figure FDA0002522585590000023
表示信道存在有界误差时第k个用户的信干噪比,
Figure FDA0002522585590000024
表示第i个发送机到第k个接收机的估计信道矩阵,Δki表示第i个发送机到第k个接收机有界误差信道矩阵,
Figure FDA0002522585590000025
表示向量或矩阵的集合;
各发射机根据本地的估计信道对发送数据进行块对角化预编码,简化SINR效用函数;具体地,
所述的块对角化预编码是指:一种迫零算法,其中,迫零限制为:
Figure FDA0002522585590000026
定义
Figure FDA0002522585590000027
其中
Figure FDA0002522585590000028
L=N(K-1);
式中:Gk表示第k个接收机对其他小区接收机的干扰信道,
Figure FDA0002522585590000029
表示第k个发射机到第K个接收机的估计信道,
Figure FDA00025225855900000210
表示L行M列的复数矩阵,L表示干扰信道的行数;
令fk=Qkvk,其中Qk是Gk零空间的正交基,vk表示待求的预编码向量;对Gk进行奇异值分解获得零奇异值对应的右奇异向量Qk
Figure FDA00025225855900000211
其中,
Figure FDA00025225855900000212
Figure FDA00025225855900000213
Figure FDA00025225855900000214
式中:Uk表示对Gk进行SVD分解获得的左特征向量,
Figure FDA00025225855900000215
表示对Gk进行SVD分解获得的右特征向量,
Figure FDA00025225855900000216
表示Uk的共轭转置,I表示单位矩阵,
Figure FDA00025225855900000217
表示M行M列的复数矩阵,Vk表示
Figure FDA0002522585590000031
的共轭转置,Vk(:,L+1:M)表示矩阵Vk第L+1到第M列的全部元素;
Figure FDA0002522585590000032
表示M行M-L列的复数矩阵;
简化的SINR效用函数是指:
Figure FDA0002522585590000033
式中:
Figure FDA0002522585590000034
表示第k个发射机到第k个接收机的估计信道,Δkk表示第k个发射机到第k个接收机的误差信道,vk表示待求的预编码向量;
Σk是对Gk进行SVD分解获得的半正定L*M阶对角矩阵;
所述步骤4包括:
定义信道存在有界误差时信干噪比的下界为
Figure FDA0002522585590000035
根据内积的性质,利用柯西不等式分析代价函数ζk,把误差参量形式的SINR下界表达式转化成基于误差界的显式表达式;
Figure FDA0002522585590000036
其中:(x)+=max{0,x},
Figure FDA0002522585590000037
式中:
Figure FDA0002522585590000038
表示
Figure FDA0002522585590000039
的第j个行向量,
Figure FDA00025225855900000310
表示归一化的
Figure FDA00025225855900000311
(·)+表示非负数运算,εki表示误差界,
Figure FDA00025225855900000312
表示Frobenius范数;
所述步骤5包括:多小区MIMO系统中最坏情况的Max-Min SINR优化问题是指:
Figure FDA00025225855900000313
subject to ||uk||2=1,||fk||2≤Pk,||Δki||F≤εki
转化后的优化问题是指:
Figure FDA0002522585590000041
式中:a、tk、bki为引入的松弛变量;
所述步骤6中的SOCP问题是指:
Figure FDA0002522585590000042
a
subject to
Figure FDA0002522585590000043
Figure FDA0002522585590000044
||uk||2=1
Figure FDA0002522585590000045
a
subject to
Figure FDA0002522585590000046
Figure FDA0002522585590000047
||fk||2≤Pk
其中:a、tk、bki为引入的松弛变量。
CN201610028295.2A 2016-01-15 2016-01-15 多小区mimo系统中块对角化辅助的鲁棒收发机设计方法 Expired - Fee Related CN106982085B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610028295.2A CN106982085B (zh) 2016-01-15 2016-01-15 多小区mimo系统中块对角化辅助的鲁棒收发机设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610028295.2A CN106982085B (zh) 2016-01-15 2016-01-15 多小区mimo系统中块对角化辅助的鲁棒收发机设计方法

Publications (2)

Publication Number Publication Date
CN106982085A CN106982085A (zh) 2017-07-25
CN106982085B true CN106982085B (zh) 2020-08-04

Family

ID=59340502

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610028295.2A Expired - Fee Related CN106982085B (zh) 2016-01-15 2016-01-15 多小区mimo系统中块对角化辅助的鲁棒收发机设计方法

Country Status (1)

Country Link
CN (1) CN106982085B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110048787A (zh) * 2019-04-04 2019-07-23 中国人民解放军63686部队 一种应用于矢量信号分析的测量滤波器选择方法
CN113411112B (zh) * 2021-06-09 2023-03-28 西安交通大学 一种分布式鲁棒多小区协作波束成形admm网络的构建方法
CN113949607B (zh) * 2021-10-22 2024-01-12 东南大学 一种智能反射面无小区系统的鲁棒性波束设计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103427889A (zh) * 2012-05-22 2013-12-04 北京邮电大学 预编码模式选择方法及其系统
CN103458521A (zh) * 2013-09-09 2013-12-18 东南大学 一种基于鲁棒性设计的mimo传输功率分配优化方法
CN103490804A (zh) * 2013-09-12 2014-01-01 江苏科技大学 基于优先度遗传模拟退火的多用户mimo系统天线选择方法
US8625560B2 (en) * 2010-12-02 2014-01-07 Samsung Electronics Co., Ltd. Method and apparatus for feeding back channel quality information in multi-user multi-input multi-output communication system
CN103684560A (zh) * 2013-12-04 2014-03-26 上海交通大学 多小区多用户系统中基于用户公平性的鲁棒预编码方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8625560B2 (en) * 2010-12-02 2014-01-07 Samsung Electronics Co., Ltd. Method and apparatus for feeding back channel quality information in multi-user multi-input multi-output communication system
CN103427889A (zh) * 2012-05-22 2013-12-04 北京邮电大学 预编码模式选择方法及其系统
CN103458521A (zh) * 2013-09-09 2013-12-18 东南大学 一种基于鲁棒性设计的mimo传输功率分配优化方法
CN103490804A (zh) * 2013-09-12 2014-01-01 江苏科技大学 基于优先度遗传模拟退火的多用户mimo系统天线选择方法
CN103684560A (zh) * 2013-12-04 2014-03-26 上海交通大学 多小区多用户系统中基于用户公平性的鲁棒预编码方法

Also Published As

Publication number Publication date
CN106982085A (zh) 2017-07-25

Similar Documents

Publication Publication Date Title
US9236922B2 (en) Method of communicating between base station and terminal based on interference alignment in multi-cell multi-user multiple-input multiple-output (MIMO) interference channel and method and apparatus of communication using interference alignment and block successive interference pre-cancellation in multi-user multiple-input multiple-output interference channel
US9379800B2 (en) Method and apparatus for channel state information feedback
US8472550B2 (en) Method for reducing interference in multi-cell multi-user wireless networks
KR101312015B1 (ko) 무선 통신 시스템 내의 무선 디바이스에 전송될 프리코딩 데이터에 대한 프리코딩 벡터를 결정하는 방법 및 장치
US20110176633A1 (en) Method and system for orthogonalized beamforming in multiple user multiple input multiple output (mu-mimo) communication systems
US8649457B2 (en) Precoding process for a transmitter of a MU-MIMO communication system
CN105049100B (zh) 一种多小区mimo系统双层预编码方法
CN102882570B (zh) 移动通信网络下设备间通信的最优收发联合处理方法
CN101867462A (zh) 一种基于最小总误码率的多基站协作线性预编码方法
CN105745893A (zh) 无线大量mimo系统中大尺度衰落系数估计
CN101242381B (zh) 多用户多输入多输出系统的线性预编码方法
CN106982085B (zh) 多小区mimo系统中块对角化辅助的鲁棒收发机设计方法
CN109818662A (zh) 全双工云接入数能一体化网络中混合波束成形方法
CN106102073A (zh) 一种异构认知无线网络下行链路干扰对齐方法
Wei et al. TDD reciprocity calibration for multi-user large-scale (massive) MIMO systems with iterative coordinate descent
CN102291727B (zh) 一种分布式协作波束成形设计方法
CN110943767A (zh) Fdd大规模mimo系统中基于信道部分互易性的预编码设计方法
CN102740325B (zh) 获得信道信息的方法、装置及波束成形优化方法、装置
US20170230203A1 (en) Method and apparatus of topological pilot decontamination for massive mimo systems
CN106850020A (zh) 非理想信道状态下mimo系统中基于组合干扰对齐方法
CN111510403B (zh) 一种基于对称lq的大规模mimo系统迭代信号检测方法
CN102291210A (zh) 一种生成预编码矩阵的方法及装置
US9048970B1 (en) Feedback for cooperative multipoint transmission systems
CN104618006B (zh) 一种分布式多用户下行mimo鲁棒波束成形方法
CN105162500B (zh) 一种多小区mimo系统协作干扰抑制方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200804

Termination date: 20220115