CN101247141A - 一种跳频信号时差估计方法 - Google Patents
一种跳频信号时差估计方法 Download PDFInfo
- Publication number
- CN101247141A CN101247141A CNA200710050704XA CN200710050704A CN101247141A CN 101247141 A CN101247141 A CN 101247141A CN A200710050704X A CNA200710050704X A CN A200710050704XA CN 200710050704 A CN200710050704 A CN 200710050704A CN 101247141 A CN101247141 A CN 101247141A
- Authority
- CN
- China
- Prior art keywords
- frequency hopping
- centerdot
- hopping signal
- time difference
- phi
- 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
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
跳频信号时差估计方法属信号处理范畴。利用两个不同位置处接收的跳频信号间的时差与这两个位置接收的跳频信号载波之间的相位差的线性对应关系,通过测量这两个不同位置处接收的跳频信号载波之间的多个跳频周期内的相位差,估计这两个不同位置处接收的跳频信号之间的时差。包括对两个不同位置处接收的跳频信号进行带通滤波;在一个跳频周期内的Hilbert变换、计算相位;确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差、相位差的中值;在K个跳频周期内的相位差中值序列;解出中值序列对应的线性方程组的最小二乘解确定跳频信号的时差估计等步骤。本方法可对跳频信号的高精度、快速的估计时差,广泛应用于跳频辐射信号源定位系统。
Description
所属技术领域
本发明属于信号处理范畴,涉及时差估计方法,尤其涉及跳频信号的高精度、快速时差估计方法。
背景技术
高精度时差估计在雷达、声纳、卫星和通信等系统的电磁波辐射源定位中具有重要的应用价值。由于不同位置处接收的信号之间的时差与信号包络之间的互相关函数的峰值位置之间存在一一对应的关系,因此可以利用该关系设计时差估计技术,通过对不同位置处接收的信号进行互相关处理,实现信号时差估计。在基于这种原理的现有时差估计技术中,典型的是广义互相关(General Cross Correlation)时差估计技术。由信号处理和参数估计的基本理论可知,这类时差估计的方差的下界与信号带宽的三次方成反比,与信噪比和信号的时间长度成反比。因此,针对窄带信号的时差估计问题,需要较大的信噪比和较长的信号时间长度,才能获得较高精度的窄带信号时差估计。
当信号是跳频信号时,接收的信号在整个信号观测时间内按跳频周期随机的改变跳频频偏。由于具有相同跳频频偏的信号仅占整个信号观测时间的一部分,使得对应该跳频频偏的信号在整个信号观测时间内的平均信息速率远小于每个跳频周期内的信息速率,即对应每个跳频频偏的信号的基带带宽实际上变得更窄了。此时如果仍利用不同位置处接收的信号之间的时差与信号包络之间的互相关函数的峰值位置之间存在的一一对应关系设计跳频信号时差估计方法,时差估计的方差的下界将难以满足实际应用的需要。
若要避免对应每个跳频频偏的信号的基带带宽变得更窄,则只能利用每个跳频周期的接收信号估计信号之间的时差。由于一个跳频周期的时间是有限的,将信号观测时间限制为一个跳频周期也将使时差估计的方差的下界难以满足实际应用的需要。
可见,如能利用较短时间内在两个不同位置处接收的跳频信号得到高精度的时差估计,将在保持高精度跳频信号时差估计性能的同时,有效的提高跳频信号时差估计的速度,从而对跳频辐射信号源定位具有重要意义。
发明内容
本发明的目的是为跳频辐射信号源的精确、快速定位提供一种高精度、快速的时差估计方法。
本发明的目的是这样达到的:
一种跳频信号时差估计方法,其特征在于:利用两个不同位置处接收的跳频信号之间的时差与这两个不同位置处接收的跳频信号载波之间的相位差存在的线性对应关系,通过测量这两个不同位置处接收的跳频信号载波之间的多个跳频周期内的相位差,估计这两个不同位置处接收的跳频信号之间的时差。其具体步骤是:
根据一个跳频周期(T)的载波频率和信号带宽,对两个不同位置处接收的跳频信号进行带通滤波;
确定两个不同位置处接收的跳频信号在一个跳频周期内的Hilbert变换;
利用Hilbert变换结果分别确定两个不同位置处接收的跳频信号在一个跳频周期内的相位;
确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差;
确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差的中值;
依次确定两个不同位置处接收的跳频信号在K(大于1的整数)个跳频周期内的相位差的中值序列;
由K个跳频周期内相位差的中值序列对应的线性方程组的最小二乘解确定跳频信号的时差估计。
所述根据一个跳频周期内的载波频率和带宽,对两个不同位置处接收的跳频信号x1(t+kT)和x2(t+kT)分别进行带通滤波,结果为:
和
其中,-T/2≤t≤T/2,a(m)为根据一个跳频周期内的载波频率和带宽设计的带通滤波器的系数,m=1,2,...,M。
所述确定两个不同位置处接收的跳频信号在一个跳频周期内的Hilbert变换分别为:
和 其中, 是单位虚数。
所述利用Hilbert变换结果分别确定两个不同位置处接收的跳频信号在一个跳频周期内的相位分别为 1(t+kT)=angle(y1(t+kT))
和 2(t+kT)=angle(y2(t+kT))
其中,angle( )表示取复数的相角。
所述确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差为:
φ(t+kT)=2(t+kT)-1(t+kT)
其中,-T/2≤t≤T/2。
所述确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差的中值为:
其中, 表示φ(t+kT)在-T/2≤t≤T/2范围内的中值。
所述依次确定两个不同位置处接收的跳频信号在K个跳频周期的相位差的中值序列为:
所述由K个跳频周期内相位差的中值序列对应的线性方程组:
的解:
本发明的积极效果:为跳频辐射信号源的精确、快速定位提供了一种高精度、快速的时差估计方法。满足了在雷达、声纳、卫星和通信等系统的电磁波辐射源定位中高精度时差估计的需要,为电磁波监测技术的发展和应用拓展了新的发展空间。根据本发明的方法进行的实验表明,将本发明提出的跳频信号时差估计方法应用于每秒250跳,驻留时间3.5毫秒,跳频频偏依次为-25kHz,225kHz,0kHz,-250kHz,250kHz,775kHz的跳频信号,在两个不同位置处接收的跳频信号时差为493纳秒,基带信噪比15dB的情况下,估计的时差为534纳秒,时间为6个跳频周期,即24毫秒。
附图说明
图1为本发明一种跳频信号时差估计方法的具体步骤流程图;
图2为本发明一种跳频信号时差估计方法的相位差模型拟合图。
具体实施方式
附图给出了本发明的一个实施例。
以下结合附图对本发明的具体实施方式进行说明。
术语与符号定义
x1(t)和x2(t) 两个不同位置处接收的跳频信号
T 跳频周期
τ 跳频信号到达两个不同位置处的时差
K 跳频周期数(大于1的整数)
fc 源信号的载频
ft和ft-τ t时刻和t-τ时刻相对于载频fc的跳频频偏
α1和α2 两个不同位置处接收的跳频信号的幅度变化
θ1和θ2 两个不同位置处接收的跳频信号的初相
v1(t)和v2(t) 加性白噪声
angle( ) 表示复数的相角
以两个不同位置处接收的跳频信号x1(t)和x2(t)之间的时差估计为例子。假设跳频信号到达两个不同位置处的时差为τ,在t时刻两个不同位置处接收的跳频信号x1(t)和x2(t)分别为:
x1(t)=α1s(t)cos(2π(fc+ft)t+θ1)+v1(t) (1)
和
x2(t)=α2s(t-τ)cos(2π(fc+ft-τ)(t-τ)+θ2)+v2(t) (2)
其中,s(t)为信号包络,fc为源信号的载频,ft和ft-τ为t时刻和t-τ时刻相对于载频fc的跳频频偏,α1和α2为两个不同位置处接收的跳频信号的幅度变化,θ1和θ2为两个不同位置处接收的跳频信号的初相,v1(t)和v2(t)为加性白噪声。
在t时刻两个不同位置处接收的跳频信号的相位差为:
φ(t)=[2π(fc+ft-τ)(t-τ)+θ2]-[2π(fc+ft)t+θ1] (3)
=-2πft-ττ+2π(ft-τ-ft)t+θ2-θ1-2πfcτ
记跳频周期为T,由于时差τ一般都是远小于跳频周期T,所以在一个跳频周期内的大部分时刻两路接收信号的跳频频偏都满足ft-τ=ft。于是,式(3)又可写成:
φ(t)=-2πft-ττ+θ2-θ1-2πfcτ (4)
与得到(3)式的道理相同,在t+kT时刻两个不同位置处接收的跳频信号的相位差为:
φ(t+kT)=-2πft+kT-ττ+θ2-θ1-2πfcτ (5)
其中,k=0,1,2,…,K-1,K是接收信号期间跳频信号的跳频次数或跳频周期数。
由(3)-(5)式可见,两个不同位置处接收的跳频信号之间的时差与这两个不同位置处接收的跳频信号载波之间的相位差之间存在线性对应关系,该线性对应关系受相对于载频的跳频频偏的影响。若已知或预先估计相对于载频的跳频频偏,则可利用两个不同位置处接收的跳频信号之间的时差与这两个不同位置处接收的跳频信号载波之间的相位差存在的线性对应关系,通过测量这两个不同位置处接收的跳频信号载波之间的多个跳频周期内的相位差,估计这两个不同位置处接收的跳频信号之间的时差。
图1示出根据本发明的设计跳频信号时差估计方法的流程图。流程开始于步骤101。在步骤102,根据一个跳频周期内的载波频率和带宽,对两个不同位置处接收的跳频信号x1(t+kT)和x2(t+kT)分别进行带通滤波,结果为:
和
其中,-T/2≤t≤T/2,a(m)为根据一个跳频周期内的载波频率和带宽设计的带通滤波器的系数,m=1,2,…,M。
在步骤1031,确定两个不同位置处接收的跳频信号在一个跳频周期内的Hilbert变换分别为:
和
其中, 是单位虚数。
在步骤1032,确定两个不同位置处接收的跳频信号在一个跳频周期内的相位分别为:
1(t+kT)=angle(y1(t+kT)) (10)
和
2(t+kT)=angle(y2(t+kT)) (11)
其中,angle( )表示取复数的相角。
在步骤1041,确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差为:
φ(t+kT)=2(t+kT)-1(t+kT) (12)
其中,-T/2≤t≤T/2。
在步骤1042,确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差的中值为:
其中, 表示φ(t+kT)在-T/2≤t≤T/2范围内的中值。
在步骤105,依次确定两个不同位置处接收的跳频信号在K个跳频周期的相位差的中值序列为:
在步骤106,由K个跳频周期内相位差的中值序列对应的线性方程组:
的解:
根据本发明的设计跳频信号时差估计方法的流程结束于步骤107。
定义相位差中值序列拟合误差为
将本发明提出的跳频信号时差估计方法应用于每秒250跳,驻留时间3.5毫秒,跳频频偏依次为-25kHz,225kHz,0kHz,-250kHz,250kHz,775kHz的跳频信号,在两个不同位置处接收的跳频信号时差为493纳秒,K=6,带内信噪比15dB的情况下,相位差中值序列拟合误差如图2所示。根据本发明的跳频信号时差估计方法估计的时差为534纳秒,时间为6个跳频周期,即24毫秒。
虽然已经参考附图对本发明的跳频信号时差估计方法以举例方式进行了描述,但是本发明不限于上述这些细节,本申请含盖权利要求范围之内的各种变型或改变。
工业应用性
将本发明提出的跳频信号时差估计方法应用于跳频辐射信号源时差定位系统,满足跳频辐射信号源时差定位系统对快速时差估计、高精度时差估计性能的要求。
Claims (8)
1、一种跳频信号时差估计方法,其特征在于:利用两个不同位置处接收的跳频信号之间的时差与这两个不同位置处接收的跳频信号载波之间的相位差存在的线性对应关系,通过测量这两个不同位置处接收的跳频信号载波之间的多个跳频周期内的相位差,估计这两个不同位置处接收的跳频信号之间的时差;其具体步骤是:
根据一个跳频周期的载波频率和信号带宽,对两个不同位置处接收的跳频信号进行带通滤波;
确定两个不同位置处接收的跳频信号在一个跳频周期内的Hilbert变换;
利用Hilbert变换结果分别确定两个不同位置处接收的跳频信号在一个跳频周期内的相位;
确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差;
确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差的中值;
依次确定两个不同位置处接收的跳频信号在K(大于1的整数)个跳频周期内的相位差的中值序列;
由K个跳频周期内相位差的中值序列对应的线性方程组的最小二乘解确定跳频信号的时差估计。
2、如权利要求1所述的跳频信号时差估计方法,其特征在于:所述根据一个跳频周期内的载波频率和带宽,对两个不同位置处接收的跳频信号分别进行带通滤波是指对跳频信号x1(t+kT)和x2(t+kT)进行带通滤波,
其滤波结果为:
和
其中,-T/2≤t≤T/2,a(m)为根据一个跳频周期内的载波频率和带宽设计的带通滤波器的系数,m=1,2,...,M。
3、如权利要求1所述的跳频信号时差估计方法,其特征在于:所述确定两个不同位置处接收的跳频信号在一个跳频周期内的Hilbert变换分别为:
和
其中, 是单位虚数。
4、如权利要求1所述的跳频信号时差估计方法,其特征在于:所述利用Hilbert变换结果分别确定两个不同位置处接收的跳频信号在一个跳频周期内的相位分别为 1(t+kT)=angle(y1(t+kT))
和 2(t+kT)=angle(y2(t+kT))
其中,angle( )表示取复数的相角。
5、如权利要求1所述的跳频信号时差估计方法,其特征在于:所述确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差为:
φ(t+kT)=2(t+kT)-1(t+kT)
其中,-T/2≤t≤T/2。
6、如权利要求1所述的跳频信号时差估计方法,其特征在于:所述确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差的中值为:
其中, 表示φ(t+kT)在-T/2≤t≤T/2范围内的中值。
7、如权利要求1所述的跳频信号时差估计方法,其特征在于:所述依次确定两个不同位置处接收的跳频信号在K个跳频周期的相位差的中值序列为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710050704XA CN101247141B (zh) | 2007-12-05 | 2007-12-05 | 一种跳频信号时差估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710050704XA CN101247141B (zh) | 2007-12-05 | 2007-12-05 | 一种跳频信号时差估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101247141A true CN101247141A (zh) | 2008-08-20 |
CN101247141B CN101247141B (zh) | 2011-07-06 |
Family
ID=39947395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200710050704XA Expired - Fee Related CN101247141B (zh) | 2007-12-05 | 2007-12-05 | 一种跳频信号时差估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101247141B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102710286A (zh) * | 2012-06-19 | 2012-10-03 | 电子科技大学 | 一种基于tod信息的快跳频同步方法 |
CN102857961A (zh) * | 2012-09-14 | 2013-01-02 | 中国人民解放军总参谋部第五十七研究所 | 一种存在频率偏移的通信信号时差测量方法 |
CN108020830A (zh) * | 2017-12-05 | 2018-05-11 | 中国人民解放军陆军工程大学 | 基于时频矩阵的跳频测距方法及系统 |
CN111510183A (zh) * | 2020-03-25 | 2020-08-07 | 北京理工大学 | 一种相干快跳频多路并行本振相位计算方法及本振 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1065092C (zh) * | 1994-02-09 | 2001-04-25 | Ntt移动通信网株式会社 | 码分多址移动通信的方法和系统 |
US6493376B1 (en) * | 1997-10-10 | 2002-12-10 | Qualcomm, Inc. | Multi-layered PN code spreading in a multi-user communications system |
JP2006133214A (ja) * | 2004-10-04 | 2006-05-25 | Topcon Corp | 時間差測定装置および測定方法並びに測距装置および測距方法 |
-
2007
- 2007-12-05 CN CN200710050704XA patent/CN101247141B/zh not_active Expired - Fee Related
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102710286A (zh) * | 2012-06-19 | 2012-10-03 | 电子科技大学 | 一种基于tod信息的快跳频同步方法 |
CN102710286B (zh) * | 2012-06-19 | 2014-04-23 | 电子科技大学 | 一种基于tod信息的快跳频同步方法 |
CN102857961A (zh) * | 2012-09-14 | 2013-01-02 | 中国人民解放军总参谋部第五十七研究所 | 一种存在频率偏移的通信信号时差测量方法 |
CN108020830A (zh) * | 2017-12-05 | 2018-05-11 | 中国人民解放军陆军工程大学 | 基于时频矩阵的跳频测距方法及系统 |
CN108020830B (zh) * | 2017-12-05 | 2019-10-18 | 中国人民解放军陆军工程大学 | 基于时频矩阵的跳频测距方法及系统 |
CN111510183A (zh) * | 2020-03-25 | 2020-08-07 | 北京理工大学 | 一种相干快跳频多路并行本振相位计算方法及本振 |
Also Published As
Publication number | Publication date |
---|---|
CN101247141B (zh) | 2011-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101807977B (zh) | 基于波形特征的空时盲自适应抗干扰方法 | |
CN101325807B (zh) | 信号波达方向估计方法 | |
CN103777190B (zh) | 一种双基地mimo雷达高速高机动目标的角度估计方法 | |
CN101645865B (zh) | 基于分数阶傅立叶变换的信道中多径时延和多普勒频移估计方法及实现该方法的系统 | |
CN104007421B (zh) | 基于全变差和压缩感知的罗兰c被动雷达toa估计方法 | |
CN102342075A (zh) | 从一个频率和时间到另一个的信道推断 | |
CN103901417A (zh) | L型阵列mimo雷达低复杂度空间目标二维角度估计方法 | |
CN103546221B (zh) | 一种宽带相干信号波达角估计方法 | |
CN103823217A (zh) | 基于双频发射的双基地mimo雷达高速运动目标参数估计方法 | |
CN101247141B (zh) | 一种跳频信号时差估计方法 | |
US10670693B2 (en) | Position determination system having a deconvolution decoder | |
CN102841344A (zh) | 一种少阵元近场宽带信号源参数估计方法 | |
CN103607361A (zh) | Alpha稳定分布噪声下时频重叠信号的参数预估方法 | |
CN104678372A (zh) | 正交频分复用雷达超分辨距离与角度值联合估计方法 | |
CN105911521A (zh) | 联合无线电信号复包络和载波相位信息的超视距目标直接定位方法 | |
CN102857961A (zh) | 一种存在频率偏移的通信信号时差测量方法 | |
CN102833016B (zh) | 一种在多径环境下针对同源信号时差的测定方法 | |
CN102131288B (zh) | 基于超分辨率uwb信号传播时延估计的室内定位方法及系统 | |
CN103018719A (zh) | 一种oth雷达发射波形的生成方法 | |
CN104360355A (zh) | 抗干扰方法和装置 | |
CN104854473A (zh) | 使用多普勒聚焦的亚奈奎斯特雷达处理 | |
US9261576B2 (en) | Navigation using range measurements to OFDM transmitters | |
CN105158751A (zh) | 一种声矢量阵快速方位估计方法 | |
CN104601512A (zh) | 一种检测相位调制信号载波频偏的方法及系统 | |
CN105549037A (zh) | 一种高精度卫星导航宽带阵列信号生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110706 Termination date: 20141205 |
|
EXPY | Termination of patent right or utility model |