预导向协方差矩阵重构的宽带自适应波束形成方法及系统技术方案

本发明专利技术属于阵列信号处理中的自适应波束形成技术领域，公开了一种基于预导向协方差矩阵重构的稳健宽带自适应波束形成方法及系统，所述方法包括：设置天线阵列；宽带信号频域处理；使用稳健波束形成准则估计预导向矩阵；通过STMV宽带谱估计得到期望信号功率和噪声功率，在新估计的导向协方差矩阵中去除期望信号部分；通过二次约束二次规划问题的求解得到进一步优化的权值；计算宽带自适应波束形成权矢量。本发明专利技术针对STMV算法中其导向协方差矩阵(STCM)误差和导向矢量误差，设计了一种稳健STMV算法，具有在具有各种误差的宽带信号情况下，波束形成输出的信干噪比更高，期望信号不失真等优点。失真等优点。失真等优点。

【技术实现步骤摘要】
预导向协方差矩阵重构的宽带自适应波束形成方法及系统


[0001]本专利技术属于阵列信号处理中的自适应波束形成
，尤其涉及一种基于预导向协方差矩阵重构的稳健宽带自适应波束形成方法及系统。

技术介绍

[0002]宽带自适应波束形成是一种信号处理技术，用于在接收端将多个天线接收到的信号进行处理，以获得目标信号并抑制干扰信号。该技术基于自适应滤波器和波束形成原理，通过对每个天线接收的信号进行加权和相位调整，从而实现对目标信号的增强和对干扰信号的抑制。它广泛应用于无线通信、雷达、无线电侦听等领域，可提高信号质量和系统性能。
[0003]当处理宽带信号时，由于信号频带宽度较大，信号波形会随着频率的变化而变化，这就使得传统的波束形成技术不能直接应用于宽带信号。与此同时，为了在宽带信号中提取出有效的信息，需要对信号进行高分辨率的采样和处理，这也会增加处理的复杂度和计算量。一种解决方法是使用子带划分的波束形成技术，这种技术将宽带信号分解成多个窄带信号，然后对每个窄带信号进行波束形成。在这种方法中，首先需要将宽带信号分成若干个相互重叠的窄带信号，然后分别对每个窄带信号进行波束形成。在这个过程中，需要确保相邻的窄带信号之间具有一定的重叠，以确保波束形成过程中没有信号损失。然后将各个窄带信号的波束形成结果进行合并，从而得到宽带信号的波束形成结果。另外一种解决方法是基于聚焦的波束成形，此方法使用时域和频域信号处理将波束聚焦到目标，提供高分辨率和低旁瓣水平，同时计算效率高。相干聚焦方法，例如空间重采样或相干信号子空间方法，可以用于在应用窄带波束形成方法之前，在宽带信号的带宽上将窄带相关谱矩阵(CSDM)聚焦到单个频率。空间重采样方法聚焦来自所有角度扇区的源，但其应用限于均匀线性阵列。相干信号子空间方法(CSM)适用于任意阵列，但通常需要DOA的初始估计来估计聚焦矩阵。之后的贝叶斯聚焦方法以及空间插值方法在一定程度上解决了初始角度估计误差的影响。这种方法基于这样的观察，即在由平面波组成的波场中，空间参数和光谱参数是相关的，因为平面波仅通过乘积取决于频率。因此，阵列流形在频率之间的变换等效于阵列的空间插值，其中新的空间坐标按期望的频率比缩放。文献中的大多数聚焦方法要么属于需要DOA先验知识的定向方法，要么属于不需要先验知识的全景方法。有几种聚焦方法试图在定向方法和全景方法之间妥协。
[0004]通过上述分析，现有技术存在的问题及缺陷为：
[0005](1)现有预导向结构的宽带自适应波束形成中，存在频域导向矢量误差以及导向协方差矩阵误差。
[0006](2)现有预导向结构的宽带自适应波束形成中，没有考虑导向协方差矩阵中存在期望信号部分造成的影响。

技术实现思路

