CN112504434A - 一种物体相对运动速度和空气声波衰减系数的测量系统及测量方法 - Google Patents

一种物体相对运动速度和空气声波衰减系数的测量系统及测量方法 Download PDF

Info

Publication number
CN112504434A
CN112504434A CN202011227342.9A CN202011227342A CN112504434A CN 112504434 A CN112504434 A CN 112504434A CN 202011227342 A CN202011227342 A CN 202011227342A CN 112504434 A CN112504434 A CN 112504434A
Authority
CN
China
Prior art keywords
signal
acoustic
attenuation coefficient
air
segmented
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
CN202011227342.9A
Other languages
English (en)
Other versions
CN112504434B (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.)
Changzhou University
Original Assignee
Changzhou 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 Changzhou University filed Critical Changzhou University
Priority to CN202011227342.9A priority Critical patent/CN112504434B/zh
Publication of CN112504434A publication Critical patent/CN112504434A/zh
Application granted granted Critical
Publication of CN112504434B publication Critical patent/CN112504434B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P3/00Measuring linear or angular speed; Measuring differences of linear or angular speeds
    • G01P3/64Devices characterised by the determination of the time taken to traverse a fixed distance
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/245Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by measuring transit time of acoustical waves

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种物体相对运动速度和空气声波衰减系数的测量系统及测量方法,所述系统以声信号发生器为一端,由声信号采集器和信号发送装置组成另一端,上述两端中一个为固定端,另一个为运动端,运动端或固定端产生的具有平稳特征的声信号,经信号发生器发送、并被安装在固定端或运动端的声信号采集器接收后,通过信号发送装置和信号接收装置被传送到信号处理与计算系统进行处理,计算得到物体的相对运动速度和空气声波衰减系数。本发明可以实现声信号的高速采集,移动端相对运动速度的跟踪测量和空气声波衰减系数的准确计算,并且使用成本低,方案多样,实施方便。

Description

