# Pricing of Embedded Options in Bank Deposits and Loans Based on Jump-Diffusion Interest Rate Model

Authors: Enlin Tang; Song Xu
Journal: Complexity
 (2021)

---

## Abstract

The marketization of interest rate is an inevitable requirement for China’s financial reform and joining the WTO to connect with the international financial market. It is also an important link to improve the marketization degree of China’s financial system. The marketization of interest rate in China is gradually advancing according to its preset mode. In the process of interest rate marketization, an unavoidable problem is that while the interest rate marketization gives the commercial banks the autonomy of capital pricing, the fluctuation of interest rate is more and more frequent. However, due to the fluctuation of interest rate, the loan as the main assets of commercial banks will be prepayed by borrowers, and the time deposit as the main liabilities of commercial banks will be withdrawn by depositors in advance; that is, embedded options are implied in asset liability items, which makes it difficult for commercial banks to accurately calculate the actual interest margin of deposits and loans and manage the interest rate risk. Therefore, it is of great significance to identify and price such embedded option value. On the basis of identifying and decomposing the embedded options in deposit and loan of commercial banks, according to the change characteristics of deposit and loan interest rate of Chinese commercial banks, this paper chooses jump-diffusion interest rate model to describe the change of benchmark interest rate of deposit and loan in China and demonstrates the advantages of this model compared with other models. Based on Monte Carlo simulation technology, the embedded options of five-year fixed deposit and ten-year prepayable loan in China are priced. On this basis, it points out that the real interest margin of commercial bank’s deposit and loan should be the nominal interest margin minus the value of deposit and loan’s embedded options. In the process of interest rate risk management, we should pay attention to the existence of embedded options and carry out effective management.


---

## Body

## 1. Introduction

Embedded options refer to the options embedded in bank assets and liabilities. When the interest rate fluctuates greatly, the customers of the bank will exercise the option implied in the assets or liabilities for the purpose of self-interest (such as prepaying the loan or withdrawing the time deposit in advance), thus bringing uncertainty to the bank’s income [1]. This uncertainty obviously constitutes a source of risk, which is called “embedded option risk.” With the advancement of interest rate marketization in China, financial derivatives with interest rate as the subject matter are increasing. At the same time, assets and liabilities with embedded option are increasingly appearing in the balance sheets of banks and other financial institutions [2, 3]. However, in China’s current financial products, including time deposits, mortgage loans, and mortgage loans, there is no pricing for the embedded options implied in deposits and loans, which leads to the pricing of deposits and loans too low. With the advancement of interest rate marketization in China, the existence of embedded options has become one of the important reasons for the interest rate risk of banks and other financial institutions. Therefore, it is of great significance for China’s commercial banks to study the pricing technology of embedded option implied in the assets and liabilities of banks and other financial institutions in terms of product development and risk management. Through the collection of option fee, banks can reduce their operating costs and refinancing risks and increase their liquidity. It is of great theoretical and practical significance to measure and control the prepayment risk of assets and the withdrawal risk of liabilities and to establish a practical and effective asset liability portfolio optimization model, which can provide decision support for the risk control of commercial banks and the optimal allocation of asset structure.


## 2. Literature Review

The sharp fluctuation of interest rate and the sharp increase in interest rate risk in the international financial market have attracted many scholars’ in-depth research on the measurement of interest rate risk and the corresponding management and control technology. Foreign scholars’ research on embedded options and interest rate risk was started earlier and made outstanding achievements. On the pricing of embedded options, Easton [4] studied the principle and method of pricing embedded options in Asian options by Monte Carlo simulation and obtained the bond price with Asian options through parameter estimation. Gilkson and Ruff [5] studied the early withdrawal problem of small retail CDs and proposed a pricing method of CDs with early withdrawal right. Glikson et al. [6] analyzed the correlation between the interest rate level and the embedded option value in CDs by using cross-sectional data. The research shows that there is a general phenomenon of early withdrawal of long-term time deposits, and the early withdrawal behavior of depositors increases with the increase in interest rate. Brooks [7] calculated the annual percentage yield (APY) and effective annual rate (EAR) of CDs to analyze the actual yield level of products. Tolstykh et al. [8] also found that the risk of embedded options may come from the change of asset and liability convexity, which will lead to unexpected changes in asset and liability prices. Turner et al. [9] studied the embedded options in deposits and analyzed the characteristics of options that are beneficial to depositors in time deposits and the time characteristics of depositors’ exercise of options. Cline and Brooks [10] analyzed the characteristics of some exotic CDs in recent years and their pricing [11].Domestic research on embedded options is more about the risk management of embedded options or the measurement of risk. For example, Wang and Zhang [12] constructed the goal programming model of interest rate risk management and thought that compared with the duration gap model, the convexity gap model had better effect in reducing interest rate risk. Luo and Wan [13] studied the interest rate risk control of bank’s assets and liabilities when the embedded option exists and proposed the control strategy of interest rate risk. Zheng and Lin [14] proposed a new standard to measure the impact of embedded options on bond interest rate risk and used the binary tree model to make an empirical analysis on the callable and resolvable financial bonds issued by China Development Bank. Ren and Li [15] selected the Vasicek interest rate model and used the method of partial differential equation to approximate the calculation formula of time deposit value. Yang [16] focused on the effective duration and effective convexity model based on embedded option and option-adjusted spread and studied its specific application in the interest rate risk management of commercial banks in China. Dai [17] used no arbitrage pricing method and OAS method to calculate the embedded options based on the decomposition of the embedded options. Yi and Liu [18] used the Monte Carlo simulation method and trigeminal tree method to reconstruct the interest rate trend, and based on this, they studied the interest rate risk of embedded options. However, there is little research on the pricing of embedded options. Zheng and Lin [19] used no arbitrage analysis method and numerical calculation method to price the embedded options and made the sensitivity analysis of the option price to various factors. Liu et al. [20] established the mathematical model of prepayable loan under the CIR interest rate model and discussed the influence of stochastic interest rate on prepayment. Zhu et al. [21] divided the embedded options of deposits into two parts such as interest rate risk and liquidity risk and used the numerical analysis method to price the fair value of the two parts. Zhu et al. [22] used the trigeminal tree analysis method to price the embedded options. Xia et al., [23] priced the deposit and loan with embedded options, and on this basis, they analyzed the impact of embedded options on the interest rate risk management of commercial banks.Through literature review, we find that in recent years, although Chinese scholars have gradually begun to study the interest rate risk of financial instruments with embedded options, the above research mainly focuses on the measurement of interest rate risk with embedded options, as well as the analysis of the risk caused by the prepayment of loans or the analysis of the characteristics of embedded options. There is a lack of research on the pricing of embedded options in deposits and loans, especially the advance execution of embedded options in deposits and loans based on the dynamic change of interest rate. As an inducement of interest rate risk, it is an inevitable choice for Chinese banks to price the embedded option scientifically in the face of external shocks after the interest rate liberalization. Generally speaking, the research on interest rate risk of financial institutions in China started relatively late, and the research on pricing model and practical application of interest rate derivatives such as embedded options is rare. This situation needs to be changed in the innovation and development of financial instruments today. In addition, in terms of research methods, the arbitrage free analysis method and differential equation method proposed by predecessors are all based on a series of harsh assumptions, and in more complex cases (such as adding the jump factor of interest rate), they may not be able to get accurate solutions. Based on its own advantages, the Monte Carlo simulation method can solve multidimensional securities pricing, which will become the development trend of this kind of financial product pricing.


## 3. Recognition of Embedded Options in Bank Deposits and Loans

### 3.1. Decomposition of Embedded Options in Bank Deposits and Loans

In the asset business of banks, loans account for the largest proportion of assets. Loans can be divided into floating rate loans and fixed rate loans [24]. For the fixed rate loan, the borrower must pay a fixed amount of interest to the bank according to the fixed interest rate under normal circumstances. Therefore, the fixed rate loan can be regarded as a redeemable bond issued by the borrower to the bank. The face value of the redeemable bond is equal to the face value of the loan, and the execution price is the loan interest rate. For banks, it is equivalent to the buyer of callable bonds, and it is short of American put options; the borrower is equivalent to the issuer of callable bonds and is also the bull of American put options. As shown in Figure 1, this is equivalent to the bank issuing an American put option to the borrower and the borrower issuing a bond to the bank.Figure 1
The decomposition of bank loan with embedded option.When banks make loan pricing, they often only pay attention to the content of ② in Figure1 and ignore the content of ① in Figure 1. In theory, the borrower should pay a certain option fee to the bank when he obtains the American option. Otherwise, banks should consider the existence of options when they determine the interest rate. However, for the consideration of laws and regulations or industry competition, banks ignore the existence of options. When the interest rate drops, the rational borrower will choose to repay in advance, and then apply to the bank for a loan with lower interest rate. The value of options will determine the motivation of borrowers to prepay. When the interest rate drops sharply and the penalty is very small, the value of options will increase, and the emergence of prepayment behavior will make banks face losses.Fixed rate loan is when the interest rate drops, the borrower prepays in order to reduce the interest expense. However, for floating rate loans, prepayment occurs when the interest rate rises. In a channel of interest rate rise, the loan interest rate will be priced according to the higher interest rate of the next period, so the borrower will choose prepayment when the funds are abundant. Whether it is a floating rate loan or a fixed rate loan, prepayment will reduce the interest income of the bank and increase the risk faced by the bank.The same is true of bank’s liability business, which mainly includes all kinds of deposit business. According to whether it can be withdrawn at any time, bank deposits can be divided into time deposits and current deposits. The current deposit is mainly for the borrower to meet the daily liquidity, which can be withdrawn at any time. The interest rate of current deposit is generally low, and its duration is short. The nature of embedded option in current deposit is not obvious, and its sensitivity to interest rate is low. Time deposit is different, and its duration and interest rate sensitivity are greater and vulnerable to changes in interest rates. With the method of financial engineering, time deposit can be decomposed into the combination of securities and options without options. If time deposit is compared with a kind of securities, then the bank can be regarded as the issuer of securities (time deposit) and options, and the depositor is equivalent to the purchaser of securities (time deposit) and options, as shown in Figure2. The value of options will determine the motivation of depositors to withdraw in advance.Figure 2
The decomposition of bank deposits with embedded option.From the above analysis, we can see that the main assets and liabilities of commercial banks have the nature of embedded options. Bank deposits and loans with embedded options can be regarded as the combination of deposits loans without embedded options and options. As shown in Table1, the value of deposits with embedded option is equal to the value of deposits without embedded option plus the value of call options; the value of loans with embedded option is equal to the value of loans without embedded option minus the value of American call options.Table 1
Decomposition of bank deposits and loans with embedded options.
Banking itemsCorresponding bond positionCorresponding option positionBank depositsCallable bonds issued by banks (short) = ordinary bonds + resale rightInterest rate call options (short)Bank loansRedeemable bonds issued by borrowers (long) = ordinary bonds – redemption rightInterest rate put options (short)


### 3.2. Analysis of the Exercise Boundary and Characteristics of Embedded Options in Deposits and Loans

#### 3.2.1. Analysis of the Implementation and Related Characteristics of Embedded Options in Deposits

Theorem 1.
For depositors, time deposits can be decomposed into a fixed debt and several American interest rate call option bulls; for banks, time deposits can be decomposed into a fixed debt and several American interest rate call option bulls.
Suppose at 0, a depositor deposits a time deposit with a maturity ofTin the bank, the continuous compound interest rate is r0,T and the deposit amount is standardized to 1. If the withdrawal in advance does not occur, the final income of the deposit is er0,TT. Assuming that the investor’s investment period is from 0 to T, the depositor’s goal is to maximize the return from 0 to T.
Suppose at timet, the central bank suddenly announced to change the level of deposit interest rate, the changed current interest rate is rt,t and the fixed deposit interest rate from T to t period is rt,T. At this point,
Ifer0,TT≥ert,tt+rt,TT−t, the income from maintaining the original time deposit exceeds the transferred deposit, and the depositor chooses to maintain the original time deposit.
Ifer0,TT≤ert,tt+rt,TT−t, the income of withdrawal and deposit in advance exceeds that of maintaining the original time deposit, and the depositor chooses to withdraw and deposit in advance, and the amount of deposit in advance is ert,tt.
To sum up, at the moment, the depositor’s final income function can be expressed as follows:(1)maxer0,TT,ert,tt+rt,TT−t=er0,TT+maxert,tt+rt,TT−t−er0,TT,0=er0,TT+maxer0,TTert,tt+rt,TT−t−r0,TT−1,0≈er0,TT+er0,TTmaxrt,tt+rt,TT−t−r0,TT,0.
In formula (1), er0,TT is a basic bond contract, and er0,TTmaxrt,tt+rt,TT−t−r0,TT,0 is an interest rate call option with an exercise price of K=r0,TT [25, 26]. The conditions for the exercise of the option are as follows:(2)rt,tt+rt,T≥r0,TT.
By transforming formula (2), we get the following results:(3)rt,tt+rt,TT≥r0,T.
In addition, after the option is exercised, the depositor obtains a fixed deposit with a term ofT−t and naturally obtains an option, that is to say, during the T−t period, the depositor can withdraw and deposit in advance, and the exercise price of the option is rt,T. Therefore, for depositors, a time deposit can be decomposed into a basic debt and several American options. However, as time goes on, the possibility of the option being executed gradually decreases, and its value gradually decreases. The value of embedded options in time deposits is reflected in a period of time shortly after the contract is signed.Corollary 1.
The necessary condition for option execution is the rise of interest rate.Proof.
The conditions of option execution are as follows:(4)rt,tt+rt,T≥r0,TT,rt,T≥r0,TT−rt,ttT−t.
Because under normal conditionsrt,t≤r0,T,(5)rt,T≥r0,TT−rt,ttT−t≥r0,TT−r0,TtT−t=r0,T.
In other words, the depositor’s immediate income is negative, which needs to be compensated by the future excess income.Corollary 2.
The more the interest rate rises, the more likely it is to withdraw in advance.Proof.
The more the interest rate rises, the less the immediate loss of early withdrawal and the more the excess income in the future, so the more likely the future income is to make up for the immediate loss, the more likely it is to withdraw in advance.Corollary 3.
The larger thet, the less likely it is to withdraw in advance.Proof.
The condition of option execution is as follows:(6)rt,T≥r0,TT−rt,ttT−t≥r0,T−rt,ttT−t+r0,T.
Therefore, the minimum executive interest rate isrt=r0,T−rt,tt/T−t+r0,T.
Because∂rt/∂t=r0,T−rt,tT/T−t2≥0, the minimum exercise interest rate level increases with the extension of time, and the possibility of option execution also decreases.


