CN116089800A - 激波流场动态压力测量信号振铃分量提取校正方法及系统 - Google Patents

激波流场动态压力测量信号振铃分量提取校正方法及系统 Download PDF

Info

Publication number
CN116089800A
CN116089800A CN202310368863.3A CN202310368863A CN116089800A CN 116089800 A CN116089800 A CN 116089800A CN 202310368863 A CN202310368863 A CN 202310368863A CN 116089800 A CN116089800 A CN 116089800A
Authority
CN
China
Prior art keywords
dynamic pressure
pressure measurement
ringing
data
flow field
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
CN202310368863.3A
Other languages
English (en)
Other versions
CN116089800B (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.)
Wuhan Institute of Technology
Original Assignee
Wuhan Institute of Technology
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 Wuhan Institute of Technology filed Critical Wuhan Institute of Technology
Priority to CN202310368863.3A priority Critical patent/CN116089800B/zh
Publication of CN116089800A publication Critical patent/CN116089800A/zh
Application granted granted Critical
Publication of CN116089800B publication Critical patent/CN116089800B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L23/00—Devices or apparatus for measuring or indicating or recording rapid changes, such as oscillations, in the pressure of steam, gas, or liquid; Indicators for determining work or energy of steam, internal-combustion, or other fluid-pressure engines from the condition of the working fluid
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M10/00—Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

本发明公开了一种激波流场动态压力测量信号振铃分量提取校正方法及系统，该方法包括以下步骤：采用经验模态分解方法对激波流场动态压力测量信号进行自适应分解，提取出激波流场动态压力测量信号的振铃分量；利用卷积神经网络来构建振铃分量畸变校正模型，数据集构建采用交错跨度式构建方法；对振铃分量中模态混叠段时间序列进行预测，实现对动态压力响应信号振铃分量的准确提取。本发明将自适应信号分解方法和卷积神经网络模型辨识方法的优点相结合，可有效减小激波流场动态压力响应信号振铃分量提取过程受频率混叠的影响，通过振铃分量的变化特征构建数据预测模型，实现激波流场动态压力响应信号振铃分量的准确提取与校正。

Description

