Optimal Control and Bifurcation Issues for Lorenz-Rössler Model

Optimal control is one of the most popular decision-making tools recently in many researches and in many areas. The Lorenz-Rossler model is one of the interesting models because of the idea of consolidation of the two models: Lorenz and ossler. This paper discusses the Lorenz-Rossler model from the bifurcation phenomena and the optimal control problem (OCP). The bifurcation property at the system equilibrium is studied and it is found that saddle-node and Hopf bifurcations can be holed under some conditions on the parameters. Also, the problem of the optimal control of Lorenz-Rossler model is discussed and it uses the Pontryagin’s Maximum Principle (PMP) to derive the optimal control inputs that achieve the optimal trajectory. Numerical examples and solutions for bifurcation cases and the optimal controlled system are carried out and shown graphically to show the effectiveness of the used procedure.


Introduction
Prediction of any system's development is a very important goal, especially in the case of chaotic systems which exist frequently in several real-life and various fields.
These systems are very important for the service of mankind. Those systems include psychology [1], secure communications [2], economics system modeling [3], medicine [4], etc. Due to the importance of these models, these deserved to be studied. Despite the recent trend of forecasting, this task is often not easy.
The chaotic models usually need to study the optimization or the OCP. The Principle of Optimality means to take the best choice (procedure) to minimize the cost (maximize the profit) of the current time (stage) and all subsequent times (stages). The OCPs in general does not have a perfect solution, so solutions are often approximate, which is another aspect of difficulty and sensitivity.
Lorenz system is a reduced version of a larger system studied earlier by Barry Saltzman [5]. This model is a system of three non-linear ordinary differential equations, is extremely sensitive to perturbation of the initial conditions, and this system has strange attractors; thus this system has a chaotic behavior for different values of parameters and different initial conditions as it is in the famous butterfly attractor which is produced for the special values of parameters of (10, 28, 8/3) respectively as these parameters shown sequentially in the model in [6]. Rössler attractor is another famous attractor which is one of the products of the works of a German biochemist Otto Eberhard Rössler. There is a similarity between Lorenz and Rössler attractors, where the latter follows an outward spiral around two fixed points. Although each variable in the system is oscillating around specified values, the oscillations are chaotic [7]. There is a merging of the two models of Lorenz and Rössler in a non-additive form, and because the new additive system losses the chaotic behavior property, some switching between variables was done and some transactions were manipulated [6] [8]. Although the chaotic systems are difficult to predict in a long-term, this paper discusses the optimal control problem of the Lorenz-Rössler through some external inputs. It is important to mention that the Lorenz-Rössler which should be studied is that which is presented in [6] and [8]. We use some procedures in the optimal control problem to determine some sets of numbers of parameters and initial conditions that achieve the optimal behavior to the goal state.
The qualitative changes in the trajectories in the phase space due to the change in one or more control parameters are called bifurcations. The bifurcating study is possible for a one-dimensional system with one parameter. But it is difficult in higher-dimensional cases, especially with several parameters. so there is little research in this area [9].
In the next section, we provide a necessary mathematical introduction. The mathematical system of the Lorenz-Rössler model is presented and a brief discussion of the stability of the system is in Section 3. In Section 4, an analytical investigation of some cases of bifurcation is discussed and many diagrams for those cases are presented. In Section 5, the optimal control problem is discussed followed by many digital examples that were made through simulations. The conclusion is presented in Section 6.

