Optimal control of stochastic system with Fractional Brownian Motion

In this paper, we introduce a class of stochastic harvesting population system with Fractional Brownian Motion (FBM), which is still unclear when the stochastic noise has the character of memorability. Stochastic optimal control problems with FBM can not be studied using classical methods, because FBM is neither a Markov pocess nor a semi-martingale. When the external environment impact on the system of FBM, the necessary and sufficient conditions for the optimization are offered through the stochastic maximum principle, Hamilton function and Itô formula in our work. To illustrate our study, we provide an example to demonstrate the obtained theoretical results, which is the expansion of certainty population system.


Introduction
In the biological population system, it is affected by a variety of external factors, which they are likely to change the population's amount. In order to control the development of biological population reasonably, it is necessary to select appropriate control variables and establish reasonable performance indicators to study the optimal control of stochastic population systems. In the paper, a nonlinear population system equation (the harvesting equation) is discussed. The typical harvesting system can be described in the form: ∂r + ∂p(r,t) ∂t = −λ(r, t, P(t))p(r, t) − u(r, t)p(r, t), p(r, 0) = p 0 , p(0, t) = A 0 β(r, t, P(t))p(r, t)dr, P(t) = A 0 p(r, t)dt.
(1. 1) where p(r, t) is the density of the population of age r at time t, A is the life expectancy, and p(A, t) = 0. λ 1 is the average mortality ratio of the population of age r at time t, β 1 is the average fertility ratio of the population of age r at time t, u 1 (r, t) is harvesting effort function, which is the control variable in the model and satisfies: 0 ≤ u(t) ≤ u max . f 1 (r, t, P(t)) + g 1 (r, t, P(t)) dB H dt is the stochastic perturbation, effecting of external environment on the population system, such as earthquakes, emigration, impacts of extra terrestrial objects, and so on.
The stochastic model has aroused concern in the recent years. Abel Cadenikls [1] used a stochastic maximum principle for systems with jumps, with applications to finance systems. Zhang [2][3][4][5][6] investigated the stability of numerical solutions for the stochastic age-dependent system. Zhang [7,8] provided the different methods for the numerical solutions in the stochastic system with age-dependent. Pei [10] focuses on asymptotic mean-square boundedness of several numerical methods applied to a class of stochastic age-dependent population equations with Poisson jumps. Emel Savku et al. [11,12] illustrated that they contributed to modern OR by hybrid (continuous-discrete) dynamics of stochastic differential equations with jumps and the optimal control. However, compared with stochastic system driven by the classical Brownian motion, FBM is a family of centered Gaussian random process indexed by the Hurst parameter H ∈ (0, 1) with continuous sample paths. Some special kinds of dynamical systems require both Wiener process and FBM to model their dynamics. Meanwhile, few has been done because classical methods to solve stochastic problems can not be used directly, since Fractional Brownian Motion (FBM) is not a semi-martingale and not a Markov process. Ma [13] developed a numerical scheme and show the convergence of the numerical approximation solution to the analytic solution for stochastic age-dependent population equations with FBM. Kloeden [14] used the multilevel Monte Carlo method introduced by Giles [15] to stochastic differential equations with Fractional Brownian Motion of Hurst parameter H > 1/2 and achieved a prescribed root mean square error of order ε with a computational effort of order ε −2 . Duncan [16] discussed the solutions semi-linear stochastic systems with FBM. Zhou [17] investigated the stability for the delayed neural networks with FBM.
On the other hand, optimal control problems have also attracted wide attention, due to their several applications in population system, economic system, finance system [11,12,[18][19][20]. Luo [18] studied optimal harvesting control problem for an age-dependent competing system of n-dimension competing species. Zhao [19] and Chen [20] talked about optimal control of different stochastic system. He [21] investigated optimal harvesting problem for age-structured species. However, an optimal control problem requires the minimization of a criterion function of the states and control inputs of the system over a set of admissible control functions [22], which creates huge troublesome. Stochastic optimal control problem driven by FBM is the bottleneck problem. In this paper, all the previous fields are combined to consider the optimal control problem of stochastic harvesting population system with FBM, the necessary and sufficient conditions for the optimization are obtained, and the example for the obtained theoretical results is illustrated. We provide below a brief summary of our results.
• We introduce the fractional Brownian noise into a class of stochastic harvesting population system and establish necessary as well as sufficient conditions of optimal control, which has not been studied before; • Using the stochastic maximum principle, Hamilton function and Itô formula to stochastic harvesting equations with Fractional Brownian Motion and study the optimal control of the system; • The example is presented, and it supports our theoretical results. The paper is divided into five sections. The assumption, notations and some basic definition are given in section 2. In section 3, we establish necessary as well as sufficient conditions of optimal control. In Section 4, an example is provided to illustrate the theoretical results. The conclusions are given in section 5.

Preliminaries of the problem
Let f (x) is continuous functions, its H order fractional derivative is defined as: Where D is differential operator, H is Hurst parameter.
In this paper, we discuss stochastic optimal control problems driven by fractional Brownian motion (fBm), and consider the following stochastic control harvesting population system with FBM interval [0,A]: Getting motivation from the above facts, we discuss the optimal control problems in the system (2.1) is Were E(·) is expectation operator.
As the standing hypotheses, we always assume that the following conditions are satisfied: 1, 2), and there exists a constant c d > 0, such that (A 5 ) All x, y ∈ R n , and there exists a constant L > 0, such that all t ∈ [0, T ] satisfied (A 6 ) Let f (t, y, u), g(t, y, u), h(t, y, u) are linear functions, we introduce f (t, y, u) = D t y + E t u + F t , g(t, y, u) = G t y + H t u + I t .
The Hamiltonian function is given by H(t, q, γ, y, u) = −uy u + q, f − λy + βy − uy + γ, g(t, y, u) . Now, we introduce the adjoint equation for our problem. The adjoint equation can be written as: Note that the couple(q, γ) is the adjoint process corresponding to the stochastic system p(r, t). The adjoint equation admits one and only one F t − adapted solution (q, γ),where L(t, y, u) = −uy u . Moreover, to ensure that the above stochastic differential equation make sense, we shall consider only those F t − predictable control processes u : u ∈ U ad that satisfy p{ T 0 |F t u t |dt < ∞, T 0 |I t u t |dt < ∞} = 1.

