CN113988144B - 一种科氏流量计固有频率获得方法 - Google Patents

一种科氏流量计固有频率获得方法 Download PDF

Info

Publication number
CN113988144B
CN113988144B CN202111608974.4A CN202111608974A CN113988144B CN 113988144 B CN113988144 B CN 113988144B CN 202111608974 A CN202111608974 A CN 202111608974A CN 113988144 B CN113988144 B CN 113988144B
Authority
CN
China
Prior art keywords
signal
index
spectrum
value
coriolis flowmeter
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
CN202111608974.4A
Other languages
English (en)
Other versions
CN113988144A (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.)
Equipment Design and Testing Technology Research Institute of China Aerodynamics Research and Development Center
Original Assignee
Equipment Design and Testing Technology Research Institute of China Aerodynamics Research and Development Center
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 Equipment Design and Testing Technology Research Institute of China Aerodynamics Research and Development Center filed Critical Equipment Design and Testing Technology Research Institute of China Aerodynamics Research and Development Center
Priority to CN202111608974.4A priority Critical patent/CN113988144B/zh
Publication of CN113988144A publication Critical patent/CN113988144A/zh
Application granted granted Critical
Publication of CN113988144B publication Critical patent/CN113988144B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02—Preprocessing
    • G06F2218/04—Denoising
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/76—Devices for measuring mass flow of a fluid or a fluent solid material
    • G01F1/78—Direct mass flowmeters
    • G01F1/80—Direct mass flowmeters operating by measuring pressure, force, momentum, or frequency of a fluid flow to which a rotational movement has been imparted
    • G01F1/84—Coriolis or gyroscopic mass flowmeters
    • G01F1/8409—Coriolis or gyroscopic mass flowmeters constructional details
    • G01F1/8436—Coriolis or gyroscopic mass flowmeters constructional details signal processing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction

Abstract

本发明公开了一种科氏流量计固有频率获得方法，涉及信号处理领域，包括：基于科氏流量计在自由衰减状态下的振动信号获得采样信号，对采样信号进行截短，得到截短信号，并对截短信号进行快速傅里叶变换，得到截短信号的频谱最大值的索引；分别先后对截短信号和采样信号的频谱进行两点插值计算，得到采样信号频谱偏移量和衰减因子初步估计值，从而得到降噪滤波器的初始参数；利用正弦信号的线性预测性质构造预测关系式，并通过降噪滤波器进一步抑制噪声的影响，得到降噪后的预测矩阵；求解矩阵得到预测系数，从而得到科氏流量计固有频率，本发明能够提高科氏流量计固有频率估计精度，减少计算量、提升固有频率获得速度。

Description

