Computation of the periodic steady-state response of nonlinear networks by extrapolation methods
Stig Skelboe · IEEE Transactions on Circuits and Systems · 1980
The problem of computing the periodic steady-state response can be formulated as solving a nonlinear equation of the formz = F(z)whereF(z)Is the solution vector for the nonlinear network after one period of integration from the initial vectorz. The convergence of the sequencey_0 , y_1 , \cdotsgenerated byY_{r+1} = F(y_r)can be accelerated by extrapolation methods. This paper presents a unified analysis of three extrapolation methods: the scalar and vector\epsilon-algorithms and the minimum polynomial extrapolation algorithm. The main result of the paper is the theorem giving conditions for quadratic convergence of the extrapolation methods. To obtain this result the methods are studied for linear problems (where F is a linear function) and the error propagation properties are investigated. For autonomous systems a function calledGsimilar toFcan be defined. In order to obtain quadratic convergence from the extrapolation methods, the derivatives of F and G must be Lipschitz continuous. The appendixes give sufficient conditions for the Lipschitz continuity. A discussion of practical problems related to the implementation of the extrapolation methods is based on the convergence theorem and the error analysis. The performance of the extrapolation methods is demonstrated and compared with other methods for steady-state analysis by four examples, two autonomous and two nonautonomous. Extrapolation methods are very easy to implement, and they are efficient for the steady-state analysis of nonlinear circuits with few reactive elements giving rise to slowly decaying transients.