运用Matlab基于LSM方法对羹式期权定价的新探究

传统期权定价方法是通过主观假定初始价格、执行价格、期限、波动率、无风险利率等条件来对期权进行定价,很少联系实际的期权市场报价对期权进行定价。本文根据股票期权市场报价,通过Matlab快速方便地求解出隐含的波动率和无风险利率,并在此基础上运用Matlab基于最／bZ．乘蒙特卡洛模拟（LSM）方法对该股票的美式期权进行定价。本文揭示了如何根据期权市场报价实现隐含波动率和无风险利率的求解,进而结合LSM方法对美式期权进行定价的一种新方法。此外,本文对LSM方法的改进技术也进行了探讨。     

The representation of American options prices under stochastic volatility and jump-diffusion dynamics

This paper considers the problem of pricing American options when the dynamics of the underlying are driven by both stochastic volatility following a square-root process as used by Heston [Rev. Financial Stud., 1993, 6, 327–343], and by a Poisson jump process as introduced by Merton [J. Financial Econ., 1976, 3, 125–144]. Probability arguments are invoked to find a representation of the solution in terms of expectations over the joint distribution of the underlying process. A combination of Fourier transform in the log stock price and Laplace transform in the volatility is then applied to find the transition probability density function of the underlying process. It turns out that the price is given by an integral dependent upon the early exercise surface, for which a corresponding integral equation is obtained. The solution generalizes in an intuitive way the structure of the solution to the corresponding European option pricing problem obtained by Scott [Math. Finance, 1997, 7(4), 413–426], but here in the case of a call option and constant interest rates.     

American options and callable bonds under stochastic interest rates and endogenous bankruptcy

A new characterization of the American-style option is proposed under a very general multifactor Markovian and diffusion framework. The efficiency of the proposed pricing solutions is shown to depend only on the use of a viable valuation method for the corresponding European-style option and for the transition density of the model’s state variables. Under a Gauss-Markov stochastic interest rates setup, these new American option pricing solutions are shown to offer a much better accuracy-efficiency trade-off than the approximations already available in the literature. This result is also used to price callable corporate bonds under an endogenous bankruptcy structural approach, by decomposing the option to call or default into a European put on the firm value plus two early exercise premium components.     

Lapse rate modeling: a rational expectation approach

The surrender option embedded in many life insurance products is a clause that allows policyholders to terminate the contract early. Pricing techniques based on the American Contingent Claim (ACC) theory are often used, though the actual policyholders' behavior is far from optimal. Inspired by many prepayment models for mortgage backed securities, this paper builds a Rational Expectation (RE) model describing the policyholders' behavior in lapsing the contract. A market model with stochastic interest rates is considered, and the pricing is carried out through numerical approximation of the corresponding two-space-dimensional parabolic partial differential equation. Extensive numerical experiments show the differences in terms of pricing and interest rate elasticity between the ACC and RE approaches as well as the sensitivity of the contract price with respect to changes in the policyholders' behavior.     

Using Richardson extrapolation techniques to price American options with alternative stochastic processes

In this paper the authors investigate the performance of the original and repeated Richardson extrapolation methods for American option pricing by implementing both the original and modified Geske?CJohnson approximation formulae. A comprehensive numerical comparison includes alternative stochastic processes of the underlying asset price. The numerical results show that whether the original or modified formula is implemented, the Richardson extrapolation techniques work very well. The repeated Richardson extrapolation strongly outperforms the original, especially when the underlying asset price follows a stochastic volatility process. Moreover, this study verifies the feasibility of the estimated error bounds of the American option prices under alternative stochastic processes by applying the repeated Richardson extrapolation method and estimating the interval of true American option values, as well as determining the number of options needed for an approximation to achieve a desired accuracy level.     

An exact and explicit solution for the valuation of American put options

In this paper, an exact and explicit solution of the well-known Black–Scholes equation for the valuation of American put options is presented for the first time. To the best of the author's knowledge, a closed-form analytical formula has never been found for the valuation of American options of finite maturity, although there have been quite a few approximate solutions and numerical approaches proposed. The closed-form exact solution presented here is written in the form of a Taylor's series expansion, which contains infinitely many terms. However, only about 30 terms are actually needed to generate a convergent numerical solution if the solution of the corresponding European option is taken as the initial guess of the solution series. The optimal exercise boundary, which is the main difficulty of the problem, is found as an explicit function of the risk-free interest rate, the volatility and the time to expiration. A key feature of our solution procedure, which is based on the homotopy-analysis method, is the optimal exercise boundary being elegantly and temporarily removed in the solution process of each order, and, consequently, the solution of a linear problem can be analytically worked out at each order, resulting in a completely analytical and exact series-expansion solution for the optimal exercise boundary and the option price of American put options.     

American option pricing under the double Heston model based on asymptotic expansion

