CN112333629A - 一种互耦未知条件下分布式阵列多目标定位方法 - Google Patents
一种互耦未知条件下分布式阵列多目标定位方法 Download PDFInfo
- Publication number
- CN112333629A CN112333629A CN202011143540.7A CN202011143540A CN112333629A CN 112333629 A CN112333629 A CN 112333629A CN 202011143540 A CN202011143540 A CN 202011143540A CN 112333629 A CN112333629 A CN 112333629A
- Authority
- CN
- China
- Prior art keywords
- mutual coupling
- signal
- observation station
- observation
- array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000008878 coupling Effects 0.000 title claims abstract description 56
- 238000010168 coupling process Methods 0.000 title claims abstract description 56
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 27
- 238000005457 optimization Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 13
- 230000003595 spectral effect Effects 0.000 claims description 6
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 238000013499 data model Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000006880 cross-coupling reaction Methods 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000035899 viability Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/023—Services making use of location information using mutual or relative location information between multiple location based services [LBS] targets or of distance thresholds
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/025—Services making use of location information using location based information parameters
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种互耦未知条件下分布式阵列多目标定位方法,每个观测站均采集得到阵列信号时域数据,求每个观测站接收数据的协方差矩阵,并计算出噪声子空间,得到仅关于目标位置参数的数学优化模型,多个观测站联合解算,获得信号源的初步直接定位结果;使用解算出来的目标位置信息,求得未知互耦系数,将求得的互耦系数补偿回数学优化模型,求得修正后的定位结果。本发明能够在阵列存在未知互耦、来波信号在信噪比较低的情况下,确保直接定位准确率与定位精度,在无需参数测量的情况下,直接解算出目标信号源的位置,计算起来更加高效。
Description
技术领域
本发明涉及无线电信号定位领域,适用于在阵列存在互耦影响的情况下,利用无线电信号数据域信息直接完成目标位置求解的直接位置确定方法。
背景技术
无线信号定位技术广泛应用于通信、导航、目标检测、雷达、安全管理等领域。在工业生产和军事应用中都发挥着重要作用。但是随着信息对抗技术的不断发展,无线电定位系统面临着巨大的威胁和挑战。无源定位技术因其不主动发射电磁信号、生存能力强、侦察作用距离远等优势,是目前研究较为广泛的定位技术。
目标直接定位方式(如图1所示),可以直接从原始信号数据域提取目标未知参数,有利于直接利用信号的波形特征,解算出目标的位置信息。并且可以有效利用各观测站接收数据间的相关性,从而进一步提高定位的准确率与精度。
当前无线电信号直接定位技术能够实现在不考虑阵列存在未知互耦影响的直接定位,然而,在天线阵列中常常因为单元间的互耦使得天线阵列的特性发生根本性的变化,从而降低直接定位的准确率与定位精度。
对阵列单元之间的互耦影响进行分析,并对其互耦做出一定的补偿,可以有效提高直接定位算法中的参数估计精度。
发明内容
为了克服现有技术的不足,本发明提供一种互耦未知条件下分布式阵列多目标定位方法,本发明提供一种在阵列存在互耦影响时的无线电信号直接定位技术,能够在阵列存在未知互耦、来波信号在信噪比较低的情况下,确保直接定位准确率与定位精度。
本发明解决其技术问题所采用的技术方案包括以下步骤:
(1)设置初始化参数;
(2)每个观测站均采集目标在采样快拍数K个时隙内辐射的无线电信号数据,并按照预设的采样信号时域模型得到目标的K个时隙的阵列信号时域数据;
(3)求每个观测站接收数据的协方差矩阵,并对协方差矩阵进行特征分解,并计算出噪声子空间;
(4)使用初始化的参数及求得的噪声子空间,得到仅关于目标位置参数的数学优化模型,多个观测站联合解算,获得信号源的初步直接定位结果;
(5)使用解算出来的目标位置信息,求得未知互耦系数;
(6)将求得的互耦系数补偿回数学优化模型,求得修正后的定位结果。
所述步骤(2)中,假设有Q个静止的信号源,L个静止的观测站,每个观测站天线阵列阵元数目均为M,且Q<M,采样的快拍数是K,第l个观测站阵列的有ml个非零互耦系数,l=1,2,…,L,且ml<M,记第q个信号源的位置坐标为q=1,2,…,Q,记第l个观测站的位置坐标为记第l个观测站的互耦系数向量为则在t时刻的观测数据模型为:
X(t)=Clal(p)s(t)+n(t) (1)
式中s(t)为来波信号,Cl为第l个观测站互耦系数的Toeplitz矩阵,al(p)为第l个阵列的导向矢量,p=[x,y]为观测站观测范围内任意一点,n(t)为均值为零的加性高斯白噪声;
al(p)=[1,β(p),…,β(p)M-1]T,β(p)=e-j2πd sin(p)/λ (3)
每个观测站均采集目标在K个时隙内辐射的无线电信号数据,并按照预设的采样信号时域模型得到目标的K个时隙的阵列信号时域数据;第l个运动观测站在第k个时隙内所接收到的采样信号时域模型Xl(k)为:
Xl(k)=Clal(p)s(k)+n(k) (4)
其中,k=1,2,…,K;
所述步骤(3)中,求第l个观测站接收数据的协方差矩阵Rl,并对其进行分解,得到:
其中,Ul,s为第l个观测站接收信号的信号子空间,∑l,s和∑l,v分别为第l个观测站接收信号的信号和噪声的特征值组成的对角阵,Ul,v为第l个观测站接收信号的噪声子空间,(·)H表示求矩阵的共轭转置。
所述步骤(4)中得到初步直接定位结果的步骤为:
根据采样信号时域模型Xl(k)和互耦系数矩阵,第l个观测站的导向矢量写为:
其中,
且αl=[μl,1,…,μl,l-1,1,αl,1,…,αl,l-1]T,αl是一个(2Pl-1)×1的非零向量,
由于接收信号的导向矢量和噪声子空间是正交的,得到(8)式:
其中Ul,v为第l个观测站的噪声子空间;
接下来通过对式(11)在观测站观测范围内进行搜索得到最大值所对应的坐标,即为定位的初步结果;
其中,det[·]表示矩阵的行列式。
所述步骤(5)中,求得估计互耦系数的步骤为:
得到估计互耦系数。
所述步骤(6)中,求得修正定位结果的步骤为:
估计出互耦系数后,将第l个观测站互耦系数cl的Toeplitz矩阵Cl,带入式(16):
在观测站观测范围内对式(16)的谱峰进行搜索后,谱峰出现的对应坐标即为直接定位结果。
本发明的有益效果在于:
1.使用直接定位的方法,在无需参数测量的情况下,直接解算出目标信号源的位置,计算起来更加高效。
2.考虑了阵列单元之间的未知互耦对直接定位算法的影响,并对其进行补偿,得到了更高的定位精度。
附图说明
图1为本发明所用直接定位技术的原理简图。
图2为本发明所用直接定位技术的流程简图。
图3为本发明假设信号源和观测站的位置分布图。
图4为本发明估计的互耦系数随信噪比变化的均方根误差(RMSE)曲线。
图5为本发明定位结果随信噪比变化的RMSE曲线。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本发明解决技术问题的步骤为:
(1)设置初始化参数。
(2)每个观测站均采集目标在采样快拍数K个时隙内辐射的无线电信号数据,并按照预设的采样信号时域模型得到目标的K个时隙的阵列信号时域数据。
(3)求每个观测站接收数据的协方差矩阵,并对协方差矩阵进行特征分解,并计算出噪声子空间。
(4)使用初始化的参数及求得的噪声子空间,经过一系列的数学推演得到仅关于目标位置参数的数学优化模型,多个观测站联合解算,获得信号源的大致位置。
(5)使用解算出来的目标位置信息,经过数学推导之后,求得未知互耦系数。
(6)将求得的互耦系数补偿回数学优化模型,求得更加精确的信号源位置。
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
如图2所示,本发明的技术方案步骤为:
步骤一:设置初始化参数;
假设有Q个静止的信号源,L个静止的观测站,如图3所示,每个观测站天线阵列阵元数目均为M,且Q<M,采样的快拍数是K,第l个观测站阵列的有ml个非零互耦系数,l=1,2,…,L,且ml<M,记第q个信号源的位置坐标为q=1,2,…,Q,记第l个观测站的位置坐标为记第l个观测站的互耦系数向量为则在t时刻的观测数据模型为:
X(t)=Clal(p)s(t)+n(t) (1)
式中s(t)为来波信号,Cl为第l个观测站互耦系数的Toeplitz矩阵,al(p)为第l个阵列的导向矢量,p=[x,y]为观测站观测范围内任意一点,n(t)为均值为零的加性高斯白噪声;
al(p)=[1,β(p),…,β(p)M-1]T,β(p)=e-j2πd sin(p)/λ (3)
每个观测站均采集目标在K个时隙内辐射的无线电信号数据,并按照预设的采样信号时域模型得到目标的K个时隙的阵列信号时域数据;第l个运动观测站在第k个时隙内所接收到的采样信号时域模型Xl(k)为:
Xl(k)=Clal(p)s(k)+n(k) (4)
其中,k=1,2,…,K;
步骤二:求第l个观测站接收数据的协方差矩阵Rl,并对其进行分解,得到:
其中,Ul,s为第l个观测站接收信号的信号子空间,∑l,s和∑l,v分别为第l个观测站接收信号的信号和噪声的特征值组成的对角阵,Ul,v为第l个观测站接收信号的噪声子空间,(·)H表示求矩阵的共轭转置;
步骤三;优化模型,求出初步直接定位结果;
根据采样信号时域模型Xl(k)和互耦系数矩阵,第l个观测站的导向矢量写为:
其中,
且αl=[μl,1,…,μl,l-1,1,αl,1,…,αl,l-1]T,αl是一个(2Pl-1)×1的非零向量,
由于接收信号的导向矢量和噪声子空间是正交的,得到(8)式:
其中Ul,v为第l个观测站的噪声子空间;
接下来通过对式(11)在观测站观测范围内进行搜索得到最大值所对应的坐标,即为定位的初步结果,det[·]表示矩阵的行列式;
步骤四:估计互耦系数;
步骤五:补偿互耦系数,修正定位结果;
估计出互耦系数后,将第l个观测站互耦系数cl的Toeplitz矩阵Cl,带入式(16):
在观测站观测范围内对式(16)的谱峰进行搜索后,谱峰出现的对应坐标即为直接定位结果。
本发明的效果通过以下仿真结果进一步说明。
观测站的个数为L=4,每个观测站为阵元个数为M=10的均匀线阵,观测站的位置坐标分别为(-100,100),(-100,-100),(100,100)和(100,-100),四个观测站的互耦系数分别为c1=[1,-0.1545+0.4755i],c2=[1,0.14142-0.14142i],c3=[1,0.43301-0.25i],c4=[1,-0.1880-0.0582i];信号源的个数为Q=2,信号源的位置坐标为(10,20),(-10,50);采样快拍数K=200。背景观测噪声为零均值的高斯白噪声且信噪比从-15dB变化到10dB。图4为估计的互耦系数随信噪比变化的均方根误差(RMSE)曲线,图5为定位结果随信噪比变化的RMSE曲线。
Claims (6)
1.一种互耦未知条件下分布式阵列多目标定位方法,其特征在于包括下述步骤:
(1)设置初始化参数;
(2)每个观测站均采集目标在采样快拍数K个时隙内辐射的无线电信号数据,并按照预设的采样信号时域模型得到目标的K个时隙的阵列信号时域数据;
(3)求每个观测站接收数据的协方差矩阵,并对协方差矩阵进行特征分解,并计算出噪声子空间;
(4)使用初始化的参数及求得的噪声子空间,得到仅关于目标位置参数的数学优化模型,多个观测站联合解算,获得信号源的初步直接定位结果;
(5)使用解算出来的目标位置信息,求得未知互耦系数;
(6)将求得的互耦系数补偿回数学优化模型,求得修正后的定位结果。
2.根据权利要求1所述的一种互耦未知条件下分布式阵列多目标定位方法,其特征在于:
所述步骤(2)中,假设有Q个静止的信号源,L个静止的观测站,每个观测站天线阵列阵元数目均为M,且Q<M,采样的快拍数是K,第个观测站阵列的有个非零互耦系数,且记第q个信号源的位置坐标为q=1,2,…,Q,记第个观测站的位置坐标为记第个观测站的互耦系数向量为则在t时刻的观测数据模型为:
al(p)=[1,β(p),…,β(p)M-1]T,β(p)=e-j2πd sin(p)/λ (3)
其中,k=1,2,…,K。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011143540.7A CN112333629A (zh) | 2020-10-23 | 2020-10-23 | 一种互耦未知条件下分布式阵列多目标定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011143540.7A CN112333629A (zh) | 2020-10-23 | 2020-10-23 | 一种互耦未知条件下分布式阵列多目标定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112333629A true CN112333629A (zh) | 2021-02-05 |
Family
ID=74312342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011143540.7A Pending CN112333629A (zh) | 2020-10-23 | 2020-10-23 | 一种互耦未知条件下分布式阵列多目标定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112333629A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114417917A (zh) * | 2021-12-31 | 2022-04-29 | 西北工业大学 | 一种未知互耦条件下直接定位方法 |
CN117202343A (zh) * | 2023-08-30 | 2023-12-08 | 中国人民解放军战略支援部队信息工程大学 | 一种针对多个宽带信号辐射源的分布式阵列协同直接定位方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2824479A1 (en) * | 2013-07-12 | 2015-01-14 | Ingenieria de Sistemas Para La Defensa de Espana, S.A. | Method and system for locating a moving vehicle |
CN107132505A (zh) * | 2017-05-19 | 2017-09-05 | 中国人民解放军信息工程大学 | 直达与非直达混合场景中的多目标直接定位方法 |
CN108919231A (zh) * | 2018-05-25 | 2018-11-30 | 长江大学 | 联合误差条件下双基地mimo雷达角度估计方法 |
-
2020
- 2020-10-23 CN CN202011143540.7A patent/CN112333629A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2824479A1 (en) * | 2013-07-12 | 2015-01-14 | Ingenieria de Sistemas Para La Defensa de Espana, S.A. | Method and system for locating a moving vehicle |
CN107132505A (zh) * | 2017-05-19 | 2017-09-05 | 中国人民解放军信息工程大学 | 直达与非直达混合场景中的多目标直接定位方法 |
CN108919231A (zh) * | 2018-05-25 | 2018-11-30 | 长江大学 | 联合误差条件下双基地mimo雷达角度估计方法 |
Non-Patent Citations (5)
Title |
---|
XIANPENG WANG等: "Robust Vehicle Localization Exploiting Two Based Stations Cooperation: A MIMO Radar Perspective", 《IEEE ACCESS》, 27 August 2018 (2018-08-27) * |
李新潮等: "阵列互耦条件下相干信源DOA估计算法研究", 《舰船电子对抗》, no. 06, 25 December 2010 (2010-12-25) * |
潘哲等: "未知互耦条件下混合信号波达方向估计", 《海军航空工程学院学报》, no. 02, 20 March 2015 (2015-03-20) * |
王伟等: "双基地MIMO雷达多目标定位及互耦参数估计", 《华中科技大学学报(自然科学版)》, no. 07, 23 July 2012 (2012-07-23) * |
郭艺夺等: "双基地MIMO雷达目标定位及互耦自校正算法", 《西安电子科技大学学报》, no. 05, 23 May 2012 (2012-05-23) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114417917A (zh) * | 2021-12-31 | 2022-04-29 | 西北工业大学 | 一种未知互耦条件下直接定位方法 |
CN117202343A (zh) * | 2023-08-30 | 2023-12-08 | 中国人民解放军战略支援部队信息工程大学 | 一种针对多个宽带信号辐射源的分布式阵列协同直接定位方法 |
CN117202343B (zh) * | 2023-08-30 | 2024-07-19 | 中国人民解放军战略支援部队信息工程大学 | 一种针对多个宽带信号辐射源的分布式阵列协同直接定位方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107132505B (zh) | 直达与非直达混合场景中的多目标直接定位方法 | |
Davis et al. | Angle estimation with adaptive arrays in external noise fields | |
CN108957387B (zh) | 一种卫星信号二维到达角估计方法及系统 | |
US11982698B2 (en) | Joint denoising and delay estimation for the extraction of pulse-width of signals in RF interference | |
CN110515037B (zh) | 非视距环境下能时频多域联合的被动定位方法 | |
CN111239677B (zh) | 基于数字阵列的多波束被动单脉冲测角方法 | |
CN111199280B (zh) | 短波信道模型误差存在下联合信号复包络和载波相位信息的多站目标源地理坐标估计方法 | |
CN108710103A (zh) | 基于稀疏阵列的强弱多目标超分辨测向与信源数估计方法 | |
CN101644760B (zh) | 一种适用于高分辨阵列的快速鲁棒的信源个数检测方法 | |
US11681006B2 (en) | Method for jointly estimating gain-phase error and direction of arrival (DOA) based on unmanned aerial vehicle (UAV) array | |
CN110426670A (zh) | 基于tls-cs的外辐射源雷达超分辨doa估计方法 | |
EP4334741A1 (en) | Method for two-dimensional and three-dimensional imaging based on collocated multiple-input multiple-output radars | |
CN113311402B (zh) | 考虑发射波形非理想正交的mimo雷达目标测向方法及系统 | |
CN108594165B (zh) | 一种基于期望最大化算法的窄带信号波达方向估计方法 | |
CN112333629A (zh) | 一种互耦未知条件下分布式阵列多目标定位方法 | |
CN115826004B (zh) | 一种基于二维角度及时差联合的三星协同直接定位方法 | |
CN111007457A (zh) | 一种基于块稀疏贝叶斯模型的辐射源直接定位方法 | |
CN114488217B (zh) | 基于深度学习的高轨卫星cei信号频率估计方法 | |
CN116299387A (zh) | 非均匀杂波下干扰正交抑制的目标智能检测方法 | |
CN115052246A (zh) | 一种未知衰减系数下基于多频率代价函数融合的宽带信号直接定位方法 | |
CN112881971B (zh) | 电磁方向性互耦效应下相干干扰源测向方法 | |
CN115575906A (zh) | 非均匀背景下干扰智能抑制的融合检测方法 | |
CN115524672A (zh) | 结构化干扰与杂波下目标鲁棒智能检测方法 | |
CN114969642A (zh) | 一种均匀线阵下的小扩展角度分布源通用信源数估计方法 | |
Ren et al. | A Fast Direct Position Determination Algorithm for LFM Signal Based on Spectrum Detection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210205 |
|
WD01 | Invention patent application deemed withdrawn after publication |