离心泵压力脉动宽带频的定量识别方法与流程

本发明涉及水力机械领域，更具体地，涉及一种离心泵压力脉动宽带频的定量识别方法。

背景技术：

在水力机械运行过程中，由于叶轮转动、旋转部件与静止部件之间的动静干涉作用、空化流动、卡门涡街、旋转失速、流动分离以及回流等因素的影响，水力机械内部瞬态压力呈现非定常脉动特性。高幅值压力脉动会引发水力机械及其系统的振动和噪声，在压力脉动频率与结构固有频率一致时甚至会引起共振。精确定量识别压力脉动的频率成分和幅值，对于评价水力机械运行稳定性具有重要意义。现有技术中的压力脉动频域的评价体系，主要着眼于转频、叶频和尾水管涡带等单一频率成分。但大量实验结果表明，水力机械内部，尤其是离心泵，往往还存在具有较高压力幅值的连续频率带，对于宽带特性频率成分，目前还没有可靠的定量识别和评价方法。

技术实现要素：

为了解决上述问题，本发明实施例提供一种克服上述问题或者至少部分地解决上述问题的离心泵压力脉动宽带频的定量识别方法。

根据本发明实施例的第一方面，提供一种离心泵压力脉动宽带频的定量识别方法，该方法包括：通过记录离心泵的瞬态压力，获取压力系数的频域数据，并提取频域数据中的每一个宽带频；对于每一个宽带频，基于筛值迭代法，采用正态分布形态对宽带频的全部极大值点进行拟合，获得第一拟合函数的中心频率；根据中心频率和高峰频率计算偏差度，并根据偏差度确认每一个宽带频对应的分布形态种类；其中，分布形态种类包括对称态分布形态或偏态分布形态，高峰频率为宽带频的幅值对应的频率；采用与分布形态种类对应的拟合方法对每一个宽带频进行拟合，获得每一个宽带频对应的目标拟合函数，并基于目标拟合函数获取宽带频的形态定量表征参数；其中，所述筛值迭代法包括如下步骤：对所述极大值点进行拟合后，若可决系数r-square小于设定精度阈值，则舍去所述全部极大值点中的目标极大值点，获得剩余极大值点；所述目标极大值点为压力系数小于拟合曲线的第一设定倍数的点；对所述剩余极大值点进行拟合，直到r-square不小于所述设定精度阈值。

根据本发明实施例的第二方面，提供了一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行程序时实现如第一方面的各种可能的实现方式中任一种可能的实现方式所提供的离心泵压力脉动宽带频的定量识别方法。

本发明实施例提供的离心泵压力脉动宽带频的定量识别方法，根据偏差度来确定宽带频对应的分布形态种类，从而能够根据分布形态种类通过合理的拟合方法对宽带频进行拟合，获得目标拟合函数，提取出宽带频的形态定量表征参数；从而能够通过形态定量表征参数对离心泵压力脉动宽带频进行定量的识别，进而可建立离心泵压力脉动频率特性的完整定量评价体系。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些图获得其他的附图。

图1为本发明实施例提供的离心泵压力脉动宽带频的定量识别方法的流程示意图；

图2为本发明另一实施例提供的离心泵压力脉动宽带频的定量识别方法的流程示意图；

图3为本发明实施例提供的无量纲化的时域图；

图4为本发明实施例提供的加hanning(汉宁)窗后的时域图；

图5为本发明实施例提供的压力系数频域图；

图6为本发明实施例提供的除去倍频的频域图；

图7为本发明实施例提供的除去白噪声后宽带频定位示意图；

图8为本发明实施例提供的采用正态分布函数对所有极大值点进行第一次拟合的示意图；

图9为本发明实施例提供的采用正态分布函数对剩余极大值点进行第二次拟合的示意图；

图10为本发明实施例提供的采用正态分布函数对高峰频率两侧所有极大值点进行过幅值点的第一次拟合的示意图；

图11为本发明实施例提供的采用正态分布函数对高峰频率两侧剩余极大值点进行过幅值点的第二次拟合的示意图；

图12为本发明实施例提供的采用正态分布函数对高峰频率两侧剩余极大值点进行过幅值点的第三次拟合的示意图；

图13为本发明实施例提供的采用指数函数对高峰频率两侧所有极大值点进行过幅值点的第一次拟合的示意图；

