Next Article in Journal
Semi-Analytical Solutions for Some Types of Nonlinear Fractional-Order Differential Equations Based on Third-Kind Chebyshev Polynomials
Next Article in Special Issue
Non-Overlapping Domain Decomposition for 1D Optimal Control Problems Governed by Time-Fractional Diffusion Equations on Coupled Domains: Optimality System and Virtual Controls
Previous Article in Journal
On Coefficient Inequalities of Starlike Functions Related to the q-Analog of Cosine Functions Defined by the Fractional q-Differential Operator
Previous Article in Special Issue
The Improved Stochastic Fractional Order Gradient Descent Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Control for Neutral Stochastic Integrodifferential Equations with Infinite Delay Driven by Poisson Jumps and Rosenblatt Process

by
Dimplekumar Chalishajar
1,*,
Ramkumar Kasinathan
2 and
Ravikumar Kasinathan
2
1
Department of Applied Mathematics, Virginia Military Institute (VMI), Lexington, VA 24450, USA
2
Department of Mathematics, PSG College of Arts & Science, Coimbatore 641014, India
*
Author to whom correspondence should be addressed.
Submission received: 21 July 2023 / Revised: 9 October 2023 / Accepted: 20 October 2023 / Published: 26 October 2023
(This article belongs to the Special Issue Optimal Control Problems for Fractional Differential Equations)

Abstract

:
In this paper, we investigate the optimal control problems for a class of neutral stochastic integrodifferential equations (NSIDEs) with infinite delay driven by Poisson jumps and the Rosenblat process in Hilbert space involving concrete-fading memory-phase space, in which we define the advanced phase space for infinite delay for the stochastic process. First, we introduce conditions that ensure the existence and uniqueness of mild solutions using stochastic analysis theory, successive approximation, and Grimmer’s resolvent operator theory. Next, we prove exponential stability, which includes mean square exponential stability, and this especially includes the exponential stability of solutions and their maps. Following that, we discuss the existence requirements of an optimal pair of systems governed by stochastic partial integrodifferential equations with infinite delay. Then, we explore examples that illustrate the potential of the main result, mainly in the heat equation, filter system, traffic signal light systems, and the biological processes in the human body. We conclude with a numerical simulation of the system studied. This work is a unique combination of the theory with practical examples and a numerical simulation.

1. Introduction

Since arbitrary fluctuations and noise are common and predictable in both natural and artificial systems, stochastic models should be studied rather than deterministic ones. Stochastic differential equations incorporate uncertainty into the mathematical representation of a specific occurrence of stochastic differential equations (SDEs). Due to their use in modeling numerous phenomena in physics, population dynamics, ecology, medicine, biology, engineering, and other areas of research, SDEs in both infinite and finite dimensions have garnered a lot of attention in recent years. For excellent information for learning about stochastic differential equations and their applications, see [1,2,3,4,5].
Since SDEs with delay are suitable for simulating physical systems with delay, such as those found in medicine, economics, epidemiology, biology, and chemistry, research on SDEs with delay has attracted a lot of attention over the past few decades; see [1,6,7,8] for a brief overview. Numerous writers have looked at the existence, uniqueness, controllability, and stability of the qualitative and quantitative characteristics of SDEs with delay solutions (see [9,10,11,12,13,14] and references therein). Researchers have been examining the existence and asymptotic behavior of some mild solutions to SDEs in Hilbert spaces using a semigroup approach [11], comparison theorem [15], Razumikhin-type theorem [16], analytic technique [17], Banach fixed point principle [18], etc.
In the literature, neutral stochastic partial differential equations (SPDEs) have been used to simulate a variety of dynamical systems [13,19,20]. Some of these equations have delayed state derivatives, which are different from SPDEs with delays that rely on the past and the present states. We merely recommend to see more information on this theory and its applications [14,21,22]. Recently, SIDEs have drawn a lot of interest with an emphasis on qualitative characteristics like optimality conditions, regularity, control issues, periodicity, and stability; we refer readers to [10,12,18]. Instead of using the continuous semigroup operators, the resolvent operator theory has been implemented as the system contains an integral term (see [23] for further details). However, the majority of academics that study exponential stability have focused only on infinite delays; see [9,18]. The solutions that are continuously dependent on the initial value have been the focus of this study. Authors in [21,22,24] discussed asymptotic and exponential stability of stochastic systems with infinite delay using the phase space given by Hale and Kato [25]. Notable is the paucity of research on the subject of defining the stochastic time delay systems involving exponential stability ([13,14,20,26]).
On the other hand, we consider the stationary Gaussian sequence ξ n n Z with variance and mean zero whose correlation function is R ( n ) : = E ( ξ 0 ξ n ) = n 2 H 2 k L ( n ) , with H 1 2 , 1 . Let f be a Hermite polynomial with rank k Let L be a slowly varying function at infinity. Also, if f admits the following expansion in Hermite polynomials
f ( x ) = j 0 c j H j ( x ) , c j = 1 j ! E ( f ( ξ 0 H ( ξ 0 ) ) ) ,
where k = min j | c j 0 1 and H j ( x ) represent the Hermite polynomial of degree j defined as H j ( x ) = ( 1 ) j e x 2 2 d j d x j e x 2 2 , then by virtue of the Non-Central Limit Theorem, this series 1 n H [ n t ] j = 1 g ( ξ j ) converges in the form of distribution in finite-dimension as n .
Z H k ( t ) = c ( H , k ) R k 0 t j = 1 k ( s y j ) + 1 2 + 1 H k d s d B ( y 1 ) d B ( y k ) .
Here, c ( H , k ) is a positive normalization constant that depends on H and k , and the multiple integral in Equation (1) is of Wiener–It o ^ -type integral of order k with respect to the standard Brownian motion B ( y k ) k R . For c > 0 , this process ( Z H k ( t ) ) t 0 is called the Hermite process and is self-similar to the Z H k ( c t ) = c H Z H k ( t ) and it has stationary increments [27].
When k = 1 , the Hermite process (1) is the fractional Brownian motion (fBm) with Hurst parameter H 1 2 , 1 [28]. If k = 2 , then the process (1) is called the Rosenblatt process, which occurs from the Non-Central Limit Theorem (see [29,30] and references therein). Self-similar processes with long-range dependence are important in practice and are employed in many areas, including finance, turbulence, internal traffic, hydrology, and econometrics, to mention a few [31,32]. The self-similar Rosenblatt process, which has stationary increments, seems to be the upper bound for long-range dependent stationary series, which is not a Gaussian process. Due to its significance in numerous fields of science and engineering, the fBm is undoubtedly the Hermite process in the class that has received the most research [12,33,34]. However, one can apply the Rosenblatt process in actual circumstances when the model’s assumption of Gaussianity is implausible. However, in contrast to the fBm, the study of the Rosenblatt process has received little attention due to its non-Gaussianity trait and complexity of its dependence structures. The study of SDEs using the Rosenblatt technique looks interesting. Because a subset of the class of self-decomposition distributions following the multivariate Rosenblatt distribution belongs to the Thorin class, it is possible to take advantage of this characteristic. Maejima and Tudor [27] showed that further properties of the Rosenblatt distribution can be discovered. Ramkumar et al. [35] recently used semigroup theory and the successive approximation method to study higher-order neutral SDEs driven by the Rosenblatt process in Hilbert space and Poisson jumps. One might consult the papers [12,14,34,36] and their references for more information on the Rosenblatt procedure. In system engineering, the optimal control issues arise frequently. The primary objective of optimum control is to identify the ideal values of the control variables for the dynamic system that maximize or minimize a specified performance index in an open-loop control. Because dynamic systems are nonlinear, finding the optimal control is a challenging and open-ended endeavor. In biomedicine, optimal control is frequently used to simulate cancer chemotherapy, and it has also recently been used in medicine and epidemiological models [37,38]. Rajivganthi and Muthukumar [39] examined the properties of Poisson jumps with fractional stochastic evolution equations and virtually automorphic solutions, as well as their optimal control. Wang et al. [40] studied the best possible control that is achieved for delayed SDEs powered by fBms. Very recently, Ramkumar et al. [41] investigated the existence of mild solutions and optimal controls for a class of fractional neutral SDEs driven by fBm and Poisson jumps in Hilbert spaces. Only a few authors have studied the optimal control results of stochastic differential/integrodifferential equations in infinite dimensional spaces (see [42,43] and references therein). Section 7 is devoted to the practical applications. We discuss the stochastic heat equation, the filter system, the traffic signal light systems, and the biological reaction on a human body to justify the practical applications of the theory studied.
More precisely, Boufoussi et al. [44] studied the global existence and uniqueness result of the mild solution for stochastic functional differential equations in Hilbert space driven by a fractional Brownian motion of the form
d x ( t ) = A x ( t ) + f ( x t ) d t + σ ( x t ) d B H ( t ) , t 0 , x ( t ) = η ( t ) , t [ r , 0 ] .
Next, Hajji et al. [45] extended the study of [44] and concerned themselves with the existence and uniqueness of mild solutions for a class of neutral functional stochastic differential equations described in the form
d x ( t ) + g ( t , x ( ρ ( t ) ) ) = A x ( t ) + f ( t , x ( ρ ( t ) ) ) d t + σ ( t ) d B Q H ( t ) , t [ 0 , T ] , t = φ ( t ) , r t 0 .
We extended the system studied by Hajji et al. [45] with numerous applications. Also, previous authors have not discussed the numerical simulation to justify their theory.
Considering the following neutral SIDEs with infinite delay driven by Poisson jumps and the Rosenblatt process given on the entire probability space ( Ω , , P ) is motivated by the explanation above:
d x ( t ) + f 1 ( t , x t ) = A x ( t ) + f 1 ( t , x t ) d t + f ( t , x t ) d t + 0 t Υ ( t s ) x ( s ) + f 1 ( s , x s ) d s d t + σ ( t , x t ) d Z H ( t ) + Z h ( t , x t , z ) N ˜ ( d t , d z ) , t 0 . x 0 ( · ) = ϕ C ζ ;
where ϕ is 0 -measurable, and the next section provides more information on the concrete fading memory-phase space C ζ . Υ ( t ) is given as a closed linear operator. Also, a closed linear operator A : D ( A ) H H is defined on a separable Hilbert space H such that D ( A ) D ( Υ ( t ) ) . Assume that B is another real separable Hilbert space. f 1 , f : [ 0 , + ) × C ζ H , σ : [ 0 , + ) × C ζ L 2 0 ( B , H ) , h : [ 0 , + ) × C ζ × Z H are appropriate functions; we see that function x t : { ( , 0 ] H , t 0 , such that x t ( θ ) = x ( t + θ ) } takes the value from the phase space C ζ . Let ( Ω , , P ) be the complete probability space. Assume that the Poisson counting measure N ( d t , d z ) induced by the Poisson point process p ( · ) in the measurable space ( Z , B ( Z ) ) . N ˜ ( d t , d z ) = N ( d t , d z ) ϕ ( d z ) d t is the compensated martingale measure. Z H ( t ) is a B -valued Rosenblatt process with parameter H 1 2 , 1 . The space of all continuous functions C ( ( , 0 ] , H ) is endowed with the norm given as
Φ = sup < θ 0 Φ ( θ ) .
We present the abstract phase space C ζ . Assume that ζ t : ( , 0 ] ( 0 , ) is a continuous function with l = 0 ζ ( s ) d s < . Assume the fading memory C ζ for ζ > 0 , defined as C ζ : = { Φ C ( ( , 0 ] , H ) : ( E Φ 2 ) 1 / 2 , is a bounded and measurable function on [ δ , 0 ] for any δ > 0 and 0 ζ ( t ) sup t δ 0 ( E Φ ( s ) 2 ) 1 / 2 d t < ; lim δ e ζ δ Φ ( δ ) exists in H } .
This space ( C ζ , · ζ ) represents a Banach space with the norm
Φ C ζ = 0 ζ ( s ) sup < δ 0 Φ ( δ ) d s , Φ C ζ .
The following lemma of phase spaces is essential for the infinite delay system.
Lemma 1.
Suppose z C ζ , then t [ 0 , T ] , z t C ζ . Moreover,
l ( E ( z ( t ) 2 ) ) 1 2 l sup 0 s t ( E z ( s ) 2 ) 1 2 + z 0 C ζ , z t C ζ ,
where l = 0 ζ ( s ) d s < .
Proof. 
It is simple to see that z t is bounded and measurable on [ a , 0 ] for a > 0 for any t [ 0 , a ] ,
z t C ζ = 0 ζ ( s ) sup δ [ s , 0 ] E | | z t ( δ ) | | d s = t ζ ( s ) sup δ [ s , 0 ] E | | z ( t + δ ) | | d s + t 0 ζ ( s ) sup δ [ s , 0 ] E | | z ( t + δ ) | | d s = t ζ ( s ) sup δ 1 [ t + s , t ] E | | z ( δ 1 ) | | d s + t 0 ζ ( s ) sup δ 1 [ t + s , t ] E | | z ( δ 1 ) | | d s t ζ ( s ) sup δ 1 [ t + s , 0 ] E | | z ( δ 1 ) | | + sup δ 1 [ 0 , t ] ( E | | z ( δ 1 ) | | 2 ) 1 / 2 d s + t 0 ζ ( s ) sup δ 1 [ 0 , t ] ( E | | z ( δ 1 ) | | 2 ) 1 / 2 d s = t ζ ( s ) sup δ 1 [ t + s , 0 ] E | | z ( δ 1 ) | | d s + 0 b ( s ) d s . sup s [ 0 , t ] ( E | | z ( s ) | | 2 ) 1 / 2 t ζ ( s ) sup δ 1 [ s , 0 ] E | | z ( δ 1 ) | | d s + l sup s [ 0 , t ] ( E | | z ( s ) | | 2 ) 1 / 2 0 ζ ( s ) sup δ 1 [ s , 0 ] E | | z ( δ 1 ) | | d s + l sup s [ 0 , t ] ( E | | z ( s ) | | 2 ) 1 / 2 = 0 ζ ( s ) sup δ 1 [ s , 0 ] | | z 0 ( δ 1 ) | | d s + l sup s [ 0 , t ] ( E | | z ( s ) | | 2 ) 1 / 2 = l sup s [ 0 , t ] ( E | | z ( s ) | | 2 ) 1 / 2 + z 0 C ζ .
Since Φ C ζ , then z t C ζ ,
z t C ζ = 0 ζ ( s ) sup δ [ s , 0 ] | | z t ( δ ) | | d s | | z t ( δ ) | | 0 ζ ( s ) d s E | | z ( t ) | | .
The proof is complete. □
The phase space C ζ defined above satisfies all the conditions of phase space given by Hall and Kato [25].
According to the author’s understanding, this research is the first to examine neutral SIDEs with infinite latency driven by Poisson jumps and the Rosenblatt process, as well as their existence and exponential stability. This paper’s key contribution is to identify prerequisites that guarantee:
(i)
The existence, uniqueness, and exponential stability of the solutions and their maps of system 2, including mean square exponential stability and nearly certain exponential stability, are explored.
(ii)
We demonstrate the outcome via stochastic methods and Grimmer’s [23] resolvent operator theory. It is important to note that system 2 was explored with a finite delay by Diop et al. [18]. They studied the existence of mild solutions with exponential stability. Due to this, our method can be viewed as an extension of the conclusion of [18] for the case of infinite delay.
(iii)
The literature does not address the best controls for neutral SIDEs with infinite delay of system (35), and this fact drives us to create the current work on this subject by both extending the existing controls and creating new ones.
(iv)
Lastly, to demonstrate the established idea, we came up with three alternative cases.
This manuscript is organized as follows: In Section 2 There are some notations and preliminary remarks. Section 3 represents the existence and uniqueness of mild solutions for neutral SIDEs with infinite delay. Section 4 ensures the conditions assuring mean square moment exponential stability of the solution z ( t ) and almost surely exponential stability of the solution z ( t ) is studied in Section 5. Section 6 displays the outcomes for optimal pairs of systems governed by the stochastic control system (35). Section 7 contains four examples to illustrate the theory. To illustrate the efficiency of the theoretical result, we study the numerical simulation in Section 8.

