Published: 09 October 2023

Extraction and diagnosis of rolling bearing fault signals based on improved wavelet transform

Zhiqing Cheng1
1Intelligent Manufacturing School, Jiangxi Technical College of Manufacturing, Nanchang, 330095, China
Views 82
Reads 25
Downloads 128


As the continuous growth of the machinery industry, the importance of rolling bearings as key connecting parts in machinery movement is also increasing. However, the extraction and diagnosis of rolling bearing fault signals are difficult, and how to use modern transform analysis methods to raise the extraction efficiency and diagnostic accuracy becomes the focus. For this, a rolling bearing fault signal extraction and diagnosis model is designed based on empirical wavelet transform. The diagnostic model is optimized by using support vector machine and quantum genetic algorithm to design a rolling bearing fault signal extraction and diagnosis model based on improved empirical wavelet transform-support vector machine. The test results show that the research method can obtain four component signals showing different anomalies when generating time domain diagrams. Only five component peaks are generated and one group is extracted as output when generating component peaks. The abnormal amplitude of envelope spectrum basically reaches 0.40×10-6 or above. The judgment accuracy of component diagnosis reaches 98.12%. The above results show that the research method has better fault signal extraction ability and better diagnostic accuracy when performing fault signal diagnosis, which can provide new technical support for rolling bearing fault signal extraction and diagnosis.

1. Introduction

In current society, the speed growth of the machinery manufacturing industry has given a huge impetus to the development of modern industrial production and social economy. However, with the use of industrial equipment continuing to increase, mechanical equipment failure is also gradually highlighted. Rolling Bearings (RBs) in industrial production as an important rotating machinery components assume a key role, the stability and reliability of its operating conditions directly affect the normal operation of machinery and equipment [1, 2]. Long-time operation and harsh working environment will make the bearings produce fatigue damage, wear, cracks and other faults, resulting in the performance of equipment degradation, increased noise, shortened life, and may even lead to safety accidents, so the extraction and diagnosis of bearing failure signals are very important [3]. The main challenge in rolling bearing fault (RBF) diagnosis is to extract and accurately classify the failure signals from a large number of combined vibration signals [4]. Currently, many manners have been utilized to RBF diagnosis, including time and frequency domain analysis and wavelet transform (WT). However, all these methods have certain limitations. Time domain analysis can only provide the basic information of the signal, and cannot analyze the features of the signal comprehensively from frequency and time-frequency domains. Although the frequency domain analysis can provide the spectral information of the signal, it cannot deal with the non-smooth signal well. As a time-frequency analysis method, WT can overcome the above shortcomings, but the traditional WT has a serious energy concentration problem, so the analysis of non-stationary signals is not as effective as it should be [5]. The latest technical contributions of support vector machine (SVM) in the field of intelligent fault diagnosis are mainly reflected in improved kernel function and feature extraction methods, multimodal data fusion, enhanced learning and semi-supervised learning methods, as well as online learning and incremental learning methods. The continuous development and innovation of these technologies provide more powerful and efficient solutions for intelligent fault diagnosis. In view of this, research attempts to utilize the advantages of empirical wavelet transform (EWT) in spectrum extraction, and innovatively combines SVM and quantum genetic algorithm (QGA) to handle nonlinear problems. A multi-method RBF signal extraction and diagnosis model based on improved wavelet transform is designed, with the aim of providing feasible technical improvements for RBF diagnosis through emerging technologies.

2. Related works

Nowadays, RB as a common connected component in various mechanical structures, the fault signal composition is more complex, and scholars all over the world have conducted a lot of research on the extraction and diagnosis of fault signals of RB. Xu et al. [6] proposed a time-based extraction method for the problem of difficult selection of non-smooth signal windows for RBFs. The process used S-transform for time domain extraction to reduce the burden of subsequent methods. The outcomes proved that the method could extract the pulse occurrence time accurately and had good function when performing RBF diagnosis. Pang et al. [7] proposed an optimization method based on variational modal extraction for the problem of many disturbing components of RBF vibration signals. The process set a new fault evaluation index, and the smoothness of the cycle was used as the fitness function. Research findings validated that the method had high analytical accuracy when performing resonant demodulation analysis. An et al. [8] put forward an LSTM-based optimization method for the problem of RBF signals affected by periodic inhomogeneity. The process used periodic sparse attention networks to enhance the time utilization efficiency, followed by feature extraction using LSTM. The outcomes proved that the method could well diagnose RBFs. Fu et al. [9] proposed an AC-FEGAN-based fault diagnosis manner for the inefficiency of RBF data labeling. The process used unsaturated loss to circumvent gradient loss and introduced residual network for feature extraction. The outcomes demonstrated that the method could raise the quality of fault samples and could perform fault diagnosis effectively. Shang et al. [10] proposed a DTW-based fault diagnosis manner for the inconspicuous RB features. The process used DTW to calculate the residual vector and designed a deep confidence network for learning rate setting. The lab findings expressed that the method had a high diagnostic accuracy.

