CN118074827A - 基于最大似然估计和线性插值的时钟漂移校正方法 - Google Patents
基于最大似然估计和线性插值的时钟漂移校正方法 Download PDFInfo
- Publication number
- CN118074827A CN118074827A CN202410353200.9A CN202410353200A CN118074827A CN 118074827 A CN118074827 A CN 118074827A CN 202410353200 A CN202410353200 A CN 202410353200A CN 118074827 A CN118074827 A CN 118074827A
- Authority
- CN
- China
- Prior art keywords
- clock drift
- cir
- linear interpolation
- maximum likelihood
- likelihood 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 37
- 238000012937 correction Methods 0.000 title claims abstract description 25
- 238000007476 Maximum Likelihood Methods 0.000 title claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 27
- 238000005070 sampling Methods 0.000 claims abstract description 24
- 230000005540 biological transmission Effects 0.000 claims abstract description 11
- 230000000694 effects Effects 0.000 description 3
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 229910052701 rubidium Inorganic materials 0.000 description 2
- IGLNJRXAVVLDKE-UHFFFAOYSA-N rubidium atom Chemical compound [Rb] IGLNJRXAVVLDKE-UHFFFAOYSA-N 0.000 description 2
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W56/00—Synchronisation arrangements
- H04W56/004—Synchronisation arrangements compensating for timing error of reception due to propagation delay
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W56/00—Synchronisation arrangements
- H04W56/0055—Synchronisation arrangements determining timing error of reception due to propagation delay
- H04W56/0065—Synchronisation arrangements determining timing error of reception due to propagation delay using measurement of signal travel time
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Error Detection And Correction (AREA)
Abstract
一种基于最大似然估计和线性插值的时钟漂移校正方法,在视距传输情况下,周期性地在一已知收发机间距离的位置处采集信道冲激响应(CIR)数据以获取时钟漂移的采样值后,对实测过程中产生的时钟漂移进行校正,本发明考虑了实际测量系统分辨率的有限性,采用最大似然估计获取时钟漂移采样值,以线性插值法对时钟漂移进行分段线性建模,实现了对时钟漂移的精确校正。
Description
技术领域
本发明涉及的是一种无线通信领域的技术,具体是一种基于最大似然估计和线性插值的时钟漂移校正方法。
背景技术
为支撑无线通信研究,需要开展信道测量来得到信道原始数据,其中一种重要的信道测量方法是时域信道测量法。在时域信道测量中,发射端和接收端常常使用独立的铷钟来提供时钟参考。由于硬件限制,收发端铷钟的频率无法达到完全一致,加以测量系统众多器件引入的相位噪声等,导致在测量过程中发生时钟漂移。为得到精确的多径分量时延测量结果,亟需准确有效的时钟漂移校正算法。
发明内容
本发明针对现有技术时间信息的记录依赖于时钟源节点,因此时间信息的分辨率受时钟源节点的分辨率所限,无法满足信道测量中的超高时间分辨率需求,提出一种基于最大似然估计和线性插值的时钟漂移校正方法,考虑了实际测量系统分辨率的有限性,采用最大似然估计获取时钟漂移采样值,以线性插值法对时钟漂移进行分段线性建模,实现了对时钟漂移的精确校正。
本发明是通过以下技术方案实现的:
本发明涉及一种基于最大似然估计和线性插值的时钟漂移校正方法,包括:
步骤1、时钟漂移采样:在视距传输情况下,周期性地在一已知收发机间距离的位置处采集信道冲激响应(CIR)数据以获取时钟漂移的采样值,具体为:
1.1基于CIR数据,粗略估计视距传输径的时延其中:Δτ为测量系统的时域分辨率,所测CIR最大值对应的索引j=argmaxk|h[k]|,h[k]为第k个CIR采样处的路径增益;
1.2采用最大似然估计获取第i次在参考位置采样时的视距传输径的精确估计 其中:测量得到的频域响应/> 为离散傅里叶变换,频域采样点fk=f1+(k-1)Δf,f1为所测频段的起始频率,Δf为测量系统的频域分辨率,为搜索区间。
1.3基于所测得的视距传输径时延的精确估计以及根据实际收发机距离所计算得出的视距传输径时延理论值,得出时钟漂移的采样值其中:τTheo为视距传输径时延的理论值。
步骤2、时钟漂移校正:对实测过程中产生的时钟漂移进行校正,具体为:
2.1基于线性插值法,对于任意时刻t测得的信道冲激响应,预测其涉及的时钟漂移其中:/>ti和tj为第i次和第j次在参考位置处采样时的时间,Δτi和Δτj为第i次和第j次采样时的时钟漂移,/>和/>为向下取整和向上取整。
2.2对时刻t处测得的信道冲激响应进行校正其中:⊙为同样大小矩阵间的哈达玛积(Hadamard product),即逐元素乘积,F=[f1,f1+Δf,f1+2Δf,…,f1+(k-1)Δf,…,f1+(K-1)Δf]为包括各个频域采样点对应频率的向量,K为测得的CIR长度,/>为离散傅里叶逆变换。
技术效果
本发明基于最大似然估计的时钟漂移采样及基于线性插值法的时钟漂移校正,突破了系统硬件带来的时间分辨率限制,可以精确估计系统时钟漂移值;通过基于线性插值法的时钟漂移校正,获得了具有精确绝对时延测量的CIR数据。
附图说明
图1为本发明系统示意图;
图2为本发明流程图;
图3为实施例效果示意图。
具体实施方式
如图1所示,为本实施例涉及一种用于无线信道测量的基于最大似然估计和线性插值的时钟漂移校正系统,包括:时钟漂移采样单元及时钟漂移校正单元,其中:时钟漂移采样单元根据在固定收发机位置采集的CIR数据,基于最大似然估计法,获得某一时刻的系统时钟漂移值,时钟漂移校正单元根据时钟漂移采样单元记录的时钟漂移值,采用线性插值法,获得系统在任意时刻的时钟漂移预测值,并对采集到的CIR数据进行时钟漂移校正,获得准确的信道测量结果。
经过具体实际实验,如图2所示,在典型的无线信道测量场景下,使用基于相关法的信道测量系统得到信道的冲激响应数据。测量中心频率为12GHz,带宽为1.536GHz,信道冲激响应采样点数为2048.为精确测量时钟漂移,发射机和接收机通过射频线缆直连。每秒记录一次信道冲激响应,以信道冲激响应中最强多径分量的时延作为计量,该数值相比第一次记录的最强多径的时延之间的差值即为时钟漂移的真实值。
为验证不同方法的准确性,取每30秒的时钟漂移记录作为参考值,使用文献中现有的线性模型法、线性插值法,以及本发明提出的基于最大似然估计和线性插值的方法,对时钟漂移进行预测,并与真实值比较,采用时钟漂移误差作为评估准则EΔτ=|ΔτP-ΔτT|,其中:ΔτP-ΔτT分别为时钟漂移误差的预测值和真实值。
如表1所示,为现有方法和本发明时钟漂移误差比较。可以明显看到,现有方法中的线性模型法性能最差,这是由于时钟漂移多符合分段线性,因此线性模型精度有限。线性插值法效果较好,但由于测量系统时域分辨率有限(本实施例为0.65ns),导致精度受限。而该发明提出的基于最大似然估计和线性插值的方法显著改善了对时钟漂移的预测精度。具体来说,线性模型法的平均误差为0.41ns,线性插值法的平均误差为0.08ns,本发明提出的方法的平局误差为0.01ns。
表1时钟漂移真实值与不同模型的预测值
与现有技术相比,本方法通过基于CIR数据的最大似然估计法,实现对时钟漂移的精准采样,估计精度高,且不受硬件限制,进而通过线性插值法对时钟漂移进行预测和校正,获得具有精确时延测量结果的信道测量数据,可以完全适用于无线信道测量中的时钟漂移校正。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。
Claims (4)
1.一种用于无线信道测量的基于最大似然估计和线性插值的时钟漂移校正方法,其特征在于,包括:
步骤1、时钟漂移采样:在视距传输情况下,周期性地在一已知收发机间距离的位置处采集信道冲激响应(CIR)数据以获取时钟漂移的采样值,具体为:
步骤2、时钟漂移校正:对实测过程中产生的时钟漂移进行校正。
2.根据权利要求1所述的时钟漂移校正方法,其特征是,所述的步骤1,具体包括:
1.1基于CIR数据,粗略估计视距传输径的时延其中:Δτ为测量系统的时域分辨率,所测CIR最大值对应的索引j=argmaxk|h[k]|,h[k]为第k个CIR采样处的路径增益;
1.2采用最大似然估计获取第i次在参考位置采样时的视距传输径的精确估计 其中:测量得到的频域响应/>为离散傅里叶变换,频域采样点fk=f1+(k-1)Δf,f1为所测频段的起始频率,Δf为测量系统的频域分辨率,为搜索区间;
1.3基于所测得的视距传输径时延的精确估计以及根据实际收发机距离所计算得出的视距传输径时延理论值,得出时钟漂移的采样值其中:τTheo为视距传输径时延的理论值。
3.根据权利要求1所述的时钟漂移校正方法,其特征是,所述的步骤2,具体包括:
2.1基于线性插值法,对于任意时刻t测得的信道冲激响应,预测其涉及的时钟漂移其中:/>ti和tj为第i次和第j次在参考位置处采样时的时间,Δτi和Δτj为第i次和第j次采样时的时钟漂移,/>和/>为向下取整和向上取整;
2.2对时刻t处测得的信道冲激响应进行校正其中:⊙为同样大小矩阵间的哈达玛积(Hadamard product),即逐元素乘积,F=[f1,f1+Δf,f1+2Δf,…,f1+(k-1)Δf,…,f1+(K-1)Δf]为包括各个频域采样点对应频率的向量,K为测得的CIR长度,/>为离散傅里叶逆变换。
4.一种实现权利要求1-3中任一所述方法的时钟漂移校正系统,其特征在于,包括:时钟漂移采样单元及时钟漂移校正单元,其中:时钟漂移采样单元根据在固定收发机位置采集的CIR数据,基于最大似然估计法,获得某一时刻的系统时钟漂移值,时钟漂移校正单元根据时钟漂移采样单元记录的时钟漂移值,采用线性插值法,获得系统在任意时刻的时钟漂移预测值,并对采集到的CIR数据进行时钟漂移校正,获得准确的信道测量结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410353200.9A CN118074827A (zh) | 2024-03-26 | 2024-03-26 | 基于最大似然估计和线性插值的时钟漂移校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410353200.9A CN118074827A (zh) | 2024-03-26 | 2024-03-26 | 基于最大似然估计和线性插值的时钟漂移校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN118074827A true CN118074827A (zh) | 2024-05-24 |
Family
ID=91109477
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410353200.9A Pending CN118074827A (zh) | 2024-03-26 | 2024-03-26 | 基于最大似然估计和线性插值的时钟漂移校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN118074827A (zh) |
-
2024
- 2024-03-26 CN CN202410353200.9A patent/CN118074827A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US12081643B2 (en) | Methods for nanosecond-scale time synchronization over a network | |
US11909852B1 (en) | Frequency and gain calibration for time synchronization in a network | |
KR100791179B1 (ko) | 위치 계산 방법 및 위치 산출 장치 | |
US20140355629A1 (en) | Communication apparatus, communication system, and communication method | |
US9768895B2 (en) | Multipath time delay estimation apparatus and method and receiver | |
US10241188B2 (en) | Method and apparatus for obtaining time of arrival TOA when mobile terminal is located | |
JP6292566B2 (ja) | 推定方法およびそれを利用した推定装置 | |
US20240204981A1 (en) | High precision timestamp detection for improved cable modem clock synchronization | |
JP2023549054A (ja) | 到達時間ベースの測距のハイブリッド方法 | |
JP2021520161A (ja) | 位相校正方法及び装置 | |
WO2015180776A1 (en) | Technique for time of arrival estimation | |
Li et al. | Improved, approximate, time-domain ML estimators of chirp signal parameters and their performance analysis | |
CN109682492B (zh) | 基于频域高斯拟合的频率估计方法 | |
CN111077371B (zh) | 一种提高相位测量精度的方法和装置 | |
CN118074827A (zh) | 基于最大似然估计和线性插值的时钟漂移校正方法 | |
WO2019034145A1 (zh) | 一种首径到达时差测量方法和装置 | |
CN117320040B (zh) | 一种基于信道状态信息的载波频偏误差长期动态估计方法 | |
KR101647946B1 (ko) | 측위 장치 및 측위 방법 | |
CN115242367B (zh) | 面向工业无线信道脉冲响应的数据误差矫正方法 | |
US20150195725A1 (en) | Apparatus and Methods for Radio Frequency Ranging | |
CN103259638A (zh) | 一种异地本振一致性误差下的基带时差估计方法 | |
CN117478469B (zh) | 一种载波频率偏移追踪方法及装置 | |
CN119148053B (zh) | 基于已知目标部分信息的定位方法、装置、设备及产品 | |
KR101009529B1 (ko) | 초 광대역 시스템의 거리 측정 방법 및 장치 | |
CN113364541A (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 |