2. Preliminary Remarks

Let H and K be two real separable Hilbert spaces and represent the norm operator in H , K . The space L ( K , H ) from K to H consists of bounded linear operators endowed with · , and system 2 is attired with a normal filtration t t 0 .
Rosenblatt Process: The basic properties of the Rosenblatt process and Wiener integral are given below. Let Z A ( t ) , t [ 0 , T ] be a one-dimensional Rosenblatt process with parameter A 1 2 , 1 on [ 0 , T ] , with T being an arbitrary fixed horizon representing the following equation [30]:
Z A ( t ) = d ( A ) 0 t 0 t y 1 y 2 t B A x ( x , y 1 ) B A x ( x , y 2 ) d B ( y 1 ) d B ( y 2 ) , for A > 1 2 ,
with B A ( t , s ) given by
B A = c A s 1 2 A s t ( x s ) A 3 2 x A 1 2 d x , for t > s ,
where c A = A ( 2 A 1 ) β ( 2 2 A , A 1 2 ) , B A ( t , s ) = 0 when t s , B ( t ) , t [ 0 , T ] is a standard Brownian motion, A = A + 1 2 and d ( A ) = 1 A + 1 A 2 ( 2 A 1 ) is a normalizing constant and β ( · , · ) is the Beta function. The covariance function of the Rosenblatt process Z A ( t ) , t [ 0 , T ] is
E Z A ( t ) Z A ( s ) = 1 2 s 2 A + t 2 A s t 2 A .
The Wiener integral can be constructed with respect to the Rosenblatt process due to its covariance structure. For more background information for this section, the reader is referred to [29,30].
Poisson Process: For t 0 , p ( t ) is a B -valued, σ -finite stationary, and the t -adapted Poisson point process has values in measurable space ( Z , B ) and an intensity measure of σ -finite v ( d z ) . We define N ( d s , d z ) as the Poisson counting measure. p ( · ) induces the Poisson counting measure and the compensatory martingale is provided by
N ˜ ( d s , d z ) = N ( d s , d z ) v ( d z ) d s .
Definition 1.
A resolvent operator R ( t ) L ( X ) , t 0 is a bounded linear operator-valued function for Equation (2) if R ( t ) attains the following:
(i) 
R ( 0 ) = I , R ( t ) N e α t t 0 , N 1 , α R .
For α > 0 , R ( t ) is exponentially stable.
(ii) 
The function t R ( t ) x is strongly continuous for each t 0 and for x X , x in Y , R ( · ) x C 1 ( [ 0 , + ) ; X ) C ( [ 0 , + ) ; Y ) and satisfies
d R ( t ) x = A R ( t ) x + 0 t Υ ( t s ) R ( s ) x d s d t = R ( t ) A x + 0 t R ( t s ) Υ ( s ) x d s d t .
The existence of solutions for Equation (3) is ensured by the following two criteria, which are adopted from Grimmer [23].
(A1)
The operator A is an infinitesimal generator of a C 0 -semigroup on X .
(A2)
For all t 0 , Υ ( t ) denotes a closed, continuous linear operator from D ( A ) to X , and Υ ( t ) belongs to L ( Y , X ) . For any y Y , the map t Υ ( t ) y is bounded, differentiable, and its derivative d Υ ( t ) y / d t is bounded and uniformly continuous on [ 0 , ) .
Remark 1.
The existence and uniqueness of R ( t ) is guaranteed by (A1) and (A2).
Consider the following deterministic integrodifferential equation:
d u ( t ) = A u ( t ) + 0 t Υ ( t s ) u ( s ) d s + a ( t ) d t , t 0 , u ( 0 ) = u 0 X
and a ( t ) : [ 0 , + ) X is a continuous function.
Lemma 2.
Let assumptions (A1) and (A2) hold good. If u is a strict solution of (3), i.e., u C 1 ( [ 0 , + ) ; X ) C ( [ 0 , + ) ; Y ) , then
u ( t ) = R ( t ) u 0 + 0 t R ( t s ) a ( s ) d s , t 0 .
Lemma 3.
If Ψ : [ 0 , T ] L 2 0 ( B , A ) satisfies 0 T Ψ ( s ) L 2 0 2 d s < then system (1) is well defined as an H -valued random variable, and
E 0 t Ψ ( s ) d Z A ( s ) 2 < C A t 2 A 1 0 t Ψ ( s ) L 2 0 2 d s .
Definition 2.
A mild solution of Equation (2) is an H -valued process z ( t ) , 0 t T satisfying the following:
(i) 
z ( t ) is t -adapted and 0 T z ( t ) 2 d t < almost surely.
(ii) 
z ( t ) is continuous for t [ 0 , T ] and satisfies
z ( t ) = R ( t ) ϕ ( 0 ) + g 1 ( 0 , ϕ ) g 1 ( t , z t ) + 0 t R ( t s ) g ( s , z s ) d s + 0 t R ( t s ) σ ( s , z s ) d Z A ( s ) + 0 t Z R ( t s ) h ( s , z s , z ) N ˜ ( d s , d z ) ,
with z 0 ( · ) = ϕ C ζ .
The following presumptions are useful to accomplish the main objective.
(H1)
R ( · ) is exponentially stable; that is, there exist two constants ϕ > 0 and M 1
R ( t ) M e ϕ t , t 0 .
(H2)
∃ a real number B 0 > 0 f 1 ( t , x ) f 1 ( t , y ) 2 B 0 x y C ζ 2 ; x , y C ζ , t 0 .
(H3)
∃ a real number B 1 > 0 , such that
( i ) f ( t , x ) f ( t , y ) 2 σ ( t , x ) σ ( t , y ) 2 B 1 x y C ζ 2 , ( i i ) Z h ( t , x , z ) h ( t , y , z ) 2 v ( d z ) d t Z h ( t , x , z ) h ( t , y , z ) 4 v ( d z ) d t 1 2 B 1 x y C ζ 2 , x , y C ζ , t 0 . ( i i i ) Z h ( t , x , z ) h ( t , y , z ) 2 v ( d z ) d t 1 2 B 1 x y C ζ 2 , x , y C ζ .
Obviously, f 1 ( t , 0 ) f ( t , 0 ) σ ( t , 0 ) h ( t , 0 , z ) 0 .
Remark 2.
To ensure that the stochastic system has a zero equilibrium solution, the condition f 1 ( t , 0 ) = f ( t , 0 ) = σ ( t , 0 ) = h ( t , 0 , z ) = 0 is taken into consideration (2). If this supposition is incorrect, it is always possible to convert the equilibrium solution for (2) into the zero equilibrium solution of another system.

