We study the pricing of American options in an incomplete market in which the dynamics of the underlying risky asset is driven by a jump diffusion process with stochastic volatility. By employing a risk-minimization criterion, we obtain the Radon-Nikodym derivative for the minimal martingale measure and consequently a linear complementarity problem (LCP) for American option price. An iterative method is then established to solve the LCP problem for American put option price. Our numerical results show that the model and numerical scheme are robust in capturing the feature of incomplete finance market, particularly the influence of market volatility on the price of American options.

1. Introduction

Since the pioneering work of Black-Scholes [1] and Merton [2] on the pricing of options, option pricing has become one of the predominant concerns in financial market. By relaxing some of the restrictive assumptions, such as the assumption of constant volatility and the geometric Brownian motion for the price of the underlying asset, many option pricing models have been proposed. Today, an increasing amount of literature has focused on stochastic volatility models (see, e.g., [3]), jump diffusion or Lévy process models [4, 5], which fit some features of the observed empirical market data. In this paper, in constructing our price model, we will take account of both the feature of stochastic volatility as in Heston's model and the feature of jump diffusion as in Bates [6]. The combination of the stochastic volatility and jump diffusion in the same model enables the model to capture the nature of the real financial market and yields a challenging and difficult mathematical problem for analysis.

Because of infrequent trading or by the presence of additional sources of uncertainty, the completeness of the market often breaks down. Thus, replicating the payoff of a contingent claim by a portfolio is possible, and various equivalent martingale measures (EMM) have been adopted. The question is that, among the equivalent martingale measures, one needs to establish certain criterions to determine the appropriate measure to be adopted. A unique equivalent martingale measure by using utility maximization has been identified by many authors including Kallsen [7], Bielecki and Jeanblanc [8], and Cvitanić et al. [9]. It was also found that, under the minimal martingale measure, there exists a unique risk-minimizing strategy hedging of contingent claims in complete market [10]. The criterion under the minimal martingale measure is thus referred to as risk-minimization criterion (see [11]). According to this, it is possible to price option with risk-minimization criterion. In this paper, we use the minimal martingale measure and consider the switch of assets prices with jump diffusion and stochastic volatility and obtain the Radon-Nikodym derivative and a linear complementarity problem for the pricing of American option.

Unlike European options, American options, which can be exercised any time before the expiring date, have to be priced usually numerically. Comparing to European options which lead to a partial integrodifferential equation, American options lead to inequalities in the form of a linear complementarity problem (LCP). In the American option pricing context, several methods have been proposed to approximate the linear complementarity problems resulting from American option pricing. These include the finite difference method, the operator splitting method, the multigrid method, and the penalty method [12]. Moreover, Tavella and Randall [13] described a stationary iterative method for pricing European options, and then Salmi and Toivanen [14] proposed a generalization of this iterative method to price American option. But they did not consider any multifactor jump diffusion models such as stochastic volatility models with jumps. Our paper aims to further develop previous numerical method to deal with multifactor jump diffusion models via the employment of the risk-minimization criterion to obtain the Radon-Nikodym derivative for the minimal martingale measure.

The rest of the paper is organized as follows. In Section 2, we discuss the pricing model with stochastic volatility and jump diffusion under the minimal martingale measure and then derive the linear complementarity problem under the risk-minimization criterion. Section 3 presents the discrete space and time approximation for the LCP problem. The method of projected successive overrelaxation (projected SOR) is used to yield our results. Numerical results and discussion are reported in Section 4.

2. The Model under Minimal Martingale Measure

Taking into account the jumps in the evolution process of the underlying asset, we introduce the following stochastic volatility (1) based on the work in Heston [3]: where is defined on the filtered complete space with 2-dimensional standard Brownian motion with correlation coefficient . The asset price has drift and is the compensated jump measure, denotes the jump measure, and is the Lévy measure, governs the volatility of the variance process, and and are the mean reversion coefficients to ensure that the variance will drift back to at the rate of .

According to the discussions in the reference by Föllmer and Schweizer [10], the optimal hedging strategy may be established with the minimal martingale measure and it has been proved that the minimal martingale measure can be determined uniquely. Using Theorem 3.1 in Heston [3], the Radon-Nikodym derivative under the minimal martingale measure is where and is the risk-free rate of interest in the discounted risky asset price process . The derivation of (2) is given in Appendix A. Correspondingly, the Brownian motion and the compensated Poisson process under the minimal martingale measure are derived as follows:

In our model, there are two dimensions of risks from the Lévy process and the Brownian motion. Thus the financial market is not complete and such a claim has an intrinsic risk. In this context, we construct the pricing model under the risk-minimization criterion by employing the minimal martingale measure method, and therefore model (1) can be written as follows:

