专利名称：一种扩展基线解模糊的相位干涉仪测向方法
技术领域：
本发明属于通信雷达技术领域，特别是无线电监测中的相位干涉仪测向方法。
背景技术：
在现有的测向体制中，相位干涉仪测向具有能被动测向，测向精度高、灵敏度高、 实时性好等优点，并且能适用于各种天线阵列形式，因此广泛应用于电子侦察领域的测向系统中。在军事领域，它可对雷达、通信、测控等辐射源进行定位，能在实施针对性电子干扰和军事打击中发挥重要作用；而在民用领域，它可用于交通管制、生命救援、无线电频谱管理等方面。在相位干涉仪测向中，相位的多值性即模糊性是影响测向是否成功的关键问题。 目前常用的解模糊方法主要有长短基线法(见文献雷达对抗原理，赵国庆；西安电子科技大学出版社，1999年)、参差基线法(见文献噪扰条件下多基线相位干涉仪解模糊算法，周亚强；陈翥；皇甫堪；电子与信息学报，2005年2月，第27卷第2期，259-261页)、虚拟基线法(见文献多基线干涉仪测向的基线设计，李建军；电子对抗，2005年第3期，8-11页)、立体基线法(见文献干涉仪测向解模糊方法，司伟建；初萍；应用科技，2007年9月，第34卷第9期，54-57页)、相关法(见文献相关运算在相位干涉仪解模糊中的应用，韩广；王斌；王成；声学技术，2010年10月，第四卷第5期，538-542页)等。长短基线法解相位模糊，方法简单，测向精度高，但要求短基线长度不能大于入射信号波长的一半，这就限制了它在射频频率高端信号中的应用。参差基线法和虚拟基线法虽然解决了波长对短基线物理实现的限制问题，但其单一的摆放形式在天线盘体积受限的情况下影响测向精度。立体基线法不受高频信号波长和摆放形式的限制，但受噪声影响较大(见文献干涉仪测向解模糊方法，司伟建；初萍；应用科技，2007年9月，第34卷第9期，54-57页)。相关法先通过测量值得到一组模糊角度值，反解相位差构成方向矢量并与各阵元输入数据做复相关，比较相关系数完成解模糊，这种方法很好地解决了高频信号对基线长度的限制且抗干扰能力强，同样，相关法也不受阵列摆放形式的限制(见文献相关运算在相位干涉仪解模糊中的应用，韩广； 王斌；王成；声学技术，2010年10月，第四卷第5期，538-542页)，可应用于圆阵测向。该方法应用于相位干涉仪测向，传统做法是选取少数阵元(比如3个)来作为测向阵列，而其余的阵元作为扩展阵元只用于解模糊。这样使得计算入射角度时对阵元信息利用程度低。另外，该方法计算相对复杂，且计算量随模糊程度指数增长。因此，如何充分利用各阵元信息来提高测向精度并且降低解模糊过程中的运算量成为问题的关键。
发明内容
针对现有技术中应用在相位干涉仪测向中的相关法存在对阵元信息利用程度低导致的测向不准确及计算量高的技术问题。因此有必要提供一种相关干涉仪测向方法。本发明公开了一种扩展基线解模糊的相关干涉仪测向方法，其具体包含以下步骤步骤1在圆阵中选取一个等边三角形作为粗测阵列，三角形的其中两条边为两条基线，求这两条基线的模糊相位差W , % ,穷举两基线的模糊相位差的得到if组模糊
相位差( 肩； 1
其中γ=,，为基线长度，]为入射信号波长，丨.丨表示向上取整，
乂」dλL J
4{K - rInkl + ψχ, - '2π 2 + ψ ., kY,k2 e2￡"- ]；
步骤2利用步骤1中的(2JT-1)2组模糊相位差在三角阵中求得￠21-1)2个方向余弦
步骤3将( ,代入圆阵阵元间相位差计算公式求得(2尤-I)2个相位差向； 步骤4求取圆阵各阵元间的实测相位差向量去；
步骤5将步骤3中得到的个相位差向量作为样本，与步骤4中的实测相位差向量#作相关运算；
步骤6选出步骤5中相似度最大的相位差向量样本，记为，计算出实测相位差向量的解模糊值。优选地，所述方法还包括
步骤7通过步骤6计算出来的实测相位差向量的解模糊值求解方向余弦的最小二乘解。优选地，所述方法还包括
步骤8通过步骤7中求得的方向余弦求解角度值。优选地，所述步骤5中相关运算的相似度函数为
M
其中，況为相位差向量长度。优选地，所述步骤6中计算出实测相位差向量的解模糊值公式为
Φ - round ((杀腿-# ) / 2!) *2ττ+ Φ
其中TOl^if (■)表示四舍五入到整数。本发明的有益效果为本发明的方法是利用三角阵作为粗测阵列来求解圆阵模糊，并通过圆阵的If个阵元来进行DOA估计，充分利用到了 i个阵元的信息，进一步提高了测向精度。同时，新算法通过求取方向余弦这个中间量，不但在解模糊过程中减小了计算量，相比于相关法对每一个模糊相位组都要进行一次入射角度值的运算，新算法只计算了一次入射角度值，大大减少了反三角函数这一运算量，提高了测向系统效率。
图1本发明的扩展基线解模糊的相位干涉仪测向方法的流程图。图2 M元均勻圆阵模型及三角阵选取模型。图3现有技术中的相关法与本发明的扩展基线法的方位角精度比较。图4现有技术中的相关法与本发明的扩展基线法的俯仰角精度比较。
具体实施例方式下面结合说明书附图详细说明本发明的具体实施方式
。本发明在相关法解模糊的思想上提供了一种扩展基线解模糊的相位干涉仪测向方法，通过本发明可以提高相关干涉仪二维测向精度，同时减少计算量，提高测向系统的测向性能。本发明主要是在圆阵中选取一个三角阵，首先通过三角阵各阵元之间的相位差求解方向余弦，然后将方向余弦代回圆阵相位差理论计算公式，通过相关运算完成圆阵的解模糊，最后求解入射角度的精确值。为描述方便，首先进行如下定义
相位干涉仪相位干涉仪测向是指根据所测得的电磁波到达各阵元间的空间相位差， 计算信号到达方向的方法。样本库样本库是信号从不同角度入射到天线阵列，在不同阵元间产生的相位差组成的向量的集合。相似度函数相似度函数是用于定量地描述两个向量之间类似程度的映射关系。 该映射的输入为待求相似度的两个向量，输出为相似程度。一种扩展基线解模糊的相位干涉仪测向方法，其流程图如图1所示，其具体包含以下步骤
步骤1在圆阵中选取一个等边三角形作为粗测阵列，三角形的其中两条边为两条基线，求这两条基线的模糊相位差M , ^ ,穷举两基线的模糊相位差的得到组模糊
相位差( ,IU ；
其中Z"= I+.，，为基线长度，。为入射信号波长，丨.丨表示向上取整，
权利要求
1. 一种扩展基线解模糊的相位干涉仪测向方法，其具体包含以下步骤步骤1在圆阵中选取一个等边三角形作为粗测阵列，三角形的其中两条边为两条基线，求这两条基线的模糊相位差^Z1, % ,穷举两基线的模糊相位差的得到(2Γ-if组模糊相位差(‘4」；其中Γ= | + 4 ,,为基线长度，，为入射信号波长，丨.丨表示向上取整， 乂」d L J= 2nki + Ψι, Φλμ = 2^ + ψ , kx,k2 e [-2Κ+1,2‘C-1]；步骤2利用步骤1中的组模糊相位差在三角阵中求得(2r-lf个方向余弦步骤3将(崦，爲)代入圆阵阵元间相位差计算公式求得(2尤-I)2个相位差向量； 步骤4求取圆阵各阵元间的实测相位差向量#；步骤5将步骤3中得到的(2Ζ-if个相位差向量作为样本，与步骤4中的实测相位差向量#作相关运算；步骤6选出步骤5中相似度最大的相位差向量样本，记为#Λ Κ，计算出实测相位差向量的解模糊值。
2.如权利要求1所述的扩展基线解模糊的相位干涉仪测向方法，其特征在所述方法还包括步骤7通过步骤6计算出来的实测相位差向量的解模糊值求解方向余弦的最小二乘解。
3.如权利要求2所述的扩展基线解模糊的相位干涉仪测向方法，其特征在所述方法还包括步骤8通过步骤7中求得的方向余弦求解角度值。
4.权利要求1所述的扩展基线解模糊的相位干涉仪测向方法，其特征在所述步骤5中相关运算的相似度函数为N,s-I其中，F为相位差向量长度。
5.权利要求1所述的扩展基线解模糊的相位干涉仪测向方法，其特征在所述步骤6中计算出实测相位差向量的解模糊值公式为#= round _ # ) / 2/rJ Φ 其中TOl^if (■)表示四舍五入到整数。
全文摘要
本发明属于通信雷达技术领域，本发明在圆阵中选取一个等边三角阵为粗测阵列，计算三角阵中两条基线的相位差，并根据最大模糊程度穷举可能的实际相位差；然后代入圆阵相位差求解公式构造相位差矢量样本库；通过相关运算求得实际圆阵相位差的模糊程度，进一步得到实际圆阵相位差矢量；最后通过求解方向余弦的最小二乘解反解出入射角度值作为精确DOA估计。应用本发明中的一种新的扩展基线解模糊相位干涉仪二维测向算法，能够获得相比于基于相关法的测向算法更高的测向精度，且所需运算量大大减小，是一种高性能的二维测角算法。
文档编号G01S3/46GK102411136SQ20111022658
公开日2012年4月11日 申请日期2011年8月9日 优先权日2011年8月9日
发明者何子述, 吴奉微, 张昕, 李亚星, 程婷, 蒋林鸿, 蒲刚, 贾可新, 邹丁秋, 郑攀 申请人:电子科技大学