Shock and Vibration
Volume 2019 (2019), Article ID 9602535, 11 pages
https://doi.org/10.1155/2019/9602535

Research Article

Short Data-Based Output-Only Identification for Time-Varying Systems with Fast Dynamic Evolution

Zhi-Sai Ma1,2 and Qian Ding1,2

¹School of Mechanical Engineering, Tianjin University, Tianjin 300350, China
²Tianjin Key Laboratory of Nonlinear Dynamics and Control, Tianjin 300350, China

Correspondence should be addressed to Qian Ding; qding@tju.edu.cn

Received 13 March 2019; Revised 23 April 2019; Accepted 9 May 2019; Published 2 June 2019

Academic Editor: Itzhak Green

Copyright © 2019 Zhi-Sai Ma and Qian Ding. 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.

Abstract.

Many engineering systems change appreciably over a relatively short time interval due to their fast evolution in the dynamics. Time-varying (TV) system’s ambient excitation is usually difficult to measure under operating conditions, and its dynamics have to be determined without measuring the excitation. Therefore, short data-based output-only identification for TV systems with fast dynamic evolution is considered in this paper. Deterministic parameter evolution methods are known to track fast dynamic evolution by postulating TV model parameters as deterministic functions of time and selecting proper functional subspaces. However, these methods require a significant number of parameters to represent complicated time-dependencies and dynamics characterized by larger numbers of degrees-of-freedom. In such cases, the ordinary least squares estimation may lead to less accurate or even unreliable estimates. A ridge regression-based deterministic parameter evolution method is proposed to overcome ill-posed problems via regularization and subsequently assessed through numerical and experimental validation. Comparative results confirm the advantages of the proposed method in terms of achievable natural frequency and power spectral density tracking, accuracy, and resolution of TV systems with fast dynamic evolution, when the response data length is relatively short.

1. Introduction

The need for time-varying (TV) system characterization is pervasive throughout the aerospace, mechanical, transportation, and manufacturing engineering practice. It is usually difficult to build the explicit dynamic model of a practical system via mechanism analysis; hence, TV system identification in a way that takes time variation explicitly into account is receiving more and more attention [1–4]. With the development of engineering applications, systems with fast varying dynamics are being widely used and their intrinsic nonstationary characteristics are increasingly inevitable. On the one hand, these systems may change appreciably over a relatively short time interval due to their fast evolution in the dynamics. On the other hand, the dynamic characteristics of TV systems may not fit into a laboratory and controlled testing may not be feasible under operating conditions. In other words, TV system’s ambient excitation is usually difficult to measure under operating conditions and its dynamics have to be determined without measuring the excitation. Therefore, short data-based output-only identification for TV systems is considered in this work in order to capture their fast varying dynamics.

Identification methods for TV systems are generally classified under the umbrella of time-frequency methods [1]. Most frequency-domain methods employ time-frequency analysis (e.g., the short time Fourier transform, Cohen’s class, wavelet-based methods, the Hilbert–Huang transform, etc.), which are based on nonparametric representations of the nonstationary signal as a simultaneous function of time and frequency. By contrast, time-domain methods employ parametric representations (e.g., time-dependent autoregressive moving average (TARMA) models), which offer a number of potential advantages, such as improved accuracy, resolution, and tracking of the TV dynamics. TARMA model-based methods may be further divided into three classes according to the type of “structure” imposed upon the evolution of the TV model parameters: unstructured parameter evolution (UPE), stochastic parameter evolution (SPE), and deterministic parameter evolution (DPE) methods [3]. The DPE methods are mainly based on explicit models of parameter evolution through projecting the parameter trajectory onto specific functional subspaces defined by deterministic basis functions (such as trigonometric, Legendre, Chebyshev, wavelets, B-splines, moving Kriging shape functions, radial basis functions, and others). By postulating TV model parameters as deterministic functions of time, the TARMA model can be transformed into a functional series (FS) TARMA model with time-invariant parameter matrix. Therefore, DPE methods are known to track fast dynamic evolution by selecting proper functional subspaces over their UPE and SPE counterparts [2–4].

