More Robust Pricing of European Options Based on Fourier Cosine Series Expansions

Fabien Le Floc'h

DOI: https://doi.org/10.48550/arXiv.2005.13248

2020-05-27

Computational Finance

Abstract:We present an alternative formula to price European options through cosine series expansions, under models with a known characteristic function such as the Heston stochastic volatility model. It is more robust across strikes and as fast as the original COS method.

What problem does this paper attempt to address?

A Novel Pricing Method for European Options Based on Fourier-Cosine Series Expansions

F. Fang,C. W. Oosterlee

DOI: https://doi.org/10.1137/080718061

IF: 2.968

2009-01-01

SIAM Journal on Scientific Computing

Abstract:Here we develop an option pricing method for European options based on the Fourier-cosine series and call it the COS method. The key insight is in the close relation of the characteristic function with the series coefficients of the Fourier-cosine expansion of the density function. In most cases, the convergence rate of the COS method is exponential and the computational complexity is linear. Its range of application covers underlying asset processes for which the characteristic function is known and various types of option contracts. We will present the method and its applications in two separate parts. The first one is this paper, where we deal with European options in particular. In a follow-up paper we will present its application to options with early-exercise features.

mathematics, applied
On the number of terms in the COS method for European option pricing

Gero Junike

2024-03-31

Abstract:The Fourier-cosine expansion (COS) method is used to price European options numerically in a very efficient way. To apply the COS method, one has to specify two parameters: a truncation range for the density of the log-returns and a number of terms N to approximate the truncated density by a cosine series. How to choose the truncation range is already known. Here, we are able to find an explicit and useful bound for N as well for pricing and for the sensitivities, i.e., the Greeks Delta and Gamma, provided the density of the log-returns is smooth. We further show that the COS method has an exponential order of convergence when the density is smooth and decays exponentially. However, when the density is smooth and has heavy tails, as in the Finite Moment Log Stable model, the COS method does not have exponential order of convergence. Numerical experiments confirm the theoretical results.

Computational Finance
The double exponential jump diffusion model for pricing European options under fuzzy environments

Lihua Zhang,Wei Zhang,Weijun Xu,Weilin Xiao

DOI: https://doi.org/10.1016/j.econmod.2012.02.005

IF: 4.7

2012-01-01

Economic Modelling

Abstract:The interest rate and volatility may have different values in the different commercial banks and financial institutions. Moreover, the fluctuations of the underlying assets are rare events, and there are not enough historical data to estimate the jump intensity in a precise sense. This paper considers European option pricing problems with the fuzzy interest rate, fuzzy drift, fuzzy volatility and fuzzy jump intensity. We present the fuzzy pricing formula of European options based on the Kou's double exponential jump diffusion model. We also obtain the crisp possibilistic mean option pricing formula in fuzzy double exponential jump diffusion model by using the crisp possibilistic mean values of the fuzzy numbers. Comparing with B-S formula, numerical analysis and empirical results show that the fuzzy double exponential jump diffusion formula and the crisp possibilistic mean option pricing formula are reasonable and can be taken as reference pricing tools for the financial investors.
Pricing currency options in a fractional Brownian motion with jumps

Wei-Lin Xiao,Wei-Guo Zhang,Xi-Li Zhang,Ying-Luo Wang

DOI: https://doi.org/10.1016/j.econmod.2010.05.010

IF: 4.7

2010-01-01

Economic Modelling

Abstract:A new framework for pricing the European currency option is developed in the case where the spot exchange rate fellows a fractional Brownian motion with jumps. An analytic formula for pricing European foreign currency options is proposed using the equivalent martingale measure and the estimation method of parameters in the pricing model is given, enabling option prices to be computed efficiently and accurately. For the purpose of understanding the pricing model, some properties of this pricing model are discussed in the latter part of this paper. Finally, the numerical simulations illustrate that our model is flexible and easy to implement.
Option pricing with Legendre polynomials

Julien Hok,Tat Lung Chan

DOI: https://doi.org/10.48550/arXiv.1610.03086

2017-03-19

