Research Article  Open Access
Damped CAPES 2D Spectral Estimation for RealValued Vibration Signals
Abstract
We propose a 2D representation in the frequencydecay factor plane of an arbitrary realworld vibration signal. The signal is expressed as the sum of a decayedattenuation sine term modulated by an amplitude function and a noise residue. We extend the combined approach of Capon estimation and amplitude and phase estimation (CAPES) to damped real vibration signals (DRCAPES). In the proposed DRCAPES method, the highresolution amplitude and phase are estimated simultaneously for both angular frequency and decay factor grids. The performance of the proposed approach is tested numerically with noisy vibration data. Results show that the DRCAPES method has an excellent frequency resolution, which helps to overcome difficulties in spectrum estimation when vibration modes are very close, and a small bias, which makes it suitable for obtaining accurate amplitude spectrums. The results also indicate that the proposed method can accurately estimate the amplitude spectrum with the use of averaging and denoising processes.
1. Introduction
The normal functioning of engineering structures in dynamic environments is often affected by large and frequent vibrations. These further induce negative effects in terms of comfort in use, fatigue durability, dynamic instability, and safety. Therefore, the structural vibration levels must be assessed under working load conditions and environmental excitation. Typically, these assessments should be conducted in the frequency domain, and thus, structural dynamics experiments have focused on spectral analysis and identification of frequencydomain dynamic characteristics, specifically the structural vibration signal. The vibration level of the structural response is usually described using the Fourier power spectrum [1, 2]; however, the signal quality could be poor because of noise pollution, short signal length, or inappropriate sampling frequency. Because the Fourier power spectrum frequency often has low resolution and the spectral value estimates are inaccurate, a more suitable and accurate spectral estimation algorithm is needed. The power density spectrum cannot be directly used to assess the vibration amplitude of a certain vibration response. In some cases, the amplitude spectrum is even more important than the power spectrum [3]. Vibrations in the physical world are damped signals and the level of vibration is closely related to the damping. Damping describes the dynamic characteristics and structural performance of the structure; therefore, it must be precisely measured. At present, highquality damping estimation is very difficult when using the real response signal. In this paper, we propose a nonparametric highprecision spectral estimation algorithm that determines the vibration signal characteristics of engineering structures in order to analyze data from structural dynamics experiments.
One classic category of spectral estimation algorithms is the set of nonparametric methods based on the fast Fourier transform (FFT), in which the frequency resolution is low. However, the FFT has wide applicability, because it has no prespectral shape [3–5]. There is also a category of highresolution parametric methods that require the observed signal to be composed of limited number of superposed sinusoidal signals contaminated by white noise. This is difficult to satisfy in structural vibration signal spectral analyses; therefore, these methods often perform poorly [6]. A third category is the set of nonparametric adaptive filterbank methods with high resolution capability, including the Capon algorithm [7] and the recently proposed amplitude and phase estimation spectral (APES) algorithm [8, 9], as well as variants and hybrids of these techniques. These methods have been used successfully for almost twenty years. They make no assumptions about the spectral structure and do not depend on the Fourier transform algorithm; as a result, they can be used for highresolution spectrum estimation for any stationary signals.
The Capon algorithm performs frequency resolution with greater efficiency than other algorithms, and the APES algorithm estimates spectral values [10, 11] with greater accuracy than most of the other algorithms. In 2000, Jakobsson and Stoica combined the two to give the CAPES algorithm, which improves the resolution and estimation accuracy of spectral values [12]. The CAPES algorithm has been used for one and twodimensional complexsignal spectral estimation in medical imaging and remote sensing [13–15], and Jakobsson et al. improved the algorithm to handle real onedimensional vibration signals without damping [16]. In 2001, Stoica and Sundin proposed new versions of the Capon and APES algorithms that can deal with complex onedimensional damping modulation signals, whose amplitude decays over time [14]. Recently, the researchers focused on the applications and fast computation algorithms of APES and CAPON method [17, 18]. However, when their algorithms are used to estimate the spectrum of the structural vibration directly the accuracy will be worse due to the imaginary part of the signals. Actually, the structural vibration signal is a real onedimensional damping modulation signal.
Based on existing nonparametric estimation algorithms, we propose a twodimensional spectral estimation algorithm to transform arbitrary real onedimensional damping vibration signals in the timedomain into a twodimensional signal in a joint domain of circular frequency and damping factor. The proposed method is expected to provide a highresolution frequency, highprecision amplitude spectrum, and highprecision damping factor estimations of the vibration signal. Through numerical examples, we verify the effectiveness of the proposed algorithm.
2. Theory
2.1. Representation of a Vibration Signal in a 2D Plane
In an degreesoffreedom (DOF) vibration system, the response to a steady forced vibration can be written as follows: where is the contribution of th mode of the vibration system, and are the th mode frequency and damping factor, and and are the amplitude and initial phase of th mode. This equation means that the response to any DOF can be represented by the decayedattenuation sine function . Based on this, we set as the discrete time series of an arbitrary DOF response in the vibrating system, , where is the length of time series. With interest frequency and damping factor , can be written as the sum of a decayedattenuation sine term modulated by the amplitude function and a residual term , as follows:
In (2), and are the amplitude function and the phase function, respectively, of the point , which is defined in a 2D plane; includes the decayed attenuation sine terms of all points except for and the random noise .
For the single DOF free response, where and are the initial displacement and the initial velocity, respectively, when (or ).
We define the relationship between the damping factor and the damping ratio as follows:
The geometrical meaning can be explained as the sine vector between the point and the original point in a 2D plane, as shown in Figure 1.
2.2. Damped CAPES Spectral Estimation for RealValued Vibration Response Data
In (2), the time series can be written as
If we let then
(i) Hankel Matrix Method. We can obtain practical samples of the time series , such as , . According to the nonparametric adaptive filterbank spectral estimation method [1], we can construct the order Hankel matrix of observation data ,
In (8), is the user parameter, also known as the finite impulse response (FIR) filter length [9], , and represents the total number of snapshots in the time series.
If we set Here, in (11) and (12), .
Then and . Equation (8) becomes where noise matrix is , We now construct the weighted secondary residuals term as the least squares estimator, as follows: where is the weight matrix and denotes the square root of Hermitian matrix . is defined as
Then, according to the nonparametric adaptive filterbank theory on spectral estimation, the weighted least squares estimator of parameter , which includes amplitude information, is given by
According to Larsson and Li’s contribution [10], the final estimation of parameter can be expressed as follows: Here, differently with the Larsson and Li’s expressions, is defined as (9) and , should be defined by
The weighting function values are as follows:
(ii) Vectorization of Hankel Matrix. According to (8), we define
The length of is .
Let
Then, according to (7),
In (27), where is the observation vector, is the observation matrix, , and , .
Spectral estimation for a signal at point has become the problem of parameter estimation of , based on , while minimizing the estimated noise. This is a typical linear estimation problem that can be solved by the following weighted least squares estimator as seen in [5]: where is the weighted values matrix. When is nonsingular, the former equation is solved by and then (21) reaches a minimum.
When the weighted matrix is set to unit matrix, we can get a twodimensional spectral estimation of the real signal in the least squares sense as follows:
According to the basic idea of Capon spectral estimation, is set to the covariance matrix of data as follows:
However, according to the idea of APES spectral estimation, is the covariance estimation of the residual , given by
Thus, we have five 2D damping spectral estimations of a real vibration signal, see in (23), (24), (31), (32), and (33).
3. Algorithm Implementation
Previous methods of estimation required the inverse of large weightedmatrices to be calculated at each 2D point . To reduce the computation time, we can utilize a matrix inversion formula to calculate the inverses in advance, leaving only a matrix multiplication in every iteration step.
Let , , , and be arbitrary matrices; according to the matrix inversion formula, we have
Corresponding to (24), let be , let be , let be , and let be ; we have and corresponding to (33), let
This gives
The large matrix inversions in (30) and (32) do not depend on the point ; therefore, these can be calculated before entering the iteration cycle.
Hence, we design the proposed real damped version of CAPES (DRCAPES) as follows:step 1:set the data as , and the user data length as ;step 2:calculate , form the order Hankel matrix and calculate ;step 3:form the 2D grid , where and are the divisors of the frequency and damping factor, respectively; , ;step 4:incrementing in a loop, calculate and according to (21) and (22); calculate according to (35); calculate according to (20); calculate the amplitude and phase of the real damping vibration signal according to (7).
4. Numerical Case Studies
Consider a 5DOF linear oscillation system with the following structural matrix:
It should be noted that, in order to simplify the problem and focus on the process of vibration signal itself, this numerical case of 5DOF oscillation system is designated as full decoupling system, with all of its physical matrices , , and being diagonal matrix. The modal parameters of the complex modal matrix, angular frequency, and damping factor can be calculated analytically as follows:
The free vibration responses are generated from the initial status vector . These responses are observed fully under an initial signaltonoise ratio (SNR) of and a sampling frequency of 1.0 Hz, as shown in Figure 2(a). The total length of the data is 800 samples and there are 40, 48, 56, 64, and 72 characteristic periods for each mode. If we define the effective length of the data to be the segment with a nonnegative SNR, these numbers are reduced to 24, 17, 21, 25, and 38, as illustrated in Figure 2(b). These characteristic periods are shorter than any other spectral estimation approaches, thereby illustrating the effectiveness of feature extraction in the proposed approach.
(a)
(b)
The response of the first DOF has the biggest initial SNR of all the responses. Therefore, was chosen to feed the proposed approach. The 3D amplitude spectrum and its corresponding contour in the ()plane are shown in Figure 3.
(a)
(b)
In Figure 3(a), the five distinct peaks lie in the plane of the natural frequency and damping factor, that is, , respectively. The bases of the peaks are well separated from each other and the widths in the frequency dimension are so narrow that each peak looks like a thin blade. In the damping factor dimension, the slopes of the ridges are gentle and the bases are wide. Each apex is located exactly at its corresponding designated value, as shown in Figure 3(b). The identified spectrum data are listed in Table 1.