WT has also been studied by some scholars. Shrivastavaet al. [11] proposed a WT-based encryption method for the secrecy of images in multimedia technology. The process used the Fourier transform for encoding and the wavelet family as a key. Research results validated that the proposed method had higher image security and better flexibility in transmission. Ma et al. [12] proposed a fault diagnosis method based on WT for the problem of parallel faults in battleship currents. The process used computationally light data drive for transient feature extraction and identification of abnormal disturbance types. It was proved that the proposed method could be effective for fault diagnosis. Krishnan and Soman [13] proposed an optimized EWT-based identification method for the problem of EEG signal identification. The process investigated time-dependent desynchronization and used variational modal decomposition incorporating EWT for multiple electrode EEG signals for analysis and extraction. The laboratory outcomes denoted that the proposed method had good recognition accuracy. Alharbey et al. [14] proposed a detection method using WT for the problem of arrhythmia detection. The process used standard deviation and Shannon entropy for feature extraction and established a safety threshold to classify the heart rate duration. The findings demonstrated that the proposed method could perform reasonable arrhythmia detection with good accuracy. Tripathyet al. [15] proposed a WT based disease detection classification method for the problem of lung disease detection. The process fixed the boundary points and extracted the frequency domain features from each disease pattern. The outcomes of the study indicated that the method had good accuracy in classifying lung diseases.

In summary, WT technology has been applied in many fields, but little research has been done on RBF signal extraction and diagnosis. To address RBF signal extraction and diagnosis, it proposed a RBF signal extraction and diagnosis method that integrated SVM and QGA into EWT to provide more reference solutions for RBF signal extraction and diagnosis by using the advantages of WT technology in attenuating modal confusion.

3. IEWT-SVM based RBF signal extraction and diagnosis model design

3.1. Design of IEWT-RBF signal processing model based on normalization processing

RB in the normal working process will produce vibration because of the surface roughness, size deviation, etc. The signal issued by this vibration is normal smooth signal [16]. When the RBF point, rolling body and the fault point produce signal fluctuations, in a certain rotation speed, the corresponding fault frequency signal will be formed. This signal is not smooth and non-linear [17, 18]. Typical RB faults can be divided into three types: inner ring, outer ring, and ball bearing faults. The inner and outer rings are both circular parts that generate vibration when faults occur, and each of them only generates one fault frequency. Ball bearings are rolling parts inside bearings that come into contact with both the inner and outer rings simultaneously, resulting in two failure frequencies when a fault occurs. There is a certain difference between the diagnosis of inner and outer ring faults and the diagnosis of ball faults. Research is conducted to design methods for diagnosing inner and outer ring faults in RBs. EWT is a non-linear signal adaptive decomposition method, which can effectively attenuate the impact of modal confusion by applying in the RBF signal processing. The process of generating fault signals from RBFs is displayed in Fig. 1.

Fig. 1Process for generating fault signals of RBs

Process for generating fault signals of RBs

As seen in Fig. 1, RBF signals are generated by multiple factors. The natural working vibration of the RB itself and the vibration caused by dimensional deviation are the most basic vibrations. When lubrication conditions change and external pressure and bearing failure occur during operation, the vibration frequency and amplitude change and the output failure signal is hidden in the vibration signal. As the RB rotates, these signals are periodic. The vibration signal is first processed by the Fourier transform, and the resulting signal spectrum is normalized as shown in Fig. 2.

As seen in Fig. 2, the normalization process takes the half-perimeter rotation as a complete processing interval, and divides the interval into consecutive N segments. ωN is set as the center of each band, then ω0 and ωn denote the upper and the lower boundary of the interval, respectively, and the length of the transition band is 2τn. The empirical wavelet function is constructed as shown in Eq. (1):