一种物体相对运动速度和空气声波衰减系数的测量系统及测 量方法
技术领域
本发明涉及一种物体相对运动速度和空气声波衰减系数的测量系统及测量方法。
背景技术
人们一般采用专用的测量工具来实现物体运动速度的测量。例如,计时测速、光电门测速、频闪照相测速、平抛运动测速、运动轨迹测速、速度计测速、旋转齿轮测速、旋转转筒测速、速度选择器测速、电流计测速、电磁流量计测速、单摆测速、脉冲信号测速、多普勒效应测速、超声波测速、旋转光源测速、激光散斑测速等,都是非常有效的测量物体运动速度的方法。空气声传播速度一般服从如下规律:
Figure BDA0002764010160000011
式中,T为温度,℃;p为空气压强,Pa;pw为相对湿度与对应温度饱和蒸汽压的乘积,Pa;c0为温度为0℃时干空气的声传播速度。
声音在传播过程中会产生衰减,声强的衰减规律服从如下公式:
Id=I0e-αS
式中,Id为在离声源距离为S处测量得到的声强,单位为Pa;I0为在声源的声强,单位为Pa;S为测点离声源的距离,单位为m;α为空气声传播衰减系数,单位为m-1
声波传播速度的测量方法已经发展的很成熟,常见的有时差法、共振干涉法、相位比较法等声速测量方法。
声音的衰减具有频散特性,也就是说声信号中不同频率的振动成分在传播过程中衰减的速率是不同的,高频成分要比低频成分衰减得快一些。材料的衰减系数是一个基本的物理参数,能够反映物体自身的内部结构和组成特性,在许多领域具有十分重要的应用价值。人们也往往利用材料的衰减系数进行材料组成、内部损伤和结构缺陷的探测。但目前声波参数的测量多数是在实验室完成的,例如,声波速度和声衰减系数的测量已经成为大学物理实验的一个必训项目,但应用于生产实际的测量需要采用专用的工具,受限于测量手段的局限性,限制了声波衰减规律的应用。
如果能够直接测量得到物体运动期间固定端与运动端之间的相对距离,则可以采用传统的衰减系数计算方法由物体运动信号源的振幅A随距离S的变化求得空气衰减系数α。继而根据运动期间时间的记录,计算得到测点与运动物体之间的相对运动速度。
发明内容
本发明要解决的技术问题是:针对现有技术中存在的空气声波参数测量实时性差、应用受限等不足,本发明提供一种实时性强、对能够在运动过程中产生平稳周期声信号的物体进行相对运动速度和空气衰减系数测量的系统及其测量方法。
本发明解决其技术问题所采用的技术方案是:一种物体相对运动速度和空气声波衰减系数测量系统,包括声信号发生器,声信号采集器,信号发送装置,信号接收装置以及信号处理与计算系统,以声信号发生器为一端,由声信号采集器和信号发送装置组成另一端,上述两端中一个为固定端,另一个为运动端。
所述的信号发送装置与信号接收装置之间为有线或无线连接,信号接收装置与信号处理与计算系统之间为有线连接,运动端或固定端产生的具有平稳特征的声信号,经信号发生器发送、并被安装在固定端或运动端的声信号采集器接收后,通过信号发送装置和信号接收装置,声信号被传送到信号处理与计算系统进行处理,继而计算得到物体的相对运动速度和空气声波衰减系数。
进一步地,所述的声信号采集器的信号采样率、信号发送装置和信号接收装置的接发信号速率均不低于44kHz,且能够采集全频段的声信号。
一种采用上述测量系统进行物体相对运动速度和空气声波衰减系数测量的方法,具有如下步骤:
S1、运动端相对于固定端以一定的速度u(t)移动的同时,声信号发生器处产生一个平稳的周期声信号y0(t),这一声信号被声信号采集器接收,并以fs的采样率采集得到采样数据点数量为M的y1(t)信号;y1(t)信号经信号发送装置和信号接收装置被传送到信号处理与计算系统中。
S2、从y1(t)信号中截取一段小于1个且不小于1/2个平稳周期的特征信号y1d(t)信号,y1d(t)信号的长度为N个数据点(N<M)。在y1(t)信号上,沿着y1(t)信号计算y1(t)与y1d(t)的互相关性系数rxy(t),rxy(t)的总长度为M-N+1。
S3、在信号段长度不小于N的限制条件下,获取rxy(t)在各个分段上的rxy(t)最大值所对应的数据点位置ti及对应rxy(ti)数值,并以ti为参照点对y1(t)信号进行分段处理,获得对应的分段y1(ti,i)信号。
S4、读取分段y1(ti,i)信号的数据点数量、峰值数值、峰谷数值、峰值位置和峰谷位置信息。
S5、计算分段y1(ti,i)信号的绝对均值、峰值指标、峭度指标、歪度指标和裕度指标数值。
S6、在信号段长、峰值数值、峰值位置、绝对均值、峰值指标、峭度指标、歪度指标和裕度指标参数组成的多维空间考察数据点的分布情况,删除含异常的数据点。
S7、在经过上述步骤S1~S6的处理,获得了分段y1(ti,i)信号及相应的信号段长、峰值数值、峰值位置、绝对均值、峰值指标、峭度指标、歪度指标和裕度指标统计参数后,根据以下两种情况,进行针对性的处理:
S7-1、在已知空气声传播速度c的前提下,利用式(1)进行空气声传播衰减系数的计算:
Figure BDA0002764010160000041
式中,α为空气声传播衰减系数,单位为m-1;K为分段y1(ti,i)信号前后连贯结对数量;Ai、Ai-1为对应第i对前后结对分段信号的峰值振幅数值;Li、Li-1为第i对前后结对分段信号的分段数据点数量;fP为信号的采样率,单位为Hz;c是空气的声传播速度,单位为m/s。Bf为与A的取值相关的系数,如果A取声信号振幅,则Bf=1;如果A为信号段的绝对均值、有效值等参数,则需要通过校正实验来确定Bf的数值。
S7-2、在未知空气声传播速度c的情况下,利用式(2)进行空气声传播速度和空气声传播衰减系数的优化搜索计算:
Figure BDA0002764010160000042
式中,β为系数,β=1α,单位为m。
S8、在计算得到空气的声传播衰减系数α之后,运动端相对于固定端的相对运动速度ui的计算公式如下:
Figure BDA0002764010160000043
式中,Ai、Ai-1:对应第i对前后结对分段信号的峰值振幅数值;Li、Li-1:第i对前后结对分段信号的分段数据点数量。定义ui相对接近运动为负值,远离运动为正值。
为了保证计算的可靠性,上述式(1)、(2)和式(3)中的A数值,也可以采用分段y1(ti,i)信号的绝对均值、有效值等替代。在这种情况下,式(1)中的系数Bf数值需要通过适当的校正实验来确定。
一种采用上述测量系统进行物体相对运动速度和空气声波衰减系数测量的方法,具有如下步骤:
(1)、所述运动端相对于固定端以一定的速度u(t)移动的同时,声信号发生器处产生一个平稳的周期性声信号y0(t),这一声信号被声信号采集器接收,并以fs的采样率采集得到采样数据点数量为M的y1(t)信号;y1(t)信号经信号发送装置和信号接收装置被传送到信号处理与计算系统中。
(2)、从y1(t)信号中截取一段小于1个且不小于1/2个平稳周期的特征信号y1d(t)信号,y1d(t)信号的长度为N个数据点(N<M),在y1(t)信号上,沿着y1(t)信号计算y1(t)与y1d(t)的互相关性系数rxy(t),rxy(t)的总长度为M-N+1。
(3)、在信号段长度不小于N的限制条件下,获取rxy(t)在各个分段上的rxy(t)最大值所对应的数据点位置ti及对应rxy(ti)数值,并以ti为参照点对y1(t)信号进行分段处理,获得对应的分段y1(ti,i)信号。
(4)、读取分段y1(ti,i)信号的数据点数量、峰值数值、峰谷数值、峰值位置和峰谷位置等信息,并进行异常数据点的删除处理。
(5)、用分段y1(ti,i)信号来构造信号总长度相同的循环分段信号,并通过FFT变换计算得到循环分段信号的主频数值fi
(6)、利用下式计算在相对运动速度为0时循环分段信号的主频f0
Figure BDA0002764010160000061
Figure BDA0002764010160000062
式中,f0为无频移(或相对运动速度为0)时循环分段信号的主频率,单位为Hz;Ai、Ai-1为对应第i对前后结对分段信号的峰值振幅数值;Δti为第i对前后结对分段信号的时间间隔差,单位为s;fi+1、fi为第i+1个和第i个循环分段信号的主频率,单位为Hz。
在计算得到每个前后结对分段信号的f0后,取平均值作为相对运动速度为0时循环分段信号的主频率。
(7)、在空气声传播速度c已知的情况下,利用式(5)来计算运动端相对于固定端的相对运动速度:
Figure BDA0002764010160000063
式中,c为空气声传播速度,单位为m/s。
(8)、在计算得到相对运动速度ui之后,空气声传播衰减系数α的计算公式如下:
Figure BDA0002764010160000064
式中,α为空气声传播衰减系数,单位为m-1;K为分段y1(ti,i)信号前后连贯结对数量;fp为声波信号采样率,单位为Hz。
本发明的有益效果是:本发明采用声信号测量和处理的途径,从实测信号中提取声波衰减系数和物体运动速度的信息,极大地简化了测量系统对软硬件设备的需求;通过结合无线或有线收发装置将信号数据传送到芯片处理系统,由芯片系统进行信号的分析与处理,可以实时获得空气衰减系数和物体运动相对速度的信息,有效地提高了测量效率和信息使用效率。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是本发明的系统结构示意图。
图2是本发明所述测量方法的流程图。
图3是实施例一中声信号发生器端固定、信号采集与接收端运动时的声信号测量实例图。
图4是实施例一中声信号发生器端运动、信号采集与接收端固定时的声信号测量实例图。
图5是实施例一中声信号发生器端固定、信号采集与接收端运动时的分段声信号图。
图6是实施例一中声信号发生器端运动、信号采集与接收端固定时的分段声信号图。
图7是实施例二中声信号发生器端运动、信号采集与接收端固定时的声信号测量实例。
图8是实施例二中声信号发生器端运动、信号采集与接收端固定时的分段声信号。
图中:1.声信号发生器,2.声信号采集器,3.信号发送装置,4.信号接收装置,5.信号处理与计算系统。
具体实施方式
现在结合附图对本发明作进一步详细的说明。这些附图均为简化的示意图,仅以示意方式说明本发明的基本结构,因此其仅显示与本发明有关的构成。
如图1所示,一种物体相对运动速度和空气声波衰减系数的测量系统,包括声信号发生器1、声信号采集器2、信号发送装置3、信号接收装置4及信号处理与计算系统5,以声信号发生器1为一端,由声信号采集器2和信号发送装置3组成另一端,这两端中一个为固定端,另一个为运动端,所述的信号发送装置3与信号接收装置4为有线或无线连接,信号接收装置4与信号处理与计算系统5为有线连接。
在声信号发生器1中产生的声信号,通过空气传播,被声信号采集器2接收,通过信号发送装置3和信号接收装置4后,声信号被传送到信号处理与计算系统5进行处理,计算得到物体的相对运动速度和空气声波衰减系数。
所述的声信号采集器2的信号采样率、信号发送装置3与信号接收装置4的接发信号速率均不低于44kHz,且能够采集全频段的声信号。
实施例一,图3是固定不动的电脑上发出的数字化声信号,通过运动着的声传感器采集声信号,并通过NIWLS9234多通道无线采集仪进行传输和接收的声信号测量实例;图4是某飞行动物发出的鸣声信号,由固定不动的声传感器接收,并采用声信号采集仪记录得到的声信号实例。
一种采用上述测量系统进行物体相对运动速度和空气声波衰减系数测量的方法,测量流程如图2所示,具有如下步骤:
S1、运动端相对于固定端以一定的速度u(t)移动的同时,声信号发生器1处产生一个平稳周期声信号y0(t),这一声信号被声信号采集器2接收,并以fs的采样率采集得到采样数据点数量为M的y1(t)信号;y1(t)信号经信号发送装置3和信号接收装置4被传送到信号处理与计算系统5中。
S2、从y1(t)信号中截取一段小于1个且不小于1/2个平稳周期的特征信号y1d(t)信号,y1d(t)信号的长度为N个数据点(其中N<M),在y1(t)信号上,沿着y1(t)信号计算y1(t)与y1d(t)的互相关性系数rxy(t),rxy(t)的总长度为M-N+1。
S3、在信号段长度不小于N的限制条件下,获取rxy(t)在各个分段上的rxy(t)最大值所对应的数据点位置ti及对应rxy(ti)数值,并以ti为参照点对y1(t)信号进行分段处理,获得对应的分段y1(ti,i)信号。其中,图5就是对图3的声信号进行信号分段的处理结果,图6是对应图4的信号进行信号分段的处理结果。
S4、读取图5和图6中每个分段y1(ti,i)信号的数据点数量(信号段长)、峰值数值、峰谷数值、峰值位置和峰谷位置等信息。
S5、计算分段y1(ti,i)信号的绝对均值、峰值指标、峭度指标、歪度指标和裕度指标等数值。
S6、在信号段长、峰值数值、峰值位置、绝对均值、峰值指标、峭度指标、歪度指标和裕度指标等参数组成的多维空间考察数据点的分布情况,删除含异常的数据点。
S7、在经过步骤S1~S6的处理,获得了分段y1(ti,i)信号及相应的信号段长、峰值数值、峰值位置、绝对均值、峰值指标、峭度指标、歪度指标和裕度指标等统计参数后,根据以下两种情况,进行针对性的处理:
S7-1、在已知空气声传播速度c的前提下,利用式(1)进行空气声传播衰减系数的计算:
Figure BDA0002764010160000091
式中,α为空气声传播衰减系数,单位为m-1;K为分段y1(ti,i)信号前后连贯结对数量;Ai、Ai-1为对应第i对前后结对分段信号的峰值振幅数值;Li、Li-1为第i对前后结对分段信号的分段数据点数量;fP为信号的采样率,单位为Hz;c是空气的声传播速度,单位为m/s。Bf为与A的取值相关的系数,如果A取声信号振幅,则Bf=1;如果A为信号段的绝对均值、有效值等参数,则需要通过校正实验来确定Bf的数值。
S7-2、在未知空气声传播速度c的情况下,利用式(2)进行空气声传播速度和空气声传播衰减系数的优化搜索计算:
Figure BDA0002764010160000101
式中,β为系数,β=1/α,单位为m。
S8、在计算得到空气的声传播衰减系数α之后,运动端相对于固定端的相对运动速度ui的计算公式如下:
Figure BDA0002764010160000102
式中,Ai、Ai-1为对应第i对前后结对分段信号的峰值振幅数值;Li、Li-1为第i对前后结对分段信号的分段数据点数量。定义ui相对接近运动为负值,远离运动为正值。
为了保证计算的可靠性,式(1)、(2)和式(3)中的A数值,也可以采用分段y1(ti,i)信号的绝对均值、有效值等替代,在这种情况下,式(1)中的系数Bf数值需要通过适当的校正实验来确定。
实施例二,图7是飞行动物果蝇发出的鸣声信号,由固定不动的声传感器接收,并采用声信号采集仪记录得到的声信号实例。
一种采用上述测量系统进行物体相对运动速度和空气声波衰减系数测量的方法,测量流程如图2所示,具有如下步骤:
(1)、所述运动端相对于固定端以一定的速度u(t)移动的同时,声信号发生器1处产生一个平稳周期声信号y0(t),这一声信号被声信号采集器2接收,并以fs的采样率采集得到采样数据点数量为M的y1(t)信号;y1(t)信号经信号发送装置3和信号接收装置4被传送到信号处理与计算系统5中;
(2)、从y1(t)信号中截取一段小于1个且不小于1/2个平稳周期的特征信号y1d(t)信号,y1d(t)信号的长度为N个数据点(N<M)。在y1(t)信号上,沿着y1(t)信号计算y1(t)与y1d(t)的互相关性系数rxy(t),rxy(t)的总长度为M-N+1;
(3)、在信号段长度不小于N的限制条件下,获取rxy(t)在各个分段上的rxy(t)最大值所对应的数据点位置ti及对应rxy(ti)数值,并以ti为参照点对y1(t)信号进行分段处理,获得对应的分段y1(ti,i)信号。其中,图8就是对应图7的信号进行信号分段的处理结果;
(4)、读取图8中每个分段y1(ti,i)信号的数据点数量(信号段长)、峰值数值、峰谷数值、峰值位置和峰谷位置等信息,并进行异常数据点的删除处理。
(5)、用分段y1(ti,i)信号来构造信号总长度相同的循环分段信号,并通过FFT变换计算得到循环分段信号的主频数值fi
(6)、利用下式计算在相对运动速度为0时循环分段信号的主频f0
Figure BDA0002764010160000111
Figure BDA0002764010160000112
式中,f0为无频移(或相对运动速度为0)时循环分段信号的主频率,单位为Hz;Ai、Ai-1为对应第i对前后结对分段信号的峰值振幅数值;Δti为第i对前后结对分段信号的时间间隔差,单位为s;fi+1、fi为第i+1个和第i个循环分段信号的主频率,单位为Hz。
在计算得到每个前后结对分段信号的f0后,取平均值作为相对运动速度为0时循环分段信号的主频率。
(7)、在空气声传播速度c已知的情况下,利用式(5)来计算运动端相对于固定端的相对运动速度:
Figure BDA0002764010160000121
式中,c为空气声传播速度,单位为m/s。
(8)、在计算得到相对运动速度ui之后,空气声传播衰减系数α的计算公式如下:
Figure BDA0002764010160000122
式中,α为空气声传播衰减系数,单位为m-1;K为分段y1(ti,i)信号前后连贯结对数量;fp为声波信号采样率,单位为Hz。
本发明根据固定端和运动端之间,在未知两端相对距离的前提下,通过测量与处理对方物体上发出的特征信号,来获取物体相对运动速度,以及所测区域空气声传播速度和空气衰减系数,从而实现声信号的高速采集,移动端相对运动速度的跟踪测量和空气声波衰减系数的准确计算。
以上述依据本发明的理想实施例为启示,通过上述的说明内容,相关工作人员完全可以在不偏离本项发明技术思想的范围内,进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容,必须要根据权利要求范围来确定其技术性范围。