图14为本发明实施例提供的采用指数函数对高峰频率两侧剩余极大值点进行过幅值点的第二次拟合的示意图；

图15为本发明实施例提供的采用指数函数对高峰频率两侧剩余极大值点进行过幅值点的第三次拟合的示意图；

图16为本发明实施例提供的宽带频定量识别示意图；

图17为本发明实施例提供的电子设备的实体结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明实施例提供一种离心泵压力脉动宽带频的定量识别方法，参见图1，该方法包括但不限于如下步骤：

步骤100、通过记录离心泵的瞬态压力，获取压力系数的频域数据，并提取频域数据中的每一个宽带频。

其中，可以采用实验或数值模拟的方法记录离心泵的瞬态压力，可以根据实际需求合理选择采样频率和记录时间。例如，采样频率可以为1024hz，记录时间可以为10s，那么，记录的数据点截取个数为2ⁿ，n为正整数；例如可以为8192。根据瞬态压力，可以获得压力系数的频域数据，具体可以为压力系数频域图。基于频域数据，可通过设定一个阈值，定位及提取频域数据中的每一个宽带频。

步骤200、对于每一个宽带频，基于筛值迭代法，采用正态分布形态对宽带频的全部极大值点进行拟合，获得第一拟合函数的中心频率。

具体地，在步骤100中提取了每一个宽带频后，对于每一个宽带频，均进行如下处理：获得宽带频的所有极大值点，采用正态分布形态进行拟合，正态分布形态函数可表示为：

其中，fi是相对频率，a0为拟合函数幅值，a拟合函数中心频率，σ是拟合函数标准差。

拟合后可以获得第一拟合函数，并可进一步获得第一拟合函数的中心频率。

步骤300，根据中心频率和高峰频率计算偏差度，并根据偏差度确认每一个宽带频对应的分布形态种类；其中，分布形态种类包括对称态分布形态或偏态分布形态，高峰频率为宽带频内最高幅值对应的频率。

具体地，可定义宽带频内最高幅值对应的频率为高峰频率。在步骤200中获得了第一拟合函数的中心频率后，可计算高峰频率与中心频率之间的偏差度，偏差度反映二者之间的偏差程度。并可进一步根据每个宽带频的偏差度确定每一个宽带频的分布形态种类。分布形态种类可以为对称态分布形态和偏态分布形态中的一个。

步骤400、采用与分布形态种类对应的拟合方法对每一个宽带频进行拟合，获得每一个宽带频对应的目标拟合函数，并基于目标拟合函数获取宽带频的形态定量表征参数。

具体地，在步骤300中获得每一个宽带频对应的分布形态种类后，可对分布形态种类为对称态分布形态的宽带频采用一类拟合方法，而对分布形态种类为偏态分布形态的宽带频采用另一类拟合方法，最终可以获得每一个宽带频对应的目标拟合函数。目标拟合函数应当是能够最大限度的包络宽带频的极大值概率的拟合函数。从而可进一步基于目标拟合函数来提取宽带频的形态定量表征参数。本发明实施例对形态定量表征参数的具体类型不作限定。

其中，上述步骤200中的筛值迭代法包括但不限于如下步骤：

步骤1、对极大值点进行拟合后，若可决系数r-square小于设定精度阈值，则舍去全部极大值点中的目标极大值点，获得剩余极大值点；目标极大值点为压力系数小于拟合曲线的第一设定倍数的点。

步骤2、对剩余极大值点进行拟合，直到r-square不小于设定精度阈值。

具体地，采用可决系数r-square来度量拟合精度。若r-square小于设定精度阈值(例如0.9)，则舍去目标极大值点，目标极大值点是压力系数小于第一设定倍数(例如0.9倍)的拟合曲线的点。并对剩余极大值点进行相应类型的分布形态(包括正态分布形态拟合或指数分布形态)的拟合，直到拟合精度不小于设定精度阈值。拟合完成后，可获得最终的拟合函数以及最终的拟合函数所对应的剩余极大值点。

因此，可采用上述筛值迭代法对步骤200中的全部极大值点进行正态分布形态的拟合，获得第一拟合函数，以及第一拟合函数对应的剩余极大值点。

本发明实施例提供的离心泵压力脉动宽带频的定量识别方法，根据偏差度来确定宽带频对应的分布形态种类，从而能够根据分布形态种类通过合理的拟合方法对宽带频进行拟合，获得目标拟合函数，提取出宽带频的形态定量表征参数；从而能够通过形态定量表征参数对离心泵压力脉动宽带频进行定量的识别，进而可建立离心泵压力脉动频率特性的完整定量评价体系。

