Research Article  Open Access
Operational Modal Identification of TimeVarying Structures via a Vector Multistage Recursive Approach in Hybrid Time and Frequency Domain
Abstract
Realtime estimation of modal parameters of timevarying structures can conduct an obvious contribution to some specific applications in structural dynamic area, such as health monitoring, damage detection, and vibration control; the recursive algorithm of modal parameter estimation supplies one of fundamentals for acquiring modal parameters in realtime. This paper presents a vector multistage recursive method of modal parameter estimation for timevarying structures in hybrid time and frequency domain, including stages of recursive estimation of timedependent power spectra, frozentime modal parameter estimation, recursive modal validation, and continuoustime estimation of modal parameters. An experimental example validates the proposed method finally.
1. Introduction
In the real world, many engineering structures, such as trafficexcited bridges, launch vehicles with varying fuel mass, airplanes in flight with varying additional aerodynamic effects, deployable and flexible geometryvariable aerospace structures, and rotating machinery, show properties changing with time. In many reallife applications, the excitation on the timevarying structures is unknown and random so that operational or outputonly methods are appropriate.
In the past decades, many timedomain parametric approaches of dynamic identification for timevarying structures were presented. Petsounis and Fassois [1, 2] presented the timedependent autoregressive moving average (TARMA) representation for the modeling of nonstationary stochastic vibration. Poulimenos and Fassois [3] surveyed and compared several approaches of TARMAbased nonstationary random vibration modeling including unstructured parameter evolution, stochastic parameter evolution, and deterministic parameter evolution. Poulimenos et al. [4, 5] estimated the modal parameters via the TARMAbased approaches and validated them with a laboratory experiment. Liu [6, 7] proposed a state spacebased approach for linear timevarying systems via decomposing a series of Hankel matrices that are assembled by output response data or additional input data with singular value decomposition (SVD). Liu and Deng [8, 9] improved the state spacebased approach for linear timevarying system through making it less sensitive to noise and validated the identification algorithm with an experiment of a moving cantilever beam.
Meanwhile, several timefrequency analysisbased nonparametric identification approaches were developed in the past decade. Ghanem and Romeo [10] presented a waveletbased identification approach, which transforms the classic governing equation of motion into a wavelet expanded form by projecting the physical responses to a series of wavelet coefficients and identified the modal parameters by solving the expandedform equation. RoshanGhias et al. [11] estimated modal parameters using smoothed pseudo WignerVille distribution (SPWVD), which represents the analytical explicit responses onto WVD plots, and estimated the natural frequency and damping ratio of a SDOF system with tracking the ridge of these plots. Meanwhile, some approaches using Hilbert transform (HT) or HilbertHuang transform (HHT) were proposed. Xu et al. [12] decomposed the responses into a series of single components with Gabor expansion and identified the modal parameters of these singlecomponent signals with HT.
The batch methods of the model parameter estimation for timevarying structures are overviewed above. Realtime estimation of modal parameters of timevarying can contribute to some specific applications, such as the health monitoring, the damage detection, and the vibration control. Furthermore, the recursive algorithms of modal parameter estimation supply necessary fundamentals for acquiring modal parameters in realtime. Some recursive methods for timevarying systems are also presented in past years. For instance, Lourens et al. [13] developed an augmented Kalman filter based method for force identification in structural dynamics. Song and Pei [14] presented a recursive approach of flutter analysis based on the timefrequency analysis. Some classic recursive algorithms are also introduced in Ljung’s book [15].
Moreover, recently the methods of mixing the time domain and frequency domain obtained good results. For instance, Calinoiu et al. [16] proposed a new mixed time and frequencydomain method for estimating the frequency. In author’s previous work, the timefrequencydomain twostage least square method of modal parameter estimation for timevarying structures is presented [17], which can estimate the frequency and mode shapes well.
This paper attempts to present a multistage recursive approach for the potential applications of realtime estimation methods of modal parameters in the hybrid time and frequency domain by choosing and reproducing the available timedomain and frequencydomain methods into the four sequential stages. The reminder of this paper is organized as follows. Section 2 proposes the conceptual design of the multistage recursive estimator. Sections 3 to 6 present the four stages of the proposed estimator, respectively. The proposed method is validated by an experiment in Section 7.
2. Design of the Multistage Recursive Estimator
Different from the common estimators of modal parameters of timevarying structures, the four stages of this proposed recursive estimator in the hybrid time and frequency domain are sequential designed. The basic idea and the procedure of the multistage recursive estimator are shown in Figure 1.
As shown in the first box in Figure 1, when the new responses at the time are measured, the response vector at the time is transferred to Stage I estimation. In this stage, the power spectra at the time are estimated. In addition, the frozentime modal parameters will be estimated by the frequencydomain method in the next stage, so the selection of the frequency band is possible, which is an advantage of the frequencydomain methods over the timedomain methods. The capability of the frequency band selection favours to achieve the parsimony of the estimation algorithm and may improve the accuracy and quality of the estimation in a limited frequency band.
In Stage II as shown in the second box in Figure 1, the frozentime modal parameters at time are estimated by the polyreference least square complex frequencydomain method [18] (pLSCF) based on the power spectra at time within the interested frequency band. After the estimation of this stage, the modal parameters at time are obtained, including the physical and mathematical ones, because the order of the parametric model is unknown and a relative high order is given. Therefore, the modal validation is necessary.
For timeinvariant structures, the modal parameters are validated by some specific criteria or the stabilization diagram. Some automatic stabilization diagram approaches have been presented [19]. The latter is commonly used in most cases of experimental modal analysis, which is achieved by varying the order of the parametric model and tracking the invariance of the physical modes with the model order. However, these criteria and diagrams are not directly suitable for timevarying structures due to the time variation. The time variation results in difficulty of comparing across multiple time instants and a massive data set due to the varying time instants as well as the varying model orders in terms of using stabilization diagram. Therefore, in Stage III as shown in the third box in Figure 1, the modal parameters at time estimated in the last stage are validated by the fuzzy clustering in a recursive way. After the validation based on the recursive fuzzy clustering, the modal parameters at time are automatically ordered by the maximal membership and the quality of the corresponding modal parameters is given by their memberships.
In Stage IV as shown in the fourth box in Figure 1, the continuoustimevarying modal parameters are estimated by a recursive least square approach via using the validated frozentime modal parameters from the last stage and the memberships of the corresponding modal parameters are used as the weights in the least square estimation.
In the proposed multistage estimator, the response data at the current and past time instants are used and once the response data at time are measured, the modal parameters at that time instant are estimated. Therefore, this method is recursive, which has the computational parsimony and is potentially suitable for the realtime acquisition of the modal estimation for timevarying structures.
To clarify the workflow of the proposed method through presenting the relations and dataflow among the different stages, Figure 2 shows the detail flowchart of the estimator.
As shown in Figure 2, there are two phases in the estimator for initialization and recursive estimation of the modal parameter estimation. In both phases, the four stages are included, but the purpose of the two phases is partly different. The phase of initialization is to find the prototypes of each mode and to achieve the convergence of Stage I by using a short onebatch data (before the given time instant ), in which Stage III is a general fuzzy clustering of the onebatch data. The phase of recursive estimation achieves the recursive estimation of modal parameter, in which the four stages are all recursive and only the data at current time instant are added through an iterative way.
3. Stage I: Recursive Pseudolinear Regressive Estimation of Power Spectra
3.1. Parametric Representation of Responses for TV Structural Systems
For multioutput linear TV structural systems, the multiresponses can be expressed by the forward vector timedependent autoregressive moving average (VTARMA) model as where is the vector of responses, and are the th autoregressive (AR) and moving average (MA) coefficient matrix, is the uncorrelated residual vector with zero mean and covariance matrix , and and denote the AR and MA orders. The AR coefficient matrix and the MA coefficient matrix are defined by Based on (1) and (2), the VTARMA model is rewritten in a general form as where is the vector form of the combined definition of AR and MA coefficients (see (2)) defined by and is defined by
3.2. Recursive Pseudolinear Regression
As shown in (3), the VTARMA model has been represented by a linear equation. However, the predict errors are not known, so they are approximated by , which are obtained at time using the estimated parameters at time . The response data matrix does not only depend on the responses but also on the ; (3) is actually a pseudolinear equation.
The pseudolinear regression (PLR) for the parameter estimation of the timeinvariant systems has been proposed [15] as which means the solution of a nonlinear equation. Therefore, the iterative solution is necessary. In the th iteration, a linear equation can be achieved as follows:
Thus linear equation (7) can be solved by
Making one new iteration, at the same time a new measurement is brought in; introducing the forgetting factor and using the matrix inverse lemma in the community of recursive regression, the recursive pseudolinear regression (RPLR) [5, 15] can be conducted as follows.
Parameter estimates update:
Prediction errors:
Adaptive gain:
“Covariance” update:
Different from the original PLR, its recursive version, RPLR, can track the timevariant features of the TV structural systems, which are represented by the timedependent parameter vector .
3.3. Power Spectra Estimation Based on the Estimates from the RPLR
Using the timedependent parameter vector of the VTARMA model as estimated above, thus, the frozentime power spectra can be calculated by [20] as follows: where and are defined by and the covariance matrix of the residual vector is estimated by where is estimated by as shown (10) and is the length of the smoothing window.
4. Stage II: pLSCFBased FrequencyDomain Estimator for FrozenTime Modal Parameters
With the whitenoise assumption, one of the fundamentals in the community of the operational modal analysis, the power spectra can be expressed by the right matrix fraction model as the same form as the representation of the frequency response functions. Thus, using the power spectra at time instant , the corresponding modal parameters at that time instant can be estimated by the polyreference least square complex frequencydomain [18] (pLSCF) method as below.
The parametric model for estimation in the frequency domain is expressed by where with the number of outputs and the complexfrequencydomain matrices and are represented by the matrix polynomials as where is the order of the matrix fraction model and the complexfrequencydomain basis function is defined as the domain form as ( is the sampling interval). Thus the unknown parameters in the matrix fraction model can be rewritten as the vector form by
The cost function of the least square estimator is commonly defined by with the frequency sampling at and the equation errors are where are the measured spectra.
The matrix form of (20) can be rewritten by with
Thus, the cost function as shown in (19) can be rewritten as where , , and .
When minimizing the cost function as (23), the normal equation of the least squares is achieved as
Substituting (24) into (25), the reduced normal equation is obtained as
Furthermore, giving the appropriate constraint on the , its nontrivial solution can be conducted by referring to [21] and can be calculated by backsubstituting into (24). The computational complexity of this stage is approximately with the reference.
After the unknown parameter vector is achieved, the modal parameter, including the modal frequency, damping ratio, and the operational reference vector can be obtained by solving a general right eigenvalue problem as with where are the operational reference vector. The modal frequency and damping ratio are contained in the poles .
Using the pLSCF method as presented above, the modal parameters at the new broughtin time instant can be estimated.
5. Stage III: Recursive Validation of the Estimated Modal Parameters
After obtaining the modal parameters at the new broughtin time instant, the validation and order determination of modal parameters are necessary. In this stage, the modal parameters are validated and ordered recursively by using a recursive fuzzy clustering.
5.1. Fuzzy Clustering
Fuzzy cluster analysis or fuzzy clustering partition is a set of data into clusters and each cluster is represented by its center, also called prototype. In other words, the aim of the fuzzy cluster analysis is to determine the prototypes and the memberships.
Let be a data set with data points. The memberships can be represented by a matrix with clusters. An element in the membership matrix , with and , denotes the degree of the th data point belonging to the th cluster relative to other clusters. Because the membership matrix, , represents a probabilistic cluster partition, there are a few constraints on : Equations (29) and (30) imply that the total membership for the th data point to all clusters is 1 and, for each cluster, at least one data point belongs to that cluster with a nonzero membership, respectively.
Let be a set of prototypes with cluster prototypes. is the prototype of the th cluster.
The fuzzy means (FCM) [22] algorithm is one of the objectivebased fuzzy clustering algorithms. Given the data set , the membership matrix , the cluster prototypes and the distance function , a scalar objective function is described by where , the exponent is the fuzzy exponent, and denotes the distance between two vectors, which does not depend on the memberships. Commonly, and is the Euclidean distance function as
The fuzzy means algorithm (FCM) is one of the most popular algorithms of fuzzy cluster analysis. There are five phases in FCM.
Phase I. Choosing the number of clusters , the fuzzy exponent , the convergence criterion , and the maximum number of iterations , the initial membership matrix is defined randomly, which satisfies the constraints as (29), (30), and (31).
Phase II. Update the cluster prototypes, , at the th iteration by
Phase III. Calculate the distances for all and using the updated cluster prototypes, , based on the distance function as defined in (33) or other distance functions.
Phase IV. Update the membership matrix, , at the th iteration by where is the null set and . If contains more than one element, for is not uniquely determined and more operations are needed [22].
Phase V. Check the termination criteria including the convergence criterion and the maximum number of iterations . If or , the iteration terminates; otherwise, repeat Phase II, Phase III, and Phase IV.
5.2. Recursive Fuzzy Clustering
When the observed data changes with the time, the recursive clustering is necessary to capture the current features of the newcoming data. An approach of recursive fuzzy clustering is introduced as follows [23].
The th cluster prototype at the time instant is defined by . The relation between the old cluster prototype and a new one can be expressed by with the following increment:
As shown in (37), the calculation of the membership at the th time instant requires the past memberships, which is against the recursive idea. Therefore, an approximate calculation of the denominator of (37) is achieved by adding the forgetting factors for the past memberships defined as follows: where is the forgetting factor and . As shown in (38), the denominator at the time instant , , can be calculated by a recursive approach. Furthermore, the current membership is defined by
As presented above, the algorithm is recursive and no convergence condition is required.
5.3. Recursive Fuzzy Validation of the Modal Parameters
Furthermore, the memberships will be considered as the weights in the next continuoustime estimation of modal parameters using a recursive least square approach. The validated modal parameters are expressed into the four subpartitions with the corresponding memberships by where is general data point, which can be the modal frequency or the damping ratio, is the corresponding membership, and is the length of the th subpartition including the data of time and before , which is updated with process of recursive estimation.
6. Stage IV: Estimation of ContinuousTime Modal Parameters
In this stage, the continuoustimerepresented modal parameters are estimated by a recursive least square approach based on the validated modal parameters of time and before . The cost function of this recursive least squares is defined by where , , and are defined in (40) and with the following equation error: where are the weights and defined by ( is the exponent for weighting), is an arbitrary component in and the predict function is expressed by a polynomial as where and the basis function are polynomials.
The estimation of the projection coefficients can be obtained by minimizing the cost function, (41), as where “” means the “argument minimization.”
According to (42) and (43), the equation errors can be rewritten by
As defined in (45), the least square problem is linear. Hence, the projection coefficients, , can be estimated [24] as follows:
The can be estimated as a recursive behavior. Define the projection coefficients of the th time instant ; the can be calculated by with
7. Experiment
In this section, experimental results validate the proposed multistage recursive estimator. The reference modal parameters are obtained by identifying a series of “frozenconfiguration” structures of corresponding to the timevarying structure.
7.1. The Structure and the Laboratory Setup
Figure 3 shows the schematic diagram of the experimental structure and its laboratory setup. The structure consists of a steel beam and a steel mass, of which the parameters are listed in Table 1. The beam is simply supported at the two ends by a heavy frame. The mass can slide along the axial direction of the beam driven by a controllable motor from the center to 800 mm away. The velocity of the mass is constantly 100 mm/s under a feedback control of the motor. Eleven piezoelectric accelerometers (sensitivity: ≈100 mV/g) measure the acceleration responses of the beam at eleven uniformly distributed positions along the axial direction of the beam, although only four are shown in Figure 3(a). A LMS SCADAS III system acquires the acceleration response signals and the force signal of the shaker. The PC records the input signals via the SCADAS system and LMS Test.Lab software and starts the motor. The measurement is triggered by the signal from the motor shown as the green arrowed connection in Figure 3(a).

