CN111442747A - 一种超声波信号处理方法 - Google Patents

一种超声波信号处理方法 Download PDF

Info

Publication number
CN111442747A
CN111442747A CN202010175175.1A CN202010175175A CN111442747A CN 111442747 A CN111442747 A CN 111442747A CN 202010175175 A CN202010175175 A CN 202010175175A CN 111442747 A CN111442747 A CN 111442747A
Authority
CN
China
Prior art keywords
curve
array
signal processing
resonance
processing method
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
CN202010175175.1A
Other languages
English (en)
Other versions
CN111442747B (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
China Nuclear Power Operation Technology Corp Ltd
Original Assignee
Research Institute of Nuclear Power Operation
China Nuclear Power Operation Technology Corp Ltd
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, China Nuclear Power Operation Technology Corp Ltd filed Critical Research Institute of Nuclear Power Operation
Priority to CN202010175175.1A priority Critical patent/CN111442747B/zh
Publication of CN111442747A publication Critical patent/CN111442747A/zh
Application granted granted Critical
Publication of CN111442747B publication Critical patent/CN111442747B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B17/00—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations
    • G01B17/02—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations for measuring thickness
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)

Abstract

本发明涉及信号处理技术领域，具体公开了一种超声波信号处理方法。该方法具体包括：1、利用超声波对被检工件进行壁厚测量，采集并存储仅包含共振波的射频信号数组；2、创建逼近离散曲线；3、以逼近离散曲线为窗函数，与采集的射频信号数组数据进行卷积计算，获得新的曲线；4、获取新曲线的周期特征点的采样位置；5、采用线性回归算法获得精确的共振波频率。本发明所述的一种超声波信号处理方法能够在保证管壁厚测量精度的同时，极大缩短壁厚信号处理时间，提高系统工作效率。

Description

