Application of the Matrix Pencil Method to Rational Fitting of Frequency-Domain Responses

This paper presents a general methodology based on the matrix pencil method (MPM) for the fitting of frequency-domain responses in order to be properly represented in time-domain analysis. By virtue of the proposed method, a rational fitting of the frequency-domain responses can be inferred, which, in turn, helps with their inclusion into time-domain calculations. The proposed technique is well suited for the estimation of any type of function including the case of those with superimposed noise. The main feature of the method is its direct solution, hence avoiding any iteration in the estimation process. Moreover, the method does not require starting poles as opposed to the vector-fitting (VF) method. This paper presents the validation of the proposed approach by fitting the frequency responses of general a priori known functions, artificially created noisy functions and those of power components, namely: the elements of a power transformer admittance matrix and the characteristic admittance and propagation functions of a single-core sheathed cable. Finally, a time-domain analysis referring to lightning transients on a transmission line is presented.

[1]  L. De Tommasi,et al.  An identification technique for macromodeling of long interconnects , 2005, SPI 2005.

[2]  B. Gustavsen,et al.  Fast Realization of the Modal Vector Fitting Method for Rational Modeling With Accurate Representation of Small Eigenvalues , 2009, IEEE Transactions on Power Delivery.

[3]  Adam Semlyen Contributions to the Theory of Calculation of Electromagnetic Transients on Transmission Lines with Frequency Dependent Parameters , 1981, IEEE Transactions on Power Apparatus and Systems.

[4]  B. Gustavsen,et al.  Computer Code for Rational Approximation of Frequency Dependent Admittance Matrices , 2002, IEEE Power Engineering Review.

[5]  T. Hughes,et al.  Signals and systems , 2006, Genome Biology.

[6]  M. Paolone,et al.  Lightning-Induced Overvoltages Transferred Through Distribution Power Transformers , 2009, IEEE Transactions on Power Delivery.

[7]  A. Semlyen,et al.  Fast and accurate switching transient calculations on transmission lines with ground return using recursive convolutions , 1975, IEEE Transactions on Power Apparatus and Systems.

[8]  A. Semlyen,et al.  Rational approximation of frequency domain responses by vector fitting , 1999 .

[9]  Z.D. Grcev,et al.  Frequency Dependent and Transient Characteristics of Substation Grounding Systems , 1997, IEEE Power Engineering Review.

[10]  R. B. Standler,et al.  A frequency-dependent finite-difference time-domain formulation for dispersive materials , 1990 .

[11]  T. Noda,et al.  Identification of a multiphase network equivalent for electromagnetic transient calculations using partitioned frequency response , 2005, IEEE Transactions on Power Delivery.

[12]  A. S. Morched,et al.  A high frequency transformer model for the EMTP , 1993 .

[13]  Taku Noda,et al.  Phase domain modeling of frequency-dependent transmission lines by means of an ARMA model , 1996 .

[14]  Tom Dhaene,et al.  An algorithm for direct identification of passive transfer matrices with positive real fractions via convex programming , 2011 .

[16]  B. Gustavsen Fast Passivity Enforcement for Pole-Residue Models by Perturbation of Residue Matrix Eigenvalues , 2008, IEEE Transactions on Power Delivery.

[17]  B. Gustavsen,et al.  Enforcing Passivity for Admittance Matrices Approximated by Rational Functions , 2001, IEEE Power Engineering Review.

[18]  Raj Mittra,et al.  Problems and Solutions Associated with Prony's Method for Processing Transient Data , 1978, IEEE Transactions on Electromagnetic Compatibility.

[19]  Jose R. Marti Accuarte Modelling of Frequency-Dependent Transmission Lines in Electromagnetic Transient Simulations , 1982 .

[20]  Hao Shi A Closed-Form Approach to the Inverse Fourier Transform and Its Applications , 2008, IEEE Transactions on Electromagnetic Compatibility.

[21]  Tapan K. Sarkar,et al.  Matrix pencil method for estimating parameters of exponentially damped/undamped sinusoids in noise , 1990, IEEE Trans. Acoust. Speech Signal Process..

[22]  Leonid Grcev,et al.  Computer analysis of transient voltages in large grounding systems , 1996 .

[23]  A. S. Morched,et al.  A universal model for accurate calculation of electromagnetic transients on overhead lines and underground cables , 1999 .

[24]  Adam Semlyen,et al.  Vector fitting by pole relocation for the state equation approximation of nonrational transfer matrices , 2000 .

[25]  T. Sarkar,et al.  Rational modeling by pencil-of-functions method , 1981, ICASSP.

[26]  Y. Hua,et al.  Generalized pencil-of-function method for extracting poles of an EM system from its transient response , 1989 .

[27]  P. Dehkhoda,et al.  A New Method for the Inclusion of Frequency Domain Responses in Time Domain Codes , 2011 .

[28]  Tapan K. Sarkar,et al.  On SVD for estimating generalized eigenvalues of singular matrix pencil in noise , 1991, IEEE Trans. Signal Process..

[29]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[30]  B. Gustavsen,et al.  Improving the pole relocating properties of vector fitting , 2006, 2006 IEEE Power Engineering Society General Meeting.

[31]  Robert H. Halstead,et al.  Matrix Computations , 2011, Encyclopedia of Parallel Computing.

[32]  Tapan K. Sarkar,et al.  Rational modeling by pencil-of-functions method , 1981, ICASSP.

[33]  F. Rachidi,et al.  Mitigation of lightning-induced overvoltages in medium Voltage distribution lines by means of periodical grounding of shielding wires and of surge arresters: modeling and experimental validation , 2004, IEEE Transactions on Power Delivery.

[34]  Jean Mahseredjian,et al.  On a new approach for the simulation of transients in power systems , 2007 .

[35]  M. Ianoz,et al.  Transient analysis of multiconductor lines above a lossy ground , 1999 .

[36]  B. Gustavsen,et al.  Wide band modeling of power transformers , 2004, IEEE Transactions on Power Delivery.

[37]  P.T.M. Vaessen Transformer model for high frequencies , 1988 .

[38]  M. Ianoz,et al.  Influence of a lossy ground on lightning-induced voltages on overhead lines , 1996 .

[39]  A. Semlyen,et al.  A robust approach for system identification in the frequency domain , 2004, IEEE Transactions on Power Delivery.

[40]  T. Sarkar,et al.  Using the matrix pencil method to estimate the parameters of a sum of complex exponentials , 1995 .

[41]  Raj Mittra,et al.  Problems and solutions associated with Prony's method for processing transient data , 1978 .

[42]  B. Gustavsen,et al.  Calculation of electromagnetic transients in transmission cables and lines taking frequency dependent effects accurately into account , 1995 .

[43]  F. Rachidi,et al.  Analysis of Transmission Lines With Arrester Termination, Considering the Frequency-Dependence of Grounding Systems , 2009, IEEE Transactions on Electromagnetic Compatibility.