发明内容
基于现有技术存在的问题,本发明的目的在于提供一种星地融合网络中联合广播和单播的波束赋形方法。该方法首先建立STIN下行单播和广播传输系统模型,提出基于层分复用(Layer Division Multiplexing,LDM)的联合波束赋形设计方案,并考虑非理想CSI的影响,将信道误差建模为范数有界误差模型,根据最差情况准则建立基于QoS约束的鲁棒波束赋形设计问题,以最小化系统的传输功率。然后采用半正定松弛(Semi-DefiniteRelaxation,SDR)的方法,结合S-Procedure及其扩展定理,将具有无穷维约束的非凸优化问题转化为具有线性矩阵不等式(Linear Matrix Inequality,LMI)的确定性优化形式。为了解决秩约束,采用一种基于罚函数的迭代算法,最终原问题被转换为一个标准半正定规划(Semi-Definite Programming,SDP)问题。
本发明采用如下技术方案解决上述技术问题:
一种星地融合网络中联合广播和单播的波束赋形方法,所述方法包括:
建立星地融合网络的下行单播和广播传输系统,即卫星和地面基站共享相同的频带,通过波束赋形以协作方式向覆盖范围内的用户提供单播或广播服务;
基于层分复用方式,将信道误差建模为范数误差模型,根据最差情况准则建立基于QoS约束的鲁棒波束赋形的最小化系统的传输功率模型;
采用半正定松弛方法、S-Procedure以及其扩展方法,将上述传输功率模型中的非凸优化问题转换为确定性优化形式,并在传输功率模型中增加卫星向用户应用的单播或广播波束赋形向量的秩约束,以及地面基站向用户应用的单播或广播波束赋形向量的秩约束;
在不考虑秩约束时,计算所述传输功率模型,获得波束赋形矩阵的初始可行解;
将罚函数替代所述秩约束,将所述初始可行解作为迭代初始值,迭代计算所述传输功率模型,获得波束赋形矩阵的优化解;
通过特征值分解的方式,对计算获得的波束赋形矩阵的优化解进行处理,得到最终的波束赋形向量。
本发明的有益效果:
本发明是一种星地融合网络中联合广播和单播鲁棒波束赋形方法。本发明针对下行广播和单播单独传输性能受限问题,首次提出了星地融合网络中的联合单播和广播传输方法,以优化传输功率为目标,提出了一种基于LDM的鲁棒波束赋形设计问题,并且采用SDR的方法,结合S-Procedure及其扩展定理,将具有无穷维约束的非凸优化问题转化为具有LMI的确定性优化形式。为了解决秩1约束,本发明采用一种基于IPF的算法,将原问题转换为一个标准SDP问题进行求解,能够快速获得波束赋形矩阵的优化解。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提出了一种星地融合网络中联合单播和广播传输波束赋形方案,同时考虑了信道状态信息(Channel State Information,CSI)在非理想的情况下的鲁棒性传输问题。
图1是本发明实施例中的一种星地融合网络中联合广播和单播的波束赋形方法流程图,如图1所示,所述方法包括:
101、建立星地融合网络的下行单播和广播传输系统,即卫星和地面基站共享相同的频带,通过波束赋形以协作方式向覆盖范围内的用户提供单播或广播服务;
在本发明实施例中,假设在STIN下行传输系统由一个地面基站和一颗卫星以及若干用户构成。在该系统中,卫星和地面基站共享相同的频带同时为系统中的用户提供单播和广播服务。卫星和地面基站通过波束赋形以协作方式向所有各自覆盖范围内的用户提供单播或广播服务。
在一些是优选实施例中,地面基站仅为其覆盖范围内的用户提供单播服务,地面基站覆盖范围外的用户由卫星提供单播服务;而地面基站和卫星都可以为所有用户提供广播服务。
图2是本发明实施例采用的星地融合网络下行传输系统模型,如图2所示,假设卫星配备有M根天线,地面基站配备有N根天线。卫星天线和地面基站天线通过回程链路由中央处理器(Central Processor,CP)协调。假设回程链路容量有限,系统在一个传输周期内共调度的用户数为U
tot,且所有用户均为单天线用户。根据用户接收单播信号的模式,这些用户可分为U
S个卫星用户和U
T个地面用户,其中U
tot=U
S+U
T。将卫星用户的用户索引集表示为
地面用户的用户索引集表示为
其中,
102、基于层分复用方式,将信道误差建模为范数误差模型,根据最差情况准则建立基于QoS约束的鲁棒波束赋形的最小化系统的传输功率模型;
在本发明实施例中,在步骤101的一些实施例的基础上,可以得到:
其中,
表示卫星到第i个卫星用户的共轭信道向量;
表示卫星向所有卫星用户应用的广播波束赋形向量;
表示广播信号,
表示卫星向第i个卫星用户应用的单播波束赋形向量;
表示卫星发送给第i个卫星用户的单播信号,并且满足
n
S,i表示卫星向第i个卫星用户信道的噪声。
其中,
表示基站发送给第j个地面用户的单播信号,并且满足
和
分别表示地面基站到第j个地面用户的信道向量和单播波束赋形向量,
表示地面基站到第j个地面用户u
T,j之间的共轭信道向量;w
T,j表示地面基站向第j个地面用户应用的单播波束赋形向量;
表示基站向所有地面用户应用的广播波束赋形向量。
表示卫星到第j个地面用户之间的信道向量。n
T,j表示地面基站向第j个地面用户信道的噪声。
由于采用了层分复用的方法,本发明可以将单播信号视为噪声,在接收端进行解码,接收端在消除解码的广播信号后解码单播信号。因此,对于卫星用户uS,i和地面用户uT,j,各自的广播与单播信干噪比(Signal to Interference plus Noise Ratios,SINR)可以分别表示为:
其中,
表示第i个卫星用户u
S,i的广播信干噪比;
表示第i个卫星用户u
S,i的单播信干噪比;
表示第j个地面用户u
T,j的广播信干噪比;
表示第j个地面用户u
T,j的单播信干噪比;
表示第i个卫星用户u
S,i的噪声功率;
表示第j个地面用户u
T,j的噪声功率。
根据香农定理,对于任意用户u∈S∪T的广播和单播速率分别表示为:
在本发明实施例中,建立基于QoS约束的鲁棒波束赋形设计问题,以最小化系统的传输功率。令
为了更符合实际场景,本发明考虑非理想CSI场景下的鲁棒波束赋形方案,基于范数有界误差模型,信道不确定性集合可以表示为:
其中,G
S,i表示卫星到第i个卫星用户u
S,i之间的信道不确定性集合;G
T,j表示卫星到第j个地面用户u
T,j之间的信道不确定性集合;H
T,j表示地面基站到第j个地面用户u
T,j之间的信道不确定性集合;Δh
T,j表示地面基站到第j个地面用户u
T,j的误差估计信道向量,h
T,j表示地面基站到第j个地面用户u
T,j的信道向量,
表示地面基站到第j个地面用户u
T,j的估计信道向量;Δg
S,i表示卫星到第i个卫星用户u
S,i之间的误差估计信道向量,g
S,i表示卫星到第i个卫星用户u
S,i之间的信道向量,
表示卫星到第i个卫星用户u
S,i之间的估计信道向量;Δg
T,j表示卫星到第j个地面用户u
T,j之间的误差估计信道向量,g
T,j表示卫星到第j个地面用户u
T,j之间的信道向量,
表示卫星到第j个地面用户u
T,j之间的估计信道向量;ε
S,i表示卫星到第i个卫星用户u
S,i之间的信道估计误差上限;ε
T,j表示卫星到第j个地面用户u
T,j之间的信道估计误差上限;χ
T,j表示地面基站到第j个地面用户u
T,j之间的信道估计误差上限。
因此,基于QoS约束的鲁棒波束赋形设计问题可以建模为:
其中,v表示卫星向用户应用的单播或广播波束赋形向量集合,
U
S表示卫星用户个数,S表示卫星用户索引集合;w表示地面基站向用户应用的单播或广播波束赋形向量集合,
U
T表示地面用户个数,T表示地面用户索引集合;v
S,k表示卫星向卫星用户k应用的单播或广播波束赋形向量;w
T,n表示地面基站向地面用户n应用的单播或广播波束赋形向量;C1表示卫星用户的单播速率约束;C2表示卫星用户的广播速率约束;C3表示地面用户的单播速率约束;C4表示地面用户的广播速率约束;
表示第i个卫星用户u
S,i的单播速率;
表示卫星用户的最小单播速率门限;
表示第i个卫星用户u
S,i的广播速率;
表示卫星用户的最小广播速率门限;
表示第j个地面用户u
T,j的单播速率;
表示地面用户的最小单播速率门限;
表示第j个地面用户u
T,j的广播速率;
表示地面用户的最小广播速率门限;Δh
T,j表示地面基站到第j个地面用户u
T,j的误差估计信道向量;Δg
S,i表示卫星到第i个卫星用户u
S,i之间的误差估计信道向量;Δg
T,j表示地面基站到第j个地面用户u
T,j之间的误差估计信道向量;G
S,i表示卫星到第i个卫星用户u
S,i之间的信道不确定性集合;G
T,j表示地面基站到第j个地面用户u
T,j之间的信道不确定性集合;H
T,j表示地面基站到第j个地面用户u
T,j之间的信道不确定性集合。
103、采用半正定松弛方法、S-Procedure以及其扩展方法,将上述传输功率模型中的非凸优化问题转换为确定性优化形式,并在传输功率模型中增加卫星向用户应用的单播或广播波束赋形向量的秩约束,以及地面基站向用户应用的单播或广播波束赋形向量的秩约束;
在本发明实施例中,考虑到上述基于QoS约束的鲁棒波束赋形的最小化系统的传输功率模型存在信道误差,而鲁棒波束赋形设计问题属于无穷维的非凸优化问题。为了解决这个难题,本发明采用S-Procedure定理及其扩展定理将原优化问题的约束转换为一个等价的形式。S-Procedure定理描述如下:
本发明首先定义
其中
m=1,2。如果存在
满足
那么
成立的充要条件为:存在τ≥0满足:
利用第一辅助变量
和第二辅助变量
对约束C1~C4依次进行转换表示为:
其中,上标H表示共轭转置矩阵;Δg
S,i表示卫星到第i个卫星用户u
S,i之间的误差估计信道向量,g
S,i表示卫星到第i个卫星用户u
S,i之间的信道向量,
表示卫星到第i个卫星用户u
S,i之间的估计信道向量;Δh
T,j表示地面基站到第j个地面用户u
T,j的误差估计信道向量,h
T,j表示地面基站到第j个地面用户u
T,j的信道向量,
表示地面基站到第j个地面用户u
T,j的估计信道向量;Δg
T,j表示卫星到第j个地面用户u
T,j之间的误差估计信道向量,g
T,j表示卫星到第j个地面用户u
T,j之间的信道向量,
表示卫星到第j个地面用户u
T,j之间的估计信道向量;
表示卫星用户的单播信干噪比阈值;
表示地面用户的单播信干噪比阈值;
表示卫星用户的广播信干噪比阈值;
表示地面用户的广播信干噪比阈值。
对于约束C1-C4,可以根据S-Procedure分别转变为如下线性矩阵不等式的形式,依次表示为:
其中,I
M×M代表M×N大小的单位矩阵,
τ
1,i表示第i个卫星用户对应的第一松弛变量;τ
2,i表示第i个卫星用户对应的第二松弛变量;τ
3,j表示第j个地面用户对应的第一松弛变量;τ
3,j表示第j个地面用户对应的第二松弛变量。
由于仍含有不确定性信道向量gT,j,约束C3和C4仍然是无穷维的。因此,根据如下引理
若D≥0,且矩阵Hq(q=1,…,6)满足:
则存在μ≥0使得上式等价于:
因此,可以将约束C3和C4重新表述为:
其中,μ1,j表示第j个地面用户对应的第三松弛变量;μ2,j表示第j个地面用户对应的第四松弛变量;且μ1,j≥0和μ2,j≥0。
由于采用了半正定松弛方法,具有线性矩阵不等式的鲁棒波束赋形设计问题的目标函数变为如下形式:
其中,tr(·)表示矩阵的迹运算。同时增加约束C5和C6:
因此,优化后的基于QoS约束的鲁棒波束赋形的最小化系统的传输功率模型表示为:
其中,
和
表示辅助变量集合,U
S表示卫星用户个数,U
T表示地面用户个数;
和
表示辅助变量;v
S,k表示卫星向卫星用户k应用的单播或广播波束赋形向量;w
T,n表示地面基站向地面用户n应用的单播或广播波束赋形向量;C1表示卫星用户的单播速率约束;C2表示卫星用户的广播速率约束;C3表示地面用户的单播速率约束;C4表示地面用户的广播速率约束;C5表示波束赋形向量的正定约束;C6表示波束赋形向量的秩约束;
表示变量集合;Y
1,i表示第i个卫星用户对应的第一松弛参数;Y
2,i表示第i个卫星用户对应的第二松弛参数;Y
3,j表示第j个地面用户对应的第一松弛参数;rank()表示矩阵的秩。
104、在不考虑秩约束时,计算所述传输功率模型,获得波束赋形矩阵的初始可行解;
在本发明实施例中,在不考虑约束C6的情况下,计算上述传输功率模型,即可以直接使用如下模型,可以得到波束赋形矩阵的初始可行解。
105、将罚函数替代所述秩约束,将所述初始可行解作为迭代初始值,迭代计算所述传输功率模型,获得波束赋形矩阵的优化解;
在本发明实施例中,由于秩约束C6的存在,鲁棒波束赋形设计问题仍然是非凸的。为了保证取得秩1的解
和
的合法性,需要在目标函数中插入迭代罚函数。首先,使用如下等式替换秩1约束C6:
tr(VS,k)-λS,k,max=0
tr(WT,n)-λT,n,max=0
其中,λ
S,k,max和λ
T,n,max分别为
和W
T,n的最大特征值。可以建立惩罚函数来改进鲁棒波束赋形设计问题的目标函数,可以表示为:
其中,αk表示卫星用户k的惩罚因子,βn表示地面用户n的惩罚因子;tr表示求矩阵的迹。
由于惩罚函数的存在,目标函数仍然是非凸的。为了解决该问题,采用一种迭代的方法来获得一个凸的目标函数。对于第t次迭代的解
和
有如下不等式
其中,v
S,k,max和w
T,n,max分别为λ
S,k,max和λ
T,n,max对应的单位特征向量;α
k表示卫星用户k的惩罚因子,β
n表示地面用户n的惩罚因子;tr表示求矩阵的迹;右上标H表示求共轭矩阵;右上标(t)表示第t次迭代;右上标(t+1)表示第t+1次迭代;
表示在第t次迭代过程中卫星向卫星用户k应用的单播或广播的最大波束赋形向量;
表示在第t次迭代过程中地面基站向地面用户n应用的单播或广播的最大波束赋形向量。
最终,基于迭代罚函数的鲁棒波束赋形设计问题的目标函数可以表述为:
s.t.C1,C2,C3,C4,C5.
其中,惩罚因子α
k和β
n的迭代过程为
其中,υ
k和κ
n为惩罚因子更新系数。上述问题为一个标准SDP问题,可以通过CVX等工具求解。
106、通过特征值分解的方式,对计算获得的波束赋形矩阵的优化解进行处理,得到最终的波束赋形向量。
为便于理解,现将本发明实施例波束赋形向量的求解过程单独表述为如下:
(1)初始化t=0,精度ι
1,ι
2,最大迭代次数T
max=20,惩罚因子
和系数υ
k,κ
n;
(2)不考虑秩1约束C6的前提下,求解具有LMI的鲁棒波束赋形优化问题,得到初始可行解
(3)将罚函数替代所述秩1约束C6,进入求解具有LMI的鲁棒波束赋形优化问题的循环过程;
(4)令t=t+1,求解基于IPF的鲁棒波束赋形优化问题,获得第t次迭代过程中的
(5)判断是否满足最大迭代次数或满足收敛条件,满足则跳出循环,否则更新惩罚因子
继续循环;
(6)通过特征值分解(Eigenvalue Decomposition,EVD)获得波束赋形向量
和
在一些实施例中,为验证本发明在星地融合网络中的联合广播和单播的波束赋形方法,本实施例中考虑一个地面基站和一个卫星的场景,假设卫星和地面信道均为频率平坦的准静态信道。设置载波频率为2GHz,带宽B为10MHz,地面基站的覆盖半径设置为500m,并假设卫星为高度1000km的LEO卫星。默认卫星天线数M=4,地面基站天线数N=4,卫星用户数U
S=4,地面用户数U
T=4。将卫星信道建模为莱斯信道,将地面信道建模为瑞利信道。假设环境噪声功率为-134dBm,用于对系统中相应参数进行归一化,因此设置噪声方差
默认卫星和地面用户的单播和广播最小速率门限满足
假设本发明中的信道估计误差的上限均相等,默认ε=ε
S,i=ε
T,j=χ
T,j=0.03。惩罚因子α
k和β
n初值为1,惩罚因子更新系数υ
k=1.5,κ
n=1.2,收敛精度ι
1=ι
2=10
-4。
如图3所示,随着用户的信道估计误差ε增大,达到收敛所需的迭代次数越多。但算法均能够在8次迭代内收敛,证明了本发明所提方法能够快速收敛。同时还可以观察到,对于较大的ε,需要较高的传输功率来满足鲁棒的QoS性能。
如图4所示,无论是在广播层还是单播层,非鲁棒算法的最小用户传输速率随着信道估计误差的增大而不断降低。相比之下,本发明所提鲁棒算法对信道估计误差不敏感,即能够保证用户的QoS需求,这证明了本发明所提方法具有很好的鲁棒性。同时,随着ε的增加,鲁棒算法的精度降低,但为了保证用户的QoS需求,使得最小用户速率随之增大,与此同时将会消耗更多的传输功率,这一结论与图3中的仿真结果相符合。
图5给出了系统传输功率与广播最小速率门限
的关系,对比了本发明所提非正交LDM方法与传统正交TDM方法在传输功率方面的性能差异。对于TDM方法,每个传输时隙T被分成两个子时隙:持续时间为T
0的子时隙用于单播传输,持续时间为T-T
0的子时隙用于广播传输。令T
0/T表示单播传输的时间比。从图5中可以看出,与TDM方法相比,LDM方案能够显著降低STIN的传输功率。不同T
0/T的TDM方法随着
的变化呈现不同的性能。具体来说,当T
0/T=0.6时,单播传输子时隙的占比大,广播传输子时隙占比小,传输功率随着
的增大快速增加;当T
0/T=0.4时,传输功率随着
的增大只有小幅增加。因此可以看出,T
0/T的值越大,传输功率随
变化的幅度越大,反之越小,从而导致对于不同的
TDM方案性能最优时的T
0/T不同。尽管如此,本发明所提方法仍然可以比T
0/T最优的TDM方法降低约6dBm的传输功率。
图6给出了在协作与非协作两种不同的传输方案下,用户单播平均速率和广播平均速率与最小速率门限
和
的关系。从图6(a)中可以看出,本发明所提协作传输方案的单播平均速率总是高于非协作传输方案,并且随着
的增加,采用协作传输带来的增益越大。这是因为
的增大将导致卫星信号传输功率增加,进而增大了卫星在单播传输过程中对地面用户的干扰。相比于非协作传输方案,本发明所提协作传输方案能够有效地缓解卫星单播信号对地面用户的干扰。从图6(b)中可以看出,虽然协作与非协作传输方案的系统广播平均速率均能够满足最小广播速率要求,但协作传输方案广播速率要比非协作传输方案更接近最小广播速率门限。这说明在广播层,协作传输方案中本文所提算法的收敛精度要高于非协作传输方案。
本发明对STIN中的下行单播和广播传输系统进行研究。在该系统中,地面基站和卫星均使用相同的时频资源,以协作的方式将专属的单播信号发送给特定用户,并将一个公共的广播信号发送给所有用户。提出基于LDM的鲁棒波束赋形设计方法,考虑非理想CSI情况下基于QoS约束的鲁棒优化问题。采用SDR方法,结合S-Procedure和IPF将具有无穷维约束的非凸优化问题转化为具有LMI的确定性凸问题迭代进行求解。仿真结果表明,本文方案具有较低的传输功耗和较好的抗干扰性能。
在本发明的描述中,需要理解的是,术语“同轴”、“底部”、“一端”、“顶部”、“中部”、“另一端”、“上”、“一侧”、“顶部”、“内”、“外”、“前部”、“中央”、“两端”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
在本发明中,除非另有明确的规定和限定,术语“安装”、“设置”、“连接”、“固定”、“旋转”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。