Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 541023, 9 pages

http://dx.doi.org/10.1155/2014/541023

Research Article

Construction of Interval Shannon Wavelet and Its Application in Solving Nonlinear Black-Scholes Equation

Liwei Liu

Mechanical Engineering Department, North China Institute of Aerospace Engineering, 133 Aimin East Road, Langfang, Hebei 065000, China

Received 17 October 2013; Revised 27 December 2013; Accepted 29 December 2013; Published 20 February 2014

Academic Editor: Xu Zhang

Copyright © 2014 Liwei Liu. 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.

Interval wavelet numerical method for nonlinear PDEs can improve the calculation precision compared with the common wavelet. A new interval Shannon wavelet is constructed with the general variational principle. Compared with the existing interval wavelet, both the gradient and the smoothness near the boundary of the approximated function are taken into account. Using the new interval Shannon wavelet, a multiscale interpolation wavelet operator was constructed in this paper, which can transform the nonlinear partial differential equations into matrix differential equations; this can be solved by the coupling technique of the wavelet precise integration method (WPIM) and the variational iteration method (VIM). At last, the famous Black-Scholes model is taken as an example to test this new method. The numerical results show that this method can decrease the boundary effect greatly and improve the numerical precision in the whole definition domain compared with Yan’s method.

1. Introduction

The nonlinear Black-Scholes equations have been increasingly attracting interest over the last two decades, since they provide more accurate values by taking into account more realistic assumptions, such as the transaction costs, risks from an unprotected portfolio, large investor’s preferences, or illiquid markets, which may have an impact on the stock price, the volatility, the drift, and the option price itself [1]. In contrast with the traditional Black-Scholes equation, it is almost impossible to find the analytical solution of the nonlinear Black-Scholes equations. Most of the numerical schemes focus on the finite difference methods [2, 3]. As the matter of fact, the wavelet precise integration method (WPIM) is a simple and effective method for linear partial differential equations proposed by Mei et al. [4]. For linear steady structural dynamic systems, their numerical results at the integration points are almost equal to those of the exact solution in machine accuracy. Yan proposed a wavelet precise integration method for the nonlinear Black-Scholes model [5] recently by combining the variational iteration method (VIM) and homotopy perturbation method (HPM) [6].

Since the definition domain of wavelet transformation is an infinite interval, the boundary effect would occur when being applied for resolving the engineering problems with bounded interval, for example, ordinary differential equations (ODEs). Consequently, it will decrease the precision and computational efficiency of the solution. Nevertheless, the boundary effect can be eliminated effectively by constructing an interval wavelet using numerical methods. There are several ways available to construct an interval wavelet. Mei et al. [7] proposed a general construction method of interval wavelet based on the restricted variational principle. According to this method, there exist some interval wavelets that have been constructed to solve PDEs in engineering [8–11].

The main purpose of this paper is to deduce a novel interval wavelet based on the restricted variational principle and construct an interval wavelet precise integration method (WPIM) for the nonlinear Black-Scholes model. According to WPIM, the nonlinear differential equation should be transformed to a system of ordinary differential equations with the multiscales wavelet interpolation operator [12], and then the nonlinear PDEs become a system of nonlinear differential equations. The matrix differential equation (MDE) can be solved almost exactly by the coupling technique of the precise integration method and the variational iteration method (VIM) [13–17]. In contrast to the traditional finite difference approximation, the numerical results obtained by the PIM are close to the computer precision and also are free from the stiff problem in solving a set of simultaneous linear time-invariant ODEs.

2. Construction of Interval Shannon Wavelet Based on Restricted Variational Principle

Mei et al. [7] proposed a kind of construction method of the interval interpolation wavelet based on the restricted variational principle. It does not matter with the representation of the wavelet function. In this section, a new interval wavelet function with this method will be constructed. It takes into account the gradient and the smoothness of the approximated function. This is helpful to improve the calculation precision of the wavelet numerical method.

The representation of Shannon wavelet [18–21] is based upon approximating the Dirac delta function as a band-limited function and is given by

The Shannon wavelet possesses many excellent numerical properties such as interpolating, relative sparse, and orthogonal properties. A perceived disadvantage of (1) is that it tends to zero quite slowly as . A direct consequence of this is that there are a large number of grid points will contribute to the derivatives calculation of approximated function. It is for this reason that Hoffman et al. [22] have suggested using the Shannon-Gabor wavelet as follows:

It has been proved that (2) can improve the localized and asymptotic behavior of the Shannon scaling function. is the width parameter (or called window size). However, the presence of the Gaussian window destroys the orthogonal properties possessed by the Shannon wavelet, effectively worsening the approximation to a Dirac delta function. So, the Shannon wavelet representation of the Dirac delta function is adopted in this paper, and it is shown that this representation ensures that the approach is identical to the weighted residual approach. The corresponding discrete formula of the Shannon scaling function is

In the following, the construction of the interval Shannon wavelet is based on the general variational principle [7]. Consider one-dimensional problem , whose domain of definition is . Assuming that the number of discrete points is , , the th discrete point of variable could be written as

If the Shannon wavelet with interpolation property is employed as the trial function , that is,

any element in the matrices and could be expressed, respectively, as

where is the Lagrange basic function.

As the definition domain of wavelet transformation is a double infinite interval, the wavelet coefficients are large near the endpoints of the bounded signal, which increases the computational error. To obtain a high computational precision, the additional condition of scalar fonctionelle in boundary could be changed to

Thus

where is a known value. According to the general variational principle, we obtain

where

Correspondingly, (7) changed to the form

The interpolation properties of and imply that and are the vectors consisting of the approximate values of in discrete points and extended discrete points outside the domain of , respectively. Typically, periodic, symmetric, and zero extension methods could be employed. Assuming that there are discrete points in , that is, , and extended points in both outsides of the domain , respectively, that is, and , the function can be extended outside of the definition domain using the symmetry method as below:

Substituting (13) into (10) results in

where

In the matrices and , and could be calculated, respectively, as

According to and the relational expression (13), (14) could be condensed in the following way:

Therefore, the function () can be represented as

where

Subsequently, as interval interpolation basic functions, the interval Shannon wavelet can be obtained from (18) as follows:

where is the support domain of the wavelet function; that is, sup .

The interval wavelet constructed in [8] is different from (20), which is represented as follows:

where

where is the amount of the external collocation points, the amount of discrete points in the definition domain is , and is the definition domain of the approximated function.

It is obvious that (21) is just the Lagrange extrapolation representation. The parameter appeared in (21) is relative to the smoothness and gradient of the function to be approximated. In theory, the increase of can improve the approximate precision. But it was pointed out in [23] that the increase of can also bring the increase of the condition number of the problem to be solved, and this can decrease the numerical precision greatly. In addition, the larger can also bring more calculation.

Table 1 shows the maximum absolute errors of the numerical approximation to curves with different gradient near the boundary corresponding to different value of . It illustrates that is not only related to the smoothness, but also to the gradient of original function near the boundary. For the same , the greater the original function’s gradient is, the bigger the error is. To the different gradient original functions in Table 1, the reasonable should be , , , and , respectively. In contrast, the interval wavelet proposed in this paper is not sensitive to the gradient of the function to be approximated, and the calculation amount is the same as the interval wavelet in [7] with . As the matter of fact, the proposed interval wavelet is constructed based on that the derivative function should be continuous, and the smaller avoids the Gibbs phenomenon appeared in the Lagrange interpolator. This is the reason why it is not sensitive to the gradient value.

Table 1: Errors comparison of the numerical approximation to curves with different gradients (, ).

()

(10⁻¹⁵)
(5)

(10⁻¹⁵)
(10)

(10⁻¹⁵)
(100)

(10⁻¹⁴)

Interval wavelet in [7]
1
3964.88397669
123954.735364
793825.449818
2336702.80902

2
5.60662627435
634.825525480
1759.70349403
599685.701097

3
0.22204460492
2.22044604925
46.9624339416
18207.2357191

4
0.22204460492
0.44408920985
0.44408920985
3343.21459405

5
0.27755575615
0.44408920985
0.44408920985
129.352084599

6
0.22204460492
0.44408920985
0.44408920985
18.0300219199

7
0.27755575615
0.44408920985
0.44408920985
0.88817841970

8
0.27755575615
0.44408920985
0.44408920985
0.88817841970

New interval wavelet
0.27755575613
0.27755575625
0.27755575615
0.02784584512

3. Black-Scholes Equation and Its Interval Wavelet Approximation

The famous Black-Scholes equation can be represented as

where () denotes the underlying asset, , denotes the expiry date, is the volatility (measures the standard deviation of the returns), and is the riskless interest rate. In general, the parameter is taken as a constant since the transaction cost is taken as zero. Obviously, the is not really a constant, and then we can obtain the nonlinear Black-Scholes equation as follows:

where denotes a nonconstant volatility. In order to solve the problem, it is necessary to perform a variable transformation as follows:

Substituting (25) into (24), the following equation can be obtained:

where

Initial condition is

And the boundary conditions are

