CN110417445B - 网络辅助全双工系统的稀疏波束设计与功率控制方法 - Google Patents

网络辅助全双工系统的稀疏波束设计与功率控制方法 Download PDF

Info

Publication number
CN110417445B
CN110417445B CN201910700869.XA CN201910700869A CN110417445B CN 110417445 B CN110417445 B CN 110417445B CN 201910700869 A CN201910700869 A CN 201910700869A CN 110417445 B CN110417445 B CN 110417445B
Authority
CN
China
Prior art keywords
remote radio
radio frequency
downlink
frequency head
representing
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.)
Active
Application number
CN201910700869.XA
Other languages
English (en)
Other versions
CN110417445A (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201910700869.XA priority Critical patent/CN110417445B/zh
Publication of CN110417445A publication Critical patent/CN110417445A/zh
Application granted granted Critical
Publication of CN110417445B publication Critical patent/CN110417445B/zh
Active 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
    • H04B7/0426Power distribution
    • 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
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W52/00Power management, e.g. TPC [Transmission Power Control], power saving or power classes
    • H04W52/04TPC
    • H04W52/38TPC being performed in particular situations
    • H04W52/42TPC being performed in particular situations in systems with time, space, frequency or polarisation diversity

Landscapes

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

Abstract

本发明公开了网络辅助全双工系统的稀疏波束设计与功率控制方法,以最大化系统和速率为目标,以满足指定的系统服务质量QoS、发射机的发射功率以及回程链路消耗为约束条件,建立了以远程无线电射频头发射机与接收机波束,以及上行用户功率为设计变量的数学优化模型。本方法通过基于迭代序贯参数凸逼近方法,将原始的非凸优化问题通过不断连续凸逼近方法逐步转化为凸优化问题,并通过不断迭代优化联合求出稀疏波束形成向量与发射功率的值。

Description

网络辅助全双工系统的稀疏波束设计与功率控制方法
技术领域
本发明涉及无线通信传输技术领域,尤其是网络辅助全双工系统的稀疏波束设计与功率控制方法。
背景技术
在过去几十年中,智能设备和视频流应用的激增导致了对上行链路和下行链路传输中更高数据速率的需求激增。目前,5G仍采用传统的时分双工和频分双工技术。然而,随着对上行下行业务不对称需求的不断增长,我们需要突破更灵活的双工模式,进一步提高频谱效率。为了使收发器同时使用相同的频率资源工作,使通信双方能够同时使用相同的频率,提出了全双工方案,其他方案,诸如混合半/全双工蜂窝网络、α-双工方案等双工模式仅限于每个小区一个或两个基站,小区边缘用户可能会获得全双工无线资源,从而导致严重干扰。
发明内容
本发明所要解决的技术问题在于,提供网络辅助全双工系统的稀疏波束设计与功率控制方法,能够使系统的和速率尽可能达到最大化。
为解决上述技术问题,本发明提供网络辅助全双工系统的稀疏波束设计与功率控制方法,包括如下步骤:
(1)发送远程无线电射频头通过信道估计获得它到所有下行用户设备和接收远程无线电射频头间的理想信道状态信息,上行用户通过信道估计获得它到所有下行用户设备和接收远程无线电射频头间的理想信道状态信息;假设系统采用基于网络辅助全双工模式的时分双工制式,且信道服从平坦衰落,即在信道相干时间内信道系数保持不变;
(2)中央处理器发送下行预编码向量到发送远程无线电射频头,发送远程无线电射频头接收该信号并将收到的信号,经过信道hD,k转发给下行用户设备,下行用户设备最终接收到的下行信息为
Figure BDA0002150745660000011
其中,
Figure BDA0002150745660000012
表示发送远程无线电射频头到所有下行用户设备的理想信道状态信息,hIUI,j,k表示第j个上行用户到第k个下行用户设备的信道且服从复高斯分布
Figure BDA0002150745660000021
Figure BDA0002150745660000022
sD,k,sU,j分别表示发送远程无线电射频头发送给下行用户设备k及上行用户j发送的数据信息;wD,k,pU,j分别表示发送远程无线电射频头发送给下行用户设备k的预编码向量及上行用户设备j的发射功率;
Figure BDA0002150745660000023
表示加性高斯白噪声;上行传输系统中,先对发送远程无线电射频头到接收远程无线电射频头的干扰功率做干扰消除处理,并对参与上行传输的用户及上行接收的远程无线电射频头做好配对;这样,接收远程无线电射频头接收到的上行用户信息为
Figure BDA0002150745660000024
其中,
Figure BDA0002150745660000025
表示干扰消除后的信道估计误差项,
Figure BDA0002150745660000026
表示加性高斯白噪声,
Figure BDA0002150745660000027
表示上行用户通过信道估计到接收远程无线电射频头间的理想信道状态信息;
(3)整个传输系统以最大化和速率为准则,以指定的上行及下行用户设备传输最小系统频谱效率,上行用户的功率及下行远程天线的预编码发射功率,下行回程容量为约束,求解传输系统的稀疏波束形成与功率控制优化问题;
建立问题模型为:
Figure BDA0002150745660000028
C1:
Figure BDA0002150745660000029
C2:
Figure BDA00021507456600000210
C3:
Figure BDA00021507456600000211
C4:RD,k≥RD,min,k
C5:RU,j≥RU,min,j
C6:
Figure BDA00021507456600000212
其中,
Figure BDA00021507456600000213
Figure BDA0002150745660000031
Figure BDA0002150745660000032
Figure BDA0002150745660000033
Figure BDA0002150745660000034
uU,j,z为接收远程无线电射频头z为解调上行用户j发射的信号而设计的接收机,接收远程无线电射频头z与上行用户采用预先配对模式,即接收远程无线电射频头z仅可解调与之配对的上行用户信号;其中,一个接收远程无线电射频头可以解调多个上行用户信号,而一个上行用户仅可以配对一个接收远程无线电射频头,CD,l,
Figure BDA0002150745660000035
RU,min,j,RD,min,k分别为预设的第l个发送远程无线电射频头的回程约束,发射功率约束,上行用户j的发射功率约束,最低速率要求,下行用户k的速率要求,
Figure BDA0002150745660000036
表示发送远程无线电射频头到所有下行用户设备的理想信道状态信息,hIUI,j,k表示第j个上行用户到第k个下行用户设备的信道且服从复高斯分布
Figure BDA0002150745660000037
Figure BDA0002150745660000038
sD,k,sU,j分别表示发送远程无线电射频头发送给下行用户设备k及上行用户j发送的数据信息;wD,k,pU,j分别表示发送远程无线电射频头发送给下行用户设备k的预编码向量及上行用户设备j的发射功率;
Figure BDA0002150745660000039
表示加性高斯白噪声,
Figure BDA00021507456600000310
表示干扰消除后的信道估计误差项,
Figure BDA00021507456600000311
表示加性高斯白噪声,
Figure BDA00021507456600000312
表示上行用户通过信道估计到接收远程无线电射频头间的理想信道状态信息;
Figure BDA00021507456600000313
为指示函数,具体定义为:
Figure BDA0002150745660000041
(4)由于指示函数为离散整数规划问题,无法直接通过优化求解处理,此处利用如下函数逼近该问题,即C7:
Figure BDA0002150745660000042
其中wD,l,k表示第l个发送远程无线电射频头到第k个下行用户设备的预编码向量,θ>>1为一个远大于1的实数;逼近后的函数由于是原问题的近似解,并不是原问题的准确解;因此,原问题可以分步骤完成;
(5)根据对数函数的单调性观察,首先将目标函数重新建模为:
Figure BDA0002150745660000043
C8:
Figure BDA0002150745660000044
C9:
Figure BDA0002150745660000045
C10:1+rD,k≥χD,k,
C11:1+rU,j≥χU,j,
C12:χU,j≥1,χD,k≥1,
C1,C2,C6,C7,
其中,χD,k,χU,j为针对下行用户设备速率以及上行用户设备速率而新引入的优化变量;
(6)在此基础上,特别地,将目标函数转化为二阶锥规划形式:
Figure BDA0002150745660000046
C13:
Figure BDA0002150745660000047
C14:
Figure BDA0002150745660000048
其中,
Figure BDA0002150745660000049
为新引入变量;
(7)根据一般的凸优化问题求解方法,利用迭代序贯参数凸逼近方法将原问题转化成凸问题求解,即步骤S1;
(8)判断原问题相邻两个解之间的差值是否大于阈值,即步骤S2;
(9)根据步骤(5)到步骤(8),重复步骤S1、步骤S2并返回下行远程天线单元与下行用户服务配对模式,即步骤S3;
(10)给定下行远程天线单元与下行用户的服务配对模式,利用利用迭代序贯参数凸逼近及块坐标下降方法求解原问题并重复步骤S1和步骤S2,即步骤S4。
优选的,步骤(1)中,发送远程无线电射频头通过信道估计获得它到所有下行用户设备和接收远程无线电射频头间的理想信道状态信息,即信道矩阵
Figure BDA0002150745660000051
Figure BDA0002150745660000052
其中,
Figure BDA0002150745660000053
表示第l个发送远程无线电射频头到第k个下行用户设备的信道向量且服从复高斯分布
Figure BDA0002150745660000054
Figure BDA0002150745660000055
表示第l个发送远程无线电射频头到第z个接收远程无线电射频头的信道向量且服从复高斯分布;上行用户通过信道估计获得它到所有下行用户设备和接收远程无线电射频头间的理想信道状态信息,即信道矩阵
Figure BDA0002150745660000056
Figure BDA0002150745660000057
其中,hIUI,j,k表示第j个上行用户到第k个下行用户设备的信道且服从复高斯分布
Figure BDA0002150745660000058
hU,j,z表示第j个上行用户到第z个接收远程无线电射频头的信道向量且服从复高斯分布
Figure BDA0002150745660000059
优选的,步骤(3)中,求解传输系统的联合稀疏波束形成与功率控制优化问题具体为:
Figure BDA00021507456600000510
C1:
Figure BDA00021507456600000511
C2:
Figure BDA00021507456600000512
C3:
Figure BDA00021507456600000513
C4:RD,k≥RD,min,k
C5:RU,j≥RU,min,j
C6:
Figure BDA00021507456600000514
其中,
Figure BDA00021507456600000515
Figure BDA0002150745660000061
Figure BDA0002150745660000062
Figure BDA0002150745660000063
Figure BDA0002150745660000064
为指示函数,具体定义为:
Figure BDA0002150745660000065
CD,l,
Figure BDA0002150745660000066
RU,min,j,RD,min,k分别为预设的第l个发送远程无线电射频头的回程约束,发射功率约束,上行用户j的发射功率约束,最低速率要求,下行用户k的速率要求。
优选的,将原问题的C7,C8,C9,C10,C11近似问题分别转化为:
C7→C15:
Figure BDA0002150745660000067
C8→C16:
Figure BDA0002150745660000068
C9→C17:
Figure BDA0002150745660000069
C10→C18:
Figure BDA00021507456600000610
C11→C19:
Figure BDA00021507456600000611
其中,
Figure BDA00021507456600000612
ωU,j为新引入的变量,
Figure BDA00021507456600000613
为上行用户j的干扰功率上界的二阶锥表示形式变量及一般表示形式变量,
Figure BDA00021507456600000614
Figure BDA00021507456600000615
(7.1)①针对C7→C15:首先新引入变量
Figure BDA0002150745660000071
同时,利用一阶泰勒近似原理,得到如下结果:
C15,C20:
Figure BDA0002150745660000072
C21:
Figure BDA0002150745660000073
其中,
Figure BDA0002150745660000074
对C20,C21分别进行泰勒展开进行凸逼近如下:
C22:
Figure BDA0002150745660000075
C23:
Figure BDA0002150745660000076
其中,
Figure BDA0002150745660000077
②针对C8→C16:利用一阶泰勒近似原理将约束C8的等式右边展开,得到C16;
③针对C9→C17:首先通过引入优化变量
Figure BDA0002150745660000078
Figure BDA0002150745660000079
βU,j,j′
Figure BDA00021507456600000710
求出C9的近似处理表达式,如下:
C24:
Figure BDA00021507456600000711
C25:
Figure BDA00021507456600000712
C26:
Figure BDA00021507456600000713
C27:
Figure BDA00021507456600000714
C28:
Figure BDA00021507456600000715
C29:
Figure BDA00021507456600000716
通过进一步凸逼近,C25,C28,C29可以转换为:
C30:
Figure BDA00021507456600000717
C31:
Figure BDA0002150745660000081
C32:
Figure BDA0002150745660000082
这样,C9可以直接转换为二阶锥规划形式C17,其中,
Figure BDA0002150745660000083
④针对C10→C18:约束C10可以变换为:
C33:
Figure BDA0002150745660000084
观察到约束C15的不等式右边为凸形式,所以,根据一阶泰勒近似原理展开,可以表示为C18;
⑤针对C11→C19:首先,将约束C11变换为:
C34:
Figure BDA0002150745660000085
进一步不等式可以处理如下:
C19,C35:
Figure BDA0002150745660000086
其中,ωU,j,为新引进变量,
Figure BDA0002150745660000087
这里,利用序贯参数凸逼近方法,C35可以进一步近似为:
C36:
Figure BDA0002150745660000088
C37:
Figure BDA0002150745660000089
其中,
Figure BDA00021507456600000810
最终,C35可以近似为:
C39:
Figure BDA0002150745660000091
优选的,迭代过程引入信任域约束方法,即:
C40:
Figure BDA0002150745660000092
C41:τ(PU,j)={|PU,j-PU,j|≤ηU,j}
τ(wD,l,k),τ(PU,j)表示信任域函数,ηD,k,l,ηU,j表示信任域阈值,当目标函数的相邻两次迭代结果的差大于阈值时,则ηD,k,l和ηU,j的值可以设为稍微大一点;反之,随着迭代步骤的进行,ηD,k,l和ηU,j的值可以相应的随之减小。
优选的,得到n+1次的迭代表达式如下:
Figure BDA0002150745660000093
优选的,给定下行远程天线单元与下行用户的服务配对模式,原问题可以转化为:
Figure BDA0002150745660000094
s.t. C1,C2,C6,C12,C13,C14,C16,
C17,C18,C19,C22,C23,C24,
C27,C30,C31,C32,C39,C40,C41
Figure BDA0002150745660000095
Figure BDA0002150745660000096
其中,
Figure BDA0002150745660000097
表示发射远程天线单元l与下行用户之间的关联集。
本发明的有益效果为:本发明通过基于迭代序贯参数凸逼近方法,在满足指定的系统服务质量QoS、发射机的发射功率以及回程链路消耗为约束条件基础上,使系统的和速率尽可能达到最大化;相比传统的时分双工模式及半正定松弛-块坐标下降方法,本发明所提方法具有更优秀的性能表现。
附图说明
图1为本发明的系统模型示意图。
图2为本发明的方法流程示意图。
图3为本发明的方法随着天线数变化时,系统与传统的云接入同时同频全双工模式以及时分双工模式的谱效提升对比示意图。
图4为本发明在不同自干扰抑制条件下,系统与现有的云接入同时同频全双工以及时分双工的性能曲线对比图示意图。
具体实施方式
如图1和图2所示,网络辅助全双工系统的稀疏波束设计与功率控制方法,包括如下步骤:
1).发送远程无线电射频头通过信道估计获得它到所有下行用户设备和接收远程无线电射频头间的理想信道状态信息,即信道矩阵
Figure BDA0002150745660000101
Figure BDA0002150745660000102
其中,
Figure BDA0002150745660000103
表示第l个发送远程无线电射频头到第k个下行用户设备的信道向量且服从复高斯分布
Figure BDA0002150745660000104
Figure BDA0002150745660000105
表示第l个发送远程无线电射频头到第z个接收远程无线电射频头的信道向量且服从复高斯分布。同理,上行用户通过信道估计获得它到所有下行用户设备和接收远程无线电射频头间的理想信道状态信息,即信道矩阵
Figure BDA0002150745660000106
Figure BDA0002150745660000107
其中,hIUI,j,k表示第j个上行用户到第k个下行用户设备的信道且服从复高斯分布
Figure BDA0002150745660000108
hU,j,z表示第j个上行用户到第z个接收远程无线电射频头的信道向量且服从复高斯分布
Figure BDA0002150745660000109
假设系统采用基于网络辅助全双工模式的时分双工制式,且信道服从平坦衰落,即在信道相干时间内信道系数保持不变;
2).中央处理器发送下行预编码向量到发送远程无线电射频头,发送远程无线电射频头接收该信号并将收到的信号,经过信道hD,k转发给下行用户设备,下行用户设备最终接收到的下行信息为
Figure BDA00021507456600001010
其中,sD,k,sU,j分别表示发送远程无线电射频头发送给下行用户设备k及上行用户j发送的数据信息,wD,k,pU,j分别表示发送远程无线电射频头发送给下行用户设备k的预编码向量及上行用户j的发射功率,
Figure BDA0002150745660000111
表示加性高斯白噪声。上行传输系统中,先对发送远程无线电射频头到接收远程无线电射频头的干扰功率做干扰消除处理,并对参与上行传输的用户及上行接收的远程无线电射频头做好配对。这样,接收远程无线电射频头接收到的上行用户信息为
Figure BDA0002150745660000112
其中,
Figure BDA0002150745660000113
表示干扰消除后的信道误差项,
Figure BDA0002150745660000114
表示加性高斯白噪声。
(3)整个传输系统以最大化和速率为准则,以指定的上行及下行用户设备传输最小系统频谱效率,上行用户的功率及下行远程天线的预编码发射功率,下行回程容量为约束,求解传输系统的稀疏波束形成与功率控制优化问题;
我们建立问题模型为:
Figure BDA0002150745660000115
C1:
Figure BDA0002150745660000116
C2:
Figure BDA0002150745660000117
C3:
Figure BDA0002150745660000118
C4:RD,k≥RD,min,k
C5:RU,j≥RU,min,j
C6:
Figure BDA0002150745660000119
其中,
Figure BDA00021507456600001110
Figure BDA00021507456600001111
Figure BDA0002150745660000121
Figure BDA0002150745660000122
Figure BDA0002150745660000123
uU,j,z为接收远程无线电射频头z为解调上行用户j发射的信号而设计的接收机。特别的,这里,接收远程无线电射频头z与上行用户采用预先配对模式,即接收远程无线电射频头z仅可解调与之配对的上行用户信号;其中,一个接收远程无线电射频头可以解调多个上行用户信号,而一个上行用户仅可以配对一个接收远程无线电射频头。
Figure BDA0002150745660000127
为指示函数,具体定义为:
Figure BDA0002150745660000124
CD,l,
Figure BDA0002150745660000125
RU,min,j,RD,min,k分别为预设的第l个发送远程无线电射频头的回程约束,发射功率约束,上行用户j的发射功率约束,最低速率要求,下行用户k的速率要求。
(4)由于指示函数为离散整数规划问题,无法直接通过优化求解处理,此处利用如下函数逼近该问题,即C7:
Figure BDA0002150745660000126
其中wD,l,k表示第l个发送远程无线电射频头到第k个下行用户设备的预编码向量,θ>>1为一个远大于1的实数。逼近后的函数由于是原问题的近似解,并不是原问题的准确解;因此,原问题可以分步骤完成;
(5)根据对数函数的单调性观察,首先将目标函数重新建模为:
Figure BDA0002150745660000131
C8:
Figure BDA0002150745660000132
C9:
Figure BDA0002150745660000133
C10:1+rD,k≥χD,k
C11:1+rU,j≥χU,j
C12:χU,j≥1,χD,k≥1,
C1,C2,C6,C7,
其中,χD,k,χU,j为针对下行用户设备速率以及上行用户设备速率而新引入的优化变量。
(6)在此基础上,特别地,将目标函数转化为二阶锥规划形式:
Figure BDA0002150745660000134
C13:
Figure BDA0002150745660000135
C14:
Figure BDA0002150745660000136
其中,
Figure BDA00021507456600001314
为新引入变量。
(7)根据一般的凸优化问题求解方法,我们将原问题的C7,C8,C9,C10,C11近似问题分别转化为:
C7→C15:
Figure BDA0002150745660000137
C8→C16:
Figure BDA0002150745660000138
C9→C17:
Figure BDA0002150745660000139
C10→C18:
Figure BDA00021507456600001310
C11→C19:
Figure BDA00021507456600001311
其中,
Figure BDA00021507456600001312
ωU,j为新引入的变量,
Figure BDA00021507456600001313
为上行用户j的干扰功率上界的二阶锥表示形式变量及一般表示形式变量,
Figure BDA0002150745660000141
Figure BDA0002150745660000142
具体转化过程,我们在步骤(8)中提供;
(8)①针对C7→C15:首先新引入变量
Figure BDA0002150745660000143
同时,利用一阶泰勒近似原理,得到如下结果:
C15,C20:
Figure BDA0002150745660000144
C21:
Figure BDA0002150745660000145
其中,
Figure BDA0002150745660000146
对C20,C21分别进行泰勒展开进行凸逼近如下:
C22:
Figure BDA0002150745660000147
C23:
Figure BDA0002150745660000148
其中,
Figure BDA0002150745660000149
②针对C8→C16:利用一阶泰勒近似原理将约束C8的等式右边展开,得到C16。
③针对C9→C17:首先通过引入优化变量
Figure BDA00021507456600001410
Figure BDA00021507456600001411
βU,j,j′
Figure BDA00021507456600001412
求出C9的近似处理表达式,如下:
C24:
Figure BDA0002150745660000151
C25:
Figure BDA0002150745660000152
C26:
Figure BDA0002150745660000153
C27:
Figure BDA0002150745660000154
C28:
Figure BDA0002150745660000155
C29:
Figure BDA0002150745660000156
通过进一步凸逼近,C25,C28,C29可以转换为:
C30:
Figure BDA0002150745660000157
C31:
Figure BDA0002150745660000158
C32:
Figure BDA0002150745660000159
这样,C9可以直接转换为二阶锥规划形式C17,其中,
Figure BDA00021507456600001510
④针对C10→C18:约束C10可以变换为:
C33:
Figure BDA00021507456600001511
观察到约束C15的不等式右边为凸形式,所以,根据一阶泰勒近似原理展开,可以表示为C18;
⑤针对C11→C19:首先,将约束C11变换为:
C34:
Figure BDA00021507456600001512
进一步不等式可以处理如下:
C19,C35:
Figure BDA00021507456600001513
其中,ωU,j,为新引进变量,
Figure BDA0002150745660000161
这里,利用序贯参数凸逼近方法,C35可以进一步近似为:
C36:
Figure BDA0002150745660000162
C37:
Figure BDA0002150745660000163
其中,
Figure BDA0002150745660000164
最终,C35可以近似为:
C39:
Figure BDA0002150745660000165
(9)由于以上步骤中采用部分近似处理,为了得到稳定解,迭代过程引入信任域约束方法,即:
C40:
Figure BDA0002150745660000166
C41:τ(PU,j)={|PU,j-PU,j|≤ηU,j}
τ(wD,l,k),τ(PU,j)表示信任域函数,ηD,k,l,ηU,j表示信任域阈值,当目标函数的相邻两次迭代结果的差大于阈值时,则ηD,k,l和ηU,j的值可以设为稍微大一点;反之,随着迭代步骤的进行,ηD,k,l和ηU,j的值可以相应的随之减小;
(10)根据步骤(5)到步骤(9),得到n+1次的迭代表达式如下:
Figure BDA0002150745660000167
s.t. C1,C2,C6,C12,C13,C14,C15,
C16,C17,C18,C19,C22,C23,C24,
C27,C30,C31,C32,C39,C40,C41
(11)给定下行远程天线单元与下行用户的服务配对模式,原问题可以转化为:
Figure BDA0002150745660000171
s.t. C1,C2,C6,C12,C13,C14,C16,
C17,C18,C19,C22,C23,C24,
C27,C30,C31,C32,C39,C40,C41
Figure BDA0002150745660000172
Figure BDA0002150745660000173
其中,
Figure BDA0002150745660000174
表示发射远程天线单元l与下行用户之间的关联集。
我们将整个实施方法总结为:
初始化
Figure BDA0002150745660000175
根据(10),迭代优化得到下行远程天线单元与下行用户的服务配对模式。
根据(11)迭代优化得到最终的实现结果:
(1)下行发射远程天线单元l与下行用户之间数据传输优化配对模式ψ(wD,l,k)。
(2)下行发射远程天线单元l的传输预编码优化设计wD,k
(3)上行用户的发射功率优化设计PU,j
图3给出了随着天线数变化时,本发明所提方案与传统的云接入同时同频全双工模式以及时分双工模式的谱效提升对比图。可以明显看到,本发明随着天线数增加,谱效有明显的提高,并且,谱效分别高于传统的云接入同时同频全双工模式以及时分双工模式5.44bps/Hz和21.45bps/Hz。图4给出了在不同上下行远程无线电射频头干扰抑制条件下,系统与现有时分双工的性能曲线对比图,明显看出,当抑制水平在-25DB以下时,算法性能明显得到提高,而传统的同时同频全双工模式,需要抑制在-30DB以下。
为了使已有的双向动态网络和全双工方案更一般化,本发明提出了一种网络辅助全双工体系结构,使远程无线电射频头能够更灵活地运行双工模式,即远程无线电射频头可以根据实际需要在双向动态网络和全双工方案模式中运行。
本发明的网络辅助全双工系统的稀疏波束设计与功率控制方法,以最大化系统和速率为目标,以满足指定的系统服务质量QoS、发射机的发射功率以及回程链路消耗为约束条件,建立了以远程无线电射频头发射机与接收机波束,以及上行用户功率为设计变量的数学优化模型。本方法通过基于迭代序贯参数凸逼近方法,将原始的非凸优化问题通过不断连续凸逼近方法逐步转化为凸优化问题,并通过不断迭代优化联合求出稀疏波束形成向量与发射功率的值。

