Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 723152, 16 pages

http://dx.doi.org/10.1155/2014/723152

Research Article

Mathematical Modeling for Lateral Displacement Induced by Wind Velocity Using Monitoring Data Obtained from Main Girder of Sutong Cable-Stayed Bridge

Gao-Xin Wang and You-Liang Ding

The Key Laboratory of Concrete and Prestressed Concrete Structures of Ministry of Education, Southeast University, Nanjing 210096, China

Received 2 April 2014; Accepted 22 May 2014; Published 19 June 2014

Academic Editor: Giuseppe Rega

Copyright © 2014 Gao-Xin Wang and You-Liang 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.

Based on the health monitoring system installed on the main span of Sutong Cable-Stayed Bridge, GPS displacement and wind field are real-time monitored and analyzed. According to analytical results, apparent nonlinear correlation with certain discreteness exists between lateral static girder displacement and lateral static wind velocity; thus time series of lateral static girder displacement are decomposed into nonlinear correlation term and discreteness term, nonlinear correlation term of which is mathematically modeled by third-order Fourier series with intervention of lateral static wind velocity and discreteness term of which is mathematically modeled by the combined models of ARMA and EGARCH. Additionally, stable power spectrum density exists in time series of lateral dynamic girder displacement, which can be well described by the fourth-order Gaussian series; thus time series of lateral dynamic girder displacement are mathematically modeled by harmonic superposition function. By comparison and verification between simulative and monitoring lateral girder displacements from September 1 to September 3, the presented mathematical models are effective to simulate time series of lateral girder displacement from main girder of Sutong Cable-Stayed Bridge.

1. Introduction

Nowadays, long span cable-stayed and suspension bridge structures are commonly constructed at home and abroad. On account of their flexible structural characteristics, displacement response from main girder of long-span bridge structure swings obviously impacted by strong aerostatic and fluctuating wind actions. According to aerostatic response analysis on Sutong Cable-Stayed Bridge by Xu et al., the lateral displacement response from main girder can approach 1.2 m under strong wind velocity 40 m/s with attack angle 0° [1]; and research results from buffeting response analysis on Golden Gate Bridge by Vincent showed that extreme buffeting amplitude from main girder can reach 1.7 m under strong wind velocity 31 m/s [2]. Such large amplitude can definitely threaten comfort and safety of the whole bridge structure. For example, severe wind vibration from main girder of Tacoma Suspension Bridge in Washington state eventually brought about collapse of the whole bridge structure under wind velocity 19 m/s [3]. Therefore, it is of great significance to research displacement response impacted by wind loads from main girder of long span bridge structures, and especially the lateral displacement response, as one fairly important part for main girder, should be specifically valued.

Theoretical exploration, numerical simulation, and wind tunnel tests for lateral displacement response have been carried out to some extent. Cheng and Xiao improved the calculation method for aerostatic stability and further concluded that instable lateral displacement was 4.24 m under critical static wind [4]; Long et al. analyzed lateral displacement response from Sidu Suspension Bridge through ANSYS finite element simulation and concluded that maximum lateral displacement at middle span was 32.26 cm, which is close to 1/1000 length of main span [5]; Yu et al. researched lateral displacement response from Xihoumen Suspension Bridge through wind tunnel tests and concluded that lateral displacement at horizontal angle 10° was larger than that of other angles [6].

However, for mechanism complexity of lateral displacement response impacted by aerostatic and fluctuating wind actions, traditional methods of theoretical deduction, numerical simulation, and wind tunnel tests are difficult to accurately reflect the actual lateral displacement response of bridge structure, on account of uncertain boundary condition, imprecise assignment of initial parameters, and inappropriate ignorance of subordinate factors. In recent years, with development of structural health monitoring technology, it is feasible to install monitoring sensors on long span bridge structures, monitoring data of which can authentically reflect bridge structural behaviors under actual environment and load actions. Although wind field of long span bridge structures has been widely monitored in recent years [7–9], lateral displacement response is rarely monitored and researched; thus real correlation regularity between lateral displacement response and wind action is still covered. Additionally, lateral displacement response under actual operation environment is also affected by other random factors, which have never been taken into account by researchers before. Therefore, lateral displacement response from main girder is necessarily researched upon monitoring data to reveal real structural behavior of long span bridges.