DPE methods have played an important role in the development of TV system identification and nonstationary signal processing. In 1970, Rao [5] approximated AR coefficients by truncated Taylor series expansion and introduced the idea of representing TV model parameters as polynomials of time. Thereafter, Mendel [6] summarized the nonstationary identification methods and classified them into UPE, SPE, and DPE, which are still in use today. Kozin [7] approximated AR coefficients by a linear combination of Legendre basis functions and applied the FS-TAR model to earthquake ground motion analysis. In 1983, Grenier [8] extended the FS-TAR model to the FS-TARMA case by examining the existence and uniqueness of a FS-TARMA model representation for a nonstationary system, which was successfully applied to speech signal modeling [9–12]. Gersh and Kitagawa [13] extended the univariate FS-TAR model to the multivariate/vector case. Niedzwiecki [14, 15] defined the idea of expanding TV model parameters onto linear combinations of basis functions as “functional series modeling” and proposed a recursive least squares for the estimation of FS-TAR model parameters. Tsatsanis and Giannakis [16, 17] proposed FS-TARMA models with exogenous input and FS time-dependent finite impulse response models and considered the issue of model structure selection. In 2000, Niedzwiecki [2] summarized the DPE methods and investigated the time and frequency characteristics of basis function estimators.

In 2006, Poulimenos and Fassois [3] systematically reviewed the UPE, SPE, and DPE methods in terms of achievable model parsimony, model parameter estimation accuracy, tracking capability of TV dynamics, and computational simplicity. Spiridonakos and Fassois extended the FS-TAR model to the FS-TARMA case [18] and proposed the adaptable FS-TARMA model characterized by basis functions that are adaptable to the signal being modeled [19]. In 2014, Spiridonakos and Fassois [4] reviewed the conventional and adaptable FS-TARMA model in terms of parameter estimation, model structure selection, and model validation and outlined the developments of the DPE methods. Recently, Yang et al. [20] have proposed an adaptable FS-TARMA modeling method by selecting moving Kriging shape functions as basis functions. Zhou et al. [21] have proposed a least squares support vector machine-based FS-TAR modeling method by selecting compactly supported radial basis functions. Bertha and Golinval [22] have focused on the identification of TV mode shapes during the FS-TARMA modeling. Li et al. [23] have focused on the problem of extracting the real modal parameters from computational ones during the adaptable FS-TAR modeling. Ma et al. [24, 25] have focused on the recursive identification of TV systems by introducing the idea of kernel methods into the FS-TARMA modeling.

Despite the foregoing progress that has been achieved so far, there are still many issues that are open and require work, for example, the feasibility of TV system identification under a relatively limited data length. DPE methods can achieve better model parsimony over their UPE and SPE counterparts, yet they still require a significant number of parameters to represent complicated time-dependencies and dynamics characterized by larger numbers of degrees-of-freedom [3]. Therefore, the aim of the paper is two-fold: (1) the existing least squares-based FS-TAR modeling is reviewed to reveal the problem of TV system identification under a relatively limited data length and (2) a ridge regression-based FS-TAR method is subsequently proposed to track fast dynamic evolution of the system being modeled by using short data.

The paper is organized as follows: the least squares-based FS-TAR modeling is presented in Section 2. Section 3 proposes a new ridge regression-based FS-TAR method to identify TV systems with fast dynamic evolution under limited data length. The proposed and existing methods are numerically and experimentally tested in Section 4 and Section 5, respectively. Finally, Section 6 gives the remarkable conclusions.

2. LS-Based FS-TAR Model Parameter Estimation

2.1. FS-TAR Model

A time-dependent autoregressive (TAR) model with designating its autoregressive (AR) order is given by [18]

where designates normalized discrete time, the discrete-time nonstationary vibration signal being modeled, an unobservable uncorrelated innovations (residual) sequence with zero mean and TV nonsingular covariance matrix , and the time-dependent AR parameter matrix. normally independently distributed random variables with the indicated mean and covariance.

Any given function can be approximated by a set of selected basis functions with arbitrary accuracy, as long as a sufficient number of basis functions is used [26]. Therefore, a TV parameter can be approximated by

where denotes a set of basis functions selected from a suitable functional subspace, the coefficient of projection, the subspace dimensionality, and the index the specific basis functions of the selected functional subspace that are included in the basis.