基于上述实施例的内容，作为一种可选实施例，提供一种通过记录离心泵的瞬态压力，获取压力系数的频域数据的方法，包括但不限于如下步骤：

步骤101、记录离心泵的瞬态压力，对瞬态压力进行无量纲化，获得压力系数。

其中，对瞬态压力进行无量纲化，压力系数表示方法为：

其中，cp是压力系数，pi是瞬态的压力，[pa]，是平均压力，[pa]，u叶片出口边的周向速度，[m/s]。

步骤102、在时间域上对压力系数加上汉宁窗函数后进行频域变换，并对频域上的压力系数进行修正，对频率进行无量纲化，获得压力系数的频域数据。

具体地，在整个时间域上对压力系数加汉宁窗函数，采用快速傅里叶变换得到压力系数频域图(即压力系数的频域数据)；将频域上的压力系数乘以修正参数(例如)进行修正。对频率进行无量纲化，相对频率表示方法为：

fc＝fp/fr

其中，fc是相对频率，fp是压力脉动的真实频率，[hz]，fr是叶轮旋转频率，[hz]。

基于上述实施例的内容，作为一种可选实施例，在上述步骤102中获得压力系数的频域数据之后，还包括如下步骤：

步骤103、去除频域数据中的转频及转频的整数倍频。换言之，在频域图中找到转频的倍数，并将其去掉。

步骤104、获取频域数据中压力系数的平均值，并将平均值作为阈值，去除频域数据中小于阈值的压力系数的点，并在去除过程中通过调整阈值的大小，去除所述宽带频的两侧的白噪声。

具体地，获得频域数据中压力系数的平均值pave，令pave作为阈值，去掉频域图中压力系数小于阈值的点。另外，若宽带频两侧还有白噪声存在，则阈值增大到1.1倍，再次去除频域数据中压力系数小于阈值的点，直到没有白噪声存在。

基于上述实施例的内容，作为一种可选实施例，提供一种步骤200中采用正态分布形态对宽带频的全部极大值点进行拟合的方法，包括但不限于如下步骤：获得宽带频的全部极大值点，并基于筛值迭代法，采用正态分布形态对全部极大值点进行拟合，获得第一拟合函数。

基于上述实施例的内容，作为一种可选实施例，提供一种根据中心频率和高峰频率计算偏差度，并根据偏差度确认每一个宽带频对应的分布形态种类的方法，包括但不限于如下步骤：

步骤301、计算中心频率与高峰频率之间的差值，根据差值及第一拟合函数的标准差计算偏差度。

其中，偏差度可以为上述差值的绝对值与一定倍数(例如6倍)的标准差之间的比值。

步骤302、若偏差度小于偏差度阈值，则确认宽带频的分布形态种类为对称态分布形态；否则，宽带频的分布形态种类为偏态分布形态。

其中，偏差度阈值可以为3％。即，若偏差度小于3％，则分布形态种类为对称态分布形态，否则为偏态分布形态。

基于上述实施例的内容，作为一种可选实施例，分布形态种类为对称态分布形态，相应地，提供一种采用与分布形态种类对应的拟合方法对每一个宽带频进行拟合，获得每一个宽带频对应的目标拟合函数的方法，包括但不限于如下步骤：

步骤401、获得第一拟合函数的幅值点，并以中心频率为界线将宽带频的全部极大值点分为左右两部分；基于筛值迭代法，采用指数分布形态分别对左右两部分的极大值点进行过第一拟合函数的幅值点的拟合，获得第二拟合函数。

具体地，若为宽带频为对称态分布形态，得到第一拟合函数的幅值点。以中心频率为界线将宽带频的所有极大值点分为左右两部分，采用指数分布形态分别对左右两边的极大值点进行过幅值点的拟合，拟合具体可采用筛值迭代法，指数分布形态函数可表示为：

其中，b和c是指数拟合的系数，f0为高峰频率。

拟合完成后可获得第二拟合函数，及第二拟合函数对应的剩余极大值点。

步骤402、若第一拟合函数对应的剩余极大值点的数量大于第二拟合函数对应的剩余极大值点的数量，则将第一拟合函数作为宽带频的目标拟合函数；否则，将第二拟合函数作为宽带频的目标拟合函数。