ψ^n(ω)=1, 1+γωnωn1-γωn+1,cosπ2β12γωn+1ω-1-γωn+1 ,1-γωn+1ωn1+γωn+1,sinπ2β12γωnω-1-γωn , 1-γωnωn1+γωn,0, otherwise,

where, ωn denotes the lower boundary of the n band; the value of γ is less than 1. The empirical scale function is constructed as shown in Eq. (2):

ϕ^nω=1, fωn1-γωn,cosπ2β12γωnω-1-γωn, 1-γωnωn1+γωn,0, otherwise,

where, γ is less than the ratio of the difference between the upper and lower boundaries of the frequency band to the sum of them. The definition of the modal component is shown in Eq. (3):


where, Wfε(n,t) means the detail coefficients, which are generated by the wavelet function and the signal inner product; Wfε(0,t) indicates the approximation coefficients, which are produced by the scale function and the signal inner product; ϕ1 denotes the empirical scale coefficients. The output vibration signal is reconstructed as shown in Eq. (4):


where, W^fε(0,ω), W^fε(n,ω) are the Fourier transform of Wfε(0,t), Wfε(n,t), respectively. To perform the spectrum differentiation, the spectrum is divided into N bands, and each band exists independently, excluding the two end points of the interval, and it also needs to determine the boundary points of N-1. The discrete scale space of the spectrum within a single interval is determined by the number of minima, and the function for calculating the number of minima is shown in Eq. (5):


where, x is the variable; t means the scale parameter; g(n,t) indicates the Gaussian kernel function. When the RB is in a harsh working environment, the factors affecting the vibration become more and the generated vibration signal is more complex, and a large number of redundant and meaningless frequency bands will be generated when the frequency band is divided, which affects the subsequent RBF signal extraction efficiency. To improve the efficiency of fault signal extraction, Otsu is used to optimize the EWT process, and the IEWT transformation process is obtained, as shown in Fig. 3.

Fig. 2Spectrum processing diagram

Spectrum processing diagram

Fig. 3IEWT transformation process

IEWT transformation process

As seen in Fig. 3, the IEWT transform processes based on Otsu are as follows. First, it performs FFT transform on the input signal, and analyzes the Fourier spectrum. Then, the scale space curve threshold is determined by the Otsu method, and the initial dividing point is derived. It decides whether to merge the components according to their mutual information and the values of adjacent components, updates the Fourier spectrum according to the merged components, and outputs the optimized component signal [19, 20]. The vibration characteristics at the beginning of the RBF are weak and easily disturbed by other environmental noise when extracted. The vibration signal function of RBF is established as shown in Eq. (6):


where, y(n) refers to the output signal; x(n) denotes the input signal; h(n) is the signal transfer function; and e(n) represents the environmental noise. The smaller the entropy recovered to x(n), the less the fault information is disturbed [21]. The value of x(n) is taken as shown in Eq. (7):


where, l represents the inverse filter length. The size of the entropy value of x(n) is measured by the target parametrization of x(n), and the objective function is established as expressed in Eq. (8):


where, o24(f(n)) is the parametric amount of x(n). The partial derivatives of the entropy value and the parametric objective function are coupled as shown in Eq. (9):


where, B is the inter correlation matrix of the input and output data obtained from the association; g denotes the inverse filter matrix obtained from the iterative calculation; A stands for the L×L autocorrelation matrix of the signal y(n).

3.2. Optimized design of IEWT RBF signal extraction and diagnosis model based on SVM and QGA

After completing the pre-processing of the vibration signal of RBF, subsequent extraction of the fault signals is required to enter the diagnosis process. The study optimizes and reconstructs the RBF diagnosis process after the introduction of MED [22]. The optimized RBF diagnosis process is shown in Fig. 4.

Fig. 4Optimized RBF diagnosis process

Optimized RBF diagnosis process

In Fig. 4, the optimized RBF diagnosis firstly separates the noise and impulse components from the fault information to obtain the modal component signal. After that, the modal component information is reconstructed using the peakness and mutual information criterion, and the components with the largest mutual information and those with too small peakness are excluded in the reconstruction. After obtaining the optimal modal components, the MED is used for noise reduction. The output vibration signal is analyzed in various aspects, and the data integration is completed. Finally, the fault diagnosis results are obtained. In the real working condition, it is not easy to sample the RBF, and the sample size of the fault signals that can be obtained is small, so the conclusions obtained by using traditional statistical analysis methods often have large deviations. SVM has good generalization ability when identifying and classifying, and has better anti-interference ability under small sample data, so it can be utilized to RBF diagnosis to improve the diagnosis accuracy [23]. The SVM works as shown in Fig. 5.