Therefore, the time-dependent AR parameter matrix can be expressed as [3, 4, 18, 20]

The TAR model can be rewritten into the following FS-TAR model by substituting equation (3) into equation (1):

with

being the time-invariant parameter matrix, and

the corresponding regression vector, where is the Kronecker product [27], is given by

Evidently, the time-varying TAR model of equation (1) is transformed into a time-invariant FS-TAR model by expanding the time-dependent model parameters onto the selected functional subspaces. As given by equation (5), the FS-TAR model is fully described by parameters (coefficients of projection).

2.2. Least Squares Estimation

Suppose the length of measured response data is , the matrix equation can be obtained by assembling each equation of form (4) for , as follows:

with , , , and .

The parameter matrix can be estimated by using the well-known least-squares estimator through minimizing the sum of squared residuals, as follows:

where denotes the Moore–Penrose pseudo inverse. The time-dependent AR parameter matrix can be subsequently estimated by using equation (3), after substituting by the obtained estimate .

However, if measured response data length is short and the inequality holds, there are fewer equations than unknown parameters and more than one satisfies equation (8). In other words, the solution is not unique and the problem is ill posed. In such cases, the ordinary least squares estimation is sensitive to errors and may lead to unreliable estimates. From the standpoint of the FS-TAR modeling, the data can be considered as “short” if the sample length () is smaller than the number of unknown parameters of the identified model ().

3. RR-Based FS-TAR Model Parameter Estimation

3.1. Ridge Regression Estimation

Ridge regression, also known as Tikhonov regularization, is a commonly used method of regularization of ill-posed problems [28, 29]. Given the regularization parameter , the least squares cost function is defined as

By using the ridge regression, we have the solution to equation (10) as

where and .

By using the matrix inversion lemma, we have

Obviously, the computational complexity of the matrix inversion in equations (11) and (12) is, respectively, and . In order to avoid the matrix inversion, we further define and as

Then, we have

The update rule for the inverse of the above growing matrix is

with and .

By combining equations (12), (13), and (16), we have

with the final estimate . Similarly, the time-dependent AR parameter matrix can be subsequently estimated by using equation (3), after substituting by the obtained estimate .

3.2. FS-TAR Model-Based Structural Dynamics

Once have been estimated, the system’s “frozen-time” poles can be obtained by solving the following general eigenvalue problem [20]:

where designates the identity matrix of order and and , respectively, the poles and eigenvector of at time instant with the mode shape . The matrix is constructed from as

The “frozen-time” natural frequencies of the system can be computed by

where designates the absolute value and the sampling period.

The response “frozen-time” power spectral density (PSD) for each time instant can be obtained as

where designates frequency in , the imaginary unit, the Hermitian transpose, and , with being the window length.

4. Numerical Validation

The performance characteristics of the proposed RR-FS-TAR method, along with the LS-FS-TAR method, are examined via Monte Carlo study focused on the identification of a TV three degrees-of-freedom structural system subject to a Gaussian, zero-mean, and uncorrelated excitation , as shown in Figure 1 [3, 30]. Numerical values of the system parameters are given in Table 1.

Figure 1: Schematic diagram of a TV structural system.

Table 1: Numerical values of the system parameters.

Property
Value

Mass (kg)

Damping (N/(m/s))

Stiffness (N/m)

Three cases are considered by choosing different values of the scaling factor . Case 1: slow dynamic evolution with , Case 2: medium dynamic evolution with , and Case 3: fast dynamic evolution with . The simulation time is, respectively, 425, 85, and 17 seconds for each case. Given a particular excitation realization, the displacement responses are obtained through the Runge–Kutta 4/5 method characterized by variable integration step and subsequently recorded at a sampling frequency 12 Hz. Therefore, the length of the samples is, respectively, 5100, 1020, and 204 for each case. The above procedure is repeated for each Monte Carlo test (number of tests: 50), with a different excitation realization used each time.

4.1. Case 1: Slow Dynamic Evolution

In this case, the displacement responses are 425s long and the sample length is . The AR order is manually selected as because the actual number of modes present is and , as given by equation (20). The functional basis spanned by trigonometric functions of the form [3, 30],

