CN113009413A - 基于正交频分复用波形的网络节点间测距方法 - Google Patents

基于正交频分复用波形的网络节点间测距方法 Download PDF

Info

Publication number
CN113009413A
CN113009413A CN202110247187.5A CN202110247187A CN113009413A CN 113009413 A CN113009413 A CN 113009413A CN 202110247187 A CN202110247187 A CN 202110247187A CN 113009413 A CN113009413 A CN 113009413A
Authority
CN
China
Prior art keywords
frequency domain
response
channel
time
noise
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
Application number
CN202110247187.5A
Other languages
English (en)
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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202110247187.5A priority Critical patent/CN113009413A/zh
Publication of CN113009413A publication Critical patent/CN113009413A/zh
Pending legal-status Critical Current

Links

Images

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
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-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
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2649Demodulators
    • H04L27/265Fourier transform demodulators, e.g. fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2657Carrier synchronisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2689Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2691Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation involving interference determination or cancellation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2689Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2695Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with channel estimation, e.g. determination of delay spread, derivative or peak tracking

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Discrete Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于正交频分复用波形的网络节点间测距方法,主要解决现有技术易受无线信道中噪声和多经干扰的问题。其实现方案是:1)接收端提取训练序列;2)根据训练序列计算初始信道频域相应;3)在初始信道频域响应对应的时域形式的信道冲击响应中去除噪声和多经干扰信号,得到处理后的信道频域响应;4)从处理后的信道频域响应中提取子载波间的相位差并计算信号飞行时间;5)利用信号飞行时间计算节点间的距离。本发明通过在信道冲击响应中去除噪声和多经干扰信号,相较于现有正交频分复用波形测距算法提高了一个量级的测距精度,且在存在噪声和多经干扰的无线信道中能稳定的实现分米级的单次测距性能,可用于通信与定位一体化系统。

Description

