CN102590598A - 基于多阈值比较的超声波信号过零点预测方法 - Google Patents

基于多阈值比较的超声波信号过零点预测方法 Download PDF

Info

Publication number
CN102590598A
CN102590598A CN2012100340292A CN201210034029A CN102590598A CN 102590598 A CN102590598 A CN 102590598A CN 2012100340292 A CN2012100340292 A CN 2012100340292A CN 201210034029 A CN201210034029 A CN 201210034029A CN 102590598 A CN102590598 A CN 102590598A
Authority
CN
China
Prior art keywords
voltage
zero crossing
signal
ultrasound wave
sample point
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
CN2012100340292A
Other languages
English (en)
Other versions
CN102590598B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201210034029.2A priority Critical patent/CN102590598B/zh
Publication of CN102590598A publication Critical patent/CN102590598A/zh
Application granted granted Critical
Publication of CN102590598B publication Critical patent/CN102590598B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Volume Flow (AREA)

Abstract

本发明公开了一种基于多阈值比较的超声波信号过零点预测方法。现有过零检测方法将信号与零电位进行电压比较,直接以测量方式确定信号过零点,存在抗干扰能力差的问题。本发明提出了一种基于多阈值比较的超声波信号过零点预测方法。该发明采用测量与计算相结合的方法,通过多阈值比较以测量方式获取样本点参数,以此进行超声波波形曲线拟合,进而采用拟合所得的曲线方程计算出超声波接收信号的过零点。本发明有效解决了传统过零检测方法抗干扰能力差的问题,提高了超声波信号过零检测的精度与可靠性。此发明适用于各种超声波检测领域,尤其适用于接受波信号信噪比低的气体超声波流量计,能够大大提高其检测精度和抗干扰能力。

Description