Fig. 5SVM classification plane

SVM classification plane

a) Initial classification plane

SVM classification plane

b) Insert optimal classification plane

In the basic case, SVM sets an optimal classification plane between two types of data point groups by analyzing the data points. It also inserts the optimal classification plane to squeeze the two types of data points towards both sides, and expands the classification interval, thus distinguishing the two types of data. The data are thus completely separated. The definition of the optimal hyperplane is shown in Eq. (10):


where, xi represents the input feature vector; yi indicates the sample category; ω means the weight vector; and b means the bias vector. The classification interval is shown in Eq. (11):


where, margin represents the classification interval value. The larger the classification interval is, the better the data points are classified, and the optimal classification plane selection is transferred to quadratic convexity, as shown in Eq. (12):

min12ω2,s.t.yiωxi+b1, i=1,2,,n,Lω,b,a=12ω2-i=1naiyiωxi+b-1,

where, ai is the introduced Lagrange multiplier. The decision function of the optimal classification plane is shown in Eq. (13):

maxi=1nai-12i=1ni=1nyiyjaiajxixj,s.t.ai0, i=1,2,,n,i=1naiyi=0,fx=sgnωx+b=sgni=1nai*yixi,yj+b*,

where, f(x) represents the optimal classification plane decision function, which is derived from the optimal classification plane solution by quadratic search for the best solution. When nonlinear data are involved, the introduction of the kernel function can map the features to a high-dimensional space, thus transforming the nonlinearly indistinguishable data into linearly distinguishable data in a high-dimensional space. The pairwise problem Lagrange is transformed as shown in Eq. (14):

minαQ(α)=k=1nαk+12k=1ni=1nykyiakaiKxkxi,s.t.i=1nakai=0,0akc, k=1,2,,n,

where, c is a constant. The decision function of the nonlinear SVM undergoes a transformation, as shown in Eq. (15):


where, yk represents the sample category; sgn is the symbolic function. The Gaussian radial basis kernel function is utilized as the kernel function, as shown in Eq. (16):


where, σ expresses the kernel parameter, and the smaller the kernel parameter, the better the learning ability of SVM. The ranking entropy is able to analyze the dynamic changes of the signal by integrating the temporal connections between data points, and the entropy value is calculated without setting the probability distribution of the signal in advance. The definition of the ranking entropy is expressed in Eq. (17):


where, Hp(m) is the definition of alignment entropy according to the Shannon entropy setting; k represents the number of reconstruction vectors; ln(m!) is the maximum value obtained by Hp(m) under the condition of pl=1/m!; the value of Hp takes the range of 0Hp1, the smaller the value is, the weaker the randomness of the time series and the lower the sequence complexity. The values of optimal kernel parameter width and penalty coefficient affect the accuracy and generalization ability of SVM. Too large a kernel parameter width or too small a penalty coefficient will lead to insufficient classification accuracy; too small a kernel parameter width or too large a penalty coefficient will lead to insufficient generalization ability of SVM [24]. Using QGA for parameter selection can improve the efficiency of optimal parameter selection. The flow of optimal parameter selection using QGA is shown in Fig. 6.

Fig. 6Optimal parameter selection process based on QGA

Optimal parameter selection process based on QGA

As seen in Fig. 6, for parameter selection, the number of populations, the maximization algebra, and the range of penalty coefficients need to be determined first. Then, the populations are initialized and binary measurements are transformed into binary string codes. The randomly generated parameter combinations are input to the SVM for processing, cross-validation training is performed, and the resulting best parameter combinations are retained. The search is continuously verified using a quantum revolving gate until the fitness function satisfies the termination condition and the resulting best optimized parameters are output. The initial population is shown in Eq. (18):


where, N represents the number of randomly generated chromosomes; Pnt denotes the first n individual that evolves to the t generation; P(t) is encoded by quantum bitwise. The quantum revolving door search is shown in Eq. (19):


