# Bayesian operational modal analysis

Bayesian operational modal analysis (BAYOMA) adopts a Bayesian system identification approach for operational modal analysis (OMA). Operational modal analysis aims at identifying the modal properties (natural frequencies, damping ratios, mode shapes, etc.) of a constructed structure using only its (output) vibration response (e.g., velocity, acceleration) measured under operating conditions. The (input) excitations to the structure are not measured but are assumed to be 'ambient' ('broadband random'). In a Bayesian context, the set of modal parameters are viewed as uncertain parameters or random variables whose probability distribution is updated from the prior distribution (before data) to the posterior distribution (after data). The peak(s) of the posterior distribution represents the most probable value(s) (MPV) suggested by the data, while the spread of the distribution around the MPV reflects the remaining uncertainty of the parameters.

## Pros and cons

In the absence of (input) loading information, the identified modal properties from OMA often have significantly larger uncertainty (or variability) than their counterparts identified using free vibration or forced vibration (known input) tests. Quantifying and calculating the identification uncertainty of the modal parameters become relevant.

The advantage of a Bayesian approach for OMA is that it provides a fundamental means via the Bayes' Theorem to process the information in the data for making statistical inference on the modal properties in a manner consistent with modeling assumptions and probability logic.

The potential disadvantage of Bayesian approach is that the theoretical formulation can be more involved and less intuitive than their non-Bayesian counterparts. Algorithms are needed for efficient computation of the statistics (e.g., mean and variance) of the modal parameters from the posterior distribution. Unlike non-Bayesian methods, the algorithms are often implicit and iterative. E.g., optimization algorithms may be involved in the determination of most probable value, which may not converge for poor quality data.

## Methods

Bayesian formulations have been developed for OMA in the time domain[1] and in the frequency domain using the spectral density matrix[2] and fast Fourier transform (FFT)[3] of ambient vibration data. Based on the formulation for FFT data, fast algorithms have been developed for computing the posterior statistics of modal parameters.[4] Recent developments based on EM algorithm[5] show promise for simpler algorithms and reduced coding effort. The fundamental precision limit of OMA has been investigated and presented as a set of uncertainty laws which can be used for planning ambient vibration tests.[6]

## Connection with maximum likelihood method

Bayesian method and maximum likelihood method (non-Bayesian) are based on different philosophical perspectives but they are mathematically connected; see, e.g.,[7] and Section 9.6 of.[4] For example,

• Assuming a uniform prior, the most probable value (MPV) of parameters in a Bayesian method is equal to the location where the likelihood function is maximized, which is the estimate in Maximum Likelihood Method
• Under a Gaussian approximation of the posterior distribution of parameters, their covariance matrix is equal to the inverse of Hessian of the negative log of likelihood function at the MPV. Generally, this covariance depends on data. However, if one assumes (hypothetically; non-Bayesian) that the data is indeed distributed as the likelihood function, then for large data size it can be shown that the covariance matrix is asymptotically equal to the inverse of the Fisher information matrix (FIM) of parameters (which has a non-Bayesian origin). This coincides with the Cramer–Rao bound in classical statistics, which gives the lower bound (in the sense of matrix inequality) of the ensemble variance of any unbiased estimator. Such lower bound can be reached by maximum-likelihood estimator for large data size.
• In the above context, for large data size the asymptotic covariance matrix of modal parameters depends on the 'true' parameter values (a non-Bayesian concept), often in an implicit manner. It turns out that by applying further assumptions such as small damping and high signal-to-noise ratio, the covariance matrix has mathematically manageable asymptotic form, which provides insights on the achievable precision limit of OMA and can be used to guide ambient vibration test planning. This is collectively referred as 'uncertainty law'.[6]

## Notes

• See monographs on non-Bayesian OMA [8][9][10] and Bayesian OMA [4]
• See OMA datasets [11]
• See Jaynes[12] and Cox[13] for Bayesian inference in general.
• See Beck[14] for Bayesian inference in structural dynamics (relevant for OMA)
• The uncertainty of the modal parameters in OMA can also be quantified and calculated in a non-Bayesian manner. See Pintelon et al.[15]

## References

1. ^ Yuen, K.V.; Katafygiotis, L.S. (2001). "Bayesian time-domain approach for modal updating using ambient data". Probabilistic Engineering Mechanics. 16 (3): 219–231. doi:10.1016/S0266-8920(01)00004-2.
2. ^ Yuen, K.V.; Katafygiotis, L.S. (2001). "Bayesian spectral density approach for modal updating using ambient data". Earthquake Engineering & Structural Dynamics. 30 (8): 1103–1123. doi:10.1002/eqe.53.
3. ^ Yuen, K.V.; Katafygiotis, L.S. (2003). "Bayesian Fast Fourier Transform approach for modal updating using ambient data". Advances in Structural Engineering. 6 (2): 81–95. doi:10.1260/136943303769013183.
4. ^ a b c Au, S.K. (2017). Operational Modal Analysis: Modeling, Inference, Uncertainty Laws. Springer.
5. ^ Li, B.; Au, S.K. (2019). "An expectation-maximization algorithm for Bayesian operational modal analysis with multiple (possibly close) modes". Mechanical Systems and Signal Processing. doi:10.1016/j.ymssp.2019.06.036.
6. ^ a b Au, S.K.; Brownjohn, J.M.W.; Mottershead, J. (2018). "Quantifying and managing uncertainty in operational modal analysis". Mechanical Systems and Signal Processing. doi:10.1016/j.ymssp.2017.09.017. hdl:10871/30384.
7. ^ Au, S.K.; Li, B. (2017). "Posterior uncertainty, asymptotic law and Cramér‐Rao bound". Mechanical Systems and Signal Processing. doi:10.1002/stc.2113.
8. ^ Van Overschee, P.; De Moor, B. (1996). Subspace Identification for Linear Systems. Boston: Kluwer Academic Publisher.
9. ^ Schipfors, M.; Fabbrocino, G. (2014). Operational Modal Analysis of Civil Engineering Structures. Springer.
10. ^ Brincker, R.; Ventura, C. (2015). Introduction to Operational Modal Analysis. John Wiley & Sons.
11. ^
12. ^ Jaynes, E.T. (2003). Probability Theory: The Logic of Science. United Kingdom: Cambridge University Press.
13. ^ Cox, R.T. (1961). The Algebra of Probable Inference. Baltimore: Johns Hopkins University Press.
14. ^ Beck, J.L. (2010). "Bayesian system identification based on probability logic". Structural Control and Health Monitoring. 17 (7): 825–847. doi:10.1002/stc.424.
15. ^ Pintelon, R.; Guillaume, P.; Schoukens, J. (2007). "Uncertainty calculation in (operational) modal analysis". Mechanical Systems and Signal Processing. 21 (6): 2359–2373. Bibcode:2007MSSP...21.2359P. doi:10.1016/j.ymssp.2006.11.007.