According to the transformation relation (25), it is easy to understand that the point is corresponding to the strike price . Obviously, the initial solution curve is smooth except around the position , where a sharp steep wave happened. So, an adaptive numerical method is suited to this problem.

According to the wavelet precision integration method, the nonlinear Black-Scholes equation can be approximately represented as

where . Let

Then the system of (30) can be expressed as the matrix format:

It is obvious that is a constant matrix according to (32). ()   is the nonlinear term in (32), which can be linearized subsectionally in terms of first-order Taylor series expansion; that is, the nonlinear terms can be taken as linear ones in the time step (,), and then (32) can be rewritten as

where

Letting and combining with (32), the expression of can be obtained as

The precise time integration scheme for the linear ordinary differential equations (33) can be expressed as

where

and then solving (33) is changed into the calculation of matrix . the exponential matrix can be calculated accurately by PIM as follows:

Let , where is a positive integer (usually take , and then ). Taking into account that is a small time interval and Δ is a very smaller value, then

which is the Taylor series expansion of exp(Δ). In order to calculate the matrix more accurately, it is necessary to factorize the matrix as

After doing times of factorization as above, a more accurate solution of can be obtained.

The calculation of the exponent matrix at different time steps is needed in solving nonlinear equations through iteration based on the precise integration method, and the algorithm of the matrix presented here can obtain all the matrices at different time steps for once.

The numerical solution of the nonlinear Black-Scholes equation is shown in Figure 1, where the riskless rate , the volatility , and strike price . With the increase of the parameter , the smoothness of the solution function near the strike price becomes better and better, and the amount of the collocation points is decreasing correspondingly. This illustrates the performance of the adaptive wavelet numerical method.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 1: Evolution of the call option price with the parameter .

Compared with the wavelet precision integration method employed in Yan’s method [5], the amount of the collocation points in the interval wavelet precision integration method proposed in this paper is decreased greatly, especially around the boundary of the definition domain. This illustrates that the interval wavelet precision integration method can eliminate the boundary effect completely. It is obviously that the decrease of the collocation points can improve the efficiency of the algorithm. The performance of Mei’s interval wavelet [4] relates to the external interpolation points (), and the choice scheme of was not given in [3]. It was pointed out in [23] that the bigger usually introduces bigger condition number which can decrease the precision greatly. This can be illustrated in Figures 1(d) and 1(e). As , the boundary effect can be eliminated effectively. On the contrary, larger almost cannot restrict it.

In addition, the interval wavelet can also improve the numerical precision compared with the common wavelet method. For convenience, we take the linear Black-Scholes equation as an example, which has exact analytical solution as follows:

where

is the call price, is the underlying asset price, is the strike price, is the riskless rate, is the maturity, is the volatility, and express the normal distribution.

All the comparisons in Table 2 are made qualitatively by comparing the calculation precision in the same time step and space mesh grid size. The first measure of error is given by

which provides a measure of the accuracy of the solution near the boundary. The second measure of error is given by

In Table 2, the parameter is taken as 1 in Mei’s method. The error of Mei’s method almost equals that of Yan’s method, and the error e ₂ of Mei’s method is smaller evidently than Yan’s method. This indicates that the interval wavelet method proposed by Mei can eliminate the boundary effect to some extent as , but it cannot decrease the maximum of the absolute error evidently. In addition, with the increasing of the parameters and , the errors of both Yan and Mei’s method become larger and larger evidently. This can be explained by the condition number theory [23]. On the contrary, the interval wavelet method proposed in this paper is not sensitive to the parameters and , and the precision is the best among all above methods.

Table 2: The numerical precision comparison.

Interval WPIM
Mei’s method
Yan’s method
Interval WPIM
Mei’s method
Yan’s method

12
10

100

200

15
10

100

200

0.3812
0.3971

0.0211
0.1195

4. Conclusion

The solution of the nonlinear Black-Scholes equation has an evident sharp steep wave near the boundary, and this can introduce evident boundary effect in the wavelet approximation such as Yan’s method. The interval wavelet approximation proposed in this paper makes sure that the derivative function near the boundary of the definition domain is continued, which is helpful to eliminate the boundary effect and improve the numerical precision in the whole definition domain. Compared with Mei’s interval wavelet method, the interval wavelet proposed in this paper is not sensitive to the amount of the levels and the time step and does not have to choose the external collocation points. So, it is a simple and robust tool for solving the Black-Scholes equation. In fact, we can construct various interval wavelet functions aiming at different problems in engineering based on the general variational principle.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (no. 41171337) and the National Key Technology R&D Program of China (no. 2012BAD35B02).

