本发明属于信号处理技术领域，特别涉及一种自适应波束形成方法，可用于信号抗干扰处理系统。
背景技术：
在通信和雷达领域，自适应天线是及其重要的部件之一。自适应天线能够自动将天线最大辐射方向对准所需信号的来波方向，将天线方向图的波瓣零陷对准干扰方向，具有极为灵活可靠的抗干扰检测能力，是空域上反电子对抗的一种有效措施，其由天线阵列和自适应处理系统构成，故又称自适应阵列。
自适应处理系统的核心是自适应波束形成，自适应波束形成是指根据一定的波束形成准则和算法实时自适应地调整权向量，使得各阵元的接收信号通过加权叠加后，输出信号的质量在所采用的准则下最优使期望信号得到最大的输出功率，零陷对准干扰方向，对干扰信号进行抑制,从而提高阵列输出的信干噪比。
自适应波束形成中的最佳滤波准则有最小均方误差mmse准则、最小方差mv准则、最大信干噪比sinr准则等。其中:
最小均方误差mmse准则使阵列输出信号与期望信号的均方误差最小，该准则不需要已知期望信号的波达方向，但会产生干扰信号。
最小方差准则mv在广义约束条件下使阵列输出方差最小，该准则必须已知期望信号的波达方向。
最大信干噪比sinr准则使有用信号功率与干扰信号及噪声功率比最大，该准则可以获得最大的信噪比，但必须已知噪声统计特性和期望信号的波达方向。
上述准则在自适应波束形成中应用广泛，但当干扰功率起伏时，自适应波束形成的零陷深度不足以完全抑制干扰，使得波束形成器的输出信号中仍含有干扰分量，抗干扰性能急剧下降。
技术实现要素：
本发明的目的在于针对上述现有技术存在的不足，提出一种基于均匀直线阵列零陷加深的自适应波束形成方法，以提高干扰功率起伏时的抗干扰性能。
本发明的技术思路是：通过采用均匀直线阵列为信号的接收模型，对天线阵列接收信号的自相关矩阵进行奇异值分解，增大干扰子空间，提高干扰功率。通过计算改进后的干扰子空间矩阵、期望信号的方向矢量，得到零陷加深后的自适应波束形成最优权系数，其实现步骤包括如下：
1、一种基于均匀直线阵列零陷加深的自适应波束形成方法，其特征在于，包括：
(1)获取均匀天线阵列的阵元数n和阵元间距d，天线阵列接收数据x(t)，获取期望信号的来波方向θ0，干扰信号个数m、干扰信号的来波方向θk，信号源频率f和发射信号波长λ；k＝1,2,...,m；
(2)利用天线阵列接收数据x(t)，计算天线阵列接收数据x(t)的自相关矩阵rxx：
rxx＝e[x(t)·(x(t))h]，
其中，e[·]表示·的数学期望，h表示共轭转置；
(3)利用步骤(1)获取阵列天线的阵元数n和阵元间距d，计算信号源在第i个阵元的时延τi，根据该时延τi，计算期望信号的方向矢量a(θ0)及均匀直线阵的导向矢量矩阵a，i＝1,2,...,n；；
(4)利用矩阵奇异值分解定理，对自相关矩阵rxx进行分解：
其中，σj表示自相关矩阵rxx特征值分解后的第j个特征值，uji表示自相关矩阵rxx特征值分解后的第j个特征向量中第i个矢量，j＝1,2,...,n；
(5)利用分解后的自相关矩阵rxx，构造新的干扰子空间：
(5a)利用分解后的自相关矩阵rxx，得到干扰子空间的m个大特征值λ及对应的特征向量u：
λ＝diag(σ1,σ2,...,σk,...σm)
其中，其中diag(σ1,σ2,...,σk,...σm)表示对角线元素σ1,σ2,...,σk,...,σm的对角矩阵，σk表示干扰子空间的第k个特征值，k＝1,2,...,m。u1i表示rxx特征值分解后的第1个特征向量中的第i个矢量，uki表示rxx特征值分解后的第k个特征向量中的第i个矢量；
(5b)利用干扰子空间的特征值λ，特征向量u，计算干扰子空间新的特征值λ1和特征向量u1：
u1＝p·u
其中，p是可以增强干扰子空间但不改变波束到达方向的参数，p＝α2·aah。α1，α2是两个数值不同的调节因子，取值范围为1.5～5，a是均匀直线阵的导向矢量矩阵；
(5c)利用新的特征值λ1和新的特征向量u1，计算新的干扰子空间矩阵rx：
rx＝u1λ1u1h，
(6)利用新的干扰子空间矩阵rx和期望信号的方向矢量a(θ0)，计算零陷加深后的自适应波束形成最优权系数ω：
其中，θ0表示期望信号的波达方向，rx-1表示rx的逆矩阵。
本发明与现有技术相比，具有如下优点：
1.本发明利用天线阵列的接收数据，对天线阵列接收信号的自相关矩阵进行奇异值分解，增大干扰子空间，提高干扰功率，使得天线方向图的零陷加深，解决了现有技术中由于干扰功率起伏、自适应波束形成的零陷深度不能完全抑制干扰的问题，提高了系统的抗干扰性能。
2.本发明直接将自相关矩阵进行奇异值分解，得到其特征值和特征向量，构建新的干扰子空间，弥补了利用阵列分解自相关矩阵时天线损失的孔径，且不会增加天线阵列排布的复杂度。
附图说明
图1为本发明的实现流程图。
图2为本发明与传统自适应调零方法的对比曲线图。
具体实施方式
以下结合附图，对本发明作进一步详细说明。
参照图1，本发明的实现步骤包括如下：
步骤1，获取均匀天线阵列的相关参数和来波信号的相关参数。
设置来波信号的相关参数，包括干扰信号个数m、第k个信号的波达方向θk，k＝1,2,...,m，期望信号的来波方向θ0、信号源频率f和发射信号波长λ；
根据天线排布情况，获取线性阵列天线的相关参数，包括线阵的阵元个数n和阵元间隔d，天线阵列接收数据x(t)。
步骤2，利用天线阵列接收数据x(t)，计算天线阵列接收数据x(t)的自相关矩阵rxx：
其中，e[·]表示·的数学期望，h表示共轭转置。
步骤3，根据阵列天线的相关参数，计算期望信号的方向矢量a(θ0)及均匀直线阵的导向矢量矩阵a。
3a)利用步骤1获取的阵列天线的阵元间隔d，计算信号源在第i个阵元的时延τi：
其中，c表示光速，θ0表示期望信号的波达方向，i＝1,2,...,n；
3b)利用信号源在第i个阵元的时延τi，计算信号源的方向矢量a(θ0)：
其中，t表示转置，τi表示信号源在第i个阵元的时延，f表示信号源频率，j表示虚数单位；
3c)利用信号源在第i个阵元的时延τi，计算均匀直线阵的导向矢量矩阵a：
a＝[a(θ1),a(θ2),...,a(θk),...,a(θm)],
其中，a(θk)表示第k个信号的方向矢量，
步骤4，利用矩阵奇异值分解定理，对自相关矩阵rxx进行分解如下：
其中，σj表示自相关矩阵rxx特征值分解后的第j个特征值，uji表示自相关矩阵rxx特征值分解后的第j个特征向量中第i个矢量，j＝1,2,...,n。
步骤5，利用分解后的自相关矩阵rxx，构造新的干扰子空间。
5a)利用分解后的自相关矩阵rxx，得到干扰子空间的m个大特征值λ及对应的特征向量u：
λ＝diag(σ1,σ2,...,σk,...σm)
其中，其中diag(σ1,σ2,...,σk,...σm)表示对角线元素σ1,σ2,...,σk,...,σm的对角矩阵，σk表示干扰子空间的第k个特征值，k＝1,2,...,m。u1i表示rxx特征值分解后的第1个特征向量中的第i个矢量，uki表示rxx特征值分解后的第k个特征向量中的第i个矢量；
5b)利用干扰子空间的特征值λ，特征向量u，计算干扰子空间新的特征值λ1和特征向量u1：
u1＝p·u
其中，p是可以增强干扰子空间但不改变波束到达方向的参数，p＝α2·aah，
α1，α2是两个数值不同的调节因子，取值范围为1.5～5，a是均匀直线阵的导向矢量矩阵；
5c)利用新的特征值λ1和新的特征向量u1，计算新的干扰子空间矩阵rx：
rx＝u1λ1u1h。
步骤6，利用新的干扰子空间矩阵rx和期望信号的方向矢量a(θ0)，计算零陷加深后的自适应波束形成最优权系数ω：
其中，θ0表示期望信号的波达方向，rx-1表示rx的逆矩阵。
本发明的效果可以通过以下仿真进一步说明：
1.仿真条件：
设置天线阵是4天线均匀线阵，天线间距是80cm；
设置一路期望信号，其中心频率是1561.098mhz，带宽是4.092mhz，方向角是-60°；设置三路干扰信号，其中心频率分别是1561.098mhz，1200mhz，1600mhz，带宽分别是4.092mhz，2mhz，6mhz，三路信号方向角分别是-20°，30°，60°；天线接收到的四路信号不相干；
设置干噪比是50db，环境中信噪比是-20db。
2.仿真内容
在上述仿真条件下，将4路信号入射到由4根天线组成的均匀线阵，用本发明和传统自适应调零算法对该均匀线阵进行自适应波束形成，结果图2，从图2可见本发明加深了干扰信号的零陷深度，有效抑制干扰信号。