The main results
This is the main result of this paper, in this section, we derive necessary conditions for a control to be optimal.
Lemma 3.1. J is Gâteaux-differentiable with differential given by To obtain Equation (3.1), we use Itô formula [24], Gronwall's inequality [25], and equivalently the formula of integration by parts. (3.1) The above equation may be rewritten as Where we denote for every u ∈ U; t ∈ [0, T ]: We note that H(t, u) is convex. Theorem 3.1. If case 1 hold, then a necessary condition for a control u * to be optimal for Problem (2.2) is that for every u * ∈ U: On the other hand, if case 2 holds, then inequality (3.2) is a sufficient condition of optimality for a control u * .
Proof: Here, we apply previous knowledge and methods to obatain the results, such as Young inequality [2], Itô integral (Lemma 2.1) and the Hölder, Burkholder-Davis-Gundy (BDG) inequalities [22]. According to (2.3), u * is an optimal control if and only if u * ∈ U: In case 1, we see that for every u ∈ U: Thus, in case 1 and in conjunction with (3.2), a necessary condition for a control u * to be optimal is Which is equivalent to (3.1).
On the other hand, in case 2, for every ∀u ∈ U: Thus, in case 2, a su1cient condition for a control u * to be optimal is that (3.6), or equivalently (3.2), holds for every u * ∈ U.

An example
For convenience, we still adopt the notation introduced in Section 2.
If u * (r, t) is optimal,the necessary condition for (4.1) is Moreover, because of the solution for state equation p(r, t) is the function of u(r, t), λ(r, t), f (r, t, p) and g(r, t, p) satisfying the assumption, we can conclude that

Conclusions
Existence and optimal control results of the stochastic model with Fractional Brownian Motion (FBM) is studied in this paper. Firstly, we introduce the fractional Brownian noise into a class of stochastic harvesting population system and establish necessary as well as sufficient conditions of optimal control, which has not been studied before. Secondly, Using the stochastic maximum principle, Hamilton function and Itô formula to stochastic harvesting equations with Fractional Brownian Motion and study the optimal control of the system. finally, the obtained theoretical results are verified by an illustrative example. As further direction, researchers are invited to investigate the optimal control problem for stochastic model by including Gâteaux-differentiable with differential.