3. Existence and Uniqueness

First, we ensure sufficiency for the existence and uniqueness of the solution for (2). The sequential approximation method and numerous stochastic analytic methods are used to accomplish this. However, in order to deal with infinite delay, we must create some novel strategies. In the sections that follow, we replace X for the Hilbert space H in (A1) and (A2). We now offer the primary finding listed below.
Theorem 1.
Let assumptions (A1), (A2), (H2), and (H3) hold good with B 0 < 1 10 . Then, system (2) has a unique mild solution.
Proof. 
We use the successive approximation technique to Equation (5).
For t 0 , x n ( t ) = ϕ ( t ) , n N , and for 0 t T ,
x n ( t ) = R ( t ) ϕ ( 0 ) + f 1 ( 0 , ϕ ) f 1 ( t , x t n ) + 0 t R ( t s ) f ( s , x s n 1 ) d s + 0 t R ( t s ) σ ( s , x s n 1 ) d Z A ( s ) + 0 t Z R ( t s ) h ( s , x s n 1 , z ) N ˜ ( d s , d z ) ;
and
x 0 ( t ) = R ( t ) ϕ ( 0 ) , 0 t T , n 1 .
Let
M T = sup 0 t T R ( t ) L ( A ) .
i.e., M T is uniformly bounded.
The proof now follows the steps below:
Step 1: { x n ( t ) } , n 0 is bounded.
Using (6)
E x n ( t ) 2 = E ϕ ( t ) C ζ 2 < , t ( , 0 ] ,
we obtain that
E sup 0 s t x 0 ( t ) M T 2 E ϕ C ζ 2 < , t [ 0 , T ] .
Then, for any n 1 , we have,
E x n ( t ) 2 = E R ( t ) ϕ ( 0 ) + f 1 ( 0 , ϕ ) f 1 ( t , x t n ) + 0 t R ( t s ) f ( s , x s n 1 ) d s + 0 t R ( t s ) σ ( s , x s n 1 ) d Z A ( s ) + 0 t Z R ( t s ) h ( s , x s n 1 , z ) N ˜ ( d s , d z ) 2 E x n ( t ) 2 = 5 i = 1 5 Δ i .
From assumption (H2), we obtain
Δ 1 = E sup 0 s t R ( t ) ϕ ( 0 ) + f 1 ( 0 , ϕ ) 2 , 2 M T 2 1 + B 0 E ϕ C ζ 2 ,
and
Δ 2 = E sup 0 s t f 1 ( s , x s n ) 2 , B 0 E sup 0 s t x s n C ζ 2 .
Next, combining (H3) and Hölder’s inequality yields
Δ 3 = M T 2 E sup 0 s t 0 t f ( s , x s n 1 ) d s 2 , M T 2 B 1 0 t E sup 0 s t x s n 1 C ζ 2 .
Now, using Lemma 2, and (H3), we obtain
Δ 4 = M T 2 E sup 0 s t 0 t σ ( s , x s n 1 ) d Z A ( s ) 2 , M T 2 C A t 2 A 0 t E sup 0 s t x s n 1 2 d s .
Finally, with (H3), Lemma 2, and Hölder’s inequality, we obtain
Δ 5 = M T 2 { E 0 t Z h ( s , x s n 1 , z ) 2 v ( d z ) d s + E 0 t Z h ( s , x s n 1 , z ) 4 v ( d z ) d s 1 2 } 3 M T 2 B 1 T 0 t E sup 0 s t x s n 1 2 d s .
Substituting (8) and (9) into (7), we obtain
E sup 0 s t x n ( s ) 2 10 M T 2 1 + B 0 E ϕ 2 + 5 B 0 E sup 0 s t x s n 2 + Q 1 0 t E sup 0 s t x s n 1 2 d s ,
where Q 1 = 5 M T 2 B 1 2 t + 2 C A t 2 A + 3 .
By using the definition of the norm · C ζ , we can write
E sup 0 s t x n ( s ) 2 Q 2 + Q 1 1 5 B 0 0 t E sup 0 s t x n 1 ( s ) 2 d s ,
where Q 2 = 10 M T 2 1 + B 0 + 5 B 0 + Q 1 T .
Therefore,
max 1 n k E sup 0 s t x n ( s ) 2 Q 2 + Q 3 0 t max 1 n k E sup 0 s t x n 1 ( s ) 2 d s , k 1 ,
where Q 3 = Q 1 1 5 B 0 . Besides that, we have
max 1 n k E x n 1 ( s ) 2 E x 0 ( s ) 2 + max 1 n k E x n ( s ) 2 M T 2 E ϕ C ζ 2 + max 1 n k E x n ( s ) 2 .
It follows that
max 1 n k E sup 0 s t x n ( s ) 2 Q 2 + Q 3 M T 2 T E ϕ C ζ 2 + Q 3 0 t max 1 n k E sup 0 s t x n ( s ) 2 d s , Q 4 + Q 3 0 t max 1 n k E sup 0 s t x n ( s ) 2 d s ,
where Q 4 = Q 2 + Q 3 M T 2 T E ϕ C ζ 2 .
Using Gronwall’s inequality,
max 1 n k E sup 0 s t x n ( s ) 2 Q 4 e Q 3 t < .
Since k is arbitrary, we have
E x n ( t ) 2 Q 4 e Q 3 t < ,
which completes the proof.
Step 2: x n , n N is a Cauchy sequence.
The successive approximations are constructed to have x n ( t ) = x n 1 on ( , 0 ] for n 1 . So we now establish
x t 0 C ζ 2 1 + M T 2 ϕ C ζ 2 , t [ 0 , T ] .
Therefore from (3.1), we have
E x t 1 x t 0 2 5 i = 1 5 Ω i .
By using (H2), we obtain
Ω 1 E R ( t ) f 1 ( 0 , ϕ ) 2 , M T 2 B 0 E ϕ C ζ 2 .
and
Ω 2 E f 1 ( t , x t 1 ) 2 B 0 E x t 1 C ζ 2 2 B 0 E x t 1 x t 0 C ζ 2 + E x t 0 C ζ 2 2 B 0 E x t 1 x t 0 C ζ 2 + 2 B 0 1 + M T 2 E ϕ C ζ 2 .
Next, using (H3) and Hölder’s inequality, we obtain
Ω 3 M T 2 E 0 t f ( s , x s 0 ) d s 2 M T 2 T 2 B 1 E 0 t x s 0 C ζ 2 d s M T 2 T 2 B 1 1 + M T 2 E ϕ C ζ 2 d s .
From Lemma 2, we have
Ω 4 M T 2 T E 0 t σ ( s , x s 0 ) d Z ( s ) M T 2 T C A t 2 A 0 T E x s 0 C ζ 2 d s M T 2 T C A t 2 A 1 + M T 2 E ϕ C ζ 2 d s .
In a similar way, using (H3), we obtain
Ω 5 = M T 2 { E 0 t Z h ( s , x s 0 , z ) 2 v ( d z ) d s + E 0 t Z h ( s , x s 0 , z ) 4 v ( d z ) d s 1 2 } 3 M T 2 B 1 T 1 + M T 2 E ϕ C ζ 2 d s .
Substituting (17)–(21) into (16), we obtain
E x 1 ( t ) x 0 ( 0 ) 2 10 B 0 E x t 1 x t 0 C ζ 2 + 5 { M T 2 B 0 + 2 B 0 1 + M T 2 + M T 2 T 2 B 1 1 + M T 2 + M T 2 T C A t 2 A 1 + M T 2 + 3 M T 2 B 1 T 1 + M T 2 } E ϕ C ζ 2 Q 5 + 10 B 0 E x t 1 x t 0 C ζ 2 ;
where
Q 5 = 5 { M T 2 B 0 + 2 B 0 1 + M T 2 + M T 2 T 2 B 1 1 + M T 2 + M T 2 T C A t 2 A 1 + M T 2 + 3 M T 2 B 1 T 1 + M T 2 } E ϕ C ζ 2 .
This yields
E x 1 ( t ) x 0 ( 0 ) 2 sup 0 s t x 1 ( t ) x 0 ( t ) 2 .
Next, we can deduce
E x 1 ( t ) x 0 ( 0 ) 2 Q 5 1 10 B 0 = Q 6 .
Similarly,
E x 2 ( t ) x 1 ( t ) 2 4 { E f 1 ( t , x t 2 ) f 1 ( t , x t 1 ) 2 + E 0 t R ( t s ) [ f ( s , x s 1 ) f ( s , x s 0 ) ] d s 2 + E 0 t R ( t s ) [ σ ( s , x s 1 ) σ ( s , x s 0 ) ] d Z A ( s ) 2 + E 0 t Z R ( t s ) [ h ( s , x s 1 , z ) h ( s , x s 0 , z ) ] N ˜ ( d s , d z ) 2 } 4 B 0 E x t 2 x t 1 C ζ 2 + 4 M T 2 B 1 T + C A t 2 A + 3 T 0 t E x s 1 x s 0 C ζ 2 d s 4 B 0 E x t 2 x t 1 C ζ 2 + 4 M T 2 B 1 T + C A t 2 A + 3 T t Q 6 .
Therefore,
E x 2 ( s ) x 1 ( s ) 2 4 M T 2 B 1 T + C A t 2 A + 3 T t 1 4 B 0 Q 6 , = Q 6 t Q 7 ,
where Q 7 = 4 M T 2 B 1 T + C A t 2 A + 3 T 1 4 B 0 .
Sequentially, we can prove that
E sup 0 s t x 3 ( s ) x 2 ( s ) 2 ( Q 7 t ) 2 2 ! Q 6 .
Using the successive iteration scheme, n 0 ,
E sup 0 s t x n + 1 ( s ) x n ( s ) 2 ( Q 7 t ) n n ! Q 6 .
Therefore,
E x m ( t ) x n ( t ) 2 Q 6 k = n m 1 ( Q 7 t ) k k ! 0 as n , m > n 0 .
This argument proves that x n ( t ) , n 0 is a Cauchy sequence in L 2 ( Ω , H ) .
Step 3: Existence and uniqueness of Equation (2).
We have x n ( t ) x ( t ) in L 2 as n . For all t [ 0 , T ] , the Borel–Cantelli lemma shows that x n ( t ) converge to x ( t ) uniformly as n , for t ( , T ] . We prove the next inequality using (H2) and (H3).
E f 1 ( t , x t n ) f 1 ( t , x t ) 2 B 0 E x t n x t C ζ 2 B 0 E sup 0 s t x s n x s 2 0 , as n .
Next,
E 0 t R ( t s ) [ f ( s , x s n ) f ( s , x s ) ] d s 2 M T 2 T B 1 0 t E x s n x s C ζ 2 d s 0 , as n .
Similarly, we have
E 0 t R ( t s ) [ σ ( s , x s n ) σ ( s , x s ) ] d Z ( s ) 2 M T 2 T C A t 2 A 0 t E x s n x s C ζ 2 d s 0 , as n ;
and
E 0 t R ( t s ) [ h ( s , x s n , z ) h ( s , x s , z ) ] N ˜ ( d s , d z ) 2 3 M T 2 B 1 T 0 t E x s n x s C ζ 2 d s 0 , as n .
Taking the limit as n yields,
x ( t ) = R ( t ) ϕ ( 0 ) + f 1 ( 0 , ϕ ) f 1 ( t , x t ) + 0 t R ( t s ) f ( s , x s ) d s + 0 t R ( t s ) σ ( s , x s ) d Z A ( s ) + 0 t Z R ( t s ) h ( s , x s , z ) N ˜ ( d s , d z ) .
Gronwall’s lemma assures the assertion which completes the proof. □
Remark 3.
The existence and uniqueness of the local solution is defined on ( , T ] , and that leads us to the global solution on ( , + ) .

