A New Optimal Homotopy Asymptotic Method for Fractional Optimal Control Problems

Solving fractional optimal control problems (FOCPs) with an approximate analytical method has been widely studied by many authors, but to guarantee the convergence of the series solution has been a challenge. We solved this by integrating the Galerkin method of optimization technique into the whole region of the governing equations for accurate optimal values of controlconvergence parameters Csj. -e arbitrary-order derivative is in the conformable fractional derivative sense. We use Euler–Lagrange equation form of necessary optimality conditions for FOCPs, and the arising fractional differential equations (FDEs) are solved by optimal homotopy asymptotic method (OHAM). -e OHAM technique speedily provides the convergent approximate analytical solution as the arbitrary order derivative approaches 1.-e convergence of the method is discussed, and its effectiveness is verified by some illustrative test examples.


Introduction
e general definition of an optimal control problem necessitates minimization of an objective function of the states and control inputs of the system over a set of admissible control functions.
e optimal control problems emerge naturally in different areas of science, engineering, and economics. Several significant works have been done in the area of optimal control of integer order derivative dynamic systems: the optimal control design for the hyperchaotic Chen system [1], an optimal control of time-varying and state time-delay systems using shooting method [2], a CasADi software for nonlinear optimal control problem [3], direct and simultaneous approach for constrained mixedinteger optimal control [4], a tailored algorithm for mixed integer nonlinear optimal control model of hybrid electric vehicles [5], and many more. e arbitrary order of fractional derivative equations (FDEs) from fractional calculus has contributed a substantial part to mathematics, engineering, and different area of sciences. e previous studies review that arbitrary order derivatives provide more exact models of several engineering control systems and show the behaviour, performance, and mathematical representation of the dynamic systems in sciences than traditional calculus, for example, the positioning of satellite in space is a model using fractional order because traditional calculus is better equipped to capture the memory-like effect observed in the system. Fractional calculus applications are but not limited to liquid-containing gas bubbles [6], neurodynamics system [7], the FOCPs with a general derivative operator [8], a new adaptive synchronization and hyperchaos control of a biological snap oscillator [9], a new mathematical model for Zika virus transmission [10], the fractional features of a harmonic oscillator with position-dependent mass [11], electrohydrodynamic flow in a cylindrical conduit [12], fractional order of HIV infection model [13], and novel solution methods for IBVPs of arbitrary-order with conformable differentiation [14]. Many areas of FDEs have been researched by many authors: the multistage ADM for solving NLP problems over a nonlinear fractional dynamical system [15], the existence of the solution to FOCPs [16], stable analytical techniques for solving FDEs [17], Caputo-Fabrizio operator sense for brucellosis model [18], and spectral collocation method of generalized fractional operator sense for fractional calculus equations [19]. Currently, to obtain exact solutions for FOCPs is still difficult and remain a challenge. But where there is not an exact solution, appropriate approximate analytical solution (AAS) is suitable and applicable. A recent approximate analytical methods are: modified ADM for fractional optimal control problems (FOCPs) [20], optimal control of a constrained fractionally damped elastic beam [21], analyses of an optimal solutions of optimization problems from fractional gradient based system using VIM [22], conformable fractional optimal control problem of an heat conduction equations using Laplace and finite Fourier sine transforms [23], spectral Galerkin approximation [24], and direct transcription methods for FOCPs [25], but the aforementioned methods lack convergence norm that will guarantee the convergence of the series solution. In 1992, Liao proposed the homotopy analysis method (HAM) [26] for solving the nonlinear problem which was later advanced to optimal homotopy asymptotic method for noninteger order [27], new fractional homotopy method for optimal control problems (OCPs) [28], and comparisons of OHAM [29]. But OHAM has never been used to solve FOCPs, which drives this research work. We extend the work of (Alizadeh Effati, 2018a), who proposed a fractional optimal control problem with an approximate analysis solution that lacks convergence norm. We provided the solution using OHAM with an optimization technique of the Galerkin method. It is hoped that the simplicity of this formulation will initiate a new interest in the area of fractional order optimal control problems. Our focus in this paper is to widen the application of the OHAM to FOCPs. e formulation, numerical scheme, and approximate analytical results for FOCPs presented in this paper are attempts to fill this gap.
We arrange the paper as follows. In Section 2, a brief introduction to the fractional calculus and the necessary optimality conditions are discussed. In Section 3, we formulate OHAM with FOCPs and show the convergence analysis of the technique. We present numerical examples and results in Section 4. Finally, in Section 5, we present the conclusion.

