Next Article in Journal
Research on the Heating Effect of a Convection Radiator Based on a Human Thermophysiological Model
Previous Article in Journal
Assessment of Soft-First-Floor Structures Reinforced by Rocking Frame Based on Seismic Resilience
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Accuracy of Frequency Estimation of Structure Vibration under Ambient Excitation: Problems, Causes, and Methods

1
State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Nanjing Hydraulic Research Institute, Nanjing 210029, China
2
College of Mechanics and Materials, Hohai University, Nanjing 210098, China
*
Authors to whom correspondence should be addressed.
Buildings 2024, 14(1), 198; https://doi.org/10.3390/buildings14010198
Submission received: 6 November 2023 / Revised: 27 December 2023 / Accepted: 10 January 2024 / Published: 12 January 2024
(This article belongs to the Section Building Structures)

Abstract

:
Accurate identification of building structure frequencies forms the basis for damage detection. The structural dynamic response signal, under ambient excitation, can be transformed into a superposition of multiple single-frequency exponentially damped sinusoids combined with random white noise. However, the peak power spectrum of the response signal tends to exhibit line splitting, compromising the precision of frequency identification. This study examines the accuracy characteristics of the single-frequency free damping vibration signal (SFFDVS) and derives the Cramer–Rao lower bound for the frequency estimator. It thoroughly analyzes the factors influencing the accuracy of SFFDVS frequency identification. The study reveals that the primary cause of spectral line splitting is the random delay inherent in SFFDVS. Based on the maximum likelihood method (MLM), this research introduces the MLM algorithm for SFFDVS and provides a simulation analysis. The findings indicate that the MLM estimation algorithm for frequency parameters effectively addresses spectral line splitting and offers robust noise resistance and recognition accuracy.

1. Introduction

Modal parameters of building structures under environmental excitation for damage detection are widely used in the fields of structural health monitoring (SHM) and post-earthquake evaluation of buildings [1,2,3]. For example, B Bhowmik et al. [4,5] proposed a novel approach of recursive canonical correlation analysis (RCCA) for robust damage detection. Application of the RCCA to the combined ambient and earthquake responses obtained from the UCLA Factor Building demonstrates the robustness of this methodology as an ideal candidate for real-time SHM.
Frequency is the most fundamental parameter of structural modes. The primary step in completing structural damage detection and SHM is to accurately identify the natural frequency of the structure as much as possible. The traditional method for identifying structural natural frequency requires both the measured input excitation and the responses [6,7]. Civil structures invariably exist within specific vibrational environments, and ambient vibration tests can determine modal parameters using only the output ambient response [8,9,10]. When a structure is subjected to ambient excitation and this excitation is presumed to be white noise, the system’s natural frequency can also be determined based on the structure’s response signals. Frequency is a crucial parameter in safety monitoring and damage identification of civil structures [11,12,13]. Numerous damage-identification methods that rely on frequency parameters have been explored [14,15,16,17,18,19]. The related research especially sheds light on the development of a new class of damage detection strategies that identify the damage in real time, using the key principles of the first-order eigen perturbation (FOEP) technique [4,20]. Frequent parameter changes are a crucial indicator for assessing structural performance in practical engineering applications. The methods of damage identification demand high accuracy in recognizing frequency parameters [21]. When applied to environmental excitation, the engineering community generally believes that the modal parameter method can precisely determine the system frequency [22,23,24]. However, accuracy in recognizing frequency parameters is often intimately linked to response signals and recognition algorithms. There is a dearth of comprehensive studies on accuracy in recognizing frequency parameters under environmental excitation.
Yang and Cao [25,26] extensively examined the frequency parameter identification algorithm for noisy complex (and real) sinusoidal signals. They introduced a high-precision frequency algorithm for such signals and conducted an error analysis. Theoretical derivations confirmed that the mean square error (MSE) of frequency parameter identification for a sinusoidal signal can approach the Cramer–Rao lower bound (CRLB) [27]. The structural response signal under environmental excitation differs significantly from the noisy sinusoidal signal model used in the communications industry. The accuracy of frequency parameter identification is influenced by noise and the response phase difference resulting from the randomness of environmental excitation, damping parameters, and spectral leakage. Aboutanios and Bernard [28] performed a Fourier analysis on complex periodic signals to estimate rough frequency. Subsequently, the identified frequency error approached the CRLB through binary interpolation of Fourier coefficients. Umesh and Tufts [29] introduced a fast maximum likelihood (FML) estimation for the parameters of multiple exponentially damped sinusoids. They proposed an iterative procedure that separates the data into its constituent components and estimates the parameters of each component individually. YU [30] suggested a state-space model method to estimate the damping factor and frequency of multi-component exponentially damped sinusoids. Simulations demonstrated that the state-space model algorithm outperformed the ESPRIT algorithm, especially at low signal-to-noise ratios (SNR).
In light of the characteristics of the structural response signal under ambient excitation, this study develops a model for the structural response signal under such excitation. The structural response can be divided into a free-damping vibration signal and a random response signal. This research deduces the CRLB for frequency parameter identification of free damping vibration signals (FDVS) and analyzes the error factors involved in frequency parameter identification. A maximum likelihood estimation algorithm for identifying frequency parameters of freely attenuated vibration signals is introduced, followed by a simulation analysis.

2. Materials and Methods

2.1. Target Model and Its Problems

