CN106371098A - 一种有阻尼自由振动超声波测距系统频率拐点判别方法 - Google Patents

一种有阻尼自由振动超声波测距系统频率拐点判别方法 Download PDF

Info

Publication number
CN106371098A
CN106371098A CN201610778428.8A CN201610778428A CN106371098A CN 106371098 A CN106371098 A CN 106371098A CN 201610778428 A CN201610778428 A CN 201610778428A CN 106371098 A CN106371098 A CN 106371098A
Authority
CN
China
Prior art keywords
point
frequency
signal
ultrasonic
data
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
CN201610778428.8A
Other languages
English (en)
Other versions
CN106371098B (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 Gongshang University
Original Assignee
Zhejiang Gongshang 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 Zhejiang Gongshang University filed Critical Zhejiang Gongshang University
Priority to CN201610778428.8A priority Critical patent/CN106371098B/zh
Publication of CN106371098A publication Critical patent/CN106371098A/zh
Application granted granted Critical
Publication of CN106371098B publication Critical patent/CN106371098B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/06Systems determining the position data of a target
    • G01S15/08Systems for measuring distance only

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种有阻尼自由振动超声波测距系统频率拐点判别方法,该方法具体为:控制器通过正弦波发生器输出正弦波,该正弦波通过功率放大电路放大后激励超声波发射器发出超声波,同时超声波接收器开始接收超声波。超声波接收器接收到微弱电信号通过整形、放大后进行A/D采样,采样后的数据由DSP数据处理器进行逐个点的FFT变换,得到每个点的幅度谱和相位谱。根据幅度谱中某一点和其相邻点的幅度值所对应的频率的差别,可以明确地得出超声波信号相邻两个点的频率变化,从而准确地判别出有阻尼自由振动的超声波测距系统的频率拐点。

Description

一种有阻尼自由振动超声波测距系统频率拐点判别方法
技术领域
本发明涉及一种测量方法,尤其涉及一种有阻尼自由振动超声波测距系统中频率拐点判别方法。
背景技术
声波是一种机械波,按振动频率可分为:20kHz频率以上的高频区、20Hz-20kHz频率的中频区以及低于20Hz的低频区。频率在20kHz以上的声波,其振动频率已经超出人耳所能感知的最高频率上限,因而称为超声波。超声波具备反射、折射、衍射以及散射等机械波的通性。超声波因其振动频率较高故而具备一些特有性质:其一是功率大,其能量比一般声波大得多,有较强的穿透能力,可以用做焊接、切削、钻孔等;其二是因其频率高,波长较短,一般的障碍物尺寸比超声波波长较长一些,故其衍射能力很差,反而具有良好的定向性,所以超声波广泛应用于工业与医学上的非接触检测上。
随着超声波测距的应用逐渐普及,对测量的精确度要求也逐渐更加严格。超声波在同一介质中的传播速度基本恒定,因此利用测量超声波在该介质中经过一段距离所用的时间,可以实现非接触式测距。超声波在不同的介质中传播的速度不同,利用这一特性,可以测定介质的特性,如测定酒精的浓度,测量锅炉中空气的温度等等。如果测定了不同介质的超声波传播特性,则可以测定不同物质的类型,如蔗糖溶液、酒精溶液、各种酒类的品质特性等等。随着社会的进步,很多场合对于超声检测的精度要求也越来越高。为了提高测距的精确度,必须精确测量超声波在介质中传播时所经过的时间(即踱越时间)。
超声波在非均匀媒介中传播会出现声波随着距离而逐渐衰减的现象,在传播的过程中也会受环境温度等物理条件的影响。在传统的超声波测距中,采用40kHz的超声波,在介质传播过程中衰减较快,长距离测量时测量误差较大,所以40kHz的超声波测距量程较小。若选择超声波频率过低,则导致其方向性不好,当选用25kHz频率的超声波时,这种频率超声波既能保证良好方向性,又能降低声衰减的影响。选取合适频率的超声波,同时利用温度标杆法,准确得出相同温度下超声波传播的速度,可以将超声波传播过程中逐渐衰减以及受环境温度影响造成的误差降到最低,甚至可以忽略不计。
超声波接收器是一种机械振动装置,接收压电振子在电信号驱动下,同时受到弹性力的阻力的作用,接收振子作有阻尼的受迫振动。此时会在贴近振子表面的层中产生压力,此压力从一个质点传到另一个质点,即声波从振子表面传播开去并形成声场。假定接收压电振子在零状态(振动位移s(t)|t=0=0,振动速度s′(t)|t=0=0)下受电信号激励,电信号产生的应变力为式中Ca为接收压电振子的电容量,d33为压电系数。其动力学方程为:
m d 2 s dt 2 + r s d s d t + K s s = F ~ s - - - ( 1 )
式中:接收压电振子的质量为m,阻尼为rs,劲度为Ks(决定于发射压电振子系统的弹性系数和几何尺寸),s是接收压电振子辐射表面的位移。
当输入电信号为带入式(1)可得:
m d 2 s dt 2 + r s d s d t + K s s = C a U m d 33 c o s ( ω t ) - - - ( 2 )
由于同时受到弹性力的阻力的作用,接收振子作有阻尼的受迫振动,其谐振角频率为:接收压电振子系统的阻尼系数其中Rs为阻尼系数,ms为系统质量,其振动状态为:
由于接收压电振子在零状态开始振动,即振动位移s(t)|t=0=0,振动速度s′(t)|t=0=0,即
A 0 = - 1 ( ω s 0 2 - ω 2 ) 2 + 4 β s 2 ω 2 ω s 0 ω s 0 2 - β s 2 - - - ( 4 )
式(3)化简可得:
以上是正弦电信号一直持续作用于换能器的振动状态,接收压电振子会经历两个振动状态:①当经过tc(起振时间)时间后,接收端振子只有超声波的受迫振动起作用,即振动角频率为ω的振幅为kA的等幅振动,此时达到位移幅值最大且稳定持续的振动状态,如果电信号不停止,接收压电振子一直处于此振动状态。②位移幅度下降,此时的振动状态是以有阻尼振动角频率为谐振频率作指数规律的衰减振动。
当输入电信号的角频率ω=ωs0时,接收压电振子的振动方程为:
s ( t ) = k 2 β s ω s 0 ( 1 - e - β s t ) s i n ( ω s 0 t ) - - - ( 7 )
对上式分析可知,此时从停止振动到达到稳定振动状态,其振动频率为谐振角频率ωs0并保持不变,变化的只是其振动的幅值。当电信号停止的时候,接收压电振子会以电信号停止时刻的振动状态为初始条件作有阻尼的自由振荡。若激励信号在N个周期后停止,在接收压电振子振动状态应为:
接收压电振子在NT时刻的振动状态:振动位移s(t)|t=NT和振动速度s′(t)|t=NT分别为:
将以上带入通解项可求得
由式(11)可知,在接收激励信号停止之后,接收压电振子会在原有振动状态下不断衰减直至能量为0,振子才停止振动,也是所谓的“拖尾”现象。
综上所述,由于超声波的传播速度受环境温度影响、超声波振子受机械阻尼作用而产生起振点延迟、“拖尾”等因素都会对测距系统的精度产生影响。
专利“基于有阻尼自由振动的超声波踱越时间测量方法”巧妙地利用超声波接收振子的有阻尼自由振动特性,通过计算超声波接收器振子受迫振动周期和振子振动于有阻尼自由振动周期的时间,可以得到超声波接收振子受超声波激励频率ω和振子振动于有阻尼自由振动频率ω'r。当超声波激励消失时,超声波接收振子会脱离超声波激励频率ω而振动于有阻尼自由振动频率ω'r。而且超声波激励频率ω和有阻尼自由振动频率ω'r都是单频点,控制器通过判别ω到ω'r的变化,来识别超声波接收振子所在位置超声波激励信号是否消失,从而准确判断出超声波传播的渡越时间,有效地提高测距精度。
但是,接收振子有阻尼自由振动的时间很短,甚至是短短几个周期,因此有阻尼自由振动频率ω'r并不容易可靠检测。
发明内容
本发明的目的是针对现有技术的不足,提供一种有阻尼自由振动超声波测距系统频率拐点判别方法。
本发明的目的是通过以下技术方案来实现的:一种有阻尼自由振动超声波测距系统频率拐点判别方法,该方法包括以下步骤:
(1)正弦波发生器输出n个波形完整且连续的正弦波信号,该正弦波信号通过功率放大电路放大后激励超声波发射器发出超声波,同时超声波接收器开始接收超声波,超声波接收器输出的微弱电信号经信号调理电路放大后,输出至A/D转换器。超声波接收器输出的微弱电信号包含了超声波接收器振子受迫振动状态和有阻尼自有振动状态两种状态的信号,超声波接收器接收的超声波频率记为ω,接收器振子的有阻尼自由振动频率记为ω'r
(2)n个正弦波输出完成后,正弦波发生器输出的标识信号S从低电平变为高电平,A/D转换器开始进行A/D采样,采样频率为fS,采样的数据按顺序存贮于数据缓冲区中,数据缓冲区的长度为X;同时,正弦波发生器输出m个波形完整但相位倒相的正弦波。此波形输出完成后,正弦波发生器不再输出信号。m为自然数,与超声波发射器的机械惯性相关,m≥1。
(3)DSP数据处理器从数据缓冲区读取A/D采样的数据进行FFT变换。从A/D采样开始的点依次取L个点D0、D1、D2···DL-1,1≤L≤2x,2x为FFT变换的数据点数,剩余的2x-L个点补0后,进行FFT变换。FFT变换后的结果就是2x点的复数,每一个点就对应着一个频率点。这个点的模值,就是该频率值下的幅度特性。相邻两点之间的频率间隔为第1个点是直流分量,它的模值是直流分量的2x倍。第2个点对应的频率为它的模值为该频率上信号幅值的2x-1倍;第3个点对应的频率为它的模值为该频率上信号幅值的2x-1倍;第4个点对应的频率为它的模值为该频率上信号幅值的2x-1倍……第2x个点对应的频率为它的模值为该频率上信号幅值的2x-1倍。因此选择合适的采样频率fS和合适的数据长度L,经过FFT变换,可以得到2x个频点的信号幅度值;选择合适的采样频率fS和合适的FFT变换数据点数2x,可以将超声波接收器接收的超声波频率ω和接收器振子的有阻尼振动频率ω'r准确对应于的整数倍,从而得到ω和ω'r两个频点信号的幅度值Aω0、Aω'r0
(4)DSP数据处理器从数据缓冲区第2点依次取L个点D1、D2、D3···DL,剩余的2x-L个点补0后,进行FFT变换,得到ω和ω'r两个频点信号的幅度值Aω1、Aω'r1,将Aω1、Aω'r1和Aω0、Aω'r0依次比较,判断这两个频点信号幅度的变化。
(5)依步骤(4)类推,DSP数据处理器从数据缓冲区第i点(i≥3)开始到X-L点依次作为起始点,读取L个点Di-1、Di、Di+1···Di+L-2(3≤i≤X-L),剩余的2x-L个点补0后,进行FFT变换,得到ω和ω'r两个频点信号的幅度值Aωi-1、Aω'ri-1,将Aωi-1、Aω'ri-1和Aωi-2、Aω'ri-2依次比较,判断这两个频点信号幅度的变化。
(6)控制器通过分别比较ω和ω'r两个频点信号的幅度值Aωi-1、Aω'ri-1和Aωi-2、Aω'ri-2,会发现Aωi逐渐减小,Aω'ri逐渐增大。控制器通过判别某点,即Aωi变为零的点和Aω'ri出现不为零的点,表明ω频率信号结束,ω'r频率信号出现,该点即为超声波接收器输出电信号的频率拐点。
本发明的有益效果是,本发明充分利用了超声波测距系统频率拐点的特性,通过A/D采样、DSP数据处理方法,可以准确、有效地解决有阻尼自由振动的超声波测距系统的频率拐点的检测问题,有效地克服了超声波接收振子的有阻尼自由振动的“拖尾”现象,准确、可靠地测量超声波的踱越时间,提高超声波测距的精度。实验表明,测量误差小于0.1mm,远远优于同类技术所能达到的技术水平。此方法对硬件电路要求不高、成本低、数据结构简单,处理速度快、应用广泛,易于普及与推广。
附图说明
图1是本发明的方法的流程框图。
图2是超声波发射波与接收信号在时间轴的对应关系图。
具体实施方式
下面根据附图详细描述本发明,本发明的目的和效果将变得更加明显。
本发明提供一种有阻尼自由振动超声波测距系统频率拐点判别方法,如图1所示,所述超声波测距系统包括控制器、正弦波发生器、功率放大电路、超声波发射器、超声波接收器、调理电路、A/D转换器、DSP数据处理器。
正弦波发生器输出的信号经功率放大器放大后激励超声波发射器输出超声波。正弦波发生器输出n个波形完整且连续的正弦波(图2中,n=15),其输出波形的个数n是可控制的,在其输出信号的同时,也输出一个标识信号S,当其输出信号时,S=0,当其不输出信号时,S=1。假定输出波形的时间为T1,不输出波形的时间为T2。为减小正弦波发生器输出信号停止后,发射器的余振对接收器的影响,正弦波发生器输出n个波形后,附加输出了m个波形完全倒相的信号周期,保证发射器不产生拖尾现象(图2中,m=1)。
超声波接收器接收超声波发射器输出的超声波,并以电信号的形式输出。超声波接收器输出的微弱电信号包含了超声波接收器振子受迫振动状态和有阻尼自有振动状态两种状态的信号,超声波接收器接收的超声波频率记为ω,接收器振子的有阻尼自由振动频率记为ω'r
接收器输出的微弱电信号经信号调理电路放大后,进入A/D转换器进行A/D采样。A/D转换器以TMS320F28335型数字信号处理器为核心,TMS320F28335型数字信号处理器具有150MHz的高速处理能力,具备32位浮点处理单元,6个DMA通道支持ADC和12位16通道ADC。采样率为fs,(图2中fs=128kHz),采样点数为X,结果存储于数据缓存区中。
DSP数据处理器从数据缓冲区中读取A/D采样的数据。假设A/D采样后,接收到信号的峰值为H,那么FFT的结果的每个点的模值就是H的2x-1倍。那么,第一个点就是直流分量,它的模值就是直流分量的2x倍。而每个点的相位就是在该频率下的信号的相位。第一个点表示直流分量(即0Hz),而最后一个点2x的再下一个点则表示采样频率fs,这中间的波形被2x-1个点平均分成2x等份,每个点的频率依次递增。假设某点n所表示的频率为:fn=(n-1)fs/2x。由上面的公式可以看出,fn所能分辨到频率为fs/2x,如果采样频率为1024Hz,采样点数为1024点,则可以分辨到1Hz。1024Hz的采样率采样1024点,刚好是1秒,换句话说就是,采样1秒时间的信号并做FFT,则结果可以分析到1Hz,如果采样2秒时间的信号并做FFT,则结果可以分析到0.5Hz。频率分辨率和采样时间是倒数关系,如果增加采样点数,即采样时间,就可以提高频率分辨力。假设FFT之后某点n用复数a+bi表示,其中,复数a+bi的实部和虚部,分别表示信号在n点幅度和相位。那么这个复数的模就是相位就是假设我们从A/D采样开始的点依次取L个点D0、D1、D2···DL-1,1≤L≤2x,2x为FFT变换的数据点数,剩余的2x-L个点补0后,进行FFT变换。FFT变换后的结果就是2x点的复数,每一个点就对应着一个频率点。这个点的模值,就是该频率值下的幅度特性。相邻两点之间的频率间隔为第1个点是直流分量,它的模值是直流分量的2x倍。第2个点对应的频率为它的模值为该频率上信号幅值的2x - 1倍;第3个点对应的频率为它的模值为该频率上信号幅值的2x - 1倍;第4个点对应的频率为它的模值为该频率上信号幅值的2x-1倍……第2x个点对应的频率为它的模值为该频率上信号幅值的2x-1倍。因此选择合适的采样频率fS和合适的数据长度L,经过FFT变换,可以得到2x个频点的信号幅度值;选择合适的采样频率fS和合适的FFT变换数据点数2x,可以将ω和ω'r准确对应于的整数倍,从而得到ω和ω'r两个频点信号的幅度值Aω0、Aω'r0。DSP数据处理器从数据缓冲区第3点开始到X-L点依次作为起始点,依次取L个点Di、Di+1、Di+2···Di+L-1(3≤i≤X-L)剩余的2x-L个点补0后,进行FFT变换,得到ω和ω'r两个频点信号的幅度值Aωi、Aω'ri,将Aωi、Aω'ri和Aωi-1、Aω'ri-1依次比较,判断这两个频点信号幅度的变化。控制器读取通过依次比较ω和ω'r两个频点信号的幅度值Aωi、Aω'ri和Aωi-1、Aω'ri-1,会发现Aωi逐渐减小,Aω'ri逐渐增大。控制器通过判别某一点,即Aωi变为零的点和Aω'ri出现不为零的点,该点即为频率的拐点。
根据以上的结果,就可以准确地判别出有阻尼自由振动的超声波测距系统的频率拐点。利用有阻尼自由振动的超声波测距系统的频率拐点判别方法可以准确地得出超声波传输的渡越时间,从而提高超声波测距系统的精度。
图1中,通过控制器来驱动正弦波发生器,正弦波发生器输出的信号经功率放大器放大后激励超声波发射器输出超声波。超声波发射器通过传输后,将超声波信号传给接收器。接收器输出的微弱电信号经信号调理电路放大后进行A/D采样,将A/D采样后得到的数据经过DSP数据处理器进行FFT数据分析,分析的结果传回到控制器里的存储器。
图2中,分别描述了正弦波发生器输出的信号、标识信号、超声波信号、超声波接收器输出信号、A/D采样信号和FFT变换后信号在时间轴上对应的关系。由图2中的(5)可见,A/D采样可以将超声波中的模拟信号转化为数字信号。由图2中的(6)可见,对A/D采样后的点分组后,逐个组进行FFT变换后,可以得到超声波信号每个点的幅度谱和相位谱,根据幅度谱中某一点和其相邻点的幅度值对应的频率差别,可以明确地得出超声波信号相邻两个点的频率变化,从而判别出频率的拐点。
实施例:本实施例有阻尼自由振动超声波测距系统频率拐点判别方法,具体包括以下步骤:
(1)正弦波发生器输出n个波形完整且连续的正弦波信号,该正弦波信号通过功率放大电路放大后激励超声波发射器发出超声波,同时超声波接收器开始接收超声波,超声波接收器输出的微弱电信号经信号调理电路放大后,输出至A/D转换器。
(2)n个正弦波输出完成后,正弦波发生器输出的标识信号S从低电平变为高电平,A/D转换器开始进行A/D采样,采样频率为fS(fs=128kHz),采样的数据按顺序存贮于数据缓冲区中,数据缓冲区的长度为X;同时,正弦波发生器输出m个波形完整但相位倒相的正弦波。此波形输出完成后,正弦波发生器不再输出信号。m为自然数,与超声波发射器的机械惯性相关,m≥1。
(3)DSP数据处理器从数据缓冲区读取A/D采样的数据进行FFT变换。从A/D采样开始的点依次取40个点D0、D1、D2···D39,2x为FFT变换的数据点数,剩余的2x-40个点补0后,进行FFT变换。FFT变换后的结果就是2x点的复数,每一个点就对应着一个频率点。这个点的模值,就是该频率值下的幅度特性。相邻两点之间的频率间隔为第1个点是直流分量,它的模值是直流分量的2x倍。第2个点对应的频率为它的模值为该频率上信号幅值的2x-1倍;第3个点对应的频率为它的模值为该频率上信号幅值的2x-1倍;第4个点对应的频率为它的模值为该频率上信号幅值的2x-1倍……第2x个点对应的频率为它的模值为该频率上信号幅值的2x-1倍。因此选择合适的采样频率fS和合适的数据长度L,经过FFT变换,可以得到2x个频点的信号幅度值;选择合适的采样频率fS和合适的FFT变换数据点数2x,可以将ω和ω'r准确对应于的整数倍,从而得到ω和ω'r两个频点信号的幅度值Aω0、Aω'r0
(4)DSP数据处理器从数据缓冲区第2点依次取40个点D1、D2、D3···D40,剩余的2x-40个点补0后,进行FFT变换,得到ω和ω'r两个频点信号的幅度值Aω1、Aω'r1,将Aω1、Aω'r1和Aω0、Aω'r0依次比较,判断这两个频点信号幅度的变化。
(5)依步骤(4)类推,DSP数据处理器从数据缓冲区第3点开始到X-39点依次作为起始点,读取40个点Di-1、Di、Di+1、Di+2···Di+38(3≤i≤X-39)剩余的2x-40个点补0后,进行FFT变换,得到ω和ω'r两个频点信号的幅度值Aωi-1、Aω'ri-1,将Aωi-1、Aω'ri-1和Aωi-2、Aω'ri-2依次比较,判断这两个频点信号幅度的变化。
(6)控制器通过分别比较ω和ω'r两个频点信号的幅度值Aωi-1、Aω'ri-1和Aωi-2、Aω'ri-2,会发现Aωi逐渐减小,Aω'ri逐渐增大。控制器通过判别某点,即Aωi变为零的点和Aω'ri出现不为零的点,表明ω频率信号结束,ω'r频率信号出现,该点即为超声波接收器输出电信号的频率拐点。
本发明专利的核心技术要点之一在于采用了FFT快速傅里叶变换,根据幅度谱中某一固定频率对应的点在不同周期内相对位置是否发生偏移,可以明确地得出超声波信号相邻两个周期的频率变化,从而判别出频率的拐点。
本发明专利的核心技术要点之二在于采用了DSP数据处理器,DSP数据处理器能够对每一个A/D采样后超声波接收器输出的数字信号进行FFT变换,处理速度快,效率高。
本发明专利的核心技术要点之三在于采用了TMS320F28335型数字信号处理器进行A/D转换,TMS320F28335型数字信号处理器是浮点型控制器,与以往的定点型控制器相比,该器件的精度高,成本低,功耗小,性能高,外设集成度高,数据存储量大,A/D转换更精确快速。
本发明充分利用了有阻尼自由振动的超声波测距系统频率拐点的特性,可以准确、有效地测量超声波的踱越时间,提高超声波测距的精度。实验表明,采用此方法测量误差小于0.1m远远优于同类技术所能达到的技术水平。

Claims (2)

1.一种有阻尼自由振动超声波测距系统频率拐点判别方法,其特征在于,该方法包括以下步骤:
(1)正弦波发生器输出n个波形完整且连续的正弦波信号,该正弦波信号通过功率放大电路放大后激励超声波发射器发出超声波,同时超声波接收器开始接收超声波,超声波接收器输出的微弱电信号经信号调理电路放大后,输出至A/D转换器。超声波接收器输出的微弱电信号包含超声波接收器振子受迫振动状态和有阻尼自有振动状态两种状态的信号,超声波接收器接收的超声波频率记为ω,接收器振子的有阻尼自由振动频率记为ω'r
(2)n个正弦波输出完成后,正弦波发生器输出的标识信号S从低电平变为高电平,A/D转换器开始进行A/D采样,采样频率为fS,采样的数据按顺序存贮于数据缓冲区中,数据缓冲区的长度为X;同时,正弦波发生器输出m个波形完整但相位倒相的正弦波。此波形输出完成后,正弦波发生器不再输出信号。m为自然数,与超声波发射器的机械惯性相关,m≥1。
(3)DSP数据处理器从数据缓冲区读取A/D采样的数据进行FFT变换。从A/D采样开始的点依次取L个点D0、D1、D2···DL-1,1≤L≤2x,2x为FFT变换的数据点数,剩余的2x-L个点补0后,进行FFT变换。FFT变换后的结果就是2x点的复数,每一个点对应一个频率点。这个点的模值,就是该频率值下的幅度特性。相邻两点之间的频率间隔为第1个点是直流分量,它的模值是直流分量的2x倍。第n个点对应的频率为它的模值为该频率上信号幅值的2x-1倍,1<n≤2x;因此选择合适的采样频率fS和合适的数据长度L,经过FFT变换,可以得到2x个频点的信号幅度值;选择合适的采样频率fS和合适的FFT变换数据点数2x,可以将超声波接收器接收的超声波频率ω和接收器振子的有阻尼振动频率ω'r准确对应于的整数倍,从而得到ω和ω'r两个频点信号的幅度值Aω0、Aω'r0
(4)DSP数据处理器从数据缓冲区第2点依次取L个点D1、D2、D3···DL,剩余的2x-L个点补0后,进行FFT变换,得到ω和ω'r两个频点信号的幅度值Aω1、Aω'r1,将Aω1、Aω'r1和Aω0、Aω'r0依次比较,判断这两个频点信号幅度的变化。
(5)依步骤(4)类推,DSP数据处理器从数据缓冲区第i点(i≥3)开始到X-L点依次作为起始点,读取L个点Di-1、Di、Di+1···Di+L-2(3≤i≤X-L),剩余的2x-L个点补0后,进行FFT变换,得到ω和ω'r两个频点信号的幅度值Aωi-1、Aω'ri-1,将Aωi-1、Aω'ri-1和Aωi-2、Aω'ri-2依次比较,判断这两个频点信号幅度的变化。
(6)控制器通过分别比较ω和ω'r两个频点信号的幅度值Aωi-1、Aω'ri-1和Aωi-2、Aω'ri-2,会发现Aωi逐渐减小,Aω'ri逐渐增大。控制器通过判别某点,即Aωi变为零的点和Aω'ri出现不为零的点,表明ω频率信号结束,ω'r频率信号出现,该点即为超声波接收器输出电信号的频率拐点。
2.根据权利要求1所述的一种有阻尼自由振动超声波测距系统频率拐点判别方法,其特征在于,所述步骤2中,A/D转换器的采样具体为:A/D转换器启动后首先清空转换结果数据缓冲区,以免上次的采样数据对当前造成不必要的误差影响;将A/D转换器转换分辨率设置为8位,采样频率为fs=128kHz,结果存储在数据缓存区中。
CN201610778428.8A 2016-08-30 2016-08-30 一种有阻尼自由振动超声波测距系统频率拐点判别方法 Expired - Fee Related CN106371098B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610778428.8A CN106371098B (zh) 2016-08-30 2016-08-30 一种有阻尼自由振动超声波测距系统频率拐点判别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610778428.8A CN106371098B (zh) 2016-08-30 2016-08-30 一种有阻尼自由振动超声波测距系统频率拐点判别方法

Publications (2)

Publication Number Publication Date
CN106371098A true CN106371098A (zh) 2017-02-01
CN106371098B CN106371098B (zh) 2018-12-14

Family

ID=57902529

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610778428.8A Expired - Fee Related CN106371098B (zh) 2016-08-30 2016-08-30 一种有阻尼自由振动超声波测距系统频率拐点判别方法

Country Status (1)

Country Link
CN (1) CN106371098B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111947769A (zh) * 2020-06-30 2020-11-17 核动力运行研究所 共振波频率确定方法及装置
CN112291676A (zh) * 2020-05-18 2021-01-29 珠海市杰理科技股份有限公司 抑制音频信号拖尾的方法及系统、芯片、电子设备
CN113740856A (zh) * 2021-07-21 2021-12-03 国创移动能源创新中心(江苏)有限公司 基于超声波波形起振特性的测距方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1233281A1 (en) * 2001-02-14 2002-08-21 ISE Ingegneria dei Sistemi Elettronici S.r.l. Apparatus and method for reducing the effects of a fluid motion on measuring distances by means of acoustic waves
CN103777202A (zh) * 2013-12-16 2014-05-07 浙江工商大学 基于有阻尼自由振动的超声波测量距离的方法
CN104101870A (zh) * 2014-05-27 2014-10-15 浙江工商大学 一种频域调制式超声波测距系统频率拐点判别方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1233281A1 (en) * 2001-02-14 2002-08-21 ISE Ingegneria dei Sistemi Elettronici S.r.l. Apparatus and method for reducing the effects of a fluid motion on measuring distances by means of acoustic waves
CN103777202A (zh) * 2013-12-16 2014-05-07 浙江工商大学 基于有阻尼自由振动的超声波测量距离的方法
CN104101870A (zh) * 2014-05-27 2014-10-15 浙江工商大学 一种频域调制式超声波测距系统频率拐点判别方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MIKE HAZAS 和 ANDY HOPPER: "Broadband ultrasonic location systems for improved indoor positioning", 《IEEE TRANSACTIONS ON MOBILE COMPUTING》 *
陈建 等: "一种高精度超声波到达时刻的检测方法", 《仪器仪表学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112291676A (zh) * 2020-05-18 2021-01-29 珠海市杰理科技股份有限公司 抑制音频信号拖尾的方法及系统、芯片、电子设备
CN112291676B (zh) * 2020-05-18 2021-10-15 珠海市杰理科技股份有限公司 抑制音频信号拖尾的方法及系统、芯片、电子设备
CN111947769A (zh) * 2020-06-30 2020-11-17 核动力运行研究所 共振波频率确定方法及装置
CN113740856A (zh) * 2021-07-21 2021-12-03 国创移动能源创新中心(江苏)有限公司 基于超声波波形起振特性的测距方法和装置
CN113740856B (zh) * 2021-07-21 2023-12-22 国创移动能源创新中心(江苏)有限公司 基于超声波波形起振特性的测距方法和装置

Also Published As

Publication number Publication date
CN106371098B (zh) 2018-12-14

Similar Documents

Publication Publication Date Title
CN105424802A (zh) 一种复合绝缘子缺陷的超声导波检测系统及其检测方法
CN101726539B (zh) 利用声表面波器件测试气体浓度的方法
CN2842393Y (zh) 材料内应力超声波测量装置
CN102043015B (zh) 长距离探测钢轨轨底缺陷的超声导波装置及方法
CN105066918A (zh) 超声水下目标测厚系统及测厚方法
CN106371098B (zh) 一种有阻尼自由振动超声波测距系统频率拐点判别方法
CN102507655A (zh) 基于压电智能骨料的钢管混凝土结构密实性监测方法
CN102338776A (zh) 应力波在木材中的传播时间测量方法以及无损检测系统
TWI225152B (en) Method for non-destructive stress wave testing of wood
CN106404247A (zh) 基于瑞利波模式的声表面波双谐振器一体化扭矩传感器
CN106950282A (zh) 一种基于纵向超声导波的玻璃钢芯棒缺陷检测方法及系统
CN103308012A (zh) 混凝土路面厚度检测系统及检测方法
CN103323525A (zh) 检测超级奥氏体不锈钢晶间腐蚀的非线性超声系统及方法
CN106802202A (zh) 一种测量各向异性材料平面应力的方法
CN110361324A (zh) 一种结合激光冲击波和兰姆波的复合材料结合力在线快速检测装置及方法
CN104764522B (zh) 一种超声功率测量方法及装置
CN107121497A (zh) 基于Duffing系统的随机共振特性的超声导波检测方法
CN111157065A (zh) 气体超声流量计超声波信号传输回路中声延时测量方法
CN109490417B (zh) 一种金属材料平面各向异性超声检测方法
CN110672714A (zh) 一种非接触式桥梁腐蚀钢筋磁感应检测方法
CN104267103A (zh) 小模数齿轮齿根残余应力的晶片式超声无损检测探头
CN100501394C (zh) 对物体进行超声波检测的方法
CN103777202A (zh) 基于有阻尼自由振动的超声波测量距离的方法
CN102519551A (zh) 一种基于伪码相关技术的超声波液位测量仪
KR100844173B1 (ko) 과실의 물성 측정을 위한 초음파 탐촉자

Legal Events

Date Code Title Description
C06 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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181214

Termination date: 20210830

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