#### 3.2.2. Analysis on the Execution and Related Characteristics of Embedded Options in Loans

Theorem 2.
For the borrower, the loan can be decomposed into a basic debt and several American interest rate put options; for the lending bank, the loan can be divided into a basic creditor’s right and several short positions of American interest rate put options.
Suppose at time 0, an investor signs a loan contract with a bank with a term ofT, the interest rate of the continuous compound interest loan is r0,T, the loan amount is standardized as 1, and the penalty is α. If there is no prepayment, the total cost he has to pay at the end of the period is er0,TT. Assuming that the borrower’s borrowing period is constant from 0 to T, the goal is to minimize the cost of the period. The credit rating of the borrower has not changed during the term of the loan.
Suppose at timet, the central bank suddenly announced to change the interest rate level, and the interest rate of the loan with a term T−t is r1,T, then if er0,TT≥er0,Tt+αer1TT−t, the cost of maintaining the loan exceeds prepayment, payment of liquidated damages, and signing a new contract, and the borrower will choose prepayment, payment of liquidated damages, and signing a new contract, and the loan amount of the new contract is er0,Tt+α; if er0,TT≤er0,Tt+αer1TT−t, the cost of maintaining the original term loan is lower than prepayment, paying liquidated damages, and signing a new contract, and the borrower maintains the original loan contract.
To sum up, at timet, the borrower’s final cost function can be expressed as follows:(7)−miner0,Tt+αer1TT−t,er0,TT=−er0,TT−miner0,Tt+αer1TT−t−er0,TT,0=−er0,TT+maxer0,TT−er0,Tt+αer1TT−t,0.
In formula (7), −er0,TTrepresents a basic debt contract with credit risk.
Andmaxer0,TT−er0,Tt+αer1TT−t,0 represents an interest rate put option with an exercise price of K=er0,TT. The conditions of option execution are as follows:(8)er0,TT≥er0,Tt+αer1TT−t,r1,T≤r0,TT−lner0,Tt+αT−t.Corollary 4.
The necessary condition of option execution is the reduction of interest rate.Proof.
The conditions of option execution are as follows:(9)er0,TT≥er0,Tt+αer1TT−t.
If the interest rate increases, thenrt,T≥r0,T, and then er0,TT≤er0,T+rt,TT−t<er0,Tt+αer1TT−t, and the option will not be executed, expressed mathematically as follows:(10)r1,T≤r0,TT−lner0,Tt+αT−t<r0,TT−r0,TtT−t=r0,T.
That is to say, for the borrower, the execution of options leads to a fine, which needs to be made up by the reduction of future interest.Corollary 5.
The larger theα, the less likely the option will be executed.Proof.
The highest interest rate of option execution is(11)r=r0,TT−lner0,Tt+αT−t,∂r∂α=−1er0,Tt+αT−t<0.
Therefore, the larger theα is, the lower the interest rate is needed before the option can be executed.Corollary 6.
The larger thet is, the less likely the prepayment is.Proof.
The highest interest rate of option execution is(12)r=r0,TT−lner0,Tt+αT−t,∂r∂t=−lner0,Tt+αT−t2<0.
After the option is executed, the borrower can reapply for a loan with a term ofT−t and naturally get another option; that is, it can prepay, pay liquidated damages, and sign a new contract during the period of T−t. The execution price is er1TT−t. Therefore, for borrowers, a term loan can be decomposed into a basic debt and several American options. With the delay of time, the possibility of the option being executed is gradually reduced, and its value is gradually reduced. The value of embedded options in loans is reflected in a period of time shortly after the contract is signed.


### 3.3. The Mechanism of the Effect of Embedded Options on the Value of Bank Deposits and Loans

Loan is the main part of commercial bank’s assets, and time deposit is the main part of commercial bank’s liabilities. Therefore, when the loan of commercial bank can be paid in advance and the time deposit can be withdrawn in advance, that is, embedded option is implied in assets and liabilities, the value of commercial bank’s assets and liabilities will be affected by the existence of embedded option. Especially when the market interest rate changes, the duration and convexity of assets and liabilities will change their original characteristics because of the existence of embedded options. This paper uses the graph to illustrate [27].For asset items whose main part is loans, especially fixed rate loans, the rise of interest rate has no effect on them. However, when the interest rate falls belowR∗, the borrowers execute the embedded options and choose to prepay, and the absolute value of asset duration becomes smaller. When the interest rate falls to a certain extent, a large area of borrowers will prepay. This makes the duration of assets further smaller and the convexity becomes negative, which is reflected in the value-interest rate curve of assets in Figure 3; that is, the curve just begins to become flat until it bends in the opposite direction, which means that the value of assets just begins to rise with the decline of interest rate, but the further decline of interest rate will lead to the borrower’s exercise of rights and the bottleneck of the rise of asset value, with the upper limit of P∗; for the liability items whose main part is time deposit, the decrease in interest rate has no effect on them. However, when the interest rate rises above R∗, the depositor executes the embedded options and chooses to withdraw in advance, which leads to the decrease in the absolute value of the duration of the liability, which is reflected in the value-interest rate curve of the liability in Figure 4; that is, the curve becomes flat. This means that the value of the liability has a lower limit as interest rates rise.Figure 3
Value-interest rate curve of assets.Figure 4
Value-interest rate curve of liabilities.In order to more clearly explain the impact of the existence of embedded options on the net assets (the difference between assets and liabilities) of commercial banks, this paper uses the mathematical expression of net assets and the value-interest rate curve of assets and liabilities to explain in more detail. Usually, in order to prevent interest rate risk exposure, commercial banks match the duration convexity of asset items and liability items to make the slope of value-interest rate curve of asset and liability the same. Here, bank assets are represented byA, liabilities by L, and net assets by E. Then, dE=dA−dL. The derivative of this method is as follows:(13)dA=dAdRAdRA+12d2AdRA2dRA2,dL=dLdRLdRL+12d2LdRL2dRL2.LetA and L have the same initial interest rate, RA=RL=R, and dRA=dRL=dR. We can get DA=−dA/AdR, DL=−dL/LdR, CA=−d2A/AdR2, and CL=−d2L/LdR2, and then the dE expression is as follows:(14)dE=−DA−KDLAdR+12CA−KCLAdR2.In equation (14), DA is the duration of assets, DL is the duration of liabilities, CA is the convexity of assets, and CL is the convexity of liabilities. K=L/A is the ratio of assets to liabilities, also known as leverage. DGAP=DA−KDL can reflect the duration matching of assets and liabilities, which is called duration gap; CGAP=CA−KCL can reflect the convexity matching of assets and liabilities, which is called convexity gap. From the expression of formula (14), we can know that the main factors affecting the net assets of commercial banks are A, DGAP=DA−KDL, CGAP=CA−KCL, and R [28–30].Next, we use the value-interest rate curve of bank assets and liabilities in Figure5 to illustrate the impact on the bank’s net assets when the asset liability items are implied with embedded options. Assume that the initial net assets of the bank are greater than zero, so the value-interest rate curve of the assets is above the value-interest rate curve of the liabilities. At the beginning, without considering the existence of embedded options, commercial banks have matched the duration and convexity of assets and liabilities. When the interest rate changes little, that is, between R∗ and R∗∗, the effect of asset liability management is very good, the net assets remain stable, and there is no interest rate risk exposure. However, when the interest rate falls below R∗∗, the depositors will not be affected, the borrowers will repay in advance, and the net assets will begin to decline, especially when the interest rate falls below R1, the net asset is E1, significantly less than the initial net assets E∗; similarly, when the interest rate rises to R∗, the lender will not be affected, and the depositor withdraws the deposit in advance, especially when the interest rate rises to R2, the net assets of the bank become E2, significantly less than the initial net assets E∗ [31].Figure 5
Value-interest rate curve of bank assets and liabilities.If we analyze formula (14), when the interest rate falls further below R∗∗, the duration gap DGAP=DA−KDL changes from 0 to negative and the convexity gap DGAP=DA−KDL also changes from 0 to negative so that the net assets become smaller; when the interest rate rises above R∗, the duration gap DGAP=DA−KDL changes from 0 to positive and the convexity gap DGAP=DA−KDL changes from 0 to negative, so the net assets become smaller. That is to say, as long as the bank’s loan can be repaid in advance and the time deposit can be withdrawn in advance, that is, when the bank’s asset liability project is implied with embedded options, the net assets will be less no matter whether the interest rate rises or falls. However, commercial banks ignore this kind of embedded option with intrinsic value, which leads to that the real interest rate difference between deposit and loan is not the difference between loan interest rate and deposit interest rate on the surface, which leads to the hidden interest loss of commercial banks, and the difficulty of interest rate risk management will increase.


## 3.1. Decomposition of Embedded Options in Bank Deposits and Loans

In the asset business of banks, loans account for the largest proportion of assets. Loans can be divided into floating rate loans and fixed rate loans [24]. For the fixed rate loan, the borrower must pay a fixed amount of interest to the bank according to the fixed interest rate under normal circumstances. Therefore, the fixed rate loan can be regarded as a redeemable bond issued by the borrower to the bank. The face value of the redeemable bond is equal to the face value of the loan, and the execution price is the loan interest rate. For banks, it is equivalent to the buyer of callable bonds, and it is short of American put options; the borrower is equivalent to the issuer of callable bonds and is also the bull of American put options. As shown in Figure 1, this is equivalent to the bank issuing an American put option to the borrower and the borrower issuing a bond to the bank.Figure 1
The decomposition of bank loan with embedded option.When banks make loan pricing, they often only pay attention to the content of ② in Figure1 and ignore the content of ① in Figure 1. In theory, the borrower should pay a certain option fee to the bank when he obtains the American option. Otherwise, banks should consider the existence of options when they determine the interest rate. However, for the consideration of laws and regulations or industry competition, banks ignore the existence of options. When the interest rate drops, the rational borrower will choose to repay in advance, and then apply to the bank for a loan with lower interest rate. The value of options will determine the motivation of borrowers to prepay. When the interest rate drops sharply and the penalty is very small, the value of options will increase, and the emergence of prepayment behavior will make banks face losses.Fixed rate loan is when the interest rate drops, the borrower prepays in order to reduce the interest expense. However, for floating rate loans, prepayment occurs when the interest rate rises. In a channel of interest rate rise, the loan interest rate will be priced according to the higher interest rate of the next period, so the borrower will choose prepayment when the funds are abundant. Whether it is a floating rate loan or a fixed rate loan, prepayment will reduce the interest income of the bank and increase the risk faced by the bank.The same is true of bank’s liability business, which mainly includes all kinds of deposit business. According to whether it can be withdrawn at any time, bank deposits can be divided into time deposits and current deposits. The current deposit is mainly for the borrower to meet the daily liquidity, which can be withdrawn at any time. The interest rate of current deposit is generally low, and its duration is short. The nature of embedded option in current deposit is not obvious, and its sensitivity to interest rate is low. Time deposit is different, and its duration and interest rate sensitivity are greater and vulnerable to changes in interest rates. With the method of financial engineering, time deposit can be decomposed into the combination of securities and options without options. If time deposit is compared with a kind of securities, then the bank can be regarded as the issuer of securities (time deposit) and options, and the depositor is equivalent to the purchaser of securities (time deposit) and options, as shown in Figure2. The value of options will determine the motivation of depositors to withdraw in advance.Figure 2
The decomposition of bank deposits with embedded option.From the above analysis, we can see that the main assets and liabilities of commercial banks have the nature of embedded options. Bank deposits and loans with embedded options can be regarded as the combination of deposits loans without embedded options and options. As shown in Table1, the value of deposits with embedded option is equal to the value of deposits without embedded option plus the value of call options; the value of loans with embedded option is equal to the value of loans without embedded option minus the value of American call options.Table 1
Decomposition of bank deposits and loans with embedded options.
Banking itemsCorresponding bond positionCorresponding option positionBank depositsCallable bonds issued by banks (short) = ordinary bonds + resale rightInterest rate call options (short)Bank loansRedeemable bonds issued by borrowers (long) = ordinary bonds – redemption rightInterest rate put options (short)


## 3.2. Analysis of the Exercise Boundary and Characteristics of Embedded Options in Deposits and Loans

### 3.2.1. Analysis of the Implementation and Related Characteristics of Embedded Options in Deposits

Theorem 1.
For depositors, time deposits can be decomposed into a fixed debt and several American interest rate call option bulls; for banks, time deposits can be decomposed into a fixed debt and several American interest rate call option bulls.
Suppose at 0, a depositor deposits a time deposit with a maturity ofTin the bank, the continuous compound interest rate is r0,T and the deposit amount is standardized to 1. If the withdrawal in advance does not occur, the final income of the deposit is er0,TT. Assuming that the investor’s investment period is from 0 to T, the depositor’s goal is to maximize the return from 0 to T.
Suppose at timet, the central bank suddenly announced to change the level of deposit interest rate, the changed current interest rate is rt,t and the fixed deposit interest rate from T to t period is rt,T. At this point,
Ifer0,TT≥ert,tt+rt,TT−t, the income from maintaining the original time deposit exceeds the transferred deposit, and the depositor chooses to maintain the original time deposit.
Ifer0,TT≤ert,tt+rt,TT−t, the income of withdrawal and deposit in advance exceeds that of maintaining the original time deposit, and the depositor chooses to withdraw and deposit in advance, and the amount of deposit in advance is ert,tt.
To sum up, at the moment, the depositor’s final income function can be expressed as follows:(1)maxer0,TT,ert,tt+rt,TT−t=er0,TT+maxert,tt+rt,TT−t−er0,TT,0=er0,TT+maxer0,TTert,tt+rt,TT−t−r0,TT−1,0≈er0,TT+er0,TTmaxrt,tt+rt,TT−t−r0,TT,0.
In formula (1), er0,TT is a basic bond contract, and er0,TTmaxrt,tt+rt,TT−t−r0,TT,0 is an interest rate call option with an exercise price of K=r0,TT [25, 26]. The conditions for the exercise of the option are as follows:(2)rt,tt+rt,T≥r0,TT.
By transforming formula (2), we get the following results:(3)rt,tt+rt,TT≥r0,T.
In addition, after the option is exercised, the depositor obtains a fixed deposit with a term ofT−t and naturally obtains an option, that is to say, during the T−t period, the depositor can withdraw and deposit in advance, and the exercise price of the option is rt,T. Therefore, for depositors, a time deposit can be decomposed into a basic debt and several American options. However, as time goes on, the possibility of the option being executed gradually decreases, and its value gradually decreases. The value of embedded options in time deposits is reflected in a period of time shortly after the contract is signed.Corollary 1.
The necessary condition for option execution is the rise of interest rate.Proof.
The conditions of option execution are as follows:(4)rt,tt+rt,T≥r0,TT,rt,T≥r0,TT−rt,ttT−t.
Because under normal conditionsrt,t≤r0,T,(5)rt,T≥r0,TT−rt,ttT−t≥r0,TT−r0,TtT−t=r0,T.
In other words, the depositor’s immediate income is negative, which needs to be compensated by the future excess income.Corollary 2.
The more the interest rate rises, the more likely it is to withdraw in advance.Proof.
The more the interest rate rises, the less the immediate loss of early withdrawal and the more the excess income in the future, so the more likely the future income is to make up for the immediate loss, the more likely it is to withdraw in advance.Corollary 3.
The larger thet, the less likely it is to withdraw in advance.Proof.
The condition of option execution is as follows:(6)rt,T≥r0,TT−rt,ttT−t≥r0,T−rt,ttT−t+r0,T.
Therefore, the minimum executive interest rate isrt=r0,T−rt,tt/T−t+r0,T.
Because∂rt/∂t=r0,T−rt,tT/T−t2≥0, the minimum exercise interest rate level increases with the extension of time, and the possibility of option execution also decreases.


### 3.2.2. Analysis on the Execution and Related Characteristics of Embedded Options in Loans

Theorem 2.
For the borrower, the loan can be decomposed into a basic debt and several American interest rate put options; for the lending bank, the loan can be divided into a basic creditor’s right and several short positions of American interest rate put options.
Suppose at time 0, an investor signs a loan contract with a bank with a term ofT, the interest rate of the continuous compound interest loan is r0,T, the loan amount is standardized as 1, and the penalty is α. If there is no prepayment, the total cost he has to pay at the end of the period is er0,TT. Assuming that the borrower’s borrowing period is constant from 0 to T, the goal is to minimize the cost of the period. The credit rating of the borrower has not changed during the term of the loan.
Suppose at timet, the central bank suddenly announced to change the interest rate level, and the interest rate of the loan with a term T−t is r1,T, then if er0,TT≥er0,Tt+αer1TT−t, the cost of maintaining the loan exceeds prepayment, payment of liquidated damages, and signing a new contract, and the borrower will choose prepayment, payment of liquidated damages, and signing a new contract, and the loan amount of the new contract is er0,Tt+α; if er0,TT≤er0,Tt+αer1TT−t, the cost of maintaining the original term loan is lower than prepayment, paying liquidated damages, and signing a new contract, and the borrower maintains the original loan contract.
To sum up, at timet, the borrower’s final cost function can be expressed as follows:(7)−miner0,Tt+αer1TT−t,er0,TT=−er0,TT−miner0,Tt+αer1TT−t−er0,TT,0=−er0,TT+maxer0,TT−er0,Tt+αer1TT−t,0.
In formula (7), −er0,TTrepresents a basic debt contract with credit risk.
Andmaxer0,TT−er0,Tt+αer1TT−t,0 represents an interest rate put option with an exercise price of K=er0,TT. The conditions of option execution are as follows:(8)er0,TT≥er0,Tt+αer1TT−t,r1,T≤r0,TT−lner0,Tt+αT−t.Corollary 4.
The necessary condition of option execution is the reduction of interest rate.Proof.
The conditions of option execution are as follows:(9)er0,TT≥er0,Tt+αer1TT−t.
If the interest rate increases, thenrt,T≥r0,T, and then er0,TT≤er0,T+rt,TT−t<er0,Tt+αer1TT−t, and the option will not be executed, expressed mathematically as follows:(10)r1,T≤r0,TT−lner0,Tt+αT−t<r0,TT−r0,TtT−t=r0,T.
That is to say, for the borrower, the execution of options leads to a fine, which needs to be made up by the reduction of future interest.Corollary 5.
The larger theα, the less likely the option will be executed.Proof.
The highest interest rate of option execution is(11)r=r0,TT−lner0,Tt+αT−t,∂r∂α=−1er0,Tt+αT−t<0.
Therefore, the larger theα is, the lower the interest rate is needed before the option can be executed.Corollary 6.
The larger thet is, the less likely the prepayment is.Proof.
The highest interest rate of option execution is(12)r=r0,TT−lner0,Tt+αT−t,∂r∂t=−lner0,Tt+αT−t2<0.
After the option is executed, the borrower can reapply for a loan with a term ofT−t and naturally get another option; that is, it can prepay, pay liquidated damages, and sign a new contract during the period of T−t. The execution price is er1TT−t. Therefore, for borrowers, a term loan can be decomposed into a basic debt and several American options. With the delay of time, the possibility of the option being executed is gradually reduced, and its value is gradually reduced. The value of embedded options in loans is reflected in a period of time shortly after the contract is signed.


## 3.2.1. Analysis of the Implementation and Related Characteristics of Embedded Options in Deposits

Theorem 1.
For depositors, time deposits can be decomposed into a fixed debt and several American interest rate call option bulls; for banks, time deposits can be decomposed into a fixed debt and several American interest rate call option bulls.
Suppose at 0, a depositor deposits a time deposit with a maturity ofTin the bank, the continuous compound interest rate is r0,T and the deposit amount is standardized to 1. If the withdrawal in advance does not occur, the final income of the deposit is er0,TT. Assuming that the investor’s investment period is from 0 to T, the depositor’s goal is to maximize the return from 0 to T.
Suppose at timet, the central bank suddenly announced to change the level of deposit interest rate, the changed current interest rate is rt,t and the fixed deposit interest rate from T to t period is rt,T. At this point,
Ifer0,TT≥ert,tt+rt,TT−t, the income from maintaining the original time deposit exceeds the transferred deposit, and the depositor chooses to maintain the original time deposit.
Ifer0,TT≤ert,tt+rt,TT−t, the income of withdrawal and deposit in advance exceeds that of maintaining the original time deposit, and the depositor chooses to withdraw and deposit in advance, and the amount of deposit in advance is ert,tt.
To sum up, at the moment, the depositor’s final income function can be expressed as follows:(1)maxer0,TT,ert,tt+rt,TT−t=er0,TT+maxert,tt+rt,TT−t−er0,TT,0=er0,TT+maxer0,TTert,tt+rt,TT−t−r0,TT−1,0≈er0,TT+er0,TTmaxrt,tt+rt,TT−t−r0,TT,0.
In formula (1), er0,TT is a basic bond contract, and er0,TTmaxrt,tt+rt,TT−t−r0,TT,0 is an interest rate call option with an exercise price of K=r0,TT [25, 26]. The conditions for the exercise of the option are as follows:(2)rt,tt+rt,T≥r0,TT.
By transforming formula (2), we get the following results:(3)rt,tt+rt,TT≥r0,T.
In addition, after the option is exercised, the depositor obtains a fixed deposit with a term ofT−t and naturally obtains an option, that is to say, during the T−t period, the depositor can withdraw and deposit in advance, and the exercise price of the option is rt,T. Therefore, for depositors, a time deposit can be decomposed into a basic debt and several American options. However, as time goes on, the possibility of the option being executed gradually decreases, and its value gradually decreases. The value of embedded options in time deposits is reflected in a period of time shortly after the contract is signed.Corollary 1.
The necessary condition for option execution is the rise of interest rate.Proof.
The conditions of option execution are as follows:(4)rt,tt+rt,T≥r0,TT,rt,T≥r0,TT−rt,ttT−t.
Because under normal conditionsrt,t≤r0,T,(5)rt,T≥r0,TT−rt,ttT−t≥r0,TT−r0,TtT−t=r0,T.
In other words, the depositor’s immediate income is negative, which needs to be compensated by the future excess income.Corollary 2.
The more the interest rate rises, the more likely it is to withdraw in advance.Proof.
The more the interest rate rises, the less the immediate loss of early withdrawal and the more the excess income in the future, so the more likely the future income is to make up for the immediate loss, the more likely it is to withdraw in advance.Corollary 3.
The larger thet, the less likely it is to withdraw in advance.Proof.
The condition of option execution is as follows:(6)rt,T≥r0,TT−rt,ttT−t≥r0,T−rt,ttT−t+r0,T.
Therefore, the minimum executive interest rate isrt=r0,T−rt,tt/T−t+r0,T.
Because∂rt/∂t=r0,T−rt,tT/T−t2≥0, the minimum exercise interest rate level increases with the extension of time, and the possibility of option execution also decreases.


## 3.2.2. Analysis on the Execution and Related Characteristics of Embedded Options in Loans

Theorem 2.
For the borrower, the loan can be decomposed into a basic debt and several American interest rate put options; for the lending bank, the loan can be divided into a basic creditor’s right and several short positions of American interest rate put options.
Suppose at time 0, an investor signs a loan contract with a bank with a term ofT, the interest rate of the continuous compound interest loan is r0,T, the loan amount is standardized as 1, and the penalty is α. If there is no prepayment, the total cost he has to pay at the end of the period is er0,TT. Assuming that the borrower’s borrowing period is constant from 0 to T, the goal is to minimize the cost of the period. The credit rating of the borrower has not changed during the term of the loan.
Suppose at timet, the central bank suddenly announced to change the interest rate level, and the interest rate of the loan with a term T−t is r1,T, then if er0,TT≥er0,Tt+αer1TT−t, the cost of maintaining the loan exceeds prepayment, payment of liquidated damages, and signing a new contract, and the borrower will choose prepayment, payment of liquidated damages, and signing a new contract, and the loan amount of the new contract is er0,Tt+α; if er0,TT≤er0,Tt+αer1TT−t, the cost of maintaining the original term loan is lower than prepayment, paying liquidated damages, and signing a new contract, and the borrower maintains the original loan contract.
To sum up, at timet, the borrower’s final cost function can be expressed as follows:(7)−miner0,Tt+αer1TT−t,er0,TT=−er0,TT−miner0,Tt+αer1TT−t−er0,TT,0=−er0,TT+maxer0,TT−er0,Tt+αer1TT−t,0.
In formula (7), −er0,TTrepresents a basic debt contract with credit risk.
Andmaxer0,TT−er0,Tt+αer1TT−t,0 represents an interest rate put option with an exercise price of K=er0,TT. The conditions of option execution are as follows:(8)er0,TT≥er0,Tt+αer1TT−t,r1,T≤r0,TT−lner0,Tt+αT−t.Corollary 4.
The necessary condition of option execution is the reduction of interest rate.Proof.
The conditions of option execution are as follows:(9)er0,TT≥er0,Tt+αer1TT−t.
If the interest rate increases, thenrt,T≥r0,T, and then er0,TT≤er0,T+rt,TT−t<er0,Tt+αer1TT−t, and the option will not be executed, expressed mathematically as follows:(10)r1,T≤r0,TT−lner0,Tt+αT−t<r0,TT−r0,TtT−t=r0,T.
That is to say, for the borrower, the execution of options leads to a fine, which needs to be made up by the reduction of future interest.Corollary 5.
The larger theα, the less likely the option will be executed.Proof.
The highest interest rate of option execution is(11)r=r0,TT−lner0,Tt+αT−t,∂r∂α=−1er0,Tt+αT−t<0.
Therefore, the larger theα is, the lower the interest rate is needed before the option can be executed.Corollary 6.
The larger thet is, the less likely the prepayment is.Proof.
The highest interest rate of option execution is(12)r=r0,TT−lner0,Tt+αT−t,∂r∂t=−lner0,Tt+αT−t2<0.
After the option is executed, the borrower can reapply for a loan with a term ofT−t and naturally get another option; that is, it can prepay, pay liquidated damages, and sign a new contract during the period of T−t. The execution price is er1TT−t. Therefore, for borrowers, a term loan can be decomposed into a basic debt and several American options. With the delay of time, the possibility of the option being executed is gradually reduced, and its value is gradually reduced. The value of embedded options in loans is reflected in a period of time shortly after the contract is signed.


## 3.3. The Mechanism of the Effect of Embedded Options on the Value of Bank Deposits and Loans

Loan is the main part of commercial bank’s assets, and time deposit is the main part of commercial bank’s liabilities. Therefore, when the loan of commercial bank can be paid in advance and the time deposit can be withdrawn in advance, that is, embedded option is implied in assets and liabilities, the value of commercial bank’s assets and liabilities will be affected by the existence of embedded option. Especially when the market interest rate changes, the duration and convexity of assets and liabilities will change their original characteristics because of the existence of embedded options. This paper uses the graph to illustrate [27].For asset items whose main part is loans, especially fixed rate loans, the rise of interest rate has no effect on them. However, when the interest rate falls belowR∗, the borrowers execute the embedded options and choose to prepay, and the absolute value of asset duration becomes smaller. When the interest rate falls to a certain extent, a large area of borrowers will prepay. This makes the duration of assets further smaller and the convexity becomes negative, which is reflected in the value-interest rate curve of assets in Figure 3; that is, the curve just begins to become flat until it bends in the opposite direction, which means that the value of assets just begins to rise with the decline of interest rate, but the further decline of interest rate will lead to the borrower’s exercise of rights and the bottleneck of the rise of asset value, with the upper limit of P∗; for the liability items whose main part is time deposit, the decrease in interest rate has no effect on them. However, when the interest rate rises above R∗, the depositor executes the embedded options and chooses to withdraw in advance, which leads to the decrease in the absolute value of the duration of the liability, which is reflected in the value-interest rate curve of the liability in Figure 4; that is, the curve becomes flat. This means that the value of the liability has a lower limit as interest rates rise.Figure 3
Value-interest rate curve of assets.Figure 4
Value-interest rate curve of liabilities.In order to more clearly explain the impact of the existence of embedded options on the net assets (the difference between assets and liabilities) of commercial banks, this paper uses the mathematical expression of net assets and the value-interest rate curve of assets and liabilities to explain in more detail. Usually, in order to prevent interest rate risk exposure, commercial banks match the duration convexity of asset items and liability items to make the slope of value-interest rate curve of asset and liability the same. Here, bank assets are represented byA, liabilities by L, and net assets by E. Then, dE=dA−dL. The derivative of this method is as follows:(13)dA=dAdRAdRA+12d2AdRA2dRA2,dL=dLdRLdRL+12d2LdRL2dRL2.LetA and L have the same initial interest rate, RA=RL=R, and dRA=dRL=dR. We can get DA=−dA/AdR, DL=−dL/LdR, CA=−d2A/AdR2, and CL=−d2L/LdR2, and then the dE expression is as follows:(14)dE=−DA−KDLAdR+12CA−KCLAdR2.In equation (14), DA is the duration of assets, DL is the duration of liabilities, CA is the convexity of assets, and CL is the convexity of liabilities. K=L/A is the ratio of assets to liabilities, also known as leverage. DGAP=DA−KDL can reflect the duration matching of assets and liabilities, which is called duration gap; CGAP=CA−KCL can reflect the convexity matching of assets and liabilities, which is called convexity gap. From the expression of formula (14), we can know that the main factors affecting the net assets of commercial banks are A, DGAP=DA−KDL, CGAP=CA−KCL, and R [28–30].Next, we use the value-interest rate curve of bank assets and liabilities in Figure5 to illustrate the impact on the bank’s net assets when the asset liability items are implied with embedded options. Assume that the initial net assets of the bank are greater than zero, so the value-interest rate curve of the assets is above the value-interest rate curve of the liabilities. At the beginning, without considering the existence of embedded options, commercial banks have matched the duration and convexity of assets and liabilities. When the interest rate changes little, that is, between R∗ and R∗∗, the effect of asset liability management is very good, the net assets remain stable, and there is no interest rate risk exposure. However, when the interest rate falls below R∗∗, the depositors will not be affected, the borrowers will repay in advance, and the net assets will begin to decline, especially when the interest rate falls below R1, the net asset is E1, significantly less than the initial net assets E∗; similarly, when the interest rate rises to R∗, the lender will not be affected, and the depositor withdraws the deposit in advance, especially when the interest rate rises to R2, the net assets of the bank become E2, significantly less than the initial net assets E∗ [31].Figure 5
Value-interest rate curve of bank assets and liabilities.If we analyze formula (14), when the interest rate falls further below R∗∗, the duration gap DGAP=DA−KDL changes from 0 to negative and the convexity gap DGAP=DA−KDL also changes from 0 to negative so that the net assets become smaller; when the interest rate rises above R∗, the duration gap DGAP=DA−KDL changes from 0 to positive and the convexity gap DGAP=DA−KDL changes from 0 to negative, so the net assets become smaller. That is to say, as long as the bank’s loan can be repaid in advance and the time deposit can be withdrawn in advance, that is, when the bank’s asset liability project is implied with embedded options, the net assets will be less no matter whether the interest rate rises or falls. However, commercial banks ignore this kind of embedded option with intrinsic value, which leads to that the real interest rate difference between deposit and loan is not the difference between loan interest rate and deposit interest rate on the surface, which leads to the hidden interest loss of commercial banks, and the difficulty of interest rate risk management will increase.


## 4. Selection of Interest Rate Model and Estimation of Model Parameters

### 4.1. Selection of Interest Rate Model

Embedded option is an interest rate derivative, so it is particularly important to describe the trend of interest rate. This paper combines the interest rate jump behavior caused by macro policy with the interest rate diffusion behavior caused by market behavior and chooses the CIR model with jump term to describe the change path of deposit and loan interest rate, in which the jump term can represent the sudden change of policy, the collapse of stock market, and the administrative control of interest rate [32].The CIR model considers that in a risk neutral world, the change of short-term interest rate follows the following order:(15)dr=αβ−rdt+σr1/2dWt.In equation (15), r represents the short-term instantaneous interest rate; α, β, and σ represent the average recovery speed of interest rate, the long-term average of interest rate, and the volatility of interest rate, respectively, and they are all constants. Wt is Brownian motion. In this paper, the jump term JdP is added to the CIR model, which is the jump-diffusion model. Its expression is as follows:(16)drt=αβ−rtdt+σtr1/2dWt+JdP.In equation (16), J represents a random variable with uniform distribution. In order to avoid the sudden depth jump causing negative interest rate, it is assumed that J follows the uniform distribution on m,n, and dP represents the Poisson distribution with jump intensity λ [33–35].


### 4.2. Estimation Method of Model Parameters

In the jump-diffusion model,Jis completely independent of parameters α, β, σ, and λ. Starting from the practicability and simplicity of parameter estimation, this paper calculates the existing interest rate data by Excel tool, obtains its parameters α, β, σ, and λ, and uses the maximum likelihood estimation method to estimate the parameters m,n in J.


#### 4.2.1. Estimation Method of Basic Parametersα, β, σ, and λ

The estimation methods of basic parametersα, β, σ, and λ are shown in Table 2.Table 2
Estimation method of basic parameters.
ParameterCalculation formulaIllustrationαy=αt+εFirst, the standard deviation of each term is calculated through the interest rate data, and then the straight line reflecting the relationship between the standard deviation and the term, and it can be obtained by linear regression. The slope of the straight line isα.ββ=1/n∑i=1nRiThe long-term average can be replaced by the average interest rate over the years, whereRi is the historical data value of the corresponding term interest rate.σσ=1/n∑i=1nRi−β2We usually use the volatility of the base market interest rate to replace the transient interest rate volatilityσ, which can be expressed by the standard deviation of interest rate in period N, where Ri and β have the same meaning as above.λλ=w/nλ is equal to the average number of events per unit time, where w is the number of interest rate jumps, n is the interval length, and its unit is month.


#### 4.2.2. Estimation Method of Parameterm,n in Uniform Distribution J

The maximum likelihood estimation (MLE) method is used to estimate the jump parameters. The MLE method is a widely used method for estimating values. It is based on the maximum likelihood principle; that is, there are several possible resultsA, B, and C, in a random experiment. If result A appears in an experiment, it can be considered that pA is larger. The principle is as follows.Let the family of probability functions of populationξ be fx;θ,θ∈Θ, where θ′=θ1,θ2,…,θk is the parameter vector to be estimated in K dimension and x1,…,xn=x is an observed value of subsample ξ1,…,ξn=ξ′, then the probability of subsample ξ falling into the field of point x′ is ∏i=1nfxi,θΔxi. It can be seen that this probability will be affected by the change of θ, which is the function of θ. The principle of the maximum likelihood method is to select the parameter value θ^ which makes the probability ∏i=1nfxi,θΔxi of the subsample in the field of observation value x1,…,xn reach the maximum as the estimation of θ. That is to say, for a fixed x1,…,xn, θ^ is selected.(17)∏i=1nfxi;θ^=SUPθ∈Θ∏i=1nfxi;θ.Based on the above principle, we discretize equation (16) by Euler as follows:(18)rt−rt−1=αβ−rt−1+εt.Then, we define the parameter vector to be estimatedθ=α,β,σ,m,n,q, where m and n are the mean distribution parameters of the jump process and are the parameter values to be determined in our method; α, β, σ, and q represent the mean recovery speed of the jump-diffusion model, the long-term mean of interest rate, and volatility and jump frequency, respectively, which can be determined in the above parameter estimation methods. Suppose when the jump occurs, i.e., dp=1, Eεt=m+n/2 and Dεt2=σ2rt−1+n−m2/12; when the jump does not occur, i.e., dp=0, Eεt=0 and Dεt2=σ2rt−1.If the probability of jump isq=λdt, the transition probability density of jump diffusion is approximately as follows:(19)frt\rt−1;θ=q12πσ2rt−1+n−m2/12exp−rt−rt−1−αβ+αrt−1−m+n/222σ2rt−1+n−m2/12+1−q12πσ2rt−1exp−rt−rt−1−αβ+αrt−122σ2rt−1.From the stationarity and Markov property of Poisson process, the likelihood function is obtained:(20)lθ=∏i=1nfrt\rt−1;θ.The likelihood functionlθ is logarithm. The results are as follows:(21)Lθ=∑i=1nlnfrt\rt−1;θ.In order to obtain the maximum value of equation (21), let the likelihood equation be as follows:(22)∂Lθ∂θj=0,j=1,2.The parameter values of MLE can be obtained by solving and verifying whether it is the maximum point [36, 37].


### 4.3. Empirical Analysis of Model Parameters

#### 4.3.1. Parameter Empirical Simulation