一种超声波信号处理方法
技术领域
本发明属于信号处理技术领域，具体涉及一种超声波信号处理方法。
背景技术
超声波在所测量的有限长材料中传播时,如果超声波的频率与材料的固有频率相等,就会发生共振现象。利用共振现象提出的超声共振测厚法是一种常用的测厚方法，具有可测下限小、精度高等优点。共振测厚的原理是T＝c/2f，T为被测工件的厚度，c为超声波在工件中传播的声速，f为共振波的频率。因此，超声共振频率的估计精度，是采用超声共振法测厚的关键。传统的方法是直接通过傅里叶变换获得其频谱最大值的频率。
高精密薄壁部件被广泛应用于航空航天和核电领域的核心部件中，实际应用过程中对这类部件的壁厚精度要求极高，通常为微米级。为了保证其尺寸精度，超声共振法则被用于高精密部件的全范围测量。以一根4m长，0.5mm厚的薄壁管为例，为了实现全范围测量，需要采集超过10万组的射频数据进行共振频率分析。同时为了保证测量精度，在进行傅里叶计算的时候需要进行补零操作。然而过长的离散数据进行傅里叶变化时，计算耗时也会随之增加。完成10万组的计算，耗时近160s，无法满足大批量管材的高速检查的需求。
发明内容
本发明的目的在于提供一种超声波信号处理方法，能够代替现有的离散傅里叶变换实现高速高精度的共振波频率计算，提高超声波对高精密部件壁厚的全范围壁厚测量数据分析的效率。
本发明的技术方案如下：一种超声波信号处理方法，该方法具体包括如下步骤：
步骤1、利用超声波对被检工件进行壁厚测量，采集并存储仅包含共振波的射频信号数组；
步骤2、创建逼近离散曲线；
步骤3、以步骤2中逼近离散曲线为窗函数，与步骤1采集的射频信号数组数据进行卷积计算，获得新的曲线；
步骤4、获取步骤3中新曲线的周期特征点的采样位置；
利用搜索算法获取曲线ArrayResonance所有拐点的横坐标，获得数组CPi；
步骤5、采用线性回归算法获得精确的共振波频率；
步骤5.1、获得共振波的周期或者半周期；
利用线性回归算法获得回归直线的斜率，即为共振波的周期或者半周期t；
Figure BDA0002410575970000021
其中，x为CPi的序号i，y为CPi数组，
Figure BDA0002410575970000022
和
Figure BDA0002410575970000023
为各自的均值；
步骤5.2、当以周期特征点计算时，共振波的共振频率为f＝1/t；若以半周期特征点计算时，共振波的共振频率为f＝1/2t。
所述的步骤2中创建逼近离散曲线的具体步骤为：
步骤2.1、在相同工况下，采集或测量被检工件的标准声速V；
步骤2.2、计算获得共振波的估计波长；
利用公式N＝T*F/V，计算获得共振波的估计波长，即采样点数，其中，T为被检工件的厚度，F为系统采样的采样频率；
步骤2.3、创建正/余弦离散曲线；
以sa＝360°/(N-1)为间距，创建正/余弦离散曲线ArrayE(n)＝sin(sa*n+θ)，其中，n为曲线的横坐标，具体为0～N-1的正整数；ArrayE(n)为n点对应的纵坐标，θ为任意角度的相位角。
所述的步骤2中创建逼近离散曲线的具体步骤为：
以sa＝360°/(N-1)为间距，创建正/余弦离散曲线ArrayE(n)＝sin(sa*n+θ)，其中，n为曲线的横坐标，具体为0～N-1的正整数；ArrayE(n)为n点对应的纵坐标，θ为任意角度的相位角；
所述的步骤5后还进一步包括步骤6为：
若t与N的差值小于1，则不用继续计算，共振波周期为f；若t与N的差值大于1，则将t作为新的N值N＝t，重复步骤3至步骤5，获得新的t值；当新获得的t与N的差值小于1，则完成计算，获得共振波准确频率f。
所述的步骤4中利用搜索算法获取曲线ArrayResonance所有拐点的横坐标，获得数组CPi的具体步骤为：
步骤4.1、从曲线ArrayResonance的第一个点开始，获取数据点ARi和ARi+1；
步骤4.2、如果ARi*ARi+1>0，则顺延至下一个点i＝i+1；
步骤4.3、如果ARi*ARi+1<＝0，则采用插值法计算0点的准确位置；
步骤4.4、完成所有点的搜索，获得曲线ArrayResonance所有拐点的横坐标，获得数组CPi。
所述的步骤4.3中，如果ARi*ARi+1<＝0，采用线性插值，则0点的准确位置在Cp0＝i+ARi/(ARi-ARi+1)。
所述的步骤1中利用超声波对被检工件进行壁厚测量，采集并存储仅包含共振波的射频信号数组的具体步骤为：
对被检工件壁厚测量超声波信号采集、存储的共振波周期不小于3个周期。
所述的步骤2中创建逼近离散曲线的具体步骤为：
以sa＝360°/(N-1)为间距，创建正/余弦离散曲线ArrayE(n)＝sin(sa*n+θ)，其中，n为曲线的横坐标，具体为0～N-1的正整数；ArrayE(n)为n点对应的纵坐标，θ为任意角度的相位角；
所述的步骤5后还进一步包括步骤6为：
若t与N的差值小于1，则不用继续计算，共振波周期为f；若t与N的差值大于1，则将t作为新的N值N＝t，重复步骤3至步骤5，获得新的t值；当完成有限次的循环重复步骤3至步骤5，获得新的t值与N的差值大于等于1，则停止循环，并获得共振波准确频率f。
利用步骤2创建的逼近离散曲线ArrayE为窗函数，对数组ArrayFre数据进行卷积计算，获得新的曲线ArrayResonance。
本发明的显著效果在于：本发明所述的一种超声波信号处理方法，能够在保证管壁厚测量精度的同时，极大缩短壁厚信号处理时间，提高系统工作效率；有的对10万个射频信号采用傅里叶变换处理的时间约为160s，而利用本发明所述的一种超声波信号处理方法处理时间约为2.5s，极大缩短了处理时间，提高了工作效率。
具体实施方式
实施例1
一种超声波信号处理方法，该方法具体包括如下步骤：
步骤1、利用超声波对被检工件进行壁厚测量，采集并存储仅包含共振波的射频信号数组；
对被检工件壁厚测量超声波信号采集、存储的共振波周期不小于3个周期，记为数组ArrayFre，同时记录系统采样的采样频率F；
步骤2、创建逼近离散曲线
步骤2.1、在相同工况下，采集或测量被检工件的标准声速V；
步骤2.2、计算获得共振波的估计波长；
利用公式N＝T*F/V，计算获得共振波的估计波长，即采样点数，其中，T为被检工件的厚度；
步骤2.3、创建正/余弦离散曲线；
以sa＝360°/(N-1)为间距，创建正/余弦离散曲线ArrayE(n)＝sin(sa*n+θ)，其中，n为曲线的横坐标，具体为0～N-1的正整数；ArrayE(n)为n点对应的纵坐标，θ为任意角度的相位角；
步骤3、以步骤2中逼近离散曲线为窗函数，与步骤1采集的射频信号数组数据进行卷积计算，获得新的曲线；
利用步骤2创建的逼近离散曲线ArrayE为窗函数，对数组ArrayFre数据进行卷积计算，获得新的曲线ArrayResonance；
步骤4、获取步骤3中新曲线的周期特征点的采样位置；
利用搜索算法获取曲线ArrayResonance所有拐点的横坐标；
以曲线ArrayResonance坐标零点为特征点，任意两个相邻零点的时间间隔为半个周期，用搜索算法可以获取所有拐点的横坐标，记为CPi，i＝0,1,2…；
步骤4.1、从曲线ArrayResonance的第一个点开始，获取数据点ARi和ARi+1；
步骤4.2、如果ARi*ARi+1>0，则顺延至下一个点i＝i+1；
步骤4.3、如果ARi*ARi+1<＝0，则采用插值法计算0点的准确位置；
如果ARi*ARi+1<＝0，采用线性插值，则0点的准确位置在Cp0＝i+ARi/(ARi-ARi+1)；
步骤4.4、完成所有点的搜索，获得曲线ArrayResonance所有拐点的横坐标，获得数组CPi；
步骤5、采用线性回归算法获得精确的共振波频率；
步骤5.1、获得共振波的周期或者半周期；
利用线性回归算法获得回归直线的斜率，即为共振波的周期或者半周期t；
Figure BDA0002410575970000061
其中，x为CPi的序号i，y为CPi数组，
Figure BDA0002410575970000062
和
Figure BDA0002410575970000063
为各自的均值；
步骤5.2、当以周期特征点计算时，共振波的共振频率为f＝1/t；若以半周期特征点计算时，共振波的共振频率为f＝1/2t。
实施例2
一种超声波信号处理方法，该方法具体包括如下步骤：
步骤1、利用超声波对被检工件进行壁厚测量，采集并存储仅包含共振波的射频信号数组；
对被检工件壁厚测量超声波信号采集、存储的共振波周期不小于3个周期，记为数组ArrayFre，同时记录系统采样的采样频率F；
步骤2、创建逼近离散曲线
以sa＝360°/(N-1)为间距，创建正/余弦离散曲线ArrayE(n)＝sin(sa*n+θ)，其中，n为曲线的横坐标，具体为0～N-1的正整数；ArrayE(n)为n点对应的纵坐标，θ为任意角度的相位角；
步骤3、以步骤2中逼近离散曲线为窗函数，与步骤1采集的射频信号数组数据进行卷积计算，获得新的曲线；
利用步骤2创建的逼近离散曲线ArrayE为窗函数，对数组ArrayFre数据进行卷积计算，获得新的曲线ArrayResonance；
步骤4、获取步骤3中新曲线的周期特征点的采样位置；
利用搜索算法获取曲线ArrayResonance所有拐点的横坐标；
以曲线ArrayResonance坐标零点为特征点，任意两个相邻零点的时间间隔为半个周期，用搜索算法可以获取所有拐点的横坐标，记为CPi，i＝0,1,2…；
步骤4.1、从曲线ArrayResonance的第一个点开始，获取数据点ARi和ARi+1；
步骤4.2、如果ARi*ARi+1>0，则顺延至下一个点i＝i+1；
步骤4.3、如果ARi*ARi+1<＝0，则采用插值法计算0点的准确位置；
如果ARi*ARi+1<＝0，采用线性插值，则0点的准确位置在Cp0＝i+ARi/(ARi-ARi+1)；
步骤4.4、完成所有点的搜索，获得曲线ArrayResonance所有拐点的横坐标，获得数组CPi；
步骤5、采用线性回归算法获得精确的共振波频率；
步骤5.1、获得共振波的周期或者半周期；
利用线性回归算法获得回归直线的斜率，即为共振波的周期或者半周期t；
Figure BDA0002410575970000071
其中，x为CPi的序号i，y为CPi数组，
Figure BDA0002410575970000072
和
Figure BDA0002410575970000073
为各自的均值；
步骤5.2、当以周期特征点计算时，共振波的共振频率为f＝1/t；若以半周期特征点计算时，共振波的共振频率为f＝1/2t；
步骤6、若t与N的差值小于1，则不用继续计算，共振波周期为f；若t与N的差值大于1，则将t作为新的N值N＝t，重复步骤3至步骤5，获得新的t值；当新获得的t与N的差值小于1，则完成计算，或者完成有限次循环后停止计算，获得共振波准确频率f。