References

J. Ankudinova and M. Ehrhardt, “On the numerical solution of nonlinear Black-Scholes equations,” Computers & Mathematics with Applications, vol. 56, no. 3, pp. 799–812, 2008.

L. Song and W. Wang, “Solution of the fractional Black-Scholes option pricing model by finite difference method,” Abstract and Applied Analysis, vol. 2013, Article ID 194286, 10 pages, 2013.

Z. Cen, A. Le, and A. Xu, “Exponential time integration and second-order difference scheme for a generalized Black-Scholes equation,” Journal of Applied Mathematics, vol. 2012, Article ID 796814, 12 pages, 2012.

S. L. Mei, Q. S. Lu, S. W. Zhang, and L. Jin, “Adaptive interval wavelet precise integration method for partial differential equations,” Applied Mathematics and Mechanics, vol. 26, no. 3, pp. 364–371, 2005.

H.-H. Yan, “Adaptive wavelet precise integration method for nonlinear Black-Scholes model based on variational iteration method,” Abstract and Applied Analysis, vol. 2013, Article ID 735919, 6 pages, 2013.

A. A. Elbeleze, A. Kılıçman, and B. M. Taib, “Homotopy perturbation method for fractional Black-Scholes European option pricing equations using Sumudu transform,” Mathematical Problems in Engineering, vol. 2013, Article ID 524852, 7 pages, 2013.

S.-L. Mei, H.-L. Lv, and Q. Ma, “Construction of interval wavelet based on restricted variational principle and its application for solving differential equations,” Mathematical Problems in Engineering, vol. 2008, Article ID 629253, 14 pages, 2008.

S.-L. Mei, “Construction of target controllable image segmentation model based on homotopy perturbation technology,” Abstract and Applied Analysis, vol. 2013, Article ID 131207, 8 pages, 2013.

R.-Y. Xing, “Wavelet-based homotopy analysis method for nonlinear matrix system and its application in Burgers equation,” Mathematical Problems in Engineering, vol. 2013, Article ID 982810, 7 pages, 2013.

S.-L. Pang, “Wavelet numerical method for nonlinear random system,” Transactions of the Chinese Society of Agricultural Machinery, vol. 38, no. 3, pp. 168–170, 2007.

S.-L. Mei, C.-J. Du, and S.-W. Zhang, “Linearized perturbation method for stochastic analysis of a rill erosion model,” Applied Mathematics and Computation, vol. 200, no. 1, pp. 289–296, 2008.

S.-L. Mei, “HAM-based adaptive multiscale meshless method for Burgers equation,” Journal of Applied Mathematics, vol. 2013, Article ID 248246, 10 pages, 2013.

J. H. He, “Variational iteration method—a kind of non-linear analytical technique: some examples,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 699–708, 1999.

J. H. He, “Variational iteration method for autonomous ordinary differential systems,” Applied Mathematics and Computation, vol. 114, no. 2-3, pp. 115–123, 2000.

J. H. He, G. C. Wu, and F. Austin, “The variational iteration method which should be followed,” Nonlinear Science Letters A, vol. 1, pp. 1–30, 2010.

J. H. He and X. H. Wu, “Variational iteration method: new development and applications,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 881–894, 2007.

J. H. He, “Asymptotic methods for solitary solutions and compactons,” Abstract and Applied Analysis, vol. 2012, Article ID 916793, 130 pages, 2012.

C. Cattani, “Harmonic wavelets towards the solution of nonlinear PDE,” Computers & Mathematics with Applications, vol. 50, no. 8-9, pp. 1191–1210, 2005.

C. Cattani, “Connection coefficients of Shannon wavelets,” Mathematical Modelling and Analysis, vol. 11, no. 2, pp. 117–132, 2006.

C. Cattani, “Shannon wavelets theory,” Mathematical Problems in Engineering, vol. 2008, Article ID 164808, 24 pages, 2008.

C. Cattani, “Shannon wavelets for the solution of integrodifferential equations,” Mathematical Problems in Engineering, vol. 2010, Article ID 408418, 22 pages, 2010.

D. K. Hoffman, G. W. Wei, D. S. Zhang, and D. J. Kouri, “Shannon-Gabor wavelet distributed approximating functional,” Chemical Physics Letters, vol. 287, no. 1-2, pp. 119–124, 1998.

S.-L. Mei and D.-H. Zhu, “Interval Shannon wavelet collocation method for fractional Fokker-Planck equation,” Advances in Mathematical Physics, vol. 2013, Article ID 821820, 12 pages, 2013.