Abstract:Here we develop an option pricing method based on Legendre series expansion of the density function. The key insight, relying on the close relation of the characteristic function with the series coefficients, allows to recover the density function rapidly and accurately. Based on this representation for the density function, approximations formulas for pricing European type options are derived. To obtain highly accurate result for European call option, the implementation involves integrating high degree Legendre polynomials against exponential function. Some numerical instabilities arise because of serious subtractive cancellations in its formulation (96) in proposition 7.1. To overcome this difficulty, we rewrite this quantity as solution of a second-order linear difference equation and solve it using a robust and stable algorithm from Olver. Derivation of the pricing method has been accompanied by an error analysis. Errors bounds have been derived and the study relies more on smoothness properties which are not provided by the payoff? functions, but rather by the density function of the underlying stochastic models. This is particularly relevant for options pricing where the payoff of the contract are generally not smooth functions. The numerical experiments on a class of models widely used in quantitative finance show exponential convergence.

Mathematical Finance,Computational Finance
From characteristic functions to multivariate distribution functions and European option prices by the damped COS method

Gero Junike,Hauke Stier

2024-06-02

Abstract:We provide a unified framework to obtain numerically certain quantities, such as the distribution function, absolute moments and prices of financial options, from the characteristic function of some (unknown) probability density function using the Fourier-cosine expansion (COS) method. The classical COS method is numerically very efficient in one-dimension, but it cannot deal very well with certain integrands in general dimensions. Therefore, we introduce the damped COS method, which can handle a large class of integrands very efficiently. We prove the convergence of the (damped) COS method and study its order of convergence. The method converges exponentially if the characteristic function decays exponentially. To apply the (damped) COS method, one has to specify two parameters: a truncation range for the multivariate density and the number of terms to approximate the truncated density by a cosine series. We provide an explicit formula for the truncation range and an implicit formula for the number of terms. Numerical experiments up to five dimensions confirm the theoretical results.

Computational Finance
An efficient pricing method for rainbow options based on two-dimensional modified sine–sine series expansions

Qing-Jiang Meng,Deng Ding

DOI: https://doi.org/10.1080/00207160.2012.749349

IF: 1.75

2013-02-14

International Journal of Computer Mathematics

Abstract:An efficient option pricing method based on the Fourier cosine series expansion was proposed by Fang and Oosterlee [A novel option pricing method based on Fourier-cosine series expansions, SIAM J. Sci. Comput. 31 (2008), pp. 826–848], and recently it has been used in two-dimensional cases by Ruijter and Oosterlee [Two-dimensional Fourier cosine series expansion method for pricing financial options, SIAM J. Sci. Comput. 34 (2012), pp. B642–B671]. In this paper, we consider a modification of the Fourier sine–sine expansion, whereby in [−1, 1] the sin(π nx) functions are replaced by sin(π(n− )x), n≥1, for pricing rainbow options. Feasible approaches for obtaining the pricing formulae are presented. A numerical comparison and a flexible and general error analysis are considered to prove the effectiveness of the method. Practical applications to basket options and correlation options under different models are also given. They indicate satisfactory convergence and efficiency as expected.

mathematics, applied
Nonparametric estimates of option prices via Hermite basis functions

Carlo Marinelli,Stefano d'Addona

DOI: https://doi.org/10.1007/s10436-023-00431-4

2023-08-12

Abstract:We consider approximate pricing formulas for European options based on approximating the logarithmic return's density of the underlying by a linear combination of rescaled Hermite polynomials. The resulting models, that can be seen as perturbations of the classical Black-Scholes one, are nonpararametric in the sense that the distribution of logarithmic returns at fixed times to maturity is only assumed to have a square-integrable density. We extensively investigate the empirical performance, defined in terms of out-of-sample relative pricing error, of this class of approximating models, depending on their order (that is, roughly speaking, the degree of the polynomial expansion) as well as on several ways to calibrate them to observed data. Empirical results suggest that such approximate pricing formulas, when compared with simple nonparametric estimates based on interpolation and extrapolation on the implied volatility curve, perform reasonably well only for options with strike price not too far apart from the strike prices of the observed sample.

Pricing of Securities
Hedging and Pricing European-type, Early-Exercise and Discrete Barrier Options using Algorithm for the Convolution of Legendre Series

Tat Lung Chan,Nicholas Hale

DOI: https://doi.org/10.48550/arXiv.1811.09257

2018-11-22

Computational Finance

Abstract:This paper applies an algorithm for the convolution of compactly supported Legendre series (the CONLeg method) (cf. Hale and Townsend 2014a), to pricing/hedging European-type, early-exercise and discrete-monitored barrier options under a Levy process. The paper employs Chebfun (cf. Trefethen et al. 2014) in computational finance and provides a quadrature-free approach by applying the Chebyshev series in financial modelling. A significant advantage of using the CONLeg method is to formulate option pricing and option Greek curves rather than individual prices/values. Moreover, the CONLeg method can yield high accuracy in option pricing and hedging when the risk-free smooth probability density function (PDF) is smooth/non-smooth. Finally, we show that our method can accurately price/hedge options deep in/out of the money and with very long/short maturities. Compared with existing techniques, the CONLeg method performs either favourably or comparably in numerical experiments.
Generalized Finite Integration Method with Laplace transform for European option pricing under Black–Scholes and Heston models