基于正交频分复用波形的网络节点间测距方法
技术领域
本发明属于测控通信领域,更进一步网络节点间的测距方法,可用于通信与定位一体化系统。
背景技术
随着无线网络的发展,人们对于基于位置的服务需求也越来越大,全球卫星定位系统目前可以达到米级的定位精度,但是卫星定位信号强度低,受环境影响大,并且无法提供更高精度的定位服务。任光亮,空间精密OFDM测距技术[J].西安电子科技大学学报,2009一文描述了一种通过提取OFDM子载波间相位差的实现更高精度的测距方法。该方法具有非常高的理论精度,但是由于信号在无线环境中易受噪声和多径等干扰的影响导致该方法中的子载波相位发生变化,在常见的室内外场景中只能获得米级的测距精度。
武汉大学在其专利号为CN201710464875.0专利文献中公开了一种多载波数字广播信号载波平滑高精度时延估计方法。其首先进行接收信号的定时同步和载波同步;然后通过对单位时间内多个到达符号的多径信息进行统计,按照最早到达路径最多频次检测原则,提取OFDM到达的首径信息;最后通过相位平滑进一步消除噪声以提高时延估计的准确性。该方法通过统计方式识别首径信息和消除测量噪声,提高了测距的稳定性。但是,该方法由于识别的首径时延单位为采样周期,因而精度较低,且处理时延大。
发明内容
本发明的目的在于针对OFDM测距技术的不足,提出一种正交频分复用波形测距中的噪声和多径抑制方法,以减小测距时延,提高OFDM测距的精度。
本发明的技术方案是,在信道频域响应中提取子载波间相位差之前,首先将信道频域响应变换为时域形式且只保留主径及其附近的信道冲击响应信息,然后将信道冲击响应再变换为信道频域响应,进行相位提取等操作直至完成测距功能。其实现步骤包括如下:
(1)对接收信号分别进行定时同步和载波同步,记录定时同步的时间信息ti,并从载波同步的信号中提取出训练序列q;
(2)获取抑制噪声和多径干扰的信道频域响应:
(2a)对训练序列q进行FFT运算,再与本地频域形式的训练序列Q进行相关运算,得到初始信道频域响应H;
(2b)对初始信道频域响应H进行IFFT运算,得到初始信道冲击响应h;
(2c)检测初始信道冲击响应h的峰值,保留峰值及其左右各k个位置上的值,其余位置置零,得到处理后的信道冲击响应r;
(2d)将处理后的信道冲击响应r进行FFT运算,得到抑制噪声和多径干扰的信道频域响应R;
(3)在抑制噪声和多径干扰的信道频域响应R的子载波之间进行相关求和运算,并将运算结果转换为相位角,即为相邻子载波之间的相位差
Figure BDA0002964519890000021
(4)将相邻子载波间的相位差
Figure BDA0002964519890000022
转换为时延信息,并与(1)记录的定时同步时间信息组合,得到精细的信号飞行时间t,再根据该飞行时间t通过TOA测距算法计算出节点间的距离:d=t·c,其中c为光速。
与现有技术相比,本发明具有以下优点:
1.本发明由于在信道冲击响应上去除主径外的噪声和多径信息,因此具有较强的抗干扰能力,提高了至少一个量级的测距精度。
2.本发明利用信道冲击响应的峰值位置补偿定时同步的误差,提高了单次测距的稳定性,从而减小了获取高精度距离值的时延。
3.本发明中由于分别从通信系统的定时同步和信道估计中获取信号飞行时间,其测距方法适合用于通信和导航的一体化设计。
附图说明
图1为本发明的实现流程框图;
图2为本发明仿真中使用的信道对子载波干扰图;
图3为用本发明对图2数据进行干扰抑制后的仿真图;
图4为用本发明在单径信道中的测距性能仿真图。
图5为用本发明在多径信道中的测距性能仿真图。
具体实施方式
下面结合附图对本发明的实施例和效果做进一步的描述:
本实施例的实施场景基于通信与定位一体化的无线自组网,自组网内各节点在分配的时隙中发送信号,接收节点在解析信号的同时提取信号飞行时间,最后利用信号飞行时间计算出发送节点到接收节点的距离。
参照图1,对本实例的实现步骤如下:
步骤1,从接收信号中获取初始信道频域响应。
1.1)接收端对接收信号进行检测,通过本地相关算法实现定时同步,并记录定时同步的时间信息ti
1.2)根据定时同步指示的信号位置对信号进行载波同步,然后提取出时域形式的训练序列q:
q=(q1,q2,…,qk,…,qN),
其中,qk是序列q中的第k个值,其表示为:
Figure BDA0002964519890000031
其中,j为虚部符号,μ为序列根植,N为序列长度,cf为N对2取模的值;
1.3)将时域形式的训练序列q进行FFT运算,再与本地频域形式的训练序列Q进行相关运算,得到初始信道频域响应H:
H=FFT(q)·Q*
其中,Q*表示Q的共轭形式。
步骤2,抑制初始信道频域响应中的噪声和多径干扰。
2.1)将初始信道频域响应H进行IFFT运算,得到初始信道冲击响应h:
h=IFFT(H);
2.2)寻找初始信道冲击响应h的峰值位置并记为m,保留m位置及其左右各k个值,将其余位置的值设为零,得到处理后的信道冲击响应r:
r=(r1,r2,…ri,…,rN),
其中,ri是序列r中的第i个值,i从1到N,N为序列长度,ri的取值为:
Figure BDA0002964519890000041
通过此操作可以去除k倍的采样周期外的噪声和多径干扰,提高的信噪比为:
Figure BDA0002964519890000042
当N=2048,且k=1时,SNR=28.3dB;
2.3)对处理后的信道冲击响应r进行FFT运算,得到抑制干扰的信道频域响应R:
R=FFT(r)。
步骤3,提取包含小数倍的信号飞行时间。
3.1)在抑制噪声和多径干扰的信道频域响应R的子载波之间进行相关求和运算,得到间隔为Δ子载波间的相位差信息W:
Figure BDA0002964519890000043
其中,N为R的长度,Ri表示R的第i个值,
Figure BDA0002964519890000044
表示R的第i+Δ个值的共轭形式;
3.2)对W进行求相位角运算,得到间隔Δ子载波间的相位差
Figure BDA0002964519890000045
为:
Figure BDA0002964519890000046
其中,angle为求相位角运算;
3.3)将子载波间的相位差
Figure BDA0002964519890000047
转换为包含小数倍的信号飞行时间tf
Figure BDA0002964519890000048
其中,tf中允许存在整数倍的信号飞行时间,并能够补偿定时同步的误差。
步骤4,计算发送节点与接收节点的距离。
根据定时同步记录的时间信息ti和包含小数倍的信号飞行时间tf,则发送节点到接收节点的距离d为:
d=(ti+tf)·c,其中,c为光速。
本发明的效果可通过以下仿真实验进一步的说明:
1.仿真实验条件:
仿真实验的硬件平台:处理器为Intel i7-8700 CPU,主频为3.2GHz,内存16GB。
仿真实验的软件平台为:Windows 10操作系统和MATLAB R2018b。
2.仿真内容与结果分析:
仿真1,使用MATLAB软件模拟包含长度为2048的训练序列经过信噪比为10dB的多径信道后的子载波相位干扰情况,结果如图2所示,其中图2(a)表示信号飞行时间产生的子载波相位,图2(b)表示图2(a)中信号被噪声和多径干扰后的子载波相位。
仿真2,使用本发明对图2(b)进行干扰抑制,结果如图3所示,从图(3)可以看出子载波的相位噪声得到了明显的抑制,并且中间约1400个子载波之间的相位差与图2(a)相似。
仿真3,为了验证本发明对噪声干扰的抑制作用,在不同信噪比下对测距结果的误差进行统计,结果如图4所示,从图4可知现有正交频分复用波形测距算法的测距误差在信噪比为0dB时有90%的概率小于4.92m,使用本发明的测距误差在信噪比为0dB时有90%的概率小于0.31m,其性能优于现有算法在信噪比为10dB的测距性能,可见本发明具有很强的降噪能力。
仿真4,在仿真3的基础上,增加包含多径干扰的莱斯信道,在不同信噪比下对测距结果的误差进行统计,结果如图5所示,其中现有正交频分复用波形测距算法的测距误差小于10m的概率均小于75%,而使用本发明的测距误差均有90%的概率小于0.60m,可见本发明的测距精度提高了一个量级,并且具有更高的稳定性。
以上仿真结果表明:本发明中通过在时域的信道冲击响应中去除噪声和多径干扰信号,使正交频分复用波形测距技术具有很强的抗干扰能力,显著提高了测距精度和稳定性。