In this paper, based on the health monitoring system installed on the main span of Sutong Cable-Stayed Bridge, GPS displacement and wind field are real-time monitored and analyzed. According to analytical results, apparent nonlinear correlation with certain discreteness exists between lateral static girder displacement and lateral static wind velocity; thus time series of lateral static girder displacement are decomposed into nonlinear correlation term and discreteness term, nonlinear correlation term of which is mathematically modeled by th-order Fourier series with intervention of lateral static wind velocity and discreteness term of which is mathematically modeled by the combined models of ARMA and EGARCH. Additionally, stable power spectrum density exists in time series of lateral dynamic girder displacement; thus time series of lateral dynamic girder displacement are mathematically modeled by harmonic superposition function. By comparison and verification between simulative and monitoring lateral displacements from September 1 to September 3, mathematical models are feasible and effective to simulate time series of lateral girder displacement from main girder of Sutong Cable-Stayed Bridge.

2. Bridge Monitoring and Sample Analysis

The bridge monitoring object for this research is the worldwide famous Sutong Cable-Stayed Bridge (in Jiangsu Province, China). Its whole structure form is single-spanned and double-hinged with the main span reaching 1088 m as shown in Figure 1, and the main girder employs flat steel box type with 36.3 m wide and 4.0 m high as shown in Figure 2. 3D ultrasonic anemometers and GPS monitoring station are installed on two flanks of midspan cross-section from main girder (resp., shown in Figures 1 and 2) to continuously acquire wind data and displacement data with sample frequency of 1 Hz.

Figure 1: Longitudinal layout of monitoring equipment on the Sutong Bridge (unit: m).

Figure 2: Transverse layout of monitoring equipment on the flat steel box girder (unit: mm). (1) GPSMS: GPS monitoring station; (2) 3DUA: 3D ultrasonic anemometer.

Specifically, the wind data from 3D ultrasonic anemometers embrace such three types as wind velocity, horizontal angle, and vertical angle in local coordinate system (Figure 3), and the girder displacement data from GPS monitoring station contains absolute locations in WGS-84 coordinate system (Figure 3), which are supposed to deduct reference locations for analysis. Until now, the storage amount of monitoring data has increased to 93 million for each measurement point. Such considerable monitoring data cannot be totally applied to actual analysis; thus the monitoring data from upstream flank in the year 2012 are specially chosen.

Figure 3: Description of the defined local and WGS-84 coordinate system. Notes: (1) local coordinate system: ; (2) WGS-84 coordinate system: ; (3) 90° of wind horizontal angle denotes the axis angle from to (within range of [0° 360°]); (4) 90° of wind vertical angle denotes axis angle from to (within range of [−90° 90°]); (5) 10.6° axis angle exists between axes and of across the bridge.

Taking the monitoring data in the whole August, for example, considering that lateral displacement effect of main girder is primarily reflected by wind load across the Sutong Bridge, hence time series of wind velocity and girder displacement are decomposed into the -axis in local coordinate system as shown in Figures 4(a) and 5(a), respectively, which can be obviously observed that either of two time series contain static variant trend in whole and dynamic stochastic fluctuation in part. Such two kinds of variation characteristics can be furthermore separated by 10-minute average process as shown in Figures 4(b), 4(c), 5(b), and 5(c), respectively.

(a) Time series of original wind velocity

(b) Time series of static wind velocity

(c) Time series of dynamic wind velocity

Figure 4: Time series of wind velocity along -axis in the whole August.

(a) Time series of girder displacement

(b) Time series of static girder displacement

(c) Time series of dynamic girder displacement

Figure 5: Time series of girder displacements along -axis in the whole August.

By comparison between Figures 4(b) and 5(b), similar variation characteristics exist between static wind velocity and static girder displacement, which can be visually described by correlation scatter plots as shown in Figure 6, indicating apparent nonlinear correlation similar to quadratic parabolic curve. Therefore, static girder displacement can be mathematically expressed by static wind velocity, with consideration of definite discreteness affected by other random factors. Moreover, time series of dynamic girder displacement depict obvious steady stochastic fluctuation, as well as no variation of its power spectrum densities by time shown in Figure 7; thus dynamic girder displacement can be mathematically simulated by harmonic superposition method.

Figure 6: Correlation scatter plots between static wind velocity and static girder displacement.

Figure 7: Power spectrum densities of dynamic girder displacement in different periods.

3. Modeling Theory and Procedure

3.1. Modeling for Static Girder Displacement

3.1.1. Modeling Theory and Method

Correlation scatter plots in Figure 6 can be acquired by combination of nonlinear correlation term and discreteness term. Furthermore, nonlinear correlation term can be mathematically modeled by the -order Fourier series; that is,