is selected. The subspace dimensionality and index are selected to approximate the TV model parameters. Figure 2 depicts the LS-based and RR-based natural frequency estimates from the 50 Monte Carlo tests and their true counterparts. As shown in Figure 2 the TV natural frequencies are adequately tracked by the proposed RR-FS-TAR method and the LS-FS-TAR method, due to the relatively larger sample length. However, the proposed RR-FS-TAR method seems to attain better overall performance in terms of tracking and accuracy as the LS-FS-TAR method-based natural frequency estimates exhibit some local tracking errors.

(a)

(b)

Figure 2: Natural frequency estimates by LS (a) and RR (b) for the TV numerical system with slow dynamic evolution (gray dots: natural frequency estimates from the 50 Monte Carlo tests; red dots: natural frequency estimates from the single test; green dashed lines: true values of natural frequencies).

4.2. Case 2: Medium Dynamic Evolution

In this case, the displacement responses are 85s long and the sample length is . The AR order is selected as . The functional basis spanned by trigonometric functions of equation (22) is here selected with the same values of subspace dimensionality and index . Figure 3 depicts the LS-based and RR-based natural frequency estimates from the 50 Monte Carlo tests and their true counterparts. Evidently, the TV natural frequencies are not adequately tracked by the LS-FS-TAR method and their estimates are quite rough. The trend of the natural frequency estimates obtained by the LS-FS-TAR method is in overall agreement with the true values, but they are of rather poor quality in terms of tracking and accuracy. By contrast, the proposed RR-FS-TAR method attains much better overall performance in terms of tracking and accuracy, even if the system changes faster and the available sample length is relatively shorter.

(a)

(b)

Figure 3: Natural frequency estimates by LS (a) and RR (b) for the TV numerical system with medium dynamic evolution (gray dots: natural frequency estimates from the 50 Monte Carlo tests; red dots: natural frequency estimates from the single test; green dashed lines: true values of natural frequencies).

4.3. Case 3: Fast Dynamic Evolution

In this case, the system changes appreciably over a relatively short time interval. The displacement responses are 17s long and the sample length is . The AR order is selected as . The functional basis spanned by trigonometric functions of equation (22) is here selected with the same values of subspace dimensionality and index . Figure 4 depicts the LS-based and RR-based natural frequency estimates from the 50 Monte Carlo tests and their true counterparts. The sample length () is much smaller than the number of unknown parameters (). In such case, the ordinary least squares estimation is sensitive to errors and may lead to unreliable estimates. This behavior is confirmed by Figure 4, where the TV natural frequencies are not tracked by the LS-FS-TAR method and their estimates are rather vaguely shown or even missing. The proposed RR-FS-TAR method is able to deal with the ill-posed problem via regularization. As shown in Figure 4, the natural frequency estimates obtained by the RR-FS-TAR method is in overall agreement with the true values, although they are rough and exhibit larger tracking errors.

(a)

(b)

Figure 4: Natural frequency estimates by LS (a) and RR (b) for the TV numerical system with fast dynamic evolution (gray dots: natural frequency estimates from the 50 Monte Carlo tests; red dots: natural frequency estimates from the single test; green dashed lines: true values of natural frequencies).

4.4. Power Spectral Density Estimates

The superior performance of the proposed RR-FS-TAR method is pronounced by the comparison of Figures 2–4 where the estimated natural frequencies are shown. In this section, we further examine the enhanced capability of the RR-FS-TAR method in terms of achievable power spectral density tracking and resolution. Figure 5 depicts the LS-based and RR-based power spectral density estimates from the single test for the TV numerical system with slow, medium, and fast dynamic evolution. The LS-based and RR-based power spectral density estimates are quite similar for the TV numerical system with slow dynamic evolution, as shown in Figure 5. However, with the decrease of the sample length, the LS-FS-TAR method obtains worse results, while the proposed RR-FS-TAR method attains much better performance in terms of tracking and resolution, even if the system changes very fast and the available sample length is quite short (). In contrast to LS-based FS-TAR estimator requiring larger sample length to represent complicated time-dependencies and dynamics, RR-based FS-TAR estimator is able to track fast dynamic evolution of the system being modeled by using short data.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5: Power spectral density estimates by LS (a, c, e) and RR (b, d, f) for the TV numerical system with slow, medium, and fast dynamic evolution (from top to bottom).