During the experiment on modal testing using ambient excitation, the input for civil structures, resulting from foundation vibrations or wind, was approximated as either white noise or an ergodic process [31,32,33]. The structure’s natural frequency can be identified based on its responses. Generally, the structural response to random excitation comprises a deterministic part (pulse or step) and a random part (assumed to have a zero mean). The dynamic equation for a single-degree-of-freedom system can be mathematically expressed as follows:
m x ¨ + c x ˙ + k x = f ( t )
where f t is a sample of the excitation random process; m is the system’s mass; k is the system’s stiffness; and c is the system’s damping coefficient.
Applying the Laplace transform to the equation above yields:
m s 2 + c s + k X s = F s + m x ˙ 0 + m s + c x 0
X s = F s m s 2 + c s + k + m m s 2 + c s + k x ˙ 0 + m s + c m s 2 + c s + k x 0
Applying the inverse Laplace transformation to the previous equation results in:
x t = 0 t F τ m h t τ d τ + x ˙ 0 h t + x 0 g t = x 1 t + x 2 t + x 3 t
where h t = 1 ω 0 e ξ ω 0 t sin 1 ξ 2 ω 0 t
g t = 1 1 ξ 2 e ξ ω 0 t cos 1 ξ 2 ω 0 t φ
φ = t g 1 ξ 1 ξ 2
The mathematical model for frequency identification is described as the combination of a free-damping vibration signal and a random white noise signal, defined as follows:
s t = i = 1 l A i e β ( t t i ) cos ω 0 t + φ i + z ( t ) t 0
where z ( t ) is a white noise signal; l is the pulse or step excitation times in random excitation; A i is the maximum amplitude caused by pulse or step excitation; t i is the time of pulse or step excitation; φ i is the phase of a response signal; and A i , t i are both random signals.
Equation (5) is sampled discretely. If the sampling frequency is f s , the sampling time length is T, the sampling length is N, and T = N / f s , then
s n = i = 1 l A i e β ( n f s t i ) cos ω 0 n f s + φ i + z n n = 0,1 , N
Based on Equation (6), the random excitation response signal for a single-degree-of-freedom system over 10 s was simulated, as depicted in Figure 1.
The derivations above indicate that the structural response signal under environmental excitation combines multiple single−frequency free attenuation vibration signals with random white noise signals. The random decrement method can isolate the system’s free attenuation vibration signals through averaging and mathematical statistics. Wavelet packet decomposition technology can divide each order of free attenuation vibration signals into single-frequency free attenuation vibration signals. Consequently, identifying the frequency of response signals under environmental excitation can be transformed into identifying the frequency parameter of single-frequency free attenuation vibration signals. However, due to the unpredictability of actual structural environmental excitations and environmental noise, there are challenges in accurately identifying the frequency of the response signal spectrum. The CRLB of the frequency parameter identification error’s mean square value for the response signal model will be derived in the subsequent sections to examine the factors leading to frequency parameter identification errors.

2.2. Theoretical Analysis and Its Causes

2.2.1. Single-Frequency Free Damping Vibration Signal

SFFDVS is described as follows:
s t = a e ξ ω 0 t cos 1 ξ 2 ω 0 t + φ ( 0 t T )
where a,   ξ , ω 0 , and φ are the amplitude, damping ratio, frequency, and phase of the signal, respectively.
When t < 0 , s t = 0 . Energy integration is performed on s(t), and the integration formula is as follows:
E = lim T T T s t 2 d t = lim T 0 T s t 2 d t = a 2 1 + ξ 2 4 ξ ω 0
Equation (8) indicates that the single-frequency free attenuation signal possesses limited energy. When SFFDVS combines with the white noise signal, the SNR across the entire domain becomes zero due to the infinite energy of the white noise signal. The genuine free attenuation signal is a finite sequence signal. Its SNR is the free attenuation signal energy ratio to the noise signal energy over the observation time, T.
Equation (7) is discretely sampled, and N sample values are obtained, forming the discrete sequence.
s n = a e ξ ω 0 n f s cos 1 ξ 2 ω 0 n f s + φ n = 0,1 , 2 , , N 1
where f s is the sampling frequency.
It is assumed that Gaussian white noise z(n) is mixed into the signal s(n); if x(n) = s(n) + z(n), the signal sequence is
x n = s ( n ) + z n n = 0,1 , 2 , , N 1
where z ( n ) is a Gaussian noise sample with zero mean and variance σ 2 .
The energy integration of s n is
E n = n = 0 N 1 a e ξ ω 0 n f s cos 1 ξ 2 ω 0 n f s 2 n = 0,1 , 2 , , N 1
Expanding the above formula yields
E n = n = 0 N 1 a e ξ ω 0 n f s cos 1 ξ 2 ω 0 n f s 2 a 2 2 n = 0 N 1 e 2 ξ ω 0 n f s a 2 2 1 e 2 ξ ω 0 N f s 1 e 2 ξ ω 0 1 f s
When the sampling length is long enough, then e 2 ξ ω 0 N f s 0 . When the system is a small damping system, based on the sampling theorem f s 2.56 ω 0 2 π , 2 ξ ω 0 1 f s 0 . If it defines 2 ξ ω 0 1 f s = x , then
E n a 2 2 1 1 e x
Based on Taylor’s expansion, e x = 1 x + 1 2 ! x 2 1 3 ! x 3 + o x 3 1 x , and by substituting into Equation (13), it can be derived that
E n a 2 2 x = a 2 f s 4 ξ ω 0
The SNR of a noisy single-degree-of-freedom free attenuation signal can be obtained as follows:
ς = a 2 f s 4 N ξ ω 0 σ 2
If the sampling time is T = N / f s , then
ς = a 2 4 T ξ ω 0 σ 2
Based on the analysis, the free attenuation vibration signal exhibits the following characteristics:
The free attenuation vibration signal is energy-limited; as t approaches infinity, a ~ approaches zero. As the sampling time extends, the SNR decreases; when the sample time, T, approaches infinity, the SNR, ς ~ , approaches zero.
The SNR of the free attenuation signal is directly proportional to the sampling frequency and the square of the initial amplitude. It is inversely proportional to the damping ratio, signal frequency, and sampling length.
The primary parameters of a single-frequency free attenuation signal encompass the initial amplitude, damping coefficient, frequency, and phase angle. The frequency and damping coefficient are the main parameters of interest in practical engineering applications.

2.2.2. CRLB of Single-Frequency Free Damping Vibration Signal

