International Journal of Engineering Research and Technology. ISSN 0974-3154, Volume 14, Number 6 (2021), pp. 502-509 © International Research Publication House. http://www.irphouse.com Stability Analysis of the Lorenz System using Hurwitz Polynomials Fabián Toledo Sánchez1, Pedro Pablo Cárdenas Alzate2 and José Rodrigo González Granada3 1,2,3 Department of Mathematics and GEDNOL, Universidad Tecnológica de Pereira, Pereira, Colombia. Abstract Among these fundamental studies is the stability analysis; this stability analysis can be performed using the results of In the qualitative study of differential equations, we can Lypaunov, who presented two methods to determine stability analyze whether or not small variations or perturbations in the that solve particular problems according to the structure of the initial conditions produce small changes in the future, this differential equations or using Hurwitz type polynomials. intuitive idea was formalized by Lyapunov in his work "The General Problem of Motion Stability". The stability analysis of If the nonlinear ordinary differential equation is expressed in a system of differential equations can be performed using the the form 푥̇ = Ax, then the problem of analyzing the stability methods proposed by Lyapunov or using criteria to obtain becomes a problem of algebraic type, since it is only enough to Hurwitz type polynomials, which provide conditions to analyze know the roots of the characteristic polynomial associated to the dynamics of the system by studying the location of the roots the matrix 퐴 which correspond to the eigenvalues and to of the characteristic polynomial associated to the system. In this observe if these have negative real part; if the above occurs it paper we present a stability analysis of the Lorenz system using is said that the system 푥̇ = Ax is asymptotically stable. If this stability criteria to obtain Hurwitz type polynomials. polynomial has the above mentioned characteristic it is said to be a Hurwitz polynomial [6], [7]. Keywords: Stability of differential equations, Hurwitz polynomials, Hurwitz criterion, Lorenz system. Therefore, a relevant fact to analyze the stability of the system 푥̇ = Ax corresponds to establish criteria to determine when the characteristic polynomial associated to the matrix 퐴 is Hurwitz, I. INTRODUCTION i.e., if all its roots have negative real part. In [8] a series of criteria are presented to obtain Hurwitz type polynomials, these The stability analysis of systems of ordinary differential criteria present some equivalences in their formulation, among equations is the object of study in various mathematical models them the Routh-Hurwitz criterion, the Lienard-Chipart applied to a varied branch of the exact sciences and control conditions, the Hermite-Biehler theorem, the stability test and theory, where the aim is to provide necessary and sufficient the Routh algorithm [9]. conditions so that in a given problem with established initial conditions, solutions close to this initial value remain close Thus, the objective of this paper is to analyze the stability of throughout the future or in a given case tend to the equilibrium the Lorenz system using the Routh-Hurwitz stability criteria solution. This notion of closeness and proximity for small and the Routh algorithm. For this purpose, four sections are variations in the initial conditions was initially worked on by presented: in the first one, the Routh-Hurwitz stability criteria outstanding mathematicians such as Lagrange and Dirichlet, and the Routh algorithm are deduced; in the second section, a but it was not until 1892 when the Russian mathematician and stability analysis of the Lorenz system is performed using the physicist Aleksandr Lyapunov laid the foundations for this criteria presented to obtain Hurwitz type polynomials and concept of proximity in his doctoral thesis entitled "The finally, the results are analyzed by means of simulations General Problem of Motion Stability" [1]. performed in MATLAB using the ode45 function for the solution of the system of nonlinear differential equations. The The study of the dynamics of these sensitive variations in the simulations validate the results obtained on the stability of the initial conditions is the fundamental work in the qualitative Lorenz system. analysis of ordinary differential equations. The importance of such sensitive variation is observed in the Lorenz system, this nonlinear model that described the movement of air masses in the atmosphere gave way to a branch of mathematics called II. HURWITZ POLYNOMIALS chaos theory, because Lorenz, studying the weather patterns, Let began to observe unusual behaviors to those that had been established that should happen. Lorenz using a numerical 풙̇ = 퐴풙 (1) simulation used 0.506 as an initial data as a less accurate be a system of linear or nonlinear ordinary differential approximation of the data 0.506127, obtaining as a result, equations, where 퐴 is a square matrix and 풙 is a vector. The according to two tabulations, two completely different climate stability of the system (1) at its equilibrium point can be scenarios [2], [3]. This small variation or perturbation in their determined by performing a study of the eigenvalues of the initial values yielded completely different climatological associated matrix 퐴. This algebraic study establishes that if the results. From the moment Lorenz formulated the problem with roots of the characteristic polynomial associated to the matrix this question, various characteristics of the system began to be A have negative real part, then we conclude that the system (1) studied, such as regions of stability, attractor basins, is asymptotically stable. bifurcations, chaos, geometric aspects, among others [4], [5]. 502 International Journal of Engineering Research and Technology. ISSN 0974-3154, Volume 14, Number 6 (2021), pp. 502-509 © International Research Publication House. http://www.irphouse.com Now, the problem of determining the stability in a system of This matrix is constructed as follows: differential equations from the study of the roots of the In the first row are the coefficients of the polynomial characteristic polynomial of 퐴 associated to the linear or (3) with odd location starting with 푎1. nonlinear system, translates into the task of finding necessary In the second row are the coefficients of the and sufficient conditions for which all the roots of this polynomial (3) with even location starting with 푎0. polynomial are located in the left half of the complex plane, i.e., The elements of each subsequent row are formed so if they have negative real part. that the component ℎ푖푗 is given by: This way of obtaining polynomials of this type with negative real part was initially proposed by the physicist James Maxwell 푎 푠푖 0 < 2푗 − 푖 ≤ 0 in 1868, who presented a solution to this problem for ℎ = { 2푗−푖 푖푗 polynomials of degree 3, then in 1877, the Canadian 0 푖푛 표푡ℎ푒푟 푐푎푠푒 mathematician Edward Routh presented an algorithm that As a result of the construction, the coefficients 푎1, 푎2, 푎3, solved the problem in a more general way providing explicit … , 푎푛 are on the main diagonal of the matrix, and all the conditions for polynomials up to degree 5 and later it was the elements of the last column are null, except the last element German mathematician Adolf Hurwitz in 1895 who presented which is 푎푛. The matrix ℋ is called Hurwitz Matrix. a properly analytical solution to the problem [6]. Routh-Hurwitz theorem. The polynomial (3), with its Definition 1: A polynomial with real coefficients is said to be positive leading coefficient (푎0 > 0), is a Hurwitz polynomial Hurwitz if all its roots have a negative real part, that is, if all its if and only if all the diagonal principal minors of the Hurwitz roots lie in ℂ−, the left half-plane of the complex plane, matrix ℋ are positive [10]. ℂ− = { 푎 + 푏푖 ∶ 푎 < 0 } (2) The principal diagonal minors of the matrix (4) are given by the following determinants, 푎1 푎3 푎5 II.I Criteria for obtaining Hurwitz polynomials 푎1 푎3 Δ1 = ∣ 푎1 ∣, Δ2 = | |, Δ3 = |푎0 푎2 푎4|, 푎0 푎2 Let 0 푎1 푎3 푛 푛−1 2 푎1 푎3 푎5 푎7 푃(휆) = 푎0휆 + 푎1휆 + ⋯ + 푎푛−2휆 + 푎푛−1휆 + 푎푛 (3) 푎0 푎2 푎4 푎6 Δ4 = | |, ⋯ , Δ푛 = 푎푛 ∙ Δ푛−1. be the characteristic polynomial associated to the matrix 퐴 of 0 푎1 푎3 푎5 the system (1). 0 푎0 푎2 푎4 In the search for algorithms to determine the location of the roots of the polynomial (3) the following question was asked: II.I.II Routh criterion how to determine if the polynomial (3) is Hurwitz? For the coefficients of the polynomial (3) the Routh This question was studied by several mathematicians and arrangement is constructed as follows, physicists, among the first who proposed to solve this type of problem was the Austrian engineer A. Stodola who at the end 푎0 푎2 푎4 푎6 ⋯ of the 19th century was interested in the problem of finding 푎1 푎3 푎5 푎7 ⋯ conditions under which all the roots of a polynomial had a negative real part; but in 1895 it was Hurwitz who presented a 푏0 푏1 푏2 푏3 ⋯ (5) 푐 푐 ⋯ solution to the previous question based on the work of Hermite. 푐0 푐1 2 3 푑 푑 ⋯ 푑0 푑1 2 3 There are several criteria to obtain Hurwitz type polynomials. ⋮ ⋮ ⋮ ⋮ ⋱ In this paper we only briefly present the Routh-Hurwitz criterion and the Routh algorithm. In the first row are the coefficients of the polynomial (3) with even location starting with 푎0. II.I.I Routh-Hurwitz criterion In the second row are the coefficients of the polynomial To present the Routh-Hurwitz criterion we first construct the (3) with odd location starting with 푎1. following matrix from the coefficients 푎0, 푎1, … , 푎푛 of the The elements of each subsequent row are formed polynomial (3): according to the following algorithm 푎 푎 푎1 푎3 푎5 푎7 ⋯ 0 0 2 푏0 = 푎2 − 푎3 , 푏1 = 푎4 − 푎5 , ⋯ 푎 푎 푎 푎 ⋯ 0 푎 푎 0 2 4 6 1 3 푎 푎 푎 ℋ = 0 1 2 3 ⋯ 0 (4) 푎1 푎3 푎2 푎 푎 ⋯ 0 0 4 6 푐0 = 푎3 − 푏1 , 푐1 = 푎5 − 푏2 , ⋯ ⋮ ⋮ ⋮ ⋮ ⋱ 0 푏0 푏1 [ 0 0 0 0 ⋯ 푎푛] 푏0 푏1 푑0 = 푏1 − 푐1 , 푑1 = 푏2 − 푐2 , ⋯ 푐0 푐1 ⋮ ⋮ 503 International Journal of Engineering Research and Technology.