Claims (7)

1.一种基于正交频分复用波形的网络节点间测距方法,其特征在于,包括如下:
(1)对接收信号分别进行定时同步和载波同步,记录定时同步的时间信息ti,并从载波同步的信号中提取出训练序列q;
(2)获取抑制噪声和多径干扰的信道频域响应:
(2a)对训练序列q进行FFT运算,再与本地频域形式的训练序列Q进行相关运算,得到初始信道频域响应H;
(2b)对初始信道频域响应H进行IFFT运算,得到初始信道冲击响应h;
(2c)检测初始信道冲击响应h的峰值,保留峰值及其左右各k个位置上的值,其余位置置零,得到处理后的信道冲击响应r;
(2d)将处理后的信道冲击响应r进行FFT运算,得到抑制噪声和多径干扰的信道频域响应R;
(3)在抑制噪声和多径干扰的信道频域响应R的子载波之间进行相关求和运算,并将运算结果转换为相位角,即为相邻子载波之间的相位差
Figure FDA0002964519880000011
(4)将相邻子载波间的相位差
Figure FDA0002964519880000012
转换为时延信息,并与(1)记录的定时同步时间信息组合,得到精细的信号飞行时间t,再根据该飞行时间t通过TOA测距算法计算出节点间的距离:d=t·c,其中c为光速。
2.根据权利要求1所述的方法,其特征在于,(1)中提取的时域形式的训练序列q,表示如下:
q=(q1,q2,…,qk,…,qN),
其中,qk是序列q中的第k个值,
Figure FDA0002964519880000013
j为虚部符号,μ为序列根植,N为序列长度,cf为N对2取模的值。
3.根据权利要求1所述的方法,其特征在于,(2a)中的初始信道频域响应H,表示为:
H=FFT(q)·Q*
其中,Q*表示Q的共轭形式。
4.根据权利要求1所述的方法,其特征在于,(2b)中的初始信道冲击响应h,表示为:
h=IFFT(H)。
5.根据权利要求1所述的方法,其特征在于,(2c)中处理后的信道冲击响应r,表示为:
r=(r1,r2,…ri,…,rN),
其中,ri是序列r中的第i个值,i从1到N,N为序列长度,ri的取值为:
Figure FDA0002964519880000021
其中,hi表示初始信道冲击响应h的第i个值,m为h的峰值位置。
6.根据权利要求1所述的方法,其特征在于,(3)中的相邻子载波之间的相位差
Figure FDA0002964519880000022
表示如下:
Figure FDA0002964519880000023
其中,N为R的长度,Δ为子载波间隔,Ri表示R的第i个值,
Figure FDA0002964519880000024
表示R的第i+Δ个值的共轭形式,angle为求相位角运算。
7.根据权利要求1所述的方法,其特征在于,(4)中计算信号飞行时间t,通过如下公式计算:
Figure FDA0002964519880000025
其中,N为序列长度,Δ为子载波间隔,ti为记录定时同步的时间信息。
CN202110247187.5A 2021-03-05 2021-03-05 基于正交频分复用波形的网络节点间测距方法 Pending CN113009413A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110247187.5A CN113009413A (zh) 2021-03-05 2021-03-05 基于正交频分复用波形的网络节点间测距方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110247187.5A CN113009413A (zh) 2021-03-05 2021-03-05 基于正交频分复用波形的网络节点间测距方法

