Research Article  Open Access
Jiling Ding, "The Hierarchical Iterative Identification Algorithm for MultiInputOutputError Systems with Autoregressive Noise", Complexity, vol. 2017, Article ID 5292894, 11 pages, 2017. https://doi.org/10.1155/2017/5292894
The Hierarchical Iterative Identification Algorithm for MultiInputOutputError Systems with Autoregressive Noise
Abstract
This paper considers the identification problem of multiinputoutputerror autoregressive systems. A hierarchical gradient based iterative (HGI) algorithm and a hierarchical least squares based iterative (HLSI) algorithm are presented by using the hierarchical identification principle. A gradient based iterative (GI) algorithm and a least squares based iterative (LSI) algorithm are presented for comparison. The simulation results indicate that the HLSI algorithm can obtain more accurate parameter estimates than the LSI algorithm, and the HGI algorithm converges faster than the GI algorithm.
1. Introduction
System identification studies mathematical models of dynamic systems by fitting experimental data to a suitable model structure [1, 2]. Many practical systems have multiple inputs and multiple outputs such as chemical processes [3, 4], automation devices [5–7], and network communication engineering [8–10]. For decades, much research has been performed on the multivariable systems [11, 12], and some typical approaches for the parameter estimation of the multivariable systems have been reported [13], such as the canonical approach [14], the iterative methods [15, 16], and the least squares methods [17]. Recently, Panda and Vijayaraghavan adopted the sequential relay feedback test to estimate the parameter of the linear multivariable systems [18]. Jafari et al. presented an iterative least squares algorithm to identify the multivariable nonlinear systems with colored noises [19].
The multivariable systems contain both parameter vectors and parameter matrices, and the systems inputs and system outputs are relevant and coupled [20–22]. For the sake of reducing the computational complexity, the hierarchical identification principle is utilized to transform a complex system into several subsystems and then to estimate the parameter vector of each subsystem [23, 24], respectively. In this literature, Schranz et al. proposed a feasible hierarchical identification process for identifying the viscoelastic model of respiratory mechanics [25]. Xu et al. developed the parameter estimation for dynamical response signals [26, 27].
The iterative methods have been widely applied in identifying the parameters of linear or nonlinear systems [28–30]. Many iterative algorithms for system identification are based on the gradient method [31]and the least squares method [32–35]. The basic idea of iterative methods is to update the parameter estimates using batch data.
This paper focuses on the parameter estimation for outputerror autoregressive (OEAR) systems using the hierarchical identification principle and the iterative identification principle and presents a hierarchical gradient based iterative (HGI) algorithm and a hierarchical least squares based iterative (HLSI) algorithm. The key is to decompose a multiinput OEAR system into two subsystems and then to identify each subsystem. The work in [36, 37] discussed the singleinput singleoutput systems, but many practical systems have multiple inputs and multiple outputs with the development of industrial technology. Compared with the work in [36, 37], this paper discusses the parameter estimation for multiinput OEAR systems and the presented HLSI algorithm can achieve higher estimation accuracy than the LSI algorithm, and the HGI algorithm also can achieve higher estimation accuracy than the GI algorithm.
The rest of this paper is organized as follows. Section 2 gives some definitions and the identification model of multiinput OEAR systems. Section 3 presents a gradient based iterative algorithm and a least squares based iterative algorithm for multiinput OEAR systems. Section 4 derives a hierarchical gradient based iterative algorithm. Section 5 derives a hierarchical least squares based iterative algorithm. Section 6 provides two illustrative examples to demonstrate the effectiveness of the proposed algorithms. Finally, concluding remarks are given in Section 7.
2. The Problem Formulation
Let us define some notation.
Symbols: meaning : an dimensional column vector whose entries are all 1 : a large positive constant, for example, : the transpose of the vector or matrix : : defined as : the maximum eigenvalue of the symmetric real matrix .
Consider the following multiinputoutputerror type models: where is the system output, , , are the system inputs, and is the colored noise with zero mean. and are polynomials in the unit backward shift operator , and Assume that the orders are known, , , and as . The colored noise can be fitted by a moving average process or an autoregressive process or an autoregressive moving average process where is the white noise with zero mean and , are polynomials in the unit backward shift operator : This paper considers the colored noise to be an autoregressive process, so the models in (1) can be taken as the multiinput OEAR systems.
Define the intermediate variables: From (4) and (7), we have The output in (1) can be written as Define the parameter vectors as and the information vectors as According to the above definitions, (8) and (9) can be written as Equation (14) is the identification model of the multiinput OEAR system.
3. The Gradient Based and Least Squares Based Iterative Algorithm
Consider the data from to and define quadratic criterion function as Let be an iteration variable and be the estimate of at iteration . Minimizing by using the negative gradient search, we can obtain where is an iterative stepsize. Because the information vector contains the unknown variables and , we use the estimates and at iteration to replace the unknown variables and ; we can obtain the gradient based iterative (GI) algorithm for estimating the parameter vector of the multiinput OEAR systems:The convergence rate of the GI algorithm is slow. To improve the convergence speed, we derive a least squares based iterative (LSI) identification algorithm. Minimizing and letting the derivative of with respect to be zero give the LSI identification algorithm for the multiinput OEAR systems:
4. The Hierarchical Gradient Based Iterative Algorithm
Define intermediate variables:Using the hierarchical identification principle, the multiinput OEAR system in (14) can be decomposed into two fictitious subsystems:Next, we identify the parameters and of each subsystem in (21) and (22), respectively. Define quadratic criterion functions asLet and be the estimates of and at iteration . Using the negative gradient search and minimizing and , we can obtainHere, and are the iterative stepsizes or convergence factors. Substituting (19) into (24) and (20) into (25), we can obtainThe parameter estimates and cannot be computed by (26) and (27), because the information vectors and contain unknown variables and , and the parameter vectors and in (26) and (27) are unknown. We solve this problem by replacing the unknown variables and with their corresponding estimates and at iteration and define the estimates and at iteration asFrom (12) and (14), we haveSubstituting and with their estimates and , we can get the estimates and at iteration :Replacing , in (26) and (27) with their estimates and , replacing in (26) with its estimate , and replacing in (27) with its estimate , we haveIn order to guarantee the convergence of and , a conservative choice isAt last, we can summarize the hierarchical gradient based iterative parameter estimation (HGI) algorithm for estimating and of the multiinput OEAR systems:The steps of computing the parameter estimates and for the multiinput OEAR systems are as follows.(1)Set the data length , let , and collect the inputoutput data .(2)Collect the inputoutput data and .(3)To initialize, let , , , , for .(4)Form and by (36) and (37), and form by (35).(5)Choose a satisfying (40) and update the estimate using (33) and using (34).(6)Read and using (41) and (42) and compute using (38) and using (39).(7)Give a small positive . If , increase by 1 and go to Step (); otherwise, obtain the parameters and and increase by 1 and go to Step ().
5. The Hierarchical Least Squares Based Iterative Algorithm
The HGI algorithm can produce higher parameter estimation accuracy compared with the GI algorithm, but it converges slowly. In order to solve this short board, we derive a hierarchical least squares based iterative algorithm for the multiinput OEAR systems.
Minimizing and letting the partial derivative of with respect to be zero and minimizing and letting the partial derivative of with respect to be zero, respectively, we can obtain the least squares estimate :Inserting (19) into (44) and (20) into (45) givesThe above estimates and are impossible to compute, since the righthand side of (46) contains the unknown parameter vector and the unknown information vectors and and the righthand side of (47) also contains the unknown parameter vector and the unknown information vectors and . We solve this difficulty by replacing , with their estimates , and replacing in (46) and in (47) with their estimates and . Then, we can summarize the hierarchical least squares based iterative (LSI) algorithm of estimating the parameter vectors and as follows:
The procedure for computing the parameter estimation and is as follows.(1)Give the data length , let , collect the inputoutput data , and give a small positive .(2)Collect the inputoutput data and .(3)To initialize, let , , for .(4)Form , , and by (51), (52), and (50), respectively.(5)Update the estimates and by (48) and (49) and read by (55).(6)Compute by (53) and by (54).(7)If , increase by 1 and go to Step (); otherwise, obtain the parameters and and increase by 1 and go to Step ().
6. Example
Example 1. Consider the following twoinput OEAR system:The inputs are taken as two persistent excitation signal sequences with zero mean and unit variance and as a white noise sequence with zero mean and variance .
Take the data length , applying the LSI algorithm and the HLSI algorithm to estimate the parameters of this example system. The parameter estimates and their errors of the LSI algorithm are shown in Table 1, the parameter estimates and their errors of HLSI algorithm are shown in Table 2, and the parameter estimation errors of the LSI and HLSI algorithms versus are shown in Figure 1.
From the simulation results in Tables 1 and 2 and Figure 1, we can draw the following conclusions. (i)The estimation errors given by the LSI algorithm and HLSI algorithm become smaller and smaller as iteration variable increases.(ii)Under the same noise variance, the estimation errors given by the HLSI algorithm are lower than that given by the LSI algorithm.(iii)The estimation accuracy of the HLSI algorithm is close to their true values; this indicates that the proposed algorithm can effectively identify the multiinput OEAR systems.


Example 2. Consider the following another twoinput OEAR system:The simulation conditions are the same as that of Example 1, and the noise variance . Take the data length . Applying the GI algorithm and the HGI algorithm to estimate the parameters of this example system, the simulation results are shown in Tables 3–5 and Figure 2.
From the simulation results in Tables 3–5 and Figure 2, we can draw the following conclusions. (i)Under the same noise variance and data length, the HGI algorithm has less estimation errors than the GI algorithm. This shows that the HGI estimation algorithm can obtain more accurate estimates than the GI algorithm.(ii)As the iteration variable increases, the HGI parameter estimates are very close to their true values.(iii)The proposed HGI algorithm requires more iterations than the HLSI algorithm to achieve almost same estimation accuracy.



7. Conclusions
Combining the iterative technique and the hierarchical identification principle, a HGI algorithm and a HLSI algorithm are derived for identifying the multiinput OEAR systems. Compared with the GI algorithm, the HGI algorithm can generate more accurate parameter estimates. Compared with the HGI algorithm, the HLSI algorithm has faster convergence speed. The proposed methods can be extended to discuss the parameter estimation of the multiinputoutput systems with colored noise [38–42] and timedelay systems [43, 44], such as network and signal processing [45–52].
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The author is grateful to her supervisor Professor Feng Ding at the Jiangnan University for his helpful suggestions and the main idea of this work comes from him and his book MultiInnovation Identification Theory and Methods, Beijing: Science Press, 2016. This work was supported by the Natural Science Foundation of Shandong Province (China, ZR2016FL08) and the Science Foundation of Jining University (China, 2016QNKJ01).
References
 L. Xu, “A proportional differential control method for a timedelay system using the Taylor expansion approximation,” Applied Mathematics and Computation, vol. 236, pp. 391–399, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 L. Xu, L. Chen, and W. Xiong, “Parameter estimation and controller design for dynamic systems from the step responses based on the Newton iteration,” Nonlinear Dynamics, vol. 79, no. 3, pp. 2155–2163, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 K. Rajarathinam, J. B. Gomm, D.L. Yu, and A. S. Abdelhadi, “PID controller tuning for a multivariable glass furnace process by genetic algorithm,” International Journal of Automation and Computing, vol. 13, no. 1, pp. 64–72, 2016. View at: Publisher Site  Google Scholar
 D. Zumoffen, M. Gonzalo, and M. Basualdo, “Improvements on multivariable control strategies tested on the Petlyuk distillation column,” Chemical Engineering Science, vol. 93, pp. 292–306, 2013. View at: Publisher Site  Google Scholar
 J. Na, M. N. Mahyuddin, G. Herrmann, X. Ren, and P. Barber, “Robust adaptive finitetime parameter estimation and control for robotic systems,” International Journal of Robust and Nonlinear Control, vol. 25, no. 16, pp. 3045–3071, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 H. Li, Y. Shi, and W. Yan, “On Neighbor Information Utilization in Distributed Receding Horizon Control for ConsensusSeeking,” IEEE Transactions on Cybernetics, vol. 46, no. 9, pp. 2019–2027, 2016. View at: Publisher Site  Google Scholar
 L. Feng, M. Wu, Q. Li et al., “Array Factor Forming for Image Reconstruction of OneDimensional Nonuniform Aperture Synthesis Radiometers,” IEEE Geoscience and Remote Sensing Letters, vol. 13, no. 2, pp. 237–241, 2016. View at: Publisher Site  Google Scholar
 G. Gu, S. Wan, and L. Qiu, “Networked stabilization for multiinput systems over quantized fading channels,” Automatica, vol. 61, pp. 1–8, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 A. Cristofaro and S. Pettinari, “Fault accommodation for multiinput linear sampleddata systems,” International Journal of Adaptive Control and Signal Processing, vol. 29, no. 7, pp. 835–854, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Ji and X. Liu, “Unified synchronization criteria for hybrid switchingimpulsive dynamical networks,” Circuits, Systems and Signal Processing, vol. 34, no. 5, pp. 1499–1517, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 L. Jia, X. Li, and M.S. Chiu, “Correlation analysis based MIMO neurofuzzy Hammerstein model with noises,” Journal of Process Control, vol. 41, pp. 76–91, 2016. View at: Publisher Site  Google Scholar
 F. Ding, F. Wang, L. Xu, and M. Wu, “Decomposition based least squares iterative identification algorithm for multivariate pseudolinear {ARMA} systems using the data filtering,” Journal of The Franklin Institute, vol. 354, no. 3, pp. 1321–1339, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 C. Yu and M. Verhaegen, “Blind multivariable {ARMA} subspace identification,” Automatica, vol. 66, pp. 3–14, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 G. Mercère and L. Bako, “Parameterization and identification of multivariable statespace systems: a canonical approach,” <