Fractional Calculus
e definition of conformable fractional derivative (CFD) preserves many properties of classical order derivatives [30]. Some features that we will adopt are as follows.
In the right case, we can similarly prove the following.

Necessary Optimality Condition.
We define optimal control problems as a minimization of a performance index (PI) subject to system dynamic constraints on the state variable z(t) and control variable u(t). e necessary optimality condition helps reduce OCPs to a system of integer order ODE. e general formulation of Caputo and Riemann-Liouville derivatives is with OCPs [31]. Using the above definition equations (1) and (2), we present an Euler-Lagrange equation form for FOCPs under some consideration as follows.
Find the optimal control u(t) that minimizes the performance index: subject to the system dynamic constraints and the initial condition where F and G are two arbitrary functions and T a α is the conformable fractional derivative [30]. Some additional terms containing z ′ s(t) may be included in equation (6) at the endpoint. When α � 1, the problem above reduces to a integer OCPs [32]. Here, we take the limits of integration from 0 to 1 and consider 0.8 ≤ α ≤ 1. According to [33], the necessary equations are obtained through the combination of equations (5) and (6) and using a Lagrange multiplier approach. e variants of the resultant is taken, which is followed by using integration by parts to revise the equation so that it does not contain variants of a derivative term, enforce the necessary terminal conditions, and set the coefficients of δλ, δz, and δu to zero. We have Where λ is the Lagrange multiplier. e details of the derivation of equations (8)-(11) are seen in [33]. e Euler-Lagrange equations for FOCP equations (6) and (7) represent the necessary conditions for the optimality of FOCP. e difference between arbitrary and integer OCP is that the arbitrary order FOCP consists of the left and the right FD. Examine that equation (8) contains left conformable derivative, whereas equation (9) contains the right conformable derivative. is demonstrates that the solution of OCP requires the insight of not only forward derivatives (FD) but also backward derivatives (BD) to account for end conditions. In integer-order OCP theories, this issue is either not addressed or not highlighted. is is broad because the BD of order 1 turns out to be negative of the FD of order 1. It can be shown that, in the lim α⟶1 , equations (8)-(10) goes back to those obtained using standard methods.