激波流场动态压力测量信号振铃分量提取校正方法及系统
技术领域
本发明属于信号处理技术领域，具体涉及一种激波流场动态压力测量信号振铃分量提取校正方法及系统。
背景技术
激波管动态压力广泛存在于爆炸测试、医疗器械、材料冲击测试、航空发动机等领域。激波管动态压力的产生是一种瞬变的动态过程，不仅测试环境复杂，而且很难控制，从而导致动态压力信号难以准确估计，严重地影响动态压力信号的测量精度。
压力传感器作为测量动态压力的元器件，对其得到的动态压力测量信号进行分析和处理至关重要。其中，对响应信号中振铃分量的提取与剔除程度与动态压力信号的重构精度息息相关。由于响应信号中包含的噪声成分较为复杂，传统的信号提取方法得到的振铃分量会出现模态混叠现象，从而导致部分振铃分量信号的缺失。
发明内容
本发明的目的在于，提供一种激波流场动态压力测量信号振铃分量提取校正方法及系统，解决传统的信号提取方法得到的振铃分量会出现模态混叠现象的问题。
本发明所采用的技术方案如下：
提供一种激波流场动态压力测量信号振铃分量提取校正方法，包括以下步骤：
S1、采用经验模态分解方法对激波流场动态压力测量信号进行自适应分解，提取出激波流场动态压力测量信号的振铃分量；
S2、从所提取出的振铃分量中选取未发生模态混叠的连续时间序列的序列数据，将未发生模态混叠的连续时间序列的序列数据分解组合为 K个等长的连续时间子序列，；其中，每个子序列的序列数据量为 d+1个，每个子序列如下所示：
式中，为未发生模态混叠的连续时间序列的第 k个序列数据；
S3、利用卷积神经网络来构建振铃分量畸变校正模型，并将子序列的 d+1个序列数据作为该振铃分量畸变校正模型的输入特征，将未发生模态混叠的连续时间序列的第 k+ d+zt个序列数据作为模型的输出特征，训练模型；其中， zt为时间跨度；
S4、采集待校正的激波流场动态压力测量信号，并进行经验模态分解，找到振铃分量中发生模态混叠的第一个数据点作为畸变点，以畸变点前第 d+zt个数据点作为起始点，构建输入特征，并输入至训练好的振铃分量畸变校正模型，得到畸变校正数据，之后利用畸变校正数据替换掉畸变点；重复该步骤，获取最终的激波流场动态压力响应信号振铃分量。
接上述技术方案，步骤S1具体包括以下步骤：
S11、识别激波流场动态压力测量信号所有的局部极小值点和局部极大值点；
S12、采用三次样条曲线分别连接所有的局部极小值点和局部最大值点，得到的下包络线和上包络线，计算上包络线和下包络线的均值：
(1)
S13、从信号中减去，得到差值信号：
(2)
如果满足本征模态函数的两个条件，则为的第一个本征模态函数分量；否则，令，重复步骤S11至步骤S13的计算过程 k’次，直到第 k’次得到的差值信号满足本征模态函数的两个条件，此时的第一个本征模态函数分量为：
(3)
S14、从信号中减去，得到残余信号：
(4)
令，重复步骤S11到步骤S14的计算过程 i次，得到第 i个本征模态函数分量为：
(5)
式中，为第 i个本征模态函数分量，为第 i次得到的残余信号；
继续上述分解过程，直到最终的残余信号成为单调函数或者只包含一个极值点，此时从中无法再分解出本征模态函数分量；综合式(4)和式(5)，动态压力测量信号表示为：
(6)
因此，动态压力测量信号被分解为 h个本征模态函数分量和一个残余分量，并且这些分量的频带由高到低变化，最终动态压力测量信号表示为：
(7)
计算经验模态分解得到的所有分量信号与动态压力测量信号之间的相关系数：
(8)
式中，和表示振动信号和第 i个模态分量对应的离散信号；和分别表示和的平均值； N为信号长度；
与动态压力测量信号之间相关系数达到设定阈值以上，就认为该分量为动态压力测量信号的振铃分量。
接上述技术方案，本征模态函数的两个条件如下：
（1）极值点与过零点数量在整个数据集上相等或最多相差一个；
（2）任意时间上，由局部极大值点估计的上包络线和局部极小值点估计的下包络线的均值为零。
接上述技术方案，设定阈值为60%。
接上述技术方案，卷积神经网络CNN包含输入层、卷积层、池化层、全连接层和输出层。
接上述技术方案，设卷积层输入为，其中A和B分别为从输入层输入的长和宽，则卷积层输出表示为：
(9)
式中，“*”为卷积运算，表示第 cn个映射， cn表示滤波核个数，为当前卷积层的第 cn个滤波核的权重矩阵，第 cn个滤波核的偏置为，为整流线性单元ReLU激活函数。
接上述技术方案，池化表示为：
(10)
式中，为乘法偏置项，为输入，表示池化操作， b为加性偏置向量，为激活函数。
接上述技术方案，全连接层表示为：
(11)
式中， P为上一个池化层的输出，为当前全连接层的输出，和分别为权值和加性偏置项，为激活函数。
接上述技术方案，选取振铃频率幅值大的未发生模态混叠的连续时间序列的序列数据训练模型。
本发明还提供一种激波流场动态压力测量信号振铃分量提取校正系统，包括：
振铃分量提取模块，用于采用经验模态分解方法对激波流场动态压力测量信号进行自适应分解，提取出激波流场动态压力测量信号的振铃分量；
振铃分量处理模块，用于从所提取出的振铃分量中选取未发生模态混叠的连续时间序列的序列数据，将未发生模态混叠的连续时间序列的序列数据分解组合为 K个等长的连续时间子序列，；其中，每个子序列的序列数据量为 d+1个，每个子序列如下所示：
式中，为未发生模态混叠的连续时间序列的第 k个序列数据；
模型构建训练模块，用于利用卷积神经网络来构建振铃分量畸变校正模型，并将子序列的 d+1个序列数据作为模型的输入特征，将未发生模态混叠的连续时间序列的第 k+d+zt个序列数据作为模型的输出特征，训练模型；其中， zt为时间跨度；
振铃分量校正模块，用于以振铃分量中发生模态混叠的第一个数据点为畸变点，以畸变点前第 d+zt个数据点为起始点，构建输入特征，并输入至振铃分量畸变校正模型，得到畸变校正数据，之后利用畸变校正数据替换掉畸变点；重复该步骤，获取最终的激波流场动态压力响应信号振铃分量。
本发明与现有技术相比，具有以下优点及有益效果：
本发明综合考虑了传统方法对振铃分量提取过程中产生的模态混叠现象，将自适应信号分解方法和卷积神经网络模型辨识方法的优点相结合，提出一种激波流场动态压力测量信号振铃分量准确提取及校正方法。该方法基于振铃分量的周期震荡衰减特性以及模态混叠影响，采用一种交错跨度式数据集构建方式，以此来扩充样本量以及增强预测的准确性；为消除由于模态混叠而导致的信号缺失等问题，利用卷积神经网络来构建振铃分量畸变校正模型，从而消除信号分解带来的模态混叠问题，实现振铃分量信号的完整提取。
附图说明
图1为本发明实施例1的激波流场动态压力测量信号振铃分量提取校正方法流程图；
图2为本发明实施例3的激波流场动态压力测量信号振铃分量提取校正方法流程图；
图3为本发明实施例5中压力传感器的响应信号图；
图4为本发明实施例5中自适应分解提取的振铃分量图；
图5为本发明实施例5中畸变校正模型的训练结果图；
图6为本发明实施例5中畸变校正模型的测试结果图；
图7为本发明实施例5中提取的完整振铃分量图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
实施例1
如图1所示，本发明实施例激波流场动态压力测量信号振铃分量提取校正方法主要包括以下步骤：
S1、采用经验模态分解方法对激波流场动态压力测量信号进行自适应分解，提取出激波流场动态压力测量信号的振铃分量；
S2、从所提取出的振铃分量中选取未发生模态混叠的连续时间序列的序列数据，将未发生模态混叠的连续时间序列的序列数据分解组合为 K个等长的连续时间子序列，；其中，每个子序列的序列数据量为 d+1个，每个子序列如下所示：
式中，为未发生模态混叠的连续时间序列的第 k个序列数据；
S3、利用卷积神经网络来构建振铃分量畸变校正模型，并将子序列的 d+1个序列数据作为该振铃分量畸变校正模型的输入特征，将未发生模态混叠的连续时间序列的第 k+ d+zt个序列数据作为模型的输出特征，训练模型；其中， zt为时间跨度；
S4、采集待校正的激波流场动态压力测量信号，并进行经验模态分解，找到振铃分量中发生模态混叠的第一个数据点作为畸变点，以畸变点前第 d+zt个数据点作为起始点，构建输入特征，并输入至训练好的振铃分量畸变校正模型，得到畸变校正数据，之后利用畸变校正数据替换掉畸变点；重复该步骤，获取最终的激波流场动态压力响应信号振铃分量。
具体地，步骤S1具体包括以下步骤：
S11、识别激波流场动态压力测量信号所有的局部极小值点和局部极大值点；
S12、采用三次样条曲线分别连接所有的局部极小值点和局部最大值点，得到的下包络线和上包络线，计算上包络线和下包络线的均值：
(1)
S13、从信号中减去，得到差值信号：
(2)
如果满足本征模态函数的两个条件，则为的第一个本征模态函数分量；否则，令，重复步骤S11至步骤S13的计算过程 k’次，直到第 k’次得到的差值信号满足本征模态函数的两个条件，此时的第一个本征模态函数分量为：
(3)
S14、从信号中减去，得到残余信号：
(4)
令，重复步骤S11到步骤S14的计算过程 i次，得到第 i个本征模态函数分量为：
(5)
式中，为第 i个本征模态函数分量，为第 i次得到的残余信号；
继续上述分解过程，直到最终的残余信号成为单调函数或者只包含一个极值点，此时从中无法再分解出本征模态函数分量；综合式(4)和式(5)，动态压力测量信号表示为：
(6)
因此，动态压力测量信号被分解为 h个本征模态函数分量和一个残余分量，并且这些分量的频带由高到低变化，最终动态压力测量信号表示为：
(7)
计算经验模态分解得到的所有分量信号与动态压力测量信号之间的相关系数：
(8)
式中，和表示振动信号和第 i个模态分量对应的离散信号；和分别表示和的平均值； N为信号长度；
与动态压力测量信号之间相关系数达到设定阈值以上，就认为该分量为动态压力测量信号的振铃分量。
该实施例主要将自适应信号分解方法和卷积神经网络模型辨识方法的优点相结合，可有效减小激波流场动态压力响应信号振铃分量提取过程受频率混叠的影响，通过振铃分量的变化特征构建数据预测模型，实现激波流场动态压力响应信号振铃分量的准确提取与校正。
实施例2
本发明的激波流场动态压力测量信号振铃分量准确提取与校正方法，包括如下步骤：
（1）动态压力响应信号自适应分解
采用经验模态分解方法对激波管动态压力响应信号（即动态压力测量信号）进行自适应分解，得到一系列分量，根据各分量信号与动态压力响应信号之间的相关性识别并提取出振铃分量。具体实现过程如下：
采用经验模态分解方法对激波管动态压力测量信号进行自适应分解，可表示为：
(1)
式中，，为分解的本征模态分量，为趋势分量。
计算各分量与动态压力测量信号之间的相关系数，剔除无关分量，提取振铃分量。
（2）振铃分量畸变校正模型建立
激波流场动态压力测量信号具有高频（振铃频率在300kHz以上）、非平稳特点，在对这种高频、非平稳信号进行EMD分解时，往往出现模态混叠问题，即部分时间区间内信号出现严重畸变，导致激波流场动态压力测量精度低，因此本发明先将激波流场动态压力测量信号进行振铃分量提取，然后分别将其中未发生模态混叠的振铃分量提取出来，作为校正样本。
针对步骤（1）中提取的振铃分量，选取振铃频率幅值较大且未发生模态混叠段构建畸变校正模型数据集，数据集构建采用交错跨度式构建方法。最后，将数据集分为训练集和测试集，构建基于卷积神经网络的振铃分量畸变校正模型。具体实现过程如下：
针对步骤（1）中提取的振铃分量，由于未发生模态混叠段数据样本量不足，无法利用神经网络进行大量训练，且经EMD分解出的振铃分量在靠近模态混叠段的数据并不能保证未受模态混叠影响。考虑到振铃分量呈周期性震荡衰减特性，对其预测采取间隔性预测不会影响其准确性，于是充分考虑振铃分量的周期震荡衰减特性以及模态混叠影响，提出一种交错跨度式数据集构建方式，以此来扩充样本量以及增强预测的准确性，增强畸变校准模型数据集，提升畸变校正效果。假设选取未发生模态混叠的连续时间序列的数据 n，将 n分解组合为K个等长的连续时间子序列，每个子序列为，，子序列数量为 d+1个，每个子序列如式（2）所示：
(2)
为保证预测数据的准确性，将的 d+1个数据作为输入特征inp，第 k+d+zt个序列数据作为输出特征otp， zt为时间跨度。定义前P个子序列的输入输出为训练集，P+1~K个子序列的输入输出作为测试集。因此，训练集和测试集输入数据分别( d+1)×P和( d+1)×（K-P），输出数据分别为1×P和1×（K-P），最后将其转化为卷积神经网络输入方式输入到神经网络来构建振铃分量畸变矫正模型。
（3）振铃分量提取及校正
利用步骤（2）中构建的振铃分量畸变校正模型，对模态混叠段时间序列进行预测，实现对动态压力响应信号振铃分量的准确提取。具体实现过程如下：
如步骤（2）中交错跨度式数据集构建方式，以发生模态混叠数据点为终点，其前 d+ zt个数据点作为起始点，按步骤（2）中数据集构建方法来重构“第一段”发生模态混叠的振铃分量，得到畸变校正数据，并替换畸变数据，完成第一段共 zt个数据点的畸变矫正；再以的最后一个数据为终点，其前 d+zt个数据点作为起始点，构建数据集并输入到振铃分量畸变校正模型中，来重构“第二段”发生模态混叠的振铃分量，得到畸变校正数据，并以此类推，实现激波管动态压力测量信号的缺失段振铃分量提取。最终得到重构的完整振铃分量。
实施例3
如图2所示，本发明实施例的激波流场动态压力测量信号振铃分量准确提取与校正方法，包括以下步骤：
（1）动态压力测量信号自适应分解；
（2）振铃分量畸变校正模型建立；
（3）振铃分量提取及校正。
1、动态压力测量信号自适应分解
采用经验模态分解对动态压力测量信号进行处理，得到一系列窄带分量，称为本征模态函数，每个本征模态函数必须满足以下两个条件：（1）极值点与过零点数量在整个数据集上相等或最多相差一个；（2）任意时间上，由局部极大值点估计的上包络线和局部极小值点估计的下包络线的均值为零。分解的基本步骤如下：
步骤（1）：识别动态压力测量信号的局部极小值点和局部极大值点；
步骤（2）：采用三次样条曲线分别连接所有的局部极小值点和局部最大值点，得到的下包络线和上包络线，计算上包络线和下包络线的均值为：
(1)
步骤(3)：从信号中减去，得到差值信号为：
(2)
如果满足本征模态函数的两个条件，则为的第一个本征模态函数分量；否则，令，重复步骤(1)至步骤(3)计算过程 k’次，直到得到的满足本征模态函数的两个条件，此时的第一个本征模态函数分量为：
(3)
步骤(4)：从动态压力测量信号中减去，得到残余信号为：
(4)
令，重复步骤(1)到步骤(4)计算过程 i次，得到第 i个本征模态函数分量为：
(5)
继续上述分解过程，直到最终的残余分量成为单调函数或者只包含一个极值点，此时从中无法再分解出本征模态函数分量。综合式(4)和式(5)，动态压力测量信号表示为：
(6)
因此，动态压力测量信号被分解为个 h个本征模态函数分量和一个残余分量，并且这些分量的频带由高到低变化。最终动态压力测量信号可表示为：
(7)
计算经验模态分解得到的所有分量信号与动态压力测量信号之间的相关系数：
(8)
式中，和表示振动信号和第 i个模态分量对应的离散信号；和分别表示和的平均值； N为信号长度。
根据相关系数意义，与动态压力测量信号之间相关系数达到60%以上，就认为该分量为动态压力测量信号的振铃分量。
2、振铃分量畸变校正模型建立
针对步骤（1）中初步提取的振铃分量，由于未发生模态混叠段数据样本量不足，无法利用神经网络进行大量训练，且经EMD分解出的振铃分量在靠近模态混叠段的数据并不能保证未受模态混叠影响。考虑到振铃分量呈周期性震荡衰减特性，对其预测采取间隔性预测不会影响其准确性，于是充分考虑振铃分量的周期震荡衰减特性以及模态混叠影响，提出一种交错跨度式数据集构建方式，以此来扩充样本量以及增强预测的准确性。为消除由于模态混叠而导致的信号缺失等问题，利用卷积神经网络来构建振铃分量畸变校正模型，从而消除信号分解带来的模态混叠问题，实现振铃分量信号的完整提取。卷积神经网络（Convolutional Neural Network，CNN）是一种深度前馈神经网络，用于处理具有网状结构的数据。CNN通过构造多个滤波器，具有强大的特征提取能力，并利用这些滤波器从输入数据层中提取具有代表性的特征。将稀疏连接与参数权重共享机制相结合，对数据维度进行时间和空间上的下采样，大大减少了训练参数的数量，有效避免了算法的过拟合，在信号预测方面具有优越的性能。
CNN包含输入层、卷积层、池化层、全连接层和输出层。其中，对于卷积层和池化层以及全连接层的设置至关重要。
对于卷积层，可根据需要来定义滤波核个数，在卷积层中，其核参数与输入数据点进行卷积。设卷积层输入为,其中A和B分别为从输入层输入的长和宽，则卷积层输出可表示为：
(9)
式中，“*”为卷积运算，表示第 cn个映射， cn表示滤波核个数，为当前卷积层的第 cn个滤波核的权重矩阵，第 cn个滤波核的偏置为。为整流线性单元ReLU激活函数。
经过卷积层后，输出特征图的维数会增加，在此基础上，使用池化层来降低输出特征图的维数。在池化层，可通过最大池化（max-pooling）或者平均池化（average-pooling）等方法来降低卷积层获得的特征图的维数。该池化过程可表示为：
(10)
式中，为乘法偏置项，为输入，表示池化操作， b为加性偏置向量，为激活函数。
与传统神经网络一样，采用全连接层的多个神经元来对其输入进行非线性拟合，所有神经元都连接到最后一个池化层的特征图的所有数据点，其操作可表示为：
(11)
式中， P为上一个池化层的输出，为当前全连接层的输出，和分别为权值和加性偏置项，为激活函数。
本方法将时间序列的输入转化为矩阵形式，单个矩阵输入为[q×1×1]，对于卷积层，可根据需要定义卷积核个数，卷积核尺寸定义为[p×1]。
在卷积层之后，输出的维数增加，因此，要使用池化层来降低输出特征图的维数，通过最大池化函数来对减少亢余特征，以进一步避免过度拟合并减少网络参数，特别的，由于卷积核尺寸为[p×1]，则池化核尺寸为[k×1]。
在多次交替卷积和池化操作之后，输出矩阵可最终由全连接层的多个神经元对输入进行非线性拟合来获得。具体数据集构建过程如下：
假设选取未发生模态混叠的连续时间序列的数据 n，将 n分解组合为K个等长的连续时间子序列，每个子序列为，,子序列数量为 d+1个，每个子序列如式12所示：
(12)
为保证预测数据的准确性，将的 d+1个数据作为输入特征inp，第 k+d+zt个序列数据作为输出特征otp， zt为时间跨度。定义前P个子序列的输入输出为训练集，P+1~K个子序列的输入输出作为测试集。因此，训练集和测试集输入数据分别( d+1)×P和( d+1)×（K-P），输出数据分别为1×P和1×（K-P），最后将其转化为卷积神经网络输入方式输入到卷积神经网络来构建振铃分量畸变矫正模型。
例如：当 k等于1， d等于4， zt等于2时，即利用未发生模态混叠的连续时间序列的第1、第2、第3、第4和第5个时序数据，来得到第7个时序数据。当 k等于2， d等于4， zt等于2时，即利用未发生模态混叠的连续时间序列的第2、第3、第4、第5和第6个时序数据，来得到第8个时序数据。
3、振铃分量提取及校正
由于子序列长度以及时间数据跨度限制，对于振铃分量模态混叠段信号无法一次性进行预测，需要根据子序列长度以及时间跨度来决定一次预测的时间序列长度。
如步骤（2）中交错跨度式数据集构建方式，以发生模态混叠数据点为终点，其前 d+ zt个数据点作为起始点，按步骤（2）中数据集构建方法来重构“第一段”发生模态混叠的振铃分量，得到畸变校正数据，并替换畸变数据，完成第一段共 zt个数据点的畸变矫正；再以的最后一个数据为终点，其前 d+zt个数据点作为起始点，构建数据集并输入到振铃分量畸变校正模型中，来重构“第二段”发生模态混叠的振铃分量，得到畸变校正数据，并以此类推，
并以此类推，最终实现激波管动态压力测量信号的缺失振铃分量提取。
(12)
最终，完成重构的完整振铃分量。
实施例4
本发明实施例的激波流场动态压力测量信号振铃分量提取校正系统，主要用于实现上述方法实施例，该系统包括：
振铃分量提取模块，用于采用经验模态分解方法对激波流场动态压力测量信号进行自适应分解，提取出激波流场动态压力测量信号的振铃分量；
振铃分量处理模块，用于从所提取出的振铃分量中选取未发生模态混叠的连续时间序列的序列数据，将未发生模态混叠的连续时间序列的序列数据分解组合为 K个等长的连续时间子序列，；其中，每个子序列的序列数据量为 d+1个，每个子序列如下所示：
式中，为未发生模态混叠的连续时间序列的第 k个序列数据；
模型构建训练模块，用于利用卷积神经网络来构建振铃分量畸变校正模型，并将子序列的 d+1个序列数据作为模型的输入特征，将未发生模态混叠的连续时间序列的第 k+d+zt个序列数据作为模型的输出特征，训练模型；其中， zt为时间跨度；
振铃分量校正模块，用于以振铃分量中发生模态混叠的第一个数据点为畸变点，以畸变点前第 d+zt个数据点为起始点，构建输入特征，并输入至振铃分量畸变校正模型，得到畸变校正数据，之后利用畸变校正数据替换掉畸变点；重复该步骤，获取最终的激波流场动态压力响应信号振铃分量。
系统的各个模块主要用于实现上述实施例的方法功能，在此不一一赘述。
实施例5
下面以ENDEVCO 8510B PR压力传感器测量由激波管系统产生的动态压力响应数据进行分析，其中数据的采样频率为5 MHz。
1、ENDEVCO 8510B PR压力传感器原始测量数据，见图3。
2、利用实施例3的步骤（1）对图3中的响应信号进行自适应分解，初步提取振铃分量，见图4。
3、利用实施例3的步骤（2）得到振铃分量畸变校正模型，训练及测试结果见图5和图6。
4、利用实施例3的步骤（3）实现完整振铃分量提取及校正，见图7。
上述方法可有效减小激波流场动态压力响应信号振铃分量提取过程受频率混叠的影响，通过振铃分量的变化特征构建数据预测模型，实现激波流场动态压力响应信号振铃分量的准确提取。
需要指出，根据实施的需要，可将本申请中描述的各个步骤/部件拆分为更多步骤/部件，也可将两个或多个步骤/部件或者步骤/部件的部分操作组合成新的步骤/部件，以实现本发明的目的。
本领域的技术人员容易理解，以上仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种激波流场动态压力测量信号振铃分量提取校正方法，其特征在于，包括以下步骤：
S1、采用经验模态分解方法对激波流场动态压力测量信号进行自适应分解，提取出激波流场动态压力测量信号的振铃分量；
S2、从所提取出的振铃分量中选取未发生模态混叠的连续时间序列的序列数据，将未发生模态混叠的连续时间序列的序列数据分解组合为K个等长的连续时间子序列，；其中，每个子序列的序列数据量为d+1个，每个子序列如下所示：
式中，为未发生模态混叠的连续时间序列的第k个序列数据；
S3、利用卷积神经网络来构建振铃分量畸变校正模型，并将子序列的d+1个序列数据作为该振铃分量畸变校正模型的输入特征，将未发生模态混叠的连续时间序列的第k+d+zt个序列数据作为模型的输出特征，训练模型；其中，zt为时间跨度；
S4、采集待校正的激波流场动态压力测量信号，并进行经验模态分解，找到振铃分量中发生模态混叠的第一个数据点作为畸变点，以畸变点前第d+zt个数据点作为起始点，构建输入特征，并输入至训练好的振铃分量畸变校正模型，得到畸变校正数据，之后利用畸变校正数据替换掉畸变点；重复该步骤，获取最终的激波流场动态压力响应信号振铃分量。
2.根据权利要求1所述的激波流场动态压力测量信号振铃分量提取校正方法，其特征在于，步骤S1具体包括以下步骤：
S11、识别激波流场动态压力测量信号所有的局部极小值点和局部极大值点；
S12、采用三次样条曲线分别连接所有的局部极小值点和局部最大值点，得到的下包络线和上包络线，计算上包络线和下包络线的均值：
(1)
S13、从信号中减去，得到差值信号：
(2)
如果满足本征模态函数的两个条件，则为的第一个本征模态函数分量；否则，令，重复步骤S11至步骤S13的计算过程k’次，直到第k’次得到的差值信号满足本征模态函数的两个条件，此时的第一个本征模态函数分量为：
(3)
S14、从信号中减去，得到残余信号：
(4)
令，重复步骤S11到步骤S14的计算过程i次，得到第i个本征模态函数分量为：
(5)
式中，为第i个本征模态函数分量，为第i次得到的残余信号；
继续上述分解过程，直到最终的残余信号成为单调函数或者只包含一个极值点，此时从中无法再分解出本征模态函数分量；综合式(4)和式(5)，动态压力测量信号表示为：
(6)
因此，动态压力测量信号被分解为h个本征模态函数分量和一个残余分量，并且这些分量的频带由高到低变化，最终动态压力测量信号表示为：
(7)
计算经验模态分解得到的所有分量信号与动态压力测量信号之间的相关系数：
(8)
式中，和表示振动信号和第i个模态分量对应的离散信号；和分别表示和的平均值；N为信号长度；
与动态压力测量信号之间相关系数达到设定阈值以上，就认为该分量为动态压力测量信号的振铃分量。
3.根据权利要求2所述的激波流场动态压力测量信号振铃分量提取校正方法，其特征在于，本征模态函数的两个条件如下：
（1）极值点与过零点数量在整个数据集上相等或最多相差一个；
（2）任意时间上，由局部极大值点估计的上包络线和局部极小值点估计的下包络线的均值为零。
4.根据权利要求2所述的激波流场动态压力测量信号振铃分量提取校正方法，其特征在于，设定阈值为60%。
5.根据权利要求1所述的激波流场动态压力测量信号振铃分量提取校正方法，其特征在于，卷积神经网络CNN包含输入层、卷积层、池化层、全连接层和输出层。
6.根据权利要求5所述的激波流场动态压力测量信号振铃分量提取校正方法，其特征在于，设卷积层输入为，其中A和B分别为从输入层输入的长和宽，则卷积层输出表示为：
(9)
式中，“*”为卷积运算，表示第cn个映射，cn表示滤波核个数，为当前卷积层的第cn个滤波核的权重矩阵，第cn个滤波核的偏置为，为整流线性单元ReLU激活函数。
7.根据权利要求5所述的激波流场动态压力测量信号振铃分量提取校正方法，其特征在于，池化表示为：
(10)
式中，为乘法偏置项，为输入，表示池化操作，b为加性偏置向量，为激活函数。
8.根据权利要求5所述的激波流场动态压力测量信号振铃分量提取校正方法，其特征在于，全连接层表示为：
(11)
式中，P为上一个池化层的输出，为当前全连接层的输出，和分别为权值和加性偏置项，为激活函数。
9.根据权利要求1所述的激波流场动态压力测量信号振铃分量提取校正方法，其特征在于，选取振铃频率幅值大的未发生模态混叠的连续时间序列的序列数据训练模型。
10.一种激波流场动态压力测量信号振铃分量提取校正系统，其特征在于，包括：
振铃分量提取模块，用于采用经验模态分解方法对激波流场动态压力测量信号进行自适应分解，提取出激波流场动态压力测量信号的振铃分量；
振铃分量处理模块，用于从所提取出的振铃分量中选取未发生模态混叠的连续时间序列的序列数据，将未发生模态混叠的连续时间序列的序列数据分解组合为K个等长的连续时间子序列，；其中，每个子序列的序列数据量为d+1个，每个子序列如下所示：
式中，为未发生模态混叠的连续时间序列的第k个序列数据；
模型构建训练模块，用于利用卷积神经网络来构建振铃分量畸变校正模型，并将子序列的d+1个序列数据作为模型的输入特征，将未发生模态混叠的连续时间序列的第k+d+ zt个序列数据作为模型的输出特征，训练模型；其中，zt为时间跨度；
振铃分量校正模块，用于以振铃分量中发生模态混叠的第一个数据点为畸变点，以畸变点前第d+zt个数据点为起始点，构建输入特征，并输入至振铃分量畸变校正模型，得到畸变校正数据，之后利用畸变校正数据替换掉畸变点；重复该步骤，获取最终的激波流场动态压力响应信号振铃分量。
CN202310368863.3A 2023-04-10 2023-04-10 激波流场动态压力测量信号振铃分量提取校正方法及系统 Active CN116089800B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310368863.3A CN116089800B (zh) 2023-04-10 2023-04-10 激波流场动态压力测量信号振铃分量提取校正方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310368863.3A CN116089800B (zh) 2023-04-10 2023-04-10 激波流场动态压力测量信号振铃分量提取校正方法及系统