Claims (11)

1.一种物体相对运动速度和空气声波衰减系数测量系统,包括声信号发生器,声信号采集器,信号发送装置,信号接收装置以及信号处理与计算系统,其特征是:以声信号发生器为一端,由声信号采集器和信号发送装置组成另一端,上述两端中一个为固定端,另一个为运动端,所述的信号发送装置与信号接收装置之间为有线或无线连接,信号接收装置与信号处理与计算系统之间为有线连接,运动端或固定端产生的具有平稳特征的声信号,经信号发生器发送、并被安装在固定端或运动端的声信号采集器接收后,通过信号发送装置和信号接收装置,声信号被传送到信号处理与计算系统进行处理而计算得到物体的相对运动速度和空气声波衰减系数。
2.如权利要求1所述的物体相对运动速度和空气声波衰减系数测量系统,其特征是:所述的声信号采集器的信号采样率、信号发送装置和信号接收装置的接发信号速率均不低于44kHz,且能够采集全频段的声信号。
3.一种采用如权利要求1~2中任一项所述的物体相对运动速度和空气声波衰减系数测量系统的测量方法,其特征是:具有如下步骤:
S1、所述运动端相对于固定端以一定的速度u(t)运动的同时,声信号发生器处产生一个平稳的周期声信号y0(t),这一声信号通过运动端和固定端之间的空气声波传送后,被声信号采集器接收,并以fs的采样率采集得到采样数据点数量为M的y1(t)信号;y1(t)信号经信号发送装置和信号接收装置被传送到信号处理与计算系统中;
S2、从y1(t)信号中截取一段小于1个且不小于1/2个平稳周期的特征信号y1d(t),y1d(t)信号的长度为N个数据点,其中N<M,在y1(t)信号上,沿着y1(t)信号计算y1(t)与y1d(t)的互相关性系数rxy(t),rxy(t)的总长度为M-N+1;
S3、在信号段长度不小于N的限制条件下,获取rxy(t)在各个分段上的rxy(t)最大值所对应的数据点位置ti及对应rxy(ti)数值,并以ti为参照点对y1(t)信号进行分段处理,获得对应的分段y1(ti,i)信号;
S4、读取分段y1(ti,i)信号的数据点数量、峰值数值、峰谷数值、峰值位置及峰谷位置信息;
S5、计算分段y1(ti,i)信号的绝对均值、峰值指标、峭度指标、歪度指标和裕度指标数值;
S6、在信号段长、峰值数值、峰值位置、绝对均值、峰值指标、峭度指标、歪度指标和裕度指标参数组成的多维空间考察数据点的分布情况,并删除含异常参数的数据点。
4.如权利要求3所述的物体相对运动速度和空气声波衰减系数测量方法,其特征是:在已知空气声传播速度c的前提下,利用式(1)计算空气声传播衰减系数:
Figure FDA0002764010150000021
式中,α为空气声传播衰减系数,单位为m-1;K为分段y1(ti,i)信号前后连贯结对数量;Ai、Ai-1为对应第i对前后结对分段信号的峰值振幅数值;Li、Li-1为第i对前后结对分段信号的分段数据点数量;fP为信号的采样率,单位为Hz;c是空气的声传播速度,单位为m/s,Bf为与A的取值方式相关的系数,如果A取声信号振幅,则Bf=1;如果A为信号段的绝对均值、有效值参数,则需要通过校正实验来确定Bf的数值。
5.如权利要求3所述的物体相对运动速度和空气声波衰减系数测量方法,其特征是:在未知空气声传播速度c的情况下,利用式(2)进行空气声传播速度和空气声传播衰减系数的优化搜索计算:
Figure FDA0002764010150000031
式中,β为系数,β=1/α,单位为m;r1、r2为约束函数的惩罚项系数。
6.如权利要求4或5所述的物体相对运动速度和空气声波衰减系数测量方法,其特征是:计算得到空气的声传播衰减系数α之后,运动端相对于固定端的相对运动速度ui的计算公式如下:
Figure FDA0002764010150000032
式中,Ai、Ai-1:对应第i对前后结对分段信号的峰值振幅数值;Li、Li-1:第i对前后结对分段信号的分段数据点数量,定义ui相对接近运动为负值,远离运动为正值。
7.如权利要求6所述的物体相对运动速度和空气声波衰减系数测量方法,其特征是:为了保证计算的可靠性,式(1)、(2)和式(3)中的A数值,也可以采用分段y1(ti,i)信号的绝对均值、有效值替代。
8.一种采用如权利要求1~2中任一项所述的物体相对运动速度和空气声波衰减系数测量系统的测量方法,其特征是:具有如下步骤:
(1)、所述运动端相对于固定端以一定的速度u(t)运动的同时,声信号发生器处产生一个平稳周期声信号y0(t),这一声信号通过运动端和固定端之间的空气声波传送后,被声信号采集器接收,并以fs的采样率采集得到采样数据点数量为M的y1(t)信号;y1(t)信号经信号发送装置和信号接收装置被传送到信号处理与计算系统中;
(2)、从y1(t)信号中截取一段小于1个且不小于1/2个平稳周期的特征信号y1d(t),y1d(t)信号的长度为N个数据点,其中N<M,在y1(t)信号上,沿着y1(t)信号计算y1(t)与y1d(t)的互相关性系数rxy(t),rxy(t)的总长度为M-N+1;
(3)、在信号段长度不小于N的限制条件下,获取rxy(t)在各个分段上的rxy(t)最大值所对应的数据点位置ti及对应rxy(ti)数值,以ti为参照点对y1(t)信号进行分段处理,获得对应的分段y1(ti,i)信号;
(4)、读取分段y1(ti,i)信号的数据点数量、峰值数值、峰谷数值、峰值位置和峰谷位置信息,并进行异常数据点的删除处理;
(5)、用分段y1(ti,i)信号来构造信号总长度相同的循环分段信号,并通过FFT变换计算得到循环分段信号的主频fi
9.如权利要求8所述的物体相对运动速度和空气声波衰减系数测量方法,其特征是:利用下式计算在相对运动速度为0时循环分段信号的主频f0
Figure FDA0002764010150000041
Figure FDA0002764010150000042
式中,f0为无频移或相对运动速度为0时循环分段信号的主频率,单位为Hz;Ai、Ai-1为对应第i对前后结对分段信号的峰值振幅数值;Δti为第i对前后结对分段信号的时间间隔差,单位为s;fi+1、fi为第i+1个和第i个循环分段信号的主频率,单位为Hz;在计算得到每个前后结对分段信号的f0后,取平均值作为相对运动速度为0时循环分段信号的主频率。
10.如权利要求8所述的物体相对运动速度和空气声波衰减系数测量方法,其特征是:在空气声传播速度c已知的情况下,利用式(5)来计算运动端相对于固定端的相对运动速度:
Figure FDA0002764010150000043
式中,c为空气声传播速度,单位为m/s。
11.如权利要求9或10所述的物体相对运动速度和空气声波衰减系数测量方法,其特征是:在计算得到相对运动速度ui之后,空气声传播衰减系数α的计算公式如下:
Figure FDA0002764010150000051
式中,α为空气声传播衰减系数,单位为m-1;K为分段y1(ti,i)信号前后连贯结对数量;fp为声波信号采样率,单位为Hz。
CN202011227342.9A 2020-11-06 2020-11-06 一种物体相对运动速度和空气声波衰减系数的测量系统及测量方法 Active CN112504434B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011227342.9A CN112504434B (zh) 2020-11-06 2020-11-06 一种物体相对运动速度和空气声波衰减系数的测量系统及测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011227342.9A CN112504434B (zh) 2020-11-06 2020-11-06 一种物体相对运动速度和空气声波衰减系数的测量系统及测量方法

Publications (2)

Publication Number Publication Date
CN112504434A true CN112504434A (zh) 2021-03-16
CN112504434B CN112504434B (zh) 2022-08-26

Family

ID=74955604

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011227342.9A Active CN112504434B (zh) 2020-11-06 2020-11-06 一种物体相对运动速度和空气声波衰减系数的测量系统及测量方法

Country Status (1)

Country Link
CN (1) CN112504434B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1534004A (en) * 1975-05-27 1978-11-29 Schlumberger Ltd Method and apparatus for determining acoustic wave parameters indicative of wave velocity in well logging
JP2006524798A (ja) * 2003-02-13 2006-11-02 ネクセンス リミテッド 様々なパラメータの高感度測定を行なうための装置、およびそのような装置で特に有用なセンサ
WO2007003058A1 (en) * 2005-07-06 2007-01-11 National Research Council Of Canada Method and system for determining material properties using ultrasonic attenuation
CN101482614A (zh) * 2008-01-07 2009-07-15 格库技术有限公司 声音传播速度建模方法、装置和系统
CN104990994A (zh) * 2015-05-28 2015-10-21 国家深海基地管理中心 基于脉冲压缩技术的沉积物中声速和衰减的宽带测量方法
CN108956781A (zh) * 2018-07-27 2018-12-07 国家海洋局第海洋研究所 沉积物声速和声衰减测量装置及方法
CN109934206A (zh) * 2019-04-08 2019-06-25 中国矿业大学(北京) 一种非平稳工况下的旋转机械故障诊断方法
CN110296327A (zh) * 2019-06-19 2019-10-01 常州大学 一种基于瞬变流频率响应分析的管道泄漏检测方法
CN110726775A (zh) * 2019-11-14 2020-01-24 四川省食品药品检验检测院(四川省药品质量研究所、四川省医疗器械检测中心) 声速与声衰减系数测量装置及方法
CN211603049U (zh) * 2019-11-14 2020-09-29 四川省食品药品检验检测院(四川省药品质量研究所、四川省医疗器械检测中心) 声速与声衰减系数测量装置

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1534004A (en) * 1975-05-27 1978-11-29 Schlumberger Ltd Method and apparatus for determining acoustic wave parameters indicative of wave velocity in well logging
JP2006524798A (ja) * 2003-02-13 2006-11-02 ネクセンス リミテッド 様々なパラメータの高感度測定を行なうための装置、およびそのような装置で特に有用なセンサ
WO2007003058A1 (en) * 2005-07-06 2007-01-11 National Research Council Of Canada Method and system for determining material properties using ultrasonic attenuation
CN101482614A (zh) * 2008-01-07 2009-07-15 格库技术有限公司 声音传播速度建模方法、装置和系统
CN104990994A (zh) * 2015-05-28 2015-10-21 国家深海基地管理中心 基于脉冲压缩技术的沉积物中声速和衰减的宽带测量方法
CN108956781A (zh) * 2018-07-27 2018-12-07 国家海洋局第海洋研究所 沉积物声速和声衰减测量装置及方法
CN109934206A (zh) * 2019-04-08 2019-06-25 中国矿业大学(北京) 一种非平稳工况下的旋转机械故障诊断方法
CN110296327A (zh) * 2019-06-19 2019-10-01 常州大学 一种基于瞬变流频率响应分析的管道泄漏检测方法
CN110726775A (zh) * 2019-11-14 2020-01-24 四川省食品药品检验检测院(四川省药品质量研究所、四川省医疗器械检测中心) 声速与声衰减系数测量装置及方法
CN211603049U (zh) * 2019-11-14 2020-09-29 四川省食品药品检验检测院(四川省药品质量研究所、四川省医疗器械检测中心) 声速与声衰减系数测量装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LIU, XU: "Velocity and attenuation dispersion relations for the effective Blot model: total-field formulation", 《NEAR SURFACE GEOPHYSICS》 *
ZAGZEBSKI, JA: "Simulation study of effects of speed of sound and attenuation on ultrasound lateral resolution", 《ULTRASOUND IN MEDICINE AND BIOLOGY》 *
伍枝平: "HIFU治疗中声衰减及声阻抗特性的测量与研究", 《硕士电子期刊出版工程科技Ⅱ辑》 *
吕苗荣等: "钻井管柱振动信号系统仿真研究", 《石油机械》 *