4. Exponential Stability

Definition 3.
For any initial value ϕ C ζ , the mild solution of 2 is said to be mean square mment exponentially stable if it is 0 -measurable and two positive real numbers α > 0 and β > 0 such that
E x ( t ) 2 α E ϕ C ζ 2 exp ( β t ) , t 0 .
The following theorem contains the major finding of this work.
Theorem 2.
Let us assume that for all conditions of Theorem 1, the following inequalities hold:
2 μ > ϕ and Q 10 1 5 B 0 ϕ < 0 .
The mild solution x ( t ) of Equation (2) and the solution x t that maps to Equation (2) are then both mean square moment exponentially stable.
Proof. 
The mean square moment exponential stability of x ( t ) : Combining Lemma 3 and (H1)–(H3), we can write
E x ( t ) 2 5 { E R ( t ) ϕ ( 0 ) + f 1 ( 0 , ϕ ) 2 + E f 1 ( t , x t ) 2 + E 0 t R ( t s ) f ( s , x s ) d s 2 + E 0 t R ( t s ) σ ( s , x s ) d Z A ( s ) 2 + E 0 t Z R ( t s ) h ( s , x s , z ) N ˜ ( d s , d z ) 2 } = 5 i = 1 5 Θ i .
Note that
Θ 1 = E R ( t ) ϕ ( 0 ) + f 1 ( 0 , ϕ ) 2 2 E R ( t ) ϕ ( 0 ) 2 + E R ( t ) f 1 ( 0 , ϕ ) 2 2 M 2 e ϕ t 1 + B 0 E ϕ C ζ 2 ,
and
Θ 2 = E f 1 ( t , x t ) 2 B 0 E x t C ζ 2 .
By Hölder’s inequality,
Θ 3 = E 0 t R ( t s ) f ( s , x s ) d s 2 M 2 0 t e 2 ϕ ( t s ) 0 t e ϕ ( t s ) E f 1 ( s , x s ) 2 d s M 2 0 t e ϕ ( t s ) d s 0 t e ϕ ( t s ) E f 1 ( s , x s ) 2 d s M 2 1 ϕ B 1 0 t e ϕ ( t s ) E x s C ζ 2 d s M 2 ϕ 1 B 1 0 t e ϕ ( t s ) E x s C ζ 2 d s .
Using Lemma 3 and (H3), we obtain
Θ 4 = E 0 t R ( t s ) σ ( s , x s ) d Z A ( s ) 2 M 2 C A t 2 A 1 ϕ B 1 0 t e ϕ ( t s ) E x s C ζ 2 d s M 2 C A t 2 A ( ϕ ) 1 B 1 0 t e ϕ ( t s ) E x s C ζ 2 d s .
Finally, using (H1) and (H3), we obtain
Θ 5 = E 0 t Z R ( t s ) h ( s , x s , z ) N ˜ ( d s , d z ) 2 M 2 { E 0 t e 2 ϕ ( t s ) Z h ( s , x ( s ) , z ) 2 v ( d z ) d s + E 0 t Z h ( s , x ( s ) , z ) 2 v ( d z ) d s 1 2 } 3 M 2 B 1 ( ϕ ) 1 0 t e ϕ ( t s ) E x s C ζ 2 d s .
Substituting (24)–(28) into (23), and using Definition 3, we obtain
E x ( t ) 2 10 M 2 e ϕ t 1 + B 0 E ϕ C ζ 2 + 5 B 0 E x s C ζ 2 + 5 M 2 ϕ 1 B 1 + M 2 C A t 2 A ( ϕ ) 1 B 1 + 5 M 2 B 1 ( ϕ ) 1 0 t e ϕ ( t s ) E x s C ζ 2 d s 10 M 2 e ϕ t 1 + B 0 E ϕ C ζ 2 + 5 B 0 E x s C ζ 2 + Q 8 0 t e ϕ ( t s ) E x s C ζ 2 d s ,
where Q 8 = 5 M 2 ϕ 1 B 1 + M 2 C A t 2 A ( ϕ ) 1 B 1 + 5 M 2 B 1 ( ϕ ) 1 .
Multiplying both sides of (29) by e ϕ t yields
e ϕ t E x ( t ) 2 10 M 2 1 + B 0 E ϕ C ζ 2 + 5 B 0 e ϕ t E x s C ζ 2 + Q 8 0 t e ϕ s E x s C ζ 2 d s .
For t > 0 , we have
e ϕ t E x ( t ) 2 10 M 2 1 + B 0 E ϕ C ζ 2 + 5 B 0 e ϕ t E x s C ζ 2 + Q 8 0 t e ϕ s E x s C ζ 2 d s ,
From Appendix A, pp. 12 of [46], we have the following estimation:
E e ϕ t x ( t ) 2 10 M 2 1 + B 0 E ϕ C ζ 2 + 5 B 0 e ( ϕ 2 μ ) t E ϕ C ζ 2 + E sup 0 s t e ϕ s x ( s ) 2 + Q 8 0 t e ( ϕ 2 μ ) s E ϕ C ζ 2 + E sup 0 s t e ϕ s x ( s ) 2 d s .
Recall that B 0 < 1 10 and ϕ 2 μ < 0 ; hence,
E sup 0 s t e ϕ s x ( s ) 2 1 1 5 B 0 10 M 2 1 + B 0 + 5 B 0 E ϕ C ζ 2 + Q 8 1 5 B 0 ( 2 μ ϕ ) 1 e ( ϕ 2 μ ) t E ϕ C ζ 2 + Q 8 1 5 B 0 0 t E sup 0 s t e ϕ s x ( s ) 2 d s = Q 9 + Q 10 0 t E sup 0 s t e ϕ s x ( s ) 2 d s ,
where
Q 9 = 10 M 2 ( 1 + B 0 ) + 5 B 0 + Q 8 2 μ ϕ 1 5 B 0 E ϕ C ζ 2 and Q 10 = Q 8 1 5 B 0 .
By virtue of Gronwall’s lemma,
E sup 0 s t e ϕ s x ( s ) 2 Q 9 e Q 10 t ,
which implies that
E x ( t ) 2 Q 9 e ( Q 10 ϕ ) t .
The result of the square moment exponential stability of solution x ( t ) is therefore satisfied given the criterion in (22). For any t 0 , we have
x t C ζ 2 e 2 μ t ϕ C ζ 2 + sup 0 s t x ( s ) 2 ,
multiplying both sides by e ϕ t ,
e ϕ t x t C ζ 2 e ( ϕ 2 μ ) t ϕ C ζ 2 + sup 0 s t e ϕ s x ( s ) 2 ϕ C ζ 2 + sup 0 s t e ϕ s x ( s ) 2 .
From (30), we have
E sup 0 s t e ϕ s x s C ζ 2 ϕ C ζ 2 + 10 M 2 1 + B 0 E ϕ C ζ 2 + 5 B 0 e ϕ t E x s C ζ 2 + Q 8 0 t e ϕ s E x s C ζ 2 d s ( 1 5 B 0 ) E sup 0 s t e ϕ s x s C ζ 2 1 + 10 M 2 1 + B 0 E ϕ C ζ 2 + Q 8 0 t e ϕ s E x s C ζ 2 d s Q 10 E ϕ C ζ 2 + Q 11 0 t E sup 0 s t e ϕ s x s C ζ 2 d s ,
where Q 10 = 1 + 10 M 2 1 + B 0 1 5 B 0 and Q 11 = Q 8 1 5 B 0 . Thus, we have
E sup 0 s t e ϕ s x s C ζ 2 Q 10 e ( Q 11 ϕ ) t .
The inequality in (31) ensures that the square moment of the solution map x t is exponentially stable.

