Journal of Function Spaces

Volume 2018, Article ID 8095728, 8 pages

https://doi.org/10.1155/2018/8095728

## Averaged Control for Fractional ODEs and Fractional Diffusion Equations

^{1}Faculty of Mathematics, University of Montenegro, George Washington st. bb, 81000 Podgorica, Montenegro^{2}Faculty of Science, Department of Physics, Division of Theoretical Physics, University of Zagreb, Bijenicka cesta 30, 10000 Zagreb, Croatia^{3}Faculty of Electrical Engineering, University of Sarajevo, Zmaja od Bosne bb, 71000 Sarajevo, Bosnia and Herzegovina

Correspondence should be addressed to Darko Mitrovic; em.ca@mokrad

Received 27 May 2018; Accepted 11 July 2018; Published 17 July 2018

Academic Editor: Gisele Mophou

Copyright © 2018 Darko Mitrovic et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We generalize results concerning averaged controllability on fractional type equations: system of fractional ODEs and the fractional diffusion equation. The proofs are accomplished by introducing appropriate Banach space in which we prove observability inequalities.

#### 1. Introduction and Notation

General problem of control theory amounts to choosing certain control which would lead a system (e.g., governed by a system of (partial) differential equations) from the given initial state to the prescribed final state. Such types of problems have obviously great potential in the sense of various applications and, thus, the control theory is very well developed. It is hard to say where the first mathematical treatment of the problem essentially originates but one can find lots of information in the standard books, e.g., [1–3].

However, in certain situation it is not possible to precisely determine the coefficients governing the process and it is natural that the coefficients depend on another (essentially stochastic) variable. In such a situation, we cannot require exact controllability of the system but so-called average controllability introduced first in [4] (to the best of our knowledge).

After that, several interesting publications appeared, e.g., [5–8] and the Ph.D. thesis [9], which will be in the basis of our approach. In the current contribution, we shall extend mentioned results concerning the averaged control (in particular, those from [4, 9]) on equations containing fractional derivatives.

In the first part of the paper, we shall deal with the following finite dimensional system of ODEs:where is the fractional derivative of order (to be defined later), is the state of the system, and are given matrices depending on the (stochastic) variable , and is the given initial state. The control problem here is to find the vector-valued function such that, for the given final state (it is independent of ), we have for the final time . Remark that here we have the Lebesgue integration with respect to . However, there are no obstacles to choose integration with respect to some probability measure as we shall do in the next part (see [6, 9]). In order to get the result, we shall adapt the method from [4]. More precisely, we shall rewrite the controllability problem into the observability one and then prove the observability estimate.

We remark that we could also consider the standard initial and final data (without the -derivative), but it is more standard and convenient for the fractional ODEs to consider the type of initial data given in (1) (see, e.g., [10]).

In the second part of the paper, we shall consider the fractional diffusion equation which has received considerable attention recently (see, e.g., [11] and references therein). In our case, it has the form (see [9]) for an appropriately chosen Hilbert space . Here, is known positive bounded function (or stochastic variable defined on a probability space ), is a bounded operator for every , the function is given, and is the fractional Laplacean (to be defined later). As before, we are looking for the control function such that, for the given (final) state , it holds that for the final time . Unlike the situation from [9] where the initial-boundary value problem for the diffusion equation is considered, here we deal with the Cauchy problem which, by means of the Fourier transform, provides us with observability estimates under appropriate assumptions on the relation between and .

The paper is organized as follows. After the Introduction, we recall the notions and notations that we are going to use. In Section 3, we deal with problem (1)-(2). In Section 4, we deal with (3)-(4).

#### 2. Notions and Notations

The fractional derivative can be defined in many ways (see, e.g., a classical book [12]) but they are all essentially equivalent. However, some definitions provide easier operations depending on a situation. Here, we shall use quite rudimentary way of defining the fractional derivative—via the notion of the Fourier multiplier operator. Let us first introduce the notion of the Fourier transform.

*Definition 1. *For a function we define the Fourier transform by

Using the Fourier transform, we can then define the Fourier multiplier operator through the following definition.

*Definition 2. *A (Fourier) multiplier operator associated with a bounded function (see, e.g., [13]) is a mapping defined bywhere is the Fourier transform while (or ) is the inverse Fourier transform.

If, for a given , the multiplier operator satisfies where is a positive constant, then its symbol is called an (Fourier) multiplier.

Finally, we can introduce the notion of the fractional derivative.

*Definition 3. *The fractional derivative , is defined as the Fourier multiplier operator with symbol : The fractional Laplacean is defined as the Fourier multiplier operator with the symbol , : Such a definition is of course not optimal since might not be unique. However, this does not affect the essence of our results and one can always fix one of the roots generated by .

Let us informally derive the integration by parts rule for the fractional derivatives. Formal derivation requires unsubstantial adaptations. We remark that we extend the involved functions by zero in domains where the functions are not defined. We have by means of the Plancherel theorem for smooth enough functions and vanishing at infinity (due to simplicity in notations, we present the rule in one-dimensional situation): where, as before, is the Fourier multiplier operator with the symbol .

Now, we shall introduce the notion of controllability and observability and provide the theorems connecting the two notions. We will take the notations from [4, 9].

*Definition 4. *We say that system (1) is controllable if for every there exists an -component control vector function such that for the prescribed final time .

