CN109581278A - 一种虚拟孔径扩展的相关干涉仪测向方法 - Google Patents

一种虚拟孔径扩展的相关干涉仪测向方法 Download PDF

Info

Publication number
CN109581278A
CN109581278A CN201811462154.7A CN201811462154A CN109581278A CN 109581278 A CN109581278 A CN 109581278A CN 201811462154 A CN201811462154 A CN 201811462154A CN 109581278 A CN109581278 A CN 109581278A
Authority
CN
China
Prior art keywords
finding
vector
database
collected
extension
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
CN201811462154.7A
Other languages
English (en)
Other versions
CN109581278B (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.)
China Research Institute of Radio Wave Propagation CRIRP
China Institute of Radio Wave Propagation CETC 22 Research Institute
Original Assignee
China Research Institute of Radio Wave Propagation CRIRP
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 China Research Institute of Radio Wave Propagation CRIRP filed Critical China Research Institute of Radio Wave Propagation CRIRP
Priority to CN201811462154.7A priority Critical patent/CN109581278B/zh
Publication of CN109581278A publication Critical patent/CN109581278A/zh
Application granted granted Critical
Publication of CN109581278B publication Critical patent/CN109581278B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Direction-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/02Direction-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/14Systems for determining direction or deviation from predetermined direction
    • G01S3/143Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Direction-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/02Direction-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/14Systems for determining direction or deviation from predetermined direction
    • G01S3/46Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems

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)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明公开了一种虚拟孔径扩展的相关干涉仪测向方法,包括如下步骤:(1)利用双通道或多通道测向接收机,记录同时采集到的IQ信号;(2)计算参考天线单元与参考天线单元之外的其他各天线单元采集到的IQ信号之间的二阶统计量复值相位差、四阶统计量复值相位差,并组合为测量矢量;(3)读取相应频点的测向数据库,计算扩展测向数据库;(4)将(2)中的测量矢量与(3)中的扩展测向数据库进行相关计算,相关值最大处作为方向的估计值。本发明所公开虚拟孔径扩展的相关干涉仪测向方法,基于测向接收机采集到的IQ信号的二阶和四阶统计信息,实现了测向天线的虚拟孔径扩展,改善了相关干涉仪测向设备在频率低端的测向性能。

Description