基于多阈值比较的超声波信号过零点预测方法
技术领域
本发明涉及一种超声波信号过零点检测方法,特别是涉及一种基于多阈值比较的超声波信号过零点预测方法。
背景技术
在超声波信号检测技术中,阈值检测技术和过零检测技术是其中两大检测技术。所谓阈值检测(阈值比较),是指用一个高出零电位许多的电压值作为比较器的触发信号,当信号的电压值高于此触发信号时,比较器输出置高,否则置低;所谓过零检测(过零比较),是指用非常接近于零电位的电压值作为比较器的触发信号,当信号的电压值高于此触发信号时,比较器输出置高,否则置低。
阈值检测相对于过零检测而言,由于其对比较电压的要求不高,因此对比较器性能的要求也不高,成本较低。但由于超声波接收信号通常存在幅值波动和波型扭曲的现象,采用阈值检测容易产生较大的误差。而由于超声波在传播过程中频率保持不变的特性,使得其在固定的声程内,在相同环境的同一介质中其过零点始终不变,因此采用过零检测技术较阈值比较技术能获得更高的精度。
不过,过零比较技术也存在其固有的弊病——处理信噪比较低信号时精度较差。在固定的声程内,在相同环境的同一介质中,超声波信号的过零点始终不变的,但由于环境影响,超声波信号线上总会耦合上不规则的噪声干扰,使得直接测量信号的过零点总会受到噪音的干扰,降低测量精度。
为了解决过零检测信噪比低的问题,一般有两种方法:一种是利用多次测量求平均的方法,简单有效。但测量次数的增加将带来整机功耗的增大,随着对产品低功耗的要求越来越高,此方法已不适用于低功耗超声波产品之中,尤其是不适用于以电池为电源的超声波仪表之中;另一种方法是购买一片高性能的过零比较器和设计一套高性能的滤波电路。此方法已牺牲成本为代价,不适用于低成本的超声波民用产品。
发明内容
本发明的目的在于提供一种基于多阈值比较的超声波信号过零点预测方法,意在用低成本有效解决超声波信号过零检测抗干扰能力差的问题。
为实现上述目的,本发明技术实施方案如下:
1、利用多阈值比较获取曲线拟合所需的样本点横纵坐标参数:
在原始超声波接收信号进入由多个电压比较器集成的多电压比较器之前,先进行电压偏置处理。通过采用加法电路把原始超声波接收信号的电压幅值抬高一个偏置电压,使其转变成电压幅值均高于零电位电压的电压偏置超声波接收信号。因此,多阈值比较所需的多个比较电压只需通过分压电阻以电阻分压方式即能产生。根据电压比较器的工作原理,电阻分压产生的比较电压Vi即为样本点的纵坐标值。当参考电压为Vref,N个分压电阻进行电阻分压时,第i个电阻值为Ri的分压电阻处的比较电压Vi通过以下公式计算:
V i = R i + R i + 1 + . . . + R n R 1 + R 2 + . . . + R n V ref
由于电压偏置超声波接收信号是一个时域信号,因此所述样本点的横坐标为一个时间参数ti。通过把多电压比较器产生的触发脉冲信号分别引入具有引脚中断功能的单片机或者专用计时芯片的计时单元中,利用触发脉冲信号上升沿触发计时单元工作即能完成对样本点横坐标参数ti的测量;
2、利用曲线拟合,计算出超声波接收信号的过零点:
采用最小二乘法曲线拟合时,至少需有四个样本点和超声波接收波波形函数为型函数才能进行曲线拟合;采用多项式等插值方法曲线拟合时,至少有二个样本点即能进行曲线拟合;以上样本点的个数由所需测量精度决定。
利用曲线拟合所得的数学表达式即能计算出超声波接收信号的过零点。
本发明具有的有益效果是:
(1)降低了对比较器和滤波电路的要求,节约了成本。
(2)采用曲线拟合,利用拟合所得的曲线方程计算超声波信号过零点,起到了软件滤波的作用,提高了超声波过零检测的精度和抗干扰能力。
本发明有效解决了传统过零检测方法抗干扰能力差的问题,提高了超声波信号过零检测的精度与可靠性。本发明适用于各种超声波检测领域,尤其适用于接受波信号信噪比低的气体超声波流量计,能够大大提高其检测精度和抗干扰能力。
附图说明
图1是本发明中采用的多阈值比较采样示意图。
图2是本发明中采用的样本点坐标示意图。
图3是本发明中采用的曲线拟合示意图。
图1中:1、原始超声波接收信号,2、电压偏置超声波接收信号,3、偏置电压,4、参考电压,5、分压电阻,6、比较电压,7、多电压比较器,8、触发脉冲,9、计时单元。
图2中:10、比较电压V1,11、比较电压V2,12、比较电压Vi,13、比较电压Vn-1,14、比较电压Vn
图3中:15、样本点,16、曲线,17、零点。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
本发明采用测量与计算相结合的方法,通过多阈值比较获取样本点,以此进行超声波波形曲线拟合,进而采用曲线拟合所得的数学表达式计算出超声波信号的过零点。
该方法主要分为以下二个步骤:
步骤1)利用多阈值比较获取曲线拟合所需的样本点横纵坐标参数。如图1所示,为了降低对多电压比较器7的性能要求,避免负电源供电,节约成本,在原始超声波接收信号1进入多电压比较器7之前,必须先进行电压偏置处理,把原始超声波接收信号1的电压幅值抬高一个偏置电压3,使其转变成电压幅值均高于零电位电压的电压偏置超声波接收信号2。因此,多阈值比较所需的多个比较电压也只需通过分压电阻5以电阻分压方式产生。根据电压比较器的工作原理,电阻分压产生的比较电压Vi即为样本点的纵坐标值。当参考电压为Vref,N个分压电阻进行电阻分压时,第i个电阻值为Ri的分压电阻处的比较电压Vi通过以下公式计算:
V i = R i + R i + 1 + . . . + R n R 1 + R 2 + . . . + R n V ref - - - ( 1 )
由于电压偏置超声波接收信号是一个时域信号,因此样本点的横坐标为一个时间参数ti。通过把多电压比较器产生的触发脉冲信号分别引入具有引脚中断功能的单片机或者专用计时芯片的计时单元9中,利用触发脉冲信号上升沿触发计时单元工作即能完成对样本点横坐标参数ti的测量。通过以上方法所得的样本点如图2所示。
步骤2)利用超声波接收波波形函数进行曲线拟合,计算出超声波接收信号的过零点。常用的曲线拟合有多项式曲线拟合、最小二乘法曲线拟合等。下面以最小二乘法曲线拟合为例,阐述曲线拟合并计算超声波接收信号过零点的过程。
通过对超声波系统的建模,获得以下超声波接收波形的时域波形函数:
y = b ( e - at · sin gt g 2 + a 2 ( - a sin ωt + g cos ωt )
- e - at · cos gt - 1 g 2 + a 2 ( g sin ωt + a cos ωt ) - - - ( 2 )
- e - at · sinh t h 2 + a 2 ( a sin ωt + h cos ωt )
- e - at · cosh t - 1 h 2 + a 2 ( h sin ωt - a cos ωt )
注:a,指数系数;b,幅值系数;g,和角频率系数;h,差角频率系数;ω,角频率,ω=2πf,f为超声波频率。
在利用获取的样本点和上述型函数进行最小二乘法进行曲线拟合时,还需分二个步骤:
步骤1:确定公式1中的a、b、g、h四个待定系数的最优初值。
首先利用换能器阻抗分析仪分别测出以下参数:
1、介质动态电阻R0
2、发射换能器动态电阻R1、动态电感L1、动态电容C1和静态电容C01
3、接收换能器动态电阻R2、动态电感L2、动态电容C2和静态电容C02
幅值系数b的最优初值取1,a、g、h三个系数的最优初值分别利用公式3、公式4和公式5计算所得,公式如下:
a = R 0 + R 1 + R 2 2 ( L 1 + L 2 ) - - - ( 3 )
g = 1 ( L 1 + L 2 ) C 1 · C 2 C 1 + C 2 + 1 ( L 1 + L 2 ) C 02 - a 2 + ω 10 - - - ( 4 )
h = 1 ( L 1 + L 2 ) C 1 · C 2 C 1 + C 2 + 1 ( L 1 + L 2 ) C 02 - a 2 - ω 10 - - - ( 5 )
步骤2:以公式2为型函数进行最小二乘法曲线拟合,获得超声波接收波形的曲线函数,计算出超声波接收信号的过零点。把已获得的样本点参数和a、b、g、h四个系数的最优值代入最小二乘法中,即可获得a、b、g、h四个系数的最优值,将其代入公式2即可获得对应以上样本点的超声波接收波形的曲线函数,通过计算即可获得超声波接收信号的过零点。如图3所示,利用样本点15拟合出曲线16,零点17即为通过曲线16计算所得的超声波信号过零点。由于采用了以所拟合曲线里各采样点距离的均方和最小为优化目标的最小二乘法曲线拟合方式,使得阈值比较采样点上噪声的影响被有效减弱,从而确保了过零点预测的准确性和可靠性。

Claims (4)

1.一种基于多阈值比较的超声波信号过零点预测方法,其特征在于该方法的步骤如下:
步骤1)利用多阈值比较获取曲线拟合所需的样本点横纵坐标参数:
在原始超声波接收信号进入多电压比较器之前,先进行电压偏置处理,把原始超声波接收信号的电压幅值抬高一个偏置电压,使其转变成电压幅值均高于零电位电压的电压偏置超声波接收信号;因此,多阈值比较所需的多个比较电压只需通过分压电阻以电阻分压方式就能产生;根据电压比较器的工作原理,电阻分压产生的比较电压Vi即为样本点的纵坐标值;当参考电压为Vref,N个分压电阻进行电阻分压时,第i个电阻值为Ri的分压电阻处的比较电压Vi通过以下公式计算:
V i = R i + R i + 1 + . . . + R n R 1 + R 2 + . . . + R n V ref
由于电压偏置超声波接收信号是一个时域信号,因此所述样本点的横坐标为一个时间参数ti;通过把多电压比较器产生的触发脉冲信号分别引入计时单元中,利用触发脉冲信号上升沿触发计时单元工作即能完成对样本点横坐标参数ti的测量;
步骤2)利用曲线拟合,计算出超声波接收信号的过零点:
采用最小二乘法曲线拟合时,至少需有四个样本点和超声波接收波波形函数为型函数才能进行曲线拟合;采用多项式等插值方法曲线拟合时,至少有二个样本点即能进行曲线拟合;以上样本点的个数由所需测量精度决定;
利用曲线拟合所得的数学表达式即能计算出超声波接收信号的过零点。
2.根据权利要求1所述的一种基于多阈值比较的超声波信号过零点预测方法,其特征在于:所述的多电压比较器由多个电压比较器集成,具有多个电压比较功能,用于对不同的电压值进行比较,以获取样本点的纵坐标值;多电压比较器中集成的电压比较器个数取决于样本点的个数。
3.根据权利要求1所述的一种基于多阈值比较的超声波信号过零点预测方法,其特征在于:所述的电压偏置处理,是指通过加法电路将原始超声波接收信号的电压幅值均抬高一个偏置电压的过程。
4.根据权利要求1所述的一种基于多阈值比较的超声波信号过零点预测方法,其特征在于:所述的计时单元,是指具有引脚中断功能的单片机计时模块或者是专用的计时芯片。
CN201210034029.2A 2012-02-15 2012-02-15 基于多阈值比较的超声波信号过零点预测方法 Active CN102590598B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210034029.2A CN102590598B (zh) 2012-02-15 2012-02-15 基于多阈值比较的超声波信号过零点预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210034029.2A CN102590598B (zh) 2012-02-15 2012-02-15 基于多阈值比较的超声波信号过零点预测方法