Publications (2)

Publication Number Publication Date
CN116089800A true CN116089800A (zh) 2023-05-09
CN116089800B CN116089800B (zh) 2023-06-27

Family

ID=86199453

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310368863.3A Active CN116089800B (zh) 2023-04-10 2023-04-10 激波流场动态压力测量信号振铃分量提取校正方法及系统

Country Status (1)

Country Link
CN (1) CN116089800B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140100798A1 (en) * 2012-10-05 2014-04-10 Siemens Energy, Inc. Turbine blade fatigue life analysis using non-contact measurement and dynamical response reconstruction techniques
CN110807524A (zh) * 2019-11-13 2020-02-18 大连民族大学 单通道信号盲源分离幅度校正方法
CN111222458A (zh) * 2020-01-06 2020-06-02 浙江工业大学 一种基于集合经验模态分解和卷积神经网络的滚动轴承故障诊断方法
DE102021000248A1 (de) * 2021-01-19 2021-03-18 Daimler Ag Verfahren zur Charakterisierung von strukturdynamischen Eigenschaften einer Struktur mittels einer Übertragungsfunktion in einem Frequenzbereich mittels eines Messsystems, sowie Messsystem
CN113051661A (zh) * 2021-02-19 2021-06-29 南京航空航天大学 基于微小容腔多动力学的高温气流动态总压智能软测量方法
US20210278330A1 (en) * 2020-03-06 2021-09-09 California Institute Of Technology Quantification of Shear Stress Field Using Tuft Visualization and Deep Learning
CN113987691A (zh) * 2021-11-25 2022-01-28 中国人民解放军国防科技大学 激波失稳的高精度混合计算方法、装置、设备和存储介质
CN114964620A (zh) * 2022-08-03 2022-08-30 中国航空工业集团公司北京长城计量测试技术研究所 机械传感器动态校准分析方法、装置、计算机设备和存储介质
US20220345831A1 (en) * 2019-10-03 2022-10-27 Ear Physics, Llc Endoscopic three-dimensional imaging systems and methods
CN115758246A (zh) * 2022-11-21 2023-03-07 东南大学 基于EMD与AlexNet的非侵入式负荷识别方法
CN115795282A (zh) * 2023-01-30 2023-03-14 武汉工程大学 激波管动态压力重构方法、装置、电子设备及存储介质

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140100798A1 (en) * 2012-10-05 2014-04-10 Siemens Energy, Inc. Turbine blade fatigue life analysis using non-contact measurement and dynamical response reconstruction techniques
US20220345831A1 (en) * 2019-10-03 2022-10-27 Ear Physics, Llc Endoscopic three-dimensional imaging systems and methods
CN110807524A (zh) * 2019-11-13 2020-02-18 大连民族大学 单通道信号盲源分离幅度校正方法
CN111222458A (zh) * 2020-01-06 2020-06-02 浙江工业大学 一种基于集合经验模态分解和卷积神经网络的滚动轴承故障诊断方法
US20210278330A1 (en) * 2020-03-06 2021-09-09 California Institute Of Technology Quantification of Shear Stress Field Using Tuft Visualization and Deep Learning
DE102021000248A1 (de) * 2021-01-19 2021-03-18 Daimler Ag Verfahren zur Charakterisierung von strukturdynamischen Eigenschaften einer Struktur mittels einer Übertragungsfunktion in einem Frequenzbereich mittels eines Messsystems, sowie Messsystem
CN113051661A (zh) * 2021-02-19 2021-06-29 南京航空航天大学 基于微小容腔多动力学的高温气流动态总压智能软测量方法
CN113987691A (zh) * 2021-11-25 2022-01-28 中国人民解放军国防科技大学 激波失稳的高精度混合计算方法、装置、设备和存储介质
CN114964620A (zh) * 2022-08-03 2022-08-30 中国航空工业集团公司北京长城计量测试技术研究所 机械传感器动态校准分析方法、装置、计算机设备和存储介质
CN115758246A (zh) * 2022-11-21 2023-03-07 东南大学 基于EMD与AlexNet的非侵入式负荷识别方法
CN115795282A (zh) * 2023-01-30 2023-03-14 武汉工程大学 激波管动态压力重构方法、装置、电子设备及存储介质

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
BOWEI QI ET AL.: "Application of EMD Combined with Deep Learning and Knowledge Graph in Bearing Fault", 《JOURNAL OF SIGNAL PROCESSING SYSTEMS》, pages 1 - 20 *
IRINA A. ET AL.: "Edge detection and machine learning for automatic flow srtuctures detection and tracking on schlieren and shadowgraph images", 《JOURNAL OF FLOW VISUALIZATION AND IMAGE PROCESSING》, pages 1 - 26 *
ZHENJIAN YAO ET AL.: "Empirical mode decomposition-adaptive least squares method for dynamic calibration of pressure sensors", 《MEASUREMENT SCIENCE AND TECHNOLOGY》, pages 1 - 10 *
ZHENJIAN YAO ET AL.: "Improved shock tube method for dynamic calibration of the sensitivity characteristic of piezoresistive pressure sensors", 《MEASUREMENT》, pages 1 - 11 *
姚贞建;王中宇;王辰辰;左思然;: "压力传感器动态特性参数不确定度评定", 北京航空航天大学学报, no. 08, pages 107 - 116 *
方晓汾;: "基于经验模态和神经网络的纯电动汽车集中式电驱系统振动信号分析与故障预测", 机械工程师, no. 07, pages 64 - 67 *

