Next Article in Journal
Information Fragmentation, Encryption and Information Flow in Complex Biological Networks
Previous Article in Journal
Ecological Function Analysis and Optimization of a Recompression S-CO2 Cycle for Gas Turbine Waste Heat Recovery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Control of Background-Based Uncertain Systems with Applications in DC Pension Plan

1
College of Mathematics and System Science, Xinjiang University, Urumqi 830046, China
2
School of Mathematics and Statistics, Wuhan University, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
Submission received: 21 March 2022 / Revised: 14 May 2022 / Accepted: 15 May 2022 / Published: 21 May 2022

Abstract

:
In this paper, we propose a new optimal control model for uncertain systems with jump. In the model, the background-state variables are incorporated, where the background-state variables are governed by an uncertain differential equation. Meanwhile, the state variables are governed by another uncertain differential equation with jump, in which both the background-state variables and the control variables are involved. Under the optimistic value criterion, using uncertain dynamic programming method, we establish the principle and the equation of optimality. As an application, the optimal investment strategy and optimal payment rate for DC pension plans are given, where the corresponding background-state variables represent the salary process. This application in DC pension plans illustrates the effectiveness of the proposed model.

1. Introduction

Since [1] first discussed the problem of stochastic control with jump, the problem of stochastic optimal control with jump has become an important branch of control theory. From the perspectives of theory and applications, especially including in finance and insurance, the involved stochastic differential equations with and without jump have been extensively studied. For example, see [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20], and the references therein.
The above-mentioned references are based on probability models. However, in reality, the finance markets usually are of model uncertainty, which means that it is difficult to determine the specific probability. Therefore, it is of great importance to study uncertainty theory with its applications in finance and insurance. For general theory and applications about uncertainty theory and optimal control of uncertain systems, we refer to [21,22,23,24,25,26,27,28,29,30,31] and the references therein. For the applications of uncertainty theory in option pricing theory and portfolio selections, see [26,27,28,29,32,33,34,35], and the references therein. For the applications of uncertainty theory in insurance, especially in pension plans, see [22,36]. In the study of optimal control of uncertain systems, the optimality criteria mainly include four criteria: expected value criterion, optimistic value criterion, pessimistic value criterion, and Hurwicz criterion. Under the optimistic value criterion, the principle of optimality and the equation of optimality for uncertain systems without jump were discussed by [26]. Recently, Ref. [27] studied the optimal control of uncertain systems with jump under the optimistic value criterion, where the state variables are governed by an uncertain differential equation with jump. Later, Ref. [37] extended those of [27] to the multidimensional setting. Nevertheless, from both the theoretical and practical point of view, the state variables are usually also affected by the environment factors except the control variables. Therefore, it is interesting and necessary to consider the optimal control of uncertain systems by incorporating the environment factors into the optimal control models under the optimistic value criterion.
In this paper, we propose a new optimal control model for uncertain systems under the optimistic value criterion. Namely, the environment factors are first understood as background variables. Then we assume that the background-state variables are governed by an uncertain differential equation, and further we assume that the state variables are also governed by another uncertain differential equation with jump in which both the background-state variable and the control variables are involved. By making use of the uncertain dynamic programming method, both the principle and the equation of optimality are established. Finally, as an application, the optimal investment strategy and the optimal payment rate for DC pension plans are discussed, where the corresponding background-state variables represent the salary process. This application in DC pension plans illustrates the effectiveness of the proposed model.
The rest of the paper is organized as follows. In Section 2, we introduce preliminaries, including basic notations of uncertainty theory. In Section 3, the optimistic value models for uncertain systems with jump are introduced and the principle of optimality is provided. Section 4 is devoted to the equation of optimality. In Section 5, as an application of the proposed optimal control model in DC pension plans, the optimal investment strategy and the optimal payment rate are obtained. Section 6 presents numerical analysis to illustrate our results. Finally, the conclusions are summarized.

2. Preliminary

2.1. Uncertainty Space

In this subsection, we collect some basic definitions of uncertainty theory which are from [21,23,24].
Let Γ be a nonempty set, and L a σ -algebra over Γ . Each element A L is called an event. A set function M defined on the σ -algebra L over Γ is called an uncertain measure if it satisfies the following four conditions:
(i)
M { Γ } = 1 ,
(ii)
M { Λ } + M { Λ c } = 1 for any event A L ,
(iii)
M { i = 1 Λ i } i = 1 M { Λ i } for every countable sequence of events Λ i .
(iv)
Let ( Γ k , L k , M k ) be uncertainty spaces for k = 1 , 2 , n . The product uncertain measure is
M i = 1 n Λ k = min 1 i n M { Λ k } .
Definition 1.
Let Γ be a nonempty set, L a σ-algebra over Γ and M an uncertain measure. Then the triplet ( Γ , L , M ) is called an uncertainty space.
Definition 2.
An uncertain variable is a function ξ from an uncertainty space ( Γ , L , M ) to the set of real numbers such that { ξ B } : = { γ Γ ξ ( γ ) B } s an event for any Borel set B.
Definition 3.
The uncertainty distribution Φ : R [ 0 , 1 ] of an uncertain variable ξ is defined by
Φ ( x ) : = M { ξ x } , x R .
The following lemma is a characterization of an uncertainty distribution, which is from [25].
Lemma 1.
A function Φ ( x ) : R [ 0 , 1 ] is an uncertainty distribution if and only if it is a monotone increasing function except Φ ( x ) = 0 and Φ ( x ) = 1 .
Definition 4.
Let ξ be an uncertain variable. Then the expected value of ξ is defined by
E [ ξ ] : = 0 + M { ξ x } d x 0 M { ξ x } d x ,
provided that at least one of the two integrals is finite.
Definition 5.
Let ξ be an uncertain variable with finite expected value E [ ξ ] . Then the variance of ξ is defined by
V [ ξ ] : = E [ ( ξ E [ ξ ] ) 2 ] .
Definition 6.
The uncertain variables ξ 1 , ξ 2 , , ξ n are said to be independent if
M i = 1 n ( ξ i B i ) = min 1 i n M { ξ i B i } ,
for any Borel sets B 1 , B 2 , · · · , B n of real numbers.

2.2. Optimistic Value and Pessimistic Value