Interestingly, the latter definition is essentially equivalent to the following one taken from [9] and adapted to system (3) (see, e.g., [2, 3] for details).

*Definition 5. *Equation (3) is controllable in average with the cost if for all there exists such that where is the Sobolev space of functions with the fractional derivative of order and Using a classical functional analytic tools such as Hahn-Banach theorem, one can also prove that controllability is equivalent to the observability. In order to introduce it, we need the notion of the adjoint problem to (3). Let be a solution to Indeed, using the Plancherel theorem, it is not difficult to check that since is actually the multiplier operator with the even symbol and thus it is an operator mapping real functions into real functions (see [14, Lemma 5]). Thus, (14) is the corresponding adjoint equation to (3). We have the following definition of observability.

*Definition 6. *We say that system (3) is observable in average if there exists such that for all it holds thatwhere is the solution to (3) and is the adjoint operator to .

In [9] one can find proof that the averaged controllability (Definition 5) is equivalent to the averaged observability (Definition 6) in a much more general situation where is replaced by a general operator defined on a general Hilbert space . Also, is taken to be from a Hilbert space while is a bounded functional from to .

The following theorem holds.

Theorem 7 (see [9]). *Problem (3) is controllable in average if and only if it is observable in average.*

#### 3. Fractional ODE

Although classical ODEs have been used for so many years as a tool for modeling of the systems with the applications in physics, biology, chemistry, medicine, and engineering, recent experiments indicate that there is a large class of complex systems with different kinetics which have a microscopic complex behavior, and their macroscopic dynamics cannot be described by classical models. We mention several contributions in different directions which is just a microsnapshot of the activities in this direction. Some examples can be found in the [15] where it has shown that fractional calculus models of viscoelastic materials are consistent with the physical principles that govern such materials. Concerning the signal processing and sampling theory in [16] authors derive an interpolation functions necessary for direct evaluation of fractional derivatives directly from the sample values. In mathematical modeling in biology several papers regarding the fractional-order differential models of biological systems with memory, such as dynamics of tumor-immune system and dynamics of HIV infection of CD4 + T cells and fractional-order predator-pray models, have been considered in [17–19]. Additional applications concern the spread of contaminants in underground water, network traffic, charge transport in amorphous semiconductors, cell diffusion process, the transmission of signals through strong magnetic fields such as those found within confined plasma, etc. A number of surveys with collections of applications can be found in [12].

In the case of system (1), the vector-valued function is the state of the system, is a -matrix governing its dynamics through the fractional ODE-system, and is a -component control vector in acting on the system through the control matrix depending on the parameter .

The matrices and are measurable and bounded with respect to . We remark that the initial state as well as the final state are required to be independent of , while it is obviously not a case with the solution of (1).

As usual in the control-type problems, we shall need the adjoint system in order to use the observability concept. The adjoint system in the case of system (1) has the form where is the adjoint matrix to the matrix .

Indeed, we have for an arbitrary function using (10) (we remind the reader again that we extend the functions by zero on the sets where they are not defined, i.e., out of the interval ): Applying the operator here, we reach (17). Following [9, Definition ] and [4, p.4 ], we have the following averaged observability inequality for the adjoint system (17): for all . We have the following theorem.

Theorem 8. *System (1) fulfills the averaged controllability property from Definition 5 if and only if the adjoint system (17) satisfies (19) and both are equivalent to the rank conditions**When these properties hold, the averaged control of minimal -norm is given by where is the solution of the adjoint system (17) corresponding to the datum minimizing the functional*

*Proof. *We have for a solution to (17) (see the derivation (18) above) according to the controllability condition from Definition 5. In other words, Definition 5 is recasted as which is the Euler-Lagrange equation associated with the minimization of the functional from (22) (see [4, p.4. ]. Now, the proof follows the steps from the proof of [4, Theorem 3].

##### 3.1. Numerical Example

At this time, we want to give a simple numerical illustration of the theory that we have presented so far. The most difficult issue here is to determine the final value (at ) for the adjoint problem (17) since it is equal to the datum minimizing the functional (22). This means that we essentially need to know explicit dependence of and from and then we can apply the Euler-Lagrange equations techniques. Here, we shall simply assume that we have a minimizing final value for the adjoint state, compute the control function using (21), and then see what final value we can get. In the further research, we shall develop a numerical method in the frame of which we are able to solve the original control problem (finding control from the equations, initial and final states).

Here, we aim to solvewhere and are as as follows: As for the adjoint system (17), we shall assume that . We then devise the following strategy.

Let matrices and corresponding to system (25) of order and initial condition be given.(1)Formulate the adjoint problem (17).(2)Solve (as will be described later) the adjoint problem and obtain the solution .(3)Use formula (21) to compute the averaged control .

So the only challenge now is to solve step (2), that is, to solve fractional-order ODE. To this end we will use Grunwald-Letnikov discretization; for further reference on modern numerical treatment of the fractional-order ODE please see [20].

*Definition 9. *Let . The operatorwhere and , is called the Grunwald-Letnikov fractional derivative of order .

This definition holds for arbitrary functions , but the convergence of the infinite sum cannot be ensured for all functions [21]. After employing Grunwald-Letnikov discretization in the case of system we have After solving with respect to and we obtain the finite difference scheme for computing step (2). where the coefficients are given by the following formulas: By assuming that we see that we reach the final state . The result is plotted on Figure 1.