CN111947769B - 共振波频率确定方法及装置 - Google Patents

共振波频率确定方法及装置 Download PDF

Info

Publication number
CN111947769B
CN111947769B CN202010622344.1A CN202010622344A CN111947769B CN 111947769 B CN111947769 B CN 111947769B CN 202010622344 A CN202010622344 A CN 202010622344A CN 111947769 B CN111947769 B CN 111947769B
Authority
CN
China
Prior art keywords
detected
frequency
discrete function
determining
wave
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
CN202010622344.1A
Other languages
English (en)
Other versions
CN111947769A (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.)
Research Institute of Nuclear Power Operation
Original Assignee
Research Institute of Nuclear Power Operation
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 Research Institute of Nuclear Power Operation filed Critical Research Institute of Nuclear Power Operation
Priority to CN202010622344.1A priority Critical patent/CN111947769B/zh
Publication of CN111947769A publication Critical patent/CN111947769A/zh
Application granted granted Critical
Publication of CN111947769B publication Critical patent/CN111947769B/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
    • G01H13/00Measuring resonant frequency

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Complex Calculations (AREA)

Abstract

本发明属于核电维修技术领域,具体涉及一种共振波频率确定方法及装置。本公开采用估算的共振波波长确定逼近离散函数,并采用该逼近离散函数对采样数据集合进行卷积,由此大大减小了需要计算的数据量,减小了计算耗时,由于逼近离散函数是根据估算的共振波波长确定的,使得采用逼近离散函数卷积处理采样数据集合得到的目标离散函数,能够较准确的拟合待测共振波的变化趋势,进而能够较准确的得到共振波的频率。

Description