Claims (8)

1.一种超声波信号处理方法，其特征在于：该方法具体包括如下步骤：
步骤1、利用超声波对被检工件进行壁厚测量，采集并存储仅包含共振波的射频信号数组；
步骤2、创建逼近离散曲线；
步骤3、以步骤2中逼近离散曲线为窗函数，与步骤1采集的射频信号数组数据进行卷积计算，获得新的曲线；
步骤4、获取步骤3中新曲线的周期特征点的采样位置；
利用搜索算法获取曲线ArrayResonance所有拐点的横坐标，获得数组CPi；
步骤5、采用线性回归算法获得精确的共振波频率；
步骤5.1、获得共振波的周期或者半周期；
利用线性回归算法获得回归直线的斜率，即为共振波的周期或者半周期t；
Figure FDA0002410575960000011
其中，x为CPi的序号i，y为CPi数组，
Figure FDA0002410575960000012
和
Figure FDA0002410575960000013
为各自的均值；
步骤5.2、当以周期特征点计算时，共振波的共振频率为f＝1/t；若以半周期特征点计算时，共振波的共振频率为f＝1/2t。
2.根据权利要求1所述的一种超声波信号处理方法，其特征在于：所述的步骤2中创建逼近离散曲线的具体步骤为：
步骤2.1、在相同工况下，采集或测量被检工件的标准声速V；
步骤2.2、计算获得共振波的估计波长；
利用公式N＝T*F/V，计算获得共振波的估计波长，即采样点数，其中，T为被检工件的厚度，F为系统采样的采样频率；
步骤2.3、创建正/余弦离散曲线；
以sa＝360°/(N-1)为间距，创建正/余弦离散曲线ArrayE(n)＝sin(sa*n+θ)，其中，n为曲线的横坐标，具体为0～N-1的正整数；ArrayE(n)为n点对应的纵坐标，θ为任意角度的相位角。
3.根据权利要求1所述的一种超声波信号处理方法，其特征在于：所述的步骤2中创建逼近离散曲线的具体步骤为：
以sa＝360°/(N-1)为间距，创建正/余弦离散曲线ArrayE(n)＝sin(sa*n+θ)，其中，n为曲线的横坐标，具体为0～N-1的正整数；ArrayE(n)为n点对应的纵坐标，θ为任意角度的相位角；
所述的步骤5后还进一步包括步骤6为：
若t与N的差值小于1，则不用继续计算，共振波周期为f；若t与N的差值大于1，则将t作为新的N值N＝t，重复步骤3至步骤5，获得新的t值；当新获得的t与N的差值小于1，则完成计算，获得共振波准确频率f。
4.根据权利要求1所述的一种超声波信号处理方法，其特征在于：所述的步骤4中利用搜索算法获取曲线ArrayResonance所有拐点的横坐标，获得数组CPi的具体步骤为：
步骤4.1、从曲线ArrayResonance的第一个点开始，获取数据点ARi和ARi+1；
步骤4.2、如果ARi*ARi+1>0，则顺延至下一个点i＝i+1；
步骤4.3、如果ARi*ARi+1<＝0，则采用插值法计算0点的准确位置；
步骤4.4、完成所有点的搜索，获得曲线ArrayResonance所有拐点的横坐标，获得数组CPi。
5.根据权利要求4所述的一种超声波信号处理方法，其特征在于：所述的步骤4.3中，如果ARi*ARi+1<＝0，采用线性插值，则0点的准确位置在Cp0＝i+ARi/(ARi-ARi+1)。
6.根据权利要求1所述的一种超声波信号处理方法，其特征在于：所述的步骤1中利用超声波对被检工件进行壁厚测量，采集并存储仅包含共振波的射频信号数组的具体步骤为：
对被检工件壁厚测量超声波信号采集、存储的共振波周期不小于3个周期。
7.根据权利要求3所述的一种超声波信号处理方法，其特征在于：所述的步骤2中创建逼近离散曲线的具体步骤为：
以sa＝360°/(N-1)为间距，创建正/余弦离散曲线ArrayE(n)＝sin(sa*n+θ)，其中，n为曲线的横坐标，具体为0～N-1的正整数；ArrayE(n)为n点对应的纵坐标，θ为任意角度的相位角；
所述的步骤5后还进一步包括步骤6为：
若t与N的差值小于1，则不用继续计算，共振波周期为f；若t与N的差值大于1，则将t作为新的N值N＝t，重复步骤3至步骤5，获得新的t值；当完成有限次的循环重复步骤3至步骤5，获得新的t值与N的差值大于等于1，则停止循环，并获得共振波准确频率f。
8.根据权利要求1所述的一种超声波信号处理方法，其特征在于：所述的步骤3具体为：利用步骤2创建的逼近离散曲线ArrayE为窗函数，对数组ArrayFre数据进行卷积计算，获得新的曲线ArrayResonance。
CN202010175175.1A 2020-03-13 2020-03-13 一种超声波信号处理方法 Active CN111442747B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010175175.1A CN111442747B (zh) 2020-03-13 2020-03-13 一种超声波信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010175175.1A CN111442747B (zh) 2020-03-13 2020-03-13 一种超声波信号处理方法