Y. Ma,C.Z. Shi,Y.C. Hon

DOI: https://doi.org/10.1016/j.enganabound.2024.105751

IF: 3.25

2024-05-06

Engineering Analysis with Boundary Elements

Abstract:In this paper, we combine a recently developed Generalized Finite Integration Method (GFIM) with Laplace transform technique for pricing options under the Black Scholes model and Heston model respectively. Instead of using traditional time-stepping process, we first perform Laplace transform on the governing equation and boundary conditions to remove the temporal derivatives. The Generalized Finite Integration Method is then exploited to handle the spatial differential operators in the transformed space. From numerical Laplace inversion algorithm, we restore the required time-dependent option price. For verification, several option pricing models governed by one-dimensional Black–Scholes equation and two-dimensional extended Heston equation are constructed to demonstrate the efficiency and feasibility of the proposed approach.

engineering, multidisciplinary,mathematics, interdisciplinary applications
On bid and ask pricing of European options via direct discretization of Choquet distorted expectations

Matteo Michielon

DOI: https://doi.org/10.1080/14697688.2024.2430309

IF: 1.986

2024-12-05

Quantitative Finance

Abstract:We present a numerical approach to compute bid and ask prices of European options within the conic finance paradigm. In particular, we show how discretizing option bid and ask pricing formulae, where the latter are expressed in terms of Choquet integrals, allows to perform efficient pricing and calibration by means of the COS method and the fast Fourier transform technique. Additionally, the approach introduced provides alternative expressions for risk-neutral European option prices by leveraging the relationship between Choquet and Riemann-Stieltjes integrals. These alternative formulations, requiring only the discretization of the cumulative distribution function of the (transformed) terminal value of the underlying asset, enhance computational efficiency under risk-neutral settings. We provide practical illustrations of the methodology when the underlying asset process follows Heston or Variance Gamma dynamics for both pricing and implied liquidity calculations.

economics,social sciences, mathematical methods,mathematics, interdisciplinary applications,business, finance
The SINC way: A fast and accurate approach to Fourier pricing

Fabio Baschetti,Giacomo Bormetti,Silvia Romagnoli,Pietro Rossi

DOI: https://doi.org/10.48550/arXiv.2009.00557

2021-05-19