5. Almost Surely Exponential Stability

Definition 4.
If the following inequality is almost certainly guaranteed, the mild solution of (2) is said to be exponentially stable if
0 > lim t sup 1 t l o g x ( t ) ϕ C ζ .
Theorem 3.
Suppose that in Theorem 2 the prerequisites are all true.
(i) 
lim t sup 1 t log x ( t ) ( Q 10 ϕ ) ϵ 2 almost surely,
(ii) 
lim t sup 1 t log x t ( Q 10 ϕ ) ϵ 2 almost surely,
which means that x ( t ) and x t are almost surely exponentially stable for any ϵ ( 0 , 1 ) .
Proof. 
Theorem 1 guarantees condition (i) of Theorem 3 n 0 .
E sup n t n + 1 x ( t ) 2 Q 9 e ( Q 10 ϕ ) ( n + 1 ) = Q 9 e e Q 10 ϕ . e ( Q 10 ϕ ) n .
For any ϵ ( 0 , 1 ) . Let I n be the interval [ n , n + 1 ] .
Q 10 = Q 10 1 5 B 0 .
Since 1 ϵ > 0 and Q 10 ϕ < 0 , by using the Markov inequality, we have
P sup t I n x ( t ) 2 > e ( Q 10 ϕ ) n ϵ E x ( t ) 2 e ( Q 10 ϕ ) n ϵ
Q 9 e Q 10 ϕ e ( Q 10 ϕ ) n ( 1 ϵ ) .
In Equation (33), the right-hand-side term is bounded by
n = 0 e ( 1 ϵ ) ( Q 10 ϕ ) n < .
Therefore, the Borel–Cantelli lemma guarantees that ∃ an integer n 0 such that n n 0 ,
sup t I n x ( t ) 2 e ( Q 10 ϕ ) n ϵ , almost surely .
Thus,
1 t log x ( t ) 2 1 n Q 10 ϕ n ϵ = ( Q 10 ϕ ) ϵ , t I n , n n 0 ,
yields,
lim t sup 1 t log x ( t ) Q 10 ϕ n ϵ 2 , almost surely ,
which shows that Theorem 3 (i) is satisfied. The argument that proves Theorem 3 (ii) follows an analogous reasoning. To be specific, we can use the assumption that solution map x t is square moment exponentially stable and reach the same conclusion by using the same logic as before. The details are left out. □

6. Optimal Control

Let control u be associated with the mild solution x u of (2). Let B be the reflexive Banach space where the values of control u are taken. B ( J , L ( B , A ) ) represents the set of operator valued functions that are uniformly bounded in J and measurable in the strong operator topology. Let L 2 ( J , B ) be the closed subspace of B containing all measurable and t -adapted, B -valued stochastic processes satisfying the condition 0 t u ( t ) B 2 d t < and endowed with the norm
u ( t ) L 2 ( J , B ) = 0 t u ( t ) B 2 d t 1 2 .
A nonempty closed bounded convex set v ( · ) is a subset of B .
Define A a d = u ( · ) L 2 ( J , B ) : u ( t ) A a d a . e , t J . We preemptively believe the following:
(H4)
B B ( J , L ( B , A ) ) .
Using Bolza problem P ˜ (see [39]), we can find an optimal pair ( x 0 , u 0 ) C ζ × A a d , such that l ( x 0 , u 0 ) l ( x u , u ) , u A a d where the cost functional
l ( x u , u ) = E 0 t J ( t , x t u , x u ( t ) , u ( t ) ) d t + E Θ ( x u ( T ) ) .
(H5)
We introduce the following hypotheses:
(1) The functional J : [ 0 , T ] × C ζ × A × B R is Borel measurable.
(2) For almost all t [ 0 , T ] , J ( t , · , · , · ) is sequentially lower semi-continuous on C ζ × A × B .
(3) For almost all t [ 0 , T ] , J ( t , x , x t , · ) is convex on B for each x C ζ .
(4) For j > 0 , μ 0 is non-negative, μ 0 L 1 ( [ 0 , T ] ; R ) ∋.
μ 0 ( t ) + d E x 2 + e E x t 2 + j E u B 2 J ( t , x ( t ) , x t , u ( t ) ) ;
where d , e 0 are constants.
Consider the optimal control for neutral SIDEs with infinite delay driven by Poisson jumps and the Rosenblatt process:
d x ( t ) + f 1 ( t , x t ) = A x ( t ) + f 1 ( t , x t ) d t + f ( t , x t ) d t + B ( t ) u ( t ) d t + 0 t Υ ( t s ) x ( s ) + f 1 ( s , x s ) d s + σ ( t , x t ) d Z A ( t ) + Z h ( t , x t , z ) N ˜ ( d t , d z ) , t J = ( 0 , T ] , t 0 , J = [ 0 , T ] . u 0 ( · ) = ϕ C ζ .
Theorem 4.
If (H1)–(H5) are satisfied, then there exists a unique mild solution of the system (35) for every u A a d of the form
x ( t ) = R ( t ) ϕ ( 0 ) + f 1 ( 0 , ϕ ) f 1 ( t , x t ) + 0 t R ( t s ) f ( s , x s ) d s + 0 t R ( t s ) B ( s ) u ( s ) d s + 0 t R ( t s ) σ ( s , x s ) d Z A ( s ) + 0 t Z R ( t s ) h ( s , x s , z ) N ˜ ( d s , d z ) , t J .
Proof. 
This theorem’s proof is comparable to that of Theorem 1, and since the successive approximation approach can be used to demonstrate the solution to system (35), it is removed. □
Theorem 5.
Assume that (H1)–(H5) are fulfilled. Also, Theorems 1 and 4 hold good. Then, the stochastic control problem (35) admits at least one optimal pair.
Proof. 
If inf l ( u ) : u A a d = + , then there is nothing to prove. Assume that inf l ( u ) : u A a d = ϵ < . Using ( H 5 ) , we have ϵ > . By definition of infimum, ∃ a minimizing sequence feasible pair ( x n , u n ) A a d { ( x , u ) : x n is a mild solution of system (35) corresponding to u n A a d } such that l ( x n , u n ) ϵ as n + . Since u n A a d , u n n 1 L 2 ( J , B ) is bounded, ∃ u 0 L 2 ( J , B ) and a subsequence is extracted from u n such that u n w u 0 in L 2 ( J , B ) . Since A a d is closed and convex, the Mazur lemma forces us to conclude that u 0 A a d . Suppose that x n and x 0 are the mild solutions of (35) corresponding to u n and u 0 , respectively; i.e.,
x n ( t ) = ϕ ( t ) , t [ r , 0 ] , R ( t ) ϕ ( 0 ) + f 1 ( 0 , ϕ ) f 1 ( t , x t ) + 0 t R ( t s ) f ( s , x s n ) d s + 0 t R ( t s ) B ( s ) u n ( s ) d s + 0 t R ( t s ) σ ( s , x s n ) d Z A ( t ) + 0 t Z R ( t s ) h ( s , x s n , z ) N ˜ ( d t , d z ) , t J .
Similarly, corresponding to u 0 , ∃ a mild solution x 0 of (35); that is,
x 0 ( t ) = ϕ ( t ) , t [ r , 0 ] , R ( t ) ϕ ( 0 ) + f 1 ( 0 , ϕ ) f 1 ( t , x t ) + 0 t R ( t s ) f ( s , x s 0 ) d s + 0 t R ( t s ) B ( s ) u 0 ( s ) d s + 0 t R ( t s ) σ ( s , x s 0 ) d Z A ( t ) + 0 t Z R ( t s ) h ( s , x s 0 , z ) N ˜ ( d t , d z ) , t J .
Using (H1)–(H5) and the Hölder inequality for t J , we have
E x n ( t ) x 0 ( t ) 2 5 B 0 E x n ( s ) x 0 ( s ) 2 + 5 M T 2 B 1 t 0 t E x n ( s ) x 0 ( s ) 2 d s + 5 M T 2 B 1 C A t 2 A 0 t E x n ( s ) x 0 ( s ) 2 d s + 10 M T 2 B 1 0 t E x n ( s ) x 0 ( s ) 2 d s + 5 M T 2 t B ( s ) u n ( s ) B ( s ) u 0 ( s ) L ( J , B ) 2 .
By Lemma 4.2 in [41], B is strongly continuous. Also, by Lebesgue’s dominated convergence theorem, we have
G 1 = 5 M T 2 B 1 t + C A t 2 A + 2 0 t E sup 0 s t x s n x s 0 d s 0 as n , G 2 = 5 M T 2 t E B ( s ) u n ( s ) B ( s ) u 0 ( s ) L ( J , B ) 2 0 as n .
For each t J , x n ( · ) , x 0 ( · ) A , we have
E x n ( t ) x 0 ( t ) 2 G 1 + G 2 1 5 B 0 , 5 B 0 < 1 .
So, let us infer that x n x 0 as n . Finally using Balder’s theorem [41] and (H5), we obtain
ϵ = lim n 0 t L ( t , x n ( t ) , x t n , u n ( t ) ) d t 0 t L ( t , x 0 ( t ) , x t 0 , u 0 ( t ) ) d t = l ( u 0 ) ϵ .
Hence, the result is followed that l attains its minimum at u 0 A a d . □
Remark 4.
The main idea of the proof here is similar to that of [41].

7. Applications

