Estimation of regions of attraction of power systems by using sum of squares programming

This paper presents an algorithm estimating the regions of attraction of power systems based on the Lyapunov function approach where a sublevel set of a Lyapunov function for a target system is used as the estimate. In particular, we focus here on the algorithm based on sum of squares (SOS) programming, which has been recently proposed, and aim to develop a simpler algorithm for the practical use. For this aim, we present an algorithm overcoming the difficulty of the SOS programming problem addressed in the existing study, i.e., the bilinear constraints, in a simpler way. In the proposed algorithm, two SOS programming problems are iteratively solved, and the number of the problems solved at each iteration is reduced to half of that in the existing algorithm. In addition, we theoretically analyze the proposed algorithm, and show the convergence under certain conditions. The performance of our algorithm is demonstrated by numerical examples.