where denotes time series of fitting static displacement for nonlinear correlation term, denotes time series of static wind velocity (Figure 4(b)), and , and are fitting parameters of th-order Fourier series.

With denotation of original static girder displacement as (Figure 5(b)), minus acquires time series of discreteness , which contains three types of stochastic characteristics (autoregression, moving average, and heteroscedasticity) and can be mathematically described by the combined models of ARMA and EGARCH. In detail, the ARMA model defines the stochastic characteristics of autoregression and moving average as follows:

where is the constant term, and , respectively, denote the orders of autoregression or moving average of , and , respectively, denote the coefficients of autoregression or moving average of with , and denotes the innovations process with time delay of . Meanwhile, the other EGARCH model defines the stochastic characteristic of heteroscedasticity as follows [9, 10]:

where denotes the conditional variance of the innovations process is the constant term, and , respectively, denote the orders of the EGARCH model, , and , respectively, denote the coefficients of the EGARCH model, is a standard, independent, and identically distributed random draw from some specified probability distribution such as Gaussian or Student’s , and

with degrees of freedom . Considering that the EGARCH model is treated as ARMA models for , thus the stationarity constraint for the EGARCH model is included by ensuring that the eigenvalues of the characteristic polynomial,

are inside the unit circle, where and are the variable and coefficients of the characteristic polynomial, respectively.

During modeling process for time series of discreteness , the autocorrelation function and partial correlation function with their lag phase are introduced for stationary test of . Specifically, the autocorrelation function can be calculated as follows [11]:

where

with denoting the amount of . And the other partial correlation function can be calculated through fitting successive autoregressive models of orders by ordinary least squares, retaining the last coefficient of each regression [12]. Besides, the AIC and BIC delimitation criteria are applied to determine model orders, with their statistical values and being, respectively, calculated as follows [13, 14]:

where denotes the optimized log-likelihood objective function (LLF) values associated with parameter estimates of the combined models, denotes the number of estimated parameters associated with each value in LLF, and denotes the sample size of the observed associated with each LLF value.

3.1.2. Detailed Procedure

Based on the theory and method above, the detailed procedure for mathematically modeling time series of static girder displacement is illustrated, taking the correlation scatter plots in the whole August in Figure 6, for example, as follows.

Step 1. Fitting Fourier series for nonlinear correlation term. By means of the MATLAB fitting tools (utilizing the third-order Fourier series (1) to fit the correlation scatter plots) [10], the mathematical model of nonlinear correlation term is straightforward and acquired as shown in Figure 8(a), together with the estimated values of Fourier parameters presented in Table 1. By substitution of estimated values together with into formula (1), time series of fitting static displacement in the whole August are acquired as shown in Figure 8(b).

Table 1: Estimated values of fitting parameters ().

Fitting parameters

Estimated values
0.1102
−0.1026
0.0218
0.01214
0.006845
−0.00618
−0.003763
0.1374

(a) The fitting Fourier series of nonlinear correlation term

(b) Time series of fitting static displacement

Figure 8: The fitting Fourier series and time series of fitting static displacement.

Step 2. Stationary test for time series of discreteness. Time series of discreteness can be acquired by minus as shown in Figure 9(a). Autocorrelation function and partial correlation function of are calculated with 50 lag phases shown in Figures 9(b) and 9(c), respectively, presenting that is slowly converging into the 95% confidence intervals as the lag phase increases, which indicates bad stationarity of for mathematical modeling. Due to this, process of first-order difference for is carried out as shown in Figure 10(a), with its and shown in Figures 10(b) and 10(c), both presenting rapid convergence into the 95% confidence intervals and verifying good stationarity for processed discreteness.

(a) Time series of original discreteness

(b) The autocorrelation function

(c) The partial correlation function

Figure 9: Time series of original discreteness with its .

(a) Time series of processed discreteness

(b) The autocorrelation function

(c) The partial correlation function

Figure 10: Time series of processed discreteness with its .

Step 3. Order determination of ARMA and EGARCH. The orders of and are relative to the convergent forms of and . That is, the and clearly present the trailing property in Figures 10(b) and 10(c), with 4th and 7th of lag phase initially converging into the 95% confidence intervals, inferring that 7 and 4 are appropriately assigned to the orders of and , respectively [11, 12]. Furthermore, the orders of and are determined utilizing the AIC and BIC delimitation criterion. In detail, the statistical values of AIC and BIC from time series of processed discreteness (Figure 10(a)) are calculated, respectively, under integer assignment of and between 0 and 8 as shown in Figure 11, presenting that the most suitable orders of and are 2 and 1, respectively, corresponding to the minimum statistical values [13, 14].