where, U(θi) is the revolving door operation; θi represents the rotation angle. The RB fault feature extraction and diagnosis model designed by this research integrates multiple methods to achieve the fusion of multimodal features, which can capture RB fault information from different dimensions and scales. By introducing QGA, representative features can be adaptively selected, and redundant information can be eliminated. Compared to traditional methods that combined wavelet transform with SVM, it has better diagnostic efficiency. When conducting fault diagnosis, research methods can provide a certain degree of interpretability for the diagnostic process and results. In practical applications, compared to general wavelet transform methods, they have better applicability and can reduce the information burden on engineering personnel.

4. Performance testing and application efficiency analysis of IEWT-SVM based RBF signal extraction and diagnosis model

4.1. Performance test of IEWT-SVM based RBF signal extraction and diagnosis model

To verify the effectiveness of the IEWT-SVM-based RBF signal extraction and diagnosis model at runtime, the study first tested the performance of the constructed diagnosis model. During the test, the data used in the experiment came from the bearing database of the electrical laboratory of Case Western Reserve University in the United States. This database is the most commonly used data set in the field of RB fault diagnosis, which belongs to the standard data set. The data of 6205-2RS JEM SKF deep groove ball bearing at the drive end of the database were selected for calculation, with the sampling frequency of 12 kHz. The sampling frequency of the simulation signal was set to 20 kHz. The damping coefficient was 0.1. The displacement constant was 5. The natural frequency was 3 kHz. The rotation frequency was 25 Hz. The impact amplitude was 1.2 m/s2. The noise amplitude was 3.0 m/s2. The harmonic amplitude was 1.5 m/s2. The attenuation rate was 800. The inner ring diameter of the experimental bearing was 25 mm. The outer ring diameter was 52 mm. The pitch diameter was 39.04 mm. The diameter of the ball is 7.94 mm. The number of rolling balls was 9. The contact angle was 90°. Two faults were set in the rolling element, with dimensions of 0.1778 mm and 0.3556 mm, respectively. When conducting application analysis, a double row roller bearing with model FAG F-80781109 TAROL 130/240-B-TVP was used, and a fault groove with a length of 5 mm, a depth of 0.7 mm, and a width of 1 mm was manufactured using electric spark. The fault experimental system used is shown in Fig. 7.

Fig. 7Fault experimental system

Fault experimental system

As shown in Fig. 7, the experimental system used was established on a high-speed train bearing experimental platform, where sensors were installed at different positions to collect data. The high-speed train bearing experimental platform was set to rotate at 1200 revolutions per minute. When conducting experiments, wavelet transform was performed on data in different states, retaining the first five modal components, and then the arrangement entropy was calculated to obtain the input feature vector set. The content of system performance testing included dividing spectrum performance testing and simulating signal time-domain diagrams. After completing parameter optimization and training, the envelope spectrum, component kurtosis, diagnostic accuracy, and component diagnostic results would be used as application analysis content. The study first tested the performance of EWT and IEWT-SVM for dividing the spectrum, and the test results are shown in Fig. 8.

Fig. 8Divide spectral images

Divide spectral images

a) EWT original partition spectrum

Divide spectral images

b) IE WT-SVM repartition spectrum

As seen in Fig. 8(a), in the same spectrum interval, IEWT divided the interval into 37 bands according to the band anomalies, and each band was short in length and contained less information. Some bands divided the same anomaly into multiple small bands, which also generated more meaningless bands. In Fig. 8(b), IEWT-SVM determined a new cut-off point according to the overall trend, and divided the interval into four frequency bands according to the new cut-off point and the anomalies, preserving the integrity of each anomalous band and effectively reducing the generation of redundant bands. The time domain diagrams of the simulated signals in the same time period were compared, and the results are shown in Fig. 9.

Fig. 9Time domain diagram of simulated signal

Time domain diagram of simulated signal

a) EWT time domain diagram of simulated signal

Time domain diagram of simulated signal

b) IE WT-SVM component signal time-domain diagram

As seen in Fig. 9(a), a complete time domain map containing all the information was obtained by IEWT in the 0-0.25 s time period, in which the overall trend and fluctuation of the signal could be roughly distinguished, but the specific composition of the signal and the amplitude of each component could not be judged intuitively. In Fig. 9(b), IEWT separated the simulated signal and obtained four component simulated signals. The harmonic, shock and noise components were separately represented, and the trend and variation of each component could be judged more intuitively. It indicated that IEWT-SVM could effectively extract the abnormalities in the signal when performing RBF signal extraction and diagnosis, and could output better diagnosis data.

