Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 165727, 9 pages

http://dx.doi.org/10.1155/2013/165727

Research Article

Option Pricing under Risk-Minimization Criterion in an Incomplete Market with the Finite Difference Method

Xinfeng Ruan,¹ Wenli Zhu,¹ Shuang Li,² and Jiexiang Huang¹

¹School of Economic Mathematics, Southwestern University of Finance and Economics, Chengdu 611130, China

²Department of Mathematics and Statistics, Curtin University, Perth, WA 6102, Australia

Received 16 January 2013; Revised 22 March 2013; Accepted 13 April 2013

Academic Editor: Guangchen Wang

Copyright © 2013 Xinfeng Ruan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We study option pricing with risk-minimization criterion in an incomplete market where the dynamics of the risky underlying asset is governed by a jump diffusion equation with stochastic volatility. We obtain the Radon-Nikodym derivative for the minimal martingale measure and a partial integro-differential equation (PIDE) of European option. The finite difference method is employed to compute the European option valuation of PIDE.

1. Introduction

Option pricing problem is one of the predominant concerns in the financial market. Since the advent of the justly celebrated Black-Scholes option pricing formula in [1], there has been an increasing amount of the literature describing the theory and its practice. In the Black-Scholes model, the appreciation rate and the volatility rate are assumed to be constants. However, more and more empirical evidence has revealed that these assumptions are not consistent with reality. Hence, many option valuation models which relax some of the restrictive assumptions in the Black-Scholes framework have been proposed and tested, such as the stochastic volatility models in [2–8], the jump diffusion or the Lévy process models in [9–12]. The model to be studied in this paper takes into account not only stochastic volatility based on Heston’s model in [8], but also the jump diffusion case. Therefore, several new problems and difficulties are exposed in our model.

Different from the Black-Scholes framework which uses jump diffusion to describe the price dynamics of underlying asset, the market of our model is incomplete; that is, it is not possible to replicate the payoff of every contingent claim by a portfolio, and there are several equivalent martingale measures. How to choose a consistent pricing measure from the set of equivalent martingale measures becomes an important problem. This means that we need to find some criterions to determine one from the set of equivalent martingale measures in some economically or mathematically motivated fashion. Föllmer and Leukert (2000), Kallsen (1999), Cvitanić et al. (2001), and Bielecki and Jeanblanc (2008) in [13–16] identified a unique equivalent martingale measure by utility maximization. Then, the option valuation under the minimal martingale measure was further developed by several researchers. Schweizer (1991) and Föllmer and Schweizer (1991) in [17, 18] found that under the minimal martingale measure, a unique risk-minimizing (or optimal) strategy hedging of contingent claims in incomplete market exists. In our paper, we name the criterion under the minimal martingale measure as a risk-minimization criterion. Thus, in the incomplete market, option pricing is approximately possible with risk-minimization criterion. As presented in this paper, our work is based on the task of Föllmer and Schweizer in [17, 18], and the purpose of this paper is to find the minimal martingale measure and the measure switch of asset prices processes with stochastic volatility and jump diffusion. By employing the minimal martingale measure, we obtain the Radon-Nikodym derivative and a partial integro-differential equation (PIDE) for the European option.

However, it is difficult to get the exact solution of the PIDE in our model. Several numerical methods have been proposed to solve the PIDE approximately. The methods include numerical integration by Chiarella and Ziogas (2009) in [19], finite elements by Matache et al. (2005) in [20], the method of lines by Meyer (1998) in [21], and the finite difference methods (FDM) including those by Carr and Hirsa (2003) in [22], d’Halluin et al. (2004) in [23], and Briani et al. (2004) in [24]. This paper employs the FDM to compute the valuation of European option, and American option. Then further studies on using other methods to compute the solution of PIDE and derive the pricing formula of another option will be included in our future study.

The rest of the paper is organized as follows. In Section 2, we present the model for the underlying market and Doob-Meyer decomposition of the risky asset. In Section 3, we investigate an explicit representation of the density process of the minimal martingale measure. In Section 4, we derive a PIDE for the European option. The European option pricing with FDM is then studied in Section 5. Numerical results are shown in Section 6, and conclusions are given in Section 7.

2. The Model

In this paper, we consider the financial market with the following two basic assets:

(i)a Bond whose price at time is given by

(ii)a Stock whose price at time is given by

Here, ,, and ; on the filtered complete space , there are , which both are 1-dimensional Brownian motions with and which is the compensated jump measure of . is the jump measure, and is the Lévy measure of the Lévy process . is given by , . Additionally, .

From above settings, the risk is at least two dimensional including the Brownian motion and Lévy process but only one risky asset in the market. It is impossible to exactly replicate the payoff of a given option by a dynamic portfolio strategy of the basic assets which is self-financing. We all know that if the dynamic portfolio strategy exists, then the initial cost of the portfolio must equal the price of the option, and the market is called complete; otherwise, an arbitrage opportunity exists, and the market is called incomplete. Therefore, the market of our model is incomplete since the perfect replication is impossible.

Compared with the Black-Scholes model, our model remedies some serious drawbacks of the Black-Scholes model, such as the constant volatility assumption and market completeness assumption. Bakshi et al. in [10] systematically analyzed the performance of the stochastic volatility model, the jump diffusion model, and the stochastic interest rate model, and concluded that a model with stochastic volatility and the jump diffusion was a better alternative to the Black-Scholes model, because the former not only performed far better but also was practically implementable. An excellent innovation of our model is combining stochastic volatility model and the jump diffusion model. Consequently, our model is better than the Black-Scholes model to portray realistic financial markets in theory and practice.

Now, we want to find a unique optimal strategy hedging of contingent claims under risk-minimization criterion. Based on [18], it is equivalent to find the minimal martingale measure from the set of equivalent martingale measures and then obtain approximate prices of contingent claims.

With 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.

3. Minimal Martingale Measure

We introduce the nations of minimal martingale measure in this section. Föllmer and Schweizer (1991) [18] noticed that the optimal hedging strategy can be computed in terms of the minimal martingale measure. Furthermore, it is uniquely determined. Hence, under the minimal martingale measure, the Radon-Nikodym derivative can be found and computed. Before that, we define the minimal martingale measure.

Definition 1 (see [18]). A local martingale measure , equivalent to the original measure , is called minimal if on and if any square-integrable -martingale which is orthogonal to remains a local martingale under .

Theorem 2 (see [18]). (i) The minimal martingale measure is uniquely determined.

(ii) exists if and only if there exists a predictable process that satisfies

Using Theorem 2, we obtain the following theorem for computing the Radon-Nikodym derivative.

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

Proof. The theory of the Girsanov transformation shows that the predictable process of bounded variation can also be computed in terms of :

Throughout this paper, we make use of the notations that defined the quadratic variation process between and and denote . Under , the predictable process of bounded variation in the Doob-Meyer decomposition of is given by

Using (6) and (9), we have

Denote , then (10) can be written as

From (4), we get

Hence,

From (11), we know that is the Doléans-Dade exponential. Thus, we obtain

where

Solving (14), we obtain in Theorem 3.

Remark 4. The Brown motions under the minimal martingale measure are

and the compensatory of is

Remark 5. Equation (2) under the minimal martingale measure is written as

To guarantee that is a martingale under the minimal martingale measure , the following corollary is necessary.

Corollary 6. Under the minimal martingale measure , is a martingale if and only if

Proof. Substituting into (18), since is a martingale, the drift term must be identical to zero. Then, we can get (19).

4. Partial Integro-Differential Equation for European Call Option

Under the minimal martingale measure , the price of the European call option at time with strike price and maturity date is given by

and .

By the fact that the discounted price of the European call option is a martingale under , we can obtain the following theorem.

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

and .

Proof. The total derivative of the discounted option price is

We make the drift term zero, since the discounted price of the European put option is a martingale. Then, we obtain in Theorem 7.

5. Finite Difference Scheme of European Call Option

5.1. The Problem

In this section, we will derive the finite difference scheme for PIDE (21). Before that, we assume that , , and . In this paper, we take Merton’s jump diffusion model with the density which is from the log-normal distribution:

Furthermore, using error function, we denote that

where , , , , , , , and are constants, and .

From (21) and (19), we get the PIDE with original parameters in the original model (2) as follows:

To represent conveniently, we denote that

Then, (25) can be written as

Now, we denote , , , and . Then we have the following problem equivalent to (27).

Problem 1. Consider

5.2. Discretization

To solve (38), we use a finite difference scheme with the following nation and approximations:

Now, we evaluate the integral term. First, by performing a change of variable , and , the and can be discretized by using the linear interpolation. Thus,

Now, by using linear interpolation, we get approximation the following:

where

Denote the following:

Then, we have

Denote the following:

Then, we have , where is the error function defined by with and .

Then,

Finally, discretization of (27) gives

It equals

where

6. Numerical Results

In this section, we adopted the basic set of appropriate parameter values of the model shown in Table 1 as our basic set of parameter values to solve the finite difference equation (38).

Table 1: Basic set of appropriate parameter values of the model.

0.23

3.46

(0.0894)²

0.0001

−0.82

1

200

1

0.039

0.14

0.47

−0.086

100

2000

30

30

Figure 1 shows the payoff with the price and the volatility. The solution of (38) is shown in Figure 2. All numerical experiments have been implemented by MATLAB R2011b software on a 2.0-GHz Intel Core PC.

Figure 1: The payoff with the price and the volatility.

Figure 2: The solution of (38).

Now, we suppose that the underlying stock price process follows (2) and there is a European call option with , , and, other parameters in Table 1.  Table 2 shows the European option prices solved by FDM with different which stands for the different accuracy.

Table 2: European option prices.

Price

8.7858
10

9.3014
20

9.6929
30

9.7006
40

7. Conclusion