(a) The statistical values of AIC

(b) The statistical values of BIC

Figure 11: The statistical values of AIC and BIC

Step 4. Parameter estimation and residual test of ordered models. Based on the specified orders above, model parameters (, , , , , and ) are further estimated to fit time series of processed discreteness, as shown in Tables 2 and 3, respectively. For testing the fitting effectiveness of estimated parameters, residuals between processed discreteness (Figure 10(a)) and its defined models with estimated parameters are analyzed using and , as shown in Figure 12, which presents that both and are consistent in the 95% confidence intervals (except is 0) and thus verifies good availability of estimated parameters for ordered models. Moreover, the standard errors between processed discreteness (Figure 10(a)) and its models with lower orders are shown in Figure 13.

Table 2: Estimating model parameters of the ARMA(7, 4) model.

The combined models
The ARMA(7, 4) model

The fitting parameters

The values
0.00
−0.49
0.94
0.75
−0.32
−0.01
0.04
0.02
0.06
−1.26
−0.48
0.68

Table 3: Estimating model parameters of the EGARCH(2, 1) model.

The combined models
The EGARCH(2, 1) model

The fitting parameters

The values
−0.76
0.27
0.65
0.28
0.06

(a) The autocorrelation function

(b) The partial correlation function

Figure 12: The correlation functions and of residuals.

(a) The standard errors when

(b) The standard errors when

Figure 13: The standard errors between processed discreteness and its models with lower orders.

Step 5. Simulation for time series of discreteness and static displacement . Based on the mathematical models of ARMA and EGARCH with estimated parameters, time series of simulative processed discreteness are shown in Figure 14(a), and through inverse calculation of first-order difference, time series of simulative original discreteness are obtained as shown in Figure 14(b). Together with time series of fitting static displacement in Figure 8(b), time series of simulative static displacement are ultimately shown in Figure 14(c), definitely similar to time series of monitoring static displacement in Figure 5(b).

(a) Time series of simulative processed discreteness

(b) Time series of simulative original discreteness

(c) Time series of simulative static displacement

Figure 14: Simulation for time series of discreteness and static displacement .

3.2. Modeling for Dynamic Girder Displacement

3.2.1. Modeling Theory and Method

Time series of dynamic girder displacement show consistent power spectrum density (Figure 7), which can be mathematically simulated by harmonic superposition method [15–17]. Primarily, function expression of power spectrum density should be confirmed to lay foundation for harmonic superposition. Considering that straightforward function fitting will ignore local feather of acute peak at 10⁻¹ Hz around of frequency, thus the power spectrum density are decomposed into two parts: part one to fit the whole trend with ignorance of acute peak; part two to specifically fit the acute peak. Furthermore, each part can be expressed by the -order Gaussian series in logarithmic form; that is,

where denotes function expression of part one, denotes function expression of part two, and , and are fitting parameters of and . Accordingly, the function expression of power spectrum density can be formulized as follows:

Based on power spectrum density above, variance within frequency band of can be calculated, utilizing the expanding characteristics of average power spectrum for stationary and stochastic process, as follows:

Then the frequency band of [] is divided into subintervals, and the whole values in the th subinterval can be represented by the corresponding value at the middle frequency ; thus formula (11) is furthermore discretized as follows:

where

As for the th subinterval, its sinusoidal function containing middle frequency and standard deviation can be expressed as follows:

where is the random variable of uniform distribution within . With superposition of from each subinterval, the time series of dynamic girder displacement can be mathematically modeled by harmonic superposition function as follows:

3.2.2. Detailed Procedure

Based on the theory and method above, the detailed procedure for mathematically modeled time series of dynamic girder displacement is illustrated, taking the correlation scatter plots from August 1 to August 8 in Figure 15(a), for example, as follows.

(a) Original power spectrum density

(b) Processed power spectrum density

Figure 15: Original and processed power spectrum densities.

Step 1. Decreasing discreteness for original power spectrum density. Considering that certain discreteness existing in original power spectrum density can cover the acute peak at 10⁻¹ Hz, adverse to fitting the -order Gaussian series, therefore, average process with double frequency scales is carried out to decrease discreteness, specifically by 10⁻³ Hz average process within frequency band [10−1.5 Hz, 10−0.5 Hz] and 10⁻² Hz average process within other frequency bands. Process result is shown in Figure 15(b), showing distinct acute peak for fitting the -order Gaussian series.

