# Explicit Solution of Reinsurance-Investment Problem for an Insurer with Dynamic Income under Vasicek Model

Authors: De-Lei Sheng
Journal: Advances in Mathematical Physics
 (2016)

---

## Abstract

Unlike traditionally used reserves models, this paper focuses on a reserve process with dynamic income to study the reinsurance-investment problem for an insurer under Vasicek stochastic interest rate model. The insurer’s dynamic income is given by the remainder after a dynamic reward budget being subtracted from the insurer’s net premium which is calculated according to expected premium principle. Applying stochastic control technique, a Hamilton-Jacobi-Bellman equation is established and the explicit solution is obtained under the objective of maximizing the insurer’s power utility of terminal wealth. Some economic interpretations of the obtained results are explained in detail. In addition, numerical analysis and several graphics are given to illustrate our results more meticulous.


---

## Body

## 1. Introduction

Existing literatures investigate the optimal reinsurance-investment problem only from the view of insurer but lack consideration to balance the interests of both insurance company and policy-holders. An effective way to balance both profits of insurer and policy-holders is reward budget. The reward budget is a dynamic item which can be subtracted from the actuarial premium income using the classic premium principle. As a consequence, insurer’s premium income is dynamic, which is given by an actuarial income under classic premium principle minus a dynamic reward budget item.Portfolio optimization is a very interesting problem, which attracted a lot of research work. Cao and Xu [1] derived the forms of proportional and excess-of-loss reinsurance contracts under the assumption that investment fund follows the logarithm-normal distribution. Cao and Zeng [2] studied the optimal reinsurance-investment problem in order to minimize the probability of ruin. However, only a handful of research works studied the portfolio optimization problem under the stochastic interest rates during the past twenty years. Grasselli [3] studies an investment problem for the HARA utility functions assuming the interest rates are described with the Cox-Ingersoll-Ross dynamics. Under the same stochastic interest rate model, Li and Wu [4] obtained an optimal policy in closed form in order to maximize a power utility for the optimal investment problem. When the interest rate is described with an ergodic Markov diffusion process, Pang [5] considered a portfolio optimization problem on an infinite time horizon by maximizing the infinite horizon expected discounted log utility of consumption. If the randomness of the interest rates and inflation are considered at the same time, Munk and Rubtsov [6] obtained a solution in closed form of stock-bond-cash portfolio problem for a risk- and ambiguity-averse investor. In addition, Korn and Kraft [7] considered just an investment problem with stochastic interest rates; however, they assumed that an investor can invest on several assets, such as savings account, stocks, and bonds, in order to maximize an insurer’s utility of terminal wealth without the usual Lipschitz assumptions. Similarly, Hainaut [8] also considered the financial market with a variety of assets, cash, stocks, and a zero coupon bond, and investigated the problem of dynamic asset allocation assuming the dynamics of the instantaneous short rates being driven by a Hull and White model. But all these previous research works did not consider the insurer having a dynamic income, they have ignored the insurers need to identify preferred customers for promotional campaigns or for customer favourable activities in a variety of forms. Such neglected problems directly lead to our study of the dynamic reward budget, which makes the income of insurance company become dynamic.This paper focuses on a reinsurance-investment problem with dynamic income for an insurer under Vasicek stochastic interest rate model. The dynamic income of an insurer is given by the remainder after a dynamic reward budget being subtracted from the insurer’s net premium which is calculated according to expected premium principle. Reward budget is an effective measure to attract more customers so that the total premium income may increase faster with timet. It is also an effective way to balance the profit of both insurance company and policy-holders. A Hamilton-Jacobi-Bellman equation is established and the explicit solution for the HJB equation is obtained by maximizing the insurer’s power utility of terminal wealth. At last, some numerical analysis and graphics are given to show the impacts of different parameters. Many interesting facts and laws hidden behind the formulas become very clear and more easy to understand.The subsequent parts of this paper are organized as follows. Section2 gives some models and a general framework will be used hereinafter. In Section 3, two Hamilton-Jacobi-Bellman equations are established for both cases of stochastic interest rate and of constant interest rate, respectively. The explicit solutions are given with the constant relative risk aversion function and strategies of two different cases are compared. In addition, the economic interpretations of the obtained results are given in the last subsection of this part. Section 4 provides some numerical analysis and graphics.


## 2. Formulation of the Model

Throughout this paper,(
Ω
,
F
,
P
,
{
F
t
}
0
≤
t
≤
T
) denotes a complete probability space satisfying the usual condition, where T
>
0 is a finite constant representing the investment time horizon; F
t stands for the information available until time t.


### 2.1. The Financial Market

We consider a financial market in which transaction amounts are so small that we can consider they have no influence on the prices, with the market being no arbitrage, frictionless, and traded continuously.Without loss of generality, we assume that the financial market is composed of two kinds of assets: cash and equities. For the sake of simplicity, we will only consider one equity asset which can indeed represent the index of a stock market.The instantaneous risk-free rater
(
t
) followed by the stochastic interest rate models of constant volatility, which is described by the Vasicek [9] model: (1)
d
r
t
=
β
-
α
r
t
d
t
+
δ
d
W
^
t
, where W
^
(
t
) is a standard Brownian motion and β
>
0, α
>
0, and δ
>
0 are constants.The priceS
0
(
t
) of the cash asset is given by(2)
d
S
0
t
S
0
t
=
r
t
d
t
,
S
0
0
=
1
.Here r
(
t
) is the stochastic interest rate.The stock market index priceS
(
t
) at time t
≥
0 can be described by the model:(3)
d
S
t
S
t
=
r
t
+
u
t
d
t
+
v
d
W
~
t
,
S
0
=
s
0
,where u
(
t
) is another positive real-valued function and the constant v
>
0 denotes the volatility rate of risk asset. ρ
2
∈
[
-
1,1
] is the correlation coefficient of two Brownian motions W
^
(
t
) and W
~
(
t
) and E
[
W
^
(
t
)
,
W
~
(
t
)
]
=
ρ
2
·
t.In model (3), u
(
t
)
>
0 is a natural hypothesis, which means that the investment on stock can make more money than investing on cash asset, which is consistent with the economic principle of high-risk high yield. In fact, only when u
(
t
)
>
0, more people are willing to invest on equities. For simplicity, we also can take the positive real-valued function u
(
t
) as a positive constant u
>
0, which does not bring substantial impact on the problem investigated in this paper.


### 2.2. The Surplus Process of Insurer

We use the classical Cramér-Lundberg model to describe the surplus process of insurer:(4)
d
R
t
=
c
d
t
-
d
∑
i
=
1
N
t
Y
i
,where N
(
t
) is a homogeneous Poisson process with intensity λ and the claim sizes {
Y
i
,
i
≥
1
} are independent and identically distributed positive random variables with E
[
Y
i
]
=
μ and E
[
Y
i
2
]
=
σ
2. According to the expected value premium principle, the pure premium rate of insurer is c
=
(
1
+
θ
)
λ
μ and θ
>
0 is the safety loading of insurer.The diffusion approximation of the compound Poisson processd
(
∑
i
=
1
N
(
t
)
Y
i
) is given by(5)
d
∑
i
=
1
N
t
Y
i
=
λ
μ
d
t
-
λ
σ
d
W
t
,where W
(
t
) is a standard Brownian motion, which is independent with W
~
(
t
). Moreover, we also assume E
[
W
(
t
)
,
W
^
(
t
)
]
=
ρ
1
·
t, ρ
1
∈
[
-
1,1
] denoting the correlation coefficient of W
(
t
) and W
^
(
t
).In order to avoid unbearable risk, the insurance company can buy proportional reinsurance to transfer part of the potential risk to reinsurance company. The insurer’s retention proportion isp
∈
[
0,1
], and the remaining 1
-
p will be borne by the reinsurer. η
>
0 is the safety loading of reinsurer:(6)
d
R
t
=
c
-
1
+
η
1
-
p
λ
μ
d
t
-
p
λ
μ
d
t
-
λ
σ
d
W
t
=
1
+
θ
λ
μ
-
1
+
η
1
-
p
λ
μ
-
p
λ
μ
d
t
+
p
λ
σ
d
W
t
.


### 2.3. The Insurer’s Wealth Process with Reward Budget

The proportion of wealth invested on stock asset is denoted byπ
(
t
) which is a progressively measurable real function, X
(
t
) with X
(
0
)
=
x
0 represents the wealth of the insurer, and then we have(7)
d
X
t
=
d
R
t
+
π
t
X
t
d
S
t
S
t
+
1
-
π
t
X
t
d
S
0
t
S
0
t
=
θ
-
η
λ
μ
+
η
p
λ
μ
+
π
t
u
t
+
r
t
X
t
d
t
+
p
λ
σ
d
W
t
+
v
π
t
X
t
d
W
~
t
.A notationd
=
(
η
-
θ
)
λ
μ
/
η
p
λ
μ
=
(
η
-
θ
)
/
η
p
=
(
1
-
θ
/
η
)
1
/
p and (
η
-
θ
)
λ
μ
=
d
η
p
λ
μ, where d is an important constant which indicates the benefit game on safety loading of insurer and reinsurer. The new notation d introduced here is interesting, existing for every combination of an insurance contract with a proportional reinsurance transaction. If η
=
θ, then d
=
0, with such proportional reinsurance transaction being more favourable for insurer. The actual reinsurance transaction must balance both benefits of insurer and reinsurer, so that 0
<
d
<
1 holds invariably.The reward budget at a fixed experience share is deducted from premium income at each timet:(8)
d
X
t
=
c
-
k
X
t
-
1
+
η
1
-
p
λ
μ
d
t
-
p
λ
μ
d
t
-
λ
σ
d
W
t
+
π
t
1
-
k
X
t
d
S
t
S
t
+
1
-
π
t
1
-
k
X
t
d
S
0
t
S
0
t
=
1
-
d
η
p
λ
μ
+
1
-
k
X
t
r
t
-
k
X
t
+
π
t
1
-
k
X
t
u
t
d
t
+
p
λ
σ
d
W
t
+
π
t
1
-
k
X
t
v
d
W
~
t
.


## 2.1. The Financial Market

We consider a financial market in which transaction amounts are so small that we can consider they have no influence on the prices, with the market being no arbitrage, frictionless, and traded continuously.Without loss of generality, we assume that the financial market is composed of two kinds of assets: cash and equities. For the sake of simplicity, we will only consider one equity asset which can indeed represent the index of a stock market.The instantaneous risk-free rater
(
t
) followed by the stochastic interest rate models of constant volatility, which is described by the Vasicek [9] model: (1)
d
r
t
=
β
-
α
r
t
d
t
+
δ
d
W
^
t
, where W
^
(
t
) is a standard Brownian motion and β
>
0, α
>
0, and δ
>
0 are constants.The priceS
0
(
t
) of the cash asset is given by(2)
d
S
0
t
S
0
t
=
r
t
d
t
,
S
0
0
=
1
.Here r
(
t
) is the stochastic interest rate.The stock market index priceS
(
t
) at time t
≥
0 can be described by the model:(3)
d
S
t
S
t
=
r
t
+
u
t
d
t
+
v
d
W
~
t
,
S
0
=
s
0
,where u
(
t
) is another positive real-valued function and the constant v
>
0 denotes the volatility rate of risk asset. ρ
2
∈
[
-
1,1
] is the correlation coefficient of two Brownian motions W
^
(
t
) and W
~
(
t
) and E
[
W
^
(
t
)
,
W
~
(
t
)
]
=
ρ
2
·
t.In model (3), u
(
t
)
>
0 is a natural hypothesis, which means that the investment on stock can make more money than investing on cash asset, which is consistent with the economic principle of high-risk high yield. In fact, only when u
(
t
)
>
0, more people are willing to invest on equities. For simplicity, we also can take the positive real-valued function u
(
t
) as a positive constant u
>
0, which does not bring substantial impact on the problem investigated in this paper.


## 2.2. The Surplus Process of Insurer

We use the classical Cramér-Lundberg model to describe the surplus process of insurer:(4)
d
R
t
=
c
d
t
-
d
∑
i
=
1
N
t
Y
i
,where N
(
t
) is a homogeneous Poisson process with intensity λ and the claim sizes {
Y
i
,
i
≥
1
} are independent and identically distributed positive random variables with E
[
Y
i
]
=
μ and E
[
Y
i
2
]
=
σ
2. According to the expected value premium principle, the pure premium rate of insurer is c
=
(
1
+
θ
)
λ
μ and θ
>
0 is the safety loading of insurer.The diffusion approximation of the compound Poisson processd
(
∑
i
=
1
N
(
t
)
Y
i
) is given by(5)
d
∑
i
=
1
N
t
Y
i
=
λ
μ
d
t
-
λ
σ
d
W
t
,where W
(
t
) is a standard Brownian motion, which is independent with W
~
(
t
). Moreover, we also assume E
[
W
(
t
)
,
W
^
(
t
)
]
=
ρ
1
·
t, ρ
1
∈
[
-
1,1
] denoting the correlation coefficient of W
(
t
) and W
^
(
t
).In order to avoid unbearable risk, the insurance company can buy proportional reinsurance to transfer part of the potential risk to reinsurance company. The insurer’s retention proportion isp
∈
[
0,1
], and the remaining 1
-
p will be borne by the reinsurer. η
>
0 is the safety loading of reinsurer:(6)
d
R
t
=
c
-
1
+
η
1
-
p
λ
μ
d
t
-
p
λ
μ
d
t
-
λ
σ
d
W
t
=
1
+
θ
λ
μ
-
1
+
η
1
-
p
λ
μ
-
p
λ
μ
d
t
+
p
λ
σ
d
W
t
.


## 2.3. The Insurer’s Wealth Process with Reward Budget

The proportion of wealth invested on stock asset is denoted byπ
(
t
) which is a progressively measurable real function, X
(
t
) with X
(
0
)
=
x
0 represents the wealth of the insurer, and then we have(7)
d
X
t
=
d
R
t
+
π
t
X
t
d
S
t
S
t
+
1
-
π
t
X
t
d
S
0
t
S
0
t
=
θ
-
η
λ
μ
+
η
p
λ
μ
+
π
t
u
t
+
r
t
X
t
d
t
+
p
λ
σ
d
W
t
+
v
π
t
X
t
d
W
~
t
.A notationd
=
(
η
-
θ
)
λ
μ
/
η
p
λ
μ
=
(
η
-
θ
)
/
η
p
=
(
1
-
θ
/
η
)
1
/
p and (
η
-
θ
)
λ
μ
=
d
η
p
λ
μ, where d is an important constant which indicates the benefit game on safety loading of insurer and reinsurer. The new notation d introduced here is interesting, existing for every combination of an insurance contract with a proportional reinsurance transaction. If η
=
θ, then d
=
0, with such proportional reinsurance transaction being more favourable for insurer. The actual reinsurance transaction must balance both benefits of insurer and reinsurer, so that 0
<
d
<
1 holds invariably.The reward budget at a fixed experience share is deducted from premium income at each timet:(8)
d
X
t
=
c
-
k
X
t
-
1
+
η
1
-
p
λ
μ
d
t
-
p
λ
μ
d
t
-
λ
σ
d
W
t
+
π
t
1
-
k
X
t
d
S
t
S
t
+
1
-
π
t
1
-
k
X
t
d
S
0
t
S
0
t
=
1
-
d
η
p
λ
μ
+
1
-
k
X
t
r
t
-
k
X
t
+
π
t
1
-
k
X
t
u
t
d
t
+
p
λ
σ
d
W
t
+
π
t
1
-
k
X
t
v
d
W
~
t
.


## 3. Main Results

In this subsection, we consider the reinsurance-investment problem with stochastic interest rate. Following Merton [10], we use an expected utility maximization criterion and assume that the insurer’s objective is to maximize the expected utility from terminal wealth, where the utility U
(
x
) is a power utility (constant relative risk aversion) function: (9)
U
x
=
x
q
,
q
∈
0,1
.This utility function plays a vital role in actuarial mathematics and insurance practice.Using the definition of Arrow-Pratt measure of absolute risk-aversion (ARA)(10)
ARA
x
=
-
U
′
′
x
U
′
x
, the coefficient of absolute risk aversion of power utility U
(
x
) is(11)
ARA
x
=
-
q
q
-
1
x
q
-
2
q
x
q
-
1
=
1
-
q
xwhich decreases with the increase of wealth.According to the definition of Arrow-Pratt measure of relative risk-aversion (RRA)(12)
RRA
x
=
x
·
ARA
x
=
-
x
U
′
′
x
U
′
x
, the coefficient of relative risk aversion of power utility U
(
x
) is(13)
RRA
x
=
x
·
ARA
x
=
-
x
q
q
-
1
x
q
-
2
q
x
q
-
1
=
1
-
qwhich does not change with the level of wealth.Understanding the economic meaning of the two expressions,(
1
-
q
) and (
1
-
q
)
/
x, is important to understand the results in this paper.A strategyΔ
=
(
p
(
t
)
,
π
(
t
)
) is said to be admissible, if ∀
t
∈
[
0
,
T
], Δ is F
(
t
) progressively measurable, and E
[
∫
0
∞
π
2
(
t
)
X
2
(
t
)
]
<
∞, p
(
t
)
∈
[
0,1
], π
(
t
)
∈
[
0
,
+
∞
). Λ denotes the set of all admissible strategies Δ and (8) has a unique (strong) solution. Suppose the insurer has a utility function U
(
x
) which is strictly concave and continuously differentiable on (
-
∞
,
+
∞
) and aims to maximize the expected utility of his/her terminal wealth.For an admissible strategyΔ
=
(
p
(
t
)
,
π
(
t
)
), the value function J
Δ
(
t
,
x
,
r
) from state (
x
,
r
) at time t is defined as(14)
J
Δ
t
,
x
,
r
=
E
U
X
T
∣
X
t
=
x
,
r
t
=
r
,and the objective function is(15)
H
t
,
x
,
r
=
max
Δ
∈
Λ
⁡
J
Δ
t
,
x
,
r
,with the boundary condition(16)
H
T
,
x
,
r
=
U
x
.The goal of the insurer is to find an optimal strategyΔ
∗
=
(
p
∗
(
t
)
,
π
∗
(
t
)
) such that J
Δ
∗
(
t
,
x
,
r
)
=
H
(
t
,
x
,
r
), where p
∗
(
t
) is called the optimal reinsurance strategy and π
∗
(
t
) is called the optimal investment strategy.


### 3.1. An Insurer with Reward Budget under the Vasicek Model

The stochastic interest rate followed by the Vasicek model is(17)
d
r
t
=
β
-
α
r
t
d
t
+
δ
d
W
^
t
. For mathematical treatment, there is no real difficulties besides representation of results being more complex when u
(
t
) is a positive real function. So we take u
(
t
) as a constant u
>
0 for the representation of results being more concise.Theorem 1.
For the optimal investment-reinsurance problem (15), assuming the utility is given by a power function U
(
x
)
=
x
q,  0
<
q
<
1, T denoting the terminal time, then the optimal strategy Δ
∗
=
(
p
∗
,
π
∗
) is given by(18)
p
∗
=
-
x
1
-
d
η
λ
μ
+
ρ
1
δ
λ
σ
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
λ
σ
2
q
-
1
,
π
∗
=
-
u
+
ρ
2
δ
v
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
v
2
1
-
k
q
-
1
.And the optimal value function is (19)
H
t
,
x
,
r
=
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
r
x
q
,in which(20)
A
=
e
2
T
α
-
1
+
k
2
q
v
2
δ
2
1
-
q
+
q
ρ
1
2
+
q
ρ
2
2
σ
2
2
α
,
B
=
e
2
T
α
T
u
2
α
2
σ
2
-
2
-
1
+
k
q
u
v
α
δ
ρ
2
σ
2
+
v
2
B
1
+
B
2
+
B
3
,
C
=
2
e
2
T
α
-
1
+
k
v
σ
C
1
+
q
C
2
α
,
E
t
=
e
2
t
α
-
1
+
k
2
q
v
2
δ
2
1
-
q
+
q
ρ
1
2
+
q
ρ
2
2
σ
2
2
α
,
F
t
=
e
2
T
α
t
u
2
α
2
σ
2
-
2
-
1
+
k
q
u
v
α
δ
ρ
2
σ
2
+
v
2
F
1
+
F
2
+
F
3
,
G
t
=
2
e
t
α
+
T
α
-
1
+
k
v
σ
-
v
α
β
σ
+
G
1
+
q
-
u
α
δ
ρ
2
σ
+
G
2
α
,where B
1, B
2, B
3, C
1, C
2, F
1, F
2, F
3, G
1, and G
2 are given by(21)
B
1
=
-
1
+
k
2
q
δ
2
1
+
q
-
1
+
ρ
1
2
+
ρ
2
2
σ
2
,
B
2
=
α
2
-
1
+
d
2
η
2
λ
μ
2
+
2
k
-
1
+
q
σ
2
,
B
3
=
2
-
1
+
k
α
σ
-
β
σ
+
q
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
C
1
=
-
v
α
β
σ
+
-
1
+
k
q
2
v
δ
2
-
1
+
ρ
1
2
+
ρ
2
2
σ
,
C
2
=
-
u
α
δ
ρ
2
σ
+
v
-
1
+
k
δ
2
σ
+
α
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
F
1
=
-
1
+
k
2
q
δ
2
1
+
q
-
1
+
ρ
1
2
+
ρ
2
2
σ
2
,
F
2
=
α
2
-
1
+
d
2
η
2
λ
μ
2
+
2
k
-
1
+
q
σ
2
,
F
3
=
2
-
1
+
k
α
σ
-
β
σ
+
q
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
G
1
=
-
1
+
k
q
2
v
δ
2
-
1
+
ρ
1
2
+
ρ
2
2
σ
,
G
2
=
v
-
1
+
k
δ
2
σ
+
α
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
.Proof.
For anyH
(
t
,
x
,
r
)
∈
C
1,2
,
2
(
[
0
,
T
]
×
R
+
×
R
+
), according to the classical results of Fleming and Soner [11], the generator can be defined as a variational operator: (22)
A
H
Δ
t
,
x
,
r
=
H
t
+
H
x
r
1
-
k
x
-
k
x
+
H
r
β
-
α
r
+
H
r
r
1
2
δ
2
+
H
x
1
-
d
p
η
λ
μ
+
H
x
r
ρ
1
δ
p
λ
σ
+
H
x
x
1
2
p
2
λ
σ
2
+
H
x
π
u
t
1
-
k
x
+
H
x
r
ρ
2
δ
v
π
1
-
k
x
+
H
x
x
1
2
v
2
π
2
1
-
k
2
x
2
, where H
t, H
x, H
r, H
x
x, H
r
r, and H
x
r denote the corresponding first- and second-order partial derivatives with respect to the corresponding variables t, x, r, respectively. The Hamilton-Jacobi-Bellman (HJB) equation for problem (15) is given by(23)
sup
Δ
∈
Λ
⁡
A
H
Δ
t
,
x
,
r
=
0with the boundary condition H
(
T
,
x
,
r
)
=
x
q.
According to the first-order necessary condition for the extreme, the following equation is obtained from (23): (24)
H
x
1
-
d
η
λ
μ
+
H
x
r
ρ
1
δ
λ
σ
+
H
x
x
p
λ
σ
2
=
0
,
H
x
u
1
-
k
x
+
H
x
r
ρ
2
δ
v
1
-
k
x
+
H
x
x
v
2
π
1
-
k
2
x
2
=
0
.
The solution of the above equation gives the optimal feedback control:(25)
p
∗
=
-
H
x
1
-
d
η
λ
μ
+
H
x
r
ρ
1
δ
λ
σ
H
x
x
λ
σ
2
,
π
∗
=
-
H
x
u
+
H
x
r
ρ
2
δ
v
H
x
x
v
2
1
-
k
x
.
Plugging (25) into (23),(26)
H
t
+
H
x
r
1
-
k
x
-
k
x
+
H
r
β
-
α
r
+
H
r
r
1
2
δ
2
-
H
x
1
-
d
η
λ
μ
+
H
x
r
ρ
1
δ
σ
2
2
H
x
x
σ
2
-
H
x
u
+
H
x
r
ρ
2
δ
v
2
2
v
2
H
x
x
=
0
.
Conjecture of the solution has the following form:(27)
H
t
,
x
,
r
=
f
t
e
h
t
r
x
q
,
H
T
,
x
,
r
=
x
q
,
0
<
q
<
1
,
f
T
=
1
,
h
T
=
0
.
The derivatives of the conjecture with respect tot, x, r, respectively, are(28)
H
t
=
f
′
t
+
f
t
h
′
t
r
e
h
t
r
x
q
,
H
r
=
f
t
e
h
t
r
x
q
h
t
,
H
r
r
=
f
t
e
h
t
r
x
q
h
t
2
,
H
x
=
f
t
e
h
t
r
q
x
q
-
1
,
H
x
x
=
f
t
e
h
t
r
q
q
-
1
x
q
-
2
,
H
x
r
=
f
t
e
h
t
r
q
x
q
-
1
h
t
.
Plugging these derivatives into (26),(29)
f
′
t
+
f
t
h
′
t
r
+
f
t
q
r
1
-
k
-
k
+
f
t
h
t
β
-
α
r
+
f
t
h
t
2
1
2
δ
2
-
1
-
d
η
λ
μ
+
h
t
ρ
1
δ
σ
2
f
t
q
2
q
-
1
σ
2
-
u
+
h
t
ρ
2
δ
v
2
f
t
q
2
q
-
1
v
2
=
0
.
Splitting (29) into two equations,(30)
h
′
t
-
α
h
t
+
q
1
-
k
=
0
,
h
T
=
0
,
f
′
t
f
t
=
k
q
-
h
t
β
-
1
2
δ
2
h
t
2
+
1
-
d
η
λ
μ
+
h
t
ρ
1
δ
σ
2
2
q
-
1
σ
2
+
u
+
h
t
ρ
2
δ
v
2
2
q
-
1
v
2
q
=
0
,
f
T
=
1
.Solving the above two ordinary differential equations, the solutions h
(
t
) and f
(
t
) are given by(31)
h
t
=
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
α
f
t
=
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
,where (32)
A
=
e
2
T
α
-
1
+
k
2
q
v
2
δ
2
1
-
q
+
q
ρ
1
2
+
q
ρ
2
2
σ
2
2
α
,
B
=
e
2
T
α
T
u
2
α
2
σ
2
-
2
-
1
+
k
q
u
v
α
δ
ρ
2
σ
2
+
v
2
B
1
+
B
2
+
B
3
,
C
=
2
e
2
T
α
-
1
+
k
v
σ
C
1
+
q
C
2
α
,
E
t
=
e
2
t
α
-
1
+
k
2
q
v
2
δ
2
1
-
q
+
q
ρ
1
2
+
q
ρ
2
2
σ
2
2
α
,
F
t
=
e
2
T
α
t
u
2
α
2
σ
2
-
2
-
1
+
k
q
u
v
α
δ
ρ
2
σ
2
+
v
2
F
1
+
F
2
+
F
3
,
G
t
=
2
e
t
α
+
T
α
-
1
+
k
v
σ
-
v
α
β
σ
+
G
1
+
q
-
u
α
δ
ρ
2
σ
+
G
2
αand in which B
1, B
2, B
3, C
1, C
2, F
1, F
2, F
3, G
1, and  G
2 are given by (33)
B
1
=
-
1
+
k
2
q
δ
2
1
+
q
-
1
+
ρ
1
2
+
ρ
2
2
σ
2
,
B
2
=
α
2
-
1
+
d
2
η
2
λ
μ
2
+
2
k
-
1
+
q
σ
2
,
B
3
=
2
-
1
+
k
α
σ
-
β
σ
+
q
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
C
1
=
-
v
α
β
σ
+
-
1
+
k
q
2
v
δ
2
-
1
+
ρ
1
2
+
ρ
2
2
σ
,
C
2
=
-
u
α
δ
ρ
2
σ
+
v
-
1
+
k
δ
2
σ
+
α
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
F
1
=
-
1
+
k
2
q
δ
2
1
+
q
-
1
+
ρ
1
2
+
ρ
2
2
σ
2
,
F
2
=
α
2
-
1
+
d
2
η
2
λ
μ
2
+
2
k
-
1
+
q
σ
2
,
F
3
=
2
-
1
+
k
α
σ
-
β
σ
+
q
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
G
1
=
-
1
+
k
q
2
v
δ
2
-
1
+
ρ
1
2
+
ρ
2
2
σ
,
G
2
=
v
-
1
+
k
δ
2
σ
+
α
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
.Plugging (31) into (27), the optimal value function is obtained as follows: (34)
H
t
,
x
,
r
=
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
r
x
q
.
Then the corresponding derivatives being plugged into (25), the optimal strategies Δ
∗ are given by (35)
p
∗
=
-
x
1
-
d
η
λ
μ
+
ρ
1
δ
λ
σ
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
λ
σ
2
q
-
1
,
π
∗
=
-
u
+
ρ
2
δ
v
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
v
2
1
-
k
q
-
1
.Theorem 2 (verification theorem).
LetH
(
t
,
x
,
r
) be a convex, twice differential solution to (23) with the boundary condition (16) such that |
H
(
0
,
x
0
,
r
0
)
|
<
∞. Then, for all t
∈
[
0
,
T
]
,
x
∈
(
0
,
∞
), (i)
H
(
t
,
x
,
r
)
≤
J
Δ
(
t
,
x
,
r
) for every admissible strategy Δ
∈
Λ.
(ii)
If there exists an admissible controlΔ
∗
=
(
p
∗
,
π
∗
) such that (36)
Δ
∗
∈
argmax
⁡
A
Δ
H
t
,
x
∗
t
,
r
∗
t then the value function H
(
t
,
x
,
r
)
=
J
Δ
∗
(
t
,
x
,
r
) and the policy Δ
∗
=
(
p
∗
,
π
∗
) is the optimal strategy corresponding to X
Δ
∗
(
t
) which is the solution of (8).Proof.
The proof is classical; refer to Fleming and Soner [11], Guan and Liang [12, 13], Gu et al. [14], and so on; we omit it here.


### 3.2. A Special Case of an Insurer with Reward Budget under a Constant Interest Rate

In this subsection, we investigate a special case that the interest rater
(
t
) is a constant r. All assumptions are the same as above except that r
(
t
)
≡
r throughout this part; thus, we have (37)
d
S
0
t
S
0
t
=
r
d
t
,
S
0
0
=
1
,
d
S
t
S
t
=
r
+
u
t
d
t
+
v
d
W
~
t
,
S
0
=
s
0
, where r, v are positive constants and u
(
t
) is a positive real-valued function.The surplus processes with claim and reinsurance are the same as (5) and (6), and W
(
t
) and W
~
(
t
) are independent standard Brownian motions. So we get the wealth process with reward budget, investing in the financial market as follows: (38)
d
X
t
=
1
-
k
X
t
r
-
k
X
t
+
1
-
d
η
p
λ
μ
+
π
t
1
-
k
X
t
u
t
d
t
+
π
t
1
-
k
X
t
v
d
W
~
t
+
p
λ
σ
d
W
t
.We consider that the utility U
(
x
) is a power utility (constant relative risk aversion) function (39)
U
x
=
x
q to maximize the expected utility of terminal wealth, where the constant q
∈
(
0,1
) is the relative risk aversion coefficient of a risk averse insurer. For an admissible strategy Δ
=
(
p
(
t
)
,
π
(
t
)
), the value function J
Δ
(
t
,
x
,
r
) from state x at time t is defined as (40)
J
Δ
t
,
x
=
E
U
X
T
∣
X
t
=
x
, and the optimal value function is given by (41)
H
t
,
x
=
max
π
,
p
⁡
J
Δ
t
,
x with the boundary condition H
(
T
,
x
,
r
)
=
x
q.Applying the classical tools of stochastic optimal control, if the optimal value functionH
(
t
,
x
)
∈
C
1,2
(
[
0
,
T
]
×
R
), then H satisfies the following Hamilton-Jacobi-Bellman (HJB) equation: (42)
H
t
+
H
x
r
-
1
+
r
k
x
+
max
p
⁡
H
x
1
-
d
p
η
λ
μ
+
H
x
x
1
2
p
2
λ
σ
2
+
max
π
⁡
H
x
π
u
t
1
-
k
x
+
H
x
x
1
2
v
2
π
2
1
-
k
2
x
2
=
0
. According to the first-order necessary condition of extreme, we obtain that (43)
p
∗
=
-
1
-
d
η
μ
H
x
σ
2
H
x
x
,
π
∗
=
-
u
t
H
x
v
2
1
-
k
x
H
x
x
.Conjecture of the solution has the following form:(44)
H
t
,
x
=
f
t
x
q
,
H
T
,
x
=
x
q
,
0
<
q
<
1
,
f
T
=
1
.The derivatives of the conjecture with respect tot, x, respectively, are(45)
H
t
=
f
′
t
x
q
,
H
x
=
f
t
q
x
q
-
1
,
H
x
x
=
f
t
q
q
-
1
x
q
-
2
.Plugging these derivatives into the HJB equation, we have(46)
f
′
t
+
f
t
q
r
-
1
+
r
k
+
1
-
d
η
μ
2
λ
q
f
t
2
σ
2
1
-
q
+
u
t
2
q
f
t
2
v
2
1
-
q
=
0
. The solution of the above ordinary differential equation is (47)
f
t
=
exp
⁡
q
∫
t
T
r
-
1
+
r
k
+
1
-
d
2
η
2
λ
μ
2
2
σ
2
1
-
q
+
u
τ
2
2
v
2
1
-
q
d
τ
.We obtain the value function(48)
H
t
,
x
=
x
q
exp
⁡
r
-
1
+
r
k
+
1
-
d
2
η
2
λ
μ
2
2
σ
2
1
-
q
q
T
-
t
+
q
∫
t
T
u
τ
2
2
v
2
1
-
q
d
τ and the optimal strategy α
∗
=
(
p
∗
,
π
∗
) as follows: (49)
p
∗
=
1
-
d
1
-
q
·
η
μ
x
σ
2
,
π
∗
=
1
1
-
k
1
-
q
·
u
t
v
2
.For simplicity, if the positive real-valued functionu
(
t
) taken as a positive constant u
>
0, it does not bring substantial impact on the problem investigated. We obtain the optimal value function (50)
H
t
,
x
=
x
q
exp
⁡
r
-
1
+
r
k
+
1
-
d
2
η
2
λ
μ
2
2
σ
2
1
-
q
+
u
2
2
v
2
1
-
q
q
T
-
t and the optimal strategy α
∗
=
(
p
∗
,
π
∗
) as follows: (51)
p
∗
=
1
-
d
η
μ
x
1
-
q
σ
2
,
π
∗
=
u
v
2
1
-
k
1
-
q
.


### 3.3. Some Economic Interpretations

The optimal strategyΔ
∗
=
(
p
∗
,
π
∗
) under a constant interest rate is (52)
p
∗
,
π
∗
=
x
1
-
q
·
1
-
d
η
μ
σ
2
,
1
1
-
k
1
-
q
·
u
v
2and the optimal strategy Δ
∗
=
(
p
∗
,
π
∗
) under the stochastic interest rates is (53)
p
∗
,
π
∗
=
x
1
-
d
η
λ
μ
+
ρ
1
δ
λ
σ
e
-
T
α
-
e
t
α
+
e
T
α
1
-
k
q
/
α
λ
σ
2
1
-
q
,
u
+
ρ
2
δ
v
e
-
T
α
-
e
t
α
+
e
T
α
1
-
k
q
/
α
v
2
1
-
k
1
-
q
.In this subsection, the results obtained under constant interest rate will serve as a reference. So as to comprehend the implicit economic meaning of abstract formula more clearly, we will make the following interpretations.Remark 3.
Just looking at the optimal strategy under the constant interest rate, it is not difficult to give the following economic interpretations.
First, the optimal reinsurance proportion(54)
p
∗
=
x
1
-
q
·
1
-
d
η
μ
σ
2
=
1
1
-
q
/
x
·
1
-
d
η
μ
σ
2 is determined by two key factors (
1
-
q
)
/
x and (
1
-
d
)
η
μ
/
σ
2.
Equation (11) shows that (
1
-
q
)
/
x is the coefficient of absolute risk aversion, which gives a measurable indicator of an insurer’s absolute risk aversion. So the expression 1
/
(
1
-
q
)
/
x represents the risk tolerance coefficient and the optimal reinsurance proportion p
∗ is in direct proportion to the risk tolerance coefficient 1
/
(
1
-
q
)
/
x. Intuitive understanding is that there would be greater risk scale in the insurer’s businesses if an insurer has greater risk tolerance (the risk tolerance coefficient is bigger), so that the insurance company having larger size of risk needs to be transferred to the reinsurance companies.
As for this expression(
1
-
d
)
η
μ
/
σ
2, it represents the rate of underwriting return. An insurer assumes a certain claim risk σ
2, the higher rate of underwriting return which would be more advantageous to the insurance company.
Second, the optimal investment proportion(55)
π
∗
=
1
1
-
q
·
1
1
-
k
·
u
v
2 is also based on three key factors (
1
-
q
), (
1
-
k
), and u
/
v
2.
Equation (13) shows that (
1
-
q
) is the coefficient of relative risk aversion, which gives a measurable indicator of an insurer’s relative risk aversion. So the expression 1
/
(
1
-
q
) represents the relative risk tolerance coefficient and the optimal investment proportion π
∗ is in direct proportion to the relative risk tolerance coefficient 1
/
(
1
-
q
). This is very natural, because the companies of high risk tolerance always prefer risk assets investment for that they is more valued high yield hidden in high risk.
As for the expressionu
/
v
2
=
(
r
+
u
)
-
r
/
v
2, it is similar to the market price of risk in financial market. (
1
-
k
) represents the total capital share an insurer can use to invest. The amount of risk investment is restricted by the insurer’s risk tolerance, so the proportion of risk investment seems smaller if the total capital share is greater.Remark 4.
Comparing the results of constant interest rate with the results of stochastic interest rate, there is one more itemx
·
ρ
1
δ
/
λ
σ
·
[
e
-
T
α
-
e
t
α
+
e
T
α
(
1
-
k
)
q
]
/
(
1
-
q
)
α of the optimal reinsurance proportion under the stochastic interest rate than the optimal reinsurance proportion under a constant interest rate:(56)
p
∗
stochastic
-
p
∗
constant
=
q
1
-
q
·
ρ
1
δ
λ
σ
·
1
-
e
-
T
-
t
α
α
1
-
k
x
.
Ifρ
1
>
0, then p
∗
s
t
o
c
h
a
s
t
i
c
>
p
∗
c
o
n
s
t
a
n
t; else if ρ
1
<
0, then p
∗
s
t
o
c
h
a
s
t
i
c
<
p
∗
c
o
n
s
t
a
n
t.
In fact,[
1
-
e
-
(
T
-
t
)
α
]
/
α is the changing expectations of the interest rate described by a Vasicek model under the risk neutral probability measure. So [
1
-
e
-
(
T
-
t
)
α
]
/
α
(
1
-
k
)
x reflects the change of insurance company’s wealth with time t. q
/
(
1
-
q
) is the characterization of risk aversion of an insurer. ρ
1
δ
/
λ
σ is a ratio of interest rate volatility rate divided by the volatility rate of claim size. It is not hard to see the crucial effects determined by the correlation of the interest rate fluctuations and the volatility of claim amount.
Similarly, there is one more itemρ
2
δ
/
v
·
[
e
-
T
α
-
e
t
α
+
e
T
α
q
]
/
(
1
-
q
)
α of the optimal investment proportion under the stochastic interest rate than the optimal investment proportion under a constant interest rate:(57)
π
∗
stochastic
-
π
∗
constant
=
q
1
-
q
·
ρ
2
δ
v
·
1
-
e
-
T
-
t
α
α
.
Ifρ
2
>
0, then π
∗
s
t
o
c
h
a
s
t
i
c
>
π
∗
c
o
n
s
t
a
n
t; else if ρ
2
<
0, then π
∗
s
t
o
c
h
a
s
t
i
c
<
π
∗
c
o
n
s
t
a
n
t.
Same as the above,[
1
-
e
-
(
T
-
t
)
α
]
/
α is the changing expectations of the interest rate described by a Vasicek model under the risk neutral probability measure. q
/
(
1
-
q
) is the characterization of risk aversion of an insurer. But δ
/
v is a ratio of interest rate volatility rate divided by the volatility rate of risky assets, so that the crucial effects are determined by the correlation of the interest rate fluctuations and the risk assets volatility.
In conclusion, the primary differences from the optimal strategy under a constant interest rate are implicit in the two expressionsx
·
ρ
1
δ
/
λ
σ
·
[
e
-
T
α
-
e
t
α
+
e
T
α
(
1
-
k
)
q
]
/
(
1
-
q
)
α and ρ
2
δ
/
v
·
[
e
-
T
α
-
e
t
α
+
e
T
α
q
]
/
(
1
-
q
)
α, which reflect the correlation of the stochastic interest rate and the volatility of claim amount and the correlation of the stochastic interest rate and the risk assets volatility, respectively. However, this kind of distinction is a synergy effect based on interest rate volatility, risky asset volatility, claim size volatility, investor’s risk aversion, and so forth.Remark 5.
In the last part of this subsection, let us focus on the actual economic meaning of the optimal value function. The optimal value function under a constant interest rate is(58)
H
t
,
x
=
x
q
exp
⁡
r
-
1
+
r
k
+
1
-
d
2
η
2
λ
μ
2
2
σ
2
1
-
q
+
u
2
2
v
2
1
-
q
q
T
-
t
=
x
q
e
1
-
k
r
q
T
-
t
e
-
k
q
T
-
t
exp
⁡
1
-
d
2
η
2
λ
μ
2
2
σ
2
1
-
q
+
u
2
2
v
2
1
-
q
q
T
-
t and the optimal value function under the stochastic interest rate is (59)
H
t
,
x
,
r
=
x
q
·
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
r
=
x
q
·
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
1
-
e
-
T
-
t
α
/
α
r
q
e
1
-
e
-
T
-
t
α
/
α
r
k
q
.
The optimal value function represents the maximum expected utility of wealth, which increases with the wealth increasing. The key factors of wealth changes mainly depend on the following four aspects: the return of risk-free investment, the share of reward budget, the rate of underwriting return, and the yields rate of risk investment. For the optimal value function under a constant interest rate,e
-
(
t
-
T
)
(
1
-
k
)
r
q represents the impact on insurer’s utility of the discounted value at time t, e
-
k
q
(
T
-
t
) indicates the loss of utility caused by the share k of reward budget, (
1
-
d
)
2
η
2
λ
μ
2
/
2
σ
2
(
1
-
q
)
q
(
T
-
t
) denotes the rate of underwriting return for the risk volatility σ
2 in insurance business, and u
2
/
2
v
2
(
1
-
q
)
q
(
T
-
t
) signifies the yields rate of risk investment for the risk volatility v
2. It is worth noting that the risk aversion factor q is included in each expression, because that the utility itself is closely related with risk preference.
For the optimal value function under the stochastic interest rate, the key factors of wealth changes mainly depend on the following aspects: the stochastic interest rate, the share of reward budget, the rate of underwriting return, and the yields rate of risk investment. We can give the similar explanations by reference to economic meanings under the constant interest rate. For example,e
e
-
T
α
-
e
t
α
+
e
T
α
q
/
α
r, e
e
-
T
α
-
e
t
α
+
e
T
α
k
q
/
α
r reflect the action effects of the stochastic interest rate and the share of reward budget, respectively. e
-
e
-
2
T
α
q
[
A
+
B
-
C
]
/
2
(
-
1
+
q
)
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
[
E
(
t
)
+
F
(
t
)
-
G
(
t
)
]
/
2
(
-
1
+
q
)
v
2
α
2
σ
2 denotes the comprehensive effect of the rate of underwriting return and the yields rate of risk investment. But it is hard to entirely separate the influence of one factor apart from the influence of other factors. Arithmetically, a combination of a variety of random factors appears in one expression in a complex way. It is noteworthy that the risk aversion factor q is an important factor of an insurer’s utility.


## 3.1. An Insurer with Reward Budget under the Vasicek Model

The stochastic interest rate followed by the Vasicek model is(17)
d
r
t
=
β
-
α
r
t
d
t
+
δ
d
W
^
t
. For mathematical treatment, there is no real difficulties besides representation of results being more complex when u
(
t
) is a positive real function. So we take u
(
t
) as a constant u
>
0 for the representation of results being more concise.Theorem 1.
For the optimal investment-reinsurance problem (15), assuming the utility is given by a power function U
(
x
)
=
x
q,  0
<
q
<
1, T denoting the terminal time, then the optimal strategy Δ
∗
=
(
p
∗
,
π
∗
) is given by(18)
p
∗
=
-
x
1
-
d
η
λ
μ
+
ρ
1
δ
λ
σ
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
λ
σ
2
q
-
1
,
π
∗
=
-
u
+
ρ
2
δ
v
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
v
2
1
-
k
q
-
1
.And the optimal value function is (19)
H
t
,
x
,
r
=
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
r
x
q
,in which(20)
A
=
e
2
T
α
-
1
+
k
2
q
v
2
δ
2
1
-
q
+
q
ρ
1
2
+
q
ρ
2
2
σ
2
2
α
,
B
=
e
2
T
α
T
u
2
α
2
σ
2
-
2
-
1
+
k
q
u
v
α
δ
ρ
2
σ
2
+
v
2
B
1
+
B
2
+
B
3
,
C
=
2
e
2
T
α
-
1
+
k
v
σ
C
1
+
q
C
2
α
,
E
t
=
e
2
t
α
-
1
+
k
2
q
v
2
δ
2
1
-
q
+
q
ρ
1
2
+
q
ρ
2
2
σ
2
2
α
,
F
t
=
e
2
T
α
t
u
2
α
2
σ
2
-
2
-
1
+
k
q
u
v
α
δ
ρ
2
σ
2
+
v
2
F
1
+
F
2
+
F
3
,
G
t
=
2
e
t
α
+
T
α
-
1
+
k
v
σ
-
v
α
β
σ
+
G
1
+
q
-
u
α
δ
ρ
2
σ
+
G
2
α
,where B
1, B
2, B
3, C
1, C
2, F
1, F
2, F
3, G
1, and G
2 are given by(21)
B
1
=
-
1
+
k
2
q
δ
2
1
+
q
-
1
+
ρ
1
2
+
ρ
2
2
σ
2
,
B
2
=
α
2
-
1
+
d
2
η
2
λ
μ
2
+
2
k
-
1
+
q
σ
2
,
B
3
=
2
-
1
+
k
α
σ
-
β
σ
+
q
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
C
1
=
-
v
α
β
σ
+
-
1
+
k
q
2
v
δ
2
-
1
+
ρ
1
2
+
ρ
2
2
σ
,
C
2
=
-
u
α
δ
ρ
2
σ
+
v
-
1
+
k
δ
2
σ
+
α
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
F
1
=
-
1
+
k
2
q
δ
2
1
+
q
-
1
+
ρ
1
2
+
ρ
2
2
σ
2
,
F
2
=
α
2
-
1
+
d
2
η
2
λ
μ
2
+
2
k
-
1
+
q
σ
2
,
F
3
=
2
-
1
+
k
α
σ
-
β
σ
+
q
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
G
1
=
-
1
+
k
q
2
v
δ
2
-
1
+
ρ
1
2
+
ρ
2
2
σ
,
G
2
=
v
-
1
+
k
δ
2
σ
+
α
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
.Proof.
For anyH
(
t
,
x
,
r
)
∈
C
1,2
,
2
(
[
0
,
T
]
×
R
+
×
R
+
), according to the classical results of Fleming and Soner [11], the generator can be defined as a variational operator: (22)
A
H
Δ
t
,
x
,
r
=
H
t
+
H
x
r
1
-
k
x
-
k
x
+
H
r
β
-
α
r
+
H
r
r
1
2
δ
2
+
H
x
1
-
d
p
η
λ
μ
+
H
x
r
ρ
1
δ
p
λ
σ
+
H
x
x
1
2
p
2
λ
σ
2
+
H
x
π
u
t
1
-
k
x
+
H
x
r
ρ
2
δ
v
π
1
-
k
x
+
H
x
x
1
2
v
2
π
2
1
-
k
2
x
2
, where H
t, H
x, H
r, H
x
x, H
r
r, and H
x
r denote the corresponding first- and second-order partial derivatives with respect to the corresponding variables t, x, r, respectively. The Hamilton-Jacobi-Bellman (HJB) equation for problem (15) is given by(23)
sup
Δ
∈
Λ
⁡
A
H
Δ
t
,
x
,
r
=
0with the boundary condition H
(
T
,
x
,
r
)
=
x
q.
According to the first-order necessary condition for the extreme, the following equation is obtained from (23): (24)
H
x
1
-
d
η
λ
μ
+
H
x
r
ρ
1
δ
λ
σ
+
H
x
x
p
λ
σ
2
=
0
,
H
x
u
1
-
k
x
+
H
x
r
ρ
2
δ
v
1
-
k
x
+
H
x
x
v
2
π
1
-
k
2
x
2
=
0
.
The solution of the above equation gives the optimal feedback control:(25)
p
∗
=
-
H
x
1
-
d
η
λ
μ
+
H
x
r
ρ
1
δ
λ
σ
H
x
x
λ
σ
2
,
π
∗
=
-
H
x
u
+
H
x
r
ρ
2
δ
v
H
x
x
v
2
1
-
k
x
.
Plugging (25) into (23),(26)
H
t
+
H
x
r
1
-
k
x
-
k
x
+
H
r
β
-
α
r
+
H
r
r
1
2
δ
2
-
H
x
1
-
d
η
λ
μ
+
H
x
r
ρ
1
δ
σ
2
2
H
x
x
σ
2
-
H
x
u
+
H
x
r
ρ
2
δ
v
2
2
v
2
H
x
x
=
0
.
Conjecture of the solution has the following form:(27)
H
t
,
x
,
r
=
f
t
e
h
t
r
x
q
,
H
T
,
x
,
r
=
x
q
,
0
<
q
<
1
,
f
T
=
1
,
h
T
=
0
.
The derivatives of the conjecture with respect tot, x, r, respectively, are(28)
H
t
=
f
′
t
+
f
t
h
′
t
r
e
h
t
r
x
q
,
H
r
=
f
t
e
h
t
r
x
q
h
t
,
H
r
r
=
f
t
e
h
t
r
x
q
h
t
2
,
H
x
=
f
t
e
h
t
r
q
x
q
-
1
,
H
x
x
=
f
t
e
h
t
r
q
q
-
1
x
q
-
2
,
H
x
r
=
f
t
e
h
t
r
q
x
q
-
1
h
t
.
Plugging these derivatives into (26),(29)
f
′
t
+
f
t
h
′
t
r
+
f
t
q
r
1
-
k
-
k
+
f
t
h
t
β
-
α
r
+
f
t
h
t
2
1
2
δ
2
-
1
-
d
η
λ
μ
+
h
t
ρ
1
δ
σ
2
f
t
q
2
q
-
1
σ
2
-
u
+
h
t
ρ
2
δ
v
2
f
t
q
2
q
-
1
v
2
=
0
.
Splitting (29) into two equations,(30)
h
′
t
-
α
h
t
+
q
1
-
k
=
0
,
h
T
=
0
,
f
′
t
f
t
=
k
q
-
h
t
β
-
1
2
δ
2
h
t
2
+
1
-
d
η
λ
μ
+
h
t
ρ
1
δ
σ
2
2
q
-
1
σ
2
+
u
+
h
t
ρ
2
δ
v
2
2
q
-
1
v
2
q
=
0
,
f
T
=
1
.Solving the above two ordinary differential equations, the solutions h
(
t
) and f
(
t
) are given by(31)
h
t
=
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
α
f
t
=
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
,where (32)
A
=
e
2
T
α
-
1
+
k
2
q
v
2
δ
2
1
-
q
+
q
ρ
1
2
+
q
ρ
2
2
σ
2
2
α
,
B
=
e
2
T
α
T
u
2
α
2
σ
2
-
2
-
1
+
k
q
u
v
α
δ
ρ
2
σ
2
+
v
2
B
1
+
B
2
+
B
3
,
C
=
2
e
2
T
α
-
1
+
k
v
σ
C
1
+
q
C
2
α
,
E
t
=
e
2
t
α
-
1
+
k
2
q
v
2
δ
2
1
-
q
+
q
ρ
1
2
+
q
ρ
2
2
σ
2
2
α
,
F
t
=
e
2
T
α
t
u
2
α
2
σ
2
-
2
-
1
+
k
q
u
v
α
δ
ρ
2
σ
2
+
v
2
F
1
+
F
2
+
F
3
,
G
t
=
2
e
t
α
+
T
α
-
1
+
k
v
σ
-
v
α
β
σ
+
G
1
+
q
-
u
α
δ
ρ
2
σ
+
G
2
αand in which B
1, B
2, B
3, C
1, C
2, F
1, F
2, F
3, G
1, and  G
2 are given by (33)
B
1
=
-
1
+
k
2
q
δ
2
1
+
q
-
1
+
ρ
1
2
+
ρ
2
2
σ
2
,
B
2
=
α
2
-
1
+
d
2
η
2
λ
μ
2
+
2
k
-
1
+
q
σ
2
,
B
3
=
2
-
1
+
k
α
σ
-
β
σ
+
q
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
C
1
=
-
v
α
β
σ
+
-
1
+
k
q
2
v
δ
2
-
1
+
ρ
1
2
+
ρ
2
2
σ
,
C
2
=
-
u
α
δ
ρ
2
σ
+
v
-
1
+
k
δ
2
σ
+
α
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
F
1
=
-
1
+
k
2
q
δ
2
1
+
q
-
1
+
ρ
1
2
+
ρ
2
2
σ
2
,
F
2
=
α
2
-
1
+
d
2
η
2
λ
μ
2
+
2
k
-
1
+
q
σ
2
,
F
3
=
2
-
1
+
k
α
σ
-
β
σ
+
q
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
,
G
1
=
-
1
+
k
q
2
v
δ
2
-
1
+
ρ
1
2
+
ρ
2
2
σ
,
G
2
=
v
-
1
+
k
δ
2
σ
+
α
-
1
+
d
δ
η
λ
μ
ρ
1
+
β
σ
.Plugging (31) into (27), the optimal value function is obtained as follows: (34)
H
t
,
x
,
r
=
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
r
x
q
.
Then the corresponding derivatives being plugged into (25), the optimal strategies Δ
∗ are given by (35)
p
∗
=
-
x
1
-
d
η
λ
μ
+
ρ
1
δ
λ
σ
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
λ
σ
2
q
-
1
,
π
∗
=
-
u
+
ρ
2
δ
v
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
v
2
1
-
k
q
-
1
.Theorem 2 (verification theorem).
LetH
(
t
,
x
,
r
) be a convex, twice differential solution to (23) with the boundary condition (16) such that |
H
(
0
,
x
0
,
r
0
)
|
<
∞. Then, for all t
∈
[
0
,
T
]
,
x
∈
(
0
,
∞
), (i)
H
(
t
,
x
,
r
)
≤
J
Δ
(
t
,
x
,
r
) for every admissible strategy Δ
∈
Λ.
(ii)
If there exists an admissible controlΔ
∗
=
(
p
∗
,
π
∗
) such that (36)
Δ
∗
∈
argmax
⁡
A
Δ
H
t
,
x
∗
t
,
r
∗
t then the value function H
(
t
,
x
,
r
)
=
J
Δ
∗
(
t
,
x
,
r
) and the policy Δ
∗
=
(
p
∗
,
π
∗
) is the optimal strategy corresponding to X
Δ
∗
(
t
) which is the solution of (8).Proof.
The proof is classical; refer to Fleming and Soner [11], Guan and Liang [12, 13], Gu et al. [14], and so on; we omit it here.


## 3.2. A Special Case of an Insurer with Reward Budget under a Constant Interest Rate

In this subsection, we investigate a special case that the interest rater
(
t
) is a constant r. All assumptions are the same as above except that r
(
t
)
≡
r throughout this part; thus, we have (37)
d
S
0
t
S
0
t
=
r
d
t
,
S
0
0
=
1
,
d
S
t
S
t
=
r
+
u
t
d
t
+
v
d
W
~
t
,
S
0
=
s
0
, where r, v are positive constants and u
(
t
) is a positive real-valued function.The surplus processes with claim and reinsurance are the same as (5) and (6), and W
(
t
) and W
~
(
t
) are independent standard Brownian motions. So we get the wealth process with reward budget, investing in the financial market as follows: (38)
d
X
t
=
1
-
k
X
t
r
-
k
X
t
+
1
-
d
η
p
λ
μ
+
π
t
1
-
k
X
t
u
t
d
t
+
π
t
1
-
k
X
t
v
d
W
~
t
+
p
λ
σ
d
W
t
.We consider that the utility U
(
x
) is a power utility (constant relative risk aversion) function (39)
U
x
=
x
q to maximize the expected utility of terminal wealth, where the constant q
∈
(
0,1
) is the relative risk aversion coefficient of a risk averse insurer. For an admissible strategy Δ
=
(
p
(
t
)
,
π
(
t
)
), the value function J
Δ
(
t
,
x
,
r
) from state x at time t is defined as (40)
J
Δ
t
,
x
=
E
U
X
T
∣
X
t
=
x
, and the optimal value function is given by (41)
H
t
,
x
=
max
π
,
p
⁡
J
Δ
t
,
x with the boundary condition H
(
T
,
x
,
r
)
=
x
q.Applying the classical tools of stochastic optimal control, if the optimal value functionH
(
t
,
x
)
∈
C
1,2
(
[
0
,
T
]
×
R
), then H satisfies the following Hamilton-Jacobi-Bellman (HJB) equation: (42)
H
t
+
H
x
r
-
1
+
r
k
x
+
max
p
⁡
H
x
1
-
d
p
η
λ
μ
+
H
x
x
1
2
p
2
λ
σ
2
+
max
π
⁡
H
x
π
u
t
1
-
k
x
+
H
x
x
1
2
v
2
π
2
1
-
k
2
x
2
=
0
. According to the first-order necessary condition of extreme, we obtain that (43)
p
∗
=
-
1
-
d
η
μ
H
x
σ
2
H
x
x
,
π
∗
=
-
u
t
H
x
v
2
1
-
k
x
H
x
x
.Conjecture of the solution has the following form:(44)
H
t
,
x
=
f
t
x
q
,
H
T
,
x
=
x
q
,
0
<
q
<
1
,
f
T
=
1
.The derivatives of the conjecture with respect tot, x, respectively, are(45)
H
t
=
f
′
t
x
q
,
H
x
=
f
t
q
x
q
-
1
,
H
x
x
=
f
t
q
q
-
1
x
q
-
2
.Plugging these derivatives into the HJB equation, we have(46)
f
′
t
+
f
t
q
r
-
1
+
r
k
+
1
-
d
η
μ
2
λ
q
f
t
2
σ
2
1
-
q
+
u
t
2
q
f
t
2
v
2
1
-
q
=
0
. The solution of the above ordinary differential equation is (47)
f
t
=
exp
⁡
q
∫
t
T
r
-
1
+
r
k
+
1
-
d
2
η
2
λ
μ
2
2
σ
2
1
-
q
+
u
τ
2
2
v
2
1
-
q
d
τ
.We obtain the value function(48)
H
t
,
x
=
x
q
exp
⁡
r
-
1
+
r
k
+
1
-
d
2
η
2
λ
μ
2
2
σ
2
1
-
q
q
T
-
t
+
q
∫
t
T
u
τ
2
2
v
2
1
-
q
d
τ and the optimal strategy α
∗
=
(
p
∗
,
π
∗
) as follows: (49)
p
∗
=
1
-
d
1
-
q
·
η
μ
x
σ
2
,
π
∗
=
1
1
-
k
1
-
q
·
u
t
v
2
.For simplicity, if the positive real-valued functionu
(
t
) taken as a positive constant u
>
0, it does not bring substantial impact on the problem investigated. We obtain the optimal value function (50)
H
t
,
x
=
x
q
exp
⁡
r
-
1
+
r
k
+
1
-
d
2
η
2
λ
μ
2
2
σ
2
1
-
q
+
u
2
2
v
2
1
-
q
q
T
-
t and the optimal strategy α
∗
=
(
p
∗
,
π
∗
) as follows: (51)
p
∗
=
1
-
d
η
μ
x
1
-
q
σ
2
,
π
∗
=
u
v
2
1
-
k
1
-
q
.


## 3.3. Some Economic Interpretations

The optimal strategyΔ
∗
=
(
p
∗
,
π
∗
) under a constant interest rate is (52)
p
∗
,
π
∗
=
x
1
-
q
·
1
-
d
η
μ
σ
2
,
1
1
-
k
1
-
q
·
u
v
2and the optimal strategy Δ
∗
=
(
p
∗
,
π
∗
) under the stochastic interest rates is (53)
p
∗
,
π
∗
=
x
1
-
d
η
λ
μ
+
ρ
1
δ
λ
σ
e
-
T
α
-
e
t
α
+
e
T
α
1
-
k
q
/
α
λ
σ
2
1
-
q
,
u
+
ρ
2
δ
v
e
-
T
α
-
e
t
α
+
e
T
α
1
-
k
q
/
α
v
2
1
-
k
1
-
q
.In this subsection, the results obtained under constant interest rate will serve as a reference. So as to comprehend the implicit economic meaning of abstract formula more clearly, we will make the following interpretations.Remark 3.
Just looking at the optimal strategy under the constant interest rate, it is not difficult to give the following economic interpretations.
First, the optimal reinsurance proportion(54)
p
∗
=
x
1
-
q
·
1
-
d
η
μ
σ
2
=
1
1
-
q
/
x
·
1
-
d
η
μ
σ
2 is determined by two key factors (
1
-
q
)
/
x and (
1
-
d
)
η
μ
/
σ
2.
Equation (11) shows that (
1
-
q
)
/
x is the coefficient of absolute risk aversion, which gives a measurable indicator of an insurer’s absolute risk aversion. So the expression 1
/
(
1
-
q
)
/
x represents the risk tolerance coefficient and the optimal reinsurance proportion p
∗ is in direct proportion to the risk tolerance coefficient 1
/
(
1
-
q
)
/
x. Intuitive understanding is that there would be greater risk scale in the insurer’s businesses if an insurer has greater risk tolerance (the risk tolerance coefficient is bigger), so that the insurance company having larger size of risk needs to be transferred to the reinsurance companies.
As for this expression(
1
-
d
)
η
μ
/
σ
2, it represents the rate of underwriting return. An insurer assumes a certain claim risk σ
2, the higher rate of underwriting return which would be more advantageous to the insurance company.
Second, the optimal investment proportion(55)
π
∗
=
1
1
-
q
·
1
1
-
k
·
u
v
2 is also based on three key factors (
1
-
q
), (
1
-
k
), and u
/
v
2.
Equation (13) shows that (
1
-
q
) is the coefficient of relative risk aversion, which gives a measurable indicator of an insurer’s relative risk aversion. So the expression 1
/
(
1
-
q
) represents the relative risk tolerance coefficient and the optimal investment proportion π
∗ is in direct proportion to the relative risk tolerance coefficient 1
/
(
1
-
q
). This is very natural, because the companies of high risk tolerance always prefer risk assets investment for that they is more valued high yield hidden in high risk.
As for the expressionu
/
v
2
=
(
r
+
u
)
-
r
/
v
2, it is similar to the market price of risk in financial market. (
1
-
k
) represents the total capital share an insurer can use to invest. The amount of risk investment is restricted by the insurer’s risk tolerance, so the proportion of risk investment seems smaller if the total capital share is greater.Remark 4.
Comparing the results of constant interest rate with the results of stochastic interest rate, there is one more itemx
·
ρ
1
δ
/
λ
σ
·
[
e
-
T
α
-
e
t
α
+
e
T
α
(
1
-
k
)
q
]
/
(
1
-
q
)
α of the optimal reinsurance proportion under the stochastic interest rate than the optimal reinsurance proportion under a constant interest rate:(56)
p
∗
stochastic
-
p
∗
constant
=
q
1
-
q
·
ρ
1
δ
λ
σ
·
1
-
e
-
T
-
t
α
α
1
-
k
x
.
Ifρ
1
>
0, then p
∗
s
t
o
c
h
a
s
t
i
c
>
p
∗
c
o
n
s
t
a
n
t; else if ρ
1
<
0, then p
∗
s
t
o
c
h
a
s
t
i
c
<
p
∗
c
o
n
s
t
a
n
t.
In fact,[
1
-
e
-
(
T
-
t
)
α
]
/
α is the changing expectations of the interest rate described by a Vasicek model under the risk neutral probability measure. So [
1
-
e
-
(
T
-
t
)
α
]
/
α
(
1
-
k
)
x reflects the change of insurance company’s wealth with time t. q
/
(
1
-
q
) is the characterization of risk aversion of an insurer. ρ
1
δ
/
λ
σ is a ratio of interest rate volatility rate divided by the volatility rate of claim size. It is not hard to see the crucial effects determined by the correlation of the interest rate fluctuations and the volatility of claim amount.
Similarly, there is one more itemρ
2
δ
/
v
·
[
e
-
T
α
-
e
t
α
+
e
T
α
q
]
/
(
1
-
q
)
α of the optimal investment proportion under the stochastic interest rate than the optimal investment proportion under a constant interest rate:(57)
π
∗
stochastic
-
π
∗
constant
=
q
1
-
q
·
ρ
2
δ
v
·
1
-
e
-
T
-
t
α
α
.
Ifρ
2
>
0, then π
∗
s
t
o
c
h
a
s
t
i
c
>
π
∗
c
o
n
s
t
a
n
t; else if ρ
2
<
0, then π
∗
s
t
o
c
h
a
s
t
i
c
<
π
∗
c
o
n
s
t
a
n
t.
Same as the above,[
1
-
e
-
(
T
-
t
)
α
]
/
α is the changing expectations of the interest rate described by a Vasicek model under the risk neutral probability measure. q
/
(
1
-
q
) is the characterization of risk aversion of an insurer. But δ
/
v is a ratio of interest rate volatility rate divided by the volatility rate of risky assets, so that the crucial effects are determined by the correlation of the interest rate fluctuations and the risk assets volatility.
In conclusion, the primary differences from the optimal strategy under a constant interest rate are implicit in the two expressionsx
·
ρ
1
δ
/
λ
σ
·
[
e
-
T
α
-
e
t
α
+
e
T
α
(
1
-
k
)
q
]
/
(
1
-
q
)
α and ρ
2
δ
/
v
·
[
e
-
T
α
-
e
t
α
+
e
T
α
q
]
/
(
1
-
q
)
α, which reflect the correlation of the stochastic interest rate and the volatility of claim amount and the correlation of the stochastic interest rate and the risk assets volatility, respectively. However, this kind of distinction is a synergy effect based on interest rate volatility, risky asset volatility, claim size volatility, investor’s risk aversion, and so forth.Remark 5.
In the last part of this subsection, let us focus on the actual economic meaning of the optimal value function. The optimal value function under a constant interest rate is(58)
H
t
,
x
=
x
q
exp
⁡
r
-
1
+
r
k
+
1
-
d
2
η
2
λ
μ
2
2
σ
2
1
-
q
+
u
2
2
v
2
1
-
q
q
T
-
t
=
x
q
e
1
-
k
r
q
T
-
t
e
-
k
q
T
-
t
exp
⁡
1
-
d
2
η
2
λ
μ
2
2
σ
2
1
-
q
+
u
2
2
v
2
1
-
q
q
T
-
t and the optimal value function under the stochastic interest rate is (59)
H
t
,
x
,
r
=
x
q
·
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
r
=
x
q
·
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
1
-
e
-
T
-
t
α
/
α
r
q
e
1
-
e
-
T
-
t
α
/
α
r
k
q
.
The optimal value function represents the maximum expected utility of wealth, which increases with the wealth increasing. The key factors of wealth changes mainly depend on the following four aspects: the return of risk-free investment, the share of reward budget, the rate of underwriting return, and the yields rate of risk investment. For the optimal value function under a constant interest rate,e
-
(
t
-
T
)
(
1
-
k
)
r
q represents the impact on insurer’s utility of the discounted value at time t, e
-
k
q
(
T
-
t
) indicates the loss of utility caused by the share k of reward budget, (
1
-
d
)
2
η
2
λ
μ
2
/
2
σ
2
(
1
-
q
)
q
(
T
-
t
) denotes the rate of underwriting return for the risk volatility σ
2 in insurance business, and u
2
/
2
v
2
(
1
-
q
)
q
(
T
-
t
) signifies the yields rate of risk investment for the risk volatility v
2. It is worth noting that the risk aversion factor q is included in each expression, because that the utility itself is closely related with risk preference.
For the optimal value function under the stochastic interest rate, the key factors of wealth changes mainly depend on the following aspects: the stochastic interest rate, the share of reward budget, the rate of underwriting return, and the yields rate of risk investment. We can give the similar explanations by reference to economic meanings under the constant interest rate. For example,e
e
-
T
α
-
e
t
α
+
e
T
α
q
/
α
r, e
e
-
T
α
-
e
t
α
+
e
T
α
k
q
/
α
r reflect the action effects of the stochastic interest rate and the share of reward budget, respectively. e
-
e
-
2
T
α
q
[
A
+
B
-
C
]
/
2
(
-
1
+
q
)
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
[
E
(
t
)
+
F
(
t
)
-
G
(
t
)
]
/
2
(
-
1
+
q
)
v
2
α
2
σ
2 denotes the comprehensive effect of the rate of underwriting return and the yields rate of risk investment. But it is hard to entirely separate the influence of one factor apart from the influence of other factors. Arithmetically, a combination of a variety of random factors appears in one expression in a complex way. It is noteworthy that the risk aversion factor q is an important factor of an insurer’s utility.


## 4. Numerical Simulations

In this section, we give some qualitative results to show the relationship between the optimal control strategies (optimal investment strategy and optimal proportional reinsurance strategy, resp.) and the parameters in our model. Unless otherwise stated, the hypothetical values of model parameters are as follows:(60)
x
=
3
,
η
=
0.1
,
λ
=
2
,
μ
=
0.2
,
b
=
0.05
,
v
=
1
,
k
=
0.1
,
T
=
10
,
δ
=
0.5
,
σ
=
2
,
α
=
5
,
β
=
0.1
,
ρ
1
=
0.3
,
ρ
2
=
0.2
,
r
=
0.08
,
q
=
0.35
,
u
=
0.15
,
d
=
0.3
.The parameteru is a risk premium which shows the excess return when the risk-free interest rate r is subtracted from the expected return of risk assets. σ
,
δ
,
v are volatility parameters, indicating the risk volatility of claim size, interest rate, or risk assets, respectively. ρ
1 plays an important role in risk interconnection between interest rate and claim size. ρ
2 plays an important role in risk collaboration between interest rate and risk asset. Such risk correlations need more special attention when you make decisions.It is known that the value function is the optimal expected utility of insurer’ wealth when the insurer takes the optimal strategy. First, let us pay attention to the optimal value function; see Figures1, 2, 3, and 4: (61)
H
t
,
x
,
r
=
e
-
e
-
2
T
α
q
A
+
B
-
C
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
e
-
2
T
α
q
E
t
+
F
t
-
G
t
/
2
-
1
+
q
v
2
α
2
σ
2
·
e
-
e
-
T
α
-
e
t
α
+
e
T
α
-
1
+
k
q
/
α
r
x
q
.Figure 1
The optimal value function varies with respect tox.Figure 2
The optimal value function varies with respect tor.Figure 3
The optimal value function varies with respect tou at low level.Figure 4
The optimal value function varies with respect tou at high level.In the case of other parameters being fixed, the influence of interest rate on the value function is not obvious, only when interest rates were significantly higher than normal value which can lead to a significant change in the value function. From the numerical illustration, we also can see that the growth trend of value function with the increase of interest rates will not appear within a short time period. Because the profit of investment on risk-free asset is so meager that it needs such a long enough time to show its observable growth trend.With interest rates at a certain level, the risk premiumu of risk investment is not too much higher than the interest rate, and the impact of risk investment income on the value function is not obvious. Only when the risk investment income is too much higher than the interest rate, the influence of risk investment income on the utility is significant, which can be seen from the impact on the value function. In addition, the long-term observation and short-term observation still have big differences.Next, we show the variation trend of optimal reinsurance proportion:(62)
p
∗
=
x
1
-
d
η
μ
σ
2
1
-
q
+
ρ
1
δ
e
-
T
α
-
e
t
α
+
e
T
α
1
-
k
q
/
α
λ
σ
1
-
q
=
x
η
μ
σ
2
·
1
-
d
1
-
q
+
ρ
1
·
δ
λ
σ
·
e
-
T
α
-
e
t
α
+
e
T
α
1
-
k
α
q
1
-
q with respect to several key parameters x, η, μ. See Figures 5 and 6.Figure 5
The optimal reinsurance proportion varies with respect tox at different η.Figure 6
The optimal reinsurance proportion varies with respect toμ at different η.As we can see from the expression of optimal reinsurance proportionp
∗ immediately, it is mainly influenced by the average claim size λ
μ, the standard deviation of claim volatility λ
σ, and the variance of claim size σ
2. Based on the possible influence of the interest rate fluctuation on the volatility of claim size, the volatility of interest rate δ, and correlation coefficient ρ
1 are also risk factors that appeared in the optimal reinsurance proportion strategy.At last, the variation tendency of optimal investment strategy with respect tou, v
(63)
π
∗
=
u
+
ρ
2
δ
v
e
-
T
α
-
e
t
α
+
e
T
α
1
-
k
q
/
α
v
2
1
-
k
1
-
q
=
u
v
2
+
ρ
2
·
δ
v
e
-
T
α
-
e
t
α
+
e
T
α
1
-
k
q
α
·
1
1
-
k
1
-
q will be illustrated; see Figures 7 and 8.Figure 7
The optimal investment strategy varies with respect tou at different k.Figure 8
The optimal investment strategy varies with respect tov at different u.With the return of risk investment being more higher than the interest rate, the share of investment on risk asset will be higher. Under the condition of having positive yields higher than the interest rate, if the fluctuation of risk investment is more smaller, the share of risk investment will increase more rapidly. If risk asset fluctuates greatly, despite having higher investment income than the risk-free asset, the growth of investment proportion is still very small.It is worth to say that the above analysis and graphics just focus on some key parameters. In fact, the influence of other parameters is also very important. The analysis based on different views can lead to different conclusions, and the work here is just serving as a modest spur to induce someone to come forward with his valuable contributions.


## 5. Conclusion

Based on the viewpoint of interaction between insurance company and policy-holder, we introduce reward budget into the insurance risk management process in this paper, which makes the insurer’s income to be dynamic. Some results are obtained under the Vasicek stochastic interest rate model; the results of which also can be generalized to some other interest rate models.



---

*Source: 1967872-2016-01-26.xml*