具体实施方式
参阅图1,本发明实施例提供一种预编码矩阵获取方法,该方法包括:
101、针对每个接收端,从所有干扰向量中选取作为干扰空间基向量的干扰向量,将剩余干扰向量表示为所述干扰空间基向量的线性组合;其中,所述剩余干扰向量为接收端的所有干扰向量中除去作为干扰空间基向量的干扰向量以外的干扰向量。
102、根据每个接收端剩余干扰向量所表示成的干扰空间基向量的线性组合及当前的信道状况,确定当前采用的预编码矩阵。
其中,所述干扰向量可以表示为信道向量与预编码向量的乘积,此时针对第k个接收端,步骤101中确定干扰空间基向量的个数为N-LK个,从L-LK个干扰向量中选择N-LK个干扰向量作为干扰空间基向量;N为第k个接收端的接收天线数,LK为第k个接收端的信号流数;L为多个接收端的信号流总数。具体的,因为第k个接收端的接收天线数为N,而第k个接收端的信号流数为LK,所以对于第k个接收端来说,其具有LK维的信号空间,N-LK维的干扰空间;由于多个接收端的信号流总数L,而第k个接收端的信号流数为LK,可知针对该第k个接收端,干扰向量的总数为L-LK个,从该L-LK个中选择N-LK个作为干扰空间基向量,那么剩余干扰向量(即需要对齐的干扰向量)为(L-LK)-(N-LK)=L-N个。
这种情况下,该步骤102的具体实现过程包括:
根据每个接收端剩余干扰向量所表示成的干扰空间基向量的线性组合,得到矩阵
根据所述矩阵和当前的信道状况,利用最优化方法,得到第一拉格朗日函数;
其中,
n=(k-1)(L-N)+n’,n’=1,2,...,L-N,α
n,0=-1,α
n,1、α
n,2、α
n,3、α
n(N-Lk)分别为第n个向量方程中第k个接收端的第n’个剩余干扰向量表示为干扰空间基向量线性组合时干扰空间基向量的系数;
中的第一列为一剩余干扰向量,其余列为干扰空间基向量;
根据每个接收端剩余干扰向量所表示成的干扰空间基向量的线性组合,得到Aw=0;根据所述Aw=0和当前的信道状况,利用最优化方法,得到第二拉格朗日函数;其中,A中的每一个子矩阵Anl的维度均为N×M,l≤L,L为数据流总数,第n个方程中的剩余干扰向量对应的Anl为αn,0与所述剩余干扰向量中的信道矩阵的乘积,第n个方程中的干扰空间基向量对应的Anl为所述干扰空间基向量中信道矩阵与所述干扰空间基向量的系数的乘积;
利用第一拉格朗日函数和第二拉格朗日函数,得到向量an和预编码矩阵w的互相关函数;其中,第一拉格朗日函数和第二拉格朗日函数的具体形式见后续实施例的详细描述。
将预设的w初始值或者向量an的初始值代入所述互相关函数,得到w当前值或者向量an的当前值,将所述w当前值或者向量an的当前值再代入所述互相关函数,直到满足预设条件为止,所述预设条件为将得到的w当前值或者向量an的当前值本次代入第一或者第二拉格朗日函数得到的值与前一次代入时得到的值的差达到第一预设门限,或者,迭代次数达到第二门限。
或者,干扰向量可以表示为信道向量与均衡向量的乘积,此时针对第k个接收端,确定干扰空间基向量的个数为M-LK个,从L-LK个干扰向量中选择M-LK个干扰向量作为干扰空间基向量;其中,M为与第k个接收端配对的发射端的发射天线数,LK为第k个接收端的信号流数;L为多个接收端的信号流总数。具体的,因为与第k个接收端配对的发射端的发射天线数为M,而第k个接收端的信号流数为LK,所以对于该发射端来说,其具有LK维的信号空间,M-LK维的干扰空间;根据多个接收端的信号流总数L,而第k个接收端的信号流数为LK,可知针对该发送端,干扰向量的总数为L-LK个,从该L-LK个中选择M-LK个作为干扰空间基向量,那么剩余干扰向量(即需要对齐的干扰向量)为(L-LK)-(M-LK)=L-M个。
这种情况下,该步骤102的具体实现过程包括:
根据每个接收端剩余干扰向量所表示成的干扰空间基向量的线性组合,得到矩阵
根据所述矩阵和当前的信道状况,利用最优化方法,得到第三拉格朗日函数;
其中,
n=(k-1)(L-M)+n’,n’=1,2,...,L-N,β
n,0=-1,β
n,1、β
n,2、β
n,3、β
n(M-Lk)分别为第n个向量方程中第k个接收端的第n’个剩余干扰向量表示为干扰空间基向量线性组合时干扰空间基向量的系数;
中的第一列为一剩余干扰向量,其余列为干扰空间基向量;
根据每个接收端剩余干扰向量所表示成的干扰空间基向量的线性组合,得到BG=0;根据所述Bw=0和当前的信道状况,利用最优化方法,得到第四拉格朗日函数;其中,B中的每一个子矩阵Bnl的维度均为M×N,l≤L,L为数据流总数,第n个方程中的剩余干扰向量对应的Bnl为βn,0与所述剩余干扰向量中的信道矩阵的乘积,第n个方程中的干扰空间基向量对应的Bnl为所述干扰空间基向量中信道矩阵与所述干扰空间基向量的系数的乘积;其中,G为均衡矩阵,是由均衡向量组成的矩阵。
利用第三拉格朗日函数和第四拉格朗日函数,得到向量an和均衡矩阵的互相关函数;
将预设的均衡矩阵G的初始值或者向量βn的初始值代入所述互相关函数,得到均衡矩阵G的当前值或者向量βn的当前值,将所述均衡矩阵G的当前值或者向量βn的当前值再代入所述互相关函数,直到满足预设收敛条件为止,所述预设收敛条件为将得到的均衡矩阵G的当前值或者向量βn的当前值本次代入第三或者第四拉格朗日函数得到的值与前一次代入时得到的值的差达到第三预设门限,或者,迭代次数达到第四门限。
其中,本发明实施例中的发射端可以是基站,接收端可以是用户设备。
本发明实施例中,针对每个接收端,都将需要对齐的干扰向量表示为干扰空间基向量的线性组合,这样达到干扰完全对齐,能够完全消除干扰,不会产生BER平层。
一、为了说明步骤101中将需要对齐的干扰向量表示为干扰空间基向量的线性组合的思想,如下先以每个用户设备只接收一个信号流为例进行描述:
第k个接收端接收的信号可以表示为:
其中,H
kkW
kx
k为接收的有用信号;
为接收的干扰信号,n
k是方差为σ
2的噪声向量。这里,H
kk表示第k个基站到第k个用户设备的信道矩阵,H
ki表示第i个基站到第k个用户设备的信道矩阵,维数为N×M;
Wk为第k个基站的预编码矩阵,xk为第k个基站发射的信号;Wi为第i个基站的预编码矩阵,xi为第i个基站发射的信号。
其中,公式(1)中:
Li为第i个用户设备的预编码向量数,也即第i个用户设备传输的信号流的数目,所有用户设备(假定用户设备一种有K个)共有个信号流,不难理解,这L个预编码向量之间必须两两独立,wim为第i个基站的第m个预编码向量,维数为M×1,xim表示发送的信号,为一标量。
首先,以Li=1的情况来介绍本发明提供的干扰对齐方案,即每个用户设备都只有一个流,则(1)式变为:
yk=Hk1w1x1+Hk2w2x2+...+Hkkwkxk+...+nk
则根据第k个用户设备的接收信号所估计出的发射信号为:
其中,
为一个维度为1×N的均衡向量,用于加权。当||g
k||
2=1时,g
k的选择并不改变噪声的功率,则获得x
k的最佳估计的一个必要条件是:对所有i,k=1,2,…,K,有:
ck为正数,
其中为了确保能最大化SINR,ck需最大化,以达到最优的检测性能。
对于第k个用户设备,当w
i确定时,要满足
i=1,2,…,K,可令a
i=H
kiw
i,则有
其中a
i=(a
i1,a
i2,…,a
iN)
T,可以得到如下含K-1个方程的线性方程组:
为了保证公式(4)有非零解,则必须使得上述K-1个方程中,至少有L-N个干扰迫零方程的系数向量是N-1个干扰迫零方程的系数向量的线性组合,即为干扰对齐的约束条件,其中,干扰迫零方程是指(4)中右边为零的方程。对于一个用户设备只接收一个信号流的情况,则:
其中,k=1,2,…,K,in|n=1,2,...,N-1∈Sk1,im|m=1,2,...,K-N∈Sk2,αkn表示组合系数,Sk1表示干扰空间基向量的集合,该集合中包含N-1个元素,Sk2表示包含需要对齐的干扰向量的集合,该集合包括K-N个元素,且
Sk1∪Sk2∪{k}={1,2,…,K}。
上述方案为了描述简便,以单流为例,即每个用户设备只接收一路信号,也即每个用户设备的信号流数为1。若第k个用户设备的信号流数为LK,则针对该用户设备,干扰向量总个数为L-LK,干扰空间的基向量的个数为N-LK,则在所有干扰向量中选择N-LK个作为干扰空间的基向量,然后将剩余干扰向量表示为干扰空间的基向量的线性组合,其中,剩余干扰向量即是需要对齐的干扰向量,其为(L-LK)-(N-LK)=L-N个。其中,N为第k个用户设备的接收天线数,L为所有用户设备的信号流总数。
为了使本发明实施例提供的干扰对齐方案更加清楚明白,如下举实例进行说明。图2为本发明实施例提供的干扰对齐示意图,假定M=6,N=3,K=5,即有5个基站,5个用户设备,每个用户设备的接收天线数都为3,每个用户设备只接收一个信号流。则对每个用户设备来说,因接收天线数为3,则接收的信号向量呈3维,除去有用信号占1维之外,还剩2维,则要想实现干扰对齐,必须使4个干扰向量中的任意3个干扰向量线性相关,这样,就需要从4个干扰向量中任意挑出2个作为干扰空间的基向量,将剩下的两个干扰向量(即需要对齐的干扰向量)表示为干扰空间的基向量的线性组合。
参见图2,对于第1个用户设备,即第1行所示,H14w4是需要对齐的干扰向量,H14w4是H12w2和H13w3的线性组合,可以表示为:
H14w4=α11H12w2+α12H13w3(6)
同理,H15w5是第二个需要对齐的干扰向量,H15w5是H12w2和H13w3的线性组合,可以表示为:
H15w5=α13H12w2+α14H13w3(7)
对第1个用户设备,被对齐的干扰向量数为L-N=2;剩余干扰向量个数为N-1=2。同样,对第二、三、四、五个用户设备,可以得到多个形如(6),(7)的方程,则得到方程组。
二、下面将每个用户设备只接收单信号流的情形推广到每个用户设备接收多个信号流,假定第k个用户设备接收的信号流数为Lk≥1,这样,由于该用户设备本身的信号流已经占用了Lk维空间(该Lk维空间称为信号空间),所以需要将干扰向量压缩到剩余的N-Lk维空间中(该N-Lk维空间称为干扰空间),即:
Hkiwim,m=1,2,…,Li,i=1,2,…,K,i≠k
换言之,对第k个用户设备,L-Lk个干扰向量必须对齐到N-Lk维的干扰空间里。假定用πk表示从L-Lk个干扰向量中取出的作为干扰空间基向量的N-Lk个向量的集合,用πk’来表示余下的(L-Lk)-(N-Lk)=L-N个需要对齐的干扰向量的集合。例如,K=3,M=10,N=15,L=18,L1=L2=L3=6,则每个用户设备的干扰向量都是L-Lk=12个,干扰空间均为N-Lk=9维,即干扰空间基向量为9个,则各用户设备的πk和πk’可以表示为:
π1={w24,w25,w26,w31,…,w36},π′1={w21,w22,w23}
π2={w14,w15,w16,w31,…,w36},π′2={w11,w12,w13}
π3={w14,w15,w16,w21,…,w26},π′3={w11,w12,w13}
采用这样的表示方法,则干扰对齐可表示为:
上述(8)式表示将πk’中的向量表示成πk中向量的线性组合,所以(8)等价为如下的向量方程:
wjm′∈π′k,k=1,2,…,K (9)
其中r=1,2,...,N-Lk,n=(k-1)(L-N)+n’,n’=1,2,...,L-N,αnr表示第r个干扰向量在第n个向量方程中的系数,显然每个方程中αnr不能全部为零。特别地,如果Lk=N,此时第k个接收端没有干扰空间,则需要将M×1的零向量放入πk。。
三、如下将详细描述步骤102的具体实现过程:
为了便于后的预编码矩阵的计算,则将公式(8)变为如下:
公式(11)对每个wjm′∈π′k成立,其中k=1,2,...,K,r=1,2,...,N-Lk,n=(k-1)(L-N)+n’,n’=1,2,...,L-N,αn,0=-1,则用矩阵形式表示公式(11)为:
其中,
为N×(N-L
k+1)的矩阵,其第一列为H
kjw
jm′,w
jm′∈π′
k,余下为π
k中的N-L
k个干扰向量,即H
kiw
im′w
im∈π
k。
公式(11)可以进一步重写为:
Aw=0
其中,
为LM×1的级联预编码向量,其中,L1为第1个用户设备接收的信号流数;LK为第K个用户设备接收的信号流数;对应L个预编码向量,而A为(L-N)NK×LM的矩阵:
其中,每一个子矩阵Anl的维度均为N×M,n=(k-1)(L-N)+n’,n’=1,2,...,L-N,l=1,2,...,L,k=1,2,...,K。由于每个w对应一个数据流l,则第n个方程中的剩余干扰向量对应的Anl为αn,0与所述剩余干扰向量中的信道矩阵的乘积,第n个方程中的干扰空间基向量对应的Anl为所述干扰空间基向量中信道矩阵与所述干扰空间基向量的系数的乘积;具体的,对于第n个方程中对应的Anl可以表示为:
其中,l→wjm ′∈π′k表示第n个方程中,第l个数据流对应的向量为剩余干扰向量,即第l个数据流对应的向量wjm′∈π′k;l→wim∈πk表示第n个方程中,第l个数据流对应的向量为干扰空间基向量,即第l个数据流对应的向量wim∈πk。
为了清楚的说明上述公式(11)的变形,如下举例说明:假定K=4,L=6,流分布为[2211],即第1个用户设备和第2个用户设备分别有两个流,第3个用户设备和第4个用户设备分别有一个流,发射天线M和接收天线N均为4,这样,L-N就等于2,所以干扰对齐方程为:
因此:
如下将干扰对齐方案转化成最优化问题:
其中,约束条件为:
这里b=[1,0,...,0]
T,
表示第k个用户设备的第i个数据流对应的预编码向量的转置与所述预编码向量的乘积为该数据流的功率P
l,其中,k=1,2,3...K;i=1,2,....L
k,该P
l为总数据流中第l个数据流的功率。比如,L=6,流分布为[2211],即第1个用户设备和第2个用户设备分别有两个流,第3个用户设备和第4个用户设备分别有一个流,则第2个用户设备的第1个数据流为总数据流中的第3个流,该第2个用户设备的第1个数据流对应的预编码向量的转置与所述预编码向量的乘积为P
3。
当J足够小(趋近于零)的时候,可以认为找到了满足(15)的最优解,也就是(8)的解。
令D为一个L×LM的块矩阵:
那么,功率限制可以表示为:
Dw=e (17)
其中,e=[P1,P2,...,PL]T。因此(15)可以等效为最小化如下第一拉格朗日(Lagrange)函数:
其中λn和μ为拉格朗日因子,将D看成已知矩阵,则可对w求复微分如下:
令(19)等于0,则得到
由(15),可以把(18)重写成如下第二拉格朗日(Lagrange)函数:
对αn求复微分如下:
令(22)等于0,则得到
公式(23)需要矩阵
满秩,否则的话,可以直接令α
n为
的零空间向量,比如取最小特征值向量,并归一化至其第一个元素为-1。
将预设的αn或者w的初始值代入公式(23)和公式(20),通过公式(23)和公式(20)的循环迭代可以依次更新αn和w,可以在满足预设收敛条件时获得最优解,预设收敛条件可以是连续两次迭代的Lc值的差值小于第一预设门限值或者迭代次数达到一个预设的第二门限值,即达到了预设的迭代次数。
特别地,如果考虑发射总功率限制,最优化问题可以变为:
其中,P
T表示发射端发送各数据流的总功率;其中,
表示第k个用户设备的第i个数据流对应的预编码向量的转置与所述预编码向量的乘积;
表示每个接收端的每个数据流对应的预编码向量的转置与所述预编码向量的乘积之和为各数据流的总功率。
不难得到,给定αn时,可以得到最优的w为:
其中,um是A//A最小特征值对应的单位特征向量。通过公式(23)和公式(25)的循环迭代可以依次更新αn和w得到最优解。
如下简单说明如何通过迭代求解预编码矩阵。比如,把更新αn和w作为一次迭代周期,可以看到,在第i次迭代过程中,对任意的α(i-1)值,由公式(20)求得wn (i)值,然后将αn (i-1)值和w(i)值代入拉格朗日函数(18),将wn (i)值代入(23)就得到了α(i)值,然后利用将αn (i)值和w(i)值进行下一个迭代周期。因此,可以得到:
这表明该算法单调不增,趋向于收敛到一个较小值,因此,可以使用上述收敛条件。
前面讨论的干扰对齐是在接收端的干扰空间进行干扰对齐,也可以在发射端的干扰空间进行干扰对齐。
具体而言,即对第k个用户设备,共有L-Lk个干扰向量:
需要将L-Lk个干扰向量对齐到M-Lk维的干扰空间中,则表示为:
这里
和
的定义与π
k和π′
k的定义类似,这里不再赘述。类似地,可以获得如下干扰对齐方程:
其中r=1,2,…,M-L
k,n=(k-1)(L-M)+n’,n’=1,2,…,L-M,β
nr表示对应于第n个向量方程中第r个干扰向量的系数,显然每个方程中不能全部为零,类似地,当L
k=M时,没有干扰空间,此时,可令N×1的零向量为
的元素,这样,干扰对齐方法依然可以使用。
得到所有用户设备的所有均衡向量集合
(k=1,2,...,K)后,可以用下式得到各用户设备的预编码向量w
k:
其中ε=0表示用迫零均衡方法,ε=1表示用最小均方误差均衡方法,σ2为噪声功率,当然也可以用其他方法通过gk求取wk,不影响本发明的实现。
图3为依据本发明上述实施例所得到的频谱效率的数字仿真结果示意图,其中M为发射天线数,N为接收天线数,K为用户设备数,L为总流数,括号内各数字表示各用户设备的信号流数,可以看出在较高信噪比下,也不会出现BER平层。
参阅图4A和图4B,本发明实施例提供一种预编码矩阵获取装置,该装置可以为独立于发射端和接收端的中心控制设备,或者是位于发射端上,该装置具体包括:
获取单元401,用于针对每个接收端,从所有干扰向量中选取作为干扰空间基向量的干扰向量;
线性组合单元402,用于将剩余干扰向量表示为所述干扰空间基向量的线性组合;其中,所述剩余干扰向量为接收端的所有干扰向量中除去作为干扰空间基向量的干扰向量以外的干扰向量;
计算单元403,用于根据每个接收端剩余干扰向量所表示成的干扰空间基向量的线性组合及当前的信道状况,确定当前采用的预编码矩阵。
其中,干扰向量可以表示为信道向量与预编码向量的乘积;
此时获取单元401用于针对第k个接收端,确定干扰空间基向量的个数为N-LK个,从L-LK个干扰向量中选择N-LK个干扰向量作为干扰空间基向量;其中,N为第k个接收端的接收天线数,LK为第k个接收端的信号流数;L为多个接收端的信号流总数。
计算单元403包括:
第一函数获取子单元4031,用于根据每个接收端剩余干扰向量所表示成的干扰空间基向量的线性组合,得到矩阵
根据所述矩阵和当前的信道状况,利用最优化方法,得到第一拉格朗日函数;其中,
n=(k-1)(L-N)+n’,n’=1,2,...,L-N,α
n,0=-1,α
n,1、α
n,2、α
n,3、α
n(N-Lk)分别为第n个向量方程中第k个接收端的第n’个剩余干扰向量表示为干扰空间基向量线性组合时基向量的系数;
中的第一列为一剩余干扰向量,其余列为干扰空间基向量;
第二函数获取子单元4032,用于根据每个接收端剩余干扰向量所表示成的干扰空间基向量的线性组合,得到Aw=0;根据所述Aw=0和当前的信道状况,利用最优化方法,得到第二拉格朗日函数;l≤L,L为数据流总数,第n个方程中的剩余干扰向量对应的Anl为αn,0与所述剩余干扰向量中的信道矩阵的乘积,第n个方程中的干扰空间基向量对应的Anl为所述干扰空间基向量中信道矩阵与所述干扰空间基向量的系数的乘积;
第一迭代函数获取子单元4033,用于利用第一拉格朗日函数和第二拉格朗日函数,得到向量an和预编码矩阵w的互相关函数;
第一迭代子单元4034,用于将预设的w初始值或者向量an的初始值代入所述互相关函数,得到w当前值或者向量an的当前值,将所述w当前值或者向量an的当前值再代入所述互相关函数,直到满足预设条件为止,所述预设条件为将得到的w当前值或者向量an的当前值本次代入第一或者第二拉格朗日函数得到的值与前一次代入时得到的值的差达到第一预设门限,或者,迭代次数达到第二门限。
或者,干扰向量表示为信道向量与均衡向量的乘积;
此时,获取单元401用于针对第k个接收端,确定干扰空间基向量的个数为M-LK个;从L-LK个干扰向量中选择M-LK个干扰向量作为干扰空间基向量;其中,M为与第k个接收端配对的发送端的发射天线数,LK为第k个接收端的信号流数;L为多个接收端的信号流总数。
计算单元403包括:
第三函数获取子单元4036,用于根据每个接收端剩余干扰向量所表示成的干扰空间基向量的线性组合,得到矩阵
根据所述矩阵和当前的信道状况,利用最优化方法,得到第三拉格朗日函数;其中,
n=(k-1)(L-M)+n’,n’=1,2,...,L-N,β
n,0=-1,β
n,1、β
n,2、β
n,3、
分别为第n个向量方程中第k个接收端的第n’个剩余干扰向量表示为干扰空间基向量线性组合时基向量的系数;
中的第一列为一剩余干扰向量,其余列为干扰空间基向量;
第四函数获取子单元4037,用于根据每个接收端剩余干扰向量所表示成的干扰空间基向量的线性组合,得到BG=0;根据所述Bw=0和当前的信道状况,利用最优化方法,得到第四拉格朗日函数;其中,B中的每一个子矩阵Bnl的维度均为M×N,l≤L,L为数据流总数,第n个方程中的剩余干扰向量对应的Bnl为βn,0与所述剩余干扰向量中的信道矩阵的乘积,第n个方程中的干扰空间基向量对应的Bnl为所述干扰空间基向量中信道矩阵与所述干扰空间基向量的系数的乘积;
第二迭代函数获取子单元4038,用于利用第三拉格朗日函数和第四拉格朗日函数,得到向量an和均衡矩阵G的互相关函数;
第二迭代子单元4039,用于将预设的均衡矩阵G的初始值或者向量βn的初始值代入所述互相关函数,得到均衡矩阵G的当前值或者向量βn的当前值,将所述均衡矩阵G的当前值或者向量βn的当前值再代入所述互相关函数,直到满足预设收敛条件为止,所述预设收敛条件为将得到的均衡矩阵G的当前值或者向量βn的当前值本次代入第三或者第四拉格朗日函数得到的值与前一次代入时得到的值的差达到第三预设门限,或者,迭代次数达到第四门限。
本发明实施例中,针对每个接收端,都将需要对齐的干扰向量表示为干扰空间基向量的线性组合,这样达到干扰完全对齐,能够完全消除干扰,不会产生BER平层。
需要说明的是,本发明上述各实施例中的每对发射端和接收端都可以拥有不同的天线配置,不影响本发明的实现。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,例如只读存储器,磁盘或光盘等。
以上对本发明实施例所提供的预编码矩阵获取方法及装置进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。