Refining the time

Scientific Reports volume 13, Article number: 5215 (2023) Cite this article

363 Accesses

Metrics details

Time–frequency ridge not only exhibits the variable process of non-stationary signal with time changing but also provides the information of signal synchronous or non-synchronous components for subsequent detection research. Consequently, the key is to decrease the error between real and estimated ridge in the time–frequency domain for accurate detection. In this article, an adaptive weighted smooth model is presented as a post-processing tool to refine the time–frequency ridge which is based on the coarse estimated time–frequency ridge using newly emerging time–frequency methods. Firstly, the coarse ridge is estimated by using multi-synchrosqueezing transform for vibration signal under variable speed conditions. Secondly, an adaptive weighted method is applied to enhance the large time–frequency energy value location of the estimated ridge. Then, the reasonable smooth regularization parameter associated with the vibration signal is constructed. Thirdly, the majorization–minimization method is developed for solving the adaptive weighted smooth model. Finally, the refined time–frequency characteristic is obtained by utilizing the stop criterion of the optimization model. Simulation and experimental signals are given to validate the performance of the proposed method by average absolute errors. Compared with other methods, the proposed method has the highest performance in refinement accuracy.

Time–frequency analysis (TFA) method is an effective tool to provide information on signal synchronous or non-synchronous components in condition monitoring and fault diagnosis under non-stationary conditions. Furthermore, the time-varying features of non-stationary signals could be characterized. TFA methods are wildly applied in radar, sonar and astronomical, biomedicine, and mechanical engineering areas1,2,3,4,5,6, etc. The conventional TFA methods are roughly divided into linear and quadratic transforms, and all of them have respective drawbacks. For example, short-time Fourier transform (STFT) and continuous wavelet transform (CWT) and so on, both of which are difficult in choosing a reasonable window parameter of TFA, which leads to time and frequency resolution in the time–frequency domain7. On the other hand, the classical quadratic transform represented by Wigner–Ville distribution (WVD), the cross-term interferences would be introduced in analyzing multi-component signal8, which decreases the readability of time–frequency, and increases the difficulty of time–frequency ridge extraction.

Mostly, the peak value search algorithm is always applied to extract the peak energy of time–frequency representation for characterizing the procedure of time-varying signal in the industry area. Nevertheless, the obtained peak ridge is a rough curve using the aforementioned time–frequency methods. Therefore, the rough curve is an approximated broken line although constructing a suitable window parameter.

To mitigate the impact of entangled background noises and interferences in analyzing time-varying signals and obtain concentrated time–frequency representation, the post-progressing tool is introduced to solve the above problems. Auger9,10 proposed a reassignment (RM) technique to concentrate the time–frequency energy into a narrow band. After that, the synchrosqueezing transform (SST)11 is proposed to squeeze the time–frequency coefficients into the instantaneous frequency (IF) trajectory along the frequency axis, the method could provide fine time–frequency readability. In other words, the blurry time–frequency representation is concentrated by using a synchrosqueezing operator when analyzing a stationary signal, as a result, an accurate time–frequency representation is obtained12. Nevertheless, the fitted time–frequency curve is heavily biased in comparison with the real IF when analyzing chirp signals or frequency-modulated signals13,14. Several years ago, Yang proposed a series of parametric time–frequency analysis methods to characterize the variety of the time-varying signal15,16,17. It is worth mentioning that author extended the conventional linear chirp kernel to a polynomial chirplet transform (PCT) by constructing a polynomial nonlinear chirplet kernel to replace the chirplet kernel in the chirplet transform. In the same way, spline-kernelled chirplet transform (SCT) is developed. (Weierstrass approximation theorem is applied to guarantee that any continuous function on a closed and bounded interval can be uniformly approximated on that interval by a polynomial to any degree of accuracy, however, the order value should be determined in advance15). Although the time–frequency trajectory of a time-varying signal is well-fitted, the time–frequency representation energy is blurry. In recent years, some useful improved techniques is proposed to process non-stationary signals, second-order STFT-based SST (FSST2)18 and high-order SST19 are developed to match amplitude modulation (AM) and frequency modulation (FM) multi-component signals20, meanwhile, the time–frequency energy is concentrated into a narrow band. However, the complexity and diversity of practical cases are difficult to determine the accurate parameters of IF17,21. Yu proposed an iterative technique to improve the time–frequency energy concentration compared with the SST method, iterative technique not only processes time-varying signals but also has been validated in the advantage of concentrating energy by computing the index of Rényi entropy22. Although the time–frequency readability is obtained by introducing a high-order synchrosqueezing operator and iterative techniques, the estimated time–frequency trajectory is broken-line.

