Based on a neural field network model with impulsive and random disturbances, a preview control method that makes full use of known future information is proposed to reduce the static error of the target signal and the transient oscillatory behavior of the controlled system when it receives random disturbance inputs. The preview controller for epileptic seizures is constructed, and the feasibility and effectiveness of clinical single-target and multi-target stimulation in epilepsy regulation are explored from a computational perspective. In addition, a performance index function is proposed to evaluate the energy consumption of controller with and without preview under different input (target) strategies. Suggestions for different strategies are given in terms of the individualized disease environment of patients. From the perspective of seizure control effectiveness and performance consumption, the results show that the preview controller has a greater advantage. The theory of preview control is applied to the control of epileptic seizures for the first time, and the conclusions of the multifaceted study provide some references for clinical trials and controller applications.
Citation: Denggui Fan, Yingxin Wang, Jiang Wu, Songan Hou, Qingyun Wang. The preview control of a corticothalamic model with disturbance[J]. Electronic Research Archive, 2024, 32(2): 812-835. doi: 10.3934/era.2024039
[1] | Wanshun Zhao, Kelin Li, Yanchao Shi . Exponential synchronization of neural networks with mixed delays under impulsive control. Electronic Research Archive, 2024, 32(9): 5287-5305. doi: 10.3934/era.2024244 |
[2] | Wei Ji . Optimal control problems with time inconsistency. Electronic Research Archive, 2023, 31(1): 492-508. doi: 10.3934/era.2023024 |
[3] | Zhaoyan Meng, Shuting Lyu, Mengqing Zhang, Xining Li, Qimin Zhang . Sufficient and necessary conditions of near-optimal controls for a stochastic listeriosis model with spatial diffusion. Electronic Research Archive, 2024, 32(5): 3059-3091. doi: 10.3934/era.2024140 |
[4] | Yang Song, Beiyan Yang, Jimin Wang . Stability analysis and security control of nonlinear singular semi-Markov jump systems. Electronic Research Archive, 2025, 33(1): 1-25. doi: 10.3934/era.2025001 |
[5] | Daliang Zhao, Yansheng Liu . Controllability of nonlinear fractional evolution systems in Banach spaces: A survey. Electronic Research Archive, 2021, 29(5): 3551-3580. doi: 10.3934/era.2021083 |
[6] | Changling Xu, Huilai Li . Two-grid methods of finite element approximation for parabolic integro-differential optimal control problems. Electronic Research Archive, 2023, 31(8): 4818-4842. doi: 10.3934/era.2023247 |
[7] | N. U. Ahmed, Saroj Biswas . Optimal strategy for removal of greenhouse gas in the atmosphere to avert global climate crisis. Electronic Research Archive, 2023, 31(12): 7452-7472. doi: 10.3934/era.2023376 |
[8] | Chengbo Yi, Jiayi Cai, Rui Guo . Synchronization of a class of nonlinear multiple neural networks with delays via a dynamic event-triggered impulsive control strategy. Electronic Research Archive, 2024, 32(7): 4581-4603. doi: 10.3934/era.2024208 |
[9] | Haijun Wang, Gege Kang, Ruifang Zhang . On optimality conditions and duality for multiobjective fractional optimization problem with vanishing constraints. Electronic Research Archive, 2024, 32(8): 5109-5126. doi: 10.3934/era.2024235 |
[10] | Xinzheng Xu, Yanyan Ding, Zhenhu Lv, Zhongnian Li, Renke Sun . Optimized pointwise convolution operation by Ghost blocks. Electronic Research Archive, 2023, 31(6): 3187-3199. doi: 10.3934/era.2023161 |
Based on a neural field network model with impulsive and random disturbances, a preview control method that makes full use of known future information is proposed to reduce the static error of the target signal and the transient oscillatory behavior of the controlled system when it receives random disturbance inputs. The preview controller for epileptic seizures is constructed, and the feasibility and effectiveness of clinical single-target and multi-target stimulation in epilepsy regulation are explored from a computational perspective. In addition, a performance index function is proposed to evaluate the energy consumption of controller with and without preview under different input (target) strategies. Suggestions for different strategies are given in terms of the individualized disease environment of patients. From the perspective of seizure control effectiveness and performance consumption, the results show that the preview controller has a greater advantage. The theory of preview control is applied to the control of epileptic seizures for the first time, and the conclusions of the multifaceted study provide some references for clinical trials and controller applications.
Impulsive differential equations are employed for the analysis of real-world phenomena that are characterized by the state of instantaneous system changes. It has been found to play a vital role in many areas, such as sampled-data control, communication networks, industrial robots, biology, and so on. Due to the widespread existence of impulsive perturbation, extensive research has been dedicated to the exploration of the stability of impulsive differential systems (see [1,2,3,4]). Meanwhile, some scholars have focused on studying the optimal impulsive control problem, which has resulted in numerous interesting findings (see [5,6,7,8,9]).
Finding the necessary optimality conditions is one of the central tasks in optimal control theory. Pontryagin and his co-authors have made milestone contributions (see [10]). As pointed out in [11]: "The mathematical significance of the maximum principle lies in that maximizing the Hamiltonian is much easier than the original control problem that is infinite-dimensional". Essentially, the Hamiltonian maximization occurs pointwise. Since Pontryagin's maximum principle was discovered, the first-order and second-order necessary conditions for optimal control problems in both finite- and infinite-dimensional spaces have been extensively researched (see [11,12]).
However, it is not always possible to find the optimal control by pointwise maximizing the Hamiltonian, and this kind of optimal control problem is often referred to as the singular control problems. For singular control problems, the first task is to discover new necessary conditions to distinguish optimal singular controls from other singular controls; one way to do this is to look for second-order conditions. It is common to seek a second-order necessary condition that requires a quadratic functional to be non-negative. Ideally, we would prefer the second-order necessary conditions to have pointwise characteristics that are similar to Pontryagin's maximum principle, which means that the pointwise mode is preferable for maximizing a particular function. The pointwise necessary optimality conditions are reviewed in [13]; the Jacobson conditions and the Goh conditions are generally considered as two types of pointwise second-order necessary optimality conditions for the optimal singular control problem. In addition, references [14,15,16] are very comprehensive references on the singular control problem. Regarding the Goh conditions, the original contribution can be seen in [17]. In short, the Goh conditions are obtained by applying Goh's transformation. Goh's transformation approaches have been designed to transform the original singular problem into a new nonsingular one; then, for the nonsingular optimal control problem, the classical Legendre-Clebsch condition may be applied. In this process, Goh's transformation may be used several times; therefore, the Goh conditions are also called the generalized Legendre-Clebsch conditions; recent research results are very abundant, for example, references [18,19,20,21] and related references. Conversely, limited attention has been given to investigating the Jacobson-type conditions.
There is an interesting story about Jacobson-type necessary conditions. Until Jacobson discovered the Jacobson-type necessary conditions, it was thought that there was no Riccati-type matrix differential equation for singular control problems. In fact, Jacobson introduced a linear matrix Riccati differential equation that is similar to the nonlinear matrix Riccati differential equation in the standard LQ problem. Thus, a "new" necessary condition for singular control was obtained in [22]. After this, Jacobson found that this new condition is different from the generalized Legendre-Clebsch condition, and Jacobson has demonstrated that these two necessary conditions are generally insufficient for optimality in [23]. Therefore, the matrix Riccati differential equation can not only solve nonsingular problems, they can also solve important optimal singular control problems; singular control and nonsingular control can even be considered under a unified framework (see Theorem 4.2 in [15,24]).
Recent studies [25,26] have established the Jacobson-type pointwise second-order necessary optimality conditions for deterministic and stochastic optimal singular control problems, respectively. Particularly, the relaxation methods for control problems have been used to solve the singular control problem of a lack of a linear structure in control problems; also, the pointwise second-order necessary conditions have been obtained in [25]. When discussing the second-order necessary conditions, only the allowable set of singular control is considered, rather than the original allowable control set, which can greatly reduce the computational expense of singular control problems. References [25,26] have also been applied to derive the pointwise second-order necessary optimality conditions for singular control problems with constraints in finite- or infinite-dimensional spaces. For further details, please refer to [13,27,28,29,30,31]. Theorem 4.3 in [25] describes the Jacobson-type second-order necessary optimality condition of singular control problems governed by pulseless controlled systems according to Pontryagin. For the definition of singular control in the classical sense and in the Pontryagin sense, see Definitions 1 and 2 in [14], as well as (4.4) and (4.5) in [25].
Regarding the importance of the impulsive system, the second-order necessary conditions for the optimal control problem governed by impulsive differential systems have been studied in [32], which focuses on impulsive systems with multi-point nonlocal and integral boundary conditions; the second-order necessary conditions of the integral form has been obtained by introducing the impulsive matrix function directly. When discussing the second-order necessary conditions for singular control, the perturbation control takes its value in full space, and not in the singular control region (see (3.13)); therefore, the conclusion is not reached by the pointwise Jacobson-type necessary conditions. In addition, references [33,34,35] consider second-order necessary conditions whereby the measure is impulsive control; the system can experience an infinite number of jumps in a finite amount of time, but this is different from what we are going to consider.
Inspired by the above discussions, we aimed to study the following optimal control problem, as governed by impulsive differential systems, which differs from the previous works.
Problem P: Let T>0 and Λ={ti∣0<t1<t2<⋯<tk<T}⊂(0,T) be given; U⊆Rm is a nonempty bounded convex open set.
minJ(u(⋅))=∫T0l(t,y(t),u(t))dt+G(y(T)), | (1.1) |
subject to
{˙y(t)=f(t,y(t),u(t)),t∈[0,T]∖Λ,y(ti+)−y(ti)=Ji(y(ti)),ti∈Λ,y(0)=y0, | (1.2) |
where u(⋅)∈Uad={u(⋅)∣u(⋅)measurable,u(t)∈U}, and l, G, f, and Ji(i=1,2,⋯,k) are the given maps.
Here, we generalize the control model of [15,25] to an impulsive controlled system. For Problem P, the difficulty lies in the introduction of the appropriate impulsive adjoint matrix differential equation, which we overcome by borrowing the method presented in [25]. The main conclusions include Theorems 3.1 and 4.6, where the former is a generalization of Theorem 4.2 in [15] to impulsive controlled systems; the latter is a similar conclusion of Theorem 4.3 in [25], but the difference is that Theorem 4.6 is a conclusion of the impulsive controlled systems, and in the classical sense, while Theorem 4.3 in [25] is a conclusion for the pulseless controlled systems, and in the Pontryagin sense.
The main novelties and contributions of this paper can be summarized as follows: (i) generalization of pulseless controlled systems to impulsive controlled systems, which is helpful to increase the applicability of the model; (ii) through the use of C([0,T]) as a dense subspace in L1([0,T]), as per the results of functional analysis, the condition of conclusion is weakened; (iii) pointwise Jacobson-type necessary conditions have been obtained, which facilitates the calculations for and distinguishes the optimal singular control from other singular controls in the classical sense (see Definition 2.2).
The outline of this paper is as follows. Some preliminaries are proposed in Section 2. In Section 3, the integral-form second-order necessary conditions are derived. In Section 4, the pointwise Jacobson-type necessary conditions are given. An example is considered to elucidate the proposed main results in Section 5, and Section 6 concludes this paper.
In this section, we will present some preliminaries, which includes basic assumptions, the definition of the singular control in the classical sense, the solvability of impulsive systems, and a lemma that has been obtained via functional analysis.
Let B⊤ denote the transposition of a matrix B. Define C1([0,T]∖Λ,Rn)={y:[0,T]⟶Rn|y is continuous differential at t∈[0,T]∖Λ} and PCl([0,T],Rn)(PCr([0,T],Rn))={y:[0,T]⟶Rn|y is continuous at t∈[0,T]∖Λ,y is left (right) continuous and exists right (left) limit at t∈Λ}, obviously endowed with the norm ‖y‖PC=sup{‖y(t+)‖,‖y(t−)‖|t∈[0,T]}, PCl([0,T],Rn); also, PCr([0,T],Rn) denotes Banach spaces.
Let us assume the following:
(A1) U⊆Rm is a nonempty bounded convex open set.
(A2) The functions denoted by F=(fl):[0,T]×Rn×U⟶Rn+1 are measurable in t and twice continuously differentiable in (y,u); for any ρ>0, there exists a constant L(ρ)>0 such that, for all y,ˆy∈Rn and u,ˆu∈U with ‖y‖,‖ˆy‖,‖u‖,‖ˆu‖≤ρ, and for all t∈[0,T] such that
{‖F(t,y,u)−F(t,ˆy,ˆu)‖≤L(ρ)(‖y−ˆy‖+‖u−ˆu‖),‖Fy(t,y,u)−Fy(t,ˆy,ˆu)‖≤L(ρ)(‖y−ˆy‖+‖u−ˆu‖),‖Fu(t,y,u)−Fu(t,ˆy,ˆu)‖≤L(ρ)(‖y−ˆy‖+‖u−ˆu‖), |
there is a constant h>0 such that
‖F(t,y,u)‖≤h(1+‖y‖), for all (t,u)∈[0,T]×U, | (2.1) |
where Fy(t,y,u)⊤ (or Fu(t,y,u)⊤) denotes the Jacobi matrix of F in y (or u).
(A3) the functions denoted by ˜Ji=(JiG):Rn⟶Rn+1 (i=1,2,⋯,k) are twice continuously differentiable in x, and, for any ρ>0, there exists a constant L(ρ)>0 such that, for all y,ˆy∈Rn with ‖y‖,‖ˆy‖≤ρ, we have
{‖˜Ji(y)−˜Ji(ˆy)‖≤L(ρ)‖y−ˆy‖,‖˜Jix(y)−˜Jiy(ˆy)‖≤L(ρ)‖y−ˆy‖, |
where ˜Jiy(y)⊤ denotes the Jacobi matrix of ˜Ji in y.
Denote H(t)=l(t,y(t),u(t))+⟨f(t,y(t),u(t)),φ(t)⟩. To simplify the notation, [t] is used to replace (t,ˉy(t),¯φ(t),ˉu(t)) when evaluating the dynamics f and the Hamiltonian H; for example, f[t]=f(t,ˉy(t),ˉu(t)) and H[t]=l(t,ˉy(t),ˉu(t))+⟨f(t,ˉy(t),ˉu(t)),¯φ(t)⟩.
Remark 2.1. H(t)=l(t,y(t),u(t))+⟨f(t,y(t),u(t)),φ(t)⟩ is often referred to as the Hamiltonian function, or, simply, the Hamiltonian. The Hamiltonian function is represented in many reports as ˜H(t)=−l(t,y(t),u(t))+⟨f(t,y(t),u(t)),φ(t)⟩; the difference between the two is that Pontryagin's maximum principle maximizes the Hamiltonian ˜H or minimizes the Hamiltonian H. Accordingly, it also leads to a difference between positivity and negativity of the optimal inequality in the maximum principle. In this article, we use H to denote the Hamiltonian; in other words, we need to minimize the Hamiltonian H.
Now, we will prove an elementary theorem that will be useful in the following sections.
Theorem 2.2. Let (A1)–(A3) hold; for any fixed u∈Uad, the control system (1.2) has a unique solution yu∈PCl([0,T],Rn) given by
yu(t)=y0+∫t0f(s,yu(s),u(s))ds+∑0<ti<tJi(yu(ti)), | (2.2) |
and there exists a constant M=M(h,T,y0,J1(0),J2(0),⋯,Jk(0)) such that
‖yu(t)‖≤Mfor all(t,u)∈[0,T]×Uad. | (2.3) |
Proof. By the qualitative theory of differential equations, it follows from (A1)–(A3) that the system of equations
{˙y(t)=f(t,y(t),u(t)),t∈[0,t1],y(0)=y0, |
has a unique solution yu∈C([0,t1],Rn) given by
yu(t)=y0+∫t0f(s,yu(s),u(s))ds,t∈[0,t1], |
and
‖yu(t)‖≤eht(ht1+‖y0‖) for all (t,u)∈[0,t1]×Uad. |
Let
y1=yu(t1)+J1(yu(t1)), |
then, one also can infer that the system of equations
{˙y(t)=f(t,y(t),u(t)),t∈(t1,t2],y(t1+)=y1, |
has a unique solution yu∈C((t1,t2],Rn) given by
yu(t)=y0+∫t0f(s,yu(s),u(s))ds+J1(yu(t1)),t∈(t1,t2], |
and
‖yu(t)‖≤eh(t−t1)(h(t2−t1)+‖y1‖) for all (t,u)∈(t1,t2]×Uad. |
Using a step by step method, together with
‖Ji(yu(ti))‖≤L(‖yu(ti)‖)‖yu(ti)‖+‖Ji(0)‖,i=1,2,⋯,k, |
it is not difficult to claim that there exists a constant M such that (2.2) and (2.3) hold. Therefore, we have completed the proof of Theorem 2.2.
To establish the main conclusion of Theorem 4.6, we now introduce the definition of singular control in the classical sense (see Definition 2 in [14], as well as (4.4) in [25]).
Definition 2.3. We refer to the elements in the following equation as singular controls in the classical sense:
¯Uad={v(⋅)∈Uad|Hu(t,ˉy(t),v(t),¯φ(t))=0,Huu(t,y(t),v(t),φ(t))≡0}. |
Remark 2.4. If the Hamiltonian H is linear in control u, then ¯Uad is a singular control set in the classical sense.
Now, we shall introduce a fundamental lemma, as derived from functional analysis, that will be utilized to establish the necessary condition for Problem P.
Lemma 2.5. Let h(t) be an n-dimensional piecewise continuous vector value function on [0,T], and suppose that
∫T0⟨h(t),a(t)⟩dt=0, |
for all n-dimensional piecewise continuous vector value functions a(t) on [0,T]; then, h(t)=0 at all continuous moments of h(t) on [0,T].
Proof. Suppose that h(t) at some continuous moments ˉt when h(ˉt)≠0 because h(t) is an n-dimensional piecewise continuous vector value function; therefore, there exists an interval Iˉt at ˉt such that
h(t)≠0,t∈Iˉt. |
In this case, given the following:
a(t)={h(t),t∈Iˉt,0,t∈[0,T]∖Iˉt, |
then
∫T0h⊤(t)h(t)dt=∫Iˉth⊤(t)h(t)dt=∫Iˉt‖h(t)‖2dt>0. |
This is a contradiction. Therefore, we have finished the proof of Lemma 2.5.
The purpose of this section is to prove Theorem 3.1, which establishes the integral form of the second-order necessary conditions for Problem P. The basic idea is that the condition of being second-order variational and non-negative is a necessary condition for an optimal control problem. We will borrow the method adopted in [25] and introduce a linear impulsive adjoint matrix to prove it.
Theorem 3.1. Let (A1)–(A3) hold and ˉu represent the optimal control of J over Uad; it is necessary that there exist functions (¯y,¯φ,¯W,¯Φ)∈PCl([0,T],Rn)×PCr([0,T],Rn)×PCr([0,T],Rn×n)×PCl([0,T],Rn×n) such that the following equations and inequality hold:
{˙ˉy(t)=f(t,ˉy(t),ˉu(t)),t∈[0,T]∖Λ,ˉy(ti+)−ˉy(ti)=Ji(ˉy(ti)),ti∈Λ,ˉy(0)=y0; | (3.1) |
{˙¯φ(t)=−fy(t,ˉy(t),ˉu(t))¯φ(t)−ly(t,ˉy(t),ˉu(t)),t∈[0,T]∖Λ,¯φ(ti−)=¯φ(ti)+Jiy(ˉy(ti))¯φ(ti),ti∈Λ,¯φ(T)=Gy(ˉy(T)); | (3.2) |
{˙¯W(t)=−fy(t,ˉy(t),ˉu(t))¯W(t)−¯W(t)fy(t,ˉy(t),ˉu(t))⊤−¯φ(t)⊤fyy(t,ˉy(t),ˉu(t))−lyy(t,ˉy(t),ˉu(t)),t∈[0,T]∖Λ,¯W(ti−)=¯W(ti)+Jiy(ˉy(ti))¯W(ti)+¯W(ti+)Jiy(ˉy(ti))⊤+Jiy(ˉy(ti))¯W(ti)Jiy(ˉy(ti))⊤+¯φ(ti)⊤Jiyy(ˉy(ti)),ti∈Λ,¯W(T)=Gyy(ˉy(T)); | (3.3) |
{˙¯Φ(t)=fy(t,ˉy(t),ˉu(t))⊤¯Φ(t),t∈[0,T]∖Λ,¯Φ(ti+)=¯Φ(ti)+Jiy(ˉy(ti))⊤¯Φ(ti),ti∈Λ,¯Φ(0)=I; | (3.4) |
12∫T0⟨Huu[t][u(t)−ˉu(t)],u(t)−ˉu(t)⟩dt+∫T0∫t0⟨[Huy[t]+¯W(t)fu[t]⊤][u(t)−ˉu(t)],¯Φ(t)¯Φ(s)−1fu[s]T[u(s)−ˉu(s)]ds⟩dt≥0for allu(⋅)∈Uad. |
Proof. Now, let (ˉy(⋅),ˉu(⋅)) be the given optimal pair and ε∈(0,1]. For an arbitrary but fixed u(⋅)∈Uad, let uε(⋅)=ˉu(⋅)+ε(u(⋅)−ˉu(⋅)). It follows from the assumption (A1) that uε(⋅)∈Uad; according to Theorem 2.2, uε determines the unique allowed state yε(⋅); then, we can get from (A2), (A3) and Theorem 2.2 (see (2.2) and (2.3)) that
‖yε(t)−ˉy(t)‖≤∫t0‖f(s,yε(s),uε(s))−f(s,ˉy(s),ˉu(s))‖ds+∑0<ti<t‖Ji(yε(ti))−Ji(ˉy(ti))‖≤L(M)∫t0(‖yε(t)−ˉy(t)‖+ε‖u(t)−ˉu(t)‖)ds+L(M)∑0<ti<t‖yε(ti)−ˉy(ti)‖. |
Using the impulse integral inequality (see [1]), we have
limε→0‖yε−ˉy‖PC=0. | (3.5) |
Let
Y(t)=limε→0Yε(t)=limε→0yε(t)−y(t)ε. |
In the same way as for (3.5), it is easy to show that
limε→0‖Yε−Y‖PC=0, | (3.6) |
and Y solves the following system of variational equations:
{˙Y(t)=fy[t]⊤Y(t)+fu[t]⊤(u(t)−ˉu(t)),t∈[0,T]∖Λ,Y(ti+)=Y(ti)+Jiy(ˉy(ti))⊤Y(ti),ti∈Λ,Y(0)=0. | (3.7) |
To obtain the first-order necessary condition for Problem P, the following proposition will be used.
Proposition 3.2. Let (A2) and (A3) hold and ¯φ∈PCr([0,T],Rn) be the solution of the impulsive adjoint equation given by (3.2). Then
∫T0⟨fu(t,ˉy(t),ˉu(t))¯φ(t),u(t)−ˉu(t)⟩dt=∫T0⟨ly(s,ˉy(s),ˉu(s)),Y(s)⟩ds+⟨Gy(ˉy(T)),Y(T)⟩. | (3.8) |
Proof. Since C([0,T]) is a dense subspace in L1([0,T]), there exist function sequences {fαy}, {lαy}⊆C([0,T] such that
fαy(⋅)⟶fy(⋅,ˉy(⋅),ˉu(⋅)) and lαy(⋅)⟶ly(⋅,ˉy(⋅),ˉu(⋅)) in L1([0,T]) as α→∞. | (3.9) |
Moreover, it follows from (A3) that the system of linear impulsive differential equations given by
{˙φα(t)=−fαy(t)φα(t)−lαy(t),t∈[0,T]∖Λ,φα(ti−)=φα(ti)+Jiy(ˉy(ti))φα(ti),ti∈Λ,φα(T)=Gy(ˉy(T)), | (3.10) |
has a unique solution φα∈PCr([0,T],Rn)⋂C1([0,T]∖Λ,Rn), and that there is a constant β>0 such that
‖φα‖PC≤β for all α. |
Hence, we have
‖φα(t)−¯φ(t)‖≤∫Tt‖lαy(s)−ly(s,ˉy(s),ˉu(s))‖ds+∫Tt‖fαy(s)−fy(s,ˉy(s),ˉu(s))‖‖φα(s)‖ds+∫Tt‖fy(s,ˉy(s),ˉu(s))‖‖φα(s)−¯φ(s)‖ds+∑t<ti<T‖Jiy(ˉy(ti))‖‖φα(ti)−¯φ(ti)‖. |
Therefore, using the same method as for (3.5), it is not difficult to show that
limα→∞‖φα−¯φ‖PC=0. | (3.11) |
Consequently, we can infer from (3.7) and (3.10) that
∫T0⟨fu(t,ˉy(t),ˉu(t))φα(t),u(t)−ˉu(t)⟩dt=∫T0⟨φα(t),fu(t,ˉy(t),ˉu(t))⊤(u(t)−ˉu(t))⟩dt=∫T0⟨φα(t),˙Y(t)−f⊤y(t,ˉy(t),ˉu(t))Y(t)⟩dt=⟨φα(T),Y(T)⟩−k∑i=1[⟨φα(ti),Y(ti+)⟩−⟨φα(ti−),Y(ti)⟩]−∫T0⟨˙φα(t)+fαy(t)φα(t),Y(t)⟩dt=⟨Gy(ˉy(T)),Y(T)⟩−∫T0⟨˙φα(t)+fαy(t)φα(t),Y(t)⟩dt=⟨Gy(ˉy(T)),Y(T)⟩+∫T0⟨lαy(t),Y(t)⟩dt. |
Let α⟶∞ in the above expression; using (3.9) and (3.11), we have (3.8). Therefore, we have finished the proof of Proposition 3.2.
Based on the above proposition, we now continue to prove Theorem 3.1.
By the optimality of ˉu, one can ascertain from the assumptions (A2) and (A3), (3.2), (3.6), (3.7), and Proposition 3.2 (see (3.8)) that
0≤limε→0J(uε)−J(ˉu)ε=limε→0∫T0⟨∫10ly(t,ˉy(t)+τ(yε(t)−ˉy(t)),ˉu(t))dτ,Yε(t)⟩dt+limε→0∫T0⟨∫10lu(t,yε(t),ˉu(t)+τε(u(t)−ˉu(t)))dτ,u(t)−ˉu(t)⟩dt+limε→0⟨∫10Gy(ˉy(T)+τ(yε(T)−ˉy(T)))dτ,Yε(T)⟩=∫T0⟨ly(t,ˉy(t),ˉu(t)),Y(t)⟩dt+∫T0⟨lu(t,ˉy(t),ˉu(t)),u(t)−ˉu(t)⟩dt+⟨Gy(ˉy(T)),Y(T)⟩=∫T0⟨lu(t,ˉy(t),ˉu(t))+fu(t,ˉy(t),ˉu(t))¯φ(t),u(t)−ˉu(t)⟩dt, |
which leads to the following optimal inequality:
∫T0⟨lu(t,ˉy(t),ˉu(t))+fu(t,ˉy(t),ˉu(t))¯φ(t),u(t)−ˉu(t)⟩dt≥0. |
Because of the arbitrariness of u∈Uad, we get
∫T0⟨lu(t,ˉy(t),ˉu(t))+fu(t,ˉy(t),ˉu(t))¯φ(t),u(t)−ˉu(t)⟩dt=0. | (3.12) |
Moreover, combining this with (A1) and Proposition 3.2, it follows that, for all continuous time points t∈[0,T] of ˉu(t), we have
Hu[t]=lu(t,ˉy(t),ˉu(t))+⟨fu(t,ˉy(t),ˉu(t)),¯φ(t)⟩=0. |
We define
¯U(t)={v∈U∣lu(t,ˉy(t),v)+⟨fu(t,ˉy(t),v),¯φ(t)⟩=0}. | (3.13) |
We refer to this as the singular control region in the classical sense, which will be used later.
Let
Zε(⋅)≡(Zε1(⋅)Zε2(⋅)⋯Zεn(⋅))⊤=Yε(t)−Y(t)ε, |
and
Zk(t)=limε→0Zεk(t). |
In the same way as for (3.5), one can also claim from (A2) and (A3) that
limε→0‖Zε−Z‖PC([0,T])=0, | (3.14) |
and Z(⋅) denotes the solution of the following system of equations:
{˙Z(t)=fy(t,ˉy(t),ˉu(t))⊤Z(t)+12⟨fuu(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],[u(t)−ˉu(t)]⟩+⟨fuy(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],Y(t)⟩+12⟨fyy(t,ˉy(t),ˉu(t))Y(t),Y(t)⟩,t∈(0,T]∖Λ,Z(ti+)=Z(ti)+J⊤iy(ˉy(ti))Z(ti)+12⟨Jiyy(ˉy(ti))Y(ti),Y(ti)⟩,ti∈Λ,Z(0)=0. | (3.15) |
where
12⟨fuu(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],[u(t)−ˉu(t)]⟩=12{⟨f1uu(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],[u(t)−ˉu(t)]⟩⟨f2uu(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],[u(t)−ˉu(t)]⟩⋮⟨fnuu(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],[u(t)−ˉu(t)]⟩}, |
12⟨fyy(t,ˉy(t),ˉu(t))Y(t),Y(t)⟩=12{⟨f1yy(t,ˉy(t),ˉu(t))Y(t),Y(t)⟩⟨f2yy(t,ˉy(t),ˉu(t))Y(t),Y(t)⟩⋮⟨fnyy(t,ˉy(t),ˉu(t))Y(t),Y(t)⟩}, |
and
⟨fuy(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],Y(t)⟩={⟨f1uy(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)]),Y(t)⟩⟨f2uy(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],Y(t)⟩⋮⟨fnuy(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],Y(t)⟩}. |
Meanwhile, one can infer from (3.7) that X:=YY⊤ is the solution to the following system of equations:
{˙X(t)=fy(t,ˉy(t),ˉu(t))⊤X(t)+X(t)fy(t,ˉy(t),ˉu(t))+fu(t,ˉy(t),ˉu(t))⊤(u(t)−ˉu(t))Y⊤(t)+Y(t)(fu(t,ˉy(t),ˉu(t))⊤(u(t)−ˉu(t)))⊤,t∈[0,T]∖Λ,X(ti+)=X(ti)+J⊤iy(ˉy(ti))X(ti)+X(ti)Jiy(ˉy(ti))+J⊤iy(ˉy(ti))X(ti)Jiy(ˉy(ti)),ti∈Λ,X(0)=0. | (3.16) |
The subsequent proposition plays a crucial role in obtaining the integral form of the second-order necessary conditions.
Proposition 3.3. Let (A2) and (A3) hold and W∈PCr([0,T],Rn×n), ¯Φ∈PCl([0,T],Rn×n) be the solution of (3.3) and (3.4). Then,
∫T0⟨¯W(t)fu(t,ˉy(t),ˉu(t))⊤(u(t)−ˉu(t)),Y(t)⟩dt=12∫T0⟨[lyy(t,ˉy(t),ˉu(t))+¯φ(t)⊤fyy(t,ˉy(t),ˉu(t))]Y(t),Y(t)⟩dt+12⟨Gyy(ˉy(T))Y(T),Y(T)⟩+12k∑i=1⟨¯φ(ti)⊤Jiyy(ˉy(ti))Y(ti)),Y(ti))⟩, | (3.17) |
and
∫T0⟨(luy(t,ˉy(t),ˉu(t))+¯φ(t)⊤fuy(t,ˉy(t),ˉu(t))+¯W(t)fu(t,ˉy(t),ˉu(t))T)[u(t)−ˉu(t)],Y(t)⟩dt=∫T0dt∫t0⟨(luy(t,ˉy(t),ˉu(t))+¯φ(t)⊤fuy(t,ˉy(t),ˉu(t))+¯W(t)fu(t,ˉy(t),ˉu(t))T)[u(t)−ˉu(t)],.¯Φ(t)¯Φ(s)−1fu(s,ˉy(s),ˉu(s))T[u(s)−ˉu(s)]⟩ds. | (3.18) |
Proof. Since C([0,T]) is a dense subspace in L1([0,T]), there exist function sequences {fαy}, {fαyy}, {fαuy}, {fαu}, {lαy}, {lαyy}, {lαuy}⊆C([0,T] such that
{fαy(⋅)⟶fy(⋅,ˉy(⋅),ˉu(⋅)),lαy(⋅)⟶ly(⋅,ˉy(⋅),ˉu(⋅)),fαyy(⋅)⟶fyy(⋅,ˉy(⋅),ˉu(⋅)),lαyy(⋅)⟶lyy(⋅,ˉy(⋅),ˉu(⋅)),fαuy(⋅)⟶fuy(⋅,ˉy(⋅),ˉu(⋅)),lαuy(⋅)⟶luy(⋅,ˉy(⋅),ˉu(⋅)),fαu⟶fu, in L1([0,T]) as α→∞. | (3.19) |
Consequently, by (A3) and (3.19), one can infer that the systems of linear impulsive matrix differential equations given by
{˙Wα(t)=−fαy(t)Wα(t)−Wα(t)fαy(t)⊤−¯φ(t)⊤fαyy(t)−lαyy(t),t∈[0,T]∖Λ,Wα(ti−)=Wα(ti+)+Jiy(ˉy(ti))Wα(ti+)+Wα(ti+)Jiy(ˉy(ti))⊤+Jiy(ˉy(ti))Wα(ti+)Jiy(ˉy(ti))⊤+¯φ(ti)⊤Jiyy(ˉy(ti)),ti∈Λ,Wα(T)=Gyy(ˉy(T)), | (3.20) |
and
{˙Φα(t)=fαy(t)⊤Φα(t),t∈[0,T]∖Λ,Φα(ti+)=Φα(ti)+Jiy(ˉy(ti))⊤Φα(ti),ti∈Λ,Φα(0)=I, | (3.21) |
each have a unique solution Wα∈PCr([0,T],Rn×n)⋂C1([0,T]∖Λ,Rn×n) and Φα∈PCl([0,T],Rn×n)⋂C1([0,T]∖Λ,Rn×n), respectively. Not only that, there exists a constant γ>0 such that
‖Wα‖PC≤γand‖Φα‖PC≤γ for all α. |
Moreover, we have
‖Wα(t)−¯W(t)‖≤∫Tt‖lαyy(s)−lyy(s,ˉy(s),ˉu(s))‖ds+∫Tt‖¯φ(s)⊤‖‖fαyy(s)−fyy(s,ˉy(s),ˉu(s))‖ds+2γ∫Tt‖fαy(s)−fy(s,ˉy(s),ˉu(s))‖ds+2∫Tt‖Wα(s)−¯W(s)‖‖f⊤y(s,ˉy(s),ˉu(s))‖ds+2∑t<ti<T(‖Jiy(ˉy(ti))‖+‖Jiy(ˉy(ti))‖2)‖Wα(ti)−¯W(ti)‖, |
and
‖Φα(t)−¯Φ(t)‖≤γ∫t0‖fαy(s)−fy(s,ˉy(s),ˉu(s))‖ds+∫t0‖fy(s,ˉy(s),ˉu(s))‖‖Φα(s)−¯Φ(s)‖ds+∑0<ti<t‖Jiy(ˉy(ti))‖‖Φα(ti)−¯Φ(ti)‖. |
In the same way as for (3.5), we obtain
limα→∞‖Wα−¯W‖PC=0andlimα→∞‖Φα−¯Φ‖PC=0. | (3.22) |
In addition, it is obvious from (3.20) that W⊤α is also a solution of (3.20). This means that
W⊤α(t)=Wα(t) for all t∈[0,T]. | (3.23) |
Since tr(AB)=tr(BA) for all k×j matrix A and j×k matrix B, we can get from (3.2), (3.7), (3.16), (3.19), (3.20), (3.22), and (3.23) that
2∫T0⟨¯W(t)fu(t,ˉy(t),ˉu(t))⊤(u(t)−ˉu(t)),Y(t)⟩dt=2limα→∞∫T0⟨Wα(t)(˙Y(t)−fy(t,ˉy(t),ˉu(t))⊤Y(t)),Y(t)⟩dt=limα→∞tr∫T0[Wα(t)(˙Y(t)−fy(t,ˉy(t),ˉu(t))⊤Y(t))Y⊤(t)+Wα(t)Y(t)(˙Y(t)⊤−Y(t)⊤fy(t,ˉy(t),ˉu(t)))]dt=limα→∞tr∫T0Wα(t)[˙X(t)−fy(t,ˉy(t),ˉu(t))⊤X(t)−X(t)fy(t,ˉy(t),ˉu(t))]dt=limα→∞{−∫T0⟨[˙Wα(t)+Wα(t)fy(t,ˉy(t),ˉu(t))⊤+fy(t,ˉy(t),ˉu(t))Wα(t)]Y(t),Y(t)⟩dt+⟨Wα(T)Y(T),Y(T)⟩+k∑i=1[⟨Wα(ti−)Y(ti),Y(ti)⟩−⟨Wα(ti)Y(ti+),Y(ti+)⟩]}=limα→∞{−∫T0⟨[˙Wα(t)+Wα(t)fαy(t)⊤+fαy(t)Wα(t)]Y(t),Y(t)⟩dt+⟨Wα(T)Y(T),Y(T)⟩+k∑i=1[⟨(Wα(ti−)−Wα(ti))Y(ti),Y(ti)⟩−⟨(Jiy(ˉy(ti))Wα(ti)+Wα(ti)Jiy(ˉy(ti))⊤)Y(ti),Y(ti)⟩−⟨Jiy(ˉy(ti))Wα(ti)Jiy(ˉy(ti))⊤Y(ti),Y(ti)⟩]}=limα→∞{∫T0⟨[lαyy(t)+¯φ(t)⊤fαyy(t)]Y(t),Y(t)⟩dt+⟨Wα(T)Y(T),Y(T)⟩+k∑i=1[⟨(Wα(ti−)−Wα(ti))Y(ti),Y(ti)⟩−⟨(Jiy(ˉy(ti))Wα(ti)+Wα(ti)Jiy(ˉy(ti))⊤)Y(ti),Y(ti)⟩−⟨Jiy(ˉy(ti))Wα(ti)Jiy(ˉy(ti))⊤Y(ti),Y(ti)⟩]}=∫T0⟨[lyy(t,ˉy(t),ˉu(t))+¯φ(t)⊤fyy(t,ˉy(t),ˉu(t))]Y(t),Y(t)⟩dt+⟨Gyy(ˉy(T))Y(T),Y(T)⟩+k∑i=1⟨¯φ(ti)⊤Jiyy(ˉy(ti))Y(ti),Y(ti)⟩, |
i.e., (3.17) holds.
Now, let us prove (3.18). By (3.4), (3.7), (3.21), and (3.22), we have
Y(t)=¯Φ(t)∫t0¯Φ(s)−1fu[s]⊤(u(s)−ˉu(s))ds=limα→0Φα(t)∫t0Φα(s)−1fu[s]⊤(u(s)−ˉu(s))ds, |
which means that (3.18) holds. Therefore, we have finished the proof of Proposition 3.3.
Based on the above propositions, we now continue to prove Theorem 3.1.
Since ˉu represents optimal control of J over Uad, together with Proposition 3.2 (see (3.8)) and (3.12), we have
∫T0⟨ly(t,ˉy(t),ˉu(t)),Y(t)⟩ds+∫T0⟨lu(t,ˉy(t),ˉu(t)),u(t)−ˉu(t)⟩dt+⟨Gy(ˉy(T)),Y(T)⟩=∫T0⟨lu(t,ˉy(t),ˉu(t))+fu(t,ˉy(t),ˉu(t))φ(t),u(t)−ˉu(t)⟩dt=0 for all u∈Uad. | (3.24) |
Taken together with (A2), (A3), and (3.24), one can get
J(uε(⋅))−J(ˉu(⋅))ε=∫T0⟨∫10(lu(t,yε(t),ˉu(t)+τε(u(t)−ˉu(t)))−lu(t,ˉy(t),ˉu(t)))dτ,u(t)−ˉu(t)⟩dt+∫T0⟨∫10(ly(t,ˉy(t)+τ(yε(t)−ˉy(t)),ˉu(t))−ly(t,ˉy(t),ˉu(t)))dτ,Yε(t)⟩dt+∫T0⟨ly(t,ˉy(t),ˉu(t)),Yε(t)−Y(t)⟩dt+⟨Gy(ˉy(T)),Yε(T)−Y(T)⟩+⟨∫10(Gy(ˉy(T)+τ(yε(T)−ˉy(T)))−Gy(ˉy(T)))dτ,Yε(T)⟩=ε∫T0⟨∫10τ∫10luu(t,yε(t),ˉu(t)+ντε(u(t)−ˉu(t)))dνdτ[u(t)−ˉu(t)],u(t)−ˉu(t)⟩dt+ε∫T0⟨∫10luy(t,ˉy(t)+τ(yε(t)−ˉy(t)),ˉu(t))(u(t)−ˉu(t))dτ,Yε(t)⟩dt+ε∫T0⟨∫10τ∫10lyy(t,ˉy(t)+ντ(yε(t)−ˉy(t)),ˉu(t))dνdτYε(t),Yε(t)⟩dt+ε∫T0⟨ly(t,ˉy(t),ˉu(t)),Zε(t)⟩dt+ε⟨Gy(ˉy(T)),Zε(T)⟩+ε⟨∫10τ∫10Gyy(ˉy(T)+ντ(yε(T)−ˉy(T)))dνdτYε(T),Yε(T)⟩. |
Then, combining this with (3.6) and (3.14), the above expression, (A2), and (A3) leads to the following:
12∫T0⟨luu(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],u(t)−ˉu(t)⟩dt+∫T0⟨luy(t,ˉy(t),ˉu(t))(u(t)−ˉu(t)),Y(t)⟩ds+12∫T0⟨lyy(t,ˉy(t),ˉu(t))Y(t),Y(t)⟩dt+12⟨Gyy(ˉy(T))Y(T),Y(T)⟩+∫T0⟨ly(t,ˉy(t),ˉu(t)),Z(t)⟩ds+⟨Gy(ˉy(T)),Z(T)⟩≥0 for all u∈Uad. |
By (3.2), we have
12∫T0⟨luu(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)],u(t)−ˉu(t)⟩dt+∫T0⟨luy(t,ˉy(t),ˉu(t))(u(t)−ˉu(t)),Y(t)⟩ds+12∫T0⟨lyy(t,ˉy(t),ˉu(t))Y(t),Y(t)⟩dt+12⟨Gyy(ˉy(T))Y(T),Y(T)⟩−∫T0⟨˙¯φ(t)+fy(t,ˉy(t),ˉu(t))¯φ(t),Z(t)⟩ds+⟨Gy(ˉy(T)),Z(T)⟩≥0 for all u∈Uad. | (3.25) |
Since C([0,T]) is a dense subspace in L1([0,T]), there exist function sequences {uα},{fαuu}⊆C([0,T]) such that
limα→∞‖uα−[u−ˉu]‖L1=0 and limα→∞‖fαuu(⋅)−fuu(⋅,ˉy(⋅),ˉu(⋅))‖L1=0. | (3.26) |
It follows immediately from (3.15), (3.19), and (3.26) that the system of equations given by
{˙Zα(t)=fαy(t)⊤Zα(t)+12uα(t)⊤fαuu(t)uα(t)+uα(t)⊤fαyu(t)Y(t)+12Y(t)⊤fαyy(t)Y(t),t∈[0,T]∖Λ,Zα(ti+)=Zα(ti)+12Y(ti))⊤Jiyy(ˉy(ti))Y(ti))+Jiy(ˉy(ti))⊤Zα(ti),ti∈Λ,Zα(0)=0, | (3.27) |
has a unique solution Zα∈PCl([0,T],Rn)⋂C1([0,T]∖Λ,Rn) and
limα→∞‖Zα−Z‖PC=0, | (3.28) |
where Z(⋅) is the solution to (3.15).
Moreover, we can infer from (3.19) and (3.26)–(3.28) that
−∫T0⟨˙¯φ(t)+fy(t,ˉy(t),ˉu(t))¯φ(t),Z(t)⟩dt+⟨Gy(ˉy(T)),Z(T)⟩=−limα→∞∫T0⟨˙¯φ(t)+fαy(t)¯φ(t),Zα(t)⟩dt+⟨Gy(ˉy(T)),Z(T)⟩=limα→∞{∫T0⟨¯φ(t),˙Zα(t)−fα⊤y(t)Zα(t)⟩dt−k∑i=1[⟨¯φ(ti−),Zα(ti)⟩−⟨¯φ(ti),Zα(ti+)⟩]}=12limα→∞∫T0⟨¯φ(t),uα(t)⊤fαuu(t)uα(t)+2uα(t)⊤fαyu(t)Y(t)+Y(t)⊤fαyy(t)Y(t)⟩dt+12k∑i=1⟨¯φ(ti),Y(ti))⊤Jiyy(ˉy(ti))Y(ti))⟩=12∫T0⟨¯φ(t),[u(t)−ˉu(t)]⊤fuu(t,ˉy(t),ˉu(t))[u(t)−ˉu(t)]⟩dt+∫T0⟨¯φ(t),(u(t)−ˉu(t))⊤(t)fuy(t,ˉy(t),ˉu(t))Y(t)⟩dt+12∫T0⟨¯φ(t),Y(t)⊤fyy(t,ˉy(t),ˉu(t))Y(t)⟩dt+12k∑i=1⟨¯φ(ti),Y(ti)⊤Jiyy(ˉy(ti))Y(ti)))⟩. | (3.29) |
Taken together with (3.29), we deduce from (3.25) that
12∫T0⟨[luu(t,ˉy(t),ˉu(t))+¯φ(t)⊤fuu(t,ˉy(t),ˉu(t))][u(t)−ˉu(t)],u(t)−ˉu(t)⟩ds+∫T0⟨[luy(t,ˉy(t),ˉu(t))+¯φ(t)⊤fuy(t,ˉy(t),ˉu(t))][u(t)−ˉu(t)],Y(t)⟩dt+12∫T0⟨[lyy(t,ˉy(t),ˉu(t))+¯φ(t)⊤fyy(t,ˉy(t),ˉu(t))]Y(t),Y(t)⟩dt+12⟨Gyy(ˉy(T))Y(T),Y(T)⟩+12k∑i=1⟨¯φ(ti)⊤Jiyy(ˉy(ti))Y(ti)),Y(ti))⟩≥0 for all u∈Uad. | (3.30) |
The following inequality follows from (3.30) and (3.17):
12∫T0⟨[luu(t,ˉy(t),ˉu(t))+¯φ(t)⊤fuu(t,ˉy(t),ˉu(t))][u(t)−ˉu(t)],u(t)−ˉu(t)⟩dt+∫T0⟨[luy(t,ˉy(t),ˉu(t))+¯φ(t)⊤fuy(t,ˉy(t),ˉu(t))+¯W(t)fu(t,ˉy(t),ˉu(t))⊤][u(t)−ˉu(t)],.Y(t)⟩dt≥0 for all u∈Uad. | (3.31) |
Then, by (3.18) and (3.31), we can show that
12∫T0⟨[luu(t,ˉy(t),ˉu(t))+¯φ(t)⊤fuu(t,ˉy(t),ˉu(t))][u(t)−ˉu(t)],u(t)−ˉu(t)⟩dt+∫T0∫t0⟨[luy(t,ˉy(t),ˉu(t))+¯φ(t)⊤fuy(t,ˉy(t),ˉu(t))+¯W(t)fu(t,ˉy(t),ˉu(t))⊤][u(t)−ˉu(t)],.¯Φ(t)¯Φ(s)−1fu(s,ˉy(s),ˉu(s))T[u(s)−ˉu(s)]ds⟩dt≥0 for all u∈Uad, | (3.32) |
Thus, we have finished the proof of Theorem 3.1.
Remark 3.4. Theorem 3.1 does not establish whether the optimal control of Problem P is singular or nonsingular; it is a unified conclusion, similar to the equations in (4.5.2) of Theorem 4.2 in [15]. Therefore, Theorem 3.1 is a generalization of Theorem 4.2 in [15] to impulsive controlled systems. Based on Theorem 3.1, singular control and nonsingular control in the classical sense can be considered under a unified framework for Problem P.
In this section, on the basis of Theorem 3.1 in the previous section, we first obtain the Legendre-Clebsch condition; then, we give a corollary for the integral form of the second-order necessary optimality conditions for optimal singular control; finally, we give the pointwise Jacobson type necessary conditions and the pointwise Legendre-Clebsch condition.
Corollary 4.1. Let (A1)–(A3) hold and ˉu denote the optimal control of J over Uad; it is necessary that there exist a pair of functions (¯y,¯φ)∈PCl([0,T],Rn)×PCr([0,T],Rn) such that (3.1), (3.2), and
Huu[t]≥0,for all the continuous time ofˉu(t),t∈[0,T], | (4.1) |
hold.
Proof. To prove Corollary 4.1, let ˉu(⋅)∈Uad; take the special control variational problem as follows:
u(t)−ˉu(t)={0,t∈[t0,ˉt),h,t∈[ˉt,ˉt+ε),0,t∈[ˉt+ε,T], | (4.2) |
where h∈Rr is a constant vector, ε is a sufficiently small positive number, ˉt is any continuous time of ˉu(t). For u(t)−ˉu(t), ¯Φ(t) satisfies (3.4); then, the solution Y(t) of the variational problem (3.7) is given by
Y(t)={0,t∈[t0,ˉt),∫tˉt¯Φ(t)¯Φ(s)−1fu[s]⊤hds,t∈[ˉt,ˉt+ε),∫ˉt+εˉt¯Φ(t)¯Φ(s)−1fu[s]⊤hds,t∈[ˉt+ε,T]. |
Then
‖Y(t)‖={0,t∈[t0,ˉt),O(ε),t∈[ˉt,T], |
where limε→0O(ε)ε=C(nonzero constant). Utilizing the continuity of the mean value theorem for integrals, we have
∫T0Y(t)⊤Hyy[t]Y(t)dt=∫ˉt+εˉtY(t)⊤Hyy[t]Y(t)dt=O(ε2)=o(ε),∫T0Y(t)⊤Huy[t][u(t)−ˉu(t)]dt=∫ˉt+εˉtY(t)⊤Huy[t]hdt=εY(ˉt)⊤Huy[ˉt]h+o(ε),∫T0[u(t)−ˉu(t)]⊤Huu[t][u(t)−ˉu(t)]dt=∫ˉt+εˉth⊤Huu[t]hdt=εh⊤Huu[ˉt]h+o(ε). | (4.3) |
Substituting (4.3) into (3.32), we have
εh⊤Huu[t]h+o(ε)≥0. |
Observe that h⊤Huu[t]h is independent of ε and ε can be arbitrarily small; we have
h⊤Huu[t]h≥0. |
Since h∈Rr is an arbitrary vector and t denotes arbitrary continuous time, (4.1) holds.
Remark 4.2. Condition (4.1) represents the Legendre-Clebsch condition for the optimal control problem Problem P. At the same time, it also shows the rationality of the conventional hypothesis Huu[t]≥0,t∈[t0,T].
Remark 4.3. For the LQ problem, where R[t]=Huu[t], when R[t]>0,∀t∈[t0,T], the problem is called a nonsingular problem; when R[t]=0,∀t∈[t0,T], the problem is called the totally singular case; when R[t]≥0,∀t∈[t0,T], the problem is called the partially singular case (see Definitions (4.4)–(4.6) in [15]). Whether or not the optimal control problem is singular, and regardless of the kind of singularity, this classification standard is often adopted.
The following is a corollary of Theorem 3.1 in the case in which Huu[t]≡0, that is, Problem P is a totally singular problem according to the definitions in [15].
Corollary 4.4. Let (A1)–(A3) hold and ˉu(⋅)∈¯Uad denotes the optimal singular control of J over Uad; it is necessary that there exist functions (¯y,¯φ,¯W,¯Φ)∈PCl([0,T],Rn)×PCr([0,T],Rn)×PCr([0,T],Rn×n)×PCl([0,T],Rn×n) that satisfy (3.1)–(3.4) and
∫T0dt∫t0⟨(¯W(t)fu[t]⊤+Huy[t])[u(t)−ˉu(t)],¯Φ(t)¯Φ(s)−1fu[s]⊤[u(s)−ˉu(s)]⟩ds≥0,∀u∈¯Uad. |
Remark 4.5. Corollary 4.4 is a similar conclusion of Theorem 4.3 in [25] for Problem P.
The following is the pointwise Jacobson-type second-order necessary optimality condition for Problem P. The conclusion is similar to that of Theorem 4.3 in [25]. Note that the set ¯U(t) (see (3.13)) of values for v differs from the set Rm of values that is described in [32]; thus, it essentially confirms the pointwise characteristic. The author of [25] has proven similar conclusions under the weaker condition, suggesting that the control region U is a Polish space. In fact, under our basic assumption (A1), we can use the method on page 93 in [15] to prove it.
Theorem 4.6. Let (A1)–(A3) hold and ˉu(⋅)∈¯Uad denote the optimal singular control of J over Uad; it is necessary that there exist functions (¯y,¯φ,¯W,¯Φ)∈PCl([0,T],Rn)×PCr([0,T],Rn)×PCr([0,T],Rn×n)×PCl([0,T],Rn×n) that satisfy (3.1)–(3.4), and, for all continuous time of ˉu(t),t∈[0,T], we have
⟨(¯W(t)fu[t]⊤+Huy[t])[v−ˉu(t)],fu[t]⊤[v−ˉu(t)]⟩≥0,∀v∈¯U(t). | (4.4) |
Proof. {Note that Definition 2.3 implies that Huu[t]≡0for allt∈[0,T], apply the same control perturbation as for (4.2);} we have
∫T0⟨(¯W(t)fu[t]⊤+Hyu[t])h,Y(t)⟩dt=∫ˉt+εˉt⟨(¯W(t)fu[t]⊤+Hyu[t])h,Y(t)⟩dt, | (4.5) |
and the dominant term in the expansion of (4.5) for sufficiently small ε is given by
(ε)2⟨(¯W(t)fu[t]⊤+Hyu[t])h,fu[t]⊤h⟩|ˉt. |
Since ˉt can be chosen as any continuous time of ˉu(t),t∈[0,T], let h=v−ˉu(t),∀v∈¯U(t); thus, (4.4) holds and we have finished the proof of Theorem 4.6.
Using the same idea as in Theorem 4.6, we can also obtain the pointwise Legendre-Clebsch necessary optimality condition corresponding to Corollary 4.1.
Corollary 4.7. Let (A1)–(A3) hold and ˉu denote the optimal control of J over Uad; it is necessary that there exist a pair of functions (¯y,¯φ)∈PCl([0,T],Rn)×PCr([0,T],Rn) such that (3.1), (3.2), and, for all continuous time of ˉu(t),t∈[0,T]
(v−ˉu(t))⊤Huu[t](v−ˉu(t))≥0,∀v∈¯U(t), | (4.6) |
hold.
Remark 4.8. Comparing Corollaries 4.7 and 4.1, it can be found that if the pointwise condition is satisfied, Huu≥0 is not required, as only (4.6) needs to be satisfied.
In this section, we will give an example to illustrate the effectiveness of Theorem 4.6.
Let
minJ(u(⋅))=y2(1), |
subject to
{˙y1(t)=u,t∈[0,1]∖0.5,˙y2(t)=−y21,t∈[0,1]∖0.5,y1(0.5+)=y1(0.5)+y1(0.5),y2(0.5+)=y2(0.5),y1(0)=0,y2(0)=0. |
Obviously, the Hamiltonian H(t,y,u,φ)=φ1u−φ2y21 satisfies the conditions for linear control, as denoted by u and u∈U=R. According to Remark (2.4), the problem is singular, i.e., ¯U(t)=R. It is not difficult to assert that ˉu≡0 denotes singular control; this is because, by ˉu≡0, we can get that ˉy1≡0,ˉy2≡0, and ¯φ1=0,t∈[0,1]; consequently, H(t,y,u,φ)=¯φ1ˉu−¯φ2ˉy21≡0, Hu≡0, and Huu≡0; by Definition 2.3, ˉu≡0,t∈[0,1] denotes singular control. The question is whether it is optimal singular control. Now, let us use Theorem 4.6 to determine that it must not be optimal singular control.
By (3.2), we have
{˙¯φ1(t)=2¯φ2ˉy1,t∈[0,1]∖0.5,˙¯φ2(t)=0,t∈[0,1]∖0.5,¯φ1(0.5−)=¯φ1(0.5),¯φ2(0.5−)=¯φ1(0.5)+¯φ2(0.5),¯φ1(1)=0,¯φ2(1)=1. |
Using (3.3), we have
˙¯W(t)=[˙¯w11(t)˙¯w12(t)˙¯w21(t)˙¯w22(t)]=−[0−2ˉy100][¯w11(t)¯w12(t)¯w21(t)¯w22(t)]−[¯w11(t)¯w12(t)¯w21(t)¯w22(t)][00−2ˉy10]+[¯φ2(t)000], |
and
[¯w11(1)¯w12(1)¯w21(1)¯w22(1)]=[0000], |
and
¯W(0.5−)=[¯w11(0.5−)¯w12(0.5−)¯w21(0.5−)¯w22(0.5−)]=[0010][¯w11(0.5)¯w12(0.5)¯w21(0.5)¯w22(0.5)]+[¯w11(0.5)¯w12(0.5)¯w21(0.5)¯w22(0.5)][0100]+[0010][¯w11(0.5)¯w12(0.5)¯w21(0.5)¯w22(0.5)][0100]. |
Substituting ˉu≡0,ˉy1≡0, and ˉy2≡0 directly into the above equations, the following results can be obtained directly
{¯φ1(t)≡0,¯φ2(t)≡1,t∈[0,1], |
and
{¯w11=t−1,t∈[0,1],¯w12=¯w21=¯w22={0,t∈[0.5,1],−0.5,t∈[0,0.5). |
By (4.4), the necessary condition for the singular control ˉu≡0 to be the optimal control scheme is given by
(t−1)v2≥0,∀v∈¯U(t). | (5.1) |
But, because of the above equations, (5.1) cannot be true for arbitrary but fixed t∈(0,1). Therefore, regarding singular control ˉu=0, according to Theorem 4.6, it must not be optimal singular control.
In this paper, we have investigated the pointwise Jacobson type necessary conditions for Problem P. By introducing an impulsive linear matrix Riccati differential equation, we have derived the integral representation of the functional second-order variational equation. On this basis, we obtained the integral form of the second-order necessary conditions and the pointwise Jacobson type necessary conditions for optimal singular control in the classical sense. Incidentally, the Legendre-Clebsch condition and the pointwise Legendre-Clebsch condition were also obtained. These conclusions have been derived under weaker conditions, thereby enriching existing conclusions. In the future, we will continue to research the pointwise Jacobson-type second-order necessary optimality conditions in the Pontryagin sense.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors are grateful to the anonymous referees for their helpful comments and valuable suggestions which have improved the quality of the manuscript. This work was supported by the National Natural Science Foundation of China (No. 12061021 and No. 11161009).
The authors declare that there is no conflict of interest.
[1] |
Y. Aghakhani, A. P. Bagshaw, C. G. Bénar, C. Hawco, F. Andermann, F. Dubeau, et al., fMRI activation during spike and wave discharges in idiopathic generalized epilepsy, Brain, 127 (2004), 1127–1144. https://doi.org/10.1093/brain/awh136 doi: 10.1093/brain/awh136
![]() |
[2] |
P. N. Taylor, G. Baier, A spatially extended model for macroscopic spike-wave discharges, J. Comput. Neurosci., 31 (2011), 679–684. https://doi.org/10.1007/s10827-011-0332-1 doi: 10.1007/s10827-011-0332-1
![]() |
[3] |
J. G. Milton, Epilepsy as a dynamic disease: a tutorial of the past with an eye to the future, Epilepsy Behav., 18 (2010), 33–34. https://doi.org/10.1016/j.yebeh.2010.03.002 doi: 10.1016/j.yebeh.2010.03.002
![]() |
[4] |
F. H. L. da Silva, W. Blanes, S. N. Kalitzin, J. Parra, P. Suffczynski, D. N. Velis, Dynamical diseases of brain systems: different routes to epileptic seizures, IEEE Trans. Biomed. Eng., 50 (2003), 540–548. https://doi.org/10.1109/TBME.2003.810703 doi: 10.1109/TBME.2003.810703
![]() |
[5] |
F. H. L. da Silva, J. P. Pijn, W. J. Wadman, Dynamics of local neuronal networks: control parameters and state bifurcations in epileptogenesis, Progr. Brain Res., 102 (1994), 359–370. https://doi.org/10.1016/S0079-6123(08)60552-X doi: 10.1016/S0079-6123(08)60552-X
![]() |
[6] |
M. C. Mackey, L. Glass, Oscillation and chaos in physiological control systems, Science, 197 (1977), 287–289. https://doi.org/10.1126/science.267326 doi: 10.1126/science.267326
![]() |
[7] |
L. Glass, M. C. Mackey, Pathological conditions resulting from instabilities in physiological control systems, Ann. N. Y. Acad. Sci., 316 (1979), 214–235. https://doi.org/10.1111/j.1749-6632.1979.tb29471.x doi: 10.1111/j.1749-6632.1979.tb29471.x
![]() |
[8] |
Z. C. Yang, D. G. Fan, Q. Y. Wang, G. M. Luan, Sharp decrease in the Laplacian matrix rank of phase-space graphs: a potential biomarker in epilepsy, Cognit. Neurodyn., 15 (2021), 649–659. https://doi.org/10.1007/s11571-020-09662-x doi: 10.1007/s11571-020-09662-x
![]() |
[9] | P. Gloor, Neurophysiological basis of generalized seizures termed centrocephalic, in The Physiopathogenesis of the Epilepsies, Springfield, (1969), 209–236. Available from: https://cir.nii.ac.jp/crid/1573950399116067840. |
[10] |
W. W. Lytton, Computer modeling of epilepsy, Nat. Rev. Neurosci., 9 (2008), 626–637. https://doi.org/10.1038/nrn2416 doi: 10.1038/nrn2416
![]() |
[11] | I. Soltesz, K. Staley, Computational Neuroscience in Epilepsy, Academic Press, San Diego, 2008. https://doi.org/10.1016/B978-0-12-373649-9.X5001-7 |
[12] |
H. R. Wilson, J. D. Cowan, Excitatory and inhibitory interactions in localized populations of model neurons, Biophys. J., 12 (1972), 1–24. https://doi.org/10.1016/S0006-3495(72)86068-5 doi: 10.1016/S0006-3495(72)86068-5
![]() |
[13] | J. G. Milton, P. H. Chu, J. D. Cowan, Spiral waves in integrate-and-fire neural network, in Advances in Neural Information Processing Systems, 5 (1992), 1001–1006. Available from: https://proceedings.neurips.cc/paper_files/paper/1992/file/07a96b1f61097ccb54be14d6a47439b0-Paper.pdf. |
[14] |
S. A. Hou, D. G. Fan, Q. Y. Wang, Regulating absence seizures by tri-phase delay stimulation applied to globus pallidus internal, Appl. Math. Mech., 43 (2022), 1399–1414. https://doi.org/10.1007/s10483-022-2896-7 doi: 10.1007/s10483-022-2896-7
![]() |
[15] |
R. P. Lesser, S. H. Kim, L. Beyderman, D. L. Miglioretti, W. R. Webber, M. Bare, et al., Brief bursts of pulse stimulation terminate afterdischarges caused by cortical stimulation, Neurology, 53 (1999), 1081–2073. https://doi.org/10.1212/WNL.53.9.2073 doi: 10.1212/WNL.53.9.2073
![]() |
[16] |
G. K. Motanedi, R. P. Lesser, D. L. Miglioretti, M. M. Yuko, B. Gordon, W. R. S. Webber, et al., Optimizing parameters for terminating cortical afterdischarges with pulse stimulation, Epilepsia, 43 (2002), 836–846. https://doi.org/10.1046/j.1528-1157.2002.24901.x doi: 10.1046/j.1528-1157.2002.24901.x
![]() |
[17] |
A. Ashourvan, S. Pequito, A. N. Khambhati, F. Mikhail, S. N. Baldassano, K. A. Davis, et al., Model-based design for seizure control by stimulation, J. Neural Eng., 17 (2020), 836–846. https://doi.org/10.1088/1741-2552/ab7a4e doi: 10.1088/1741-2552/ab7a4e
![]() |
[18] | J. Guckenheimer, P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, Springer-Verlag, New York, 1983. https://doi.org/10.1007/978-1-4612-1140-2 |
[19] | S. H. Strogatz, Nonlinear Dynamics and Chaos, Addison-Wesley, New York, 1994. Available from: http://hdl.handle.net/1813/97. |
[20] |
H. S. Haghighi, A. H. Markazi, Control of epileptic seizures by electrical stimulation: a model-based study, Biomed. Phys. Eng. Express, 7 (2021), 065009. https://doi.org/10.1088/2057-1976/ac240d doi: 10.1088/2057-1976/ac240d
![]() |
[21] |
P. Suffczynski, S. Kalitzin, F. H. L. da Silva, Dynamics of non-convulsive epileptic phenomena modeled by a bistable neuronal network, Neuroscience, 126 (2004), 467–-484. https://doi.org/10.1016/j.neuroscience.2004.03.014 doi: 10.1016/j.neuroscience.2004.03.014
![]() |
[22] |
P. N. Taylor, Y. J. Wang, M. Goodfellow, J. Dauwels, F. Moeller, U. Stephani, et al., A computational study of stimulus driven epileptic seizure abatement, PLoS One, 9 (2014), e114316. https://doi.org/10.1371/journal.pone.0114316 doi: 10.1371/journal.pone.0114316
![]() |
[23] |
M. Breakspear, J. A. Roberts, J. R. Terry, S. Rodrigues, N. Mahant, P. A. Robinson, A unifying explanation of primary generalized seizures through nonlinear brain modeling and bifurcation analysis, Cereb. Cortex, 16 (2006), 1296–1313. https://doi.org/10.1093/cercor/bhj072 doi: 10.1093/cercor/bhj072
![]() |
[24] |
M. Cetin, Model-based robust suppression of epileptic seizures without sensory measurements, Cognit. Neurodyn., 14 (2020), 51–67. https://doi.org/10.1007/s11571-019-09555-8 doi: 10.1007/s11571-019-09555-8
![]() |
[25] |
D. G. Fan, L. Y. Zhang, Q. Y. Wang, Transition dynamics and adaptive synchronization of time-delay interconnected corticothalamic systems via nonlinear control, Nonlinear Dyn., 94 (2018), 2807–2825. https://doi.org/10.1007/s11071-018-4526-1 doi: 10.1007/s11071-018-4526-1
![]() |
[26] |
M. A. Kramer, B. A. Lopour, H. E. Kirsch, A. J. Szeri, Bifurcation control of a seizing human cortex, Phys. Rev. E: Stat. Nonlinear Biol. Soft Matter Phys., 73 (2006), 041928. https://doi.org/10.1103/PhysRevE.73.041928 doi: 10.1103/PhysRevE.73.041928
![]() |
[27] | P. N. Taylor, J. Thomas, N. Sinha, J. Dauwels, M. Kaiser, T. Thesen, et al., Optimal control based seizure abatement using patient derived connectivity, Front. Neurosci., 9 (2015). https://doi.org/10.3389/fnins.2015.00202 |
[28] |
O. Zakary, M. Rachik, I. Elmouki, On the analysis of a multi-regions discrete SIR epidemic model: an optimal control approach, Int. J. Dyn. Control, 5 (2017), 917–930. https://doi.org/10.1007/s40435-016-0233-2 doi: 10.1007/s40435-016-0233-2
![]() |
[29] |
E. H. Essoufi, A. Zafrar, Boundary optimal control of time-space SIR model with nonlinear Robin boundary condition, Int. J. Dyn. Control, 10 (2022), 1279–1290. https://doi.org/10.1007/s40435-021-00886-1 doi: 10.1007/s40435-021-00886-1
![]() |
[30] |
S. Mollah, S. Biswas, Optimal control for the complication of Type 2 diabetes: the role of awareness programs by media and treatment, Int. J. Dynam. Control, 11 (2023), 877–891. https://doi.org/10.1007/s40435-022-01013-4 doi: 10.1007/s40435-022-01013-4
![]() |
[31] |
N. Birla, A. Swarup, Optimal preview control: a review, Optim. Control. Appl. Methods, 36 (2015), 241–268. https://doi.org/10.1002/oca.2106 doi: 10.1002/oca.2106
![]() |
[32] |
X. Yu, F. C. Liao, L. Li, New results on observer-based robust preview tracking control for Lipschitz nonlinear systems, J. Vib. Control, 27 (2020), 2081–2096. https://doi.org/10.1177/1077546320953650 doi: 10.1177/1077546320953650
![]() |
[33] |
A. Khalil, N. Fezans, Gust load alleviation for flexible aircraft using discrete-time preview control, Aeronaut. J., 125 (2020), 341–364. https://doi.org/10.1017/aer.2020.85 doi: 10.1017/aer.2020.85
![]() |
[34] |
L. Li, F. C. Liao, Design of a preview controller for discrete-time systems based on LMI, Math. Probl. Eng., 2015 (2015), 179126. https://doi.org/10.1155/2015/179126 doi: 10.1155/2015/179126
![]() |
[35] |
E. Sitnikova, Thalamo-cortical mechanisms of sleep spindles and spike–wave discharges in rat model of absence epilepsy (a review), Epilepsy Res., 89 (2010), 17–26. https://doi.org/10.1016/j.eplepsyres.2009.09.005 doi: 10.1016/j.eplepsyres.2009.09.005
![]() |
[36] |
B. M. Bouwman, P. Suffczynski, F. H. L. da Silva, E. Maris, C. V. van Rijn, GABAergic mechanisms in absence epilepsy: a computational model of absence epilepsy simulating spike and wave discharges after vigabatrin in WAG/Rij rats, Eur. J. Neurosci., 25 (2007), 2783–2790. https://doi.org/10.1111/j.1460-9568.2007.05533.x doi: 10.1111/j.1460-9568.2007.05533.x
![]() |
[37] |
J. C. Young, A. Paolini, M. Pedersen, G. Jackson, Genetic absence epilepsy: effective connectivity from piriform cortex to mediodorsal thalamus, Epilepsy Behav., 97 (2019), 219–228. https://doi.org/10.1016/j.yebeh.2019.05.042 doi: 10.1016/j.yebeh.2019.05.042
![]() |
[38] |
F. Marten, S. Rodrigues, P. Suffczynski, M. Richardson, J. R. Terry, Derivation and analysis of an ordinary differential equation mean-field model for studying clinically recorded epilepsy dynamics, Phys. Rev. E: Stat. Nonlinear Biol. Soft Matter Phys., 79 (2009), 021911. https://doi.org/10.1103/PHYSREVE.79.021911 doi: 10.1103/PHYSREVE.79.021911
![]() |
[39] |
Y. J. Wang, M. Goodfellow, P. N. Taylor, G. Baier, Phase space approach for modeling of epileptic dynamics, Phys. Rev. E: Stat. Nonlinear Biol. Soft Matter Phys., 85 (2012), 061918. https://doi.org/10.1103/PHYSREVE.85.061918 doi: 10.1103/PHYSREVE.85.061918
![]() |
[40] |
P. N. Taylor, M. Goodfellow, Y. J. Wang, G. Baier, Towards a large-scale model of patient-specific epileptic spike-wave discharges, Biol. Cybern., 107 (2013), 83–94. https://doi.org/10.1007/s00422-012-0534-2 doi: 10.1007/s00422-012-0534-2
![]() |
[41] | X. Yu, Preview Control for Several Classes of Nonlinear Systems (in Chinese), Ph.D thesis, University of Science and Technology Beijing, 2019. |
[42] | T. E. T. Tsuchiya, F. C. Liao, The Latest Automatic Control Technology - Digital Preview and Predictive Control, Beijing Science and Technology Press, Beijing, 1994. |