Claims (6)

1.网络辅助全双工系统的稀疏波束设计与功率控制方法,其特征在于,包括如下步骤:
(1)发送远程无线电射频头通过信道估计获得它到所有下行用户设备和接收远程无线电射频头间的理想信道状态信息,上行用户通过信道估计获得它到所有下行用户设备和接收远程无线电射频头间的理想信道状态信息;假设系统采用基于网络辅助全双工模式的时分双工制式,且信道服从平坦衰落,即在信道相干时间内信道系数保持不变;
(2)中央处理器发送下行预编码向量到发送远程无线电射频头,发送远程无线电射频头接收信号并将收到的信号,经过信道hD,k转发给下行用户设备,下行用户设备最终接收到的下行信息为
Figure FDA0003032278300000011
其中,
Figure FDA0003032278300000012
表示发送远程无线电射频头到所有下行用户设备的理想信道状态信息,hIUI,j,k表示第j个上行用户到第k个下行用户设备的信道且服从复高斯分布
Figure FDA0003032278300000013
Figure FDA0003032278300000014
sD,k,sU,j分别表示发送远程无线电射频头发送给下行用户设备k及上行用户j发送的数据信息;wD,k,pU,j分别表示发送远程无线电射频头发送给下行用户设备k的预编码向量及上行用户设备j的发射功率;
Figure FDA0003032278300000015
表示加性高斯白噪声;上行传输系统中,先对发送远程无线电射频头到接收远程无线电射频头的干扰功率做干扰消除处理,并对参与上行传输的用户及上行接收的远程无线电射频头做好配对;这样,接收远程无线电射频头接收到的上行用户信息为
Figure FDA0003032278300000016
其中,
Figure FDA0003032278300000017
表示干扰消除后的信道估计误差项,
Figure FDA0003032278300000018
表示加性高斯白噪声,
Figure FDA0003032278300000019
表示上行用户通过信道估计到接收远程无线电射频头间的理想信道状态信息;
(3)整个传输系统以最大化和速率为准则,以指定的上行及下行用户设备传输最小系统频谱效率,上行用户的功率及下行远程天线的预编码发射功率,下行回程容量为约束,求解传输系统的稀疏波束形成与功率控制优化问题;
建立问题模型为:
Figure FDA0003032278300000021
C1:
Figure FDA0003032278300000022
C2:
Figure FDA0003032278300000023
C3:
Figure FDA0003032278300000024
C4:RD,k≥RD,min,k
C5:RU,j≥RU,min,j
C6:
Figure FDA0003032278300000025
其中,
Figure FDA0003032278300000026
Figure FDA0003032278300000027
Figure FDA0003032278300000028
Figure FDA0003032278300000029
Figure FDA00030322783000000210
uU,j,z为接收远程无线电射频头z为解调上行用户j发射的信号而设计的接收机,接收远程无线电射频头z与上行用户采用预先配对模式,即接收远程无线电射频头z仅解调与之配对的上行用户信号;其中,一个接收远程无线电射频头解调多个上行用户信号,而一个上行用户仅配对一个接收远程无线电射频头,CD,l,
Figure FDA00030322783000000211
RU,min,j,RD,min,k分别为预设的第l个发送远程无线电射频头的回程约束,发射功率约束,上行用户j的发射功率约束,最低速率要求,下行用户k的速率要求,
Figure FDA00030322783000000212
表示发送远程无线电射频头到所有下行用户设备的理想信道状态信息,hIUI,j,k表示第j个上行用户到第k个下行用户设备的信道且服从复高斯分布
Figure FDA0003032278300000031
Figure FDA0003032278300000032
sD,k,sU,j分别表示发送远程无线电射频头发送给下行用户设备k及上行用户j发送的数据信息;wD,k,pU,j分别表示发送远程无线电射频头发送给下行用户设备k的预编码向量及上行用户设备j的发射功率;
Figure FDA0003032278300000033
表示加性高斯白噪声,
Figure FDA0003032278300000034
表示干扰消除后的信道估计误差项,
Figure FDA0003032278300000035
表示加性高斯白噪声,
Figure FDA0003032278300000036
表示上行用户通过信道估计到接收远程无线电射频头间的理想信道状态信息;
Figure FDA0003032278300000037
为指示函数,具体定义为:
Figure FDA0003032278300000038
(4)由于指示函数为离散整数规划问题,无法直接通过优化求解处理,此处利用如下函数逼近该问题,即C7:
Figure FDA0003032278300000039
其中wD,l,k表示第l个发送远程无线电射频头到第k个下行用户设备的预编码向量,θ>>1为一个远大于1的实数;逼近后的函数由于是原问题的近似解,并不是原问题的准确解;因此,原问题分步骤完成;
(5)根据对数函数的单调性观察,首先将目标函数重新建模为:
Figure FDA00030322783000000310
C8:
Figure FDA00030322783000000311
C9:
Figure FDA00030322783000000312
C10:1+rD,k≥χD,k,
C11:1+rU,j≥χU,j,
C12:χU,j≥1,χD,k≥1,
C1,C2,C6,C7,
其中,χD,k,χU,j为针对下行用户设备速率以及上行用户设备速率而新引入的优化变量;
(6)在此基础上,将目标函数转化为二阶锥规划形式:
Figure FDA0003032278300000041
C13:
Figure FDA0003032278300000042
C14:
Figure FDA0003032278300000043
其中,
Figure FDA0003032278300000044
为新引入变量;
(7)根据一般的凸优化问题求解方法,利用迭代序贯参数凸逼近方法将原问题转化成凸问题求解,即步骤S1;
(8)判断原问题相邻两个解之间的差值是否大于阈值,即步骤S2;
(9)根据步骤(5)到步骤(8),重复步骤S1、步骤S2并返回下行远程天线单元与下行用户服务配对模式,即步骤S3;
(10)给定下行远程天线单元与下行用户的服务配对模式,利用迭代序贯参数凸逼近及块坐标下降方法求解原问题并重复步骤S1和步骤S2,即步骤S4。
2.如权利要求1所述的网络辅助全双工系统的稀疏波束设计与功率控制方法,其特征在于,步骤(1)中,发送远程无线电射频头通过信道估计获得它到所有下行用户设备和接收远程无线电射频头间的理想信道状态信息,即
Figure FDA0003032278300000045
其中,
Figure FDA0003032278300000046
表示第l个发送远程无线电射频头到第k个下行用户设备的信道向量且服从复高斯分布
Figure FDA0003032278300000047
Figure FDA0003032278300000048
表示第l个发送远程无线电射频头到第z个接收远程无线电射频头的信道向量且服从复高斯分布;hU,j,z表示第j个上行用户到第z个接收远程无线电射频头的信道向量且服从复高斯分布
Figure FDA0003032278300000049
3.如权利要求1所述的网络辅助全双工系统的稀疏波束设计与功率控制方法,其特征在于,将原问题的C7、C8、C9、C10、C11近似问题分别转化为:
C7→C15:
Figure FDA00030322783000000410
C8→C16:
Figure FDA0003032278300000051
C9→C17:
Figure FDA0003032278300000052
C10→C18:
Figure FDA0003032278300000053
C11→C19:
Figure FDA0003032278300000054
其中,
Figure FDA0003032278300000055
Figure FDA0003032278300000056
ωU,j为新引入的变量,
Figure FDA0003032278300000057
为上行用户j的干扰功率上界的二阶锥表示形式变量及一般表示形式变量,
Figure FDA0003032278300000058
Figure FDA0003032278300000059
(7.1)①针对C7→C15:首先新引入变量
Figure FDA00030322783000000510
Figure FDA00030322783000000511
同时,利用一阶泰勒近似原理,得到如下结果:
Figure FDA00030322783000000512
其中,
Figure FDA00030322783000000513
对C20,C21分别进行泰勒展开进行凸逼近如下:
Figure FDA00030322783000000514
其中,
Figure FDA00030322783000000515
②针对C8→C16:利用一阶泰勒近似原理将约束C8的等式右边展开,得到C16;
③针对C9→C17:首先通过引入优化变量
Figure FDA0003032278300000061
Figure FDA0003032278300000062
βU,j,j′
Figure FDA0003032278300000063
求出C9的近似处理表达式,如下:
C24:
Figure FDA0003032278300000064
C25:
Figure FDA0003032278300000065
C26:
Figure FDA0003032278300000066
C27:
Figure FDA0003032278300000067
C28:
Figure FDA0003032278300000068
C29:
Figure FDA0003032278300000069
通过进一步凸逼近,C25,C28,C29转换为:
C30:
Figure FDA00030322783000000610
C31:
Figure FDA00030322783000000611
C32:
Figure FDA00030322783000000612
这样,C9直接转换为二阶锥规划形式C17,其中,
Figure FDA00030322783000000613
④针对C10→C18:约束C10变换为:
C33:
Figure FDA00030322783000000614
观察到约束C15的不等式右边为凸形式,所以,根据一阶泰勒近似原理展开,表示为C18;
⑤针对C11→C19:首先,将约束C11变换为:
C34:
Figure FDA00030322783000000615
进一步不等式处理如下:
C19,C35:
Figure FDA0003032278300000071
其中,ωU,j,为新引进变量,
Figure FDA0003032278300000072
这里,利用序贯参数凸逼近方法,C35进一步近似为:
C36:
Figure FDA0003032278300000073
C37:
Figure FDA0003032278300000074
其中,
Figure FDA0003032278300000075
最终,C35近似为:
C39:
Figure FDA0003032278300000076
4.如权利要求3所述的网络辅助全双工系统的稀疏波束设计与功率控制方法,其特征在于,迭代过程引入信任域约束方法,即:
Figure FDA0003032278300000077
τ(wD,l,k),τ(PU,j)表示信任域函数,ηD,k,l,ηU,j表示信任域阈值。
5.如权利要求4所述的网络辅助全双工系统的稀疏波束设计与功率控制方法,其特征在于,得到n+1次的迭代表达式如下:
Figure FDA0003032278300000078
6.如权利要求4所述的网络辅助全双工系统的稀疏波束设计与功率控制方法,其特征在于,给定下行远程天线单元与下行用户的服务配对模式,原问题转化为:
Figure FDA0003032278300000081
s.t.C1,C2,C6,C12,C13,C14,C16,
C17,C18,C19,C22,C23,C24,
C27,C30,C31,C32,C39,C40,C41
Figure FDA0003032278300000082
Figure FDA0003032278300000083
其中,
Figure FDA0003032278300000084
表示发射远程天线单元l与下行用户之间的关联集。
CN201910700869.XA 2019-07-31 2019-07-31 网络辅助全双工系统的稀疏波束设计与功率控制方法 Active CN110417445B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910700869.XA CN110417445B (zh) 2019-07-31 2019-07-31 网络辅助全双工系统的稀疏波束设计与功率控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910700869.XA CN110417445B (zh) 2019-07-31 2019-07-31 网络辅助全双工系统的稀疏波束设计与功率控制方法

