Numerical approximation of \((n-1)\)-dimensional stable manifolds in large systems such as the power system (Q1376280)

The authors consider a nonlinear system whose dynamics is described by the equation \(\dot x = f(x)\) which has an equilibrium point at the origin \(x=0\). It is supposed that the linearized matrix at the equilibrium has one eigenvalue with a positive real part and all other eigenvalues have negative real parts. From the authors' conclusions: ``This paper proposes a simple method for computing quadratic approximations of the stable manifolds on the stability boundary. Geometric algorithms are developed for testing whether an initial condition that is ``near'' a saddle, in fact, lies within the region of attraction or outside. The algorithm is valid under certain Morse-Smale like assumptions and is easy to implement in large-system models. Power system examples used for illustrating the algorithms provide promising results on stability monitoring type analysis. The basic algorithm presented here can be readily extended for computing the stable manifolds of unstable periodic orbits and can also be modified to accommodate differential-algebraic models''.