Also Published As

Publication number Publication date
CN116089800B (zh) 2023-06-27

Similar Documents

Publication Publication Date Title
CN109033719B (zh) 一种风力涡轮机轴承故障诊断方法
CN104038181B (zh) 一种基于nlms算法的自适应滤波器的构建方法
Park A time series data prediction scheme using bilinear recurrent neural network
CN116316591A (zh) 基于混合双向门控循环的短期光伏功率预测方法及系统
CN106374879A (zh) 一种基于cse有效的fir滤波器优化方法
CN111365624A (zh) 一种输卤管道泄漏检测的智能终端与方法
CN116089800B (zh) 激波流场动态压力测量信号振铃分量提取校正方法及系统
CN110415180A (zh) 一种基于小波卷积神经网络的sar图像去噪方法
Hou et al. Deep neural network denoising model based on sparse representation algorithm for ECG signal
CN114152440A (zh) 一种滚动轴承故障特征提取方法
CN116465623B (zh) 一种基于稀疏Transformer的齿轮箱寿命预测方法
Sakow et al. Design and optimisation of regression-type small phase shift FIR filters and FIR-based differentiators with optimal local response in LS-sense
CN109239006B (zh) 一种基于湿度补偿模型的物质识别方法、装置及存储介质
CN115753067B (zh) 一种基于数字孪生的强噪声下旋转机械故障诊断方法
CN110969238A (zh) 一种校准用电数据的方法及装置
CN113378485A (zh) 一种基于meemd与sg滤波的海杂波去噪方法及装置
CN110889197A (zh) 基于神经网络的自适应前馈式主动降噪方法、计算机可读存储介质、电子设备
CN111582132B (zh) 基于改进eemd与pcnn的气体泄漏信号降噪方法
CN112634391B (zh) 一种基于压缩感知的灰度图像深度重建及故障诊断系统
CN115438301A (zh) 一种基于ICEEMDAN二次分解耦合informer模型的设备运行趋势预测方法
CN112761843A (zh) 一种共轨喷油器故障诊断方法
CN110705373A (zh) 一种基于小波包分析和残差网络的槽电压信号分类方法
CN114235396B (zh) 一种齿轮减速器故障特征提取方法
CN114330850B (zh) 一种用于气候预测的异常相对倾向生成方法与系统
CN117404853B (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