一种科氏流量计固有频率获得方法
技术领域
本发明涉及信号处理领域，具体地，涉及一种科氏流量计固有频率获得方法。
背景技术
数字驱动首先通过自激信号启振流量管，当振动幅值到达一定时，停止激励，进入零驱动状态，此时流量管处于自由衰减状态；然后采用参数估计算法获取流量管固有频率；最后根据频率合成驱动信号，并结合非线性幅值控制算法，驱动流量管稳幅振动。
该方案的关键点在于停止激励后，能够快速有效地获取流量管固有频率。当流量管处于自由衰减状态时，振动模型为单频衰减实信号，如式(1)所示。
Figure 100002_DEST_PATH_IMAGE001
（1）
式中：
Figure 100002_DEST_PATH_IMAGE002
、
Figure 100002_DEST_PATH_IMAGE003
、
Figure 100002_DEST_PATH_IMAGE004
和
Figure 100002_DEST_PATH_IMAGE005
分别表示采样信号初幅值、频率、初相位和衰减因子；
Figure 100002_DEST_PATH_IMAGE006
表示采样时刻点，
Figure 100002_DEST_PATH_IMAGE007
常为偶数，表示信号长度；
Figure 100002_DEST_PATH_IMAGE008
是均值为0，方差为
Figure 100002_DEST_PATH_IMAGE009
的加性高斯白噪声。
根据对信号的不同处理方式，科氏流量计固有频率获得方法主要可分为时域法和频域法两大类。时域法是在时域对信号进行计算得到固有频率，如线性预测法，STMB法等，此类方法通过构造矩阵估计固有频率，其思路简单、分辨率高，在中高信噪比条件下的固有频率估计精度高，但受噪声影响大，在低信噪比条件下的固有频率估计精度较低，且涉及到迭代求解矩阵，方法计算量较大、实时性较差。频域法将时域信号转换到频域进行频谱分析，可借助硬件实现，其计算速度快，且具有更强的抗噪性，但受频谱泄漏的影响，在中高信噪比或者信号频率较低条件下的固有频率估计精度较低，且频率分辨率较低，不适合固有频率特别低的科氏流量计。典型方法有迭代插值法、加窗插值法和频谱泄漏校正法等，这些方法在一定程度上抑制了频谱泄漏的影响，提高了固有频率估计精度，但仍有待提高。
针对稳幅振动状态下的科氏流量计参数估计问题，发明专利202011468681.6公布了一种时频结合的含噪信号参数估计新算法。该方法首先采用频域法中的快速傅里叶算法对采样信号进行处理，得到降噪滤波器参数，并通过时域法中的线性预测法构造了滤波后的预测矩阵，通过求解矩阵得到了信号参数。但是上述方法存在以下技术问题：其模型简单无法处理复杂问题，其计算量大。
发明内容
为提高科氏流量计固有频率的估计精度，并减少计算量、提升计算速度，结合频域法抗噪性好、计算量小和时域法参数估计精度高的特点，本发明提供了一种科氏流量计固有频率获得方法。
为实现上述目的，本发明提供了一种科氏流量计固有频率获得方法，所述方法包括：
采集科氏流量计在自由衰减状态下的振动信号获得采样信号；
对采样信号进行截短获得截短信号；
获得所述截短信号的频谱，提取所述截短信号频谱最大值的索引，获得第一索引；
在所述第一索引前后进行插值，分别获得第一插值点频谱值和第二插值点频谱值；
求取所述第一索引对应频谱值、所述第一插值点频谱值和所述第二插值点频谱值中的最大值获得最大频谱值，获得所述最大频谱值对应的第二索引；
基于所述第二索引获得所述采样信号频谱最大值对应的第三索引；
在所述第三索引的前后进行插值，分别获得第三插值点频谱值和第四插值点频谱值；
基于所述第三插值点频谱值和所述第四插值点频谱值，计算所述采样信号的频谱偏移量估计值和衰减因子的初步估计值；
基于所述频谱偏移量和衰减因子的初步估计值，计算获得滤波器参数，基于所述滤波器参数获得滤波器；
构造所述采样信号的预测关系式，将所述预测关系式带入所述滤波器获得预测矩阵；
求解所述预测矩阵获得所述科氏流量计的固有频率。
其中，本方法首先为抑制噪声的影响、提高计算速度，对采样信号进行截短，得到截短信号，并对截短信号进行快速傅里叶变换(Fast Fourier Transform, FFT)，得到截短信号的频谱最大值的索引；其次，分别先后对截短信号和采样信号的频谱进行两点插值计算，得到采样信号频谱偏移量和衰减因子初步估计值，从而得到降噪滤波器的初始参数；然后，利用正弦信号的线性预测性质构造预测关系式，并通过降噪滤波器进一步抑制噪声的影响，得到降噪后的预测矩阵；最后，利用最小二乘法求解矩阵，得到预测系数，从而得到科氏流量计的固有频率。
其中，本发明中的科氏流量计固有频率获得方法能够减少计算量，提高了科氏流量计的固有频率估计精度，具有结合频域法抗噪性好、计算量小和时域法参数估计精度高的特点。
优选的，本方法中所述采样信号为：
Figure 100002_DEST_PATH_IMAGE010
式中：
Figure 611447DEST_PATH_IMAGE002
、
Figure 263008DEST_PATH_IMAGE003
、
Figure 777166DEST_PATH_IMAGE004
和
Figure 906796DEST_PATH_IMAGE005
分别表示所述采样信号初幅值、频率、初相位和衰减因子；
Figure 100002_DEST_PATH_IMAGE011
表示采样时刻点，
Figure 675163DEST_PATH_IMAGE007
表示信号长度；
Figure 232046DEST_PATH_IMAGE008
是均值为0，方差为
Figure 233500DEST_PATH_IMAGE009
的加性高斯白噪声。
优选的，本方法中所述截短信号的长度为所述采样信号长度的一半。其中，这样设计的目的是减少计算量。
优选的，所述方法具体为：对采样信号进行截短，得到长度为采样信号长度一半的截短信号
Figure 100002_DEST_PATH_IMAGE012
，
Figure 100002_DEST_PATH_IMAGE013
，对所述截短信号进行快速傅里叶变换，提取所述截短信号频谱最大值的索引，获得第一索引。
优选的，本方法中所述截短信号频谱
Figure 100002_DEST_PATH_IMAGE014
的计算方式为：
Figure 100002_DEST_PATH_IMAGE015
Figure 100002_DEST_PATH_IMAGE016
其中，
Figure 100002_DEST_PATH_IMAGE017
表示截短信号，
Figure 100002_DEST_PATH_IMAGE018
表示截短信号频谱，
Figure 100002_DEST_PATH_IMAGE019
表示截短信号频谱索引，
Figure 100002_DEST_PATH_IMAGE020
为自然底数，
Figure 100002_DEST_PATH_IMAGE021
为虚数单位，
Figure 100002_DEST_PATH_IMAGE022
为第一索引，
Figure 100002_DEST_PATH_IMAGE023
表示取复数序列
Figure 445782DEST_PATH_IMAGE018
的模值，
Figure 100002_DEST_PATH_IMAGE024
表示提取序列
Figure 583502DEST_PATH_IMAGE023
中最大值的索引。
优选的，本方法中所述第一索引与其前后两个插值点的间隔均为0.5，所述第三索引与其前后两个插值点的间隔均为0.5。选择插值点间隔为0.5，是因为频谱偏移量
Figure 100002_DEST_PATH_IMAGE025
的取值边界是0.5(
Figure 100002_DEST_PATH_IMAGE026
)，利用插值点间隔0.5计算出来的频谱偏移量更加准确。
优选的，本方法中所述第一插值点频谱值为
Figure 100002_DEST_PATH_IMAGE027
，所述第二插值点频谱值为
Figure 100002_DEST_PATH_IMAGE028
，所述第二索引为
Figure 100002_DEST_PATH_IMAGE029
，所述第三索引为
Figure 100002_DEST_PATH_IMAGE030
，其中：
Figure 100002_DEST_PATH_IMAGE031
；
Figure 100002_DEST_PATH_IMAGE032
；
Figure 100002_DEST_PATH_IMAGE033
；
Figure 100002_DEST_PATH_IMAGE034
。
优选的，本方法中所述第三插值点频谱值为
Figure 100002_DEST_PATH_IMAGE035
，所述第四插值点频谱值为
Figure 100002_DEST_PATH_IMAGE036
，所述采样信号的频谱偏移量估计值为
Figure 100002_DEST_PATH_IMAGE037
，所述采样信号的衰减因子初步估计值为
Figure 100002_DEST_PATH_IMAGE038
，其中：
Figure 100002_DEST_PATH_IMAGE039
；
Figure 100002_DEST_PATH_IMAGE040
；
Figure 100002_DEST_PATH_IMAGE041
；
其中，
Figure 100002_DEST_PATH_IMAGE042
和
Figure 100002_DEST_PATH_IMAGE043
分别表示取复数
Figure 100002_DEST_PATH_IMAGE044
的实部和虚部。
优选的，本方法中所述滤波器为
Figure 100002_DEST_PATH_IMAGE045
，其中：
Figure 100002_DEST_PATH_IMAGE046
Figure 100002_DEST_PATH_IMAGE047
表示复变量，
Figure 100002_DEST_PATH_IMAGE048
和
Figure 100002_DEST_PATH_IMAGE049
为滤波器参数，计算方式为：
Figure 100002_DEST_PATH_IMAGE050
。
优选的，本方法中所述预测关系式为：
Figure 100002_DEST_PATH_IMAGE051
其中，
Figure 100002_DEST_PATH_IMAGE052
，
Figure 100002_DEST_PATH_IMAGE053
，
Figure 100002_DEST_PATH_IMAGE054
,
Figure 100002_DEST_PATH_IMAGE055
,
Figure 100002_DEST_PATH_IMAGE056
为单位冲击函数；
所述预测矩阵为：
Figure 100002_DEST_PATH_IMAGE057
其中，
Figure 100002_DEST_PATH_IMAGE058
和
Figure DEST_PATH_IMAGE059
分别表示采样信号
Figure DEST_PATH_IMAGE060
和冲击信号
Figure DEST_PATH_IMAGE061
经过滤波器后的滤波信号；
求解所述预测矩阵获得所述科氏流量计的固有频率
Figure DEST_PATH_IMAGE062
，其中：
Figure DEST_PATH_IMAGE063
。
本发明提供的一个或多个技术方案，至少具有如下技术效果或优点：
与发明专利202011468681.6公布的方法(以下称为传统方法)相比，本发明存在两个主要不同点：
1) 本方法与传统方法针对的信号模型不一样，传统方法是针对的稳幅振动状态下的科氏流量计振动信号，为幅值恒定信号，只包含幅值、频率和初相位3个参数，参数少、模型简单，而本发明针对的是自由衰减状态下的科氏流量计振动信号，为衰减信号，含有幅值、频率、初相位和衰减因子4个参数，模型更加复杂，增加了参数估计方法的设计难度，也无法直接使用传统方法。
2)在利用频域法求解滤波器参数时，传统方法对N点采样信号进行补零，得到2N点补零信号，再对补零信号进行处理，增加了计算量，但本发明方法没有直接采用这种方式，而是设计了更加巧妙的滤波器参数求解方法，通过对N点采样信号进行截短，得到N/2点截短信号，再对截短信号进行处理，利用较少的信号点得到了滤波器参数，降低了算法的计算量，提升了算法的实时性。
与其他针对衰减信号的参数估计算法相比：本发明中的科氏流量计固有频率获得方法，具有结合频域法抗噪性好、计算量小和时域法参数估计精度高的特点，能够减少计算量，提高了科氏流量计固有频率的估计精度。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本发明的一部分，并不构成对本发明实施例的限定；
图1为科氏流量计固有频率获得方法的流程示意图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在相互不冲突的情况下，本发明的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述范围内的其他方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
请参考图1，图1为科氏流量计固有频率获得方法的流程示意图，本发明提供了一种科氏流量计固有频率获得方法，所述方法包括：
采集科氏流量计在自由衰减状态下的振动信号获得采样信号；
对采样信号进行截短获得截短信号；
获得所述截短信号的频谱，提取所述截短信号频谱最大值的索引，获得第一索引；
在所述第一索引前后进行插值，分别获得第一插值点频谱值和第二插值点频谱值；
求取所述第一索引对应频谱值、所述第一插值点频谱值和所述第二插值点频谱值中的最大值获得最大频谱值，获得所述最大频谱值对应的第二索引；
基于所述第二索引获得所述采样信号频谱最大值对应的第三索引；
在所述第三索引的前后进行插值，分别获得第三插值点频谱值和第四插值点频谱值；
基于所述第三插值点频谱值和所述第四插值点频谱值，计算所述采样信号的频谱偏移量估计值和衰减因子的初步估计值；
基于所述频谱偏移量和衰减因子的初步估计值，计算获得滤波器参数，基于所述滤波器参数获得滤波器；
构造所述采样信号的预测关系式，将所述预测关系式带入所述滤波器获得预测矩阵；
求解所述预测矩阵获得所述科氏流量计的固有频率。
其中，本实施例中，振动信号的获得方式为：科氏流量计利用自激函数启振后，流量管幅值达到一定时，停止自激，使流量管处于自由衰减状态，并采集流量管此时的振动信号。
在本实施例中，本方法具体包括：
对采样信号进行截短，得到长度为采样信号长度一半的截短信号
Figure DEST_PATH_IMAGE064
，并进行快速傅里叶变换(Fast FourierTransform, FFT)，提取频谱最大值的索引。
Figure DEST_PATH_IMAGE065
(2)
Figure 560555DEST_PATH_IMAGE016
(3)
其中，
Figure 49305DEST_PATH_IMAGE017
表示截短信号，
Figure 520737DEST_PATH_IMAGE018
表示截短信号频谱，
Figure DEST_PATH_IMAGE066
表示截短信号频谱索引，
Figure 496653DEST_PATH_IMAGE020
为自然底数，
Figure 660918DEST_PATH_IMAGE021
为虚数单位，
Figure 636964DEST_PATH_IMAGE022
表示截短信号频谱最大值的索引，
Figure 646508DEST_PATH_IMAGE023
表示取复数序列
Figure 24400DEST_PATH_IMAGE018
的模值，
Figure 828408DEST_PATH_IMAGE024
表示提取序列
Figure 39553DEST_PATH_IMAGE023
中最大值的索引。
在截短信号
Figure 852788DEST_PATH_IMAGE017
频谱索引
Figure 85186DEST_PATH_IMAGE022
左边进行插值，插值间隔为0.5，得到插值点频谱值
Figure DEST_PATH_IMAGE067
。
Figure DEST_PATH_IMAGE068
(4)
在截短信号
Figure 794516DEST_PATH_IMAGE017
频谱索引
Figure 479576DEST_PATH_IMAGE022
右边进行插值，插值间隔为0.5，得到插值点频谱值
Figure 345769DEST_PATH_IMAGE028
。
Figure DEST_PATH_IMAGE069
(5)
通过比较索引
Figure 167095DEST_PATH_IMAGE022
处的频谱值
Figure DEST_PATH_IMAGE070
和其左右两边的插值点频谱值
Figure DEST_PATH_IMAGE071
、
Figure DEST_PATH_IMAGE072
，求取三个频谱值中最大值的索引
Figure DEST_PATH_IMAGE073
。
Figure DEST_PATH_IMAGE074
(6)
从而得到N点采样信号频谱最大值的索引
Figure DEST_PATH_IMAGE075
。
Figure DEST_PATH_IMAGE076
(7)
得到采样信号频谱最大值的索引后，在采样信号频谱索引
Figure 470163DEST_PATH_IMAGE075
左边插值，间隔为0.5，得到插值点频谱值
Figure 642518DEST_PATH_IMAGE035
。
Figure DEST_PATH_IMAGE077
(8)
在采样信号频谱索引
Figure 266397DEST_PATH_IMAGE075
右边插值，间隔为0.5，得到插值点频谱值
Figure 722655DEST_PATH_IMAGE036
。
Figure 304946DEST_PATH_IMAGE040
(9)
利用插值点频谱值，计算频谱偏移量和衰减因子的初步估计值
Figure DEST_PATH_IMAGE078
。
Figure DEST_PATH_IMAGE079
(10)
式中：
Figure 167860DEST_PATH_IMAGE037
和
Figure 126589DEST_PATH_IMAGE078
分别表示频谱偏移量
Figure 670309DEST_PATH_IMAGE025
估计值和衰减因子
Figure DEST_PATH_IMAGE080
的初步估计值，
Figure DEST_PATH_IMAGE081
和
Figure DEST_PATH_IMAGE082
分别表示取复数
Figure DEST_PATH_IMAGE083
的实部和虚部。
设计传递函数为
Figure 361185DEST_PATH_IMAGE045
的降噪滤波器：
Figure DEST_PATH_IMAGE084
(11)
式中：
Figure 226241DEST_PATH_IMAGE047
表示复变量，滤波器参数
Figure 988661DEST_PATH_IMAGE048
和
Figure 904664DEST_PATH_IMAGE049
由频谱偏移量和衰减因子的初步估计值的计算结果确定。
Figure DEST_PATH_IMAGE085
(12)
根据正弦函数的预测性质，构造预测关系式：
Figure DEST_PATH_IMAGE086
(13)
式中：
Figure 32020DEST_PATH_IMAGE052
，
Figure 416997DEST_PATH_IMAGE053
，
Figure 717528DEST_PATH_IMAGE054
,
Figure 753617DEST_PATH_IMAGE055
，
Figure 848612DEST_PATH_IMAGE056
为单位冲击函数。
将式(13)代入式(11)，可得预测矩阵：
Figure 704572DEST_PATH_IMAGE057
(14)
式中：
Figure 808795DEST_PATH_IMAGE058
和
Figure 417499DEST_PATH_IMAGE059
分别表示采样信号
Figure 948975DEST_PATH_IMAGE060
和冲击信号
Figure 292232DEST_PATH_IMAGE061
经过滤波器后的滤波信号。
利用最小二乘法计算式(14)，得到预测系数，从而得到所述科氏流量计的固有频率。
Figure 200145DEST_PATH_IMAGE063
(15)
得到固有频率后，合成驱动信号，驱动流量管稳幅振动。
尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (10)

