CN107656238B - 一种基于双查找表方法的高频地波雷达测向新方法 - Google Patents
一种基于双查找表方法的高频地波雷达测向新方法 Download PDFInfo
- Publication number
- CN107656238B CN107656238B CN201710819351.9A CN201710819351A CN107656238B CN 107656238 B CN107656238 B CN 107656238B CN 201710819351 A CN201710819351 A CN 201710819351A CN 107656238 B CN107656238 B CN 107656238B
- Authority
- CN
- China
- Prior art keywords
- ratio
- beams
- output
- region
- value
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000002159 abnormal effect Effects 0.000 claims abstract description 28
- 238000010586 diagram Methods 0.000 claims abstract description 17
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 8
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 8
- 230000002547 anomalous effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 6
- 238000001514 detection method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process 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
- G01S3/16—Systems for determining direction or deviation from predetermined direction using amplitude comparison of signals derived sequentially from receiving antennas or antenna systems having differently-oriented directivity characteristics or from an antenna system having periodically-varied orientation of directivity characteristic
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;最后根据是否落入异常区来判定查找相应的表格。本发明的一种基于双查找表方法的高频地波雷达测向新方法,采用查表法进行测向,仅需与已经做好的归一化幅度表进行比对,即可获知测向角度,方法简单,计算量成数量级的形式减少,且在精度能够保持的情况下,可以大幅度节约时间成本。
Description
技术领域
本发明涉及雷达技术领域,特别是涉及一种基于双查找表方法的高频地波雷达测向新方法。
背景技术
高频地波雷达有着全天候工作、超视距探测等优点,在监测专属经济区、保护海洋环境等方面有这广泛应用。但该种雷达波长较长,对目标方位不敏感,导致目标测向精度较低。
现有方法包括数字波束合成法(DBF)、比幅测向方法等。DBF法中,来波方向上的导向矢量为其中α为目标方位角,λ为波长,d为阵元间距,N为阵元个数。经过数字波束合成(DBF)处理后得到输出y=w1 TBw。其中w1为复加权矢量,T表示转置,B为切比雪夫加权矩阵,可以增强方向图主瓣并抑制旁瓣。w1=w1(θ),θ可以以一定扫描精度遍历雷达测向的角度范围,当得到的输出y最大时,θ即为测得的目标方位角α。
DBF法测向的原理虽然简单,但该方法运算量较大,每个来波都得进行一次全部可能角度的计算,即一次目标测向可能得进行上千次运算才能得到目标方位角,且测向精度非常依赖扫描精度:扫描精度不够,测向就会发生很大偏差;精度提高了,运算次数也成倍上升;
比幅测向法中原本不同波束中心指向的方向图形状并不完全一样,用高斯函数近似的方向图始终与真实方向图存在一定偏差,这意味着传统比幅测向方法的假设与实际情况并不完全一致,这导致了传统比幅测向方法得出的测向结果始终存在误差。
发明内容
本发明的目的在于提供一种基于双查找表方法的高频地波雷达测向新方法。
一种基于双查找表方法的高频地波雷达测向新方法,其具体制作步骤为:
1.生成表格1和表格2,而表格2主要是针对来波方向落入异常区(解释在本段末)中设计生成的。主要步骤如下:
1.1设定雷达扫描起始角度θ0、雷达扫描末端角度θn-1和波束步长step的值,波束步长step的大小设定,会影响到波束个数n的大小,即n=(θn-1-θ0)/step,也会对实际的测向精度产生影响。
1.2确定波束位置和个数n后,与原始通道数据进行波束合成,并根据y=w1(θi)*B*w(α),i=0,1,2…n-1,计算得到波束合成的输出,并得到了雷达方向分布图,如下;
1.3根据雷达方向图分布比较输出结果,得到每一个来波方向的输出最大值ymax1和次大值ymax2(此时对应二者已定波束为相邻波束),用前一波束的输出比上后一波束的输出,即可得到当前假定来波的比值ratio=ymax1/ymax2or ratio=ymax2/ymax1;
1.4以扫描精度不断改变来波角度α,重复上述操作即可得到完整的比值表,该比值表即为表1;
1.5从初始扫描位置向左(或向右)移动半个波束步长,即雷达扫描起始角度变为θ0'=θ0+step/2,波束步长step和波束个数n不变,重复以上操作中的1.2—1.4,这样得到新的比值表,即表2;
2.得到比值表1的异常区R:
2.1每四个相邻波束间均会存在这一小异常区Ri,i=1,2,…n-2,以i-1,i,i+1,i+2这四个波束为例,得到第i-1个和第i+1个波束的交点对应的角度值β,查找表1得到角度值β对应的第i-1和第i个波束对应的比值ratioi;得到第i个和第i+2个波束的交点对应的角度值β',查找表1得到角度值β'对应的第i和第i+1个波束对应的比值ratioi+1,这样ratio<ratioi∪ratio>ratioi+1则组成了这一个小异常区Ri;
2.2而所有波束间的小异常区将会组成一个大异常区,使之所有的来波都不会落入这些小的异常区Ri,i=1,2,…n-2内,于是比值表1的异常区R=R1∪R2∪R3∪…∪Rn-2;
3.对于未知来波α,对表1的所有波束进行波束合成,找出所有输出的最大值和次大值,得到二者比值ratio1',这一比值亦是前一波束输出除以后一波束输出的结果,判断ratio1'是否落入阈值区R内,若没有则直接查找表1即可得出未知来波方向的角度值α;若落入阈值区R内,则重新对表2的所有波束进行波束合成,找出所有输出的最大值和次大值ratio2',查找表2即可得出未知来波方向的角度值α。
有益效果在于:本发明在真实方向图上,得到已知的不同相邻波束与来波之间的幅度比值,排序后得到最大值和次大值,并记录成表,测得所有角度的来波方向上的幅度比值,便可依据最为接近的幅度比值查到未知来波的方向角度,基于比幅法的思想给出算法的框架,而采用已知波束的实际方向图,避免因高斯函数作为近似方向图带来的误差;其归一化幅度比值表一经生成,在阵列间距不变的情况下,是不用调整的,无需大量计算,即不同方向来波可以即时与比值表进行比对即可准确获知目标方位;
增加了方法上的多样性,在大多数数据点采用查表法测向的情况下,对于不能采用查表法则采用较为耗时的DBF法;
本发明采用查表法进行测向,仅需与已经做好的归一化幅度表进行比对,即可获知测向角度,方法简单,计算量成数量级的形式减少,且在精度能够保持的情况下,可以大幅度节约时间成本。
附图说明
图1为本发明工作流程示意图;
图2为本发明雷达真实方向分布示意图;
图3为本发明异常区示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1-3,一种基于双查找表方法的高频地波雷达测向新方法,其具体步骤为:
实施例一
1.生成表格1和表格2,而表格2主要是针对来波方向落入异常区(解释在本段末)中设计生成的,主要步骤如下:
1.1设定雷达扫描起始角度θ0、雷达扫描末端角度θn-1和波束步长step的值,波束步长step的大小设定,会影响到波束个数n的大小,即n=(θn-1-θ0)/step,也会对实际的测向精度产生影响;
1.2确定波束位置和个数n后,与原始通道数据进行波束合成,并根据y=w1(θi)*B*w(α),i=0,1,2…n-1,计算得到波束合成的输出,并得到了雷达方向分布图;
1.3根据雷达方向图分布比较输出结果,得到每一个来波方向的输出最大值ymax1和次大值ymax2(此时对应二者已定波束为相邻波束),用前一波束的输出比上后一波束的输出,即可得到当前假定来波的比值ratio=ymax1/ymax2or ratio=ymax2/ymax1;
1.4以扫描精度不断改变来波角度α,重复上述操作即可得到完整的比值表,该比值表即为表1;
1.5从初始扫描位置向左(或向右)移动半个波束步长,即雷达扫描起始角度变为θ0'=θ0+step/2,波束步长step和波束个数n不变,重复以上操作中的1.2—1.4。这样得到新的比值表,即表2;
2.得到比值表1的异常区R:
2.1每四个相邻波束间均会存在这一小异常区Ri,i=1,2,…n-2,以i-1,i,i+1,i+2这四个波束为例,得到第i-1个和第i+1个波束的交点对应的角度值β,查找表1得到角度值β对应的第i-1和第i个波束对应的比值ratioi;得到第i个和第i+2个波束的交点对应的角度值β',查找表1得到角度值β'对应的第i和第i+1个波束对应的比值ratioi+1,这样ratio<ratioi∪ratio>ratioi+1则组成了这一个小异常区Ri;
2.2而所有波束间的小异常区将会组成一个大异常区,使之所有的来波都不会落入这些小的异常区Ri,i=1,2,…n-2内,于是比值表1的异常区R=R1∪R2∪R3∪…∪Rn-2;
3.对于未知来波α,对表1的所有波束进行波束合成,找出所有输出的最大值和次大值,得到二者比值ratio1',这一比值亦是前一波束输出除以后一波束输出的结果,判断ratio1'是否落入阈值区R内,若没有则直接查找表1即可得出未知来波方向的角度值α;若落入阈值区R内,则重新对表2的所有波束进行波束合成,找出所有输出的最大值和次大值ratio2',查找表2即可得出未知来波方向的角度值α。
本发明基于双查找表方法的高频地波雷达测向新方法具有以下优点:
1、本发明给出了查表法的思想,对所有来波角度进行扫描,得到每个角度与此相邻波束间的归一化幅度比值,将其记录成表格形式;对新的来波进行测向时仅需查询归一化幅度比值表,即可获知来波方向;
2、避免因高斯函数作为近似方向图带来的误差;其归一化幅度比值表一经生成,在阵列间距不变的情况下,是不用调整的,无需大量计算,即不同方向来波可以即时与比值表进行比对即可准确获知目标方位;
3、参考冗余的思想,设计出两套归一化幅度比值表可供选择,能有效避免因落入异常区的来波方向测得的结果发生异常这一情况。
4、采用查表法进行测向,仅需与已经做好的归一化幅度表进行比对,即可获知测向角度,方法简单,计算量成数量级的形式减少,且在精度能够保持的情况下,可以大幅度节约时间成本。
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其效物界定。
Claims (1)
1.一种基于双查找表方法的高频地波雷达测向新方法,其特征在于:所述基于双查找表方法的高频地波雷达测向新方法步骤为:
1.生成表格1和表格2,而表格2主要是针对来波方向落入异常区中设计生成的,主要步骤如下:
1.1设定雷达扫描起始角度θ0、雷达扫描末端角度θn-1和波束步长step的值,波束步长step的大小设定,会影响到波束个数n的大小,即n=(θn-1-θ0)/step,也会对实际的测向精度产生影响;
1.2确定波束位置和个数n后,与原始通道数据进行波束合成,并根据y=w1(θi)*B*w(α),i=0,1,2…n-1,计算得到波束合成的输出,并得到了雷达方向分布图;
1.3根据雷达方向图分布比较输出结果,得到每一个来波方向的输出最大值ymax1和次大值ymax2,此时对应二者已定波束为相邻波束,用前一波束的输出比上后一波束的输出,即可得到当前假定来波的比值ratio=ymax1/ymax2 or ratio=ymax2/ymax1;
1.4以扫描精度不断改变来波角度α,重复上述操作即可得到完整的比值表,该比值表即为表格1;
1.5从初始扫描位置向左,或向右,移动半个波束步长,即雷达扫描起始角度变为θ0'=θ0+step/2,波束步长step和波束个数n不变,重复以上操作中的1.2—1.4;这样得到新的比值表,即表格2;
对于异常区的解释,波束2和波束3的指向分别为-35度和-15度,在-35度和-15之间,绝大部分区域的最大值和次大值都是波束2和波束3,但是在波束2的指向线和阈值线1之间,最大值为波束2,次大值为波束1,即次大值为波束2和波束3之外的波束,波束2和阈值线1之间的区域就是异常区,因为这块区域在生成表格时用波束2和波束3作比值存入表格,而查找表时最大值次大值不是波束2和波束3,因此会引起测向错误,同理图中波束3和阈值线2之间也是异常区,如果将波束平移一半步长,波束异常区就会变成正常区域,因为使用两套波束,生成两个表格,如果比幅值正常,则用表格1,如果比幅值落入异常区,则用表格2;
2.得到比值表格1的异常区R:
2.1每四个相邻波束间均会存在这一小异常区Ri,i=1,2,…n-2,以i-1,i,i+1,i+2这四个波束为例,得到第i-1个和第i+1个波束的交点对应的角度值β,查找表格1得到角度值β对应的第i-1和第i个波束对应的比值ratioi;得到第i个和第i+2个波束的交点对应的角度值β',查找表格1得到角度值β'对应的第i和第i+1个波束对应的比值ratioi+1,这样ratio<ratioi∪ratio>ratioi+1则组成了这一个小异常区Ri;
2.2而所有波束间的小异常区将会组成一个大异常区,使之所有的来波都不会落入这些小的异常区Ri,i=1,2,…n-2内,于是比值表格1的异常区R=R1∪R2∪R3∪…∪Rn-2;
3.对于未知来波α,对表格1的所有波束进行波束合成,找出所有输出的最大值和次大值,得到二者比值ratio1',这一比值亦是前一波束输出除以后一波束输出的结果,判断ratio1'是否落入阈值区R内,若没有则直接查找表格1即可得出未知来波方向的角度值α;若落入阈值区R内,则重新对表格2的所有波束进行波束合成,找出所有输出的最大值和次大值ratio2',查找表格2即可得出未知来波方向的角度值α。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710819351.9A CN107656238B (zh) | 2017-09-12 | 2017-09-12 | 一种基于双查找表方法的高频地波雷达测向新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710819351.9A CN107656238B (zh) | 2017-09-12 | 2017-09-12 | 一种基于双查找表方法的高频地波雷达测向新方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107656238A CN107656238A (zh) | 2018-02-02 |
CN107656238B true CN107656238B (zh) | 2020-10-30 |
Family
ID=61129512
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710819351.9A Expired - Fee Related CN107656238B (zh) | 2017-09-12 | 2017-09-12 | 一种基于双查找表方法的高频地波雷达测向新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107656238B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507638A (zh) * | 2018-11-30 | 2019-03-22 | 南京长峰航天电子科技有限公司 | 一种比幅测向方法及*** |
CN110794362B (zh) * | 2019-09-30 | 2022-04-12 | 西安空间无线电技术研究所 | 一种短脉冲高功率微波快速测向***和方法 |
CN111220954B (zh) * | 2019-12-05 | 2022-07-22 | 上海无线电设备研究所 | 一种基于自校幅度归一化的雷达角误差修正方法 |
CN111142061B (zh) * | 2019-12-24 | 2022-06-14 | 泰凌微电子(上海)股份有限公司 | 一种信号到达角估计方法、装置及计算机可读存储介质 |
CN112782645B (zh) * | 2020-12-30 | 2021-11-19 | 无锡国芯微电子***有限公司 | 一种四臂螺旋天线数据拟合测角方法 |
JP2022165747A (ja) | 2021-04-20 | 2022-11-01 | 三菱重工業株式会社 | 方向探知装置、強度差テーブルの取得方法、方向探知方法及び方向探知プログラム |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102171721A (zh) * | 2008-08-26 | 2011-08-31 | 高通股份有限公司 | 用于边缘定向图像缩放的双查找表设计 |
CN103048645A (zh) * | 2013-01-24 | 2013-04-17 | 南京莱斯信息技术股份有限公司 | 雷达回波尾迹的显示***及显示方法 |
CN103714175A (zh) * | 2014-01-07 | 2014-04-09 | 国家卫星海洋应用中心 | 一种低分辨率卫星数据产品空间查询方法及装置 |
CN104122527A (zh) * | 2014-07-14 | 2014-10-29 | 中国人民解放军国防科学技术大学 | 一种基于查表法的圆阵相位干涉仪宽带瞬时测向方法 |
CN106066473A (zh) * | 2016-05-23 | 2016-11-02 | 西安电子科技大学 | 正交波形下mimo雷达目标的多波束联合测角和点迹融合方法 |
CN106772305A (zh) * | 2017-01-23 | 2017-05-31 | 西安电子科技大学 | 一种非正交波形下集中式mimo雷达的目标点迹融合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101675112B1 (ko) * | 2010-01-21 | 2016-11-22 | 삼성전자주식회사 | 거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치 |
-
2017
- 2017-09-12 CN CN201710819351.9A patent/CN107656238B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102171721A (zh) * | 2008-08-26 | 2011-08-31 | 高通股份有限公司 | 用于边缘定向图像缩放的双查找表设计 |
CN103048645A (zh) * | 2013-01-24 | 2013-04-17 | 南京莱斯信息技术股份有限公司 | 雷达回波尾迹的显示***及显示方法 |
CN103714175A (zh) * | 2014-01-07 | 2014-04-09 | 国家卫星海洋应用中心 | 一种低分辨率卫星数据产品空间查询方法及装置 |
CN104122527A (zh) * | 2014-07-14 | 2014-10-29 | 中国人民解放军国防科学技术大学 | 一种基于查表法的圆阵相位干涉仪宽带瞬时测向方法 |
CN106066473A (zh) * | 2016-05-23 | 2016-11-02 | 西安电子科技大学 | 正交波形下mimo雷达目标的多波束联合测角和点迹融合方法 |
CN106772305A (zh) * | 2017-01-23 | 2017-05-31 | 西安电子科技大学 | 一种非正交波形下集中式mimo雷达的目标点迹融合方法 |
Non-Patent Citations (2)
Title |
---|
ship detection in strong clutter environment based on adaptive regression;Ji Yonggang et al.;《2014 International Conference on Computer Science and Electronic Technology (ICCSET 2014)》;20160504;全文 * |
一种阵列式高频地波雷达比幅测向新方法;张玲 等;《中国海洋大学学报》;20170228;第47卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107656238A (zh) | 2018-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107656238B (zh) | 一种基于双查找表方法的高频地波雷达测向新方法 | |
US7369833B2 (en) | Method and apparatus for improving the directivity of an antenna | |
US8354960B2 (en) | Method for low sidelobe operation of a phased array antenna having failed antenna elements | |
CN109946664B (zh) | 一种主瓣干扰下的阵列雷达导引头单脉冲测角方法 | |
CN110673086A (zh) | 一种基于数字阵列雷达的二维角度超分辨方法 | |
CN110197112B (zh) | 一种基于协方差修正的波束域Root-MUSIC方法 | |
CN105137409B (zh) | 基于幅相约束的目标信号稳健空时自适应处理方法 | |
CN110967673B (zh) | 一种多域联合抗主瓣干扰方法 | |
CN112596022B (zh) | 低轨星载多波束正六边形相控阵天线的波达角估计方法 | |
CN114236465B (zh) | 一种搜索式二维高精度比幅测向方法及装置 | |
Dropkin et al. | Superresolution for scanning antenna | |
CN107505592B (zh) | 一种基于多波束雷达粗测向的通信接入方法 | |
CN112698267A (zh) | 一种干扰源测试定位方法 | |
CN115825944B (zh) | 基于外辐射源雷达的单快拍多目标来波方向估计方法 | |
CN109541557B (zh) | 一种雷达天线保护通道生成方法 | |
CN112363108A (zh) | 信号子空间加权超分辨的波达方向检测方法及*** | |
RU2304789C1 (ru) | Способ радиолокационного сопровождения траектории объекта | |
US7372402B2 (en) | Method for enhancing the measuring accuracy in an antenna array | |
CN103792525A (zh) | 一种分布式宽带相控阵雷达阵列基线长度及带宽优化方法 | |
CN110018466B (zh) | 一种基于分子阵预处理的music波束形成方法 | |
Kohler et al. | Wideband direction of arrival estimation using monopulse in Rotman lens beamspace | |
US11994603B1 (en) | Method of sparse array oriented approach for DOA estimation of grating lobe target filtering | |
RU2306580C1 (ru) | Способ измерения угловых координат объекта в процессе захвата и сопровождения траектории в стробах | |
CN110456312B (zh) | 一种基于弧形等相位面的波束展宽方法 | |
Raia et al. | Performance of a hybrid non-contiguous planar array: effect of phase shifter quantization |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201030 |