Abstract:The goal of this paper is to investigate the method outlined by one of us (PR) in Cherubini et al. (2009) to compute option prices. We name it the SINC approach. While the COS method by Fang and Osterlee (2009) leverages the Fourier-cosine expansion of truncated densities, the SINC approach builds on the Shannon Sampling Theorem revisited for functions with bounded support. We provide several results which were missing in the early derivation: i) a rigorous proof of the convergence of the SINC formula to the correct option price when the support grows and the number of Fourier frequencies increases; ii) ready to implement formulas for put, Cash-or-Nothing, and Asset-or-Nothing options; iii) a systematic comparison with the COS formula for several log-price models; iv) a numerical challenge against alternative Fast Fourier specifications, such as Carr and Madan (1999) and Lewis (2000); v) an extensive pricing exercise under the rough Heston model of Jaisson and Rosenbaum (2015); vi) formulas to evaluate numerically the moments of a truncated density. The advantages of the SINC approach are numerous. When compared to benchmark methodologies, SINC provides the most accurate and fast pricing computation. The method naturally lends itself to price all options in a smile concurrently by means of Fast Fourier techniques, boosting fast calibration. Pricing requires to resort only to odd moments in the Fourier space. A previous version of this manuscript circulated with the title `Rough Heston: The SINC way'.

Pricing of Securities,Computational Finance
Pricing European Options in Realistic Markets

Martin Schaden

DOI: https://doi.org/10.48550/arXiv.physics/0210025

2002-10-06

Physics and Society

Abstract:We investigate the relation between the fair price for European-style vanilla options and the distribution of short-term returns on the underlying asset ignoring transaction and other costs. We compute the risk-neutral probability density conditional on the total variance of the asset's returns when the option expires. If the asset's future price has finite expectation, the option's fair value satisfies a parabolic partial differential equation of the Black-Scholes type in which the variance of the asset's returns rather than a trading time is the evolution parameter. By immunizing the portfolio against large-scale price fluctuations of the asset, the valuation of options is extended to the realistic case\cite{St99} of assets whose short-term returns have finite variance but very large, or even infinite, higher moments. A dynamic Delta-hedged portfolio that is statically insured against exceptionally large fluctuations includes at least two different options on the asset. The fair value of an option in this case is determined by a universal drift function that is common to all options on the asset. This drift is interpreted as the premium for an investment exposed to risk due to exceptionally large variations of the asset's price. It affects the option valuation like an effective cost-of-carry for the underlying in the Black-Scholes world would. The derived pricing formula for options in realistic markets is arbitrage free by construction. A simple model with constant drift qualitatively reproduces the often observed volatility -skew and -term structure.
Stochastic Expansion for the Pricing of Asian and Basket Options

Fabien Le Floc'h

2024-08-11

Abstract:We present closed analytical approximations for the pricing of basket options, also applicable to Asian options with discrete averaging under the Black-Scholes model with time-dependent parameters. The formulae are obtained by using a stochastic Taylor expansion around a log-normal proxy model and are found to be highly accurate for Asian options in practice as well as for vanilla options with discrete dividends.

Pricing of Securities,Computational Finance,Mathematical Finance
A Novel Fourth-Order Finite Difference Scheme for European Option Pricing in the Time-Fractional Black–Scholes Model

Xin Cai,Yihong Wang

DOI: https://doi.org/10.3390/math12213343

IF: 2.4

2024-10-25

Mathematics

Abstract:This paper addresses the valuation of European options, which involves the complex and unpredictable dynamics of fractal market fluctuations. These are modeled using the α-order time-fractional Black–Scholes equation, where the Caputo fractional derivative is applied with the parameter α ranging from 0 to 1. We introduce a novel, high-order numerical scheme specifically crafted to efficiently tackle the time-fractional Black–Scholes equation. The spatial discretization is handled by a tailored finite point scheme that leverages exponential basis functions, complemented by an L1-discretization technique for temporal progression. We have conducted a thorough investigation into the stability and convergence of our approach, confirming its unconditional stability and fourth-order spatial accuracy, along with (2−α)-order temporal accuracy. To substantiate our theoretical results and showcase the precision of our method, we present numerical examples that include solutions with known exact values. We then apply our methodology to price three types of European options within the framework of the time-fractional Black–Scholes model: (i) a European double barrier knock-out call option; (ii) a standard European call option; and (iii) a European put option. These case studies not only enhance our comprehension of the fractional derivative's order on option pricing but also stimulate discussion on how different model parameters affect option values within the fractional framework.

mathematics
European Option Pricing With a Fast Fourier Transform Algorithm for Big Data Analysis

Shuang Xiao,Shi-Hua Ma,Guo Li,Samar K. Mukhopadhyay

DOI: https://doi.org/10.1109/tii.2015.2500885

IF: 12.3

2016-06-01

IEEE Transactions on Industrial Informatics

Abstract:Several empirical studies show that, under multiple risks, markets exhibit many new properties, such as volatility smile and cluster fueled by the explosion of transaction data. This paper attempts to capture these newly developed features using the valuation of European options as a vehicle. Statistical analysis performed on the data collected from the currency option market clearly shows the coexistence of mean reversion, jumps, volatility smile, and leptokurtosis and fat tail. We characterize the dynamics of the underlying asset in this kind of environment by establishing a coupled stochastic differential equation model with triple characteristics of mean reversion, nonaffine stochastic volatility, and mixed-exponential jumps. Moreover, we propose a characteristic function method to derive the closed-form pricing formula. We also present a fast Fourier transform (FFT) algorithm-based numerical solution method. Finally, extensive numerical experiments are conducted to validate both the modeling methodology and the numerical algorithm. Results demonstrate that the model behaves well in capturing the properties observed in the market, and the FFT numerical algorithm is both accurate and efficient in addressing large amount of data.

automation & control systems,computer science, interdisciplinary applications,engineering, industrial
Efficient pricing of path‐dependent interest rate derivatives

Allan Jonathan da Silva,Jack Baczynski,José V. M. Vicente

DOI: https://doi.org/10.1002/asmb.2877

2024-06-02

Applied Stochastic Models in Business and Industry

Abstract:Interest rate derivative pricing is a critical aspect of fixed‐income markets, where efficient methods are essential. This study introduces a novel approach to pricing path‐dependent interest rate derivatives within a broad class of affine jumps. The study's particular setting is the Fourier‐cosine series (COS) method adaptation, which offers an accurate and computationally efficient method for pricing interest rate derivatives. The Fourier‐cosine series approach can be used to compute probability density functions and option pricing with a linear computing complexity and exponential convergence rate. The lack of a quick and precise pricing technique for Asian interest rate options in diverse fixed‐income market scenarios is a research gap that is being addressed. This approach closes this gap by providing quasi‐closed and closed‐form equations for a range of density and characteristic functions, resulting in precise pricing. The results demonstrate the versatility of the COS method in interest rate markets. Similar to what has been previously reported for stock options, the numerical findings demonstrate the extreme precision and computing speed of the pricing and hedging estimations provided here. This method is an innovative approach to interest rate derivative pricing, offering researchers and practitioners a powerful tool for efficiently calculating prices and calibrating options across strikes and maturities.

statistics & probability,mathematics, interdisciplinary applications,operations research & management science
A fractional Black-Scholes model with stochastic volatility and European option pricing

Xin-Jiang He,Sha Lin

DOI: https://doi.org/10.1016/j.eswa.2021.114983

IF: 8.5

2021-09-01

Expert Systems with Applications

Abstract:<p>In this paper, we introduce the stochastic volatility into the FMLS (finite moment log-stable) model to capture the effect of both jumps and stochastic volatility. However, this additional stochastic source adds another degree of complexity in seeking for analytical formula when pricing European options, as the involved FPDE (fractional partial differential equation) system governing option prices is now of three dimensions. Albeit difficult, we have still managed to present an analytical solution expressed in terms of Fourier cosine series, after a two-step solution procedure is developed for the target FPDE system. This solution is different from the most literature as it is truly explicit, involving no Fourier inversion. It is also shown through the numerical experiments that it converges very rapidly and has potential to be applied in practice.</p>

computer science, artificial intelligence,engineering, electrical & electronic,operations research & management science
A closed-form pricing formula for European options with market liquidity risk

Puneet Pasricha,Song-Ping Zhu,Xin-Jiang He

DOI: https://doi.org/10.1016/j.eswa.2021.116128

IF: 8.5

2022-03-01

Expert Systems with Applications

Abstract:In this paper, the impact of liquidity on the underlying asset is taken into account when pricing European options through a discounting factor which depends on two factors, i.e., market liquidity risk modeled as a mean-reverting stochastic process and the sensitivity of the underlying to market liquidity. A closed-form pricing formula for liquidity-adjusted European options is derived in the form of an infinite series using Karhunen–Loève expansion for the Ornstein–Uhlenbeck process. The convergence of the series solution is theoretically proved to guarantee closedness so that market practitioners can adopt the new formula when they need to account in market liquidity risk. The speed of convergence is demonstrated through numerical experiments. Finally, the accuracy of the newly derived formula is shown by comparing option prices calculated with our formula and those obtained from Monte-Carlo simulation, and various properties of our formula are also investigated.

computer science, artificial intelligence,engineering, electrical & electronic,operations research & management science

More Robust Pricing of European Options Based on Fourier Cosine Series Expansions

A Novel Pricing Method for European Options Based on Fourier-Cosine Series Expansions

On the number of terms in the COS method for European option pricing

The double exponential jump diffusion model for pricing European options under fuzzy environments

Pricing currency options in a fractional Brownian motion with jumps

Option pricing with Legendre polynomials

From characteristic functions to multivariate distribution functions and European option prices by the damped COS method

An efficient pricing method for rainbow options based on two-dimensional modified sine–sine series expansions

Nonparametric estimates of option prices via Hermite basis functions

Hedging and Pricing European-type, Early-Exercise and Discrete Barrier Options using Algorithm for the Convolution of Legendre Series

Generalized Finite Integration Method with Laplace transform for European option pricing under Black–Scholes and Heston models

On bid and ask pricing of European options via direct discretization of Choquet distorted expectations

The SINC way: A fast and accurate approach to Fourier pricing

Pricing European Options in Realistic Markets

Stochastic Expansion for the Pricing of Asian and Basket Options

A Novel Fourth-Order Finite Difference Scheme for European Option Pricing in the Time-Fractional Black–Scholes Model

European Option Pricing With a Fast Fourier Transform Algorithm for Big Data Analysis

Efficient pricing of path‐dependent interest rate derivatives

A fractional Black-Scholes model with stochastic volatility and European option pricing

A closed-form pricing formula for European options with market liquidity risk