CN114978262A - 一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法 - Google Patents
一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法 Download PDFInfo
- Publication number
- CN114978262A CN114978262A CN202210538796.0A CN202210538796A CN114978262A CN 114978262 A CN114978262 A CN 114978262A CN 202210538796 A CN202210538796 A CN 202210538796A CN 114978262 A CN114978262 A CN 114978262A
- Authority
- CN
- China
- Prior art keywords
- vector
- user
- precoding
- formula
- beam satellite
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000004891 communication Methods 0.000 title claims abstract description 27
- 239000013598 vector Substances 0.000 claims abstract description 158
- 238000005457 optimization Methods 0.000 claims abstract description 62
- 239000011159 matrix material Substances 0.000 claims abstract description 47
- 238000013461 design Methods 0.000 claims abstract description 19
- 230000005540 biological transmission Effects 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000009826 distribution Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000013139 quantization Methods 0.000 claims description 9
- 238000009827 uniform distribution Methods 0.000 claims description 9
- 230000009977 dual effect Effects 0.000 claims description 8
- 230000014509 gene expression Effects 0.000 claims description 8
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 6
- 230000004907 flux Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 4
- 238000012804 iterative process Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18513—Transmission in a satellite or space-based system
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Abstract
本发明公开了一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,包括:构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型;获取关于信道自相关矩阵的统计信道模型;构建多波束卫星系统和速率最大化的鲁棒混合预编码优化设计问题;在部分连接结构下,分别将多波束卫星系统和速率最大化的鲁棒数字预编码和模拟预编码的优化设计问题等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题;结合对应的等价优化问题的拉格朗日函数及其KKT条件,通过迭代算法获得最优的数字预编码和模拟预编码矢量;本发明能够显著提高系统的发射功率效率,降低了硬件成本和运行成本,同时提高了系统对定位角度估计误差的鲁棒性。
Description
技术领域
本发明涉及多波束卫星通信鲁棒混合预编码方法的技术领域,特别是涉及一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法。
背景技术
近年来,人们对更高质量的通信服务和全球无缝接入的需求迅速增加,特别是在后5G及未来的6G网络,这使得从太空提供互联网服务的概念引起了学术界和工业界的深刻关注。值得注意的是,卫星通信系统在提供无缝覆盖和扩展地面通信网络方面具有优势,并可实现更低的延迟、更低的生产成本、更大的覆盖范围,提供无处不在的全球无线接入。但是由于目前低轨卫星(LEO)趋于小型化、低生产成本和有限的收发器功率,硬件成本和功耗成为卫星预编码器设计的关键考虑因素。
近几年,具有高能效和低实现成本的混合预编码技术得到了研究界的广泛关注。在混合预编码结构中,混合收发器通常将大量低成本的模拟移相器(PS)和功率放大器与低维数字信号处理单元相结合,使用少数昂贵的射频链路(RF),这种混合预编码器的设计降低了系统生产成本和功率消耗,符合当前高能效LEO卫星通信系统的设计需求。此外,在毫米波频段,由于发送天线尺寸小,因此可以将许多天线阵元封装在一个紧凑的区域内,有利于发送波形在空间域上有着良好的方向自由度和三维覆盖范围,并在具有严重路径损耗的毫米波信道提供足够的预编码增益。然而,毫米波卫星下行链路的混合预编码器的设计并不是一件容易的事,因为模拟预编码器和数字预编码器的优化问题通常是非凸的,求解最优的预编码矢量比较困难。
发明内容
有鉴于此,本发明的目的在于提供一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,以提高多波束卫星通信系统的传输效率,并降低硬件成本和系统功耗。
为了实现上述目的,本发明采用如下技术方案:
一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,所述方法包括:
步骤S1、基于多波束卫星对各用户的位置角度估计误差以及由于卫星姿态轨道控制而导致的公共角度误差,构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型;再基于该多波束卫星下行信道矢量模型构建用户侧的接收信号模型;
步骤S2、对多波束卫星系统的和速率关于定位角度误差变量求数学期望,得到第k个用户的遍历传输速率,获取关于信道自相关矩阵的统计信道模型;
步骤S3、基于部分连接的混合预编码结构,关于模拟预编码矢量和数字预编码矢量构建多波束卫星系统和速率最大化的鲁棒混合预编码优化设计问题,约束条件为多波束卫星的每个天线阵元的功率值均小于某一门限值,其中,该优化设计问题定义为第一优化问题;
步骤S4、将步骤S3中构建的第一优化问题等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题,其中,该功率最小化问题,定义为第二优化问题;
步骤S5、固定模拟预编码矩阵,结合第二优化问题的拉格朗日函数及其KKT条件,将混合预编码结构中数字预编码的最优结构建模为一个广义特征值问题的解,计算最优数字预编码矢量;
步骤S6、固定最优数字预编码矢量,结合关于模拟预编码矢量的等价转化问题的拉格朗日函数及其对应的KKT条件,将模拟预编码的最优方向建模为一个广义特征值问题的解,计算获得最优模拟预编码矢量;其中,该等价转化后的关于模拟预编码矢量的设计问题定义为第三优化问题;
步骤S7、对模拟预编码矢量的相位进行离散化操作,同时保证多波束卫星系统和速率最大化;
步骤S8、利用次梯度投影技术及拉格朗日函数与其对偶函数的强对偶性确定拉格朗日乘子;
步骤S9、循环迭代步骤S4、步骤S5、步骤S6、步骤S7和步骤S8直至系统和速率性能收敛,获得最优数字预编码矢量和具有离散相位的模拟预编码矢量。
进一步的,在所述步骤S1中,假设卫星对用户存在角度定位误差,引入定位角度误差变量后,多波束卫星到第k个用户的下行信道参数矢量hk为:
在公式中,t表示时间,f表示系统频率,vk表示平面阵天线相应矢量,表示由于多波束卫星移动导致的多普勒频偏,τk,l表示多波束卫星到第k个用户的第l条传播径的传播延时,表示传播延时τk,l的最小值,即表示第k个用户的下行信道增益,且满足:
进一步的,在所述步骤S1中,构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型,其还包括构建如下角度关系:
在公式中,θk为用户k相对于低轨道卫星平面阵天线y轴的准确角度,为用户k相对于低轨道卫星平面阵天线x轴的准确角度,和表示卫星对第k个用户的位置角度估计误差,和表示由于卫星姿态轨道控制而导致的公共角度检测误差,表示为用户k相对于低轨道卫星平面阵天线的y轴的估计角度,表示为用户k相对于低轨道卫星平面阵天线x轴的估计角度;
进一步的,所述步骤S2具体包括:
公式中,SINRk表示第k个用户的信干噪比,并且定义为:
公式中,平面阵天线相应矢量vk的第n项表示为:
其中,信道自相关矩阵Rk中的第m行第n列元素[Rk]m,n表示为:
其中,a和b分别为第一中间变量和第二中间变量;
则得到:
其中各中间变量分别定义为:
公式中,A、B和Z分别为第三中间变量、第四中间变量和第五中间变量;
其中,各个中间变量定义为:
其中,P、Q、D、F、C、E分别为第六中间变量、第七中间变量、第八中间变量、第九中间变量、第十中间变量和第十一中间变量。
进一步的,所述步骤S3包括:
P1:
因此,第一优化问题表述如下:
在公式中,Pn是第n个天线的最大可用功率。
进一步的,所述步骤S4具体包括:
通过下面等价以使QoS约束更容易处理:
因此优化问题P2重新表述为第二优化问题:
进一步的,所述步骤S5具体包括:
固定模拟预编码矩阵WRF,则第二优化问题P2相对于数字预编码矢量wD,k的拉格朗日函数表示为:
根据KKT条件,最优数字预编码矢量需要满足:
因此得到如下关于数字预编码矢量的广义特征值问题:
再令
γk看作矩阵对(SD,k,ND,k)的最大广义特征值,uD,k是对应于γk的广义特征向量,即,
γk=max.generalized eigenvector(SD,k,ND,k),
uD,k=max.generalized eigenvector(SD,k,ND,k);
因此一旦确定了γk和uD,k,则能够计算第k个用户的数字预编码的功率值ρD,k;
定义:
数字预编码的功率矢量ρD=[ρD,1,…,ρD,K]T计算为:
因此,最优数字预编码矢量wD,k能够通过所提出的解结构精确计算。
进一步的,所述步骤S6具体包括:
固定步骤S5得到的最优数字预编码矢量wD,k,为了方便对模拟预编码矢量处理,首先需要将用户k的混合预编码矢量WRFwD,k转换为以下等价形式:
因此,模拟预编码矩阵WRF的设计简化为对矢量wA的设计;
在获得最佳模拟预编码矢量wA后,利用一个量化操作来保证该离散相位约束;
在给定数字预编码矢量wD,k以及拉格朗日乘数vk和qn的前提下,等价转化后的关于模拟预编码矢量的设计问题表述为:
上述的第三优化问题P3的拉格朗日函数表示为
根据拉格朗日函数的KKT条件,最优模拟预编码矢量wA满足如下广义特征值问题:
在部分连接的结构中,每个RF链路都连接到一个天线子集,其中对应于第j个RF链路的模拟预编码矢量wA,j被认为是独立于其他射频链路的,因此对第i个RF链路进行广义特征值分解,而不是直接计算wA;
令
对应于第i个RF链路的最优模拟预编码矢量由下式给出:
进一步的,所述步骤S7具体包括:
根据等价优化问题P2的拉格朗日函数和其对偶函数,以及两者之间的强对偶关系,得到:
在该公式中,Ptotal是系统允许的最大发射功率;
在实际应用中,目标SINR的约束门限值γk需要尽可能地取最大值,因此求解如下凸优化问题确定vk的值:
在所述步骤S9中,在发送天线单天线功率约束下,通过采用一个简单的终止策略来停止单天线功率约束的迭代过程,即:
在该公式中,εn表示最大可超过第n个天线阵元额定功率的功率值;
当系统的和速率性能收敛时,即得到最优模拟预编码器和最优数字预编码器。
本发明的有益效果是:
1、本发明通过对系统和速率求数学期望并结合对单天线阵元的功率约束,有效减少多波束卫星用户定位角度不准确性带来的不利影响。此外,利用具有低分辨率相移器的混合预编码结构,能够在保证多波束卫星传输速率的同时,提高卫星的发射功率效率,降低硬件实现成本;
2、本发明将卫星对用户定位时存在的角度不确定性以及由于卫星姿态抖动而导致的检测角度误差融合至信道建模中,并且根据角度误差变量服从某一概率分布这一特点,建立了能够实现混合预编码性能对定位角度误差鲁棒性的统计信道模型;
3、本发明为多波束低轨卫星下行通信场景建立了一种搭载低分辨率移相器的部分连接混合预编码结构。为了保证每个天线的功率放大器工作在线性区域内,避免信号失真,本发明采用了一种更加实际的单天线功率约束。
4、本发明在单天线阵元的功率约束下,实现卫星通信系统的遍历和速率最大化,分别建立关于数字预编码矢量和模拟预编码矢量的优化问题,并根据各自拉格朗日函数及其KKT条件,将最优鲁棒数字预编码和模拟预编码的结构建模为对应广义特征值问题的解。同时,本发明开发出了一种高效的迭代算法求解最优数字预编码矢量和模拟预编码矢量,进而获得最优的混合预编码器分连接结构下的混合预编码器的示意图。
附图说明
图1为本实施例1中提出的一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法的流程示意图;
图2为本实施例1中卫星天线面阵对用户存在定位角度误差的示意图;
图3为本实施例1采用的部分连接结构下的混合预编码器的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
参见图1-图3,本实施例提供一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,该方法的具体流程如图1所示,该方法具体包括如下的步骤:
步骤S1、基于多波束卫星对各用户的位置角度估计误差以及由于卫星姿态轨道控制而导致的公共角度误差,构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型;再基于该多波束卫星下行信道矢量模型构建用户侧的接收信号模型;
具体的说,在本实施例中,该步骤S1包括:
针对一多波束卫星下行传输场景,在该场景中,LEO卫星同时为K单天线用户提供服务;且在卫星侧配备均匀平面阵列(UPA),由NT=MxMy天线组成,其中Mx和My是在分别为x和y轴上的天线数。射频(RF)链路数为NRF大于或等于用户数,即NRF≥K,每条射频链路与带有天线的子天线阵列连接。
如图2所示,假设卫星对用户存在角度定位误差,引入定位角度误差变量后,多波束卫星到第k个用户的下行信道参数矢量hk为:
在公式中,t表示时间,f表示系统频率,vk表示平面阵天线相应矢量,表示由于多波束卫星移动导致的多普勒频偏,τk,l表示多波束卫星到第k个用户的第l条传播径的传播延时,表示传播延时τk,l的最小值,即表示第k个用户的下行信道增益,且满足:
具体的说,在本实施例中,构建多波束卫星下行信道矢量模型还包括:
构建如下角度关系:
在公式中,θk为用户k相对于低轨道卫星平面阵天线y轴的准确角度,为用户k相对于低轨道卫星平面阵天线x轴的准确角度,和表示卫星对第k个用户的位置角度估计误差,和表示由于卫星姿态轨道控制而导致的公共角度检测误差,表示为用户k相对于低轨道卫星平面阵天线的y轴的估计角度,表示为用户k相对于低轨道卫星平面阵天线x轴的估计角度。
由于模拟预编码器配备了移相器和功率放大器,因此模拟预编码器可以同时对输入的信号进行幅度和相位的变化。本实施例定义了模拟预编码器的第(i,j)个元素其中,连续相位θi,j在[0,2π]区间内取值,ai,j为信号的幅值。在实际卫星系统中,θi,j通常通过离散相位来近似以降低硬件成本,因此可以通过一组低分辨率相移器实现,并且离散相位的取值属于集合其中B是控制量化相位的精度的比特数。
具体的说,在本实施例中,定义sk为从LEO卫星发送到第k个用户的传输信号,那么第k个用户的接收信号可表示为
在部分连接架构中,模拟预编码权值矩阵WRF可以表示为
步骤S2、对多波束卫星系统的和速率关于定位角度误差变量求数学期望,得到第k个用户的遍历传输速率,获取关于信道自相关矩阵的统计信道模型;
具体的说,在本实施例中,该步骤S2包括:
公式中,SINRk表示第k个用户的信干噪比,并且定义为:
公式中,平面阵天线相应矢量vk的第n项可以表示为:
其中,信道自相关矩阵Rk中的第m行第n列元素[Rk]m,n可以表示为:
其中,a和b分别为第一中间变量和第二中间变量。
则可以得到:
其中各中间变量分别定义为:
公式中,A、B和Z分别为第三中间变量、第四中间变量和第五中间变量;
其中,各个中间变量定义为:
其中,P、Q、D、F、C、E分别为第六中间变量、第七中间变量、第八中间变量、第九中间变量、第十中间变量和第十一中间变量。
步骤S3、基于部分连接的混合预编码结构,关于模拟预编码矢量和数字预编码矢量构建多波束卫星系统和速率最大化的鲁棒混合预编码优化设计问题,约束条件为多波束卫星的每个天线阵元的功率值均小于某一门限值,其中,该优化设计问题定义为第一优化问题;
因此,第一优化问题可以表述如下:
在公式中,Pn是第n个天线的最大可用功率。
步骤S4、将步骤S3中构建的第一优化问题等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题,其中,该功率最小化问题,定义为第二优化问题;
进一步,通过下面等价以使QoS约束更容易处理:
步骤S5、固定模拟预编码矩阵,结合第二优化问题的拉格朗日函数及其KKT条件,将混合预编码结构中数字预编码的最优结构建模为一个广义特征值问题的解,计算最优数字预编码矢量;
具体的说,在本实施例中,固定模拟预编码矩阵WRF,则第二优化问题P2相对于数字预编码矢量wD,k的拉格朗日函数可以表示为:
根据KKT条件,最优数字预编码矢量需要满足:
因此可以得到如下关于数字预编码矢量的广义特征值问题:
再令
将最优数字预编码矢量表示为其中用户k的数字预编码矢量的功率ρD,k和单位范数数字预编码矢量的方向uD,k可以通过对矩阵SD,k和ND,k进行广义特征值分解来计算。特别地,γk可以看作矩阵对(SD,k,ND,k)的最大广义特征值,uD,k是对应于γk的广义特征向量,即,
γk=max.generalized eigenvector(SD,k,ND,k),
uD,k=max.generalized eigenvector(SD,k,ND,k);
因此一旦确定了γk和uD,k,则可计算第k个用户的数字预编码的功率值ρD,k。
定义:
数字预编码的功率矢量ρD=[ρD,1,…,ρD,K]T可以计算为:
因此,最优数字预编码矢量wD,k可以通过所提出的解结构精确计算。
步骤S6、固定最优数字预编码矢量,结合关于模拟预编码矢量的等价转化问题的拉格朗日函数及其对应的KKT条件,将模拟预编码的最优方向建模为一个广义特征值问题的解,计算获得最优模拟预编码矢量;其中,该等价转化后的关于模拟预编码矢量的设计问题定义为第三优化问题;
具体的说,在本实施例中,固定步骤S5得到的最优数字预编码矢量wD,k,为了方便对模拟预编码矢量处理,首先需要将用户k的混合预编码矢量WRFwD,k转换为以下等价形式:
因此,模拟预编码矩阵WRF的设计可以简化为对矢量wA的设计。
在给定数字预编码矢量wD,k以及拉格朗日乘数vk和qn的前提下,等价转化后的关于模拟预编码矢量的设计问题表述为:
进一步的,问题P3的拉格朗日函数可以表示为
根据拉格朗日函数的KKT条件,最优模拟预编码矢量wA满足如下广义特征值问题:
在部分连接的结构中,每个RF链路都连接到一个天线子集,其中对应于第j个RF链路的模拟预编码矢量wA,j可以被认为是独立于其他射频链路的,因此对第i个RF链路进行广义特征值分解,而不是直接计算wA。令
对应于第i个RF链路的最优模拟预编码矢量由下式给出:
步骤S7、对模拟预编码矢量的相位进行离散化操作,同时保证多波束卫星系统和速率最大化;
步骤S8、利用次梯度投影技术及拉格朗日函数与其对偶函数的强对偶性确定拉格朗日乘子;具体的说,在本实施例中,为了确定拉格朗日乘子和的最优值,采用次梯度投影技术并结合KKT条件qn≥0,和qnPn=0,得到关于的迭代表达式Q为:
根据等价优化问题P2的拉格朗日函数和其对偶函数,以及两者之间的强对偶关系,可以得到:
在该公式中,Ptotal是系统允许的最大发射功率。
在实际应用中,目标SINR的约束门限值γk需要尽可能地取最大值,因此求解如下凸优化问题确定vk的值:
步骤S9、循环迭代步骤S4、步骤S5、步骤S6、步骤S7和步骤S8直至系统和速率性能收敛,获得最优数字预编码矢量和具有离散相位的模拟预编码矢量。
具体的说,在本实施例中,在发送天线单天线功率约束下,可以采用一个简单的终止策略来停止单天线功率约束的迭代过程,即:
在该公式中,εn表示最大可超过第n个天线阵元额定功率的功率值。当系统的和速率性能收敛时,即得到最优模拟预编码器和最优数字预编码器。
本发明未详述之处,均为本领域技术人员的公知技术。
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。
Claims (10)
1.一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,其特征在于,所述方法包括:
步骤S1、基于多波束卫星对各用户的位置角度估计误差以及由于卫星姿态轨道控制而导致的公共角度误差,构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型;再基于该多波束卫星下行信道矢量模型构建用户侧的接收信号模型;
步骤S2、对多波束卫星系统的和速率关于定位角度误差变量求数学期望,得到第k个用户的遍历传输速率,获取关于信道自相关矩阵的统计信道模型;
步骤S3、基于部分连接的混合预编码结构,关于模拟预编码矢量和数字预编码矢量构建多波束卫星系统和速率最大化的鲁棒混合预编码优化设计问题,约束条件为多波束卫星的每个天线阵元的功率值均小于某一门限值,其中,该优化设计问题定义为第一优化问题;
步骤S4、将步骤S3中构建的第一优化问题等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题,其中,该功率最小化问题,定义为第二优化问题;
步骤S5、固定模拟预编码矩阵,结合第二优化问题的拉格朗日函数及其KKT条件,将混合预编码结构中数字预编码的最优结构建模为一个广义特征值问题的解,计算最优数字预编码矢量;
步骤S6、固定最优数字预编码矢量,结合关于模拟预编码矢量的等价转化问题的拉格朗日函数及其对应的KKT条件,将模拟预编码的最优方向建模为一个广义特征值问题的解,计算获得最优模拟预编码矢量;其中,该等价转化后的关于模拟预编码矢量的设计问题定义为第三优化问题;
步骤S7、对模拟预编码矢量的相位进行离散化操作,同时保证多波束卫星系统和速率最大化;
步骤S8、利用次梯度投影技术及拉格朗日函数与其对偶函数的强对偶性确定拉格朗日乘子;
步骤S9、循环迭代步骤S4、步骤S5、步骤S6、步骤S7和步骤S8直至系统和速率性能收敛,获得最优数字预编码矢量和具有离散相位的模拟预编码矢量。
2.根据权利要求1所述的一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,其特征在于,在所述步骤S1中,假设卫星对用户存在角度定位误差,引入定位角度误差变量后,多波束卫星到第k个用户的下行信道参数矢量hk为:
在公式中,t表示时间,f表示系统频率,vk表示平面阵天线相应矢量,表示由于多波束卫星移动导致的多普勒频偏,τk,l表示多波束卫星到第k个用户的第l条传播径的传播延时,表示传播延时τk,l的最小值,即 表示第k个用户的下行信道增益,且满足:
3.根据权利要求2所述的一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,其特征在于,在所述步骤S1中,构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型,其还包括构建如下角度关系:
在公式中,θk为用户k相对于低轨道卫星平面阵天线y轴的准确角度,为用户k相对于低轨道卫星平面阵天线x轴的准确角度,和表示卫星对第k个用户的位置角度估计误差,和表示由于卫星姿态轨道控制而导致的公共角度检测误差,表示为用户k相对于低轨道卫星平面阵天线的y轴的估计角度,表示为用户k相对于低轨道卫星平面阵天线x轴的估计角度;
4.根据权利要求3所述的一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,其特征在于,所述步骤S2具体包括:
公式中,SINRk表示第k个用户的信干噪比,并且定义为:
公式中,平面阵天线相应矢量vk的第n项表示为:
其中,信道自相关矩阵Rk中的第m行第n列元素[Rk]m,n表示为:
其中,a和b分别为第一中间变量和第二中间变量;
则得到:
其中各中间变量分别定义为:
公式中,A、B和Z分别为第三中间变量、第四中间变量和第五中间变量;
其中,各个中间变量定义为:
其中,P、Q、D、F、C、E分别为第六中间变量、第七中间变量、第八中间变量、第九中间变量、第十中间变量和第十一中间变量。
7.根据权利要求6所述的一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,其特征在于,所述步骤S5具体包括:
固定模拟预编码矩阵WRF,则第二优化问题P2相对于数字预编码矢量wD,k的拉格朗日函数表示为:
根据KKT条件,最优数字预编码矢量需要满足:
因此得到如下关于数字预编码矢量的广义特征值问题:
再令
γk看作矩阵对(SD,k,ND,k)的最大广义特征值,uD,k是对应于γk的广义特征向量,即,
γk=max.generalized eigenvector(SD,k,ND,k),
uD,k=max.generalized eigenvector(SD,k,ND,k);
因此一旦确定了γk和uD,k,则能够计算第k个用户的数字预编码的功率值ρD,k;
定义:
数字预编码的功率矢量ρD=[ρD,1,…,ρD,K]T计算为:
因此,最优数字预编码矢量wD,k能够通过所提出的解结构精确计算。
8.根据权利要求7所述的一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,其特征在于,所述步骤S6具体包括:
固定步骤S5得到的最优数字预编码矢量wD,k,为了方便对模拟预编码矢量处理,首先需要将用户k的混合预编码矢量WRFwD,k转换为以下等价形式:
因此,模拟预编码矩阵WRF的设计简化为对矢量wA的设计;
在获得最佳模拟预编码矢量wA后,利用一个量化操作来保证该离散相位约束;
在给定数字预编码矢量wD,k以及拉格朗日乘数vk和qn的前提下,等价转化后的关于模拟预编码矢量的设计问题表述为:
上述的第三优化问题P3的拉格朗日函数表示为
根据拉格朗日函数的KKT条件,最优模拟预编码矢量wA满足如下广义特征值问题:
在部分连接的结构中,每个RF链路都连接到一个天线子集,其中对应于第j个RF链路的模拟预编码矢量wA,j被认为是独立于其他射频链路的,因此对第i个RF链路进行广义特征值分解,而不是直接计算wA;
令
对应于第i个RF链路的最优模拟预编码矢量由下式给出:
10.根据权利要求9所述的一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法,其特征在于,在所述步骤S8中,为了确定拉格朗日乘子和的最优值,采用次梯度投影技术并结合KKT条件qn≥0,和qnPn=0,得到关于的迭代表达式Q为:
根据等价优化问题P2的拉格朗日函数和其对偶函数,以及两者之间的强对偶关系,得到:
在该公式中,Ptotal是系统允许的最大发射功率;
在实际应用中,目标SINR的约束门限值γk需要尽可能地取最大值,因此求解如下凸优化问题确定vk的值:
在所述步骤S9中,在发送天线单天线功率约束下,通过采用一个简单的终止策略来停止单天线功率约束的迭代过程,即:
在该公式中,εn表示最大可超过第n个天线阵元额定功率的功率值;
当系统的和速率性能收敛时,即得到最优模拟预编码器和最优数字预编码器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210538796.0A CN114978262B (zh) | 2022-05-17 | 2022-05-17 | 一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210538796.0A CN114978262B (zh) | 2022-05-17 | 2022-05-17 | 一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114978262A true CN114978262A (zh) | 2022-08-30 |
CN114978262B CN114978262B (zh) | 2024-01-26 |
Family
ID=82982428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210538796.0A Active CN114978262B (zh) | 2022-05-17 | 2022-05-17 | 一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114978262B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115459834A (zh) * | 2022-09-19 | 2022-12-09 | 浙江大学 | 一种低轨卫星物联网的低分辨率均衡方法 |
CN115792981A (zh) * | 2023-02-06 | 2023-03-14 | 深圳大学 | 一种基于阵列天线的可见卫星探测方法 |
CN116488686A (zh) * | 2023-04-14 | 2023-07-25 | 中山大学 | 一种低轨卫星多组多播混合预编码方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101080846A (zh) * | 2004-06-03 | 2007-11-28 | 美商内数位科技公司 | 具智能天线的卫星通信用户装置及其方法 |
CN105191313A (zh) * | 2013-01-04 | 2015-12-23 | 三星电子株式会社 | 使用考虑相移的上采样的可伸缩视频编码方法和设备以及可伸缩视频解码方法和设备 |
CN105978664A (zh) * | 2016-06-24 | 2016-09-28 | 中国科学院国家空间科学中心 | 一种用于遥感卫星的高效超高速有效载荷数据传输系统 |
-
2022
- 2022-05-17 CN CN202210538796.0A patent/CN114978262B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101080846A (zh) * | 2004-06-03 | 2007-11-28 | 美商内数位科技公司 | 具智能天线的卫星通信用户装置及其方法 |
CN105191313A (zh) * | 2013-01-04 | 2015-12-23 | 三星电子株式会社 | 使用考虑相移的上采样的可伸缩视频编码方法和设备以及可伸缩视频解码方法和设备 |
CN105978664A (zh) * | 2016-06-24 | 2016-09-28 | 中国科学院国家空间科学中心 | 一种用于遥感卫星的高效超高速有效载荷数据传输系统 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115459834A (zh) * | 2022-09-19 | 2022-12-09 | 浙江大学 | 一种低轨卫星物联网的低分辨率均衡方法 |
CN115792981A (zh) * | 2023-02-06 | 2023-03-14 | 深圳大学 | 一种基于阵列天线的可见卫星探测方法 |
CN116488686A (zh) * | 2023-04-14 | 2023-07-25 | 中山大学 | 一种低轨卫星多组多播混合预编码方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN114978262B (zh) | 2024-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Méndez-Rial et al. | Hybrid MIMO architectures for millimeter wave communications: Phase shifters or switches? | |
CN111294095B (zh) | 基于统计csi的irs辅助大规模mimo无线传输方法 | |
CN114978262B (zh) | 一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法 | |
Satyanarayana et al. | Hybrid beamforming design for full-duplex millimeter wave communication | |
Eltayeb et al. | Opportunistic beam training with hybrid analog/digital codebooks for mmWave systems | |
CN113765553B (zh) | 一种基于机器学习的多波束卫星通信系统鲁棒预编码方法 | |
CN110011712B (zh) | 一种面向毫米波大规模多输入多输出的混合预编码方法 | |
CN110535500B (zh) | 一种基于深度学习的毫米波mimo混合波束成形优化方法 | |
You et al. | Hybrid analog/digital precoding for downlink massive MIMO LEO satellite communications | |
CN110943768B (zh) | 一种毫米波大规模mimo系统的混合预编码码本联合设计方法 | |
CN110311717B (zh) | 基于方向调制的稳健的混合波束成形设计方法 | |
Sun et al. | A robust secure hybrid analog and digital receive beamforming scheme for efficient interference reduction | |
Alluhaibi et al. | Capacity maximisation for hybrid digital-to-analog beamforming mm-wave systems | |
CN115085774A (zh) | 一种基于克拉美罗界的通感融合混合波束赋形方法 | |
Li et al. | Bidirectional positioning assisted hybrid beamforming for massive MIMO systems | |
CN113572503A (zh) | 一种低复杂度的基于gp的改进混合波束成形方法 | |
Yu et al. | Low-cost intelligent reflecting surface aided Terahertz multiuser massive MIMO: design and analysis | |
Song et al. | Efficient hybrid beamforming for relay assisted millimeter-wave multi-user massive MIMO | |
CN116545482A (zh) | 一种ris辅助采用低精度dac的多用户mimo下行传输方法 | |
CN114844537B (zh) | 深度学习辅助鲁棒大规模mimo收发联合方法 | |
Majidzadeh et al. | Rate maximization for partially connected hybrid beamforming in single-user MIMO systems | |
CN112636800B (zh) | 基于毫米波大规模mimo多用户场景的混合预编码方法 | |
Arvola et al. | Two-layer precoding for dimensionality reduction in massive MIMO | |
CN115001545A (zh) | 一种毫米波大规模mimo系统中的混合预编码方法 | |
Pu et al. | Low complexity hybrid beamforming based on orthogonal constraint and phase extraction |
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 |