CN109541531B - 一种降低频差对时间延迟估计影响的方法 - Google Patents
一种降低频差对时间延迟估计影响的方法 Download PDFInfo
- Publication number
- CN109541531B CN109541531B CN201811398832.8A CN201811398832A CN109541531B CN 109541531 B CN109541531 B CN 109541531B CN 201811398832 A CN201811398832 A CN 201811398832A CN 109541531 B CN109541531 B CN 109541531B
- Authority
- CN
- China
- Prior art keywords
- phase
- frequency
- frequency difference
- signals
- time delay
- 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.)
- Active
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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Complex Calculations (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提出了一种降低频差对时间延迟估计影响的方法,包括以下步骤:首先,获取时域信号y(n)的相位角,记作phase_y(n),并进行初始赋值x0(n)=x(n),k0=0,τ0=0;然后,进行循环迭代运算,利用信号的时域相位进行比对,进行频差补偿,计算出不受频率分辨率限制的频差;最后,采用广义相关算法对频差补偿后的信号做相关,获取信号之间的时延估计值。本发明基于傅里叶变换频移特性进行频率补偿,可以对于远小于频率分辨率的频率误差进行有效补偿,不受频率分辨率对频差修正的限制,极大地提高了时延估计的精度。
Description
技术领域
本发明涉及测试技术领域,特别涉及一种降低频差对时间延迟估计影响的方法。
背景技术
TDOA(Time Difference Of Arrival)定位技术,同DOA(Direction Of Arrival)定位技术一样,是一种经典的无源定位方式,由于其量测的TDOA精度较DOA精度高,得到了广泛的研究和应用,其定位精度很大程度上取决于时延估计的精度。
在进行时延估计时,一般都是假设各接收机之间不存在频率差,即本振是完全相同的,认为中频信号之间不存在频率误差的影响。而在实际应用中,两台接收机难以保证本振频率的一致,因此信号之间存在着频率差,即使是远小于频率分辨率的频率误差,也可能对时延估计精度产生致命的影响,尤其降低了窄带信号和非幅度调制信号的时延估计精度,导致信号的时延估计误差很大甚至时延估计的失效。
现有的频差修正方法,主要通过频域相关函数对频差做准确估计和补偿,对补偿后的信号做时延估计。当频差小于离散傅立叶变换的频率分辨率时,无法估计两路信号的频差,进而不能做频差补偿并基于频差补偿进行时延估计。
发明内容
本发明提出了一种降低频差对时间延迟估计影响的方法,利用信号的时域相位进行比对,计算出不受频率分辨率限制的频差,基于频差补偿,有效的提高信号时延的估计精度。
本发明的技术方案是这样实现的:
一种降低频差对时间延迟估计影响的方法,两台接收机分别接收来自同一个辐射源发射的信号,接收的信号分别是x(n)和y(n),x(n)和y(n)的关系表示为:
y(n)=Ax(n-τ)e-jΔw(n-τ)
其中,A表示归一化的幅度比值,Δw表示本振不一致造成的两信号之间的频差,τ表示时间延迟;
包括以下步骤:
首先,获取时域信号y(n)的相位角,记作phase_y(n),并进行初始赋值x0(n)=x(n),k0=0,τ0=0;
然后,进行循环迭代运算,利用信号的时域相位进行比对,进行频差补偿,计算出不受频率分辨率限制的频差;
最后,采用广义相关算法对频差补偿后的信号做相关,获取信号之间的时延估计值。
可选地,所述迭代运算的过程包括以下步骤:
步骤(b),获取时域信号xi(n)的相位角,表示为phase_xi(n);
步骤(c),计算y(n)的相位角和xi(n)的相位角之间的相位差:phase_errori(n)=phase_y(n)-phase_xi(n);
步骤(d),剔除phase_errori(n)中奇异点,表示为phase_ERRi(n);
步骤(e),对phase_ERRi(n)进行线性拟合,得到拟合直线的斜率ki;
步骤(g),计算xi(n)和y(n)的时差τi;
其中,Num表示循环迭代的次数,ki表示相位差的拟合斜率。
本发明的有益效果是:
基于傅里叶变换频移特性进行频率补偿,可以对于远小于频率分辨率的频率误差进行有效补偿,不受频率分辨率对频差修正的限制,极大地提高了时延估计的精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一种降低频差对时间延迟估计影响的方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
两台接收机A和B分别接收来自同一个辐射源发射的信号,接收的信号分别是x(n)和y(n)。在实际工程应用中,接收机一般以IQ数据的形式给出接收信号,因此,在处理过程中,主要针对复数数据进行处理。
如果不考虑噪声的影响,x(n)和y(n)的关系可以表示为:
y(n)=Ax(n-τ)e-jΔw(n-τ)
其中,A表示归一化的幅度比值,Δw表示本振不一致造成的两信号之间的频差,τ表示时间延迟。为了对时间延迟进行有效的估计,需要对Δw进行有效的补偿,从而消除频差的影响。
为了消除频差的影响,如图1所示,本发明提出了一种降低频差对时间延迟估计影响的方法,包括以下步骤:
首先,获取时域信号y(n)的相位角,记作phase_y(n),并进行初始赋值x0(n)=x(n),k0=0,τ0=0。x(n)表示原始数据序列,赋值给x0(n),下次更新赋值给x1(n),继续迭代,赋值给xi(n),xi(n)表示对序列x(n)进行第i次修正后得到的序列,i=0,1,2,….。
然后,进行循环迭代运算,利用信号的时域相位进行比对,进行频差校正,计算出不受频率分辨率限制的频差。
最后,采用广义相关算法对频差校正后的信号做相关,获取信号之间的时延估计值。
所述迭代运算的过程包括以下步骤:
步骤(a),对信号xi(n)的频谱进行校正,校正公式为:反傅里叶变换得到新的时域信号xi(n)。xi(n)表示对序列x(n)进行第i次修正后得到的序列;i=0,1,2,…。Fi(w)表示序列xi(n)的傅里叶变换。
步骤(b),获取新的时域信号xi(n)的相位角,表示为phase_xi(n)。
步骤(c),计算时域信号y(n)的相位角和时域信号xi(n)的相位角之间的相位差:phase_errori(n)=phase_y(n)-phase_xi(n)。
步骤(d),剔除phase_errori(n)中奇异点,此时表示为phase_ERRi(n)。
步骤(e),对phase_ERRi(n)进行线性拟合,得到拟合直线的斜率ki。
步骤(g),计算xi(n)和y(n)的时差τi。
其中,Num表示循环迭代的次数,ki表示相位差的拟合斜率。每一次迭代都对x(n)做一次纠正,采用预先设置的Num次迭代基本可以消除信号之间的频差。
上述迭代过程步骤(e)得到的线性拟合表达式是:phase_ERRi(n)=kin+b,其中b表示拟合直线与y轴的交点,ki表示斜率,该斜率对应两接收机之间的频差,由相位拟合曲线得到,因而不受频率分辨率的影响。
本发明与目前频差修正方法相比,主要基于傅里叶变换频移特性进行频率校正,可以对于远小于频率分辨率的频率误差进行有效补偿,不受频率分辨率对频差修正的限制,极大地提高了时延估计的精度。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种降低频差对时间延迟估计影响的方法,其特征在于,两台接收机分别接收来自同一个辐射源发射的信号,接收的信号分别是x(n)和y(n),x(n)和y(n)的关系表示为:
y(n)=Ax(n-τ)e-jΔw(n-τ)
其中,A表示归一化的幅度比值,Δw表示本振不一致造成的两信号之间的频差,τ表示时间延迟;
包括以下步骤:
首先,获取时域信号y(n)的相位角,记作phase_y(n),并进行初始赋值x0(n)=x(n),k0=0,τ0=0;
然后,进行循环迭代运算,利用信号的时域相位进行比对,进行频差补偿,计算出不受频率分辨率限制的频差;
最后,采用广义相关算法对频差补偿后的信号做相关,获取信号之间的时延估计值;
所述迭代运算的过程包括以下步骤:
步骤(a),对信号xi(n)的频谱进行校正,校正公式为:反傅里叶变换得到新的时域信号xi(n);xi(n)表示对序列x(n)进行第i次修正后得到的序列;i=0,1,2,…;Fi(w)表示序列xi(n)的傅里叶变换;
步骤(b),获取新的时域信号xi(n)的相位角,表示为phase_xi(n);
步骤(c),计算时域信号y(n)的相位角和时域信号xi(n)的相位角之间的相位差:phase_errori(n)=phase_y(n)-phase_xi(n);
步骤(d),剔除phase_errori(n)中奇异点,表示为phase_ERRi(n);
步骤(e),对phase_ERRi(n)进行线性拟合,得到拟合直线的斜率ki;
步骤(g),计算xi(n)和y(n)的时差τi;
其中,Num表示循环迭代的次数,ki表示相位差的拟合斜率;每一次迭代都对x(n)做一次纠正;
迭代过程步骤(e)得到的线性拟合表达式是:phase_ERRi(n)=kin+b,其中b表示拟合直线与y轴的交点,ki表示斜率,该斜率对应两接收机之间的频差,由相位拟合曲线得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811398832.8A CN109541531B (zh) | 2018-11-22 | 2018-11-22 | 一种降低频差对时间延迟估计影响的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811398832.8A CN109541531B (zh) | 2018-11-22 | 2018-11-22 | 一种降低频差对时间延迟估计影响的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109541531A CN109541531A (zh) | 2019-03-29 |
CN109541531B true CN109541531B (zh) | 2020-07-31 |
Family
ID=65849374
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811398832.8A Active CN109541531B (zh) | 2018-11-22 | 2018-11-22 | 一种降低频差对时间延迟估计影响的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109541531B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110275133B (zh) * | 2019-06-27 | 2021-06-29 | 清华大学 | 一种基于虚拟到达时差的非共视脉冲信号无源定位方法 |
CN110261821B (zh) * | 2019-07-18 | 2021-12-10 | 中电科思仪科技股份有限公司 | 一种基于信息解调相关的时差估计方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007004307A1 (en) * | 2005-07-06 | 2007-01-11 | National Institute Of Information And Communications Technology Incorporated Administrative Agency | Delay estimation apparatus and method |
CN101692163A (zh) * | 2009-09-24 | 2010-04-07 | 中国计量科学研究院 | 一种频率标准远程校准方法及其*** |
CN104237871A (zh) * | 2013-06-08 | 2014-12-24 | 中国科学院声学研究所 | 一种基于相位补偿的时延差估计方法 |
US8948320B2 (en) * | 2010-07-30 | 2015-02-03 | National Instruments Corporation | Frame and symbol timing recovery for unbursted packetized transmissions using constant-amplitude continuous-phase frequency-modulation |
CN105309018A (zh) * | 2013-10-31 | 2016-02-03 | 华为技术有限公司 | 一种传输方法和装置 |
CN107645770A (zh) * | 2016-07-13 | 2018-01-30 | 华为技术有限公司 | 一种相位校准方法及装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8174444B2 (en) * | 2009-09-26 | 2012-05-08 | Rincon Research Corporation | Method of correlating known image data of moving transmitters with measured radio signals |
CN103809192A (zh) * | 2014-02-25 | 2014-05-21 | 浙江理工大学 | 一种gnss接收机的动态校正算法 |
US9297884B1 (en) * | 2014-04-30 | 2016-03-29 | The United States Of America As Represented By The Director, National Security Agency | Systems and methods of frequency estimation for geolocation |
US9310464B2 (en) * | 2014-07-09 | 2016-04-12 | Deere & Company | Determining location of a receiver with a multi-subcarrier signal |
CN106488557B (zh) * | 2016-10-28 | 2019-05-03 | 中国人民解放军信息工程大学 | 基于WiFi信号的LDCC-PDF分级时延估计方法 |
CN107870317B (zh) * | 2017-11-10 | 2020-02-14 | 中国电子科技集团公司第四十一研究所 | 一种基于时差计算的窄带信号tdoa定位方法及装置 |
-
2018
- 2018-11-22 CN CN201811398832.8A patent/CN109541531B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007004307A1 (en) * | 2005-07-06 | 2007-01-11 | National Institute Of Information And Communications Technology Incorporated Administrative Agency | Delay estimation apparatus and method |
CN101692163A (zh) * | 2009-09-24 | 2010-04-07 | 中国计量科学研究院 | 一种频率标准远程校准方法及其*** |
US8948320B2 (en) * | 2010-07-30 | 2015-02-03 | National Instruments Corporation | Frame and symbol timing recovery for unbursted packetized transmissions using constant-amplitude continuous-phase frequency-modulation |
CN104237871A (zh) * | 2013-06-08 | 2014-12-24 | 中国科学院声学研究所 | 一种基于相位补偿的时延差估计方法 |
CN105309018A (zh) * | 2013-10-31 | 2016-02-03 | 华为技术有限公司 | 一种传输方法和装置 |
CN107645770A (zh) * | 2016-07-13 | 2018-01-30 | 华为技术有限公司 | 一种相位校准方法及装置 |
Non-Patent Citations (2)
Title |
---|
《Relations between fractional operations and time-frequency distributions, and their application》;Soo-Chang Pei、Jian-Jiun Ding;《IEEE Transactions on Signal Processing》;20010831;第49卷(第8期);1638-1655 * |
《一种基于频差补偿的相位谱时延估计方法》;邱天爽、尤国红、沙岚、赵小平、高阳;《大连理工大学学报》;20120131;第52卷(第1期);90-94 * |
Also Published As
Publication number | Publication date |
---|---|
CN109541531A (zh) | 2019-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
So et al. | Simple formulae for bias and mean square error computation [DSP tips and tricks] | |
US9449594B2 (en) | Adaptive phase difference based noise reduction for automatic speech recognition (ASR) | |
CN103760522B (zh) | 用于时差估计与多站时钟误差校准的方法及*** | |
CN109541531B (zh) | 一种降低频差对时间延迟估计影响的方法 | |
CN103399307B (zh) | 一种校正通道误差的方法及装置 | |
CN107085140B (zh) | 基于改进的SmartDFT算法的非平衡***频率估计方法 | |
CN107024855A (zh) | 基于双数字到时间转换器(dtc)的差分相关双采样dtc校准 | |
CN104142425B (zh) | 一种正弦信号频率估计的相位匹配方法 | |
US10371557B2 (en) | Equaliser for the intermediate frequency channel of an FMCW fill level sensor | |
JP2017215208A (ja) | パルスレーダ | |
Jin et al. | Robust cyclic beamforming against cycle frequency error in Gaussian and impulsive noise environments | |
Zhu et al. | Noise influence on estimation of signal parameter from the phase difference of discrete Fourier transforms | |
CN107966676A (zh) | 复杂噪声环境下阵列天线角度和信源个数联合估计方法 | |
JP5992129B2 (ja) | キャリブレーション装置 | |
Wang et al. | Active sampling rate calibration method for acoustic sensor networks | |
WO2016165416A1 (zh) | 一种信噪比确定方法及装置 | |
Chakraborty | A new adaptive filter for estimating and tracking the delay and the amplitude of a sinusoid | |
JP2004150825A (ja) | スペクトル分析装置およびスペクトル分析方法 | |
CN105681240B (zh) | 一种适应于低信噪比环境下的高精度频偏估计方法 | |
CN107679014B (zh) | 一种基于内点罚函数和相邻相关数据段的频偏估计方法 | |
JP2009031425A (ja) | 雑音推定装置と方法およびプログラム | |
CN103905348A (zh) | 基于相关函数线性预测和泰勒分解的双阶段频率估计方法 | |
WO2016062033A1 (zh) | 频率补偿处理方法及装置 | |
CN108540284B (zh) | 一种连续变量量子密钥分发后处理外差探测相位补偿方法 | |
Díaz et al. | Channel phase calibration based on Savitzky-Golay filter in time-domain for OFDM systems |
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 | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province Patentee after: CLP kesiyi Technology Co.,Ltd. Address before: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province Patentee before: CHINA ELECTRONIC TECHNOLOGY INSTRUMENTS Co.,Ltd. |