CN107576964B - 线性变频信号的回波时间测量方法 - Google Patents

线性变频信号的回波时间测量方法 Download PDF

Info

Publication number
CN107576964B
CN107576964B CN201710740712.0A CN201710740712A CN107576964B CN 107576964 B CN107576964 B CN 107576964B CN 201710740712 A CN201710740712 A CN 201710740712A CN 107576964 B CN107576964 B CN 107576964B
Authority
CN
China
Prior art keywords
echo
pulse
time
frequency
width
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
Application number
CN201710740712.0A
Other languages
English (en)
Other versions
CN107576964A (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.)
Xian University of Technology
Original Assignee
Xian University of Technology
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 Xian University of Technology filed Critical Xian University of Technology
Priority to CN201710740712.0A priority Critical patent/CN107576964B/zh
Publication of CN107576964A publication Critical patent/CN107576964A/zh
Application granted granted Critical
Publication of CN107576964B publication Critical patent/CN107576964B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明线性变频信号的回波时间测量方法,包括以下步骤:步骤1,超声波传感器发射周期先线性增加后线性减小的脉冲串,记录脉冲串的频率转折时刻t0;步骤2,经过延时以消除超声波传感器的余振影响;步骤3,超声波传感器同时接收超声波回波信号,经过放大、滤波电路处理后得到有效回波信号;步骤4,将有效回波信号经过模拟电压比较器整形为数字回波脉冲序列,并传送给处理器,计算出各个回波脉冲的宽度;步骤5,找到回波脉冲的宽度变化转折点m,即为回波脉冲频率转折点,记录该点的时刻为t[m];步骤6,计算回波时间techo,techo=t[m]‑t0。本发明解决了现有技术中存在的测量回波信号精度低、实时性差、计算量大的问题。

Description

