一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法

1.本发明涉及多波束卫星通信鲁棒混合预编码方法的技术领域，特别是涉及一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法。


背景技术：

2.近年来，人们对更高质量的通信服务和全球无缝接入的需求迅速增加，特别是在后5g及未来的6g网络，这使得从太空提供互联网服务的概念引起了学术界和工业界的深刻关注。值得注意的是，卫星通信系统在提供无缝覆盖和扩展地面通信网络方面具有优势，并可实现更低的延迟、更低的生产成本、更大的覆盖范围，提供无处不在的全球无线接入。但是由于目前低轨卫星(leo)趋于小型化、低生产成本和有限的收发器功率，硬件成本和功耗成为卫星预编码器设计的关键考虑因素。
3.近几年，具有高能效和低实现成本的混合预编码技术得到了研究界的广泛关注。在混合预编码结构中，混合收发器通常将大量低成本的模拟移相器(ps)和功率放大器与低维数字信号处理单元相结合，使用少数昂贵的射频链路(rf)，这种混合预编码器的设计降低了系统生产成本和功率消耗，符合当前高能效leo卫星通信系统的设计需求。此外，在毫米波频段，由于发送天线尺寸小，因此可以将许多天线阵元封装在一个紧凑的区域内，有利于发送波形在空间域上有着良好的方向自由度和三维覆盖范围，并在具有严重路径损耗的毫米波信道提供足够的预编码增益。然而，毫米波卫星下行链路的混合预编码器的设计并不是一件容易的事，因为模拟预编码器和数字预编码器的优化问题通常是非凸的，求解最优的预编码矢量比较困难。


技术实现要素：