1.一种科氏流量计固有频率获得方法，其特征在于，所述方法包括：
采集科氏流量计在自由衰减状态下的振动信号获得采样信号；
对采样信号进行截短获得截短信号；
获得所述截短信号的频谱，提取所述截短信号频谱最大值的索引，获得第一索引；
在所述第一索引前后进行插值，分别获得第一插值点频谱值和第二插值点频谱值；
求取所述第一索引对应频谱值、所述第一插值点频谱值和所述第二插值点频谱值中的最大值获得最大频谱值，获得所述最大频谱值对应的第二索引；
基于所述第二索引获得所述采样信号频谱最大值对应的第三索引；
在所述第三索引的前后进行插值，分别获得第三插值点频谱值和第四插值点频谱值；
基于所述第三插值点频谱值和所述第四插值点频谱值，计算所述采样信号的频谱偏移量估计值和衰减因子的初步估计值；
基于所述频谱偏移量和衰减因子的初步估计值，计算获得滤波器参数，基于所述滤波器参数获得滤波器；
构造所述采样信号的预测关系式，将所述预测关系式带入所述滤波器获得预测矩阵；
求解所述预测矩阵获得所述科氏流量计的固有频率。
2.根据权利要求1所述的科氏流量计固有频率获得方法，其特征在于，所述采样信号为：
Figure DEST_PATH_IMAGE001
式中：
Figure DEST_PATH_IMAGE002
、
Figure DEST_PATH_IMAGE003
、
Figure DEST_PATH_IMAGE004
和
Figure DEST_PATH_IMAGE005
分别表示所述采样信号初幅值、频率、初相位和衰减因子；
Figure DEST_PATH_IMAGE006
表示采样时刻点，
Figure DEST_PATH_IMAGE007
表示信号长度；
Figure DEST_PATH_IMAGE008
是均值为0，方差为
Figure DEST_PATH_IMAGE009
的加性高斯白噪声。
3.根据权利要求1所述的科氏流量计固有频率获得方法，其特征在于，所述截短信号的长度为所述采样信号长度的一半。
4.根据权利要求3所述的科氏流量计固有频率获得方法，其特征在于，所述方法具体为：对采样信号进行截短，得到长度为采样信号长度一半的截短信号
Figure DEST_PATH_IMAGE010
，
Figure DEST_PATH_IMAGE011
，对所述截短信号进行快速傅里叶变换，提取所述截短信号频谱最大值的索引，获得第一索引。
5.根据权利要求4所述的科氏流量计固有频率获得方法，其特征在于，所述截短信号频谱
Figure DEST_PATH_IMAGE012
的计算方式为：
Figure DEST_PATH_IMAGE013
Figure DEST_PATH_IMAGE014
其中，
Figure 483181DEST_PATH_IMAGE010
表示截短信号，
Figure 134742DEST_PATH_IMAGE012
表示截短信号频谱，
Figure DEST_PATH_IMAGE015
表示截短信号频谱索引，
Figure DEST_PATH_IMAGE016
为自然底数，
Figure DEST_PATH_IMAGE017
为虚数单位，
Figure DEST_PATH_IMAGE018
为第一索引，
Figure DEST_PATH_IMAGE019
表示取复数序列
Figure DEST_PATH_IMAGE020
的模值，
Figure DEST_PATH_IMAGE021
表示提取序列
Figure 6490DEST_PATH_IMAGE019
中最大值的索引。
6.根据权利要求1所述的科氏流量计固有频率获得方法，其特征在于，所述第一索引与其前后两个插值点的间隔均为0.5，所述第三索引与其前后两个插值点的间隔均为0.5。
7.根据权利要求5所述的科氏流量计固有频率获得方法，其特征在于，所述第一插值点频谱值为
Figure DEST_PATH_IMAGE022
，所述第二插值点频谱值为
Figure DEST_PATH_IMAGE023
，所述第二索引为
Figure DEST_PATH_IMAGE024
，所述第三索引为
Figure DEST_PATH_IMAGE025
，其中：
Figure DEST_PATH_IMAGE026
；
Figure DEST_PATH_IMAGE027
；
Figure DEST_PATH_IMAGE028
；
Figure DEST_PATH_IMAGE029
。
8.根据权利要求7所述的科氏流量计固有频率获得方法，其特征在于，所述第三插值点频谱值为
Figure DEST_PATH_IMAGE030
，所述第四插值点频谱值为
Figure DEST_PATH_IMAGE031
，所述采样信号的频谱偏移量估计值为
Figure DEST_PATH_IMAGE032
，所述采样信号的衰减因子初步估计值为
Figure DEST_PATH_IMAGE033
，其中：
Figure DEST_PATH_IMAGE034
；
Figure DEST_PATH_IMAGE035
；
Figure DEST_PATH_IMAGE036
；
其中，
Figure DEST_PATH_IMAGE037
和
Figure DEST_PATH_IMAGE038
分别表示取复数
Figure DEST_PATH_IMAGE039
的实部和虚部。
9.根据权利要求8所述的科氏流量计固有频率获得方法，其特征在于，所述滤波器为
Figure DEST_PATH_IMAGE040
，其中：
Figure DEST_PATH_IMAGE041
Figure DEST_PATH_IMAGE042
表示复变量，
Figure DEST_PATH_IMAGE043
和
Figure DEST_PATH_IMAGE044
为滤波器参数，计算方式为：
Figure DEST_PATH_IMAGE045
。
10.根据权利要求9所述的科氏流量计固有频率获得方法，其特征在于，所述预测关系式为：
Figure DEST_PATH_IMAGE046
其中，
Figure DEST_PATH_IMAGE047
，
Figure DEST_PATH_IMAGE048
，
Figure DEST_PATH_IMAGE049
,
Figure DEST_PATH_IMAGE050
,
Figure DEST_PATH_IMAGE051
为单位冲击函数；
所述预测矩阵为：
Figure DEST_PATH_IMAGE052
其中，
Figure DEST_PATH_IMAGE053
和
Figure DEST_PATH_IMAGE054
分别表示采样信号
Figure DEST_PATH_IMAGE055
和冲击信号
Figure DEST_PATH_IMAGE056
经过滤波器后的滤波信号；
求解所述预测矩阵获得所述科氏流量计的固有频率
Figure DEST_PATH_IMAGE057
，其中：
Figure DEST_PATH_IMAGE058
。
CN202111608974.4A 2021-12-27 2021-12-27 一种科氏流量计固有频率获得方法 Active CN113988144B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111608974.4A CN113988144B (zh) 2021-12-27 2021-12-27 一种科氏流量计固有频率获得方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111608974.4A CN113988144B (zh) 2021-12-27 2021-12-27 一种科氏流量计固有频率获得方法