The smoothing technique is wildly applied in scientific research and industry areas. The sampled data are always affected by vibration, electromagnetic interference, transmission path, quantization error, and so on; consequently, the obtained data is mutational, with spikes, and jump23,24,25. Therefore, it is important to confirm the obtained data is reliable and available before signal processing. Aimed at solving the problem of broken-line for extracting time–frequency trajectory and then to refine the rough curve for obtaining a more accurate curve. Firstly, Yang applied by PCT or SCT method to obtain IF trajectory16,26, secondly, searching the peak values of time–frequency representation and then fitted it, and finally, smoothed the rough curve by the least square method (LSM) to obtain more accuracy estimated curve. Nevertheless, if the feature matrix is non-invertible or ill-conditioned, the analytical solution of LSM cannot be obtained. Non-invertible means that the data is linear correlation and redundancy. For an ill-conditioned matrix, the obtained analytical solution is sensitive to little change in a coefficient matrix or constant term. Therefore, the regularization term is added to the optimal function to avoid the aforementioned problems. The most famous methods named ridge regression, least absolute shrinkage, and selection operator (LASSO). The regularization term of ridge regression is L2-norm which is differentiable. Nevertheless, the super-parameter selection is a great important problem in a ridge regression model. In 2017, Chen proposed a method that formulates an optimal demodulation problem to construct a time–frequency filter bank for obtaining a narrow-band signal27. The author applied a ridge regression model to smooth the time–frequency curve, the smaller penalty parameter is constructed to ensure a smoother time–frequency trajectory28. L1-norm is applied in the LASSO model, which could select an argument and squeeze the coefficient of the negligible argument in zero value. Therefore, the LASSO model also called the smooth model in optimization fields and it is a perfect tool to de-noise vibration signals. L1-norm could not differentiable at zero point and the obtained solution is not analytical. Sometimes, the regularization parameter is always set constant value instead of changing with the signal, and L1-norm and L2-norm are applied to avoid non-invertible and sensitivity to little change of coefficient matrix or constant term in the LSM method.

Therefore, in this article, an adaptive weighted smooth model (AWMM) is proposed to solve the aforementioned problems. The regularization parameter associated with the vibration signal is constructed, which does not depend on any prior knowledge of the tested signal. Furthermore, the prior regularization parameter can be determined by the signal itself. Majorization–minimization (MM) method is introduced to solve the problem of non-differentiable at zero point. Based on the estimated coarse time–frequency ridge by the multi-synchrosqueezing transform (MSST) method22, the ridge is smoothed and then achieves high accuracy using AWMM. The proposed model not only could eliminate the unrelated components of the estimated coarse IF but also provide the refined IF accurately.

This article is organized as follows: the theoretical background of MSST and AWMM is displayed in "Method". The completed refine procedure of the vibration sigal is shown in "Numerical simulations". In "Experiment investigation", the performance of the proposed method is validated by simulation and experimental signals. Finally, the conclusion is shown in "Conclusion".

Inspired by the formula of IF smooth construction in27, the constructed model could be further to be improved, because the key penalty parameter of the model is difficult to be determined. In this section, a signal-driven technique is introduced to solve the above problem and the model is enhanced to improve the IF accuracy in linear and nonlinear time-varying conditions. The common optimal models are used to eliminate the unconcern components of signals and make errors between the estimated and the actual values decrease, for example, LASSO and ridge regression et al. To convenient to express the above two methods, the former called L1-based optimal function and the latter named L2-based optimal function.

The smooth model is constructed as follows:

where the \(\tilde{f}\) is the calculated coarse IF of the signal, the estimated IF could be a nonlinear curve, thus \(\tilde{f} = [\tilde{f}(t_{0} ),\tilde{f}(t_{1} ),...,\tilde{f}(t_{N - 1} )]\), and the \(f\) is the corresponding refined IF, \(f = [f(t_{0} ),f(t_{1} ),...,f(t_{N - 1} )]\). The constructed model can refer to27,29. To decrease the end effects caused by the difference operation, the second-order difference matrix is given as \({\mathbf{D}} = \left[ {\begin{array}{*{20}c} { - 1} & {2} & {\begin{array}{*{20}c} { - 1} & {} & {} \\ \end{array} } & {} \\ {} & { - 1} & {\begin{array}{*{20}c} {2} & { - 1} & {} \\ \end{array} } & {} \\ {} & {} & {\begin{array}{*{20}c} {} & {...} & {} \\ \end{array} } & {} \\ {} & {} & {\begin{array}{*{20}c} {} & { - 1} & 2 \\ \end{array} } & { - 1} \\ \end{array} } \right]\) , the size of the matrix \((N - 2) \times N\), and the N has defined as the length of \(f\) and \(\lambda\) is the regularization parameter. It is important to set a suitable \(\lambda\) initially, in the subsequent section, the rule of the determined parameter would be given. The penalty term of the proposed model is to let the coefficients of the signal approximate zero or equal to zero and further eliminate the unrelated components of the signal. Sometimes, the regularization parameter is always set constant value instead of changing with the signal and it is obvious that the same parameter corresponding to each point is unsuitable. Therefore, an adaptive weighted technique is introduced to address the above problem, and then the initial regularization parameter is determined by the signal. The corresponding formula of the initial value is set

and the formula of adaptive weighted is given

where the j is the iteration count and the value is from 1 to J, when \(j = 1\), the weight matrix could be a unit matrix I, \(W = diag(w_{1} ,w_{2} ,...,w_{N - 2} )\), the size of this matrix is \((N - 2) \times (N - 2)\). Equation (1) could be rewritten as Eq. (4).

To express the penalty term briefly, the \(W^{\prime} = \lambda_{0} WD\), thus, the refined IF \(f\) could be calculated as

Considering that the conventional penalty function is non-differentiable at zero point, the majorization–minimization (MM) algorithm is applied to realize non-differential elimination at zero point, the pivotal of the MM algorithm is to seek a majorizer \(G(f,u)\) of \(F(f)\), and the majorizer \(g(f,u)\) of \(\varphi (f)\), (\(\varphi (f) = \left\| {W^{\prime}f} \right\|_{1}\)), they must meet the following formula:

where the two variables satisfy the condition \(f,u \in R\), after calculating Eqs. (7) and (8), the obtained equations are as follows

in which the unknown variables m and b could be solved by the method of undetermined coefficients (MUC), the obtained equations are shown

Therefore, the majorizer function is detailed to represent in Eq. (13)

It is noted that the L1 norm of the optimal model can be defined as

therefore, Eq. (13) can be revised as

where the diagonal matrix is represented \([\Lambda (u)] = diag(\frac{{\varphi^{\prime}(u_{n} )}}{{u_{n} }})\) and scalar \(b(u) = \sum\nolimits_{n = 0}^{N - 1} {[\varphi (u_{n} ) - \frac{1}{2}\varphi^{\prime}(u_{n} )} ]\). When we considered the penalty term of Eq. (4), Eq. (15) would be revised

when the equation \(f = u\), the majorizer function \(G(f,u)\) is given

this is a minimized problem and its analytical solution could be obtained

The proposed model not only could eliminate the unrelated components of the estimated coarse IF but also provide the refined IF accurately. The important parameters of the model can be determined adaptively based on the signal itself.