线性变频信号的回波时间测量方法
技术领域
本发明属于超声波测距技术领域,涉及一种线性变频信号的回波时间测量方法。
背景技术
超声波脉冲回波测距技术是目前广泛应用的非接触测距方法[1],目前在脉冲回波测距技术中有阀门法、互相关函数法、频域相位法和时域法[2]、FSK(频移键控)等检测方法[3];由于超声波回波信号传输过程衰减严重,回波前沿难检测,阀门法精度无法保证[2],互相关函数法、频域相位法和时域方法需要对回波信号采样,需要高速的模数转换器和大量存储器,分析方法计算量大,实用性较低[2],FSK(频移键控)是检测超声波频率发生跳变的时刻来确定超声波在介质中的渡越时间的方法,由于超声波传感器在发射和接收过程需要较长的起振时间[5],而脉冲回波法发送脉冲个数有限,在较短时间内频率不能稳定,回波频率切换点不易找到。
[1]潘仲明简盈王跃科.基于两步相关法的大量程超声波测距技术[J].电子测量与仪表学报2006 20(5):73-76。
[2]程晓畅苏绍景王跃科潘仲明祝琴.超声回波信号调制及其包络相关时延估计算法[J].传感技术学报2006 19(6)2571-2577。
[3]孟升卫曹义尹洪涛.基于FSK的超声波测距技术研究[J].电子测量技术200932(5):150-156。
[4]徐斌.一种考虑起振延迟的低频超声波短距高精度测量校准方法[J].传感技术学报2013 26(5):666-669。
[5]章小兵宋爱国唐鸿儒.基于单接收头的超声波多目标测距[J].传感技术学报2007 20(5):1167-1170。
发明内容
本发明的目的是提供一种线性变频信号的回波时间测量方法,解决了现有技术中存在的测量回波信号精度低、实时性差、计算量大的问题。
本发明所采用的技术方案是,线性变频的回波时间测量方法,包括以下步骤:
步骤1,超声波传感器发射周期先线性增加后线性减小的脉冲串,记录发射脉冲串的频率转折时刻t0
步骤2,经过延时以消除超声波传感器的余振影响;
步骤3,超声波传感器在发射脉冲串的同时接收超声波回波信号,超声波回波信号经过放大、滤波电路处理后得到有效回波信号;
步骤4,将有效回波信号经过模拟电压比较器整形为数字回波脉冲序列,并传送给处理器,计算出各个回波脉冲的宽度;
步骤5,找到回波脉冲的宽度变化转折点m,即为回波脉冲频率转折点,记录该点的时刻为t[m];
步骤6,计算回波时间techo,techo=t[m]-t0
步骤1具体为:
超声波传感器发送周期以
Figure BDA0001389049810000021
为起始,周期渐变的步长为Δt的,周期先线性增加后线性减小的脉冲串,当发送宽度为
Figure BDA0001389049810000022
的脉冲时,记录当时的时刻t0,即为发射脉冲串的频率转折时刻;
其中,fs为超声波传感器的中心频率,n为脉冲串中脉冲的个数。
步骤1中,n取值为奇数,n的取值范围为7-15。
步骤4具体为:
按照回波脉冲的次序记下每个回波脉冲的上升边沿时刻,并记录到上升边沿时刻数组t[i]中,i为回波脉冲的序数;
根据t[i]确定各个回波脉冲的宽度,记录到回波脉冲宽度数组w[i]中,w[i]=t[i]-t[i-1],i为回波脉冲的序数。
步骤5具体为:
根据回波脉冲宽度数组w[i],找到回波脉冲宽度变化转折点m:即在接收到的回波脉冲序列中,存在一段包含n个脉冲的连续脉冲串,回波脉冲满足先连续递增
Figure BDA0001389049810000031
个脉冲达到m点,再从m点连续递减
Figure BDA0001389049810000032
个脉冲的规则,将其中心位置m点确定为回波脉冲频率转折点,记录该点的时刻为t[m]。
本发明的有益效果是,线性变频的回波时间测量方法,通过在超声波传感器中心频率附近,发送逐脉冲频率先线性递减再线性递增的脉冲串,有助于传感器跟随激励信号变化,回波信号输出比较稳定;再将回波信号幅度放大并整形成数字脉冲序列,测量每个脉冲的宽度(频率),以确定回波信号的频率转折点,发射信号的频率转折点和回波信号的频率转折点之时间差为有效的回波时间;本发明不需要采样模拟回波信号,降低了系统对模数转换器、存储器等超声波检测距装置的硬件要求,成本较低,不需要复杂的算法,可靠性强,容易实现,具有较高的实用性。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
线性变频的回波时间测量方法,包括以下步骤:
步骤1,超声波传感器发射周期先线性增加后线性减小的脉冲串,记录发射脉冲串的频率转折时刻t0
步骤1具体为:
超声波传感器发送周期以
Figure BDA0001389049810000041
为起始,周期渐变的步长为Δt的,周期先线性增加后线性减小的脉冲串,即脉冲串中的每个脉冲的宽度分别为
Figure BDA0001389049810000042
Figure BDA0001389049810000043
当发送宽度为
Figure BDA0001389049810000044
的脉冲时,即脉冲串中心位置的脉冲,记录当时的时刻t0,即为发射脉冲串的频率转折时刻;
其中,fs为超声波传感器的中心频率,n为脉冲串中脉冲的个数,n取值为奇数,n的取值范围为7-15,Δt可根据超声波传感器的频带宽度取
Figure BDA0001389049810000045
的0.1-0.25%;
发送这种周期(频率)渐变的脉冲串,便于超声波传感器跟随激励信号变化,其超声波回波信号输出比较稳定。
步骤2,经过一段时间的延时(大约毫秒级),以消除超声波传感器的余振影响,延时时间依据具体的超声波传感器的特性而定。
步骤3,超声波传感器在发射脉冲串的同时接收超声波回波信号,超声波回波信号经过放大、滤波电路处理后得到有效回波信号,有效回波信号的幅值较大,干扰成分较低。
步骤4,将有效回波信号经过模拟电压比较器整形为数字回波脉冲序列,传送给处理器,每个数字脉冲的宽度对应回波信号的周期,每个数字脉冲的位置对应具体回波信号的位置;计算出各个回波脉冲的宽度;
步骤4具体为:
按照回波脉冲的次序记下每个回波脉冲的上升边沿时刻,并记录到上升边沿时刻数组t[i]中,i为回波脉冲的序数;
根据t[i]确定各个回波脉冲的宽度,记录到回波脉冲宽度数组w[i]中,w[i]=t[i]-t[i-1],i为回波脉冲的序数。
步骤5,找到回波脉冲的宽度变化转折点m,即为回波脉冲频率转折点,记录该点的时刻为t[m];
步骤5具体为:
根据回波脉冲宽度数组w[i],找到回波脉冲宽度变化转折点m:
若m满足:
Figure BDA0001389049810000051
w[m]>w[m+1],
w[m+1]>w[m+2],
Figure BDA0001389049810000052
根据回波脉冲宽度数组w[i],找到回波脉冲宽度变化转折点m:即在接收到的回波脉冲序列中,存在一段包含n个脉冲的连续脉冲串,回波脉冲满足先连续递增
Figure BDA0001389049810000053
个脉冲达到m点,再从m点连续递减
Figure BDA0001389049810000054
个脉冲的规则,将其中心位置m点确定为回波脉冲频率转折点,记录该点的时刻为t[m];采用这种方式确定回波脉冲频率转折点而不能根据发射的脉冲串的脉冲宽度严格比较以确定回波脉冲频率转折点,这是根据超声波传感器的特性决定的,超声波发射和接收短暂的脉冲串,频率不能稳定,但是频率变化趋势能够反映出来。
步骤6,计算回波时间techo,techo=t[m]-t0
实施例:采用工作频率为125Khz的收发一体化超声波传感器,以STM32F405RGT6处理器为核心实现本发明。
超声波发射脉冲串采用定时器TIM2控制,脉冲发送规则为:以超声波传感器中心频率125Khz为起点,频率分5个脉冲逐步降到124Khz,再从124Khz分4个脉冲逐步升到125Khz,发送共9个脉冲,经过放大,发射脉冲幅值大约为200Vpp,当发射脉冲的频率为124Khz时,使定时器TIM3从零开始计数。
接收的回波信号经过总增益大约为50db三级放大滤波,经过电压比较器转换成数字脉冲信号,输入至处理器的定时器TIM3输入端,利用定时器TIM3的捕获功能记录每个数字回波脉冲上升边沿的时刻。
发射超声波脉冲串后,经过大约1ms延时,去除超声波传感器的余振的影响,打开定时器TIM3的捕获中断,在每个数字回波脉冲上升沿,产生定时器TIM3的捕获中断,将TIM3当前的计数值依次保存到数组t[40]中,经过大约20ms延时,接收过程完成。
通过程序计算收到的各个回波脉冲宽度,保存到数组w[40]中:w[i]=t[i]-t[i-1];
根据收到的回波脉冲宽度数组w[i],在其中查找回波脉冲频率转折点,即查找回波脉冲频率连续逐渐增加5次,又连续降低4次的点,若找到的点为m,可以确定回波时间techo为t[m]。
实验表明,该测试装置在恒定室温环境下,在目标距离为1米范围以内,回波测量误差小于±5us,在目标距离为2米范围,回波测量误差小于±10us,若再大于该范围误差明显增加;考虑声波传播速度为340m/s,可以折算出距离测量误差在目标距离为1米范围以内,大约为±0.85mm,在目标距离为2米范围,回波测量误差小于±1.7mm,测量精度明显高于传统方法。
通过上述方式,线性变频的回波时间测量方法,通过在超声波传感器中心频率附近,发送逐脉冲频率先线性递减再线性递增的脉冲串,有助于传感器跟随激励信号变化,回波信号输出比较稳定;再将回波信号幅度放大并整形成数字脉冲序列,测量每个脉冲的宽度(频率),以确定回波信号的频率转折点,发射信号的频率转折点和回波信号的频率转折点之差为有效的回波时间;本发明不需要采样模拟回波信号,降低了系统对模数转换器、存储器等超声波检测距装置的硬件要求,成本较低,不需要复杂的算法,可靠性强,容易实现,具有较高的实用性。解决了利用回波前沿检测回波信号精度低的问题;解决了基于互相关函数法、频域相位法和时域方法需要高速的模数转换器和大量数据存储,分析方法计算量大,实时性差,实用性较低的问题,解决了FSK(频移键控)检测回波时间频率切换点难以确定的问题。