The value of the American put option at time , under strike price and maturity date , is where , representing the stopping times in , is to be determined to maximize the discounted option payoff.

Under the minimal martingale measure , we obtain the following partial integrodifferential equation (see Appendix B for the detailed derivation): where

Under Merton’s jump diffusion model, we assume , where is the jump intensity and is distribution of the jump size , and assume that the stochastic volatility is the log-normal distribution with the density Let ; thus the notations under the log-normal density functions above are ,, and .

Now we take the time transform and introduce an operator by

Then, the value of the American put option satisfies the following linear complementarity problem: with the boundary conditions and the initial condition

3. Numerical Discretization

We use the uniform grid , , and on the domain . Denoting and using the central difference scheme, we have

Through a long derivation, we obtain from (10) the following set of discrete equations and inequalities in matrix form: where ,, in which is a block tridiagonal matrix due to the first six terms in (10), and is a full matrix contributed from the two integral terms of (10). The matrix can be expressed in the form in which ,  , , and .

To calculate the matrix due to the two integral terms in (10), let and rewrite the integrals by

The above integrals can then be discretized by using linear interpolation, namely, where

Let , and . Then (19) becomes where in which ,  , , and ,. Hence .

In a similar way, we obtain , where in which ,  ,  , , and .

Hence, the integral part of the operator is .

Now for American options we obtain a semidiscrete LCP. The Crank-Nicolson method [15] is then used for time discretization. Let ; we obtain which can be denoted by where and .

The projected SOR algorithm PSOR() can then be used to solve the LCP problem. The algorithm is as follows.

For , , where the relaxation parameter ), wherein .

4. Numerical Results and Discussion

In this section, we present a set of results to demonstrate the feature of the model and the robustness of the model for the determination of American put option price and then investigate the influence of market volatility. The basic parameter values used in the computation are given in Table 1. The numerical algorithm is implemented in computer by using the MATLAB 2012 software package. Figure 1 shows the payoff as the function of stock price. The numerical solutions obtained by solving the LCP problem (24) are showed in Figure 2. We also compute the value of American option price under various values of volatility. The results are plotted in Figure 3.

From the numerical results, we can see that volatility has no influence on the payoff of American option, but it plays an important role in American option pricing. From Figures 2 and 3, we can see that, when the stock price is relatively low, the influence of volatility is not significant, while, as the stock price increases, the influence of volatility becomes more and more significant. It can also be found that, at the same level of stock price, the larger the volatility is, the more the influence it plays on American option price is, which is in agreement with the market reality. Thus, the results show that the model and numerical scheme are robust in capturing the feature of incomplete finance market, particularly the influence of market volatility on the price of American options.


A. Derivation of the Radon-Nikodym Derivative under the Minimal Martingale Measure

This appendix derives the Radon-Nikodym derivative under the minimal martingale measure and then consequently obtains (5) for American option pricing model.

According to the Doob-Meyer decomposition, the discounted risky asset price process, , is a special semimartingale and can be written as with where is the martingale part of and is the predictable process of finite variation.

The minimal martingale measure is unique [10] if there exists a predictable process that satisfies

By using (A.4), we get the theorem below for calculating the Radon-Nikodym derivative.

Theorem 1. The Radon-Nikodym derivative under the minimal martingale measure is

Proof. According to the Girsanov transformation theory, the predictable process of bounded variation can be computed in terms of , where is the quadratic variation process between and . Under the minimal martingale measure , the predicable process of bounded variation in the Doob-Meyer decomposition of is given by
Using (A.4) and (A.7), we have
From (A.2), we get Let ; then (A.8) can be written as . Hence, Because is the Doléans-Dade exponential, we obtain where Solving equations (A.11), we obtain as given in Theorem 2.

B. Derivation of the Partial Integrodifferential Equation (7)

Here we derive the partial integrodifferential equation (PIDE) (7), which is used in Section 2 to construct the linear complementarity problem (LCP) for the pricing of American put option.

Under the measure of minimal martingale , the price of American put option , having strike price and maturity date , is as follows: where , representing the stopping times in , is to be determined to maximize the discounted option payoff.

To derive the LCP problem for the American put option, we first need to derive the following theorem for fixed , which is for European put option.

Theorem 2. The price of the European option satisfies the following PIDE:

Proof. Using the Itô formula of semimartingales, we get By substituting (5) into the equation above, we obtain By combining the like terms of , we get
Because the discounted price of the European style option is a martingale, we set the drift term to zero, and consequently we obtain the equation as given in Theorem 2. By rearranging and simplifying (B.2), we obtain (7) in Section 2. Consequently, based on the operator derived from this theorem, the LCP problem (11) is constructed in Section 2 for solving the American option model.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.