The Formulation of OHAM with FOCPs
We exemplify the basic concept of the OHAM by considering the following general FDEs system: with initial conditions where T a α is the conformable fractional derivative, l k is a linear operator, N k is a nonlinear operator, t is an independent variable, x k (t) is unknown function, φ is the problem domain, and g k (t) is a known function. According to OHAM, one can construct a homotopy map H k (ϕ i (t, p): φ × [0, 1] ⟶ φ which satisfies that equation (12) can be constructed using OHAM as [26] ( where embedding parameter (ℓ) is 0 ≤ ℓ ≤ 1, auxiliary function H k (ℓ) ∀ℓ ≠ 0, unknown function (5 k (t, ℓ)), and H(0) � 0. Definitely, when ℓ � 0 and ℓ � 1, it holds that 5 k (t, 0) � ψ k,0 (t) and 5 k (t, 1) � ψ k (t), respectively. us, as ℓ moves from 0 to 1, the solution 5 k (t, ℓ) approaches from ψ k,0 (t) to ψ k (t), where initial guess ψ k,0 (t) satisfies the linear operator generated from equation (14) for ℓ � 0 as H k (ℓ) is given as where C s j can be known later. We get approximate solution by expanding 5 k (t, ℓ, C j ) in Taylor's series in terms of ℓ using equation (14) with collections of the coefficient like power of ℓ, gives the governing equations ψ i,0 (t) in a linear form in equation (15). en, the 1 st problem is given as where N k,m (ψ 0 (t), ψ k,1 (t), . . . , ψ k,m (t)) is the coefficient of ℓ m , produced by expanding N k (5 k (t, ℓ, C j )) in series relating to ℓ e convergence of series solution equation (21) relies on C s j . If it is convergent at ℓ � 1, it gives solution to equation (12) as Using equation (22) in equation (12), we have expression for the residual error as If International Journal of Differential Equations then ψ k (t, C j ) is the exact solution. Usually, such a case does not occur. We adopt Galerkin method to find the optimal values C s j as follows: Minimize the functional where the values of a and b depend on the given problem.
With these known C s j , the approximate analytical solution equation (12) is well known.
We calculate the correctness of OHAM by (2) Error norm L ∞

Convergence Analysis of the Method
is governed by equation (22) under the definitions equations (19) and (20), it achieves the solution of equations (12) and (13).
If the C j is chosen properly, then equation (36) leads to the solution of equations (12) and (13)

Example 1. Find the control u(t) that minimizes the quadratic performance index
subject to system dynamics and the initial condition e exact solution for the system for α � 1 is given as where Following equations (8)-(11) procedure, we have Using OHAM, the solutions for fractional-order are acquired. We choose the linear and nonlinear operators in the following forms: Using homotopy in equation (14), where Substituting 5 1 (t, ℓ), 5 2 (t, ℓ), and H k (ℓ) into equations (47) and (48), and equating the coefficient of likes power of ℓ, it gives linear FDEs as follows: Using the operator Lemmas 1 and 2 on the above equations with initial condition gives Summing up the solution from equations (58)-(63), the 3 rd -order approximate analytical method generated by OHAM, for ℓ � 1, is We determine C 1 and C 2 by using the procedure mentioned in equations (24)- (26).
(69) Figure 1 shows the approximate analytical solutions of the state variable z(t) and control variable u(t) at different values of α using OHAM after the 3 rd order of approximate solution. Table 1 shows the control-convergence parameter at different values of fractional-order for both state and control variables, Table 2 shows comparisons and absolute errors between the proposed method and the exact solution at integer order α � 1, and Table 3 shows the optimal values of the cost function and the comparisons between OHAM and VIM. It is observed that OHAM converges to optimal values faster in CPU times and iterations.
(91) Figure 2 shows the approximate analytical solutions of the states variables z 1 (t) and z 2 (t) and control variable u(t) at different values of α using OHAM after the 3 rd order of approximate solution. Tables 4 and 5 show the controlconvergence parameter at different values of fractional-order for both states and control variables, Table 6 shows .577999908 x 10 − 10 0.3913043497 0.9 − 6.661323216 x 10 − 10 0.3612134421 0.8 − 6.755673906 x 10 − 10 0.3125152370 Table 6: Comparisons of (OHAM and exact, α � 1).  Table 7 shows absolute errors between the proposed method and the exact solution at integer order at α � 1, and Table 8 shows the optimal values of the cost function and the comparisons between OHAM and VIM. It is observed that OHAM converges to optimal values faster in CPU times and iterations.

Conclusion
In this research article, we examine the effectiveness of the OHAM to generate an accurate approximate analytical solution to FOCPs. We derive the necessary optimal condition and then solve the system of nonlinear FDEs with OHAM. e technique helps provide an accurate approximate analytical solution to the problems. e fractional derivative operator is in the conformable fractional derivative operator sense. We show that Galerkin method helps to determine the optimal values of the control-convergence parameter. It can be seen that the OHAM results corresponded with the exact solution at α � 1, which shows the effectiveness of the proposed method to provide an accurate approximate analytical solution. OHAM is used to solve this type of problems for the first time and gave optimal values to the convergence norm challenge for accurate series solution. No existing approximate solution is known to have considered such a nonlinear constrained programming problem with OHAM ahead of this one. OHAM technique is very easy to implement and efficient for solving FOCPs.

Data Availability
No data were used to support this study.  International Journal of Differential Equations 9