The parameter estimation model of SFFDVS is given by Equation (9). For convenience of calculation, normalization is conducted. This study considers β = ξ ω 0 / f s and ω 1 = 1 ξ 2 ω 0 / f s ; thus,
x n = a e β n cos ω 1 n + φ + z ( n )
where a , ω 1 , β , a n d   φ are unknown, and θ = a , ω 1 , β , φ T .
From the parameter estimation model, the following can be derived:
z n = x n s ( n )   n = 0 ,   1 ,   2 , , N 1
The probability density function is presented as follows:
f z n , θ = 1 2 π σ e x n s ( n ) 2 σ 2 2
The joint probability density function is described as follows:
f x , θ = i = 0 N 1 f z i , θ = 1 2 π σ 2 N 2 e x p 1 2 σ 2 n = 0 N 1 x n s ( n ) 2
and
L x , θ = l n f x , θ = N 2 ln 2 π σ 2 1 2 σ 2 n = 0 N 1 x n s ( n ) 2
The unbiased CR bounds consist of the diagonal elements of the inverse of the Fisher information matrix, J. When z ( n ) represents a Gaussian noise sample with a zero mean and a variance of σ 2 , the Fisher information matrix is as follows [34]:
J i , j = 1 σ 2 n = 0 N 1 s ( n ) θ i · s ( n ) θ j
The following can be derived based on Equation (22):
J 1,1 = 1 σ 2 n = 0 N 1 s ( n ) θ 1 · s ( n ) θ 1 = 1 4 σ 2 e i 2 φ 1 x a + e i 2 φ 1 x b + 2 1 e 2 β
where x a = e 2 β + i 2 ω 1 , x b = e 2 β i 2 ω 1
J 1,2 = 1 σ 2 n = 0 N 1 s ( n ) θ 1 · s ( n ) θ 2 i a 4 σ 2 e i 2 φ x a 1 x a 2 e i 2 φ x b 1 x b 2
J 1,3 = 1 σ 2 n = 0 N 1 s ( n ) θ 1 · s ( n ) θ 3 a e i 2 φ x b 4 σ 2 1 x b 2 a e i 2 φ x a 4 σ 2 1 x a 2 a 1 + 2 β 8 σ 2 β 2
J 1,4 = 1 σ 2 n = 0 N 1 s ( n ) θ 1 · s ( n ) θ 4 = a 2 σ 2 n = 0 N 1 e 2 β n sin 2 ω 1 n + 2 φ a i 4 σ 2 e i 2 φ 1 x b e i 2 φ 1 x a
J 2,2 = 1 σ 2 n = 0 N 1 s ( n ) θ 2 · s ( n ) θ 2 a 2 e 2 β 1 + e 2 β 2 σ 2 e 2 β 1 3 a 2 8 σ 2 e i 2 φ x a 1 + x a 1 x a 3 + e i 2 φ x b 1 + x b 1 x b 3
J 2,3 = 1 σ 2 n = 0 N 1 s ( n ) θ 2 · s ( n ) θ 3 i a 2 4 σ 2 e i 2 φ x b 1 + x b 1 x b 3 e i 2 φ x a 1 + x a 1 x a 3
J 2,4 = 1 σ 2 n = 0 N 1 s ( n ) θ 2 · s ( n ) θ 4 = a 2 e 2 β 2 σ 2 e 2 β 1 2 a 2 e i 2 φ x b 4 σ 2 1 x b 2 a 2 e i 2 φ x a 4 σ 2 1 x a 2
J 3,3 = 1 σ 2 n = 0 N 1 s n θ 3 · s n θ 3 a 2 e 2 β 1 + e 2 β 2 σ 2 e 2 β 1 3 + a 2 4 σ 2 e i 2 φ x a 1 + x a 1 x a 3 + e i 2 φ x b 1 + x b 1 x b 3
J 3,4 = 1 σ 2 n = 0 N 1 s n θ 3 · s n θ 4 i a 2 4 σ 2 e i 2 φ x a 1 x a 2 e i 2 φ x b 1 x b 2
J 4,4 = 1 σ 2 n = 0 N 1 s ( n ) θ 4 · s ( n ) θ 4 a 2 4 σ 2 β
The Fisher information matrix is symmetric; thus:
J J 1,1 J 1,2 J 1,2 J 2,2 J 1,3 J 1,4 J 2,3 J 2,4 J 1,3 J 2,3 J 1,4 J 2,4 J 3,3 J 3,4 J 3,4 J 4,4
Based on Fisher’s matrix, the CRLB of the parameters can be estimated.
E ω ^ 1 ω 1 2 = E θ ^ 2 θ 2 2 J 1 θ 22
E β ^ β 2 = E θ ^ 3 θ 3 2 J 1 θ 33
After normalization processing with ω 0 = ω 1 × f s / 1 ξ 2 , the following is derived:
E ω ^ 0 ω 0 2 = f s 2 · E ω ^ 1 ω 1 2 / 1 ξ 2 f s 2 J 1 ( θ ) 22 1 ξ 2
Hence, obtaining the analytic formula for the inverse matrix of the Fisher matrix proved challenging. A simulation analysis was conducted to better understand the identification accuracy of the parameters under consideration.

2.2.3. Characteristics of CRLB of Single-Frequency Free Damping Vibration Signal