Step 2. Fitting Gaussian series for two parts of processed power spectrum density. The processed power spectrum density can be divided into two parts: the whole trend with ignorance of acute peak as shown in Figure 16(a); the specific acute peak as shown in Figure 16(b). By means of the MATLAB fitting tools (utilizing the fourth-order Gaussian series (9) to fit two parts) [10], the mathematical models and of two parts are, respectively, shown in Figures 15(a) and 15(b), together with estimated parameters presented in Table 4. The combination of and is shown in Figure 16(c).

Table 4: Estimated parameter values of the fourth-order Gaussian series.

Fitting curves
Estimated parameter values of the fourth-order Gaussian Series

−5.36
0.45
−4.86
−0.73
−4.98
−2.98
−0.75
−0.22
1.69
0.38
1.73
0.36

0.58
0.85
0.14
0.70
−0.96
−0.96
−1.20
−1.10
0.01
0.05
0.04
0.15

(a) The first part and its fitting curve

(b) The second part and its fitting curve

(c) The combination of and

Figure 16: Two parts of processed power spectrum densities and their fitting curves.

Step 3. Harmonic superposition for fitting power spectrum density is divided into 25000 subintervals within frequency bands [10⁻⁵ Hz, 10−0.5 Hz]. In each subinterval, one middle frequency and its corresponding value are existent and then substituted into (15) for summation; thus time series of dynamic girder displacement can be acquired as shown in Figure 17(a) (simulated for 86400 s). By comparing its power spectrum density shown in Figure 17(b) with the monitoring one shown in Figure 15(a), good similarity of the whole trend and the acute peak verifies appropriateness of mathematical models for dynamic girder displacement.

(a) Time series of simulative dynamic displacement

(b) Power spectrum density of

Figure 17: Time series of simulative dynamic displacement and its power spectrum density.

4. Model Test and Evaluation

According to mathematical modeling process above, time series of lateral girder displacement are expressed by combination of third-order Fourier series , ARMA(7,4), EGARCH, and harmonic superposition function . For verifying feasibility and effectiveness of the whole mathematical models, time series from September 1 to September 3 are simulated with intervention of monitoring static wind (Figure 18(a)) and then compared with monitoring ones during same period (Figure 18(b)).

(a) Time series of simulative lateral displacement

(b) Time series of monitoring lateral displacement

Figure 18: Time series of simulative and monitoring lateral girder displacement.

According to mathematical modeling theory and procedure above, comparison is divided into two parts: (1) time series of static girder displacement; (2) time series of dynamic girder displacement. As for the first part, its simulative and monitoring results are shown in Figure 19(a) and linear fitting curve of correlation scatter plots is shown in Figure 19(b), presenting consistent variation tendency in Figure 19(a) and approximating to in Figure 19(b) (where and , respectively, denote simulative and monitoring results), which verifies good feasibility and effectiveness of mathematical models for static girder displacement. As for the second part, power spectrum densities of simulative and monitoring results are shown in Figure 20, presenting uniform variation tendency of both whole trends and acute peaks, which verifies good feasibility and effectiveness of mathematical models for dynamic girder displacement. Therefore, mathematical models above can be reasonably utilized to simulate time series of lateral girder displacement from main girder of Sutong Cable-Stayed Bridge.

(a) Simulative and monitoring results

(b) Correlation scatter plots with fitting curve

Figure 19: Simulative and monitoring results and their correlation scatter plots with fitting curve.

Figure 20: Power spectrum densities of simulative and monitoring results.

5. Conclusions

Based on monitoring data from main girder of Sutong Cable-Stayed Bridge, time series of lateral girder displacement effect are mathematically modeled by methods of fitting Fourier series and Gaussian series, combined models of ARMA(7,4) and EGARCH(2,1), and harmonic superposition function. And conclusions can be drawn as follows.