即在第一拟合函数对应的剩余极大值点的数量大于第二拟合函数对应的剩余极大值点的数量时，采用正态分布形态拟合宽带频；否则，采用指数分布形态拟合宽带频。

基于上述实施例的内容，作为一种可选实施例，分布形态种类为偏态分布形态；相应地，提供一种采用与分布形态种类对应的拟合方法对每一个宽带频进行拟合，获得每一个宽带频对应的目标拟合函数的方法，包括但不限于如下步骤：

步骤401’、获得宽带频的幅值点，并以高峰频率为界线将宽带频的全部极大值点分为左右两部分。

具体地，若宽带频为偏态分布形态，则取宽带频的幅值点，并以高峰频率为界线将宽带频的所有极大值点分为左右两部分。

步骤402’、基于筛值迭代法，采用正态分布形态分别对左右两部分的极大值点进行过宽带频的幅值点的拟合，获得第三拟合函数；以及，基于筛值迭代法，采用指数分布形态分别对左右两部分的极大值点进行过宽带频的幅值点的拟合，获得第四拟合函数。

具体地，首先采用正态分布形态分别对左右两边的极大值点进行过幅值点的拟合，拟合采用筛值迭代法，获得第三拟合函数。然后采用指数分布形态分别对左右两边的极大值点进行过幅值点的拟合，拟合采用筛值迭代法，获得第四拟合函数。应当说明的是，本发明实施例对采用正态分布形态进行拟合和采用指数分布形态进行拟合的先后顺序不作限定。

步骤403’、若第三拟合函数对应的剩余极大值点的数量大于第四拟合函数对应的剩余极大值点的数量，则将第三拟合函数作为宽带频的目标拟合函数；否则，将第四拟合函数作为宽带频的目标拟合函数。

具体地，由于拟合采用的是筛值迭代法，因此拟合完成后的拟合函数对应有不同数量的极大值点。通过比较第三拟合函数对应的剩余极大值点的数量与第四拟合函数对应的剩余极大值点的数量，获得比较后的结果，从而根据比较后的结果来采用相应的分布形态拟合宽带频。即在第三拟合函数对应的剩余极大值点的数量大于第四拟合函数对应的剩余极大值点的数量时，采用正态分布形态拟合宽带频；否则，采用指数分布形态拟合宽带频。

基于上述实施例的内容，作为一种可选实施例，形态定量表征参数包括：带宽、带宽系数、带宽中心频率和宽带频幅值中的至少一个；相应地，提供一种步骤400中基于目标拟合函数获取宽带频的形态定量表征参数的方法，包括但不限于如下步骤：

作压力系数为第二设定倍数的宽带频的幅值的直线，获得直线与目标拟合函数的两个交点及两个交点之间的距离；带宽为两个交点之间的距离；带宽系数为带宽与宽带频的幅值之间的比值；若分布形态种类为对称态分布，则带宽中心频率为目标拟合函数的中心频率；若分布形态种类为偏态分布，则带宽中心频率为宽带频的高峰频率；所述宽带频幅值为所述目标拟合函数的最大值。

其中，第二设定倍数可以为0.5倍，本发明实施例对此不作限定。作压力系数等于0.5倍幅值的直线，得到直线与拟合函数的两个交点及两点间的距离，将该距离定义为宽带频的带宽，将该距离与幅值的比值定义为带宽系数，将对称态分布的中心频率或偏态分布的高峰频率定义为宽带中心频率；将目标拟合函数的最大值最为宽带频幅值。

为了对上述实施例进行说明，本发明另一实施例提供一种离心泵压力脉动宽带频的定量识别方法，参见图2，该方法包括但不限于如下步骤：

s1、记录水泵的瞬态压力，采样频率不低于1024hz，记录时间为10s，截取的数据点个数为2ⁿ，n＝13，对瞬态压力进行无量纲化，得到压力系数；

s2、在整个时间域上对压力系数加hanning(汉宁)窗函数，并进行频域变换，并对频域的压力系数进行修正，对频率进行无量纲化，得到相对频率；

s3、找到频域中所有转频及其整数倍频，并将其去除；

s4、得到频域中压力系数的平均值pave，令pave作为阈值，去掉频域图中压力系数小于阈值的点，若宽带频两侧还有白噪声存在，则阈值增大到1.1倍，直到没有白噪声存在；