共振波频率确定方法及装置
技术领域
本发明属于核电维修技术领域,具体涉及一种共振波频率确定方法及装置。
背景技术
超声波在所测量的有限长材料中传播时,如果超声波的频率与材料的固有频率相等,就会发生共振现象。利用共振现象提出的超声共振测厚法是一种常用的测厚方法,具有可测下限小、精度高等优点。共振测厚的原理是T=c/2f,T为被测工件的厚度,c为超声波在工件中传播的声速,f为共振波的频率。因此,超声共振频率的估计精度,是采用超声共振法测厚的关键。传统的方法是直接通过傅里叶变换获得其频谱最大值的频率。
高精密薄壁部件被广泛应用于航空航天和核电领域的核心部件中,实际应用过程中对这类部件的壁厚精度要求极高,通常为微米级。为了保证其尺寸精度,超声共振法则被用于高精密部件的全范围测量。以一根4m长,0.5mm厚的薄壁管为例,为了实现全范围测量,需要采集超过10万组的射频数据进行共振频率分析。同时为了保证测量精度,在进行傅里叶计算的时候需要进行补零操作。然而过长的离散数据进行傅里叶变化时,计算耗时也会随之增加。完成10万组的计算,耗时近160s,无法满足大批量管材的高速检查的需求。
发明内容
为克服相关技术中存在的问题,提供了一种共振波频率确定方法及装置。
根据本公开实施例的一方面,提供一种共振波频率确定方法,所述方法包括:
步骤100,获取待测共振波的采样数据集合,所述采样数据集合包括多个采样数据,各相邻采样数据的采集间隔时长相同;
步骤101,根据检测被检工件的超声波的频率、所述超声波的标准传播速度以及预估的被检工件的厚度,确定待测共振波的预估波长值;
步骤102,根据所述预估波长值,确定逼近离散函数;
步骤103,采用所述逼近离散函数作为窗口函数,对所述采样数据集合进行卷积处理,得到目标离散函数;
步骤104,根据所述目标离散函数,确定拐点数据集合,所述拐点数据集合包括所述目标离散函数各拐点对应的横坐标;
步骤105,根据所述拐点数据集合以及线性回归模型,确定待测共振波的频率。
在一种可能的实现方式中,所述步骤102包括:
将式1确定为逼近离散函数:
Figure BDA0002563450940000021
其中,A(n)为逼近离散函数,N为预估波长值,θ为任意角度的相位角,n∈(0,N-1)的正整数。
在一种可能的实现方式中,所述步骤104,包括:
获取目标离散函数的各横坐标Xi对应的纵坐标Yi,Xi为正整数,i为正整数;
判断Yi*Yi+1与0的大小关系;
在Yi*Yi+1≤0的情况下,根据Yi、Yi+1,采用插值法确定目标离散函数的拐点对应的横坐标。
在一种可能的实现方式中,所述方法还包括:
在确定的待测共振波的频率的倒数与所述预估波长值之差小于或等于预设阈值的情况下,将确定的待测共振波的频率作为最终的计算结果;
在确定的待测共振波的频率的倒数与所述预估波长值之差大于预设阈值的情况下,将确定的待测共振波的频率的倒数作为预估波长值,重复所述步骤102至步骤105的操作,直至所述预估波长值之差小于或等于预设阈值。
根据本公开实施例的另一方面,提供一种共振波频率确定装置,所述装置包括:
获取模块,用于获取待测共振波的采样数据集合,所述采样数据集合包括多个采样数据,各相邻采样数据的采集间隔时长相同;
第一确定模块,用于根据检测被检工件的超声波的频率、所述超声波的标准传播速度以及预估的被检工件的厚度,确定待测共振波的预估波长值;
第二确定模块,用于根据所述预估波长值,确定逼近离散函数;
卷积模块,用于采用所述逼近离散函数作为窗口函数,对所述采样数据集合进行卷积处理,得到目标离散函数;
第三确定模块,用于根据所述目标离散函数,确定拐点数据集合,所述拐点数据集合包括所述目标离散函数各拐点的横坐标;
第四确定模块,用于根据所述拐点数据集合以及线性回归模型,确定待测共振波的频率。
在一种可能的实现方式中,所述第二确定模块包括:
第一确定子模块,用于将式1确定为逼近离散函数:
Figure BDA0002563450940000031
其中,A(n)为逼近离散函数,N为预估波长值,θ为任意角度的相位角,n∈(0,N-1)的正整数。
在一种可能的实现方式中,所述第三确定模块包括:
获取子模块,用于获取目标离散函数的各横坐标Xi对应的纵坐标Yi,Xi为正整数,i为正整数;
判断子模块,用于判断Yi*Yi+1与0的大小关系;
确定子模块,用于在Yi*Yi+1≤0的情况下,根据Yi、Yi+1,采用插值法确定目标离散函数的拐点对应的横坐标。
在一种可能的实现方式中,所述装置还包括:
第一验证模块,用于在确定的待测共振波的频率的倒数与所述预估波长值之差小于或等于预设阈值的情况下,将确定的待测共振波的频率作为最终的计算结果;
第二验证模块,用于在确定的待测共振波的频率的倒数与所述预估波长值之差大于预设阈值的情况下,将确定的待测共振波的频率的倒数作为预估波长值,重复所述第二确定模块、卷积模块、第三确定模块、第四确定模块的操作,直至所述预估波长值之差小于或等于预设阈值。
根据本公开实施例的另一方面,提供一种无线控制装置,包括:
处理器;
用于存储处理器可执行指令的存储器;
其中,所述处理器被配置为:
执行上述方法。
根据本公开实施例的另一方面,提供一种非临时性计算机可读存储介质,当所述存储介质中的指令由处理器执行时,使得处理器能够执行上述方法。
本发明的有益效果在于:本公开采用估算的共振波波长确定逼近离散函数,并采用该逼近离散函数对采样数据集合进行卷积,由此大大减小了需要计算的数据量,减小了计算耗时,由于逼近离散函数是根据估算的共振波波长确定的,使得采用逼近离散函数卷积处理采样数据集合得到的目标离散函数,能够较准确的拟合待测共振波的变化趋势,进而能够较准确的得到共振波的频率。
附图说明
图1是根据一示例性实施例示出的一种共振波频率确定方法的流程图。
图2是根据一示例性实施例示出的待测共振波信号的示意图。
图3是根据一示例性实施例示出的逼近离散函数的示意图。
图4是根据一示例性实施例示出的目标离散函数的示意图。
图5是根据一示例性实施例示出的一种共振波频率确定装置的框图。
图6是根据一示例性实施例示出的一种共振波频率确定装置的框图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
图1是根据一示例性实施例示出的一种共振波频率确定方法的流程图。该方法可以应用于台式电脑或服务器等终端设备,本公开实施例对终端设备的类型不做限定。如图1所示,该方法可以包括:
步骤100,获取待测共振波的采样数据集合,所述采样数据集合包括多个采样数据,相邻采样数据的采集间隔时长相同。
举例来讲,可以采用预设频率的超声波对被检工件进行壁厚测量,采集并存储仅包含共振波的射频信号的采样数据形成采样数据集合,图2是根据一示例性实施例示出的待测共振波信号的示意图,如图2所示,该共振波的射频信号所包含的共振波周期可以大于3个周期共振波周期。
步骤101,根据检测被检工件的超声波的频率、所述超声波的标准传播速度以及预估的被检工件的厚度,确定待测共振波的预估波长值。
举例来讲,可以根据N=T*F/V确定待测共振波的预设波长值N,
其中,T为预估的被检工件的厚度,F为检测被检工件的超声波的频率,V为检测被检工件的超声波在相同工况下的标准传播速度(该标准传播速度可以例如为经验值,或历史平均值)。
步骤102,根据所述预估波长值,确定逼近离散函数;
举例来讲,所述步骤102可以包括:
将式1确定为逼近离散函数:
Figure BDA0002563450940000061
其中,A(n)为逼近离散函数,N为预估波长值,θ为任意角度的相位角,n∈(0,N-1),且n为正整数。可以例如将预设波长值N处理为正整数,以便于后续计算。图3是根据一示例性实施例示出的逼近离散函数的示意图。如图3所示,该逼近离散函数为正弦离散函数。
在另一个示例中,该逼近离散函数也可以为余弦离散函数,例如,
Figure BDA0002563450940000062
步骤103,采用所述逼近离散函数作为窗口函数,对所述采样数据集合进行卷积处理,得到目标离散函数;
举例来讲,该采样数据集合可以为{b1、b2、……bM},其中,M为正整数,且M>N,上述式1中逼近离散函数的纵坐标集合可以为{a1、a2、……aN},可以将窗口函数的移动步长设置为1,对采样数据集合进行卷积处理,得到目标离散函数,其中,目标离散函数的离散点集合可以为{(X1,Y1),(X2,Y2),……(XM-N+1,YM-N+1)},其中,{X1,X2,……XM-N+1}为自然数序列,
Yi=a1×bi+a2×bi+1+……+aN×bi+N,i为正整数,i∈(1,M-N+1)。
步骤104,根据所述目标离散函数,确定拐点数据集合,所述拐点数据集合包括所述目标离散函数各拐点对应的横坐标;图4是根据一示例性实施例示出的目标离散函数的示意图,如图4所示,拐点可以为各数值为0的纵坐标对应的横坐标。
接上例,可以获取目标离散函数的各横坐标Xi对应的纵坐标Yi,Xi为正整数,i为正整数;
可以判断Yi*Yi+1与0的大小关系;并在Yi*Yi+1≤0的情况下,根据Yi、Yi+1,采用插值法确定目标离散函数的拐点对应的横坐标。以线性插值法为例,可以确定拐点对应的横坐标Ck=i+Yi/Yi-Yi+1。由此可以获取拐点数据集合{C1,C2,……Cp},k∈(1,p),k为正整数,p为正整数。
步骤105,根据所述拐点数据集合以及线性回归模型,确定待测共振波的频率。接上例,可以根据式2确定共振波的半周期t。
Figure BDA0002563450940000071
其中,
Figure BDA0002563450940000072
为自然数组{1、2、3……p}的均值,/>
Figure BDA0002563450940000073
为数据集合{C1,C2,……Cp}的均值。
最后可以确定待测共振波的频率f=1/2t。
在一种可能的实现方式中,所述方法还可以包括:
在确定的待测共振波的频率的倒数与所述预估波长值之差小于或等于预设阈值(例如1)的情况下,将确定的待测共振波的频率作为最终的计算结果;
在确定的待测共振波的频率的倒数与所述预估波长值之差大于预设阈值的情况下,将确定的待测共振波的频率的倒数作为预估波长值,重复所述步骤102至步骤105的操作,直至所述预估波长值之差小于或等于预设阈值。
这样,既可以减小了计算耗时,且计算精度会随着数据长度增加而进一步增加。
在一种应用示例中,相关技术针对10万个射频信号的处理时间约为160秒,采用本公开的方法针对同样10万个射频信号的处理时长约为2.5秒,可见,本公开实施例的方法大大缩短了计算时间。
在一种可能的实现方式中,提供一种共振波频率确定装置,所述装置包括:
获取模块,用于获取待测共振波的采样数据集合,所述采样数据集合包括多个采样数据,相邻采样数据的采集间隔时长相同;
第一确定模块,用于根据检测被检工件的超声波的频率、所述超声波的标准传播速度以及预估的被检工件的厚度,确定待测共振波的预估波长值;
第二确定模块,用于根据所述预估波长值,确定逼近离散函数;
卷积模块,用于采用所述逼近离散函数作为窗口函数,对所述采样数据集合进行卷积处理,得到目标离散函数;
第三确定模块,用于根据所述目标离散函数,确定拐点数据集合,所述拐点数据集合包括所述目标离散函数各拐点的横坐标;
第四确定模块,用于根据所述拐点数据集合以及线性回归模型,确定待测共振波的频率。
在一种可能的实现方式中,所述第二确定模块包括:
第一确定子模块,用于将式1确定为逼近离散函数:
Figure BDA0002563450940000081
其中,A(n)为逼近离散函数,N为预估波长值,θ为任意角度的相位角,n∈(0,N-1)的正整数。
在一种可能的实现方式中,所述第三确定模块包括:
获取子模块,用于获取目标离散函数的各横坐标Xi对应的纵坐标Yi,Xi为正整数,i为正整数;
判断子模块,用于判断Yi*Yi+1与0的大小关系;
确定子模块,用于在Yi*Yi+1≤0的情况下,根据Yi、Yi+1,采用插值法确定目标离散函数的拐点对应的横坐标。
在一种可能的实现方式中,所述装置还包括:
第一验证模块,用于在确定的待测共振波的频率的倒数与所述预估波长值之差小于或等于预设阈值的情况下,将确定的待测共振波的频率作为最终的计算结果;
第二验证模块,用于在确定的待测共振波的频率的倒数与所述预估波长值之差大于预设阈值的情况下,将确定的待测共振波的频率的倒数作为预估波长值,重复所述第二确定模块、卷积模块、第三确定模块、第四确定模块的操作,直至所述预估波长值之差小于或等于预设阈值。
图5是根据一示例性实施例示出的一种共振波频率确定装置的框图。例如,装置800可以是移动电话,计算机,数字广播终端,消息收发设备,游戏控制台,平板设备,医疗设备,健身设备,个人数字助理等。
参照图5,装置800可以包括以下一个或多个组件:处理组件802,存储器804,电源组件806,多媒体组件808,音频组件810,输入/输出(I/O)的接口812,传感器组件814,以及通信组件816。
处理组件802通常控制装置800的整体操作,诸如与显示,电话呼叫,数据通信,相机操作和记录操作相关联的操作。处理组件802可以包括一个或多个处理器820来执行指令,以完成上述的方法的全部或部分步骤。此外,处理组件802可以包括一个或多个模块,便于处理组件802和其他组件之间的交互。例如,处理组件802可以包括多媒体模块,以方便多媒体组件808和处理组件802之间的交互。
存储器804被配置为存储各种类型的数据以支持在装置800的操作。这些数据的示例包括用于在装置800上操作的任何应用程序或方法的指令,联系人数据,电话簿数据,消息,图片,视频等。存储器804可以由任何类型的易失性或非易失性存储设备或者它们的组合实现,如静态随机存取存储器(SRAM),电可擦除可编程只读存储器(EEPROM),可擦除可编程只读存储器(EPROM),可编程只读存储器(PROM),只读存储器(ROM),磁存储器,快闪存储器,磁盘或光盘。
电源组件806为装置800的各种组件提供电力。电源组件806可以包括电源管理系统,一个或多个电源,及其他与为装置800生成、管理和分配电力相关联的组件。
多媒体组件808包括在所述装置800和用户之间的提供一个输出接口的屏幕。在一些实施例中,屏幕可以包括液晶显示器(LCD)和触摸面板(TP)。如果屏幕包括触摸面板,屏幕可以被实现为触摸屏,以接收来自用户的输入信号。触摸面板包括一个或多个触摸传感器以感测触摸、滑动和触摸面板上的手势。所述触摸传感器可以不仅感测触摸或滑动动作的边界,而且还检测与所述触摸或滑动操作相关的持续时间和压力。在一些实施例中,多媒体组件808包括一个前置摄像头和/或后置摄像头。当装置800处于操作模式,如拍摄模式或视频模式时,前置摄像头和/或后置摄像头可以接收外部的多媒体数据。每个前置摄像头和后置摄像头可以是一个固定的光学透镜系统或具有焦距和光学变焦能力。
音频组件810被配置为输出和/或输入音频信号。例如,音频组件810包括一个麦克风(MIC),当装置800处于操作模式,如呼叫模式、记录模式和语音识别模式时,麦克风被配置为接收外部音频信号。所接收的音频信号可以被进一步存储在存储器804或经由通信组件816发送。在一些实施例中,音频组件810还包括一个扬声器,用于输出音频信号。
I/O接口812为处理组件802和外围接口模块之间提供接口,上述外围接口模块可以是键盘,点击轮,按钮等。这些按钮可包括但不限于:主页按钮、音量按钮、启动按钮和锁定按钮。
传感器组件814包括一个或多个传感器,用于为装置800提供各个方面的状态评估。例如,传感器组件814可以检测到装置800的打开/关闭状态,组件的相对定位,例如所述组件为装置800的显示器和小键盘,传感器组件814还可以检测装置800或装置800一个组件的位置改变,用户与装置800接触的存在或不存在,装置800方位或加速/减速和装置800的温度变化。传感器组件814可以包括接近传感器,被配置用来在没有任何的物理接触时检测附近物体的存在。传感器组件814还可以包括光传感器,如CMOS或CCD图像传感器,用于在成像应用中使用。在一些实施例中,该传感器组件814还可以包括加速度传感器,陀螺仪传感器,磁传感器,压力传感器或温度传感器。
通信组件816被配置为便于装置800和其他设备之间有线或无线方式的通信。装置800可以接入基于通信标准的无线网络,如WiFi,2G或3G,或它们的组合。在一个示例性实施例中,通信组件816经由广播信道接收来自外部广播管理系统的广播信号或广播相关信息。在一个示例性实施例中,所述通信组件816还包括近场通信(NFC)模块,以促进短程通信。例如,在NFC模块可基于射频识别(RFID)技术,红外数据协会(IrDA)技术,超宽带(UWB)技术,蓝牙(BT)技术和其他技术来实现。
在示例性实施例中,装置800可以被一个或多个应用专用集成电路(ASIC)、数字信号处理器(DSP)、数字信号处理设备(DSPD)、可编程逻辑器件(PLD)、现场可编程门阵列(FPGA)、控制器、微控制器、微处理器或其他电子元件实现,用于执行上述方法。
在示例性实施例中,还提供了一种包括指令的非临时性计算机可读存储介质,例如包括指令的存储器804,上述指令可由装置800的处理器820执行以完成上述方法。例如,所述非临时性计算机可读存储介质可以是ROM、随机存取存储器(RAM)、CD-ROM、磁带、软盘和光数据存储设备等。
图6是根据一示例性实施例示出的一种共振波频率确定装置的框图。例如,装置1900可以被提供为一服务器。参照图6,装置1900包括处理组件1922,其进一步包括一个或多个处理器,以及由存储器1932所代表的存储器资源,用于存储可由处理组件1922的执行的指令,例如应用程序。存储器1932中存储的应用程序可以包括一个或一个以上的每一个对应于一组指令的模块。此外,处理组件1922被配置为执行指令,以执行上述方法。
装置1900还可以包括一个电源组件1926被配置为执行装置1900的电源管理,一个有线或无线网络接口1950被配置为将装置1900连接到网络,和一个输入输出(I/O)接口1958。装置1900可以操作基于存储在存储器1932的操作系统,例如Windows ServerTM,MacOS XTM,UnixTM,LinuxTM,FreeBSDTM或类似。
在示例性实施例中,还提供了一种非易失性计算机可读存储介质,例如包括计算机程序指令的存储器1932,上述计算机程序指令可由装置1900的处理组件1922执行以完成上述方法。
本公开可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质,其上载有用于使处理器实现本公开的各个方面的计算机可读程序指令。
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身,诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如,通过光纤电缆的光脉冲)、或者通过电线传输的电信号。
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备,或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令,并转发该计算机可读程序指令,以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码,所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等,以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中,通过利用计算机可读程序指令的状态信息来个性化定制电子电路,例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA),该电子电路可以执行计算机可读程序指令,从而实现本公开的各个方面。
这里参照根据本公开实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本公开的各个方面。应当理解,流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合,都可以由计算机可读程序指令实现。
这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器,从而生产出一种机器,使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时,产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中,这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作,从而,存储有指令的计算机可读介质则包括一个制造品,其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。
也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上,使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤,以产生计算机实现的过程,从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。
附图中的流程图和框图显示了根据本公开的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分,所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

