Option Pricing with Stochastic Volatility, Equity Premium, and Interest Rates

Nicole Hao,Echo Li,Diep Luong-Le
2024-08-28
Abstract:This paper presents a new model for options pricing. The Black-Scholes-Merton (BSM) model plays an important role in financial options pricing. However, the BSM model assumes that the risk-free interest rate, volatility, and equity premium are constant, which is unrealistic in the real market. To address this, our paper considers the time-varying characteristics of those parameters. Our model integrates elements of the BSM model, the Heston (1993) model for stochastic variance, the Vasicek model (1977) for stochastic interest rates, and the Campbell and Viceira model (1999, 2001) for stochastic equity premium. We derive a linear second-order parabolic PDE and extend our model to encompass fixed-strike Asian options, yielding a new PDE. In the absence of closed-form solutions for any options from our new model, we utilize finite difference methods to approximate prices for European call and up-and-out barrier options, and outline the numerical implementation for fixed-strike Asian call options.
Mathematical Finance,Analysis of PDEs,Computational Finance
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is the limitations of existing option - pricing models (such as the Black - Scholes - Merton model) in the actual market. Specifically, the Black - Scholes - Merton model assumes that the risk - free interest rate, volatility, and equity premium are constants, which is not true in reality. In order to more accurately reflect market conditions, the author introduced a new option - pricing model that takes into account the time - varying characteristics of these parameters. ### Detailed description of the main problem 1. **Limitations of existing models**: - **Black - Scholes - Merton (BSM) model**: It is unrealistic in the actual market to assume that the risk - free interest rate, volatility, and equity premium are constants. 2. **Introduction of time - varying parameters**: - To improve the BSM model, the author integrated the following models: - The **Heston (1993)** model for stochastic volatility. - The **Vasicek (1977)** model for stochastic risk - free interest rate. - The **Campbell and Viceira (1999, 2001)** model for stochastic equity premium. 3. **Objectives of the new model**: - **More accurate pricing**: By introducing time - varying risk - free interest rate, volatility, and equity premium, the model is made closer to the actual market situation. - **Expanded application scope**: It is not only applicable to traditional European options but also extended to complex financial instruments such as fixed - strike Asian options. 4. **Mathematical model**: - The author derived a linear second - order parabolic partial differential equation (PDE) and further extended it to cover fixed - strike Asian options, obtaining a new PDE. - In the absence of a closed - form solution, the finite - difference method is used to approximately solve the prices of European call options and knock - out barrier options, and the numerical implementation of fixed - strike Asian call options is outlined. ### Key features of the new model - **Stochastic differential equations (SDEs)**: \[ \begin{cases} dS(t) = (\mu + X(t)+R(t))S(t)dt+\sqrt{\sigma_s(t)}S(t)dW_1(t) \\ dX(t)=-\kappa_x X(t)dt+\sigma_x(\rho_x dW_1(t)+\sqrt{1 - \rho_x^2}dW_2(t)) \\ d\sigma_s(t)=\kappa_s(\sigma-\sigma_s(t))dt+\eta\sqrt{\sigma_s(t)}(\rho_s dW_1(t)+\sqrt{1-\rho_s^2}dW_3(t)) \\ dR(t)=\kappa_r(r - R(t))dt+\sigma_r(\rho_r dW_1(t)+\sqrt{1-\rho_r^2}dW_4(t)) \end{cases} \] where \(S(t)\) is the price of the underlying asset, \(X(t)\) is the deviation of the equity premium from its mean, \(\sigma_s(t)\) is the volatility, \(R(t)\) is the risk - free interest rate, and \(W_1, W_2, W_3, W_4\) are independent Brownian motions. - **Partial differential equation (PDE)**: For European derivatives, the derived PDE is: \[ V_t = R(V - S V_S - X V_X-\sigma_s V_{\sigma_s}-R V_R)-\frac{1}{2}\sigma_s S^2 V_{SS}-\frac{1}{2}\sigma_x^2 V_{XX}-\frac{1}{2}\eta^2\sigma_s V_{\sigma_s\sig}