Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 2371826, 6 pages

http://dx.doi.org/10.1155/2016/2371826

## Sparsity-Homotopy Perturbation Inversion Method with Wavelets and Applications to Black-Scholes Model and Todaro Model

^{1}School of Finance, Harbin University of Commerce, Harbin 150028, China^{2}School of Management, Harbin University of Science and Technology, Harbin 150080, China

Received 17 April 2016; Accepted 28 June 2016

Academic Editor: Thomas Schuster

Copyright © 2016 Yixin Dou and Zhihao Wang. 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

Sparsity regularization method plays an important role in reconstructing parameters. Compared with traditional regularization methods, sparsity regularization method has the ability to obtain better performance for reconstructing sparse parameters. However, sparsity regularization method does not have the ability to reconstruct smooth parameters. For overcoming this difficulty, we combine a sparsity regularization method with a wavelet method in order to transform smooth parameters into sparse parameters. We use a sparsity-homotopy perturbation inversion method to improve the accuracy and stability and apply the proposed method to reconstruct parameters for a Black-Scholes option pricing model and a Todaro model. Numerical experiments show that the proposed method is convergent and stable.

#### 1. Introduction

The reconstruction of parameters plays an important role in financial mathematics, such as the reconstructions of a volatility and a policy parameter [1, 2] and, in other fields, imaging enhancement techniques [3, 4], seismic signals [5, 6], and electrocardiogram signals (ECG) [7, 8]. With the development of economy and financial mathematics, the reconstructions of a volatility and a policy parameter have been widely used in many real applications. In general, the reconstruction of parameters is ill-posed. In other words, the small noisy level of measurement data can lead to the large error of reconstruction [9]. In order to overcome the ill-posedness, some regularization methods are developed. The most popular method is Tikhonov regularization method, which is composed of a fitting term and a penalty term in norm. The aims of those two terms are to match measurement data and to suppress noises, respectively.

The numerical methods for Tikhonov regularization method have been conducted, such as a well-known Landweber method [10], a Gauss-Newton method [11], a regularizing Newton-Kaczmarz method [12], and a multiscale smoothing method [13]. These methods have the ability to reconstruct smooth parameters in the case of sufficient measurement data. When measurement data are limited, smooth parameters are very difficult to be reconstructed. In the fields of economy and finance, reconstructed parameters are smooth and measurement data are limited. Hence, we need to use wavelet and curvelet transformations from smooth parameters to sparse parameters (i.e., the number of nonzero elements of parameter is very limited). Sparsity regularization methods are used to reconstruct sparse parameters. Compared with Tikhonov regularization method, sparsity regularization methods are not differentiable, and hence some specific techniques were developed to overcome this difficulty, such as Bregman iterations [14–17]. For reducing the computational time, a homotopy perturbation inversion method has been widely used in real applications, such as nonlinear systems, optimal control, and heat transfer equation [18–20].

In this paper, we combine a homotopy perturbation inversion method with a sparsity regularization method, in order to improve the accuracy and stability. After introducing the proposed method, we reconstruct two economic parameters based on a Black-Scholes option pricing model and a Todaro model. From the numerical experiments, a sparsity-homotopy perturbation inversion method with wavelets is convergent and stable.

#### 2. Sparsity Regularization Method

The reconstruction of parameters is ill-posed, and hence we should employ the regularization method to reconstruct parameters. The Tikhonov regularization method can reconstruct smooth parameters while measurement data are sufficient. However, measurement data are limited in real applications. We combine a sparsity regularization method with a wavelet method, in order to improve the performance of reconstructing smooth parameters for limited measurement data.

##### 2.1. Tikhonov Regularization Method

Many inverse problems can be formulated as a nonlinear operator equation where denote a nonlinear operator, reconstructed parameter, and measurement data, respectively. We assume measurement data contaminated by noise where stand for the real measurement data and noisy level in the norm.

Ill-posedness means that the small noisy level included in measurement data may lead to the large error of reconstruction. We apply the regularization method in order to overcome ill-posedness. It is very important to balance a fitting term and a penalty term. Next, we introduce a Tikhonov regularization functional as follows: where is a fitting term, is a penalty term, and is a regularized parameter balancing the fitting term and the penalty term. For reconstruction, we should minimize the functional (3):

The minimizer satisfies the Euler equation: where is the Fréchet derivative. The Landweber method is widely used for solving (5). The Landweber method can be written as follows: where denotes the iteration number. Equation (6) is a well-known Landweber method that is stable; however, the convergent rate is slow and the accuracy is low.

##### 2.2. Sparsity Regularization Method with Wavelets

The Tikhonov regularization method can obtain better reconstruction for smoothing parameters; however, its performance for reconstructing sparse parameters is worse. Hence, the classical Tikhonov regularization method is modified as the following sparsity regularization method: where the norm of means the number of nonzeros in . The minimization problem (7) has the ability to obtain the best reconstruction for sparse parameters. However, the penalty term of (7) is nondifferentiable and minimization problem (7) is a NP problem. For speeding up minimization problem (7), we rewrite (7) as where the norm of means **. In functional (8), the norm replaces the norm in order to overcome a NP problem. This improvement can reduce the cost time of (7) dramatically.

Because the penalty term of (8) is nondifferentiable, we modify problem (8) as follows: where replaces , and the auxiliary parameter is a positive real number. As the Tikhonov regularization method, we can apply the Landweber method to minimize functional (9):

For a smooth parameter, we combine a sparsity regularization method with a wavelet method. A wavelet transformation converts a smooth parameter into a sparse parameter. A wavelet transformation is written as where stand for an inverse transformation matrix and the sparse representation of , respectively. Taking (11) into (9) and (10), we arrive at

#### 3. Sparsity-Homotopy Perturbation Inversion Method

Because the Landweber method (13) is very slow, it is not realistic to reconstruct large scale real applications. We use a homotopy perturbation inversion method to modify the classical Landweber method. For convenience, we omit and replace with .

Setting a homotopy mapping where is an embedding parameter and is an initial guess value. Hence,

We write as the power series of , and obtain the approximation of functional (12):

We expand in (14) as a Taylor series near : Following the power of , we can get Following formula (20), the parameter for the noisy measurement data is reconstructed by the first two terms: where is the iteration number.

When the parameter is reconstructed by the first term, Equation (22) is a well-known Landweber method. Equation (21) is a modified version of (22), which is called a homotopy perturbation inversion method. The convergent rate is faster and the accuracy is higher than a Landweber method. By using a wavelet transformation, we can obtain

#### 4. Numerical Experiments

##### 4.1. Reconstructing Sparse Volatility

The reconstruction of a volatility is investigated based on the Black-Scholes (B-S) option pricing model. We use a finite difference method to solve the forward problem. In order to improve the accuracy and convergent rate of a Landweber method, we design a homotopy perturbation inversion method to reconstruct a volatility. The B-S formula is widely used in the field of a derivative pricing, when the price changes of a derivative satisfy the standard geometric Brown motion. The boundary conditions of the different values vary with the different types of derivatives. When boundary conditions are given, a derivative pricing can be obtained by solving the B-S formula of the derivative pricing model. When the derivative is taken as option, we define the forward problem as the determination of the option pricing. For simplicity, we take the European call option as example.

The relationship between the European call and put options shows that the reconstructed volatility should be the same by using the call options market quotes or the put options market quotes. For the European call option in the time interval , let be the European call option pricing, and satisfies where is stock prices with time, is strike price, is interest rate, is dividend, is availability period, is time, and is volatility.

In the B-S formula, except that the volatility is a free variable, the other parameters and variables are decided by the prevailing market conditions or contracts. Hence, the volatility is a very important parameter. For an option, we can infer the volatility from the B-S formula.

The inverse problem is defined as follows: reconstructing the local volatility from the measurement data and .

We define a nonlinear vector-valued function , namely, . For testing the performance of sparse regularization and Tikhonov regularization methods, we set and set stock prices , interest rate , and strike price . We add 0%, 0.5%, and 1% Gaussian random noises to the measurement data in order to test the stability. The exact volatility is provided as follows:

We transform the smooth volatility into the sparse volatility where is a sparse representation. Functional (12) is rewritten as

The results of the proposed method (I) and Tikhonov method (II) are listed in Table 1.