Definition 7.
Let ξ be an uncertain variable, and α ( 0 , 1 ] . Then ξ sup ( α ) : = sup { r M { ξ r } α } is called the α-optimistic value to ξ; and ξ inf ( α ) : = inf { r M { ξ r } α } is called the α-pessimistic value to ξ.
The following lemma is about properties of optimistic value, which is from [21,24].
Lemma 2.
Assume that ξ and η are independent uncertain variables and α ( 0 , 1 ] . Then we have
(i) 
if λ 0 , then ( λ ξ ) sup ( α ) = λ ξ sup ( α ) , and ( λ ξ ) inf ( α ) = λ ξ inf ( α ) ;
(ii) 
λ < 0 , then ( λ ξ ) sup ( α ) = λ ξ inf ( α ) , and ( λ ξ ) inf ( α ) = λ ξ sup ( α ) ;
(iii) 
( ξ + η ) sup ( α ) = ξ sup ( α ) + η sup ( α ) , ( ξ + η ) inf ( α ) = ξ inf ( α ) + η inf ( α ) .
Definition 8.
An uncertain process C t is said to be a Liu process if
(i) 
C 0 = 0 and almost all sample paths are Lipschitz continuous,
(ii) 
C t has stationary and independent increments,
(iii) 
every increment C s + t C s is a normal distributed uncertain variable with expected value 0 and variance t 2 , whose uncertainty distribution is
Φ ( x ) : = 1 + exp π x 3 t 1 , x R .
Let C t be a Liu process, and ξ : = Δ C t : = C t + Δ t C t . Then for any 0 < α < 1 , α -optimistic and α -pessimistic values of ξ are
ξ sup ( α ) = 3 Δ t π ln 1 α α
and
ξ inf ( α ) = 3 Δ t π ln 1 α α ,
respectively, for example, see Example 1.7 of [30] or (1) and (2) of [27].
The following definitions are about optimal control with jump of uncertainty theory, which are from [38].
Definition 9.
An uncertain process V t is said to be a V jump process with parametrs r 1 and r 2 ( 0 < r 1 < r 2 < 1 ) for t 0 if
(i) 
V 0 = 0 ,
(ii) 
V t has stationary and independent increments,
(iii) 
every increment V s + t V s is a Z jump uncertain variable Z ( r 1 , r 2 , t ) , whose uncertainty distribution is
Ψ ( x ) : = 0 , if x < 0 , 2 r 1 t x , if 0 x t 2 , r 2 + 2 ( 1 r 2 ) t ( x t 2 ) , if t 2 x < t , 1 , if x t .
Let V t be a V jump uncertain process, and η = Δ V t = V t + Δ t V t . Then for any α ( 0 , 1 ) , it follows from the definition of α -optimistic value and α -pessimistic value that
η sup ( α ) = ( 1 α 2 ( 1 r 2 ) ) Δ t , if 0 < α < 1 r 2 , Δ t 2 , if 1 r 2 α < 1 r 1 , 1 α 2 r 1 Δ t , if 1 r 1 α < 1 ,
and
η inf ( α ) = α 2 r 1 Δ t , if 0 < α r 1 , Δ t 2 , if r 1 < α r 2 , ( 1 1 α 2 ( 1 r 2 ) ) Δ t , if r 2 < α < 1 ,
respectively.
Definition 10.
Suppose that C t is an Liu process, f and g are two functions. Then
d X t = f ( t , X t ) d t + g ( t , X t ) d C t
is called an uncertain differential equation. A solution is a uncertain process X t that satisfies (6) identically in t.
Remark 1.
The uncertain differential Equation (6) means the solution X t meets the uncertain integral equation
X t = X 0 + 0 t f ( t , X t ) d s + 0 t g ( t , X t ) d C s .
Definition 11.
Suppose that C t is a Liu process, V t is an uncertain V jump process, and g 1 , g 2 and g 3 are some given functions. Then
d X t = g 1 ( t , X t ) d t + g 2 ( t , ( X t ) d C t + g 3 ( t , X t ) d V t
is called an uncertain differential equation with jump. A solution is an uncertain process X t that satisfies (7) identically in t.

3. Optimistic Value Model under Background-State of Uncertain Optimal Control with Jump

In the problem of uncertain optimal control, we should determine some optimization criteria to optimize objective function of involving uncertain variables and convert the uncertain objective to its definite equivalent goal. In the uncertain optimal control, there are many criteria, for example: expected value, optimistic value, pessimistic value and Hurwicz criterion. Under [27], they discussed optimal control problem of uncertain dynamic systems with jump under the optimistic value criterion. In this paper, we involve the background-state variables and discuss optimal control problem under the optimistic value criterion for this kind of systems, where both the background-state variables and the control variables are involved.
Assume that C t is a Liu process, V t is an uncertain V-jump process with parameters r 1 and r 2 ( 0 < r 1 < r 2 < 1 ) , where C t and V t are independent. The confidence level α ( 0 , 1 ) . For any 0 < t < T , an optimistic value model of uncertain optimal control problem with jump as follows
J ( t , x , l ) = sup D t D t T f ( s , X s , L s , D s ) d s + G ( T , X T , L T ) sup ( α ) subject to d X s = ν ( s , X s , L s , D s ) d s + λ ( s , X s , L s , D s ) d C s + χ ( s , X s , L s , D s ) d V s , d L s = m ( L s ) d s + n ( L s ) d C s , X t = x , L t = l ,
where X s is the state variable, L s is called the background-state variable, D s is the control variable and it subject to a constraint set D . The function f : [ 0 , T ] × R × R × R R is the objective function, and the function G : [ 0 , T ] × R × R R is the terminal reward. [ t T f ( s , X s , L s , D s ) d s + G ( T , X T , L T ) ] sup ( α ) denotes the α -optimistic value to the uncertain variable in middle bracket. ν , λ , χ are three functions of time s, state X s , background-state L s and control D s . Furthermore, m, n are two functions of L s . All functions mentioned are continuous. Now we give the following principle of optimality to solve the proposed model.
Theorem 1.
For any ( t , x , l ) [ 0 , T ) × R × R , Δ t > 0 with t + Δ t < T , we have
J ( t , x , l ) = sup D t D { f ( t , x , l , D t ) Δ t + J ( t + Δ t , x + Δ X t , l + Δ L t ) + o ( Δ t ) } sup ( α ) ,
where x + Δ X t = X t + Δ t , l + Δ L t = L t + Δ t .
Proof. 
We use J ¯ ( t , x , l ) to denote the right side of (9). It follows from the definition of J ( t , x , l ) that
J ( t , x , l ) { t t + Δ t f ( s , X s , L s , D s | [ t , t + Δ t ) ) d s + t + Δ t T f ( s , X s , L s , D s | [ t + Δ t , T ] ) d s + G ( T , X T , L T ) } sup ( α ) ,
for any D s , where D s | [ t , t + Δ t ) and D s | [ t + Δ t , T ] are the values of decision variable D s restricted on [ t , t + Δ t ) and [ t + Δ t , T ] , respectively.
For any Δ t > 0 , by using Taylor series expansion, we get
t t + Δ t f ( s , X s , L s , D s | [ t , t + Δ t ) ) d s = f ( t , x , l , D ( t , x , l ) ) Δ t + o ( Δ t ) .
Thus
J ( t , x , l ) f ( t , x , l , D t ) Δ t + o ( Δ t ) + t + Δ t T f ( s , X s , L s , D s | [ t + Δ t , T ] ) d s + G ( T , X T , L T ) sup ( α ) .
Taking the supremum with respect to D s | [ t + Δ t , T ] in (10), then we get
sup D t D t + Δ t T f ( s , X s , L s , D s | [ t + Δ t , T ] ) d s + G ( T , X T , L T ) sup ( α ) = J ( t + Δ t , x + Δ X t , l + Δ L t ) .
Then the right side of (10) becomes that
sup D t D f ( t , x , l , D t ) Δ t + J ( t + Δ t , x + Δ X t , l + Δ L t ) + o ( Δ t ) sup ( α ) = J ¯ ( t , x , l ) .
Now we get J ( t , x , l ) J ¯ ( t , x , l ) .
On the other hand, for all D s , we have
t T f ( s , X s , L s , D s ) d s + G ( T , X T , L T ) sup ( α ) = f ( t , x , l , D t ) Δ t + o ( Δ t ) + t + Δ t T f ( s , X s , L s , D s | [ t + Δ t , T ] ) d s + G ( T , X T , L T ) sup ( α ) f ( t , x , l , D t ) Δ t + o ( Δ t ) + [ J ( t + Δ t , x + Δ X t , l + Δ L t ) ] J ¯ ( t , x , l ) .
Hence, J ( t , x , l ) J ¯ ( t , x , l ) , and then J ( t , x , l ) = J ¯ ( t , x , l ) . Theorem 1 is proved. □

4. Optimality Condition

We derive the following equation of optimality by the principle of optimality above.
Theorem 2.
Let J ( t , x , l ) be twice continuously differentiable on [ 0 , T ] × R × R . Then we have
J t = sup D t D f + ν J x + m ( L t ) J l + 3 π ln 1 α α ( λ J x + n ( L t ) J l ) + k | χ J x | ,
where J t , J x and J l are the partial derivatives of the function J ( t , x , l ) in t, x and l, respectively, and f, ν, λ, χ, J t , J x , J l , J t t , J x x , J l l , J t x , J t l , J x l denote f ( t , x , l , D t ) , ν ( t , x , l , D t ) , λ ( t , x , l , D t ) , χ ( t , x , l , D t ) , J t ( t , x , l ) , J x ( t , x , l ) , J l ( t , x , l ) , J t t ( t , x , l ) , J x x ( t , x , l ) , J l l ( t , x , l ) , J t x ( t , x , l ) , J t l ( t , x , l ) , J x l ( t , x , l ) , respectively, and
(1) 
when χ J x 0 , (i) if 0 < α < 1 r 2 , then k = 1 α 2 ( 1 r 2 ) , (ii) if 1 r 2 α < 1 r 1 , then k = 1 2 , (iii) if 1 r 1 α < 1 , then k = 1 α 2 r 1 ;
(2) 
when χ J x < 0 , (i) if 0 < α < r 1 , then k = α 2 r 1 , (ii) if r 1 < α r 2 , then k = 1 2 , (iii) if r 2 < α < 1 , then k = 1 1 α 2 ( 1 r 2 ) .
Remark 2.
Suppose state variable X t does not depend on L t , the J in (8) is a function of ( t , x ) , therefore Theorem 2 goes back to the classic setting, for example, see the model in [37].
Proof. 
For any Δ t > 0 , note that there exist constants K 1 > 0 and K 2 > 0 , such that Δ X t K 1 Δ t and Δ L t K 2 Δ t , by Taylor series expansion, we have that
J ( t + Δ t , x + Δ X t , l + Δ L t ) = J ( t , x , l ) + J t ( t , x , l ) Δ t + J x ( t , x , l ) Δ X t + J l ( t , x , l ) Δ L t + 1 2 J t t ( t , x , l ) Δ t 2 + 1 2 J x x ( t , x , l ) Δ X t 2 + 1 2 J l l ( t , x , l ) Δ L t 2 + J t x ( t , x , l ) Δ t Δ X t + J t l ( t , x , l ) Δ t Δ L t + J x l ( t , x , l ) Δ X t Δ L t + o ( Δ t 2 ) ,
where the infinitesimal o ( Δ t 2 ) satisfies that
lim Δ t 0 o ( Δ t 2 ) Δ t 2 = 0 .
Note that Δ X t = ν ( t , X t , L t , D t ) Δ t + λ ( t , X t , L t , D t ) Δ C t + χ ( t , X t , L t , D t ) Δ V t , Δ L t = m ( L t ) Δ t + n ( L t ) Δ C t . Substituting (12) into (9) yields that
J ( t , x , l ) = sup D t D { f ( t , x , l , D t ) Δ t + [ J ( t , x , l ) + J t ( t , x , l ) Δ t + J x ( t , x , l ) Δ X t + 1 2 J t t ( t , x , l ) Δ t 2 + 1 2 J x x ( t , x , l ) Δ X t 2 + 1 2 J l l ( t , x , l ) Δ L t 2 + J l ( t , x , l ) Δ L t + J t x ( t , x , l ) Δ t Δ X t + J t l ( t , x , l ) Δ t Δ L t + J x l ( t , x , l ) Δ X t Δ L t + o ( Δ t 2 ) ] sup ( α ) } .
Then we have
0 = sup D t D { f ( t , x , l , D t ) Δ t + J t ( t , x , l ) Δ t + [ J x ( t , x , l ) Δ X t + J l ( t , x , l ) Δ L t + 1 2 J t t ( t , x , l ) Δ t 2 + 1 2 J x x ( t , x , l ) Δ X t 2 + 1 2 J l l ( t , x , l ) Δ L t 2 + J t x ( t , x , l ) Δ t Δ X t + J t l ( t , x , l ) Δ t Δ L t + J x l ( t , x , l ) Δ X t Δ L t ] sup ( α ) + o ( Δ t 2 ) } .
Let uncertain variable ξ = Δ C t , η = Δ V t , then it follows from the uncertain differential equation in model (8), that
Δ X t = ν Δ t + λ ξ + χ η , Δ L t = m ( L t ) Δ t + n ( L t ) ξ ,
where ν , λ , χ denote ν ( t , X t , L t , D t ) , λ ( t , X t , L t , D t ) , χ ( t , X t , L t , D t ) , respectively.
Substituting (15) into (14) yields that
0 = sup D t D { f ( t , x , l , D t ) Δ t + J t Δ t + ν J x Δ t + m ( L t ) J l Δ t + [ ( λ J x + n ( L t ) J l + ν λ J x x Δ t + m ( L t ) n ( L t ) J l l Δ t + λ J t x Δ t + n ( L t ) J t l Δ t + ν n ( L t ) J x l Δ t + λ m ( L t ) J x l Δ t ) ξ + 1 2 λ 2 J x x + 1 2 n ( L t ) 2 J l l + λ n ( L t ) J x l ξ 2 + χ J x + ν χ J x x Δ t + χ J t x Δ t + m ( L t ) χ J x l Δ t η + 1 2 χ 2 J x x η 2 + λ χ J x x + n ( L t ) χ J x l ξ η ] sup ( α ) + o ( Δ t 2 ) } = sup D t D f Δ t + J t Δ t + ν J x Δ t + m ( L t ) J l Δ t + a ξ + b ξ 2 + p η + q η 2 + f ξ η sup ( α ) + o ( Δ t 2 ) ,
where
a = λ [ J x + ν J x x Δ t + J t x Δ t + m ( L t ) J x l Δ t ] + n ( L t ) [ J l + m ( L t ) J l l Δ t + J t l Δ t + ν J x l Δ t ] , b = λ 2 1 2 J x x + 1 2 n ( L t ) 2 J l l + λ n ( L t ) J x l , p = χ [ J x + ν J x x Δ t + J t x Δ t + m ( L t ) J x l Δ t ] , q = χ 2 1 2 J x x , f = λ χ J x x + n ( L t ) χ J x l .
Let Q ( ξ , η ) = a ξ + b ξ 2 + p η + q η 2 + f ξ η , thus the equation becomes
0 = sup D t D f Δ t + J t Δ t + ν J x Δ t + m ( L t ) J l Δ t + [ Q ] sup ( α ) + o ( Δ t 2 ) ,
where [ Q ( ξ , η ) ] sup ( α ) is denoted by [ Q ] sup ( α ) .
Since | b ξ 2 | 1 2 | b | ( ξ 2 + ξ 2 ) , | q η 2 | 1 2 | q | ( η 2 + η 2 ) , | f ξ η | 1 2 | f | ( ξ 2 + η 2 ) , we have
a ξ + p η ( | b | + 1 2 | f | ) ξ 2 + ( 1 2 | f | + | q | ) η 2 a ξ + b ξ 2 + p η + q η 2 + f ξ η a ξ + p η + ( | b | + 1 2 | f | ) ξ 2 + ( 1 2 | f | + | q | ) η 2
It follows from the independence of ξ and η that
[ a ξ + p η ( | b | + 1 2 | f | ) ξ 2 + ( 1 2 | f | + | q | ) η 2 ] sup ( α ) [ a ξ + b ξ 2 + p η + q η 2 + f ξ η ] sup ( α ) a ξ + p η + ( | b | + 1 2 | f | ) ξ 2 + ( 1 2 | f | + | q | ) η 2 sup ( α )
According to Theorem 4 in Sheng and Zhu (2013), for any ε > 0 small enough, we have
a ξ + p η ( | b | + 1 2 | f | ) ξ 2 sup ( α ) 3 π ln 1 α ε α + ε | a | Δ t ( 3 π ln 2 ε ε ) 2 ( | b | + 1 2 | f | ) Δ t 2
a ξ + p η + ( | b | + 1 2 | f | ) ξ 2 sup ( α ) 3 π ln 1 α + ε α ε | a | Δ t + ( 3 π ln 2 ε ε ) 2 ( | b | + 1 2 | f | ) Δ t 2
According to Theorem 5.1 in Deng et al. (2018), we have
(1)
if p 0 , then
p η ( | q | + 1 2 | f | ) η 2 sup ( α ) ( 1 α 2 ( 1 r 2 ) ) p Δ t ( | q | + 1 2 | f | ) Δ t 2 , if 0 < α < 1 r 2 , 1 2 p Δ t ( | q | + 1 2 | f | ) Δ t 2 , if 1 r 2 α < 1 r 1 , 1 α 2 r 1 p Δ t ( | q | + 1 2 | f | ) Δ t 2 , if 1 r 1 α < 1 ,
p η + ( | q | + 1 2 | f | ) η 2 sup ( α ) ( 1 α 2 ( 1 r 2 ) ) p Δ t + ( | q | + 1 2 | f | ) Δ t 2 , if 0 < α < 1 r 2 , 1 2 p Δ t + ( | q | + 1 2 | f | ) Δ t 2 , if 1 r 2 α < 1 r 1 , 1 α 2 r 1 p Δ t + ( | q | + 1 2 | f | ) Δ t 2 , if 1 r 1 α < 1 ,
(2)
if p < 0 , then
p η ( | q | + 1 2 | f | ) η 2 sup ( α ) α 2 r 1 p Δ t ( | q | + 1 2 | f | ) Δ t 2 , if 0 < α r 1 , 1 2 p Δ t ( | q | + 1 2 | f | ) Δ t 2 , if r 1 < α r 2 , ( 1 1 α 2 ( 1 r 2 ) ) p Δ t ( | q | + 1 2 | f | ) Δ t 2 , if r 2 < α < r 1 ,
p η + ( | q | + 1 2 | f | ) η 2 sup ( α ) α 2 r 1 p Δ t + ( | q | + 1 2 | f | ) Δ t 2 , if 0 < α r 1 , 1 2 p Δ t + ( | q | + 1 2 | f | ) Δ t 2 , if r 1 < α r 2 , ( 1 1 α 2 ( 1 r 2 ) ) p Δ t + ( | q | + 1 2 | f | ) Δ t 2 , if r 2 < α < r 1 ,
Therefore,
(1)
if p 0 , then
[ Q ] sup ( α ) [ 3 π ln 1 α ε α + ε | a | + ( 1 α 2 ( 1 r 2 ) ) p ] Δ t h Δ t 2 , if 0 < α < 1 r 2 , [ 3 π ln 1 α ε α + ε | a | + 1 2 p ] Δ t h Δ t 2 , if 1 r 2 α < 1 r 1 , [ 3 π ln 1 α ε α + ε | a | + 1 α 2 r 1 p ] Δ t h Δ t 2 , if 1 r 1 α < 1 ,
[ Q ] sup ( α ) [ 3 π ln 1 α + ε α ε | a | + ( 1 α 2 ( 1 r 2 ) ) p ] Δ t + h Δ t 2 , if 0 < α < 1 r 2 , [ 3 π ln 1 α + ε α ε | a | + 1 2 p ] Δ t + h Δ t 2 , if 1 r 2 α < 1 r 1 , [ 3 π ln 1 α + ε α ε | a | + 1 α 2 r 1 p ] Δ t + h Δ t 2 , if 1 r 1 α < 1 ,
(2)
if p < 0 , then
[ Q ] sup ( α ) [ 3 π ln 1 α ε α + ε | a | + α 2 r 1 p ] Δ t h Δ t 2 , if 0 < α < 1 r 2 , [ 3 π ln 1 α ε α + ε | a | + 1 2 p ] Δ t h Δ t 2 , if 1 r 2 α < 1 r 1 , [ 3 π ln 1 α ε α + ε | a | + ( 1 1 α 2 ( 1 r 2 ) ) p ] Δ t h Δ t 2 , if 1 r 1 α < 1 ,
[ Q ] sup ( α ) [ 3 π ln 1 α + ε α ε | a | + α 2 r 1 p ] Δ t + h Δ t 2 , if 0 < α < 1 r 2 , [ 3 π ln 1 α + ε α ε | a | + 1 2 p ] Δ t + h Δ t 2 , if 1 r 2 α < 1 r 1 , [ 3 π ln 1 α + ε α ε | a | + ( 1 1 α 2 ( 1 r 2 ) ) p ] Δ t + h Δ t 2 , if 1 r 1 α < 1 ,
where h = [ ( 3 π ln 2 ε ε ) 2 ( | b | + 1 2 | f | ) + | q | + 1 2 f ] .
(1)
If p 0 , then by (17), (27), for Δ t > 0 , if 0 < α < 1 r 2 , there exists a control u = u ε , Δ t such that
ε Δ t f Δ t + J t Δ t + ν J x Δ t + m ( L t ) J l Δ t + [ Q ] sup ( α ) + o ( Δ t ) f Δ t + J t Δ t + ν J x Δ t + m ( L t ) J l Δ t + 3 π ln 1 α + ε α ε | a | Δ t + 1 α 2 ( 1 r 2 ) p Δ t + h Δ t 2 + o ( Δ t ) .
Since
a = λ [ J x + ν J x x Δ t + J t x Δ t + m ( L t ) J x l Δ t ] + n ( L t ) [ J l + m ( L t ) J l l Δ t + J t l Δ t + ν J x l Δ t ] , p = χ [ J x + ν J x x Δ t + J t x Δ t + m ( L t ) J x l Δ t ] .
Substituting them into (30) and dividing both sides of the above inequality by Δ t , we get
ε f + J t + ν J x + m ( L t ) J l + 3 π ln 1 α + ε α ε | λ J x + n ( L t ) J l + h 1 ( Δ t ) | + 1 α 2 ( 1 r 2 ) χ J x + h 2 ( Δ t ) + h 3 ( ε , Δ t ) J t + sup D t D { f + ν J x + m ( L t ) J l + 3 π ln 1 α + ε α ε | λ J x + n ( L t ) J l + h 1 ( Δ t ) | + 1 α 2 ( 1 r 2 ) χ J x + h 2 ( Δ t ) } + h 3 ( ε , Δ t ) ,
where
h 1 ( Δ t ) = λ ( ν J x x Δ t + J t x Δ t + m ( L t ) J x l Δ t ) + n ( L t ) ( m ( L t ) J l l Δ t + J t l Δ t + ν J x l Δ t ) , h 2 ( Δ t ) = χ ( ν J x x Δ t + J t x Δ t + m ( L t ) J x l Δ t ) , h 3 ( ε , Δ t ) = h Δ t ,
and h 1 ( Δ t ) 0 , h 2 ( Δ t ) 0 , h 3 ( ε , Δ t ) 0 as Δ t 0 . Letting Δ t 0 and then ε 0 results in
0 J t + sup D t D { f + ν J x + m ( L t ) J l + 3 π ln 1 α α | λ J x + n ( L t ) J l | + 1 α 2 ( 1 r 2 ) χ J x } ,
if χ J x 0 and 0 < α < 1 r 2 .
In the same way, by (17) and (26), we can get
0 J t + sup D t D { f + ν J x + m ( L t ) J l + 3 π ln 1 α α | λ J x + n ( L t ) J l | + 1 α 2 ( 1 r 2 ) χ J x } ,
if χ J x 0 and 0 < α < 1 r 2 .
Combining (32) and (33), we obtain
J t = sup D t D f + ν J x + m ( L t ) J l + 3 π ln 1 α α | λ J x + n ( L t ) J l | + 1 α 2 ( 1 r 2 ) χ J x ,
if χ J x 0 and 0 < α < 1 r 2 .
According to (17), (26) and (27), using the similar techniques, we have
J t = sup D t D f + ν J x + m ( L t ) J l + 3 π ln 1 α α | λ J x + n ( L t ) J l | + 1 2 χ J x ,
if χ J x 0 , 1 r 2 α < 1 r 1 and
J t = sup D t D f + ν J x + m ( L t ) J l + 3 π ln 1 α α | λ J x + n ( L t ) J l | + 1 α 2 r 1 χ J x ,
if χ J x 0 and 1 r 1 α < 1 .
(2)
If p < 0 , similar to the above method, we use (17), (29), (33) to derive the equation of optimality for χ J x < 0 . Therefore, Theorem 2 is proved. □

5. An Optimal Control Problem of DC Pension Fund

In recent years, pension fund management has become a popular and significant subject to retirees because it plays an essential role in the financial market and in the social security system. The dynamic control-theoretical framework was first applied to a pension fund problem by [4] by assuming that the pension fund can be invested in a risk-free asset and a risky asset whose return follows random jump-diffusion processes. Ref. [27] assumed risk asset returns follow an uncertain process with jump and made use of optimistic value criterion to optimize objective function of involving uncertain variables. We assume that the contribution of pension is related to the salary factors of members, then the DC pension plan control problem may be solved by the optimistic value model of uncertain optimal control with jump.

5.1. Finance Market

We assume that the financial market consists of two assets, a risk-free asset (i.e., the bank account or bond), and a single risk asset (i.e., stock).
The price of the risk-free asset S 0 ( t ) at the time t evolves according to the following uncertain process
d S 0 ( t ) S 0 ( t ) = r d t ,
where r is a constant and represents the risk-free interest rate.
The price of the risk asset S ( t ) at the time t evolves according to the following uncertain process with jump
d S ( t ) S ( t ) = μ d t + σ 1 d C ( t ) + σ 2 d V ( t ) ,
where μ is the appreciation rate of the risk asset and σ 1 , σ 2 is the volatility rate, μ , σ 1 , and σ 2 are all positive constants, and C ( t ) is a canonical process, V ( t ) is a V-jump uncertain process. In general, we assume that μ > r .
The salary level is denoted by L ( t ) at time t. We assume that L ( t ) follows a uncertain growth given by
d L ( t ) L ( t ) = μ L d t + σ L d C ( t ) ,
where μ L is the expected rate of return on salary, σ L is the salary volatility caused by the fluctuation of risk asset.
We assume that the pension contribution rate is θ L ( t ) , where 0 < θ < 1 is a constant.

5.2. Wealth Process

Assuming that the plan managers can invest in both the risk-free and the risky assets described by (37) and (38), respectively, and use the fund to pay retirement benefits. Let x 0 denote the initial wealth of this fund, ω ( t ) denote the investment proportion that the plan managers invest in the risky asset at time t, and X ( t ) denote the wealth of the pension fund at time t after adapting the investment strategy ω ( t ) , B ( t ) is the pension payment rate at time t. Then the fund’s value follows the dynamics
d X ( t ) = [ 1 ω ( t ) ] X ( t ) d S 0 ( t ) S 0 ( t ) + ω ( t ) X ( t ) d S ( t ) S ( t ) + θ L ( t ) [ μ L d t + σ L d C ( t ) ] B ( t ) d t , X ( 0 ) = x 0 .
Using (37)–(39), we can easily rewrite (40) as
d X ( t ) = [ X ( t ) r + ( μ r ) ω ( t ) X ( t ) + θ μ L L ( t ) B ( t ) ] d t + [ σ 1 ω ( t ) X ( t ) + σ L θ L ( t ) ] d C ( t ) + [ ω ( t ) X ( t ) σ 2 ] d V ( t ) , X 0 = x 0 .

5.3. Optimization Model

Our goal is to seek the optimal investment strategy ω ( t ) and payment rate B ( t ) to minimize the accumulated losses, thus we establish the following optimal model of pension fund.
J ( t , x , l ) = min ω ( t ) , B ( t ) t e ρ s α 1 ( B ( s ) b m ) 2 + α 2 ( ω ( s ) X ( s ) x p ) 2 d s inf ( α ) subject to d X ( t ) = [ X ( t ) r + ( μ r ) ω ( t ) X ( t ) + θ μ L L ( t ) B ( t ) ] d t + [ σ 1 ω ( t ) X ( t ) + σ L θ L ( t ) ] d C ( t ) + [ ω ( t ) X ( t ) σ 2 ] d V ( t ) , d L ( t ) = L ( t ) [ μ L d t + σ L d C ( t ) ] , X t = x , L ( t ) = l ,
where, α 1 > 0 , α 2 > 0 , and α 1 + α 2 = 1 . α ( 0 , 1 ) denotes a given confidence level, ρ > 0 denotes the discount rate. b m denotes the constant target contribution rate and x p denotes the constant target funding level.
It follows from Lemma 2 that model (42) is equivalent to the following model (43).
J ( t , x , l ) = max ω ( t ) , B ( t ) t e ρ s α 1 ( B ( s ) b m ) 2 + α 2 ( ω ( s ) X ( s ) x p ) 2 d s sup ( α ) subject to d X ( t ) = [ X ( t ) r + ( μ r ) ω ( t ) X ( t ) + θ μ L L ( t ) B ( t ) ] d t + [ σ 1 ω ( t ) X ( t ) + σ L θ L ( t ) ] d C ( t ) + [ ω ( t ) X ( t ) σ 2 ] d V ( t ) , d L ( t ) = L ( t ) [ μ L d t + σ L d C ( t ) ] , X t = x , L ( t ) = l .

5.4. The Solution to the Model

By applying the equation of optimality (11), we get
J t = max ω ( t ) , B ( t ) { e ρ t [ α 1 ( B ( t ) b m ) 2 + α 2 ( ω ( t ) x x p ) 2 ] + x r + ( μ r ) ω ( t ) x + θ μ L l B ( t ) J x + μ L l J l + 3 π ln 1 α α | [ σ 1 ω ( t ) x + σ L θ l ] J x + σ L l J l | + k ω ( t ) x σ 2 J x } = max ω ( t ) , B ( t ) H ( ω ( t ) , B ( t ) ) ,
where H ( ω ( t ) , B ( t ) ) represents the term in the brackets.
Now we solve the (44)
(1)
If ( σ 1 ω ( t ) x + σ L θ l ) J x + σ L l J l 0 , we differentiate the expression in brackets with respect to ω ( t ) and B ( t ) to find that
H ( ω ( t ) , B ( t ) ) ω ( t ) = 2 e ρ t α 2 ( ω ( t ) x x p ) x + ( μ r ) x J x + 3 π ln 1 α α σ 1 x J x + k x σ 2 J x = 0 ,
H ( ω ( t ) , B ( t ) ) B ( t ) = 2 e ρ t α 1 ( B ( t ) b m ) J x = 0 .
Solving Equations (45) and (46), we get
ω ( t ) = 1 x [ x p + 1 2 α 2 k ˜ e ρ t J x ] , B ( t ) = b m 1 2 α 1 e ρ t J x ,
where k ˜ = μ r + 3 π ln 1 α α σ 1 + k σ 2 .
Substituting them into (44) implies
J t = a ˜ e ρ t J x 2 + ( x r + θ c ˜ l + b ˜ ) J x + c ˜ l J l ,
where a ˜ = 1 4 α 1 + 1 4 α 2 k ˜ 2 , b ˜ = k ˜ x p b m , c ˜ = μ L + 3 π ln 1 α α σ L .
Multiplying both sides of equation by e ρ t
e ρ t J t = a ˜ ( e ρ t J x ) 2 + ( x r + θ c ˜ l + b ˜ ) e ρ t J x + c ˜ l e ρ t J l .
Next we solve the partial differential Equation (48).
Supposing J ( t , x , l ) = e ρ t Q ( x , l ) , then differentiating both sides with respect to t, x, and l, then J t = ρ e ρ t Q ( x , l ) , J x = e ρ t Q x , J l = e ρ t Q l . Substituting them into (48) yields
ρ Q ( x , l ) = a ˜ Q x 2 + ( x r + θ c ˜ l + b ˜ ) Q x + c ˜ l Q l .
Assuming Q ( x , l ) = A ( x 2 + 2 h x l + b l 2 + a x + g l + u ) , then Q x = 2 A x + 2 A h l + A a , Q l = 2 A b l + 2 A h x + A g ) . Substituting them into (49) yields
( ρ A 4 A 2 a ˜ 2 A r ) x 2 + ( ρ A a 4 A 2 a a ˜ A a r 2 A b ˜ ) x + ( ρ A b 4 A 2 h 2 a ˜ 2 A h θ c ˜ 2 A b c ˜ ) l 2 + ( ρ A g 4 A 2 h a a ˜ A a θ c ˜ 2 A h b ˜ A g c ˜ ) l + ( 2 ρ A h 8 A 2 h a ˜ 2 A h r 2 A θ c ˜ 2 A h c ˜ ) x l + ( ρ A u A 2 a 2 a ˜ A a b ˜ ) = 0 .
Decomposing Equation (50) obtains
ρ A 4 A 2 a ˜ 2 A r = 0 , ρ A a 4 A 2 a a ˜ A a r 2 A b ˜ = 0 , ρ A b 4 A 2 h 2 a ˜ 2 A h θ c ˜ 2 A b c ˜ = 0 , ρ A g 4 A 2 h a a ˜ A a θ c ˜ 2 A h b ˜ A g c ˜ = 0 , 2 ρ A h 8 A 2 h a ˜ 2 A h r 2 A θ c ˜ 2 A h c ˜ = 0 , ρ A u A 2 a 2 a ˜ A a b ˜ = 0 .
By solving Equation (51), we get
A = ρ 2 r 4 a ˜ , a = 2 b ˜ r , h = θ c ˜ r c ˜ , b = θ c ˜ r c ˜ 2 , g = 2 θ b ˜ c ˜ r ( r c ˜ ) , u = b ˜ 2 r 2 .
Thus,
J x = e ρ t Q x , Q x = 2 A x + 2 A h l + A a = ( ρ 2 r ) 2 a ˜ x + θ c ˜ ( ρ 2 r ) 2 a ˜ ( r c ˜ ) l + b ˜ ( ρ 2 r ) 2 a ˜ r .
Then,
J x = e ρ t ( ρ 2 r ) 2 a ˜ x + θ c ˜ ( ρ 2 r ) 2 a ˜ ( r c ˜ ) l + b ˜ ( ρ 2 r ) 2 a ˜ r .
So the optimal investment rate and the payment rate are determined, respectively, by
ω * ( t ) = 1 x [ x p + 1 2 α 2 k ˜ e ρ t J x ] = 1 x x p + 1 2 α 2 k ˜ ( ρ 2 r ) 2 a ˜ x + θ c ˜ ( ρ 2 r ) 2 a ˜ ( r c ˜ ) l + b ˜ ( ρ 2 r ) 2 a ˜ r ,
B * ( t ) = b m 1 2 α 1 e ρ t J x = b m 1 2 α 1 ( ρ 2 r ) 2 a ˜ x + θ c ˜ ( ρ 2 r ) 2 a ˜ ( r c ˜ ) l + b ˜ ( ρ 2 r ) 2 a ˜ r .
(2)
If ( σ 1 ω ( t ) x + σ L θ l ) J x + σ L l J l < 0 , then applying the similar method to the above processes, we can get results (55) and (56), where k ˜ = μ r 3 π ln 1 α α σ 1 + k σ 2 .
Theorem 3.
For the optimization model (43), the optimal investment strategy is given by
ω * ( t ) = 1 x x p + 1 2 α 2 k ˜ ( ρ 2 r ) 2 a ˜ x + θ c ˜ ( ρ 2 r ) 2 a ˜ ( r c ˜ ) l + b ˜ ( ρ 2 r ) 2 a ˜ r ,
the payment rate is given by
B * ( t ) = b m 1 2 α 1 ( ρ 2 r ) 2 a ˜ x + θ c ˜ ( ρ 2 r ) 2 a ˜ ( r c ˜ ) l + b ˜ ( ρ 2 r ) 2 a ˜ r ,
where
a ˜ = 1 4 α 1 + 1 4 α 2 k ˜ 2 , b ˜ = k ˜ x p b m , c ˜ = μ L + 3 π ln 1 α α σ L .
If ( σ 1 ω ( t ) x + σ L θ l ) J x + σ L l J l 0 ,
k ˜ = μ r + 3 π ln 1 α α σ 1 + k σ 2 ,
if ( σ 1 ω ( t ) x + σ L θ l ) J x + σ L l J l < 0 ,
k ˜ = μ r 3 π ln 1 α α σ 1 + k σ 2 .
Remark 3.
The optimal payment rate B * ( t ) , the optimal investment proportion ω * ( t ) and the optimal value J ( t , x , l ) depend on the parameters x p , b m , α 1 , α 2 , ρ , r , θ , l , μ , μ L , σ 1 , σ 2 , σ L and the fund amount x.