Four simulation groups were established to examine the effects of sampling frequency, phase angle, damping coefficient, and noise on recognition accuracy, each assessing the impact of distinct parameters.
(1)
Influence of sample frequency on frequency identification accuracy 
Assuming that a = 1.0 ,   β = 0.05 ω 1 f s , ω 1 = 62.8 f s ( f = 10   H z ) , φ = 0 , the noise signal variance σ = 0.2 , and the data sampling length N = 20,001 , then f s 30   H z ,   500   H z . The Fisher matrix was determined, and subsequently, the CRLB of the frequency was estimated using Equation (36) by substituting these parameters into Equations (23)–(33). The computed results are listed in Table 1, with a corresponding visual representation in Figure 2.
The CRLB of a real sinusoidal signal is
E ω ^ 0 ω 0 2 12 ζ N ( N 2 1 )
where ζ is the SNR ( ζ = a 2 2 σ 2 ); N is the length of the data.
The CRLB is the result of normalized sampling frequency, which can be obtained by considering the influence of sampling frequency in the same way as Equation (36).
E ω ^ 01 ω 01 2 12 f s 2 ζ N ( N 2 1 )
The free attenuation signal differs significantly from the sinusoidal signal. When the data length of the sinusoidal signal remains constant, and the sampling theorem is met, a lower sampling frequency results in a smaller frequency identification variance. However, simulation results for free attenuation signals indicate that free fading signals’ frequency and time domain resolutions align. As the time domain resolution is enhanced, so is the frequency domain resolution. The simulation results reveal that:
(1)
Elevating the sampling frequency can decrease the lower limit of frequency identification variance. Enhancing the sampling frequency improves the frequency identification accuracy at low sampling frequencies.
(2)
Continuously increasing the sampling frequency does not indefinitely enhance recognition accuracy. When the sampling frequency is 40 times the recognition frequency, further increases in the sampling frequency minimally impact recognition accuracy.
(3)
Free-damping vibration signals differ markedly from sinusoidal signals. A lower sampling frequency yields higher frequency identification accuracy with consistent data length. A higher sampling frequency increases frequency identification accuracy with a sufficiently long signal length. Thus, a higher sampling frequency can be employed for signal acquisition to boost frequency identification accuracy in actual modal tests.
(4)
According to the influence of the sampling frequency of the free-damping vibration signals on the accuracy of frequency identification, a real-time eigen perturbation strategy can be further adopted to adaptively adjust the sampling frequency in real-time modal testing, thereby improving the accuracy of frequency identification [5].
(2)
Influence of damping coefficient on frequency identification accuracy 
Assuming that a = 1.0 ,   β = ξ ω 1 f s , ω 1 = 62.8 f s ( f = 10   H z ) , φ = 0 , noise signal variance σ = 0.2 , data sampling length N = 20,001 , and f s = 500   H z , then ξ 0.01 ,   0.1 . The Fisher matrix was computed, and subsequently, the CRLB of the frequency was determined using Equation (36) by substituting these parameters into Equations (23)–(33). Table 2 details the outcomes, whereas Figure 3 represents them.
The simulation findings show that the frequency identification error also rises as the signal damping coefficient increases. Hence, in systems with a larger damping coefficient, the impact of this coefficient must be accounted for during system frequency identification.
(3)
Influence of noise on frequency identification accuracy 
If a = 1.0 ,   β = 0.05 ω 1 f s , ω 1 = 62.8 f s ( f = 10   H z ) , φ = 0 , data sampling length N = 20,001 , and f s = 500   H z , then   σ 0.01 ,   0.2 . The Fisher matrix was determined, and then the CRLB of the frequency was estimated using Equation (36) by substituting these parameters into Equations (23)–(33). The calculation results are listed in Table 3, and the graph is shown in Figure 4.
The simulation results indicate that the CRLB of frequency parameter identification is proportional to a 2 σ 2 . As noise increases, the accuracy of frequency identification decreases. Consequently, the initial step in frequency parameter identification should involve signal denoising. Enhancing the SNR can notably enhance the frequency identification accuracy.
(4)
Influence of the phase of the signal on frequency identification accuracy 
Assuming that a = 1.0 ,   β = 0.05 ω 1 f s , ω 1 = 62.8 f s ( f = 10   H z ) , data sampling length N = 20,001 , f s = 500   H z , and   σ = 0.2 , then   φ 0 ,   360 . By substituting the above parameters into Equations (23)–(33), the Fisher matrix was calculated, and then the CRLB of the frequency was estimated from Equation (36). The calculation results are listed in Table 4 and illustrated in Figure 5.
The simulation results indicate the following: (1) The phase angle significantly affects the lower variance of frequency parameter identification. The variance is lowest when the phase angle is approximately 54° and 234° and highest at approximately 144° and 324°. (2) The difference between the maximum and minimum values of the CRLB is 0.0001 for various phase angles.
Since the phase angle’s influence is less pronounced than the noise, sampling frequency, and damping coefficient, the randomness of excitation in actual sampled signals renders the phase angle uncertain. Consequently, the superposition of free-damping signals with identical frequencies but differing phase angles can lead to the spectral line-splitting phenomenon. This uncertainty can adversely affect the accuracy of frequency identification.

2.3. Algorithm Derivation and Its Methods

