The value of a European option satisfies the Black-Scholes equation with appropriately specified final and boundary conditions.We transform the problem to an initial boundary value problem in dimensionless form.There ...The value of a European option satisfies the Black-Scholes equation with appropriately specified final and boundary conditions.We transform the problem to an initial boundary value problem in dimensionless form.There are two parameters in the coefficients of the resulting linear parabolic partial differential equation.For a range of values of these parameters,the solution of the problem has a boundary or an initial layer.The initial function has a discontinuity in the first-order derivative,which leads to the appearance of an interior layer.We construct analytically the asymptotic solution of the equation in a finite domain.Based on the asymptotic solution we can determine the size of the artificial boundary such that the required solution in a finite domain in x and at the final time is not affected by the boundary.Also,we study computationally the behaviour in the maximum norm of the errors in numerical solutions in cases such that one of the parameters varies from finite (or pretty large) to small values,while the other parameter is fixed and takes either finite (or pretty large) or small values. Crank-Nicolson explicit and implicit schemes using centered or upwind approximations to the derivative are studied.We present numerical computations,which determine experimentally the parameter-uniform rates of convergence.We note that this rate is rather weak,due probably to mixed sources of error such as initial and boundary layers and the discontinuity in the derivative of the solution.展开更多
We introduce a weak Galerkin finite element method for the valuation of American options governed by the Black-Scholes equation. In order to implement, we need to solve the optimal exercise boundary and then introduce...We introduce a weak Galerkin finite element method for the valuation of American options governed by the Black-Scholes equation. In order to implement, we need to solve the optimal exercise boundary and then introduce an artificial boundary to make the computational domain bounded. For the optimal exercise boundary, which satisfies a nonlinear Volterra integral equation, it is resolved by a higher-order collocation method based on graded meshes. With the computed optimal exercise boundary, the front-fixing technique is employed to transform the free boundary problem to a one- dimensional parabolic problem in a half infinite area. For the other spatial domain boundary, a perfectly matched layer is used to truncate the unbounded domain and carry out the computation. Finally, the resulting initial-boundary value problems are solved by weak Galerkin finite element method, and numerical examples are provided to illustrate the efficiency of the method.展开更多
In this paper we are concerned with the pricing of lookback options with American type constrains. Based on the differential linear complementary formula associated with the pricing problem, an implicit difference sch...In this paper we are concerned with the pricing of lookback options with American type constrains. Based on the differential linear complementary formula associated with the pricing problem, an implicit difference scheme is constructed and analyzed. We show that there exists a unique difference solution which is unconditionally stable. Using the notion of viscosity solutions, we also prove that the finite difference solution converges uniformly to the viscosity solution of the continuous problem. Furthermore, by means of the variational inequality analysis method, the O(△t + △x^2)-order error estimate is derived in the discrete L2-norm provided that the continuous problem is sufficiently regular. In addition, a numerical example is provided to illustrate the theoretical results.展开更多
基金partially supported by the Russian Foundation for Basic Research(grant No.07-01-00729)the Singapore Academic Research Funds R-146-000-064-112 and R-146-000-099-112the Boole Centre for Research in Informatics at the National University of Ireland,Cork and by the Mathematics Applications Consortium for Science and Industry in Ireland(MACSI)under the Science Foundation Ireland Mathematics Initiative.
文摘The value of a European option satisfies the Black-Scholes equation with appropriately specified final and boundary conditions.We transform the problem to an initial boundary value problem in dimensionless form.There are two parameters in the coefficients of the resulting linear parabolic partial differential equation.For a range of values of these parameters,the solution of the problem has a boundary or an initial layer.The initial function has a discontinuity in the first-order derivative,which leads to the appearance of an interior layer.We construct analytically the asymptotic solution of the equation in a finite domain.Based on the asymptotic solution we can determine the size of the artificial boundary such that the required solution in a finite domain in x and at the final time is not affected by the boundary.Also,we study computationally the behaviour in the maximum norm of the errors in numerical solutions in cases such that one of the parameters varies from finite (or pretty large) to small values,while the other parameter is fixed and takes either finite (or pretty large) or small values. Crank-Nicolson explicit and implicit schemes using centered or upwind approximations to the derivative are studied.We present numerical computations,which determine experimentally the parameter-uniform rates of convergence.We note that this rate is rather weak,due probably to mixed sources of error such as initial and boundary layers and the discontinuity in the derivative of the solution.
基金Acknowledgements The authors would like to thank the anonymous referees for their careflll reading of the manuscript and their valuable comments. The authors also wish to thank the High Performance Computing Center of Jilin University and C, omputing Center of ,lilin Province for essential support. This work was supported by the National Natural Science Foundation of China Grant Nos. 11271157, 11371171), the Open Project Program of the State Key Lab of CAD&CG (A1302) of Zhejiang University, the Scientific Research Foundation for bleturned Scholars, Ministry of Education of China. and UIBE (11QD17).
文摘We introduce a weak Galerkin finite element method for the valuation of American options governed by the Black-Scholes equation. In order to implement, we need to solve the optimal exercise boundary and then introduce an artificial boundary to make the computational domain bounded. For the optimal exercise boundary, which satisfies a nonlinear Volterra integral equation, it is resolved by a higher-order collocation method based on graded meshes. With the computed optimal exercise boundary, the front-fixing technique is employed to transform the free boundary problem to a one- dimensional parabolic problem in a half infinite area. For the other spatial domain boundary, a perfectly matched layer is used to truncate the unbounded domain and carry out the computation. Finally, the resulting initial-boundary value problems are solved by weak Galerkin finite element method, and numerical examples are provided to illustrate the efficiency of the method.
基金supported in part by the National Basic Research Program(2007CB814906)the National Natural Science Foundation of China(10771031,10471019,10471103,and 10771158)+1 种基金Social Science Foundation of the Ministry of Education of China(Numerical methods for convertible bonds,06JA630047)Tianjin Natural Science Foundation(07JCYBJC14300)and Tianjin University of Finance and Economics
文摘In this paper we are concerned with the pricing of lookback options with American type constrains. Based on the differential linear complementary formula associated with the pricing problem, an implicit difference scheme is constructed and analyzed. We show that there exists a unique difference solution which is unconditionally stable. Using the notion of viscosity solutions, we also prove that the finite difference solution converges uniformly to the viscosity solution of the continuous problem. Furthermore, by means of the variational inequality analysis method, the O(△t + △x^2)-order error estimate is derived in the discrete L2-norm provided that the continuous problem is sufficiently regular. In addition, a numerical example is provided to illustrate the theoretical results.