6. Numerical Analysis

In this section, we provide a numerical analysis to characterize the dynamic behavior of the optimal investment strategy and the optimal payment rate. We fix the parameter values according to the modeling background. r = 0.02 , μ = 0.05 , σ 1 = 0.15 , σ 2 = 0.2 , μ L = 0.03 , σ L = 0.06 , θ = 0.02 , α = 0.95 , α 1 = 0.6 , α 2 = 1 α 1 , x 0 = 0 , x p = 5 , p m = 5 , ρ = 0.05 , r 1 = 0.2 , r 2 = 0.7 , when α = 0.2 , k = 1 α / ( 2 ( 1 r 2 ) ) , when α = 0.5 , k = 0.5 , when α = 0.8 , k = ( 1 α ) / ( 2 r 1 ) .
Figure 1 shows the effect of parameter α on the optimal investment proportion ω * ( t ) and the optimal payment rate B * ( t ) . From Figure 1a,b, we can see that both ω * ( t ) and B * ( t ) increase when α increases, with all other parameters being fixed. The confidence level α reflects the risk preference of a pension fund manager. Larger α means that the pension fund manager is risk averse. In other words, the pension fund manager would like more prudently to run the fund to achieve his/ her expected management targets. Figure 1 says that for a prudent manager, he/she can invest more in the financial market so that can make more profit, while he/she can also pay more money to the retirees.
Effects of model main parameters μ , σ 1 and σ 2 on ω * ( t ) and B * ( t ) are shown in Figure 2. The graphs in Figure 2a–f plot the values of the optimal investment proportion ω * ( t ) and the optimal payment rate B * ( t ) with respect to the wealth x at time 0, when the parameters μ , σ 1 and σ 2 influencing the stock’s price change. From Figure 2a,b, we can find that the values of ω * ( t ) and B * ( t ) increase as μ increases. This is consistent with the intuition that when the return rate μ of the stock becomes higher, the pension fund manager would naturally like to invest more in the stock to make more profit. At the same time, the pension fund manager is able to pay more to the retirees who participate in the plan. From Figure 2c–f, we can see that ω * ( t ) and B * ( t ) decrease as σ 1 and σ 2 increases, respectively. The corresponding economic explanation is as follows. Higher values of σ 1 and σ 2 represent higher uncertainty of the fluctuation of the price of the stock. In other words, the higher the values of σ 1 and σ 2 are, the more risky the stock is. Therefore, a risk averse pension fund manager would most likely to reduce the amount invested in the stock, and has to lower the payment rate to the pension members because of the reduction of the profit from the stock market.
In Figure 3, we examine the sensitivity of ω * ( t ) and B * ( t ) with respect to the parameters influencing salary level. From Figure 3a–d, we can see that both ω * ( t ) and B * ( t ) are insensitive to μ L and σ L .