Claims (3)

1.线性变频的回波时间测量方法,其特征在于,包括以下步骤:
步骤1,超声波传感器发射周期先线性增加后线性减小的脉冲串,记录发射脉冲串的频率转折时刻t0
步骤2,经过延时以消除所述超声波传感器的余振影响;
步骤3,所述超声波传感器在发射脉冲串的同时接收超声波回波信号,所述超声波回波信号经过放大、滤波电路处理后得到有效回波信号;
步骤4,将所述有效回波信号经过模拟电压比较器整形为数字回波脉冲序列,并传送给处理器,计算出各个回波脉冲的宽度;所述步骤4具体为:
按照回波脉冲的次序记下每个回波脉冲的上升边沿时刻,并记录到上升边沿时刻数组t[i]中,i为回波脉冲的序数;
根据t[i]确定各个回波脉冲的宽度,记录到回波脉冲宽度数组w[i]中,w[i]=t[i]-t[i-1],i为回波脉冲的序数;
步骤5,找到回波脉冲的宽度变化转折点m,即为回波脉冲频率转折点,记录该点的时刻为t[m];所述步骤5具体为:
根据回波脉冲宽度数组w[i],找到回波脉冲宽度变化转折点m:即在接收到的回波脉冲序列中,存在一段包含n个脉冲的连续脉冲串,回波脉冲满足先连续递增
Figure FDA0002405984720000011
个脉冲达到m点,再从m点连续递减
Figure FDA0002405984720000012
个脉冲的规则,将其中心位置m点确定为回波脉冲频率转折点,记录该点的时刻为t[m];
步骤6,计算回波时间techo,techo=t[m]-t0
2.根据权利要求1所述的线性变频的回波时间测量方法,其特征在于,所述步骤1具体为:
超声波传感器发送周期以
Figure FDA0002405984720000021
为起始,周期渐变的步长为Δt的,周期先线性增加后线性减小的脉冲串,当发送宽度为
Figure FDA0002405984720000022
的脉冲时,记录当时的时刻t0,即为发射脉冲串的频率转折时刻;
其中,fs为超声波传感器的中心频率,n为脉冲串中脉冲的个数。
3.根据权利要求2所述的线性变频的回波时间测量方法,其特征在于,所述步骤1中,n取值为奇数,n的取值范围为7-15。
CN201710740712.0A 2017-08-25 2017-08-25 线性变频信号的回波时间测量方法 Active CN107576964B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710740712.0A CN107576964B (zh) 2017-08-25 2017-08-25 线性变频信号的回波时间测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710740712.0A CN107576964B (zh) 2017-08-25 2017-08-25 线性变频信号的回波时间测量方法

Publications (2)

Publication Number Publication Date
CN107576964A CN107576964A (zh) 2018-01-12
CN107576964B true CN107576964B (zh) 2020-05-22

Family

ID=61035258

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710740712.0A Active CN107576964B (zh) 2017-08-25 2017-08-25 线性变频信号的回波时间测量方法

Country Status (1)

Country Link
CN (1) CN107576964B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112344888B (zh) * 2020-10-23 2021-10-22 华南理工大学 一种基于lfm信号互相关相位特性的超声波测厚方法
CN114088972B (zh) * 2021-11-21 2024-04-05 吉林大学 一种基于相位双频法的超声波风速风向测量系统及方法
CN113970744B (zh) * 2021-12-26 2022-04-05 湖南天羿领航科技有限公司 基于线性调频超声波的测距方法、系统及介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101029932A (zh) * 2007-01-24 2007-09-05 哈尔滨工业大学 基于调制域测量的超声波回波前沿检测方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000046931A (ja) * 1998-07-24 2000-02-18 Mitsubishi Electric Corp レーダ装置
CN100558299C (zh) * 2007-04-30 2009-11-11 陆尧胜 一种产程监护装置
CN101089653B (zh) * 2007-07-20 2011-03-09 西安理工大学 近程调频连续波fmcw雷达抗干扰方法
US7924216B2 (en) * 2008-04-30 2011-04-12 Rosemount Tank Radar Ab Method of determining a disturbance echo profile for a radar level gauge system
CN101852854B (zh) * 2010-06-07 2012-10-31 华南理工大学 一种水下多波束测探系统及其探测方法
CN102004256B (zh) * 2010-09-09 2012-07-04 北京航空航天大学 基于空间谱全息存储的激光干涉测距系统
CN102707288B (zh) * 2012-06-13 2013-11-06 西安理工大学 超声波回波时间的检测方法
CN103345145B (zh) * 2013-05-27 2016-03-30 北京空间飞行器总体设计部 一种利用激光进行星载时钟测量的方法
CN103454643B (zh) * 2013-09-10 2015-05-27 浙江工商大学 恒定声压fsk超声波渡越时间精确测量方法
CN204178554U (zh) * 2014-10-27 2015-02-25 北京握奇智能科技有限公司 一种自由流收费的激光扫描车型识别系统
CN106383346B (zh) * 2016-09-14 2019-04-02 深圳天眼激光科技有限公司 脉冲激光扫描回波接收电路、接收方法及脉冲式激光扫描仪

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101029932A (zh) * 2007-01-24 2007-09-05 哈尔滨工业大学 基于调制域测量的超声波回波前沿检测方法