With risk-minimization criterion, we employ the minimal martingale measure to solve the pricing problem in an incomplete market. Then, we obtain the Radon-Nikodym derivative under the minimal martingale measure and a partial integro-differential equation (PIDE) for the European option. Since it is difficult to get the exact solution of PIDE in our model, a FDM scheme is proposed to compute the solution approximately. Finally, we complete a European call option price, and it is shown that our method is stable and locally accurate.

References

F. Black and M. Scholes, “The pricing of options and corporate liabilities,” Journal of Political Economy, vol. 81, no. 3, pp. 637–659, 1973.

J. Hull and A. White, “The pricing of options on assets with stochastic volatilities,” The Journal of Finance, vol. 42, no. 2, pp. 281–300, 1987.

L. Scott, “Option pricing when the variance changes randomly: theory, estimation, and an application,” The Journal of Financial and Quantitative Analysis, vol. 22, no. 4, pp. 419–438, 1987.

J. Wiggins, “Option values under stochastic volatility: theory and empirical estimates,” Journal of Financial Economics, vol. 19, no. 2, pp. 351–372, 1987.

S. Melino and S. M. Turnbull, “Pricing foreign currency options with stochastic volatility,” Journal of Econometrics, vol. 45, no. 1-2, pp. 239–265, 1990.

S. Melino and S. M. Turnbull, “Misspecification and the pricing and hedging of long-term foreign currency options,” Journal of International Money and Finance, vol. 14, no. 3, pp. 373–393, 1995.

E. M. Stein and J. C. Stein, “Stock price distributions with stochastic volatility: an analytic approach,” The Review of Financial Studies, vol. 4, no. 4, pp. 727–752, 1991.

S. Heston, “A closed-form solution for options with stochastic volatility with applications to bond and currency options,” The Review of Financial Studies, vol. 6, no. 2, pp. 327–343, 1993.

R. Merton, “Option pricing when underlying stock returns are discontinuous,” Journal of Financial Economics, vol. 3, no. 1-2, pp. 125–144, 1976.

G. Bakshi, C. Cao, and Z. Chen, “Empirical performance of alternative option pricing models,” The Journal of Finance, vol. 52, no. 5, pp. 2003–2049, 1997.

F. E. Benth and T. Meyer-Brandis, “The density process of the minimal entropy martingale measure in a stochastic volatility model with jumps,” Finance and Stochastics, vol. 9, no. 4, pp. 563–575, 2005.

D. Bates, “The crash of '87: was it expected? The evidence from options markets,” The Journal of Finance, vol. 46, no. 3, pp. 1009–1044, 1991.

H. Föllmer and P. Leukert, “Efficient hedging: cost versus shortfall risk,” Finance and Stochastics, vol. 4, no. 2, pp. 117–146, 2000.

J. Kallsen, “A utility maximization approach to hedging in incomplete markets,” Mathematical Methods of Operations Research, vol. 50, no. 2, pp. 321–338, 1999, Financial optimization

J. Cvitanić, W. Schachermayer, and H. Wang, “Utility maximization in incomplete markets with random endowment,” Finance and Stochastics, vol. 5, no. 2, pp. 259–272, 2001.

T. Bielecki and M. Jeanblanc, Indifference Prices in Indifference Pricing, Theory and Applications, Financial Engineering, Edited by R. Carmona, Princeton University Press, 2008.

M. Schweizer, “Option hedging for semimartingales,” Stochastic Processes and their Applications, vol. 37, no. 2, pp. 339–363, 1991.

H. F. Föllmer and M. Schweizer, “Hedging of contingent claims under incomplete information,” in Applied Stochastic Analysis, M. H. A. Davis and R. J. Elliott, Eds., vol. 5, pp. 389–414, Stochastics Monographs, London, UK, 1991.

C. Chiarella and A. Ziogas, “American call options under jump-diffusion processes: a fourier transform approach,” Applied Mathematical Finance, vol. 16, no. 1-2, pp. 37–79, 2009.

A. M. Matache, C. Schwab, and T. P. Wihler, “Fast numerical solution of parabolic integrodifferential equations with applications in finance,” SIAM Journal on Scientific Computing, vol. 27, no. 2, pp. 369–393, 2005.

G. H. Meyer, “The numerical valuation of options with underlying jumps,” Acta Mathematica, vol. 67, no. 1, pp. 69–82, 1998.

P. Carr and A. Hirsa, “Why be backward? forward equations for American options,” Risk, vol. 16, no. 1, pp. 103–107, 2003.

Y. d'Halluin, P. A. Forsyth, and G. Labahn, “A penalty method for American options with jump diffusion processes,” Numerische Mathematik, vol. 97, no. 2, pp. 321–352, 2004.

M. Briani, C. La Chioma, and R. Natalini, “Convergence of numerical schemes for viscosity solutions to integro-differential degenerate parabolic problems arising in financial theory,” Numerische Mathematik, vol. 98, no. 4, pp. 607–646, 2004.