7. Conclusions

In this paper, we have proposed a new optimal control model for uncertain systems. Unlike the classic optimal control model for uncertain systems, the proposed new optimal control model takes into account environmental factors. Under the optimistic value criterion, we have established the principle and equation of optimality. As an example, an application to DC pension plans has been given to illustrate the proposed optimal control model for uncertain systems. Numerical studies are also given to show the sensitivity analysis of the optimal solution to the model parameters. For further topics, it would be interesting to consider the principle and equation of optimality under other optimality criteria such as the Hurwitz criterion.

Author Contributions

Conceptualization, W.L.; Supervision, Y.H.; Writing—original draft, W.W. and X.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China grant number 11961064.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Merton, R. Optimum consumption and portfolio rules in a continuous-time model. J. Econ. Theory 1971, 3, 373–413. [Google Scholar] [CrossRef] [Green Version]
  2. Fleming, W.; Rishel, R. Deterministic and Stochastic Optimal Control; Springer: New York, NY, USA, 1975. [Google Scholar]
  3. Whittle, P. Optimization over Time: Dynamic Programming and Stochastic Control; Wiley: New York, NY, USA, 1983. [Google Scholar]
  4. Boulier, J.; Trussant, E.; Florens, D. A dynamic model for pension funds management. In Proceedings of the Fifth AFIR International Colloquium, Brussels, Belgium, 7–8 September 1995; pp. 361–384. [Google Scholar]
  5. Yong, J.; Zhou, X. Stochastic Controls: Hamiltonian Systems and HJB Equations; Springer: New York, NY, USA, 1999. [Google Scholar]
  6. Taksar, M. Optimal risk and dividend distribution control models for an insurance company. Math. Methods Oper. Res. 2000, 51, 1–42. [Google Scholar] [CrossRef]
  7. Yang, H.; Zhang, L. Optimal investment for insurer with jump-diffusion risk process. Insur. Math. Econ. 2005, 37, 615–634. [Google Scholar] [CrossRef]
  8. Schmidli, H. Stochastic Control in Discrete Time. In Stochastic Control in Insurance; Springer: London, UK, 2008. [Google Scholar]
  9. Bjork, T.; Murgoci, A. A general theory of markovian time inconsistent stochastic control problems. SSRN Electron. J. 2010, 336, 1694759. [Google Scholar] [CrossRef] [Green Version]
  10. Fleming, W.H.; Soner, H.M. Controlled Markov Process and Viscosity Solutions; Springer: New York, NY, USA, 2010. [Google Scholar]
  11. Josa-Fombellida, R.; Rinco-Zapatero, J. Stochastic pension funding when the benefit and the risky asset follow jump diffusion processes. Eur. J. Oper. Res. 2012, 220, 404–413. [Google Scholar] [CrossRef]
  12. He, L.; Liang, Z. Optimal investment strategy for the DC plan with the return of premiums clauses in a mean variance framework. Insur. Math. Econ. 2013, 53, 643–649. [Google Scholar] [CrossRef]
  13. Guan, G.; Liang, Z. Optimal management of DC pension plan in a stochastic interest rate and stochastic volatility framework. Insur. Math. Econ. 2014, 57, 58–66. [Google Scholar] [CrossRef]
  14. Sun, J.; Li, Z.; Zeng, Y. Precommitment and equilibrium investment strategies for defined contribution pension plans under a jump-diffusion model. Insur. Math. Econ. 2016, 67, 158–172. [Google Scholar] [CrossRef]
  15. Nkeki, C.; Mcaleer, M. Optimal investment and optimal additional voluntary contribution rate of a DC pension fund in a jump diffusion environment. Ann. Financ. Econ. 2018, 12, 1750017. [Google Scholar] [CrossRef]
  16. Zhao, P.; Zhou, B.; Wang, J. Non–Gaussian Closed Form Solutions for Geometric Average Asian Options in the Framework of Non-Extensive Statistical Mechanics. Entropy 2018, 20, 71. [Google Scholar] [CrossRef] [Green Version]
  17. Mwanakatwe, P.; Song, L.; Hagenimana, E.; Wang, X. Management strategies for a defined contribution pension fund under the hybrid stochastic volatility model. Comput. Appl. Math. 2019, 38, 45. [Google Scholar] [CrossRef]
  18. Lima, L. Nonlinear Stochastic Equation within an It? Prescription for Modelling of Financial Market. Entropy 2019, 21, 530. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Lindgren, J. Efficient Markets and Contingent Claims Valuation: An Information Theoretic Approach. Entropy 2020, 22, 1283. [Google Scholar] [CrossRef] [PubMed]
  20. Zhang, X.; Guo, J. Optimal defined contribution pension management with salary and risky assets following jump diffusion processes. East Asian J. Appl. Math. 2020, 10, 22–39. [Google Scholar]
  21. Liu, B. Uncertainty Theory, 2nd ed.; Springer: Berlin, Germany, 2007. [Google Scholar]
  22. Liu, B. Fuzzy process, hybrid process and uncertain process. J. Uncertain Syst. 2008, 2, 3–16. [Google Scholar]
  23. Liu, B. Some research problems in uncertainty theory. J. Uncertain Syst. 2009, 3, 3–10. [Google Scholar]
  24. Liu, B. Uncertainty Theory: A Branch of Mathematics for Modeling Human Uncertainty; Springer: Berlin, Germany, 2010. [Google Scholar]
  25. Liu, B. Uncertainty Theory, 5th ed.; Springer: Berlin, Germany, 2020. [Google Scholar]
  26. Sheng, L.; Zhu, Y. Optimistic value model of uncertain optimal control. Int. J. Uncertain. Fuzziness Knowl.-Based Syst. 2013, 21 (Suppl. 01), 75–87. [Google Scholar] [CrossRef]
  27. Deng, L.; Chen, Y. Optimal control of uncertain systems with jump under optimistic value criterion. Eur. J. Control 2017, 38, 7–15. [Google Scholar] [CrossRef]
  28. Zhu, Y. Uncertain optimal control with application to a portfolio selection model. Cybern. Syst. Int. J. 2010, 41, 535–547. [Google Scholar] [CrossRef]
  29. Zhu, Y. Uncertain fractional differential equations and an interest rate model. Math. Methods Appl. Sci. 2015, 38, 3359–3368. [Google Scholar] [CrossRef]
  30. Zhu, Y. Uncertain Optimal Control; Springer Nature: Singapore, 2019. [Google Scholar]
  31. Deng, L.; Shen, J.; Chen, Y. Hurwicz model of uncertain optimal control with jump. Math. Methods Appl. Sci. 2020, 43, 10054–10069. [Google Scholar] [CrossRef]
  32. Huang, X.; Di, H. Uncertain portfolio selection with background risk. Appl. Math. Comput. 2016, 276, 284–296. [Google Scholar] [CrossRef]
  33. Li, B.; Zhu, Y.; Sun, Y.; Aw, G.; Teo, K. Multi-period portfolio selection problem under uncertain environment with bankruptcy constraint. Appl. Math. Model. 2018, 56, 539–550. [Google Scholar] [CrossRef] [Green Version]
  34. Lu, Z.; Yan, H.; Zhu, Y. European option pricing model based on uncertain fractional differential equation. Fuzzy Optim. Decis. Mak. 2019, 18, 199–217. [Google Scholar] [CrossRef]
  35. Jin, T.; Sun, Y.; Zhu, Y. Time integral about solution of an uncertain fractional order differential equation and application to zero-coupon bond model. Appl. Math. Comput. 2020, 372, 124991. [Google Scholar] [CrossRef]
  36. Liu, Y.; Yang, M.; Zhai, J.; Bai, M. Portfolio selection of the defined contribution pension fund with uncertain return and salary: A multi-period mean-variance model. J. Intell. Fuzzy Syst. 2018, 34, 2363–2371. [Google Scholar] [CrossRef]
  37. Deng, L.; You, Z.; Chen, Y. Optimistic value model of multidimensional uncertain optimal control with jump. Eur. J. Control 2018, 39, 1–7. [Google Scholar] [CrossRef]
  38. Deng, L.; Zhu, Y. Uncertain optimal control with jump. ICIC Express Lett. Part B Appl. 2012, 3, 419–424. [Google Scholar]
