CN109901104A - 一种通过估计时差解干涉仪测向模糊的方法 - Google Patents
一种通过估计时差解干涉仪测向模糊的方法 Download PDFInfo
- Publication number
- CN109901104A CN109901104A CN201910311268.XA CN201910311268A CN109901104A CN 109901104 A CN109901104 A CN 109901104A CN 201910311268 A CN201910311268 A CN 201910311268A CN 109901104 A CN109901104 A CN 109901104A
- Authority
- CN
- China
- Prior art keywords
- signal
- fuzzy
- phase difference
- array element
- estimation
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 7
- 238000003491 array Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000010129 solution processing Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
Landscapes
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明公开了一种通过估计时差解干涉仪测向模糊的方法,包括估计阵元间信号的相位差和时差、计算阵元间解模糊相位差和计算信号到达方向,通过估计时差去除干涉仪测相位差的模糊,从而实现干涉仪测向解模糊。本发明通过两个阵元阵列,构成一条基线,进行了测向模拟仿真实验,在相位差模糊的情况下,通过时差测量实现了干涉仪测向解模糊处理,典型测向精度优于1°。
Description
技术领域
本发明涉及通信技术领域,尤其涉及一种通过估计时差解干涉仪测向模糊的方法。
背景技术
干涉仪测向是无线电测向的一种通用技术,通常需要多条基线实现干涉仪解模糊处理。可以利用多基线,采用参差基线法或虚拟基线法解模糊实现高精度干涉仪测向。
由于相位差测量存在2π角度模糊问题,导致干涉仪在阵元间距大于半波长时出现测向模糊,现有干涉仪测向***需要3个或以上阵元,典型情况4阵元构成多基线干涉仪,长基线或虚拟长基线实现高精度的方向估计,短基线或虚拟短基线实现解模糊处理,才能实现高精度干涉仪测向。
发明内容
为了解决上述问题,本发明通过在干涉仪阵列中额外利用信号达到时差信息,在两阵元间距大于半波长时可实现高精度干涉仪测向的解模糊处理,具体的,本发明提出一种通过估计时差解干涉仪测向模糊的方法,包括以下步骤:
S1.估计阵元间信号的相位差和时差:采集信号样本,通过相关法得到视在相位差估计值和时差估计值其中
S2.计算阵元间解模糊相位差:阵元接收信号相位差公式如下,
其中φ为相位差,D为阵元间距,θ为辐射源信号入射角,λ为波长,f为信号频率,τ为时差;更大的阵元间距D将获得更高的测角精度,但当时,视在相位差值存在角度模糊问题,由时差估计值可得到相位差模糊数估计值:
其中[]表示四舍五入取整,从而得到无模糊的相位差估计值
S3.计算信号到达方向:根据以及无模糊的相位差估计值可得到无模糊的方向估计值:
进一步的,所述步骤S1中,对于脉冲信号,可以利用脉冲到达时间相减得到时差估计值
进一步的,所述步骤S2中,在信号频率f未知的情况下,采集到信号样本后,还需要根据信号样本估计频率。
本发明的有益效果在于:本发明通过两个阵元阵列,构成一条基线,进行了测向模拟仿真实验,在相位差模糊的情况下,通过时差测量实现了干涉仪测向解模糊处理,典型测向精度优于1°。
附图说明
图1是干涉仪构型几何及阵元间信号关系图;
图2是信噪比20dB条件下的含噪信号样本;
图3是信噪比20dB条件下的100次角度估计结果。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图说明本发明的具体实施方式。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明,即所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
本实施例的具体参数如下:
频率f:300MHz(对应的信号波长1m);
信号入射角θ:55°;
阵元间距D:4m;
采样率(复采样)fs:4GSPS。
设入射信号为脉冲串,脉宽0.2μs,脉冲重复周期0.5μs。如图1所示是干涉仪构型几何及阵元间信号关系图。
本实施例的具体步骤如下:
1)采集信号并估计阵元间信号的视在相位差和时差
采集信号样本,取实部如图2所示,信噪比为20dB,注意到阵元1样本较阵元2样本有延迟,直接通过相关可得到角弧度1.6573rad。
直接通过相关法估计时差,得到相关延迟为44个采样样本,延迟为3.3个信号周期。对于脉冲信号,也可以直接用脉冲到达时间相减得到时差估计。
2)计算阵元间的解模糊相位差
信号真实相位差在3.3个周期左右,模糊数为3,考虑相位差估计值,可得到实际相位的角弧度估计值
3)计算信号到达方向
换算为角度值与实际角度相差0.3129°。
信噪比20dB输入,100次测量估计结果如图3所示。
测量结果的统计均值为55.0665°,角度估计均方根误差0.6370°。
需要说明的是,在某些频率未知的情况下,首先需要进行频率估计,然后利用上面的3个步骤得到角度估计值。
以上所述仅是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (3)
1.一种通过估计时差解干涉仪测向模糊的方法,其特征在于,包括以下步骤:
S1.估计阵元间信号的相位差和时差:采集信号样本,通过相关法得到视在相位差估计值和时差估计值其中
S2.计算阵元间解模糊相位差:阵元接收信号相位差公式如下,
其中φ为相位差,D为阵元间距,θ为辐射源信号入射角,λ为波长,f为信号频率,τ为时差;更大的阵元间距D将获得更高的测角精度,但当时,相位差计算值φo∈[-π,+π]存在2π模糊问题,由时差估计值可得到相位差模糊数估计值:
其中[]表示四舍五入取整,从而得到无模糊的相位差估计值
S3.计算信号到达方向:根据以及无模糊的相位差估计值可得到无模糊的方向估计值:
2.根据权利要求1所述的一种通过估计时差解干涉仪测向模糊的方法,其特征在于,所述步骤S1中,对于脉冲信号,可以利用脉冲到达时间相减得到时差估计值
3.根据权利要求1所述的一种通过估计时差解干涉仪测向模糊的方法,其特征在于,所述步骤S2中,在信号频率f未知的情况下,采集到信号样本后,还需要根据信号样本估计频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910311268.XA CN109901104A (zh) | 2019-04-18 | 2019-04-18 | 一种通过估计时差解干涉仪测向模糊的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910311268.XA CN109901104A (zh) | 2019-04-18 | 2019-04-18 | 一种通过估计时差解干涉仪测向模糊的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109901104A true CN109901104A (zh) | 2019-06-18 |
Family
ID=66955059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910311268.XA Pending CN109901104A (zh) | 2019-04-18 | 2019-04-18 | 一种通过估计时差解干涉仪测向模糊的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109901104A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110187303A (zh) * | 2019-06-24 | 2019-08-30 | 中国电子科技集团公司第二十九研究所 | 一种单脉冲比幅比相联合测向方法 |
CN110687498A (zh) * | 2019-11-06 | 2020-01-14 | 天津津航计算技术研究所 | 动基线干涉仪 |
CN111693936A (zh) * | 2020-06-19 | 2020-09-22 | 中国电子科技集团公司第二十九研究所 | 一种基于射频模拟接收***的测向方法及*** |
CN111693934A (zh) * | 2020-06-19 | 2020-09-22 | 中国电子科技集团公司第二十九研究所 | 一种基于数字干涉***的测向方法及*** |
CN112485756A (zh) * | 2020-11-20 | 2021-03-12 | 中国电子科技集团公司第二十九研究所 | 基于dbf阵列天线的干涉仪信号处理方法 |
CN114428225A (zh) * | 2022-04-06 | 2022-05-03 | 中国电子科技集团公司第二十九研究所 | 一种基于准匹配滤波的多辐射源到达角度测量方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102411136A (zh) * | 2011-08-09 | 2012-04-11 | 电子科技大学 | 一种扩展基线解模糊的相位干涉仪测向方法 |
CN102540138A (zh) * | 2011-11-25 | 2012-07-04 | 华中科技大学 | 一种多基线相位搜索式二维空间谱测向方法 |
CN104635201A (zh) * | 2013-11-06 | 2015-05-20 | 郁涛 | 一种基于相位差分判别的无模糊测向方法 |
CN105676174A (zh) * | 2016-01-25 | 2016-06-15 | 西安电子科技大学 | 基于相位差分最小二乘法的单站无源快速定位方法 |
CN106772239A (zh) * | 2017-02-20 | 2017-05-31 | 上海微小卫星工程中心 | 旋转基线干涉仪定位***的模拟源测试***及其使用方法 |
CN108254718A (zh) * | 2017-12-04 | 2018-07-06 | 西安电子科技大学 | 基于旋转干涉仪的多目标角度解模糊方法 |
-
2019
- 2019-04-18 CN CN201910311268.XA patent/CN109901104A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102411136A (zh) * | 2011-08-09 | 2012-04-11 | 电子科技大学 | 一种扩展基线解模糊的相位干涉仪测向方法 |
CN102540138A (zh) * | 2011-11-25 | 2012-07-04 | 华中科技大学 | 一种多基线相位搜索式二维空间谱测向方法 |
CN104635201A (zh) * | 2013-11-06 | 2015-05-20 | 郁涛 | 一种基于相位差分判别的无模糊测向方法 |
CN105676174A (zh) * | 2016-01-25 | 2016-06-15 | 西安电子科技大学 | 基于相位差分最小二乘法的单站无源快速定位方法 |
CN106772239A (zh) * | 2017-02-20 | 2017-05-31 | 上海微小卫星工程中心 | 旋转基线干涉仪定位***的模拟源测试***及其使用方法 |
CN108254718A (zh) * | 2017-12-04 | 2018-07-06 | 西安电子科技大学 | 基于旋转干涉仪的多目标角度解模糊方法 |
Non-Patent Citations (3)
Title |
---|
元亮: "基于时差解相位干涉仪模糊的测向方法研究", 《舰船电子对抗》 * |
狄慧: "联合到达时间估计的长基线测向相位解模糊算法研究", 《电子学报》 * |
王国林: "长基线干涉仪测向关键技术研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110187303A (zh) * | 2019-06-24 | 2019-08-30 | 中国电子科技集团公司第二十九研究所 | 一种单脉冲比幅比相联合测向方法 |
CN110687498A (zh) * | 2019-11-06 | 2020-01-14 | 天津津航计算技术研究所 | 动基线干涉仪 |
CN110687498B (zh) * | 2019-11-06 | 2023-04-28 | 天津津航计算技术研究所 | 动基线干涉仪 |
CN111693936A (zh) * | 2020-06-19 | 2020-09-22 | 中国电子科技集团公司第二十九研究所 | 一种基于射频模拟接收***的测向方法及*** |
CN111693934A (zh) * | 2020-06-19 | 2020-09-22 | 中国电子科技集团公司第二十九研究所 | 一种基于数字干涉***的测向方法及*** |
CN111693934B (zh) * | 2020-06-19 | 2022-08-09 | 中国电子科技集团公司第二十九研究所 | 一种基于数字干涉***的测向方法及*** |
CN112485756A (zh) * | 2020-11-20 | 2021-03-12 | 中国电子科技集团公司第二十九研究所 | 基于dbf阵列天线的干涉仪信号处理方法 |
CN112485756B (zh) * | 2020-11-20 | 2022-07-05 | 中国电子科技集团公司第二十九研究所 | 基于dbf阵列天线的干涉仪信号处理方法 |
CN114428225A (zh) * | 2022-04-06 | 2022-05-03 | 中国电子科技集团公司第二十九研究所 | 一种基于准匹配滤波的多辐射源到达角度测量方法及装置 |
CN114428225B (zh) * | 2022-04-06 | 2022-06-14 | 中国电子科技集团公司第二十九研究所 | 一种基于准匹配滤波的多辐射源到达角度测量方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109901104A (zh) | 一种通过估计时差解干涉仪测向模糊的方法 | |
JP6778336B2 (ja) | Rf信号の到来角度決定方法及びシステム | |
CN106231670B (zh) | 一种无接触感知定位方法 | |
CN102411136B (zh) | 一种扩展基线解模糊的相位干涉仪测向方法 | |
CN105759241B (zh) | 基于时差频差的直接定位方法 | |
Naseri et al. | Cooperative simultaneous localization and mapping by exploiting multipath propagation | |
CN112073895A (zh) | 一种基于csi高精度人员定位跟踪方法 | |
Liu et al. | Azimuth and elevation estimation with rotating long-baseline interferometers | |
US11762085B2 (en) | Device, system and method for localization of a target in a scene | |
CN108398670B (zh) | 一种基于旋转干涉仪的脉冲信号测向方法及装置 | |
CN105044667A (zh) | 一种运动目标的双星跟踪方法、装置和*** | |
Keunecke et al. | Deriving 2D TOA/TDOA IEEE 802.11 g/n/ac location accuracy from an experimentally verified fading channel model | |
CN106569180B (zh) | 一种基于Prony方法的方位估计算法 | |
RU2711341C1 (ru) | Способ двухмерного пеленгования | |
CN111741429B (zh) | 基于信号分布Wasserstein距离度量的Wi-Fi室内定位方法 | |
KR101762510B1 (ko) | 위치 추정 정확도 향상을 위한 이동형 기준 노드 위치 결정 방법 | |
Qin et al. | A ML-based direct localization method for multiple sources with moving arrays | |
Sanchez et al. | From theory to practice: Implementations issues for techniques based on AoA and TDoA | |
Searle | An examination of bias in SODA interferometry | |
Webb et al. | Radar imaging of three-dimensional surfaces using limited data | |
Witzgall | Two-sensor tracking of maneuvering transmitters | |
CN108226850A (zh) | 一种基于抛物线拟合的单脉冲鉴相方法 | |
Kaune | Performance analysis of passive emitter tracking using TDOA, AOA and FDOA measurements | |
Liu et al. | Displacement Measurement of Weak Targets with Imaging Radar | |
Sippel | Holographic 3D Indoor Localization |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190618 |