一种基于机器学习的多波束卫星通信系统鲁棒预编码方法与流程

1.本发明涉及卫星通信的技术领域，尤其涉及一种基于机器学习的多波束卫星通信系统鲁棒预编码方法。


背景技术：

2.目前，多波束卫星系统已经显示出其在5g网络以及未来6g网络实现无处不在的全球无线接入的巨大潜力。其中下行链路的预编码器设计在卫星通信中起着至关重要的作用，而现有的预编码器设计方法通常是基于完全已知的信道状态信息和总功率约束下进行的。由于卫星的高速移动性以及卫星姿态抖动，理想的信道状态信息通常很难在低轨道多波束卫星上获得，现有的预编码器设计方法会导致卫星下行预编码的性能下降，进而导致对用户的服务质量下降。另外，伴随着卫星通常趋向小型化以及信息传输处理过程中能量消耗的急剧增长，使得卫星通信系统中的功率消耗在系统设计时成为需要着重考虑的因素。因此如何实现低功耗和信息鲁棒的传输是当前卫星通信系统设计的趋势。
3.近些年来，人工智能技术飞速发展，并且在通信领域等到广泛的应用。该技术通过离线训练神经网络，并直接在线预测结果的方式，能够有效降低算法实现复杂度。在实际的卫星通信系统中，信号传输距离较远，导致传输延时较高，因此在系统设计过程中，信号传输的实时性是非常重要的。如何结合机器学习技术，降低预编码方法的实现复杂度和系统处理时间，进而提升卫星通信系统的信号传输的实时性，降低星载设备的功率消耗是当前卫星通信系统设计中需要解决的问题。


技术实现要素：