4.2. Analysis of the application effectiveness of IEWT-SVM based RBF signal extraction and diagnosis model

The envelope spectrum of the wheel-to-wheel bearing containing the inner ring spalling fault was extracted, and the envelope spectrum generated by using ordinary spatial noise reduction was compared with that generated by IEWT-SVM, and the results are shown in Fig. 10.

Fig. 10Envelope spectrum after noise reduction

Envelope spectrum after noise reduction

a) Normal space noise reduction

Envelope spectrum after noise reduction


In Fig. 10(a), the overall amplitude of the envelope spectrum generated by normal spatial noise reduction was low in the 500 Hz length, with the highest value below 0.25 (m⋅s-2). The amplitude fluctuation was obvious, but the obtained spectrum fluctuated frequently, leading to confusion between some normal high amplitude values and abnormal amplitude values around 0.01 (m⋅s-2). In Fig. 10(b), the overall amplitude of the envelope spectrum generated by IEWT-SVM was also low. After the noise reduction optimization process, some of the highest amplitude values reached 0.88×10-6 and the abnormal amplitude values basically reached above 0.40×10-6, which had obvious high and low differences with the normal amplitude values and could distinguish the abnormal amplitude values from the normal amplitude values more effectively. It was also necessary to generate the peak value to derive the reconstructed signal when judging the fault, and the ensemble empirical mode decomposition (EEMD) was compared with the component peak situation generated by IEWT-SVM, and the results are shown in Fig. 11.

Fig. 11Ensemble empirical mode decomposition (EEMD) vs with the component peak situation generated by IEWT-SVM

Ensemble empirical mode decomposition (EEMD) vs  with the component peak situation generated by IEWT-SVM


Ensemble empirical mode decomposition (EEMD) vs  with the component peak situation generated by IEWT-SVM


In Fig. 11(a), EEMD generated 8 component quantitative peaks with the maximum peak of 4.4533 and the minimum peak of 2.6853, and finally output 3 components that were larger than the average value of 3.3766. As seen in Fig. 11(b), IEWT-SVM generated 5 component quantitative peaks, with the maximum peak of 22.1611 and the minimum peak of 2.7686. Finally, it output 1 component quantitative peak with the maximum peak, which could reduce the pressure of subsequent data processing and improve efficiency. The kernel function width and penalty coefficients were iteratively searched for using SVM-based GA and IEWT-SVM, and the accuracy of RBF diagnosis was derived using the obtained kernel function width and penalty coefficients, and the outcomes are shown in Fig. 12.

Fig. 12Iterative optimization results

Iterative optimization results

a) GA

Iterative optimization results


As seen in Fig. 12(a), the GA reached the highest accuracy after 50 iterations out of 200 iterations, with a short rising stagnation between the 20th and 30th iteration intervals, with an initial accuracy of 90.4 % and a final stable maximum accuracy of 94.6 %. As seen in Fig. 12(b), the IEWT-SVM reached the highest accuracy after 24 iterations out of 200 iterations, with a short rising stagnation between the 3rd and 10th iteration intervals and a final stable maximum accuracy of 96.7 %. It indicated that IEWT-SVM had better iteration speed and judgment accuracy. The results were shown in Fig. 13, using GA, Particle Swarm Optimization (PSO) method and IEWT-SVM for component diagnosis under the generated optimal kernel function width and penalty coefficients, and three sample labels were set for diagnosis.

Fig. 13Component diagnosis results

Component diagnosis results

a) GA

Component diagnosis results

b) PSO

Component diagnosis results


When conducting component diagnosis, all three methods classified faulty and faultless samples. From Fig. 13(a), during component diagnosis by GA, there were no misjudgments for the classification that originally belonged to the first class label results, but some samples of the second and third class labels were mistakenly classified into the first class label. In the classification results, there were 4 errors in selecting faultless samples, and the overall judgment accuracy was 90.33 %. From Fig. 13(b), when PSO performed component diagnosis, there was a certain amount of misjudgment in the results of various label samples, but the overall number was less than GA. In the classification results, there were three errors in selecting faultless samples, with an overall accuracy rate of 94.33 %. From Fig. 13(c), IEWT-SVM mistakenly classified a small number of samples belonging to the first class label into the third class label during component judgment, and mistakenly classified a small number of samples belonging to the second class label into the first and third class labels. The classification that originally belonged to the third class label results did not show any misjudgment, and the overall judgment accuracy was 98.12 %. When using research methods for diagnosis, fault samples were more accurately detected, and the proportion of fault samples accurately classified into labels was higher. It denoted that the research method had better model sensitivity and detection accuracy, and the probability of misjudgment was lower. When label classification was performed on faulty samples, only faulty samples were classified, and all faultless samples were excluded, indicating good specificity of the research method.