5. Experimental Validation

5.1. Experimental System

The experimental system [25, 31] is composed of the TV structure, the exciter system (an exciter and a power amplifier), force and motion transducers (an impedance head and 15 accelerometers), the measurement and acquisition module (LMS SCADAS III system), control systems (two Faulhaber DC motors and corresponding motion controllers), and boundary conditions, as shown in Figure 6. The TV structure is a coupled system consisting of a simply supported beam and a moving mass. The mass, driven by two DC motors, can slide on the beam. In this section, two cases are considered by adjusting the moving speed of the mass, as follows: Case 1: the speed is 0.130 m/s and Case 2: the speed is 0.325 m/s.

(a)

(b)

Figure 6: Schematic diagram of the experimental system and photo of its setup.

During the test, the mass starts at 0.5 m and stops at 1.8 m. Therefore, the durations for the mass to move over 1.3 m are, respectively, 10 and 4 seconds. The shaker is switched off and removed, and the TV structure is only excited by the motion of the moving mass. Obviously, the interaction between the moving mass and the beam is difficult to measure. Therefore, the excitation of the experimental system is unobservable and the system’s dynamic characteristics have to be determined by exclusively using measured vibration responses. The vibration responses of the simply supported beam are measured by accelerometers, placed at 15 uniformly distributed positions along the beam from left to right, as shown in Figure 6. The acceleration responses are measured at a frequency 512 Hz and subsequently resampled at 256 Hz. Therefore, the length of the samples is, respectively, 2560 and 1024 for the two cases. The above procedure is repeated for each Monte Carlo test (number of tests: 30), and the TV structure undergoes the same variation in each test.

5.2. Identification Results and Discussions

The proposed RR-FS-TAR method and the LS-FS-TAR method are experimentally compared and assessed based on the nonstationary acceleration measurements of the laboratory TV structure. The AR order is manually selected as , and the functional basis spanned by trigonometric functions of equation (22) is selected with the subspace dimensionality and index . Obviously, the sample length ( or ) is much smaller than the number of unknown parameters () and the least squares estimation may lead to unreliable estimates.

Figure 7 depicts the LS-based and RR-based natural frequency estimates from the 30 Monte Carlo tests when . The first four baseline natural frequencies of the TV structure are obtained by using the “frozen-time” technique [31] and denoted by green dashed lines in Figure 7. The gray dots denote the natural frequency estimates from the 30 tests, and the red dots denote the natural frequency estimates from the single test. Evidently, the TV natural frequencies are not adequately tracked by the LS-FS-TAR method and their estimates are quite rough, while the proposed RR-FS-TAR method attains much better tracking performance and the corresponding natural frequency estimates exhibit less false modes. The LS-based and RR-based power spectral density estimates from the single test for the TV structure are also shown in Figure 7. Compared to the LS-based power spectral density estimates, the RR-based estimates appear smoother and are of much better quality in terms of tracking and resolution.

(a)

(b)

(c)

(d)

Figure 7: Natural frequency and power spectral density estimates by LS (a, c) and RR (b, d) for the laboratory TV structure when (gray dots: natural frequency estimates from the 30 Monte Carlo tests; red dots: natural frequency estimates from the single test; green dashed lines: baseline values of natural frequencies).

Similarly, Figure 8 depicts the LS-based and RR-based natural frequency estimates from the 30 Monte Carlo tests and power spectral density estimates from the single test when . On the one hand, with the decrease of the sample length from 2560 to 1024, both the LS-FS-TAR and RR-FS-TAR methods obtain worse results. Compared to the identification results in Figure 7, natural frequency estimates exhibit larger overall scatter and power spectral density estimates appear more rough, as shown in Figure 8. On the other hand, the enhanced capability of the proposed RR-FS-TAR method is further pronounced to track fast dynamic evolution of the system based on short data. Compared to the LS-based identification results, the RR-based natural frequency estimates exhibit lower overall scatter, less false modes, and smaller tracking errors. Furthermore, the RR-based power spectral density estimates are quite smooth and the ridges are in good agreement with the baseline values, even if the structure changes very fast and the available sample length is relatively short.

(a)

(b)

(c)

(d)

