# The Study of Mean-Variance Risky Asset Management with State-Dependent Risk Aversion under Regime Switching Market

Authors: Shuang Li; Yu Yang; Yanli Zhou; Yonghong Wu; Xiangyu Ge
Journal: Journal of Function Spaces
 (2021)

---

## Abstract

How do investors require a distribution of the wealth among multiple risky assets while facing the risk of the uncontrollable payment for random liabilities? To cope with this problem, firstly, this paper explores the approach of asset-liability management under the state-dependent risk aversion with only risky assets, which has been considered under a continuous-time Markov regime-switching setting. Next, based on this realistic modelling, an extended Hamilton-Jacob-Bellman (HJB) system has been necessarily established for solving the optimization problem of asset-liability management. It has been derived closed-form analytical expressions applied in the time-inconsistent investment with optimal control theory to see that happens to the optimal value of the function. Ultimately, numerical examples presented with comparisons of the analytical results under different market conditions are exposed to analyse numerically the developed mean variance asset liability management strategy. We find that our proposed model can explain the financial phenomena more effectively and accurately.


---

## Body

## 1. Introduction

Portfolio optimization selection problem, well known as an essential topic in financial markets, has been done in deep researches by many scholars after the first reported by Markowitz [1]. The most frequently used method of optimal asset allocation strategies is HJB equation, i.e., the Hamilton-Jacobi-Bellman equation (see Detemple and Fernando [2], Björk et al. [3]). In the analysis of portfolio optimization, utility function and several system parameters are given to find the optimal values of the control parameters to realise the final utility maximization. Previous researches in this area are classified for the endogenous habit formation [2], the classic constant relative risk aversion (CRRA) by Yu and Yuan [4], the hyperbolic discounting [3], and the utilities like the mean-variance utility proposed by Li et al. [5]. In recent paper by Li et al. [6], the analytical solution portfolio optimization problem involving stochastic short-term interest rates is provided, which can be controlled by the mean-variance utility function with state dependent risk aversion (SDRA). The paper [6] uses the Nash equilibrium for the subgame strategy to concrete analytical expressions of value function and control policy of equilibrium and figure out under the condition of the stochastic short-term interest rates, how do investors with “natural risk aversion” achieve optimal control policies by simplifying financial settings.Under the framework of mean variance equilibrium asset liability management with SDRA, some extended models have been constructed, such as the mean-variance asset-liability management problem by regime-switching models, as well as mean-variance models with only risky assets (see Bening and Koroley [7]; the asset-only models to asset-liability models have been greatly expanded by Yao et al. [8, 9]). A geometric method raised by Leippold et al. [10] is supposed to apply into the multiperiod mean-variance asset-liability management model by taking the implied mean-variance of liability frontier into consideration. A study by Chiu and Li [11] reported that the influence of the rebalancing frequency is quantified to determine the allocation of optimal initial funds. The work of extension into a continuous-time setting has been developed with the aid of a stochastic linear quadratic control approach. Based on the assumptions used in Leippold et al. [10], analytical results have been derived in a complete market with discussing the impact of liability on the optimal funding ratio. To construct more realistic models, more focus has been put on studying the asset-liability management under the market behavior in the face of many restriction conditions, for example, an uncertain investment horizon (see Li and Ng [12]; Li and Yao [13]), regime-switching to describe phenomena between “Bullish” and “Bearish” markets (see Elliott et al. [14]; Wei et al. [15]; Wu and Li [16]; Wu and Chen [17]; Yu [18]), the choice of optimal portfolio selection for assets with transaction costs without short sales (see Li et al. [19]), portfolio selection under partial information (see Xiong and Zhou [20]), bankruptcy control (see Li and Li [21]), jump-diffusion in financial markets (see Lim [22]; Zeng and Li [23]), and stochastic volatility and stochastic interest rates (see Lim [22]; Lim and Zhou [24]). Also, various studies of assets and liabilities management problem have been carried out in some particular field with application in insurance and pension fund, including Drijver [25] for pension funds Hilli et al. [26] for a Finnish pension company, and Gerstner et al. [27] and Chiu and Wong [28] for life insurance policies.Among them, regime-switching models have become popular in finance and related fields, which is expected to describe the characteristics of different markets (called “Bullish” versus“Bearish”). A limited number of regimes have been applied to represent the various patterns of the market states. According to diverse financial markets as the change market pattern occurs, indices for instance the interest rate, appreciation rates, and volatilities of stock and liability may be different. Boyle and Draviam’s study [29] is an interesting example of regime-switching modelling applied in option pricing achieved by [29], followed by Elliott and Siu [30], who have embedded the regime-switching modelling into the bond valuation, the concept of which has been put forward in the portfolio selection problem by Zhou and Yin [31], Chen et al. [32], and Chen and Yang [33]. The research studied in [32, 33] involves both the asset-liability feature and Markovian regime-switching modelling. As we all know, the models with only risk assets are valuable to be studied. Yao et al. [8] were working on the research of the continuous-time mean-variance model for only risky assets. It is rare for risk-free asset in reality; as a matter of fact, a relatively long investment is considered, corresponding to the stochastic nature real interest rates and the inflation risk (see Viceira [34]). The previous method of a nominal risk-free asset incorporated into the market will simplify the process of selecting portfolio but degenerate the GMV strategy to a bank deposit strategy with zero risks, which is not favourable to investors. Besides, the empirical evidence in the study by DeMiguel et al. [35] shows that the static global minimum-variance (GMV) strategy with only risky assets (derived by Markowitz [1]) tends to be better in performance out-of-sample among all estimated optimal strategies. Then in general, the properties of the time-consistent MV strategies have been shown in a market only with risky assets by Chi [36] on the analysis of Yao et al. [8] and Zhang et al. [37].In this paper, on the basis of the work of Björk et al. [3], it is determinate to make a further realistic financial model, and it makes sense to select a regime-switching market with only risky assets. Afterwards, the general expansion of the HJB equation will be reached according to the control theory with time inconsistency by Björk and Murgoci [38]. Finally, it proceeds the numerical illustrations to show our extended results and state the relationships with previous researches. The rest of the paper is completed as follows: the setting of the financial market will be explained in Section 2, with the developed structure of mean-variance asset-liability management with state-dependent risk aversion in a regime-switching market with only risky assets. Also, the HJB equation is generalized to the general situation. In Section 3, three different cases with derived solutions will be illustrated in details. More numerical examples are presented in Section 4 with corresponding figures and illustrations, and a conclusion is given in Section 5.


## 2. Model Formulation

In a given probability space filtered,Ω,P,F,It0≤t≤T, let Wt=W1t,W2t,⋯,Wmt′ be a standard m-dimensional Brownian motion with definition of Ω,P,F over the period of 0,T. Since the involution of individual investments has been found, a few number of investors will not make much effect on the whole market. The mode of the market dynamics is described by a Markov chain process αt. For that sense, the processes of Wt and αt are independent of each other. It=σWs,αs;0≤s≤t could be augmented in the case of all the P-null sets in F, where F=IT. Some finite T is used to denote the range of investment time. All random variables taken into consideration, in this paper, are defined within this filtered probability space. Assuming that there are d regimes for the market state, it means that the Markov chain αt gets one of the values from the set of 1,2,3,⋯,d every time. By assumption, a generator Q=qkjd×d in the Markov chain with the stationary transition probabilities such that pkjt=ℙαs+t=j∣αs=k, where s,t≥0,k,j=1,2,⋯,d, qkj=d/dtpkjtt=0 and ∑j=1dqkj=0,qjj=−∑j=1dqkj<0,qkj>0. A financial market with continuous-time under the standard assumptions has been considered. Concretely speaking, the market assumptions in this paper are listed here with permission for continuous trading, no transaction cost or tax in trading, and infinitely divisible assets.


### 2.1. Financial Market

#### 2.1.1. Assets

Suppose that an investor decides to allocate his wealth amongn+1 risky assets. The prices of these risky assets meet the following requirements of stochastic differential equations (SDE) driven by the geometric Brownion motion (GBM) (1):
(1)dPit=Pitbit,αtdt+∑h=1mσiht,αtdWht,0≤t≤T,Pi0=pi>0,i=0,1,2,⋯,n,where Pit means the initial prices of the risky assets, pi,i=0,1,2,⋯,n; αt is defined as volatility factor; and b0t,αt,b1t,αt,⋯,bnt,αt and σiht,αtn+1×m refer to the appreciation rate vector and the volatility matrix of these assets, respectively, with assumption of positive continuous bounded deterministic functions of time t. As mentioned above, the GBM vector Wt=W1t,W2t,⋯,Wmt′ is supposed to have a detailed description of all the random factors influencing the prices of risky assets.


#### 2.1.2. Wealth Process

It is assumed that an investor endowed with an initial wealthX0 at time 0 is intended to invest in the market dynamically through the period of 0,T. Here, Xt stands for the total wealth at time t for an investor, and uit denotes the amount investment in asset i and Nit for the total of units of asset i in an investor’s portfolio, i=0,1,2⋯,n. The sum of investment in the 0th asset after the deduction of liability is described as Xt−∑i=1nuit. Therefore, under the conditions above, the wealth held by the investor at time t, Xt is shown as follows:
(2)dXt=∑i=0nNitdPit=∑i=0nuitdPitPit.To further simplifydXt in equation (2), the SDE will be represented as
(3)dXt=b0t,αtXt+B′t,αtutdt+Xtσ0′t,αt+u′tσt,αtdWt,X0=X0,α0=k0,where
(4)ut=u1t,u2t,⋯,unt′n×1,Bt,k=b1t,k−b0t,k,b2t,k−b0t,k⋯,bnt,k−b0t,k′n×1,σit,k=σi1t,k,σi2t,k,⋯,σimt,k′m×1,σt,k=σ1t,k−σ0t,k,σ2t,k−σ0t,k,⋯,σnt,k−σ0t,k′n×m.We assume that all the functions are measurable and uniformly bounded in0,T. Here, LF2t,T;ℝn+1 is denoted as the set of all ℝn+1-valued and measurable stochastic processes fs,αs are adjusted to Iss≥t on 0,T such that
(5)E∫tTfs,αs2ds<+∞.


#### 2.1.3. Liability Process

In fact, the investor in the financial market is exposed to the uncontrollable liability, with value process by the following SDE:(6)dLt=μt,αtdt+ρ′t,αtdWt,L0=l0,where Lt is the stochastic liability process and l0 is defined as the initial value of the liability. Besides, μt,αt and ρt,αt=ρ1t,αt,ρ2t,αt,⋯,ρmt,αt′ are expressed as the appreciation and volatility in liability, respectively, on the assumption of stochastic functions at time t with Markov process αt. In addition, generally, the liability is functioned as the real liability excluding the random income of the investor. As a result, it turns out to be negative liabilities; the random income of the investor can be more than the real liability.Remark 1.
It is clearly to see the correlation between liability value and risky assets in the dynamic processes bym-dimensional geometric Brownian motion Wt. Since the investment portfolio has n+1 risky assets and one liability, it leads to m≥n+2 in the asset liability management model.


#### 2.1.4. Surplus Process

LetSt=Xt−Lt represent the current surplus of our fortune. By substituting the wealth and liability processes, we have
(7)dSt=b0t,αtSt−μt,αt+B′t,αtutdt+Stσ0′t,αt−ρ′t,αt+u′tσt,αtdWt,dLt=μt,αtdt+ρ′t,αtdWt,S0=s0=x0−l0,L0=l0,α0=k0.


### 2.2. Mean-Variance Risky Asset Management (MVRAM)

First,U0,T is denoted as the set of all avaliable strategies u0t,u1t,⋯,unt over 0,T. Naturally, the MVRAM problem will give emphasis on finding optimal admissible strategy to maximize mean-variance utility at terminal time T. So the objective function of Jt,s,l,k,u and the equilibrium value function of Vt,s,l,k are described mathematically as follows:
(8)Jt,s,k,u=Et,s,kSuT−γs,k2Vart,s,kSuT,(9)Vt,s,k=maxu·∈U0,TEt,s,kST−γs,k2Vart,s,kST=Jt,s,k,u∗.where Et,s,k·=E·∣Sut,k=s,Lut,k=l,αut=k, in which Sut,k,Lut,k,αut successively represent the surplus process, liability, and market dynamics obtained by using the control strategy u=u0t,u1t,⋯,unt, and rs,k means the risk aversion coefficient depending on s and k. As a result, Jt,s,k,u can be defined as
(10)Jt,s,k,u=Et,s,kFs,k,SuT+Gs,k,Et,s,kSuT,where
(11)Fs,k,y=y−γs,k2y2,Gs,k,y=γs,k2y2,where in here y represents Et,s,kSuT.Second, letAbe infinitesimal generator, for any fixed u∈U; the controlled infinitesimal generator Au corresponded to
(12)AuWt,s,k=∂Wt,s,k∂t+b0s−μ+B′u∂Wt,s,k∂s+∑j=1dqkjWt,s,j+12sσ0′−ρ′+u′σ∂Wt,s,k∂s2sσ0−ρ+σ′u.Based on the analysis of Björk Bjrk and Murgoci [38], with the definition of equilibrium control in equation (9) and the infinitesimal generator Au in equation (12), the HJB equation will be extended as follow, as well as the verification of theorem.Theorem 1 (verification theorem).
It is assumed thatV,f,g is a solution to the following extended HJB system with the supremum of control law u^ in the equation. Then, u^ is subject to an equilibrium control law, and V is supposed to the corresponding value function.
(13)supu∈UAuVt,s,k−Auft,s,k,s,k+Aufs,kt,s,k−AuG∘gt,s,k+Hugt,s,k=0,0≤t≤T,Au∧fm,ht,s,k=0,0≤t≤T,Au∧gt,s,k=0,0≤t≤T,VT,s,k=s,fm,hT,s,k=s−γm,h2s2,gT,s,k=s,with the following conditions and definitions:
(14)fm,ht,s,k=ft,s,k,m,h,G∘gt,s,k=Gs,k,gt,s,k,Hugt,s,k=Gys,k,y×Augt,s,k,Gys,k,y=∂Gs,k,y∂y.
Moreover,f and g have the following probabilistic representations:
(15)ft,s,k,m,h=Et,s,kFm,h,Su∧T,gt,s,k=Et,s,kSu∧T,Vt,s,k=ft,s,k+Gs,k,gt,s,k.Proof.
On the basis of the HJB equation in [2] and the objective function [3], it can be derived as
(16)Vt,St,kt=supu∈UJt,St,kt,Uand consequently have
(17)Jt,St,kt,U=Et,St,ktFs,k,STU+Gs,k,Et,St,ktSTU,where the forms of the function F and function G are described in (11).
Whenl>t, we have
(18)Jl,Sl,kl,U=El,Sl,klFsl,kl,STU+Gsl,kl,El,Sl,klSTU,El,Sl,klFsl,kl,STU=fUl,Sl,kl,sl,kl,El,Sl,klSTU=gUl,Sl,kl.
Hence, equation (18) above can simply be represented as
(19)Jl,Sl,kl,U=fUl,Sl,kl,sl,kl+GUsl,kl,gUl,Sl,kl.
So the expectations of the equation can be shown(20)Et,St,ktJl,Sl,kl,U=Et,St,ktfUl,Sl,kl,sl,kl+Et,St,ktGsl,kl,gUl,Sl,kl,and substituting this result into the definition of (17), we then have
(21)Et,St,ktJl,Sl,kl,U=Jt,St,kt,U+Et,St,ktfUl,Sl,kl,sl,kl−Et,St,ktFs,k,STU+Et,St,ktGsl,kl,gUl,Sl,kl−Gs,k,Et,St,ktSTU.
After the process of iteration, we obtain(22)Et,St,ktFsl,kl,STU=Et,St,ktEl,Sl,klFsl,kl,STU=Et,St,ktfUl,Sl,kl,sl,kl,(23)Et,St,ktSTU=Et,St,ktEl,Sl,klSTU=Et,St,ktgUl,Sl,kl.
By substituting the results of (22) and (23) back into the equation of (21), we can get
(24)Et,St,ktJl,Sl,kl,U−Et,St,ktfUl,Sl,kl,sl,kl+Et,St,ktfUl,Sl,kl,sl,kl−Jt,St,kt,U−Et,St,ktGsl,kl,gUl,Sl,kl+Gs,k,Et,St,ktgUl,Sl,kl=0.
Then(25)supu∈UEt,St,ktJl,Sl,kl,U−Et,St,ktfUl,Sl,kl,sl,kl+Et,St,ktfUl,Sl,kl,sl,kl−Jt,St,kt,U−Et,St,ktGsl,kl,gUl,Sl,kl+Gst,kt,Et,St,ktgUl,Sl,kl=0.
Through our proposed problem (16) with the definition of the control law in the classic work, it can be found out that the control U coincides with the equilibrium law u^ in l,T, and we formulate the following results:
(26)Jl,Sl,kl,u^=Vl,Sl,kl,fUl,Sl,kl,s,k=fl,Sl,kl,s,k,gUl,Sl,kl=gl,Sl,kl.
Thus, the optimization problem of (25) can be solved as
(27)supu∈UEt,St,ktVl,Sl,kl−Et,St,ktfl,Sl,kl,sl,kl+Et,St,ktfl,Sl,kl,sl,kl−Vt,St,kt−Et,St,ktGsl,kl,gl,Sl,kl+Gyt,Et,St,ktgUl,Sl,kl=0.
Here, by using the operator denotations of similarity, we have(28)Et,St,ktVl,Sl,kl−Vt,St,kt=AuV,Et,St,ktfl,Sl,kl,sl,kl=Auf,Et,St,ktfl,Sl,kl=Aufs,k,Et,St,ktGsl,kl,gl,Sl,kl=AuG,Gst,kt,Et,St,ktgUl,Sl,kl=Hug.
The derivation of the extended HJB equation with stochastic volatility will be given as,(29)supu∈UAuVt,s,k−Auft,s,k,s,k+Aufs,kt,s,k−AuG∗gt,s,k+Hugt,s,k=0,0≤t≤T.


## 2.1. Financial Market

### 2.1.1. Assets

Suppose that an investor decides to allocate his wealth amongn+1 risky assets. The prices of these risky assets meet the following requirements of stochastic differential equations (SDE) driven by the geometric Brownion motion (GBM) (1):
(1)dPit=Pitbit,αtdt+∑h=1mσiht,αtdWht,0≤t≤T,Pi0=pi>0,i=0,1,2,⋯,n,where Pit means the initial prices of the risky assets, pi,i=0,1,2,⋯,n; αt is defined as volatility factor; and b0t,αt,b1t,αt,⋯,bnt,αt and σiht,αtn+1×m refer to the appreciation rate vector and the volatility matrix of these assets, respectively, with assumption of positive continuous bounded deterministic functions of time t. As mentioned above, the GBM vector Wt=W1t,W2t,⋯,Wmt′ is supposed to have a detailed description of all the random factors influencing the prices of risky assets.


### 2.1.2. Wealth Process

It is assumed that an investor endowed with an initial wealthX0 at time 0 is intended to invest in the market dynamically through the period of 0,T. Here, Xt stands for the total wealth at time t for an investor, and uit denotes the amount investment in asset i and Nit for the total of units of asset i in an investor’s portfolio, i=0,1,2⋯,n. The sum of investment in the 0th asset after the deduction of liability is described as Xt−∑i=1nuit. Therefore, under the conditions above, the wealth held by the investor at time t, Xt is shown as follows:
(2)dXt=∑i=0nNitdPit=∑i=0nuitdPitPit.To further simplifydXt in equation (2), the SDE will be represented as
(3)dXt=b0t,αtXt+B′t,αtutdt+Xtσ0′t,αt+u′tσt,αtdWt,X0=X0,α0=k0,where
(4)ut=u1t,u2t,⋯,unt′n×1,Bt,k=b1t,k−b0t,k,b2t,k−b0t,k⋯,bnt,k−b0t,k′n×1,σit,k=σi1t,k,σi2t,k,⋯,σimt,k′m×1,σt,k=σ1t,k−σ0t,k,σ2t,k−σ0t,k,⋯,σnt,k−σ0t,k′n×m.We assume that all the functions are measurable and uniformly bounded in0,T. Here, LF2t,T;ℝn+1 is denoted as the set of all ℝn+1-valued and measurable stochastic processes fs,αs are adjusted to Iss≥t on 0,T such that
(5)E∫tTfs,αs2ds<+∞.


### 2.1.3. Liability Process

In fact, the investor in the financial market is exposed to the uncontrollable liability, with value process by the following SDE:(6)dLt=μt,αtdt+ρ′t,αtdWt,L0=l0,where Lt is the stochastic liability process and l0 is defined as the initial value of the liability. Besides, μt,αt and ρt,αt=ρ1t,αt,ρ2t,αt,⋯,ρmt,αt′ are expressed as the appreciation and volatility in liability, respectively, on the assumption of stochastic functions at time t with Markov process αt. In addition, generally, the liability is functioned as the real liability excluding the random income of the investor. As a result, it turns out to be negative liabilities; the random income of the investor can be more than the real liability.Remark 1.
It is clearly to see the correlation between liability value and risky assets in the dynamic processes bym-dimensional geometric Brownian motion Wt. Since the investment portfolio has n+1 risky assets and one liability, it leads to m≥n+2 in the asset liability management model.


### 2.1.4. Surplus Process

LetSt=Xt−Lt represent the current surplus of our fortune. By substituting the wealth and liability processes, we have
(7)dSt=b0t,αtSt−μt,αt+B′t,αtutdt+Stσ0′t,αt−ρ′t,αt+u′tσt,αtdWt,dLt=μt,αtdt+ρ′t,αtdWt,S0=s0=x0−l0,L0=l0,α0=k0.


## 2.1.1. Assets

Suppose that an investor decides to allocate his wealth amongn+1 risky assets. The prices of these risky assets meet the following requirements of stochastic differential equations (SDE) driven by the geometric Brownion motion (GBM) (1):
(1)dPit=Pitbit,αtdt+∑h=1mσiht,αtdWht,0≤t≤T,Pi0=pi>0,i=0,1,2,⋯,n,where Pit means the initial prices of the risky assets, pi,i=0,1,2,⋯,n; αt is defined as volatility factor; and b0t,αt,b1t,αt,⋯,bnt,αt and σiht,αtn+1×m refer to the appreciation rate vector and the volatility matrix of these assets, respectively, with assumption of positive continuous bounded deterministic functions of time t. As mentioned above, the GBM vector Wt=W1t,W2t,⋯,Wmt′ is supposed to have a detailed description of all the random factors influencing the prices of risky assets.


## 2.1.2. Wealth Process

It is assumed that an investor endowed with an initial wealthX0 at time 0 is intended to invest in the market dynamically through the period of 0,T. Here, Xt stands for the total wealth at time t for an investor, and uit denotes the amount investment in asset i and Nit for the total of units of asset i in an investor’s portfolio, i=0,1,2⋯,n. The sum of investment in the 0th asset after the deduction of liability is described as Xt−∑i=1nuit. Therefore, under the conditions above, the wealth held by the investor at time t, Xt is shown as follows:
(2)dXt=∑i=0nNitdPit=∑i=0nuitdPitPit.To further simplifydXt in equation (2), the SDE will be represented as
(3)dXt=b0t,αtXt+B′t,αtutdt+Xtσ0′t,αt+u′tσt,αtdWt,X0=X0,α0=k0,where
(4)ut=u1t,u2t,⋯,unt′n×1,Bt,k=b1t,k−b0t,k,b2t,k−b0t,k⋯,bnt,k−b0t,k′n×1,σit,k=σi1t,k,σi2t,k,⋯,σimt,k′m×1,σt,k=σ1t,k−σ0t,k,σ2t,k−σ0t,k,⋯,σnt,k−σ0t,k′n×m.We assume that all the functions are measurable and uniformly bounded in0,T. Here, LF2t,T;ℝn+1 is denoted as the set of all ℝn+1-valued and measurable stochastic processes fs,αs are adjusted to Iss≥t on 0,T such that
(5)E∫tTfs,αs2ds<+∞.


## 2.1.3. Liability Process

In fact, the investor in the financial market is exposed to the uncontrollable liability, with value process by the following SDE:(6)dLt=μt,αtdt+ρ′t,αtdWt,L0=l0,where Lt is the stochastic liability process and l0 is defined as the initial value of the liability. Besides, μt,αt and ρt,αt=ρ1t,αt,ρ2t,αt,⋯,ρmt,αt′ are expressed as the appreciation and volatility in liability, respectively, on the assumption of stochastic functions at time t with Markov process αt. In addition, generally, the liability is functioned as the real liability excluding the random income of the investor. As a result, it turns out to be negative liabilities; the random income of the investor can be more than the real liability.Remark 1.
It is clearly to see the correlation between liability value and risky assets in the dynamic processes bym-dimensional geometric Brownian motion Wt. Since the investment portfolio has n+1 risky assets and one liability, it leads to m≥n+2 in the asset liability management model.


## 2.1.4. Surplus Process

LetSt=Xt−Lt represent the current surplus of our fortune. By substituting the wealth and liability processes, we have
(7)dSt=b0t,αtSt−μt,αt+B′t,αtutdt+Stσ0′t,αt−ρ′t,αt+u′tσt,αtdWt,dLt=μt,αtdt+ρ′t,αtdWt,S0=s0=x0−l0,L0=l0,α0=k0.


## 2.2. Mean-Variance Risky Asset Management (MVRAM)

First,U0,T is denoted as the set of all avaliable strategies u0t,u1t,⋯,unt over 0,T. Naturally, the MVRAM problem will give emphasis on finding optimal admissible strategy to maximize mean-variance utility at terminal time T. So the objective function of Jt,s,l,k,u and the equilibrium value function of Vt,s,l,k are described mathematically as follows:
(8)Jt,s,k,u=Et,s,kSuT−γs,k2Vart,s,kSuT,(9)Vt,s,k=maxu·∈U0,TEt,s,kST−γs,k2Vart,s,kST=Jt,s,k,u∗.where Et,s,k·=E·∣Sut,k=s,Lut,k=l,αut=k, in which Sut,k,Lut,k,αut successively represent the surplus process, liability, and market dynamics obtained by using the control strategy u=u0t,u1t,⋯,unt, and rs,k means the risk aversion coefficient depending on s and k. As a result, Jt,s,k,u can be defined as
(10)Jt,s,k,u=Et,s,kFs,k,SuT+Gs,k,Et,s,kSuT,where
(11)Fs,k,y=y−γs,k2y2,Gs,k,y=γs,k2y2,where in here y represents Et,s,kSuT.Second, letAbe infinitesimal generator, for any fixed u∈U; the controlled infinitesimal generator Au corresponded to
(12)AuWt,s,k=∂Wt,s,k∂t+b0s−μ+B′u∂Wt,s,k∂s+∑j=1dqkjWt,s,j+12sσ0′−ρ′+u′σ∂Wt,s,k∂s2sσ0−ρ+σ′u.Based on the analysis of Björk Bjrk and Murgoci [38], with the definition of equilibrium control in equation (9) and the infinitesimal generator Au in equation (12), the HJB equation will be extended as follow, as well as the verification of theorem.Theorem 1 (verification theorem).
It is assumed thatV,f,g is a solution to the following extended HJB system with the supremum of control law u^ in the equation. Then, u^ is subject to an equilibrium control law, and V is supposed to the corresponding value function.
(13)supu∈UAuVt,s,k−Auft,s,k,s,k+Aufs,kt,s,k−AuG∘gt,s,k+Hugt,s,k=0,0≤t≤T,Au∧fm,ht,s,k=0,0≤t≤T,Au∧gt,s,k=0,0≤t≤T,VT,s,k=s,fm,hT,s,k=s−γm,h2s2,gT,s,k=s,with the following conditions and definitions:
(14)fm,ht,s,k=ft,s,k,m,h,G∘gt,s,k=Gs,k,gt,s,k,Hugt,s,k=Gys,k,y×Augt,s,k,Gys,k,y=∂Gs,k,y∂y.
Moreover,f and g have the following probabilistic representations:
(15)ft,s,k,m,h=Et,s,kFm,h,Su∧T,gt,s,k=Et,s,kSu∧T,Vt,s,k=ft,s,k+Gs,k,gt,s,k.Proof.
On the basis of the HJB equation in [2] and the objective function [3], it can be derived as
(16)Vt,St,kt=supu∈UJt,St,kt,Uand consequently have
(17)Jt,St,kt,U=Et,St,ktFs,k,STU+Gs,k,Et,St,ktSTU,where the forms of the function F and function G are described in (11).
Whenl>t, we have
(18)Jl,Sl,kl,U=El,Sl,klFsl,kl,STU+Gsl,kl,El,Sl,klSTU,El,Sl,klFsl,kl,STU=fUl,Sl,kl,sl,kl,El,Sl,klSTU=gUl,Sl,kl.
Hence, equation (18) above can simply be represented as
(19)Jl,Sl,kl,U=fUl,Sl,kl,sl,kl+GUsl,kl,gUl,Sl,kl.
So the expectations of the equation can be shown(20)Et,St,ktJl,Sl,kl,U=Et,St,ktfUl,Sl,kl,sl,kl+Et,St,ktGsl,kl,gUl,Sl,kl,and substituting this result into the definition of (17), we then have
(21)Et,St,ktJl,Sl,kl,U=Jt,St,kt,U+Et,St,ktfUl,Sl,kl,sl,kl−Et,St,ktFs,k,STU+Et,St,ktGsl,kl,gUl,Sl,kl−Gs,k,Et,St,ktSTU.
After the process of iteration, we obtain(22)Et,St,ktFsl,kl,STU=Et,St,ktEl,Sl,klFsl,kl,STU=Et,St,ktfUl,Sl,kl,sl,kl,(23)Et,St,ktSTU=Et,St,ktEl,Sl,klSTU=Et,St,ktgUl,Sl,kl.
By substituting the results of (22) and (23) back into the equation of (21), we can get
(24)Et,St,ktJl,Sl,kl,U−Et,St,ktfUl,Sl,kl,sl,kl+Et,St,ktfUl,Sl,kl,sl,kl−Jt,St,kt,U−Et,St,ktGsl,kl,gUl,Sl,kl+Gs,k,Et,St,ktgUl,Sl,kl=0.
Then(25)supu∈UEt,St,ktJl,Sl,kl,U−Et,St,ktfUl,Sl,kl,sl,kl+Et,St,ktfUl,Sl,kl,sl,kl−Jt,St,kt,U−Et,St,ktGsl,kl,gUl,Sl,kl+Gst,kt,Et,St,ktgUl,Sl,kl=0.
Through our proposed problem (16) with the definition of the control law in the classic work, it can be found out that the control U coincides with the equilibrium law u^ in l,T, and we formulate the following results:
(26)Jl,Sl,kl,u^=Vl,Sl,kl,fUl,Sl,kl,s,k=fl,Sl,kl,s,k,gUl,Sl,kl=gl,Sl,kl.
Thus, the optimization problem of (25) can be solved as
(27)supu∈UEt,St,ktVl,Sl,kl−Et,St,ktfl,Sl,kl,sl,kl+Et,St,ktfl,Sl,kl,sl,kl−Vt,St,kt−Et,St,ktGsl,kl,gl,Sl,kl+Gyt,Et,St,ktgUl,Sl,kl=0.
Here, by using the operator denotations of similarity, we have(28)Et,St,ktVl,Sl,kl−Vt,St,kt=AuV,Et,St,ktfl,Sl,kl,sl,kl=Auf,Et,St,ktfl,Sl,kl=Aufs,k,Et,St,ktGsl,kl,gl,Sl,kl=AuG,Gst,kt,Et,St,ktgUl,Sl,kl=Hug.
The derivation of the extended HJB equation with stochastic volatility will be given as,(29)supu∈UAuVt,s,k−Auft,s,k,s,k+Aufs,kt,s,k−AuG∗gt,s,k+Hugt,s,k=0,0≤t≤T.


## 3. Solution Scheme

### 3.1. The Case with a Generatedγs,k

Theorem 2.
Under the general form of states dependent risk aversion, the optimal control strategy of MVRAM among risky assets is(30)u^t,s,k=−σσ´−1fss,k+γs,kggsfsss,k+γs,kggssB+sσσ0−σρ,where f.s,k and g represent the function f.s,kt,s,k and gt,s,k by partial equations in (32).Proof.
By using the definition of the infinitesimal generator, we simplify and thus have(31)AuVt,s,k−Auft,s,k,s,k+Aufs,kt,s,k−AuG∘gt,s,k+Hugt,s,k=Aufs,kt,s,k+Hugt,s,k.
The resulting extended HJB equations can be rewritten as(32)supu∈Ufts,kt,s,k+γs,kgt,s,kgtt,s,k+μlfls,kt,s,k+γs,kgt,s,kglt,s,k+b0s−μ+B′ufss,kt,s,k+γs,kgt,s,kgst,s,k+12sσ0′−ρ′+u′σfsss,kt,s,k+γs,kgt,s,kgsst,s,ksσ0′+σ0′−ρ′l+u′σ′+∑j=1dqkjfs,kt,s,j+γs,kgt,s,kgt,s,j=0,ftm,ht,s,k+b0s−μ+B′u^fsm,ht,s,k+12sσ0′+−ρ′+u∧′σfssm,ht,s,ksσ0′−ρ′+u∧′σ′+∑j=1dqkjfm,ht,s,j=0,gtt,s,k+b0s−μ+B′u^gst,s,k+12sσ0′−ρ′+u∧′σgsst,s,ksσ0′−ρ′+u∧′σ′+∑j=1dqkjgt,s,j=0,fT,s,k,m,h=s−γm,h2s2,gT,s,k=s.
Adding up all the terms related tou in (32), we have
(33)12u′σσ′ufsss,k+γs,kggss+fss,k+γs,kggsB′u+sσ0′−ρ′′σufsss,k+γs,kggss=12ε1u′σσ′u+2ε1ε2B′σσ´−1σ+sσ0′−ρ′ε1σ′u=12ε1σ′u+1ε1ε2σ′σσ′−1B+sσ0′−ρ′ε1´σ′u+1ε1ε2B′σσ′−1σ+sσ0′−ρ′ε1′−12ε1ε2σ′σσ′−1B+sσ0′−ρ′ε1´´ε2σ′σσ′−1B+sσ0′−ρ′ε1,where
(34)ε1=fsss,k+γs,kggss,ε2=fss,k+γs,kggs.
Therefore, the first-order condition forut corresponding to optimal strategy can be described as
(35)u^t,s,k=−σσ′−1ε2ε1B+sσσ0−σρ=−σσ′−1fss,k+γs,kggsfsss,k+γs,kggssB+sσσ0−σρ.


### 3.2. The Case with a Natural Choiceγs,k

Here, we haveγs,k=γk/s as a special form of γs,k. Then, equation (30) has been changed into the following form:
(36)u^t,s,k=−σσ′−1sfss,k+γkggssfsss,k+γkggssB+sσσ0−σρ.Then, we make use of the following natural choice ofγs,k:
(37)gt,s,k=at,ks+nt,k,ft,s,k,m,h=at,ks+nt,k−γh2mAt,ks2+2Dt,ks+Nt,k.By differentiation, we have(38)gt=ȧs+ṅ,gs=a,gss=0,ft=ȧs+ṅ−γh2mȦs2+2Ḋs+Ṅ,fs=a−γhmAs+D,fss=−γhmA.Substituting the above expressions intou^t,s,k, we have
(39)u^=k1s+k2,where
(40)k1=σσ´−1a−γkA+γka2γkAB−σσ0,k2=σσ´−1a·n−DAB+σρ.Then, we can simplify(41)b0s−μ+B′u^=b0+B′k1s+B′k2−μ=M1s+M2,sσ0−ρ+σ′u^=σ0+σ′k1s+σ′k2−ρ=H1s+H2.By substituting these expressions into (32), we also have the following alternative expressions for Au∧fm,h and Au∧g:
(42)Au∧fm,h=−γh2mȦ+2AM1+H1′H1A+∑j=1dqkjAjs2+ȧ−γh2m2Ḋ+M1a−γhmD+M2−γhmA−γhmAH1′H2+∑j=1dqkjaj−γh2m2Djs+ṅ−γh2mṄ+M2a−γhmD−γh2mAH2′H2+∑j=1dqkjnj−γh2mNj=0,Au∧g=ȧ+M1a+∑j=1dqkjajs+ṅ+M2a+∑j=1dqkjnj=0.From equation (42), we can have the following system of ordinary differential equations (ODE):
(43)Ȧ+2M1A+H1′H1A+∑j=1dqkjAj=0,ȧ−γhmḊ+M1a−γhmD−γhmM2A−γhmAH1′H2+∑j=1dqkjaj−γhmDj=0,ṅ−γh2mṄ+M2a−γhmD−γh2mAH2′H2+∑j=1dqkjnj−γh2mNj=0,ȧ+M1a+∑j=1dqkjaj=0,ṅ+M2a+∑j=1dqkjnj=0.By simplifying and substituting the expressions ofM1,M2,H1,H2, and k1,k2 into the above equation (43), we have the following system of ODEs:
(44)Ȧ+2b0−δ2−δ1A+a+γa22δ1γ2A+∑j=1dqkjAj=0,Ḋ+b0−δ2−δ1D+a21+γanδ1γA+δ3−μA+∑j=1dqkjDj=0,Ṅ−δ1Aa2n2−D2+2δ3−μD+∑j=1dqkjNj=0,ȧ+b0+a−γA+γa2γAδ1−δ2a+∑j=1dqkjaj=0,ṅ+a2δ1An−aDAδ1+δ3−μa+∑j=1dqkjnj=0.with the terminal conditions
(45)AT,k=aT,k=1,DT,k=NT,k=nT,k=0,where
(46)δ1=B′σσ′−1B,δ2=B′σσ′−1σσ0,δ3=B′σσ′−1σρ.Following the process of simplification, the equilibrium control can be represented as (39)
(47)u^t,s,k=k1t,ks+k2,where
(48)k1=σσ′−1a−γkA+γka2γkAB−σσ0,k2=σσ′−1a·n−DAB+σρ,and the equilibrium value function of correspondence is given as
(49)Vt,s,k=γk2a2t,k−At,k+at,ks+γk2sn2t,k−Nt,k+γkat,knt,k−Dt,k+nt,k,where A·,D·,N·,a·, and n· satisfy the ODE system in (44).


### 3.3. The Case without Liability

By lettinglt=0, the asset liability problem will be tackled as a portfolio selection problem. We have St=Xt, Vt,s,k=Jt,s,k,u^·, and then, the portfolio optimal control u^ turns to be
(50)u^t,x,k=σσ′−1a−γA+γa2γAB−σσ0x+BAa·n−D,which corresponds to equilibrium value function defined as
(51)Vt,x,k=γ2a2−A+ax+γa·n−D+n+γ2xn2−N,where A,a,D,n,N satisfy the following system of ODEs:
(52)Ȧ+2b0−δ2−δ1A+a+γa22δ1γ2A+∑j=1dqkjAj=0,Ḋ+b0−δ2−δ1D+a21+γanδ1γA+∑j=1dqkjDj=0,Ṅ+a2n2−D2Aδ1+∑j=1dqkjNj=0,ȧ+b0+a−γA+γa2γAδ1−δ2a+∑j=1dqkjaj=0,ṅ+a2δ1An−aDAδ1+∑j=1dqkjnj=0,with the terminal conditions
(53)AT,k=aT,k=1,DT,k=NT,k=nT,k=0.The ODEs are still complicated, which requires to make a further restriction such that the first asset is risk-free, namely,b0t,k=rt,σ0t,k=0,⋯,0′, and we have δ2=0, nor does liability. Hence, the portfolio optimal control u^ is shown as follows:
(54)u^t,x,k=σσ′−1a−γA+γa2γABx+BAa·n−D.Then, it comes to the equilibrium value function as(55)Vt,x,k=γ2a2−A+ax+γa·n−D+n+γ2xn2−N,where A,a,D,n,N satisfy the following system of ODEs:
(56)Ȧ+2r−δ1A+a+γa22δ1γ2A+∑j=1dqkjAj=0,Ḋ+r−δ1D+a21+γanδ1γA+∑j=1dqkjDj=0,Ṅ+a2n2−D2Aδ1+∑j=1dqkjNj=0,ȧ+r+a−γA+γa2γAδ1a+∑j=1dqkjaj=0,ṅ+a2δ1An−aDAδ1+∑j=1dqkjnj=0,with the terminal conditions
(57)AT,k=aT,k=1,DT,k=NT,k=nT,k=0.


## 3.1. The Case with a Generatedγs,k

Theorem 2.
Under the general form of states dependent risk aversion, the optimal control strategy of MVRAM among risky assets is(30)u^t,s,k=−σσ´−1fss,k+γs,kggsfsss,k+γs,kggssB+sσσ0−σρ,where f.s,k and g represent the function f.s,kt,s,k and gt,s,k by partial equations in (32).Proof.
By using the definition of the infinitesimal generator, we simplify and thus have(31)AuVt,s,k−Auft,s,k,s,k+Aufs,kt,s,k−AuG∘gt,s,k+Hugt,s,k=Aufs,kt,s,k+Hugt,s,k.
The resulting extended HJB equations can be rewritten as(32)supu∈Ufts,kt,s,k+γs,kgt,s,kgtt,s,k+μlfls,kt,s,k+γs,kgt,s,kglt,s,k+b0s−μ+B′ufss,kt,s,k+γs,kgt,s,kgst,s,k+12sσ0′−ρ′+u′σfsss,kt,s,k+γs,kgt,s,kgsst,s,ksσ0′+σ0′−ρ′l+u′σ′+∑j=1dqkjfs,kt,s,j+γs,kgt,s,kgt,s,j=0,ftm,ht,s,k+b0s−μ+B′u^fsm,ht,s,k+12sσ0′+−ρ′+u∧′σfssm,ht,s,ksσ0′−ρ′+u∧′σ′+∑j=1dqkjfm,ht,s,j=0,gtt,s,k+b0s−μ+B′u^gst,s,k+12sσ0′−ρ′+u∧′σgsst,s,ksσ0′−ρ′+u∧′σ′+∑j=1dqkjgt,s,j=0,fT,s,k,m,h=s−γm,h2s2,gT,s,k=s.
Adding up all the terms related tou in (32), we have
(33)12u′σσ′ufsss,k+γs,kggss+fss,k+γs,kggsB′u+sσ0′−ρ′′σufsss,k+γs,kggss=12ε1u′σσ′u+2ε1ε2B′σσ´−1σ+sσ0′−ρ′ε1σ′u=12ε1σ′u+1ε1ε2σ′σσ′−1B+sσ0′−ρ′ε1´σ′u+1ε1ε2B′σσ′−1σ+sσ0′−ρ′ε1′−12ε1ε2σ′σσ′−1B+sσ0′−ρ′ε1´´ε2σ′σσ′−1B+sσ0′−ρ′ε1,where
(34)ε1=fsss,k+γs,kggss,ε2=fss,k+γs,kggs.
Therefore, the first-order condition forut corresponding to optimal strategy can be described as
(35)u^t,s,k=−σσ′−1ε2ε1B+sσσ0−σρ=−σσ′−1fss,k+γs,kggsfsss,k+γs,kggssB+sσσ0−σρ.


## 3.2. The Case with a Natural Choiceγs,k

Here, we haveγs,k=γk/s as a special form of γs,k. Then, equation (30) has been changed into the following form:
(36)u^t,s,k=−σσ′−1sfss,k+γkggssfsss,k+γkggssB+sσσ0−σρ.Then, we make use of the following natural choice ofγs,k:
(37)gt,s,k=at,ks+nt,k,ft,s,k,m,h=at,ks+nt,k−γh2mAt,ks2+2Dt,ks+Nt,k.By differentiation, we have(38)gt=ȧs+ṅ,gs=a,gss=0,ft=ȧs+ṅ−γh2mȦs2+2Ḋs+Ṅ,fs=a−γhmAs+D,fss=−γhmA.Substituting the above expressions intou^t,s,k, we have
(39)u^=k1s+k2,where
(40)k1=σσ´−1a−γkA+γka2γkAB−σσ0,k2=σσ´−1a·n−DAB+σρ.Then, we can simplify(41)b0s−μ+B′u^=b0+B′k1s+B′k2−μ=M1s+M2,sσ0−ρ+σ′u^=σ0+σ′k1s+σ′k2−ρ=H1s+H2.By substituting these expressions into (32), we also have the following alternative expressions for Au∧fm,h and Au∧g:
(42)Au∧fm,h=−γh2mȦ+2AM1+H1′H1A+∑j=1dqkjAjs2+ȧ−γh2m2Ḋ+M1a−γhmD+M2−γhmA−γhmAH1′H2+∑j=1dqkjaj−γh2m2Djs+ṅ−γh2mṄ+M2a−γhmD−γh2mAH2′H2+∑j=1dqkjnj−γh2mNj=0,Au∧g=ȧ+M1a+∑j=1dqkjajs+ṅ+M2a+∑j=1dqkjnj=0.From equation (42), we can have the following system of ordinary differential equations (ODE):
(43)Ȧ+2M1A+H1′H1A+∑j=1dqkjAj=0,ȧ−γhmḊ+M1a−γhmD−γhmM2A−γhmAH1′H2+∑j=1dqkjaj−γhmDj=0,ṅ−γh2mṄ+M2a−γhmD−γh2mAH2′H2+∑j=1dqkjnj−γh2mNj=0,ȧ+M1a+∑j=1dqkjaj=0,ṅ+M2a+∑j=1dqkjnj=0.By simplifying and substituting the expressions ofM1,M2,H1,H2, and k1,k2 into the above equation (43), we have the following system of ODEs:
(44)Ȧ+2b0−δ2−δ1A+a+γa22δ1γ2A+∑j=1dqkjAj=0,Ḋ+b0−δ2−δ1D+a21+γanδ1γA+δ3−μA+∑j=1dqkjDj=0,Ṅ−δ1Aa2n2−D2+2δ3−μD+∑j=1dqkjNj=0,ȧ+b0+a−γA+γa2γAδ1−δ2a+∑j=1dqkjaj=0,ṅ+a2δ1An−aDAδ1+δ3−μa+∑j=1dqkjnj=0.with the terminal conditions
(45)AT,k=aT,k=1,DT,k=NT,k=nT,k=0,where
(46)δ1=B′σσ′−1B,δ2=B′σσ′−1σσ0,δ3=B′σσ′−1σρ.Following the process of simplification, the equilibrium control can be represented as (39)
(47)u^t,s,k=k1t,ks+k2,where
(48)k1=σσ′−1a−γkA+γka2γkAB−σσ0,k2=σσ′−1a·n−DAB+σρ,and the equilibrium value function of correspondence is given as
(49)Vt,s,k=γk2a2t,k−At,k+at,ks+γk2sn2t,k−Nt,k+γkat,knt,k−Dt,k+nt,k,where A·,D·,N·,a·, and n· satisfy the ODE system in (44).


## 3.3. The Case without Liability

By lettinglt=0, the asset liability problem will be tackled as a portfolio selection problem. We have St=Xt, Vt,s,k=Jt,s,k,u^·, and then, the portfolio optimal control u^ turns to be
(50)u^t,x,k=σσ′−1a−γA+γa2γAB−σσ0x+BAa·n−D,which corresponds to equilibrium value function defined as
(51)Vt,x,k=γ2a2−A+ax+γa·n−D+n+γ2xn2−N,where A,a,D,n,N satisfy the following system of ODEs:
(52)Ȧ+2b0−δ2−δ1A+a+γa22δ1γ2A+∑j=1dqkjAj=0,Ḋ+b0−δ2−δ1D+a21+γanδ1γA+∑j=1dqkjDj=0,Ṅ+a2n2−D2Aδ1+∑j=1dqkjNj=0,ȧ+b0+a−γA+γa2γAδ1−δ2a+∑j=1dqkjaj=0,ṅ+a2δ1An−aDAδ1+∑j=1dqkjnj=0,with the terminal conditions
(53)AT,k=aT,k=1,DT,k=NT,k=nT,k=0.The ODEs are still complicated, which requires to make a further restriction such that the first asset is risk-free, namely,b0t,k=rt,σ0t,k=0,⋯,0′, and we have δ2=0, nor does liability. Hence, the portfolio optimal control u^ is shown as follows:
(54)u^t,x,k=σσ′−1a−γA+γa2γABx+BAa·n−D.Then, it comes to the equilibrium value function as(55)Vt,x,k=γ2a2−A+ax+γa·n−D+n+γ2xn2−N,where A,a,D,n,N satisfy the following system of ODEs:
(56)Ȧ+2r−δ1A+a+γa22δ1γ2A+∑j=1dqkjAj=0,Ḋ+r−δ1D+a21+γanδ1γA+∑j=1dqkjDj=0,Ṅ+a2n2−D2Aδ1+∑j=1dqkjNj=0,ȧ+r+a−γA+γa2γAδ1a+∑j=1dqkjaj=0,ṅ+a2δ1An−aDAδ1+∑j=1dqkjnj=0,with the terminal conditions
(57)AT,k=aT,k=1,DT,k=NT,k=nT,k=0.


## 4. Numerical Example

According to the research studied by Chen et al. [32] and Wei et al. [15], the market state is either “Bullish” or “Bearish” with the assumption of d=2, where Regime 1 refers to “Bullish” and Regime 2 to “Bearish.” In this section, it is exposed to illustrate the results numerically obtained in Section 3. We present the optimal control strategy and the equilibrium value in three different situations: the first situation is state-dependent risk aversion in a regime switching market with one bond and one risky asset without liability; the second situation is different from the first with liability, while the last one is in the same situation but with two risky assets.


### 4.1. The Case with Liability

Similarly, all the related constant parameters in situation two can been seen in Table1, which are specified for the illustrative purpose.Table 1
All the related constant parameters in the case with liability.
Parameter (symbol)Regime 1 (i=1)Regime 2 (i=2)Exit time (Ti)1010Risk-free interest rate (ri)0.040.04Risky asset appreciation rate (bi)0.20.05Risky asset volatility (σi)0.30.07Liability appreciation rate (μi)0.080.04Liability volatility (ρi)0.30.1Differentiation of transition probabilities (qi)0.30.7Risk aversion (γi)0.50.9The corresponding ODE system becomes(58)Ȧ1+2r1−δ11A1+a1+γ1a122δ11γ12A1+q11A1+q12A2=0,Ȧ2+2r2−δ12A2+a2+γ2a222δ12γ22A2+q21A1+q22A2=0,Ḋ1+r1−δ11D1+δ31−μ1A1+a121+γ1a1n1δ11γ1A1+q11D1+q12D2=0,Ḋ2+r2−δ12D2+δ32−μ2A2+a221+γ2a2n2δ12γ2A2+q21D1+q22D2=0,Ṅ1+2δ31−μ1D1+a12n12−D12A1δ11+q11N1+q12N2=0,Ṅ2+2δ32−μ2D2+a22n22−D22A2δ12+q21N1+q22N2=0,ȧ1+r1+a1−γ1A1+γ1a12γ1A1δ11a1+q11a1+q12a2=0,ȧ2+r2+a2−γ2A2+γ2a22γ2A2δ12a2+q21a1+q22a2=0,ṅ1+a12δ11A1n1−a1D1A1δ11+δ31−μ1a1+q11n1+q12n2=0,ṅ2+a22δ12A2n2−a2D2A2δ12+δ32−μ2a2+q21n1+q22n2=0,with the terminal conditions
(59)A1T=A2T=a1T=a2T=1,D1T=D2T=N1T=N2T=n1T=n2T=0.where
(60)δ11=b1−r12σ12,δ12=b2−r22σ22,δ31=b1−r1ρ1σ1,δ32=b2−r2ρ2σ2.


### 4.2. The Case with a Natural Choiceγs,k

Again, all the related constant parameters in situation three have been represented in Table2. All parameters are specified for the illustrative purpose.Table 2
All the related constant parameters in the case with a natural choiceγs,k.
Parameter (symbol)Regime 1 (i=1)Regime 2 (i=2)Exit time (Ti)1010The first risky asset appreciation rate (b0i)0.30.06Other risky asset appreciation rate (b1i)0.20.05The first risky asset volatility (σ0i)0.30.07Other risky asset volatility (σ1i)0.20.09Liability appreciation rate (μi)0.080.04Liability volatility (ρi)0.30.1Differentiation of transition probabilities (qi)0.30.7Risk aversion (γi)0.50.9The resulting ODE system are shown as follows:(61)Ȧ1+2b01−δ21−δ11A1+a1+γ1a122δ11γ12A1+q11A1+q12A2=0,Ȧ2+2b02−δ22−δ12A2+a2+γ2a222δ12γ22A2+q21A1+q22A2=0,Ḋ1+b01−δ21−δ11D1+δ31−μ1A1+a121+γ1a1n1δ11γ1A1+q11D1+q12D2=0,Ḋ2+b02−δ22−δ12D2+δ32−μ2A2+a221+γ2a2n2δ12γ2A2+q21D1+q22D2=0,Ṅ1+2δ31−μ1D1+a12n12−D12A1δ11+q11N1+q12N2=0,Ṅ2+2δ32−μ2D2+a22n22−D22A2δ12+q21N1+q22N2=0,ȧ1+b01−δ21+a1−γ1A1+γ1a12γ1A1δ11a1+q11a1+q12a2=0,ȧ2+b02−δ22+a2−γ2A2+γ2a22γ2A2δ12a2+q21a1+q22a2=0,ṅ1+a12δ11A1n1−a1D1A1δ11+δ31−μ1a1+q11n1+q12n2=0,ṅ2+a22δ12A2n2−a2D2A2δ12+δ32−μ2a2+q21n1+q22n2=0,with the terminal conditions
(62)A1T=A2T=a1T=a2T=1,D1T=D2T=N1T=N2T=n1T=n2T=0.where
(63)δ11=b11−b012σ11−σ012,δ12=b12−b022σ12−σ022,δ21=b11−b01σ01σ11−σ01,δ22=b12−b02σ02σ12−σ02,δ31=b11−b01ρ1σ11−σ01,δ32=b12−b02ρ2σ12−σ02.


### 4.3. The Case without Liability

All the related constant parameters in situation one have been displayed in Table3, which are specified for the illustrative purpose.Table 3
All the related constant parameters in the case without liability.
Parameter (symbol)Regime 1 (i=1)Regime 2 (i=2)Exit time (Ti)1010Risk-free interest rate (ri)0.040.04Risky asset appreciation rate (bi)0.20.05Risky asset volatility (σi)0.30.07Liability appreciation rate (μi)0.080.04Liability volatility (ρi)0.30.1Differentiation of transition probabilities (qi)0.30.7Risk aversion (γi)0.50.9The ODE system of (56) can then be simplified in the form of the following expressions:
(64)Ȧ1+2r1−δ11A1+a1+γ1a122δ11γ12A1+q11A1+q12A2=0,Ȧ2+2r2−δ12A2+a2+γ2a222δ12γ22A2+q21A1+q22A2=0,Ḋ1+r1−δ11D1+a121+γ1a1n1δ11γ1A1+q11D1+q12D2=0,Ḋ2+r2−δ12D2+a221+γ2a2n2δ12γ2A2+q21D1+q22D2=0,Ṅ1+a12n12−D12A1δ11+q11N1+q12N2=0,Ṅ2+a22n22−D22A2δ12+q21N1+q22N2=0,ȧ1+r1+a1−γ1A1+γ1a12γ1A1δ11a1+q11a1+q12a2=0,ȧ2+r2+a2−γ2A2+γ2a22γ2A2δ12a2+q21a1+q22a2=0,ṅ1+a12δ11A1n1−a1D1A1δ11+q11n1+q12n2=0,ṅ2+a22δ12A2n2−a2D2A2δ12+q21n1+q22n2=0,with the terminal conditions
(65)A1T=A2T=a1T=a2T=1,D1T=D2T=N1T=N2T=n1T=n2T=0.where
(66)δ11=b1−r12σ12,δ12=b2−r22σ22.From Figures1–6, in the “Bearish” market, we could see the optimal control strategy u is increasing as time goes by along with the increase of equilibrium value. They capture the optimal strategy of the utility and keep it for the whole time, because the investor pays attention to the utility functions during the entire time. It is reasonable for the time consistent investor to sacrifice parts of the utility or happiness to secure sufficient budget in order to avoid unpredictable deficits. On the other hand, from Figures 1–6 for the “Bullish” market, the optimal control strategy is growing with time while the equilibrium value decreases. At the beginning of the investment in the “Bullish” market, the investor is confident enough to employ the strategy which would lead to a optimize equilibrium value. Over time, the investor will have less investment time to invest to maximize their current utility as well as the benefit from the “Bullish” market because of the state-dependent risk aversion.Figure 1
Optimal Controlu.Figure 2
Equilibrium valueV.Figure 3
Optimal controlu.Figure 4
Equilibrium valueV.Figure 5
Optimal controlu.Figure 6
Equilibrium valueV.In Figures7 and 8, they compare the results of optimal control strategy and the equilibrium value among three different situations under the “Bullish” market.Figure 7
Comparison of optimal controlu in Bull market.Figure 8
Comparison of equilibrium value in Bull market.In Figure7, with the comparison of optimal control strategy (u) in Bull market among three different situations, we find that there is no significant difference between the situations with and without liability. When the investors come across the “Bullish” market, they are willing to invest the risky asset instead of bearing the liability. However, in the situation with all risky assets, the optimal control strategy provides a way to invest in one of risky assets, which highly depends on the parameters of the two risky assets. In Figure 8, the comparison analysis on the equilibrium values in Bull market among three different situations has been conducted. The figure shows decreasing evidence along with levelling off, of which the reason has been explained above. For the situation with all risky assets, the equilibrium value has a more stable trend which also depends on the selected two risky assets.In the last picture, the comparison has been made about the optimal control strategy by varying the risk aversion coefficient. The effect of the risk aversion has been analysed in Figure9, which presents the optimal control strategy (u) under three different risk aversions. Obviously, as the risk aversion increases, it can be seen that the investor is less likely to invest in risky assets, which is consistent with common sense.Figure 9
Optimal controlu under different risk aversions.


## 4.1. The Case with Liability

Similarly, all the related constant parameters in situation two can been seen in Table1, which are specified for the illustrative purpose.Table 1
All the related constant parameters in the case with liability.
Parameter (symbol)Regime 1 (i=1)Regime 2 (i=2)Exit time (Ti)1010Risk-free interest rate (ri)0.040.04Risky asset appreciation rate (bi)0.20.05Risky asset volatility (σi)0.30.07Liability appreciation rate (μi)0.080.04Liability volatility (ρi)0.30.1Differentiation of transition probabilities (qi)0.30.7Risk aversion (γi)0.50.9The corresponding ODE system becomes(58)Ȧ1+2r1−δ11A1+a1+γ1a122δ11γ12A1+q11A1+q12A2=0,Ȧ2+2r2−δ12A2+a2+γ2a222δ12γ22A2+q21A1+q22A2=0,Ḋ1+r1−δ11D1+δ31−μ1A1+a121+γ1a1n1δ11γ1A1+q11D1+q12D2=0,Ḋ2+r2−δ12D2+δ32−μ2A2+a221+γ2a2n2δ12γ2A2+q21D1+q22D2=0,Ṅ1+2δ31−μ1D1+a12n12−D12A1δ11+q11N1+q12N2=0,Ṅ2+2δ32−μ2D2+a22n22−D22A2δ12+q21N1+q22N2=0,ȧ1+r1+a1−γ1A1+γ1a12γ1A1δ11a1+q11a1+q12a2=0,ȧ2+r2+a2−γ2A2+γ2a22γ2A2δ12a2+q21a1+q22a2=0,ṅ1+a12δ11A1n1−a1D1A1δ11+δ31−μ1a1+q11n1+q12n2=0,ṅ2+a22δ12A2n2−a2D2A2δ12+δ32−μ2a2+q21n1+q22n2=0,with the terminal conditions
(59)A1T=A2T=a1T=a2T=1,D1T=D2T=N1T=N2T=n1T=n2T=0.where
(60)δ11=b1−r12σ12,δ12=b2−r22σ22,δ31=b1−r1ρ1σ1,δ32=b2−r2ρ2σ2.


## 4.2. The Case with a Natural Choiceγs,k

Again, all the related constant parameters in situation three have been represented in Table2. All parameters are specified for the illustrative purpose.Table 2
All the related constant parameters in the case with a natural choiceγs,k.
Parameter (symbol)Regime 1 (i=1)Regime 2 (i=2)Exit time (Ti)1010The first risky asset appreciation rate (b0i)0.30.06Other risky asset appreciation rate (b1i)0.20.05The first risky asset volatility (σ0i)0.30.07Other risky asset volatility (σ1i)0.20.09Liability appreciation rate (μi)0.080.04Liability volatility (ρi)0.30.1Differentiation of transition probabilities (qi)0.30.7Risk aversion (γi)0.50.9The resulting ODE system are shown as follows:(61)Ȧ1+2b01−δ21−δ11A1+a1+γ1a122δ11γ12A1+q11A1+q12A2=0,Ȧ2+2b02−δ22−δ12A2+a2+γ2a222δ12γ22A2+q21A1+q22A2=0,Ḋ1+b01−δ21−δ11D1+δ31−μ1A1+a121+γ1a1n1δ11γ1A1+q11D1+q12D2=0,Ḋ2+b02−δ22−δ12D2+δ32−μ2A2+a221+γ2a2n2δ12γ2A2+q21D1+q22D2=0,Ṅ1+2δ31−μ1D1+a12n12−D12A1δ11+q11N1+q12N2=0,Ṅ2+2δ32−μ2D2+a22n22−D22A2δ12+q21N1+q22N2=0,ȧ1+b01−δ21+a1−γ1A1+γ1a12γ1A1δ11a1+q11a1+q12a2=0,ȧ2+b02−δ22+a2−γ2A2+γ2a22γ2A2δ12a2+q21a1+q22a2=0,ṅ1+a12δ11A1n1−a1D1A1δ11+δ31−μ1a1+q11n1+q12n2=0,ṅ2+a22δ12A2n2−a2D2A2δ12+δ32−μ2a2+q21n1+q22n2=0,with the terminal conditions
(62)A1T=A2T=a1T=a2T=1,D1T=D2T=N1T=N2T=n1T=n2T=0.where
(63)δ11=b11−b012σ11−σ012,δ12=b12−b022σ12−σ022,δ21=b11−b01σ01σ11−σ01,δ22=b12−b02σ02σ12−σ02,δ31=b11−b01ρ1σ11−σ01,δ32=b12−b02ρ2σ12−σ02.


## 4.3. The Case without Liability

All the related constant parameters in situation one have been displayed in Table3, which are specified for the illustrative purpose.Table 3
All the related constant parameters in the case without liability.
Parameter (symbol)Regime 1 (i=1)Regime 2 (i=2)Exit time (Ti)1010Risk-free interest rate (ri)0.040.04Risky asset appreciation rate (bi)0.20.05Risky asset volatility (σi)0.30.07Liability appreciation rate (μi)0.080.04Liability volatility (ρi)0.30.1Differentiation of transition probabilities (qi)0.30.7Risk aversion (γi)0.50.9The ODE system of (56) can then be simplified in the form of the following expressions:
(64)Ȧ1+2r1−δ11A1+a1+γ1a122δ11γ12A1+q11A1+q12A2=0,Ȧ2+2r2−δ12A2+a2+γ2a222δ12γ22A2+q21A1+q22A2=0,Ḋ1+r1−δ11D1+a121+γ1a1n1δ11γ1A1+q11D1+q12D2=0,Ḋ2+r2−δ12D2+a221+γ2a2n2δ12γ2A2+q21D1+q22D2=0,Ṅ1+a12n12−D12A1δ11+q11N1+q12N2=0,Ṅ2+a22n22−D22A2δ12+q21N1+q22N2=0,ȧ1+r1+a1−γ1A1+γ1a12γ1A1δ11a1+q11a1+q12a2=0,ȧ2+r2+a2−γ2A2+γ2a22γ2A2δ12a2+q21a1+q22a2=0,ṅ1+a12δ11A1n1−a1D1A1δ11+q11n1+q12n2=0,ṅ2+a22δ12A2n2−a2D2A2δ12+q21n1+q22n2=0,with the terminal conditions
(65)A1T=A2T=a1T=a2T=1,D1T=D2T=N1T=N2T=n1T=n2T=0.where
(66)δ11=b1−r12σ12,δ12=b2−r22σ22.From Figures1–6, in the “Bearish” market, we could see the optimal control strategy u is increasing as time goes by along with the increase of equilibrium value. They capture the optimal strategy of the utility and keep it for the whole time, because the investor pays attention to the utility functions during the entire time. It is reasonable for the time consistent investor to sacrifice parts of the utility or happiness to secure sufficient budget in order to avoid unpredictable deficits. On the other hand, from Figures 1–6 for the “Bullish” market, the optimal control strategy is growing with time while the equilibrium value decreases. At the beginning of the investment in the “Bullish” market, the investor is confident enough to employ the strategy which would lead to a optimize equilibrium value. Over time, the investor will have less investment time to invest to maximize their current utility as well as the benefit from the “Bullish” market because of the state-dependent risk aversion.Figure 1
Optimal Controlu.Figure 2
Equilibrium valueV.Figure 3
Optimal controlu.Figure 4
Equilibrium valueV.Figure 5
Optimal controlu.Figure 6
Equilibrium valueV.In Figures7 and 8, they compare the results of optimal control strategy and the equilibrium value among three different situations under the “Bullish” market.Figure 7
Comparison of optimal controlu in Bull market.Figure 8
Comparison of equilibrium value in Bull market.In Figure7, with the comparison of optimal control strategy (u) in Bull market among three different situations, we find that there is no significant difference between the situations with and without liability. When the investors come across the “Bullish” market, they are willing to invest the risky asset instead of bearing the liability. However, in the situation with all risky assets, the optimal control strategy provides a way to invest in one of risky assets, which highly depends on the parameters of the two risky assets. In Figure 8, the comparison analysis on the equilibrium values in Bull market among three different situations has been conducted. The figure shows decreasing evidence along with levelling off, of which the reason has been explained above. For the situation with all risky assets, the equilibrium value has a more stable trend which also depends on the selected two risky assets.In the last picture, the comparison has been made about the optimal control strategy by varying the risk aversion coefficient. The effect of the risk aversion has been analysed in Figure9, which presents the optimal control strategy (u) under three different risk aversions. Obviously, as the risk aversion increases, it can be seen that the investor is less likely to invest in risky assets, which is consistent with common sense.Figure 9
Optimal controlu under different risk aversions.


## 5. Conclusions

In this paper, the mean variance model of asset-liability management has been discussed in the case of state-dependent risk aversion with only risky assets. Based on the continuous-time Markov regime-switching, this paper derives an analytical optimal control expressions theoretically in a more realistic financial market and then makes numerical analysis on a series of special cases. From the numerical results, this paper reveals the feasibility and application of introducing factors, such as regime-switching, liability, and risky asset in a mean-variance optimization framework, and also shows the relationships between a set of risk aversions and the optimal controller.



---

*Source: 5476781-2021-11-20.xml*