Publications (1)

Publication Number Publication Date
CN113009413A true CN113009413A (zh) 2021-06-22

Family

ID=76407385

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110247187.5A Pending CN113009413A (zh) 2021-03-05 2021-03-05 基于正交频分复用波形的网络节点间测距方法

Country Status (1)

Country Link
CN (1) CN113009413A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114222242A (zh) * 2021-11-08 2022-03-22 中国电子科技集团公司第五十四研究所 基于数字相控阵体制的通信定位一体化方法及系统
WO2023015556A1 (zh) * 2021-08-13 2023-02-16 Oppo广东移动通信有限公司 距离确定方法、装置、设备及存储介质
CN117538854A (zh) * 2024-01-09 2024-02-09 腾讯科技(深圳)有限公司 测距方法、装置、计算机设备和计算机可读存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002063866A2 (en) * 2001-02-02 2002-08-15 Rosum Corporation Position location using broadcast digital television signals
CN106789791A (zh) * 2017-02-28 2017-05-31 北京科技大学 基于共轭对称训练序列的移动通信系统载波频偏估计方法
CN107302514A (zh) * 2017-06-19 2017-10-27 武汉大学 一种多载波数字广播信号载波平滑高精度时延估计方法
CN108040028A (zh) * 2017-12-22 2018-05-15 中国人民解放军国防科技大学 一种基于本地序列互相关检测的ofdm系统抗干扰信号检测与同步方法
CN110691056A (zh) * 2019-11-27 2020-01-14 中国人民解放军国防科技大学 一种无线通信系统的同步方法、装置、设备及存储介质
CN111884978A (zh) * 2020-07-30 2020-11-03 电子科技大学 一种基于ofdm抗脉冲噪声的符号同步方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002063866A2 (en) * 2001-02-02 2002-08-15 Rosum Corporation Position location using broadcast digital television signals
CN106789791A (zh) * 2017-02-28 2017-05-31 北京科技大学 基于共轭对称训练序列的移动通信系统载波频偏估计方法
CN107302514A (zh) * 2017-06-19 2017-10-27 武汉大学 一种多载波数字广播信号载波平滑高精度时延估计方法
CN108040028A (zh) * 2017-12-22 2018-05-15 中国人民解放军国防科技大学 一种基于本地序列互相关检测的ofdm系统抗干扰信号检测与同步方法
CN110691056A (zh) * 2019-11-27 2020-01-14 中国人民解放军国防科技大学 一种无线通信系统的同步方法、装置、设备及存储介质
CN111884978A (zh) * 2020-07-30 2020-11-03 电子科技大学 一种基于ofdm抗脉冲噪声的符号同步方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
任光亮等: "空间精密OFDM测距技术", 《西安电子科技大学学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023015556A1 (zh) * 2021-08-13 2023-02-16 Oppo广东移动通信有限公司 距离确定方法、装置、设备及存储介质
CN114222242A (zh) * 2021-11-08 2022-03-22 中国电子科技集团公司第五十四研究所 基于数字相控阵体制的通信定位一体化方法及系统
CN114222242B (zh) * 2021-11-08 2024-01-26 中国电子科技集团公司第五十四研究所 基于数字相控阵体制的通信定位一体化方法及系统
CN117538854A (zh) * 2024-01-09 2024-02-09 腾讯科技(深圳)有限公司 测距方法、装置、计算机设备和计算机可读存储介质
CN117538854B (zh) * 2024-01-09 2024-04-09 腾讯科技(深圳)有限公司 测距方法、装置、计算机设备和计算机可读存储介质