Figure 1. (a) Effect of α on the optimal investment proportion ω * ( t ) ; (b) Effect of α on the optimal payment rate B * ( t ) .
Figure 1. (a) Effect of α on the optimal investment proportion ω * ( t ) ; (b) Effect of α on the optimal payment rate B * ( t ) .
Entropy 24 00734 g001
Figure 2. (a) Effect of μ on the optimal investment proportion ω * ( t ) ; (b) Effect of μ on the optimal payment rate B * ( t ) ; (c) Effect of σ 1 on the optimal investment proportion ω * ( t ) ; (d) Effect of σ 1 on the optimal payment rate B * ( t ) ; (e) Effect of σ 2 on the optimal investment proportion ω * ( t ) ; (f) Effect of σ 2 on the optimal payment rate B * ( t ) .
Figure 2. (a) Effect of μ on the optimal investment proportion ω * ( t ) ; (b) Effect of μ on the optimal payment rate B * ( t ) ; (c) Effect of σ 1 on the optimal investment proportion ω * ( t ) ; (d) Effect of σ 1 on the optimal payment rate B * ( t ) ; (e) Effect of σ 2 on the optimal investment proportion ω * ( t ) ; (f) Effect of σ 2 on the optimal payment rate B * ( t ) .
Entropy 24 00734 g002aEntropy 24 00734 g002b
Figure 3. (a) Effect of μ L on the optimal investment proportion ω * ( t ) ; (b) Effect of μ L the optimal payment rate B * ( t ) ; (c) Effect of σ L on the optimal investment proportion ω * ( t ) ; (d) Effect of σ L the optimal payment rate B * ( t ) .
Figure 3. (a) Effect of μ L on the optimal investment proportion ω * ( t ) ; (b) Effect of μ L the optimal payment rate B * ( t ) ; (c) Effect of σ L on the optimal investment proportion ω * ( t ) ; (d) Effect of σ L the optimal payment rate B * ( t ) .
Entropy 24 00734 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, W.; Wu, W.; Tang, X.; Hu, Y. Optimal Control of Background-Based Uncertain Systems with Applications in DC Pension Plan. Entropy 2022, 24, 734. https://0-doi-org.brum.beds.ac.uk/10.3390/e24050734

AMA Style

Liu W, Wu W, Tang X, Hu Y. Optimal Control of Background-Based Uncertain Systems with Applications in DC Pension Plan. Entropy. 2022; 24(5):734. https://0-doi-org.brum.beds.ac.uk/10.3390/e24050734

Chicago/Turabian Style

Liu, Wei, Wanying Wu, Xiaoyi Tang, and Yijun Hu. 2022. "Optimal Control of Background-Based Uncertain Systems with Applications in DC Pension Plan" Entropy 24, no. 5: 734. https://0-doi-org.brum.beds.ac.uk/10.3390/e24050734

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop