Numerical simulation of a Finite Moment Log Stable model for a European call option

Compared to the classical Black-Scholes model for pricing options, the Finite Moment Log Stable (FMLS) model can more accurately capture the dynamics of the stock prices including large movements or jumps over small time steps. In this paper, the FMLS model is written as a fractional partial differential equation and we will present a new numerical scheme for solving this model. We construct an implicit numerical scheme with second order accuracy for the FMLS and consider the stability and convergence of the scheme. In order to reduce the storage space and computational cost, we use a fast bi-conjugate gradient stabilized method (FBi-CGSTAB) to solve the discrete scheme. A numerical example is presented to show the efficiency of the numerical method and to demonstrate the order of convergence of the implicit numerical scheme. Finally, as an application, we use the above numerical technique to price a European call option. Furthermore, by comparing the FMLS model with the classical B-S model, the characteristics of the FMLS model are also analyzed.     

A new approximation for the generalized fractional derivative and its application to generalized fractional diffusion equation

In this paper, we derive a fourth order approximation for the generalized fractional derivative that is characterized by a scale function z(t) and a weight function w(t) . Combining the new approximation with compact finite difference method, we develop a numerical scheme for a generalized fractional diffusion problem. The stability and convergence of the numerical scheme are proved by the energy method, and it is shown that the temporal and spatial convergence orders are both 4. Several numerical experiments are provided to illustrate the efficiency of our scheme.     

High-order computational methods for option valuation under multifactor models

Many of the different numerical techniques in the partial differential equations framework for solving option pricing problems have employed only standard second-order discretization schemes. A higher-order discretization has the advantage of producing low size matrix systems for computing sufficiently accurate option prices and this paper proposes new computational schemes yielding high-order convergence rates for the solution of multi-factor option problems. These new schemes employ Galerkin finite element discretizations with quadratic basis functions for the approximation of the spatial derivatives in the pricing equations for stochastic volatility and two-asset option problems and time integration of the resulting semi-discrete systems requires the computation of a single matrix exponential. The computations indicate that this combination of high-order finite elements and exponential time integration leads to efficient algorithms for multi-factor problems. Highly accurate European prices are obtained with relatively coarse meshes and high-order convergence rates are also observed for options with the American early exercise feature. Various numerical examples are provided for illustrating the accuracy of the option prices for Heston’s and Bates stochastic volatility models and for two-asset problems under Merton’s jump-diffusion model.     

Finite Difference/Collocation Method for Two-Dimensional Sub-Diffusion Equation with Generalized Time Fractional Derivative

In this paper, we propose a finite difference/collocation method for two-dimensional time fractional diffusion equation with generalized fractional operator. The main purpose of this paper is to design a high order numerical scheme for the new generalized time fractional diffusion equation. First, a finite difference approximation formula is derived for the generalized time fractional derivative, which is verified with order $2-\alpha$ $(0<\alpha<1)$. Then, collocation method is introduced for the two-dimensional space approximation. Unconditional stability of the scheme is proved. To make the method more efficient, the alternating direction implicit method is introduced to reduce the computational cost. At last, numerical experiments are carried out to verify the effectiveness of the scheme.     

The Numerical Simulation of Space-Time Variable Fractional Order Diffusion Equation

Many physical processes appear to exhibit fractional order behavior that may vary with time or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. Numerical methods and analysis of stability and convergence of numerical scheme for the variable fractional order partial differential equations are quite limited and difficult to derive. This motivates us to develop efficient numerical methods as well as stability and convergence of the implicit numerical methods for the space-time variable fractional order diffusion equation on a finite domain. It is worth mentioning that here we use the Coimbra-definition variable time fractional derivative which is more efficient from the numerical standpoint and is preferable for modeling dynamical systems. An implicit Euler approximation is proposed and then the stability and convergence of the numerical scheme are investigated. Finally, numerical examples are provided to show that the implicit Euler approximation is computationally efficient.     

A series of high‐order quasi‐compact schemes for space fractional diffusion equations based on the superconvergent approximations for fractional derivatives

Based on the superconvergent approximation at some point (depending on the fractional order α, but not belonging to the mesh points) for Grünwald discretization to fractional derivative, we develop a series of high‐order quasi‐compact schemes for space fractional diffusion equations. Because of the quasi‐compactness of the derived schemes, no points beyond the domain are used for all the high‐order schemes including second‐order, third‐order, fourth‐order, and even higher‐order schemes; moreover, the algebraic equations for all the high‐order schemes have the completely same matrix structure. The stability and convergence analysis for some typical schemes are made; the techniques of treating the fractional derivatives with nonhomogeneous boundaries are introduced; and extensive numerical experiments are performed to confirm the theoretical analysis or verify the convergence orders. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1345–1381, 2015     

Finite difference scheme for time-space fractional diffusion equation with fractional boundary conditions

In this paper, the finite difference scheme is developed for the time-space fractional diffusion equation with Dirichlet and fractional boundary conditions. The time and space fractional derivatives are considered in the senses of Caputo and Riemann-Liouville, respectively. The stability and convergence of the proposed numerical scheme are strictly proved, and the convergence order is O(τ^2−α+h²). Numerical experiments are performed to confirm the accuracy and efficiency of our scheme.     

A Fourier method and an extrapolation technique for Stokes’ first problem for a heated generalized second grade fluid with fractional derivative

In this article, we consider Stokes’ first problem for a heated generalized second grade fluid with fractional derivative (SFP-HGSGF). Implicit and explicit numerical approximation schemes for the SFP-HGSGF are presented. The stability and convergence of the numerical schemes are discussed using a Fourier method. In addition, the solvability of the implicit numerical approximation scheme is also analyzed. A Richardson extrapolation technique for improving the order of convergence of the implicit scheme is proposed. Finally, a numerical test is given. The numerical results demonstrate the good performance of our theoretical analysis.     

High order finite difference methods on non-uniform meshes for space fractional operators

In the past decades, the finite difference methods for space fractional operators develop rapidly; to the best of our knowledge, all the existing finite difference schemes, including the first and high order ones, just work on uniform meshes. The nonlocal property of space fractional operator makes it difficult to design the finite difference scheme on non-uniform meshes. This paper provides a basic strategy to derive the first and high order discretization schemes on non-uniform meshes for fractional operators. And the obtained first and second schemes on non-uniform meshes are used to solve space fractional diffusion equations. The error estimates and stability analysis are detailedly performed; and extensive numerical experiments confirm the theoretical analysis or verify the convergence orders.     

Analytical and numerical solutions of a two‐dimensional multi‐term time‐fractional Oldroyd‐B model

In this paper, we consider a two‐dimensional multi‐term time‐fractional Oldroyd‐B equation on a rectangular domain. Its analytical solution is obtained by the method of separation of variables. We employ the finite difference method with a discretization of the Caputo time‐fractional derivative to obtain an implicit difference approximation for the equation. Stability and convergence of the approximation scheme are established in the L_∞ ‐norm. Two examples are given to illustrate the theoretical analysis and analytical solution. The results indicate that the present numerical method is effective for this general two‐dimensional multi‐term time‐fractional Oldroyd‐B model.     

Bivariate hierarchical Hermite spline quasi-interpolation

Spline quasi-interpolation (QI) is a general and powerful approach for the construction of low cost and accurate approximations of a given function. In order to provide an efficient adaptive approximation scheme in the bivariate setting, we consider quasi-interpolation in hierarchical spline spaces. In particular, we study and experiment the features of the hierarchical extension of the tensor-product formulation of the Hermite BS quasi-interpolation scheme. The convergence properties of this hierarchical operator, suitably defined in terms of truncated hierarchical B-spline bases, are analyzed. A selection of numerical examples is presented to compare the performances of the hierarchical and tensor-product versions of the scheme.     

A fast second‐order difference scheme for the space–time fractional equation

In this paper, a fast second‐order accurate difference scheme is proposed for solving the space–time fractional equation. The temporal Caputo derivative is approximated by ?L2 ‐1_σ formula which employs the sum‐of‐exponential approximation to the kernel function appeared in Caputo derivative. The second‐order linear spline approximation is applied to the spatial Riemann–Liouville derivative. At each time step, a fast algorithm, the preconditioned conjugate gradient normal residual method with a circulant preconditioner (PCGNR), is used to solve the resulting system that reduces the storage and computational cost significantly. The unique solvability and unconditional convergence of the difference scheme are shown by the discrete energy method. Numerical examples are given to verify numerical accuracy and efficiency of the difference schemes.     

美式期权定价问题的数值方法

本文研究美式股票看跌期权定价问题的数值方法。通过将问题转化为等价的变分不等式方程，分别建立了半离散和全离散有限元逼近格式。并给出了有限元解的收敛性和稳定性分析。数值实验表明本文算法是一个高效和收敛的算法。     

A fully discrete spectral method for fractional Cattaneo equation based on Caputo–Fabrizo derivative

Recently Caputo and Fabrizio introduced a new derivative with fractional order without singular kernel. The derivative can be used to describe the material heterogeneities and the fluctuations of different scales. In this article, we derived a new discretization of Caputo–Fabrizio derivative of order α (1 < α < 2) and applied it into the Cattaneo equation. A fully discrete scheme based on finite difference method in time and Legendre spectral approximation in space is proposed. The stability and convergence of the fully discrete scheme are rigorously established. The convergence rate of the fully discrete scheme in H¹ norm is O(τ² + N^1?m), where τ, N and m are the time‐step size, polynomial degree and regularity in the space variable of the exact solution, respectively. Furthermore, the accuracy and applicability of the scheme are confirmed by numerical examples to support the theoretical results.     

Fast evaluation and high accuracy finite element approximation for the time fractional subdiffusion equation