(1)Scatter plots between lateral static wind velocity and lateral static displacement present apparent nonlinear correlation, which is similar to quadratic parabolic curve, and time series of lateral dynamic displacement contain obvious stable power spectrum density with no variation by time.(2)Time series of lateral static displacement can be decomposed into nonlinear correlation term and discreteness term. Moreover, nonlinear correlation term can be mathematically modeled by third-order Fourier series with intervention of lateral static wind velocity, and discreteness term can be mathematically modeled by the combined models of ARMA and EGARCH.(3)Through decreasing discreteness in double frequency scales and division in double frequency bands, power spectrum density of lateral dynamic displacement can be mathematically modeled by the fourth-order Gaussian series, and time series of lateral dynamic displacement can be further mathematically modeled by harmonic superposition function.(4)By the comparison between simulative and monitoring lateral displacement effect from September 1 to September 3, mathematical models are feasible and effective to simulate time series of lateral girder displacement from main girder of Sutong Cable-Stayed Bridge.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the National Science and Technology Support Program (no. 2014BAG07B01), the National Natural Science Foundation (no. 51178100), and Key Program of Ministry of Transport (no. 2011318223190 and no. 2013319223120).

References

F. Xu, A. Chen, and Z. Zhang, “Aerostatic wind effects on the Sutong Bridge,” in Proceedings of the 3rd International Conference on Intelligent System Design and Engineering Applications (ISDEA 2013), pp. 247–256, January 2013.

M. Gu and H. F. Xiang, “Analysis of buffeting response and its control of Yangpu Bridge,” Journal of Tongji University, vol. 21, no. 3, pp. 307–314, 1993 (Chinese).

D. Green and W. G. Unruh, “The failure of the Tacoma bridge: a physical model,” American Journal of Physics, vol. 74, no. 8, pp. 706–716, 2006.

J. Cheng and R.-C. Xiao, “A simplified method for lateral response analysis of suspension bridges under wind loads,” Communications in Numerical Methods in Engineering, vol. 22, no. 8, pp. 861–874, 2006.

X.-H. Long, L. Li, and L. Hu, “Time domain analysis of buffeting responses of sidu river suspension bridge,” Engineering Mechanics, vol. 27, no. 1, pp. 113–117, 2010 (Chinese).

M. Yu, H.-L. Liao, M.-S. Li, C.-M. Ma, and M. Liu, “Field measurement and wind tunnel test of buffeting response of long-span bridge under skew wind,” Journal of Experiments in Fluid Mechanics, vol. 27, no. 3, pp. 51–55, 2013 (Chinese).

H. Wang, A. Li, J. Niu, Z. Zong, and J. Li, “Long-term monitoring of wind characteristics at Sutong Bridge site,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 115, pp. 39–47, 2013.

F. T. Lombardo, J. A. Main, and E. Simiu, “Automated extraction and classification of thunderstorm and non-thunderstorm wind data for extreme-value analysis,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 97, no. 3-4, pp. 120–131, 2009.

H. Wang, A.-Q. Li, C.-K. Jiao, and X.-P. Li, “Characteristics of strong winds at the Runyang Suspension Bridge based on field tests from 2005 to 2008,” Journal of Zhejiang University: Science A, vol. 11, no. 7, pp. 465–476, 2010.

MATLAB, MATLAB online manual, Natick, Mass, USA, 2010

E. Castaño and J. Martínez, “Use of the crosscorrelation function in the identification of ARMA models,” Revista Colombiana de Estadistica, vol. 31, no. 2, pp. 293–310, 2008.

T. I. Lin and H. J. Ho, “A simplified approach to inverting the autocovariance matrix of a general ARMA(p, q) process,” Statistics and Probability Letters, vol. 78, no. 1, pp. 36–41, 2008.

B. Clarke, J. Clarke, and C. W. Yu, “Statistical problem classes and their links to information theory,” Econometric Reviews, vol. 33, no. 1–4, pp. 337–371, 2014.

J. C. Escanciano, I. N. Lobato, and L. Zhu, “Automatic specification testing for vector autoregressions and multivariate nonlinear time series models,” Journal of Business & Economic Statistics, vol. 31, no. 4, pp. 426–437, 2013.

J. A. Munoz, J. R. Espinoza, and C. R. Baier, “Decoupled and modular harmonic compensation for multilevel statcoms,” IEEE Transactions on Industrial Electronics, vol. 61, no. 6, pp. 2743–2753, 2014.

S. Yan and W. Zheng, “Wind load simulation by superposition of harmonic,” Journal of Shenyang Architectural and Civil Engineering Institute, vol. 21, no. 1, pp. 1–4, 2005 (Chinese).

C.-X. Li and C.-Z. Liu, “RBF-neural-network-based harmony superposition method,” Journal of Vibration and Shock, vol. 29, no. 1, pp. 112–116, 2010.