一种被动雷达空间谱估计与合成波束的联合测向方法与流程

本发明涉及一种联合测向方法，特别是一种被动雷达空间谱估计与合成波束的联合测向方法。

背景技术：

目前典型的微波被动导引头常用测向方法有比幅法测向、干涉仪测向、空间谱估计。上述测向方法均基于对电场分布的测量，从而计算出电波到达方向。其中比幅法测向是基于对目标幅度的测量，干涉仪测向是基于对目标相位差的测量，而空间谱估计是基于对目标信号电磁特征的多次测量。空间谱估计是利用天线阵列实现超分辨测向的一种测向体制，但其计算量大，在进行目标大空域搜索时，难以满足测向实时性的要求。

技术实现要素：

本发明目的在于提供一种被动雷达空间谱估计与合成波束的联合测向方法，解决微波被动导引头大空域搜索的实时性问题。

一种被动雷达空间谱估计与合成波束的联合测向方法的具体步骤为：

第一步构建被动雷达空间谱估计与合成波束的联合测向系统

被动雷达空间谱估计与合成波束的联合测向系统，包括：宽带接收天线阵、接收通道、信号处理单元、信号提取模块、信号采集解算模块、空域扫描模块和空间谱估计模块；所述信号提取模块置于接受通道内，所述信号采集解算模块、空域扫描模块和空间谱估计模块置于信号处理单元内。

所述宽带接收天线阵包含n个天线单元，以宽带接收天线阵中心为坐标原点建立坐标系，以第i个天线为例，则第i个天线坐标表示为(xi，yi，zi)；宽带接收天线阵与接收通道用微波电缆连接；接收通道的两个输出端口与信号处理单元的射频电缆和控制电缆连接；信号处理单元的通信接口与弹上惯导设备连接。

信号提取模块的功能为：提取宽带接收天线阵接收信号的幅度和相位；

信号采集解算模块的功能为：采集并解算接收通道的幅度和相位；

空域扫描模块的功能为：进行数字波束合成并进行空域扫描；

空间谱估计模块的功能为：进行空间谱估计计算，获取信号的入射方向。

第二步信号提取模块提取宽带接收天线阵的幅度和相位

信号提取模块根据宽带接收天线阵接收目标辐射源辐射的微波信号，经由接收通道下变频，频率为ω，由此得到宽带接收天线阵接收到的信号频率f0及各天线的幅度和相位。

第三步信号采集解算模块采集并解算接收通道的幅度和相位

信号处理单元采集接收通道输出的模拟信号，并根据计算得到通道的幅度和相位信息。第i个通道的幅度为ei，相位为φi。

第四步空域扫描模块进行数字波束合成并进行空域扫描

空域扫描模块进行数字波束合成和空域扫描处理，设来波方位角为α，来波俯仰角为β，则第i个通道与宽带接收天线阵中心的相位差为：

ψi＝2πf0(xisinβ+yisinαcosβ+zicosαcosβ)/c。

式中：c为光速。

又根据预先的通道功率补偿值，第i个通道预先的通道功率补偿值表示为ei'，得到一组向量如下：

a＝[e1'ψ1,e2'ψ2,…,ei'ψi,…,en'ψn]^t。

根据得到的通道的幅度和相位信息，得到一组向量如下：

b＝[e1φ1,e2φ2,…,eiφi,…,enφn]^t。

分别改变α、β的取值，使α、β的取值在全视场范围内以一定步长扫描，计算空域扫描值：

c＝ab^t

选取c值最大点对应的α、β的取值，记为αm、βm。

第五步空间谱估计模块进行空间谱估计计算并获取信号的入射方向

空间谱估计模块根据得到的通道的幅度和相位信息，得到数据矩阵x如下，其中n为通道数，m为数据点数，

由x得到数据协方差矩阵r，

r＝xx^h；

对r进行特征值分解得到r的噪声子空间导向矢量un；

用a与un采用music算法对方位角αm、俯仰角βm附近空域进行搜索，搜索范围按实际情况确定，公式为：

式中：p为谱峰值。

找出极大值点对应的角度即信号的入射方向，完成测向。

至此，完成了被动雷达空间谱估计与合成波束联合测向。

本发明采用合成波束进行空域粗扫，空间谱估计在粗扫的结果上进行精扫，比单纯采用谱估计进行全空域扫描的计算量降低了40％以上，能够在1ms内实现大空域搜索功能。

附图说明

图1一种被动雷达空间谱估计与合成波束联合测向方法所述系统组成示意图；

图2一种被动雷达空间谱估计与合成波束联合测向方法所述宽带接收天线阵示意图。

1.宽带接收天线阵2.接收通道3.信号处理单元4.信号提取模块5.信号采集解算模块6.空域扫描模块7.空间谱估计模块

具体实施方式

一种被动雷达空间谱估计与合成波束的联合测向方法的具体步骤为：

第一步构建被动雷达空间谱估计与合成波束的联合测向系统

被动雷达空间谱估计与合成波束的联合测向系统，包括：宽带接收天线阵1、接收通道2、信号处理单元3、信号提取模块、信号采集解算模块5、空域扫描模块6和空间谱估计模块7；所述信号提取模块置于接受通道内，所述信号采集解算模块5、空域扫描模块6和空间谱估计模块7置于信号处理单元3内。

所述宽带接收天线阵1包含n个天线单元，以宽带接收天线阵1中心为坐标原点建立坐标系，以第i个天线为例，则第i个天线坐标表示为xi，yi，zi；宽带接收天线阵1与接收通道2用微波电缆连接；接收通道2的两个输出端口与信号处理单元3的射频电缆和控制电缆连接；信号处理单元3的通信接口与弹上惯导设备连接。

信号提取模块的功能为：提取宽带接收天线阵1接收信号的幅度和相位；

信号采集解算模块5的功能为：采集并解算接收通道2的幅度和相位；

空域扫描模块6的功能为：进行数字波束合成并进行空域扫描；

空间谱估计模块7的功能为：进行空间谱估计计算，获取信号的入射方向。

第二步信号提取模块4提取宽带接收天线阵1的幅度和相位

信号提取模块根据宽带接收天线阵1接收目标辐射源辐射的微波信号，经由接收通道2下变频，频率为ω，由此得到宽带接收天线阵1接收到的信号频率f0及各天线的幅度和相位。

第三步信号采集解算模块5采集并解算接收通道2的幅度和相位

信号处理单元3采集接收通道2输出的模拟信号，并根据计算得到通道的幅度和相位信息。第i个通道的幅度为ei，相位为φi。

第四步空域扫描模块6进行数字波束合成并进行空域扫描

空域扫描模块6进行数字波束合成和空域扫描处理，设来波方位角为α，来波俯仰角为β，则第i个通道与宽带接收天线阵1中心的相位差为：

ψi＝2πf0(xisinβ+yisinαcosβ+zicosαcosβ)/c。

式中：c为光速。

又根据预先的通道功率补偿值，第i个通道预先的通道功率补偿值表示为ei'，得到一组向量如下：

a＝[e1'ψ1,e2'ψ2,…,ei'ψi,…,en'ψn]^t。

根据得到的通道的幅度和相位信息，得到一组向量如下：

b＝[e1φ1,e2φ2,…,eiφi,…,enφn]^t。

分别改变α、β的取值，使α、β的取值在全视场范围内以一定步长扫描，计算空域扫描值：

c＝ab^t

选取c值最大点对应的α、β的取值，记为αm、βm。

第五步空间谱估计模块7进行空间谱估计计算并获取信号的入射方向

空间谱估计模块7根据得到的通道的幅度和相位信息，得到数据矩阵x如下，其中n为通道数，m为数据点数，

由x得到数据协方差矩阵r，

r＝xx^h；

对r进行特征值分解得到r的噪声子空间导向矢量un；

用a与un采用music算法对方位角αm、俯仰角βm附近空域进行搜索，搜索范围按实际情况确定，公式为：

式中：p为谱峰值。

找出极大值点对应的角度即信号的入射方向，完成测向。

至此，完成了被动雷达空间谱估计与合成波束联合测向。