(a)
(b)
Figure 3(b) shows the real configuration of the structure and the laboratory setup. For the measurement of the baseline models, the input force is measured by an impedance head (sensitivity: 2.248 V/kN). For the measurement of the timevarying case, the input force is assumed as unknown.
7.2. Baseline Parametric Identification of the FixedConfiguration Structure via the General TimeInvariant Approach
The timevarying structure in this paper can also be considered as a linear parametervarying (LPV) structure. The structural dynamic characteristics are functions of the position of the mass. While other parameters are fixed, the position of the mass can be continuously timevarying. In other words, when the mass stays at a fixed position, the structural dynamic characteristics, such as modal parameters, can be obtained by the available LTI system identification approaches, which are the “frozentime” or frozenconfiguration reference for the timevarying case [4].
The frozenconfiguration structure is the baseline model for the timevarying structure. Moving the mass from the center of the beam to 800 mm away with a spacing of 10 mm, 81 baseline models can be obtained. The eleven acceleration signals and the force signal are measured with a sampling frequency of 512 Hz and record length of 2 s. Figure 4 shows the averaged FRFs (the bandwidth is 0–160 Hz, the average number is 5, and the FRFs are estimated by H1 method [25]) between the excitation and Outputs 1, 4, 7, and 11 for the 81 baseline models.
(a)
(b)
(c)
(d)
In the FRFs shown in Figure 4, there are five peaks with respect to frequency, which indicate the five modes in the bandwidth (0–160 Hz).
In order to identify the modal parameters of the baseline models, the least squares complex exponential method (LSCE) [25] is used. Figures 5(a) and 5(b) show the modal frequency and damping ratio of the baseline models, respectively.
(a)
(b)
7.3. TimeVarying Testing
In this section, the structure is timevarying due to the mass sliding continuously. The setup is the same as that for the “fixedconfiguration” structure. The sampling frequency is 512 Hz and the record length is 8 s. For the sake of comparison between the estimation results of the timevarying structure and the reference values of the baseline models, the relationship between the position of the mass and the time variable is achieved by using the known constant velocity of the mass.
The excitation generated by the shaker is Gaussian white noise. The acceleration is measured at 11 uniformly distributed locations, which forms the basic data set for the modal parameter estimation of the timevarying structure.
Considering the uniformity, Figures 6(a), 6(b), 6(c), and 6(d) show the responses of Outputs 1, 4, 7, and 11 of the 11 outputs.
(a)
(b)
(c)
(d)
7.4. Recursive Pseudolinear Regressive Estimation of Power Spectra
Once the response data are measured, the power spectra at time can be obtained. For the sake of presenting the results as the sequence of the stages, the whole batch of the results is shown in the following illustrations rather than as the sequence of the time instant by time instant, although the estimation is recursively processed by time.
Assume that the interested frequency band is from 70 Hz to 150 Hz; the cross power spectra for Outputs 1, 4, 7, and 11 referred to Output 7 are shown in Figure 7.
(a)
(b)
(c)
(d)
7.5. pLSCFBased FrequencyDomain Estimator for FrozenTime Modal Parameters
Based on the pLSCF estimator in the frequency domain and the power spectra (cross power spectra for Outputs 1–11 referred to Output 7), the modal parameters at time are estimated as shown in Figure 8.
(a) Resonance frequency
(b) Damping ratio
7.6. Recursive Validation of the Estimated Modal Parameters
This section validates the estimated modal frequency as shown in Figure 8 using the recursive fuzzy clustering introduced in Section 5. The clustering results are shown in Figure 9, in which the letters “A” and “B” denote the two clusters of the modal frequency, the color of these letters indicates the memberships of the modal frequency and the black dots show the prototypes of these clusters. The prototypes of the clusters move from left to right due to recursive behavior of the recursive fuzzy clustering algorithm. In this case, only the modal frequency is used for validation and the modal frequency before 0.5 s ( s) is considered as the start data for initializing the fuzzy clustering. In the initialization phase, the parameters for fuzzy clustering are and .
As shown in Figure 9, modal parameters are clustered and the memberships generally measure the quality of the corresponding modal parameters. The latter discriminates the lowquality estimates through weighting the data with the memberships in the next stage, as presented in (42).
7.7. Estimation of ContinuousTime Modal Parameters
Based on the validated and ordered modal parameters from the last stage, the continuoustime modal parameters are estimated by using the approach presented in Section 6. The results of the continuoustime modal frequency and damping ratio are shown in Figures 10(a) and 10(b). In these recursive least squares, is 7 and .
(a) Resonance frequency
(b) Damping ratio
As shown in Figure 10, the continuoustime estimation of the modal frequency coincides with the results of the baseline models and that of the damping ratio is at the same level as the results of the baseline models.
8. Conclusion
This paper proposes a vector multistage recursive model parameter estimator for timevarying structures in the hybrid time and frequency domain. The research reproduces the available approaches of power spectra estimation, frequencydomain estimator for timeinvariant structures, fuzzy clustering, and so forth, to present a recursive way of modal parameter estimation. This estimator has the capability of selecting the frequency band, which is helpful to achieve the parsimony of the estimation algorithm and may improve the accuracy and quality of the estimation in a limited frequency band. The recursive fuzzy clustering in the proposed estimator achieves the automatically recursive validation of the modal parameters at the current time instant. The final stage can obtain the continuoustime modal parameters through recursive least squares. The contribution of this paper is mainly the combination in the workflow itself and the necessary reproduction of each element in the proposed method.
In case the response data at time are measured, the modal parameters at that time instant are estimated, so it is recursive, which has the computational parsimony, and is potentially suitable for the realtime acquisition of the modal estimation for timevarying structures and is further possibly used in the applications such as the mode tracking, monitoring, and adaptive control of vibration.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors acknowledge the support for the work presented in this paper from the Beijing Institute of Technology Foundation for Basic Research (Grant no. 20120142009) and National Natural Science Foundation of China (Grant no. 11402022).
References
 K. A. Petsounis and S. D. Fassois, “Nonstationary functional series TARMA vibration modelling and analysis in a planar manipulator,” Journal of Sound and Vibration, vol. 231, no. 5, pp. 1355–1376, 2000. View at: Publisher Site  Google Scholar
 K. A. Petsounis and S. D. Fassois, “Parametric timedomain methods for the identification of vibrating structuresa critical comparison and assessment,” Mechanical Systems and Signal Processing, vol. 15, no. 6, pp. 1031–1060, 2001. View at: Publisher Site  Google Scholar
 A. G. Poulimenos and S. D. Fassois, “Parametric timedomain methods for nonstationary random vibration modelling and analysis—a critical survey and comparison,” Mechanical Systems and Signal Processing, vol. 20, no. 4, pp. 763–816, 2006. View at: Publisher Site  Google Scholar
 A. G. Poulimenos and S. D. Fassois, “Outputonly stochastic identification of a timevarying structure via functional series TARMA models,” Mechanical Systems and Signal Processing, vol. 23, no. 4, pp. 1180–1204, 2009. View at: Publisher Site  Google Scholar
 M. D. Spiridonakos and S. D. Fassois, “Parametric identification of a timevarying structure based on vector vibration response measurements,” Mechanical Systems and Signal Processing, vol. 23, no. 6, pp. 2029–2048, 2009. View at: Publisher Site  Google Scholar
 K. Liu, “Identification of linear timevarying systems,” Journal of Sound and Vibration, vol. 206, no. 4, pp. 487–505, 1997. View at: Publisher Site  Google Scholar  MathSciNet
 K. Liu, “Extension of modal analysis to linear timevarying systems,” Journal of Sound and Vibration, vol. 226, no. 1, pp. 149–167, 1999. View at: Publisher Site  Google Scholar
 K. Liu and L. Deng, “Experimental verification of an algorithm for identification of linear timevarying systems,” Journal of Sound and Vibration, vol. 279, no. 3–5, pp. 1170–1180, 2005. View at: Publisher Site  Google Scholar
 K. Liu and L. Deng, “Identification of pseudonatural frequencies of an axially moving cantilever beam using a subspacebased algorithm,” Mechanical Systems and Signal Processing, vol. 20, no. 1, pp. 94–113, 2006. View at: Publisher Site  Google Scholar
 R. Ghanem and F. Romeo, “A waveletbased approach for the identification of linear timevarying dynamical systems,” Journal of Sound and Vibration, vol. 234, no. 4, pp. 555–576, 2000. View at: Publisher Site  Google Scholar
 A. RoshanGhias, M. B. Shamsollahi, M. Mobed, and M. Behzad, “Estimation of modal parameters using bilinear joint timefrequency distributions,” Mechanical Systems and Signal Processing, vol. 21, no. 5, pp. 2125–2136, 2007. View at: Publisher Site  Google Scholar
 X. Xu, Z. Zhang, H. Hua, and Z. Chen, “Identification of timevarying modal parameters using linear timefrequency representation,” Chinese Journal of Mechanical Engineering, vol. 16, no. 4, pp. 445–448, 2003. View at: Publisher Site  Google Scholar
 E. Lourens, E. Reynders, G. De Roeck, G. Degrande, and G. Lombaert, “An augmented Kalman filter for force identification in structural dynamics,” Mechanical Systems and Signal Processing, vol. 27, no. 1, pp. 446–460, 2012. View at: Publisher Site  Google Scholar
 S.B. Song and C.M. Pei, “A recursive method of timefrequency analysis for the signal processing of flutter test with progression variable speed,” Chinese Journal of Aeronautics, vol. 18, no. 3, pp. 213–217, 2005. View at: Publisher Site  Google Scholar
 L. Ljung, System Identification: Theory for the User, Prentice Hall, Upper Saddle River, NJ, USA, 1999.
 N. Calinoiu, S. Bachir, and C. Duvanaud, “Mixed time and frequencydomain identification of radio frequency power amplifiers,” in Proceedings of the 40th European Microwave Conference (EuMC ’10), pp. 1018–1021, Paris, France, September 2010. View at: Google Scholar
 S.D. Zhou, W. Heylen, P. Sas, and L. Liu, “Parametric modal identification of timevarying structures and the validation approach of modal parameters,” Mechanical Systems and Signal Processing, vol. 47, no. 12, pp. 94–119, 2014. View at: Publisher Site  Google Scholar
 H. van der Auweraer, P. Guillaume, P. Verboven, and S. Vanlanduit, “Application of a faststabilizing frequency domain parameter estimation method,” Journal of Dynamic Systems, Measurement and Control, Transactions of the ASME, vol. 123, no. 4, pp. 651–658, 2001. View at: Publisher Site  Google Scholar
 H. van der Auweraer and B. Peeters, “Discriminating physical poles from mathematical poles in high order systems: use and automation of the stabilization diagram,” in Proceedings of the 21st IEEE Instrumentation and Measurement Technology Conference (IMTC ’04), vol. 3, pp. 2193–2198, Cornu, Italy, May 2004. View at: Publisher Site  Google Scholar
 G. C. Reinsel, Elements of Multivariate Time Series Analysis, Springer, Berlin, Germany, 1993.
 B. Cauberghe, P. Guillaume, P. Verboven, S. Vanlanduit, and E. Parloo, “On the influence of the parameter constraint on the stability of the poles and the discrimination capabilities of the stabilisation diagrams,” Mechanical Systems and Signal Processing, vol. 19, no. 5, pp. 989–1014, 2005. View at: Publisher Site  Google Scholar
 F. Höppner, F. Klawonn, R. Kruse, and T. Runkler, Fuzzy Cluster Analysis: Methods for Classification, Data Analysis and Image Recognition, John Wiley & Sons, 1999.
 D. Dovžan and I. Škrjanc, “Recursive fuzzy cmeans clustering for recursive fuzzy identification of timevarying processes,” ISA Transactions, vol. 50, no. 2, pp. 159–169, 2011. View at: Publisher Site  Google Scholar
 S. Van Huffel and J. Vandewalle, The Total Least Squares Problem: Computational Aspects and Analysis, Society for Industrial and Applied Mathematics, 1991. View at: Publisher Site  MathSciNet
 W. Heylen, S. Lammens, and P. Sas, Modal Analysis Theory and Testing, Katholieke Universiteit Leuven, Leuven, Belgium, 2007.
Copyright
Copyright © 2015 SiDa Zhou 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.