Differential quadrature method is employed by numerous researchers due to its numerical accuracy and computational efficiency,and is mentioned as potential alternative of conventional numerical methods.In this paper,a...Differential quadrature method is employed by numerous researchers due to its numerical accuracy and computational efficiency,and is mentioned as potential alternative of conventional numerical methods.In this paper,a differential quadrature based numerical scheme is developed for solving volterra partial integro-differential equation of second order having a weakly singular kernel.The scheme uses cubic trigonometric B-spline functions to determine the weighting coefficients in the differential quadrature approximation of the second order spatial derivative.The advantage of this approximation is that it reduces the problem to a first order time dependent integro-differential equation(IDE).The proposed scheme is obtained in the form of an algebraic system by reducing the time dependent IDE through unconditionally stable Euler backward method as time integrator.The scheme is validated using a homogeneous and two nonhomogeneous test problems.Conditioning of the system matrix and numerical convergence of the method are analyzed for spatial and temporal domain discretization parameters.Comparison of results of the present approach with Sinc collocation method and quasi-wavelet method are also made.展开更多
This article studies the development of two numerical techniques for solving convection-diffusion type partial integro-differential equation(PIDE)with a weakly singular kernel.Cubic trigonometric B-spline(CTBS)functio...This article studies the development of two numerical techniques for solving convection-diffusion type partial integro-differential equation(PIDE)with a weakly singular kernel.Cubic trigonometric B-spline(CTBS)functions are used for interpolation in both methods.The first method is CTBS based collocation method which reduces the PIDE to an algebraic tridiagonal system of linear equations.The other method is CTBS based differential quadrature method which converts the PIDE to a system of ODEs by computing spatial derivatives as weighted sum of function values.An efficient tridiagonal solver is used for the solution of the linear system obtained in the first method as well as for determination of weighting coefficients in the second method.An explicit scheme is employed as time integrator to solve the system of ODEs obtained in the second method.The methods are tested with three nonhomogeneous problems for their validation.Stability,computational efficiency and numerical convergence of the methods are analyzed.Comparison of errors in approximations produced by the present methods versus different values of discretization parameters and convection-diffusion coefficients are made.Convection and diffusion dominant cases are discussed in terms of Peclet number.The results are also compared with cubic B-spline collocation method.展开更多
This paper develops and analyzes a moving mesh finite difference method for solving partial integro-differential equations. First, the time-dependent mapping of the coordinate transformation is approximated by a a pie...This paper develops and analyzes a moving mesh finite difference method for solving partial integro-differential equations. First, the time-dependent mapping of the coordinate transformation is approximated by a a piecewise linear function in time. Then, piecewise quadratic polynomial in space and an efficient method to discretize the memory term of the equation is designed using the moving mesh approach. In each time slice, a simple piecewise constant approximation of the integrand is used, and thus a quadrature is constructed for the memory term. The central finite difference scheme for space and the backward Euler scheme for time are used. The paper proves that the accumulation of the quadrature error is uniformly bounded and that the convergence of the method is second order in space and first order in time. Numerical experiments are carried out to confirm the theoretical predictions.展开更多
A framework for the optimal sparse-control of the probability density function of a jump-diffusion process is presented. This framework is based on the partial integro-differential Fokker-Planck (FP) equation that gov...A framework for the optimal sparse-control of the probability density function of a jump-diffusion process is presented. This framework is based on the partial integro-differential Fokker-Planck (FP) equation that governs the time evolution of the probability density function of this process. In the stochastic process and, correspondingly, in the FP model the control function enters as a time-dependent coefficient. The objectives of the control are to minimize a discrete-in-time, resp. continuous-in-time, tracking functionals and its L2- and L1-costs, where the latter is considered to promote control sparsity. An efficient proximal scheme for solving these optimal control problems is considered. Results of numerical experiments are presented to validate the theoretical results and the computational effectiveness of the proposed control framework.展开更多
We consider pricing options in a jump-diffusion model which requires solving a partial integro-differential equation.Discretizing the spatial direction with a fourth order compact scheme leads to a linear system of or...We consider pricing options in a jump-diffusion model which requires solving a partial integro-differential equation.Discretizing the spatial direction with a fourth order compact scheme leads to a linear system of ordinary differential equations.For the temporal direction,we utilize the favorable boundary value methods owing to their advantageous stability properties.In addition,the resulting large sparse system can be solved rapidly by the GMRES method with a circulant Strang-type preconditioner.Numerical results demonstrate the high order accuracy of our scheme and the efficiency of the preconditioned GMRES method.展开更多
This article establishes a universal robust limit theorem under a sublinear expectation framework.Under moment and consistency conditions,we show that,forα∈(1,2),the i.i.d.sequence{(1/√∑_(i=1)^(n)X_(i),1/n∑_(i=1)...This article establishes a universal robust limit theorem under a sublinear expectation framework.Under moment and consistency conditions,we show that,forα∈(1,2),the i.i.d.sequence{(1/√∑_(i=1)^(n)X_(i),1/n∑_(i=1)^(n)X_(i)Y_(i),1/α√n∑_(i=1)^(n)X_(i))}_(n=1)^(∞)converges in distribution to L_(1),where L_(t=(ε_(t),η_(t),ζ_(t))),t∈[0,1],is a multidimensional nonlinear Lévy process with an uncertainty■set as a set of Lévy triplets.This nonlinear Lévy process is characterized by a fully nonlinear and possibly degenerate partial integro-differential equation(PIDE){δ_(t)u(t,x,y,z)-sup_(F_(μ),q,Q)∈■{∫_(R^(d)δλu(t,x,y,z)(dλ)with.To construct the limit process,we develop a novel weak convergence approach based on the notions of tightness and weak compactness on a sublinear expectation space.We further prove a new type of Lévy-Khintchine representation formula to characterize.As a byproduct,we also provide a probabilistic approach to prove the existence of the above fully nonlinear degenerate PIDE.展开更多
A non-local abstract Cauchy problem with a singular integral is studied, which is a closed system of two evolution equations for a real-valued function and a function-valued function. By proposing an appropriate Banac...A non-local abstract Cauchy problem with a singular integral is studied, which is a closed system of two evolution equations for a real-valued function and a function-valued function. By proposing an appropriate Banach space, the well-posedness of the evolution system is proved under some boundedness and smoothness conditions on the coefficient functions. Furthermore, an isomorphism is established to extend the result to a partial integro-differential equation with a singular convolution kernel, which is a generalized form of the stationary Wigner equation. Our investigation considerably improves the understanding of the open problem concerning the well-posedness of the stationary Wigner equation with in ow boundary conditions.展开更多
In this work we investigate the pricing of swing options in a model where the underlying asset follows a jump diffusion process.We focus on the derivation of the partial integro-differential equation(PIDE)which will b...In this work we investigate the pricing of swing options in a model where the underlying asset follows a jump diffusion process.We focus on the derivation of the partial integro-differential equation(PIDE)which will be applied to swing contracts and construct a novel pay-off function from a tree-based pay-off matrix that can be used as initial condition in the PIDE formulation.For valuing swing type derivatives we develop a theta implicit-explicit finite difference scheme to discretize the PIDE using a Gaussian quadrature method for the integral part.Based on known results for the classical theta-method the existence and uniqueness of solution to the new implicit-explicit finite difference method is proven.Various numerical examples illustrate the usability of the proposed method and allow us to analyse the sensitivity of swing options with respect to model parameters.In particular the effects of number of exercise rights,jump intensities and dividend yields will be investigated in depth.展开更多
In this paper,we construct and analyze a Crank-Nicolson fitted finite volume scheme for pricing European options under regime-switching Kou’s jumpdiffusion model which is governed by a system of partial integro-diffe...In this paper,we construct and analyze a Crank-Nicolson fitted finite volume scheme for pricing European options under regime-switching Kou’s jumpdiffusion model which is governed by a system of partial integro-differential equations(PIDEs).We show that this scheme is consistent,stable and monotone as the mesh sizes in space and time approach zero,hence it ensures the convergence to the solution of continuous problem.Finally,numerical experiments are performed to demonstrate the efficiency,accuracy and robustness of the proposed method.展开更多
文摘Differential quadrature method is employed by numerous researchers due to its numerical accuracy and computational efficiency,and is mentioned as potential alternative of conventional numerical methods.In this paper,a differential quadrature based numerical scheme is developed for solving volterra partial integro-differential equation of second order having a weakly singular kernel.The scheme uses cubic trigonometric B-spline functions to determine the weighting coefficients in the differential quadrature approximation of the second order spatial derivative.The advantage of this approximation is that it reduces the problem to a first order time dependent integro-differential equation(IDE).The proposed scheme is obtained in the form of an algebraic system by reducing the time dependent IDE through unconditionally stable Euler backward method as time integrator.The scheme is validated using a homogeneous and two nonhomogeneous test problems.Conditioning of the system matrix and numerical convergence of the method are analyzed for spatial and temporal domain discretization parameters.Comparison of results of the present approach with Sinc collocation method and quasi-wavelet method are also made.
文摘This article studies the development of two numerical techniques for solving convection-diffusion type partial integro-differential equation(PIDE)with a weakly singular kernel.Cubic trigonometric B-spline(CTBS)functions are used for interpolation in both methods.The first method is CTBS based collocation method which reduces the PIDE to an algebraic tridiagonal system of linear equations.The other method is CTBS based differential quadrature method which converts the PIDE to a system of ODEs by computing spatial derivatives as weighted sum of function values.An efficient tridiagonal solver is used for the solution of the linear system obtained in the first method as well as for determination of weighting coefficients in the second method.An explicit scheme is employed as time integrator to solve the system of ODEs obtained in the second method.The methods are tested with three nonhomogeneous problems for their validation.Stability,computational efficiency and numerical convergence of the methods are analyzed.Comparison of errors in approximations produced by the present methods versus different values of discretization parameters and convection-diffusion coefficients are made.Convection and diffusion dominant cases are discussed in terms of Peclet number.The results are also compared with cubic B-spline collocation method.
基金partly supported by SRF for ROCS, SEMsupported by a grant from the "project 211 (phase Ⅲ)" of the Southwestern University of Finance and Economics
文摘This paper develops and analyzes a moving mesh finite difference method for solving partial integro-differential equations. First, the time-dependent mapping of the coordinate transformation is approximated by a a piecewise linear function in time. Then, piecewise quadratic polynomial in space and an efficient method to discretize the memory term of the equation is designed using the moving mesh approach. In each time slice, a simple piecewise constant approximation of the integrand is used, and thus a quadrature is constructed for the memory term. The central finite difference scheme for space and the backward Euler scheme for time are used. The paper proves that the accumulation of the quadrature error is uniformly bounded and that the convergence of the method is second order in space and first order in time. Numerical experiments are carried out to confirm the theoretical predictions.
文摘A framework for the optimal sparse-control of the probability density function of a jump-diffusion process is presented. This framework is based on the partial integro-differential Fokker-Planck (FP) equation that governs the time evolution of the probability density function of this process. In the stochastic process and, correspondingly, in the FP model the control function enters as a time-dependent coefficient. The objectives of the control are to minimize a discrete-in-time, resp. continuous-in-time, tracking functionals and its L2- and L1-costs, where the latter is considered to promote control sparsity. An efficient proximal scheme for solving these optimal control problems is considered. Results of numerical experiments are presented to validate the theoretical results and the computational effectiveness of the proposed control framework.
基金supported by the research grant UL020/08-Y2/MAT/JXQ01/FST,RG063/08-09S/SHW/FST from University of Macao,and the research grant from FDCT of Macao.
文摘We consider pricing options in a jump-diffusion model which requires solving a partial integro-differential equation.Discretizing the spatial direction with a fourth order compact scheme leads to a linear system of ordinary differential equations.For the temporal direction,we utilize the favorable boundary value methods owing to their advantageous stability properties.In addition,the resulting large sparse system can be solved rapidly by the GMRES method with a circulant Strang-type preconditioner.Numerical results demonstrate the high order accuracy of our scheme and the efficiency of the preconditioned GMRES method.
基金supported by the National Key R&D Program of China(Grant No.2018YFA0703900)the National Natural Science Foundation of China(Grant No.11671231)+2 种基金the Qilu Young Scholars Program of Shandong Universitysupported by the Tian Yuan Projection of the National Natural Science Foundation of China(Grant Nos.11526205,11626247)the National Basic Research Program of China(973 Program)(Grant No.2007CB814900(Financial Risk)).
文摘This article establishes a universal robust limit theorem under a sublinear expectation framework.Under moment and consistency conditions,we show that,forα∈(1,2),the i.i.d.sequence{(1/√∑_(i=1)^(n)X_(i),1/n∑_(i=1)^(n)X_(i)Y_(i),1/α√n∑_(i=1)^(n)X_(i))}_(n=1)^(∞)converges in distribution to L_(1),where L_(t=(ε_(t),η_(t),ζ_(t))),t∈[0,1],is a multidimensional nonlinear Lévy process with an uncertainty■set as a set of Lévy triplets.This nonlinear Lévy process is characterized by a fully nonlinear and possibly degenerate partial integro-differential equation(PIDE){δ_(t)u(t,x,y,z)-sup_(F_(μ),q,Q)∈■{∫_(R^(d)δλu(t,x,y,z)(dλ)with.To construct the limit process,we develop a novel weak convergence approach based on the notions of tightness and weak compactness on a sublinear expectation space.We further prove a new type of Lévy-Khintchine representation formula to characterize.As a byproduct,we also provide a probabilistic approach to prove the existence of the above fully nonlinear degenerate PIDE.
基金National Natural Science Foundation of China (Grant Nos. 1167103& 91630130, 91434201, 11421101).
文摘A non-local abstract Cauchy problem with a singular integral is studied, which is a closed system of two evolution equations for a real-valued function and a function-valued function. By proposing an appropriate Banach space, the well-posedness of the evolution system is proved under some boundedness and smoothness conditions on the coefficient functions. Furthermore, an isomorphism is established to extend the result to a partial integro-differential equation with a singular convolution kernel, which is a generalized form of the stationary Wigner equation. Our investigation considerably improves the understanding of the open problem concerning the well-posedness of the stationary Wigner equation with in ow boundary conditions.
文摘In this work we investigate the pricing of swing options in a model where the underlying asset follows a jump diffusion process.We focus on the derivation of the partial integro-differential equation(PIDE)which will be applied to swing contracts and construct a novel pay-off function from a tree-based pay-off matrix that can be used as initial condition in the PIDE formulation.For valuing swing type derivatives we develop a theta implicit-explicit finite difference scheme to discretize the PIDE using a Gaussian quadrature method for the integral part.Based on known results for the classical theta-method the existence and uniqueness of solution to the new implicit-explicit finite difference method is proven.Various numerical examples illustrate the usability of the proposed method and allow us to analyse the sensitivity of swing options with respect to model parameters.In particular the effects of number of exercise rights,jump intensities and dividend yields will be investigated in depth.
基金supported by the National Natural Science Foundation of China(Nos.11971354,and 11701221)the Special Basic Cooperative Research Programs of Yunnan Provincial Undergraduate Universities’Association(No.2019FH001-079)the Fundamental Research Funds for the Central Universities(No.22120210555).
文摘In this paper,we construct and analyze a Crank-Nicolson fitted finite volume scheme for pricing European options under regime-switching Kou’s jumpdiffusion model which is governed by a system of partial integro-differential equations(PIDEs).We show that this scheme is consistent,stable and monotone as the mesh sizes in space and time approach zero,hence it ensures the convergence to the solution of continuous problem.Finally,numerical experiments are performed to demonstrate the efficiency,accuracy and robustness of the proposed method.