Mathematical Introduction
It is known that the optimal control problem requires: A mathematical form for the system to be controlled, description of the constraints, determination of the goal to be accomplished, usually it is an additional boundary condition and de-Open Journal of Optimization termination for the performance measure [10] [11] [12]. Firstly, we consider the simplest example to find the shortest length between two specified locations Consider the case of no constraints on the variables, and for simplicity, we assume that the minimal curve is given as the graph of a smooth function ( ), The problem now can be illustrated as minimizing the integral which is given by where J is called a functional or the objective function, that must be minimized with respect to t. In a simple case of one dependent variable ( ) x t and no constraints, and conditionally that the optimal curve ( ) x t exist and unique, the general problem now is to find the optimal curve that minimizes the functional where g is a continuous function in all its variables and has a continuous first and second orders partial derivatives with respect to all its variables. Moreover, 0 t and f t are fixed. If we consider a small variation in the carve, then where  is a small parameter and ( ) . See Figure 1.
Thus, the functional in (2.2) can be written as The necessary condition in order to be extremum function is [13] Under the assumption that x and all its derivatives are continuous, with some mathematical processes and considering that ( ) ( ) Equation (  Let us now consider the n first ordinary differential equations (ODEs) as the following forms , and suppose that the optimal state exist and unique, that is the vector of n twice differentiable function. In this case, the functional (2.2) can take the following form Under the same conditions in the case of one dependent variable, the function ( ) g x which makes the integral (2.9) an extremum must satisfy the n simultaneous E-L equations, which are given by with the boundary conditions See [11] and [13] for more details. Now, in the case when there are constraints on the states and control variables. The functional that need to maximize (minimize) takes the following form ( ) where 0 0 λ ≥ and one can get 0 1 λ = for maximization [12]. Theorem: Assume is the optimal function that maximizes the objec- is the corresponding trajectory, then * u must satisfy the following conditions [11] [14] ( ) ( ) This system consists of 2n nonlinear differential equations with n initial con- For more details about this theorem and its proof see [11] and [14]. Notable: An additional equation is required if f t is indeterminate.

Lorenz-Rössler Mathematical
The Lorenz-Rössler system is a three-dimensional system with five parameters.
This system is described by the following equations as presented in [6] and [8] ( ) ( ) where 1 2 , x x and 3 x are the state variables of the system, 1 2 3 1 2 , , , , a a a b b and 3 b are the system parameters. Clearly, the zero-state is not a solution of the system (3.1) because the system is not homogeneous, and with a few mathematical calculations we can be sure that this system has the following possible equilibrium states ( ) , , It is easy to show that the system (3.1) under some conditions, is unstable at least at one of its steady-states, so be it 1 E . The Jacobian matrix W of the model (3.1) is given by x a According the linear stability analysis and theory of linear differential equations, we strive to find the eigenvalues of 1 W . The determinant equation of 1 W is given by the following equation:

Bifurcation of Lorenz-Rössler System
In this section, we discuss the bifurcation phenomenon of the considered system. At the first equilibrium point, 1 E depends upon the following characteristic Equation (3.8), that can be rewritten as the following:

Optimal Control Problem in a Period of Time
In the case of constraints on the control variables, Pontryagin maximum principleis considered as a design tool to get the best possible trajectory for a dynamical system by providing a necessary condition that must hold for an optimum, but not (in general) sufficient conditions [11]. So, the basic problem is to find the controllers that maximize/minimize the functional J in (2.12), which is called a co-state function and considered as a type of L-m. We use the PMP to find the best possible controllers with respect to a choice of specified measure. Open Journal of Optimization Subject to: The controlled system of (3.1) that is given by  i e are the controlling inputs that be determined by the PMP with respect to the optimality measure for the system (3.1) near its stead-states.
The Hamiltonian equations are given by: with the following boundary conditions:

Numerical Simulation
In the following, some numerical solutions of the system in Equations (5.12), -The optimal control to the stationary state Figure 5.
-The optimal control to the stationary state ( ) is shown in Figure 6 where  ( ) ( )  Figure 5 and Figure 6 indicate that, according to the assumed values of the parameters, the optimal controlled state 1 2 3 , , x x x converge with time to the assumed goal level 0, 0, 5, respectively in Figure 5 and to 0.450, 0.675, 0.325 respectively in Figure 6. Also, in each case the co-state variables 1 2 3 , , γ γ γ disappear with time. The assumed goal levels are represented by the dotted lines. All results indicate the possibility of the optimal control of the Lorenz-Rössler system, and the PMP has shown excellent results in achieving the optimal behavior of the system.

Conclusion
Many studies can be implemented on the Lorenz-Rössler model, but in this paper, we have focused on the issues of the bifurcations and the optimal control problem of the system. The bifurcation analysis of the system at the equilibrium state ( ) ciple is considered to solve the optimal control problem. The optimal control inputs were analytically derived and it is found that they are functions of the co-state variables which disappear when the system arrives at the ideal state.
Analytical methods are used to solve the necessary conditions, while the non-linear differential equations of the optimal controlled system are solved numerically by the math software Maple, and then some illustrative solutions are shown graphically.