Example 1: Consider the following neutral SIDEs driven by Poisson jumps and the Rosenblatt process with infinite delay:
t x ( t , ζ ) + e μ θ ( t ) f 1 ^ ( t , x ( t + θ ( t ) , ζ ) ) = 2 ζ 2 x ( t , ζ ) + e μ θ ( t ) f 1 ^ ( t , x ( t + θ ( t ) , ζ ) ) d t + 0 t B ( s , ζ ) u ( s , ζ ) ^ d s + 0 t b ( t s ) 2 ζ 2 x ( t , ζ ) + e μ θ ( t ) f 1 ^ ( t , x ( t + θ ( t ) , ζ ) ) d t + 0 e 2 μ s f ^ ( t , s , x ( t + s , ζ ) ) d s d t + e μ θ ( t ) σ ( t , x ( t + θ ( t ) , ζ ) ) d Z A ( t ) + e μ θ ( t ) h ^ ( t , x ( t + θ ( t ) , ζ ) , z ) N ˜ ( d t , d z ) , t 0 , ζ [ 0 , π ] , x ( t , 0 ) = x ( t , π ) = 0 , t 0 , x ( θ , ζ ) = x 0 ( θ , ζ ) , θ ( , 0 ] , ζ [ 0 , π ] .
Here, Z A ( t ) denotes a Q-Rosenblatt process, h : R + × R R and f 1 , σ : R + × R R , f : R + × ( , 0 ] × R R are continuous functions, x 0 C ζ and b C 1 ( R + , R ) .
Let A = L 2 ( 0 , π ) with the norm · . Define A : D ( A ) A A by A = 2 ζ 2 with domain D ( A ) = A 2 ( 0 , π ) A 0 1 ( 0 , π ) . Obviously, A generates a strongly continuous semigroup ( S ( t ) ) t 0 on A . So, (A1) is guaranteed. The operator Υ : D ( A ) A A is Υ ( t ) ( x ) = b ( t ) A x .
The resolvent operator ( R ( t ) ) t 0 decays exponentially; i.e.,
R ( t ) M e ϕ t .
Let us assume the following:
(i)
f 1 ^ ( t , 0 ) = f ^ ( t , s , 0 ) = σ ^ ( t , 0 ) = h ^ ( t , 0 , z ) = 0 ; t 0 , s 0 .
(ii)
There exist real numbers B 0 ^ ( 0 , 1 10 ) , B 1 ^ > 0
B 0 ^ x 1 x 2 2 f 1 ^ ( t , x 1 ) f 1 ^ ( t , x 2 ) 2 . B 1 ^ x 1 x 2 2 f ^ ( t , x 1 ) f ^ ( t , x 2 ) 2 σ ^ ( t , x 1 ) σ ^ ( t , x 2 ) 2 z h ^ ( t , x 1 ) h ^ ( t , x 2 ) 2 .
(iii)
∃ an integrable function r : ( , 0 ] [ 0 , + ) such that
f ^ ( t , s , x 1 ) f ^ ( t , s , x 2 ) 2 r ( s ) x 1 x 2 2 ; x 1 , x 2 R .
For ϕ C ζ L 2 0 ( R , A ) , ζ [ 0 , π ] ,
f 1 ( t , ϕ 1 ) ( ζ ) = e μ θ ( t ) f 1 ^ ( t , x ( t + θ ( t ) , ζ ) ) σ ( t , ϕ 1 ) ( ζ ) = e μ θ ( t ) σ ^ ( t , x ( t + θ ( t ) , ζ ) ) h ( t , x ( t + θ ( t ) , ζ ) , z ) = e μ θ ( t ) h ^ ( t , x ( t + θ ( t ) , ζ ) , z ) ,
which involve a variable delay term and the following equation contains a distributed delay term:
f ( t , ϕ 1 ) ( ζ ) = 0 e 2 μ s f ( t , s , ϕ 1 ( s ) ( ζ ) ) d s .
Set x ( t ) ζ = x ( t , ζ ) and ϕ ( θ ) ( ζ ) = x 0 ( θ , ζ ) θ ( , 0 ] . Then, (37) represents (2).
Using condition (ii) with the norm · C ζ yields
B 0 ^ ϕ 1 ϕ 2 C ζ 2 f 1 ( t , ϕ 1 ) f 1 ( t , ϕ 2 ) 2 , B 1 ^ ϕ 1 ϕ 2 C ζ 2 f ^ ( t , ϕ 1 ) f ^ ( t , ϕ 1 ) 2 σ ^ ( t , ϕ 1 ) σ ^ ( t , ϕ 1 ) 2 z h ^ ( t , ϕ 1 ) h ^ ( t , ϕ 1 ) 2 ,
using condition (iii) with Hölder’s inequality,
B 1 ^ ϕ 1 ϕ 2 C ζ 2 f 1 ( t , ϕ 1 ) f 1 ( t , ϕ 2 ) 2 , t 0 and ϕ 1 , ϕ 2 C ζ .
Thus, (37) follows the presumptions of Theorem 1 with B 0 = B 0 ^ and B 1 = B 1 ^ . Theorems 2 and 3 assure the mild solution of (37) with mean square surely and almost surely exponentially stable operators with the following conditions:
5 M 2 B 1 1 ϕ 2 + 2 ϕ < 1 5 B 0 and ϕ < 2 μ .
Let us consider the cost functional as
l ( x , u ) ( ζ ) = 0 π 0 1 E x ( t , ζ ) 2 d t d ζ + 0 π 0 1 u ( t , ζ ) 2 d t d ζ .
Here, B u ( t , ζ ) = 0 t B ( s , ζ ) u ( s , ζ ) ^ d s . Now,
0 π 0 1 B u ( t , ζ ) 2 d t d ζ 0 π 0 t B ( s , ζ ) u ( s , ζ ) ^ d s 2 d t d ζ M 2 0 π 0 t 0 1 u ( t , ζ ) 2 d s d t d ζ M 2 × m e s ( π ) × u ( t , ζ ) L ( [ 0 , 1 ] , [ 0 , π ] ) 2 ,
where M 2 = max B ( t , ζ ) 2 : t [ 0 , 1 ] . Consequently, it may be said that B is a bounded linear operator in L 2 ( [ 0 , 1 ] , [ 0 , π ] ) . Additionally, if Theorem 4’s presumptions are met, at least one optimal pair of ( x ( t , ζ ) , u ( t , ζ ) ) exists.
Example 2: Filter System. An advanced filter is a system that applies mathematical operations to a digitized sign that has been examined in order to tone down or enhance some of the prepared signal’s highlights. We represented a filter design for our framework, which is depicted in Figure 1 driven by the strategies studied in [47]. A simple block diagram pattern is shown in Figure 1 to help improve the feasibility of an arrangement with the fewest input sources available.
  • Modulator Product (MP)-1 acquires the input A, τ α ( σ ) , which represents the output as A τ α ( σ ) .
  • (MP)-2 acquires x ( ι ) , E, and represents E ( ι , x ( ι ) ) .
  • (MP)-3 acquires v ( ι ) , B, and represents B v ( ι ) .
  • (MP)-4 acquires x 0 and R α ( σ ) at time σ = 0 , and represents R α ( σ ) x 0 .
  • The integral of τ α ( σ ) [ E ( σ , x ( σ ) ) + B v ( σ ) ] is executed by integrators over σ .
  • Inputs τ α ( σ ) and E are multiplied and come together with an integrator output over ( 0 , σ ) .
  • τ α ( σ ) and B are multiplied and come together over ( 0 , σ ) .
Finally, all integrator signals are sent to a digital network which attains the output x ( t ) .
x ( σ ) = R α ( σ ) x 0 + 0 σ τ α ( σ ι ) E ( ι , x ( ι ) ) d ι + 0 σ τ α ( σ ι ) B v ( ι ) d ι , σ V .
Example 3: Traffic Signal Light System. If ( Ω , , { t } t 0 , P ) are complete probability space, then { t } , t 0 are natural filters that meet the usual requirement (i.e., they are right continuous and 0 contains all P -null sets). The family of piecewise right continuous functions P C ( [ r , 0 ] , R n ) is defined as ϕ : [ r , 0 ] R n with respect to ϕ = s u p r u 0 | ϕ ( u ) | , for r > 0 . Let P L 0 p ( [ r , 0 ] , R n ) denote the family for all 0 -measurable, P C ( [ r , 0 ] , R n ) valued random variables satisfying s u p r u 0 E | ϕ ( u ) | p < . Let E [ . ] denote the corresponding expectation operator with respect to the probability measure P . Assume that B t = B ( t ) = ( B 1 ( t ) , B 2 ( t ) , , B m ( t ) ) T is an m-dimensional Brownian motion that is specified on the complete probability space. Let τ ( t ) be a right continuous fBm on the probability space ( Ω , , P ) taking values in a finite state space S = { r 1 , r 2 , . . . , r m } with the generator Q = ( q i j ) m × m given by
P ( τ ( t + Δ t ) = r j | τ ( t ) = r i ) = q i j Δ t + o ( Δ t ) , if i j , 1 + q i i Δ t + o ( Δ t ) , if i = j ,
where Δ t > 0 . Here, q i j 0 is the transition rate from r i to r j if r i r j , while q i i = j i q i j = q i . Define ξ 0 = 0 ,   ξ k = i n f { t > ξ k 1 ; τ ( t ) τ ( ξ k 1 ) } and r ( n ) = τ ( ξ n ) . From fBm theory, we know that { r ( n ) , n Z } is an fBm. Its transition probability is r i j ( 1 ) = ( 1 δ i j ) q i j q i . Here, we assume that Rossenblatt process { τ ( t ) , t 0 } is independent of the Brownian motion and the fBm.
Electronic Control System: If the volume of the traffic on the route reaches A, the red light is turned on. The green light is on if the vehicle flows below A. With the generator, we define a two-state Rosenblatt chain { τ ( t ) , t 0 } .
Q = q 11 q 12 q 21 q 22
{the vehicle flow exceeds A} = {the red light on} = { τ ( t ) = r 1 } .
{the vehicle flow does not exceed A} = {the green light on} = { τ ( t ) = r 2 } .
The electronic control system will receive either the red light or the green light signal. However, the sign transport latency varies depending on the sign brightness. In addition, noise will disturb the electrical control system, making impulsive signals (such as a quick voltage change) inevitable. The electronic control system can be characterized by the following two systems, depending on whether the red light (or the green light) is on:
d x ( t ) = f ( t , x ( t ) , x ( t r i ) ) d t + g ( t , x ( t ) , x ( t r i ) ) d B t , x ( t k ) = I k ( t k r i , x ( t k r i ) ) , k = 1 , 2 , , t = t k ;
where i = 1 , 2 . Although the movement of the vehicle is unpredictable, the Rosenblatt chain τ ( t ) governs the random switching of the electronic control system during the aforementioned two systems.
Example 4: Biological impact of an infinite delay on the human body: After overnight fasting, a patient is called to the hospital for the routine yearly body check-up to measure blood sugar and other diseases caused due to the age factor. Last year, he was diagnozed with the normal lab results, but this time he detected diabetes. The patient was curious to know when exactly the diabetes entered the body with the particular month and date, but the doctor showed insufficiency in answering his question. There is no medical instrument which gives us this information so far. This is the practical example of an infinite time delay problem which we have studied in system (2). During this time, the patient took sweets in his diet very frequently and that resulted in boosting up sugar levels, which is a study of stochastic process.
Remark 5.
Ramkumar et al. [48] studied the river-dam-contamination problem as one of the applications of the stochastic delay system with optimal control. By taking f 1 = 0 , system (35) is converted to [48]. Thus, our system is in the more generalized form.

8. Numerical Simulation

In this section, we provide numerical simulations of the application of Equation (38). This example includes several interesting and tricky parts that relate to approximating a numerical solution. If it merely contained a parabolic partial differential equation (PDE), we could solve it using a variety of software packages and the techniques detailed in the cited [38,39] sources. If there was merely a neutral delay term, we could apply the techniques in [40]. All of these elements are represented in this equation along with an integral term.
ζ ϑ ( ζ , θ ) + ζ 2 + e ζ ϑ ( ζ φ , θ ) 2 18 = 2 θ 2 ϑ ( ζ , θ ) + ζ 2 + e ζ ϑ ( ζ φ , θ ) 2 18 + 0 ζ Λ ˜ ( ζ ζ ) ϑ ( ζ , θ ) + ζ 2 + e ζ ϑ ( ζ φ , θ ) 2 18 + c ( ζ , θ ) + e ζ ϑ ( ζ φ , θ ) / 2 + ϑ ( ζ , sin ζ ϑ ( ζ , θ ) ) / 2 9 + ζ e ζ 3 + e ζ d Z A ( ζ ) d ζ + e ζ , θ ( t ) cos θ ν ( ζ ) ( 1 η 2 ) N ˜ ( d t , d z ) , ζ ( 0 , 1 ] , θ [ 0 , π ] , ϑ ( ζ , 0 ) = 0 , ϑ ( ζ , θ ) = φ ( ζ , θ ) , θ [ 0 , π ] ,   ζ ( , 0 ) , η ( 0 , 1 ) .
Now we can write
e ζ , θ ( t ) cos θ ν ( ζ ) ( 1 η 2 ) N ˜ ( d t , d z ) = e ζ , θ ( t ) cos θ ν ( ζ ) ( 1 η 2 ) ϕ d ζ d θ .
When θ = π 3 , η = 0.5 and ϕ = 0.5 , then by using [49], we have
e ζ , θ ( t ) cos θ ν ( ζ ) ( 1 η 2 ) ϕ d ζ d θ ( 0.4193 ) 2 E ν 2 = 0.1758 E ν 2 ,
and we know that
E ν 2 C .
Since (37) with infinite delay and mixed fBm does not yet have an established software suite that can handle all of these components, we created the one using Matlab. The code contains all of the details of the simulation; however, the following are the main methods. The time derivative was made apparent by using a forward finite difference. Each spatial θ derivative was approximated by a centered difference. The integral term was approximated at each discretization point using the trapezoid rule, which is implemented by the Matlab function trapz.m. Difference derivatives on the mesh were also used to roughly estimate the delayed derivative terms.
The simulations provided (Figure 1 and Figure 2) used the following functions and settings. For each timestep, we used n = 2 and a total of 400 spatial points, with 20 points in each spatial dimension. We used 5000 timesteps, so d ι = 0.0002 . The fBm parameter h = 0.7 . The functions are ϕ ( ζ , θ ) = ι 2 + i n θ i 2 , Λ ( s ) = cos ( s ) , and E ν 2 0.5 . Figure 2 shows the function at the beginning, when ι = 0 , and halfway through the experiment. At about the halfway point of the experiment and at the end, Figure 3 shows the simulation.

9. Conclusions

For a class of neutral SIDEs driven by Poisson jumps and Rosenblatt processes utilizing concrete-fading memory-phase space C ζ defined for the stochastic process, we looked into the optimal control difficulties. Grimmer’s resolvent operator theory, successive approximation, and stochastic analysis theory have all been used to examine the existence and uniqueness of mild solutions. The exponential stability, including mean square exponential stability, is then demonstrated. The best-performing pair of systems demonstrates the solution of exponential stability and their maps and is most likely controlled by stochastic partial integrodifferential equations with infinite delay. In order to support the hypothesis, we created a numerical framework. The code includes all of the simulation information for the Matlab finite difference approach. In this study, theoretical proof and numerical estimates are combined in an original way. The following are the research’s potential future applications:
  • The same concept can be expanded utilizing Riemann–Liouville (R-L) and Caputo derivatives for the fractional order/hybrid fractional order system with deviating inputs. Consider the Hilfer fractional system with non-instantaneous impulses and state-dependent delay.
  • The new work with the numerical simulation will be the system’s trajectory controllability.
  • The Measure of Noncompactness method can be used in place of the technique employed in this paper. Additionally, the same system and several types of fractional order SIDEs can be studied using monotone operator theory.
  • The approach of "Integral Contractor with Regularity" can be used to weaken Lipschitz continuity in the nonlinear operators.

Author Contributions

Methodology: D.C. and R.K. (Ramkumar Kasinathan); software: D.C. and R.K. (Ravikumar Kasinathan); validation, R.K. (Ramkumar Kasinathan) and R.K. (Ravikumar Kasinathan); formal analysis: R.K. (Ramkumar Kasinathan) and R.K. (Ravikumar Kasinathan); investigation: D.C.; resources, R.K. (Ramkumar Kasinathan) and R.K. (Ravikumar Kasinathan); data accuration, D.C., R.K. (Ramkumar Kasinathan) and R.K. (Ravikumar Kasinathan); writing—R.K. (Ravikumar Kasinathan); original draft preparation, R.K. (Ramkumar Kasinathan); writing—review and editing; D.C.; visualization, R.K. (Ramkumar Kasinathan) and R.K. (Ravikumar Kasinathan); supervision, D.C.; project administration, R.K. (Ramkumar Kasinathan). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

No animal or human body was involved in conducting this research.

Data Availability Statement

The authors did not use any scientific data yet to be determined for this research.

Acknowledgments

