本发明涉及阵列信号处理技术领域，提供一种鲁棒自适应波束形成的优化方法及系统。本发明构建波束形成输出的最大化信干噪比函数，引入采样矩阵，并等效为以优化权重向量为目标的优化问题，求优化问题，获得最优权重向量，实现最大化信干噪比SINR。本发明通过构建最大化信干噪比函数，并引入采样矩阵代替干扰加噪声的协方差矩阵Ri+n，简化计算过程，降低计算量，能够在短时间内求出最优导向矢量解；通过将最大化信干噪比函数在数学上等效为以优化权重向量w为目标的优化问题，使波束形成器的鲁棒性更强，具有更强的抗干扰能力，实现了最大化信号与干扰加噪声比的目标。
1.一种鲁棒自适应波束形成的优化方法，其特征在于：具体步骤如下：步骤S1、构建波束形成输出的最大化信干噪比函数，用于求解最大化信干噪比的问题；
步骤S2、在最大化信干噪比函数中引入采样矩阵 代替干扰加噪声的协方差矩阵Ri+n；
步骤S3、将所述最大化信干噪比函数等效为与采样矩阵 相关的、以优化权重向量w为目标的优化问题；
步骤S4、求解步骤S3中的优化问题，获得最优权重向量w，实现最大化信干噪比SINR；
在步骤S1中，所述最大化信干噪比函数表达式为：其中的w为阵列天线优化设计的波束形成器权重向量，σs为期望信号的功率，a是每个天H
线波达方向角组成的导向矢量，Ri+n为干扰加噪声的协方差矩阵，() 指矩阵的共轭转置；
在步骤S2中，使用采样矩阵 代替干扰加噪声的协方差矩阵Ri+n；
所述采样矩阵 表达式为：其中，T为样本数据的个数，x(k)是阵列在时刻k时的观测矢量，x(k)＝s(k)+i(k)+n(k)，其中s(k)代表期望信号，i(k)代表干扰信号，n(k)为噪声信号；
所述优化问题最优解为 “*”指最优解，w 指w的最优值，其中 由求逆取得；其中，采样矩阵为 权重向量为w、导向矢量为a。
2.根据权利要求1所述的鲁棒自适应波束形成的优化方法，其特征在于：所述步骤S4具体包括：
步骤S4.1、引入阵列输出功率并进行简化，使得求解最大化信干噪比的问题等效于求解最大化阵列输出功率的问题；
步骤S4.2、根据阵列输出功率，构建最大化阵列输出功率函数，用于求解最大化阵列输出功率的问题，并将最大化阵列输出功率函数表示为以优化导向矢量a为目标的优化问题；
步骤S4.3、解步骤S4.2中的优化问题得到最优导向矢量a ，通过a 与 计算出步骤S3*
中的优化问题中的最优权重向量w，实现最大化信干噪比SINR。
3.根据权利要求2所述的鲁棒自适应波束形成的优化方法，其特征在于：在步骤S4.1中，阵列输出功率的表达式为：
其中y(k)指波束形成器在时刻k的输出信号；y(k)＝wx(k)；
对阵列输出功率进行简化，得到阵列输出功率使得求解最大化信干噪比的问题等效于求解最大化阵列输出功率的问题；
N(1‑η1)≤||a||≤N(1+η2)H 2 2
或|aia|≤γi||a||其中，表示采样矩阵，a表示导向矢量，η1为第一用户参数，η2为第二用户参数，N表示*
天线单元个数，a为最优导向矢量，a0是期望信号的导向矢量；
aCa≥Δ1使期望信号的导向矢量a0从期望信号源所在的角扇区Θ中分离出来，也包含了所有干扰信号的导向矢量ai的线性组合，其中：H
C＝∫Θd(θ)d(θ)dθ，C表示阵列几何形状矩阵，其结构由天线阵列几何形状定义；Θ为期望信号源所在的角扇区，Θ＝[θmin，θmax]，d(θ)为当前的角扇区导向矢量，θ为波达方向角；
为约束的基准线，a Ca≥Δ1是为导向矢量a增加不确定集，使用双边的范数约束考虑了阵列响应向量a的增益扰动；其中η1为第一用户参数，η2为第二用户参数，N为天线单元个数；
H 2 2 H 2 2|a0a|≥γ0||a||中，其中γ0为函数在θ0处对应的基准线；|aia|≤γi||a||中，ai是H 2 2 H 2
干扰信号的导向矢量，而γi则为函数在θi处对应的基准线；|aia|≤γi||a||是|a0a|≥2
|a0 a|≥γ0||a|| 来源cos∠(a0，a)≥γ0，为球约束中的一个圆锥状约束；|ai a| ≤
γi||a||来源cos∠(ai，a)≤γi，是|a0a|≥γi||a||的补集，形状为一个球中缺少代表a0的那份圆锥。
4.根据权利要求3所述的鲁棒自适应波束形成的优化方法，其特征在于：步骤S4.3的具体步骤如下：
步骤S4.3.1、输入观测矢量矩阵x(k)，计算采样矩阵步骤S4.3.2、输入角扇区导向矢量d(θ)，计算阵列几何形状矩阵C与约束的基准线Δ1；
步骤S4.3.3、输入采样矩阵 阵列几何形状矩阵C、约束的基准线Δ1，期望信号的导向矢量a0或干扰信号的导向矢量ai、基准线γ、第一用户参数η1、第二用户参数η2、天线单元个数N，利用最大化阵列输出功率函数求解导向矢量a，计得导向矢量a；
步骤S4.3.4、若求得导向矢量a秩为一，该导向矢量a即为最优导向矢量a ；若求得导向*
矢量a秩不为一，对导向矢量a进行降秩，求得秩为一的最优导向矢量a；
步骤S4.3.5、输入最优导向矢量a与采样矩阵 计算出最优权重向量w ；
步骤S4.3.6、输入最优权重向量w与观测矢量矩阵x(k)，计算出波束形成器输出信号y(k)，得到的波束形成器输出信号y(k)拥有最大化信干噪比、最大化阵列输出功率。
5.根据权利要求4所述的鲁棒自适应波束形成的优化方法，其特征在于：步骤S4.2中的优化问题转换为如下第一问题：H
N(1‑η1)≤||a||≤N(1+η2)H 2 2
其中约束|a0a|≥γ0||a||变换为a(γ0I‑a0a0)a≤0，I是单位阵；
a(γ0I‑a0a0)a≤0非凸的二次约束二次规划问题P1等价于如下第二问题：s.t tr(A1Y)≥Δ1N(1‑η1)≤tr(A2Y)≤N(1+η2)tr(A3Y)≤0
其中tr()表示矩阵的迹，rank(Y)＝1表示矩阵秩为一，Y＝aa， A1＝C，A2＝I，H
将第二问题进行半正定松弛并转化为凸问题P2：s.t tr(A1Y)≥Δ1N(1‑η1)≤tr(A2Y)≤N(1+η2)tr(A3Y)≤0
解得最优解Y，若Y秩不为一，通过矩阵降秩算法进行降秩，使得Y秩为一。
6.根据权利要求4所述的鲁棒自适应波束形成的优化方法，其特征在于：步骤S4.2中的优化问题转换为如下第一问题：H
N(1‑η1)≤||a||≤N(1+η2)H 2 2
其中约束|aia|≤γi||a||变换为a(γiI‑aiai)a≥0，I是单位阵；
a(γiI‑aiai)a≥0非凸的二次约束二次规划问题P1等价于如下第二问题：s.t tr(A1Y)≥Δ1N(1‑η1)≤tr(A2Y)≤N(1+η2)tr(A3Y)≥0
其中tr()表示矩阵的迹，rank(Y)＝1表示矩阵秩为一，Y＝aa， A1＝C，A2＝I，A3＝γiI‑aiaiH；
将第二问题等价于如下第三问题：s.t tr(A1Y)≥Δ1N(1‑η1)≤tr(A2Y)≤N(1+η2)tr(A3Y)≥0
解得最优解Y，若Y秩不为一，通过矩阵降秩算法进行降秩，使得Y秩为一。
一种鲁棒自适应波束形成的优化方法及系统
技术领域
[0001] 本发明涉及阵列信号处理技术领域，具体涉及一种鲁棒自适应波束形成的优化方法及系统。
背景技术
[0002] 目前阵列信号处理已广泛应用于雷达、声纳、通信、麦克风阵列语音和音频处理等众多领域，在阵列处理中，因为传统的自适应波束形成技术，如Capon波束形成的优化方法，
对估计的和实际的信号导向矢量之间的微小差异，指向和天线校准误差，以及其他的失配
情况，具有较弱的抗干扰能力。这催生了鲁棒自适应波束形成技术出现，它是一种可以显著
提高阵列输出信干噪比和主瓣宽度和旁瓣电平等性能的方法，而且鲁棒自适应波束形成技
术在失配情况下非常有效。因此研究具有鲁棒性或稳健性的自适应波束形成就具有重大意
义。
[0003] 众所周知，因为实际和假设的导向矢量之间的失配，最小方差无失真响应(MVDR)波束形成器的性能显著下降。为了提升波束形成器的性能，设计时对先验信息的要求越少
越好，因为先验知识通常是不精确和不完整的。所以在估计信号导向矢量时加入各种约束，
可提高波束形成的性能。一般使用波束形成器输出信干噪比和波束形成器输出功率作为评
估性能，越高越好。
[0004] 中国发明专利公开(公告)号CN108462521B(公开日为2018‑08‑28)，公开了一种自适应阵列天线的抗干扰实现方法，旨在提供一种计算量小，实时性强，抗干扰效果更好的抗
干扰方法。该发明通过下述技术方案予以实现：在天线阵元数为M的均匀圆阵天线中，数据
预处理函数将M路数字信号进行N阶时域抽头得到连续的空时二维快拍数据，同时根据上位
机提供的期望信号来向得到空时二维导向矢量；最小方差无失真响应MVDR算法简化实现函
数利用空时二维导向矢量计算导向矢量剩余矩阵，再结合空时二维快拍数据根据其内置的
波束合成函数和权值迭代函数同时进行空时二维波束合成和空时二维权值迭代；阵列输出
函数将波束合成后的数字信号经数模DA转换器转换成模拟信号，模拟信号经上变频函数上
变频输出送给后端处理机。现有技术在实现上述效果的同时，存在抗干扰能力较差与计算
速度较低的缺陷。
发明内容
[0005] 本发明的目的是克服现有技术的上述缺陷，提供一种基于导向矢量估计的、加入了新的约束的、最小方差无失真响应的、鲁棒自适应波束形成的优化方法。
[0006] 为了实现上述目的，本发明提供了一种鲁棒自适应波束形成的优化方法及系统。
[0007] 一种鲁棒自适应波束形成的优化方法，具体步骤如下：
[0008] 步骤S1、构建波束形成输出的最大化信干噪比函数，用于求解最大化信干噪比的问题；
[0009] 步骤S2、在最大化信干噪比函数中引入采样矩阵 代替干扰加噪声的协方差矩阵Ri+n；
[0010] 步骤S3、将所述最大化信干噪比函数等效为与采样矩阵 相关的、以优化权重向量w为目标的优化问题；
[0011] 步骤S4、求解步骤S3中的优化问题，获得最优权重向量w*，实现最大化信干噪比SINR。
[0012] 作为优选方案，所述步骤S4具体包括：
[0013] 步骤S4.1、引入阵列输出功率并进行简化，使得求解最大化信干噪比的问题等效于求解最大化阵列输出功率的问题；
[0014] 步骤S4.2、根据阵列输出功率，构建最大化阵列输出功率函数，用于求解最大化阵列输出功率的问题，并将最大化阵列输出功率函数表示为以优化导向矢量a为目标的优化
问题；
[0015] 步骤S4.3、解步骤S4.2中的优化问题得到最优导向矢量a*，通过a*与 计算出步*
骤S3中的优化问题中的最优权重向量w，实现最大化信干噪比SINR。
[0016] 作为优选方案，在步骤S1中，所述最大化信干噪比函数表达式为：
[0017]
[0018] 其中的w为阵列天线优化设计的波束形成器权重向量，σs为期望信号的功率，a是H
每个天线波达方向角组成的导向矢量，Ri+n为干扰加噪声的协方差矩阵，() 指矩阵的共轭
转置。
[0019] 作为优选方案，在步骤S2中，使用采样矩阵 代替干扰加噪声的协方差矩阵Ri+n；
[0020] 所述采样矩阵 表达式为：
[0021]
[0022] 其中，T为样本数据的个数，x(k)是阵列在时刻k时的观测矢量，x(k)＝s(k)+i(k)+n(k)，其中s(k)代表期望信号，i(k)代表干扰信号，n(k)为噪声信号。
[0023] 作为优选方案，在步骤S3中，所述优化问题表达式为：
[0024]
[0025] s.t wHa＝1
[0026] 所述优化问题最优解为 “*”指最优解，w*指w的最优值，其中由 求逆取得；其中，采样矩阵 权重向量w、导向矢量a。
[0027] 作为优选方案，在步骤S4.1中，阵列输出功率的表达式为：
[0028] E[|y(k)|2]＝E[|w*Hx(k)|2]
[0029] 其中y(k)指波束形成器在时刻k的输出信号；y(k)＝wHx(k)；
[0030] 对阵列输出功率进行简化，得到阵列输出功率使得求解最大化信干噪比的问题等效于
求解最大化阵列输出功率的问题；
[0031] 所述步骤S4.2中的优化问题表达式为：
[0032]
[0033] s.t aHCa≥Δ1
[0034] N(1‑η1)≤||a||2≤N(1+η2)
[0035] |a0Ha|2≥γ0||a||2
[0036] 或|aiHa|2≤γi||a||2
[0037] 其中，采样矩阵 导向矢量a、波达方向角θ、为第一用户参数η1、第二用户参数η2，*
天线单元个数N的优化问题，最优导向矢量a；
[0038] aHCa≥Δ1使期望信号的导向矢量a0从期望信号源所在的角扇区Θ中分离出来，也H
包含了所有干扰信号的导向矢量ai的线性组合，其中C＝∫Θd(θ)d(θ)dθ，其结构由天线阵列
几何形状定义，Θ为期望信号源所在的角扇区，Θ＝[θmin，θmax]，d(θ)为当前的角扇区导向
H
矢量；Δ1＝minθ∈Θd(θ)Cd(θ)，为约束的基准线；
[0039] aHCa≥Δ1是为导向矢量a增加不确定集，使用双边的范数约束考虑了阵列响应向量a的增益扰动；其中η1为第一用户参数，η2为第二用户参数，N为天线单元个数；
[0040] |a0Ha|2≥γ0||a||2中，a0是期望信号的导向矢量，其中γ0为函数在θ0处对应的基H 2 2
准线；|aia|≤γi||a|| 中，ai是干扰信号的导向矢量，而γi则为函数在θi处对应的基准
H 2 2 H 2 2
线；|aia|≤γi||a||是|a0a|≥γi||a||的补集；
[0041] |a0Ha|2≥γ0||a||2来源cos∠(a0,a)≥γ0，为球约束中的一个圆锥状约束；|aiHa2 2 H 2 2
|≤γi||a||来源cos∠(ai,a)≤γi，是|a0a|≥γi||a||的补集，形状为一个球中缺少
代表a0的那份圆锥。
[0042] 作为优选方案，步骤S4.3的具体步骤如下：
[0043] 步骤S4.3.1、输入观测矢量矩阵x(k)，计算采样矩阵 ；
[0044] 步骤S4.3.2、输入角扇区导向矢量d(θ)，计算阵列几何形状矩阵C与约束的基准线Δ1；
[0045] 步骤S4.3.3、输入采样矩阵 阵列几何形状矩阵C、约束的基准线Δ1，期望信号的导向矢量a0或干扰信号的导向矢量ai、基准线γ、第一用户参数η1、第二用户参数η2、天线单
元个数N，利用最大化阵列输出功率函数求解导向矢量a，计得导向矢量a；
[0046] 步骤S4.3.4、若求得导向矢量a秩不为一，该导向矢量a即为最优导向矢量a*；若求*
得导向矢量a秩不为一，对导向矢量a进行降秩，求得秩为一的最优导向矢量a；
[0047] 步骤S4.3.5、输入最优导向矢量a*与采样矩阵 计算出最优权重向量w*；
[0048] 步骤S4.3.6、输入最优权重向量w*与观测矢量矩阵x(k)，计算出波束形成器输出信号y(k)，得到的波束形成器输出信号y(k)拥有最大化信干噪比、最大化阵列输出功率。
[0049] 作为优选方案，步骤S4.2中的优化问题转换为如下问题：
[0050]
[0051] s.t aHCa≥Δ1
[0052] N(1‑η1)≤||a||2≤N(1+η2)
[0053] |a0Ha|2≥γ0||a||2
[0054] 其中约束|a0Ha|2≥γ0||a||2变换为aH(γ0I‑a0a0H)a≤0，I是单位阵；
[0055] 上述问题变为非凸的二次约束二次规划问题：
[0056]
[0057] s.t aHCa≥Δ1
[0058] N(1‑η1)≤||a||2≤N(1+η2)
[0059] aH(γ0I‑a0a0H)a≤0
[0060] 上述问题等价于如下问题：
[0061]
[0062] s.t tr(A1Y)≥Δ1
[0063] N(1‑η1)≤tr(A2Y)≤N(1+η2)
[0064] tr(A3Y)≤0
[0065] rank(Y)＝1
[0066] Y≥0H
[0067] 其中tr()表示矩阵的迹，rank(Y)＝1表示矩阵秩为一，Y＝aa ， A1＝C，H
A2＝I，A3＝γ0I‑a0a0；
[0068] 将问题进行半正定松弛，将问题转化为凸问题：
[0069]
[0070] s.t tr(A1Y)≥Δ1
[0071] N(1‑η1)≤tr(A2Y)≤N(1+η2)
[0072] tr(A3Y)≤0
[0073] Y≥0
[0074] 解得最优解Y*，若Y*秩不为一，通过矩阵降秩算法进行降秩，使得Y*秩为一。
[0075] 作为优选方案，步骤S4.2中的优化问题转换为如下问题：
[0076]
[0077] s.t aHCa≥Δ1
[0078] N(1‑η1)≤||a||2≤N(1+η2)
[0079] |aiHa|2≤γi||a||2
[0080] 其中约束|aiHa|2≤γi||a||2变换为aH(γiI‑aiaiH)a≥0，I是单位阵；
[0081] 上述问题变为非凸的二次约束二次规划问题：
[0082]
[0083] s.t aHCa≥Δ1
[0084] N(1‑η1)≤||a||2≤N(1+η2)
[0085] aH(γiI‑aiaiH)a≥0
[0086] 上述问题等价于如下问题：
[0087]
[0088] s.t tr(A1Y)≥Δ1
[0089] N(1‑η1)≤tr(A2Y)≤N(1+η2)
[0090] tr(A3Y)≥0
[0091] rank(Y)＝1
[0092] Y≥0H
[0093] 其中tr()表示矩阵的迹，rank(Y)＝1表示矩阵秩为一，Y＝aa ， A1＝C，A2H
＝I，A3＝γiI‑aiai；
[0094] 将上述问题等价于如下问题：
[0095]
[0096] s.t tr(A1Y)≥Δ1
[0097] N(1‑η1)≤tr(A2Y)≤N(1+η2)
[0098] tr(A3Y)≥0
[0099] Y≥0
[0100] 解得最优解Y*，若Y*秩不为一，通过矩阵降秩算法进行降秩，使得Y*秩为一。
[0101] 在估计导向矢量的问题中加入了新的约束，可以提高鲁棒自适应波束形成器的性能，使波束形成器的鲁棒性更强，具有更强的抗干扰能力。求解过程计算量小，能够在短时
间内求出最优导向矢量解，实现了最大化信号与干扰加噪声比的目标。
[0102] 此外，本发明提供了一种鲁棒自适应波束形成的系统，包括最大化信干噪比函数、最大化阵列输出功率函数；在最大化信干噪比函数中，求解最大化信干噪比；在最大化阵列
输出功率函数，求解最大化阵列输出功率；简化后，求解最大化信干噪比的问题等效于求解
最大化阵列输出功率的问题。
[0103] 与现有技术相比，本发明的有益效果在于：
[0104] 本发明提供了一种鲁棒自适应波束形成的优化方法，通过构建最大化信干噪比函数，并引入采样矩阵 代替干扰加噪声的协方差矩阵Ri+n，简化计算过程，降低计算量，能够
在短时间内求出最优导向矢量解；通过将最大化信干噪比函数在数学上等效为以优化权重
向量w为目标的优化问题，使波束形成器的鲁棒性更强，具有更强的抗干扰能力，实现了最
大化信号与干扰加噪声比的目标。
[0105] 本发明提供了一种鲁棒自适应波束形成的系统，利用最大化信干噪比函数、最大化阵列输出功率函数，简化后，将求解最大化信干噪比的问题等效于求解最大化阵列输出
功率的问题，提高鲁棒自适应波束形成器的性能，使波束形成器的鲁棒性更强，具有更强的
抗干扰能力，求解过程计算量小，能够在短时间内求出最优导向矢量解，实现了最大化信干
噪比的目标。
附图说明
[0106] 图1是本发明实施例的鲁棒自适应波束形成的优化方法及系统的线形阵列天线结构示意图。
[0107] 图2是本发明实施例的鲁棒自适应波束形成的优化方法及系统的期望信号的导向矢量示意图。
[0108] 图3是本发明实施例的鲁棒自适应波束形成的优化方法及系统的流程示意图。
[0109] 其中，1‑天线；2‑期望信号。
具体实施方式
[0110] 下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
[0111] 在本发明的描述中，需要说明的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为
基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗
示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对
本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗
示相对重要性。
[0112] 在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可
以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是
两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本
发明中的具体含义。
[0113] 此外，在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。
[0114] 具体的，如图1至3所示，本实施例的鲁棒自适应波束形成的优化方法及系统。
[0115] 一种鲁棒自适应波束形成的优化方法，具体步骤如下：
[0116] 步骤S1、构建波束形成输出的最大化信干噪比函数，用于求解最大化信干噪比的问题；
[0117] 步骤S2、在最大化信干噪比函数中引入采样矩阵 代替干扰加噪声的协方差矩阵Ri+n；
[0118] 步骤S3、将所述最大化信干噪比函数等效为与采样矩阵 相关的、以优化权重向量w为目标的优化问题；
[0119] 步骤S4、求解步骤S3中的优化问题，获得最优权重向量w*，实现最大化信干噪比SINR。最大化SINR被称为Capon波束形成。
[0120] 进一步的，步骤S4具体包括：
[0121] 步骤S4.1、引入阵列输出功率并进行简化，使得求解最大化信干噪比的问题等效于求解最大化阵列输出功率的问题；通过天线阵列的几何形状与期望信号的参数，对导向
矢量a进行估计；
[0122] 步骤S4.2、根据阵列输出功率，构建最大化阵列输出功率函数，用于求解最大化阵列输出功率的问题，并将最大化阵列输出功率函数表示为以优化导向矢量a为目标的优化
问题；
* *
[0123] 步骤S4.3、解步骤S4.2中的优化问题得到最优导向矢量a ，通过a 与 计算出步*
骤S3中的优化问题中的最优权重向量w，实现最大化信干噪比SINR。
[0124] 进一步的，在步骤S1中，所述最大化信干噪比函数表达式为：
[0125]
[0126] 其中的w为阵列天线优化设计的波束形成器权重向量，σs为期望信号的功率，a是H
每个天线波达方向角组成的导向矢量，Ri+n为干扰加噪声的协方差矩阵，() 指矩阵的共轭
转置。
[0127] 进一步的，在步骤S2中，使用采样矩阵 代替干扰加噪声的协方差矩阵Ri+n；
[0128] 所述采样矩阵 表达式为：
[0129]
[0130] 其中，T为样本数据的个数，x(k)是阵列在时刻k时的观测矢量，x(k)＝s(k)+i(k)+n(k)，其中s(k)代表期望信号，i(k)代表干扰信号，n(k)为噪声信号。
[0131] 进一步的，在步骤S3中，所述步骤S3中的优化问题表达式为：
[0132]
[0133] s.t wHa＝1
[0134] 步骤S3中的优化问题最优解为 “*”指最优解，w*指w的最优值，其中 由 求逆取得；其中，采样矩阵 权重向量w、导向矢量a。
[0135] 进一步的，在步骤S4.1中，阵列输出功率的表达式为：
[0136] E[|y(k)|2]＝E[|w*Hx(k)|2]
[0137] 其中y(k)指波束形成器在时刻k的输出信号；y(k)＝wHx(k)。
[0138] 对阵列输出功率进行简化，得到阵列输出功率使得求解最大化信干噪比的问题等效于
求解最大化阵列输出功率的问题。
[0139] 进一步的，所述步骤S4.2中的优化问题表达式为：
[0140]
[0141] s.t aHCa≥Δ1
[0142] N(1‑η1)≤||a||2≤N(1+η2)
[0143] |a0Ha|2≥γ0||a||2
[0144] 或|aiHa|2≤γi||a||2
[0145] 其中，采样矩阵 导向矢量a、波达方向角θ、为第一用户参数η1、第二用户参数η2，*
天线单元个数N的优化问题，最优导向矢量a；
[0146] aHCa≥Δ1使期望信号的导向矢量a0从期望信号源所在的角扇区Θ中分离出来，也H
包含了所有干扰信号的导向矢量ai的线性组合，其中C＝∫Θd(θ)d(θ)dθ，其结构由天线阵列
几何形状定义，Θ为期望信号源所在的角扇区，Θ＝[θmin,θmax]，d(θ)为当前的角扇区导向
H
矢量；Δ1＝minθ∈Θd(θ)Cd(θ)，为约束的基准线；
[0147] aHCa≥Δ1是为导向矢量a增加不确定集，使用双边的范数约束考虑了阵列响应向量a的增益扰动；其中η1为第一用户参数，η2为第二用户参数，N为天线单元个数；
[0148] |a0Ha|2≥γ0||a||2中，a0是期望信号的导向矢量，其中γ0为函数在θ0处对应的基H 2 2
准线；|aia|≤γi||a|| 中，ai是干扰信号的导向矢量，而γi则为函数在θi处对应的基准
H 2 2 H 2 2
线；|aia|≤γi||a||是|a0a|≥γi||a||的补集；
[0149] |a0Ha|2≥γ0||a||2来源cos∠(a0,a)≥γ0，为球约束中的一个圆锥状约束；|aiHa2 2 H 2 2
|≤γi||a||来源cos∠(ai,a)≤γi，是|a0a|≥γi||a||的补集，形状为一个球中缺少
代表a0的那份圆锥。
[0150] 进一步的，步骤S4.3的具体步骤如下：
[0151] 步骤S4.3.1、输入观测矢量矩阵x(k)，计算采样矩阵
[0152] 步骤S4.3.2、输入角扇区导向矢量d(θ)，计算阵列几何形状矩阵C与约束的基准线Δ1；
[0153] 步骤S4.3.3、输入采样矩阵 阵列几何形状矩阵C、约束的基准线Δ1，期望信号的导向矢量a0或干扰信号的导向矢量ai、基准线γ、第一用户参数η1、第二用户参数η2、天线单
元个数N，利用最大化阵列输出功率函数求解导向矢量a，计得导向矢量a；
[0154] 步骤S4.3.4、若求得导向矢量a秩不为一，该导向矢量a即为最优导向矢量a*；若求*
得导向矢量a秩不为一，对导向矢量a进行降秩，求得秩为一的最优导向矢量a；
* *
[0155] 步骤S4.3.5、输入最优导向矢量a与采样矩阵 计算出最优权重向量w ；
[0156] 步骤S4.3.6、输入最优权重向量w*与观测矢量矩阵x(k)，计算出波束形成器输出信号y(k)，得到的波束形成器输出信号y(k)拥有最大化信干噪比、最大化阵列输出功率。
[0157] 进一步的，步骤S4.2中的优化问题转换为如下问题：
[0158]
[0159] s.t aHCa≥Δ1
[0160] N(1‑η1)≤||a||2≤N(1+η2)
[0161] |a0Ha|2≥γ0||a||2
[0162] 其中约束|a0Ha|2≥γ0||a||2变换为aH(γ0I‑a0a0H)a≤0，I是单位阵；
[0163] 上述问题变为非凸的二次约束二次规划问题：
[0164]
[0165] s.t aHCa≥Δ1
[0166] N(1‑η1)≤||a||2≤N(1+η2)
[0167] aH(γ0I‑a0a0H)a≤0
[0168] 上述问题等价于如下问题：
[0169]
[0170] s.t tr(A1Y)≥Δ1
[0171] N(1‑η1)≤tr(A2Y)≤N(1+η2)
[0172] tr(A3Y)≤0
[0173] rank(Y)＝1
[0174] Y≥0H
[0175] 其中tr()表示矩阵的迹，rank(Y)＝1表示矩阵秩为一，Y＝aa ， A1＝C，H
A2＝I，A3＝γ0I‑a0a0；
[0176] 将问题进行半正定松弛，将问题转化为凸问题：
[0177]
[0178] s.t tr(A1Y)≥Δ1
[0179] N(1‑η1)≤tr(A2Y)≤N(1+η2)
[0180] tr(A3Y)≤0
[0181] Y≥0
[0182] 解得最优解Y*，若Y*秩不为一，通过矩阵降秩算法进行降秩，使得Y*秩为一。
[0183] 进一步的，步骤S4.2中的优化问题转换为如下问题：
[0184]
[0185] s.t aHCa≥Δ1
[0186] N(1‑η1)≤||a||2≤N(1+η2)
[0187] |aiHa|2≤γi||α||2
[0188] 其中约束|aiHα|2≤γi||α||2变换为αH(γiI‑aiaiH)a≥0，I是单位阵；
[0189] 上述问题变为非凸的二次约束二次规划问题：
[0190]
[0191] s.t aHCa≥Δ1
[0192] N(1‑η1)≤||a||2≤N(1+η2)
[0193] aH(γiI‑aiaiH)a≥0
[0194] 上述问题等价于如下问题：
[0195]
[0196] s.t tr(A1Y)≥Δ1
[0197] N(1‑η1)≤tr(A2Y)≤N(1+η2)
[0198] tr(A3Y)≥0
[0199] rank(Y)＝1
[0200] Y≥0
[0201] 其中tr()表示矩阵的迹，rank(Y)＝1表示矩阵秩为一，Y＝aaH， A1＝C，A2H
＝I，A3＝γiI‑aiai；
[0202] 将上述问题等价于如下问题：
[0203]
[0204] s.t tr(A1Y)≥Δ1
[0205] N(1‑η1)≤tr(A2Y)≤N(1+η2)
[0206] tr(A3Y)≥0
[0207] Y≥0
[0208] 解得最优解Y*，若Y*秩不为一，通过矩阵降秩算法进行降秩，使得Y*秩为一。
[0209] 本实施例采用新的基于导向矢量估计的最小方差无失真响应的鲁棒自适应波束形成的优化方法。这种新的算法主要在估计导向矢量的问题中加入了新的约束，提高鲁棒
自适应波束形成器的性能，使波束形成器的鲁棒性更强，求解过程计算量小，能够在短时间
内求出最优导向矢量解，实现了最大化信号与干扰加噪声比的目标；考虑到期望信号的导
向矢量不能精确估计的情况，在估计导向矢量的问题中加入了新的约束，通过求解导向矢
量，以达到最大化信干噪比的目的。
[0210] 此外，本实施例提供了一种鲁棒自适应波束形成的系统，包括最大化信干噪比函数、最大化阵列输出功率函数；在最大化信干噪比函数中，求解最大化信干噪比；在最大化
阵列输出功率函数，求解最大化阵列输出功率；简化后，求解最大化信干噪比的问题等效于
求解最大化阵列输出功率的问题。
[0211] 本实施例的具体实施过程为：
[0212] 步骤S1、构建波束形成输出的最大化信干噪比函数，用于求解最大化信干噪比的问题；
[0213] 步骤S2、在最大化信干噪比函数中引入采样矩阵 代替干扰加噪声的协方差矩阵Ri+n；
[0214] 步骤S3、将所述最大化信干噪比函数等效为与采样矩阵 相关的、以优化权重向量w为目标的优化问题；
[0215] 步骤S4、求解步骤S3中的优化问题，获得最优权重向量w*，实现最大化信干噪比SINR。
[0216] 综上，本发明实施例提供一种鲁棒自适应波束形成的优化方法及系统：
[0217] 本发明实施例提供了一种鲁棒自适应波束形成的优化方法，通过构建最大化信干噪比函数，并引入采样矩阵 代替干扰加噪声的协方差矩阵Ri+n，简化计算过程，降低计算
量，能够在短时间内求出最优导向矢量解；通过将最大化信干噪比函数在数学上等效为以
优化权重向量w为目标的优化问题，使波束形成器的鲁棒性更强，具有更强的抗干扰能力，
实现了最大化信号与干扰加噪声比的目标。
[0218] 本发明实施例提供了一种鲁棒自适应波束形成的系统，利用最大化信干噪比函数、最大化阵列输出功率函数，简化后，将求解最大化信干噪比的问题等效于求解最大化阵
列输出功率的问题，提高鲁棒自适应波束形成器的性能，使波束形成器的鲁棒性更强，具有
更强的抗干扰能力，求解过程计算量小，能够在短时间内求出最优导向矢量解，实现了最大
化信干噪比的目标。
[0219] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和替换，这些改进和替换
也应视为本发明的保护范围。