Maximum likelihood estimation seeks to estimate parameters based on known observations, without prior knowledge of the unknown quantity under consideration. In applying the maximum likelihood estimation method, the parameters to be estimated are considered constant, though they remain unknown, while the known observation data are treated as random sequences. Let x = ( x 1 , x 2 , , x N ) represent N observations of random variables and let f x | θ , θ Θ represent the joint conditional probability density function of observed samples with specific parameters, where θ is the potential range of values. The joint conditional probability density function is as follows:
f θ = f x | θ i = 1 n f x i | θ
Assuming that the joint conditional probability density function exists and is constrained, the maximum likelihood estimation is the estimation value θ ^ that optimizes the likelihood function f x | θ , θ Θ . Given a sampling frequency of 1 Hz, the probability density function for a single-frequency free damping vibration signal is presented as follows:
f x , θ = 1 2 π σ 2 N 2 e x p 1 2 σ 2 n = 0 N 1 x n a e β n c o s ω 1 n + φ 2
where β = ξ ω 0 ,   ω 1 = 1 ξ 2 ω 0 .
When the identification parameter θ ^ results in f x , θ achieving its maximum value, this parameter represents the maximum likelihood estimation. Taking the logarithm of f(x,θ) yields the following expression:
ln f x , θ = N 2 l n 2 π σ 2 1 2 σ 2 n = 0 N 1 x n a e β n c o s ω 1 n + φ 2
and
L x , θ = n = 0 N 1 x n a e β n c o s ω 1 n + φ 2
Expanding Equation (42) yields
L x , θ = n = 0 N 1 x 2 ( n ) + a 2 e 2 β n c o s 2 ω 1 n + φ 2 a x n e β n c o s ω 1 n + φ
By ignoring the known time series, n = 0 N 1 x 2 ( n ) , from Equation (43), we can derive
L 0 x , θ = a 2 4 β 2 a n = 0 N 1 x n e β n c o s ω 1 n + φ
When a and β are unknown values, then
L 0 x , ω = n = 0 N 1 x n e β n c o s ω n + φ
When the parameter ω ^ causes f x , ω 1 to achieve its maximum, ln f x , ω ^ reaches its maximum, whereas L x , θ and L 0 x , θ achieve their minima and L 0 x , ω 1 attains its maximum. Maximum likelihood estimation is equivalent to the best square approximation estimation.
Based on Equation (45), L 0 x , ω has the greatest value when identifying the parameter ω ^ ω 1 , where ω ^ is the maximum likelihood estimation parameter.
Taking the partial derivative of L 0 x , ω , we can obtain
L x x , ω = L 0 x , θ ω = 2 a n = 0 N 1 x n n e β n s i n ω n + φ
When Equation (40) attains the maximum value, L x x , ω = 0 . When ω ^ ω 1 , we can obtain
n = 0 N 1 x n n e β n s i n ω ^ n + φ = 0  
Likewise, if a , β , a n d   φ in the parameter   θ = a , ω 1 , β , φ T are deduced similarly to the parameter ω 1 , then
L x x , a = L 0 x , θ a = a 2 β 2 n = 0 N 1 x n e β n c o s ω 1 n + φ = 0
When a ^ a , then
a ^ = 4 β n = 0 N 1 x n e β n c o s ω 1 n + φ
When β ^ β , the following can be derived:
n = 0 N 1 x n n e β ^ n c o s ω 1 n + φ a 8 β ^ 2
The following can be obtained from Equations (49) and (50):
β ^ = n = 0 N 1 x n e β n c o s ω ^ n + φ 2 n = 0 N 1 x n n e β n c o s ω ^ n + φ
When φ ^ φ , then
φ ^ = tan 1 n = 0 N 1 x n e β n s i n ω ^ n n = 0 N 1 x n e β n c o s ω ^ n
Based on the above deduction, the frequency identification of a free-damping vibration signal can be categorized into the following steps:
(1)
The frequency and damping coefficient parameters are initially estimated using the self-power spectrum;
(2)
The initial estimated frequency and damping coefficient parameters are substituted into Equation (52), and the corresponding calculations is performed to obtain phase angle parameters;
(3)
The initial estimated damping coefficient and phase angle estimation parameter are substituted into Equation (45), and the new estimated frequency corresponds to its maximum value;
(4)
The iteration is repeated from the second to the third step to continuously improve recognition accuracy.
Figure 6 depicts the flowchart of frequency parameter identification using the maximum likelihood method.

3. Results

3.1. Simulation Signal Model

For the simulation calculation model (Equation (10)), assuming that a = 1.0 ,   β = 0.05 ω 1 f s , ω 1 = 62.8 f s ( f = 10   H z ) , φ = 0 , data sampling length N = 20,001 , f s = 500   H z , and σ 0.01 ,   0.2 , the SNR ranges between 30 and −30 dB. Figure 7 illustrates the time domain curves at different SNRs.

3.2. Error Analysis

A rough estimation of the damping coefficient and frequency is required during the estimation process using the maximum likelihood algorithm. The rough estimation values are simulated within a specific error range to evaluate the impact of the accuracy of these rough estimates on frequency estimation. Assuming that the error in the rough estimate of the frequency ranges from −20% to 20%, meaning the rough estimate frequency lies between 8 and 12 Hz, the effect on recognition results is presented in Table 5. Based on the simulation findings, when the rough frequency estimate ranges from −20% to 20%, the recognition error can be diminished to less than 0.26% after four iterations. Consequently, the maximum likelihood frequency estimation algorithm is not significantly influenced by the initial estimate and can rapidly converge to the actual frequency value through numerous iterations.
Similarly, the impact of the initial estimate of the damping coefficient on frequency identification accuracy is also simulated. Assuming that the initial error in estimating the damping coefficient ranges from −50% to 50%, meaning the damping coefficient varies between 0.025 and 0.075, the maximum likelihood method estimates the frequency for different damping coefficients. These simulation results are documented in Table 6. When the damping coefficient ranges from −50% to 50%, the most considerable discrepancy in identification frequency is only 0.034 Hz, and the maximum error in frequency identification is 0.34%. Hence, in the practical algorithm, the initial estimate of the damping coefficient is determined using the power spectrum half-power bandwidth method, and the damping coefficient does not require iterative calculations. The maximum likelihood frequency estimation algorithm is not sensitive to changes in the damping coefficient.
The identification errors in the initial amplitude A and damping coefficient ξ of a free damping signal exert minimal influence on frequency parameter identification. However, the phase angle significantly impacts frequency parameter identification. Therefore, only the phase angle requires iterative solutions when employing the maximum likelihood method for iterative identification. The identification results for the initial amplitude and damping coefficient within a specific range do not compromise the accuracy of frequency parameter identification.

3.3. Comparative Analysis of Recognition Results of Different Methods