The proposed approach performed well in identifying both the frequency and damping factor. Under the noise levels given above, the frequency errors were almost zero. The damping factor errors were also smaller with the proposed method than with any other approach.
4.1. Frequency Resolution of 3D Amplitude Spectrum
The traditional CAPES method is known for its high frequency resolution and accurate spectral value estimation. To investigate the frequency resolution of the proposed DRCAPES approach, we conducted a case study based on the oscillation system described above.
The effect of closed peaks was simulated by changing the physical parameters , , or . In this case, we simply fixed the second mode to 0.06 Hz and gradually adjusted the first mode frequency from 0.05 Hz to 0.06 Hz while maintaining and . A new 5DOF oscillation system with a different was obtained for each iteration and the free vibration responses of the system were calculated under the initial SNRs mentioned above.
Figure 4 shows the amplitude spectra of four typical steps when the first mode approaches the second. The other modes were not identified and are not illustrated in the figures. As a comparison, the global amplitude spectrum identified by an FFT method was also found and the nine spectral lines corresponding to the same frequency window are plotted in the figures. It can be seen that the frequency resolution of DRCAPES is much better than that of the FFT method; from Figures 4(a)–4(c), the two peaks of the FFT spectrum approach one another quickly and in Figure 4(d) the peaks have overlapped completely, whereas the two 3D peaks given by the DRCAPES method remain separate and distinct, with the shape and height of the peaks basically invariant.
(a)
(b)
(c)
(d)
To further illustrate the high resolution of the proposed approach, we define a mode closeness measure and frequency resolution in a physical sense as follows.
Let be a measure of the closeness of two adjacent modes (or adjacent peaks) in a spectrum, in which and is the th mode or peak in spectrum. Larger values imply closer distances between two spectrum lines, modes, or peaks.
When adjacent peaks are close to each other, the physicalsense frequency resolution can be defined as the logarithm of a ratio of the minimum height of two adjacent peaks and the saddle value between them, as follows:
Here, is the saddle value of two adjacent hills in the spectrum and , are the values of the two adjacent peaks.
Figure 5 illustrates how the results of the two processes converge in terms of the mode closeness and frequency resolution. As the mode closeness increases, the frequency resolutions of both the DRCAPES method and the FFT method decrease. The resolution of the FFT method dropped close to zero at a relatively small value of mode closeness. Compared to FFT, the proposed DRCAPES has a higher frequency resolution, even when the two peaks are extremely close to each other.
4.2. Bias and Accuracy of the Amplitude Spectrum Given by the DRCAPES Method
The bias and accuracy of the amplitude spectrum identified by the DRCAPES method can be revealed by numerical simulations. Four simulations were performed under a fixed observation noise level; that is, the initial SNR was 16.6 dB for the first mode component of the response corresponding to the 1st DOF. Each simulation was run with the same parameters, except for the damping factor, which was set as , 0.0015, 0.0025, and 0.003. In each simulation, the following processes were repeated 1000 times: generation of random response, simultaneous spectrum estimation with the proposed method, and the FFT method.
Figure 6 shows the scatterplots of the simulation results. Each subplot corresponds to a designated damping factor, that is, , 0.0015, 0.0025, and 0.003. The mean and bias of these estimates by these two methods, indicated as , , , and , and the true values of the amplitude (designated value) are also plotted.
(a)
(b)
(c)
(d)
The results show that the spectrum bias of the DRCAPES method is smaller than that of the FFT method. Regardless of the designated damping factor, the bias of the DRCAPES method stabilized to a certain level below the true amplitude, whereas the FFT estimation was biased on either side depending on the designated damping factor. Additionally, as shown in Figure 6, the extent of scatter of the FFT result is smaller than that of the proposed method; this may mean that the accuracy of the proposed DRCAPES method can be further improved.
4.3. Influence of Noise on the Amplitude Spectrum
In the above discussion, the noise is fixed at a very low level. Further examination investigated the influence of noise level on the quality of the estimated amplitude spectrum. Similar to the numerical simulation case studies above, we simulated four different damping factors. In each simulation, the noise level was divided into 100 steps with an initial SNR interval of dB. At each noise level, the simulation process included the generation of a random response and estimation of the amplitude spectrum using the proposed method and the FFT method. The mean and standard deviation of the results given by both methods at different SNRs are illustrated in Figure 7. For brevity, we only present three of the typical damping factors in this paper; that is, 0.003, 0.002, and 0.001.
(a) HZ,
(b) HZ,
(c) HZ,
(d) HZ,
(e) HZ,
(f) HZ,
Figures 7(a), 7(c), and 7(e) show that the proposed DRCAPES method produces less bias than the FFT method and that damping factor effects were separated from the amplitude effects in the amplitude spectrum estimation. The proposed method is not affected by the damping factor, whereas the FFT method performed poorly with a swingintolerant bias when the damping factor was changed.
The curves for both methods gently plateau when the initial SNR is set to a specific value, that is, −6.0 dB for Figure 7(a) and −10 dB for both Figures 7(c) and 7(e). Before this plateau, the mean of the estimated amplitude dropped sharply from a value that is far higher than the true amplitude or designated value. The plateau in the curves proves that an optimum workable SNR zone exists for each individual freevibration response, which implies that a denoising process is needed if the initial SNR of the raw response data is below this workable zone.
Figures 7(b), 7(d), and 7(f) show that the observation noise level has a strong influence on the statistical qualities of the amplitude spectrum estimation, that is, the standard deviation σ. As the initial SNR was increased, the standard deviation of both methods decreased. Under most noise levels, the standard deviation of the proposed method was significantly greater than that of the FFT method. The reason for this may be that the spectrum estimation of the proposed method is more complicated than that of the FFT method, where variance is propagated. A higher SNR means a smaller standard deviation and statistically better spectrum estimation; this conclusion indicates that we should focus on improving the quality of the amplitude spectrum.
5. Conclusion
The vibration response of structures in the real physical world is evaluated and modulated by amplitude damping. By representing the free vibration response in two dimensions (frequency and damping factor), we extended the hybrid CAPES approach, which is a highresolution amplitude spectrum estimation method, to estimate the amplitude spectrum of the structural vibration response. Compared to the FFT method, the proposed DRCAPES has a much better resolution in the frequency domain and is, therefore, more suitable for overcoming the difficulties in spectrum estimation when vibration modes are close together. The downward bias of the amplitude spectrum estimated by DRCAPES is smaller than that given by the FFT method and is not influenced by the damping factor; therefore, the proposed method is more accurate than FFTlike approaches when determining amplitude. However, with higher SNRs, secondary operations, such as averaging and denoising, are needed to improve the amplitude spectrum estimate.
Notations
Damping factor, decay  
:  Damping ratio 
:  Logarithmic decrement 
:  Damping coefficient 
:  Damped natural angular frequency 
:  Natural angular frequency 
:  Lumped mass at the th degree of freedom (modal) 
:  Modal stiffness at the th degree of freedom 
:  Initial displacement 
:  Total number of snapshots in the time series 
:  Weight matrix for Capon spectral estimation method 
:  Vectorization of Hankel matrix 
:  Angular frequency vector 
:  Modal phase matrix 
:  Frequency resolution 
:  Defined as (9) 
:  Logarithm of frequency resolution 
:  Sinusoidal component of amplitude 
:  Modal stiffness at the th degree of freedom 
:  Diagnosis mass matrix (lumped mass matrix) 
:  Modal stiffness matrix (diagnosis) 
:  Modal damping matrix (diagnosis) 
:  Mass matrix 
:  Damping matrix 
:  Stiffness matrix 
:  Angular frequency (measured) 
:  Initial velocity 
:  Hankel matrix of observation data 
:  Least squares estimator matrix of the weighted secondary residuals term 
:  Weight matrix for APES spectral estimation method 
:  Observation matrix 
:  Damping factor vector 
:  Initial signaltonoise ratio of each contribution in each response 
:  Defined by (21) 
:  Defined by (12) 
:  Cosine component of amplitude 
:  Exponential function 
:  Number of degrees of freedom 
:  Response of steady forced vibration 
:  Amplitude function 
:  Phase function 
:  Time series 
:  Natural frequency 
:  Residual term 
:  Length of the time series 
:  Finite impulse response (FIR) filter length 
:  Solution to the weighted least squares parameter and its estimator 
:  Weighted least squares parameter and its estimator 
:  2D spectral estimation of real signal in the sense of least squares 
:  Modal amplitude matrix 
:  Closeness of two adjacent modes (or adjacent peaks) in spectrum 
:  Defined by (22) 
:  Defined by (11) 
:  Noise matrix. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research is supported by the State HighTech Research and Development Plans (863), Grant no. 2014AA110402; State Meteorological Administration Special Funds of Meteorological Industry Research, Grant no. 201306102; the Project of National Key Technology R&D Program in the 12th Five Year Plan of China (Grant no. 2012BAJ11B01); the Fundamental Research Funds for the Central Universities.
References
 J. Penzien, “Dynamics of structures,” Computers and Structures, Inc., 2003. View at: Google Scholar
 D. J. Ewins, “Modal Testing, Theory, Practice, and Application,” 2000. View at: Google Scholar
 B. P. Lathi, Linear Systems and Signals, Oxford University Press, New York, NY, USA, 2nd edition, 2004.
 P. Stoica and R. Moses, Introduction to Spectral Analysis, Prentice Hall, Upper Saddle River, NJ, USA, 1997.
 D. P. Kalogeras, An Overview of the Spectral Estimation Methods, vol. 12 of Mathematics and Computers in Science and Engineering Series, 2013.
 D. G. Manolakis, V. K. Ingle, and S. M. Kogon, Statistical and Adaptive Signal Processing, McGrawHill, 2000.
 J. Capon and S. Haykin, “Maximumlikelihood spectral estimation,” in Nonlinear Methods of Spectral Analysis, Springer, New York, NY, USA, 1983. View at: Google Scholar
 J. Li and P. Stoica, “An adaptive filtering approach to spectral estimation and SAR imaging,” IEEE Transactions on Signal Processing, vol. 44, no. 6, pp. 1469–1484, 1996. View at: Publisher Site  Google Scholar
 P. Stoica, H. Li, and J. Li, “New derivation of the APES filter,” IEEE Signal Processing Letters, vol. 6, no. 8, pp. 205–206, 1999. View at: Publisher Site  Google Scholar
 E. G. Larsson and J. Li, HighResolution Nonparametric Spectral Analysis: Theory and Applications, 2003.
 H. Li, J. Li, and P. Stoica, “Performance analysis of forwardbackward matchedfilterbank spectral estimators,” IEEE Transactions on Signal Processing, vol. 46, no. 7, pp. 1954–1966, 1998. View at: Publisher Site  Google Scholar
 A. Jakobsson and P. Stoica, “Combining Capon and APES for estimation of spectral lines,” Circuits, Systems, and Signal Processing, vol. 19, no. 2, pp. 159–169, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. G. Larsson, G. Liu, P. Stoica, and J. Li, “Highresolution SAR imaging with angular diversity,” IEEE Transactions on Aerospace and Electronic Systems, vol. 37, no. 4, pp. 1359–1372, 2001. View at: Publisher Site  Google Scholar
 P. Stoica and T. Sundin, “Nonparametric NMR spectroscopy,” Journal of Magnetic Resonance, vol. 152, no. 1, pp. 57–69, 2001. View at: Publisher Site  Google Scholar
 Z.S. Liu, R. Wu, and J. Li, “Complex ISAR imaging of maneuvering targets via the Capon estimator,” IEEE Transactions on Signal Processing, vol. 47, no. 5, pp. 1262–1271, 1999. View at: Publisher Site  Google Scholar
 A. Jakobsson, T. Ekman, and P. Stoica, “CAPON and APES spectrum estimation for realvalued signals,” in Proceedings of the IEEE Digital Signal Processing Workshop, August 1998. View at: Google Scholar
 K. S. Ryu, M. Sepehrifar, and A. A. Kishk, “High accuracy peak location and amplitude spectral estimation via tuning APES method,” Digital Signal Processing, vol. 20, no. 2, pp. 552–560, 2010. View at: Publisher Site  Google Scholar
 G. O. Glentis and A. Jakobsson, “Computationally efficient damped Capon and APES spectral estimation,” in Proceedings of the 21st European Signal Processing Conference (EUSIPCO '13), 2013. View at: Google Scholar
Copyright
Copyright © 2014 Danhui Dan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.