This paper focuses on pricing American put options under the double Heston model proposed by Christoffersen et al. By introducing an explicit exercise rule, we obtain the asymptotic expansion of the solution to the partial differential equation for pricing American put options. We calculate American option price by the sum of the European option price and the early exercise premium. The early exercise premium is calculated by the difference between the American and European option prices based on asymptotic expansions. The European option price is obtained by the efficient COS method. Based on the obtained American option price, the double Heston model is calibrated by minimizing the distance between model and market prices, which yields an optimization problem that is solved by a differential evolution algorithm combined with the Matlab function fmincon.m. Numerical results show that the pricing approach is fast and accurate. Empirical results show that the double Heston model has better performance in pricing short-maturity American put options and capturing the volatility term structure of American put options than the Heston model.     

Pricing Options under Stochastic Interest Rates: A New Approach

We will generalize the Black-Scholes option pricing formula by incorporating stochastic interest rates. Although the existing literature has obtained some formulae for stock options under stochastic interest rates, the closed-form solutions have been known only under the Gaussian (Merton type) interest rate processes. We will show that an explicit solution, which is an extended Black-Scholes formula under stochastic interest rates in certain asymptotic sense, can be obtained by extending the asymptotic expansion approach when the interest rate volatility is small. This method, called the small-disturbance asymptotics for Itô processes, has recently been developed by Kunitomo and Takahashi (1995, 1998) and Takahashi (1997). We found that the extended Black-Scholes formula is decomposed into the original Black-Scholes formula under the deterministic interest rates and the adjustment term driven by the volatility of interest rates. We will illustrate the numerical accuracy of our new formula by using the Cox–Ingersoll–Ross model for the interest rates.     

A generalization of reset options pricing formulae with stochastic interest rates

A generalization of reset call options with predetermined dates is derived in the case of time-dependent volatility and time-dependent interest rate by applying martingale method and change of nume?aire or change of probability measure. An analytical pricing formula for the reset call option is also obtained when the interest rate follows an extended Vasicek’s model. Numerical results show that the correlated coefficient between the stock price and interest rate is almost unacted on the price of reset call option with short maturity and Monte Carlo method is inefficient. Monte Carlo method should be only used if there is no closed-formed solution for option pricing.     

A perturbative approach to Bermudan options pricing with applications

In this paper we address the problem of the valuation of Bermudan option derivatives in the framework of multi-factor interest rate models. We propose a solution in which the exercise decision entails a properly defined series expansion. The method allows for the fast computation of both a lower and an upper bound for the option price, and a tight control of its accuracy, for a generic Markovian interest rate model. In particular, we show detailed computations in the case of the Bond Market Model. As examples we consider the case of a zero coupon Bermudan option and a coupon bearing Bermudan option; in order to demonstrate the wide applicability of the proposed methodology we also consider the case of a last generation payoff, a Bermudan option on a CMS spread bond.     

A quasi-analytical interpolation method for pricing American options under general multi-dimensional diffusion processes

We present a quasi-analytical method for pricing multi-dimensional American options based on interpolating two arbitrage bounds, along the lines of Johnson in J Financ Quant Anal 18(1):141–148 (1983). Our method allows for the close examination of the interpolation parameter on a rigorous theoretical footing instead of empirical regression. The method can be adapted to general diffusion processes as long as quick and accurate pricing methods exist for the corresponding European and perpetual American options. The American option price is shown to be approximately equal to an interpolation of two European option prices with the interpolation weight proportional to a perpetual American option. In the Black-Scholes model, our method achieves the same efficiency as the quadratic approximation of Barone-Adesi and Whaley in J Financ 42:301–320 (1987), with our method being generally more accurate for out-of-the-money and long-maturity options. When applied to Heston’s stochastic volatility model, our method is shown to be extremely efficient and fairly accurate.     

Adding and subtracting Black-Scholes: A new approach to approximating derivative prices in continuous-time models

We develop a new approach to approximating asset prices in the context of continuous-time models. For any pricing model that lacks a closed-form solution, we provide a closed-form approximate solution, which relies on the expansion of the intractable model around an “auxiliary” one. We derive an expression for the difference between the true (but unknown) price and the auxiliary one, which we approximate in closed-form, and use to create increasingly improved refinements to the initial mispricing induced by the auxiliary model. The approach is intuitive, simple to implement, and leads to fast and extremely accurate approximations. We illustrate this method in a variety of contexts including option pricing with stochastic volatility, computation of Greeks, and the term structure of interest rates.     

American options under stochastic volatility: control variates,maturity randomization & multiscale asymptotics

American options are actively traded worldwide on exchanges, thus making their accurate and efficient pricing an important problem. As most financial markets exhibit randomly varying volatility, in this paper we introduce an approximation of an American option price under stochastic volatility models. We achieve this by using the maturity randomization method known as Canadization. The volatility process is characterized by fast and slow-scale fluctuating factors. In particular, we study the case of an American put with a single underlying asset and use perturbative expansion techniques to approximate its price as well as the optimal exercise boundary up to the first order. We then use the approximate optimal exercise boundary formula to price an American put via Monte Carlo. We also develop efficient control variates for our simulation method using martingales resulting from the approximate price formula. A numerical study is conducted to demonstrate that the proposed method performs better than the least squares regression method popular in the financial industry, in typical settings where values of the scaling parameters are small. Further, it is empirically observed that in the regimes where the scaling parameter value is equal to unity, fast and slow-scale approximations are equally accurate.     

