^{1}

^{*}

^{1}

^{*}

^{1}

^{*}

In this paper, Homotopy perturbation method is used to find the approximate solution of the optimal control of linear systems. In this method the initial approximations are freely chosen, and a Homotopy is constructed with an embedding parameter , which is considered as a “small parameter”. Some examples are given in order to find the approximate solution and verify the efficiency of the proposed method.

Optimal control problems arise in a wide variety of disciplines. optimal control theory has also been used with great success in areas as diverse as economics to biomedicine [

In this paper we solve the optimal control problem by combine perturbation method. To this end, there are quite a few fundamentally diverse approaches, some of which can be found in [18,19]. The homotopy method is a powerful numerical method for solving nonlinear algebraic and functional equations. The main advantage over classical methods is that the method enjoys global convergence. However, it is not used as widely as these, mainly owing to being poorly covered in the Russian literature.

The Belgian mathematician Lahaye was the first to use the homotopy method for the numerical solution of equations. He considered the case of a single equation. He used discrete continuation by the Newton method. Later, Lahaye [

Consider the following optimal control problem

where, , and, are the time invariant given matrices. The control function is an admissible control if it is piecewise continues in t for each t in the given interval. It is assumed the control is bounded, that is, a closed, bounded, subset of exists, such that the control function takes its values form. The input can be derived by minimizing the quadratic performance index, where and are symmetric positive semi-definite and is symmetric positive definite. By using Pontryaging’s maximum principle, the optimal control law, can be achieved for system (1.1) (see [

Non-linear techniques for solving linear and non-linear problems have been dominated by the perturbation methods, which have found wide applications in engineering. But, like other non-linear analytical techniques, perturbation methods have their own limitations, Firstly, almost all perturbation methods are based on small parameters so that the approximate solutions can be expressed in a series of small parameters. This so called small parameter assumption greatly restricts applications of perturbation techniques, as is well known, an hefty gigantic of linear and non-linear problems have no small parameters at all. Secondly, the determination of small parameters seems to be a special art requiring special techniques. An appropriate choice of small parameters leads to ideal results, however, an unsuitable choice of small parameters results in bad effects. In 1997, Liu [

One may consider the following nonlinear differential equation (see [31-36])

with natural boundary conditions or tangentiality conditions as:

where A is a general differential operator, B is a boundary operator, is a known analytic function and is the boundary of the domain.

The operator can, generally, be divided into two parts L and N, where L is Linear, while N is nonlinear, so that (2.1) may written as:

By homotopy perturbation technique, we construct a homotopy which satisfies

or

where is an embedding parameter, and is an initial approximate solution of Equation (2.1).

Obviously from Equation (2.5)

By changing continuously from zero to unity the Equations (2.6) and (2.7) show that will change from to. In topology, this changing is called deformation, and are called homotopy functions.

In this method, using the homotopy parameter, we assume that the solution of Equation (2.5) is a power series of:

Letting results in the approximate solution of Equation (2.1) as:

Series (2.9) is convergent for most cases, the convergent rate depends upon the nonlinear operator A(v).

In this section, we apply the homotopy perturbation method to solve the optimal control system (1.1).

Consider Hamiltonian of the control system (1.1) as:

where is known as the costate variable. By Pontryagin’s maximum principle, the optimal control must satisfy the following equation:

where is a solution of the adjoint equation

with the terminal condition

Thus, from Equation (3.2), the optimal control law is

From control system (1.1) and adjoint Equation (3.3) one have:

Implementing the optimal control as a closed loop if the solution to the adjoint Equation (3.3) is assumed like Equation (3.4) as a linear function of the states in the form( see [

By using Equations (3.3), (3.6) and (3.7), we have

where the first equality follows from Equation (3.7) and the second one from Equation (3.6). Hence

Since the above equation must hold for all nonzero, must satisfy the following matrix Riccati equation

Considering Equations (3.5) and (3.7), we can see that the optimal control law is given as

and can be computed using the following relation

where, and

with conditions, and

In this section, we present some examples to show the reliability and efficiency of the method described in the previous section. In the following examples, we assume.

Example 4.1. Consider a single-input scalar system as follows (see [

According to system (1.1), we have and by using (3.12), we have

thus

So

by using (2.8), let and so,

from equating the terms with identical power of p,

where are considered as initial approximations, and imposing boundary condition, so

by using (4.3), we have

From (2.9), we have

Figures 1 and 2 show the approximated value of and, respectively.

Example 4.2. Consider a single-input scalar system as follows (see [

According to system (1.1), we have and, by using (3.12), we have

thus

So

by using (2.8), we have

from equating the terms with identical power of p,

where are considered as initial approximations, Setting and imposing boundary condition, so

by using (4.3), we have

From (2.9), we have

.

Figures 3 and 4 show the approximated value of x(t) and u(t), respectively.

Example 4.3. Consider a single-input scalar system as follows:

According to system (1.1), we have

and and by using (3.12), we have

thus

So

by using (2.8), we have

from equating the terms with identical power of p,

where, are considered as initial approximations, Setting and imposing boundary condition, the approximate and exact value for is the

exact value for is and the approximate value for and obtained from this algorithm in this following form

In this paper, we solve the optimal control problems us-

ing Homotopy perturbation method. Embedding parameter can be taken into account as a perturbation parameter. Full advantage of the traditional perturbation techniques can be taken by the novel method. The initial approximation can be freely chosen with unknown constants, which can be identified via various methods [

At last, Homotopy perturbation method is applicable method which calculates the approximate solution of linear and nonlinear problems, particularly optimal control problems.