Also Published As

Publication number Publication date
CN112504434B (zh) 2022-08-26

Similar Documents

Publication Publication Date Title
CN106842128B (zh) 运动目标的声学跟踪方法及装置
Qin et al. Rolling bearing fault diagnosis with adaptive harmonic kurtosis and improved bat algorithm
CN114397475B (zh) 一种适用于超声波水表的水流流速测量方法
CN106643982B (zh) 一种基于声波共振频率的液位测量方法
CN110596247B (zh) 一种温度变化环境下的超声结构健康监测方法
CN103645248B (zh) 一种基于超声相速度的高温合金晶粒度评价方法
CN109459745A (zh) 一种利用辐射噪声估计运动声源速度的方法
CN111157065A (zh) 气体超声流量计超声波信号传输回路中声延时测量方法
CN104748704A (zh) 薄壁结构超声共振测厚频谱分析内插校正方法
CN108888918B (zh) 一种用于复杂路径下多目标运动速度测量系统及方法
CN105717198A (zh) 一种面向结构冲击定位的单频-重估计music方法
CN112504434B (zh) 一种物体相对运动速度和空气声波衰减系数的测量系统及测量方法
CN102670252B (zh) 颅内压无创测量方法及系统
CN111189930A (zh) 基于管体的裂纹检测方法、系统、设备及存储介质
CN109696660B (zh) 一种检测自由声场传声筒幅值灵敏度与相位精确测量方法
Zheng et al. A new characteristic peaks group judgement method for the accurate measurement of time‐of‐flight in the ultrasonic gas flowmeter
Zhang Flow measurement of natural gas in pipeline based on 1d-convolutional neural network
CN103399172A (zh) 一种超声波风速测量方法
CN112577558B (zh) 基于云端加边缘计算的超声波流量计量系统及边缘设备
CN109782018B (zh) 一种超声波风速风向测量实验仪及测量方法
CN102749134B (zh) 一种水声声管测量声速的方法
CN110529200B (zh) 一种嵌入式汽轮机叶片健康监测系统
Yao et al. Time-of-flight Extraction Method Based on Cross-correlation for Stress Measurement
CN110794208A (zh) 一种基于数字相位推算的宽带快速调制域测频方法
CN204143486U (zh) 一种基于碰撞声音识别的籽粒计数传感器

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