In this article, an efficient algorithm for the evaluation of the Caputo fractional derivative and the superconvergence property of fully discrete finite element approximation for the time fractional subdiffusion equation are considered. First, the space semidiscrete finite element approximation scheme for the constant coefficient problem is derived and supercloseness result is proved. The time discretization is based on the L1‐type formula, whereas the space discretization is done using, the fully discrete scheme is developed. Under some regularity assumptions, the superconvergence estimate is proposed and analyzed. Then, extension to the case of variable coefficients is also discussed. To reduce the computational cost, the fast evaluation scheme of the Caputo fractional derivative to solve the fractional diffusion equations is designed. Finally, numerical experiments are presented to support the theoretical results.     

Numerical method for generalized time fractional KdV-type equation

In this article, an efficient numerical method for linearized and nonlinear generalized time-fractional KdV-type equations is proposed by combining the finite difference scheme and Petrov–Galerkin spectral method. The scale and weight functions involved in generalized fractional derivative cause too much difficulty in discretization and numerical analysis. Fortunately, motivated by finite difference method for fractional differential equation on graded mesh, the stability and convergence of the constructed method are established rigorously. It is proved that the full discretization schemes of generalized time-fractional KdV-type equation is unconditionally stable in linear case. While for nonlinear case, it is stable under a CFL condition and for not small ϵ, coefficient of the high-order spatial differential term. In addition, the full discretization schemes with respect to linear and nonlinear cases respectively converge to the associated exact solutions with orders and , where τ, α, N and m accordingly indicate the time step size, the order of the fractional derivative, polynomial degree, and regularity of the exact solution. Numerical experiments are carried out to support the theoretical results.     

A General Approximation Scheme for Attractors of Abstract Parabolic Problems

We consider semilinear problems of the form u′ = Au + f(u), where A generates an exponentially decaying compact analytic C ₀-semigroup in a Banach space E, f:E ^α → E is differentiable globally Lipschitz and bounded (E ^α = D((?A)^α) with the graph norm). Under a very general approximation scheme, we prove that attractors for such problems behave upper semicontinuously. If all equilibrium points are hyperbolic, then there is an odd number of them. If, in addition, all global solutions converge as t → ±∞, then the attractors behave lower semicontinuously. This general approximation scheme includes finite element method, projection and finite difference methods. The main assumption on the approximation is the compact convergence of resolvents, which may be applied to many other problems not related to discretization.     

The upwind finite difference fractional steps method for combinatorial system of dynamics of fluids in porous media and its application

For combinatorial system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as implicit-explicit difference scheme, calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L ² norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources. Keywords: combinatorial system, multilayer dynamics of fluids in porous media, two-class upwind finite difference fractional steps method, convergence, numerical simulation of energy sources.     

Testing robustness in calibration of stochastic volatility models

Many numerical aspects are involved in parameter estimation of stochastic volatility models. We investigate a model for stochastic volatility suggested by Hobson and Rogers [Complete models with stochastic volatility, Mathematical Finance 8 (1998) 27] and we focus on its calibration performance with respect to numerical methodology.In recent financial literature there are many papers dealing with stochastic volatility models and their capability in capturing European option prices; in Figà-Talamanca and Guerra [Towards a coherent volatility pricing model: An empirical comparison, Financial Modelling, Phisyca-Verlag, 2000] a comparison between some of the most significant models is done. The model proposed by Hobson and Rogers seems to describe quite well the dynamics of volatility.In Figà-Talamanca and Guerra [Fitting the smile by a complete model, submitted] a deep investigation of the Hobson and Rogers model was put forward, introducing different ways of parameters' estimation. In this paper we test the robustness of the numerical procedures involved in calibration: the quadrature formula to compute the integral in the definition of some state variables, called offsets, that represent the weight of the historical log-returns, the discretization schemes adopted to solve the stochastic differential equation for volatility and the number of simulations in the Monte Carlo procedure introduced to obtain the option price.The main results can be summarized as follows. The choice of a high order of convergence scheme is not fully justified because the option prices computed via calibration method are not sensitive to the use of a scheme with 2.0 order of convergence or greater. The refining of the approximation rule for the integral, on the contrary, allows to compute option prices that are often closer to market prices. In conclusion, a number of 10 000 simulations seems to be sufficient to compute the option price and a higher number can only slow down the numerical procedure.     

Fast calculation based on a spatial two-grid finite element algorithm for a nonlinear space–time fractional diffusion model

In this article, a spatial two-grid finite element (TGFE) algorithm is used to solve a two-dimensional nonlinear space–time fractional diffusion model and improve the computational efficiency. First, the second-order backward difference scheme is used to formulate the time approximation, where the time-fractional derivative is approximated by the weighted and shifted Grünwald difference operator. In order to reduce the computation time of the standard FE method, a TGFE algorithm is developed. The specific algorithm is to iteratively solve a nonlinear system on the coarse grid and then to solve a linear system on the fine grid. We prove the scheme stability of the TGFE algorithm and derive a priori error estimate with the convergence result O(Δt² + h^{r + 1 − η} + H^{2r + 2 − 2η}) . Finally, through a two-dimensional numerical calculation, we improve the computational efficiency and reduce the computation time by the TGFE algorithm.