4.发明目的：为了克服现有技术中存在的不足，本发明提供一种基于机器学习的多波束卫星通信系统鲁棒预编码方法，该发明能够有效减少多波束卫星用户定位角度不准确性带来的不利影响，以及提高多波束卫星通信系统的实时性和传输性能。
5.技术方案：为了实现上述发明目的，本发明提供了一种基于机器学习的多波束卫星通信系统鲁棒预编码方法，包括以下步骤，
6.步骤1，基于多波束卫星对各用户的位置角度估计误差以及由于卫星姿态轨道控制而导致的公共角度误差，构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型；
7.步骤2，对多波束卫星系统的和速率关于定位角度误差变量求数学期望，得到第k个用户的遍历传输速率，获取关于信道自相关矩阵的统计信道模型；
8.步骤3，构建多波束卫星系统和速率最大化的鲁棒预编码优化设计问题，约束条件为多波束卫星的每个天线阵元的功率值均小于某一门限值；
9.步骤4，将多波束卫星系统和速率最大化的鲁棒预编码优化设计问题等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题；
10.步骤5，结合等价优化问题的拉格朗日函数及其kkt条件，将鲁棒预编码的最优结
构建模为其广义特征值问题的解，获得最优预编码矢量，即原始问题的最优解；
11.步骤6，结合机器学习的方法基于信道自相关矩阵预测优化问题所需的拉格朗日乘子。
12.进一步的，在本发明中：所述步骤1中构建多波束卫星下行信道矢量模型还包括，构建如下角度关系：
[0013][0014][0015]
其中，θ
k
为用户k相对于低轨道卫星平面阵天线y轴的准确角度，为用户 k相对于低轨道卫星平面阵天线
x
轴的准确角度，和表示卫星对第k个用户的位置角度估计误差，和表示由于卫星姿态轨道控制而导致的公共角度检测误差，表示为用户k相对于低轨道卫星平面阵天线的y轴的估计角度，表示为用户k相对于低轨道卫星平面阵天线
x
轴的估计角度。
[0016]
进一步的，在本发明中：所述步骤1中，平面阵天线相应矢量满足以下关系：
[0017][0018]
其中，表示平面阵天线在
x
轴上的阵元响应矢量，表示平面阵天线在y轴上的阵元响应矢量，且二者分别满足：
[0019][0020][0021]
进一步的，在本发明中：所述步骤2中第k个用户的多波束卫星通信系统的遍历和速率为：
[0022][0023]
其中，且k为总用户数，表示计算数学期望符号，表示噪声功率值，w
k
为第k个用户的预编码矢量，w
i
为第i个用户的预编码矢量， h
k
为多波束卫星到第k个用户的下行信道参数矢量，h
i
为多波束卫星到第i个用户的下行信道参数矢量，上标(
·
)
h
表示对矢量或者矩阵进行共轭转置操作，表示对矢量或者矩阵进行共轭转置操作，表示矢量w
i
和h
i
的共轭转置；
[0024]
sinr
k
为第k个用户的信干噪比，sinr
k
满足：
[0025][0026]
其中，w
k
为第k个用户的预编码矢量，上标(
·
)
h
表示共轭转置，h
k
为多波束卫星到
第k个用户关于用户角度信息的信道相应矢量。
[0027]
进一步的，在本发明中：所述步骤3中由于系统遍历和速率没有闭式表达式，使得直接处理比较困难，因此这里根据杰森不等式确定的上限为
[0028][0029]
所述系统和速率上限最大化的鲁棒预编码优化设计问题p1建模表示为：
[0030][0031][0032]
其中p
n
为第
n
个天线阵元的最大允许的发射功率，集合满足集合满足
[0033]
进一步的，在本发明中：所述步骤4中，单天线功率约束下系统和速率最大的优化问题p1等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题p2，
[0034][0035][0036][0037]
用户信干噪比保障约束e{sinr
k
}的表达式可以近似为：
[0038][0039]
其中，目标sinr的约束门限值为其中，目标sinr的约束门限值为为问题p1的能够实现最大和速率，表达式的推导需要对角度误差变量求数学期望。
[0040]
进一步的，在本发明中：所述步骤4中，基于多波束卫星下行信道参数模型，将信道自相关矩阵化简为：
[0041][0042]
其中，平面阵天线相应矢量的第
n
项可以表示为：
[0043]
[0044]
n＝0,
…
,m
x
m
y
‑
1,
[0045]
其中，信道自相关矩阵中的第m行第n列元素可以表示为：
[0046][0047]
其中，表示信道功率值，[
·
]
m,n
表示取矩阵的第
m
行第
n
列的元素，(
·
)
n
表示取矢量的第n个元素，表示对矢量中的第m个元素和第n 个元素的共轭的乘积求数学期望；
[0048]
且：
[0049][0050]
m＝0,
…
,m
x
m
y
‑
1,
[0051]
n＝0,
…
,m
x
m
y
‑
1,
[0052]
其中，a和b分别为第一中间变量和第二中间变量；
[0053]
若卫星对第k个用户的位置角度估计误差和服从均匀分布，即：
[0054][0055][0056]
其中，θ
l
和θ
u
分别表示角度误差取值的上限和下限，和分别表示角度误差取值的上限和下限；
[0057]
上式表示在区间u[θ
l
,θ
u
]服从均匀分布，在区间服从均匀分布，进一步定义：
[0058]
δ
θ
＝θ
u
‑
θ
l
[0059][0060]
其中，δ
θ
和分别表示角度误差和的上下限的差值；
[0061]
则可以得到：
[0062][0063]
将上述表达式中的和进行积分，得到：
[0064][0065]
其中各中间变量分别定义为
[0066][0067][0068][0069]
其中，a、b和z分别为第三中间变量、第四中间变量和第五中间变量；
[0070]
若服从均值为μ
θ,k
、方差为的高斯分布，服从均值为方差为的高斯分布，则和的概率密度函数表示为：
[0071][0072][0073]
其中，表示的概率密度函数，表示的概率密度函数；
[0074]
此时可以表示为：
[0075][0076]
对上述的计算式进行不定积分得到：
[0077][0078]
其中，各个中间变量定义为
[0079][0080][0081][0082][0083][0084][0085]
其中，p、q、d、f、c、e分别为第六中间变量、第七中间变量、第八中间变量、第九中间变量、第十中间变量和第十一中间变量
[0086]
进一步的，在本发明中：所述步骤5中等价优化问题p2的拉格朗日函数为：
[0087][0088]
其中，为对应于sinr约束的拉格朗日乘子，对应于n
t
个单天线功率约束的拉格朗日乘子，e
n
表示第n个元素为1、其他元素为0的矢量，其维度由前后相乘矩阵决定，表示对矢量e
n
进行转置，p
n
表示第n个天线阵元的功率上限值；
[0089]
为了简化计算，令则：
[0090][0091]
其中，对角矩阵对拉格朗日函数求导得到：
[0092][0093]
根据kkt条件，最优的预编码矢量需要满足：
[0094][0095]
因此可以得到：
[0096][0097]
令矩阵并且假设最优预编码矢量量为最优预编码矢量，则可以将最优预编码矢量的方向视为对应矩阵对(s
k
,n
k
)最大广义特征值问题的特征向量，最大广义特征值为γ
k
，和γ
k
分别为：
[0098]
γ
k
＝max.generalized eigenvalue(s
k
,n
k
)
[0099][0100]
当拉格朗日乘子和确定后，最优预编码器的方向和γ
k
可以唯一确定，并且当等价优化问题p3取得最优解时，其sinr约束满足e{sinr
k
}＝γ
k
，据此条件可以求得最优的功率分配，即：
[0101][0102]
其中，ρ为功率分配矢量，且ρ＝[ρ1,
…
,ρ
k
]
t
，矩阵f的维度为k
×
k，其第 k行第i列的元素为：
[0103][0104]
此时可以计算得到最优预编码矢量
[0105]
进一步的，在本发明中：所述步骤5中，为了确定拉格朗日乘子和的最优值，采用次梯度投影技术并结合kkt条件q
n
≥0,和 q
n
p
n
＝0,得到关于的迭代表达式q为：
[0106][0107]
其中，q
i
表示为第i迭代过程中的对应于单天线功率约束的拉格朗日乘子矩阵t
i
表示迭代步长，且t
i
＝1/i；
[0108]
根据等价优化问题p2的拉格朗日函数和其对偶函数，以及两者之间的强对偶关系，可以得到：
[0109][0110]
其中，p
t
为系统实际发送功率；将上述计算式进行重新排列，可以关于的不动点方程式：
[0111][0112]
其中，最优拉格朗日乘子矩阵
[0113]
在实际应用中，目标sinr的约束门限值γ
k
需要尽可能地取最大值，结合下述的凸优化问题确定v
k
的值：
[0114][0115][0116][0117]
其中，p
total
表示系统最大传输功率；为了降低上述优化问题的求解复杂度，将关于v
k
的不动点方程近似为：
[0118][0119]
其中计算辅助因子β
k
＝||h
k
||2，上述优化问题通过注水算法进行求解。
[0120]
进一步的，在本发明中：所述步骤6还包括以下步骤，
[0121]
步骤6
‑
1，构建卷积神经网络；
[0122]
步骤6
‑
2，收集训练样本并对卷积神经网络进行训练；
[0123]
步骤6
‑
3，训练后的神经网络以信道自相关参数矩阵的实部和虚部作为输入，拉格朗日乘子和作为输出；
[0124]
步骤6
‑
4，输入信道自相关参数矩阵通过神经网络直接获得拉格朗日乘子。
[0125]
有益效果：本发明与现有技术相比，其有益效果是：通过对系统和速率求数学期望并结合对单天线阵元的功率约束，有效减少多波束卫星用户定位角度不准确性带来的不利影响，在人工智能技术的辅助下，与不考虑卫星定位存在角度误差的传统方法相比，提高多波束卫星通信系统的实时性和传输性能；
[0126]
(1)本发明将卫星对用户定位时存在的角度不确定性以及由于卫星姿态抖动而导致的检测角度误差融合至信道建模中，并且根据角度误差变量服从某一概率分布这一特点，建立了能够实现预编码性能对定位角度误差鲁棒性的统计信道模型；
[0127]
(2)本发明在单天线阵元的功率约束下，实现卫星通信系统的遍历和速率最大化，根据优化问题的拉格朗日函数及其kkt条件，通过将最优鲁棒预编码的结构建模为广义特征值问题的解。同时，本发明开发出了一种高效的迭代算法求解最优预编码方向及功率分配，进而获得最优的预编码矢量
[0128]
(3)本发明结合机器学习技术，通过离线训练的方式，实现直接对求解过程中拉格朗日乘子进行预测，取代求解过程中的算法迭代过程，满足卫星通信系统中对计算复杂度以及系统实时性的要求，降低了系统实现复杂度和系统功率消耗。
附图说明
[0129]
图1为本发明提出的基于机器学习的多波束卫星通信系统鲁棒预编码方法的整体流程示意图；
[0130]
图2为本发明中配备平面阵天线的多波束卫星系统结构示意图；
[0131]
图3为本发明中卫星侧天线面阵对用户存在定位角度误差的示意图；
[0132]
图4为本发明中用于预测拉格朗日乘子的卷积神经网络结构示意图。
具体实施方式
[0133]
下面结合附图对本发明的技术方案做进一步的详细说明：
[0134]
本发明可以用许多不同的形式实现，而不应当认为限于这里所述的实施例。相反，提供这些实施例以便使本公开透彻且完整，并且将向本领域技术人员充分表达本发明的范围。
[0135]
如图1所示，为本发明提出的一种基于机器学习的多波束卫星通信系统鲁棒预编码方法的整体流程示意图，该方法包括以下步骤：
[0136]
步骤1，基于多波束卫星对各用户的位置角度估计误差以及由于卫星姿态轨道控制而导致的公共角度误差，构建包含用户位置定位不确定性的多波束卫星下行信道矢量模型；
[0137]
参照图2的示意，图2为多波束卫星系统示意图，假设在多波束卫星下行传输场景中，同时服务k个单天线用户，并且多波束卫星配备的平面阵天线的天线阵元数n
t
为：
[0138]
n
t
＝m
x
×
m
y
[0139]
其中，m
x
和m
y
分别为面阵天线x轴和y轴上的天线阵元数，x轴表示水平方向，y轴表示垂直方向。
[0140]
参照图3的示意，假设卫星对用户存在角度定位误差，引入定位角度误差变量后，多波束卫星到第k个用户的下行信道参数矢量h
k
为:
[0141][0142]
其中，t表示时间，f表示系统频率，表示平面阵天线相应矢量，表示由于多波束卫星移动导致的多普勒频偏，τ
k,l
表示多波束卫星到第k个用户的第l条传播径的传播延时，表示传播延时τ
k,l
的最小值，即的最小值，即表示第k个用户的下行信道增益，且满足：
[0143][0144]
其中，l
k
表示多波束卫星到第k个用户的总传播径数，表示到第k个用户的第l条传播径的复数增益，表示在第l条传播径上由于第k个用户移动而导致的多普勒频移，表示用户侧传输延时且，且满足：
[0145][0146]
在多波束卫星通信场景中，信号沿直射径传播，因此卫星的第k个用户的下行信道增益可以被视为服从莱斯因子为k
k
、功率值为的莱斯分布。
[0147]
具体的，构建多波束卫星下行信道矢量模型还包括，构建如下角度关系：
[0148][0149][0150]
其中，θ
k
为用户k相对于低轨道卫星平面阵天线y轴的准确角度，为用户k相对于低轨道卫星平面阵天线x轴的准确角度，和表示卫星对第k个用户的位置角度估计误差，和表示由于卫星姿态轨道控制而导致的公共角度检测误差，表示为用户k相对于低轨道卫星平面阵天线的y轴的估计角度，表示为用户k相对于低轨道卫星平面阵天线x轴的估计角度。进一步将平面阵天线相应矢量满足以下关系：
[0151][0152]
其中，表示平面阵天线在x轴上的阵元响应矢量，表示平面阵天线在y轴上的阵元响应矢量，且二者分别满足：
[0153][0154][0155]
步骤2，对多波束卫星系统的和速率关于定位角度误差变量求数学期望，得到第k个用户的遍历传输速率，获取关于信道自相关矩阵的统计信道模型；
[0156]
具体的，第k个用户的多波束卫星通信系统的遍历和速率为：
[0157][0158]
其中，k为总用户数，且k为总用户数，且表示计算数学期望符号，，表示噪声功率值，w
k
为第k个用户的预编码矢量，w
i
为第i个用户的预编码矢量，h
k
为多波束卫星到第k个用户的下行信道参数矢量，h
i
为多波束卫星到第i个用户的下行信道参数矢量，上标(
·
)
h
表示对矢量或者矩阵进行共轭转置操作，表示矢量w
i
和h
i
的共轭转置；
[0159]
sinr
k
为第k个用户的信干噪比，且sinr
k
满足：
[0160][0161]
步骤3，构建多波束卫星系统和速率最大化的鲁棒预编码优化设计问题，约束条件为多波束卫星的每个天线阵元的功率值均小于某一门限值；
[0162]
其中，多波束卫星系统和速率最大化的鲁棒预编码优化设计问题p1的优化目标为
[0163]
具体的，由于系统遍历和速率没有闭式表达式，使得直接处理比较困难，因此这里根据杰森不等式确定的上限并表示为
[0164][0165]
则多波束卫星系统和速率最大化的鲁棒预编码优化设计问题p1可建模表示为：
[0166][0167][0168]
其中，p
n
为第n个天线阵元的最大允许的发射功率，集合集合
[0169]
步骤4，将多波束卫星系统和速率最大化的鲁棒预编码优化设计问题等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题；
[0170]
具体的，为了能够得到最优鲁棒预编码的解结构，将步骤3的多波束卫星系统和速率最大化的鲁棒预编码优化设计问题p1等价转化为用户信干噪比保障和单天线功率约束下的功率最小化问题p2，即：
[0171][0172][0173][0174]
其中，表示问题p1取得最优解时系统能够实现的最大和速率值。因此当问题p2实现与问题p1相同的最大和速率时，能够保证系统的发送功率等于或小于问题p1所需要的发送功率值，问题p2可以实现与问题p1相同或更优的最优预编码矢量。
[0175]
为了求解问题p2，将最大和速率通过杰森不等式进行近似为：
[0176][0177]
即得到以下的等价关系：
[0178][0179]
其中，γ
k
为第k个用户的目标信干噪比值。
[0180]
此时将优化目标转化为约束条件分别为用户信干噪比保障约束以及单天线功率约束其中，p
n
为第n个天线阵元的最大允许的发射功率，且此时，问题p2可以表示为：
[0181][0182][0183][0184]
其中，用户信干噪比保障约束e{sinr
k
}可近似为：
[0185][0186]
其中目标sinr的约束门限值γ
k
为为为问题p1的能够实现最大和速率，为信道自相关矩阵，且表达式的推导需要对角度误差变量求数学期望。
[0187]
进一步的，假设多波束卫星只能获得各服务用户的不完美的关于离开角 aod的信道信息，不失一般性，且卫星对第k个用户的位置角度估计误差和同时服从某一特定的概率分布，例如均匀分布或高斯分布等，卫星姿态轨道控制而导致的公共角度检测误差和的数值通常在10
‑3rad量级，因此在信道自相关矩阵的计算过程中可以忽略和基于多波束卫星下行信道参数模型，将信道自相关矩阵化简为：
[0188][0189]
其中，平面阵天线相应矢量的第n项可以表示为：
[0190]
[0191]
n＝0,
…
,m
x
m
y
‑
1,
[0192]
其中，信道自相关矩阵中的第m行第n列元素可以表示为：
[0193][0194]
其中，表示信道功率值，[
·
]
m,n
表示取矩阵的第m行第n列的元素，(
·
)
n
表示取矢量的第n个元素。表示对矢量中的第m个元素和第 n个元素的共轭的乘积求数学期望；
[0195]
且：
[0196][0197]
m＝0,
…
,m
x
m
y
‑
1,
[0198]
n＝0,
…
,m
x
m
y
‑
1,
[0199]
其中，a和b分别为第一中间变量和第二中间变量；
[0200]
若卫星对第k个用户的位置角度估计误差和服从均匀分布，即：
[0201][0202][0203]
其中，θ
l
和θ
u
分别表示角度误差取值的上限和下限，和分别表示角度误差取值的上限和下限；
[0204]
上式表示在区间u[θ
l
,θ
u
]服从均匀分布，在区间服从均匀分布，进一步定义：
[0205]
δ
θ
＝θ
u
‑
θ
l
[0206][0207]
其中，δ
θ
和分别表示角度误差和的上下限的差值；
[0208]
则可以得到：
[0209][0210]
将上述表达式中的和进行积分，得到：
[0211][0212]
其中，
[0213][0214][0215][0216]
其中，a、b和z分别为第三中间变量、第四中间变量和第五中间变量。
[0217]
若服从均值为μ
θ,k
、方差为的高斯分布，服从均值为方差为的高斯分布，则和的概率密度函数表示为：
[0218][0219][0220]
其中，表示的概率密度函数，表示的概率密度函数。
[0221]
此时可以表示为：
[0222][0223]
对上述的计算式进行不定积分得到：
[0224][0225]
其中：
[0226][0227][0228][0229][0230][0231][0232]
其中，p、q、d、f、c、e分别为第六中间变量、第七中间变量、第八中间变量、第九中间变量、第十中间变量和第十一中间变量。
[0233]
可以理解的是，本发明中使用的中间变量均无具体含义。
[0234]
步骤5，结合等价优化问题p2的拉格朗日函数及其kkt条件，将鲁棒预编码的最优结构建模为其广义特征值问题的解，获得最优预编码矢量，即原始问题p1的最优解；
[0235]
具体的，等价优化问题p2的拉格朗日函数为：
[0236][0237]
其中，v
k
表示第k个用户的sinr约束的拉格朗日乘子，q
n
对应于第n个单天线功率约束的拉格朗日乘子，e
n
表示第n个元素为1、其他元素为0的矢量，其维度由前后相乘矩阵决定；(
·
)
t
表示对一矢量进行转置操作，则表示对矢量e
n
进行转置，p
n
表示第n个天线阵元的功率上限值。
[0238]
为了简化计算，令则：
[0239][0240]
其中，对拉格朗日函数求导得到：
[0241][0242]
根据kkt条件，最优的预编码矢量需要满足：
[0243][0244]
因此可以得到：
[0245][0246]
令矩阵并且假设最优预编码矢量量为最优预编码矢量，则可以将最优预编码矢量的方向视为对应矩阵对(s
k
,n
k
)最大广义特征值问题的特征向量，最大广义特征值为γ
k
，和γ
k
分别为：
[0247]
γ
k
＝max.generalized eigenvalue(s
k
,n
k
)
[0248][0249]
为了表示方便，本实施例利用表示脚标k的取值范围为从1到k，中脚标n的取值范围为从1到n
t
，当拉格朗日乘子和确定后，最优预编码器的方向和γ
k
可以唯一确定，并且当等价优化问题p3取得最优解时，其sinr约束将取等号，即此时e{sinr
k
}＝γ
k
,据此条件可以求得最优的功率分配，即：
[0250][0251]
其中，ρ为功率分配矢量，且ρ＝[ρ1,
…
,ρ
k
]
t
，矩阵f的维度为k
×
k，其第k 行第i列的元素[f]
k,i
为：
[0252][0253]
此时可计算得到最优预编码矢量
[0254]
进一步的，为了确定拉格朗日乘子和的最优值，本实施例采用次梯度投影技术并结合kkt条件q
n
≥0,和q
n
p
n
＝0,得到关于的迭代表达式q满足：
[0255][0256]
其中，q
i
表示为第i迭代过程中的对应于单天线功率约束的拉格朗日乘子矩阵t
i
表示迭代步长，且t
i
＝1/i。
[0257]
根据等价优化问题p3的拉格朗日函数和其对偶函数，以及两者之间的强对偶关系，可以得到：
[0258][0259]
其中，p
t
为系统实际发送功率。
[0260]
将上述计算式进行重新排列，可以关于的不动点方程式：
[0261][0262]
其中，最优拉格朗日乘子矩阵
[0263]
在实际应用中，目标sinr的约束门限值γ
k
需要尽可能地取最大值，结合下述的凸优化问题确定v
k
的值：
[0264][0265][0266][0267]
其中，p
total
表示系统最大传输功率。
[0268]
为了降低上述优化问题的求解复杂度，将关于v
k
的不动点方程近似为：
[0269][0270]
其中计算辅助因子β
k
＝||h
k
||2，且上述优化问题通过注水算法进行求解。
[0271]
上述优化问题通过注水算法进行求解。
[0272]
步骤6，结合机器学习的方法基于信道自相关矩阵预测优化问题所需的拉格朗日乘子。
[0273]
具体的，所述步骤6还包括以下步骤，
[0274]
步骤6
‑
1，构建卷积神经网络；参照图4的示意，为卷积神经网络的结构示意图，卷积神经网络包括输入层、卷积层、批量标准化层、激活层、平坦层、全连接层和输出层。
[0275]
步骤6
‑
2，收集训练样本并对卷积神经网络进行训练；训练样本包括信道自相关参数矩阵和拉格朗日乘子和通过离线训练的方式对卷积神经网络进行训练。
[0276]
步骤6
‑
3，训练后的神经网络以信道自相关参数矩阵的实部和虚部作为输入，拉格朗日乘子和作为输出；
[0277]
步骤6
‑
4，输入信道自相关参数矩阵通过神经网络，可以直接获得拉格朗日乘子，避免算法实现过程中的迭代过程，降低算法实现复杂度，直接计算满足原始问题的最优预编码矢量。
[0278]
应说明的是，以上所述实施例仅表达了本发明的部分实施方式，其描述并不能理解为对本发明专利范围的限制。应当指出的是，对于本技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干改进，这些均应落入本发明的保护范围。