Also Published As

Publication number Publication date
CN107576964A (zh) 2018-01-12

Similar Documents

Publication Publication Date Title
CN103542901B (zh) 一种流量计
US20170052253A1 (en) Method and System for Ultrasound Time-of-Flight Measurement
US8648612B2 (en) Calibration of a distance measuring device
CN107576964B (zh) 线性变频信号的回波时间测量方法
CN107860430B (zh) 一种基于时差法的超声波气体流量计时间差测量方法
WO1997014936A1 (en) Digital speed determination in ultrasonic flow measurements
US11885874B2 (en) Acoustic distance measuring circuit and method for low frequency modulated (LFM) chirp signals
CN109579950B (zh) 一种气体超声波流量计的防错波检测装置
JP2007327956A (ja) 距離計測方法及び距離計測装置
CN201177666Y (zh) 一种根据相位差比较的超声波测距装置
EP3903124A1 (en) Ultrasonic echo processing in presence of doppler shift
CN111693136A (zh) 一种采用回波信号自相关相位谱的声表面波谐振器频率估计算法
CN111238587A (zh) 一种超声波换能器传输时间校准方法
EP2457111A1 (en) Measurement method and apparatus
CN102967334B (zh) 利用对信号包络线处理测量流体流量的系统及方法
US8639462B2 (en) Method and system for determining the time-of-flight of a signal
KR100739506B1 (ko) 정합필터의 간략한 계산을 사용한 초음파 거리 정밀측정방법
CN110988853A (zh) 基于有效峰值的超声波到达时间计算方法
CN113970744B (zh) 基于线性调频超声波的测距方法、系统及介质
CN116046116A (zh) 一种基于Lobe算法的超声波水表流量校准方法
JP2003279396A (ja) 超音波流量計
KR20170035555A (ko) 초음파 유량계의 신호처리시스템
US11725967B2 (en) Gas volume determination in fluid
RU2612201C1 (ru) Способ определения дистанции гидролокатором
CN106248991A (zh) 一种低信噪比下的adcp解测速模糊方法

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