Claims (4)

1.一种共振波频率确定装置,其特征在于,所述装置包括:
获取模块,用于获取待测共振波的采样数据集合,所述采样数据集合包括多个采样数据,各相邻采样数据的采集间隔时长相同;
获取模块包括:采用预设频率的超声波对被检工件进行壁厚测量,采集并存储包含共振波的射频信号的采样数据形成采样数据集合;
第一确定模块,用于根据检测被检工件的超声波的频率、所述超声波的标准传播速度以及预估的被检工件的厚度,确定待测共振波的预估波长值;
第一确定模块包括:根据N=T*F/V确定待测共振波的预设波长值N,T为预估的被检工件的厚度,F为检测被检工件的超声波的频率,V为检测被检工件的超声波在相同工况下的标准传播速度;
第二确定模块,用于根据所述预估波长值,确定逼近离散函数;
卷积模块,用于采用所述逼近离散函数作为窗口函数,对所述采样数据集合进行卷积处理,得到目标离散函数;
第三确定模块,用于根据所述目标离散函数,确定拐点数据集合,所述拐点数据集合包括所述目标离散函数各拐点的横坐标;
第四确定模块,用于根据所述拐点数据集合以及线性回归模型,确定待测共振波的频率。
2.根据权利要求1所述的装置,其特征在于,所述第二确定模块包括:
第一确定子模块,用于将式(1)确定为逼近离散函数:
Figure FDA0003951198190000011
其中,A(n)为逼近离散函数,N为预估波长值,θ为任意角度的相位角,n∈(0,N-1)的正整数。
3.根据权利要求1所述的装置,其特征在于,所述第三确定模块包括:
获取子模块,用于获取目标离散函数的各横坐标Xi对应的纵坐标Yi,Xi为正整数,i为正整数;
判断子模块,用于判断Yi*Yi+1与0的大小关系;
确定子模块,用于在Yi*Yi+1≤0的情况下,根据Yi、Yi+1,采用插值法确定目标离散函数的拐点对应的横坐标。
4.根据权利要求1所述的装置,其特征在于,所述装置还包括:
第一验证模块,用于在确定的待测共振波的频率的倒数与所述预估波长值之差小于或等于预设阈值的情况下,将确定的待测共振波的频率作为最终的计算结果;
第二验证模块,用于在确定的待测共振波的频率的倒数与所述预估波长值之差大于预设阈值的情况下,将确定的待测共振波的频率的倒数作为预估波长值,重复所述第二确定模块、卷积模块、第三确定模块、第四确定模块的操作,直至所述预估波长值之差小于或等于预设阈值。
CN202010622344.1A 2020-06-30 2020-06-30 共振波频率确定方法及装置 Active CN111947769B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010622344.1A CN111947769B (zh) 2020-06-30 2020-06-30 共振波频率确定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010622344.1A CN111947769B (zh) 2020-06-30 2020-06-30 共振波频率确定方法及装置

Publications (2)

Publication Number Publication Date
CN111947769A CN111947769A (zh) 2020-11-17
CN111947769B true CN111947769B (zh) 2023-03-28

Family

ID=73337824

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010622344.1A Active CN111947769B (zh) 2020-06-30 2020-06-30 共振波频率确定方法及装置

Country Status (1)

Country Link
CN (1) CN111947769B (zh)

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4305294A (en) * 1980-01-10 1981-12-15 Rockwell International Corporation Ultrasonic apparatus and method for measuring wall thickness
US7600442B2 (en) * 2004-11-16 2009-10-13 H&B System Co., Ltd. Ultrasonic probing method and apparatus therefor utilizing resonance phenomenon
CN103148815B (zh) * 2013-01-30 2015-12-23 大连理工大学 基于声压反射系数自相关函数的薄层厚度超声检测方法
CN104748704A (zh) * 2013-12-31 2015-07-01 中核武汉核电运行技术股份有限公司 薄壁结构超声共振测厚频谱分析内插校正方法
CN104034287B (zh) * 2014-06-05 2016-06-22 大连理工大学 一种弹性各向异性金属基体热障涂层厚度超声测量方法
CN106371098B (zh) * 2016-08-30 2018-12-14 浙江工商大学 一种有阻尼自由振动超声波测距系统频率拐点判别方法
CN109684650B (zh) * 2017-10-18 2021-08-27 大唐移动通信设备有限公司 一种确定中心频率的方法和装置
CN108225545A (zh) * 2017-12-13 2018-06-29 珠海格力节能环保制冷技术研究中心有限公司 机械谐振频率测量方法、装置、存储介质及测量仪
CN108004388B (zh) * 2017-12-22 2019-05-03 上海海事大学 一种确定超声振动时效激振频率的方法
CN111442747B (zh) * 2020-03-13 2021-11-30 中核武汉核电运行技术股份有限公司 一种超声波信号处理方法

Also Published As

Publication number Publication date
CN111947769A (zh) 2020-11-17

Similar Documents

Publication Publication Date Title
TWI781359B (zh) 人臉和人手關聯檢測方法及裝置、電子設備和電腦可讀儲存媒體
CN111551191B (zh) 传感器外参数标定方法及装置、电子设备和存储介质
CN107692997B (zh) 心率检测方法及装置
KR20200131305A (ko) 키포인트 검출 방법, 장치, 전자 기기 및 기억 매체
WO2018157630A1 (zh) 关联用户推荐方法及装置
CN111539410B (zh) 字符识别方法及装置、电子设备和存储介质
CN111612070A (zh) 基于场景图的图像描述生成方法及装置
CN107508573B (zh) 晶振振荡频率校正方法及装置
CN111947903B (zh) 振动异常定位方法及装置
CN111640114B (zh) 图像处理方法及装置
CN111931844A (zh) 图像处理方法及装置、电子设备和存储介质
WO2021103994A1 (zh) 用于信息推荐的模型训练方法、装置、电子设备以及介质
CN112082656A (zh) 温度确定方法及装置、电子设备、存储介质
CN108984628B (zh) 内容描述生成模型的损失值获取方法及装置
CN111947769B (zh) 共振波频率确定方法及装置
CN115512116B (zh) 图像分割模型优化方法、装置、电子设备及可读存储介质
CN110121115B (zh) 精彩视频片段的确定方法及装置
CN114547073B (zh) 时序数据的聚合查询方法、装置及存储介质
CN116310633A (zh) 一种关键点检测模型训练方法及关键点检测方法
WO2020164261A1 (zh) 液质检测结果的质量控制方法及装置
CN112802032A (zh) 图像分割网络的训练和图像处理方法、装置、设备及介质
CN113850777B (zh) 数据处理装置、方法、电子设备及存储介质
CN111949924B (zh) 渗透率的确定方法及装置、电子设备和存储介质
CN116626766B (zh) 水质监测设备状态检测方法、装置、电子设备和存储介质
CN109388782B (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