Similar Documents

Publication Publication Date Title
CN113009413A (zh) 基于正交频分复用波形的网络节点间测距方法
CN110224968B (zh) 一种ofdm通信系统中的帧定时同步方法和装置
US11695672B2 (en) Communication system determining time of arrival using matching pursuit
CN107257324B (zh) 一种ofdm系统中的时频联合同步方法及装置
CN107911329B (zh) 一种信号分析仪ofdm信号解调方法
WO2010060732A1 (en) Frequency offset estimation
WO2017174004A1 (zh) 载波同步方法和装置
CN114268528B (zh) 一种帧捕获同步方法及装置
WO2017174003A1 (zh) 定时同步方法和装置
CN113259291B (zh) 利用水声连续信号动态多普勒跟踪实现的相位补偿方法
CN114615122B (zh) 通信信号的频偏确定方法及装置
CN108989259B (zh) 无线综测仪窄带物理上行共享信道的时偏估计方法及系统
CN113438730A (zh) 一种基于gfdm信号的无线定位方法
CN112738000B (zh) 一种pss分块互相关检测方法
CN101621493B (zh) Ofdm的频率偏移估计的判决方法
JP3538104B2 (ja) マルチキャリア信号検波装置
CN105187351B (zh) 一种多径信道下的ofdm定时同步检测方法
CN107276654B (zh) 信号处理方法和系统
CN107276926B (zh) 信道估计方法和装置
CN101102114A (zh) 基于二值全通序列保护间隔填充的频域信道估计方法
CN113315734A (zh) 用于卫星的载波同步方法、装置、电子设备及可读介质
US8983007B2 (en) Receiver apparatus, reception method, communication system, and communication method
CN107294889B (zh) 载波同步方法和装置
CN107276925B (zh) 信道估计方法和装置
CN107277913B (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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20210622