一种虚拟孔径扩展的相关干涉仪测向方法
技术领域
本发明属于频谱管理中的测向技术领域,特别涉及该领域中的一种虚拟孔径扩展的相关干涉仪测向方法。
背景技术
对辐射源进行测向是频谱管理领域的常见任务之一,通常由固定监测站或者移动监测站(主要为车载式)中的测向天线与接收机来实现。
目前测向技术体制主要是沃森瓦特方法和相关干涉仪方法,其中沃森瓦特方法主要适用于HF频段和VHF频段的低端,其测向精度受到天线单元的方向图的不圆度的限制,而且对于电波反射效应、设备误差等因素敏感;相关干涉仪方法适用于HF及以上频段,其对电波反射效应、设备误差等因素的鲁棒性较高,测向精度主要受到阵列孔径的限制。
在移动监测站中,由于天线体积和数量受到了平台的各种限制,VHF/UHF频段通常由一副测向天线覆盖整个频段。如果使用相关干涉仪方法,对于频率中高端而言,天线孔径较为合适,而对于频率低端而言,天线孔径偏小,相应的测向性能因此受到影响。理论上可以考虑引入沃森瓦特方法提高频率低端的测向性能,但是,因为移动平台对频率低端的电波具有明显的反射效应,所以沃森瓦特方法实际效果不佳。也就是说,如何提升频率低端的测向性能成为亟待解决的问题。
发明内容
本发明所要解决的技术问题就是提供一种可以解决现有相关干涉仪测向技术中频率低端测向性能受限问题的虚拟孔径扩展的相关干涉仪测向方法。
本发明采用如下技术方案:
一种虚拟孔径扩展的相关干涉仪测向方法,其改进之处在于,包括如下步骤:
(1)利用双通道或多通道测向接收机,记录同时采集到的IQ信号;
(2)计算参考天线单元与参考天线单元之外的其他各天线单元采集到的IQ信号之间的二阶统计量复值相位差、四阶统计量复值相位差,并组合为测量矢量;
具体为:将(1)记录的参考天线单元采集到的IQ信号记为s0(t),与其同时采集到的第n个天线单元采集到的IQ信号记为sn(t),n=1,2,…,N-1,其中,N为天线单元数量;
二阶统计量相位差记为计算公式为:
其中,E{·}代表数学期望,|·|代表绝对值,(·)*代表复数共轭,j代表虚部单位;
四阶统计量相位差记为计算公式为:
其中,δ代表预设的门限值;
测量矢量记为v,计算公式为:
测量矢量的长度为2N-2;
(3)读取相应频点的测向数据库,计算扩展测向数据库;
具体为:扩展测向数据库中的每一个矢量的长度为测向数据库中的每一个矢量的长度的两倍;
扩展测向数据库中的每一个矢量的前二分之一的元素与测向数据库中的每一个矢量的元素相同;
扩展测向数据库中的每一个矢量的后二分之一的元素为测向数据库中的每一个矢量的元素的平方;
(4)将(2)中的测量矢量与(3)中的扩展测向数据库进行相关计算,相关值最大处作为方向的估计值。
本发明的有益效果是:
本发明所公开虚拟孔径扩展的相关干涉仪测向方法,基于测向接收机采集到的IQ信号的二阶和四阶统计信息,实现了测向天线的虚拟孔径扩展,提供了一种相关干涉仪测向方法,改善了相关干涉仪测向设备在频率低端的测向性能。
本发明所公开虚拟孔径扩展的相关干涉仪测向方法,不需要对现有的相关干涉仪测向设备重新进行厂内校准、建立测向数据库的过程,仅需对测向设备的算法内核进行更新,易于被推广使用,具有较好的实用价值。
附图说明
图1是本发明实施例1所公开测向方法的流程示意图;
图2是本发明实施例1所公开测向方法步骤(1)的流程示意图;
图3是本发明实施例1所公开测向方法步骤(2)的流程示意图;
图4是本发明实施例1所公开测向方法步骤(3)的流程示意图;
图5是本发明实施例1所公开测向方法步骤(4)的仿真图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例1,假设某测向设备的基本配置如下:
假设测向天线有九个天线单元(编号为0,1,……,8),其中第0个单元为参考天线单元;
假设接收机有两个接收通道,其中一个接通参考天线单元,另外通过天线选择开关的形式接通其他天线单元;
假设连续采样点数为4096;
如图1所示,本实施例所公开的虚拟孔径扩展的相关干涉仪测向方法,包括如下步骤:
步骤一、利用双通道测向接收机,记录同时采集到的IQ信号;
步骤二、计算参考天线单元与参考天线单元之外的其他各天线单元接收到的IQ信号之间的二阶统计量复值相位差、四阶统计量复值相位差,并组合为测量矢量;
步骤三、读取相应频点的测向数据库,计算扩展测向数据库;
步骤四、将步骤二中的测量矢量与步骤三中的扩展测向数据库进行相关计算,相关值最大处作为方向的估计值;
如图2所示,步骤一通过依次接入天线单元的形式,记录8组同时采集到的IQ信号(共计65536个IQ点),即
如图3所示,步骤二具体包括:计算参考天线单元与参考天线单元之外的其他各天线单元接收到的IQ信号之间的二阶统计量复值相位差、四阶统计量复值相位差,并组成矢量,即:
二阶统计量相位差记为计算公式为:
其中,E{·}代表数学期望,|·|代表绝对值,(·)*代表复数共轭,j代表虚部单位;
四阶统计量复值相位差记为计算公式为:
其中,δ代表预设的门限值;
测量矢量记为v,组合方式为:
其中,(·)T代表矩阵转置,测量矢量的长度为16;
如图4所示,步骤三具体包括:读取相应频点的测向数据库,计算扩展测向数据库,其中:
扩展测向数据库中的每一个矢量的长度为16;
扩展测向数据库中的每一个矢量的前8个元素与测向数据库中的每一个矢量的元素相同;
扩展测向数据库中的每一个矢量的后8个元素为测向数据库中的每一个矢量的元素的平方;
如图5所示,步骤四将步骤二中的测量矢量与步骤三中的扩展测向数据库进行相关计算,相关值最大处作为方向的估计值,其中信号中心频率为50MHz,位于VHF/UHF频段的频率低端,入射角度为100度。

Claims (1)

1.一种虚拟孔径扩展的相关干涉仪测向方法,其特征在于,包括如下步骤:
(1)利用双通道或多通道测向接收机,记录同时采集到的IQ信号;
(2)计算参考天线单元与参考天线单元之外的其他各天线单元采集到的IQ信号之间的二阶统计量复值相位差、四阶统计量复值相位差,并组合为测量矢量;
具体为:将(1)记录的参考天线单元采集到的IQ信号记为s0(t),与其同时采集到的第n个天线单元采集到的IQ信号记为sn(t),n=1,2,…,N-1,其中,N为天线单元数量;
二阶统计量相位差记为计算公式为:
其中,E{·}代表数学期望,|·|代表绝对值,(·)*代表复数共轭,j代表虚部单位;
四阶统计量相位差记为计算公式为:
其中,δ代表预设的门限值;
测量矢量记为v,计算公式为:
测量矢量的长度为2N-2;
(3)读取相应频点的测向数据库,计算扩展测向数据库;
具体为:扩展测向数据库中的每一个矢量的长度为测向数据库中的每一个矢量的长度的两倍;
扩展测向数据库中的每一个矢量的前二分之一的元素与测向数据库中的每一个矢量的元素相同;
扩展测向数据库中的每一个矢量的后二分之一的元素为测向数据库中的每一个矢量的元素的平方;
(4)将(2)中的测量矢量与(3)中的扩展测向数据库进行相关计算,相关值最大处作为方向的估计值。
CN201811462154.7A 2018-12-03 2018-12-03 一种虚拟孔径扩展的相关干涉仪测向方法 Active CN109581278B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811462154.7A CN109581278B (zh) 2018-12-03 2018-12-03 一种虚拟孔径扩展的相关干涉仪测向方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811462154.7A CN109581278B (zh) 2018-12-03 2018-12-03 一种虚拟孔径扩展的相关干涉仪测向方法

Publications (2)

Publication Number Publication Date
CN109581278A true CN109581278A (zh) 2019-04-05
CN109581278B CN109581278B (zh) 2022-03-18

Family

ID=65926623

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811462154.7A Active CN109581278B (zh) 2018-12-03 2018-12-03 一种虚拟孔径扩展的相关干涉仪测向方法

Country Status (1)

Country Link
CN (1) CN109581278B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114563756A (zh) * 2022-01-12 2022-05-31 中国电子科技集团公司第十研究所 双通道相关干涉仪测向样本线性插值方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070273576A1 (en) * 2006-05-27 2007-11-29 Bae Systems Information And Electronic Systems Integration, Inc. Multipath resolving correlation interferometer direction finding
CN102411136A (zh) * 2011-08-09 2012-04-11 电子科技大学 一种扩展基线解模糊的相位干涉仪测向方法
CN102419432A (zh) * 2011-08-25 2012-04-18 电子科技大学 一种基于虚拟基线的圆阵相位干涉仪二维测向方法
WO2012134456A1 (en) * 2011-03-29 2012-10-04 Son Nguyen Directional radio-signal-detection apparatus and methods of use
CN102749612A (zh) * 2012-07-25 2012-10-24 中国人民解放军信息工程大学 一种相关干涉仪的快速测向方法
JP2017036990A (ja) * 2015-08-10 2017-02-16 三菱電機株式会社 到来方向推定装置
CN107422295A (zh) * 2017-08-30 2017-12-01 浙江大学 基于互质阵列虚拟域等价信号原子范数表示的波达方向估计方法
CN107490780A (zh) * 2017-06-01 2017-12-19 同方电子科技有限公司 一种可抑制均匀分布的相位误差的测向方法
CN108051772A (zh) * 2017-10-25 2018-05-18 西南电子技术研究所(中国电子科技集团公司第十研究所) 幅相联合测量来波方位信息的方法
CN108254719A (zh) * 2017-12-19 2018-07-06 杭州电子科技大学 一种适用于电磁频谱监测系统的信号源定位方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070273576A1 (en) * 2006-05-27 2007-11-29 Bae Systems Information And Electronic Systems Integration, Inc. Multipath resolving correlation interferometer direction finding
WO2012134456A1 (en) * 2011-03-29 2012-10-04 Son Nguyen Directional radio-signal-detection apparatus and methods of use
CN102411136A (zh) * 2011-08-09 2012-04-11 电子科技大学 一种扩展基线解模糊的相位干涉仪测向方法
CN102419432A (zh) * 2011-08-25 2012-04-18 电子科技大学 一种基于虚拟基线的圆阵相位干涉仪二维测向方法
CN102749612A (zh) * 2012-07-25 2012-10-24 中国人民解放军信息工程大学 一种相关干涉仪的快速测向方法
JP2017036990A (ja) * 2015-08-10 2017-02-16 三菱電機株式会社 到来方向推定装置
CN107490780A (zh) * 2017-06-01 2017-12-19 同方电子科技有限公司 一种可抑制均匀分布的相位误差的测向方法
CN107422295A (zh) * 2017-08-30 2017-12-01 浙江大学 基于互质阵列虚拟域等价信号原子范数表示的波达方向估计方法
CN108051772A (zh) * 2017-10-25 2018-05-18 西南电子技术研究所(中国电子科技集团公司第十研究所) 幅相联合测量来波方位信息的方法
CN108254719A (zh) * 2017-12-19 2018-07-06 杭州电子科技大学 一种适用于电磁频谱监测系统的信号源定位方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
HW WEI 等: "Performance analysis and comparison of correlative interferometers for direction finding", 《IEEE 10TH INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING PROCEEDINGS》 *
JUNG-HOON LEE 等: "Interferometer direction-finding system with improved DF accuracy using two different array configurations", 《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》 *
XIAOMING GOU 等: "Biquaternion noncircular MUSIC", 《MULTIDIMENSIONAL SYSTEMS AND SIGNAL PROCESSING》 *
刘剑 等: "利用高阶信息的干涉仪测向算法", 《航空兵器》 *
吴奉微 等: "基于虚拟阵列变换的干涉仪测向算法", 《现代雷达》 *
张小飞 等: "基于多普勒频移及变化率的地面单载波干扰源定位新方法研究及仿真", 《第22届全国电磁兼容学术会议论文选》 *
苟晓鸣 等: "基于压缩感知的相关干涉仪测向算法", 《电波科学学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114563756A (zh) * 2022-01-12 2022-05-31 中国电子科技集团公司第十研究所 双通道相关干涉仪测向样本线性插值方法