In this section, linear and non-linear simulated signals are used to demonstrate the capability of the AWMM to smooth the time–frequency curves. We focus on the comparisons between the AWMM method and other common smooth techniques in addressing linear and nonlinear signals. The comparisons mainly focus on the smooth accuracy between the real IF curve and the post-processed curve. Due to the mean absolute error (MAE) does not appear positive and negative off settings in the assessment of the error of estimated and real values, this index is introduced in this article to measure the performance of the proposed method. The absolute is a mathematical function that makes a number positive. The obtained MAE value is less than 1. Especially, the MAE value will no longer be calculated if the calculation results of the comparison method are too different. The sampling frequency is 100 Hz. It is necessary to compare similar results to test the performance of the mentioned methods. The refined and real curves differ greatly, and the calculated index value of MAE is meaningless. Considering that the experimental data provide the coarse curves, we develop the comparison cases of time–frequency analysis methods in simulation parts. We use traditional and enhanced time–frequency analysis methods to check the performance of AWSM, such as CWT and SST. CWT is a multi-resolution analysis method, which can process stationary and non-stationary signal well. SST can concentrate time–frequency energy into a limit band for separating signal's components.

Herein, a linear simulated signal is modeled and its corresponding IF is given

where the time duration is 4 s. The obtained results of the L2-based and L1-based optimal functions, the proposed method, and polynomial curve fitting-based LSM are exhibited in Fig. S1, all results are presented in the Supplementary information document, namely "All computed compared results.pdf".

The smooth result generated by the proposed method is given in Fig. 1a, which matches the point from 0.45 to 3.58 s, the obtained fitting region is the biggest for all the mentioned methods. That is 3.13 s. The performance of the proposed method is verified by the calculated fitting region. The real and estimated curves is displayed in Fig. 1b, which the red presents coarse IF and blue is defined as real IF. The coarse IF is extracted from time–frequency plane using MSST method. As a consequence, the matched region results of the refined line by using the above methods are given, the AWMM method could match most of the IF trajectory. The accuracy level of the above method is testified by calculating the index of MAE, the calculated result is 0.0411, which is smaller than the MAE of Fig. 1b. The MAE of real and estimated curves is 0.0791. To a certain extent, the proposed method could improve the accuracy of the refined curve. Furthermore, the calculated the index values of all methods are demonstrated in Table S1, which is in the Supplementary information document, namely "All computed compared results.pdf". The results of SST and CWT methods are displayed in Fig. S2, meanwhile, the corresponding smooth results is filled in Table S2, which is in the Supplementary information document, namely "All computed compared results.pdf".

Simulated signal. (a) Obtained result by using the proposed method, (b) the estimated and real curves.

The sampled time is 6.5 s and the simulated non-linear signal is as follows:

This test is to consider the performance of the AWMM method in smoothing the sinusoidal signal. In the same way, the estimated coarse curve is smoothed by the penalty function, the proposed method, and the LSM method, and the corresponding results are given in Fig. S2, all results are presented in the Supplementary information document, namely "All computed compared results.pdf". Similarly, the red line is the real IF curve and the blue one is the smooth curve. In the non-linear time–frequency ridge refinement case, the enlarged locations are the peak and trough of the rough curve. The calculated result by using the AWMM method is shown in Fig. 2a. No matter the peak or trough location, the refined curve is fitted accurately. Compare with Fig. 2b, that consists of real and estimated curves, the calculated MAE value is 0.0553, which is bigger than the value of the AWMM method. Most points are matched with the red line and the MAEs of the aforementioned methods are calculated as in Table S2, which presents in the Supplementary information document. The minimum value belongs to the method proposed in this section and provides the refinement IF with the highest accuracy in nonlinear time-varying conditions. On the other hand, the performance of the MSST method is validated by comparing SST and CWT methods, the results are shown in Fig. S4 and Table S4, they are presented in the Supplementary information document.

Simulated signal. (a) Obtained result by using the proposed method, (b) the estimated and real curves.

The calculated the index values of all methods are demonstrated in Table S2 and S2, from the tables, although the MAE value of the proposed model is not the smallest in linear case, the most curve trajectory is tracked by comparing with LSM and L2-based model. Both of LSM, L2-based model and AWSM method all have very close value.

In nonlinear case, AWSM has the smallest MAE value by comparing with the other methods, furthermore, the nonlinear operation environments are frequently in practical application. Therefore, the performance of the proposed model can be checked by the two cases.