Publications (2)

Publication Number Publication Date
CN102590598A true CN102590598A (zh) 2012-07-18
CN102590598B CN102590598B (zh) 2014-01-29

Family

ID=46479522

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210034029.2A Active CN102590598B (zh) 2012-02-15 2012-02-15 基于多阈值比较的超声波信号过零点预测方法

Country Status (1)

Country Link
CN (1) CN102590598B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932037A (zh) * 2015-12-30 2017-07-07 浙江大学 一种时差式超声波流量计接收信号波峰序列位置判定方法
CN107367633A (zh) * 2016-05-11 2017-11-21 佛山市顺德区美的电热电器制造有限公司 用于电压过零检测电路的过零处理方法、装置及设备
CN109597057A (zh) * 2018-12-06 2019-04-09 苏州镭图光电科技有限公司 一种激光回波测距方法及测距装置
CN110376424A (zh) * 2018-04-12 2019-10-25 佛山市顺德区美的电热电器制造有限公司 电磁加热系统及其电压过零信号处理方法、装置
CN110531348A (zh) * 2019-08-02 2019-12-03 陈小年 雷达测距方法、装置、计算机设备和存储介质
CN113740856A (zh) * 2021-07-21 2021-12-03 国创移动能源创新中心(江苏)有限公司 基于超声波波形起振特性的测距方法和装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1645885A1 (de) * 2004-10-08 2006-04-12 Abb Research Ltd. Verfahren zur Bestimmung des Stromnulldurchgangs in einem Umrichter
CN101542662A (zh) * 2006-12-01 2009-09-23 Abb技术有限公司 用于预测在电流通路中电流的未来行为的方法和设备
CN101833031A (zh) * 2010-05-24 2010-09-15 中国北车股份有限公司大连电力牵引研发中心 交流信号过零检测电路
CN102243263A (zh) * 2011-04-08 2011-11-16 苏州路之遥科技股份有限公司 一种交流电流值检测及电流过零点检测电路

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1645885A1 (de) * 2004-10-08 2006-04-12 Abb Research Ltd. Verfahren zur Bestimmung des Stromnulldurchgangs in einem Umrichter
CN101542662A (zh) * 2006-12-01 2009-09-23 Abb技术有限公司 用于预测在电流通路中电流的未来行为的方法和设备
CN101833031A (zh) * 2010-05-24 2010-09-15 中国北车股份有限公司大连电力牵引研发中心 交流信号过零检测电路
CN102243263A (zh) * 2011-04-08 2011-11-16 苏州路之遥科技股份有限公司 一种交流电流值检测及电流过零点检测电路

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932037A (zh) * 2015-12-30 2017-07-07 浙江大学 一种时差式超声波流量计接收信号波峰序列位置判定方法
CN106932037B (zh) * 2015-12-30 2019-04-02 浙江大学 一种时差式超声波流量计接收信号波峰序列位置判定方法
CN107367633A (zh) * 2016-05-11 2017-11-21 佛山市顺德区美的电热电器制造有限公司 用于电压过零检测电路的过零处理方法、装置及设备
CN107367633B (zh) * 2016-05-11 2019-08-27 佛山市顺德区美的电热电器制造有限公司 用于电压过零检测电路的过零处理方法、装置及设备
CN110376424A (zh) * 2018-04-12 2019-10-25 佛山市顺德区美的电热电器制造有限公司 电磁加热系统及其电压过零信号处理方法、装置
CN110376424B (zh) * 2018-04-12 2021-05-25 佛山市顺德区美的电热电器制造有限公司 电磁加热系统及其电压过零信号处理方法、装置
CN109597057A (zh) * 2018-12-06 2019-04-09 苏州镭图光电科技有限公司 一种激光回波测距方法及测距装置
CN110531348A (zh) * 2019-08-02 2019-12-03 陈小年 雷达测距方法、装置、计算机设备和存储介质
CN110531348B (zh) * 2019-08-02 2021-09-10 陈小年 雷达测距方法、装置、计算机设备和存储介质
CN113740856A (zh) * 2021-07-21 2021-12-03 国创移动能源创新中心(江苏)有限公司 基于超声波波形起振特性的测距方法和装置
CN113740856B (zh) * 2021-07-21 2023-12-22 国创移动能源创新中心(江苏)有限公司 基于超声波波形起振特性的测距方法和装置