Publications (2)

Publication Number Publication Date
CN111442747A true CN111442747A (zh) 2020-07-24
CN111442747B CN111442747B (zh) 2021-11-30

Family

ID=71648846

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010175175.1A Active CN111442747B (zh) 2020-03-13 2020-03-13 一种超声波信号处理方法

Country Status (1)

Country Link
CN (1) CN111442747B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111947769A (zh) * 2020-06-30 2020-11-17 核动力运行研究所 共振波频率确定方法及装置

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1673592A2 (en) * 2003-10-14 2006-06-28 Corporation Lord Magnetostrictive sensor for measuring distances
DE102005005965A1 (de) * 2005-02-09 2006-08-17 SONOTEC Dr. zur Horst-Meyer & Münch oHG Ultraschallkompaktsensoreinrichtung zur Erkennung von Kolben in hydraulischen Systemen
CN102183228A (zh) * 2011-01-05 2011-09-14 南京航空航天大学 一种油膜厚度的超声测量方法
CN103048006A (zh) * 2012-12-17 2013-04-17 西安石油大学 一种通用的超声波发射装置及方法
US20140052389A1 (en) * 2011-01-24 2014-02-20 Vega Grieshaber Kg Phase-Based Tracking
CN104034287A (zh) * 2014-06-05 2014-09-10 大连理工大学 一种弹性各向异性金属基体热障涂层厚度超声测量方法
CN104748704A (zh) * 2013-12-31 2015-07-01 中核武汉核电运行技术股份有限公司 薄壁结构超声共振测厚频谱分析内插校正方法
CN105157810A (zh) * 2015-05-12 2015-12-16 南阳理工学院 一种全自动声速测量仪及其测量方法
CN106442332A (zh) * 2016-10-19 2017-02-22 高飞 非线性光声谐振谱传感与成像方法及装置
CN106999151A (zh) * 2014-11-25 2017-08-01 皇家飞利浦有限公司 超声系统和方法
CN107064941A (zh) * 2017-04-28 2017-08-18 湖北工业大学 一种三段式脉冲激励与特征提取的超声测距方法
CN110159253A (zh) * 2019-06-20 2019-08-23 太平洋远景石油技术(北京)有限公司 一种超声波成像测井方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1673592A2 (en) * 2003-10-14 2006-06-28 Corporation Lord Magnetostrictive sensor for measuring distances
DE102005005965A1 (de) * 2005-02-09 2006-08-17 SONOTEC Dr. zur Horst-Meyer & Münch oHG Ultraschallkompaktsensoreinrichtung zur Erkennung von Kolben in hydraulischen Systemen
CN102183228A (zh) * 2011-01-05 2011-09-14 南京航空航天大学 一种油膜厚度的超声测量方法
US20140052389A1 (en) * 2011-01-24 2014-02-20 Vega Grieshaber Kg Phase-Based Tracking
CN103048006A (zh) * 2012-12-17 2013-04-17 西安石油大学 一种通用的超声波发射装置及方法
CN104748704A (zh) * 2013-12-31 2015-07-01 中核武汉核电运行技术股份有限公司 薄壁结构超声共振测厚频谱分析内插校正方法
CN104034287A (zh) * 2014-06-05 2014-09-10 大连理工大学 一种弹性各向异性金属基体热障涂层厚度超声测量方法
CN106999151A (zh) * 2014-11-25 2017-08-01 皇家飞利浦有限公司 超声系统和方法
CN105157810A (zh) * 2015-05-12 2015-12-16 南阳理工学院 一种全自动声速测量仪及其测量方法
CN106442332A (zh) * 2016-10-19 2017-02-22 高飞 非线性光声谐振谱传感与成像方法及装置
CN107064941A (zh) * 2017-04-28 2017-08-18 湖北工业大学 一种三段式脉冲激励与特征提取的超声测距方法
CN110159253A (zh) * 2019-06-20 2019-08-23 太平洋远景石油技术(北京)有限公司 一种超声波成像测井方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
PAN QI等: "《Study on quantitative eddy current array detection technology of cracks at free segment in heat transfer tubes of nuclear power plants》", 《45TH ANNUAL REVIEW OF PROGRESS IN QUANTITATIVE NONDESTRUCTIVE EVALUATION》 *
ZHANG YICHENG等: "《Model based reliability analysis of PA ultrasonic testing for weld of hydro turbine runner》", 《PROCEDIA ENGINEERING》 *
张俊，等: "《基于Ｂｏｏｔｓｔｒａｐ方法的自动超声检测可靠性分析》", 《武汉大学学报（工学版）》 *
张益成，等: "《相控阵超声检测计算模型的建立》", 《武汉大学学报(工学版)》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111947769A (zh) * 2020-06-30 2020-11-17 核动力运行研究所 共振波频率确定方法及装置

