CN102944360A - Imbalanced signal extracting method of dual-rotor system with slight speed difference - Google Patents
Imbalanced signal extracting method of dual-rotor system with slight speed difference Download PDFInfo
- Publication number
- CN102944360A CN102944360A CN2012104076284A CN201210407628A CN102944360A CN 102944360 A CN102944360 A CN 102944360A CN 2012104076284 A CN2012104076284 A CN 2012104076284A CN 201210407628 A CN201210407628 A CN 201210407628A CN 102944360 A CN102944360 A CN 102944360A
- Authority
- CN
- China
- Prior art keywords
- cos
- centerdot
- omega
- sin
- lambda
- 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
Links
Images
Landscapes
- Testing Of Balance (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention provides an imbalanced signal extracting method of a dual-rotor system with slight speed difference, and aims to solve the problem of difficultly in extracting imbalanced signals of the system in the prior art. With regard to un-damped real-value sinusoid signals, an improved Prony method (IPM for short) can be adopted to rather precisely estimate the frequency of each sinusoid signal. For the characteristics of beat vibration signals in the dual-rotor system with slight speed difference4, a band-pass tracking filtering method is adopted to reduce noise and then samples are extracted from filtered data, the IPM is adopted to realize frequency identification of two adjacent frequency signals; and then a least square method is used for identifying the amplitude and the phase of the beat vibration signals.
Description
Technical field
The invention belongs to the mechanical measurement technique field, particularly a kind of unbalanced signal extracting method of Dual Rotor System with Little Rotation Speed Difference.
Background technology
Dual Rotor System with Little Rotation Speed Difference is widely used in the industries such as oil, chemical industry, pharmacy, two synthetic signals that shake of clapping of power frequency component that frequency is close in the vibration signal of this type systematic, and this extracts to unbalanced signal and brings larger difficulty.
Unbalanced signal extracting method for this type of rotor-support-foundation system has two classes at present, one class is not understand shooting method, clap the peak and clap paddy from time domain waveform analysis, then according to the bat principle of shaking, amplitude through simply calculating two rotor oscillations and phase place, the method needs a bat above Measuring Time of cycle of shaking at least, and needs could determine by complicated algorithm which rotor two amplitudes respectively belong to, and the imbalance of one-shot measurement goes heavily rate to only have about 70%.Equations of The Second Kind is to separate shooting method, isolates two power frequency components by signal processing method, thereby extracts separately amplitude and phase place.Traditional solution shooting method generally adopts relevant function method, and the method need to be done digital integration in the cycle is shaken in whole bat could obtain degree of precision, the same shortcoming that needs Measuring Time long that exists, and accurately obtain to clap also difficult of the cycle of shaking..
It is poor that above-mentioned two methods all need to measure accurately the speed of two rotors, less fast difference measurements error will produce larger impact to the signal extraction precision, and internal rotor generally is enclosed in the outer rotor and speed probe can't be installed in the real system, adds the small speed of actual embedded measurement system Measurement accuracy poor also being difficult for and realizes.Although adopt fft analysis can obtain signal frequency from vibration signal, the data sample in needing for a long time just can obtain higher frequency resolution, and the method real-time is too poor and be not suitable for engineering and use.
Summary of the invention
The present invention proposes a kind of unbalanced signal extracting method of Dual Rotor System with Little Rotation Speed Difference, the problem that is difficult to extract with the unbalanced signal that solves in the prior art for this system.
The Prony method can be used than the short data sample and realize the high precision spectrum estimation, be widely used in Fault Signal Analyses in HV Transmission identification, for the real-valued sinusoidal signal of non-decay, adopt improved Prony method (Improved Prony Method is called for short IPM) can estimate more accurately the frequency of each sinusoidal signal.The present invention is directed to bat in the poor birotor of the dead slow speed signal characteristic that shakes, adopt first the logical tracking filter method noise reduction of band, then sample drawn from filtered data, adopt the IPM method to realize the frequency identification of two near by frequency signals, adopt at last least square method to identify amplitude and the phase place of clapping the signal that shakes.
Technical scheme of the present invention is that a kind of unbalanced signal extracting method of Dual Rotor System with Little Rotation Speed Difference may further comprise the steps:
A1 makes that the sampling period is T, and the vibration signal take T as the cycle to described Dual Rotor System with Little Rotation Speed Difference carries out equal interval sampling, obtains one group of data sequence { D (kT) }, k=1, and 2 ... K, K are a positive integer;
A2, getting centre frequency is f
0, quality factor q>2 adopt the digital band pass tracking filter that { D (kT) } carried out filtering, obtain { y (kT) }, k=1, and 2 ... K, wherein, the pulsed transfer function of digital band pass tracking filter is
A3 with a default rule, makes data to { y (kT) } and extracts, the acquisition data sequence x[n] }, n=0 in the formula, 1 ... N-1 is time series, and N is a positive integer;
A4 establishes { x[n] } and contains consisting of of noisy sinusoidal signal
x[n]=a
1cos(ω
1n+φ
1)+a
2cos(ω
2n+φ
2)+e(n)(1)
N=0 in the formula, 1 ... N-1 is time series, a
1, a
2Be unbalanced signal amplitude, ω
1, ω
2Be unbalanced signal frequency and φ
1, φ
2Be the unbalanced signal phase place,
The capable difference equation matrix of structure N-4 is
X[n]=[x[n]+x[n+4],x[n+1]+x[n+3],x[n+2]](2)
The Prony method shows, the homogeneous solution of above-mentioned equation has been specified e
I ωIn 4 rank multinomial coefficients, because sampled data all is real number, e
I ωIn 4 rank polynomial expressions can be rewritten into formula about cos (2 ω), thereby can solve with the Chebyshev polynomial expression.
By matrix X structure variance matrix G,
M=N-4 in the formula, all summations are limited to up and down from 0 to M-1, and the minimal eigenvalue of variance matrix G is by its secular equation polynomial expression
c
1λ
3+c
2λ
2+c
3λ+c
4=0(4)
Coefficient ask for, this coefficient c is
Then minimal eigenvalue is
Wherein S, T computing formula are,
Minimal eigenvalue λ
0The characteristic of correspondence vector is
Obtain the minimal eigenvalue λ of G matrix by matrix decomposition
0With its characteristic of correspondence vector
Then the total least square solution of matrix X is
Known cos (2 θ)=2 (cos (θ)) by the Chebyshev polynomial expression
2-1(10)
According to the Prony Method And Principle, the frequency of any in homogeneous solution v and 2 unknown signalings all satisfies equation
2cos(2ω
k)v
1+2cos(ω
k)v
2+v
3=0(11)
Then the frequency estimation of unknown signaling is
Obtain the frequencies omega of unbalanced signal in the Dual Rotor System with Little Rotation Speed Difference vibration signal
1And ω
2
A5 utilizes least square method to estimate ω
1And ω
2Amplitude and phase place, formula (1) can be write as
x[n]=a
1cos(ω
1n)cos(φ
1)-a
1sin(ω
1n)sin(φ
1)
+a
2cos(ω
2n)cos(φ
2)-a
2sin(ω
2n)sin(φ
2)+e(n)(13)
If y
1=a
1Cos (φ
1), y
2=a
1Sin (φ
1), y
3=a
2Cos (φ
2), y
4=a
2Sin (φ
2), treat that estimated parameter is y
1, y
2, y
3, y
4, according to each time observed reading x[n] (n=1,2 ... N), can draw the system of linear equations (14) of Ay=x, ask its least square solution can pick out each parameter
The frequencies omega of unbalanced signal then
1And ω
2The estimated value of amplitude and phase place is
φ
1=arctan(y
2/y
1),φ
2=arctan(y
4/y
3)(15)。
The present invention is directed to bat in the poor birotor of the dead slow speed signal characteristic that shakes, adopt first the logical tracking filter method noise reduction of band, then sample drawn from filtered data, adopt the IPM method to realize the frequency identification of two near by frequency signals, at last before the filtering data as sample (because of wave filter influential to unbalanced signal), adopt amplitude and the phase place of two unbalanced signals of least square method identification.
Description of drawings
Fig. 1 is that the present invention is used for the process flow diagram that faint unbalanced signal extracts.
Embodiment
The present invention is a kind of unbalanced signal extracting method for Dual Rotor System with Little Rotation Speed Difference, comprises following content:
Digital band pass tracking filter link, the pulsed transfer function of wave filter is
Each coefficient is as follows in the formula:
Wherein Q is the wave filter quality factor, f
0Centered by frequency, T is the sampling period.
This wave filter is unattenuated in centre frequency place amplitude, without phase shift, but the Q value is when larger, near the phase shift variations the centre frequency is violent, and transient process is long.For this reason, the Q value should not be selected excessive, is more than 2 times of inner and outer rotors difference on the frequency but generally should make filter passbands wide.
Parameter identification link based on two signals of IPM:
If contain consisting of of a noisy K sinusoidal signal
x[n]=a
1cos(ω
1n+φ
1)+a
2cos(ω
2n+φ
2)+e(n)(1)
N=0 in the formula, 1 ... N-1 is time series, amplitude a
1, a
2, frequencies omega
1, ω
2And phase
1, φ
2Be unknown parameter.
The capable difference equation matrix of structure N-4 is as follows:
X[n]=[x[n]+x[n+4],x[n+1]+x[n+3],x[n+2]](2)
Known by the Prony Method And Principle, the homogeneous solution of above-mentioned DIFFERENCE EQUATIONS is e
I ωIn the polynomial coefficient in 4 rank, when sampled data is real number, about e
I ω4 rank polynomial expressions can use cos (2 ω) to substitute, thereby can utilize the Chebyshev polynomial solving.
For asking the total least square solution of equation (2), first by matrix X structure modified variance matrix G,
M=N-4 in the formula, all summations are limited to from 0 to M-1 up and down.Because the G matrix is symmetrical matrix, only need the upper triangle of calculating or lower triangle to get final product.
The minimal eigenvalue of G can pass through its secular equation polynomial expression
c
1λ
3+c
2λ
2+c
3λ+c
4=0(4)
Coefficient ask for.Coefficient c is
Then minimal eigenvalue is
Wherein S, T computing formula are
λ
0The characteristic of correspondence vector is
Obtain the minimal eigenvalue λ of G matrix by matrix decomposition
0With its characteristic of correspondence vector
Then the total least square solution of X is
Known by the Chebyshev polynomial expression,
cos(2ω)=2(cos(ω))
2-1(10)
The frequency of any in homogeneous solution v and 2 unknown signalings all satisfies equation
2cos(2ω
k)v
1+2cos(ω
k)v
2+v
3=0(11)
Then the frequency estimation of unknown signaling is
After signal frequency estimates, utilize least square method (LSM) can estimate amplitude and phase place.Formula (1) can be write as
x[n]=a
1cos(ω
1n)cos(φ
1)-a
1sin(ω
1n)sin(φ
1)
+a
2cos(ω
2n)cos(φ
2)-a
2sin(ω
2n)sin(φ
2)+e(n)(13)
If y
1=a
1Cos (φ
1), y
2=a
1Sin (φ
1), y
3=a
2Cos (φ
2), y
4=a
2Sin (φ
2), treat that estimated parameter is y
1, y
2, y
3, y
4, according to each time observed reading x[n] (n=1,2 ... N), can draw the system of linear equations (14) of Ay=x, ask its least square solution can pick out each parameter
Then the estimated value of amplitude and phase place is
φ
1=arctan(y
2/y
1),φ
2=arctan(y
4/y
3) (15)
The below extracts with certain Dual Rotor System with Little Rotation Speed Difference (the inner and outer rotors frequency is respectively 28.283Hz and 28.1Hz) unbalanced signal specific implementation method of the present invention is described as an example.
Process flow diagram divided for 5 steps finished as shown in Figure 1, and the concrete operations in per step are as follows:
1, getting the sampling period is 1000Hz, and equal interval sampling obtains 1000 data;
2, get quality factor q=20, centre frequency is 28.1Hz, adopts digital band pass tracking filter algorithm to carry out filtering operation;
3, from the 401st of data after the filtering, per 16 data extract one, extract altogether 20 data as sample;
4, with formula (2)~(12) in the data substitution above-mentioned IP M method, obtain the frequencies omega of two unbalanced signals
1, ω
2
5, with ω
1, ω
2As given value, adopt 1000 data in (1) as sample, according to formula (14) structure system of linear equations, obtain y by the least square solution of system of equations
1~ y
4, obtain at last amplitude and the phase place of two signals according to formula (15).
The constructive simulation signal is as follows:
x(t)=a
1cos(2πf
1t+φ
1)+a
2cos(2πf
2t+φ
2)+e(n)(16)
F wherein
1=28.1, f
2=28.283(be rotating speed be respectively 1566 and 1577rpm), a
1=1, a
2=0.5, φ
1=30, φ
2=230, e (n) is white noise signal.
When signal to noise ratio snr=50, shown in 10 measurement result following tables 1.
Table 1 is 10 measurement results when SNR=50
Sequence number | f 1 | f 2 | a 1 | a 2 | ω 1 | φ 2 |
1 | 26.10042594 | 26.28008422 | 1.0026 | 0.5069 | 30.3252 | 50.7706 |
2 | 26.09998165 | 26.2807692 | 1.0004 | 0.5051 | 30.3741 | 50.8167 |
3 | 26.09995846 | 26.28186484 | 1.002 | 0.5032 | 30.1599 | 50.2225 |
4 | 26.10016344 | 26.28014628 | 1.006 | 0.5084 | 30.3457 | 50.405 |
5 | 26.10058556 | 26.2820633 | 1.0049 | 0.5045 | 30.0185 | 49.8685 |
6 | 26.10121111 | 26.27879443 | 1.0137 | 0.5157 | 30.3361 | 50.1702 |
7 | 26.10032799 | 26.28130369 | 1.0049 | 0.5057 | 30.1575 | 50.1219 |
8 | 26.10048734 | 26.28268253 | 1.0029 | 0.5025 | 29.9647 | 49.8353 |
9 | 26.10057793 | 26.28099092 | 1.0064 | 0.5075 | 30.1613 | 50.0876 |
10 | 26.10033701 | 26.28045194 | 1.0062 | 0.5078 | 30.2561 | 50.2589 |
Claims (1)
1. the unbalanced signal extracting method of a Dual Rotor System with Little Rotation Speed Difference is characterized in that, may further comprise the steps:
A1 makes that the sampling period is T, and the vibration signal take T as the cycle to described Dual Rotor System with Little Rotation Speed Difference carries out equal interval sampling, obtains one group of data sequence { D (kT) }, k=1, and 2 ... K, K are a positive integer;
A2, getting centre frequency is f
0, quality factor q>2 adopt the digital band pass tracking filter that { D (kT) } carried out filtering, obtain { y (kT) }, k=1, and 2 ... K, wherein, the pulsed transfer function of digital band pass tracking filter is
A3 with a default rule, makes data to { y (kT) } and extracts, the acquisition data sequence x[n] }, n=0 in the formula, 1 ... N-1 is time series, and N is a positive integer;
A4 establishes { x[n] } and contains consisting of of noisy sinusoidal signal
x[n]=a
1cos(ω
1n+φ
1)+a
2cos(ω
2n+φ
2)+e(n)(1)
N=0 in the formula, 1 ... N-1 is time series, a
1, a
2Be unbalanced signal amplitude, ω
1, ω
2Be unbalanced signal frequency and φ
1, φ
2Be the unbalanced signal phase place,
The capable difference equation matrix of structure N-4 is
X[n]=[x[n]+x[n+4],x[n+1]+x[n+3],x[n+2]](2)
The Prony method shows, the homogeneous solution of above-mentioned equation has been specified e
I ωIn 4 rank multinomial coefficients, because sampled data all is real number, e
I ωIn 4 rank polynomial expressions can be rewritten into formula about cos (2 ω), thereby can solve with the Chebyshev polynomial expression,
By matrix X structure variance matrix G,
M=N-4 in the formula, all summations are limited to up and down from 0 to M-1, and the minimal eigenvalue of variance matrix G is by its secular equation polynomial expression
c
1λ
3+c
2λ
2+c
3λ+c
4=0(4)
Coefficient ask for, this coefficient c is
Then minimal eigenvalue is
Wherein S, T computing formula are,
Minimal eigenvalue λ
0The characteristic of correspondence vector is
Obtain the minimal eigenvalue λ of G matrix by matrix decomposition
0With its characteristic of correspondence vector
Then the total least square solution of matrix X is
Known cos (2 ω)=2 (cos (ω)) by the Chebyshev polynomial expression
2-1(10)
The frequency of any in homogeneous solution v and 2 unknown signalings all satisfies equation
2cos(2ω
k)v
1+2cos(ω
k)v
2+v
3=0(11)
Then the frequency estimation of unknown signaling is
Obtain the frequencies omega of unbalanced signal in the Dual Rotor System with Little Rotation Speed Difference vibration signal
1And ω
2
A5 utilizes least square method to estimate ω
1And ω
2Amplitude and phase place, formula (1) can be write as
x[n]=a
1cos(ω
1n)cos(φ
1)-a
1sin(ω
1n)sin(φ
1)
+a
2cos(ω
2n)cos(φ
2)-a
2sin(ω
2n)sin(φ
2)+e(n)(13)
If y
1=a
1Cos (φ
1), y
2=a
1Sin (φ
1), y
3=a
2Cos (φ
2), y
4=a
2Sin (φ
2), treat that estimated parameter is y
1, y
2, y
3, y
4, according to each time observed reading x[n] (n=1,2 ... N), can draw the system of linear equations (14) of Ay=x, ask its least square solution can pick out each parameter,
The frequencies omega of unbalanced signal then
1And ω
2The estimated value of amplitude and phase place is
φ
1=arctan(y
2/y
1),φ
2=arctan(y
4/y
3)(15)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210407628.4A CN102944360B (en) | 2012-10-23 | 2012-10-23 | A kind of unbalanced signal extracting method of Dual Rotor System with Little Rotation Speed Difference |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210407628.4A CN102944360B (en) | 2012-10-23 | 2012-10-23 | A kind of unbalanced signal extracting method of Dual Rotor System with Little Rotation Speed Difference |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102944360A true CN102944360A (en) | 2013-02-27 |
CN102944360B CN102944360B (en) | 2016-05-04 |
Family
ID=47727322
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210407628.4A Active CN102944360B (en) | 2012-10-23 | 2012-10-23 | A kind of unbalanced signal extracting method of Dual Rotor System with Little Rotation Speed Difference |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102944360B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103780356A (en) * | 2014-02-21 | 2014-05-07 | 上海师范大学 | Design method for two-level precodes of cognitive MIMO communication system |
CN104200118A (en) * | 2014-09-15 | 2014-12-10 | 吉林大学 | Automatic balancing machine vibration signal processing method |
CN114623977A (en) * | 2022-03-12 | 2022-06-14 | 北京化工大学 | Paddle fan coaxial counter-rotating structure automatic balance control method based on micro speed difference |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0399234A (en) * | 1989-09-13 | 1991-04-24 | Ube Ind Ltd | Diagnosing method for abnormality of rotary machine |
JP2010043891A (en) * | 2008-08-11 | 2010-02-25 | Ihi Corp | Device and method of measuring rotation balance of high-speed rotating body |
CN102353500A (en) * | 2011-06-03 | 2012-02-15 | 上海师范大学 | Extraction method of unbalanced signal for dynamic balance measurement |
CN102520246A (en) * | 2011-12-05 | 2012-06-27 | 西安交通大学 | Constant frequency phasor extraction method |
-
2012
- 2012-10-23 CN CN201210407628.4A patent/CN102944360B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0399234A (en) * | 1989-09-13 | 1991-04-24 | Ube Ind Ltd | Diagnosing method for abnormality of rotary machine |
JP2010043891A (en) * | 2008-08-11 | 2010-02-25 | Ihi Corp | Device and method of measuring rotation balance of high-speed rotating body |
CN102353500A (en) * | 2011-06-03 | 2012-02-15 | 上海师范大学 | Extraction method of unbalanced signal for dynamic balance measurement |
CN102520246A (en) * | 2011-12-05 | 2012-06-27 | 西安交通大学 | Constant frequency phasor extraction method |
Non-Patent Citations (3)
Title |
---|
K.S.SHIM等: "Use of Prony analysis to extract sync information of low frequency oscillation from measured data", 《EUROPEAN TRANSACTIONS ON ELECTRICAL POWER》 * |
ZIQIANG ZHANG等: "Optimization and realization of a rotor dynamic balance measuring algorithm", 《PROCEEDINGS OF 3RD INTERNATIONAL CONGRESS ON IMAGE AND SIGNAL PROCESSING(CISP 2010)》 * |
孙晓明等: "采用改进自适应Prony方法的电力故障信号的分析与处理", 《中国电机工程学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103780356A (en) * | 2014-02-21 | 2014-05-07 | 上海师范大学 | Design method for two-level precodes of cognitive MIMO communication system |
CN103780356B (en) * | 2014-02-21 | 2017-03-15 | 上海师范大学 | A kind of method for designing of the two-stage precoding of cognitive MIMO communication system |
CN104200118A (en) * | 2014-09-15 | 2014-12-10 | 吉林大学 | Automatic balancing machine vibration signal processing method |
CN114623977A (en) * | 2022-03-12 | 2022-06-14 | 北京化工大学 | Paddle fan coaxial counter-rotating structure automatic balance control method based on micro speed difference |
Also Published As
Publication number | Publication date |
---|---|
CN102944360B (en) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107505135B (en) | Rolling bearing composite fault extraction method and system | |
CN105092241B (en) | A kind of gear local fault diagnosis method and system | |
CN104375111B (en) | The method that quick high accuracy refinement correction is carried out to intensive spectrum | |
CN104215456B (en) | Plane clustering and frequency-domain compressed sensing reconstruction based mechanical fault diagnosis method | |
Li et al. | Research on test bench bearing fault diagnosis of improved EEMD based on improved adaptive resonance technology | |
CN107843740B (en) | A kind of rotating speed measurement method of fusion vibration and voice signal spectrum signature | |
CN105547698A (en) | Fault diagnosis method and apparatus for rolling bearing | |
CN104502099A (en) | Cyclic frequency extraction method for characteristic components of transient conditions of gearbox | |
CN102353500B (en) | Extraction method of unbalanced signal for dynamic balance measurement | |
Rodopoulos et al. | A parametric approach for the estimation of the instantaneous speed of rotating machinery | |
CN104200118A (en) | Automatic balancing machine vibration signal processing method | |
CN109459131A (en) | A kind of the time-frequency characteristics extracting method and device of rotating machinery multi-channel Vibration Signal | |
CN104034412A (en) | Rotary machine fault feature extraction method based on fractional order holographic principle | |
Lin et al. | A review and strategy for the diagnosis of speed-varying machinery | |
CN108398260B (en) | Method for quickly evaluating instantaneous angular speed of gearbox based on mixed probability method | |
CN102944360B (en) | A kind of unbalanced signal extracting method of Dual Rotor System with Little Rotation Speed Difference | |
Zhao et al. | Adaptive scaling demodulation transform: Algorithm and applications | |
CN111307426A (en) | Rotating machinery fault feature extraction method based on FrFT-EWT principle | |
CN104677486B (en) | The aero-engine vibration signal Method for Phase Difference Measurement reconstructed based on tacho-pulse | |
CN112345247B (en) | Fault diagnosis method and device for rolling bearing | |
CN103913271B (en) | Method for extracting dynamic unbalance signals of rotor at non-stable rotational speed | |
CN105973999B (en) | The faint fractional harmonic characteristic recognition method of rotor crack based on enhancing phase Waterfall plot | |
CN106027454A (en) | Frequency offset estimation method of QAM (Quadrature Amplitude Modulation) signal based on frequency compensation | |
CN114486252B (en) | Rolling bearing fault diagnosis method of vector mode maximum envelope | |
Song et al. | Multispectral balanced automatic fault diagnosis for rolling bearings under variable speed conditions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |