# Pricing Convertible Bonds with Credit Risk under Regime Switching and Numerical Solutions

Authors: Wei-Guo Zhang; Ping-Kang Liao
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

This paper discusses the convertible bonds pricing problem with regime switching and credit risk in the convertible bond market. We derive a Black-Scholes-type partial differential equation of convertible bonds and propose a convertible bond pricing model with boundary conditions. We explore the impact of dilution effect and debt leverage on the value of the convertible bond and also give an adjustment method. Furthermore, we present two numerical solutions for the convertible bond pricing model and prove their consistency. Finally, the pricing results by comparing the finite difference method with the trinomial tree show that the strength of the effect of regime switching on the convertible bond depends on the generator matrix or the regime switching strength.


---

## Body

## 1. Introduction

Convertible bond is a kind of the most important financing instruments, so the convertible bond market occupies an important position in the international financial market. American convertible bond market is the largest market in the world, and it has issued more than 400 trillion dollars in total from 1980 to 2011. Hong Kong is an international financial center; there are about 70 trillion yuan of convertible bonds in 2012. Although the amount of convertible bonds issued in developing countries is much less than that in developed countries, the convertible bonds markets of some countries are developing rapidly. For example, China issued more than 60 trillion yuan of new convertible bonds in 2010, which is almost three times the level of four years ago. Because of the importance of the convertible bonds in the financial market, the convertible bond pricing problem is a hot topic.Taking the corporate value as a basic variable, Ingersoll [1] constructed a structural model for pricing convertible bonds by deriving a Black-Scholes-type partial differential equation based on Black-Scholes’ theory. Following the work of Ingersoll [1], Brennan and Schwartz [2] explored the valuation of the convertible bonds with dividends and callable provision. However, the structural approach has a shortcoming which is the fact that the data of the corporate value is difficult to measure and observe. To overcome this shortcoming, McConnell and Schwartz [3] firstly selected the stock price as the basic variable, which becomes a mainstream later. With further research, the clauses of convertible bond also are in-depth studied. Kimura and Shinohara [4] and Yang et al. [5] explored the effect of reset clause on noncallable convertible bonds; they also derived an exact solution on the valuation of the convertible bonds with and without dilution effect. The numerical solutions are also adapted to price the convertible bonds. Brennan and Schwartz [6] took stochastic interest rate into account firstly and built the so-called two-factor model for convertible bond pricing that is popular since it is available. Then, a three-factor model was established by Davis and Lischka [7] to value the convertible bonds with stochastic credit risk. To solve the complex multifactor models, characteristics/finite elements and trinomial tree model are applied to the valuation of the convertible bonds such as Barone-Adesi et al. [8] and Xu [9]. Some other researches consider special convertible bonds or other factors of the convertible bonds. Yagi and Sawaki [10] proposed a valuation model of callable-puttable reverse convertible bonds, which are issued by a company to be exchanged for the shares of another company. Instead of the popular geometric Brownian motion model, Labuschagne and Offwood [11] valued the convertible bonds with the CGMY stock price process. Lee and Yang [12] presented the unexplained portion of the valuation model of convertible bonds based on MCDM.Credit risk is another important factor that affects the value of convertible bonds. As early as the structural approach is available, the credit risk is considered by comparing the corporate value and the convertible bond value. However, this idea is not popular for the shortcoming of the structural approach. There are two main ways to deal with the credit risk in the pricing models now. One way measures the credit risk by the credit spread of bonds. McConnell and Schwartz [3] applied this idea to the valuation of the convertible bonds. Tsiveriotis and Fernandes [13] innovatively defined the “cash only part of the convertible bonds,” whose discount rate contains credit spread and is different from the rest of the convertible bonds. This idea is also applied in binomial tree, trinomial tree, and multifactor model of the convertible bond pricing with credit risk, such as the work of Bardhan et al. [14]. The other way measures the credit risk by the default intensity. The default intensity to price convertible bonds with credit risk was introduced to measure the credit risk in the early work of Duffie and Singleton [15] and Takahashi et al. [16]. However, these works are not reasonable enough for they cannot measure the credit risk accurately so that they need to be improved. Thus, Ayache et al. [17] proposed a new model for the convertible bonds pricing with credit risk (AFV model). Different from previous work, they adopt default intensity to measure the credit risk and build a model by deriving the PDE of the convertible bonds. This idea is also introduced to trinomial tree and multifactor model, such as the work of Chambers and Liu [18]. Milanov et al. [19] set a binomial tree model for the valuation of convertible bond and prove that their model converges in continuous time to the AFV model.Since the early 20th century, the global economy interchanges between the boom and the bust frequently. Similarly, the security market also interchanges between the “bull market” and the “bear market” frequently. The stock price processes have different characteristics in different states of economy. Considering the regime switching, some options pricing models with regime switching are proposed. The options pricing with regime switching can be traced back to the work of Naik [20], who establishes a formula of option pricing with only two states. Then, Buffington and Elliot [21] built pricing models of European options and American options with regime switching that the number of states is uncertain but finite. The valuation of exotic options was also studied by Boyle and Draviam [22] with finite difference method. Yuen and Yang [23] proposed a modified trinomial tree model to price complex options with regime switching. The valuation of the currency options with stochastic volatility and stochastic interest rate under regime switching was studied by Siu et al. [24] later. Goutte [25] considered general regime switching stochastic volatility models where both the asset and the volatility dynamics depend on the values of a Markov jump process and obtain pricing and hedging formulae by risk minimization approach. As a kind of medium-term or long-term securities, the value of the convertible bonds may be affected by the regime switching. The regime switching is also considered in the convertible bond pricing. Song et al. [26] developed a valuation model for a perpetual convertible bond by the valuation model for the perpetual American option. However, they do not take credit risk, dilution effect, and debt leverage into account, and the perpetual convertible is not the most common kind of convertible bond.The value of the convertible bond is affected by credit risk, market environment, dilution, and so forth. As we know, the financial leverage of the corporate has also effect on the return of shareholders, which will affect the convertible bonds finally. This paper will discuss the pricing convertible bonds with credit risk under regime switching, which take the dilution effect and the debt leverage into account. The rest of this paper is organized as follows. Section2 derives the Black-Scholes-type partial differential equation of the convertible bonds and proposes the theoretical model. Section 3 considers the debt leverage and gives a modification of the historical volatility. Section 4 introduces main details of the numerical solutions and proves their consistency. Section 5 presents a numerical analysis for the features of the convertible bonds. Finally, the main conclusions are given in Section 6.


## 2. Theoretical Model of the Convertible Bond Pricing

In this section, we derive the Black-Scholes-type partial differential equation of the convertible bonds and try to build up the pricing model. Firstly, we give the dynamic process of the stock priceS and the process of the bond value B. Here we assume that the states of the economy are modeled by a continuous-time, finite-state, observable Markov chain X
=
{
X
t
:
t
≥
0
} on a complete probability space (
Ω
,
F
,
P
) with a finite-state space χ
=
{
e
1
,
e
2
,
…
,
e
L
}, where unit vector e
k
=
(
0
,
…
,
0,1
,
0
,
…
0
)
′
∈
R
L represents the kth regime or the kth state. Let A
=
{
a
i
j
}
i
,
j
=
1,2
,
…
,
L be the generator matrix of the Markov chain  X, which satisfies the following formula:
(1)
∑
j
=
1
L
a
i
j
=
0
,
i
=
1,2
,
…
,
L
.
Then, the semimartingale decomposition of continuous-time Markov chain X is given by
(2)
d
X
t
=
A
X
t
d
t
+
d
M
t
,
where M
=
{
M
t
:
t
≥
0
} is an R
L-valued martingale.Letr
=
{
r
t
∣
t
≥
0
} be the process of risk-free rate. In the market with regime switching, the risk-free rate only depends on the current state of the market, so the risk-free rate at time t can be expressed by inner product form as
(3)
r
t
=
r
(
X
t
)
=
〈
r
,
X
t
〉
,
where r
=
(
r
1
,
r
2
,
…
,
r
L
)
∈
R
L. Similarly, let μ
=
{
μ
t
:
t
≥
0
}, q
=
{
q
t
:
t
≥
0
}, and σ
=
{
σ
t
:
t
≥
0
} be the appreciation rate, the continuous dividend rate, and the volatility of underlying stock, respectively. They can be expressed as
(4)
μ
t
=
μ
(
X
t
)
=
〈
μ
,
X
t
〉
,
q
t
=
q
(
X
t
)
=
〈
q
,
X
t
〉
,
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
σ
t
=
σ
(
X
t
)
=
〈
σ
,
X
t
〉
.
Here μ
=
(
μ
1
,
μ
2
,
…
,
μ
L
)
∈
R
L, q
=
(
q
1
,
q
2
,
…
,
q
L
)
∈
R
L, and σ
=
(
σ
1
,
σ
2
,
…
,
σ
L
)
∈
R
L.Continuous processB
=
{
B
t
:
t
≥
0
} denotes the process of bond value without risk as
(5)
d
B
t
=
r
t
B
t
d
t
.
Following the assumptions at the first of this part, the stock price process with initial value S
0 satisfies the equation
(6)
d
S
t
=
(
μ
t
-
q
t
)
S
t
d
t
+
σ
t
S
t
d
W
t
,
where W
=
{
W
t
}
t
∈
[
0
,
T
] is a standard Brownian motion.We consider a convertible bond with maturityT and face value F. Without losing generality, we assume that the convertible bond has the following clauses. (
1
) The convertible bond pays a fixed coupon amount c
i at time t
c
i  (
i
=
1,2
,
…
,
M
). (
2
) The convertible bond can be converted to κ share after time T
con  (
0
≤
T
con
≤
T
). (
3
) The convertible bond is callable by the issuers at an interval price B
c after time T
call  (
T
con
≤
T
call
≤
T
). (
4
) The convertible bond is puttable by the holders at a price B
p after time T
put  (
T
con
≤
T
put
≤
T
). Let V
k
(
k
=
1,2
,
…
,
L
) be the price of convertible bond at the kth regime, and set V
=
(
V
1
,
V
2
,
…
,
V
L
).We now consider the valuation of the convertible bonds with credit risk under regime switching. When the company goes bankrupt, the holders of the convertible bonds can get some compensation, which is usually less thanF. What is more, if all assets of the company are more valuable than all debt or the company is merged or reorganized, the stock of the company is still valuable. Therefore, we make the following assumptions.(1)
The default probability approximately is equal toλ
t
d
t for d
t, where λ
t
=
λ
(
X
t
)
=
〈
λ
,
X
t
〉 and λ
=
(
λ
1
,
λ
2
,
…
,
λ
L
).
(2)
Upon default, the holders can choose to receive the compensationR
t
F per convertible bond or convert to κ
(
1
-
η
t
) shares. Here R
t
=
R
(
X
t
)
=
〈
R
,
X
t
〉 is the recovery factor and η
t
=
η
(
X
t
)
=
〈
η
,
X
t
〉 is the jump ratio of underlying stock. Set R
=
(
R
1
,
R
2
,
…
,
R
L
) and η
=
(
η
1
,
η
2
,
…
,
η
L
).In the following discussion, we try to establish the pricing model of convertible bonds by constructing a hedging portfolio based on no arbitrage principle. We assume that the hedging portfolio∏ includes a convertible bond V
t and -
β
t shares S
t at time t. When there is no coupon payment in [
t
,
t
+
d
t
], the value of hedging portfolio at time t is
(7)
∏
t
=
V
t
-
β
t
S
t
.
Absent of default, the change in value of the hedging portfolio after time d
t is
(8)
d
V
-
β
t
d
S
.
Upon default, the loss of the hedging portfolio after time d
t is
(9)
(
V
t
-
β
t
S
t
)
-
(
max
⁡
{
κ
(
1
-
η
t
)
S
t
,
R
t
B
}
-
β
t
(
1
-
η
t
)
S
t
)
=
V
t
-
β
t
η
t
S
t
-
max
⁡
{
κ
(
1
-
η
t
)
S
t
,
R
t
B
}
.
As the default probability is λ
t
d
t and no default probability is 1
-
λ
t
d
t in time d
t, the expected change in value of the hedging portfolio after time d
t is the weighted average of the hedging portfolio value change upon default and absent of default as
(10)
d
∏
t
=
(
1
-
λ
t
d
t
)
(
d
V
-
β
t
d
S
)
-
λ
t
d
t
(
V
t
-
β
t
η
t
S
t
-
max
⁡
{
κ
(
1
-
η
t
)
S
t
,
R
t
B
}
)
.
Applying Ito formula with regime switching to (10), we have
(11)
d
∏
t
=
[
∂
V
∂
t
+
(
μ
t
-
q
t
)
S
t
∂
V
∂
S
+
1
2
σ
t
2
S
t
2
∂
2
V
∂
S
2
-
β
μ
t
S
t
]
d
t
+
λ
t
d
t
(
V
t
-
β
t
η
t
S
t
)
+
λ
t
d
t
max
⁡
{
κ
(
1
-
η
t
)
S
t
,
R
t
B
}
+
(
μ
t
S
t
∂
V
∂
S
-
β
t
μ
t
S
t
)
d
W
t
+
〈
V
,
d
X
t
〉
.
Following the no arbitrage principle, we have
(12)
d
∏
t
=
r
t
∏
t
d
t
.
That is,
(13)
r
t
(
V
t
-
β
t
S
t
)
d
t
=
[
∂
V
∂
t
+
(
μ
t
-
q
t
)
S
t
∂
V
∂
S
+
1
2
σ
t
2
S
t
2
∂
2
V
∂
S
2
-
β
μ
t
S
t
+
λ
t
(
V
t
-
β
t
η
t
S
t
)
+
λ
t
max
⁡
⁡
{
κ
(
1
-
η
t
)
S
t
,
R
t
B
}
+
〈
V
,
A
X
t
〉
1
2
]
d
t
+
(
μ
t
S
t
∂
V
∂
S
-
β
t
μ
t
S
t
)
d
W
t
+
〈
V
,
d
M
t
〉
.
Similar to Buffington and Elliot [21], all the terms in the right of (13) must be equal to the left of equation; that is
(14)
μ
t
S
t
∂
V
∂
S
-
β
t
μ
t
S
t
=
0
or
∂
V
∂
S
=
β
t
.
Consider
(15)
r
t
(
V
t
-
β
t
S
t
)
d
t
=
∂
V
∂
t
+
(
μ
t
-
q
t
)
S
t
∂
V
∂
S
+
1
2
σ
t
2
S
t
2
∂
2
V
∂
S
2
-
β
μ
t
S
t
+
λ
t
(
V
t
-
β
t
η
t
S
t
)
+
λ
t
max
⁡
{
κ
(
1
-
η
t
)
S
t
,
R
t
B
}
+
〈
V
,
A
X
t
〉
.
Substituting (14) into (15), we obtain
(16)
∂
V
∂
t
+
(
r
t
-
q
t
+
λ
t
η
t
)
S
t
∂
V
∂
S
+
1
2
σ
t
2
S
t
2
∂
2
V
∂
S
2
-
(
r
t
+
λ
t
)
V
t
+
〈
V
,
A
X
〉
+
λ
t
max
⁡
{
κ
(
1
-
η
t
)
S
t
,
R
t
B
}
=
0
.
Equation (16) is partial difference equation or the so-called Black-Scholes-type partial differential equation of the convertible bond price with credit risk under regime switching.
(17)
r
k
=
〈
r
,
e
k
〉
,
q
k
=
〈
q
,
e
k
〉
,
σ
k
=
〈
σ
,
e
k
〉
,
λ
k
=
〈
λ
,
e
k
〉
,
η
k
=
〈
η
,
e
k
〉
,
R
k
=
〈
R
,
e
k
〉
,
V
k
=
V
k
,
t
=
〈
V
,
e
k
〉
,
k
=
1,2
,
…
,
L
.
Equation (16) can be rewritten as
(18)
∂
V
k
∂
t
+
(
r
k
-
q
k
+
λ
k
η
k
)
S
t
∂
V
k
∂
S
+
1
2
σ
k
2
S
t
2
∂
2
V
k
∂
S
2
-
(
r
k
+
λ
k
)
V
k
+
∑
l
=
1
N
a
k
l
V
l
+
λ
k
max
⁡
{
κ
S
t
(
1
-
η
k
)
,
R
k
B
}
=
0
,
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
k
=
1,2
,
…
,
L
.We must notice that (18) is only applicable during two coupon payments. When the coupon payment is paid, the price per convertible bond will be reduced by the amount of coupon at once. Therefore, the price of the convertible bonds at time t
c
i satisfies
(19)
V
t
c
i
+
=
V
t
c
i
-
-
c
i
,
i
=
1,2
,
…
,
M
-
1
.When all investors convert their convertible bonds, some liabilities become the owner’s equity and the number of shares increases. We assume that the company issuesm stocks and n convertible bonds. If all convertible bonds are converted at time t, owner’s equity value is m
S
t
+
n
F totally. As the number of stocks becomes m
+
n
κ, the value per stock is (
m
S
t
+
n
F
)
/
(
m
+
n
κ
), which is less than S
t when S
t
>
F
/
κ. This is the so-called dilution effect in convertible bonds. So the investor will get κ
(
m
S
t
+
n
F
)
/
(
m
+
n
κ
) rather than κ
S
t per convertible bond when the convertible bonds are converted. Thus, the dilution effect can be adjusted by changing the payment of the convertible bonds.A rational investor seeks to maximize the value of the convertible bond at any time. The value of a convertible bond must be greater than or equal to its conversion value during the conversion time. Thus, we get the conversion provision:(20)
V
k
,
t
≥
(
m
S
t
+
n
F
)
κ
m
+
n
κ
,
t
∈
[
T
con
,
T
]
.
At last time T, the company will give a return which is more than or equal to face return to call the convertible bonds. Thus, at time T, the value of convertible bonds must be
(21)
V
k
,
T
=
max
⁡
{
(
m
S
T
+
n
B
)
κ
m
+
n
κ
,
F
+
c
M
}
.
Thus, we get the terminal condition of the partial difference equation (18). If the convertible bonds are puttable at time t, the investor can get at least B
p for par convertible bond. Thus, we get the put provision as
(22)
V
k
,
t
≥
B
p
,
t
∈
[
T
put
,
T
]
.
Similarly, the value of a convertible bond must be less than the maximum of conversion value and B
c when the convertible bonds are callable at time t, or the convertible bonds will be called by the company. This call provision can be written as
(23)
V
k
,
t
≤
max
⁡
{
(
m
S
t
+
n
F
)
κ
m
+
n
κ
,
B
c
}
,
t
∈
[
T
call
,
T
]
.Considering constraints of the convertible bonds including conversion provision, put provision, and call provision, we get the pricing model of convertible bonds with credit risk under regime switching as follows:(24)
∂
V
k
∂
t
+
(
r
k
-
q
k
+
λ
k
η
k
)
S
t
∂
V
k
∂
S
+
1
2
σ
k
2
S
t
2
∂
2
V
k
∂
S
2
-
(
r
k
+
λ
k
)
V
k
+
∑
l
=
1
N
a
k
l
V
l
+
λ
k
max
⁡
{
κ
S
t
(
1
-
η
k
)
,
R
k
F
}
=
0
V
k
,
T
=
max
⁡
{
(
m
S
T
+
n
B
)
κ
m
+
n
κ
,
F
+
c
M
}
V
k
,
t
≥
(
m
S
t
+
n
F
)
κ
m
+
n
κ
,
t
∈
[
T
con
,
T
]
V
k
,
t
≥
B
p
,
t
∈
[
T
put
,
T
]
V
k
,
t
≤
max
⁡
{
(
m
S
t
+
n
F
)
κ
m
+
n
κ
,
B
c
}
,
t
∈
[
T
call
,
T
]
V
k
,
t
c
i
+
=
V
k
,
t
c
i
-
-
c
i
,
i
=
1,2
,
…
,
M
-
1
.In fact, the model (24) is a boundary value problem. We can solve this boundary value problem by using the theory of partial difference equation. But the boundary problem is difficult that it has no explicit solution. What is more, if we set λ
1
=
λ
2
=
⋯
=
λ
L
=
0, the model (24) becomes the pricing model of the convertible bonds without credit risk. Hence, we can treat the pricing model without credit risk as a special case of the model (24).


## 3. Dilution Effect and Debt Leverage

As King [27] showed, dilution and leverage have effect on the value of convertible bond. Thus, the dilution effect and debt leverage will be considered in the pricing model in this section. The dilution effect of the convertible bonds stems from the issuance of new shares when the convertible bonds are converted. In the model (24), the dilution effect has been taken into account. Thus, we only consider the debt leverage in the following discussion. As we know, the debt leverage, or financial leverage, is often used as an indicator of financial risk. So once the debt leverage changes, the risk of shareholders’ returns may be changed. Hence, the volatility of the shares returns will be affected by debt leverage. This inspires us to modify the volatility to reflect the effect of the debt leverage on the convertible bond value. Thus, the volatility in the pricing model of the convertible bonds should be the volatility after the convertible bond issues or modification of historical volatility before convertible bond issues.MM theory about the impact of the financial leverage on the capital cost and the risk provides us with a good idea to deal with the debt leverage. For simplicity, we add some necessary assumptions based on the Black-Scholes theory as follows.(1)
No convertible bond issuance cost. The convertible bond issues at par value.
(2)
There are no new bonds, the convertible bonds and the warrants issuing during the existence of the convertible bonds. The level of debt is stable during the existence of the convertible bonds.
(3)
The operational efficiency risk is unchanged during the existence of the convertible bonds, and the return and risk of the capital are the same as before.
(4)
The return rate and volatility of the stocks are the same as that of the ownership.Following the assumptions and MM theory, we have the modification of the historical volatility of the stocks as Proposition1.Proposition 1.
There arem stocks and n convertible bonds issued by the company right now. Let E, D
0, and D
1 be the market value of the company owner’ equity, the market value of corporate liabilities, and the total face value of the convertible bonds. Let σ
k
,
0 and σ
k
,
1(k
=
1,2
,
…
,
L) be the volatility of the stocks at kth regime before and after the convertible bonds issue. Let ρ be income tax rate of enterprise. Then
(25)
σ
k
,
1
=
E
+
(
D
0
+
D
1
)
(
1
-
ρ
)
E
+
D
0
(
1
-
ρ
)
σ
k
,
0
=
m
S
0
+
(
D
0
+
m
F
)
(
1
-
ρ
)
m
S
0
+
D
0
(
1
-
ρ
)
σ
k
,
0
.Proof.
Letr
k
,
0 and r
k
,
1 be the returns of owner’s equity at the kth regime. Let r
k
,
u and σ
k
,
u be the corporate returns and the volatility at the kth regime without any debt. Let r
k
,
d be the average rate of the corporate debt at the kth regime. Because there is not any income except the convertible bonds when the convertible bonds issue, the owner’s equity is unchanged. Then the relationship between the rate of shareholders returns with and without debt at the kth regime satisfies
(26)
r
k
,
0
=
r
k
,
u
+
D
0
E
(
r
k
,
u
-
r
k
,
d
)
(
1
-
ρ
)
=
E
+
D
0
(
1
-
ρ
)
E
r
k
,
u
-
D
0
(
1
-
ρ
)
E
r
k
,
d
,
k
=
1,2
,
…
,
L
.
Because r
k
,
d is constant when the debt is known, taking the variance of both sides of (26), we have
(27)
σ
k
,
0
2
=
Var
⁡
[
r
k
,
0
]
=
(
E
+
D
0
(
1
-
ρ
)
E
)
2
Var
⁡
[
r
k
,
u
]
=
(
E
+
D
0
(
1
-
ρ
)
E
σ
k
,
u
)
2
.
Similarly, we can derive the relationship between the volatility before and after the convertible bonds issue as follows:
(28)
σ
k
,
1
2
=
Var
⁡
[
r
k
,
1
]
=
(
E
+
(
D
0
+
D
1
)
(
1
-
ρ
)
E
σ
k
,
u
)
2
=
(
E
+
(
D
0
+
D
1
)
(
1
-
ρ
)
E
+
D
0
(
1
-
ρ
)
σ
k
,
0
)
2
.
Comparing (27) and (28), we deduce
(29)
σ
k
,
1
=
E
+
D
0
+
D
1
E
+
D
0
σ
k
,
0
=
m
S
0
+
(
D
0
+
n
F
)
(
1
-
ρ
)
m
S
0
+
D
0
(
1
-
ρ
)
σ
k
,
0
.
This completes the proof.Following Proposition1, we can easily get modified volatility of the model (24) when debt leverage is considered. But it is worth noting that the method is only applicable for historical volatility, not for implied volatility and other volatilities, which are the current volatilities after the convertible bonds issue so that the debt leverage is concluded and the volatilities are not necessary to be modified.


## 4. Numerical Solutions

The pricing model of the convertible bonds with credit risk under regime switching introduced in Section2 is complex and hard to achieve an exact solution, so we need the numerical solution method to solve this model. The main numerical solution methods of convertible bonds pricing include the finite difference method, the finite element method, the binomial or the trinomial tree, and the Monte Carlo simulation. The finite difference method is a basic one that solves directly the partial difference equation with the boundary conditions. The tree model is another basic one that solves the pricing problem based on the process of stock price. What is more, without regime switching, it is known that the tree model is consistent with the finite difference method. Here we will choose the finite difference method and the tree model to solve the pricing problem. In this section we will give the details of the finite difference and the trinomial tree to solve this pricing model and explore their consistency.


### 4.1. The Finite Difference Scheme

The finite difference method is the most basic numerical solution of convertible bond pricing by solving the partial difference equation. We can improve the accuracy and get satisfactory result by increasing time steps and price steps. For convergence effect and convergence rate of the finite difference method, we apply the Crank-Nicolson method to solve the model (18) and get approximate solution. Divide the area [
t
,
T
]
×
[
0
,
S
max
⁡
] into M
×
N small identical rectangles. Set Δ
t
=
(
T
-
t
)
/
M and Δ
S
=
S
max
⁡
/
N. Let V
k
,
i
,
j be the price of the convertible bond under the kth regime, time i
Δ
t, and stock pricej
Δ
t. To get the Euler explicit scheme, we write the difference approximation as
(30)
∂
V
k
∂
t
=
V
k
,
i
+
1
,
j
-
V
k
,
i
,
j
Δ
t
,
∂
V
k
∂
S
=
V
k
,
i
,
j
+
1
-
V
k
,
i
,
j
-
1
2
Δ
S
,
i
i
i
i
i
i
i
i
i
i
i
i
∂
V
k
2
∂
2
S
=
V
k
,
i
,
j
+
1
+
V
k
,
i
,
j
-
1
-
2
V
k
,
i
,
j
Δ
S
2
.
Substituting (30) in (18), we have
(31)
V
k
,
i
+
1
,
j
-
V
k
,
i
,
j
Δ
t
+
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
S
V
k
,
i
,
j
+
1
-
V
k
,
i
,
j
-
1
2
Δ
S
+
1
2
σ
k
2
j
2
Δ
S
2
V
k
,
i
,
j
+
1
+
V
k
,
i
,
j
-
1
-
2
V
k
,
i
,
j
Δ
S
2
-
(
r
k
+
λ
k
)
V
k
,
i
,
j
+
λ
k
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
j
Δ
S
)
+
∑
l
=
1
L
a
k
l
V
l
,
i
,
j
=
0
.
To get the Euler implicit scheme, we write the difference approximation as
(32)
∂
V
k
∂
t
=
V
k
,
i
+
1
,
j
-
V
k
,
i
,
j
Δ
t
,
∂
V
k
∂
S
=
V
k
,
i
+
1
,
j
+
1
-
V
k
,
i
+
1
,
j
-
1
2
Δ
S
,
i
i
i
i
i
i
i
i
i
i
i
i
∂
V
k
2
∂
2
S
=
V
k
,
i
+
1
,
j
+
1
+
V
k
,
i
+
1
,
j
-
1
-
2
V
k
,
i
+
1
,
j
Δ
S
2
.
Substituting (32) in (18), we get
(33)
V
k
,
i
+
1
,
j
-
V
k
,
i
,
j
Δ
t
+
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
S
V
k
,
i
+
1
,
j
+
1
-
V
k
,
i
+
1
,
j
-
1
2
Δ
S
+
1
2
σ
k
2
j
2
Δ
S
2
V
k
,
i
+
1
,
j
+
1
+
V
k
,
i
+
1
,
j
-
1
-
2
V
k
,
i
+
1
,
j
Δ
S
2
-
(
r
k
+
λ
k
)
V
k
,
i
,
j
+
λ
k
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
j
Δ
S
)
+
∑
l
=
1
L
a
k
l
V
l
,
i
,
j
=
0
.
By (31) and (33), we obtain
(34)
V
k
,
i
+
1
,
j
-
V
k
,
i
,
j
Δ
t
+
(
r
k
-
q
k
+
λ
k
η
k
)
×
j
V
k
,
i
,
j
+
1
-
V
k
,
i
,
j
-
1
+
V
k
,
i
+
1
,
j
+
1
-
V
k
,
i
+
1
,
j
-
1
4
-
(
r
k
+
λ
k
)
V
k
,
i
,
j
+
∑
l
=
1
L
a
k
l
V
l
,
i
,
j
+
1
4
σ
k
2
j
2
(
V
k
,
i
+
1
,
j
+
1
+
V
k
,
i
+
1
,
j
-
1
-
2
V
k
,
i
+
1
,
j
+
V
k
,
i
,
j
+
1
+
V
k
,
i
,
j
-
1
-
2
V
k
,
i
,
j
)
+
λ
k
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
j
Δ
S
)
=
0
.
Multiplying Δ
t at both sides of (34), transposing, and simplifying, we get
(35)
[
1
4
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
t
-
1
4
σ
k
2
j
2
Δ
t
]
V
k
,
i
,
j
-
1
+
[
1
+
1
4
(
r
k
+
λ
k
)
j
Δ
t
+
1
2
σ
k
2
j
2
Δ
t
]
V
k
,
i
,
j
+
[
-
1
4
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
t
-
1
4
σ
k
2
j
2
Δ
t
]
V
k
,
i
,
j
+
1
-
∑
l
=
1
L
a
k
l
V
l
,
i
,
j
=
[
-
1
4
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
t
+
1
4
σ
k
2
j
2
Δ
t
]
V
k
,
i
+
1
,
j
-
1
+
[
1
-
1
2
σ
k
2
j
2
Δ
t
]
V
k
,
i
+
1
,
j
+
[
1
4
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
t
+
1
4
σ
k
2
j
2
Δ
t
]
V
k
,
i
+
1
,
j
+
1
+
λ
k
Δ
t
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
j
Δ
S
)
.
Solving (35) from time (
N
-
1
)
Δ
t to time 0 recursively and taking into account the constraints of the convertible bonds, we can achieve the price of the convertible bonds at any time and any stock price. Furthermore, there are always an upper bound and a lower bound in the finite difference model. Here we set
(36)
V
k
,
i
,
N
+
1
=
κ
N
Δ
S
,
V
k
,
i
,
1
=
0
,
k
=
1,2
,
…
,
L
,
i
=
1,2
,
…
,
M
.
Early studies have shown that the Crank-Nicolson scheme is second-order convergence which refers to time. Therefore, the convertible bond price calculated by the finite difference converges to the theoretical value when the interval Δ
t tends to infinity small.


### 4.2. The Trinomial Tree Model

The trinomial tree model is another important numerical solution and is widely used in derivatives pricing. The trinomial tree method has advantage in computational efficiency to the finite difference method. In the following section we are going to apply the trinomial tree to the valuation of convertible bonds. We divide interval[
t
,
T
] into M small identical intervals and set Δ
t
=
(
T
-
t
)
/
M. When the market state is in the kth regime, the default probability is
(37)
w
k
=
1
-
e
-
λ
k
Δ
t
=
λ
k
Δ
t
+
o
(
Δ
t
)
≈
λ
k
Δ
t
,
k
=
1,2
,
…
,
L
.
Following the trinomial tree of Boyle [28], we set
(38)
u
k
=
e
α
k
σ
k
Δ
t
,
d
k
=
1
u
k
=
e
-
α
k
σ
k
Δ
t
,
m
k
=
1
,
α
k
>
0
,
k
=
1,2
,
…
,
L
.
Let π
k
,
u
,
π
k
,
m
,
π
k
,
d  (
k
=
1,2
,
…
,
L
) denote the probability when the stock price increases, remains, and decreases, respectively. Considering the expected return and volatility of stock price under the risk-neutral market and applying the risk-free arbitrage principle, these jump probabilities satisfy
(39)
π
k
,
u
+
π
k
,
m
+
π
k
,
d
=
1
,
π
k
,
u
u
k

+

π
k
,
m

+

π
k
,
d
d
k

=

e
(
r
k
+
λ
k
η
k
-
q
k
)
Δ
t
,
π
k
,
u
u
k
2

+

π
k
,
m

+

π
k
,
d
d
k
2

=

e
2
(
r
k
+
λ
k
η
k
-
q
k
)
Δ
t
+
σ
k
2
Δ
t
,
k

=

1,2
,
…
,
L
.
Solving (39), we have
(40)
π
k
,
u
=
(
U
k
-
Q
k
)
-
u
k
(
Q
k
-
1
)
(
u
k
-
d
k
)
(
1
-
d
k
)
,
π
k
,
d
=
(
U
k
-
Q
k
)
-
d
k
(
Q
k
-
1
)
(
u
k
-
d
k
)
(
u
k
-
1
)
,
π
k
,
m
=
1
-
π
k
,
u
-
π
k
,
d
,
k
=
1,2
,
…
,
L
,
where U
k
=
e
(
r
k
+
λ
k
η
k
-
q
k
)
Δ
t
,
Q
k
=
e
2
(
r
k
+
λ
k
η
k
-
q
k
)
Δ
t
+
σ
k
2. Then, we get the recurrence relation equations as
(41)
V
(
k
,
S
,
t
)
=
∑
l
=
1
L
p
k
l
e
-
r
l
Δ
t
{
π
l
,
u
[
(
1
-
w
l
)
V
(
l
,
u
l
S
,
t
+
Δ
t
)
+
w
l
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
u
l
S
)
]
+
π
l
,
u
[
(
1
-
w
l
)
V
(
l
,
m
l
S
,
t
+
Δ
t
)
+
w
l
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
m
l
S
)
]
+
π
l
,
u
[
(
1
-
w
l
)
V
(
l
,
d
l
S
,
t
+
Δ
t
)
+
w
l
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
d
l
S
)
]
}
=
∑
l
=
1
L
p
k
l
e
-
r
l
Δ
t
{
(
1
-
w
l
)
[
π
l
,
u
V
(
l
,
u
l
S
,
t
+
Δ
t
)
+
π
l
,
m
V
(
l
,
m
l
S
,
t
+
Δ
t
)
+
π
l
,
d
V
(
l
,
d
l
S
,
t
+
Δ
t
)
]
+
w
l
[
π
l
,
u
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
u
l
S
)
+
π
l
,
m
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
m
l
S
)
+
π
l
,
d
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
d
l
S
)
]
}
.
We can set the parameter α
k to decrease nodes of the trinomial tree. Thus, we set
(42)
α
k
=
σ
0
σ
k
,
k
=
1,2
,
…
,
L
,
σ
0
>
max
⁡
1
≤
k
≤
L
(
σ
k
)
.
That is,
(43)
u
k
=
e
σ
0
Δ
t
,
d
k
=
e
-
σ
0
Δ
t
,
m
k
=
1
,
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
k
=
1,2
,
…
,
L
.
Yuen and Yang [23] proposed that σ
0 should be taken as
(44)
σ
0
=
max
⁡
1
≤
k
≤
L
(
σ
k
)
+
(
1.5
-
1
)
mean
1
≤
k
≤
L
(
σ
k
)
.
Thus, the number of the nodes of the trinomial tree decreases so much that we can save many storage spaces and much computation time.


### 4.3. The Consistency of the Finite Difference and the Trinomial Tree

As we know, without regime switching, the tree model is consistent with the finite difference method, so we are going to explore their consistency with regime switching. Starting with the trinomial tree, we can achieve the relation between the trinomial tree model and the partial differential equations as the following proposition.Proposition 2.
When the intervalΔ
t tends to zero, (41) converges to (18).Proof.
For anyb
>
0, applying the multivariate Taylor’s series, we have
(45)
V
(
k
,
b
S
,
t
)
=
V
(
k
,
S
+
(
b
-
1
)
S
,
t
+
Δ
t
)
=
V
(
k
,
S
,
t
)
+
(
b
-
1
)
S
∂
V
∂
S
+
Δ
t
∂
V
∂
t
+
1
2
(
b
-
1
)
2
S
2
∂
V
∂
S
+
o
(
Δ
t
)
.
Then, we obtain
(46)
π
k
,
u
V
(
k
,
u
k
S
,
t
+
Δ
t
)
+
π
k
,
m
V
(
k
,
m
k
S
,
t
+
Δ
t
)
+
π
k
,
d
V
(
k
,
d
k
S
,
t
+
Δ
t
)
=
(
π
k
,
u
+
π
k
,
m
+
π
k
,
d
)
(
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
)
+
(
(
π
k
,
u
u
k
+
π
k
,
m
m
k
+
π
k
,
d
d
k
)
-
(
π
k
,
u
+
π
k
,
m
+
π
k
,
d
)
)
S
∂
V
∂
S
+
1
2
(
(
π
k
,
u
+
π
k
,
m
+
π
k
,
d
)
1
2
(
π
k
,
u
u
k
2
+
π
k
,
m
m
k
2
+
π
k
,
d
d
k
2
)
-
2
(
π
k
,
u
u
k
+
π
k
,
m
m
k
+
π
k
,
d
d
k
)
+
(
π
k
,
u
+
π
k
,
m
+
π
k
,
d
)
1
2
)
S
2
∂
V
∂
S
+
o
(
Δ
t
)
.
=
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
e
(
r
k
+
λ
k
η
k
)
Δ
t
-
1
)
S
∂
V
∂
S
+
1
2
(
e
2
(
r
k
+
λ
k
η
k
)
Δ
t
+
σ
k
2
Δ
t
-
2
e
(
r
k
+
λ
k
η
k
)
Δ
t
+
1
)
×
S
2
∂
V
∂
S
+
o
(
Δ
t
)
=
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
k
2
S
2
Δ
t
∂
V
∂
S
+
o
(
Δ
t
)
.
Multiplying ∏
h
=
1
L
e
r
l
Δ
t both sides of (41), we have
(47)
(
∏
h
=
1
L
e
r
l
Δ
t
)
V
(
k
,
S
,
t
)
=
∑
l
=
1
L
p
k
l
(
∏
h
≠
l
e
r
h
Δ
t
)
×
{
(
1
-
w
l
)
[
π
l
,
u
V
(
l
,
u
l
S
,
t
+
Δ
t
)
IIIIIIIIIIIIIIIIIII
+
π
l
,
m
V
(
l
,
m
l
S
,
t
+
Δ
t
)
IIIIIIIIIIIIIIIIIII
+
π
l
,
d
V
(
l
,
d
l
S
,
t
+
Δ
t
)
]
+
w
l
[
π
l
,
u
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
u
l
S
)
+
π
l
,
m
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
m
l
S
)
+
π
l
,
d
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
d
l
S
)
]
}
.
Applying Taylor’s series, we have
(48)
(
∏
h
=
1
L
(
1
+
r
h
Δ
t
+
o
(
Δ
t
)
)
)
V
(
k
,
S
,
t
)
=
(
1
+
∑
h
=
1
L
r
h
Δ
t
+
o
(
Δ
t
)
)
V
(
k
,
S
,
t
)
=
(
1
+
a
k
k
Δ
t
+
o
(
Δ
t
)
)
(
∏
h
≠
k
(
1
+
r
h
Δ
t
+
o
(
Δ
t
)
)
)
×
{
∂
V
∂
t
(
1
-
λ
k
Δ
t
+
o
(
Δ
t
)
)
×
[
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
Δ
t
∂
V
∂
S
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
+
1
2
σ
k
2
S
2
Δ
t
∂
V
∂
S
+
o
(
Δ
t
)
]
iiiiiiiii
+
(
λ
k
Δ
t
+
o
(
Δ
t
)
)
×
[
π
k
,
u
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
u
k
S
)
+
π
k
,
m
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
m
k
S
)
+
π
k
,
d
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
d
k
S
)
]
∂
V
∂
t
}
+
∑
l
≠
k
(
a
k
l
Δ
t
+
o
(
Δ
t
)
)
(
∏
h
≠
l
(
1
+
r
h
Δ
t
+
o
(
Δ
t
)
)
)
×
{
∂
V
∂
S
(
1
-
λ
l
Δ
t
+
o
(
Δ
t
)
)
×
[
∂
V
∂
S
V
(
l
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
l
+
λ
l
η
l
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
l
2
S
2
Δ
t
∂
V
∂
S
+
o
(
Δ
t
)
]
+
(
λ
l
Δ
t
+
o
(
Δ
t
)
)
×
[
π
l
,
u
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
u
l
S
)
iiiiiiiiiiiiiiiiii
+
π
l
,
m
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
m
l
S
)
iiiiiiiiiiiiiiiiii
+
π
l
,
d
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
d
l
S
)
]
∂
V
∂
S
}
=
(
1
+
a
k
k
Δ
t
+
o
(
Δ
t
)
)
(
1
+
∑
h
≠
k
r
h
Δ
t
+
o
(
Δ
t
)
)
×
{
∂
V
∂
t
(
1
-
λ
k
Δ
t
+
o
(
Δ
t
)
)
×
[
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
k
2
S
2
Δ
t
∂
V
∂
S
+
o
(
Δ
t
)
]
+
(
λ
k
Δ
t
+
o
(
Δ
t
)
)
×
[
π
k
,
u
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
u
k
S
)
+
π
k
,
m
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
m
k
S
)
+
π
k
,
d
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
d
k
S
)
]
∂
V
∂
t
}
+
∑
l
≠
k
(
a
k
l
Δ
t
+
o
(
Δ
t
)
)
(
1
+
∑
h
≠
l
r
h
Δ
t
+
o
(
Δ
t
)
)
×
{
∂
V
∂
S
(
1
-
λ
l
Δ
t
+
o
(
Δ
t
)
)
×
[
∂
V
∂
S
V
(
l
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
l
+
λ
l
η
l
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
l
2
S
2
Δ
t
∂
V
∂
S
+
o
(
Δ
t
)
]
+
(
λ
l
Δ
t
+
o
(
Δ
t
)
)
×
[
π
l
,
u
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
u
l
S
)
iiiiiiiiiiiii
+
π
l
,
m
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
m
l
S
)
iiiiiiiiiiiii
+
π
l
,
d
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
d
l
S
)
]
∂
V
∂
t
}
=
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
k
2
S
2
Δ
t
∂
V
∂
S
-
λ
k
V
k
Δ
t
+
∑
h
≠
k
r
h
V
k
Δ
t
+
∑
l
=
1
L
a
k
l
V
k
Δ
t
+
o
(
Δ
t
)
+
λ
k
Δ
t
[
π
k
,
u
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
u
k
S
)
iiiiiiiiiiiiiiiiiii
+
π
k
,
m
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
m
k
S
)
iiiiiiiiiiiiiiiiiii
+
π
k
,
d
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
d
k
S
)
]
.
Eliminating the similar items of (48), we have
(49)
r
k
V
k
Δ
t
=
Δ
t
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
k
2
S
2
Δ
t
∂
V
∂
S
-
λ
k
V
k
Δ
t
+
∑
l
=
1
L
a
k
l
V
k
Δ
t
+
o
(
Δ
t
)
+
λ
k
Δ
t
[
π
k
,
u
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
u
k
S
)
+
π
k
,
m
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
m
k
S
)
+
π
k
,
d
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
d
k
S
)
]
.
Let both sides of (49) be divided by Δ
t. Then, we have
(50)
λ
k
V
k
=
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
∂
V
∂
S
+
1
2
σ
k
2
S
2
∂
V
∂
S
-
λ
k
V
k
+
∑
l
=
1
L
a
k
l
V
k
+
o
(
1
)
+
λ
k
[
π
k
,
u
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
u
k
S
)
+
π
k
,
m
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
m
k
S
)
+
π
k
,
d
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
d
k
S
)
]
.
Let Δ
t
→
0; (50) converges to
(51)
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
∂
V
∂
S
+
1
2
σ
k
2
S
2
∂
V
∂
S
-
(
r
k
+
λ
k
)
V
k
+
∑
l
=
1
L
a
k
l
V
k
+
λ
k
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
S
)
=
0
.
This completes the proof.Proposition2 indicates that the trinomial tree model of the convertible bonds with credit risk under regime switching is consistent with the partial differential equations of the convertible bonds. That is, the trinomial tree model of the convertible bonds is consistent with the finite difference model of the convertible bonds.


## 4.1. The Finite Difference Scheme

The finite difference method is the most basic numerical solution of convertible bond pricing by solving the partial difference equation. We can improve the accuracy and get satisfactory result by increasing time steps and price steps. For convergence effect and convergence rate of the finite difference method, we apply the Crank-Nicolson method to solve the model (18) and get approximate solution. Divide the area [
t
,
T
]
×
[
0
,
S
max
⁡
] into M
×
N small identical rectangles. Set Δ
t
=
(
T
-
t
)
/
M and Δ
S
=
S
max
⁡
/
N. Let V
k
,
i
,
j be the price of the convertible bond under the kth regime, time i
Δ
t, and stock pricej
Δ
t. To get the Euler explicit scheme, we write the difference approximation as
(30)
∂
V
k
∂
t
=
V
k
,
i
+
1
,
j
-
V
k
,
i
,
j
Δ
t
,
∂
V
k
∂
S
=
V
k
,
i
,
j
+
1
-
V
k
,
i
,
j
-
1
2
Δ
S
,
i
i
i
i
i
i
i
i
i
i
i
i
∂
V
k
2
∂
2
S
=
V
k
,
i
,
j
+
1
+
V
k
,
i
,
j
-
1
-
2
V
k
,
i
,
j
Δ
S
2
.
Substituting (30) in (18), we have
(31)
V
k
,
i
+
1
,
j
-
V
k
,
i
,
j
Δ
t
+
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
S
V
k
,
i
,
j
+
1
-
V
k
,
i
,
j
-
1
2
Δ
S
+
1
2
σ
k
2
j
2
Δ
S
2
V
k
,
i
,
j
+
1
+
V
k
,
i
,
j
-
1
-
2
V
k
,
i
,
j
Δ
S
2
-
(
r
k
+
λ
k
)
V
k
,
i
,
j
+
λ
k
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
j
Δ
S
)
+
∑
l
=
1
L
a
k
l
V
l
,
i
,
j
=
0
.
To get the Euler implicit scheme, we write the difference approximation as
(32)
∂
V
k
∂
t
=
V
k
,
i
+
1
,
j
-
V
k
,
i
,
j
Δ
t
,
∂
V
k
∂
S
=
V
k
,
i
+
1
,
j
+
1
-
V
k
,
i
+
1
,
j
-
1
2
Δ
S
,
i
i
i
i
i
i
i
i
i
i
i
i
∂
V
k
2
∂
2
S
=
V
k
,
i
+
1
,
j
+
1
+
V
k
,
i
+
1
,
j
-
1
-
2
V
k
,
i
+
1
,
j
Δ
S
2
.
Substituting (32) in (18), we get
(33)
V
k
,
i
+
1
,
j
-
V
k
,
i
,
j
Δ
t
+
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
S
V
k
,
i
+
1
,
j
+
1
-
V
k
,
i
+
1
,
j
-
1
2
Δ
S
+
1
2
σ
k
2
j
2
Δ
S
2
V
k
,
i
+
1
,
j
+
1
+
V
k
,
i
+
1
,
j
-
1
-
2
V
k
,
i
+
1
,
j
Δ
S
2
-
(
r
k
+
λ
k
)
V
k
,
i
,
j
+
λ
k
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
j
Δ
S
)
+
∑
l
=
1
L
a
k
l
V
l
,
i
,
j
=
0
.
By (31) and (33), we obtain
(34)
V
k
,
i
+
1
,
j
-
V
k
,
i
,
j
Δ
t
+
(
r
k
-
q
k
+
λ
k
η
k
)
×
j
V
k
,
i
,
j
+
1
-
V
k
,
i
,
j
-
1
+
V
k
,
i
+
1
,
j
+
1
-
V
k
,
i
+
1
,
j
-
1
4
-
(
r
k
+
λ
k
)
V
k
,
i
,
j
+
∑
l
=
1
L
a
k
l
V
l
,
i
,
j
+
1
4
σ
k
2
j
2
(
V
k
,
i
+
1
,
j
+
1
+
V
k
,
i
+
1
,
j
-
1
-
2
V
k
,
i
+
1
,
j
+
V
k
,
i
,
j
+
1
+
V
k
,
i
,
j
-
1
-
2
V
k
,
i
,
j
)
+
λ
k
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
j
Δ
S
)
=
0
.
Multiplying Δ
t at both sides of (34), transposing, and simplifying, we get
(35)
[
1
4
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
t
-
1
4
σ
k
2
j
2
Δ
t
]
V
k
,
i
,
j
-
1
+
[
1
+
1
4
(
r
k
+
λ
k
)
j
Δ
t
+
1
2
σ
k
2
j
2
Δ
t
]
V
k
,
i
,
j
+
[
-
1
4
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
t
-
1
4
σ
k
2
j
2
Δ
t
]
V
k
,
i
,
j
+
1
-
∑
l
=
1
L
a
k
l
V
l
,
i
,
j
=
[
-
1
4
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
t
+
1
4
σ
k
2
j
2
Δ
t
]
V
k
,
i
+
1
,
j
-
1
+
[
1
-
1
2
σ
k
2
j
2
Δ
t
]
V
k
,
i
+
1
,
j
+
[
1
4
(
r
k
-
q
k
+
λ
k
η
k
)
j
Δ
t
+
1
4
σ
k
2
j
2
Δ
t
]
V
k
,
i
+
1
,
j
+
1
+
λ
k
Δ
t
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
j
Δ
S
)
.
Solving (35) from time (
N
-
1
)
Δ
t to time 0 recursively and taking into account the constraints of the convertible bonds, we can achieve the price of the convertible bonds at any time and any stock price. Furthermore, there are always an upper bound and a lower bound in the finite difference model. Here we set
(36)
V
k
,
i
,
N
+
1
=
κ
N
Δ
S
,
V
k
,
i
,
1
=
0
,
k
=
1,2
,
…
,
L
,
i
=
1,2
,
…
,
M
.
Early studies have shown that the Crank-Nicolson scheme is second-order convergence which refers to time. Therefore, the convertible bond price calculated by the finite difference converges to the theoretical value when the interval Δ
t tends to infinity small.


## 4.2. The Trinomial Tree Model

The trinomial tree model is another important numerical solution and is widely used in derivatives pricing. The trinomial tree method has advantage in computational efficiency to the finite difference method. In the following section we are going to apply the trinomial tree to the valuation of convertible bonds. We divide interval[
t
,
T
] into M small identical intervals and set Δ
t
=
(
T
-
t
)
/
M. When the market state is in the kth regime, the default probability is
(37)
w
k
=
1
-
e
-
λ
k
Δ
t
=
λ
k
Δ
t
+
o
(
Δ
t
)
≈
λ
k
Δ
t
,
k
=
1,2
,
…
,
L
.
Following the trinomial tree of Boyle [28], we set
(38)
u
k
=
e
α
k
σ
k
Δ
t
,
d
k
=
1
u
k
=
e
-
α
k
σ
k
Δ
t
,
m
k
=
1
,
α
k
>
0
,
k
=
1,2
,
…
,
L
.
Let π
k
,
u
,
π
k
,
m
,
π
k
,
d  (
k
=
1,2
,
…
,
L
) denote the probability when the stock price increases, remains, and decreases, respectively. Considering the expected return and volatility of stock price under the risk-neutral market and applying the risk-free arbitrage principle, these jump probabilities satisfy
(39)
π
k
,
u
+
π
k
,
m
+
π
k
,
d
=
1
,
π
k
,
u
u
k

+

π
k
,
m

+

π
k
,
d
d
k

=

e
(
r
k
+
λ
k
η
k
-
q
k
)
Δ
t
,
π
k
,
u
u
k
2

+

π
k
,
m

+

π
k
,
d
d
k
2

=

e
2
(
r
k
+
λ
k
η
k
-
q
k
)
Δ
t
+
σ
k
2
Δ
t
,
k

=

1,2
,
…
,
L
.
Solving (39), we have
(40)
π
k
,
u
=
(
U
k
-
Q
k
)
-
u
k
(
Q
k
-
1
)
(
u
k
-
d
k
)
(
1
-
d
k
)
,
π
k
,
d
=
(
U
k
-
Q
k
)
-
d
k
(
Q
k
-
1
)
(
u
k
-
d
k
)
(
u
k
-
1
)
,
π
k
,
m
=
1
-
π
k
,
u
-
π
k
,
d
,
k
=
1,2
,
…
,
L
,
where U
k
=
e
(
r
k
+
λ
k
η
k
-
q
k
)
Δ
t
,
Q
k
=
e
2
(
r
k
+
λ
k
η
k
-
q
k
)
Δ
t
+
σ
k
2. Then, we get the recurrence relation equations as
(41)
V
(
k
,
S
,
t
)
=
∑
l
=
1
L
p
k
l
e
-
r
l
Δ
t
{
π
l
,
u
[
(
1
-
w
l
)
V
(
l
,
u
l
S
,
t
+
Δ
t
)
+
w
l
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
u
l
S
)
]
+
π
l
,
u
[
(
1
-
w
l
)
V
(
l
,
m
l
S
,
t
+
Δ
t
)
+
w
l
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
m
l
S
)
]
+
π
l
,
u
[
(
1
-
w
l
)
V
(
l
,
d
l
S
,
t
+
Δ
t
)
+
w
l
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
d
l
S
)
]
}
=
∑
l
=
1
L
p
k
l
e
-
r
l
Δ
t
{
(
1
-
w
l
)
[
π
l
,
u
V
(
l
,
u
l
S
,
t
+
Δ
t
)
+
π
l
,
m
V
(
l
,
m
l
S
,
t
+
Δ
t
)
+
π
l
,
d
V
(
l
,
d
l
S
,
t
+
Δ
t
)
]
+
w
l
[
π
l
,
u
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
u
l
S
)
+
π
l
,
m
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
m
l
S
)
+
π
l
,
d
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
d
l
S
)
]
}
.
We can set the parameter α
k to decrease nodes of the trinomial tree. Thus, we set
(42)
α
k
=
σ
0
σ
k
,
k
=
1,2
,
…
,
L
,
σ
0
>
max
⁡
1
≤
k
≤
L
(
σ
k
)
.
That is,
(43)
u
k
=
e
σ
0
Δ
t
,
d
k
=
e
-
σ
0
Δ
t
,
m
k
=
1
,
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
k
=
1,2
,
…
,
L
.
Yuen and Yang [23] proposed that σ
0 should be taken as
(44)
σ
0
=
max
⁡
1
≤
k
≤
L
(
σ
k
)
+
(
1.5
-
1
)
mean
1
≤
k
≤
L
(
σ
k
)
.
Thus, the number of the nodes of the trinomial tree decreases so much that we can save many storage spaces and much computation time.


## 4.3. The Consistency of the Finite Difference and the Trinomial Tree

As we know, without regime switching, the tree model is consistent with the finite difference method, so we are going to explore their consistency with regime switching. Starting with the trinomial tree, we can achieve the relation between the trinomial tree model and the partial differential equations as the following proposition.Proposition 2.
When the intervalΔ
t tends to zero, (41) converges to (18).Proof.
For anyb
>
0, applying the multivariate Taylor’s series, we have
(45)
V
(
k
,
b
S
,
t
)
=
V
(
k
,
S
+
(
b
-
1
)
S
,
t
+
Δ
t
)
=
V
(
k
,
S
,
t
)
+
(
b
-
1
)
S
∂
V
∂
S
+
Δ
t
∂
V
∂
t
+
1
2
(
b
-
1
)
2
S
2
∂
V
∂
S
+
o
(
Δ
t
)
.
Then, we obtain
(46)
π
k
,
u
V
(
k
,
u
k
S
,
t
+
Δ
t
)
+
π
k
,
m
V
(
k
,
m
k
S
,
t
+
Δ
t
)
+
π
k
,
d
V
(
k
,
d
k
S
,
t
+
Δ
t
)
=
(
π
k
,
u
+
π
k
,
m
+
π
k
,
d
)
(
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
)
+
(
(
π
k
,
u
u
k
+
π
k
,
m
m
k
+
π
k
,
d
d
k
)
-
(
π
k
,
u
+
π
k
,
m
+
π
k
,
d
)
)
S
∂
V
∂
S
+
1
2
(
(
π
k
,
u
+
π
k
,
m
+
π
k
,
d
)
1
2
(
π
k
,
u
u
k
2
+
π
k
,
m
m
k
2
+
π
k
,
d
d
k
2
)
-
2
(
π
k
,
u
u
k
+
π
k
,
m
m
k
+
π
k
,
d
d
k
)
+
(
π
k
,
u
+
π
k
,
m
+
π
k
,
d
)
1
2
)
S
2
∂
V
∂
S
+
o
(
Δ
t
)
.
=
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
e
(
r
k
+
λ
k
η
k
)
Δ
t
-
1
)
S
∂
V
∂
S
+
1
2
(
e
2
(
r
k
+
λ
k
η
k
)
Δ
t
+
σ
k
2
Δ
t
-
2
e
(
r
k
+
λ
k
η
k
)
Δ
t
+
1
)
×
S
2
∂
V
∂
S
+
o
(
Δ
t
)
=
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
k
2
S
2
Δ
t
∂
V
∂
S
+
o
(
Δ
t
)
.
Multiplying ∏
h
=
1
L
e
r
l
Δ
t both sides of (41), we have
(47)
(
∏
h
=
1
L
e
r
l
Δ
t
)
V
(
k
,
S
,
t
)
=
∑
l
=
1
L
p
k
l
(
∏
h
≠
l
e
r
h
Δ
t
)
×
{
(
1
-
w
l
)
[
π
l
,
u
V
(
l
,
u
l
S
,
t
+
Δ
t
)
IIIIIIIIIIIIIIIIIII
+
π
l
,
m
V
(
l
,
m
l
S
,
t
+
Δ
t
)
IIIIIIIIIIIIIIIIIII
+
π
l
,
d
V
(
l
,
d
l
S
,
t
+
Δ
t
)
]
+
w
l
[
π
l
,
u
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
u
l
S
)
+
π
l
,
m
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
m
l
S
)
+
π
l
,
d
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
d
l
S
)
]
}
.
Applying Taylor’s series, we have
(48)
(
∏
h
=
1
L
(
1
+
r
h
Δ
t
+
o
(
Δ
t
)
)
)
V
(
k
,
S
,
t
)
=
(
1
+
∑
h
=
1
L
r
h
Δ
t
+
o
(
Δ
t
)
)
V
(
k
,
S
,
t
)
=
(
1
+
a
k
k
Δ
t
+
o
(
Δ
t
)
)
(
∏
h
≠
k
(
1
+
r
h
Δ
t
+
o
(
Δ
t
)
)
)
×
{
∂
V
∂
t
(
1
-
λ
k
Δ
t
+
o
(
Δ
t
)
)
×
[
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
Δ
t
∂
V
∂
S
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
+
1
2
σ
k
2
S
2
Δ
t
∂
V
∂
S
+
o
(
Δ
t
)
]
iiiiiiiii
+
(
λ
k
Δ
t
+
o
(
Δ
t
)
)
×
[
π
k
,
u
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
u
k
S
)
+
π
k
,
m
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
m
k
S
)
+
π
k
,
d
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
d
k
S
)
]
∂
V
∂
t
}
+
∑
l
≠
k
(
a
k
l
Δ
t
+
o
(
Δ
t
)
)
(
∏
h
≠
l
(
1
+
r
h
Δ
t
+
o
(
Δ
t
)
)
)
×
{
∂
V
∂
S
(
1
-
λ
l
Δ
t
+
o
(
Δ
t
)
)
×
[
∂
V
∂
S
V
(
l
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
l
+
λ
l
η
l
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
l
2
S
2
Δ
t
∂
V
∂
S
+
o
(
Δ
t
)
]
+
(
λ
l
Δ
t
+
o
(
Δ
t
)
)
×
[
π
l
,
u
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
u
l
S
)
iiiiiiiiiiiiiiiiii
+
π
l
,
m
max
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
m
l
S
)
iiiiiiiiiiiiiiiiii
+
π
l
,
d
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
d
l
S
)
]
∂
V
∂
S
}
=
(
1
+
a
k
k
Δ
t
+
o
(
Δ
t
)
)
(
1
+
∑
h
≠
k
r
h
Δ
t
+
o
(
Δ
t
)
)
×
{
∂
V
∂
t
(
1
-
λ
k
Δ
t
+
o
(
Δ
t
)
)
×
[
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
k
2
S
2
Δ
t
∂
V
∂
S
+
o
(
Δ
t
)
]
+
(
λ
k
Δ
t
+
o
(
Δ
t
)
)
×
[
π
k
,
u
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
u
k
S
)
+
π
k
,
m
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
m
k
S
)
+
π
k
,
d
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
d
k
S
)
]
∂
V
∂
t
}
+
∑
l
≠
k
(
a
k
l
Δ
t
+
o
(
Δ
t
)
)
(
1
+
∑
h
≠
l
r
h
Δ
t
+
o
(
Δ
t
)
)
×
{
∂
V
∂
S
(
1
-
λ
l
Δ
t
+
o
(
Δ
t
)
)
×
[
∂
V
∂
S
V
(
l
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
l
+
λ
l
η
l
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
l
2
S
2
Δ
t
∂
V
∂
S
+
o
(
Δ
t
)
]
+
(
λ
l
Δ
t
+
o
(
Δ
t
)
)
×
[
π
l
,
u
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
u
l
S
)
iiiiiiiiiiiii
+
π
l
,
m
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
m
l
S
)
iiiiiiiiiiiii
+
π
l
,
d
max
⁡
⁡
(
R
l
B
,
κ
(
1
-
η
l
)
d
l
S
)
]
∂
V
∂
t
}
=
V
(
k
,
S
,
t
)
+
Δ
t
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
k
2
S
2
Δ
t
∂
V
∂
S
-
λ
k
V
k
Δ
t
+
∑
h
≠
k
r
h
V
k
Δ
t
+
∑
l
=
1
L
a
k
l
V
k
Δ
t
+
o
(
Δ
t
)
+
λ
k
Δ
t
[
π
k
,
u
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
u
k
S
)
iiiiiiiiiiiiiiiiiii
+
π
k
,
m
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
m
k
S
)
iiiiiiiiiiiiiiiiiii
+
π
k
,
d
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
d
k
S
)
]
.
Eliminating the similar items of (48), we have
(49)
r
k
V
k
Δ
t
=
Δ
t
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
Δ
t
∂
V
∂
S
+
1
2
σ
k
2
S
2
Δ
t
∂
V
∂
S
-
λ
k
V
k
Δ
t
+
∑
l
=
1
L
a
k
l
V
k
Δ
t
+
o
(
Δ
t
)
+
λ
k
Δ
t
[
π
k
,
u
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
u
k
S
)
+
π
k
,
m
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
m
k
S
)
+
π
k
,
d
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
d
k
S
)
]
.
Let both sides of (49) be divided by Δ
t. Then, we have
(50)
λ
k
V
k
=
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
∂
V
∂
S
+
1
2
σ
k
2
S
2
∂
V
∂
S
-
λ
k
V
k
+
∑
l
=
1
L
a
k
l
V
k
+
o
(
1
)
+
λ
k
[
π
k
,
u
max
⁡
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
u
k
S
)
+
π
k
,
m
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
m
k
S
)
+
π
k
,
d
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
d
k
S
)
]
.
Let Δ
t
→
0; (50) converges to
(51)
∂
V
∂
t
+
(
r
k
+
λ
k
η
k
)
S
∂
V
∂
S
+
1
2
σ
k
2
S
2
∂
V
∂
S
-
(
r
k
+
λ
k
)
V
k
+
∑
l
=
1
L
a
k
l
V
k
+
λ
k
max
⁡
(
R
k
B
,
κ
(
1
-
η
k
)
S
)
=
0
.
This completes the proof.Proposition2 indicates that the trinomial tree model of the convertible bonds with credit risk under regime switching is consistent with the partial differential equations of the convertible bonds. That is, the trinomial tree model of the convertible bonds is consistent with the finite difference model of the convertible bonds.


## 5. Numerical Example and Analysis

In order to explore the effect of credit risk and regime switching on the value of the convertible bonds and show the difference of the pricing efficiency between the trinomial tree method and the finite difference method, we will introduce a numerical example in this section. The information of the convertible bonds used for the numerical example is given in Table1.Table 1
Data for numerical example.
Face value
100
Coupon payments
2% annually
Conversion ratio
10
Maturity
5 years
Conversion period
1–5 years
Put period
2–5 years
Put price
103
Call period
2–5 years
Call Price
110
Amount of share
100 million
Amount of convertible bonds
1 million
Debt before convertible bonds
400 millionFirst of all, we will compare the convergence, pricing effectiveness, and efficiency of the numerical solutions for the trinomial tree method and the finite difference method. For simplifying, we assume that there are only two states in the market. The parameters of the stock price process and the default used for the numerical example are given in Table2.Table 2
Parameters of numerical example.
Regime
Risk-free rate
Historical volatility
Continuous dividends
Hazard rate
Recovery factor
Jump ratio
Regime 1
0.06
0.20
0.02
0.01
0.75
0.50
Regime 2
0.08
0.35
0.03
0.03
0.50
0.75We assume that the current stock price is 10, the corporate income tax rate is 25%, and there is no call provision. Let the generator matrix be(52)
A
=
(
-
0.5
0.5
0.5
-
0.5
)
.
For the finite difference method, we set that the upper bound of the stock price is equal to 30 and the number of the stock price steps is equal to 300. For the trinomial tree, we set the parameters the same as Yuen and Yang [23]. Table 3 and Figure 1 present the numerical results for the current price and computation time of the convertible bonds for different time steps and different numerical solutions by MATLAB 7.0.Table 3
Comparison of the finite difference and the trinomial tree.
Time steps
200
400
800
1600
3200
6400
128000
Finite difference
Regime 1
108.7974
108.7925
108.7835
108.7754
108.7666
108.7606
108.7569
Regime 2
110.1492
110.1557
110.1538
110.1500
110.1437
110.1392
110.1363
Computation time (s)
119.0600
253.1584
504.5852
1113.8
2275.4
4983.7
10823
Trinomial tree
Regime 1
108.9136
108.7876
108.7509
1087970
108.7581
108.7841
108.7598
Regime 2
110.2378
110.1481
110.1272
110.1639
110.1367
110.1578
110.1389
Computation time (s)
0.8580
3.1980
11.4817
46.3323
192.5520
763.3285
2932.1Comparison of the finite difference and the trinomial tree on convergence.
(a)
(b)From Table3 and Figure 1, the finite difference method is monotonically converging to the theoretical value; the trinomial tree method is also convergent but not monotonous. Although the numerical results of the trinomial tree and the finite difference are different, the difference between two numerical solutions is so small that two numerical solutions are consistent, which is the same as conclusions in Section 4.3. In terms of convergence, we need more time steps to make the numerical solution converge. However, the results tend to be stable when the number of time steps tends to infinity and the value changes a little when the number of time steps is over 2000. Hence, the numerical results are acceptable when the number of time steps is over 2000. In terms of the computation time, the trinomial tree needs less time than the finite difference when the number of time steps is the same and not large. However, the trinomial tree may need more time than the finite difference when the number of time steps is large enough. What is more, the numerical result of the finite difference is more stable than that of the trinomial tree. But we should also notice that the finite difference method also achieves the convertible bond price at any time and any stock price. So it is convenient to use the finite difference to study the effect of time and stock price on the convertible bond and calculate the optimal conversion price. We can choose one numerical method as we need.Next, we will study the relation of stock, time, and convertible bond price. Following Table1, Table 2, and generator matrix (50), the numerical result of convertible bond pricing at different stock prices and time by the finite difference method is graphed as Figure 2.Convertible bonds value for different times and stock prices.
(a)
(b)Figure2 shows that the trends of different regimes are the same but only different in value. Once time is given, the convertible bond price increases with the stock price increasing for the conversion value of the convertible bond. That is, the higher the stock price is, the higher the convertible bond price is. Similarly, once the stock price is given, the convertibles bond price changes wavily with the time increasing, which mainly is due to the discrete coupons payment. That is, between two coupon payments, the convertible bond price increases with the time. The reason is that the present value of the coupon increases when the time of coupon payment is near. And when the coupon is paid, the value of a convertible bond will jump with the same amount of coupon. But the height of the waves decreases with time. These indicate that the nearer the time is, the lower the convertible bond price is, which is the same as that without regime switching. The reason is that the value of option part becomes smaller and more coupons are paid when the maturity is coming.At last, we will study the effect of regime switching on the convertible bond price. Applying the trinomial tree method, Table1, and Table 2, let the number of time steps be 1000. We assume that there are only two states in the market. Set
(53)
A
=
(
-
a
1
a
1
a
2
-
a
2
)
.
By changing the parameters a
1 and a
2 from 0 to 10 continuously, we obtain the convertible bond prices at different generator matrixes and graph the results as Figure 3.Convertible bonds value for different generator matrix.
(a)
(b)Figure3 shows that different generator matrixes correspond to different convertible bond prices, though another condition is the same. The larger the parameters a
1 and a
2 are, the stronger the regime switching is. The regime switching strength presents the probability and frequency of state transition. So the stronger the regime switching is, the closer the convertible bond prices at different states are. In fact, when the parameters a
1 and a
2 tend to infinity together, the convertible bond prices tend to be consistent. We also find that the strength of the effect of regime switching depends on the generator matrix or the regime switching strength. Based on the theory of Markov chain, the larger the element of the generator matrix is, the more frequent one regime changes to the other. With the effect of other regimes, the value of the convertible bonds will be affected by that of other regimes. It is easy to know from this that the bigger the element of generator matrix is, the closer the value of the two regimes is. For example, the bigger a
1is, the closer the value of regime 2 is to that of regime 1. Here we assume there are only two states in the market, but these conclusions are also applicable for more than two states.In addition, we can use the finite difference scheme to calculate the optimal conversion price of the convertible bond at any time. If the stock price is more than the optimal conversion price, the value of the convertible bonds is not more than their conversion value, so that the investors will always choose to convert the convertible bonds. We can calculate the optimal conversion price by seeking the least stock price in which the value of the convertible bonds is equal to their conversion price. We take the numerical example with Table1, Table 2, and (50). The optimal conversion prices of callable and noncallable convertible bonds are presented in Figure 4.The optimal conversion prices of shares for callable and noncallable convertible bonds.
(a)
(b)Figure4 indicates that the optimal conversion prices of shares at different regimes are always different except on expiry date or callable date. The reason is that the convertible bond prices at different regimes with the same conditions are always different except on expiry date or during call period, but the conversion values of convertible bonds are always the same. The regime switching has effect not only on the value of convertible bonds but also on the conversion behavior. Compared with the callable convertible bonds, the optimal conversion value of shares for noncallable convertible bonds is much lower. Within the call period, the optimal conversion price of shares is equal to the call price. If the convertible bond has a call clause, the investors may choose to convert the convertible bonds early. The unlimited call clause generally limits the value of the convertible bonds and the optimal conversion price of shares, so the issuers will set up the call clause and the put clause with some requirements or even cancel them.


## 6. Conclusions

The financial market exhibits different states in different stages, and different states interchange, such as the interconversion of “bull market” and “bear market” in the stock market. This paper has discussed the convertible bond pricing with credit risk under regime switching. We have obtained the pricing model of the convertible bonds by deriving the Black-Scholes-type partial differential equation of the convertible bonds. Considering the dilution effect and the debt leverage, we have given an operational modification on the pricing model, which offers a new idea in the convertible bond pricing. Our proposed two numerical solutions are convergent and consistent. The regime switching has effect on the value of the convertible bond; its strength depends on the generator matrix. Using this model, we can also find the optimal conversion strategy.The proposed pricing model is one-factor model. Multifactor models including two-factor model can be studied in the future. In addition, there exists not only regime switching in the market but also other market features, such as long memory, jump, and fuzziness. Usually, these market features are mixed, not separated. Pricing convertible bonds under mixed environments will be an interesting topic in the future.



---

*Source: 381943-2014-05-19.xml*