4.有鉴于此，本发明的目的在于提供一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法，以提高多波束卫星通信系统的传输效率，并降低硬件成本和系统功耗。
5.为了实现上述目的，本发明采用如下技术方案：
6.一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法，所述方法包括：
7.步骤s1、基于多波束卫星对各用户的位置角度估计误差以及由于卫星姿态轨道控制而导致的公共角度误差，构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型；再基于该多波束卫星下行信道矢量模型构建用户侧的接收信号模型；
8.步骤s2、对多波束卫星系统的和速率关于定位角度误差变量求数学期望，得到第k个用户的遍历传输速率，获取关于信道自相关矩阵的统计信道模型；
9.步骤s3、基于部分连接的混合预编码结构，关于模拟预编码矢量和数字预编码矢量构建多波束卫星系统和速率最大化的鲁棒混合预编码优化设计问题，约束条件为多波束卫星的每个天线阵元的功率值均小于某一门限值，其中，该优化设计问题定义为第一优化
问题；
10.步骤s4、将步骤s3中构建的第一优化问题等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题，其中，该功率最小化问题，定义为第二优化问题；
11.步骤s5、固定模拟预编码矩阵，结合第二优化问题的拉格朗日函数及其kkt条件，将混合预编码结构中数字预编码的最优结构建模为一个广义特征值问题的解，计算最优数字预编码矢量；
12.步骤s6、固定最优数字预编码矢量，结合关于模拟预编码矢量的等价转化问题的拉格朗日函数及其对应的kkt条件，将模拟预编码的最优方向建模为一个广义特征值问题的解，计算获得最优模拟预编码矢量；其中，该等价转化后的关于模拟预编码矢量的设计问题定义为第三优化问题；
13.步骤s7、对模拟预编码矢量的相位进行离散化操作，同时保证多波束卫星系统和速率最大化；
14.步骤s8、利用次梯度投影技术及拉格朗日函数与其对偶函数的强对偶性确定拉格朗日乘子；
15.步骤s9、循环迭代步骤s4、步骤s5、步骤s6、步骤s7和步骤s8直至系统和速率性能收敛，获得最优数字预编码矢量和具有离散相位的模拟预编码矢量。
16.进一步的，在所述步骤s1中，假设卫星对用户存在角度定位误差，引入定位角度误差变量后，多波束卫星到第k个用户的下行信道参数矢量hk为:
[0017][0018]
在公式中，t表示时间，f表示系统频率，vk表示平面阵天线相应矢量，表示由于多波束卫星移动导致的多普勒频偏，τ
k,l
表示多波束卫星到第k个用户的第l条传播径的传播延时，表示传播延时τ
k,l
的最小值，即表示第k个用户的下行信道增益，且满足：
[0019][0020]
其中，lk表示多波束卫星到第k个用户的总传播径数，表示到第k个用户的第l条传播径的复数增益，表示在第l条传播径上由于第k个用户移动而导致的多普勒频移，表示用户侧传输延时，且满足：
[0021]
进一步的，在所述步骤s1中，构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型，其还包括构建如下角度关系：
[0022]
[0023]
在公式中，θk为用户k相对于低轨道卫星平面阵天线y轴的准确角度，为用户k相对于低轨道卫星平面阵天线x轴的准确角度，和表示卫星对第k个用户的位置角度估计误差，和表示由于卫星姿态轨道控制而导致的公共角度检测误差，表示为用户k相对于低轨道卫星平面阵天线的y轴的估计角度，表示为用户k相对于低轨道卫星平面阵天线x轴的估计角度；
[0024]
平面阵天线相应矢量满足以下关系：
[0025][0026]
在公式中，表示平面阵天线在x轴上的阵元响应矢量，表示平面阵天线在y轴上的阵元响应矢量，且二者分别满足：
[0027][0028][0029]
进一步的，所述步骤s2具体包括：
[0030]
第k个用户的遍历传输速率定义为：
[0031][0032]
公式中，sinrk表示第k个用户的信干噪比，并且定义为：
[0033][0034]
其中，表示第k个用户的数字预编码矢量，上标(
·
)h表示共轭转置，hk为多波束卫星到第k个用户关于用户角度信息的信道相应矢量；
[0035]
由于系统遍历和速率没有闭式表达式，使得直接处理rk比较困难，因此这里根据杰森不等式确定的上限为：
[0036][0037]
再令基于多波束卫星下行信道参数模型，信道自相关矩阵rk化简为：
[0038]
[0039]
公式中，平面阵天线相应矢量vk的第n项表示为：
[0040][0041]
其中，信道自相关矩阵rk中的第m行第n列元素[rk]
m,n
表示为：
[0042][0043]
其中，表示信道功率值，[
·
]
m,n
表示取矩阵的第m行第n列的元素，(
·
)n表示取矢量的第n个元素，表示对矢量vk中的第m个元素(vk)m和第n个元素的共轭的乘积求数学期望；且：
[0044][0045]
其中，a和b分别为第一中间变量和第二中间变量；
[0046]
若卫星对第k个用户的位置角度估计误差和服从均匀分布，即：
[0047][0048]
公式中，θ
l
和θu分别表示角度误差取值的上限和下限，和分别表示角度误差取值的上限和下限；
[0049]
上式表示在区间u[θ
l
,θu]服从均匀分布，在区间服从均匀分布；
[0050]
进一步定义：δ
θ
＝θ
u-θ
l
，公式中，δ
θ
和δ
φ
分别表示角度误差和的上下限的差值；
[0051]
则得到：
[0052][0053]
将上述表达式中的和进行积分，得到：
[0054][0055]
其中各中间变量分别定义为：
[0056][0057][0058][0059]
公式中，a、b和z分别为第三中间变量、第四中间变量和第五中间变量；
[0060]
若服从均值为μ
θ,k
、方差为的高斯分布，服从均值为方差为的高斯分布，则和的概率密度函数表示为：
[0061][0062][0063]
公式中，表示的概率密度函数，表示的概率密度函数；
[0064]
此时表示为：
[0065][0066]
对上述的计算式进行不定积分得到：
[0067][0068]
其中，各个中间变量定义为：
[0069][0070][0071][0072][0073][0074][0075]
其中，p、q、d、f、c、e分别为第六中间变量、第七中间变量、第八中间变量、第九中间变量、第十中间变量和第十一中间变量。
[0076]
进一步的，所述步骤s3包括：
[0077]
优化目标设置为最大化和速率上限并受以下约束：
[0078]
1、每个天线上的天线功率限制，即其
中为天线索引集；
[0079]
2、每个ps上的离散相位，即其中rf链路索引集合为
[0080]
p1:
[0081]
因此，第一优化问题表述如下：
[0082]
在公式中，pn是第n个天线的最大可用功率。
[0083]
进一步的，所述步骤s4具体包括：
[0084]
令表示第一优化问题p1的取得最优解时各用户能够实现最大传输速率，则将该第一优化问题p1等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题p2，其具体表示为：
[0085]
p2:
[0086]
通过下面等价以使qos约束更容易处理：
[0087][0088]
其中同时近似为：
[0089][0090]
假设和均独立服从高斯分布，即和
[0091]
因此优化问题p2重新表述为第二优化问题：
[0092]
p2:
[0093]
其中进一步写为
[0094][0095]
进一步的，所述步骤s5具体包括：
[0096]
固定模拟预编码矩阵w
rf
，则第二优化问题p2相对于数字预编码矢量w
d,k
的拉格朗日函数表示为：
[0097][0098]
在公式中，是对应于qos约束的拉格朗日乘数，对应于n
t
个天线功率约束，
[0099]
再将写成以下紧凑形式：
[0100][0101]
在该公式中，对角矩阵
[0102]
则的导数表示为：
[0103][0104]
根据kkt条件，最优数字预编码矢量需要满足：
[0105][0106]
因此得到如下关于数字预编码矢量的广义特征值问题：
[0107][0108]
再令
[0109][0110]
将最优数字预编码矢量表示为其中用户k的数字预编码矢量的功率ρ
d,k
和单位范数数字预编码矢量的方向u
d,k
通过对矩阵s
d,k
和n
d,k
进行广义特征值分解来计算；
[0111]
γk看作矩阵对(s
d,k
,n
d,k
)的最大广义特征值，u
d,k
是对应于γk的广义特征向量，即，
[0112]
γk＝max.generalized eigenvector(s
d,k
,n
d,k
),
[0113]ud,k
＝max.generalized eigenvector(s
d,k
,n
d,k
)；
[0114]
注意到问题p2的约束在最优解处等式成立，即
[0115][0116]
因此一旦确定了γk和u
d,k
，则能够计算第k个用户的数字预编码的功率值ρ
d,k
；
[0117]
定义：
[0118][0119]
该公式中，矩阵
[0120]
那么等价写为：
[0121][0122]
数字预编码的功率矢量ρd＝[ρ
d,1
,
…
,ρ
d,k
]
t
计算为：
[0123][0124]
因此，最优数字预编码矢量w
d,k
能够通过所提出的解结构精确计算。
[0125]
进一步的，所述步骤s6具体包括：
[0126]
固定步骤s5得到的最优数字预编码矢量w
d,k
，为了方便对模拟预编码矢量处理，首先需要将用户k的混合预编码矢量w
rfwd,k
转换为以下等价形式：
[0127][0128]
在该公式中，对角矩阵定义为：
[0129]
[0130]
矢量是由对应于n
rf
个射频链路的模拟预编码矢量堆叠而成；
[0131]
因此，模拟预编码矩阵w
rf
的设计简化为对矢量wa的设计；
[0132]
省略第二优化问题p2中对于相移器的离散相位约束
[0133]
在获得最佳模拟预编码矢量wa后，利用一个量化操作来保证该离散相位约束；
[0134]
在给定数字预编码矢量w
d,k
以及拉格朗日乘数vk和qn的前提下，等价转化后的关于模拟预编码矢量的设计问题表述为：
[0135]
p3:
[0136]
不等式等价表示为:
[0137][0138]
上述的第三优化问题p3的拉格朗日函数表示为
[0139][0140]
在该公式中，定义为：写成：
[0141][0142]
根据迹求偏导的计算结构：拉格朗日函数关于模拟预编码矢量wa的偏导表示为：
[0143][0144]
根据拉格朗日函数的kkt条件，最优模拟预编码矢量wa满足如下广义特征值问题：
[0145][0146]
在部分连接的结构中，每个rf链路都连接到一个天线子集，其中对应于第j个rf链路的模拟预编码矢量w
a,j
被认为是独立于其他射频链路的，因此对第i个rf链路进行广义特征值分解，而不是直接计算wa；
[0147]
令
[0148][0149]
定义矩阵和分别取自矩阵sa和na的第((i-1)n
t
1)到第(in
t
1)列，第((i-1)n
t
1)到第(in
t
1)行；矢量w
a,i
满足：
[0150][0151]
对应于第i个rf链路的最优模拟预编码矢量由下式给出：
[0152][0153]
进一步的，所述步骤s7具体包括：
[0154]
对模拟预编码矢量的相位进行离散化操作，即将的连续相位映射到可行的离散相位集中，相位量化操作为：
[0155][0156]
在该公式中，表示量化到离散相位集中最近的离散点，模拟预编码矢量的最佳离散相位计算为：
[0157]
进一步的，在所述步骤s8中，为了确定拉格朗日乘子和的最优值，采用次梯度投影技术并结合kkt条件qn≥0,和qnpn＝0,得到关于的迭代表达式q为：
[0158][0159]
在该公式中，混合预编码矢量ti为迭代步长且ti＝1/i，取最大值的运算符是按元素定义的；
[0160]
根据等价优化问题p2的拉格朗日函数和其对偶函数，以及两者之间的强对偶关
系，得到：
[0161][0162]
在该公式中，p
total
是系统允许的最大发射功率；
[0163]
将上述计算式进行重新排列，得到关于的不动点方程式：
[0164][0165]
在该公式中，最优拉格朗日乘子矩阵
[0166]
在实际应用中，目标sinr的约束门限值γk需要尽可能地取最大值，因此求解如下凸优化问题确定vk的值：
[0167][0168]
在所述步骤s9中，在发送天线单天线功率约束下，通过采用一个简单的终止策略来停止单天线功率约束的迭代过程，即：
[0169][0170]
在该公式中，εn表示最大可超过第n个天线阵元额定功率的功率值；
[0171]
当系统的和速率性能收敛时，即得到最优模拟预编码器和最优数字预编码器。
[0172]
本发明的有益效果是：
[0173]
1、本发明通过对系统和速率求数学期望并结合对单天线阵元的功率约束，有效减少多波束卫星用户定位角度不准确性带来的不利影响。此外，利用具有低分辨率相移器的混合预编码结构，能够在保证多波束卫星传输速率的同时，提高卫星的发射功率效率，降低硬件实现成本；
[0174]
2、本发明将卫星对用户定位时存在的角度不确定性以及由于卫星姿态抖动而导致的检测角度误差融合至信道建模中，并且根据角度误差变量服从某一概率分布这一特点，建立了能够实现混合预编码性能对定位角度误差鲁棒性的统计信道模型；
[0175]
3、本发明为多波束低轨卫星下行通信场景建立了一种搭载低分辨率移相器的部分连接混合预编码结构。为了保证每个天线的功率放大器工作在线性区域内，避免信号失真，本发明采用了一种更加实际的单天线功率约束。
[0176]
4、本发明在单天线阵元的功率约束下，实现卫星通信系统的遍历和速率最大化，分别建立关于数字预编码矢量和模拟预编码矢量的优化问题，并根据各自拉格朗日函数及其kkt条件，将最优鲁棒数字预编码和模拟预编码的结构建模为对应广义特征值问题的解。同时，本发明开发出了一种高效的迭代算法求解最优数字预编码矢量和模拟预编码矢量，进而获得最优的混合预编码器分连接结构下的混合预编码器的示意图。
附图说明
[0177]
图1为本实施例1中提出的一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法的流程示意图；
[0178]
图2为本实施例1中卫星天线面阵对用户存在定位角度误差的示意图；
[0179]
图3为本实施例1采用的部分连接结构下的混合预编码器的示意图。
具体实施方式
[0180]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0181]
实施例1
[0182]
参见图1-图3，本实施例提供一种基于低分辨率相移器的多波束卫星通信鲁棒混合预编码方法，该方法的具体流程如图1所示，该方法具体包括如下的步骤：
[0183]
步骤s1、基于多波束卫星对各用户的位置角度估计误差以及由于卫星姿态轨道控制而导致的公共角度误差，构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型；再基于该多波束卫星下行信道矢量模型构建用户侧的接收信号模型；
[0184]
具体的说，在本实施例中，该步骤s1包括：
[0185]
针对一多波束卫星下行传输场景，在该场景中，leo卫星同时为k单天线用户提供服务；且在卫星侧配备均匀平面阵列(upa)，由n
t
＝m
xmy
天线组成，其中m
x
和my是在分别为x和y轴上的天线数。射频(rf)链路数为n
rf
大于或等于用户数，即n
rf
≥k，每条射频链路与带有天线的子天线阵列连接。
[0186]
如图2所示，假设卫星对用户存在角度定位误差，引入定位角度误差变量后，多波束卫星到第k个用户的下行信道参数矢量hk为:
[0187][0188]
在公式中，t表示时间，f表示系统频率，vk表示平面阵天线相应矢量，表示由于多波束卫星移动导致的多普勒频偏，τ
k,l
表示多波束卫星到第k个用户的第l条传播径的传播延时，表示传播延时τ
k,l
的最小值，即表示第k个用户的下行信道增益，且满足：
[0189][0190]
其中，lk表示多波束卫星到第k个用户的总传播径数，表示到第k个用户的第l条传播径的复数增益，表示在第l条传播径上由于第k个用户移动而导致的多普勒频移，表示用户侧传输延时，且满足：
[0191]
在多波束卫星通信场景中，信号沿直射径传播，因此卫星的第k个用户的下行信道增益可以被视为服从莱斯因子为功率值为的莱斯分布。
[0192]
具体的说，在本实施例中，构建多波束卫星下行信道矢量模型还包括：
[0193]
构建如下角度关系：
[0194][0195]
在公式中，θk为用户k相对于低轨道卫星平面阵天线y轴的准确角度，为用户k相对于低轨道卫星平面阵天线x轴的准确角度，和表示卫星对第k个用户的位置角度估计误差，和表示由于卫星姿态轨道控制而导致的公共角度检测误差，表示为用户k相对于低轨道卫星平面阵天线的y轴的估计角度，表示为用户k相对于低轨道卫星平面阵天线x轴的估计角度。
[0196]
进一步将平面阵天线相应矢量满足以下关系：
[0197][0198]
在公式中，表示平面阵天线在x轴上的阵元响应矢量，表示平面阵天线在y轴上的阵元响应矢量，且二者分别满足：
[0199][0200]
如图3所示，根据混合预编码技术，k个用户的传输符号首先由数字基带预编码器进行处理，之后通过n
rf
个rf链路上变频到射频域，然后信号由维度为n
t
×nrf
的模拟预编码器w
rf
处理得到适合传输的发送波形。
[0201]
由于模拟预编码器配备了移相器和功率放大器，因此模拟预编码器可以同时对输入的信号进行幅度和相位的变化。本实施例定义了模拟预编码器的第(i,j)个元素其中，连续相位θ
i,j
在[0,2π]区间内取值，a
i,j
为信号的幅值。在
实际卫星系统中，θ
i,j
通常通过离散相位来近似以降低硬件成本，因此可以通过一组低分辨率相移器实现，并且离散相位的取值属于集合其中b是控制量化相位的精度的比特数。
[0202]
为了能够表示模拟预编码器中的低分辨率ps，这里定义了其中集合是每个ps的可能取值集，可以看出，比较大的比特数b可以实现更高分辨率的相移器，进而实现更高的系统性能，但代价是更高的硬件成本和功耗。
[0203]
具体的说，在本实施例中，定义sk为从leo卫星发送到第k个用户的传输信号，那么第k个用户的接收信号可表示为
[0204][0205]
在公式中，为leo卫星和第k个用户之间的信道参数矢量。是用户索引集；表示零均值和方差为的加性高斯白噪声(awgn)，信号sk的平均功率被假定为1，即表示第k个用户的数字预编码权值矢量；
[0206]
在部分连接架构中，模拟预编码权值矩阵w
rf
可以表示为
[0207][0208]
其中是对应第i个射频链路的模拟预编码权值矢量，i＝1,
…
,n
rf
。
[0209]
步骤s2、对多波束卫星系统的和速率关于定位角度误差变量求数学期望，得到第k个用户的遍历传输速率，获取关于信道自相关矩阵的统计信道模型；
[0210]
具体的说，在本实施例中，该步骤s2包括：
[0211]
第k个用户的遍历传输速率定义为：
[0212][0213]
公式中，sinrk表示第k个用户的信干噪比，并且定义为：
[0214][0215]
其中，表示第k个用户的数字预编码矢量，上标(
·
)h表示共轭转置，hk为多波束卫星到第k个用户关于用户角度信息的信道相应矢量。
[0216]
由于系统遍历和速率没有闭式表达式，使得直接处理比较困难，因此这里根据杰森不等式确定的上限为：
[0217][0218]
再令基于多波束卫星下行信道参数模型，信道自相关矩阵rk化简为：
[0219][0220]
公式中，平面阵天线相应矢量vk的第n项可以表示为：
[0221][0222]
其中，信道自相关矩阵rk中的第m行第n列元素[rk]
m,n
可以表示为：
[0223][0224]
其中，表示信道功率值，[
·
]
m,n
表示取矩阵的第m行第n列的元素，(
·
)n表示取矢量的第n个元素，表示对矢量vk中的第m个元素(vk)m和第n个元素的共轭的乘积求数学期望；且：
[0225][0226]
其中，a和b分别为第一中间变量和第二中间变量。
[0227]
若卫星对第k个用户的位置角度估计误差和服从均匀分布，即：
[0228][0229]
公式中，θ
l
和θu分别表示角度误差取值的上限和下限，和分别表示角度误差取值的上限和下限；
[0230]
上式表示在区间u[θ
l
,θu]服从均匀分布，在区间服从均匀分布；
[0231]
进一步定义：δ
θ
＝θ
u-θ
l
，
[0232]
公式中，δ
θ
和δ
φ
分别表示角度误差和的上下限的差值；
[0233]
则可以得到：
[0234][0235]
将上述表达式中的和进行积分，得到：
[0236]
[0237]
其中各中间变量分别定义为：
[0238][0239][0240][0241]
公式中，a、b和z分别为第三中间变量、第四中间变量和第五中间变量；
[0242]
若服从均值为μ
θ,k
、方差为的高斯分布，服从均值为方差为的高斯分布，则和的概率密度函数表示为：
[0243][0244][0245]
公式中，表示的概率密度函数，表示的概率密度函数；
[0246]
此时可以表示为：
[0247][0248]
对上述的计算式进行不定积分得到：
[0249][0250]
其中，各个中间变量定义为：
[0251][0252][0253][0254][0255][0256][0257]
其中，p、q、d、f、c、e分别为第六中间变量、第七中间变量、第八中间变量、第九中间变量、第十中间变量和第十一中间变量。
[0258]
步骤s3、基于部分连接的混合预编码结构，关于模拟预编码矢量和数字预编码矢量构建多波束卫星系统和速率最大化的鲁棒混合预编码优化设计问题，约束条件为多波束卫星的每个天线阵元的功率值均小于某一门限值，其中，该优化设计问题定义为第一优化问题；
[0259]
具体的说，在本实施例中，为了使相应的优化问题易于处理并与实际实现场景保持一致，优化目标设置为最大化和速率上限并受以下约束：
[0260]
1、每个天线上的天线功率限制，即其中为天线索引集；
[0261]
2、每个ps上的离散相位，即其中rf链路索引集合为
[0262]
因此，第一优化问题可以表述如下：
[0263]
p1:
[0264]
在公式中，pn是第n个天线的最大可用功率。
[0265]
步骤s4、将步骤s3中构建的第一优化问题等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题，其中，该功率最小化问题，定义为第二优化问题；
[0266]
具体的说，在本实施例中，令表示第一优化问题p1的取得最优解时各用户能够实现最大传输速率，则将该第一优化问题p1等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题p2，其具体表示为：
[0267]
p2:
[0268]
进一步，通过下面等价以使qos约束更容易处理：
[0269][0270]
其中同时e{sinrk}可以近似为：
[0271][0272]
在本实施例中假设和均独立服从高斯分布，即和因此优化问题p2可以重新表述第二优化问题：
[0273]
p2:
[0274]
其中可以进一步写为
[0275][0276]
步骤s5、固定模拟预编码矩阵，结合第二优化问题的拉格朗日函数及其kkt条件，将混合预编码结构中数字预编码的最优结构建模为一个广义特征值问题的解，计算最优数字预编码矢量；
[0277]
具体的说，在本实施例中，固定模拟预编码矩阵w
rf
，则第二优化问题p2相对于数字预编码矢量w
d,k
的拉格朗日函数可以表示为：
[0278][0279]
在公式中，是对应于qos约束的拉格朗日乘数，对应于n
t
个天线功率约束，
[0280]
为方便起见，可以写成以下紧凑形式：
[0281][0282]
在该公式中，对角矩阵
[0283]
则的导数可以表示为：
[0284][0285]
根据kkt条件，最优数字预编码矢量需要满足：
[0286][0287]
因此可以得到如下关于数字预编码矢量的广义特征值问题：
[0288][0289]
再令
[0290][0291]
将最优数字预编码矢量表示为其中用户k的数字预编码矢量的功率ρ
d,k
和单位范数数字预编码矢量的方向u
d,k
可以通过对矩阵s
d,k
和n
d,k
进行广义特征值分解来计算。特别地，γk可以看作矩阵对(s
d,k
,n
d,k
)的最大广义特征值，u
d,k
是对应于γk的广义特征向量，即，
[0292]
γk＝max.generalized eigenvector(s
d,k
,n
d,k
),
[0293]ud,k
＝max.generalized eigenvector(s
d,k
,n
d,k
)；
[0294]
注意到问题p2的约束在最优解处等式成立，即
[0295][0296]
因此一旦确定了γk和u
d,k
，则可计算第k个用户的数字预编码的功率值ρ
d,k
。
[0297]
定义：
[0298][0299]
该公式中，矩阵
[0300]
那么可以等价写为：
[0301][0302]
数字预编码的功率矢量ρd＝[ρ
d,1
,
…
,ρ
d,k
]
t
可以计算为：
[0303][0304]
因此，最优数字预编码矢量w
d,k
可以通过所提出的解结构精确计算。
[0305]
步骤s6、固定最优数字预编码矢量，结合关于模拟预编码矢量的等价转化问题的拉格朗日函数及其对应的kkt条件，将模拟预编码的最优方向建模为一个广义特征值问题的解，计算获得最优模拟预编码矢量；其中，该等价转化后的关于模拟预编码矢量的设计问题定义为第三优化问题；
[0306]
具体的说，在本实施例中，固定步骤s5得到的最优数字预编码矢量w
d,k
，为了方便对模拟预编码矢量处理，首先需要将用户k的混合预编码矢量w
rfwd,k
转换为以下等价形式：
[0307][0308]
在该公式中，对角矩阵定义为：
[0309][0310]
矢量是由对应于n
rf
个射频链路的模拟预编码矢量堆叠而成。
[0311]
因此，模拟预编码矩阵w
rf
的设计可以简化为对矢量wa的设计。
[0312]
暂时省略第二优化问题p2中对于相移器的离散相位约束在获得最佳模拟预编码矢量wa后，利用一个量化操作来保证该离散相位约束。
[0313]
在给定数字预编码矢量w
d,k
以及拉格朗日乘数vk和qn的前提下，等价转化后的关于模拟预编码矢量的设计问题表述为：
[0314]
p3:
[0315]
类似地，不等式可以等价表示为:
[0316][0317]
进一步的，问题p3的拉格朗日函数可以表示为
[0318][0319]
在该公式中，定义为：可以写成：
[0320][0321]
根据迹求偏导的计算结构：拉格朗日函数关于模拟预编码矢量wa的偏导表示为：
[0322][0323]
根据拉格朗日函数的kkt条件，最优模拟预编码矢量wa满足如下广义特征值问题：
[0324][0325]
在部分连接的结构中，每个rf链路都连接到一个天线子集，其中对应于第j个rf链路的模拟预编码矢量w
a,j
可以被认为是独立于其他射频链路的，因此对第i个rf链路进行广义特征值分解，而不是直接计算wa。令
[0326][0327]
定义矩阵和分别取自矩阵sa和na的第((i-1)n
t
1)到第(in
t
1)列，第((i-1)n
t
1)到第(in
t
1)行。矢量w
a,i
满足：
[0328][0329]
对应于第i个rf链路的最优模拟预编码矢量由下式给出：
[0330][0331]
步骤s7、对模拟预编码矢量的相位进行离散化操作，同时保证多波束卫星系统和速率最大化；
[0332]
具体的说，在本实施例中，对模拟预编码矢量的相位进行离散化操作，即将的连续相位映射到可行的离散相位集中，相位量化操作为：
[0333][0334]
在该公式中，表示量化到离散相位集中最近的离散点。模拟预编码矢量的最
佳离散相位可以计算为：
[0335]
步骤s8、利用次梯度投影技术及拉格朗日函数与其对偶函数的强对偶性确定拉格朗日乘子；具体的说，在本实施例中，为了确定拉格朗日乘子和的最优值，采用次梯度投影技术并结合kkt条件qn≥0,和qnpn＝0,得到关于的迭代表达式q为：
[0336][0337]
在该公式中，混合预编码矢量ti为迭代步长且ti＝1/i，取最大值的运算符是按元素定义的。
[0338]
根据等价优化问题p2的拉格朗日函数和其对偶函数，以及两者之间的强对偶关系，可以得到：
[0339][0340]
在该公式中，p
total
是系统允许的最大发射功率。
[0341]
将上述计算式进行重新排列，可以关于的不动点方程式：
[0342][0343]
在该公式中，最优拉格朗日乘子矩阵
[0344]
在实际应用中，目标sinr的约束门限值γk需要尽可能地取最大值，因此求解如下凸优化问题确定vk的值：
[0345][0346]
步骤s9、循环迭代步骤s4、步骤s5、步骤s6、步骤s7和步骤s8直至系统和速率性能收敛，获得最优数字预编码矢量和具有离散相位的模拟预编码矢量。
[0347]
具体的说，在本实施例中，在发送天线单天线功率约束下，可以采用一个简单的终止策略来停止单天线功率约束的迭代过程，即：
[0348][0349]
在该公式中，εn表示最大可超过第n个天线阵元额定功率的功率值。当系统的和速率性能收敛时，即得到最优模拟预编码器和最优数字预编码器。
[0350]
本发明未详述之处，均为本领域技术人员的公知技术。
[0351]
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。