The frequency parameters of freely attenuated vibration signals at varying SNRs were identified using the periodogram, dichotomy interpolation, and maximum likelihood methods. One hundred groups of random simulation identifications were conducted for signals with the same SNR. The frequency identification error is depicted in Figure 8. The maximum likelihood method significantly enhanced the identification accuracy. This improvement was particularly noticeable when the SNR is low, leading to a notable increase in the frequency recognition correctness rate.
The variances in recognition results for the three algorithms were calculated based on the recognition results for each group of signals with identical SNRs. The CRLB bounds for these results were also computed under the respective conditions. Table 7 lists the findings, whereas Figure 9 displays the relationship curves between frequency identification variance and signal-to-noise ratio for the various identification methods.
From the error analysis, several conclusions can be drawn:
  • The recognition accuracy of the dichotomy interpolation method marginally surpasses that of the periodogram method.
  • The maximum likelihood method’s recognition accuracy is significantly superior to the periodogram and the dichotomy interpolation methods across all SNR conditions. Approximately an order of magnitude enhances its accuracy.
  • As the SNR decreases, the variance in the maximum likelihood recognition approaches the lower bound of the CRLB variance, indicating that the maximum likelihood recognition algorithm is highly noise-resistant.
  • Enhancing the frequency resolution and the number of data calculation bits of the maximum likelihood method can further boost recognition accuracy.
  • The maximum likelihood method requires an iterative solution, making its computational demand considerably greater than that of the periodogram method.

4. Conclusions

The study focuses on the free-damping signal mixed with white noise to enhance the frequency parameter identification accuracy of the structural response signal under ambient excitation. The primary conclusions are:
1.
The free-damping vibration signal possesses limited energy. As t approaches infinity, a ~ approaches zero. With an extended sampling time, the SNR decreases; when the sample time T approaches infinity, the SNR ς ~ approaches zero.
2.
Based on the CRBL for the free damping signal and the simulation results, the following observations were made:
(a)
An increase in the sampling frequency reduces the lower bound of frequency identification variance. Specifically, when the sampling frequency is low, raising it improves frequency identification accuracy. However, increasing the sampling frequency does not always enhance recognition accuracy. After the sampling frequency reaches 40 times the recognition frequency, further increases in the sampling frequency minimally affect the recognition accuracy.
(b)
A higher signal damping coefficient results in greater frequency identification error. For systems with large damping coefficients, the effect of the damping coefficient should be factored in during system frequency identification.
(c)
The CRB for frequency parameter identification is proportional to a 2 σ 2 . Increased noise leads to decreased frequency identification accuracy.
(d)
The influence of the phase angle is less significant than the noise, sampling frequency, and damping coefficient. However, phase angles can vary due to the unpredictability of excitation in actual sampled signals. Consequently, the superposition of free-damping signals with identical frequencies but differing phase angles can split spectral lines. This variability can reduce the accuracy of actual frequency identification.
3.
Compared to the periodogram method and the dichotomy interpolation method, the recognition accuracy of the maximum likelihood method is substantially superior under various SNR conditions, improving accuracy by roughly an order of magnitude.

Author Contributions