In this section, the proposed method is further tested by a rolling bearing under non-stationary conditions, such as linear time-varying, and nonlinear time-varying. The collected signals are from the Guilin University of Electronic Technology lab and the types of experimental bearing are ER-12K and ER-16K. The experiments were conducted on the machinery fault simulator test rig of SpectraQuest Co, which is shown in Fig. 3. Two accelerometers are installed on the rolling bearing in vertical and parallel directions, respectively.

The MFS-MG test rig.

In this subsection, the linear time-varying vibration signal is collected to testify to the performance of the proposed smooth method. The sampling frequency is set at 25.6 kHz and the sampled signal length is 12.8 s. To improve the computation efficiency, we select 153,600 samples as the tested signal. Meanwhile, the key-phase signal is recorded by a tachometer, and then the real IF is calculated for comparison, furthermore, due to the calculation method and other reasons, the real IF obtained is not smooth. The time–frequency representation conducted by the MSST method, as is shown in Fig. 4a, and the corresponding rough IF curve is displayed in Fig. 4b, the broken line is presented by magnifying the IF trajectory. Figure S5 presented the IF refinement results of L2-based and L1-based optimal functions, and the LSM method. The obtained results are presented in the Supplementary information document, namely "All computed compared results.pdf". In Fig. 5a, the green line is the estimated line and the blue line is the refined IF curve. The green line is surrounded by the green line and it's a smooth line. From Fig. 5b, the refined curve is close to the real IF and track the variable tendency. The compared results are displayed in Fig. S6, and the MAEs of the aforementioned methods are calculated as in Table S5, which could refer to the Supplementary information document, namely "All computed compared results.pdf".

Linear time-varying vibration signal. (a) Obtained time–frequency representation of the signal by using MSST, (b) corresponding coarse estimated IF.

Results of the refined IF. (a) The estimated and refined result, (b) the real and refined result.

The IF of the vibration signal is an important indicator for the condition monitoring of rotating machinery, especially in complex operating conditions. In this section, the real IF of the collected signal is up and down 2 Hz fluctuations, and the baseline is 38 Hz. The sampling frequency is 12.8 kHz and the signal length is 13.28 s. The generated by the MSST method is shown in Fig. 6a and it's the estimated IF is given in Fig. 6b. Fig. S7 shows the IF refinement results of the vibration signal by using L2-based and L1-based optimal functions, and the LSM method. The refined and real results are exhibited in Fig. S8. Both of them are presented in the Supplementary information document, namely "All computed compared results.pdf".

Non-linear time-varying vibration signal. (a) Obtained time–frequency representation of the signal by using MSST, (b) corresponding coarse estimated IF.

The green line is the estimated line and the blue line is the refined IF curve. In Fig. 7a, the broken line not only is smoothed but also infinitely close to the estimated line using the AWMM method. Similarly, the blue line is the real IF curve and the red line is the estimated IF curve. From Fig. 7b, the fitting effects are more accurate than the mentioned methods, which are provided by the proposed method. the MAEs of the aforementioned methods are calculated as in Table S6, which could refer to the Supplementary information document, namely "All computed compared results.pdf".

Results of the refined IF. (a) The estimated and refined result, (b) the real and refined result.

In this article, an adaptive weighted smooth model for smoothing ridge and improving estimation accuracy is developed. An adaptive weighted method is utilized to enhance the large energy value location of the estimated ridge. The regularization parameter is determined by the signal automatically. Meanwhile, the MM-based iterative method is employed to solve the construction convex model. Based on the estimated coarse time–frequency ridge by MSST computation, the ridge is smoothed to achieve high accuracy using AWMM. Thereafter, the index of MAE values is adopted to check the performance of the proposed method. The numerical and physical experiments are performed and the results show that the proposed method is more accurate than the commonly used polynomial curve fitting-based LSM method and L2-based norm regularization method. Moreover, the proposed method is superior to L1-based norm regularization with the same regularization parameter. Nevertheless, the proposed method has the main drawback to process fast time-varying signals. Future work can mainly consider developing the general refined method and expanding the proposed method application in multiple work conditions.

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

Wu, E. Q. et al. Fault diagnosis of rotating machinery using Gaussian process and EEMD-treelet. Mech. Mach. Theory. 174, 104919 (2022).

MATH Google Scholar

Xin, Y., Li, S. & Wang, J. A new fault feature extraction method for non-stationary signal based on advanced synchrosqueezing transform. J. Vib. Eng. Technol. 7, 291–299 (2019).

Article Google Scholar

Boashash, B. & Aissa-El-Bey, A. Robust multisensor time-frequency signal processing: A tutorial review with illustrations of performance enhancement in selected application areas. Digit. Signal Process. 77, 153–186 (2018).

Article MathSciNet Google Scholar

Park, C. et al. Time-frequency analysis of EEG asymmetry using bivariate empirical mode decomposition. IEEE Trans. Neural Syst. Rehabil. Eng. 19, 366–373 (2011).

Article PubMed Google Scholar

Cen, J., Yang, Z. H., Liu, X., Xiong, J. B. & Chen, H. H. A review of data-driven machinery fault diagnosis using machine learning algorithms. J. Vib. Eng. Technol. 9, 1–27 (2022).

Google Scholar

Malla, C. & Panigrahi, I. Review of condition monitoring of rolling element bearing using vibration analysis and other techniques. J. Vib. Eng. Technol. 7, 407–414 (2019).

Article Google Scholar

Feng, Z. P., Chen, X. W. & Wang, T. Y. Time-varying demodulation analysis for rolling element bearing fault diagnosis under variable speed conditions. J. Sound Vib. 400, 71–85 (2017).

Article ADS Google Scholar

Djurović, I. & Stanković, L. An algorithm for the Wigner distribution based instantaneous frequency estimation in a high noise environment. Signal Process. 84, 631–643 (2004).

Article MATH Google Scholar

Auger, F. et al. Time-frequecny reassignment and synchrosqueezing: An overview. IEEE Signal Process. 30, 32–41 (2013).

Article ADS Google Scholar

Auger, F. & Flandrin, P. Improving the readability of time-frequency and time-scale representations by the reassignment method. IEEE Trans. Signal Process. 43, 1068–1089 (1995).

Article ADS Google Scholar

Daubechies, I., Lu, J. F. & Wu, H.-T. Synchrosqueezed wavelet transforms: An empirical mode decomposition-like tool. Appl. Comput. Harmon. Anal. 3, 243–261 (2011).

Article MathSciNet MATH Google Scholar

Lv, Y., Ma, Y. B., Yuan, R. & Lv, S. T. Velocity synchronous chirplet extracting transform: An effective tool for fault diagnosis of variable-speed rotational machinery. IEEE Sens. J. 22, 13201–13211 (2022).

Article ADS Google Scholar

Thakur, G. & Wu, H.-T. Synchrosqueezing-based recovery of instantaneous frequency from nonuniform samples. SIAM J. Math. Anal. 43, 2078–2095 (2011).

Article MathSciNet MATH Google Scholar

Huang, Z., Zhang, J., Zhao, T. & Sun, Y. Synchrosqueezing S-transform and its application in seismic spectral decomposition. IEEE Trans. Geosci. Remote Sens. 54, 817–825 (2016).

Article ADS Google Scholar

Yang, Y., Peng, Z. K., Zhang, W. M. & Meng, G. Parametric time-frequency analysis methods and their engineering applications: A review of recent advances. Mech. Syst. Signal Process. 119, 182–221 (2019).

Article ADS Google Scholar

Ma, Y. B., Lv, Y., Yuan, R. & Ge, M. Synchro spline-kernelled chirplet extracting transform: A useful tool for characterizing time-varying features under noisy environments and applications to bearing fault diagnosis. Measurement 181, 109574 (2021).

Article Google Scholar

Gelman, L. & Petrunin, I. Time-frequency higher-order spectra with adjustment to the instantaneous frequency variation. Int. J. Adapt. Control Signal Process. 24, 178–187 (2008).

MATH Google Scholar

Ma, Y. B., Lv, Y., Yuan, R. & Song, G. B. Matching synchroextracting transform for mechanical fault diagnosis under variable-speed conditions. IEEE Trans. Instrum. Meas. 71, 3505612 (2022).

Google Scholar

Lv, S. T., Lv, Y., Yuan, R. & Li, H. High-order synchroextracting transform for characterizing signals with strong AM-FM features and its application in mechanical fault diagnosis. Mech. Syst. Signal Process. 172, 108959 (2022).

Article Google Scholar

Oberlin, T., Meignen, S. & Perrier, V. Second-order synchrosqueezing transform or invertible reassignment? Towards ideal time-frequency representations. IEEE Trans. Signal Process. 63, 1335–1344 (2015).

Article ADS MathSciNet MATH Google Scholar

Feng, Z. P. & Chen, X. W. Adaptive iterative generalized demodulation for nonstationary complex signal analysis: Principle and application in rotating machinery fault diagnosis. Mech. Syst. Signal Proces. 110, 1–27 (2018).

Article ADS Google Scholar

Yu, G., Wang, Z. H. & Zhao, P. Multisynchrosqueezing transform. IEEE Trans. Ind. Electron. 66, 5441–5455 (2018).

Article Google Scholar

Ren, M. J., Cheung, C. F. & Kong, L. B. A task specific uncertainty analysis method for least-squares-based form characterization of ultra-precision freedom surfaces. Meas. Sci. Technol. 23, 54005 (2012).

Article ADS Google Scholar

Pennecchi, F. & Malengo, A. A. Weighted total least squares algorithm for any fitting model with correlated variables. Metrologia 50, 654–662 (2013).

Article ADS Google Scholar

Poshusta, R. D. Spreadsheet error analysis of least-squares fitted models using Monte Carlo simulation. Comput. Phys. 5, 248–252 (1991).

Article ADS Google Scholar

Peng, D. K., Smith, W. A., Randall, R. B., Peng, Z. X. & Mechefske, C. K. Speed estimation in planetary gearboxes: A method for reducing impulsive noise. Mech. Syst. Signal Process. 159, 107786 (2021).

Article Google Scholar

Chen, S. Q., Peng, Z. K., Yang, Y., Dong, X. J. & Zhang, W. M. Intrinsic chirp component decomposition by using Fourier series representation. Signal Process. 137, 319–327 (2017).

Article Google Scholar

McNeill, S. Decomposing a signal into short-time narrow-banded modes. J. Sound Vib. 373, 325–339 (2016).

Article ADS Google Scholar

Tu, X. T., Swärd, J., Jakobsson, A. & Li, F. C. Estimating nonlinear chirp modes exploiting sparsity. Signal Process. 183, 107952 (2021).

Article Google Scholar

Download references

The authors are grateful to the support of National Natural Science Foundation of China (No. U1909217), the Zhejiang Natural Science Foundation of China (No. LD21E050001) and the Wenzhou Major Science and Technology Innovation Project of China (Nos. ZG2021019, ZG2021027).

College of Mechanical and Electrical Engineering, Wenzhou University, Wenzhou, 325035, People's Republic of China

Yi Liu & Jiawei Xiang

School of Mechanical and Electrical Engineering, Guilin University of Electronic Technology, Guilin, 541004, People's Republic of China

Yi Liu & Zhansi Jiang

School of Mathematics and Computer Science, Northwest Minzu University, Lanzhou, 730000, People's Republic of China

Hang Xiang

Pingyang Institute of Intelligent Manufacturing, Wenzhou University, Wenzhou, People's Republic of China

Jiawei Xiang

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Y.L.: methodology, investigation, writing—original draft. H.X.: formal analysis, validation. Z.J.: data curation, investigation. J.X.: supervision, project administration, investigation, resources, writing—review and editing.

Correspondence to Jiawei Xiang.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Liu, Y., Xiang, H., Jiang, Z. et al. Refining the time–frequency characteristic of non-stationary signal for improving time–frequency representation under variable speeds. Sci Rep 13, 5215 (2023). https://doi.org/10.1038/s41598-023-32333-w

Download citation

Received: 16 December 2022

Accepted: 26 March 2023

Published: 30 March 2023

DOI: https://doi.org/10.1038/s41598-023-32333-w

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.