Figure 8: Natural frequency and power spectral density estimates by LS (a, c) and RR (b, d) for the laboratory TV structure when (gray dots: natural frequency estimates from the 30 Monte Carlo tests; red dots: natural frequency estimates from the single test; green dashed lines: baseline values of natural frequencies).

For practical systems with larger numbers of degrees-of-freedom and complicated TV dynamics, the number of unknown parameters () of the FS-TAR model is usually quite large. In the meantime, the measured response data length () is limited, when systems change appreciably over a relatively short time interval. In such cases, the ordinary least squares estimation may lead to less accurate or even unreliable estimates, and the proposed ridge regression-based FS-TAR method is suggested to identify these systems with fast dynamic evolution based on short data.

6. Conclusions

The FS-TAR model can achieve good representation parsimony, yet it requires a significant number of parameters to represent complicated time-dependencies and dynamics characterized by larger numbers of degrees-of-freedom. Therefore, relatively longer data are usually needed to overcome ill-posed problems and ensure the reliability of least square estimation. However, in many practical cases, systems change appreciably over a relatively short time interval and their dynamics have to be determined based on limited measurement data. In this work, a short data-based output-only identification method for TV systems with fast dynamic evolution was proposed via regularization of ill-posed problems. The proposed ridge regression-based FS-TAR method was employed to identify numerical and experimental systems with different dynamic evolution and compared against the least squares-based FS-TAR method via Monte Carlo study. Comparison results confirm the advantages of the proposed method in terms of achievable natural frequency and power spectral density tracking, accuracy, and resolution of TV systems with fast dynamic evolution, when the available measurement data length is relatively short.

Data Availability

The data used to support the findings of this study are available from the authors upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors acknowledge the support from the National Natural Science Foundation of China (grant no. 11802201) and China Postdoctoral Science Foundation (grant no. 2017M621075).

References

L. Garibaldi and S. Fassois, “MSSP special issue on the identification of time varying structures and systems,” Mechanical Systems and Signal Processing, vol. 47, no. 1-2, pp. 1-2, 2014.

M. Niedzwiecki, Identification of Time-varying Processes, John Wiley & Sons, New York, NY, USA, 2000.

A. G. Poulimenos and S. D. Fassois, “Parametric time-domain methods for non-stationary random vibration modelling and analysis—a critical survey and comparison,” Mechanical Systems and Signal Processing, vol. 20, no. 4, pp. 763–816, 2006.

M. D. Spiridonakos and S. D. Fassois, “Non-stationary random vibration modelling and analysis via functional series time-dependent ARMA (FS-TARMA) models—a critical survey,” Mechanical Systems and Signal Processing, vol. 47, no. 1-2, pp. 175–224, 2014.

T. S. Rao, “The fitting of non-stationary time-series models with time-dependent parameters,” Journal of the Royal Statistical Society: Series B (Methodological), vol. 32, no. 2, pp. 312–322, 1970.

J. M. Mendel, Discrete Techniques of Parameter Estimation: The Equation Error Formulation, Marcel Dekker, New York, NY, USA, 1973.

F. Kozin, “Estimation and modeling of nonstationary time series,” in Proceedings of the Symposium of Applications of Computer Mathod in Engineering, pp. 603–612, Los Angeles, CA, USA, August 1977.

Y. Grenier, “Time-dependent ARMA modeling of nonstationary signals,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 31, no. 4, pp. 899–911, 1983.

Y. Grenier, “Estimation of non-stationary moving average models,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP ‘83), pp. 268–271, Boston, MA, USA, April 1983.

Y. Grenier, “Time-frequency analysis using time-dependent ARMA models,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP ‘84), pp. 270–273, San Diego, CA, USA, January 1984.

M. Chevalier, C. Chollet, and Y. Grenier, “Speech analysis and restitution using time-dependent autoregressive models,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP ‘85), pp. 501–504, Tampa, FL, USA, 1985.

Y. Grenier, “Nonstationary ARMA models via simultaneous AR and MA estimation,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP ‘86), pp. 2339–2342, Tokyo, Japan, 1986.

W. Gersh and G. Kitagawa, “A time-varying multivariate autoregressive modeling of econometric time series,” in Proceedings of the Business and Statics ASA Annual Meeting, pp. 399–404, Toronto, Canada, 1983.

