EM-aided fast posterior covariance computation in Bayesian FFT method

Wei Zhu,Binbin Li
DOI: https://doi.org/10.1016/j.ymssp.2024.111211
IF: 8.4
2024-04-01
Mechanical Systems and Signal Processing
Abstract:In operational modal analysis, the Bayesian FFT method approximates the posterior distribution of modal parameters with a multivariate normal distribution, whose covariance is given as the inverse of the Hessian matrix of the negative log-likelihood function (NLLF). The current algorithm derives the Hessian matrix in an element-by-element way, which is tedious to computer-coding and inefficient in challenging cases, e.g., the closely-spaced modes. An expectation–maximization (EM) based approach is developed in this paper to provide a competitive alternative in terms of computational efficiency and mathematical elegancy. It takes advantage of the relationship between the NLLF and the complete-data NLLF, expressing the Hessian matrix of NLLF in terms of the gradient and Hessian matrix of the complete-data NLLF. Since the complete-data NLLF has a simpler structure of modal parameters, the derivation of its gradient and Hessian matrix is easier, especially with the usage of complex derivatives and matrix calculus. A nullspace projection method is also proposed to overcome the singularity of the Hessian matrix due to mode shape norm constraints. Synthetic, laboratory, and field data are applied to verify the performance of the proposed method. It shows that the proposed EM method yields identical results as the current algorithm, with a speedup of tens or even hundreds. The highly efficient evaluation of the Hessian matrix offers a further opportunity to accelerate the Bayesian FFT method, e.g., using the Newton method.
engineering, mechanical
What problem does this paper attempt to address?