CN105631149A - 一种基于相关性求时差的方法 - Google Patents

一种基于相关性求时差的方法 Download PDF

Info

Publication number
CN105631149A
CN105631149A CN201610001516.7A CN201610001516A CN105631149A CN 105631149 A CN105631149 A CN 105631149A CN 201610001516 A CN201610001516 A CN 201610001516A CN 105631149 A CN105631149 A CN 105631149A
Authority
CN
China
Prior art keywords
time difference
signal
signals
correlation
dependency
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
Application number
CN201610001516.7A
Other languages
English (en)
Other versions
CN105631149B (zh
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.)
Southwest Minzu University
Original Assignee
Southwest Minzu 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 Southwest Minzu University filed Critical Southwest Minzu University
Priority to CN201610001516.7A priority Critical patent/CN105631149B/zh
Publication of CN105631149A publication Critical patent/CN105631149A/zh
Application granted granted Critical
Publication of CN105631149B publication Critical patent/CN105631149B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/12—Timing analysis or timing optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Unknown Time Intervals (AREA)

Abstract

本发明提出一种基于相关性求时差的方法。传统的基于相关性求时差的方法是采用互相关函数计算两信号的相关性，得到的时差的准确度与信号采样率有关，采样率越高，所得时差的准确度就越高，然而过高的采样率会增加系统的开销。针对该问题，本发明使两信号的采样周期不相等且均能被它们之间的差值整除，并采用新的式子计算两信号的相关性和时差。采用本发明的方法，可以在不增加采样率的条件下提高所求时差的准确度。

Description

一种基于相关性求时差的方法
技术领域
本发明属于信号处理领域，特别涉及一种基于相关性求时差的方法。
背景技术
互相关函数是一种计算两信号相关性的常用方法，其计算式子为，其中S1和S2为两信号的采样序列，m和n均为整数。互相关函数常被用来分析两信号的时差，如：文献《基于信号互相关的无源时差定位技术仿真分析》（《电子信息对抗技术》2010年第2期，P49-54）在多站无源定位中，利用互相关函数提取各站之间的信号时差；文献《互相关时差分析技术及其应用——以胜利油田三维地震资料连片处理为例》（《石油物探》2010年第1期，P23-29）在三维地震资料连片处理中，利用互相关函数定量描述两个地震信号的时间延迟；文献《基于时间互相关的超声测距信号获取方法》（《仪表技术与传感器》2014年第6期，P126,127,130）在超声测距中，利用互相关函数得到超声发射与接收信号间的时差；文献《数字式时差法超声流量计的设计与实现》（《自动化仪表》2014年第9期，P80-83）在超声流量计的设计与实现中，利用互相关函数计算顺逆流传播时间差。
基于互相关函数得到的两信号时差的准确度与信号采样率有关，采样率越高，所得时差的准确度就越高。然而过高的采样率会增加系统的开销。
发明内容
为了提高两信号时差计算的准确度，与目前采用增加采样率的方法不同，本发明内容是使这两信号的采样周期不相等且均能被它们之间的差值整除，并采用新的式子计算这两信号的相关性和时差。
本发明的技术方案为：
步骤1：使两信号的采样周期分别为T1和T2，T1不等于T2，且T1与T2之间差值的绝对值△T能够整除T1和T2。对这两信号进行采样，分别得到采样序列S1和S2。
步骤2：通过下式计算采样序列S1和S2的相关性：
其中，m和n均为整数。
步骤3：M为得到的所有相关性值中的最大相关峰所对应的m值，通过下式计算这两信号时差：
该时差值的正和负分别表示采样序列S1早于和晚于S2。
本发明具有如下技术效果或优点：在基于相关性求时差的方法中，本发明可以降低因采样间隔而带来的误差，在理论上，最大地可使该项误差由T1/2（或T2/2）降到为△T/2。
具体实施方式
下面结合具体实施例来描述本发明。
设两信号中的信号一为sin(10*2π*t)，信号二为sin(10*2π*(t+0.05))，即两信号的频率均为10赫兹，信号一的初始时刻为0秒，信号二的初始时刻为0.05秒，信号一与信号二的时间差为-0.05秒。采用本发明计算这两信号时差的过程如下：
步骤1：对这两信号分别设置1/30秒和1/25秒这两个不同的采样周期，其中这两个采样周期均能被它们差值的绝对值1/150（即1/25-1/30）秒整除。对这两信号采样，采样时长均为1秒，得到信号一的采样序列S1如下表所示（由于数据量较大，省略号处数据未显示）：
信号二的采样序列S1如下表所示（由于数据量较大，省略号处数据未显示）：
步骤2：依据本说明书的“发明内容”中的“步骤2”的式子进行计算，其中S1和S2分别为这两信号的采样序列，T1=1/30秒，T2=1/25秒，△T=1/150秒，m和n为整数。计算结果就是这两信号的相关性值，如下表所示（由于数据量较大，省略号处数据已未显示）：
步骤3：得到在所有相关性值中的最大相关峰所对应的m值为-8，计算两信号时差为-8*(T2-T1)=-8*(1/25-1/30)=-0.0533秒。
该值与两信号真实时间差-0.05秒对比，误差为0.0033秒，该误差优于采用传统互相关函数求两信号时差时的误差（当对这两信号用1/30秒的周期进行采样时，采用传统互相关函数所求得的两信号时差的误差为0.0167秒）。

Claims (3)

1.一种基于相关性求时差的方法，包括如下步骤：
步骤1：对两信号采样；
步骤2：计算两信号采样序列的相关性；
步骤3：根据相关性计算两信号时差；
其特征在于，在步骤1中，两信号的采样周期不相等，而且这两信号的采样周期均能被这两信号的采样周期的差值整除。
2.如权利要求1所述的一种基于相关性求时差的方法，其特征是，所述步骤2采用式子计算两信号的相关性，其中S1和S2分别表示两信号的采样序列，T1和T2分别表示两信号的采样周期，△T为T1与T2的差值的绝对值，m和n均为整数。
3.如权利要求2所述的一种基于相关性求时差的方法，其特征是，所述步骤3采用式子计算两信号时差，其中M为最大相关峰对应的m值。
CN201610001516.7A 2016-01-05 2016-01-05 一种基于相关性求时差的方法 Expired - Fee Related CN105631149B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610001516.7A CN105631149B (zh) 2016-01-05 2016-01-05 一种基于相关性求时差的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610001516.7A CN105631149B (zh) 2016-01-05 2016-01-05 一种基于相关性求时差的方法

Publications (2)

Publication Number Publication Date
CN105631149A true CN105631149A (zh) 2016-06-01
CN105631149B CN105631149B (zh) 2019-05-03

Family

ID=56046078

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610001516.7A Expired - Fee Related CN105631149B (zh) 2016-01-05 2016-01-05 一种基于相关性求时差的方法

Country Status (1)

Country Link
CN (1) CN105631149B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107621638A (zh) * 2017-08-01 2018-01-23 昆明理工大学 一种基于平移信号峰值求两个脉冲信号间的时间差的方法
CN108809438A (zh) * 2017-04-27 2018-11-13 深圳市中兴微电子技术有限公司 一种时差补偿方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090174976A1 (en) * 2006-02-28 2009-07-09 Kabushiki Kaisha Toshiba Current differential relay device, signal processing method thereof and power transmission line protection system
CN103592619A (zh) * 2013-10-15 2014-02-19 清华大学 基于单路欠采样的时频差估计方法
CN104605888A (zh) * 2015-01-22 2015-05-13 飞依诺科技(苏州)有限公司 超声信号的时延估计方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090174976A1 (en) * 2006-02-28 2009-07-09 Kabushiki Kaisha Toshiba Current differential relay device, signal processing method thereof and power transmission line protection system
CN103592619A (zh) * 2013-10-15 2014-02-19 清华大学 基于单路欠采样的时频差估计方法
CN104605888A (zh) * 2015-01-22 2015-05-13 飞依诺科技(苏州)有限公司 超声信号的时延估计方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
何先龙等: "基于多重互相关函数计算地探信号到时差", 《地球物理学进展》 *
陈慧等: "基于信号互相关的无源时差定位技术仿真分析", 《电子信息对抗技术》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108809438A (zh) * 2017-04-27 2018-11-13 深圳市中兴微电子技术有限公司 一种时差补偿方法和装置
CN108809438B (zh) * 2017-04-27 2020-03-06 深圳市中兴微电子技术有限公司 一种时差补偿方法和装置
CN107621638A (zh) * 2017-08-01 2018-01-23 昆明理工大学 一种基于平移信号峰值求两个脉冲信号间的时间差的方法
CN107621638B (zh) * 2017-08-01 2020-11-17 昆明理工大学 一种基于平移信号峰值求两个脉冲信号间的时间差的方法

Also Published As

Publication number Publication date
CN105631149B (zh) 2019-05-03

Similar Documents

Publication Publication Date Title
CN107907878B (zh) 高精度获取fmcw雷达距离测量值的方法
CN102565857B (zh) 自动剩余动校正方法
CN108535773B (zh) 一种微地震事件检测方法及系统
CN104991440A (zh) 高精度irig-b(ac)码解调方法及装置
CN105307095B (zh) 一种基于fft的高分辨率音频频率测量方法
CN102508031A (zh) 一种基于傅里叶级数的局部放电脉冲相角测量方法
CN103852785B (zh) 地层各向异性的评价方法
CN105677958B (zh) 一种基于互相关函数求时差的方法
CN105631149A (zh) 一种基于相关性求时差的方法
CN103543331B (zh) 一种计算电信号谐波和间谐波的方法
CN106330342A (zh) 一种低计算复杂度的水声通信多普勒因子估计方法
CN110048795A (zh) 一种地震仪采集数据时钟的方法及装置
CN109283581A (zh) 一种基于深度域地震波频散分析的储层含气性评价方法
CN106324681B (zh) 一种时频连续小波变换的井震标定方法
CN109856672B (zh) 基于深度波数谱的瞬变波包提取方法、存储介质与终端
CN104111481A (zh) 同步时钟相位差测量系统和方法
CN103821499A (zh) 用于油井动液面深度检测的声音信号处理方法
CN107064942B (zh) 一种高精度多普勒信号模拟方法
CN108490490B (zh) 裂缝带表征方法、装置、电子设备及计算机存储介质
CN105158799B (zh) 基于相位准则的物理探测回波信号识别与提取方法
CN105917244B (zh) 将脉冲源定位在色散介质中的方法
CN104849551B (zh) 一种谐相角分析方法
KR101446439B1 (ko) Fmcw 레이더의 고정밀 주파수 추정 방법 및 고정밀 주파수 추정 장치
CN107036706B (zh) 一种套管振动井口监听检测设备
CN108777582B (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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190503

Termination date: 20220105

CF01 Termination of patent right due to non-payment of annual fee