5. Conclusions

Aiming at the problems of difficult extraction of fault information and high error rate in RBF signal extraction and diagnosis, this research used normalization processing, QGA and SVM for information processing and analysis, and designed an improved RBF signal extraction and diagnosis model. The study first used normalization processing to pre-process the signal spectrum of RBFs and obtained fault information by MED noise reduction. Then, the fault data samples were classified by SVM. QGA was introduced to derive the best calculation parameters, and finally the research method was compared experimentally with EWT, EEMD and GA. The results showed that when performing the spectrum division performance test, EWT generated 37 frequency bands and the research method generated 4 frequency bands. When performing the simulation signal time domain map generation, EWT generated 1 unclassified time domain map containing all information and the research method generated 4 time domain maps showing different disturbance factors, respectively. When performing the noise reduction envelope spectrum analysis, the amplitude of ordinary spatial noise reduction could not clearly distinguish the anomalous. When performing the noise reduction envelope analysis, the normal spatial noise reduction amplitude could not distinguish the abnormal amplitude significantly, while the abnormal amplitude of the research method basically reached 0.40×10-6, which was significantly higher than the normal amplitude. When performing the component diagnosis, the overall judgment accuracy of GA and PSO was 90.33 % and 94.33 %, respectively, which was lower than the 98.12 % of the research method. The above outcomes illustrated that the proposed method had higher extraction efficiency of RBF signals, strong fault information processing capability, and also higher diagnostic accuracy. However, the study was tested in an ideal environment without considering the variability of some actual harsh working conditions, so further testing of the effect in harsh working conditions is needed to enrich the experimental results.


  • Y. Li et al., “Rolling bearing fault diagnosis based on quantum LS-SVM,” EPJ Quantum Technology, Vol. 9, No. 1, pp. 1–15, Dec. 2022,
  • M. Cui, Y. Wang, X. Lin, and M. Zhong, “Fault diagnosis of rolling bearings based on an improved stack autoencoder and support vector machine,” IEEE Sensors Journal, Vol. 21, No. 4, pp. 4927–4937, 2020.
  • H. Wang, F. Wu, and L. Zhang, “Application of variational mode decomposition optimized with improved whale optimization algorithm in bearing failure diagnosis,” Alexandria Engineering Journal, Vol. 60, No. 5, pp. 4689–4699, Oct. 2021,
  • J. Zhang, J. Wu, B. Hu, and J. Tang, “Intelligent fault diagnosis of rolling bearings using variational mode decomposition and self-organizing feature map,” Journal of Vibration and Control, Vol. 26, pp. 1886–1897, 2020.
  • Y. Qiao, L. Chu, X. Chen, C. Guo, and X. Xu, “The control strategy for 4WD hybrid vehicle based on wavelet transform,” SAE International Journal of Advances and Current Practices in Mobility, Vol. 4, No. 1, pp. 182–190, Apr. 2021,
  • Y. Xu, L. Wang, A. Hu, and G. Yu, “Time-extracting S-transform algorithm and its application in rolling bearing fault diagnosis,” Science China Technological Sciences, Vol. 65, No. 4, pp. 932–942, Apr. 2022,
  • B. Pang, M. Nazari, Z. Sun, J. Li, and G. Tang, “An optimized variational mode extraction method for rolling bearing fault diagnosis,” Structural Health Monitoring, Vol. 21, No. 2, pp. 558–570, 2022.
  • Y. An, K. Zhang, Q. Liu, Y. Chai, and X. Huang, “Rolling bearing fault diagnosis method base on periodic sparse attention and LSTM,” IEEE Sensors Journal, Vol. 22, No. 12, pp. 12044–12053, 2022.
  • W. Fu, X. Jiang, C. C. Tan, B. Li, and B. Chen, “Rolling bearing fault diagnosis in limited data scenarios using feature enhanced generative adversarial networks,” IEEE Sensors Journal, Vol. 22, No. 9, pp. 8749–8759, 2022.
  • S. Zhiwu, L. Xia, L. Wanxiang, G. Maosheng, and Y. Yan, “A rolling bearing fault diagnosis method based on fastDTW and an AGBDBN,” Insight-Non-Destructive Testing and Condition Monitoring, Vol. 62, No. 8, pp. 457–463, Aug. 2020,
  • A. Shrivastava, J. B. Sharma, and S. D. Purohit, “Image Encryption Based on Fractional Wavelet Transform, Arnold Transform with Double Random Phases in the HSV Color Domain,” Recent Advances in Computer Science and Communications, Vol. 15, No. 1, pp. 5–13, Jan. 2022,
  • Y. Ma, A. Maqsood, D. Oslebo, and K. Corzine, “Wavelet transform data-driven machine learning-based real-time fault detection for naval DC pulsating loads,” IEEE Transactions on Transportation Electrification, Vol. 8, No. 2, pp. 1956–1965, 2022.
  • K. K. Krishnan and K. P. Soman, “Comparison of variational mode decomposition and empirical wavelet transform methods on EEG signals for motor imaginary applications,” International Journal of Biomedical Engineering and Technology, Vol. 38, No. 3, pp. 267–285, 2022.
  • R. A. Alharbey, S. Alsubhi, K. Daqrouq, and A. Alkhateeb, “The continuous wavelet transform using for natural ECG signal arrhythmias detection by statistical parameters,” Alexandria Engineering Journal, Vol. 61, No. 12, pp. 9243–9248, Dec. 2022,
  • R. K. Tripathy, S. Dash, A. Rath, G. Panda, and R. B. Pachori, “Automated detection of pulmonary diseases from lung sound signals using fixed-boundary-based empirical wavelet transform,” IEEE Sensors Letters, Vol. 6, No. 5, pp. 1–4, 2022.
  • B. Chen, N. Zeng, X. Cao, S. Zhou, W. He, and S. Tian, “Unsupervised learning-driven intelligent fault diagnosis algorithm for high-end bearing,” SCIENTIA SINICA Technologica, Vol. 52, No. 1, pp. 165–179, Jan. 2022,
  • H. Alqahtani and A. Ray, “Forecasting and detection of fatigue cracks in polycrystalline alloys with ultrasonic testing via discrete wavelet transform,” Journal of Nondestructive Evaluation, Diagnostics and Prognostics of Engineering Systems, Vol. 4, No. 3, pp. 1–15, Aug. 2021,
  • X. Zhang, W. Zhang, W. Sun, X. Sun, and S. Kumar Jha, “A robust 3-D medical watermarking based on wavelet transform for data protection,” Computer Systems Science and Engineering, Vol. 41, No. 3, pp. 1043–1056, 2022,
  • W. Qiao, Z. Li, W. Liu, and E. Liu, “Fastest-growing source prediction of US electricity production based on a novel hybrid model using wavelet transform,” International Journal of Energy Research, Vol. 46, No. 2, pp. 1766–1788, 2022.
  • M. Barma and U. M. Modibbo, “Multiobjective mathematical optimization model for municipal solid waste management with economic analysis of reuse/recycling recovered waste materials,” Journal of Computational and Cognitive Engineering, Vol. 1, No. 3, pp. 122–137, Jan. 2022,
  • S. Choudhuri, S. Adeniye, and A. Sen, “Distribution alignment using complement entropy objective and adaptive consensus-based label refinement for partial domain adaptation,” Artificial Intelligence and Applications, Vol. 1, No. 1, pp. 43–51, 2023,
  • Gheyath Othman and Diyar Qader Zeebaree, “The applications of discrete wavelet transform in image processing: A review,” Journal of Soft Computing and Data Mining, Vol. 1, No. 2, pp. 31–43, Dec. 2020.
  • Y. Guo, Z. Mustafaoglu, and D. Koundal, “Spam detection using bidirectional transformers and machine learning classifier algorithms,” Journal of Computational and Cognitive Engineering, Vol. 2, No. 1, pp. 5–9, Apr. 2022,
  • A. Islam, F. Othman, N. Sakib, and H. M. H. Babu, “Prevention of shoulder-surfing attack using shifting condition with the digraph substitution rules,” Artificial Intelligence and Applications, Vol. 1, No. 1, pp. 58–68, 2023.

About this article

06 June 2023
07 September 2023
09 October 2023
rolling bearing
fault sign

The authors have not disclosed any funding.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.