s5、采用s4中方法调整阈值，定位所有的宽带频，并将每一个宽带频单独提取；

s6、得到宽带频所有极大值点，采用正态分布函数进行拟合，采用可决系数r-square度量拟合精度，若r-square小于0.9，则舍去压力系数小于0.9倍拟合曲线的点，并对剩余的极大值点进行正态分布函数拟合，直到拟合精度不小于0.9(以下简称筛值迭代法)；

s7、得到拟合函数的中心频率a＝11.8和标准差σ＝0.1984，定义宽带频最高幅值对应的频率为高峰频率f0＝11.75，计算高峰频率与中心频率的差值，定义偏差度为上述差值的绝对值与6倍标准差的比值，若偏差度小于3％，则为对称态分布，否则为偏态分布；

s8、若为对称态分布形态，得到拟合函数的幅值点，并以中心频率为界线将宽带频的所有极大值点分为左右两部分，采用指数分布形态分别对左右两边的极大值点进行过幅值点的拟合，拟合时采用筛值迭代法；

s9、若s6中剩余极大值点的个数多于s8中，则采用正态分布形态拟合宽带频，否则采用指数分布形态；

s10、若为偏态分布形态，取宽带频的幅值点(11.75,0.0017)，并以高峰频率为界线将宽带频的所有极大值点分为左右两部分，采用正态分布形态分别对左右两边的极大值点进行过幅值点的拟合，拟合时采用筛值迭代法；

s11、若为偏态分布形态，取宽带频的幅值点(11.75,0.0017)，并以高峰频率为界线将宽带频的所有极大值点分为左右两部分，采用指数分布形态分别对左右两边的极大值点进行过幅值点的拟合，拟合时采用筛值迭代法；

s12、若s10中剩余极大值点的个数多于s11中，则采用正态分布形态拟合宽带频，否则采用指数分布形态；

s13、若宽带频为对称态分布形态，采用s9中拟合方法，若为偏态分布形态，则采用s12中拟合方法，作压力系数等于0.5倍幅值的直线，得到直线与拟合函数的两个交点及两点间的距离，将该距离定义为宽带频的带宽，将该距离与幅值的比值定义为带宽系数，将对称态分布的中心频率和偏态分布的高峰频率定义为宽带中心频率，将拟合函数的最大值定义为宽带频幅值。

以下提供一个结合双吸离心泵压力室内的压力脉动数据的示例对本发明实施例提供的方法进行进一步说明：

如图3所示，采样频率为1024hz，记录的数据点截取个数为8192，对压力进行无量纲化，得到压力脉动时域图。

如图4所示，为防止能量泄露，在整个时间域上加hanning(汉宁)窗函数。

如图5所示，经过快速傅里叶变换得到压力系数频域图，将频域的压力系数乘以进行修正，对频率进行无量纲化，找到所有高幅值的倍频，确定白噪声的幅值。

如图6所示，去除转频及其整数倍频。

如图7所示，得到频域中压力系数的平均值pave，令pave作为阈值，去掉频域图中压力系数小于阈值的点，得到第一个宽带频的相对频率范围为11.0-12.6。

如图8所示，找到宽带频所有极大值点，并进行正态分布函数第一次拟合，r-square＝0.6987；如图9所示，去掉所有小于0.9倍拟合曲线的极大值点，并进行正态分布函数第二次拟合，r-square＝0.9155；正态分布中心频率a＝11.8，宽带频最高幅值点对应的频率f0＝11.75，标准差σ＝0.1984，偏差度为4.2％，判定该宽带频为偏态分布形态。

如图10所示，取宽带频的幅值点(11.75,0.0017)，并以高峰频率为界线将宽带频的所有极大值点分为左右两部分，采用正态分布形态分别对左右两边极大值点进行过幅值点的函数拟合，左边r-square＝0.1611，右边r-square＝0.4655；如图11所示，去掉所有小于0.9倍拟合曲线的极大值点，并进行正态分布函数第二次拟合，左边r-square＝0.7962，右边r-square＝0.8752；如图12所示，去掉所有小于0.9倍拟合曲线的极大值点，并进行正态分布函数第三次拟合，左边r-square＝0.9211，右边r-square＝0.9434。

如图13所示，取宽带频的幅值点(11.75,0.0017)，并以高峰频率为界线将宽带频的所有极大值点分为左右两部分，采用指数分布形态分别对左右两边极大值点进行过幅值点的函数拟合，左右r-square＝0.5552，右边r-square＝0.7457；如图14所示，去掉所有小于0.9倍拟合曲线的极大值点，并进行指数函数第二次拟合，左右r-square＝0.0.8936，右边r-square＝0.9254；如图15所示，去掉所有小于0.9倍拟合曲线的极大值点，并进行指数函数第三次拟合，左边r-square＝0.9487，右边r-square＝0.9481。

正态分布函数拟合剩余极大值点个数为58，指数函数拟合剩余极大值点个数为26，最终选择如图12所示的正态分布函数对宽带频进行拟合。

如图16所示，做直线cp＝0.5a0，得到直线与拟合函数的交点及两点间的距离，将该距离定义为宽带频的带宽w＝0.4115，将该距离与幅值的比值定义为带宽系数wc＝242.06，将偏态分布的高峰频率定义为宽带中心频率fmid＝11.75，宽带频幅值a0＝0.0017。

综上，本发明实施例提供的离心泵压力脉动宽带频的定量识别方法，包括，记录压力脉动系数时域信号；得到无量纲频域信号；消除基频及其倍频；确定阈值，消除白噪声，并依次定位宽带频段；基于筛值迭代法，采用正态分布形态函数拟合宽带频的极大值点，计算宽带频与拟合函数的偏差度，若小于3％判定为满足对称态分布形态，否则为偏态分布形态；对于对称态分布形态的宽带频，基于筛值迭代法，分别对中心频率两边的极大值点采用指数分布形态拟合，比较正态分布拟合和指数拟合剩余极大值点的数量，并选择剩余极大值点多的拟合方法；对于偏态分布形态的宽带频，基于筛值迭代法，分别对高峰频率两边的极大值点采用正态分布形态拟合与指数分布形态拟合，选择剩余极大值点多的拟合方法；确定带宽、带宽系数、宽带中心频率和宽带频幅值。该方法对离心泵压力脉动宽带频进行定量的识别，进而可建立离心泵压力脉动频率特性的完整定量评价体系。

本发明实施例提供了一种电子设备，如图17所示，该设备包括：处理器(processor)501、通信接口(communicationsinterface)502、存储器(memory)503和通信总线504，其中，处理器501，通信接口502，存储器503通过通信总线504完成相互间的通信。处理器501可以调用存储器503上并可在处理器501上运行的计算机程序，以执行上述各实施例提供的离心泵压力脉动宽带频的定量识别方法，例如包括：通过记录离心泵的瞬态压力，获取压力系数的频域数据，并提取频域数据中的每一个宽带频；对于每一个宽带频，基于筛值迭代法，采用正态分布形态对宽带频的全部极大值点进行拟合，获得第一拟合函数的中心频率；根据中心频率和高峰频率计算偏差度，并根据偏差度确认每一个宽带频对应的分布形态种类；其中，分布形态种类包括对称态分布形态或偏态分布形态，高峰频率为宽带频的幅值对应的频率；采用与分布形态种类对应的拟合方法对每一个宽带频进行拟合，获得每一个宽带频对应的目标拟合函数，并基于目标拟合函数获取宽带频的定量形态定量表征参数。

此外，上述的存储器503中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-onlymemory)、随机存取存储器(ram，randomaccessmemory)、磁碟或者光盘等各种可以存储程序代码的介质。

本发明实施例还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各实施例提供离心泵压力脉动宽带频的定量识别方法，例如包括：通过记录离心泵的瞬态压力，获取压力系数的频域数据，并提取频域数据中的每一个宽带频；对于每一个宽带频，基于筛值迭代法，采用正态分布形态对宽带频的全部极大值点进行拟合，获得第一拟合函数的中心频率；根据中心频率和高峰频率计算偏差度，并根据偏差度确认每一个宽带频对应的分布形态种类；其中，分布形态种类包括对称态分布形态或偏态分布形态，高峰频率为宽带频的幅值对应的频率；采用与分布形态种类对应的拟合方法对每一个宽带频进行拟合，获得每一个宽带频对应的目标拟合函数，并基于目标拟合函数获取宽带频的定量形态定量表征参数。

以上所描述的电子设备等实施例仅仅是示意性的，其中作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如rom/ram、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分方法。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。