The authors would like to thank the referees and the editor for their careful comments and valuable suggestions to improve this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mao, X. Stochastic Differential Equations and Applications; Horwood Publishing Limited: Chichester, UK, 2007. [Google Scholar]
  2. Oksendal, B. Stochastic differential equations. In Stochastic Differential Equations; Springer: Berlin/Heidelberg, Germany, 2003; pp. 65–84. [Google Scholar]
  3. Prato, D.G.; Zabczyk, J. Stochastic Equations in Infinite Dimensions; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  4. Evans, L.C. An Introduction to Stochastic Differential Equations; University of California: Berkeley, CA, USA, 2013. [Google Scholar]
  5. Ravikumar, K.; Ramkumar, K.; Chalishajar, D. Existence and Stability Results for Second-Order Neutral Stochastic Differential Equations With Random Impulses and Poisson Jumps. Eur. J. Math. Anal. 2021, 1, 1–18. [Google Scholar] [CrossRef]
  6. Helge, H.; Ksendal, B.; Ube, J.; Zhang, T. Stochastic partial differential equations. In Stochastic Partial Differential Equations; Birkhäuser: Boston, MA, USA, 1996; pp. 141–191. [Google Scholar]
  7. Kostikov, Y.A.; Romanenkov, A.M. Approximation of the multidimensional optimal control problem for the heat equation (applicable to computational fluid dynamics (CFD)). Civ. Eng. J. 2020, 6, 743–768. [Google Scholar] [CrossRef]
  8. Trung, T.T. Smart city and modelling of its unorganized flows using cell machines. Civ. Eng. J. 2020, 6, 954–960. [Google Scholar] [CrossRef]
  9. Dieye, M.; Diop, M.A.; Ezzinbi, K. On exponential stability of mild solutions for some stochastic partial integrodifferential equations. Stat. Probab. Lett. 2017, 123, 61–76. [Google Scholar] [CrossRef]
  10. Dieye, M.; Diop, M.A.; Ezzinbi, K. Almost sure asymptotic stability for some stochastic partial functional integrodifferential equations on Hilbert spaces. Cogent Math. Stat. 2019, 6, 1602928. [Google Scholar] [CrossRef]
  11. Taniguchi, T.; Liu, K.; Truman, A. Existence, uniqueness, and asymptotic behavior of mild solutions to stochastic functional differential equations in Hilbert spaces. J. Differ. Equ. 2002, 181, 72–91. [Google Scholar] [CrossRef]
  12. Ramkumar, K.; Ravikumar, K.; Anguraj, A. Existence and Exponential Stability for Neutral Impulsive Stochastic Integrodifferential Equations with Fractional Brownian Motion Driven by Poisson Jumps. J. Vib. Test. Syst. Dyn. 2020, 4, 311–324. [Google Scholar] [CrossRef]
  13. Anguraj, A.; Ravikumar, K. Existence and stability results of impulsive stochastic partial neutral functional differential equations with infinite delays and Poisson jumps. Discontinuity Nonlinearity Complex. 2020, 9, 245–255. [Google Scholar] [CrossRef]
  14. Anguraj, A.; Ramkumar, K.; Elsayed, E.M. Existence, uniqueness and stability of impulsive stochastic partial neutral functional differential equations with infinite delays driven by a fractional Brownian motion. Discontinuity Nonlinearity Complex. 2020, 9, 327–337. [Google Scholar] [CrossRef]
  15. Govindan, T.E. Stability of mild solutions of stochastic evolution equations with variable delay. Stochastics Anal. Appl. 2003, 5, 1059–1077. [Google Scholar] [CrossRef]
  16. Kai, L.; Shi, Y. Razumikhin-type theorems of infinite dimensional stochastic functional differential equations. In Proceedings of the IFIP Conference on System Modeling and Optimization, Turin, Italy, 18–22 July 2005; pp. 237–247. [Google Scholar]
  17. Taniguchi, T. Almost sure exponential stability for stochastic partial functional differential equations. Stoch. Anal. Appl. 1998, 16, 965–975. [Google Scholar] [CrossRef]
  18. Diop, M.A.; Ezzinbi, K.; Lo, M. Exponential stability for some stochastic neutral partial functional integrodifferential equations with delays and Poisson jumps. Semigroup Forum 2014, 88, 595–609. [Google Scholar] [CrossRef]
  19. Chen, H.; Zhu, C.; Zhang, Y. A note on exponential stability for impulsive neutral stochastic partial functional differential equations. Appl. Math. Comput. 2014, 227, 139–147. [Google Scholar] [CrossRef]
  20. Cui, J.; Yan, L. Successive approximation of neutral stochastic evolution equations with infinite delay and Poisson jumps. Appl. Math. Comput. 2012, 218, 6776–6784. [Google Scholar] [CrossRef]
  21. Mao, W.; Hu, L.; Mao, X. Neutral stochastic functional differential equations with Lévy jumps under the local Lipschitz condition. Adv. Differ. Equ. 2017, 2017, 1–24. [Google Scholar] [CrossRef]
  22. Yue, C. Neutral stochastic functional differential equations with infinite delay and Poisson jumps in the Cg space. Appl. Math. Comput. 2014, 237, 595–604. [Google Scholar] [CrossRef]
  23. Grimmer, R.C. Resolvent operators for integral equations in a Banach space. Trans. Am. Math. Soc. 1982, 273, 333–349. [Google Scholar] [CrossRef]
  24. Ren, Y.; Xia, N. Existence, uniqueness and stability of the solutions to neutral stochastic functional differential equations with infinite delay. Appl. Math. Comput. 2009, 210, 72–79. [Google Scholar] [CrossRef]
  25. Hale, K.J.; Kato, J. Phase space for retarded equations with infinite delay. Funck. Ekvacioj 1978, 21, 11–41. [Google Scholar]
  26. Wu, F.; Yin, G.; Mei, H. Stochastic functional differential equations with infinite delay: Existence and uniqueness of solutions, solution maps, Markov properties, and ergodicity. J. Differ. Equ. 2017, 262, 1226–1252. [Google Scholar] [CrossRef]
  27. Maejima, M.; Tudor, C.A. Self-similar processes with stationary increments in the second Wiener chaos. Probab. Math. Stat. 2012, 32, 167–186. [Google Scholar]
  28. Pipiras, V.; Taqqu, M.S. Regularization and integral representations of Hermite processes. Stat. Probab. Lett. 2010, 80, 2014–2023. [Google Scholar] [CrossRef]
  29. Taqqu, M. Weak convergence to fractional Brownian motion and to the Rosenblatt process. Adv. Appl. Probab. 1975, 7, 249. [Google Scholar] [CrossRef]
  30. Tudor, C.A. Analysis of the Rosenblatt process. ESAIM Probab. Stat. 2008, 12, 230–257. [Google Scholar] [CrossRef]
  31. Abry, P.; Pipiras, V. Wavelet-based synthesis of the Rosenblatt process. Signal Process. 2006, 86, 2326–2339. [Google Scholar] [CrossRef]
  32. Leonenko, N.N.; Anh, V.V. Rate of convergence to the Rosenblatt distribution for additive functionals of stochastic processes with long-range dependence. J. Appl. Math. Stoch. Anal. 2001, 14, 27–46. [Google Scholar] [CrossRef]
  33. Boufoussi, B.; Hajji, S. Neutral stochastic functional differential equations driven by a fractional Brownian motion in a Hilbert space. Stat. Probab. Lett. 2012, 82, 1549–1558. [Google Scholar] [CrossRef]
  34. Caraballo, T.; Diop, M.A. Neutral stochastic delay partial functional integro-differential equations driven by a fractional Brownian motion. Front. Math. China 2013, 8, 745–760. [Google Scholar] [CrossRef]
  35. Ramkumar, K.; Ravikumar, K.; Anguraj, A.; Hamdy, M. Ahmed, Well-posedness results for higher-order neutral stochastic differential equations driven by Poisson jumps and Rosenblatt process. Filomat 2021, 35, 353–365. [Google Scholar] [CrossRef]
  36. Caraballo, T.; Garrido-Atienza, M.J.; Taniguchi, T. The existence and exponential behavior of solutions to stochastic delay evolution equations with a fractional Brownian motion. Nonlin. Anal. Theory Methods Appl. 2011, 74, 3671–3684. [Google Scholar] [CrossRef]
  37. Urszula, L.; Schättler, H. Antiangiogenic therapy in cancer treatment as an optimal control problem. SIAM J. Control. Optim. 2007, 46, 1052–1079. [Google Scholar]
  38. Area, I.; Ndairou, F.; Nieto, J.J.; Silva, C.J.; Torres, D.F. Ebola model and optimal control with vaccination constraints. J. Ind Manag Optim. 2018, 14, 427–446. [Google Scholar] [CrossRef]
  39. Rajivganthi, C.; Muthukumar, P. Almost automorphic solutions for fractional stochastic differential equations and its optimal control. Optim. Control. Appl. Methods 2016, 37, 663–681. [Google Scholar] [CrossRef]
  40. Wang, Q.; Chen, F.; Huang, F. Maximum principle for optimal control problem of stochastic delay differential equations driven by fractional Brownian motions. Optim. Control Appl. Methods 2016, 37, 90–107. [Google Scholar] [CrossRef]
  41. Ramkumar, K.; Ravikumar, K.; Varshini, S. Fractional neutral stochastic differential equations with Caputo fractional derivative: Fractional Brownian motion, Poisson jumps, and optimal control. Stoch. Anal. Appl. 2021, 39, 157–176. [Google Scholar] [CrossRef]
  42. Diop, M.A.; Ezzinbi, K.; Guindo, P.D.A. Optimal controls problems for some impulsive stochastic integro-differential equations with state-dependent delay. Stochastics 2022, 94, 1–35. [Google Scholar] [CrossRef]
  43. Diop, M.A.; Guindo, P.D.A.; Fall, M.; Diakhby, A. Optimal controls for stochastic functional integrodifferential equations. Electron. J. Math. Anal. Appl. 2021, 9, 241–261. [Google Scholar]
  44. Hajji, S.; Lakhel, E. Existence and uniqueness of mild solutions to neutral SFDE driven by a fractional Brownian motion with non-Lipschitz coefficients. arXiv 2013, arXiv:1312.6147. [Google Scholar]
  45. Boufoussi, B.; Hajji, S.; Lakhel, E.H. Functional differential equations in Hilbert spaces driven by a fractional Brownian motion. Afr. Mat. 2012, 23, 173–194. [Google Scholar] [CrossRef]
  46. Zouine, A.; Bouzahir, H.; Vargas, A.N. Stability for stochastic neutral integro-differential equations with infinite delay and Poisson jumps. RMS Res. Math. Stat. 2021, 8, 1979733. [Google Scholar] [CrossRef]
  47. Zahoor, S.; Naseem, S. Design and implementation of an efficient FIR digital filter. Cogent Eng. 2017, 4, 1323373. [Google Scholar] [CrossRef]
  48. Ramkumar, K.; Ravikumar, K.; Chalishajar, D. Existence trajectory and optimal control of Clarke subdifferential stochastic integrodifferential inclusions suffered by non-instantaneous impulses and deviated arguments. Result Control. Optim. 2023, 13, 100295. [Google Scholar] [CrossRef]
  49. Mattuvarkuzhali, C.; Balasubramaniam, P. pth Moment stability of fractional stochastic differential inclusions via resolvent operators driven by the Rosenblatt process and Poisson jumps with impulses. Stochastics 2019, 92, 1157–1174. [Google Scholar] [CrossRef]
Figure 1. Filter System:block diagram: the feasibility with the minimum input source arrangement. All letters in the figure correspond to the Definition 2.
Figure 1. Filter System:block diagram: the feasibility with the minimum input source arrangement. All letters in the figure correspond to the Definition 2.
Fractalfract 07 00783 g001
Figure 2. Graph of ϑ ( ζ = 0 , θ ) and ϑ ( ζ = 0.24925 , θ ) . Function at the start of the simulation, when ι = 0 , and a third of the way through. Blue indicates lower values of ϑ , yellow indicates higher values of ϑ .
Figure 2. Graph of ϑ ( ζ = 0 , θ ) and ϑ ( ζ = 0.24925 , θ ) . Function at the start of the simulation, when ι = 0 , and a third of the way through. Blue indicates lower values of ϑ , yellow indicates higher values of ϑ .
Fractalfract 07 00783 g002
Figure 3. ϑ ( ζ = 0.4995 , θ ) and ϑ ( ζ = 1 , θ ) . Simulation at approximately the two-thirds point of the simulation. Blue indicates lower values of ϑ , yellow indicates higher values of ϑ .
Figure 3. ϑ ( ζ = 0.4995 , θ ) and ϑ ( ζ = 1 , θ ) . Simulation at approximately the two-thirds point of the simulation. Blue indicates lower values of ϑ , yellow indicates higher values of ϑ .
Fractalfract 07 00783 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chalishajar, D.; Kasinathan, R.; Kasinathan, R. Optimal Control for Neutral Stochastic Integrodifferential Equations with Infinite Delay Driven by Poisson Jumps and Rosenblatt Process. Fractal Fract. 2023, 7, 783. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7110783

AMA Style

Chalishajar D, Kasinathan R, Kasinathan R. Optimal Control for Neutral Stochastic Integrodifferential Equations with Infinite Delay Driven by Poisson Jumps and Rosenblatt Process. Fractal and Fractional. 2023; 7(11):783. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7110783

Chicago/Turabian Style

Chalishajar, Dimplekumar, Ramkumar Kasinathan, and Ravikumar Kasinathan. 2023. "Optimal Control for Neutral Stochastic Integrodifferential Equations with Infinite Delay Driven by Poisson Jumps and Rosenblatt Process" Fractal and Fractional 7, no. 11: 783. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7110783

Article Metrics

Back to TopTop