Next Article in Journal
Market Misreaction? Leverage and Mergers and Acquisitions
Previous Article in Journal
A 3-Dimensional Frame of Reference for Prevention of Risk in Supply Chain
Previous Article in Special Issue
Mathematical Foundations for Balancing the Payment System in the Trade Credit Market
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Control Strategies for the Premium Policy of an Insurance Firm with Jump Diffusion Assets and Stochastic Interest Rate

1
Laboratoire de Mathématiques Appliquées, Université de Biskra, PB 145, Biskra 07000, Algeria
2
Departament de Matemàtiques i Informàtica, Universitat de Barcelona, Gran Via 585, 08007 Barcelona, Spain
*
Author to whom correspondence should be addressed.
J. Risk Financial Manag. 2022, 15(3), 143; https://doi.org/10.3390/jrfm15030143
Submission received: 17 November 2021 / Revised: 1 March 2022 / Accepted: 14 March 2022 / Published: 17 March 2022
(This article belongs to the Special Issue Mathematical Finance with Applications)

Abstract

:
In this paper, we present a stochastic optimal control model to optimize an insurance firm problem in the case where its cash-balance process is assumed to be described by a stochastic differential equation driven by Teugels martingales. Noticing that the insurance firm is able to control its cash-balance dynamics by regulating the underlying premium rate, the aim of the policy maker is to select an appropriate premium in order to minimize the total deviation of the state process to some pre-set target level. As a part of stochastic maximum principle approach, a verification theorem is used to fulfill this achievement.

1. Introduction

An insurance is a contract, represented by a policy, used as a method of protection against losses, whether big or small. This means that the insured can receive or reimburse some financial amount to offset his or her losses from an insurance company. Due the fact that an insurance premium is the amount of money that a person or a company ought to pay for an insurance policy, one can perceive it in two different ways. In one side, it can be considered as an income by the insurance company. On the other side, it can also be considered as a liability in that the insurer must provide coverage for claims being made against the policy. In the present paper, we focus on the first case, where the policy maker looks forward to maximize the terminal wealth of its firm’s cash-balance under a demand law.
The main problem in optimal control theory is to characterize an optimal control process. There are two main approaches, the Pontryagin’s maximum principle and the Bellman’s dynamic programming. We note that optimal control theory has been used for example in Hipp and Vogt (2003) in order to determine an optimal dynamic unlimited excess of loss reinsurance strategy to minimize the infinite time ruin probability, see also Højgaard and Taksar (1998), where the authors applied a proportional reinsurance policy for diffusion models in order to find a policy that maximizes a given return function before the time of ruin. In Cairns (2000), Cairns studied the optimization problem of stochastic pension fund models in continuous time. The mean-variance portfolio management for an insurance company was studied by Josa-Fombellida and Rincón-Zapatero (2008) by using dynamic programming techniques and also by Xie et al. (2008) using the general stochastic linear quadratic control technique. In Moore and Young (2006), Moore and Young obtained optimal dynamic consumption, investment, and insurance strategies, using the dynamic programming principle and a Markov chain approximation method. Ngwira and Gerrad in Ngwira and Gerrard (2007) showed that the optimal contribution and asset allocation policies have similar forms as in the pure diffusion case, but with a modification due to the effect of jumps. In Huang et al. (2010), Huang et al. explicitly derived the insurance company’s optimal premium strategy and the associated optimal cost function. We refer the reader to a list of recent papers and the references therein for more details in this subject, Asmussen et al. (2019); Li and Young (2021); Lin et al. (2020); Supian and Mamat (2021).
Motivated by the above results, in this paper, we solve an optimal premium policy problem of an insurance firm. The main tool used in proving our main results is the Pontryagin maximum principle. More precisely, the sufficient condition of optimality. Noting that, in pretty much all of the previous papers, the authors dealt with the problem of optimal insurance in continuous-time models, we impose here to work with a quite general semi-martingale framework assuming that the liability process is driven by both a Brownian motion and a family of pairwise orthogonal martingales associated with a Lévy process. This kind of models comes naturally from the fact that in many real cases, the continuity of the trajectories condition cannot be satisfied. Indeed, the empirical distribution of the cash balance process tends to deviate from normal distributions, either due to inspected dusters or huge profits, many successive incidents or even because of the lack of continuity in the real world of applications.
The paper is organized as follows. In Section 2, we introduce the preliminaries and formulate our control problem. In Section 3, we derive a sufficient condition of optimality for stochastic differential equations driven by a Brownian motion and a family of Teugels martingales. In Section 4, we apply our result to study problems of optimal premium rate of an insurance firm with jump diffusion assets with stochastic interest rate.

2. Preliminaries and Problem Statement

2.1. Preliminaries

Let T > 0 and Ω , F , F , P be a complete filtered probability space supporting a one-dimensional standard Brownian motion W and a one-dimensional Lévy process L with triplet ( γ , σ 2 , ν ) , defined on [ 0 , T ] , independent each other. Recall that γ R , σ 2 > 0 and ν are Lévy measures. We denote by E the expectation with respect to P . Moreover, F : = { F t , t [ 0 , T ] } denotes the complete right-continuous natural filtration generated by processes W and L .
Recall briefly the L 2 theory of Lévy processes as it is presented in Nualart and Schoutens (2000). Assume there exists α > 0 such that for every ε > 0 ,
( ε , ε ) c e α z ν d z < .
Note that this guarantee that ν has moments of all orders.
Denote by Δ L t : = L t L t the jumps of the Lévy process and define the power jump processes L 1 t : = L t and
L i t : = 0 < s t Δ L s i , i 2 .
If we define
Y i t : = L i t E L i t , i 1 ,
then the family of Teugels martingales { H ( i ) ( t ) } i = 1 is defined by
H i t : = j = 1 i a i j Y j t
where coefficients a i j correspond to the orthonormalization of polynomials 1 , x , x 2 , … with respect to measure μ d x = x 2 ν d x + σ 2 δ 0 d x . Recall that { H i t } i = 1 is a family of strongly orthogonal martingales such that for any i and j, H i , H j t = δ i j t and H i , H j H i , H j is a martingale.
Throughout this article, we will use the following spaces:
  • l 2 ( R ) : the Hilbert space of real-valued sequences x = x n n 0 with norm
    i = 1 x i 2 1 2 < .
  • P 2 R : the Hilbert space of l 2 ( R ) -valued functions f i i 0 such that
    i = 1 0 T | f i t | 2 d t 1 2 < .
  • l F 2 0 , T ; R : the Hilbert space of P 2 R -valued F -predictable processes f i t i 0 such that
    E 0 T i = 1 | f i t | 2 d t 1 2 < .
  • L F 2 0 , T ; R : the Hilbert space of real-valued F -adapted processes f t such that
    E 0 T f t 2 d t 1 2 < .
  • S F 2 0 , T ; R : the Banach space of real-valued, F -adapted and càdlàg processes f t such that
    E sup 0 t T f t 2 1 2 < .
  • L 2 Ω , F , P , R : the Hilbert space of real-valued square integrable random variables defined on Ω , F , P .

2.2. Description of the Model

In this paper, we handle a cash management problem of an insurance firm (also called cash-balance problem). The insurer aims to manage the firm’s operating cash in order to meet demand. We consider that the insurance firm can constantly modify its premium policy while investing in certain financial accounts. Its concern is to pay due benefits, but at the same time to stabilize the insurance schemes by avoiding a large deviation on the cash balance. We address an optimal premium policy problem of an insurance firm under stochastic interest rate when the liability process, also called the payment function, B t , is modeled by the stochastic differential equation:
d B t = b t + v t d t + a t d W t + i = 1 π i t d H i t , B ( 0 ) = 0 ,
where b is a locally bounded function that denotes the liability rate, that is, the expected liability per unit time due to premium loading, v is the premium rate (premium policy) and a and π i i = 1 are locally bounded functions that denote the volatility rates measuring the liability risks belonging respectively to the Brownian and Teugels martingale components.
For simplicity, from now on, we will write π ( t ) : = π i ( t ) i = 1 , H ( t ) : = H i t i = 1 and
π ( t ) d H ( t ) : = i = 1 π i t d H i t .
Assume moreover that the cash balance process of the insurer X · is described by the formula
X t = e Δ t X 0 0 t e Δ s d B s ,
where X 0 = x 0 represents the initial reserve and Δ ( · ) represents the stochastic interest rate, that follows the dynamics
d Δ t = δ t d t + α t d W t , t 0 , T , Δ 0 = 0 ,
with δ and α locally bounded functions.
Note that X ( t ) is the difference between the initial capital and the net expenses up to time t.
Now, the Itô formula applied to the process X · leads to the following controlled Stochastic Differential Equation (SDE):
d X t = f 1 t , X t , v t d t + σ 1 t , X t , v t d W t + π ( t ) d H ( t ) , X 0 = x ,
where
f 1 t , X t , v t = δ t + 1 2 α 2 t X t + α t a t + b t + v t ,
and
σ 1 t , X t , v t = X t α t + a t .
Herein, the process v · stands for the control variable. We require that process v · is adapted, with càdlàg trajectories and taking values in U , a non-empty convex subset of R , such that the fourth-power condition
E 0 T v t 4 d t <
holds. Furthermore, we assume there exists a positive constant c 0 such that the SDE (3) has a unique solution satisfying the terminal constraint
E X T = c 0 .
This last equality means that the insurance firm is looking for some regulatory requirement c 0 described by the average value of its cash balance process at the terminal time T.
Herein, a control variable is said to be an admissible control if and only if it satisfies all the above conditions.
Let us point out that the aim of the policy maker is to simultaneously minimize the deviation between the firm’s cash-balance process and its dynamic benchmark over the set of all admissible controls, which will be denoted by U , the cost of the premium policy over the whole time interval 0 , T , and the terminal variance of the cash-balance process under some given constraint. Therefore, it is quite natural that the cost functional takes the following form:
J v · = E 0 T e β t g 1 t , X t , v t d t + e β T φ 1 X T ,
with
g 1 t , x , v = 1 2 R t x A t 2 + N t v 2 ( t ) ,
and
φ 1 x = 1 2 M x c 0 2 .
Here, β 0 is a discounting factor, A t is some dynamic pre-set target, representing the dynamic benchmark of X , processes R · , N · and the positive constant M are the weighting factors which make the cost functional (5) more general and flexible to control the preference of the policy-maker. Furthermore, we suppose that process A converges to c 0 as t goes to T .
Now, we can formulate the firm’s optimal premium problem as
Problem A:
To find v ^ U such that v ^ minimizes the cost function (5) subject to (3) and the state constraint (4).
To deal with this problem, we need to impose the following assumption on the previous coefficients:
Assumption 1.
Functions R · 0 , N · > 0 , N 1 · and A · are all deterministic and bounded on the time interval 0 , T .
Let us now reformulate the above control problem (3)–(5) as a generalized stochastic recursive optimal control problem with state constraint by introducing the following backward stochastic differential equation (BSDE), in term of processes Y ( · ) , Z ( · ) and K ( · ) = { K i ( · ) } i = 1 :
d Y t = g 1 t , X t , v t β Y t d t Z t d W t K ( t ) d H ( t ) , Y T = φ 1 X T .
Notice that (3) together with (8) form a semi-coupled Forward-Backward Stochastic Differential Equation (FBSDE) driven by both the Teugels martingales and an independent Brownian motion. Lemma 2.1 in Meng and Tang (2009) shows that under (H1), (3) admits a unique solution X t for each v · U . Moreover, we have
sup 0 t T E X 4 t < + .
The proof of this estimate is placed in Lemma 1 Section 3.
As a consequence, the terminal condition of (8) is square-integrable. Then, for the foregoing v t and X t , thanks to Theorem 3.1 in Bahlali et al. (2003), the Backward Stochastic Differential Equation (BSDE) (8) admits a unique solution Y · , Z · , K · under (H1). That is, for any v · U , the semi-coupled FBSDE consisting of (3) and (8) admits a unique solution:
X · , Y · , Z · , K · .
Obviously, by using the dual technique to the BSDE (8) one can get Y 0 = J v · . Then, we can reformulate Problem A in the following way:
Problem B:
To find v ^ U such that
J v ^ · = E ( Y ( 0 ) )
subject to (3), (4) and (8).
In the next section, we are going to prove a sufficient condition of optimality for the above problem, but in a more general form by assuming that coefficients are not necessarily linear with respect to the state variables.

3. Sufficient Condition of Optimality

3.1. Problem Formulation

Motivated by the above optimal premium policy of an insurance firm, we are now going to focus on the following control problem where the state process is described by the following controlled stochastic differential equation driven by both a Brownian motion and a family of Teugels martingales,
d X t = f t , X t , v t d t + σ t , X t , v t d W t + i = 1 π i t , X t , v t d H i t , X 0 = x ,
and the general stochastic differential utility is given by the BSDE
d Y t = g t , X t , v t d t Z t d W t i = 1 K i t d H i t , Y T = φ X T .
where X t : = X t , Y t , Z t , K t , and
f : 0 , T × Ω × R × U R , σ : 0 , T × Ω × R × U R , π : 0 , T × Ω × R × U P 2 R , g : 0 , T × Ω × R × R × R × P 2 R × U R , φ : 0 , T × Ω × R R ,
are progressively measurable functions.
In the sequel, for notational simplicity, we shall use the shorthand notation
π t , X t , v t d H t and K t d H t
instead of
i = 1 π i t , X t , v t d H i t and i = 1 K i t d H i t
respectively, where K t = K i t i = 1 and π t = π i t i = 1 .
Let us now introduce the following basic assumption on coefficients that will be needed in the sequel:
Assumption 2.
(i)
The function g is F progressively measurable for all y , z , k R × R × P 2 R and for any v · U ,
g t , X t , v t = G t , Y t , Z t , K t + R t X t A t 2 + N t v 2 t ,
with
E 0 T g s , 0 , 0 , 0 , 0 , v t 2 d s < .
(ii)
For every x R , φ L 2 , f , σ and π are F progressively measurable and for any v · U ,
E 0 T f s , 0 , v t 4 + σ s , 0 , v t 4 + π s , 0 , v t P 2 R 4 d s < .
(iii)
Functions G, f, σ and π are continuously differentiable with bounded derivatives with respect to x , y , z , k and v.
Note that Lemma 2.1 in Meng and Tang (2009) shows that under the Assumption (H2) the SDE (10) admits a unique solution that belongs to S F 2 0 , T , R . On the other hand, since the function g is uniformly Lipschitz with respect to y, z and k, by using the assumptions (H2), one can easily check that the BSDE (11) satisfies all the conditions in Theorem 3.1 in Bahlali et al. (2003), and hence, it has a unique solution that belongs to S F 2 0 , T , R × L F 2 0 , T , R × l F 2 0 , T , R .
The following key lemma gives an L 4 estimate of the solution of the state equation under Assumption (H2). Since Assumption (H1) implies Assumption (H2), the result is still valid under Assumption (H1).
Lemma 1.
Assume that X · is a solution of the state Equation (10). Then, under Assumption (H2), the following estimate holds true:
sup 0 t T E X 4 t < + .
Proof. 
First, we introduce the stopping times
R k = inf t : X t k ,
assuming k great enough to include X ( 0 ) = x in [ k , k ] . Then, applying Itô’s formula to X ( t R k ) 4 we obtain
X 4 t R k = X 4 0 + 4 0 t R k X 3 s f s , X s , v ( s ) d s + 4 0 t R k X 3 s σ s , X s , v ( s ) d W s + 4 0 t R k X 3 s π s , X s , v ( s ) d H s + 6 0 t R k X 2 s σ 2 s , X s , v ( s ) d s + 6 0 t R k X s 2 | | π s , X s , v ( s ) | | 2 d s + 6 i , j = 1 0 t R k X 2 s ( π i π j ) s , X s , v ( s ) d M i , j ( t ) + I ( t R k )
where
I ( t ) : = 0 < s t X 4 s X 4 s 4 X 3 s Δ X s 12 X 2 s Δ X s 2 .
The stochastic integrals with respect to W , H and M i , j in the right hand side of the above equality are in fact martingales with zero expectation. Then, taking expectations and later absolute values, one finds
E X t R k 4 E X 0 4 + 4 E 0 t R k X s 3 f s , X s , v ( s ) d s + 6 E 0 t R k X s 2 σ s , X s , v ( s ) 2 d s + 6 E 0 t R k X s 2 | | π s , X s , v ( s ) | | 2 d s + E ( | I ( t R k ) | ) .
In view of Young’s inequality, one can easily check that
E X t R k 4 E X 0 4 + 9 E 0 t R k X s 4 d s + E ( | I ( t R k ) | ) + E 0 t R k f s , X s , v ( s ) 4 d s + 3 E 0 t R k σ s , X s , v ( s ) 4 d s + 3 E 0 t R k | | π s , X s , v ( s ) | | 4 d s .
From now on, C will denote a constant that changes from line to line.
Using Lipschitz condition in (H2), we have
E X t R k 4 E X 0 4 + E I ( t R k ) + C E 0 t R k X s 4 d s + C E 0 t R k f s , 0 , v ( s ) 4 + σ s , 0 , v ( s ) 4 + | | π s , 0 , v ( s ) | | 4 d s .
We infer, using Gronwall’s inequality, that
E X t R k 4 e C ( t R k ) E X 0 4 + E I ( t R k ) + C E 0 t R k f s , 0 , v s 4 + σ s , 0 , v s 4 + | | π s , 0 , v s | | 4 d s .
This implies, using Fatou’s lemma and taking the supremum over the interval 0 , T , that
sup 0 t T E X t 4 C E X 0 4 + E I ( T ) + E 0 T f t , 0 , v t 4 + σ t , 0 , v t 4 + | | π t , 0 , v t | | 4 d s ,
Now, we turn out to estimate E ( | I ( T ) | ) . For each k , it is clear that the stopped process X ˜ = X 1 0 , R k is bounded by k, and it is a semi-martingale as a product of two semi-martingales.
Note that the estimate
0 < s t X ˜ 4 s X ˜ 4 s 4 X ˜ 3 s Δ X ˜ s 12 X ˜ 2 s Δ X ˜ s 2
C X ˜ , X ˜ t
holds true P a . s . This is a consequence of the fact that X ˜ · takes values in the compact interval k , k and the fact that for h x = x 4 , it is easy to see that
h x h y y x h x y x 2 h x C y x 2 .
Thus,
0 < s T | X ˜ 4 s X ˜ 4 s 4 X ˜ 3 s Δ X s 12 X ˜ 2 s Δ X ˜ s 2 |
C 0 < s T Δ X ˜ s 2 C X ˜ , X ˜ T < , P a . s .
Since the last inequality is valid for X ˜ , for each k, it also remains valid for X. Combining (14) with (13), we get the desired result and this completes the proof. □

3.2. Verification Theorem

In this subsection, we study Problem B with a more general state process. We establish a sufficient stochastic maximum principle for stochastic control of forward-backward SDEs driven by Brownian motion and Teugels martingales where the control domain is assumed to be convex.
First, by combining (10) and (11), we get the following controlled semi-coupled FBSDE driven by Brownian motion and Teugels martingales:
X t = X 0 + 0 t f s , X s , v s d s + 0 t σ s , X s , v s d W s + 0 t π s , X s , v s d H s , Y t = φ X T + t T g s , X s , v s d s t T Z s d W s t T K s d H s ,
We introduce the following cost functional:
J v · = E 0 T g s , X s , v s d s + φ X T ,
The optimal control problem is to minimize the cost functional J · over the set of all admissible controls.
To deal with the above control problem, we first define the Hamiltonian function
H : 0 , T × Ω × R × R × R × P 2 R × U × R × R × R × l 2 R R
by
H t , x , y , z , k , v , p , q , λ , ρ = g t , x , y , z , k , v λ + f t , x , v p + σ t , x , v q + π t , x , v ρ .
We will use the shorthand notation
H x t = H x t , x , y , z , k , v , p , q , λ , ρ ,
and similarly, with H y t , H z t , H k t and H v t . Similar notations are used for l t where l = f , σ , π and g .
The adjoint equations are described by the following stochastic Hamiltonian system:
d λ t = H y t d t + H z t d W t + H k t d H t , d p t = H x t d t + q t d W t + ρ t d H t , λ 0 = 1 , p T = λ T φ X T .
In addition to (H2), we need the following assumption:
Assumption 3.
If v ^ · U and
X ^ · , Y ^ · , Z ^ · , K ^ · , p ^ · , q ^ · , λ ^ · , ρ ^ ·
are the corresponding solutions of (10), (11) and (18), respectively, we assume
(i)
Functions x φ x and
t , x , y , z , k , v H t , x , y , z , k , v , p ^ t , q ^ t , λ ^ t , ρ ^ t
are convex.
(ii)
Function H satisfies
H t , X ^ t , Y ^ t , Z ^ t , K ^ t , v t , p ^ t , q ^ t , λ ^ t , ρ ^ t H t , X ^ t , Y ^ t , Z ^ t , K ^ t , v ^ t , p ^ t , q ^ t , λ ^ t , ρ ^ t
for any v · U and for almost all t , w 0 , T × Ω .
Then, we have the following sufficient condition for an optimal control of Problem B.
Theorem 1
(Sufficient maximum principle). Let v ^ · U with the corresponding solutions
X ^ · , Y ^ · , Z ^ · , K ^ · , p ^ · , q ^ · , λ ^ · , ρ ^ ·
of (10), (11) and (18) respectively. Assume (H3). Then, v ^ · is an optimal control for Problem B.
Proof. 
Let X v · , Y v · , Z v · , K v · , v · be an admissible solution of (15). It follows from the definition of the cost functional (16) that
J v ^ · J v · = E Y ^ 0 Y v 0 .
From the forward component of SDE (18), the right hand side of the above equality can be rewritten
E Y ^ 0 Y v 0 = E Y ^ 0 Y v 0 λ ^ 0 .
Applying Itô’s formula to ( Y ^ t Y v t ) λ ^ t from t = 0 to t = T and using the fact that H i , H j t = δ i j t and H i , H j H i , H j is a martingale, we obtain
E Y ^ 0 Y v 0 λ ^ 0 = E φ X ^ T φ X v T λ ^ T E 0 T Y ^ t Y v t d λ ^ t E 0 T λ ^ t d Y ^ t Y v t E 0 T H z t Z ^ t Z v t d t E 0 T H k t K ^ t K v t d t .
Since φ is convex, one can get
E φ X ^ T φ X v T λ ^ T E X ^ T X v T φ X v T λ ^ T
We remark that p ^ T = φ X ^ T λ ^ T , then
E Y ^ 0 Y v 0 λ ^ 0 E X ^ T X v T p ^ T E 0 T H y t Y ^ t Y v t d t + E 0 T ( λ ^ t g ^ t g v t ) d t E 0 T H z t Z ^ t Z v t d t E [ 0 T H k t K ^ t K v t d t ] .
On the other hand, Itô’s formula applied to X ^ t X t p ^ t gives us
E X ^ T X v T p ^ T = E 0 T X ^ t X v t d p ^ t + E 0 T p ^ t d X ^ t X v t + E 0 T σ ^ t σ t q ^ t d t + E 0 T π ^ t π t ρ ^ t d t .
Substituting (21) into (20), it follows immediately that
E Y ^ 0 Y v 0 E 0 T H x t X ^ t X v t d t + 0 T f ^ t f t p ^ t + σ ^ t σ t q ^ t + π ^ t π t ρ ^ t + λ ^ t g ^ t g t H y t Y ^ t Y v t H z t Z ^ t Z v t H k t K ^ t K v t d t .
Then
E Y ^ 0 Y v 0 E 0 T H ^ t H t H x t X ^ t X v t H y t Y ^ t Y v t H z t Z ^ t Z v t H k t K ^ t K v t d t .
By virtue of the convexity property of the Hamiltonian H with respect to x , y , z , k for almost all t , w 0 , T × Ω , one can get
E Y ^ 0 Y v 0 E 0 T H x t X ^ t X v t + H y t Y ^ t Y v t + H z t Z ^ t Z v t + H k t K ^ t K v t + H v t v ^ t v t H x t X ^ t X v t H y t Y ^ t Y v t H z t Z ^ t Z v t H k K ^ t K v t d t = E 0 T H v t v ^ t v t d t .
By invoking the necessary condition of optimality (19), we conclude that
E Y ^ 0 Y v 0 0 ,
which implies,
J v ^ · J v · 0 , v · U .
This finishes the proof. □

4. Applications

4.1. Optimal Premium Problem

In this subsection, firstly, we use the Lagrangian method to treat the terminal state constraint, and after that, we apply the sufficient maximum principle to deal with the resulting unconstrained optimization problem. Throughout this subsection, we assume that α = 0 in (2), which means that the insurance firm only invests in a money account with compounded interest rate δ t , and hence, Δ t = 0 t δ s d s . We recall that this kind of problem in the continuous Brownian case was solved in Huang et al. (2010).
Then, SDE (3) becomes
d X t = δ t X t + b t + v t d t + a t d W t + π t d H t , X 0 = x .
By using the Lagrangian multiplier method, the cost function (5) becomes
J v ( · ) = E 0 T e β t g 1 t , X t , v t d t + e β T φ 1 X T + θ X T c 0
where θ , the Lagrange multiplier, is some constant to be determined. Then, we can reformulate Problem A, as the following one: we look forward to find v ^ U , which minimizes the cost function (23) subject to (22) and (4).
For this end, let us firstly define the Hamiltonian function
H : 0 , T × R × U × R × R × P 2 R R ,
by
H t , X , v , p , q , ρ = δ t X + b t + v p + a t q + π t ρ + 1 2 e β t R t X A t 2 + N t v 2 .
Then, the adjoint equation can be rewritten in a Hamiltonian form as
d p t = δ t p t + X t A t R t e β t d t q t d W t ρ t d H t , p T = θ + M e β T X T c 0 ,
where X · stands for the optimal process under the optimal premium policy v · and satisfies the Equation (22).
Let ϕ and ψ be the solutions of
ϕ t + 2 δ t ϕ t N 1 t e β t ϕ 2 t + R t e β t = 0 , ϕ T = M e β T ,
and
ψ t + δ t N 1 t e β t ϕ t ψ t + b t ϕ t A t R t e β t = 0 , ψ T = θ c 0 M e β T .
We now are in position to state and prove the main result of this subsection.
Theorem 2.
Under Assumption (H1), the optimal premium policy is given by
v ^ t = N 1 t e β t ϕ t X t + ψ t ,
where X · satisfies (22) and ϕ and ψ are the solutions of (26) and (27), respectively. Moreover, the optimal cost functional is given by
J v ^ = 1 2 0 T e β t R t A 2 t d t + M e β T c 0 2 + 1 2 ϕ 0 x 2 + ψ 0 x c 0 θ + 1 2 0 T ϕ t a 2 t + ϕ t π 2 t + ψ t 2 b t N 1 t e β t ψ t d t .
Proof. 
We shall divide the proof into several steps.
Step 1: We start by proving the existence of an optimal premium policy. Since for each t, (5) is quadratic with respect to X t , v t , X T , and the weight N t of v 2 t is larger than 0, then there exists an optimal premium policy v ^ · , which solves Problem A. Indeed, from Theorem 6.1 in Meng and Tang (2009), the optimal premium policy v ^ · satisfies
0 = H v t t , X t , v t , p t , q t , ρ t = v t δ t X t + b t + v t p t + a t q t + π t ρ t + 1 2 e β t R t X t A t 2 + N t v 2 t = p t + e β t N t v t ,
and is given by
v ^ t = N 1 t e β t p t .
Substituting (29) into (22) and combining it with (25), we obtain the generalized Hamiltonian system
d X t = δ t X t N 1 t e β t p t + b t d t + a t d W t + π t d H t , d p t = δ t p t + X t A t R t e β t d t q t d W t ρ t d H t , X 0 = x , p T = θ + M e β T X T c 0 ,
which is a coupled FBSDE system. Using a similar argument to the ones used in Baghery et al. (2014), we can verify that FBSDE (30) admits a unique solution under (H1). Note that the optimal premium policy exists and is unique. Furthermore the relation (29) implies that it is linear with respect to p .
Step 2: In this step, we are going to prove that the optimal premium policy v ^ · is in fact a linear feedback of the optimal process X · . First of all, according to the terminal condition of p · in (30), it is quite natural to suggest that
p t = ϕ t X t + ψ t ,
with ϕ T = M e β T and ψ T = θ c 0 M e β T . Then, we apply Itô’s formula to p · to obtain
d p t = ϕ t + δ t ϕ t X t + ψ t + b t N 1 t e β t p t ϕ t d t + ϕ t δ t d W t + ϕ t π t d H t .
Comparing their generator terms with those of the BSDE in (30) we get (26) and (27).
It is well known that (26) is a standard Riccati differential equation that admits a unique solution under (H1), so does (27) and
ψ t = t T b s ϕ s A s R s e β s Λ t s d s + θ c 0 M e β T Λ t T ,
where
Λ t s = exp { t s δ r N 1 r e β r ϕ r d r } .
Step 3: In this step, we are going to find the value of θ . Combining (29) and (31), we get
d X t = δ t N 1 t e β t ϕ t X t + b t N 1 t e β t ψ t d t + a t d W t + π t d H t , X 0 = x .
Define X ˜ t = E X t , then
X ˜ t = δ t N 1 t e β t ϕ t X ˜ t + b t N 1 t e β t ψ t , X ^ 0 = x
Solving (33) and keeping in mind the terminal constraint (4), we easily derive
c 0 = x Λ 0 T + 0 T b t N 1 t e β t ψ t Λ t T d t .
Note that ϕ in (26) does not depend on θ and ψ in (27) is linear with respect to θ . Inserting (27) into (34), we get the equation
θ = 0 T b t N 1 t F t Q t Λ t T d t + x Λ 0 T c 0 0 T N 1 t e β t Λ t 2 T d t
where
F t = t T e β t b s ϕ s A s R s e β t s Λ t s d s
and
Q t = c 0 M exp { t T δ s N 1 s e β s ϕ s β d s } .
Step 4: We now proceed to determine the optimal cost functional. Substituting (29) into (5) we get
J v ^ · = 1 2 E 0 T e β t R t + N 1 t e 2 β t ϕ 2 t X 2 t 2 A t R t N 1 t e 2 β t ϕ t ψ t X t + N 1 t e 2 β t ψ 2 t d t + M e β T X 2 T 2 c 0 M e β T X T + 1 2 0 T e β t R t A 2 t d t + M e β T c 0 2 .
From Itô’s formula,
d ϕ t X 2 t + 2 ψ t X t = ϕ t a 2 t + 2 ψ t b t N 1 t e β t ψ t + 2 e β t A t R t N 1 t e 2 β t ϕ t ψ t X t e β t R t + N 1 t e 2 β t ϕ 2 t X 2 t d t + 2 σ t ϕ t X t + ψ t d W t + 2 π t ϕ t X t + ψ t d H t + ϕ t i , j π i t π j t d H i , H j t .
Integrating from 0 to T, taking expectations on both sides of the above equality and using the fact that H i , H j t H i , H j t is an F martingale and H i , H j t = δ i j t , we get
E ϕ T X 2 T + 2 ψ T X T = 0 T ϕ t a 2 t + π 2 t + 2 ψ t b t N 1 t e β t ψ t d t + ϕ 0 x 2 + 2 ψ 0 x + 2 E 0 T e β t A t R t N 1 t e 2 β t ϕ t ψ t X t d t E 0 T e β t R t + N 1 t e 2 β t ϕ 2 t X 2 t d t .
On the other hand, it follows from (26) and (27) that
E ϕ T X 2 T + 2 ψ T X T = E M e β T X 2 T 2 c 0 M e β T X T + 2 θ X T .
Inserting the above two equalities into (35) we obtain the optimal cost functional (28). This gives the desired result. □
Specifically, in case we rule out the terminal constraint (4), that is, θ = 0 , Theorem 2 solves Problem A without the constraint (4).
Note also that ϕ , ψ and θ do not depend on a and π . This leads to the fact that the optimal cost function J u increases with respect to σ and π (describing the liability risk). This implies that the more uncertainty in the liability process, the higher the operation costs of our optimal premium policy.

4.2. Optimal Premium Problem under Stochastic Interest Rate

In the current subsection, we want to discuss the optimal premium problem (Problem B), in the case where the interest rate is allowed to be stochastic. More precisely, we shall consider two different cases. In the first one, we assume that the payment function and the stochastic interest rate are given by the same Brownian motion and in the second case, we assume that they are given by different and independent Brownian motions.

4.2.1. First Case

In this case, we proceed to solve Problem B, assuming that the utility function is that of (9) and the cash-balance process satisfies (3), where the stochastic interest rate is given by (2). We also point out that we solve a mean variance problem in this section as a particular case of Problem B.
Since the coefficients in (3), (6) and (7) are linear, we shall assume
d X t = a 1 t X t + a 2 t v t + a 3 t d t + b 1 t X t + b 2 t v t + b 3 t d W t + π t d H t , X 0 = x ,
and
d Y t = 1 2 c 1 t X 2 t + c 2 t X t + c 3 t v 2 t + c 4 t β Y t d t Z t d W t K t d H t Y T = 1 2 X T a 2
where a is real constant and a 1 , a 2 , a 3 , b 1 , b 2 0 , b 3 , c 1 , c 2 , c 3 , c 4 are deterministic functions satisfying some properties. Note that particular choices of the coefficients of (36) and (37) generate the SDE (10) and the BSDE (11).
Here, the Hamiltonian (17) gets the form
H = 1 2 c 1 t X 2 t + c 2 t X t + c 3 t v 2 t + c 4 t β Y t λ t + a 1 t X t + a 2 t v t + a 3 t p t + b 1 t X t + b 2 t v t + b 3 t q t + π t ρ t ,
λ t = e β t and p · , q · , ρ · satisfies the adjoint backward equation
d p t = c 1 t X t λ t + a 1 t p t + b 1 t q t + 1 2 c 2 t λ t d t + q t d W t + ρ t d H t , p T = θ + λ t φ X T .
By using the sufficient condition (19) given in Section 3, we get
λ t c 3 t v t + a 2 t p t + b 2 t q t = 0 .
To solve (38), we try a process p of the form
p t = ϕ t X t + ψ t ,
where ϕ and ψ are derivable functions with continuous derivatives.
Applying Itô’s formula to (40) and using (36), we get
d p t = ϕ t a 1 t + ϕ t X t + ϕ t a 2 t v t + ϕ t a 3 t + ψ t d t + ϕ t b 1 t X t + b 2 t v t + b 3 t d W t + ϕ t π t d H t
Comparing with (38), we obtain
ρ t = ϕ t π t , q t = ϕ t b 1 t X t + b 2 t v t + b 3 t
and
c 1 t X t λ t + a 1 t p t + b 1 t q t + 1 2 c 2 t λ t = ϕ t a 1 t + ϕ t X t + ϕ t a 2 t v t + ϕ t a 3 t + ψ t .
Substituting (42) into (39), we obtain
v ^ t = a 2 t + b 1 t b 2 t ϕ t X ^ t + a 2 t ψ t + b 2 t b 3 t ϕ t λ t c 3 t + b 2 2 t ϕ t .
On the other hand, from (43) we have
v ^ t = c 1 t λ t + 2 a 1 t ϕ t + b 1 2 t ϕ t + ϕ t X ^ t a 2 t + b 1 t b 2 t ϕ t a 1 t ψ t + b 1 t b 3 t ϕ t + 1 2 c 2 t λ t + a 3 t ϕ t + ψ t a 2 t + b 1 t b 2 t ϕ t .
Combining (44) and (45), we obtain
ϕ 2 t a 2 t + b 1 t b 2 t 2 b 2 2 t 2 a 1 t + b 1 2 t ϕ t ϕ t b 2 2 t c 3 t 2 a 1 t + b 1 2 t + c 1 t b 2 2 t ϕ t λ t c 3 t ϕ t λ t = c 1 t c 3 t λ 2 t , ϕ T = λ T .
and
ψ t ϕ t a 2 t a 2 t + b 1 t b 2 t a 1 t b 2 2 t c 3 t a 1 t λ t ψ t c 3 t λ t + ϕ t b 2 2 t = ϕ 2 t a 2 t + b 1 t b 2 t b 2 t b 3 t b 2 2 t b 1 t b 3 t + a 3 t + c 3 t b 1 t b 3 t + a 3 t ϕ t λ t + 1 2 c 2 t b 2 2 t ϕ t + c 3 t c 2 t λ t λ t , ψ T = θ λ T c 0 .
Let us emphasize that, under the boundedness of its coefficients, the first order differential Equation (46) admits a unique solution. Moreover, a simple computation shows the (47) has an explicit solution, which is given by
ψ t = θ λ T c 0 exp Λ t T exp Λ t T t T e Λ t s G t s d s ,
where
Λ t T = t T ϕ s a 2 s + b 1 s b 2 s a 2 s a 1 s b 2 2 s c 3 s λ s + b 2 2 s ϕ s d s t T c 3 s a 1 s λ s c 3 s λ s + ϕ s b 2 2 s d s ,
and
G t s = t T ϕ 2 s b 2 s b 3 s a 2 s + b 1 s b 2 s c 3 s λ s + ϕ s b 2 2 s d s + t T ϕ 2 s b 2 2 s b 1 s b 3 s + a 3 s c 3 s λ s + ϕ s b 2 2 s d s + t T c 3 s b 1 s b 3 s + a 3 s λ s ϕ s c 3 s λ s + ϕ s b 2 2 s d s . + t T 1 2 c 2 s b 2 2 s ϕ s + c 3 s c 2 s λ s λ s c 3 s λ s + ϕ s b 2 2 s d s
Thus, we have the following theorem:
Theorem 3.
Let X v t be the cash balance satisfying (36). Then, the optimal premium policy of Problem A is
v ^ t = a 2 t + b 1 t b 2 t ϕ t X ^ t + a 2 t ψ t + b 2 t b 3 t ϕ t c 3 t λ t + b 2 2 t ϕ t ,
where λ ( t ) = e β t and ϕ and ψ are the solutions of (46) and (47), respectively.
This is clearly the unique solution under Assumption (H1). Notice that we can get explicit solutions in some particular cases. For example, if g 1 = 0 we can see the problem as a mean-variance one, where our objective is to find v t such that it minimizes
V X T = E X T E X T 2 ,
under the terminal constraint condition (4).
By the Lagrangian multiplier method, the problem can be reduced to minimize the following equivalent problem:
J v · = E 1 2 X v T c 0 2 + θ X v T c 0 .
In this case, the Hamiltonian (17) gets the form
H = a 1 t X t + a 2 t v t + a 3 t p t + b 1 t X t + b 2 t v t + b 3 t q t + π t ρ t .
Hence, the adjoint equation takes the following form
d p t = a 1 t p t + b 1 t q t d t + q t d W t + ρ t d H t p T = θ + X T c 0 .
Note that in this case λ t = 1 , t 0 , T .
In the following corollary, we solve Problem (48) assuming that X v t satisfies (36).
Corollary 1.
Let X v t be the cash balance satisfying (36). Then, the optimal premium policy of Problem (48) is
v ^ t = a 2 t + b 1 t b 2 t ϕ t X ^ t + a 2 t ψ t + b 2 t b 3 t ϕ t b 2 2 t ϕ t ,
where ϕ and ψ are given, respectively, by
ϕ t = exp { t T a 2 s + b 1 s b 2 s 2 2 a 1 s + b 1 2 s b 2 2 s b 2 2 s d s }
and
ψ t = θ c 0 exp Λ t T exp Λ t T t T exp Λ t s ϕ s a 3 s b 2 s a 2 s b 3 s b 2 s d s ,
with
Λ t T = t T a 2 s a 2 s + b 1 s b 2 s a 1 s b 2 2 s b 2 2 s d s .

4.2.2. Second Case

In this case, we assume the liability of the surplus process and the interest rate are given by two different and independent Brownian motions. Let W ¯ : = W ¯ t t 0 be another Brownian motion defined in Ω , F , P and independent of W and L . Assume
F t = σ W ( s ) , W ¯ ( s ) , L ( s ) , 0 s t N .
We consider that the interest rate is given by the stochastic differential equation
d Δ t = δ t d t + α t d W ¯ t , t 0 , T , Δ 0 = 0 .
Itô formula applied to the process X · described by (1) in Section 2, leads to
d X t = X t δ t + 1 2 α 2 t + b t + v t d t + X t α t d W ¯ t + a t d W t + π t d H t , X 0 = x .
Define a process V such that
X t α t d W ¯ t + a t d W t = X 2 t α 2 t + σ 2 t 1 2 d V t ,
that is,
V t = 0 t X 2 s α 2 s + a 2 s 1 2 X s α s d W ¯ s + a s d W s .
Clearly, the process V is a continuous martingale with quadratic variation V t = t , and so, it must be a standard Brownian motion, see for example Theorem 6.1 in Chung and Williams (1990). Thus, we can write Equation (50) by
d X t = X t δ t + 1 2 α 2 t + b t + v t d t + X 2 t α 2 t + a 2 t 1 2 d V t + π t d H t , X 0 = x .
We shall derive the solution of Problem A. Here, the Hamiltonian (17) gets the form
H = X t δ t + 1 2 α 2 t + b t + v t p + X 2 t α 2 t + a 2 t 1 2 q + π t ρ + 1 2 e β t R t X t A t 2 + N t v 2 t .
The adjoint process p · , q · , ρ · satisfies the BSDE
d p t = a 4 t p t + a 5 t + R t X t A t e β t d t q t d V t ρ t d H t , p T = θ + M e β T X T c 0 .
where a 4 t = δ t + 1 2 α 2 t and a 5 t = α 2 t X t q t X 2 t α 2 t + a 2 t . Notice that in this case, λ t = 1 ,   t 0 , T .
By the same technique used in Section 4.1, we infer that the optimal premium policy is given by
v ^ t = N 1 t e β t ϕ t X t + ψ t ,
where X · satisfies (51) and ϕ and ψ are, respectively, the solutions of the equations
ϕ t + 2 δ t + α 2 t ϕ t N 1 t e β t ϕ 2 t + R t e β t = 0 ϕ T = M e β T ,
and
ψ t + δ t + 1 2 α 2 t N 1 t e β t ϕ t ψ t + b t ϕ t A t R t e β t = 0 , ψ T = θ c 0 M e β T .
Moreover, the optimal cost functional is given by
J v ^ · = 1 2 ( 0 T e β t R t A 2 t d t + M e β T c 0 2 ) + 1 2 ϕ 0 x 2 + ψ 0 x c 0 θ + 1 2 0 T ϕ t a 2 t + ϕ t α 2 t + ψ t 2 b t N 1 t e β t ψ t d t .

Author Contributions

All authors have equally contributed to the paper in all aspects. All authors have read and agreed to the published version of the manuscript.

Funding

Research of Josep Vives is partially financed by the Spanish MICIN grant PID2020-118339GB-100 (2021-2024).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Asmussen, Søren, Bent Jesper Christensen, and Julie Thøgersen. 2019. Nash equilibrium premium strategies for push–pull competition in a frictional non-life insurance market. Insurance: Mathematics and Economics 87: 92–100. [Google Scholar] [CrossRef]
  2. Baghery, Fouzia, Nabil Khelfallah, Brahim Mezerdi, and Isabelle Turpin. 2014. Fully coupled forward backward stochastic differential equations driven by Lévy processes and application to differential games. Random Operators and Stochastic Equations 22: 151–61. [Google Scholar] [CrossRef]
  3. Bahlali, Khaled, M’hamed Eddahbi, and El Hassan Essaky. 2003. BSDEs associated with Lévy processes and application to PDIE. Journal of Applied Mathematics and Stochastic Analysis 16: 1–17. [Google Scholar] [CrossRef] [Green Version]
  4. Cairns, Andrew. 2000. Some notes on the dynamics and optimal control of stochastic pension fund models in continuous time. Astin Bulletin 30: 19–55. [Google Scholar] [CrossRef] [Green Version]
  5. Chung, Kai Lai, and Ruth J. Williams. 1990. Introduction to Stochastic Integration, 2nd ed. Basel: Birkhaüser. [Google Scholar]
  6. Hipp, Christian, and Michael Vogt. 2003. Optimal dynamic XL-reinsurance. Astin Bulletin 33: 193–207. [Google Scholar] [CrossRef] [Green Version]
  7. Højgaard, Bjarne, and Michael Taksar. 1998. Optimal proportional reinsurance policies for diffusion models. Scandinavian Actuarial Journal 2: 166–80. [Google Scholar] [CrossRef]
  8. Huang, Jianhui, Guangchen Wang, and Zhen Wu. 2010. Optimal premium policy of an insurance firm: Full and partial information. Insurance: Mathematics and Economics 47: 208–15. [Google Scholar] [CrossRef]
  9. Josa-Fombellida, Ricardo, and Juan Pablo Rincón-Zapatero. 2008. Mean-variance portfolio and contribution selection in stochastic pension funding. European Journal of Operational Research 187: 120–37. [Google Scholar] [CrossRef] [Green Version]
  10. Kalfin, Sukono Firman, Sudrajat Supian, and Mustafa Mamat. 2021. Mitigation and models for determining premiums for natural disaster insurance due to excessive rainfall. Journal of Physics: Conference Series 1722: 012058. [Google Scholar] [CrossRef]
  11. Li, Danping, and Virginia R. Young. 2021. Bowley solution of a mean–variance game in insurance. Insurance: Mathematics and Economics 98: 35–43. [Google Scholar] [CrossRef]
  12. Lin, Jiaxin, Juliang Zhang, and Tai Chiu Edwin Cheng. 2020. Optimal pricing and return policy and the value of freight insurance for a retailer facing heterogeneous consumers with uncertain product values. International Journal of Production Economics 229: 107767. [Google Scholar] [CrossRef]
  13. Meng, QingXin, and MaoNing Tang. 2009. Necessary and sufficient conditions for optimal control of stochastic systems associated with Lévy processes. Science in China Series F: Information Sciences 52: 1982–92. [Google Scholar] [CrossRef]
  14. Moore, Kristen S., and Virginia R. Young. 2006. Optimal insurance in a continuous-time model. Insurance: Mathematics and Economics 39: 47–68. [Google Scholar] [CrossRef]
  15. Ngwira, Bernard, and Russell Gerrard. 2007. Stochastic pension fund control in the presence of Poisson jumps. Insurance: Mathematics and Economics 40: 283–92. [Google Scholar] [CrossRef]
  16. Nualart, David, and Wim Schoutens. 2000. Chaotic and predictable representations for Lévy processes. Stochastic Processes and their Applications 90: 109–22. [Google Scholar] [CrossRef] [Green Version]
  17. Xie, Shuxiang, Zhongfei Li, and Shouyang Wang. 2008. Continuous-time portfolio selection with liability: Mean–variance model and stochastic LQ approach. Insurance: Mathematics and Economics 42: 943–53. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guerdouh, D.; Khelfallah, N.; Vives, J. Optimal Control Strategies for the Premium Policy of an Insurance Firm with Jump Diffusion Assets and Stochastic Interest Rate. J. Risk Financial Manag. 2022, 15, 143. https://doi.org/10.3390/jrfm15030143

AMA Style

Guerdouh D, Khelfallah N, Vives J. Optimal Control Strategies for the Premium Policy of an Insurance Firm with Jump Diffusion Assets and Stochastic Interest Rate. Journal of Risk and Financial Management. 2022; 15(3):143. https://doi.org/10.3390/jrfm15030143

Chicago/Turabian Style

Guerdouh, Dalila, Nabil Khelfallah, and Josep Vives. 2022. "Optimal Control Strategies for the Premium Policy of an Insurance Firm with Jump Diffusion Assets and Stochastic Interest Rate" Journal of Risk and Financial Management 15, no. 3: 143. https://doi.org/10.3390/jrfm15030143

APA Style

Guerdouh, D., Khelfallah, N., & Vives, J. (2022). Optimal Control Strategies for the Premium Policy of an Insurance Firm with Jump Diffusion Assets and Stochastic Interest Rate. Journal of Risk and Financial Management, 15(3), 143. https://doi.org/10.3390/jrfm15030143

Article Metrics

Back to TopTop