Computation of Gram Matrix and Its Partial Derivative Using Precise Integration Method for Linear Time-Invariant Systems

Gram matrix is an important tool in system analysis and design as it provides a description of the input-output behavior for system; its partial derivative matrix is often required in some numerical algorithms. It is essential to study computation of these matrices. Analytical methods only work in some special circumstances; for example, the system matrix is diagonal matrix or Jordan matrix. In most cases, numerical integration method is needed, but there are two problems when compute using traditional numerical integration method. One is low accuracy: as high accuracy requires extremely small integration step, it will result in large amount of computation; and another is stability and stiffness issues caused by the dependence on the property of system matrix. In order to overcome these problems, this paper proposes an efficient numerical method based on the key idea of precise integration method (PIM) for the Gram matrix and its partial derivative of linear time-invariant systems. Since matrix inverse operation is not required in this method, it can be used with high precision no matter the system is normal or singular. The specific calculation algorithm and block diagram are also given. Finally, numerical examples are given to demonstrate the correctness and validity of this method.

[1]  Liubin Yan,et al.  A numerical method of calculating first and second derivatives of dynamic response based on Gauss precise time step integration method , 2010 .

[2]  Francis T.K. Au,et al.  Assessment and improvement of precise time step integration method , 2006 .

[3]  Noël Tanguy,et al.  Efficient Gram-matrix computation for irrational resonant systems using Kautz models , 2007, Signal Process..

[4]  Francis T.K. Au,et al.  On the precise integration methods based on Padé approximations , 2009 .

[5]  A.V. Garcia,et al.  On the Use of Gram Matrix in Observability Analysis , 2008, IEEE Transactions on Power Systems.

[6]  Ariovaldo V. Garcia,et al.  Active and Reactive Observability Restoration Subproblems in Generalized Power System State Estimation , 2008 .

[7]  P. Agathoklis,et al.  The computation of Gram matrix via the solution of Lyapunov equation , 1991, [1991] Proceedings of the 30th IEEE Conference on Decision and Control.

[8]  V. Sreeram,et al.  The Computation of Gram Matrix via Impulse-Response Grammians , 1991, 1991 American Control Conference.

[9]  Ian Postlethwaite,et al.  Multivariable Feedback Control: Analysis and Design , 1996 .

[10]  G. Korres.  A Gram Matrix-Based Method for Observability Restoration , 2011, IEEE Transactions on Power Systems.

[11]  A.V. Garcia,et al.  A Numerical Method for Finding Spanning Trees in Power System State Estimation , 2006, 2006 International Conference on Power System Technology.

[12]  Ariovaldo V. Garcia,et al.  Identifying critical sets in state estimation using Gram matrix , 2009, 2009 IEEE Bucharest PowerTech.

[13]  Nikolaos M. Manousakis,et al.  State estimation and observability analysis for phasor measurement unit measured systems , 2012 .

[14]  W. Zhong,et al.  A Precise Time Step Integration Method , 1994 .

[15]  Frederic Ward Williams,et al.  Accurate high-speed computation of non-stationary random structural response , 1997 .

[16]  Noël Tanguy,et al.  Orthogonal decomposition of derivatives and antiderivatives for easy evaluation of extended Gram matrix , 2006, Signal Process..

[17]  Ariovaldo V. Garcia,et al.  Power system observability analysis based on gram matrix and minimum norm solution , 2009, 2009 IEEE Power & Energy Society General Meeting.

[18]  W. Zhong,et al.  On precise integration method , 2004 .

[19]  Diane L. Peters.  Coupling and Controllability in Optimal Design and Control. , 2010 .

[20]  V. Sreeram,et al.  Evaluation of Gram matrix diagonal elements using inners technique , 1991 .

[21]  V. K. Jain,et al.  Identification of linear systems through a Grammian technique , 1970 .

[22]  Weidong Li,et al.  An Energy Conservation Algorithm for Nonlinear Dynamic Equation , 2012, J. Appl. Math..

[23]  G. Korres.  A numerical algorithm for observability analysis based on Gram matrix , 2010 .

[24]  Frederic Ward Williams,et al.  A high precision direct integration scheme for structures subjected to transient dynamic loading , 1995 .

[25]  Victor Sreeram,et al.  On the computation of the Gram matrix in time domain and its application , 1993, IEEE Trans. Autom. Control..