Publications (2)

Publication Number Publication Date
CN113988144A CN113988144A (zh) 2022-01-28
CN113988144B true CN113988144B (zh) 2022-03-08

Family

ID=79734489

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111608974.4A Active CN113988144B (zh) 2021-12-27 2021-12-27 一种科氏流量计固有频率获得方法

Country Status (1)

Country Link
CN (1) CN113988144B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114584445B (zh) * 2022-03-08 2024-01-12 上海季丰电子股份有限公司 信号的频率估计方法、装置、电子设备和介质
CN115270893B (zh) * 2022-09-26 2022-12-06 中国空气动力研究与发展中心设备设计与测试技术研究所 一种高精度科氏流量计数字信号处理方法
CN115265691B (zh) * 2022-09-26 2022-12-09 中国空气动力研究与发展中心设备设计与测试技术研究所 一种科氏流量计振动频率跟踪方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5767665A (en) * 1994-09-13 1998-06-16 Fuji Electric Co. Ltd. Phase difference measuring apparatus and mass flowmeter thereof
CN102004186A (zh) * 2010-11-08 2011-04-06 北京航空航天大学 一种消除频谱泄漏的高准确度正弦信号测量方法与装置
CN105973349A (zh) * 2016-04-19 2016-09-28 温州科技职业学院 一种科氏流量计固有频率的在线预测方法
CN112964929A (zh) * 2021-01-14 2021-06-15 中国空气动力研究与发展中心设备设计与测试技术研究所 含噪多频衰减信号参数估计新算法
CN113156206A (zh) * 2020-12-07 2021-07-23 中国空气动力研究与发展中心设备设计与测试技术研究所 时频结合的含噪信号参数估计新算法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7404336B2 (en) * 2000-03-23 2008-07-29 Invensys Systems, Inc. Correcting for two-phase flow in a digital flowmeter
CN102128656B (zh) * 2011-02-25 2013-09-04 合肥工业大学 一种微弯型科氏质量流量计数字信号处理方法和系统
CN113341220B (zh) * 2021-08-05 2021-11-02 中国空气动力研究与发展中心设备设计与测试技术研究所 含噪多频衰减实信号频率估计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5767665A (en) * 1994-09-13 1998-06-16 Fuji Electric Co. Ltd. Phase difference measuring apparatus and mass flowmeter thereof
CN102004186A (zh) * 2010-11-08 2011-04-06 北京航空航天大学 一种消除频谱泄漏的高准确度正弦信号测量方法与装置
CN105973349A (zh) * 2016-04-19 2016-09-28 温州科技职业学院 一种科氏流量计固有频率的在线预测方法
CN113156206A (zh) * 2020-12-07 2021-07-23 中国空气动力研究与发展中心设备设计与测试技术研究所 时频结合的含噪信号参数估计新算法
CN112964929A (zh) * 2021-01-14 2021-06-15 中国空气动力研究与发展中心设备设计与测试技术研究所 含噪多频衰减信号参数估计新算法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Accurate frequency estimation of multiple complex and real sinusoids based on iterative interpolation;Ze-Long Mou 等;《Digital Signal Processing》;20210727;第117卷;第1-12页 *
COMPUTATIONAL AND EXPERIMENTAL STUDY OF FREQUENCIES AND NATURAL MODE OF WELDED SHELL OF CORIOLIS FLOWMETER;Yaushev, A.A. 等;《https://www.researchgate.net/publication/326474944》;20181231;第45-51页 *
含噪多频信号频率估计算法;陈鹏 等;《振动与冲击》;20210802;第40卷(第14期);第138-143页 *
基于自适应陷波器的科氏流量计信号频率估计方法研究;游海;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20170315(第03期);第C030-259页 *
科氏流量计振动信号参数估计的最小二乘法;杨辉跃 等;《计算机工程与应用》;20161231;第52卷(第3期);第252-255页 *
频谱泄漏校正的多频实信号频率估计算法;陈鹏 等;《振动与冲击》;20210826;第40卷(第15期);第107-113页 *
高精度的科氏流量计频率估计方法及应用;沈廷鳌 等;《仪表技术与传感器》;20191130(第11期);第37-44页 *