M. Niedzwiecki, “Functional series modeling approach to identification of nonstationary stochastic systems,” IEEE Transactions on Automatic Control, vol. 33, no. 10, pp. 955–961, 1988.

M. Niedzwiecki, “Recursive functional series modeling estimators for identification of time-varying plants-more bad news than good?” IEEE Transactions on Automatic Control, vol. 35, no. 5, pp. 610–616, 1990.

M. K. Tsatsanis and G. B. Giannakis, “Time-varying system identification and model validation using wavelets,” IEEE Transactions on Signal Processing, vol. 41, no. 12, pp. 3512–3523, 1993.

M. K. Tsatsanis and G. B. Giannakis, “Modelling and equalization of rapidly fading channels,” International Journal of Adaptive Control and Signal Processing, vol. 10, no. 23, pp. 159–176, 1996.

M. D. Spiridonakos and S. D. Fassois, “Parametric identification of a time-varying structure based on vector vibration response measurements,” Mechanical Systems and Signal Processing, vol. 23, no. 6, pp. 2029–2048, 2009.

M. D. Spiridonakos and S. D. Fassois, “Adaptable functional series TARMA models for non-stationary signal representation and their application to mechanical random vibration modeling,” Signal Processing, vol. 96, pp. 63–79, 2014.

W. Yang, L. Liu, S.-D. Zhou, and Z.-S. Ma, “Moving kriging shape function modeling of vector TARMA models for modal identification of linear time-varying structural systems,” Journal of Sound and Vibration, vol. 354, pp. 254–277, 2015.

S.-D. Zhou, Y.-C. Ma, L. Liu, J. Kang, Z.-S. Ma, and L. Yu, “Output-only modal parameter estimator of linear time-varying structural systems based on vector TAR model and least squares support vector machine,” Mechanical Systems and Signal Processing, vol. 98, pp. 722–755, 2018.

M. Bertha and J.-C. Golinval, “Identification of non-stationary dynamical systems using multivariate ARMA models,” Mechanical Systems and Signal Processing, vol. 88, pp. 166–179, 2017.

W. Li, V.-H. Vu, Z. Liu, M. Thomas, and B. Hazel, “Application of adaptable functional series vector time-dependent autoregressive model for extraction of real modal parameters for identification of time-varying systems,” Measurement, vol. 103, pp. 143–156, 2017.

Z.-S. Ma, L. Liu, S.-D. Zhou, and L. Yu, “Output-only modal parameter recursive estimation of time-varying structures via a kernel ridge regression FS-TARMA approach,” Shock and Vibration, vol. 2017, Article ID 8176593, 14 pages, 2017.

Z.-S. Ma, L. Liu, S.-D. Zhou et al., “Parametric output-only identification of time-varying structures using a kernel recursive extended least squares TARMA approach,” Mechanical Systems and Signal Processing, vol. 98, pp. 684–701, 2018.

G. G. Walter, Wavelets and Other Orthogonal Systems with Applications, Taylor & Francis, Boca Raton, FL, USA, 1994.

H. Lütkepohl, Handbook of Matrices, John Wiley & Sons, New York, NY, USA, 1996.

A. N. Tikhonov, A. V. Goncharsky, V. V. Stepanov, and A. G. Yagola, Numerical Methods for the Solution of Ill-Posed Problems, Springer Science & Business Media, Dordrecht, Netherlands, 1995.

A. E. Hoerl and R. W. Kennard, “Ridge regression: biased estimation for nonorthogonal problems,” Technometrics, vol. 12, no. 1, pp. 55–67, 1970.

A. Poulimenos, M. Spiridonakos, and S. Fassois, “Parametric time-domain methods for non-stationary random vibration identification and analysis: an overview and comparison,” in Proceedings of the ISMA2006 International Conference on Noise and Vibration Engineering, pp. 2885–2905, Leuven, Belgium, September 2006.

Z.-S. Ma, L. Liu, S.-D. Zhou, and W. Yang, “Modal parameter estimation of the coupled moving-mass and beam time-varying system,” in Proceedings of the ISMA2014 International Conference on Noise and Vibration Engineering, pp. 587–596, Leuven, Belgium, 2014.