Also Published As

Publication number Publication date
CN109581278B (zh) 2022-03-18

Similar Documents

Publication Publication Date Title
CN111044971B (zh) 一种二维干涉仪测向系统
KR100416026B1 (ko) 무선 통신 시스템에서 측정된 무선 주파수 신호 전파특성을 비교하는 방법 및 시스템
KR100441068B1 (ko) 무선 통신 시스템의 시그너춰 구역에 대해 무선 주파수시그너춰를 생성하는 방법 및 시스템
CN106324559B (zh) 一种大基线四元阵宽带信号测向系统及方法
CN110351658B (zh) 一种基于卷积神经网络的室内定位方法
CN104515909B (zh) 一种基于相关法的大天线方向图测量方法
CN103630878B (zh) 一种基于迭代自适应方法的频带融合方法及装置
CN101487888B (zh) 一种空间谱谱峰搜索方法
CN111381226B (zh) 一种基于多频段接收信号强度的测距增强方法
CN111413675A (zh) 一种信号采集设备及其采集方法
CN106950528A (zh) 一种基于线性调频信号的波达方向估计方法
CN104144495A (zh) 一种基于方向传感器与wlan网络的指纹定位方法
CN111352083B (zh) 一种高频地波雷达多接收通道增益自动校准方法及装置
CN103338024A (zh) 天线组阵中时延的互补卡尔曼滤波装置与方法
CN106569180B (zh) 一种基于Prony方法的方位估计算法
CN109581278A (zh) 一种虚拟孔径扩展的相关干涉仪测向方法
CN102540153A (zh) 一种基于站间直达波干扰的阵列幅相误差校准方法
CN106707280B (zh) 基于配准和曲线拟合的合成孔径雷达有效基线估计方法
CN113960523A (zh) 基于fpga的通用超宽带校正测向方法及系统
CN112880633A (zh) 基于伯格算法的海面高度测量方法
Comsa et al. Wireless localization using time difference of arrival in narrow-band multipath systems
CN113702901B (zh) 一种基于信号相关性增强的互相关时延估计方法
CN114035149A (zh) 一种提升干涉仪测向系统灵敏度的方法
CN112327052B (zh) 一种快速高精度测频方法和系统
CN100399730C (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