CN109407501B - 一种基于相关信号处理的时间间隔测量方法 - Google Patents
一种基于相关信号处理的时间间隔测量方法 Download PDFInfo
- Publication number
- CN109407501B CN109407501B CN201811583685.1A CN201811583685A CN109407501B CN 109407501 B CN109407501 B CN 109407501B CN 201811583685 A CN201811583685 A CN 201811583685A CN 109407501 B CN109407501 B CN 109407501B
- Authority
- CN
- China
- Prior art keywords
- signal
- sampling
- time interval
- phase
- measured
- 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
- G04—HOROLOGY
- G04F—TIME-INTERVAL MEASURING
- G04F10/00—Apparatus for measuring unknown time intervals by electric means
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Unknown Time Intervals (AREA)
Abstract
本发明公开一种基于相关信号处理的时间间隔测量方法,该方法包括:对第一待测电信号和第二待测电信号进行整形处理,分别得到第一整形信号和第二整形信号;对第一整形信号和第二整形信号进行数学采样处理,分别得到第一采样信号和第二采样信号;对第一采样信号和第二采样信号进行重建处理,分别得到第一重建信号和第二重建信号;对第一重建信号和第二重建信号进行数学相关运算,得到相关函数;对相关函数进行相位密度估计计算,获得相位估计值,其中相位估计值为使得相关函数的值最大的相位值;及基于相位估计值获得时间间隔测量值。本发明的时间间隔测量方法通过时间拉伸与数字信号处理,以实现皮秒及的时间间隔测量精度。
Description
技术领域
本发明涉及时间间隔测量方法。更具体地,涉及一种基于相关信号处理的时间间隔测量方法。
背景技术
时间间隔测量主要用于精确表征两个发生事件之间的时间间隔,是时间计量、测试领域的重要研究问题之一,已广泛应用于现代科学技术等多个领域,包括精密时间频率传递、雷达、无线电导航定位、通信、激光测距、光子物理等。
时间间隔测量一般将两个事件转换为方便处理的两个电脉冲信号,通过逻辑门或者模拟电路等对两个电脉冲信号进行特定处理后得到两个电脉冲之间的时间差。常见的时间间隔测量方法包括:电子计数法、扩展法,模拟内插法、延迟线内插法、抽头延迟线法、查分延迟线法、游标法、时间-幅度转化法和时间数字转换法等,其都是通过电路对输入的电脉冲信号进行直接延迟、锁存和转换等纯硬件处理方法,系统测量精度有限,不能满足高精度的要求。
目前,普通的时间间隔测量仪器采用脉冲填充法,成本较低,但测量误差还在纳秒量级,不能满足激光测距、卫星导航定位、粒子飞行探测、频率基准等方面的要求。而高精度的测量方法中,如基于模拟时间扩展的计数法、基于AD变换器的模拟时间-幅度转换法、基于延迟线的时间-数字变换器(TDC)法和基于冲击振荡器的频率游标法等,测量分辨率都已达到了皮秒良机,但其复杂的电路设计和昂贵的造价却限制了其应用。
因此,需要提供一种既能满足时间间隔测量精度要求又能降低电路设计复杂度和造价的时间间隔测量方法。
发明内容
本发明的目的在于提供一种基于相关信号处理的时间间隔测量方法,通过时间拉伸与数字信号处理,以实现皮秒及的时间间隔测量精度,该时间间隔测量方法既能满足时间间隔测量精度要求又能降低电路设计复杂度和造价。
为达到上述目的,本发明采用下述技术方案:
一种基于相关信号处理的时间间隔测量方法,该方法包括:
对第一待测电信号和第二待测电信号进行整形处理,分别得到第一整形信号和第二整形信号;
对第一整形信号和第二整形信号进行数学采样处理,分别得到第一采样信号和第二采样信号;
对第一采样信号和第二采样信号进行重建处理,分别得到第一重建信号和第二重建信号;
对第一重建信号和第二重建信号进行数学相关运算,得到相关函数;
对相关函数进行相位密度估计计算,获得相位估计值,其中相位估计值为使得相关函数的值最大的相位值;及
基于相位估计值获得时间间隔测量值。
优选地,对第一待测电信号和第二待测电信号进行整形处理包括:
基于公式s1(t)=s(t)对第一待测信号s(t)进行整形处理,得到第一整形信号s1(t);
基于公式s2(t)=s(t-Θ)对第二待测信号s(t-Θ)进行整形处理,得到第二整形信号s2(t);
其中,第一整形信号和第二整形信号具有相同的幅度和持续时间,Θ为第一待测信号和第二待测信号之间的实际时间间隔。
进一步优选地,对第一整形信号和第二整形信号进行数学采样处理包括:
基于公式x1(nTs)=s1(nTs)+w1(nTs)对第一整形信号s1(t)进行采样处理,得到第一采样信号x1(nTs);
基于公式x2(nTs)=s2(nTs)+w2(nTs)对第二整形信号s2(t)进行采样处理,得到第二采样信号x2(nTs);
其中,采样处理的采样频率为fs,Ts为采样间隔,w1(nTs)为第一白噪声信号,w2(nTs)为第二白噪声信号。
进一步优选地,对第一采样信号和第二采样信号进行重建处理包括:
其中,W为第一待测电信号s(t)和第二待测电信号s(t-Θ)的信号带宽,sinc(x)=sin(x)/x,f0为第一待测电信号s1(t)和第二待测电信号s2(t)的中心频率。
进一步优选地,对第一重建信号和第二重建信号进行数学相关运算包括:
进一步优选地,对相关函数进行相位密度估计计算包括:
进一步优选地,基于相位估计值获得时间间隔测量值包括:
进一步优选地,基于公式对第一待测电信号s(t)和第二待测电信号s(t-Θ)的频率响应的频率区域进行限制,其中M为非负整数,fs为采样频率,第一待测电信号s(t)和第二待测电信号s(t-Θ)的带宽小于fs/2且其中心频率在f(2M+1)s/4附近。
本发明的有益效果如下:
本发明中一种基于相关信号处理的时间间隔测量方法,通过时间拉伸与数字信号处理,能够实现皮秒及的时间间隔测量精度,既能满足时间间隔测量精度要求又能降低电路设计复杂度和造价,同时避免了系统固有时间间隔测量偏移和分辨力下降的隐患。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明。
图1示出一种基于相关信号处理的时间间隔测量方法的方法流程图。
图2示出基于相关信号处理的时间间隔测量方法原理图。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
本发明中提供一种基于相关信号处理的时间间隔测量方法,通过巧妙的时间拉伸与数字信号处理,能够实现皮秒级的时间间隔测量精度。基于相关信号处理的理论的时间间隔测量方法主要是利用两个待测电脉冲信号(或其转换形式)之间独有的互相关特性,通过数学算法进行相关信号处理运算,得到两个待测电脉冲信号之间的相位差,此相位差即为两脉冲信号之间的精确时间间隔。
如图1所述,该方法包括:
对第一待测电信号和第二待测电信号进行整形处理,分别得到第一整形信号和第二整形信号;
对第一整形信号和第二整形信号进行数学采样处理,分别得到第一采样信号和第二采样信号;
对第一采样信号和第二采样信号进行重建处理,分别得到第一重建信号和第二重建信号;
对第一重建信号和第二重建信号进行数学相关运算,得到相关函数;
对相关函数进行相位密度估计计算,获得相位估计值,其中相位估计值为使得相关函数的值最大的相位值;及
基于相位估计值获得时间间隔测量值。
下面结合图2对该方法进行详细说明。
第一步两待测电脉冲信号整形模型
本方法基于系统通过特定电路已对两待测电信号整形拉伸为幅度和持续时间完全相同的待处理信号,分别用s1(t)和s2(t)表示如下,其中s(t)等同于为第一路整形后的信号,Θ表示两个时间之间的真实时间间隔。
s1(t)=s(t) (1)
s2(t)=s(t-Θ) (2)
第二步信号采样
对第一步输出的两路信号s1(t)和s2(t)进行数学采样,通过频率fs采样后的信号为x1(nTs)和x2(nTs):
x1(nTs)=s1(nTs)+w1(nTs) (3)
x2(nTs)=s2(nTs)+w2(nTs) (4)
式中Ts为采样间隔,w1(nTs)和w2(nTs)为白噪声信号。
第三步信号重构
其中W为s1(t)和s2(t)信号带宽,sinc(x)=sin(x)/x;f0为s1(t)和s2(t)信号中心频率。
第四步重建信号的相关信号运算
第五步相位精密估计计算
对第四步相关计算后的数据进行相位提取,其估计值为使得相关函数值最大的相位值:
第六步时间间隔测量值获取
另外,在一种基于相关信号处理的时间间隔测量方法中,为保证被采样信号的频谱不混叠,需要整形信号进行带宽限制,s1(t)和s2(t)频率响应特性应该限制在式(9)所示的频率区域。其中M是一个非负整数,fs为采样频率,s1(t)和s2(t)信号的带宽应小于fs/2,并且中心频率应当在f(2M+1)s/4附近。
如此才能保证完全正确的恢复原始输入信号。
同时,针对系统设计的不同中心频率和带宽的整形信号,系统采样率需要通过调整实验最终才能确定,并且采样频率优选具有小数频率的特殊数字,避免由于整数倍采样造成采样信号相位不能遍历原始信号相位的问题,造成系统固有时间间隔测量偏移和分辨力下降的隐患。
应注意的是,本发明的说明书和权利要求书及上述附图中的属于“第一”、“第二”等是用于区别不同的对象,而不是用于描述特定顺序。此外,术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法或设备固有的气体步骤或单元。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (4)
1.一种基于相关信号处理的时间间隔测量方法,其特征在于,所述方法包括:
对第一待测电信号和第二待测电信号进行整形处理,分别得到第一整形信号和第二整形信号;
对所述第一整形信号和所述第二整形信号进行数学采样处理,分别得到第一采样信号和第二采样信号;
对所述第一采样信号和所述第二采样信号进行重建处理,分别得到第一重建信号和第二重建信号;
对所述第一重建信号和所述第二重建信号进行数学相关运算,得到相关函数;
对所述相关函数进行相位密度估计计算,获得相位估计值,其中所述相位估计值为使得所述相关函数的值最大的相位值;及
基于所述相位估计值获得时间间隔测量值;
所述对第一待测电信号和第二待测电信号进行整形处理包括:
基于公式s1(t)=s(t)对所述第一待测电信号s(t)进行整形处理,得到所述第一整形信号s1(t);
基于公式s2(t)=s(t-Θ)对所述第二待测电信号s(t-Θ)进行整形处理,得到所述第二整形信号s2(t);
其中,所述第一整形信号和所述第二整形信号具有相同的幅度和持续时间,Θ为所述第一待测电信号和所述第二待测电信号之间的实际时间间隔;
所述对所述第一整形信号和所述第二整形信号进行数学采样处理包括:
基于公式x1(nTs)=s1(nTs)+w1(nTs)对所述第一整形信号s1(t)进行采样处理,得到所述第一采样信号x1(nTs);
基于公式x2(nTs)=s2(nTs)+w2(nTs)对所述第二整形信号s2(t)进行采样处理,得到所述第二采样信号x2(nTs);
其中,所述采样处理的采样频率为fs,Ts为采样间隔,w1(nTs)为第一白噪声信号,w2(nTs)为第二白噪声信号;
所述对所述第一采样信号和所述第二采样信号进行重建处理包括:
其中,W为所述第一待测电信号s(t)和所述第二待测电信号s(t-Θ)的信号带宽,sin c(x)=sin(x)/x,f0为所述第一待测电信号s1(t)和所述第二待测电信号s2(t)的中心频率;
所述对所述第一重建信号和所述第二重建信号进行数学相关运算包括:
所述对所述相关函数进行相位密度估计计算包括:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811583685.1A CN109407501B (zh) | 2018-12-24 | 2018-12-24 | 一种基于相关信号处理的时间间隔测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811583685.1A CN109407501B (zh) | 2018-12-24 | 2018-12-24 | 一种基于相关信号处理的时间间隔测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109407501A CN109407501A (zh) | 2019-03-01 |
CN109407501B true CN109407501B (zh) | 2020-10-27 |
Family
ID=65460787
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811583685.1A Active CN109407501B (zh) | 2018-12-24 | 2018-12-24 | 一种基于相关信号处理的时间间隔测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109407501B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110865531B (zh) * | 2019-11-20 | 2022-08-19 | 北京无线电计量测试研究所 | 基于非线性回归的时间间隔测量方法及系统 |
CN111010235B (zh) * | 2019-12-19 | 2021-07-06 | 北京无线电计量测试研究所 | 一种收发机和光学时频传递装置 |
CN114637184B (zh) * | 2022-03-14 | 2023-12-08 | 浙江师范大学 | 一种基于相关信号处理的时间间隔快速测量方法 |
CN116736680B (zh) * | 2023-08-16 | 2023-11-28 | 中国科学技术大学 | 一种带通信号的高精度时间信息提取方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1114773A (ja) * | 1997-06-24 | 1999-01-22 | Nec Corp | タイムインターバル計測方式及び距離測定装置 |
DE102004057465A1 (de) * | 2003-12-23 | 2005-07-28 | Riegl Laser Measurement Systems Gmbh | Verfahren zur Messung von Zeitintervallen |
CN1862234A (zh) * | 2006-04-18 | 2006-11-15 | 西安电子科技大学 | 一种基于时-空关系的时间间隔与频率的测量方法 |
CN202372798U (zh) * | 2011-11-16 | 2012-08-08 | 贵州航天计量测试技术研究所 | 一种ns级时间间隔测量装置 |
CN102830610A (zh) * | 2012-09-17 | 2012-12-19 | 中国电子科技集团公司第四十一研究所 | 一种基于tdc-gp21的精密时间自动测量电路及方法 |
CN103401557A (zh) * | 2013-08-12 | 2013-11-20 | 龙芯中科技术有限公司 | 时间数字转换器和时间间隔测量方法 |
CN103676622A (zh) * | 2013-10-28 | 2014-03-26 | 中国电子科技集团公司第四十一研究所 | 一种高精度的正负时间间隔测量方法及装置 |
-
2018
- 2018-12-24 CN CN201811583685.1A patent/CN109407501B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1114773A (ja) * | 1997-06-24 | 1999-01-22 | Nec Corp | タイムインターバル計測方式及び距離測定装置 |
DE102004057465A1 (de) * | 2003-12-23 | 2005-07-28 | Riegl Laser Measurement Systems Gmbh | Verfahren zur Messung von Zeitintervallen |
CN1862234A (zh) * | 2006-04-18 | 2006-11-15 | 西安电子科技大学 | 一种基于时-空关系的时间间隔与频率的测量方法 |
CN202372798U (zh) * | 2011-11-16 | 2012-08-08 | 贵州航天计量测试技术研究所 | 一种ns级时间间隔测量装置 |
CN102830610A (zh) * | 2012-09-17 | 2012-12-19 | 中国电子科技集团公司第四十一研究所 | 一种基于tdc-gp21的精密时间自动测量电路及方法 |
CN103401557A (zh) * | 2013-08-12 | 2013-11-20 | 龙芯中科技术有限公司 | 时间数字转换器和时间间隔测量方法 |
CN103676622A (zh) * | 2013-10-28 | 2014-03-26 | 中国电子科技集团公司第四十一研究所 | 一种高精度的正负时间间隔测量方法及装置 |
Non-Patent Citations (4)
Title |
---|
Accurate Time Interval Measurement Method Based on Vernier Caliper Principle;Yuanyuan Song,等;《2010 International Conference on Image Analysis and Signal Processing》;20100601;553-555 * |
基于相位处理的时间间隔测量系统设计;段亚东;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20170315(第3期);C042-3499 * |
脉冲激光测距中高精度时间间隔的测量;宋建辉,等;《光学精密工程》;20090531;第17卷(第5期);1046-1050 * |
高精度时间间隔测量系统;邢燕;《电子测量技术》;20100531;第33卷(第5期);1-3,9 * |
Also Published As
Publication number | Publication date |
---|---|
CN109407501A (zh) | 2019-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109407501B (zh) | 一种基于相关信号处理的时间间隔测量方法 | |
CN104007316B (zh) | 一种欠采样速率下的高精度频率测量方法及其测量仪 | |
CN103941087B (zh) | 欠采样速率下的高频余弦信号的频率测量方法及其装置 | |
CN105656485B (zh) | 一种多通道时间交错adc测量校准方法和装置 | |
CN103837740A (zh) | 一种高精度数字瞬时测频方法及装置 | |
CN110208601B (zh) | 基于fpga的瞬时测频方法及其数字接收机 | |
CN105629061A (zh) | 一种基于高稳定度宽基准脉冲的精密频率测量装置 | |
CN105182069A (zh) | 一种异频架构下的高分辨群量子化相位处理方法 | |
CN103983849A (zh) | 一种实时高精度的电力谐波分析方法 | |
Chen et al. | Robust precise time difference estimation based on digital zero-crossing detection algorithm | |
Su et al. | Digital Instantaneous Frequency Measurement of a Real Sinusoid Based on Three Sub‐Nyquist Sampling Channels | |
CN112964931B (zh) | 基于双通道欠采样的非理想多阻尼谐波信号参数测量方法 | |
US10924130B1 (en) | High resolution digital trigger detector | |
CN103941280B (zh) | 基于冲激响应不变法的数字核脉冲高斯成形方法 | |
CN105680858A (zh) | 一种估计tiadc并行采集系统时间偏移误差的方法 | |
CN104991119A (zh) | 一种消除伪峰、谱泄漏效应的互素谱分析方法及其装置 | |
CN116781079A (zh) | 一种基于参考通道的tiadc时间失配误差校准电路 | |
CN116360235A (zh) | 一种基于SerDes的TDC实现装置 | |
CN114047683B (zh) | 一种基于正交采样内插的时间间隔测量方法和装置 | |
CN110808929A (zh) | 相减策略的实复转换式信噪比估计算法 | |
Zhu et al. | A high-precision time interval measurement method using phase-estimation algorithm | |
CN109581518A (zh) | 一种无需同步电流采集传输和频域测量sip数据采集和处理方法 | |
CN112129983B (zh) | 一种基于等时间间隔等效取样的波形恢复数据处理方法 | |
CN115361108B (zh) | 一种超宽带高精度接收时间戳获取方法及装置 | |
CN104618288B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |