CN111610487B - 均匀线阵虚拟合成短时大孔径高精度测向方法 - Google Patents
均匀线阵虚拟合成短时大孔径高精度测向方法 Download PDFInfo
- Publication number
- CN111610487B CN111610487B CN202010436801.8A CN202010436801A CN111610487B CN 111610487 B CN111610487 B CN 111610487B CN 202010436801 A CN202010436801 A CN 202010436801A CN 111610487 B CN111610487 B CN 111610487B
- Authority
- CN
- China
- Prior art keywords
- antenna array
- time
- speed switch
- array element
- short
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 9
- 238000005070 sampling Methods 0.000 claims abstract description 26
- 238000006243 chemical reaction Methods 0.000 claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims description 17
- 125000004122 cyclic group Chemical group 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 10
- 230000005855 radiation Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 5
- 230000002194 synthesizing effect Effects 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了均匀线阵虚拟合成短时大孔径高精度测向方法,具体步骤如下:利用高速开关逻辑驱动器,来驱动一个高速开关,以Δt为高速开关转换时间,轮巡切换接通各个天线阵元,循环的获取到各个天线阵元的接收数据;对获取的各个天线阵元的接收数据进行相位补偿校准和幅度均衡控制,来形成虚拟短时大孔径,并获取到基于虚拟短时大孔径得到的取样数据;对获取的基于虚拟短时大孔径得到的取样数据进行下变频和数字化处理;本发明是基于均匀线阵,利用高速开关轮巡切换接通各个天线阵元,经过相位补偿校准和幅度均衡控制,虚拟形成短时大孔径,并由虚拟短时大孔径得到的取样数据,通过信号处理来实现高精度测向。
Description
技术领域
本发明涉及电子侦察技术领域,是对辐射源信号的测向方法,具体为均匀线阵虚拟合成短时大孔径高精度测向方法。
背景技术
辐射源信号的波达角估计是利用电子侦察接收到的辐射源信号,通过信号处理来估计辐射源信号的到达方向;目前,辐射源信号的波达角估计方法主要有振幅法测向、相位法测向和空间谱估计测向等方式。
振幅法测向是利用接收到的信号幅度信息进行测向,其测向精度与阵列的方向图有关,当阵元间距一定时,阵元数目越大,阵列的孔径也越大,方向图主瓣越窄,测向精度越高,这种方法的测向精度相对较低,现代测向***的应用较少;相位法测向是利用天线阵接收到的相位差信息进行测向,典型的有多基线相位干涉仪测向,其利用多基线解决不模糊测向范围及测向精度间的矛盾,短基线决定不模糊测向范围,而测向精度与最长基线长度成反比,也就是阵列的基线越长,测向精度越高,这种方法测向精度相对较高,现代工程方面的应用较多;空间谱估计测向是利用现代谱估计技术进行测向,其测向精度与阵列孔径相关,阵列孔径越大,测向精度越高,这类方法理论上可实现极高的测向精度,但现代工程方面的应用尚有许多问题待解决,即工程应用较少;
上述辐射源信号测向方法的共同特点是测向精度与天线阵列的孔径大小密切相关,要获得高精度测向结果,必须保证阵列的孔径足够大,这也导致很多现代雷达电子***的天线做得越来越庞大;
而在天线孔径一定的情况下,为了提高测向***的测向精度,来提出一种均匀线阵虚拟合成短时大孔径高精度测向方法,该方法是基于均匀线阵,利用高速开关轮巡切换接通各个天线阵元,经过相位补偿校准和幅度均衡控制,虚拟形成短时大孔径,并由虚拟短时大孔径得到的取样数据,通过信号处理来实现高精度测向;
为了解决上述缺陷,现提供一种技术方案。
发明内容
本发明的目的在于提供均匀线阵虚拟合成短时大孔径高精度测向方法,本发明是基于均匀线阵,利用高速开关轮巡切换接通各个天线阵元,经过相位补偿校准和幅度均衡控制,虚拟形成短时大孔径,并由虚拟短时大孔径得到的取样数据,通过信号处理来实现高精度测向。
本发明的目的可以通过以下技术方案实现:
均匀线阵虚拟合成短时大孔径高精度测向方法,具体步骤如下:
步骤一:利用高速开关逻辑驱动器,来驱动一个高速开关,以Δt为高速开关转换时间,轮巡切换接通各个天线阵元,循环的获取到各个天线阵元的接收数据;
步骤二:对获取的各个天线阵元的接收数据进行相位补偿校准和幅度均衡控制,来形成虚拟短时大孔径,并获取到基于虚拟短时大孔径得到的取样数据,且相位补偿校准为现有的相位补偿的校准方式,幅度均衡控制也可通过现有的幅度均衡器来设置操作,主要是将接收数据来形成虚拟短时大孔径,并获取其中的取样数据;
步骤三:对获取的基于虚拟短时大孔径得到的取样数据进行下变频和数字化处理,且取样数据的下变频和数字化处理,以及后续的广义傅里叶变换方式均为现有技术,并通过广义傅里叶变换,得到目标信号的空间谱,再依据谱峰搜索,完成对目标信号到达方向的估计,实现目标信号的高精度测向,即利用广义傅里叶变换后的空间谱,来对其到达方位的一维搜索,实现对到达方向的估计和高精度测向。
本发明的有益效果:
本发明是利用高速开关逻辑驱动器,来驱动一个高速开关,以极短的时间间隔,从均匀线阵的第一个天线阵元开始,顺序的切换接通各个天线阵元,循环获取各个天线阵元的接收数据,再对接收数据进行相位补偿校准和幅度均衡控制,形成虚拟短时大孔径得到的取样数据,基于取样数据,经广义傅里叶变换方式的信号处理,来实现高精度测向;即本发明与传统的均匀线阵测向方法相比较,能够形成更大的虚拟孔径;且本发明的测向精度高,在相同阵列条件下,比传统的均匀线阵测向方法的测向误差降低至少50%。
附图说明
为了便于本领域技术人员理解,下面结合附图对本发明作进一步的说明;
图1为本发明的天线阵列布阵示意图;
图2为本发明的虚拟短时大孔径原理图;
图3为本发明的空间谱估计结果示意图;
图4为本发明的测向误差与信噪比关系示意图。
具体实施方式
数学模型:
假定观测站的位置固定,测向天线为一维均匀线阵,其阵元数为L,阵元间距为d,测向阵列模型如图1所示,观测目标满足远场条件;
其中,s(t,x)为时间t和位移x的函数,当x固定时,即阵元空间不动,s(t,x)为时间的函数,ωs决定单位时间内的电波相位变化的弧度数;当t固定时,即时间不动,s(t,x)为空间位移的函数,Ω0决定单位长度上的电波相位变化的弧度数;ω0为信号角频率,而为空间角频率;
如果接收机是在运动中接收信号,则位置x是时间的函数,接收信号模型可表示为:
s(t,x)=Aexp[jωst+jΩ0x(t)+Φ0]
其中,x(t)表示t时刻的接收天线阵元所在的位置。
虚拟合成短时大孔径高精度测向原理:
步骤一:利用高速开关逻辑驱动器,来驱动一个高速开关,以Δt为高速开关转换时间,轮巡切换接通各个天线阵元,循环的获取到各个天线阵元的接收数据,其具体方式如下:
利用高速开关逻辑驱动器,来驱动一个高速开关从均匀线阵的第一个天线阵元开始,顺序的切换接通各个天线阵元,获取各个天线阵元的接收数据;当所有天线阵元取样结束后,控制该高速开关继续切换到第一个天线阵元,重复进行顺序取样,如此循环反复;
当第一个天线阵元的位置为x=0,高速开关按1#,2#,…,L#天线阵元的顺序轮巡切换接通各个天线阵元,高速开关转换时间为Δt,tn=nΔt时刻的高速开关位于x(tn)处,而x(tn)的表示形式如下:
x(tn)=d mod(n/L)
其中,n表示tn时刻为第n次取样;
即第n次取样时刻的各个天线阵元的接收数据表示如下:
s(n)=Aexp[jωsnΔt+jΩ0x(tn)+Φ0]
当0≤n<L时
s(n)=Aexp[jωsnΔt+jΩ0dn+Φ0]
=Aexp[jωsnΔt+jω0n+Φ0]
当kL≤n<(k+1)L-1时
s(n)=Aexp[jωsnΔt+jΩ0d(n-kL)+Φ0]
=Aexp[jωsnΔt+jω0n+Φ0-jkω0L]
则其与无限长的线性阵列相比较,在n=kL时刻产生一个相位突变,其值为kω0L,其中,k=[n/L]为取整函数;
步骤二:对获取的各个天线阵元的接收数据进行相位补偿校准和幅度均衡控制,来形成虚拟短时大孔径,并获取到基于虚拟短时大孔径得到的取样数据,且相位补偿校准为现有的相位补偿的校准方式,幅度均衡控制也可通过现有的幅度均衡器来设置操作,主要是将接收数据来形成虚拟短时大孔径,并获取其中的取样数据;
即虚拟短时大孔径原理图为图2所示,基于虚拟短时大孔径得到的取样数据,再将其传输至接收机,通过信号处理来实现高精度测向;
步骤三:对获取的基于虚拟短时大孔径得到的取样数据进行下变频和数字化处理,且取样数据的下变频和数字化处理,以及后续的广义傅里叶变换方式均为现有技术,并通过广义傅里叶变换,得到目标信号的空间谱,再依据谱峰搜索,完成对目标信号到达方向的估计,实现目标信号的高精度测向,其具体方式如下:
对s(n)进行广义傅里叶变换,得到目标信号的空间谱:
验证实验:
如图3所示,其为SNR=0dB时的空间谱,如图4所示为不同信噪比条件下的测向误差,则当信噪比SNR>-10dB时,本发明的测向误差低于0.05°;
将本发明与传统方法做性能比较,基于上述的实验条件,选择典型的基于均匀线阵的相位法测向和空间谱测向进行对比实验,实验中的信号入射角在(-60,60)范围内选择7个不同的角度,在信噪比SNR=-10dB时的三种不同方法的测向误差实验结果如表1所示:
表1-三种不同方法的测向误差实验结果表
实验结果表明,均匀线阵虚拟合成短时大孔径高精度测向方法,当信噪比SNR>-10dB时,其测向误差低于0.05°,而均匀线阵的传统测向方法的测向误差高于0.1°;即本发明通过虚拟合成短时大孔径,可实现高精度测向,在相同阵列条件下,比传统的均匀线阵测向方法的测向误差降低至少50%。
以上内容仅仅是对本发明结构所作的举例和说明,所属本技术领域的技术人员对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,只要不偏离发明的结构或者超越本权利要求书所定义的范围,均应属于本发明的保护范围。
Claims (2)
1.均匀线阵虚拟合成短时大孔径高精度测向方法,其特征在于,具体步骤如下:
步骤一:利用高速开关逻辑驱动器,来驱动一个高速开关,以Δt为高速开关转换时间,轮巡切换接通各个天线阵元,循环的获取到各个天线阵元的接收数据;
步骤二:对获取的各个天线阵元的接收数据进行相位补偿校准和幅度均衡控制,来形成虚拟短时大孔径,并获取到基于虚拟短时大孔径得到的取样数据;
步骤三:对获取的基于虚拟短时大孔径得到的取样数据进行下变频和数字化处理,并通过当Ω=Ω0时,即利用广义傅里叶变换后的空间谱,来对其到达方位的一维搜索,得到目标信号的空间谱,再依据谱峰搜索,完成对目标信号到达方向的估计,实现目标信号的高精度测向;
式中s(n)=Aexp[jωsnΔt+jΩ0dn+Φ0]=Aexp[jωsnΔt+jΩ0n+Φ0];
j为虚数部;ωs决定单位时间内的电波相位变化的弧度数;
Ω0决定单位长度上的电波相位变化的弧度数;
Ω为长度上的电波相位变化的弧度数;
x(tn)=d mod(n/L),其中L为阵元数,d为阵元间距,tn为第n次取样的时刻。
2.根据权利要求1所述的均匀线阵虚拟合成短时大孔径高精度测向方法,其特征在于,所述步骤一中的循环的获取到各个天线阵元的接收数据,其具体方式如下:
利用高速开关逻辑驱动器,来驱动一个高速开关从均匀线阵的第一个天线阵元开始,顺序的切换接通各个天线阵元,循环获取各个天线阵元的接收数据;当所有天线阵元取样结束后,控制该高速开关继续切换到第一个天线阵元,重复进行顺序取样,如此循环反复;
当第一个天线阵元的位置为x=0,高速开关按1#,2#,…,L#天线阵元的顺序轮巡切换接通各个天线阵元,高速开关转换时间为Δt,tn=nΔt时刻的高速开关位于x(tn)处,
第n次取样时刻的各个天线阵元的接收数据表示如下:
s(n)=Aexp[jωsnΔt+jΩ0x(tn)+Φ0]
当0≤n<L时
s(n)=Aexp[jωsnΔt+jΩ0dn+Φ0]
=Aexp[jωsnΔt+jΩ0n+Φ0]
当kL≤n<(k+1)L-1时
s(n)=Aexp[jωsnΔt+jΩ0d(n-kL)+Φ0]
=Aexp[jωsnΔt+jω0n+Φ0-jkω0L]
则其与无限长的线性阵列相比较,在n=kL时刻产生一个相位突变,其值为kω0L,其中,k=[n/L]为取整函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010436801.8A CN111610487B (zh) | 2020-05-21 | 2020-05-21 | 均匀线阵虚拟合成短时大孔径高精度测向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010436801.8A CN111610487B (zh) | 2020-05-21 | 2020-05-21 | 均匀线阵虚拟合成短时大孔径高精度测向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111610487A CN111610487A (zh) | 2020-09-01 |
CN111610487B true CN111610487B (zh) | 2022-08-09 |
Family
ID=72198947
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010436801.8A Active CN111610487B (zh) | 2020-05-21 | 2020-05-21 | 均匀线阵虚拟合成短时大孔径高精度测向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111610487B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112821905B (zh) * | 2020-12-30 | 2022-01-14 | 北京航空航天大学杭州创新研究院 | 一种基于evs的ads-b接收机数字增益补偿***及其补偿方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007039774A1 (en) * | 2005-10-06 | 2007-04-12 | Roke Manor Research Limited | Array calibration with signals of opportunity |
CN111903218B (zh) * | 2010-03-23 | 2012-05-16 | 中国电子科技集团公司第五十四研究所 | 基于中国剩余定理的干涉仪测向方法 |
CN106772294A (zh) * | 2016-12-30 | 2017-05-31 | 南京长峰航天电子科技有限公司 | 一种特殊辐射信号模拟装置及实现方法 |
KR20180113267A (ko) * | 2017-04-06 | 2018-10-16 | 서울대학교산학협력단 | 균일 선형 배열 안테나 기반의 가상 수신신호 생성을 이용한 도래각 추정 방법 및 장치 |
CN110515033A (zh) * | 2019-08-09 | 2019-11-29 | 南京航空航天大学 | 一种基于Toeplitz矩阵恢复的欠通道测向***及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11456809B2 (en) * | 2018-05-31 | 2022-09-27 | Electronics And Telecommunications Research Institute | Position estimation method for estimating position of interference signal source and position estimation system for performing the method |
-
2020
- 2020-05-21 CN CN202010436801.8A patent/CN111610487B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007039774A1 (en) * | 2005-10-06 | 2007-04-12 | Roke Manor Research Limited | Array calibration with signals of opportunity |
CN111903218B (zh) * | 2010-03-23 | 2012-05-16 | 中国电子科技集团公司第五十四研究所 | 基于中国剩余定理的干涉仪测向方法 |
CN106772294A (zh) * | 2016-12-30 | 2017-05-31 | 南京长峰航天电子科技有限公司 | 一种特殊辐射信号模拟装置及实现方法 |
KR20180113267A (ko) * | 2017-04-06 | 2018-10-16 | 서울대학교산학협력단 | 균일 선형 배열 안테나 기반의 가상 수신신호 생성을 이용한 도래각 추정 방법 및 장치 |
CN110515033A (zh) * | 2019-08-09 | 2019-11-29 | 南京航空航天大学 | 一种基于Toeplitz矩阵恢复的欠通道测向***及方法 |
Non-Patent Citations (4)
Title |
---|
Direction of Arrival (DOA) Estimation Using Electrically Small Tuned Dipole Antennas;Seunghyeon 等;《IEEE Transactions on Antennas and Propagation》;20061113;第54卷(第11期);第3292-3301页 * |
Direction of arrival estimation with arbitrary virtual antenna arrays using low cost inertial measurement units;Muhammad Atif Yaqoob 等;《2013 IEEE International Conference on Communications Workshops (ICC)》;20131031;第2164-7038页 * |
智能天线波达方向估计算法的研究;董增军;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20090115(第01期);第I136-23页 * |
米波数字阵列雷达低仰角测高方法研究;朱伟;《中国优秀博硕士学位论文全文数据库(博士)信息科技辑》;20131115(第11期);第I136-56页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111610487A (zh) | 2020-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109633650B (zh) | 一种车载毫米波雷达多通道阵列天线幅相校正装置及方法 | |
US7692575B2 (en) | Radar target detection method and radar apparatus using the same | |
CN112887002B (zh) | 一种智能反射面辅助通信的波束域信道角度估计方法 | |
CN110429964B (zh) | 一种基于二维相控天线阵列的快速精确波束跟踪方法 | |
WO2004088347A1 (fr) | Procede d'estimation de la direction de l'arrivee d'un faisceau de fixation dans l'espace | |
CN110149126B (zh) | 一种3d-mimo***的波束赋形方法及波束赋形装置 | |
KR101807499B1 (ko) | 상호간섭효과를 보상한 원형 배열 안테나 기반의 도래각 추정 장치 및 추정 방법 | |
CN111610487B (zh) | 均匀线阵虚拟合成短时大孔径高精度测向方法 | |
CN113189592A (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
JP3629195B2 (ja) | 到来方向推定装置及び到来方向推定方法 | |
CN112858994B (zh) | 一种基于均匀圆阵的比幅测向方法 | |
CN113644941A (zh) | 一种基于大规模mimo接收阵列结构的模糊相位快速消除方法 | |
CN114079520A (zh) | 一种信号估计方法、设备、装置及存储介质 | |
JP5722026B2 (ja) | 到来方向推定方式 | |
CN107359922A (zh) | 一种波束扫描与doa相结合的邻居发现与精确对准方法 | |
JP4967218B2 (ja) | レーダ装置 | |
CN111190135A (zh) | 一种适用于任意阵列的校准方法 | |
JP3138728B2 (ja) | アレーアンテナの較正方法 | |
JP2005207837A (ja) | 方向探知機 | |
CN112731279B (zh) | 一种基于混合天线子阵列的到达角估计方法 | |
CN111934737B (zh) | 一种宽带数字阵列的时延补偿方法 | |
US11598843B2 (en) | Apparatus and method for analyzing radar signals | |
WO2008129114A1 (en) | An improved method for channel sounding and an apparatus using the method | |
JP3113837B2 (ja) | バトラーマトリクス形式のnポートモードフォーマのバイアスエラー減少方法及び装置 | |
WO2010026503A1 (en) | Angle of arrival (basis selection) and location estimation system |
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 |