Also Published As

Publication number Publication date
CN102590598B (zh) 2014-01-29

Similar Documents

Publication Publication Date Title
CN102590598B (zh) 基于多阈值比较的超声波信号过零点预测方法
CN101701984B (zh) 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法
CN101701982B (zh) 基于加窗插值fft的电力系统谐波检测方法
CN101813725B (zh) 一种低频率信号的相位差测量方法
CN104808056B (zh) 一种基于比较器转换的频率特性测试方法与装置
CN106768109A (zh) 基于回波上升段峰值拟合和基于回波能量点定位的气体超声波流量计信号处理方法
CN106483374A (zh) 一种基于Nuttall双窗全相位FFT的谐波间谐波检测方法
CN102565574A (zh) 船舶电能质量监测系统
CN103197141A (zh) 一种测量电力系统信号频率及谐波参数的方法
CN109579950B (zh) 一种气体超声波流量计的防错波检测装置
CN103235180A (zh) 一种电网间谐波测量方法
CN103353550A (zh) 一种测量电力系统信号频率及谐波参数的方法
CN103063909B (zh) 一种基于功率谱的线性调频信号参数估值方法
CN102809687B (zh) 一种交流电频率的数字化测量方法
CN108318852A (zh) 一种智能电能表方波影响试验方法
CN105890685A (zh) 一种基于累积相位差的超声波流量测量装置
CN103399204A (zh) 一种基于Rife-Vincent(II)窗插值FFT的谐波与间谐波检测方法
CN102095929B (zh) 一种快速测量交流电信号频率的方法
CN103472333A (zh) 风电并网电能质量综合性能检测方法
CN105425128A (zh) 一种变压器局部放电超声波检测及精确定位装置与方法
CN101718816B (zh) 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法
CN103983849A (zh) 一种实时高精度的电力谐波分析方法
CN103822702A (zh) 一种用于超声相控阵换能器声场测量数据采集同步触发的脉冲信号转换装置
CN103795411A (zh) 基于五项最大旁瓣衰减窗三谱线插值测试sfdr的方法
CN102841343A (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
C14 Grant of patent or utility model
GR01 Patent grant