Also Published As

Publication number Publication date
CN111442747B (zh) 2021-11-30

Similar Documents

Publication Publication Date Title
CN101825660B (zh) 欠采样下的正弦信号频率的高效测量方法及实施装置
CN108921014B (zh) 一种基于改进噪声包络信号识别的螺旋桨轴频搜索方法
CN108444653A (zh) 一种面阵型超声波气体泄漏点定位检测装置及方法
CN107328993B (zh) 频率测量装置和频率测量方法
CN105841645B (zh) 一种基于电磁超声的测厚方法
CN111442747B (zh) 一种超声波信号处理方法
CN103236825A (zh) 一种用于高精度数据采集系统的数据校正方法
CN103399204A (zh) 一种基于Rife-Vincent（II）窗插值FFT的谐波与间谐波检测方法
CN106078508B (zh) 一种砂轮磨损在线监测方法
CN102269803B (zh) 基于时间延迟的离散频谱低频成分的校正方法
CN113343173A (zh) 一种布里渊频移提取方法
CN109117816A (zh) 基于六阶样条插值小波的信号奇异点检测方法
Yang et al. Study on ways to restrain end effect of Hilbert-Huang transform
CN104748704A (zh) 薄壁结构超声共振测厚频谱分析内插校正方法
CN112394223B (zh) 一种信号分量频率和初相位的联合估计方法
CN111964586A (zh) 一种基于随机噪声修正的白光干涉信号处理方法
CN107171748B (zh) 欠采样的多阵列协同测频与直接定位方法
CN114964401A (zh) 一种基于fft的科氏流量计信号处理方法
CN110297199B (zh) 一种基于全相位fft的铯光泵磁力仪频率测量方法及系统
CN114061451A (zh) 超精密位置探测光电信号数据拟合方法及其装置
CN111486804B (zh) 用于精密部件厚度测量的信号处理方法及测量方法
Dan et al. Comparative analysis of frequency estimation methods
Yue et al. Modified algorithm of sinusoid signal frequency estimation based on Quinn and Aboutanios iterative algorithms
CN117169590B (zh) 一种基于软件变采样率的电力谐波分析的方法和装置
CN114142853B (zh) 一种基于插值dft信号同步的数字锁相放大处理方法

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