Publications (2)

Publication Number Publication Date
CN110417445A CN110417445A (zh) 2019-11-05
CN110417445B true CN110417445B (zh) 2021-06-11

Family

ID=68364436

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910700869.XA Active CN110417445B (zh) 2019-07-31 2019-07-31 网络辅助全双工系统的稀疏波束设计与功率控制方法

Country Status (1)

Country Link
CN (1) CN110417445B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112188605B (zh) * 2020-09-29 2022-11-01 东南大学 无线携能通信场景下网络辅助全双工系统收发机设计方法
CN114885423A (zh) * 2022-03-29 2022-08-09 网络通信与安全紫金山实验室 一种网络辅助全双工系统能效优化方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104993859A (zh) * 2015-08-05 2015-10-21 中国电子科技集团公司第五十四研究所 一种适用于时间异步环境下的分布式波束赋形方法
CN105191163A (zh) * 2013-03-28 2015-12-23 华为技术有限公司 用于稀疏波束形成设计的系统和方法
CN108718228A (zh) * 2018-05-10 2018-10-30 东南大学 网络辅助的灵活双工无线通信系统及其方法
CN109818662A (zh) * 2019-03-12 2019-05-28 电子科技大学 全双工云接入数能一体化网络中混合波束成形方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9407302B2 (en) * 2012-12-03 2016-08-02 Intel Corporation Communication device, mobile terminal, method for requesting information and method for providing information
US10367677B2 (en) * 2016-05-13 2019-07-30 Telefonaktiebolaget Lm Ericsson (Publ) Network architecture, methods, and devices for a wireless communications network

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105191163A (zh) * 2013-03-28 2015-12-23 华为技术有限公司 用于稀疏波束形成设计的系统和方法
CN104993859A (zh) * 2015-08-05 2015-10-21 中国电子科技集团公司第五十四研究所 一种适用于时间异步环境下的分布式波束赋形方法
CN108718228A (zh) * 2018-05-10 2018-10-30 东南大学 网络辅助的灵活双工无线通信系统及其方法
CN109818662A (zh) * 2019-03-12 2019-05-28 电子科技大学 全双工云接入数能一体化网络中混合波束成形方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Sparse Beamforming for an Ultradensely Distributed Antenna System With Interlaced Clustering;Xinjiang Xia 等;《IEEE Access》;20190131;第7卷;第15069-15085页 *