Also Published As

Publication number Publication date
CN113988144A (zh) 2022-01-28

Similar Documents

Publication Publication Date Title
CN113988144B (zh) 一种科氏流量计固有频率获得方法
CN114061678B (zh) 一种科氏流量计数字驱动方法
CN110617964A (zh) 用于滚动轴承故障诊断的同步压缩变换阶比分析法
CN110763462B (zh) 一种基于同步压缩算子的时变振动信号故障诊断方法
CN109285561B (zh) 一种基于自适应窗长的船舶螺旋桨空化噪声调制谱特征保真增强方法
CN105092241B (zh) 一种齿轮局部故障诊断方法及系统
CN110852201A (zh) 一种基于多脉冲包络谱匹配的脉冲信号检测方法
CN105307095B (zh) 一种基于fft的高分辨率音频频率测量方法
CN109063668B (zh) 一种基于峰值保留降采样的冲击信号包络解调方法
CN109581052A (zh) 一种迭代插值的实复转换频率估计方法
CN111912521A (zh) 一种非平稳信号的频率检测方法和存储介质
CN110347970B (zh) 分数阶同步提取广义s变换时频分解与重构方法
CN115265691B (zh) 一种科氏流量计振动频率跟踪方法及系统
CN111856401A (zh) 一种基于互谱相位拟合的时延估计方法
CN108390698A (zh) 一种基于插值fft算法的电力线载波参数测量方法
Zhu et al. Noise influence on estimation of signal parameter from the phase difference of discrete Fourier transforms
CN111948454B (zh) 一种抗衰减直流分量干扰的同步相量测量方法及系统
CN111623968B (zh) 一种基于自适应形态滤波的故障特征提取方法
CN109948487B (zh) 基于时频谱相关性分析的旋转机械故障特征提取方法
CN115014765B (zh) 一种通过声信号提取滚动轴承保持架故障特征的方法
CN116383629A (zh) 一种变转速滚动轴承故障诊断方法
CN115655719A (zh) 一种轴承振动信号分阶段降噪方法及轴承故障识别方法
CN112559973B (zh) 自适应多分量线性调频信号参数估计方法
CN108363853B (zh) 一种基于多传感器相关去噪的发动机转速测量方法
CN110779723B (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