CN111239677A - 基于数字阵列的多波束被动单脉冲测角方法 - Google Patents
基于数字阵列的多波束被动单脉冲测角方法 Download PDFInfo
- Publication number
- CN111239677A CN111239677A CN202010003975.5A CN202010003975A CN111239677A CN 111239677 A CN111239677 A CN 111239677A CN 202010003975 A CN202010003975 A CN 202010003975A CN 111239677 A CN111239677 A CN 111239677A
- Authority
- CN
- China
- Prior art keywords
- target
- array
- data
- value
- angle
- 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
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000001514 detection method Methods 0.000 claims abstract description 9
- 239000013598 vector Substances 0.000 claims description 22
- 238000005070 sampling Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 2
- 230000003044 adaptive effect Effects 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 10
- 238000002474 experimental method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000007123 defense Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000003471 anti-radiation Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000004088 simulation Methods 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)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于数字阵列的多波束被动单脉冲测角方法,以提高传统多波束比幅测向方法的稳健性。具体步骤包括:步骤1,电子侦察设备将多阵元数据进行数字同时多波束处理,完成指定空域的信号截获;步骤2,对上述多个波束的数据进行门限检测,确定目标对应的波束序号和时刻序号;步骤3,采用被动自适应单脉冲方法对存在目标的数据进行处理,得到目标的角度值。本发明提高了传统多波束比幅测向方法的精度和稳健性,提升了电子侦察设备的测向能力。
Description
技术领域
本发明涉及电子侦察无源测向技术,具体涉及一种基于数字阵列的多波束被动单脉冲测角方法。
背景技术
无源测向技术是一种不需自主地发射电磁波,完全依靠被动接收辐射源的电磁波,测量出目标方位的技术。无源测向技术因为不自主地发射电磁波,所以隐身性能较好,同时还具有探测距离远、安全性能高、隐蔽接收电磁波等优点,是现代防空预警体系的重要组成部分。尤其是在现代隐身技术、超低空突防技术、反辐射导弹技术和电子干扰技术等高速发展的情况下,更加突显出被动测向技术的优势,在航空航海、侦察预警、测量测绘等方面得到了广泛的应用。
传统的无源测向方法包括比幅测向法,相位干涉仪测向法、空间谱估计测向法等。比幅测向法为全方向振幅单脉冲测向技术,其测向原理是通过比较两个或多个独立的方向性天线所接收的来自远场的电磁波信号强弱来估计获取目标的方位信息,然而该方法易受通道幅度特性不一致的影响,导致比幅测向法的测向精度不高。相位干涉仪测向法是通过阵列天线接收来自远场同一信号的相位差确定辐射源的入射方向,然而该方法对通道相位特性敏感,且缺乏同时多信号测量能力。空间谱估计测向法是利用阵元在天线中的不同位置来获得所接收到信号的各种参数,结合相应的数学运算,对来波信号进行估计,从而得到辐射源的位置信息。然而该方法计算量高,且易受阵列非理想因素影响,导致空间谱估计在电子战中的实用化还面临着诸多难题。
发明内容
本发明的目的在于提供一种基于数字阵列的多波束被动单脉冲测角方法,以提高传统多波束比幅测向方法的稳健性,进而提升电子侦察设备的测向能力。
实现本发明目的的技术解决方案为:一种基于数字阵列的多波束被动单脉冲测角方法,步骤如下:
步骤1,电子侦察设备将多阵元数据进行数字同时多波束处理,完成指定空域的信号截获,转入步骤2;
步骤2,对上述多个波束的数据进行门限检测,确定目标对应的波束序号和时刻序号,转入步骤3;
步骤3,采用被动自适应单脉冲方法对存在目标的数据进行处理,得到目标的角度值。
本发明与现有技术相比,其显著优点在于:
(1)本发明基于数字阵列,采用数字同时多波束的方法,相比于传统多波束比幅测向方法提高了空域覆盖范围。
(2)本发明通过被动单脉冲方法精测量目标角度,相比于传统多波束比幅测向方法需要事先存储鉴角曲线参数的问题,本方法自适应程度较高。
附图说明
图1为本发明基于数字阵列的多波束被动单脉冲测角方法的流程图。
图2为实施例中角度估计均方根误差随目标输入信噪比变化曲线图。
具体实施方式
下面结合附图对本发明作进一步详细描述。
结合图1,本发明所述的一种基于数字阵列的多波束被动单脉冲测角方法,包括以下步骤:
步骤1,电子侦察设备将多阵元数据进行数字同时多波束处理,完成指定空域的信号截获,具体如下:
1a)电子侦察设备对N元线阵的接收回波采样,得到回波数据矢量xk,其中xk为N×1维的矢量,下标k表示采样时刻,k=1,…,K,K为时间段内总的采样点;
1b)在侦察***要求的空域覆盖范围内均匀采样M点,得到一组波束指向点集[θ1,…,θM];
1c)利用波束指向点集,计算波束指向权矢量wm
其中,n为阵元序号,m为波束序号,d为阵元间距,λ为信号波长,T表示矩阵的转置;
1d)对接收回波数据进行数字同时多波束处理,则第m个波束在k时刻的数据ym,k为
其中,H表示矩阵的共轭转置。
转入步骤2。
步骤2,对上述多个波束的数据进行门限检测,确定目标对应的波束序号和时刻序号,具体如下:
2a)计算每一个接收波束数据的输出功率值,即
pm,k=(ym,k)2
其中,pm,k为第m个波束在k时刻的接收数据功率值;
2b)将上述计算的输出功率值与预设的门限值η进行比较,即
其中,η为设定的门限值,H1表示检测统计量大于门限值,H0表示检测统计量小于门限值;
2c)通过上述检测,假定有Q个数据的功率值超过门限,这些数据对应的波束序号与时刻序号集为{(mq,kq)}q=1,…Q。
转入步骤3。
步骤3,采用被动自适应单脉冲方法对存在目标的数据进行处理,得到目标的角度值,具体如下:
3a)令q=1,得到对应的波束序号与时间序号为(mq1,kq1),根据波束序号与波束指向的对应关系,得到相应的波束指向与时间序号为(θq1,kq1);
其中,sΣ为阵列和导向矢量,fq1=sinθq1为空间频率,sΔ为阵列差导向矢量,N为阵元数,j为虚数符号,diag()表示将矢量对角化,linspace(a,b,c)表示在范围(a,b)内均匀采样c点;
3c)计算和波束权矢量wΣ和差波束权矢量wΔ
3d)得到权矢量后,利用下式估计角度
其中,fq1,t为目标实际的空间频率值,Faa为主波束中心角度对应的斜率,Fa为主波束中心角度处的自适应差和比输出,Δu为目标相对于天线主波束的偏差,自适应差和比输出和斜率的具体形式为
其中Real()表示取复数的实部,μa为偏差校准系数;
3e)得到目标的空间频率值后,计算目标的实际角度估计值θq1,t
θq1,t=asin(fq1,t)
其中asin()表示求取反正弦值;
3f)令q=2,…,Q,重复执行之步骤3a~3e,得到所有的目标角度估计值[θq1,t,…,θqQ,t]。
下面结合仿真实验对本发明的效果做进一步说明。
实验中接收阵列为等距线阵,阵元数为10,阵元间距为0.1m。目标信号波长为0.2m,目标方位角为1.5°,目标信噪比SNR从5dB变化到20dB。阵元幅度误差为5%,服从高斯分布。实验中将传统多波束比幅测向与数字阵列多波束被动单脉冲测向方法进行对比分析。实验中以均方根误差为准则衡量两种方法性能。采用蒙特卡罗实验计算均方根误差,实验次数为100次。
实验结果如图2所示。由图2可以看出数字阵列多波束被动单脉冲测向方法和传统多波束比幅测向方法的均方根误差均随目标信噪比的增大而减小,且数字阵列多波束被动单脉冲测向方法性能优于传统多波束比幅测向方法。
Claims (4)
1.一种基于数字阵列的多波束被动单脉冲测角方法,其特征在于,步骤如下:
步骤1,电子侦察设备将多阵元数据进行数字同时多波束处理,完成指定空域的信号截获,转入步骤2;
步骤2,对上述多个波束的数据进行门限检测,确定目标对应的波束序号和时刻序号,转入步骤3;
步骤3,采用被动自适应单脉冲方法对存在目标的数据进行处理,得到目标的角度值。
2.根据权利要求1所述的基于数字阵列的多波束被动单脉冲测角方法,其特征在于,所述步骤1,电子侦察设备将多阵元数据进行数字同时多波束处理,完成指定空域的信号截获,具体如下:
1a)电子侦察设备对N元线阵的接收回波采样,得到回波数据矢量xk,其中xk为N×1维的矢量,下标k表示采样时刻,k=1,…,K,K为时间段内总的采样点;
1b)在***要求的空域覆盖范围内均匀采样M点,得到一组波束指向点集[θ1,…,θM];
1c)利用波束指向点集,计算波束指向权矢量wm
其中,n为阵元序号,m为波束序号,d为阵元间距,λ为信号波长,T表示矩阵的转置;
1d)进行数字同时多波束处理,则第m个波束在k时刻的数据ym,k为
其中,H表示矩阵的共轭转置。
4.根据权利要求3所述的基于数字阵列的多波束被动单脉冲测角方法,其特征在于,所述步骤3,采用被动自适应单脉冲方法对存在目标的数据进行处理,得到目标的角度值,具体如下:
3a)令q=1,得到对应的波束序号与时间序号为(mq1,kq1),根据波束序号与波束指向的对应关系,得到相应的波束指向与时间序号为(θq1,kq1);
其中,sΣ为阵列和导向矢量,fq1=sinθq1为空间频率,sΔ为阵列差导向矢量,λ为信号波长,d为阵元间距,N为阵元数,j为虚数符号,diag()表示将矢量对角化,linspace(a,b,c)表示在范围(a,b)内均匀采样c点;
3c)计算和波束权矢量wΣ和差波束权矢量wΔ
3d)得到权矢量后,利用下式估计角度
其中,fq1,t为目标实际的空间频率值,Faa为主波束中心角度对应的斜率,Fa为主波束中心角度处的自适应差和比输出,Δu为目标相对于天线主波束的偏差,自适应差和比输出和斜率的具体形式为
其中Real()表示取复数的实部,μa为偏差校准系数;
3e)得到目标的空间频率值后,计算目标的实际角度估计值θq1,t
θq1,t=asin(fq1,t)
其中asin()表示求取反正弦值;
3f)令q=2,…,Q,重复执行之步骤3a~3e,得到所有的目标角度估计值[θq1,t,…,θqQ,t]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010003975.5A CN111239677B (zh) | 2020-01-03 | 2020-01-03 | 基于数字阵列的多波束被动单脉冲测角方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010003975.5A CN111239677B (zh) | 2020-01-03 | 2020-01-03 | 基于数字阵列的多波束被动单脉冲测角方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111239677A true CN111239677A (zh) | 2020-06-05 |
CN111239677B CN111239677B (zh) | 2023-10-31 |
Family
ID=70872311
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010003975.5A Active CN111239677B (zh) | 2020-01-03 | 2020-01-03 | 基于数字阵列的多波束被动单脉冲测角方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111239677B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112187315A (zh) * | 2020-08-25 | 2021-01-05 | 南京理工大学 | 一种大规模扩频通信数字阵列同时多用户快速角度估计方法 |
CN112881975A (zh) * | 2021-01-08 | 2021-06-01 | 电子科技大学 | 基于子阵特征矩阵联合对角化的单脉冲和差波束测角方法 |
CN113820653A (zh) * | 2021-08-04 | 2021-12-21 | 西安电子科技大学 | 基于动态和差波束的米波雷达低仰角目标doa估计方法 |
CN114428225A (zh) * | 2022-04-06 | 2022-05-03 | 中国电子科技集团公司第二十九研究所 | 一种基于准匹配滤波的多辐射源到达角度测量方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005121581A (ja) * | 2003-10-20 | 2005-05-12 | Mitsubishi Electric Corp | レーダ装置 |
JP2013174498A (ja) * | 2012-02-24 | 2013-09-05 | Mitsubishi Electric Corp | モノパルス測角装置およびモノパルス測角方法 |
CN104166136A (zh) * | 2014-07-11 | 2014-11-26 | 河海大学 | 一种基于干扰子空间跟踪的高效自适应单脉冲测角方法 |
CN109669178A (zh) * | 2018-12-25 | 2019-04-23 | 中国航天科工集团八五研究所 | 一种星载三阵元单脉冲两维测向方法 |
CN110058193A (zh) * | 2019-04-22 | 2019-07-26 | 西安空间无线电技术研究所 | 一种基于单接收通道的数字多波束测角方法及*** |
-
2020
- 2020-01-03 CN CN202010003975.5A patent/CN111239677B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005121581A (ja) * | 2003-10-20 | 2005-05-12 | Mitsubishi Electric Corp | レーダ装置 |
JP2013174498A (ja) * | 2012-02-24 | 2013-09-05 | Mitsubishi Electric Corp | モノパルス測角装置およびモノパルス測角方法 |
CN104166136A (zh) * | 2014-07-11 | 2014-11-26 | 河海大学 | 一种基于干扰子空间跟踪的高效自适应单脉冲测角方法 |
CN109669178A (zh) * | 2018-12-25 | 2019-04-23 | 中国航天科工集团八五研究所 | 一种星载三阵元单脉冲两维测向方法 |
CN110058193A (zh) * | 2019-04-22 | 2019-07-26 | 西安空间无线电技术研究所 | 一种基于单接收通道的数字多波束测角方法及*** |
Non-Patent Citations (2)
Title |
---|
曾操 等: "相控阵子阵级和差多波束测角方法" * |
陈亮 等: "改进贝叶斯聚焦的宽带自适应单脉冲算法" * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112187315A (zh) * | 2020-08-25 | 2021-01-05 | 南京理工大学 | 一种大规模扩频通信数字阵列同时多用户快速角度估计方法 |
CN112881975A (zh) * | 2021-01-08 | 2021-06-01 | 电子科技大学 | 基于子阵特征矩阵联合对角化的单脉冲和差波束测角方法 |
CN112881975B (zh) * | 2021-01-08 | 2023-09-08 | 电子科技大学 | 基于子阵特征矩阵联合对角化的单脉冲和差波束测角方法 |
CN113820653A (zh) * | 2021-08-04 | 2021-12-21 | 西安电子科技大学 | 基于动态和差波束的米波雷达低仰角目标doa估计方法 |
CN114428225A (zh) * | 2022-04-06 | 2022-05-03 | 中国电子科技集团公司第二十九研究所 | 一种基于准匹配滤波的多辐射源到达角度测量方法及装置 |
CN114428225B (zh) * | 2022-04-06 | 2022-06-14 | 中国电子科技集团公司第二十九研究所 | 一种基于准匹配滤波的多辐射源到达角度测量方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111239677B (zh) | 2023-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10571544B2 (en) | Direction finding using signal power | |
CN111239677B (zh) | 基于数字阵列的多波束被动单脉冲测角方法 | |
CN109738854B (zh) | 一种天线阵列来波方向的到达角估计方法 | |
Davis et al. | Angle estimation with adaptive arrays in external noise fields | |
CN108549059B (zh) | 一种复杂地形条件下的低空目标仰角估计方法 | |
CN102544755B (zh) | 一种基于强散射点的均匀线阵校准方法 | |
CN109946664B (zh) | 一种主瓣干扰下的阵列雷达导引头单脉冲测角方法 | |
CN103412286B (zh) | 基于mimo雷达的发射极化优化doa估计方法 | |
CN110673086A (zh) | 一种基于数字阵列雷达的二维角度超分辨方法 | |
CN107390197B (zh) | 基于特征空间的雷达自适应和差波束测角方法 | |
CN109765519B (zh) | 一种模数混合天线阵列的角度估计方法 | |
CN104502904B (zh) | 一种鱼雷自导波束锐化方法 | |
CN110426670B (zh) | 基于tls-cs的外辐射源雷达超分辨doa估计方法 | |
CN110196417B (zh) | 基于发射能量集中的双基地mimo雷达角度估计方法 | |
KR101807499B1 (ko) | 상호간섭효과를 보상한 원형 배열 안테나 기반의 도래각 추정 장치 및 추정 방법 | |
CN113126047A (zh) | 基于远场点源的自适应通道校正方法 | |
CN109212466B (zh) | 一种基于量子蜻蜓演化机制的宽带测向方法 | |
CN113835063A (zh) | 一种无人机阵列幅相误差与信号doa联合估计方法 | |
CN111693979B (zh) | 一种基于泰勒展开的数字阵列单脉冲测角方法 | |
Wu et al. | Array pattern synthesis and robust beamforming for a complex sonar system | |
CN112415469B (zh) | 一种两维数字阵列雷达快速干扰测向方法 | |
Zhou et al. | Research on interferometer direction finding technology based on digital beam forming | |
Wang et al. | A practical calibration method of linear UHF Yagi arrays for ship target detection application | |
Xiong et al. | Space time adaptive processing for airborne MIMO radar based on space time sampling matrix | |
CN112333629A (zh) | 一种互耦未知条件下分布式阵列多目标定位方法 |
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 |