A fast and accurate numerical method for the symmetric Lévy processes based on the Fourier transform

来源 :第八届工业与应用数学国际大会 | 被引量 : 0次 | 上传用户:dongfa
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  We propose a fast and accurate numerical method based on Fourier transform to solve Kolmogorov forward equations of symmetric scalar Lévy processes.The method is based on the accurate numerical formulas for Fourier transform proposed by Ooura.These formulas are combined with nonuniform fast Fourier transform(FFT)and fractional FFT to speed up the numerical computations.Moreover,we propose a formula for numerical indefinite integration on equispaced grids as a component of the method.The proposed integration formula is based on the sinc-Gauss sampling formula,which is a function approximation formula.This integration formula is also combined with the FFT.Therefore,all steps of the proposed method are executed using the FFT and its variants.The proposed method realizes exponential convergence and allows us to be free from some special treatments for a non-smooth initial condition and numerical time integration.
其他文献
We construct a new monotone finite volume method for diffusion equations on star-shaped polygonal meshes.A distinct feature of the new scheme is that the discrete stencil of normal flux on a cell-edge
会议
Ramsey graphs theory deals with regularity and coloring of graphs.There are many interesting applications of Ramsey graphs theory,such as in the fields of communications,information retrieval,and deci
会议
Let F,G,and H are simple graphs.We say F→(G,H)if for every 2-coloring of the edges of F there exist a monochromatic G or H in F.The Ramsey number r(G,H)is defined as min V(F)— F→(G,H)and the restricte
会议
By using the finite projective geometry method,the separating properties of linear constant-weight codes are presented.An algorithm is given for computing the cardinality of separating coordinate posi
会议
Modular homogeneous two-weight codes are useful in graph and combination theory,and the graph of a modular homogeneous two-weight code is strongly regular.We will construct several classes of modular
会议
Fractional power series solutions for nonlinear fractional ordinary differential equations are considered by using a fast and efficient algorithm for the Adomian polynomials.
会议
Covariance(ensemble)inflation is a process typically required in ensemble data assimilation,where an algorithm is used to periodically increase the ensemble variance in order to prevent ensemble colla
会议
A quaternion rational surface is a surface generated from two rational space curves by quaternion multiplication.The goal of this talk is to demonstrate how to apply syzygies to analyze quaternion rat
会议
A second order accurate finite difference scheme is presented for the 2-D and 3-D Cahn-Hilliard equation,and an error analysis with an improved convergence constant is provided.The unique solvability
会议
We consider two bio-economic ratio-dependent predator-prey models for prey harvesting,one without and the other with additional food supply to the predators.We analyze the stability of the equilibrium
会议