CN116973834A - 一种强干扰下基于协同波束形成置信度的测向方法 - Google Patents

一种强干扰下基于协同波束形成置信度的测向方法 Download PDF

Info

Publication number
CN116973834A
CN116973834A CN202310955692.4A CN202310955692A CN116973834A CN 116973834 A CN116973834 A CN 116973834A CN 202310955692 A CN202310955692 A CN 202310955692A CN 116973834 A CN116973834 A CN 116973834A
Authority
CN
China
Prior art keywords
beam forming
confidence
finding
determining
directions
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.)
Granted
Application number
CN202310955692.4A
Other languages
English (en)
Other versions
CN116973834B (zh
Inventor
张兴虎
庄洪林
万群
刘宁
汤四龙
刘彦红
于波
宋强
常凯
刘翼鹏
彭翔宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Tongfang Industrial Co Ltd
Institute of Systems Engineering of PLA Academy of Military Sciences
Original Assignee
University of Electronic Science and Technology of China
Tongfang Industrial Co Ltd
Institute of Systems Engineering of PLA Academy of Military Sciences
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China, Tongfang Industrial Co Ltd, Institute of Systems Engineering of PLA Academy of Military Sciences filed Critical University of Electronic Science and Technology of China
Priority to CN202310955692.4A priority Critical patent/CN116973834B/zh
Publication of CN116973834A publication Critical patent/CN116973834A/zh
Application granted granted Critical
Publication of CN116973834B publication Critical patent/CN116973834B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明涉及一种强干扰下基于协同波束形成置信度的测向方法,属于无线电测向技术领域。本发明针对测向设备因附近存在强干扰导致对弱信号难以可靠测向的问题,对天线阵列不同时刻、不同方向的快拍信号进行协同波束形成,利用协同波束形成置信度,实现强干扰下对弱信号测向的目的。本发明可用于测向设备与附近的用频设备工作于相同的时域、频域的场合,可为无线电干扰源监测提供更宽频段、更大范围、更连续的测向信息。

Description

一种强干扰下基于协同波束形成置信度的测向方法
技术领域
本发明属于无线电测向技术领域,涉及一种强干扰下基于协同波束形成置信度的测向方法,具体涉及一种对天线阵列不同时刻、不同方向的快拍信号进行协同波束形成,利用协同波束形成置信度进行测向的方法。
背景技术
无线电信号的来波方向是无线电频谱信息中的重要参数,对合法的频谱资源使用者、无意的无线电干扰源、非法的频谱资源占用者以及战场环境中的无线电装备的监测与定位具有重要意义。目前,已有很多估计无线电信号来波方向的方法,特别是常规波束形成、卡彭(Capon)波束形成等基于方位功率轮廓的空间谱测向方法,在测向设备上得到了广泛的应用。
测向站的测向设备在开放的无线电信号环境中工作时,各种无线电信号在同一时域、频域混叠,使得距离测向站较远的用频设备的来波信号淹没在距离测向站较近的用频设备的来波信号和噪声中,导致基于方位功率轮廓的空间谱测向方法对距离测向站较远的用频设备的来波信号测向性能恶化、测向功能失效。因此,为了提高测向站适应复杂无线电信号环境的能力、对弱信号的测向能力、扩大测向站有效覆盖的地理范围,在无线电测向技术领域出现了引入方位功率轮廓之外的空间谱特征进行测向的一些方法。
多信号分类(MUSIC)方法是一种典型的伪空间谱测向方法,就不是直接利用方位功率轮廓信息进行测向,而是利用阵列快拍数据中信号来波方向向量与阵列快拍数据所在信号空间对应的正交子空间之间的正交关系进行来波方向估计。但是,这类伪空间谱测向方法的测向性能依赖利用阵列快拍数据估计信号空间、正交子空间的精度,适用于高信噪比的场合,对信噪比较低、附近方向存在强干扰的信号而言,依然存在难以可靠测向的问题。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是如何提供一种强干扰下基于协同波束形成置信度的测向方法,以解决测向设备因附近存在强干扰导致对弱信号难以可靠测向的问题。
(二)技术方案
为了解决上述技术问题,本发明提出一种强干扰下基于协同波束形成置信度的测向方法,该方法包括如下步骤:
S1、首先设置测向设备的天线个数、搜索的方向个数、搜索的方向集合、与搜索的方向集合对应的方向向量集合,测向设备的所有天线每次采集的阵列快拍信号的快拍数;
S2、其次确定测向设备的所有天线在时刻1采集的阵列快拍信号组成的矩阵、在与时刻1不同的时刻2采集的阵列快拍信号组成的矩阵;
S3、然后对应搜索的方向集合中的方向,确定时刻1的波束形成器的输入和权向量、时刻2的波束形成器的输入和权向量;
S4、接着对应搜索的方向集合中的方向,确定时刻1的波束形成器的输出、时刻2的波束形成器的输出,进而确定协同波束形成置信度;
S5、最后,由协同波束形成置信度集合中的最大值对应的搜索方向确定强干扰下基于协同波束形成置信度的测向结果。
(三)有益效果
本发明提出一种强干扰下基于协同波束形成置信度的测向方法,本发明的有益效果是:使用本发明提出的一种强干扰下基于协同波束形成置信度的测向方法,可在时刻1的阵列快拍信号中同时、同频混叠了来自一个来波方向的强干扰信号和另一个来波方向的弱信号,并且弱信号功率小于噪声功率,时刻2的阵列快拍信号中存在与弱信号波形相同的信号,这些信号的来波方向都未知的情况下,对时刻1的阵列快拍信号和时刻2的阵列快拍信号进行协同波束形成,由协同波束形成置信度即可确定弱信号的来波方向。因此,本发明可用于测向设备与附近的用频设备工作于相同的时域、频域的场合,可为无线电干扰源监测提供更宽频段、更大范围、更连续的测向信息。
具体实施方式
为使本发明的目的、内容和优点更加清楚,下面结合实施例,对本发明的具体实施方式作进一步详细描述。
本发明针对测向设备因附近存在强干扰导致对弱信号难以可靠测向的问题,对天线阵列不同时刻、不同方向的快拍信号进行协同波束形成,利用协同波束形成置信度,实现强干扰下对弱信号测向的目的。
本发明的技术方案为:
一种强干扰下基于协同波束形成置信度的测向方法,该方法包括如下步骤:
S1、首先设置测向设备的天线个数、搜索的方向个数、搜索的方向集合、与搜索的方向集合对应的方向向量集合,测向设备的所有天线每次采集的阵列快拍信号的快拍数;
S2、其次确定测向设备的所有天线在时刻1采集的阵列快拍信号组成的矩阵、在与时刻1不同的时刻2采集的阵列快拍信号组成的矩阵;
S3、然后对应搜索的方向集合中的方向,确定时刻1的波束形成器的输入和权向量、时刻2的波束形成器的输入和权向量;
S4、接着对应搜索的方向集合中的方向,确定时刻1的波束形成器的输出、时刻2的波束形成器的输出,进而确定协同波束形成置信度;
S5、最后,由协同波束形成置信度集合中的最大值对应的搜索方向确定强干扰下基于协同波束形成置信度的测向结果。
本发明具体包括以下步骤:
S1、设置测向设备的天线个数为M、搜索的方向个数为N、搜索方向为θn,搜索的方向集合为{θ1,θ2,...,θN}、与搜索方向θn对应的方向向量为a(θn),n=1,2,...,N,所有方向向量组成的集合为{a(θ1),a(θ2),......,a(θN)},测向设备的所有天线每次采集的阵列快拍信号的快拍数为L;
S2、确定测向设备的所有天线在时刻1采集的阵列快拍信号,组成M×L阶矩阵X;在与时刻1不同的时刻2采集的阵列快拍信号,组成M×L阶矩阵Y;
S3、对应搜索的方向集合中的搜索方向θn,确定时刻1的波束形成器的输入为X,权向量为
其中,w0n)=(XXH)-1a(θn),H表示向量或矩阵的共轭转置,-1表示矩阵的逆,n=1,2,...,N;
对应搜索的方向集合中的搜索方向θn,确定时刻2的波束形成器的输入为Y,权向量为
w2n)=(XYH)-1w1n)
S4、对应搜索的方向集合中的搜索方向θn,确定时刻1的波束形成器的输出为
时刻2的波束形成器的输出为
对应搜索的方向集合中的搜索方向θn,确定协同波束形成置信度为
其中,|| ||为向量的范数,n=1,2,...,N;
S5、确定协同波束形成置信度的集合{g(θ1),g(θ2),...,g(θN)}中的最大值对应的搜索方向为强干扰下基于协同波束形成置信度的测向结果。
实施例1:
在本例中,设置测向设备天线个数M=8、搜索的方向个数N=181和搜索的方向集合{-90,-89,...,89,90}度、与搜索的方向集合一一对应的方向向量组成的集合为{a(-90),a(-89),......,a(89),a(90)},采集阵列快拍信号的快拍数L=32,64,...,4096,本例中L=4096。在测向设备与附近用频设备同时、同频工作的情况下,相对于测向设备,在测向设备的所有天线在时刻1采集的阵列快拍信号中,附近用频设备的信号来波方向为30.76度,信噪比为30dB,弱信号的来波方向为20.68度,信噪比为-5dB;在测向设备的所有天线在时刻2采集的阵列快拍信号中,与弱信号具有相同波形的信号来波方向为10.29度,信噪比为10dB。
利用常规波束形成、卡彭(Capon)波束形成等测向方法,只能准确的确定测向设备附近的用频设备的信号来波方向,弱信号来波方向的测向均方根误差较大,都超过9度;而采用本发明方法,弱信号来波方向的测向均方根误差为0.32度,实现了测向设备与附近的用频设备工作于相同的时域、频域的情况下对弱信号测向的目的。
本发明的有益效果是:使用本发明提出的一种强干扰下基于协同波束形成置信度的测向方法,可在时刻1的阵列快拍信号中同时、同频混叠了来自一个来波方向的强干扰信号和另一个来波方向的弱信号,并且弱信号功率小于噪声功率,时刻2的阵列快拍信号中存在与弱信号波形相同的信号,这些信号的来波方向都未知的情况下,对时刻1的阵列快拍信号和时刻2的阵列快拍信号进行协同波束形成,由协同波束形成置信度即可确定弱信号的来波方向。因此,本发明可用于测向设备与附近的用频设备工作于相同的时域、频域的场合,可为无线电干扰源监测提供更宽频段、更大范围、更连续的测向信息。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。