The choice of interest rate data follows the principle of trading volume, stability, and relevance. Considering the timeliness of the data, this paper only selects the benchmark interest rate of deposits and loans published by the people’s Bank of China for the last 20 years, that is, 240 months from January 1, 2001, to January 1, 2021.(1) Long-Term Return Average of Interest Rateα. By linear regression of the data in Table 3, we can get the slope of the curve of volatility and maturity, which is the long-term return mean of interest rate α. The results of EViews regression are shown in Figure 6. The linear function is as follows: (23)y=0.0005x+0.0062.Table 3
Average interest rate of time deposit of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearβ0.00520.02250.02560.02830.03470.04080.0473Figure 6
EViews regression results.According to equation (23), the parameter α=0.0005.(2) Long-Term Average Interest Rateβ. Using the parameter estimation method in Table 2, the results of β are shown in Table 3.(3) Volatilityσ. Using the parameter estimation method in Table 2, the results of σ are shown in Table 4.Table 4
Time deposit interest rate volatility of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearσ0.001890.005910.006340.006930.007510.007970.00826(4) Jump Strengthλ. According to the statistics of the change times of time deposit interest rate in recent ten years, the jump times of time deposit interest rate of each term are shown in Table 5.Table 5
Jump times of time deposit interest rate of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearHop count14222526262525Since the monthly data are 240, the correspondingλ value can be obtained by dividing the number of jumps by 240. The results of λ are shown in Table 6.Table 6
Jump intensity of time deposit interest rate of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearλ0.05830.09170.1040.1080.1070.1040.104(5) Estimation of Parameterm,nin Mean Distribution J. Taking the 5-year fixed deposit as an example, the maximum likelihood method is used to calculate with MATLAB. The deviation less than 10−7 is taken as the standard to estimate the parameters, and the result shows that m=−0.016 and n=0.016.After estimating all the model parameters, we can find that the interest rate of five-year fixed deposit follows a stochastic process:(24)drt=0.00050.0473−rtdt+0.00826rtdwt+Jdp.The model of equation (24) is discretized, and the results are as follows:(25)Rt+1=Rt+0.000500.0473−RtΔt+0.00826Rtdwt+Jdp.In equation (25), J follows the mean distribution of (−0.016, 0.016) and dP represents the Poisson distribution of jump intensity λ=0.104.Using the same parameter estimation method and the same interval of loan interest rate, we estimate the jump-diffusion model of 10-year loan interest rate as follows:(26)Rt+1=Rt+0.000540.0692−RtΔt+0.00947Rtdwt+Jdp.In equation (26), J follows the mean distribution of (−0.014, 0.014) and dP represents the Poisson distribution of jump intensity λ=0.104.


#### 4.3.2. Empirical Test of the Model

(1) The Process of Empirical Test. Taking the five-year deposit interest rate as an example, this paper tests the model from the aspects of mean reversion, volatility, and jump and compares the CIR model with jump term and the CIR model without jump term with the selected interest rate sample data to analyze the standard error of the two models. The specific process is as follows:Firstly, the sample path of the underlying variable in the validity period is simulated under the risk neutral measure. The discrete stochastic model has been given in the previous paper:(27)Rt+1=Rt+0.000500.0473−RtΔt+0.00826Rtdwt+Jdp.In equation (27), Rt is the existing five-year fixed deposit interest rate of 0.0275; Δt is 1/12; Brownian motion dwt can be expressed by ΔtZ, where Z∼N0,1 normal distribution can realize the sample sampling of standard normal distribution in MATLAB; J and dP can be completed by the identification sampling of commands in MATLAB.Secondly, the corresponding function values on the sample path are estimated. According to the properties of the underlying variables, the absolute value of the deviation between the simulated data and the real data on each sample path is estimated:(28)di=Rim−Ris.In equation (28), Rim is the simulated value of interest rate and Ris is the real data of interest rate.Finally, the estimated value of the underlying variable is estimated. The estimated value of the underlying variable is obtained by averaging the corresponding function values on the simulated sample path:(29)d=1N∑i=1Ndi.In this paper, the Monte Carlo simulation method is used to simulate the interest rate model for many times to test the mean reversion of the model, the rationality of the limit value, and the error comparison between the diffusion model and the jump-diffusion model.Figure 7
A simulated interest rate path of the jump-diffusion model.(2) Empirical Test ResultFrom the two results in Figure8, we can see that the corresponding error value of the jump-diffusion model is 0.0417 and the corresponding error value of the interest rate model without jump term is 0.0462. Therefore, it can be seen that the jump-diffusion model is more in line with the trend of China’s benchmark deposit and loan interest rate.Figure 8
Comparison of simulated value error under two interest rate models: (a) in presence of jump term; (b) without jump term.
(a)(b)


## 4.1. Selection of Interest Rate Model

Embedded option is an interest rate derivative, so it is particularly important to describe the trend of interest rate. This paper combines the interest rate jump behavior caused by macro policy with the interest rate diffusion behavior caused by market behavior and chooses the CIR model with jump term to describe the change path of deposit and loan interest rate, in which the jump term can represent the sudden change of policy, the collapse of stock market, and the administrative control of interest rate [32].The CIR model considers that in a risk neutral world, the change of short-term interest rate follows the following order:(15)dr=αβ−rdt+σr1/2dWt.In equation (15), r represents the short-term instantaneous interest rate; α, β, and σ represent the average recovery speed of interest rate, the long-term average of interest rate, and the volatility of interest rate, respectively, and they are all constants. Wt is Brownian motion. In this paper, the jump term JdP is added to the CIR model, which is the jump-diffusion model. Its expression is as follows:(16)drt=αβ−rtdt+σtr1/2dWt+JdP.In equation (16), J represents a random variable with uniform distribution. In order to avoid the sudden depth jump causing negative interest rate, it is assumed that J follows the uniform distribution on m,n, and dP represents the Poisson distribution with jump intensity λ [33–35].


## 4.2. Estimation Method of Model Parameters

In the jump-diffusion model,Jis completely independent of parameters α, β, σ, and λ. Starting from the practicability and simplicity of parameter estimation, this paper calculates the existing interest rate data by Excel tool, obtains its parameters α, β, σ, and λ, and uses the maximum likelihood estimation method to estimate the parameters m,n in J.


### 4.2.1. Estimation Method of Basic Parametersα, β, σ, and λ

The estimation methods of basic parametersα, β, σ, and λ are shown in Table 2.Table 2
Estimation method of basic parameters.
ParameterCalculation formulaIllustrationαy=αt+εFirst, the standard deviation of each term is calculated through the interest rate data, and then the straight line reflecting the relationship between the standard deviation and the term, and it can be obtained by linear regression. The slope of the straight line isα.ββ=1/n∑i=1nRiThe long-term average can be replaced by the average interest rate over the years, whereRi is the historical data value of the corresponding term interest rate.σσ=1/n∑i=1nRi−β2We usually use the volatility of the base market interest rate to replace the transient interest rate volatilityσ, which can be expressed by the standard deviation of interest rate in period N, where Ri and β have the same meaning as above.λλ=w/nλ is equal to the average number of events per unit time, where w is the number of interest rate jumps, n is the interval length, and its unit is month.


### 4.2.2. Estimation Method of Parameterm,n in Uniform Distribution J

The maximum likelihood estimation (MLE) method is used to estimate the jump parameters. The MLE method is a widely used method for estimating values. It is based on the maximum likelihood principle; that is, there are several possible resultsA, B, and C, in a random experiment. If result A appears in an experiment, it can be considered that pA is larger. The principle is as follows.Let the family of probability functions of populationξ be fx;θ,θ∈Θ, where θ′=θ1,θ2,…,θk is the parameter vector to be estimated in K dimension and x1,…,xn=x is an observed value of subsample ξ1,…,ξn=ξ′, then the probability of subsample ξ falling into the field of point x′ is ∏i=1nfxi,θΔxi. It can be seen that this probability will be affected by the change of θ, which is the function of θ. The principle of the maximum likelihood method is to select the parameter value θ^ which makes the probability ∏i=1nfxi,θΔxi of the subsample in the field of observation value x1,…,xn reach the maximum as the estimation of θ. That is to say, for a fixed x1,…,xn, θ^ is selected.(17)∏i=1nfxi;θ^=SUPθ∈Θ∏i=1nfxi;θ.Based on the above principle, we discretize equation (16) by Euler as follows:(18)rt−rt−1=αβ−rt−1+εt.Then, we define the parameter vector to be estimatedθ=α,β,σ,m,n,q, where m and n are the mean distribution parameters of the jump process and are the parameter values to be determined in our method; α, β, σ, and q represent the mean recovery speed of the jump-diffusion model, the long-term mean of interest rate, and volatility and jump frequency, respectively, which can be determined in the above parameter estimation methods. Suppose when the jump occurs, i.e., dp=1, Eεt=m+n/2 and Dεt2=σ2rt−1+n−m2/12; when the jump does not occur, i.e., dp=0, Eεt=0 and Dεt2=σ2rt−1.If the probability of jump isq=λdt, the transition probability density of jump diffusion is approximately as follows:(19)frt\rt−1;θ=q12πσ2rt−1+n−m2/12exp−rt−rt−1−αβ+αrt−1−m+n/222σ2rt−1+n−m2/12+1−q12πσ2rt−1exp−rt−rt−1−αβ+αrt−122σ2rt−1.From the stationarity and Markov property of Poisson process, the likelihood function is obtained:(20)lθ=∏i=1nfrt\rt−1;θ.The likelihood functionlθ is logarithm. The results are as follows:(21)Lθ=∑i=1nlnfrt\rt−1;θ.In order to obtain the maximum value of equation (21), let the likelihood equation be as follows:(22)∂Lθ∂θj=0,j=1,2.The parameter values of MLE can be obtained by solving and verifying whether it is the maximum point [36, 37].


## 4.2.1. Estimation Method of Basic Parametersα, β, σ, and λ

The estimation methods of basic parametersα, β, σ, and λ are shown in Table 2.Table 2
Estimation method of basic parameters.
ParameterCalculation formulaIllustrationαy=αt+εFirst, the standard deviation of each term is calculated through the interest rate data, and then the straight line reflecting the relationship between the standard deviation and the term, and it can be obtained by linear regression. The slope of the straight line isα.ββ=1/n∑i=1nRiThe long-term average can be replaced by the average interest rate over the years, whereRi is the historical data value of the corresponding term interest rate.σσ=1/n∑i=1nRi−β2We usually use the volatility of the base market interest rate to replace the transient interest rate volatilityσ, which can be expressed by the standard deviation of interest rate in period N, where Ri and β have the same meaning as above.λλ=w/nλ is equal to the average number of events per unit time, where w is the number of interest rate jumps, n is the interval length, and its unit is month.


## 4.2.2. Estimation Method of Parameterm,n in Uniform Distribution J

The maximum likelihood estimation (MLE) method is used to estimate the jump parameters. The MLE method is a widely used method for estimating values. It is based on the maximum likelihood principle; that is, there are several possible resultsA, B, and C, in a random experiment. If result A appears in an experiment, it can be considered that pA is larger. The principle is as follows.Let the family of probability functions of populationξ be fx;θ,θ∈Θ, where θ′=θ1,θ2,…,θk is the parameter vector to be estimated in K dimension and x1,…,xn=x is an observed value of subsample ξ1,…,ξn=ξ′, then the probability of subsample ξ falling into the field of point x′ is ∏i=1nfxi,θΔxi. It can be seen that this probability will be affected by the change of θ, which is the function of θ. The principle of the maximum likelihood method is to select the parameter value θ^ which makes the probability ∏i=1nfxi,θΔxi of the subsample in the field of observation value x1,…,xn reach the maximum as the estimation of θ. That is to say, for a fixed x1,…,xn, θ^ is selected.(17)∏i=1nfxi;θ^=SUPθ∈Θ∏i=1nfxi;θ.Based on the above principle, we discretize equation (16) by Euler as follows:(18)rt−rt−1=αβ−rt−1+εt.Then, we define the parameter vector to be estimatedθ=α,β,σ,m,n,q, where m and n are the mean distribution parameters of the jump process and are the parameter values to be determined in our method; α, β, σ, and q represent the mean recovery speed of the jump-diffusion model, the long-term mean of interest rate, and volatility and jump frequency, respectively, which can be determined in the above parameter estimation methods. Suppose when the jump occurs, i.e., dp=1, Eεt=m+n/2 and Dεt2=σ2rt−1+n−m2/12; when the jump does not occur, i.e., dp=0, Eεt=0 and Dεt2=σ2rt−1.If the probability of jump isq=λdt, the transition probability density of jump diffusion is approximately as follows:(19)frt\rt−1;θ=q12πσ2rt−1+n−m2/12exp−rt−rt−1−αβ+αrt−1−m+n/222σ2rt−1+n−m2/12+1−q12πσ2rt−1exp−rt−rt−1−αβ+αrt−122σ2rt−1.From the stationarity and Markov property of Poisson process, the likelihood function is obtained:(20)lθ=∏i=1nfrt\rt−1;θ.The likelihood functionlθ is logarithm. The results are as follows:(21)Lθ=∑i=1nlnfrt\rt−1;θ.In order to obtain the maximum value of equation (21), let the likelihood equation be as follows:(22)∂Lθ∂θj=0,j=1,2.The parameter values of MLE can be obtained by solving and verifying whether it is the maximum point [36, 37].


## 4.3. Empirical Analysis of Model Parameters

### 4.3.1. Parameter Empirical Simulation

The choice of interest rate data follows the principle of trading volume, stability, and relevance. Considering the timeliness of the data, this paper only selects the benchmark interest rate of deposits and loans published by the people’s Bank of China for the last 20 years, that is, 240 months from January 1, 2001, to January 1, 2021.(1) Long-Term Return Average of Interest Rateα. By linear regression of the data in Table 3, we can get the slope of the curve of volatility and maturity, which is the long-term return mean of interest rate α. The results of EViews regression are shown in Figure 6. The linear function is as follows: (23)y=0.0005x+0.0062.Table 3
Average interest rate of time deposit of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearβ0.00520.02250.02560.02830.03470.04080.0473Figure 6
EViews regression results.According to equation (23), the parameter α=0.0005.(2) Long-Term Average Interest Rateβ. Using the parameter estimation method in Table 2, the results of β are shown in Table 3.(3) Volatilityσ. Using the parameter estimation method in Table 2, the results of σ are shown in Table 4.Table 4
Time deposit interest rate volatility of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearσ0.001890.005910.006340.006930.007510.007970.00826(4) Jump Strengthλ. According to the statistics of the change times of time deposit interest rate in recent ten years, the jump times of time deposit interest rate of each term are shown in Table 5.Table 5
Jump times of time deposit interest rate of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearHop count14222526262525Since the monthly data are 240, the correspondingλ value can be obtained by dividing the number of jumps by 240. The results of λ are shown in Table 6.Table 6
Jump intensity of time deposit interest rate of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearλ0.05830.09170.1040.1080.1070.1040.104(5) Estimation of Parameterm,nin Mean Distribution J. Taking the 5-year fixed deposit as an example, the maximum likelihood method is used to calculate with MATLAB. The deviation less than 10−7 is taken as the standard to estimate the parameters, and the result shows that m=−0.016 and n=0.016.After estimating all the model parameters, we can find that the interest rate of five-year fixed deposit follows a stochastic process:(24)drt=0.00050.0473−rtdt+0.00826rtdwt+Jdp.The model of equation (24) is discretized, and the results are as follows:(25)Rt+1=Rt+0.000500.0473−RtΔt+0.00826Rtdwt+Jdp.In equation (25), J follows the mean distribution of (−0.016, 0.016) and dP represents the Poisson distribution of jump intensity λ=0.104.Using the same parameter estimation method and the same interval of loan interest rate, we estimate the jump-diffusion model of 10-year loan interest rate as follows:(26)Rt+1=Rt+0.000540.0692−RtΔt+0.00947Rtdwt+Jdp.In equation (26), J follows the mean distribution of (−0.014, 0.014) and dP represents the Poisson distribution of jump intensity λ=0.104.


### 4.3.2. Empirical Test of the Model

(1) The Process of Empirical Test. Taking the five-year deposit interest rate as an example, this paper tests the model from the aspects of mean reversion, volatility, and jump and compares the CIR model with jump term and the CIR model without jump term with the selected interest rate sample data to analyze the standard error of the two models. The specific process is as follows:Firstly, the sample path of the underlying variable in the validity period is simulated under the risk neutral measure. The discrete stochastic model has been given in the previous paper:(27)Rt+1=Rt+0.000500.0473−RtΔt+0.00826Rtdwt+Jdp.In equation (27), Rt is the existing five-year fixed deposit interest rate of 0.0275; Δt is 1/12; Brownian motion dwt can be expressed by ΔtZ, where Z∼N0,1 normal distribution can realize the sample sampling of standard normal distribution in MATLAB; J and dP can be completed by the identification sampling of commands in MATLAB.Secondly, the corresponding function values on the sample path are estimated. According to the properties of the underlying variables, the absolute value of the deviation between the simulated data and the real data on each sample path is estimated:(28)di=Rim−Ris.In equation (28), Rim is the simulated value of interest rate and Ris is the real data of interest rate.Finally, the estimated value of the underlying variable is estimated. The estimated value of the underlying variable is obtained by averaging the corresponding function values on the simulated sample path:(29)d=1N∑i=1Ndi.In this paper, the Monte Carlo simulation method is used to simulate the interest rate model for many times to test the mean reversion of the model, the rationality of the limit value, and the error comparison between the diffusion model and the jump-diffusion model.Figure 7
A simulated interest rate path of the jump-diffusion model.(2) Empirical Test ResultFrom the two results in Figure8, we can see that the corresponding error value of the jump-diffusion model is 0.0417 and the corresponding error value of the interest rate model without jump term is 0.0462. Therefore, it can be seen that the jump-diffusion model is more in line with the trend of China’s benchmark deposit and loan interest rate.Figure 8
Comparison of simulated value error under two interest rate models: (a) in presence of jump term; (b) without jump term.
(a)(b)


## 4.3.1. Parameter Empirical Simulation

The choice of interest rate data follows the principle of trading volume, stability, and relevance. Considering the timeliness of the data, this paper only selects the benchmark interest rate of deposits and loans published by the people’s Bank of China for the last 20 years, that is, 240 months from January 1, 2001, to January 1, 2021.(1) Long-Term Return Average of Interest Rateα. By linear regression of the data in Table 3, we can get the slope of the curve of volatility and maturity, which is the long-term return mean of interest rate α. The results of EViews regression are shown in Figure 6. The linear function is as follows: (23)y=0.0005x+0.0062.Table 3
Average interest rate of time deposit of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearβ0.00520.02250.02560.02830.03470.04080.0473Figure 6
EViews regression results.According to equation (23), the parameter α=0.0005.(2) Long-Term Average Interest Rateβ. Using the parameter estimation method in Table 2, the results of β are shown in Table 3.(3) Volatilityσ. Using the parameter estimation method in Table 2, the results of σ are shown in Table 4.Table 4
Time deposit interest rate volatility of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearσ0.001890.005910.006340.006930.007510.007970.00826(4) Jump Strengthλ. According to the statistics of the change times of time deposit interest rate in recent ten years, the jump times of time deposit interest rate of each term are shown in Table 5.Table 5
Jump times of time deposit interest rate of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearHop count14222526262525Since the monthly data are 240, the correspondingλ value can be obtained by dividing the number of jumps by 240. The results of λ are shown in Table 6.Table 6
Jump intensity of time deposit interest rate of each term.
Interest rate varietiesCurrent3 months6 months1 year2 year3 year5 yearλ0.05830.09170.1040.1080.1070.1040.104(5) Estimation of Parameterm,nin Mean Distribution J. Taking the 5-year fixed deposit as an example, the maximum likelihood method is used to calculate with MATLAB. The deviation less than 10−7 is taken as the standard to estimate the parameters, and the result shows that m=−0.016 and n=0.016.After estimating all the model parameters, we can find that the interest rate of five-year fixed deposit follows a stochastic process:(24)drt=0.00050.0473−rtdt+0.00826rtdwt+Jdp.The model of equation (24) is discretized, and the results are as follows:(25)Rt+1=Rt+0.000500.0473−RtΔt+0.00826Rtdwt+Jdp.In equation (25), J follows the mean distribution of (−0.016, 0.016) and dP represents the Poisson distribution of jump intensity λ=0.104.Using the same parameter estimation method and the same interval of loan interest rate, we estimate the jump-diffusion model of 10-year loan interest rate as follows:(26)Rt+1=Rt+0.000540.0692−RtΔt+0.00947Rtdwt+Jdp.In equation (26), J follows the mean distribution of (−0.014, 0.014) and dP represents the Poisson distribution of jump intensity λ=0.104.


## 4.3.2. Empirical Test of the Model

(1) The Process of Empirical Test. Taking the five-year deposit interest rate as an example, this paper tests the model from the aspects of mean reversion, volatility, and jump and compares the CIR model with jump term and the CIR model without jump term with the selected interest rate sample data to analyze the standard error of the two models. The specific process is as follows:Firstly, the sample path of the underlying variable in the validity period is simulated under the risk neutral measure. The discrete stochastic model has been given in the previous paper:(27)Rt+1=Rt+0.000500.0473−RtΔt+0.00826Rtdwt+Jdp.In equation (27), Rt is the existing five-year fixed deposit interest rate of 0.0275; Δt is 1/12; Brownian motion dwt can be expressed by ΔtZ, where Z∼N0,1 normal distribution can realize the sample sampling of standard normal distribution in MATLAB; J and dP can be completed by the identification sampling of commands in MATLAB.Secondly, the corresponding function values on the sample path are estimated. According to the properties of the underlying variables, the absolute value of the deviation between the simulated data and the real data on each sample path is estimated:(28)di=Rim−Ris.In equation (28), Rim is the simulated value of interest rate and Ris is the real data of interest rate.Finally, the estimated value of the underlying variable is estimated. The estimated value of the underlying variable is obtained by averaging the corresponding function values on the simulated sample path:(29)d=1N∑i=1Ndi.In this paper, the Monte Carlo simulation method is used to simulate the interest rate model for many times to test the mean reversion of the model, the rationality of the limit value, and the error comparison between the diffusion model and the jump-diffusion model.Figure 7
A simulated interest rate path of the jump-diffusion model.(2) Empirical Test ResultFrom the two results in Figure8, we can see that the corresponding error value of the jump-diffusion model is 0.0417 and the corresponding error value of the interest rate model without jump term is 0.0462. Therefore, it can be seen that the jump-diffusion model is more in line with the trend of China’s benchmark deposit and loan interest rate.Figure 8
Comparison of simulated value error under two interest rate models: (a) in presence of jump term; (b) without jump term.
(a)(b)


## 5. Pricing of Embedded Options of Deposits and Loans Based on Monte Carlo Simulation

### 5.1. Application of Monte Carlo Simulation in Option Pricing

#### 5.1.1. Analysis of the Optimal Execution Strategy of Embedded Options

Before analyzing and pricing the embedded option, it is regarded as a standard American option, and the embedded option has the characteristics of early execution. The standard of early execution is that the value obtained by executing the option is greater than or equal to the value obtained by not executing in advance. Because the embedded option has a strong path dependence, it will have a different value because of the different interest rate changes, so there may be many opportunities for early execution in the duration of the option, but only one time is optimal, making the value of the option maximum, and the value at this time is the real value of the option. Therefore, the steps to find the optimal execution strategy are as follows:(1)Based on the change of interest rate, calculate the corresponding interest income or interest expense at each time point. For deposits, when the interest rate changes, the interest income at this time point is the sum of the demand interest before the deposit transfer point and the time deposit income after the deposit transfer point; for loans, when the interest rate changes, the interest expense at this time point is the sum of the interest expense before refinancing, the penalty for breach of contract, and the interest expense after refinancing.(2)
Optimal Execution Strategy. A point on the path of interest rate change makes the interest income maximum or the interest expense minimum, then this point is the optimal execution strategy point, and the American option value is the option value of this point.The pricing of embedded option is to find the option value of the optimal execution strategy and discount the expected return at that time to the current time. Therefore, we can know that for a path of interest rate change, we can find the best execution strategy at each execution time from the back to the front. Through repeated deduction, we finally get the best execution strategy and the corresponding discount value of income on each path.


#### 5.1.2. Application of Monte Carlo Simulation

(1) Monte Carlo Simulation Technology. The embedded option has the characteristics of multidimension and complexity, so the pricing method of numerical analysis is suitable. The numerical analysis model can effectively simulate various variables in reality, and this method has become an indispensable tool for embedded option valuation. The implementation of numerical analysis includes network analysis and Monte Carlo simulation. Network analysis is not suitable for pricing embedded options. Firstly, many fundamental characteristics of embedded options cannot be described by network methods, such as the jump-diffusion model describing the trend of interest rate change, which cannot be simulated by network methods. Secondly, for some high-dimensional weighted bonds with more underlying variables, it is impossible to price them. Thirdly, the discrete form of low dimensional equivalency is its fatal defect, that is to say, its default interest rate will only change once in the above equal time interval, while the actual situation is that the embedded option can be executed at any time before the maturity date of the option.When the dimension of financial derivative securities is relatively high, Monte Carlo simulation technology will become the only effective numerical analysis method; that is, it can effectively solve the “dimension disaster” problem. The essence of Monte Carlo simulation is a numerical analysis method to obtain the actual estimation of the mathematical expectation of random distribution by approximately replacing the overall expectation of random distribution with the sample mean of random sampling. Sometimes, when the pricing problem of complex options cannot be solved by analytical or numerical methods, Monte Carlo simulation can obtain the interest rate path and the final profit and loss of options by inputting assumed parameters and thousands of computer simulations, then obtain the distribution of option value, and finally use the average result to calculate the price of complex options [38–40].(2) Application of Monte Carlo Simulation Technology in Option Pricing. According to the nonarbitrage asset pricing theory, the price of financial derivative securities is the discounted expected value of its profit and loss function, and the expected value is the expected value under the risk neutral probability measure. According to the theory of probability statistics and stochastic analysis, these expected values can also be expressed in the form of integral in a specific region. Monte Carlo simulation is to estimate the price of financial derivative securities by simulating these mathematical expectations or integrals. The steps are as follows:SupposeSt is the price of the underlying asset at time t, r is the risk-free market interest rate, σ is the volatility of the underlying asset return, T is the maturity of the option, and X is the exercise price of the option. Under the background of Black–Scholes, the underlying asset price St satisfies the geometric Brownian motion:(30)dSt=rStdt+σStdWt.In equation (30), Wt is the geometric Brownian motion under the risk neutral probability measure, then the return function at the maturity date of the option can be expressed as follows:(31)C=e−rTES¯−X+.S¯ in equation (31) is the average price of the underlying asset in 0,T time.The specific process of Monte Carlo simulation to realize option value calculation is as follows:(1)
Simulaten paths of underlying asset prices and then get the average prices of each path S1¯,S2¯,…,Sn¯;(2)
LetY=ES¯−X+, calculate the period weight of the average price on each path and get the value Y1,Y2,…,Yn;(3)
LetCj=e−rTYj,n=1,2,…,n and average the value of C to get C¯; that is, C¯=1/n∑i=1nCj; this is the option value obtained by simulation.


### 5.2. Pricing Steps of Embedded Options for Commercial Banks’ Deposits and Loans

#### 5.2.1. Pricing Steps of Embedded Options in Deposits

Using Monte Carlo simulation to price the embedded options in commercial bank deposits, the steps are as follows:(1) Simulate the Sample Path of the Underlying Variable within the Validity Period under the Risk Neutral Measure. The underlying variable here is deposit interest rate, and the interest rate model describing the change of deposit interest rate is the jump-diffusion model: (32)drt=αβ−rtdt+σtr1/2dWt+JdP.Equation (32) is discretized and simulated. Take Δt=1 month to simulate the possible path of interest rate after n years (12∗n months) in the current initial interest rate state. By using this simulation method repeatedly, thousands of interest rate paths can be simulated to meet the needs of Monte Carlo simulation.(2) Calculation of Embedded Option Value on Simulation Path of Each Sample. Each interest rate sample path may involve the early withdrawal of deposits, so there may be multiple options on one interest rate path. However, based on the optimal strategy problem of American option execution, there is only one optimal execution point, and then the discount of the option value of the optimal execution point is the value of the embedded option on the interest rate path: (33)C=exp−rtmaxexph0,T⋅t⋅expRt,T⋅T−t−expR0,T⋅T,0.In equation (33), C is the value of the option, r is the risk-free interest rate, R0,T is the interest rate of the previous contract, h0,T is the current interest rate, Rt,T is the interest rate after the change of time t, and T is the deposit term.(3) Average the Corresponding Option Values on the Simulated Sample Path. The average value of the corresponding options on the simulated sample path is obtained, and the estimated value of the embedded options in the time deposit is obtained:(34)C¯=1n∑i=1nCi.


#### 5.2.2. Pricing Steps of Embedded Options in Loans

The pricing process of embedded options in loans is similar to the calculation process of embedded option value in time deposits; except that the penalty for prepayment is to pay liquidated damages, that is, to charge a certain proportion of liquidated damages from the prepayment loans, and the proportion coefficient is expressed byα. From 2019, CCB’s prepayment policy for housing loans is to charge 3% of prepayment for less than one year, 2% for one to two years, and 1% for more than two years.Although the calculation of embedded options in loans is one more penaltyα than that in deposits, the Monte Carlo simulation method can still be used, and the calculation steps are as follows:(1) Simulate the Change Path of Loan Interest Rate. The 10-year loan interest rate model is shown in formula (26).(2) Calculation of Embedded Option Value on Simulation Path of Each Sample. Due to the existence of penalty α, the calculation of embedded options in loans is slightly different from that in deposits: (35)C=exp−rtmaxexpr0,T⋅T−α+expr0,t⋅t⋅exprt,T⋅T−t,0.In equation (35), C is the embedded option value, r is the risk-free interest rate, α is the penalty when prepaying the loan, t is the time when the option is exercised, T is the loan term, r0,T is the interest rate agreed in the previous loan contract, and rt,T is the interest rate in the new loan contract signed at t.(3) Average the Corresponding Option Values on the Simulated Sample Path. By averaging the corresponding options on the simulated sample path, we can get the estimated value of the embedded options in loans:(36)C¯=1n∑i=1nCi.


### 5.3. Pricing of Embedded Options Based on Monte Carlo Simulation

#### 5.3.1. An Empirical Study of Embedded Option Pricing in Deposits

We take the latest five-year fixed deposit interest rate of 0.0275 on January 1, 2021, as the initial interest rate and the one-year fixed deposit interest rate of 0.015 as the discount rate, and the deposit amount is standardized to 1, and the embedded option value is calculated in the five-year fixed deposit through simulation.We have analyzed the selection of the interest rate model and the principle and steps of the Monte Carlo simulation method. In order to make the simulation results more accurate, we need to use MATLAB software. We know from the characteristics of Monte Carlo simulation that the accuracy of the simulation results has a great relationship with the simulation times. The more the simulation times are, the more accurate the results will be. Therefore, 100000 times of simulation are carried out here, and the simulation results are shown in Figure9.It can be clearly seen from the convergence chart of option value in Figure9 that the value of embedded option in 5-year fixed deposit of commercial banks is 0.0012. In other words, the real deposit interest rate of commercial banks is 0.0012 more than the nominal deposit interest rate.Figure 9
Simulation results of embedded option value in 5-year fixed deposit.


#### 5.3.2. Empirical Analysis of Embedded Option Pricing in Loans

We take the latest 10-year loan interest rate of 0.049 on January 1, 2021, as the initial interest rate and the one-year fixed deposit interest rate of 0.015 as the discount rate, and the loan amount is standardized to 1. Assuming that the customer needs to pay 1% penalty for prepayment, we calculate the embedded option value in the 10-year loan through simulation and simulate 100000 knots. The results are as follows.It can be clearly seen from the convergence chart of option value in Figure10 that the value of embedded option in 10-year loan of commercial bank is 0.0028. In other words, the 10-year real loan interest rate of commercial banks is 0.0028 less than the nominal loan interest rate.Figure 10
Simulation results of embedded option value in 10-year loan.


## 5.1. Application of Monte Carlo Simulation in Option Pricing

### 5.1.1. Analysis of the Optimal Execution Strategy of Embedded Options

Before analyzing and pricing the embedded option, it is regarded as a standard American option, and the embedded option has the characteristics of early execution. The standard of early execution is that the value obtained by executing the option is greater than or equal to the value obtained by not executing in advance. Because the embedded option has a strong path dependence, it will have a different value because of the different interest rate changes, so there may be many opportunities for early execution in the duration of the option, but only one time is optimal, making the value of the option maximum, and the value at this time is the real value of the option. Therefore, the steps to find the optimal execution strategy are as follows:(1)Based on the change of interest rate, calculate the corresponding interest income or interest expense at each time point. For deposits, when the interest rate changes, the interest income at this time point is the sum of the demand interest before the deposit transfer point and the time deposit income after the deposit transfer point; for loans, when the interest rate changes, the interest expense at this time point is the sum of the interest expense before refinancing, the penalty for breach of contract, and the interest expense after refinancing.(2)
Optimal Execution Strategy. A point on the path of interest rate change makes the interest income maximum or the interest expense minimum, then this point is the optimal execution strategy point, and the American option value is the option value of this point.The pricing of embedded option is to find the option value of the optimal execution strategy and discount the expected return at that time to the current time. Therefore, we can know that for a path of interest rate change, we can find the best execution strategy at each execution time from the back to the front. Through repeated deduction, we finally get the best execution strategy and the corresponding discount value of income on each path.


### 5.1.2. Application of Monte Carlo Simulation

(1) Monte Carlo Simulation Technology. The embedded option has the characteristics of multidimension and complexity, so the pricing method of numerical analysis is suitable. The numerical analysis model can effectively simulate various variables in reality, and this method has become an indispensable tool for embedded option valuation. The implementation of numerical analysis includes network analysis and Monte Carlo simulation. Network analysis is not suitable for pricing embedded options. Firstly, many fundamental characteristics of embedded options cannot be described by network methods, such as the jump-diffusion model describing the trend of interest rate change, which cannot be simulated by network methods. Secondly, for some high-dimensional weighted bonds with more underlying variables, it is impossible to price them. Thirdly, the discrete form of low dimensional equivalency is its fatal defect, that is to say, its default interest rate will only change once in the above equal time interval, while the actual situation is that the embedded option can be executed at any time before the maturity date of the option.When the dimension of financial derivative securities is relatively high, Monte Carlo simulation technology will become the only effective numerical analysis method; that is, it can effectively solve the “dimension disaster” problem. The essence of Monte Carlo simulation is a numerical analysis method to obtain the actual estimation of the mathematical expectation of random distribution by approximately replacing the overall expectation of random distribution with the sample mean of random sampling. Sometimes, when the pricing problem of complex options cannot be solved by analytical or numerical methods, Monte Carlo simulation can obtain the interest rate path and the final profit and loss of options by inputting assumed parameters and thousands of computer simulations, then obtain the distribution of option value, and finally use the average result to calculate the price of complex options [38–40].(2) Application of Monte Carlo Simulation Technology in Option Pricing. According to the nonarbitrage asset pricing theory, the price of financial derivative securities is the discounted expected value of its profit and loss function, and the expected value is the expected value under the risk neutral probability measure. According to the theory of probability statistics and stochastic analysis, these expected values can also be expressed in the form of integral in a specific region. Monte Carlo simulation is to estimate the price of financial derivative securities by simulating these mathematical expectations or integrals. The steps are as follows:SupposeSt is the price of the underlying asset at time t, r is the risk-free market interest rate, σ is the volatility of the underlying asset return, T is the maturity of the option, and X is the exercise price of the option. Under the background of Black–Scholes, the underlying asset price St satisfies the geometric Brownian motion:(30)dSt=rStdt+σStdWt.In equation (30), Wt is the geometric Brownian motion under the risk neutral probability measure, then the return function at the maturity date of the option can be expressed as follows:(31)C=e−rTES¯−X+.S¯ in equation (31) is the average price of the underlying asset in 0,T time.The specific process of Monte Carlo simulation to realize option value calculation is as follows:(1)
Simulaten paths of underlying asset prices and then get the average prices of each path S1¯,S2¯,…,Sn¯;(2)
LetY=ES¯−X+, calculate the period weight of the average price on each path and get the value Y1,Y2,…,Yn;(3)
LetCj=e−rTYj,n=1,2,…,n and average the value of C to get C¯; that is, C¯=1/n∑i=1nCj; this is the option value obtained by simulation.


## 5.1.1. Analysis of the Optimal Execution Strategy of Embedded Options

Before analyzing and pricing the embedded option, it is regarded as a standard American option, and the embedded option has the characteristics of early execution. The standard of early execution is that the value obtained by executing the option is greater than or equal to the value obtained by not executing in advance. Because the embedded option has a strong path dependence, it will have a different value because of the different interest rate changes, so there may be many opportunities for early execution in the duration of the option, but only one time is optimal, making the value of the option maximum, and the value at this time is the real value of the option. Therefore, the steps to find the optimal execution strategy are as follows:(1)Based on the change of interest rate, calculate the corresponding interest income or interest expense at each time point. For deposits, when the interest rate changes, the interest income at this time point is the sum of the demand interest before the deposit transfer point and the time deposit income after the deposit transfer point; for loans, when the interest rate changes, the interest expense at this time point is the sum of the interest expense before refinancing, the penalty for breach of contract, and the interest expense after refinancing.(2)
Optimal Execution Strategy. A point on the path of interest rate change makes the interest income maximum or the interest expense minimum, then this point is the optimal execution strategy point, and the American option value is the option value of this point.The pricing of embedded option is to find the option value of the optimal execution strategy and discount the expected return at that time to the current time. Therefore, we can know that for a path of interest rate change, we can find the best execution strategy at each execution time from the back to the front. Through repeated deduction, we finally get the best execution strategy and the corresponding discount value of income on each path.


## 5.1.2. Application of Monte Carlo Simulation

(1) Monte Carlo Simulation Technology. The embedded option has the characteristics of multidimension and complexity, so the pricing method of numerical analysis is suitable. The numerical analysis model can effectively simulate various variables in reality, and this method has become an indispensable tool for embedded option valuation. The implementation of numerical analysis includes network analysis and Monte Carlo simulation. Network analysis is not suitable for pricing embedded options. Firstly, many fundamental characteristics of embedded options cannot be described by network methods, such as the jump-diffusion model describing the trend of interest rate change, which cannot be simulated by network methods. Secondly, for some high-dimensional weighted bonds with more underlying variables, it is impossible to price them. Thirdly, the discrete form of low dimensional equivalency is its fatal defect, that is to say, its default interest rate will only change once in the above equal time interval, while the actual situation is that the embedded option can be executed at any time before the maturity date of the option.When the dimension of financial derivative securities is relatively high, Monte Carlo simulation technology will become the only effective numerical analysis method; that is, it can effectively solve the “dimension disaster” problem. The essence of Monte Carlo simulation is a numerical analysis method to obtain the actual estimation of the mathematical expectation of random distribution by approximately replacing the overall expectation of random distribution with the sample mean of random sampling. Sometimes, when the pricing problem of complex options cannot be solved by analytical or numerical methods, Monte Carlo simulation can obtain the interest rate path and the final profit and loss of options by inputting assumed parameters and thousands of computer simulations, then obtain the distribution of option value, and finally use the average result to calculate the price of complex options [38–40].(2) Application of Monte Carlo Simulation Technology in Option Pricing. According to the nonarbitrage asset pricing theory, the price of financial derivative securities is the discounted expected value of its profit and loss function, and the expected value is the expected value under the risk neutral probability measure. According to the theory of probability statistics and stochastic analysis, these expected values can also be expressed in the form of integral in a specific region. Monte Carlo simulation is to estimate the price of financial derivative securities by simulating these mathematical expectations or integrals. The steps are as follows:SupposeSt is the price of the underlying asset at time t, r is the risk-free market interest rate, σ is the volatility of the underlying asset return, T is the maturity of the option, and X is the exercise price of the option. Under the background of Black–Scholes, the underlying asset price St satisfies the geometric Brownian motion:(30)dSt=rStdt+σStdWt.In equation (30), Wt is the geometric Brownian motion under the risk neutral probability measure, then the return function at the maturity date of the option can be expressed as follows:(31)C=e−rTES¯−X+.S¯ in equation (31) is the average price of the underlying asset in 0,T time.The specific process of Monte Carlo simulation to realize option value calculation is as follows:(1)
Simulaten paths of underlying asset prices and then get the average prices of each path S1¯,S2¯,…,Sn¯;(2)
LetY=ES¯−X+, calculate the period weight of the average price on each path and get the value Y1,Y2,…,Yn;(3)
LetCj=e−rTYj,n=1,2,…,n and average the value of C to get C¯; that is, C¯=1/n∑i=1nCj; this is the option value obtained by simulation.


## 5.2. Pricing Steps of Embedded Options for Commercial Banks’ Deposits and Loans

### 5.2.1. Pricing Steps of Embedded Options in Deposits

Using Monte Carlo simulation to price the embedded options in commercial bank deposits, the steps are as follows:(1) Simulate the Sample Path of the Underlying Variable within the Validity Period under the Risk Neutral Measure. The underlying variable here is deposit interest rate, and the interest rate model describing the change of deposit interest rate is the jump-diffusion model: (32)drt=αβ−rtdt+σtr1/2dWt+JdP.Equation (32) is discretized and simulated. Take Δt=1 month to simulate the possible path of interest rate after n years (12∗n months) in the current initial interest rate state. By using this simulation method repeatedly, thousands of interest rate paths can be simulated to meet the needs of Monte Carlo simulation.(2) Calculation of Embedded Option Value on Simulation Path of Each Sample. Each interest rate sample path may involve the early withdrawal of deposits, so there may be multiple options on one interest rate path. However, based on the optimal strategy problem of American option execution, there is only one optimal execution point, and then the discount of the option value of the optimal execution point is the value of the embedded option on the interest rate path: (33)C=exp−rtmaxexph0,T⋅t⋅expRt,T⋅T−t−expR0,T⋅T,0.In equation (33), C is the value of the option, r is the risk-free interest rate, R0,T is the interest rate of the previous contract, h0,T is the current interest rate, Rt,T is the interest rate after the change of time t, and T is the deposit term.(3) Average the Corresponding Option Values on the Simulated Sample Path. The average value of the corresponding options on the simulated sample path is obtained, and the estimated value of the embedded options in the time deposit is obtained:(34)C¯=1n∑i=1nCi.


### 5.2.2. Pricing Steps of Embedded Options in Loans

The pricing process of embedded options in loans is similar to the calculation process of embedded option value in time deposits; except that the penalty for prepayment is to pay liquidated damages, that is, to charge a certain proportion of liquidated damages from the prepayment loans, and the proportion coefficient is expressed byα. From 2019, CCB’s prepayment policy for housing loans is to charge 3% of prepayment for less than one year, 2% for one to two years, and 1% for more than two years.Although the calculation of embedded options in loans is one more penaltyα than that in deposits, the Monte Carlo simulation method can still be used, and the calculation steps are as follows:(1) Simulate the Change Path of Loan Interest Rate. The 10-year loan interest rate model is shown in formula (26).(2) Calculation of Embedded Option Value on Simulation Path of Each Sample. Due to the existence of penalty α, the calculation of embedded options in loans is slightly different from that in deposits: (35)C=exp−rtmaxexpr0,T⋅T−α+expr0,t⋅t⋅exprt,T⋅T−t,0.In equation (35), C is the embedded option value, r is the risk-free interest rate, α is the penalty when prepaying the loan, t is the time when the option is exercised, T is the loan term, r0,T is the interest rate agreed in the previous loan contract, and rt,T is the interest rate in the new loan contract signed at t.(3) Average the Corresponding Option Values on the Simulated Sample Path. By averaging the corresponding options on the simulated sample path, we can get the estimated value of the embedded options in loans:(36)C¯=1n∑i=1nCi.


## 5.2.1. Pricing Steps of Embedded Options in Deposits

Using Monte Carlo simulation to price the embedded options in commercial bank deposits, the steps are as follows:(1) Simulate the Sample Path of the Underlying Variable within the Validity Period under the Risk Neutral Measure. The underlying variable here is deposit interest rate, and the interest rate model describing the change of deposit interest rate is the jump-diffusion model: (32)drt=αβ−rtdt+σtr1/2dWt+JdP.Equation (32) is discretized and simulated. Take Δt=1 month to simulate the possible path of interest rate after n years (12∗n months) in the current initial interest rate state. By using this simulation method repeatedly, thousands of interest rate paths can be simulated to meet the needs of Monte Carlo simulation.(2) Calculation of Embedded Option Value on Simulation Path of Each Sample. Each interest rate sample path may involve the early withdrawal of deposits, so there may be multiple options on one interest rate path. However, based on the optimal strategy problem of American option execution, there is only one optimal execution point, and then the discount of the option value of the optimal execution point is the value of the embedded option on the interest rate path: (33)C=exp−rtmaxexph0,T⋅t⋅expRt,T⋅T−t−expR0,T⋅T,0.In equation (33), C is the value of the option, r is the risk-free interest rate, R0,T is the interest rate of the previous contract, h0,T is the current interest rate, Rt,T is the interest rate after the change of time t, and T is the deposit term.(3) Average the Corresponding Option Values on the Simulated Sample Path. The average value of the corresponding options on the simulated sample path is obtained, and the estimated value of the embedded options in the time deposit is obtained:(34)C¯=1n∑i=1nCi.


## 5.2.2. Pricing Steps of Embedded Options in Loans

The pricing process of embedded options in loans is similar to the calculation process of embedded option value in time deposits; except that the penalty for prepayment is to pay liquidated damages, that is, to charge a certain proportion of liquidated damages from the prepayment loans, and the proportion coefficient is expressed byα. From 2019, CCB’s prepayment policy for housing loans is to charge 3% of prepayment for less than one year, 2% for one to two years, and 1% for more than two years.Although the calculation of embedded options in loans is one more penaltyα than that in deposits, the Monte Carlo simulation method can still be used, and the calculation steps are as follows:(1) Simulate the Change Path of Loan Interest Rate. The 10-year loan interest rate model is shown in formula (26).(2) Calculation of Embedded Option Value on Simulation Path of Each Sample. Due to the existence of penalty α, the calculation of embedded options in loans is slightly different from that in deposits: (35)C=exp−rtmaxexpr0,T⋅T−α+expr0,t⋅t⋅exprt,T⋅T−t,0.In equation (35), C is the embedded option value, r is the risk-free interest rate, α is the penalty when prepaying the loan, t is the time when the option is exercised, T is the loan term, r0,T is the interest rate agreed in the previous loan contract, and rt,T is the interest rate in the new loan contract signed at t.(3) Average the Corresponding Option Values on the Simulated Sample Path. By averaging the corresponding options on the simulated sample path, we can get the estimated value of the embedded options in loans:(36)C¯=1n∑i=1nCi.


## 5.3. Pricing of Embedded Options Based on Monte Carlo Simulation

### 5.3.1. An Empirical Study of Embedded Option Pricing in Deposits

We take the latest five-year fixed deposit interest rate of 0.0275 on January 1, 2021, as the initial interest rate and the one-year fixed deposit interest rate of 0.015 as the discount rate, and the deposit amount is standardized to 1, and the embedded option value is calculated in the five-year fixed deposit through simulation.We have analyzed the selection of the interest rate model and the principle and steps of the Monte Carlo simulation method. In order to make the simulation results more accurate, we need to use MATLAB software. We know from the characteristics of Monte Carlo simulation that the accuracy of the simulation results has a great relationship with the simulation times. The more the simulation times are, the more accurate the results will be. Therefore, 100000 times of simulation are carried out here, and the simulation results are shown in Figure9.It can be clearly seen from the convergence chart of option value in Figure9 that the value of embedded option in 5-year fixed deposit of commercial banks is 0.0012. In other words, the real deposit interest rate of commercial banks is 0.0012 more than the nominal deposit interest rate.Figure 9
Simulation results of embedded option value in 5-year fixed deposit.


### 5.3.2. Empirical Analysis of Embedded Option Pricing in Loans

We take the latest 10-year loan interest rate of 0.049 on January 1, 2021, as the initial interest rate and the one-year fixed deposit interest rate of 0.015 as the discount rate, and the loan amount is standardized to 1. Assuming that the customer needs to pay 1% penalty for prepayment, we calculate the embedded option value in the 10-year loan through simulation and simulate 100000 knots. The results are as follows.It can be clearly seen from the convergence chart of option value in Figure10 that the value of embedded option in 10-year loan of commercial bank is 0.0028. In other words, the 10-year real loan interest rate of commercial banks is 0.0028 less than the nominal loan interest rate.Figure 10
Simulation results of embedded option value in 10-year loan.


## 5.3.1. An Empirical Study of Embedded Option Pricing in Deposits

We take the latest five-year fixed deposit interest rate of 0.0275 on January 1, 2021, as the initial interest rate and the one-year fixed deposit interest rate of 0.015 as the discount rate, and the deposit amount is standardized to 1, and the embedded option value is calculated in the five-year fixed deposit through simulation.We have analyzed the selection of the interest rate model and the principle and steps of the Monte Carlo simulation method. In order to make the simulation results more accurate, we need to use MATLAB software. We know from the characteristics of Monte Carlo simulation that the accuracy of the simulation results has a great relationship with the simulation times. The more the simulation times are, the more accurate the results will be. Therefore, 100000 times of simulation are carried out here, and the simulation results are shown in Figure9.It can be clearly seen from the convergence chart of option value in Figure9 that the value of embedded option in 5-year fixed deposit of commercial banks is 0.0012. In other words, the real deposit interest rate of commercial banks is 0.0012 more than the nominal deposit interest rate.Figure 9
Simulation results of embedded option value in 5-year fixed deposit.


## 5.3.2. Empirical Analysis of Embedded Option Pricing in Loans

We take the latest 10-year loan interest rate of 0.049 on January 1, 2021, as the initial interest rate and the one-year fixed deposit interest rate of 0.015 as the discount rate, and the loan amount is standardized to 1. Assuming that the customer needs to pay 1% penalty for prepayment, we calculate the embedded option value in the 10-year loan through simulation and simulate 100000 knots. The results are as follows.It can be clearly seen from the convergence chart of option value in Figure10 that the value of embedded option in 10-year loan of commercial bank is 0.0028. In other words, the 10-year real loan interest rate of commercial banks is 0.0028 less than the nominal loan interest rate.Figure 10
Simulation results of embedded option value in 10-year loan.


## 6. Conclusion

For depositors, time deposits can be divided into a fixed claim and several American interest rate call options. The necessary condition for the execution of the embedded option in deposit is that the interest rate rises, and the more the interest rate rises, the greater the possibility of the option execution. In addition, the possibility of option execution decreases with the passage of time, generally only in the early implementation. For the borrower, the loan can be decomposed into a basic debt and several American interest rate put options. The necessary condition for the execution of embedded option in loan is that the interest rate falls, and the more the interest rate falls, the more likely the option is to be executed. The greater the penalty, the less likely the option is to be executed. The largert is, the less likely the option is to be executed, which is generally only executed in the early stage.Interest rate is the underlying asset of embedded option; based on the in-depth study of the characteristics and movement law of interest rate, this paper obtains an interest rate model suitable for Chinese market characteristics; that is, adding jump factor to the CIR model to describe the jump of interest rate caused by government macrocontrol or market mutation, which is the jump-diffusion model. Through the Monte Carlo simulation test, the jump-diffusion model has smaller error value than the CIR model, which is more suitable to describe China’s interest rate market.Based on the jump-diffusion model of interest rate, this paper uses the Monte Carlo simulation method to analyze and price the embedded options in 5-year fixed deposit and 10-year loan. The results show that the embedded options of 5-year fixed deposit and 10-year loan are 0.0012 and 0.0028, respectively. Bank’s real deposit and loan spread = (loan interest rate – American interest rate put option price) − (deposit interest rate + American interest rate call option price). Therefore, we can see that the existence of embedded options significantly reduces the bank’s real deposit and loan spread.



---

*Source: 9975536-2021-05-17.xml*