Also Published As

Publication number Publication date
CN110417445A (zh) 2019-11-05

Similar Documents

Publication Publication Date Title
CN111314935B (zh) 基于noma-mec系统的下行传输时延最小化的方法
CN111970035A (zh) 一种mimo全双工双向安全通信系统预编码矩阵的优化方法
CN110417445B (zh) 网络辅助全双工系统的稀疏波束设计与功率控制方法
CN110149127B (zh) 一种基于noma技术的d2d通信系统预编码向量优化方法
CN113825159B (zh) 基于智能反射面的无线携能通信系统鲁棒资源分配方法
CN110299934A (zh) 一种无线携能全双工中继系统的保密传输方法
WO2023165183A1 (zh) 一种多输入多输出收发机设计方法及装置
CN109348500B (zh) 一种硬件损伤条件下满足双向swipt中继系统的资源分配方法
CN109068382B (zh) 一种基于时延QoS的NOMA跨层功率分配方法
CN116155327A (zh) 一种低精度双向无蜂窝大规模mimo系统功率分配方法
CN114885423A (zh) 一种网络辅助全双工系统能效优化方法及系统
CN110191476B (zh) 一种基于可重构天线阵列的非正交多址接入方法
CN116760448A (zh) 一种基于mimo-noma的星地融合网络资源高效分配方法
CN108832978B (zh) 一种包含直传链路的多用户mimo中继系统的联合预编码方法
CN112188605B (zh) 无线携能通信场景下网络辅助全双工系统收发机设计方法
CN108449737B (zh) 分布式天线系统中基于d2d的下行链路高能效功率分配方法
CN107332614A (zh) 可见光通信非正交多址接入技术稳健波束成形的优化方法
CN112332882B (zh) 一种基于毫米波全双工中继通信的鲁棒混合收发机设计方法
CN108923831B (zh) 一种发射信号的预编码方法和装置
TW202145005A (zh) 多輸入輸出系統之參數估測方法
CN107517088B (zh) 星地一体化系统基于混合信道衰落的干扰分析方法
CN108540192B (zh) 中断概率约束下的全双工多天线中继发射波束赋形方法
CN110601736A (zh) 一种多天线全双工认知无线电能量捕获与信息传输方法
CN116321236A (zh) Ris辅助的安全无蜂窝大规模mimo系统能效优化方法
CN112929062B (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