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

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

Info

Publication number
CN105631149B
CN105631149B CN201610001516.7A CN201610001516A CN105631149B CN 105631149 B CN105631149 B CN 105631149B CN 201610001516 A CN201610001516 A CN 201610001516A CN 105631149 B CN105631149 B CN 105631149B
Authority
CN
China
Prior art keywords
time difference
signals
correlation
sampling
accuracy
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.)
Expired - Fee Related
Application number
CN201610001516.7A
Other languages
English (en)
Other versions
CN105631149A (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

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

Claims (1)

1.一种基于相关性求时差的方法，包括如下步骤：
步骤1：对两信号采样；
步骤2：计算两信号采样序列的相关性；
步骤3：根据相关性计算两信号时差；
其特征在于，在步骤1中，两信号的采样周期不相等，而且这两信号的采样周期均能被这两信号的采样周期的差值整除；步骤2采用式子计算两信号的相关性，其中S 1和S 2分别表示两信号的采样序列，T 1和T 2分别表示两信号的采样周期，△T为T 1与T 2的差值的绝对值，m和n均为整数；步骤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 CN105631149A (zh) 2016-06-01
CN105631149B true 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)

Families Citing this family (2)

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

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103592619A (zh) * 2013-10-15 2014-02-19 清华大学 基于单路欠采样的时频差估计方法
CN104605888A (zh) * 2015-01-22 2015-05-13 飞依诺科技(苏州)有限公司 超声信号的时延估计方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4908018B2 (ja) * 2006-02-28 2012-04-04 株式会社東芝 電流差動継電装置とその信号処理方法、および送電線保護システム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
基于信号互相关的无源时差定位技术仿真分析;陈慧等;《电子信息对抗技术》;20100330;第25卷(第2期);全文
基于多重互相关函数计算地探信号到时差;何先龙等;《地球物理学进展》;20090831;第24卷(第4期);全文

Also Published As

Publication number Publication date
CN105631149A (zh) 2016-06-01

Similar Documents

Publication Publication Date Title
CN111505716B (zh) 一种基于时间同步抽取广义Chirplet变换的地震时频分析方法
CN107607934B (zh) 一种时差、频差、频差变化率联合估计方法
CN113406698B (zh) 一种基于纵横波解耦的双相介质弹性波逆时偏移成像方法
CN105783974B (zh) 一种线性调频信号的检测、参数估计方法及系统
WO2017167191A1 (zh) 地震数据处理方法和装置
CN109188362B (zh) 一种麦克风阵列声源定位信号处理方法
CN103018555B (zh) 一种高精度的电力参数软件同步采样方法
CN109884694B (zh) 一种基于挤压加窗傅里叶变换的高铁震源地震信号时频分析方法
CN105588883B (zh) 三维岩石力学参数获取方法和系统
CN103675850B (zh) 一种基于分段分裂基傅里叶变换的双频gps信号捕获方法
CN105631149B (zh) 一种基于相关性求时差的方法
CN102298097A (zh) 一种估计云闪雷电脉冲信号到达时间差的方法
CN110140326A (zh) 对设备的性能数据进行采样的方法和装置
CN104614764B (zh) 一种基于扰动弹性阻抗反演的裂缝储层流体识别方法
CN105677958B (zh) 一种基于互相关函数求时差的方法
CN109117816A (zh) 基于六阶样条插值小波的信号奇异点检测方法
CN110247867B (zh) 水声多普勒估计方法及装置、水声通信方法及系统
CN104456089B (zh) 一种流体管道多泄漏点定位方法
CN111060961A (zh) 基于多信息约束反演的品质因子确定方法、装置及系统
CN104597491A (zh) 一种基于方位弹性阻抗差异的裂缝储层参数反演方法
CN107783135B (zh) 一种三元矢量阵被动测距方法
CN103175986A (zh) 一种油气水三相流气液相流速测量方法
CN106330342A (zh) 一种低计算复杂度的水声通信多普勒因子估计方法
CN115113243A (zh) 一种联合频率插值估计的导航信号捕获方法
CN108398664A (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190503

Termination date: 20220105