1.本发明涉及目标定位技术领域，尤其涉及一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法。
背景技术：2.众所周知，目标定位技术已广泛应用于无线通信、智慧城市、自动驾驶、导航遥测、紧急救助、安全管理等诸多工业与信息技术领域，其同时也是目标监测、态势感知、战略预警等国防安全领域中不可或缺的支撑技术。依据观测平台的数量可以将目标定位系统划分为单平台定位系统和多平台定位系统两大类，其中单平台定位系统无需同步和信息传输，具有更高的灵活性。
3.近些年来，直接定位(dpd—direct position determination)技术受到了国内外学者的广泛青睐[a.amar,a.j.weiss.localization of narrowband radio emittersbased on doppler frequency shifts[j].ieee transactions on signal processing,2008, 56(11):5500-5508.][t.tirer,a.j.weiss.high resolution localization of narrowbandradio emitters based on doppler frequency shifts[j].signal processing,2017,141(12): 288-298.]，已经成为一个重要的研究热点。人们通常熟知的无线定位技术需要两个计算步骤(即两步定位)：第1步是从传感器所截获到的目标信号中提取用于定位的中间参数，其中涉及空域、时域、频域以及能量域等参数；第2步则利用这些中间参数通过数据融合解算出目标的位置信息。与之不同的是，直接定位技术省略了中间环节，该技术是从原始采样信号中依据某种准则直接获得目标位置信息。依据信息论可知，从原始数据到最终估计结果之间，每增加一步处理就会引入一定的不确定性，从而导致信息损失。因此，直接定位能够比两步定位获得更高的估计精度。根据现有的研究成果可知，直接定位技术具有以下4点优势： (1)直接定位技术无需估计中间参数，因此不用对中间参数估计误差的统计特性进行精确评估；(2)直接定位技术可以利用不同观测阵列的采集信号来自相同目标这一先验信息，更有利于提高定位精度；(3)直接定位技术更便于利用信号波形信息，从而大幅提升目标定位精度；(4)直接定位技术能够避免“目标
─
量测”数据关联问题。
[0004]
由于直接定位方法是在信号域完成目标定位，因此其计算量相对较大，不利于实时定位。另一方面，基于单运动平台的无线定位难以获得较大的阵列孔径，而阵列孔径大小直接影响着定位性能。
技术实现要素：[0005]
本发明针对现有直接定位方法存在的实时定位计算量大、定位性能差的问题，针对小阵列孔径单运动观测平台定位场景，提出一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，其不仅能够提高在小阵列孔径条件下的定位精度，还能够极大降低实
时定位的计算量。
[0006]
为了实现上述目的，本发明首先，在定位区域内划分两维网格点，依次在每个网格点处放置信号源，并利用一个具有相同航迹和速度的运动阵列接收该信号。然后，针对每个网格点处的信号，该运动阵列利用多个时隙构造被动合成孔径阵列输出，并基于此形成被动合成孔径阵列输出协方差矩阵，通过特征值分解得到首1特征矩阵样本集。随后，针对待定位窄带信号源，利用相同航迹和速度的运动阵列以被动合成孔径的方式接收该信号，并计算其阵列输出协方差矩阵，得到首1特征矩阵。接着，将待定位窄带信号源对应的首1特征矩阵与前面构造的首 1特征矩阵样本集中的每个矩阵进行复相关运算，取相关值最大的9个点进行曲面拟合，从而形成二次曲面。最后，求解该二次曲面最大值对应的坐标，并将其作为最终定位结果。本发明的基于单运动阵列被动合成孔径的相关干涉仪直接定位方法的具体实施步骤如下：
[0007]
步骤1：在感兴趣的定位区域内均匀划分若干个离散点，共形成n
x
×
ny个两维网格点；
[0008]
步骤2：依次在步骤1中每个两维网格点对应的坐标处放置一个窄带信号源，并利用一个移动阵列天线在k个短时隙以被动合成孔径的方式接收辐射信号，针对每个两维网格点，所述移动阵列天线按照相同的航迹与速度行驶；
[0009]
步骤3：依次针对步骤1中的每个两维网格点，分别在k个时隙构造被动合成孔径阵列输出；
[0010]
步骤4：针对步骤1中的所有两维网格点，分别在k个时隙构造被动合成孔径阵列输出协方差矩阵，并形成首1特征矩阵样本集
[0011]
步骤5：针对待定位窄带信号源，利用与步骤2相同的移动阵列天线在k个短时隙以被动合成孔径的方式接收辐射信号，该移动阵列天线按照与步骤2相同的航迹和速度行驶；
[0012]
步骤6：针对待定位窄带信号源，分别在k个时隙构造被动合成孔径阵列输出；
[0013]
步骤7：针对待定位窄带信号源，分别在k个时隙构造被动合成孔径阵列输出协方差矩阵，并形成首1特征矩阵
[0014]
步骤8：将步骤7中的首1特征矩阵与步骤4中的首1特征矩阵样本集中的每个矩阵做复相关运算，得到相关值
[0015]
步骤9：取相关值最大值附近的9个点进行二次曲面拟合；
[0016]
步骤10：基于最小二乘二次曲面拟合参数给出最终定位结果。
[0017]
进一步地，所述步骤1中，在感兴趣的定位区域内均匀划分若干个离散点，共形成n
x
×
ny个两维网格点；假设定位区域在x轴坐标的区间范围为[x
min
,x
max
]，并以
△
x作为取样间隔，共有n
x
个离散点，其中第n1个离散点的x轴坐标为定位区域在y轴坐标的区间范围为[y
min
,y
max
]，并以
△
y作为取样间隔，共有ny个离散点，其中第n2个离散点的y轴坐标为因此，代表一个两维网格点，共有n
x
×
ny个两维网格点。
[0018]
进一步地，所述步骤2中，依次在步骤1中每个两维网格点对应的坐标处放置一个
窄带信号源，并利用一个移动阵列天线在k个短时隙以被动合成孔径的方式接收辐射信号，针对每个两维网格点，该运动阵列按照相同的航迹与速度行驶。
[0019]
进一步地，所述步骤3中，针对网格点将第k个短时隙内的运动阵列输出进行数字下变频处理，得到如下输出
[0020][0021]
式中表示在第k个短时隙内的阵列流形向量；表示信号源位置向量；sk(t)表示在第k个短时隙内的信号复包络；表示在第k 个短时隙内的阵列加性噪声；表示在第k个短时隙内的多普勒频移，其表达式为
[0022][0023]
其中pk表示运动阵列在第k个短时隙内的位置向量；表示运动阵列在第 k个短时隙内的速度向量；c表示信号传播速度；fo表示信号中心频率；
[0024]
为了得到被动合成孔径阵列输出，在第k个短时隙内对每个阵元共产生m 个延时输出，延时量分别为τ1,τ2,
…
,τm，其中第m个延时的阵列输出可以表示为
[0025][0026]
式中表示第m个延时的阵列加性噪声；
[0027]
最后将m个延时的阵列输出进行合并，可得第k个短时隙内的被动合成孔径阵列输出，其表达式为
[0028][0029]
式中表示在第k 个短时隙内的被动合成孔径阵列加性噪声；表示在第k个短时隙内的被动合成孔径阵列流形向量，其中
[0030]
进一步地，所述步骤4中，针对网格点构造第k个短时隙的被动合成孔径阵列输出协方差矩阵如下式所示
[0031][0032]
式中t
l
表示第l个采样时间；l表示信号采样数量；表示第k个短时隙内第l个采样时间的被动合成孔径阵列输出；
[0033]
对矩阵进行特征值分解，其最大特征值对应的首1特征向量(即第1 个元素等于1的特征向量)为于是形成如下首1特征矩阵样本集
[0034][0035]
进一步地，所述步骤5中，针对待定位窄带信号源，利用与步骤2相同的移动阵列天线在k个短时隙以被动合成孔径的方式接收辐射信号，该运动阵列按照与步骤2相同的航迹和速度行驶。
[0036]
进一步地，所述步骤6中，针对待定位窄带信号源，将第k个短时隙内的运动阵列输出进行数字下变频处理，得到如下输出
[0037][0038]
式中表示在第k个短时隙内针对待定位信号源的阵列流形向量；表示待定位信号源的位置向量；表示在第k个短时隙内待定位信号复包络；表示在第k个短时隙内针对待定位信号源的阵列加性噪声；表示在第k个短时隙内针对待定位信号源的多普勒频移，其表达式为
[0039][0040]
其中pk表示运动阵列在第k个短时隙内的位置向量；表示运动阵列在第 k个短时隙内的速度向量；c表示信号传播速度；fo表示信号中心频率；
[0041]
为了得到被动合成孔径阵列输出，在第k个短时隙内对每个阵元共产生m 个延时输出，延时量分别为τ1,τ2,
…
,τm，其中第m个延时的阵列输出可以表示为
[0042][0043]
式中表示针对待定位信号源的第m个延时的阵列加性噪声；
[0044]
最后将m个延时的阵列输出进行合并，可得第k个短时隙内的针对待定位信号源的被动合成孔径阵列输出，其表达式为
[0045][0046]
式中表示第k个短时隙内的针对待定位信号源的被动合成孔径阵列加性噪声；表示第 k个短时隙内的针对待定位信号源的被动合成孔径阵列流形向量，其中
[0047]
进一步地，所述步骤7中，针对待定位窄带信号源，构造第k个短时隙的被动合成孔径阵列输出协方差矩阵如下式所示
[0048][0049]
其中t
l
表示第l个采样时间；l表示信号采样数量；表示第k个短时隙内第l个采样时间针对待定位信号源的被动合成孔径阵列输出；
[0050]
对矩阵进行特征值分解，其最大特征值对应的首1特征向量(即第1个元素等于1的特征向量)为于是形成如下首1特征矩阵
[0051][0052]
进一步地，所述步骤8中，将步骤7中的首1特征矩阵与步骤4中的首1 特征矩阵样本集中的每个矩阵做复相关运算，得到相关值其表达式为
[0053][0054]
进一步地，所述步骤9中，取相关值最大值附近的9个点进行曲面拟合，令相关值最大值对应的序号n1和n2分别为n
x
和ny，则取，则取9个点进行二次曲面拟合，并按照下式构造矩阵z和向量ρ
[0055][0056]
式中
[0057][0058]
于是最小二乘二次曲面拟合参数为
[0059]
θ
ls
＝(z
t
z)-1zt
ρ。
[0060]
进一步地，所述步骤10中，基于最小二乘二次曲面拟合参数给出最终定位结果(记为)，其表达式为
[0061][0062]
式中《θ
ls
》2、《θ
ls
》3、《θ
ls
》4、《θ
ls
》5、《θ
ls
》6分别表示最小二乘二次曲面拟合参数θ
ls
中的第2、3、4、5、6个元素。
[0063]
与现有技术相比，本发明具有的有益效果：
[0064]
本发明针对小阵列孔径单运动观测平台定位场景，提出了一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法。由于直接定位方法是在信号域完成目标定位，因此其计算量相对较大，不利于实时定位。针对该问题，本发明首次将传统的相关干涉仪测向技术引入到直接定位问题中，可在提高定位精度的前提下，有效降低实时计算的复杂度。另一方面，基于单运动平台的无线定位难以获得较大的阵列孔径，而阵列孔径大小直接影响着定位性能，对此，本发明融入了被动合成孔径的思想，其中基于多普勒信息将信号的时间增益转化为空间增益，从而得到等效的大阵列孔径，能在小阵列孔径的条件下有效提高定位精度。综上所述，本发明首次将传统的相关干涉仪测向技术引入到直接定位问题中，能够同时具有相关干涉仪测向体制的优势(实时计算量低)和直接定位体制的优势(可避免信息损失)，并且其中融入了被动合成孔径的思想，从而进一步提高了在小阵列孔径条件下的无线定位精度。
附图说明
[0065]
图1是本发明实施例基于单运动阵列被动合成孔径的相关干涉仪直接定位方法原理框图；
[0066]
图2是定位均方根误差随着信噪比的变化曲线；
[0067]
图3是定位均方根误差随着每个时隙内的信号采样点数的变化曲线；
[0068]
图4是定位均方根误差随着信号源频率的变化曲线。
具体实施方式
[0069]
下面结合附图和具体的实施例对本发明做进一步的解释说明：
[0070]
如图1所示，一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，包括：
[0071]
步骤1：在感兴趣的定位区域内均匀划分若干个离散点，共形成n
x
×
ny个两维网格点；
[0072]
步骤2：依次在步骤1中每个两维网格点对应的坐标处放置一个窄带信号源，并利用一个移动阵列天线在k个短时隙以被动合成孔径的方式接收辐射信号，针对每个两维网格点，所述移动阵列天线按照相同的航迹与速度行驶；
[0073]
步骤3：依次针对步骤1中的每个两维网格点，分别在k个时隙构造被动合成孔径阵列输出；
[0074]
步骤4：针对步骤1中的所有两维网格点，分别在k个时隙构造被动合成孔径阵列输出协方差矩阵，并形成首1特征矩阵样本集
[0075]
步骤5：针对待定位窄带信号源，利用与步骤2相同的移动阵列天线在k个短时隙以
被动合成孔径的方式接收辐射信号，该移动阵列天线按照与步骤2相同的航迹和速度行驶；
[0076]
步骤6：针对待定位窄带信号源，分别在k个时隙构造被动合成孔径阵列输出；
[0077]
步骤7：针对待定位窄带信号源，分别在k个时隙构造被动合成孔径阵列输出协方差矩阵，并形成首1特征矩阵
[0078]
步骤8：将步骤7中的首1特征矩阵与步骤4中的首1特征矩阵样本集中的每个矩阵做复相关运算，得到相关值
[0079]
步骤9：取相关值最大值附近的9个点进行二次曲面拟合；
[0080]
步骤10：基于最小二乘二次曲面拟合参数给出最终定位结果。
[0081]
进一步地，所述步骤1中，在感兴趣的定位区域内均匀划分若干个离散点，共形成n
x
×
ny个两维网格点；假设定位区域在x轴坐标的区间范围为[x
min
,x
max
]，并以
△
x作为取样间隔，共有n
x
个离散点，其中第n1个离散点的x轴坐标为定位区域在y轴坐标的区间范围为[y
min
,y
max
]，并以
△
y作为取样间隔，共有ny个离散点，其中第n2个离散点的y轴坐标为因此，代表一个两维网格点，共有n
x
×
ny个两维网格点。
[0082]
进一步地，所述步骤2中，依次在步骤1中每个两维网格点对应的坐标处放置一个窄带信号源，并利用一个移动阵列天线在k个短时隙以被动合成孔径的方式接收辐射信号，针对每个两维网格点，该运动阵列按照相同的航迹与速度行驶。
[0083]
进一步地，所述步骤3中，针对网格点将第k个短时隙内的运动阵列输出进行数字下变频处理，得到如下输出
[0084][0085]
式中表示在第k个短时隙内的阵列流形向量；表示信号源位置向量；sk(t)表示在第k个短时隙内的信号复包络；表示在第k 个短时隙内的阵列加性噪声；表示在第k个短时隙内的多普勒频移，其表达式为
[0086][0087]
其中pk表示运动阵列在第k个短时隙内的位置向量；表示运动阵列在第 k个短时隙内的速度向量；c表示信号传播速度；fo表示信号中心频率；
[0088]
为了得到被动合成孔径阵列输出，在第k个短时隙内对每个阵元共产生m 个延时输出，延时量分别为τ1,τ2,
…
,τm，其中第m个延时的阵列输出可以表示为
[0089][0090]
式中表示第m个延时的阵列加性噪声；
[0091]
最后将m个延时的阵列输出进行合并，可得第k个短时隙内的被动合成孔径阵列输出，其表达式为
[0092][0093]
式中表示在第k 个短时隙内的被动合成孔径阵列加性噪声；表示在第k个短时隙内的被动合成孔径阵列流形向量，其中
[0094]
进一步地，所述步骤4中，针对网格点构造第k个短时隙的被动合成孔径阵列输出协方差矩阵如下式所示
[0095][0096]
式中t
l
表示第l个采样时间；l表示信号采样数量；表示第k个短时隙内第l个采样时间的被动合成孔径阵列输出；
[0097]
对矩阵进行特征值分解，其最大特征值对应的首1特征向量(即第1 个元素等于1的特征向量)为于是形成如下首1特征矩阵样本集
[0098][0099]
进一步地，所述步骤5中，针对待定位窄带信号源，利用与步骤2相同的移动阵列天线在k个短时隙以被动合成孔径的方式接收辐射信号，该运动阵列按照与步骤2相同的航迹和速度行驶。
[0100]
进一步地，所述步骤6中，针对待定位窄带信号源，将第k个短时隙内的运动阵列输出进行数字下变频处理，得到如下输出
[0101][0102]
式中表示在第k个短时隙内针对待定位信号源的阵列流形向量；表示待定位信号源的位置向量；表示在第k个短时隙内待定位信号复包络；表示在第k个短时隙内针对待定位信号源的阵列加性噪声；表示在第k个短时隙内针对待定位信号源的多普勒频移，其表达式为
[0103][0104]
其中pk表示运动阵列在第k个短时隙内的位置向量；表示运动阵列在第 k个短时隙内的速度向量；c表示信号传播速度；fo表示信号中心频率；
[0105]
为了得到被动合成孔径阵列输出，在第k个短时隙内对每个阵元共产生m 个延时输出，延时量分别为τ1,τ2,
…
,τm，其中第m个延时的阵列输出可以表示为
[0106][0107]
式中表示针对待定位信号源的第m个延时的阵列加性噪声；
[0108]
最后将m个延时的阵列输出进行合并，可得第k个短时隙内的针对待定位信号源的被动合成孔径阵列输出，其表达式为
[0109][0110]
式中表示第k个短时隙内的针对待定位信号源的被动合成孔径阵列加性噪声；表示第 k个短时隙内的针对待定位信号源的被动合成孔径阵列流形向量，其中
[0111]
进一步地，所述步骤7中，针对待定位窄带信号源，构造第k个短时隙的被动合成孔径阵列输出协方差矩阵如下式所示
[0112][0113]
其中t
l
表示第l个采样时间；l表示信号采样数量；表示第k个短时隙内第l个采样时间针对待定位信号源的被动合成孔径阵列输出；
[0114]
对矩阵进行特征值分解，其最大特征值对应的首1特征向量(即第1个元素等于1的特征向量)为于是形成如下首1特征矩阵
[0115][0116]
进一步地，所述步骤8中，将步骤7中的首1特征矩阵与步骤4中的首1 特征矩阵样本集中的每个矩阵做复相关运算，得到相关值其表达式为
[0117][0118]
进一步地，所述步骤9中，取相关值最大值附近的9个点进行曲面拟合，令相关值最大值对应的序号n1和n2分别为n
x
和ny，则取，则取9个点进行二次曲面拟合，并按照下式构造矩阵z和向量ρ
[0119][0120]
式中
[0121][0122]
于是最小二乘二次曲面拟合参数为
[0123]
θ
ls
＝(z
t
z)-1zt
ρ。
[0124]
进一步地，所述步骤10中，基于最小二乘二次曲面拟合参数给出最终定位结果(记为)，其表达式为
[0125][0126]
式中《θ
ls
》2、《θ
ls
》3、《θ
ls
》4、《θ
ls
》5、《θ
ls
》6分别表示最小二乘二次曲面拟合参数θ
ls
中的第2、3、4、5、6个元素。
[0127]
为验证本发明效果，考虑基于单运动观测平台的定位场景，在区域 [-2km 2km]
×
[-2km 2km]内以0.02km为间隔划分两维网格点，运动阵列为8 元均匀圆阵，半径为3m，信号为窄带调频信号，采样频率为50khz，运动阵列的初始位置向量为[7.5 8.2]
t
(km)，运动速度向量为[-0.1
ꢀ‑
0.08]
t
(km/s)，运动阵列利用k＝11个时隙对信号源进行定位，每个时隙之间的时间为4s。待定位信号源的位置向量为[0.1
ꢀ‑
0.1]
t
(km)。首先将信号源频率固定为200mhz，每个时隙内的信号采样点数固定为500，图2给出了定位均方根误差随着信噪比的变化曲线；然后将信号源频率固定为200mhz，信噪比固定为-2db，图3给出了定位均方根误差随着每个时隙内的信号采样点数的变化曲线；最后将信噪比固定为-2db，每个时隙内的信号采样点数固定为500，图4给出了定位均方根误差随着信号源频率的变化曲线。
[0128]
从图2至图4中可以看出：相比于未进行被动合成孔径的两步定位方法，新方法通过被动合成孔径有效提升了定位精度。此外，新方法同时具有相关干涉仪测向体制的优势(实时计算量低)和直接定位体制的优势(可避免信息损失)。
[0129]
综上，本发明针对小阵列孔径单运动观测平台定位场景，提出了一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法。由于直接定位方法是在信号域完成目标定位，因此其计算量相对较大，不利于实时定位。针对该问题，本发明首次将传统的相关干涉仪测向技术引入到直接定位问题中，可在提高定位精度的前提下，有效降低实时计算的复杂度。另一方面，基于单运动平台的无线定位难以获得较大的阵列孔径，而阵列孔径大小直
接影响着定位性能，对此，本发明融入了被动合成孔径的思想，其中基于多普勒信息将信号的时间增益转化为空间增益，从而得到等效的大阵列孔径，能在小阵列孔径的条件下有效提高定位精度。综上所述，本发明首次将传统的相关干涉仪测向技术引入到直接定位问题中，能够同时具有相关干涉仪测向体制的优势(实时计算量低)和直接定位体制的优势 (可避免信息损失)，并且其中融入了被动合成孔径的思想，从而进一步提高了在小阵列孔径条件下的无线定位精度。
[0130]
以上所示仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
技术特征：1.一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，其特征在于，包括：步骤1：在定位区域内均匀划分若干个离散点，共形成n
x
×
n
y
个两维网格点；步骤2：依次在步骤1中每个两维网格点对应的坐标处放置一个窄带信号源，并利用一个移动阵列天线在k个短时隙以被动合成孔径的方式接收辐射信号，针对每个两维网格点，所述移动阵列天线按照相同的航迹与速度行驶；步骤3：依次针对步骤1中的每个两维网格点，分别在k个时隙构造被动合成孔径阵列输出；步骤4：针对步骤1中的所有两维网格点，分别在k个时隙构造被动合成孔径阵列输出协方差矩阵，并形成首1特征矩阵样本集步骤5：针对待定位窄带信号源，利用与步骤2相同的移动阵列天线在k个短时隙以被动合成孔径的方式接收辐射信号，该移动阵列天线按照与步骤2相同的航迹和速度行驶；步骤6：针对待定位窄带信号源，分别在k个时隙构造被动合成孔径阵列输出；步骤7：针对待定位窄带信号源，分别在k个时隙构造被动合成孔径阵列输出协方差矩阵，并形成首1特征矩阵步骤8：将步骤7中的首1特征矩阵与步骤4中的首1特征矩阵样本集中的每个矩阵做复相关运算，得到相关值步骤9：取相关值最大值附近的9个点进行二次曲面拟合；步骤10：基于最小二乘二次曲面拟合参数给出最终定位结果。2.根据权利要求1所述的一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，其特征在于，所述步骤3中，针对两维网格点将第k个短时隙内的运动阵列输出进行数字下变频处理，得到如下输出式中表示在第k个短时隙内的阵列流形向量；表示信号源位置向量；s
k
(t)表示在第k个短时隙内的信号复包络；表示在第k个短时隙内的阵列加性噪声；表示在第k个短时隙内的多普勒频移，其表达式为其中p
k
表示运动阵列在第k个短时隙内的位置向量；表示运动阵列在第k个短时隙内的速度向量；c表示信号传播速度；f
o
表示信号中心频率；为了得到被动合成孔径阵列输出，在第k个短时隙内对每个阵元共产生m个延时输出，延时量分别为τ1,τ2,
…
,τ
m
，其中第m个延时的阵列输出表示为式中表示在第k个短时隙内第m个延时的阵列加性噪声；
最后将m个延时的阵列输出进行合并，得到第k个短时隙内的被动合成孔径阵列输出，其表达式为式中表示在第k个短时隙内的被动合成孔径阵列加性噪声；表示在第k个短时隙内的被动合成孔径阵列流形向量，其中3.根据权利要求1所述的一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，其特征在于，所述步骤4中，针对二维网格点构造第k个短时隙的被动合成孔径阵列输出协方差矩阵如下式所示式中t
l
表示第l个采样时间；l表示信号采样数量；表示第k个短时隙内第l个采样时间的被动合成孔径阵列输出；对矩阵进行特征值分解，其最大特征值对应的首1特征向量为于是形成如下首1特征矩阵样本集4.根据权利要求1所述的一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，其特征在于，所述步骤6中，针对待定位窄带信号源，将第k个短时隙内的运动阵列输出进行数字下变频处理，得到如下输出式中表示在第k个短时隙内针对待定位信号源的阵列流形向量；表示待定位信号源的位置向量；表示在第k个短时隙内待定位信号复包络；表示在第k个短时隙内针对待定位信号源的阵列加性噪声；表示在第k个短时隙内针对待定位信号源的多普勒频移，其表达式为其中p
k
表示运动阵列在第k个短时隙内的位置向量；表示运动阵列在第k个短时隙内的速度向量；c表示信号传播速度；f
o
表示信号中心频率；为了得到被动合成孔径阵列输出，在第k个短时隙内对每个阵元共产生m个延时输出，
延时量分别为τ1,τ2,
…
,τ
m
，其中第m个延时的阵列输出可以表示为式中表示针对待定位信号源的第m个延时的阵列加性噪声；最后将m个延时的阵列输出进行合并，得到第k个短时隙内的针对待定位信号源的被动合成孔径阵列输出，其表达式为式中表示第k个短时隙内的针对待定位信号源的被动合成孔径阵列加性噪声；表示第k个短时隙内的针对待定位信号源的被动合成孔径阵列流形向量，其中5.根据权利要求1所述的一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，其特征在于，所述步骤7中，针对待定位窄带信号源，构造第k个短时隙的被动合成孔径阵列输出协方差矩阵如下式所示其中t
l
表示第l个采样时间；l表示信号采样数量；表示第k个短时隙内第l个采样时间针对待定位信号源的被动合成孔径阵列输出；对矩阵进行特征值分解，其最大特征值对应的首1特征向量为于是形成如下首1特征矩阵6.根据权利要求1所述的一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，其特征在于，所述步骤8中，将步骤7中的首1特征矩阵与步骤4中的首1特征矩阵样本集中的每个矩阵做复相关运算，得到相关值其表达式为7.根据权利要求1所述的一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，其特征在于，所述步骤9中，令相关值最大值对应的序号n1和n2分别为n
x
和n
y
，取，取9个点进行二次曲面拟合，并按照下式构造矩阵z和向量ρ
式中于是最小二乘二次曲面拟合参数为θ
ls
＝(z
t
z)-1
z
t
ρ。8.根据权利要求7所述的一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，其特征在于，所述步骤10中，基于最小二乘二次曲面拟合参数给出最终定位结果其表达式为式中<θ
ls
>2、<θ
ls
>3、<θ
ls
>4、<θ
ls
>5、<θ
ls
>6分别表示最小二乘二次曲面拟合参数θ
ls
中的第2、3、4、5、6个元素。
技术总结本发明公开了一种基于单运动阵列被动合成孔径的相关干涉仪直接定位方法，包括：在定位区域内划分两维网格点，在各点放置信号源，并利用具有相同航迹和速度的运动阵列接收各信号；针对各信号，运动阵列利用多个时隙构造被动合成孔径阵列输出及对应协方差矩阵，通过特征值分解得到首1特征矩阵样本集；针对待定位窄带信号源，利用相同航迹和速度的运动阵列以被动合成孔径的方式接收各信号，并计算其阵列输出协方差矩阵，得到首1特征矩阵；将所述首1特征矩阵与矩阵样本集中的每个矩阵进行复相关运算，取相关值最大9点进行曲面拟合，形成二次曲面，进而得到最终定位结果。本发明提高了在小阵列孔径条件下的定位精度，并极大降低实时定位的计算量。时定位的计算量。时定位的计算量。
技术研发人员：王鼎 高卫港 尹洁昕 张莉 杨宾 郑娜娥 王建辉 胡德秀 唐涛 任修坤
受保护的技术使用者：中国人民解放军战略支援部队信息工程大学
技术研发日：2022.06.25
技术公布日：2022/11/1