Claims (10)

1.一种强干扰下基于协同波束形成置信度的测向方法,其特征在于,该方法包括如下步骤:
S1、首先设置测向设备的天线个数、搜索的方向个数、搜索的方向集合、与搜索的方向集合对应的方向向量集合,测向设备的所有天线每次采集的阵列快拍信号的快拍数;
S2、其次确定测向设备的所有天线在时刻1采集的阵列快拍信号组成的矩阵、在与时刻1不同的时刻2采集的阵列快拍信号组成的矩阵;
S3、然后对应搜索的方向集合中的方向,确定时刻1的波束形成器的输入和权向量、时刻2的波束形成器的输入和权向量;
S4、接着对应搜索的方向集合中的方向,确定时刻1的波束形成器的输出、时刻2的波束形成器的输出,进而确定协同波束形成置信度;
S5、最后,由协同波束形成置信度集合中的最大值对应的搜索方向确定强干扰下基于协同波束形成置信度的测向结果。
2.如权利要求1所述的强干扰下基于协同波束形成置信度的测向方法,其特征在于,所述步骤S1具体包括:设置测向设备的天线个数为M、搜索的方向个数为N、搜索方向为θn,搜索的方向集合为{θ12,…,θN}、与搜索方向θn对应的方向向量为a(θn),n=1,2,…,N,所有方向向量组成的集合为{α(θ1),a(θ2),……,a(θN)},测向设备的所有天线每次采集的阵列快拍信号的快拍数为L。
3.如权利要求2所述的强干扰下基于协同波束形成置信度的测向方法,其特征在于,所述步骤S1中,搜索的方向个数N=181和搜索的方向集合{-90,-89,…,89,90}度、与搜索的方向集合一一对应的方向向量组成的集合为{a(-90),a(-89),……,a(89),a(90)}。
4.如权利要求2所述的强干扰下基于协同波束形成置信度的测向方法,其特征在于,所述步骤S1中,采集阵列快拍信号的快拍数L=32,64,…,4096。
5.如权利要求2所述的强干扰下基于协同波束形成置信度的测向方法,其特征在于,所述步骤S2具体包括:确定测向设备的所有天线在时刻1采集的阵列快拍信号,组成M×L阶矩阵X;在与时刻1不同的时刻2采集的阵列快拍信号,组成M×L阶矩阵Y。
6.如权利要求5所述的强干扰下基于协同波束形成置信度的测向方法,其特征在于,所述步骤S3中,确定时刻1的波束形成器的输入和权向量具体包括:
对应搜索的方向集合中的搜索方向θn,确定时刻1的波束形成器的输入为X,权向量为
其中,w0n)=(XXH)-1a(θn),H表示向量或矩阵的共轭转置,-1表示矩阵的逆,n=1,2,…,N。
7.如权利要求6所述的强干扰下基于协同波束形成置信度的测向方法,其特征在于,所述步骤S3中,确定时刻2的波束形成器的输入和权向量具体包括:
对应搜索的方向集合中的搜索方向θn,确定时刻2的波束形成器的输入为Y,权向量为
w2n)=(XYH)-1w1n)。
8.如权利要求7所述的强干扰下基于协同波束形成置信度的测向方法,其特征在于,所述步骤S4中,对应搜索的方向集合中的方向,确定时刻1的波束形成器的输出、时刻2的波束形成器的输出具体包括:
对应搜索的方向集合中的搜索方向θn,确定时刻1的波束形成器的输出为
时刻2的波束形成器的输出为
9.如权利要求8所述的强干扰下基于协同波束形成置信度的测向方法,其特征在于,所述步骤S4中,确定协同波束形成置信度具体包括:
对应搜索的方向集合中的搜索方向θn,确定协同波束形成置信度为
其中,|| ||为向量的范数,n=1,2,…,N。
10.如权利要求9所述的强干扰下基于协同波束形成置信度的测向方法,其特征在于,所述步骤S5具体包括:确定协同波束形成置信度的集合{g(θ1),g(θ2),…,g(θN)}中的最大值对应的搜索方向为强干扰下基于协同波束形成置信度的测向结果。
CN202310955692.4A 2023-08-01 2023-08-01 一种强干扰下基于协同波束形成置信度的测向方法 Active CN116973834B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310955692.4A CN116973834B (zh) 2023-08-01 2023-08-01 一种强干扰下基于协同波束形成置信度的测向方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310955692.4A CN116973834B (zh) 2023-08-01 2023-08-01 一种强干扰下基于协同波束形成置信度的测向方法

Publications (2)

Publication Number Publication Date
CN116973834A true CN116973834A (zh) 2023-10-31
CN116973834B CN116973834B (zh) 2024-01-30

Family

ID=88481118

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310955692.4A Active CN116973834B (zh) 2023-08-01 2023-08-01 一种强干扰下基于协同波束形成置信度的测向方法

Country Status (1)

Country Link
CN (1) CN116973834B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030231547A1 (en) * 2002-06-18 2003-12-18 Tsih Yang Method and apparatus for passive acoustic imaging using a horizontal line array
CN104980201A (zh) * 2014-04-07 2015-10-14 想象技术有限公司 波束成形矩阵的重排序
CN107045117A (zh) * 2016-09-23 2017-08-15 南京工程学院 基于Capon波束形成定位方法及装置
CN110161489A (zh) * 2019-05-21 2019-08-23 西安电子科技大学 一种基于伪框架的强弱信号测向方法
CN112949100A (zh) * 2020-11-06 2021-06-11 中国人民解放军空军工程大学 一种机载雷达抗主瓣干扰的方法
CN113466899A (zh) * 2021-08-13 2021-10-01 电子科技大学 高信噪比环境下基于小快拍数的导航接收机波束形成方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030231547A1 (en) * 2002-06-18 2003-12-18 Tsih Yang Method and apparatus for passive acoustic imaging using a horizontal line array
CN104980201A (zh) * 2014-04-07 2015-10-14 想象技术有限公司 波束成形矩阵的重排序
CN107045117A (zh) * 2016-09-23 2017-08-15 南京工程学院 基于Capon波束形成定位方法及装置
CN110161489A (zh) * 2019-05-21 2019-08-23 西安电子科技大学 一种基于伪框架的强弱信号测向方法
CN112949100A (zh) * 2020-11-06 2021-06-11 中国人民解放军空军工程大学 一种机载雷达抗主瓣干扰的方法
CN113466899A (zh) * 2021-08-13 2021-10-01 电子科技大学 高信噪比环境下基于小快拍数的导航接收机波束形成方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
苏成晓;罗景青;: "一种强干扰下均匀圆阵弱信号测向新方法", 现代防御技术, no. 05 *

Also Published As

Publication number Publication date
CN116973834B (zh) 2024-01-30

Similar Documents

Publication Publication Date Title
CN106526530B (zh) 基于传播算子的2-l型阵列二维doa估计算法
CN109407055B (zh) 基于多径利用的波束形成方法
CN112051540B (zh) 一种快速高精度测向方法
CN113391281B (zh) 一种信号失配时的极化雷达加权扩展目标检测方法与系统
JPH11231033A (ja) 到来波推定方法
CN105334435B (zh) 一种基于任意阵形的自适应局部放电超声监测方法
Liu et al. Discrimination of mainlobe deceptive target with meter-wave FDA-MIMO radar
CN109143230A (zh) 一种应用于s波段穿墙雷达的二维波达方向估计方法
CN116973834B (zh) 一种强干扰下基于协同波束形成置信度的测向方法
Muhamed et al. Comparison of conventional subspace based DOA estimation algorithms with those employing property-restoral techniques: simulation and measurements
Fang et al. DOA estimation for sources with large power differences
CN109655093A (zh) 传感器阵列的幅相误差有源校正方法
CN113075613B (zh) 一种卫星移动通信系统中的干扰测向方法
Shirai et al. An Experimental Investigation of the MUSIC-based Wireless Position Location using LCX antenna at 5GHz band
CN116973835B (zh) 一种多信号混叠的伪空间谱累积测向方法
CN115792792B (zh) 极低样本下基于协方差矩阵迭代估计的测向方法与系统
CN113985346B (zh) 一种复杂电磁环境下多径目标doa估计方法
CN114296040B (zh) 一种极化切换雷达抗主瓣干扰优化方法
Jhang et al. Computationally efficient DOA estimation for massive uniform linear array
CN116819430B (zh) 一种强辐射源背景下的同频信号的测向方法
Chang et al. Polyphase Decomposition Based Space Frequency Adaptive Processing in Antenna Array GNSS Receivers
Du et al. AE-STAP Algorithm for Space-Time Anti-Jamming
Hill Multiple signal DF using superresolution: a practical assessment
CN117590319A (zh) 一种随机被动匹配滤波的弱信号高分辨测向方法
He DOA Estimation Approach Based on Parallel Stochastic Resonance System and Signal Correlation Calculation

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
GR01 Patent grant
GR01 Patent grant