[0007]针对现有技术存在的问题，本专利技术提供了一种基于预导向协方差矩阵重构的稳健
宽带自适应波束形成方法及系统。
[0008]本专利技术是这样实现的，一种基于预导向协方差矩阵重构的稳健宽带自适应波束形成方法，所述基于预导向协方差矩阵重构的稳健宽带自适应波束形成方法包括：
[0009]步骤一，设置天线阵列；
[0010]步骤二，宽带信号频域处理；
[0011]步骤三，预导向矩阵估计，使用稳健波束形成准则估计预导向矩阵；
[0012]步骤四，导向协方差矩阵(STCM)重构，通过导向最小方差(STMV)宽带谱估计得到期望信号功率和噪声功率，在新估计的导向协方差矩阵中去除期望信号部分；
[0013]步骤五，优化用于波束形成的导向矢量，通过二次约束二次规划问题的求解得到进一步优化的权值；
[0014]步骤六，计算宽带自适应波束形成权矢量。
[0015]进一步，所述步骤一，设置一个一维均匀线阵，共有M个阵元，阵列间距d＝λ2，设M＝10个全向阵元组成的均匀线阵，干扰是线性调制频率(LFM)信号，并且期望的信号是由滤波的高斯噪声产生的，其中每个频率的噪声功率可能不同，两个相邻传感器之间的距离是最高频率波长的一半，信号带宽默认为30MHz，信号载波中心频率范围默认为25GHz
‑
27GHz，信号传播速度为真空光速，就入射角而言，期望的信号从θ0＝0
°
,一个INR＝30dB的强干扰来自θ1＝40
°
,一个INR＝10dB的弱干扰来自θ2＝
‑
40
°
；期望信号位于角扇区Θ
s
＝[θ0‑3°
,θ0+3
°
]，干扰信号位于Θ
i
＝[θ1‑3°
,θ1+3
°
]∪[θ2‑3°
,θ2+3
°
]，纯噪声扇区为Θ
n
。
[0016]进一步，所述步骤二，包括：
[0017](1)得到阵列接受时域信号X(n)，第m个阵元的时域输出为：
[0018]X
m
(t)＝S0(t
‑
τ
m
(θ0))+N
m
(t)
[0019]其中S0(t)、τ
m
(θ0)、N
m
(t)分别代表期望信号波形，第m个阵元对于对应方向θ0的传播延迟，以及第m个传感器处的噪声加干扰；
[0020](2)将时域快拍X(t)进行L点快速傅里叶变换，生成频域快照，对于第l个FFT单元、第m个阵元和第n个频域快照，可以写出：
[0021][0022]x
m,l
(n)、s
0,l
(n)、n
m,l
(n)分别代表在频率f
l
＝(l
‑
1)/LT
s
处的傅里叶变换，T
s
为采样周期，中心频率信号的阵列导向矢量模型可以写作：
[0023][0024]第n个频域快拍按照信号结构写为：
[0025][0026](3)互谱密度矩阵(CSDM)计算：
[0027][0028](4)导向协方差矩阵(STCM)计算：
[0029]对于任意角度导向协方差矩阵为：
[0030][0031]STCM的结构可以写为：
[0032][0033]其中为期望信号总功率，为第k个干扰信号在第l个频率处的功率，为高斯白噪声总功率，其中T
l
(θ)为第l个频率的预导向矩阵形式为：
[0034][0035]进一步，所述步骤三，采用RCB的方法校正每个频点的期望信号导向矢量，校正预导向矩阵，包括：
[0036](1)在每个频点上解决期望信号导向矢量估计问题：
[0037][0038]其中为先验知识，为假设的期望信号导向矢量；a
l
为要估计的真实期望信号导向矢量；优化问题的最优解的形式为
[0039][0040]其中λ
l
为对应频率处期望信号导向矢量使用由拉格朗日乘子法获得的对角加载因子，λ...

【技术保护点】

【技术特征摘要】
1.一种基于预导向协方差矩阵重构的稳健宽带自适应波束形成方法，其特征在于，通过使用稳健波束形成准则估计预导向矩阵，提高了波束形成的稳健性；通过STMV宽带谱估计得到期望信号功率和噪声功率，在新估计的导向协方差矩阵中去除期望信号部分，进一步优化了预导向矩阵；通过二次约束二次规划问题的求解，得到进一步优化的权值，提高了波束形成的性能；通过计算宽带自适应波束形成权矢量，实现了对宽带信号的有效处理，提高了波束形成的整体性能。2.根据权利要求1所述的基于预导向协方差矩阵重构的稳健宽带自适应波束形成方法，其特征在于，包括：步骤一，设置天线阵列；步骤二，宽带信号频域处理；步骤三，预导向矩阵估计，使用稳健波束形成准则估计预导向矩阵；步骤四，STCM重构，通过STMV宽带谱估计得到期望信号功率和噪声功率，在新估计的导向协方差矩阵中去除期望信号部分；步骤五，优化用于波束形成的导向矢量，通过二次约束二次规划问题的求解得到进一步优化的权值；步骤六，计算宽带自适应波束形成权矢量。3.根据权利要求2所述的基于预导向协方差矩阵重构的稳健宽带自适应波束形成方法，其特征在于，所述步骤一，设置一个一维均匀线阵，共有M个阵元，阵列间距d＝λ/2，设M＝10个全向阵元组成的均匀线阵，干扰是线性调制频率(LFM)信号，并且期望的信号是由滤波的高斯噪声产生的，其中每个频率的噪声功率可能不同，两个相邻传感器之间的距离是最高频率波长的一半，信号带宽默认为30MHz，信号载波中心频率范围默认为25GHz
‑
27GHz，信号传播速度为真空光速，就入射角而言，期望的信号从θ0＝0
°
,一个INR＝30dB的强干扰来自θ1＝40
°
,一个INR＝10dB的弱干扰来自θ2＝
‑
40
°
；期望信号位于角扇区Θ
s
＝[θ0‑3°
,θ0+3
°
]，干扰信号位于Θ
i
＝[θ1‑3°
,θ1+3
°
]∪[θ2‑3°
,θ2+3
°
]，纯噪声扇区为Θ
n
。4.根据权利要求2所述的基于预导向协方差矩阵重构的稳健宽带自适应波束形成方法，其特征在于，所述步骤二，包括：(1)得到阵列接受时域信号X(n)，第m个阵元的时域输出为：X
m
(t)＝S0(t
‑
τ
m
(θ0))+N
m
(t)其中S0(t)、τ
m
(θ0)、N
m
(t)分别代表期望信号波形，第m个阵元对于对应方向θ0的传播延迟，以及第m个传感器处的噪声加干扰；(2)将时域快拍X(t)进行L点快速傅里叶变换，生成频域快照，对于第l个FFT单元、第m个阵元和第n个频域快照，可以写出：x
m,l
(n)、s
0,l
(n)、n
m,l
(n)分别代表在频率f
l
＝(l
‑
1)/LT
s
处的傅里叶变换，T
s
为采样周期，中心频率信号的阵列导向矢量模型可以写作：第n个频域快拍按照信号结构写为：
(3)互谱密度矩阵(CS...

【专利技术属性】
技术研发人员：张捷珂，郑植，滕云龙，
申请(专利权)人：电子科技大学长三角研究院湖州，
类型：发明
国别省市：