Conceptualization, C.D., X.C. and L.T.; methodology, C.D., X.C. and J.W.; formal analysis, C.D. and J.W.; investigation, J.W. and F.P.; writing—original draft preparation, C.D. and J.W.; writing—review and editing, L.T., X.C. and F.P. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge the financial support of the National Key Research and Development Program of China (SN: 2021YFC3090103); the Water Conservancy Technology Demonstration Project (SN: SF-202309); and the Special Fund Project of Basic Scientific Research Business Fee for Central Public Welfare Scientific Research Institutes (SN: Y422009, Y422010).

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the author C.D. ([email protected]) on request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Masi, A.; Santarsiero, G.; Chiauzzi, L. Development of a seismic risk mitigation methodology for public buildings applied to the hospitals of Basilicata region (Southern Italy). Soil Dyn. Earthq. Eng. 2014, 65, 30–42. [Google Scholar] [CrossRef]
  2. Guo, Y.; Guo, S. Structural Damage Detection And Image Based On The Zero Migration Record Of The Synthetical Pseudo-Spectrum. Eng. Mech. 2014, 31, 29–36. [Google Scholar] [CrossRef]
  3. Gromyko, P.V.; Fedin, K.V.; Seleznev, V.S.; Kolesnikov, Y.I.; Ngomayezwe, L. Monitoring of the hydraulic units operation of the Sayano-Shushenskaya hydroelectric power plant using remote seismic observations. Earthq. Eng. Struct. Dyn. 2023, 52, 335–349. [Google Scholar] [CrossRef]
  4. Bhowmik, B.; Tripura, T.; Hazra, B.; Pakrashi, V. Robust linear and nonlinear structural damage detection using recursive canonical correlation analysis. Mech. Syst. Signal Process. 2019, 136, 106499. [Google Scholar] [CrossRef]
  5. Bhowmik, B.; Tripura, T.; Hazra, B.; Pakrashi, V. Real time structural modal identification using recursive canonical correlation analysis and application towards online structural damage detection. J. Sound Vib. 2019, 468, 115101. [Google Scholar] [CrossRef]
  6. Yuan, P.-P.; Zhang, J.; Feng, J.-Q.; Wang, H.-H.; Ren, W.-X.; Wang, C. An improved time-frequency analysis method for structural instantaneous frequency identification based on generalized S-transform and synchroextracting transform. Eng. Struct. 2022, 252, 113657. [Google Scholar] [CrossRef]
  7. Civera, M.; Calamai, G.; Zanotti Fragonara, L. Experimental modal analysis of structural systems by using the fast relaxed vector fitting method. Struct. Control Health Monit. 2021, 28, e2695. [Google Scholar] [CrossRef]
  8. Xu, G.; Guo, T.; Li, A.; Zhang, R.; Zhu, R.; Liu, F. Influence of ambient temperature on seismic performance of elastomeric isolation structural system. Structures 2023, 55, 1763–1773. [Google Scholar] [CrossRef]
  9. Çetindemir, O.; Zülfikar, A.C.; Apaydın, N.M. Dynamic Modal Identification of a Long-Span Suspension Bridge Using Measured Acceleration Data during a Storm. J. Struct. Eng. 2023, 149, 11982. [Google Scholar] [CrossRef]
  10. MacLeod, E.; Arjomandi, K. Dynamic bridge weigh-in-motion using estimated modal parameters from ambient vibration tests. Eng. Struct. 2023, 289, 116254. [Google Scholar] [CrossRef]
  11. Yi, Q.; Wilcox, P.; Hughes, R. Modelling and evaluation of carbon fibre composite structures using high-frequency eddy current imaging. Compos. Part B Eng. 2023, 248, 110343. [Google Scholar] [CrossRef]
  12. Khademi, P.; Mousavi, M.; Dackermann, U.; Gandomi, A.H. Time–frequency analysis of ultrasonic signals for quality assessment of bonded concrete. Constr. Build. Mater. 2023, 403, 133062. [Google Scholar] [CrossRef]
  13. He, W.; Lai, W.W.-L.; Sui, X.; Giannopoulos, A. Delamination characterization in thin asphalt pavement structure using dispersive GPR data. Constr. Build. Mater. 2023, 402, 132834. [Google Scholar] [CrossRef]
  14. Damadipour, M.; Tarinejad, R. Structural system identification based on combining weighted transmissibility and wavelet transform. Struct. Control Health Monit. 2021, 29, e2868. [Google Scholar] [CrossRef]
  15. Ahi, M.; Ahmadian, H. Nonlinear model updating of frictional structures through frequency-energy analysis. Nonlinear Dy-Namics 2022, 110, 95–116. [Google Scholar] [CrossRef]
  16. Hou, J.; Cao, S.; Hu, H.; Zhou, Z.; Wan, C.; Noori, M.; Li, P.; Luo, Y. Vortex-Induced Vibration Recognition for Long-Span Bridges Based on Transfer Component Analysis. Buildings 2023, 13, 2012. [Google Scholar] [CrossRef]
  17. Ding, Z.; Li, J.; Hao, H. Non-probabilistic method to consider uncertainties in structural damage identification based on Hybrid Jaya and Tree Seeds Algorithm. Eng. Struct. 2020, 220, 110925. [Google Scholar] [CrossRef]
  18. Sadeghi, F.; Yu, Y.; Zhu, X.; Li, J. Damage identification of steel-concrete composite beams based on modal strain energy changes through general regression neural network. Eng. Struct. 2021, 244, 112824. [Google Scholar] [CrossRef]
  19. Zhao, B.; Lei, D.; Fu, J.; Yang, L.; Xu, W. Experimental study on micro-damage identification in reinforced concrete beam with wavelet packet and DIC method. Constr. Build. Mater. 2019, 210, 338–346. [Google Scholar] [CrossRef]
  20. Kato, T. Perturbation Theory for Linear Operators; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013; Volume 132, p. 132. [Google Scholar]
  21. Wang, H.; Barone, G.; Smith, A. A novel multi-level data fusion and anomaly detection approach for infrastructure damage identification and localisation. Eng. Struct. 2023, 292, 116473. [Google Scholar] [CrossRef]
  22. Soleymani, A.; Jahangir, H.; Rashidi, M.; Mojtahedi, F.F.; Bahrami, M.; Javanmardi, A. Damage Identification in Reinforced Concrete Beams Using Wavelet Transform of Modal Excitation Responses. Buildings 2023, 13, 1955. [Google Scholar] [CrossRef]
  23. Eiras, J.N.; Payan, C.; Rakotonarivo, S.; Garnier, V. Experimental modal analysis and finite element model updating for structural health monitoring of reinforced concrete radioactive waste packages. Constr. Build. Mater. 2018, 180, 531–543. [Google Scholar] [CrossRef]
  24. Feng, D.; Feng, M.Q. Computer vision for SHM of civil infrastructure: From dynamic response measurement to damage detection—A review. Eng. Struct. 2018, 156, 105–117. [Google Scholar] [CrossRef]
  25. Yang, C. Study on Frequency Estimators under Noise Environment. Ph.D. Thesis, South China University of Technology, Guangzhou, China, 2010. [Google Scholar]
  26. Cao, Y. Research on Frequency Estimator for a Noisy Real Sinusoid. Ph.D. Thesis, South China University of Technology, Guangzhou, China, 2012. [Google Scholar]
  27. Xin, Y.; Hao, H.; Li, J.; Wang, Z.-C.; Wan, H.-P.; Ren, W.-X. Bayesian based nonlinear model updating using instantaneous characteristics of structural dynamic responses. Eng. Struct. 2019, 183, 459–474. [Google Scholar] [CrossRef]
  28. Aboutanios, E.; Mulgrew, B. Iterative frequency estimation by interpolation on Fourier coefficients. IEEE Trans. Signal Process. 2005, 53, 1237–1242. [Google Scholar] [CrossRef]
  29. Umesh, S.; Tufts, D. Estimation of parameters of exponentially damped sinusoids using fast maximum likelihood estimation with application to NMR spectroscopy data. IEEE Trans. Signal Process. 1996, 44, 2245–2259. [Google Scholar] [CrossRef]
  30. Zhang, J.; Swain, A.K.; Nguang, S.K. Parameter estimation of exponentially damped sinusoids based on state-space model. J. Jilin Univ. 2019, 49, 2083–2088. [Google Scholar]
  31. Chen, G.-W.; Omenzetter, P.; Beskhyroun, S. Operational modal analysis of an eleven-span concrete bridge subjected to weak ambient excitations. Eng. Struct. 2017, 151, 839–860. [Google Scholar] [CrossRef]
  32. Gentile, C.; Saisi, A. Operational modal testing of historic structures at different levels of excitation. Constr. Build. Mater. 2013, 48, 1273–1285. [Google Scholar] [CrossRef]
  33. Giordano, E.; Mendes, N.; Masciotta, M.G.; Clementi, F.; Sadeghi, N.H.; Silva, R.A.; Oliveira, D.V. Expeditious damage index for arched structures based on dynamic identification testing. Constr. Build. Mater. 2020, 265, 120236. [Google Scholar] [CrossRef]
  34. Kay, S.M. Fundamentals of Statistical Signal Processing: Estimation Theory; Prentice Hall Ptr: Upper Saddle River, NJ, USA, 1993. [Google Scholar]
