CN113378699B - Signal processing method, device, equipment and storage medium - Google Patents
Signal processing method, device, equipment and storage medium Download PDFInfo
- Publication number
- CN113378699B CN113378699B CN202110639082.4A CN202110639082A CN113378699B CN 113378699 B CN113378699 B CN 113378699B CN 202110639082 A CN202110639082 A CN 202110639082A CN 113378699 B CN113378699 B CN 113378699B
- Authority
- CN
- China
- Prior art keywords
- signal
- vibration signal
- period
- calculating
- parameter
- 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
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000005070 sampling Methods 0.000 claims description 113
- 238000004458 analytical method Methods 0.000 claims description 49
- 238000004422 calculation algorithm Methods 0.000 claims description 44
- 238000012545 processing Methods 0.000 claims description 28
- 238000004590 computer program Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 10
- 125000004122 cyclic group Chemical group 0.000 abstract 4
- 238000010586 diagram Methods 0.000 description 18
- 230000000737 periodic effect Effects 0.000 description 18
- 230000008569 process Effects 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 239000003638 chemical reducing agent Substances 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/12—Classification; Matching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本申请实施例提供一种信号处理方法、装置、设备及存储介质。该方法包括：获取目标设备的振动信号；根据振动信号，确定振动信号的持续时长，以及振动信号的周期，其中，持续时长包括振动信号的至少两个周期；根据振动信号，计算持续时长内每个周期对应的信号有效值；根据每个周期对应的信号有效值，计算振动信号的周期平稳参数。根据本申请实施例，可以自动计算振动信号的周期平稳参数，定量分析振动信号的周期平稳性，进而提高周期平稳性的确定效率和准确性。
The embodiments of the present application provide a signal processing method, apparatus, device and storage medium. The method includes: obtaining a vibration signal of a target device; determining the duration of the vibration signal and the period of the vibration signal according to the vibration signal, wherein the duration includes at least two periods of the vibration signal; calculating the effective value of the signal corresponding to each period within the duration according to the vibration signal; and calculating the cyclic stationary parameters of the vibration signal according to the effective value of the signal corresponding to each period. According to the embodiments of the present application, the cyclic stationary parameters of the vibration signal can be automatically calculated, and the cyclic stationary nature of the vibration signal can be quantitatively analyzed, thereby improving the efficiency and accuracy of determining the cyclic stationary nature.
Description
Technical Field
The present application relates to the field of fault diagnosis technologies, and in particular, to a signal processing method, apparatus, device, and storage medium.
Background
The fault diagnosis technology is a technology for detecting abnormal conditions of equipment by monitoring state parameters of the equipment, analyzing and diagnosing fault reasons after the abnormal conditions are detected, and aims to find potential faults of the equipment so as to prevent equipment accidents from happening.
At present, vibration signals of equipment are usually collected to analyze fault conditions, but a plurality of signal-based analysis methods exist, each method has an application range and needs to be selected according to the periodic stability of the signals, so the periodic stability analysis of the signals is very important for fault diagnosis. However, the periodic stability of the signal is usually determined according to experience of a device manager, and the efficiency and the accuracy are low.
Disclosure of Invention
The embodiment of the application provides a signal processing method, a device, equipment and a storage medium, which can automatically calculate the cyclostationary parameter of a vibration signal and quantitatively analyze the cyclostationary of the vibration signal.
In a first aspect, an embodiment of the present application provides a signal processing method, including:
Acquiring a vibration signal of a target device;
determining the duration of the vibration signal and the period of the vibration signal according to the vibration signal, wherein the duration comprises at least two periods of the vibration signal;
calculating a signal effective value corresponding to each period in the duration according to the vibration signal;
And calculating the cyclostationary parameter of the vibration signal according to the signal effective value corresponding to each period.
In some implementations of the first aspect, the vibration signal is collected by a sampling device, the vibration signal including a plurality of sampling points, determining a period of the vibration signal from the vibration signal, comprising:
determining a plurality of period sampling points according to the sampling points of the vibration signals;
calculating a signal autocorrelation coefficient corresponding to each period sampling point according to the instantaneous value of the sampling point of the vibration signal;
Under the condition that at least two signal autocorrelation coefficients meet the preset proportion condition, calculating the period of the vibration signal according to the sampling frequency of the sampling device and the period sampling point number corresponding to the smallest autocorrelation parameter in the at least two signal autocorrelation coefficients.
In some implementations of the first aspect, calculating the period of the vibration signal according to a sampling frequency of the sampling device and a period sampling point number corresponding to a smallest autocorrelation parameter of the at least two signal autocorrelation coefficients includes:
And calculating the quotient of the periodic sampling point number corresponding to the minimum autocorrelation parameter and the sampling frequency, and taking the quotient as the period of the vibration signal.
In some implementations of the first aspect, calculating a signal effective value corresponding to each period in the duration from the vibration signal includes:
and calculating the signal effective value corresponding to each period in the duration according to the instantaneous value of the sampling point of the vibration signal in each period in the duration.
In some implementations of the first aspect, calculating the cyclostationary parameter of the vibration signal according to the signal effective value corresponding to each period includes:
calculating the average value of the signal effective values according to the signal effective values corresponding to each period;
And calculating variance according to the effective value and the average value of the signal corresponding to each period, and taking the variance as a cyclostationary parameter of the vibration signal.
In some implementations of the first aspect, the method further includes:
determining an analysis algorithm corresponding to the vibration signal according to the cyclostationary parameter;
the vibration signal is analyzed according to an analysis algorithm.
In some implementations of the first aspect, determining an analysis algorithm corresponding to the vibration signal according to the cyclostationary parameter includes:
under the condition that the cyclostationary parameter is larger than or equal to a preset parameter threshold value, determining a frequency domain analysis algorithm as an analysis algorithm corresponding to the vibration signal;
and under the condition that the cyclostationary parameter is smaller than a preset parameter threshold, determining a time domain analysis algorithm and/or a time-frequency domain analysis algorithm as an analysis algorithm corresponding to the vibration signal.
In a second aspect, an embodiment of the present application provides a signal processing apparatus, including:
the acquisition module is used for acquiring a vibration signal of the target equipment;
the determining module is used for determining duration time of the vibration signal and period of the vibration signal according to the vibration signal, wherein the duration time comprises at least two periods of the vibration signal;
the calculating module is used for calculating the signal effective value corresponding to each period in the duration according to the vibration signal;
And the calculation module is also used for calculating the cyclostationary parameter of the vibration signal according to the signal effective value corresponding to each period.
In some implementations of the second aspect, the vibration signal is collected by a sampling device, the vibration signal includes a plurality of sampling points, and the determining module includes:
a determining unit, configured to determine a plurality of period sampling points according to the sampling points of the vibration signal;
the first calculation unit is used for calculating a signal autocorrelation coefficient corresponding to the sampling point of each period according to the instantaneous value of the sampling point of the vibration signal;
The first calculating unit is further configured to calculate a period of the vibration signal according to a sampling frequency of the sampling device and a period sampling point number corresponding to a smallest autocorrelation parameter of the at least two signal autocorrelation coefficients when the at least two signal autocorrelation coefficients satisfy a preset ratio condition.
In some implementations of the second aspect, the first computing unit is specifically configured to:
And calculating the quotient of the periodic sampling point number corresponding to the minimum autocorrelation parameter and the sampling frequency, and taking the quotient as the period of the vibration signal.
In some implementations of the second aspect, the computing module includes:
and the second calculating unit is used for calculating the signal effective value corresponding to each period in the duration according to the instantaneous value of the sampling point of the vibration signal in each period in the duration.
In some implementations of the second aspect, the computing module includes:
a third calculating unit, configured to calculate an average value of the signal effective values according to the signal effective values corresponding to each period;
The third calculation unit is further configured to calculate a variance according to the signal effective value and the average value corresponding to each period, and use the variance as a cyclostationary parameter of the vibration signal.
In some implementations of the second aspect, the determining module is further configured to determine an analysis algorithm corresponding to the vibration signal according to the cyclostationary parameter;
the device also comprises an analysis module for analyzing the vibration signal according to an analysis algorithm.
In some implementations of the second aspect, the determining module includes:
A fourth determining unit, configured to determine that the frequency domain analysis algorithm is an analysis algorithm corresponding to the vibration signal when the cyclostationary parameter is greater than or equal to a preset parameter threshold;
And the fourth determining unit is further used for determining a time domain analysis algorithm and/or a time-frequency domain analysis algorithm as an analysis algorithm corresponding to the vibration signal under the condition that the cyclostationary parameter is smaller than a preset parameter threshold.
In a third aspect, an embodiment of the present application provides a signal processing apparatus, where the apparatus includes a processor and a memory storing computer program instructions, and where the processor implements the signal processing method described in the first aspect or any of the realizations of the first aspect when executing the computer program instructions.
In a fourth aspect, embodiments of the present application provide a computer readable storage medium having stored thereon computer program instructions which, when executed by a processor, implement the signal processing method described in the first aspect or any of the realizable modes of the first aspect.
The signal processing method, the device, the equipment and the storage medium provided by the embodiment of the application can firstly determine the duration of the vibration signal and the period of the vibration signal according to the vibration signal of the target equipment, then calculate the signal effective value corresponding to each period in the duration according to the vibration signal, and then calculate the cyclostationary parameter of the vibration signal according to the signal effective value corresponding to each period, thereby quantitatively analyzing the cyclostationary property of the vibration signal through the cyclostationary parameter and improving the determination efficiency and accuracy of the cyclostationary property.
Drawings
In order to more clearly illustrate the technical solution of the embodiments of the present application, the drawings that are needed to be used in the embodiments of the present application will be briefly described, and other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic diagram of a signal processing system according to an embodiment of the present application;
fig. 2 is a schematic flow chart of a signal processing method according to an embodiment of the present application;
FIG. 3 is a schematic diagram of a vibration signal provided by an embodiment of the present application;
FIG. 4 is a schematic diagram of a signal autocorrelation coefficient according to an embodiment of the present application;
FIG. 5 is a schematic diagram of signal effective values according to an embodiment of the present application;
FIG. 6 is a schematic diagram of another vibration signal provided by an embodiment of the present application;
FIG. 7 is a schematic diagram of another signal autocorrelation coefficients according to an embodiment of the present application;
FIG. 8 is a schematic diagram of another signal effective value provided by an embodiment of the present application;
fig. 9 is a schematic structural diagram of a signal processing device according to an embodiment of the present application;
Fig. 10 is a schematic structural diagram of a signal processing apparatus according to an embodiment of the present application.
Detailed Description
Features and exemplary embodiments of various aspects of the present application will be described in detail below, and in order to make the objects, technical solutions and advantages of the present application more apparent, the present application will be described in further detail below with reference to the accompanying drawings and examples. It should be understood that the detailed description and specific examples, while indicating the application, are intended for purposes of illustration only and are not intended to limit the scope of the application. It will be apparent to one skilled in the art that the present application may be practiced without some of these specific details. The following description of the embodiments is merely intended to provide a better understanding of the application by showing examples of the application.
It is noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising" does not exclude the presence of additional identical elements in a process, method, article, or apparatus that comprises the element.
Aiming at the problems in the background art, the embodiment of the application provides a signal processing method, a device, equipment and a storage medium, which can firstly determine the duration of a vibration signal and the period of the vibration signal according to the vibration signal of target equipment, then calculate the signal effective value corresponding to each period in the duration according to the vibration signal, and then calculate the cyclostationary parameter of the vibration signal according to the signal effective value corresponding to each period, thereby quantitatively analyzing the cyclostationary property of the vibration signal through the cyclostationary parameter and improving the determination efficiency and accuracy of the cyclostationary property.
The signal processing method, device, equipment and storage medium provided by the embodiment of the application are described in detail below through specific embodiments and application scenes thereof with reference to the accompanying drawings.
Fig. 1 is a schematic architecture diagram of a signal processing system according to an embodiment of the present application, as shown in fig. 1, where the signal processing system may include an electronic device 110 and a target device 120, and the embodiment of the present application is not limited herein.
The electronic device 110 may be a mobile electronic device or a non-mobile electronic device. For example, the mobile electronic device may be a cell phone, tablet computer, notebook computer, palm top computer or Ultra mobile personal computer (Ultra-Mobile Personal Computer, UMPC) or the like, and the non-mobile electronic device may be a server, network attached storage (Network Attached Storage, NAS) or personal computer (Personal Computer, PC) or the like. The target device 120 may be a rotating device such as a speed reducer, a generator, or a motor.
As one example, the signal processing system may be applied to equipment failure diagnosis scenarios in the fields of energy, electricity, traffic, manufacturing, and the like. As shown in fig. 1, the electronic device 110 may acquire a vibration signal, such as a voltage signal or a current signal, of the target device 120, then determine a duration of the vibration signal and a period of the vibration signal according to the vibration signal, then calculate an effective value of a signal corresponding to each period in the duration, that is, an effective value of an instantaneous value of the signal in each period, according to the vibration signal, and calculate a cyclostationary parameter of the vibration signal according to the effective value of the signal corresponding to each period, for characterizing cyclostationary performance of the vibration signal.
The signal processing method provided by the embodiment of the application will be described below. The main execution body of the signal processing method may be the electronic device 110 in the signal processing system shown in fig. 1.
Fig. 2 is a flow chart of a signal processing method according to an embodiment of the present application, as shown in fig. 2, the signal processing method may include the following steps:
s210, acquiring a vibration signal of the target device.
Specifically, in a state where the target device is operating, the vibration signal of the target device acquired by the sampling device may be acquired. The sampling device may be a sensor, and the vibration signal may be a voltage signal or a current signal, etc., including a plurality of sampling points.
S220, determining the duration of the vibration signal and the period of the vibration signal according to the vibration signal.
In one embodiment, the duration of the vibration signal may be calculated from the number of sampling points of the vibration signal and the sampling frequency of the sampling device. For example, the quotient obtained by dividing the sampling point number by the sampling frequency is the duration.
Meanwhile, a plurality of period sampling points can be determined according to the sampling points of the vibration signal. Wherein the periodic sampling points represent the number of sampling points in one candidate period, and the plurality of periodic sampling points represent the number of sampling points in a plurality of different candidate periods.
And calculating a signal autocorrelation coefficient corresponding to the sampling point of each period according to the instantaneous value of the sampling point of the vibration signal. Wherein the signal autocorrelation coefficients represent the degree of autocorrelation of the vibration signal.
In case at least two signal autocorrelation coefficients meet a preset scaling condition, for example, at least two signal autocorrelation coefficients meet an increasing scaling. The period of the vibration signal can be rapidly and accurately calculated according to the sampling frequency of the sampling device and the period sampling point number corresponding to the smallest autocorrelation parameter in the at least two signal autocorrelation coefficients. For example, a quotient of the number of periodic sampling points corresponding to the smallest autocorrelation parameter divided by the sampling frequency may be calculated as the period of the vibration signal.
As an example, the sampling frequency of the sampling device is m, the sampling length of the vibration signal is n, i.e. the vibration signal comprises n equally spaced sampling points, the instantaneous value series defining the sampling points is x (n), n=0, 1, 2. The number of periodic sampling points is k=1, 2. The signal autocorrelation coefficient corresponding to each period sampling point can be calculated according to the instantaneous value of the sampling point and the signal autocorrelation formula. Illustratively, the signal autocorrelation formula may be as follows:
If the ratio between the autocorrelation coefficients of the signals corresponding to k=k1, k2, k3, and k4 is 1:2:3:4, respectively, the quotient of k1 divided by m, i.e., m/k1, can be calculated, and the period of the vibration signal is m/k1.
S230, calculating a signal effective value corresponding to each period in the duration according to the vibration signal.
In one embodiment, the effective value of the signal corresponding to each period in the duration, that is, the effective value of the instantaneous value of the signal in each period, can be quickly and accurately calculated according to the instantaneous value of the sampling point of the vibration signal in each period in the duration.
As one example, the signal effective value for each period may be calculated from the instantaneous value of the sampling point and the signal effective value formula for each period. Illustratively, the signal effective value formula may be as follows:
Wherein Rms (T) represents the signal effective value of the T-th period in the duration, N 1 represents the sequence number of the last sampling point in the T-th period in the duration, N 0 represents the sequence number of the first sampling point in the T-th period in the duration, N 1-N0 represents the number of sampling points in the T-th period in the duration, and T represents the number of periods in the duration.
S240, calculating the cyclostationary parameter of the vibration signal according to the signal effective value corresponding to each period.
In one embodiment, an average value of the signal effective values may be calculated according to the signal effective values corresponding to each period, and a variance may be calculated according to the signal effective values and the average value corresponding to each period, with the variance being used as a cyclostationary parameter of the vibration signal. Thus, the difference of sampling points in each period can be accurately represented through variance, and the period stability of the vibration signal can be intuitively reflected.
As one example, the variance may be calculated from the signal effective value and the average value corresponding to each period, and a variance formula. Illustratively, the variance formula may be as follows:
Where Stability represents a cyclostationary parameter, rms (T) represents the signal effective value of the T-th period in the duration, and μ represents the average value of the signal effective values of the T periods.
In the embodiment of the application, the duration of the vibration signal and the period of the vibration signal can be determined according to the vibration signal of the target equipment, then the signal effective value corresponding to each period in the duration is calculated according to the vibration signal, and then the cyclostationary parameter of the vibration signal is calculated according to the signal effective value corresponding to each period, so that the cyclostationary property of the vibration signal is quantitatively analyzed through the cyclostationary parameter, and the determination efficiency and accuracy of the cyclostationary property are improved.
In one embodiment, an analysis algorithm corresponding to the vibration signal may be determined based on the cyclostationary parameter, and the vibration signal may be analyzed based on the analysis algorithm. The frequency domain analysis algorithm may be determined to be an analysis algorithm corresponding to the vibration signal when the cyclostationary parameter is greater than or equal to a preset parameter threshold, and the time domain analysis algorithm and/or the time-frequency domain analysis algorithm may be determined to be an analysis algorithm corresponding to the vibration signal when the cyclostationary parameter is less than the preset parameter threshold. The frequency domain analysis algorithm can be a cepstrum analysis algorithm or a refined spectrum analysis algorithm, the time domain analysis algorithm can be a classical time domain analysis method, and the time-frequency domain analysis algorithm can be a wavelet transformation analysis algorithm. Therefore, an analysis algorithm suitable for the vibration signal can be determined based on the cyclostationary parameter, and the analysis effect is improved.
The following describes a signal processing method provided by the embodiment of the present application with reference to a simulation example, which is specifically as follows:
as shown in fig. 3, the vibration signal is x (n) =sin (0.125 n), the sampling frequency is 4, the period of the vibration signal is 4s, and the number of sampling points in the period is 16.
And (3) simulating the vibration signal, determining the duration of the vibration signal and a plurality of periodic sampling points, and calculating a signal autocorrelation coefficient corresponding to each periodic sampling point according to the instantaneous value of the sampling point and the formula (1). The autocorrelation coefficients of the signal corresponding to the sampling points of each period may be as shown in fig. 4. In fig. 4, the signal autocorrelation coefficient is 8 when the periodic sampling point is 16, the signal autocorrelation coefficient is 16 when the periodic sampling point is 32, the signal autocorrelation coefficient is 24 when the periodic sampling point is 48, the signal autocorrelation coefficient is 32 when the periodic sampling point is 64, and the ratio of the signal autocorrelation coefficients is 1:2:3:4. Therefore, it is possible to determine 16, which is the minimum number of period samples, as the target number of period samples of the vibration signal, and 4, which is the period of the vibration signal, obtained by dividing 16 by 4, is the period of the vibration signal, that is, 4s.
The signal effective value for each period in the duration is then calculated from the instantaneous value of the sampling point in each period and equation (2). The signal valid value corresponding to each period may be as shown in fig. 5. It is known that the signal effective values of 4 periods of the vibration signal are equal, and the variance calculated according to the formula (3) is 0, that is, the cyclostationary parameter is 0. The effectiveness of the signal processing method provided by the embodiment of the application can be verified through simulation signals.
The signal processing method provided by the embodiment of the application is described below with reference to a specific example, which is specifically as follows:
The target equipment is a speed reducer, a speed sensor is arranged on the shell of the speed reducer, the sampling frequency of the speed sensor is 2048Hz, the collected vibration signal of the speed reducer is shown in fig. 6, and a plurality of burrs exist in the sampling signal.
Processing the vibration signal, determining the duration of the vibration signal and a plurality of periodic sampling points, and calculating a signal autocorrelation coefficient corresponding to each periodic sampling point according to the instantaneous value of the sampling point and the formula (1). The autocorrelation coefficients of the signal corresponding to the sampling points of each period may be as shown in fig. 7. In fig. 7, the signal autocorrelation coefficient is 3539 when the period sampling point is 408, 7135 when the period sampling point is 814, 10280 when the period sampling point is 1222, 12560 when the period sampling point is 1630, 16970 when the period sampling point is 2038, and the ratio of the signal autocorrelation coefficients is about 1:2:3:4:5. Therefore, it is possible to determine 408, which is the minimum number of period samples, as the target number of period samples of the vibration signal, and to divide 408 by 2048 to obtain 0.2 as the period of the vibration signal, that is, the period of the vibration signal is 0.2s.
The signal effective value for each period in the duration is then calculated from the instantaneous value of the sampling point in each period and equation (2). The signal valid value corresponding to each period may be as shown in fig. 8. It is known that the effective values of the signals of the 1 st to 5 th periods of the vibration signal are 0.1133, 0.1416, 0.2761, 0.1798, 0.2244, respectively, and the variance calculated according to the formula (3) is 0.0034, that is, the cyclostationary parameter is 0.0034. The validity of the signal processing method provided by the embodiment of the application can be verified through a specific example.
Based on the signal processing method provided by the embodiment of the present application, the embodiment of the present application further provides a signal processing apparatus, as shown in fig. 9, a signal processing apparatus 900 may include:
an acquisition module 910, configured to acquire a vibration signal of the target device.
A determining module 920, configured to determine a duration of the vibration signal and a period of the vibration signal according to the vibration signal, where the duration includes at least two periods of the vibration signal.
And the calculating module 930 is configured to calculate a signal effective value corresponding to each period in the duration according to the vibration signal.
The calculating module 930 is further configured to calculate a cyclostationary parameter of the vibration signal according to the signal effective value corresponding to each period.
In some embodiments, the vibration signal is collected by a sampling device, the vibration signal comprising a plurality of sampling points, the determining module 920 comprising:
and the determining unit is used for determining a plurality of period sampling points according to the sampling points of the vibration signal.
And the first calculation unit is used for calculating the signal autocorrelation coefficient corresponding to the sampling point of each period according to the instantaneous value of the sampling point of the vibration signal.
The first calculating unit is further configured to calculate a period of the vibration signal according to a sampling frequency of the sampling device and a period sampling point number corresponding to a smallest autocorrelation parameter of the at least two signal autocorrelation coefficients when the at least two signal autocorrelation coefficients satisfy a preset ratio condition.
In some embodiments, the first computing unit is specifically configured to:
And calculating the quotient of the periodic sampling point number corresponding to the minimum autocorrelation parameter and the sampling frequency, and taking the quotient as the period of the vibration signal.
In some embodiments, the computing module 930 includes:
and the second calculating unit is used for calculating the signal effective value corresponding to each period in the duration according to the instantaneous value of the sampling point of the vibration signal in each period in the duration.
In some embodiments, the computing module 930 includes:
And the third calculating unit is used for calculating the average value of the signal effective values according to the signal effective values corresponding to each period.
The third calculation unit is further configured to calculate a variance according to the signal effective value and the average value corresponding to each period, and use the variance as a cyclostationary parameter of the vibration signal.
In some embodiments, the determining module 920 is further configured to determine an analysis algorithm corresponding to the vibration signal according to the cyclostationary parameter.
The signal processing device 900 further comprises an analysis module for analyzing the vibration signal according to an analysis algorithm.
In some embodiments, the determination module 920 includes:
and the fourth determining unit is used for determining that the frequency domain analysis algorithm is the analysis algorithm corresponding to the vibration signal under the condition that the cyclostationary parameter is larger than or equal to the preset parameter threshold value.
And the fourth determining unit is further used for determining a time domain analysis algorithm and/or a time-frequency domain analysis algorithm as an analysis algorithm corresponding to the vibration signal under the condition that the cyclostationary parameter is smaller than a preset parameter threshold.
It can be understood that each module/unit in the signal processing apparatus 900 shown in fig. 9 has a function of implementing each step in the signal processing method provided in the embodiment of the present application, and can achieve the corresponding technical effects, which are not described herein for brevity.
Fig. 10 is a schematic structural diagram of a signal processing apparatus according to an embodiment of the present application. As shown in fig. 10, the signal processing device may include a processor 1001 and a memory 1002 storing computer program instructions.
In particular, the processor 1001 may include a central processing unit (Central Processing Unit, CPU), or an Application SPECIFIC INTEGRATED Circuit (ASIC), or may be configured as one or more integrated circuits that implement embodiments of the present application.
Memory 1002 may include mass storage for data or instructions. By way of example, and not limitation, memory 1002 may include a hard disk drive (HARD DISK DRIVE, HDD), floppy disk drive, flash memory, optical disk, magneto-optical disk, magnetic tape, or universal serial bus (Universal Serial Bus, USB) drive, or a combination of two or more of the foregoing. The memory 1002 may include removable or non-removable (or fixed) media, where appropriate. The memory 1002 may be internal or external to the signal processing device, where appropriate. In a particular embodiment, the memory 1002 is a non-volatile solid state memory. In a particular embodiment, the memory 1002 includes Read Only Memory (ROM). The ROM may be mask programmed ROM, programmable ROM (PROM), erasable PROM (EPROM), electrically Erasable PROM (EEPROM), electrically rewritable ROM (EAROM), or flash memory, or a combination of two or more of these, where appropriate.
The processor 1001 may implement the signal processing method provided by the embodiment of the present application by reading and executing the computer program instructions stored in the memory 1002, and achieve the corresponding technical effects achieved by executing the method according to the embodiment of the present application, which is not described herein for brevity.
In one example, the signal processing device may also include a communication interface 1003 and a bus 1010. As shown in fig. 10, the processor 1001, the memory 1002, and the communication interface 1003 are connected to each other by a bus 1010, and perform communication with each other.
The communication interface 1003 is mainly used for implementing communication among the modules, devices, units and/or apparatuses in the embodiment of the application.
Bus 1010 includes hardware, software, or both that couple the components of the signal processing device to each other. By way of example, and not limitation, the buses may include an accelerated graphics Port (ACCELERATED GRAPHICS Port, AGP) or other graphics Bus, an enhanced industry Standard architecture (Extended Industry Standard Architecture, EISA) Bus, a Front Side Bus (FSB), a HyperTransport (HT) interconnect, an industry Standard architecture (Industry Standard Architecture, ISA) Bus, an Infiniband interconnect, a Low Pin Count (LPC) Bus, a memory Bus, a Micro Channel Architecture (MCA) Bus, a Peripheral Component Interconnect (PCI) Bus, a PCI-Express (PCI-X) Bus, a Serial Advanced Technology Attachment (SATA) Bus, a video electronics standards Association local (VLB) Bus, or other suitable Bus, or a combination of two or more of these. Bus 1010 may include one or more buses, where appropriate. Although embodiments of the application have been described and illustrated with respect to a particular bus, the application contemplates any suitable bus or interconnect.
The signal processing device can execute the signal processing method in the embodiment of the application, thereby realizing the corresponding technical effects of the signal processing method provided by the embodiment of the application.
In addition, the embodiment of the application also provides a computer readable storage medium, wherein the computer readable storage medium is stored with computer program instructions, and the signal processing method provided by the embodiment of the application is realized when the computer program instructions are executed by a processor.
It should be clear that, all embodiments in this specification are described in a progressive manner, and the same or similar parts of all embodiments are referred to each other, so that for brevity, no further description is provided. The present application is not limited to the specific configurations and processes described above and shown in the drawings. For the sake of brevity, a detailed description of known methods is omitted here. In the above embodiments, several specific steps are described and shown as examples. The method processes of the present application are not limited to the specific steps described and shown, but various changes, modifications and additions, or the order between steps may be made by those skilled in the art after appreciating the spirit of the present application.
The functional blocks shown in the above-described structural block diagrams may be implemented in hardware, software, firmware, or a combination thereof. When implemented in hardware, it may be, for example, an electronic Circuit, application SPECIFIC INTEGRATED Circuit (ASIC), appropriate firmware, plug-in, function card, or the like. When implemented in software, the elements of the application are the programs or code segments used to perform the required tasks. The program or code segments may be stored in a machine readable medium or transmitted over transmission media or communication links by a data signal carried in a carrier wave. A "machine-readable medium" may include any medium that can store or transfer information. Examples of machine-readable media include electronic circuitry, semiconductor Memory devices, read-Only Memory (ROM), flash Memory, erasable ROM (EROM), floppy disks, CD-ROMs, optical disks, hard disks, fiber optic media, radio Frequency (RF) links, and the like. The code segments may be downloaded via computer networks such as the internet, intranets, etc.
It should also be noted that the exemplary embodiments mentioned in this disclosure describe some methods or systems based on a series of steps or devices. The present application is not limited to the order of the above-described steps, that is, the steps may be performed in the order mentioned in the embodiments, or may be performed in a different order from the order in the embodiments, or several steps may be performed simultaneously.
Aspects of the present disclosure are described above with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the disclosure. It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, enable the implementation of the functions/acts specified in the flowchart and/or block diagram block or blocks. Such a processor may be, but is not limited to being, a general purpose processor, a special purpose processor, an application specific processor, or a field programmable logic circuit. It will also be understood that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware which performs the specified functions or acts, or combinations of special purpose hardware and computer instructions.
In the foregoing, only the specific embodiments of the present application are described, and it will be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working processes of the systems, modules and units described above may refer to the corresponding processes in the foregoing method embodiments, which are not repeated herein. It should be understood that the scope of the present application is not limited thereto, and any equivalent modifications or substitutions can be easily made by those skilled in the art within the technical scope of the present application, and they should be included in the scope of the present application.
Claims (8)
1. A method of signal processing, the method comprising:
Acquiring a vibration signal of a target device;
determining the duration of the vibration signal and the period of the vibration signal according to the vibration signal, wherein the duration comprises at least two periods of the vibration signal;
determining a plurality of period sampling points according to the sampling points of the vibration signals;
calculating a signal autocorrelation coefficient corresponding to each period sampling point according to the instantaneous value of the sampling point of the vibration signal;
Calculating the period of the vibration signal according to the sampling frequency of the sampling device and the period sampling point number corresponding to the smallest autocorrelation parameter in the at least two signal autocorrelation coefficients under the condition that the at least two signal autocorrelation coefficients meet the preset proportion condition, specifically, calculating the quotient of the period sampling point number corresponding to the smallest autocorrelation parameter divided by the sampling frequency, and taking the quotient as the period of the vibration signal;
Calculating a signal effective value corresponding to each period in the duration according to the vibration signal;
and calculating the cyclostationary parameter of the vibration signal according to the signal effective value corresponding to each period.
2. The method of claim 1, wherein calculating a signal effective value corresponding to each period in the duration from the vibration signal comprises:
And calculating a signal effective value corresponding to each period in the duration according to the instantaneous value of the sampling point of the vibration signal in each period in the duration.
3. The method according to any one of claims 1-2, wherein calculating the cyclostationary parameter of the vibration signal based on the signal effective value corresponding to each period comprises:
calculating the average value of the signal effective values according to the signal effective values corresponding to each period;
And calculating variance according to the signal effective value corresponding to each period and the average value, and taking the variance as a cyclostationary parameter of the vibration signal.
4. The method according to claim 1, wherein the method further comprises:
determining an analysis algorithm corresponding to the vibration signal according to the cyclostationary parameter;
And analyzing the vibration signal according to the analysis algorithm.
5. The method of claim 4, wherein said determining an analysis algorithm corresponding to said vibration signal based on said cyclostationary parameter comprises:
under the condition that the cyclostationary parameter is larger than or equal to a preset parameter threshold, determining a frequency domain analysis algorithm as an analysis algorithm corresponding to the vibration signal;
And under the condition that the cyclostationary parameter is smaller than a preset parameter threshold, determining a time domain analysis algorithm and/or a time-frequency domain analysis algorithm as an analysis algorithm corresponding to the vibration signal.
6. A signal processing apparatus, the apparatus comprising:
the acquisition module is used for acquiring a vibration signal of the target equipment;
The system comprises a determining module, a sampling point counting module and a processing module, wherein the determining module is used for determining the duration of the vibration signal and the period of the vibration signal according to the vibration signal, wherein the duration comprises at least two periods of the vibration signal;
The system comprises a sampling device, a calculation module, a signal self-correlation coefficient calculation module, a signal effective value calculation module and a signal effective value calculation module, wherein the sampling device is used for calculating a sampling point of a vibration signal according to an instantaneous value of the sampling point of the vibration signal, the sampling point corresponds to the signal self-correlation coefficient of each period;
the calculation module is further used for calculating the cyclostationary parameter of the vibration signal according to the signal effective value corresponding to each period.
7. A signal processing device comprising a processor and a memory storing computer program instructions, the processor implementing the signal processing method according to any one of claims 1-5 when executing the computer program instructions.
8. A computer readable storage medium, characterized in that the computer readable storage medium has stored thereon computer program instructions, which when executed by a processor, implement the signal processing method according to any of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110639082.4A CN113378699B (en) | 2021-06-08 | 2021-06-08 | Signal processing method, device, equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110639082.4A CN113378699B (en) | 2021-06-08 | 2021-06-08 | Signal processing method, device, equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113378699A CN113378699A (en) | 2021-09-10 |
| CN113378699B true CN113378699B (en) | 2024-12-06 |
Family
ID=77572859
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110639082.4A Active CN113378699B (en) | 2021-06-08 | 2021-06-08 | Signal processing method, device, equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113378699B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107976304A (en) * | 2016-10-21 | 2018-05-01 | 计算系统有限公司 | The mechanical disorder prediction analyzed based on the periodical information to signal |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010271185A (en) * | 2009-05-21 | 2010-12-02 | Mitsubishi Electric Corp | Device for monitoring abnormal condition |
| CN101749256B (en) * | 2010-01-08 | 2013-11-13 | 浙江大学 | Large axial flow fan unbalance recognition method based on auto-correlation |
| CN101936812B (en) * | 2010-08-25 | 2011-10-05 | 重庆大学 | Detection method of armored vehicle engine working state |
| US9791422B2 (en) * | 2013-07-02 | 2017-10-17 | Computational Systems, Inc. | Analysis of periodic information in a signal |
| CN105865793B (en) * | 2016-03-25 | 2018-04-20 | 西北工业大学 | A kind of method for improving more rotor vibration monitoring of aero-engine precision |
| JP6768489B2 (en) * | 2016-12-22 | 2020-10-14 | 三菱電機株式会社 | Signal processing device and signal processing method |
| CN107461342B (en) * | 2017-08-16 | 2019-05-31 | 江苏大学 | A kind of centrifugal pump fault inline diagnosis method and system |
| CN111947957A (en) * | 2020-08-17 | 2020-11-17 | 浙江中控技术股份有限公司 | Vibration signal sampling method and device, storage medium and electronic equipment |
-
2021
- 2021-06-08 CN CN202110639082.4A patent/CN113378699B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107976304A (en) * | 2016-10-21 | 2018-05-01 | 计算系统有限公司 | The mechanical disorder prediction analyzed based on the periodical information to signal |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113378699A (en) | 2021-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108896973B (en) | Radar data calibration method, trace point quality evaluation method and storage medium | |
| CN106850511B (en) | Method and device for identifying access attack | |
| KR101426863B1 (en) | A method for recognizing radar intra-pulse modulation type using features | |
| CN113378699B (en) | Signal processing method, device, equipment and storage medium | |
| KR101348635B1 (en) | Diagnosis apparatus and methof for broken rotor bar in induction motors | |
| CN117872166B (en) | Method and device for detecting thermal runaway of energy storage battery and electronic equipment | |
| CN116904569B (en) | Signal processing method, device, electronic equipment, medium and product | |
| CN106936561B (en) | Side channel attack protection capability assessment method and system | |
| CN108596215B (en) | Multi-modal signal analysis and separation method, device, equipment and storage medium | |
| CN117520887A (en) | Methods, devices, equipment and storage media for determining components at risk of hydraulic turbine failure | |
| CN116628481A (en) | Electronic countermeasure information source identification method, system, computer and readable storage medium | |
| CN116775409A (en) | Power consumption evaluation method, device, electronic equipment and computer readable storage medium | |
| CN112566013B (en) | Target equipment positioning method, device, equipment and computer storage medium | |
| CN105959035B (en) | A kind of direct sequence signal intercepts and captures detection method | |
| CN114092219A (en) | Model verification method, apparatus, electronic device and storage medium | |
| CN109870621B (en) | Operation judgment method and device | |
| JPWO2018142620A1 (en) | Packet format estimation device and packet format estimation program | |
| CN114035996A (en) | Problem locating method, apparatus, device and readable storage medium | |
| CN115575979A (en) | Integrity risk probability determination method and device and electronic equipment | |
| CN114547133B (en) | Multi-dimensional dataset-based conversational attribution analysis method, device and equipment | |
| CN117007083B (en) | Vehicle ranging capability assessment method, system, equipment and medium thereof | |
| JPH1075218A (en) | Method and device for identifying characteristic boundary of data groups | |
| Guo et al. | Adaptive signal decomposition methods for vibration signals of rotating machinery | |
| CN109472051B (en) | Hardware Trojan horse detection method based on short-time Fourier transform | |
| CN119621389A (en) | Frozen frame generation method and related device |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20240927 Address after: Room 506, Building 1.0B, 5F Entrepreneurship Park, No. 118 Science Avenue, High tech Zone, Hefei City, Anhui Province, China 230088 Applicant after: Hefei Xinli Technology Co.,Ltd. Country or region after: China Address before: 518054 14 / F, Weisheng technology building, 9966 Shennan Avenue, Maling community, Yuehai street, Nanshan District, Shenzhen City, Guangdong Province Applicant before: Xilenma (Shenzhen) Technology Co.,Ltd. Country or region before: China |
|
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20241104 Address after: Room 7-8020, Aokai City Plaza, No. 1777 Zhonghui Avenue, Qianzhou Street, Huishan District, Jiangsu Province (Huishan Station area of the urban railway) Applicant after: Wuxi Wanguo Tianyi Technology Co.,Ltd. Country or region after: China Address before: Room 506, Building 1.0B, 5F Entrepreneurship Park, No. 118 Science Avenue, High tech Zone, Hefei City, Anhui Province, China 230088 Applicant before: Hefei Xinli Technology Co.,Ltd. Country or region before: China |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20250317 Address after: 7th Floor, 688 Zhenghe Avenue, Huishan Economic Development Zone, Wuxi City, Jiangsu Province, China 214000 Patentee after: Wuxi Lingsheng Technology Co.,Ltd. Country or region after: China Address before: Room 7-8020, Aokai City Plaza, No. 1777 Zhonghui Avenue, Qianzhou Street, Huishan District, Jiangsu Province (Huishan Station area of the urban railway) Patentee before: Wuxi Wanguo Tianyi Technology Co.,Ltd. Country or region before: China |