The price of risk and ambiguity in an intertemporal general equilibrium model of asset prices

We consider a version of the intertemporal general equilibrium model of Cox et?al. (Econometrica 53:363–384, 1985) with a single production process and two correlated state variables. It is assumed that only one of them, Y ₂, has shocks correlated with those of the economy’s output rate and, simultaneously, that the representative agent is ambiguous about its stochastic process. This implies that changes in Y ₂ should be hedged and its uncertainty priced, with this price containing risk and ambiguity components. Ambiguity impacts asset pricing through two channels: the price of uncertainty associated with the ambiguous state variable, Y ₂, and the interest rate. With ambiguity, the equilibrium price of uncertainty associated with Y ₂ and the equilibrium interest rate can increase or decrease, depending on: (i) the correlations between the shocks in Y ₂ and those in the output rate and in the other state variable; (ii) the diffusion functions of the stochastic processes for Y ₂ and for the output rate; and (iii) the gradient of the value function with respect to Y ₂. As applications of our generic setting, we deduct the model of Longstaff and Schwartz (J Financ 47:1259–1282, 1992) for interest-rate-sensitive contingent claim pricing and the variance-risk price specification in the option pricing model of Heston (Rev Financ Stud 6:327–343, 1993). Additionally, it is obtained a variance-uncertainty price specification that can be used to obtain a closed-form solution for option pricing with ambiguity about stochastic variance.     

A closed-form solution for options with stochastic volatility with applications to bond and currency options

I use a new technique to derive a closed-form solution for theprice of a European call option on an asset with stochasticvolatility. The model allows arbitrary correlation between volatilityand spot asset returns. I introduce stochastic interest ratesand show how to apply the model to bond options and foreigncurrency options. Simulations show that correlation betweenvolatility and the spot asset's price is important for explainingreturn skewness and strike-price biases in the Black-Scholes(1973) model. The solution technique is based on characteristicfunctions and can be applied to other problems     

Fourier transformation and the pricing of average-rate derivatives

In this article we propose a method to compute the density of the arithmetic average of a Markov process. This approach is then applied to the pricing of average rate options (Asian options). It is demonstrated that as long as a closed form formula is available for the discount bond price when the underlying process is treated as the riskless interest rate, analytical formulas for the density function of the arithmetic average and the Asian option price can be derived. This includes the affine class of term structure models. The Cox et al. (1985) square root interest rate process is used as an example. When the underlying process follows a geometric Brownian motion, a very efficient numerical method is proposed for computing the density function of the average. Extensions of the techniques to the cases of multiple state variables are also discussed.      

Pricing American options under stochastic volatility and stochastic interest rates

We introduce a new analytical approach to price American options. Using an explicit and intuitive proxy for the exercise rule, we derive tractable pricing formulas using a short-maturity asymptotic expansion. Depending on model parameters, this method can accurately price options with time-to-maturity up to several years. The main advantage of our approach over existing methods lies in its straightforward extension to models with stochastic volatility and stochastic interest rates. We exploit this advantage by providing an analysis of the impact of volatility mean-reversion, volatility of volatility, and correlations on the American put price.     

Option pricing and hedging under a stochastic volatility Lévy process model

In this paper, we discuss a stochastic volatility model with a Lévy driving process and then apply the model to option pricing and hedging. The stochastic volatility in our model is defined by the continuous Markov chain. The risk-neutral measure is obtained by applying the Esscher transform. The option price using this model is computed by the Fourier transform method. We obtain the closed-form solution for the hedge ratio by applying locally risk-minimizing hedging.     

A market-induced mechanism for stock pinning

Abstract

We propose a model to describe stock pinning on option expiration dates. We argue that if the open interest on a particular contract is unusually large, delta-hedging in aggregate by floor market-makers can impact the stock price and drive it to the strike price of the option. We derive a stochastic differential equation for the stock price which has a singular drift that accounts for the price-impact of delta-hedging. According to this model, the stock price has a finite probability of pinning at a strike. We calculate analytically and numerically this probability in terms of the volatility of the stock, the time-to-maturity, the open interest for the option under consideration and a ‘price elasticity’ constant that models price impact.     

Option Pricing under Stochastic Interest Rates: An Empirical Investigation

Using daily data of the Nikkei 225 index, call option prices and call money rates of the Japanese financial market,a comparison is made of the pricing performance of stock option pricing modelsunder several stochastic interest rate processes proposedby the existing term structure literature.The results show that (1) one option pricing modelunder a specific stochastic interest ratedoes not significantly outperformanother option pricing model under an alternative stochasticinterest rate, and (2) incorporating stochastic interest ratesinto stock option pricing does not contribute to the performanceimprovement of the original Black–Scholes pricing formula.