Figure 1. The random excitation response signal of single−degree−of−freedom system: (a) time domain signal; (b) power spectral density(psd) and the frequency.
Figure 1. The random excitation response signal of single−degree−of−freedom system: (a) time domain signal; (b) power spectral density(psd) and the frequency.
Buildings 14 00198 g001
Figure 2. The relation between the CRLB and sampling frequency.
Figure 2. The relation between the CRLB and sampling frequency.
Buildings 14 00198 g002
Figure 3. The relation between CRLB and damping coefficient.
Figure 3. The relation between CRLB and damping coefficient.
Buildings 14 00198 g003
Figure 4. The relation between CRLB and SNR.
Figure 4. The relation between CRLB and SNR.
Buildings 14 00198 g004
Figure 5. The relation between CRLB and phase.
Figure 5. The relation between CRLB and phase.
Buildings 14 00198 g005
Figure 6. Flowchart of maximum likelihood algorithm.
Figure 6. Flowchart of maximum likelihood algorithm.
Buildings 14 00198 g006
Figure 7. Time domain curves of different SNRs.
Figure 7. Time domain curves of different SNRs.
Buildings 14 00198 g007
Figure 8. The box plot of the frequency identification error at different SNRs: (a) SNR = 30 dB; (b) SNR = 5 dB; (c) SNR = −10 dB; (d) SNR = −30 dB.The red dots represent discrete data points outside the interval.
Figure 8. The box plot of the frequency identification error at different SNRs: (a) SNR = 30 dB; (b) SNR = 5 dB; (c) SNR = −10 dB; (d) SNR = −30 dB.The red dots represent discrete data points outside the interval.
Buildings 14 00198 g008
Figure 9. Relation of frequency identification variance and SNRs of different identification methods.
Figure 9. Relation of frequency identification variance and SNRs of different identification methods.
Buildings 14 00198 g009
Table 1. CRLB of frequency identification corresponding to different sampling frequencies.
Table 1. CRLB of frequency identification corresponding to different sampling frequencies.
f s (Hz)3050100150200250300400500
CRLB0.01740.01020.00510.00340.00250.00200.00170.00130.0010
Table 2. Different damping coefficients correspond to the CRLB.
Table 2. Different damping coefficients correspond to the CRLB.
ξ 0.010.020.030.040.050.060.080.10
CRLB8.0465 × 10−66.4468 × 10−62.1812 × 10−45.1882 × 10−40.00100.00180.00420.0084
Table 3. The CRLB of frequency identification at different SNRs.
Table 3. The CRLB of frequency identification at different SNRs.
σ 0.01000.01640.02700.03470.04460.05730.07360.12130.2000
σ 2 / a 2 0.00010.00030.00070.00120.00200.00330.00540.01470.0400
SNRs (dB)30201050−5−10−20−30
CRLB2.545 × 10−66.844 × 10−61.855 × 10−53.064 × 10−55.062 × 10−58.355 × 10−51.378 × 10−43.744 × 10−41.018 × 10−3
Table 4. The CRLB of frequency identification at different phases.
Table 4. The CRLB of frequency identification at different phases.
Phase (°) 0365490126144180198360
CRLB0.001020.000960.000950.000990.001050.001050.001020.000990.00102
Table 5. Identification results of the four iterations for initial estimation of different frequencies (Hz).
Table 5. Identification results of the four iterations for initial estimation of different frequencies (Hz).
Initial Estimation8.08.49.09.610.010.411.011.412.0
Results of one iteration9.7029.7289.7909.8949.99010.08010.18010.22210.264
Second iteration result9.8489.8629.8929.9449.99210.04010.08010.10410.124
Triple iteration result9.9209.9269.9409.9649.98810.01010.03010.04410.054
Results of four iterations9.9749.9769.9789.9869.9909.99410.00010.00210.006
Table 6. Frequency identification results of initial estimation of different damping coefficients.
Table 6. Frequency identification results of initial estimation of different damping coefficients.
Initial Estimation0.0250.0350.0400.0450.0500.0550.0600.0700.075
Frequency estimation (Hz)10.00710.0089.9999.9979.9939.9899.9879.9779.973
Table 7. Frequency parameter identification variance of different identification methods.
Table 7. Frequency parameter identification variance of different identification methods.
σ 0.01000.01640.02700.03470.04460.05730.07360.12130.2000
SNR dB30201050−5−10−20−30
perodogram0.00230.00330.00530.00600.01150.01260.01450.03010.1032
About0.00140.00240.00450.00540.01010.01210.01390.02970.1013
ML1.771 × 10−42.548 × 10−43.638 × 10−44.528 × 10−47.058 × 10−49.818 × 10−41.100 × 10−32.200 × 10−34.600 × 10−3
CRLB2.545 × 10−66.844 × 10−61.855 × 10−53.064 × 10−55.062 × 10−58.355 × 10−51.378 × 10−43.744 × 10−41.018 × 10−3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Deng, C.; Wen, J.; Tang, L.; Cai, X.; Peng, F. The Accuracy of Frequency Estimation of Structure Vibration under Ambient Excitation: Problems, Causes, and Methods. Buildings 2024, 14, 198. https://doi.org/10.3390/buildings14010198

AMA Style

Deng C, Wen J, Tang L, Cai X, Peng F. The Accuracy of Frequency Estimation of Structure Vibration under Ambient Excitation: Problems, Causes, and Methods. Buildings. 2024; 14(1):198. https://doi.org/10.3390/buildings14010198

Chicago/Turabian Style

Deng, Chang, Jiaqi Wen, Lei Tang, Xin Cai, and Feng Peng. 2024. "The Accuracy of Frequency Estimation of Structure Vibration under Ambient Excitation: Problems, Causes, and Methods" Buildings 14, no. 1: 198. https://doi.org/10.3390/buildings14010198

APA Style

Deng, C., Wen, J., Tang, L., Cai, X., & Peng, F. (2024). The Accuracy of Frequency Estimation of Structure Vibration under Ambient Excitation: Problems, Causes, and Methods. Buildings, 14(1), 198. https://doi.org/10.3390/buildings14010198

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop