CN114460539A - 基于相位差分干涉处理的被动合成孔径辐射源定位方法 - Google Patents
基于相位差分干涉处理的被动合成孔径辐射源定位方法 Download PDFInfo
- Publication number
- CN114460539A CN114460539A CN202210135040.1A CN202210135040A CN114460539A CN 114460539 A CN114460539 A CN 114460539A CN 202210135040 A CN202210135040 A CN 202210135040A CN 114460539 A CN114460539 A CN 114460539A
- Authority
- CN
- China
- Prior art keywords
- distance
- signal
- radiation source
- positioning
- phase difference
- 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
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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/0246—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves involving frequency difference of arrival or Doppler measurements
-
- 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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
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
本发明公开了基于相位差分干涉处理的被动合成孔径辐射源定位方法,利用干涉仪天线(两副并排放置的接收天线)接收辐射源信号,将两路接收信号进行相位差分,低通滤波后提取差分信号相位进行线性拟合,记录零相位所对应的时刻,结合定位平台飞行速度获得方位向距离;提取干涉仪中一副天线的接收信号零相位时刻所对频率,即为残余频偏。将单个天线的接收信号去调制,利用不同调频率的匹配滤波器搜索,估计多普勒调频率,获得距离向距离,结合方位向距离,实现最终定位。本发明能够精确估计出残余频偏与方位向距离,且能够将被动合成孔径辐射源定位由二维搜索减至一维搜索,大大降低运算复杂度。
Description
技术领域
本发明涉及一种基于相位差分干涉处理的被动合成孔径辐射源定位方法,属于被动辐射源侦察定位领域。
背景技术
经典被动合成孔径辐射源定位方法通过设置不同调频率与初始频偏的匹配滤波器进行相关运算完成二维搜索,提取最大相关值所对应的调频率与频偏解算距离向与方位向距离实现辐射源定位。由于多普勒信号近似为线性调频信号,模糊函数为山脊形,对频偏不敏感。因此在定位过程中,接收信号残余频偏难以准确估计,导致方位向定位误差较大,限制了被动合成孔径辐射源定位系统的实际应用性能。
发明内容
为解决经典被动合成孔径辐射源定位方法在实际应用过程中所面临的接收信号残余频偏难以准确估计,方位向定位误差过大的问题,本发明公开了基于相位差分干涉处理的被动合成孔径辐射源定位方法。利用干涉仪天线与相位差分运算消除接收信号中的残余频偏,提取相位差分信号的0相位时刻及对应多普勒信号频率实现残余频偏的精确估计与方位向高精度定位。
本发明的技术方案具体包括如下步骤:
基于相位差分干涉处理的被动合成孔径辐射源定位方法,包括以下步骤:
步骤一:接收天线为两规格相同、平行于飞行轨迹放置的干涉仪天线;
步骤二:定位平台以v匀速直线飞越目标辐射源上空,两接收天线均为正侧视,持续接收信号,分别记为x1(t)与x2(t):
步骤三:计算两通道接收信号x1(t)与x2(t)间的相位差值:
步骤四:对相位差分运算后的信号x(t)低通滤波、计算相位、进行线性拟合,提取0相位时刻tp,获得方位向距离Ra的精确估计值;
Ra=vtp (4)
步骤五:提取单通道接收信号xs(t)tp时刻的信号频率f1,为残余频偏的精确估计值;
步骤六:提取[tp-Ts/2,tp+Ts/2]内定位平台的三维坐标(px(t),py(t),pz(t)),Ts为合成孔径时间,px(t)、py(t)、pz(t)分别为x、y、z轴的定位平台瞬时位置坐标;设置距离向距离搜索范围r1,r2,r3,…,rL,并生成距离向距离搜索匹配滤波器mi(t),i=1...L,L为匹配滤波器的个数;
步骤七:对xs(t)进行去调制预处理,截取[tp-Ts/2,tp+Ts/2]时段内的信号xa(t),与本地匹配滤波器m1(t),m2(t),m3(t),…,mL(t)相乘后进行求和运算;
Mi=∫xa(t)×mi(t),i=1,2,3,…,L (5)
步骤八:搜索匹配滤波运算结果M1,M1,…,ML中的最大值Mb,其对应距离rb即为距离向距离估计值R0,结合前文所得方位向距离估计值Ra实现目标辐射源定位。
其中,所述步骤五中单通道接收信号xs(t)的获取方法为从接收信号x1(t)与x2(t)中任选一信号作为xs(t):
xs(t)=x1(t)或xs(t)=x2(t) (6)。
进一步,所述步骤六中距离向距离搜索匹配滤波器m1(t),m2(t),m3(t),…,mL(t)的生成方法为利用距离向距离搜索值r1,r2,r3,…,rL与定位平台的三维坐标(px(t),py(t),pz(t)),计算过顶时段[tp-Ts/2,tp+Ts/2]内的斜距,除以光速获得时延,依据几何关系与残余频偏f1生成本地距离搜索匹配滤波器组:
其中,fr为接收机下变频频率;h为无人机飞行高度。
本发明公开了基于相位差分干涉处理的被动合成孔径辐射源定位方法,用于接收信号残余频偏的精确估计与方位向高精度定位。本发明解决了现有被动合成孔径辐射源定位方法难以精准估计接收信号残余频偏从而导致定位失准的问题。
本发明利用干涉仪天线(两副并排放置的接收天线)接收辐射源信号,将两路接收信号进行相位差分,低通滤波后提取差分信号相位进行线性拟合,记录零相位所对应的时刻,结合定位平台飞行速度获得方位向距离;提取干涉仪中一副天线的接收信号零相位时刻所对频率,即为残余频偏。将单个天线的接收信号去调制,利用不同调频率的匹配滤波器搜索,估计多普勒调频率,获得距离向距离,结合方位向距离,实现最终定位。与传统被动合成孔径辐射源定位方法相比,基于相位差分干涉处理的被动合成孔径辐射源定位方法能够精确估计出残余频偏与方位向距离,且能够将被动合成孔径辐射源定位由二维搜索减至一维搜索,大大降低运算复杂度。本发明可应用于被动辐射源定位领域。
附图说明
图1为现有技术残余频偏对方位向定位结果影响示意图;
图2为现有技术线性调频信号模糊函数示意图;
图3为本发明基于相位差分干涉处理的被动合成孔径辐射源定位方法处理流程示意图;
图4为实施例定位过程简化几何模型示意图;
图5为实施例干涉仪接收天线结构示意图;
图6为实施例相位差分结果幅度/相位曲线示意图;
图7为实施例噪声环境下,相位差分结果线性拟合示意图;
图8为实施例距离向距离搜索示意图。
具体实施方式
下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
参照图1,在传统被动合成孔径通信辐射源定位方法中,残余频偏会影响方位向定位结果,即方位向匹配时刻会随残余频偏变化而变化。载波多普勒信号为线性调频信号,模糊函数为山脊形(如图2),对残余不敏感,导致接收信号中的残余频偏难以精确估计,进而影响方位向定位精度。
针对该问题,本发明提供了一种残余频偏与方位向距离精确估计的方法及系统,其流程示意图如图3所示,该方法包括以下步骤:
步骤一:为验证方法的可行性,选择无人机作为定位平台,定位过程的简化几何模型如图4所示。地面目标辐射源为通信辐射源,载频f0为5GHz,调制方式为BPSK。无人机所携带的接收机采样率为2MHz,接收天线为干涉仪天线(如图5所示),即两规格相同、平行于飞行轨迹的正侧视接收天线,间距l为0.4m,景中心斜距即距离向距离R0为10km,合成孔径时间为1s。
步骤二:无人机定位平台以v=15m/s匀速直线飞越目标辐射源上空,两通道的接收信号x1(t)与x2(t)可以表示为:
步骤三:计算两通道接收信号x1(t)与x2(t)间的相位差值,如公式(3)所示。
步骤四:相位差分信号x(t)近似为一个正弦信号。本场景中该信号频率近似为0.01Hz,幅度/相位曲线如图6所示。定位平台过顶时,辐射源到两接收天线距离相同,差分信号x(t)相位为0。可以此为依据,对相位差分信号x(t)低通滤波、计算相位、进行线性拟合(如图7所示),提取0相位时刻tp估计方位向距离Ra。
Ra=vtp (4)
步骤五:相位差分信号x(t)的0相位时刻tp对应的载波多普勒信号频率为0Hz,此时提取两接收通道中的一路接收信号(x1(t)或x2(t))的频率f1,即为残余频偏的精确估计值。
步骤六:提取[tp-0.5,tp+0.5]间无人机定位平台的三维坐标(px(t),py(t),pz(t)),px(t)、py(t)、pz(t)分别为x、y、z轴的定位平台瞬时位置坐标。并生成距离向距离搜索匹配滤波器mi(t),i=1...L,L为匹配滤波器的个数;
其中,fr为接收机下变频频率;h为无人机飞行高度。
步骤七:提取两接收通道中的一路接收信号(x1(t)或x2(t))的载波多普勒信号进行匹配滤波处理搜索距离向距离。以目标辐射源信号为BPSK调制信号为例,利用平方谱去除数据调制,截取[tp-0.5,tp+0.5]时段内的信号xa(t)(xa(t)=x1 2(t)或x2 2(t),t∈[tp-0.5tp+0.5]),与本地匹配滤波器m1(t),m2(t),m3(t),…,mL(t)相乘后进行求和运算。
Mi=∫xa(t)×mi(t),i=1,2,3,…,L (6)
步骤八:搜索匹配滤波运算结果M1,M1,…,ML中的最大值Mb,其对应距离rb即为距离向距离估计值R0,结合前文所得方位向距离估计值Ra实现目标辐射源定位(距离向距离搜索如图7所示)。
Claims (3)
1.基于相位差分干涉处理的被动合成孔径辐射源定位方法,其特征在于,包括以下步骤:
步骤一:接收天线为两规格相同、平行于飞行轨迹放置的干涉仪天线;
步骤二:定位平台以v匀速直线飞越目标辐射源上空,两接收天线均为正侧视,持续接收信号,分别记为x1(t)与x2(t):
步骤三:计算两通道接收信号x1(t)与x2(t)间的相位差值:
步骤四:对相位差分运算后的信号x(t)低通滤波、计算相位、进行线性拟合,提取0相位时刻tp,获得方位向距离Ra的精确估计值;
Ra=vtp (4)
步骤五:提取单通道接收信号xs(t)tp时刻的信号频率f1,为残余频偏的精确估计值;
步骤六:提取[tp-Ts/2,tp+Ts/2]内定位平台的三维坐标(px(t),py(t),pz(t)),Ts为合成孔径时间,px(t)、py(t)、pz(t)分别为x、y、z轴的定位平台瞬时位置坐标;设置距离向距离搜索范围,并生成距离向距离搜索匹配滤波器mi(t),i=1...L,L为匹配滤波器的个数;
步骤七:对xs(t)进行去调制预处理,截取[tp-Ts/2,tp+Ts/2]时段内的信号xa(t),与本地匹配滤波器mi(t)相乘后进行求和运算;
Mi=∫xa(t)×mi(t),i=1,2,3,…,L (5)
步骤八:搜索匹配滤波运算结果Mi中的最大值Mb,其对应距离rb即为距离向距离估计值R0,结合前文所得方位向距离估计值Ra实现目标辐射源定位。
2.根据权利要求1所述的基于相位差分干涉处理的被动合成孔径辐射源定位方法,其特征在于,所述步骤五中单通道接收信号xs(t)的获取方法为从接收信号x1(t)与x2(t)中任选一信号作为xs(t):
xs(t)=x1(t)或xs(t)=x2(t) (6)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210135040.1A CN114460539B (zh) | 2022-02-14 | 2022-02-14 | 基于相位差分干涉处理的被动合成孔径辐射源定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210135040.1A CN114460539B (zh) | 2022-02-14 | 2022-02-14 | 基于相位差分干涉处理的被动合成孔径辐射源定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114460539A true CN114460539A (zh) | 2022-05-10 |
CN114460539B CN114460539B (zh) | 2022-10-21 |
Family
ID=81413195
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210135040.1A Active CN114460539B (zh) | 2022-02-14 | 2022-02-14 | 基于相位差分干涉处理的被动合成孔径辐射源定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114460539B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115877318A (zh) * | 2022-12-07 | 2023-03-31 | 北京理工大学 | 一种基于多孔径交叉定位的辐射源定位方法 |
CN117572401A (zh) * | 2023-11-24 | 2024-02-20 | 北京理工大学 | 一种基于测量子孔径多普勒相位差的无源定位方法 |
CN117741561A (zh) * | 2024-01-17 | 2024-03-22 | 中国科学院空天信息创新研究院 | 一种辐射源目标三维无源定位装置和方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1304006A (zh) * | 1999-11-23 | 2001-07-18 | 吴光夏 | 多通道阀装置 |
US20100085242A1 (en) * | 2008-10-07 | 2010-04-08 | Sungkyunkwan University Foundation For Corporate Collaboration | Method of sensor network localization through reconstruction of radiation pattern |
US20110057839A1 (en) * | 2009-09-10 | 2011-03-10 | Seifer Arnold D | Method and apparatus for passive geolocation using synthetic-aperture processing |
CN104080167A (zh) * | 2014-06-27 | 2014-10-01 | 东南大学 | 基于硬件特性的无线目标被动定位方法 |
CN105099499A (zh) * | 2015-07-15 | 2015-11-25 | 北京理工大学 | 类噪声Chirp扩频信号设计及快速捕获方法 |
CN107037422A (zh) * | 2017-05-11 | 2017-08-11 | 西北大学 | 一种面向多重应用的被动式定位方法 |
US20170332341A1 (en) * | 2013-12-02 | 2017-11-16 | At&T Intellectual Property I, L.P. | Method and apparatus for performing a passive indoor localization of a mobile endpoint device |
CN108089155A (zh) * | 2017-12-28 | 2018-05-29 | 西北工业大学 | 一种深海环境下单水听器声源被动定位方法 |
CN108680941A (zh) * | 2018-04-27 | 2018-10-19 | 东南大学 | 一种基于单站点的飞行器定位系统及方法 |
CN109597073A (zh) * | 2018-12-20 | 2019-04-09 | 中国科学院声学研究所南海研究站 | 一种无人机载微型近距极化干涉合成孔径雷达系统 |
CN111416635A (zh) * | 2020-03-19 | 2020-07-14 | 上海感悟通信科技有限公司 | 一种测角方法和测角设备 |
CN111896955A (zh) * | 2020-08-06 | 2020-11-06 | 武汉大学 | 一种船载sar交轨干涉处理方法 |
CN113534045A (zh) * | 2021-07-02 | 2021-10-22 | 中国船舶重工集团公司第七二四研究所 | 一种基于匀速圆周运动的单站无源定位方法 |
-
2022
- 2022-02-14 CN CN202210135040.1A patent/CN114460539B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1304006A (zh) * | 1999-11-23 | 2001-07-18 | 吴光夏 | 多通道阀装置 |
US20100085242A1 (en) * | 2008-10-07 | 2010-04-08 | Sungkyunkwan University Foundation For Corporate Collaboration | Method of sensor network localization through reconstruction of radiation pattern |
US20110057839A1 (en) * | 2009-09-10 | 2011-03-10 | Seifer Arnold D | Method and apparatus for passive geolocation using synthetic-aperture processing |
US20170332341A1 (en) * | 2013-12-02 | 2017-11-16 | At&T Intellectual Property I, L.P. | Method and apparatus for performing a passive indoor localization of a mobile endpoint device |
CN104080167A (zh) * | 2014-06-27 | 2014-10-01 | 东南大学 | 基于硬件特性的无线目标被动定位方法 |
CN105099499A (zh) * | 2015-07-15 | 2015-11-25 | 北京理工大学 | 类噪声Chirp扩频信号设计及快速捕获方法 |
CN107037422A (zh) * | 2017-05-11 | 2017-08-11 | 西北大学 | 一种面向多重应用的被动式定位方法 |
CN108089155A (zh) * | 2017-12-28 | 2018-05-29 | 西北工业大学 | 一种深海环境下单水听器声源被动定位方法 |
CN108680941A (zh) * | 2018-04-27 | 2018-10-19 | 东南大学 | 一种基于单站点的飞行器定位系统及方法 |
CN109597073A (zh) * | 2018-12-20 | 2019-04-09 | 中国科学院声学研究所南海研究站 | 一种无人机载微型近距极化干涉合成孔径雷达系统 |
CN111416635A (zh) * | 2020-03-19 | 2020-07-14 | 上海感悟通信科技有限公司 | 一种测角方法和测角设备 |
CN111896955A (zh) * | 2020-08-06 | 2020-11-06 | 武汉大学 | 一种船载sar交轨干涉处理方法 |
CN113534045A (zh) * | 2021-07-02 | 2021-10-22 | 中国船舶重工集团公司第七二四研究所 | 一种基于匀速圆周运动的单站无源定位方法 |
Non-Patent Citations (8)
Title |
---|
LITING ZHANG: "Emitter Localization Algorithm Based on Passive Synthetic Aperture", 《IEEE XPLORE》 * |
W.Q. WANG: "Time and phase synchronisation via direct-path signal for bistatic synthetic aperture radar systems", 《IET RADAR SONAR NAVIG》 * |
张莉婷: "基于被动合成孔径的单星无源高精度定位方法", 《航天电子对抗》 * |
朱伟强等: "基于相位差变化率测量的单站定位方法", 《系统工程与电子技术》 * |
王涛等: "星载双天线干涉SAR系统总体技术研究", 《航天器工程》 * |
王贞松等: "星载合成孔径雷达有源相控阵天线研究", 《电子与信息学报》 * |
蒋一等: "基于ADPCA的三通道合成孔径雷达目标定位技术研究", 《舰船电子工程》 * |
黄瑶等: "一种双通道SAR动目标检测和参数估计方法", 《系统工程与电子技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115877318A (zh) * | 2022-12-07 | 2023-03-31 | 北京理工大学 | 一种基于多孔径交叉定位的辐射源定位方法 |
CN115877318B (zh) * | 2022-12-07 | 2024-04-16 | 北京理工大学 | 一种基于多孔径交叉定位的辐射源定位方法 |
CN117572401A (zh) * | 2023-11-24 | 2024-02-20 | 北京理工大学 | 一种基于测量子孔径多普勒相位差的无源定位方法 |
CN117741561A (zh) * | 2024-01-17 | 2024-03-22 | 中国科学院空天信息创新研究院 | 一种辐射源目标三维无源定位装置和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114460539B (zh) | 2022-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114460539B (zh) | 基于相位差分干涉处理的被动合成孔径辐射源定位方法 | |
JP4683918B2 (ja) | フィルタ中心決定のための方法および装置 | |
US8174444B2 (en) | Method of correlating known image data of moving transmitters with measured radio signals | |
CN101017202B (zh) | 一种雷达高度表及采用该表对飞行器位置的测量方法 | |
CA1279391C (en) | Passive ranging method and apparatus | |
CN107271997B (zh) | 机载多通道cssar地面运动目标运动参数估计方法 | |
JP2011149962A (ja) | 機体座標において目標に対する干渉角を決定するための方法および装置 | |
CN111736187B (zh) | 基于被动合成孔径的高精度高灵敏单星gnss定位方法 | |
JP2002526782A (ja) | 干渉計的合成開口レーダー高度計 | |
CN114624649B (zh) | 对残余频偏不敏感的机载被动合成孔径辐射源定位方法 | |
CN111090109B (zh) | 星载频差提取对频差快速变化的补偿方法 | |
US6680691B2 (en) | Methods and apparatus for accurate phase detection | |
CN112666517A (zh) | 一种基于时差测量的小型无人机信号定位系统及方法 | |
CN110161503A (zh) | 基于三维等距圆模型的近距宽域sar高分辨成像方法 | |
JP2005525562A (ja) | レーダ距離アンビギュオウスの解法のための方法および装置 | |
CN104515974B (zh) | 微波着陆机载设备角度、测距数据处理方法 | |
US6639545B1 (en) | Methods and apparatus to determine a target location in body coordinates | |
CN102937709B (zh) | 基于直达波的双基地sar无同步回波聚焦成像方法 | |
CN112595328A (zh) | 视觉辅助稀疏无线电测量的月球导航定位方法 | |
CN108287353B (zh) | 天基无人机通讯卫星定位跟踪方法 | |
CN110736988B (zh) | 双基地pfa运动目标参数估计和成像方法 | |
CN118112567B (zh) | 一种多通道sar运动误差关联成像处理方法 | |
Dang et al. | A fast two-dimensional velocity estimation method for multi-channel UWB SAR | |
CN117233697B (zh) | 一种针对多个窄带信号运动辐射源的分布式阵列协同直接定位方法 | |
Guo et al. | Short baseline direct position determination based on precise synchronization |
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 |