A drift‐free simulation method for pricing commodity derivatives

Having in view the pricing of commodity derivatives in Libor Market Model (LMM) setting, we first analyze the set of basic rates we need to formulate the model by using the spanning tree concept taken from graph theory. Next, we present an efficient procedure for Monte Carlo simulation of the dynamics of the rates associated to LMM, avoiding the presence of the rates dependent drifts (drift‐free simulation) and the presence of negative deflated bond prices and negative forward rates. The method is based upon a new parameterization of the martingales introduced by Glasserman and Zhao and it is extended to a Cross‐Market Model for commodities. Finally, a particular example of commodity derivative (spread option) pricing problem is considered. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling the evolution of credit spreads using the Cox process within the HJM framework: A CDS option pricing model

In this paper a simulation approach for defaultable yield curves is developed within the Heath et al. (1992) framework. The default event is modelled using the Cox process where the stochastic intensity represents the credit spread. The forward credit spread volatility function is affected by the entire credit spread term structure. The paper provides the defaultable bond and credit default swap option price in a probability setting equipped with a subfiltration structure. The Euler–Maruyama stochastic integral approximation and the Monte Carlo method are applied to develop a numerical scheme for pricing. Finally, the antithetic variable technique is used to reduce the variance of credit default swap option prices.     

基于贝叶斯MCMC算法的美式期权定价

鉴于美式期权的定价具有后向迭代搜索特征,本文结合Longstaff和Schwartz提出的美式期权定价的最小二乘模拟方法,研究基于马尔科夫链蒙特卡洛算法对回归方程系数的估计,实现对美式期权的双重模拟定价.通过对无红利美式看跌股票期权定价进行大量实证模拟,从期权价值定价误差等方面同著名的最小二乘蒙特卡洛模拟方法进行对比分析,结果表明基于MCMC回归算法给出的美式期权定价具有更高的精确度.模拟实证结果表明本文提出的对美式期权定价方法具有较好的可行性、有效性与广泛的适用性.该方法的不足之处就是类似于一般的蒙特卡洛方法,会使得求解的计算量有所加大.     

资产跳跃下CM策略多期收入保证价格模拟

恒定混合策略(CM策略)多期收入保证价格是保本基金发行方采取设置止损的CM\linebreak策略作为投资策略时收取保 本费的理论依据, 其中标的资产由复合泊松过程和维纳过程共同驱动, 这一定价问题内嵌奇异期权, 蒙特卡罗模拟方法擅长处理这种高维数量金融问题. 基于风险中性测度推导出多期收入保证价格的现值表达式, 用条件蒙特卡罗推导出这一现值表达式的模拟公式. 在给定参数下分别用普通蒙特卡罗和条件蒙特卡罗计算CM策略多期收入保证价格的数值解, 结果显示两种蒙特卡罗方法均能有效计算其数值解, 之后通过给定显著性水平下的置信区间长度评价两种方法的精确度, 结果显示条件蒙特卡罗比普通蒙特卡罗有很大改进. 接着运用条件蒙特卡罗模拟研究多期收入保证价格对不同参数范围的变化情况.     

Pricing Variance Swaps in a Hybrid Model of Stochastic Volatility and Interest Rate with Regime-Switching

In this paper, we consider the problem of pricing discretely-sampled variance swaps based on a hybrid model of stochastic volatility and stochastic interest rate with regime-switching. Our modeling framework extends the Heston stochastic volatility model by including the Cox-Ingersoll-Ross (CIR) stochastic interest rate model. In addition, certain model parameters in our model switch according to a continuous-time observable Markov chain process. This enables our model to capture several macroeconomic issues such as alternating business cycles. A semi-closed form pricing formula for variance swaps is derived. The pricing formula is assessed through numerical implementation, where we validate our pricing formula against the Monte Carlo simulation. The impact of incorporating regime-switching for pricing variance swaps is also discussed, where variance swaps prices with and without regime-switching effects are examined in our model. We also explore the economic consequence for the prices of variance swaps by allowing the Heston-CIR model to switch across three different regimes.     

A linear algebraic method for pricing temporary life annuities and insurance policies

We recast the valuation of annuities and life insurance contracts under mortality and interest rates, both of which are stochastic, as a problem of solving a system of linear equations with random perturbations. A sequence of uniform approximations is developed which allows for fast and accurate computation of expected values. Our reformulation of the valuation problem provides a general framework which can be employed to find insurance premiums and annuity values covering a wide class of stochastic models for mortality and interest rate processes. The proposed approach provides a computationally efficient alternative to Monte Carlo based valuation in pricing mortality-linked contingent claims.     

MC方差减小技术在算术平均亚式外汇期权定价中的应用

建立了利率和汇率波动率均为随机情形下算术平均亚式外汇期权的定价模型.由于其定价问题求解十分困难,运用蒙特卡罗(Monte Carlo)方法并结合控制变量方差减小技术进行模拟,有效地减小了模拟方差,得到了期权定价问题的数值结果.     

上证50ETF期权定价有效性的研究:基于B-S-M模型和蒙特卡罗模拟

上证50ETF期权作为中国资本市场上股票期权的第一个试点产品,其定价问题尤为重要。本文分别运用B-S-M期权定价模型和蒙特卡罗模拟方法对其定价进行实证研究,分析结果表明:1)IGARCH模型比传统的GARCH模型更能较好地拟合上证50ETF的波动率;2)当模拟次数为1000时,蒙特卡罗方法的效率一致地高于B-S-M模型,并且除了对偶变量技术的拟蒙特卡罗其他模型的精确度也都高于B-S-M模型;3)B-S-M模型和蒙特卡罗模拟方法都可以较为准确地、有效地模拟出上证50ETF期权价格。这些研究将为今后期权定价模型的发展和完善提供必要的参考和指引。     

基于LSMC模型的煤炭资源价值定价研究

煤炭资源价值定价可以抽象为一种美式期权定价问题.最小二乘蒙特卡洛模拟(LSMC)方法是解决美式期权定价问题的一个有效途径.详尽地分析了Cortazar等人的基于资源价格、利率和便利收益随机变动的三因素定价模型,利用向量Ito定理提出了三因素模型中价格、利率和便利收益变量的递推公式.对LSMC方法原理进行了细致的阐述,总结出实现LSMC方法的完整过程,并在Matlab环境下编制了LSMC算法实现程序,进行算例计算.算例结果表明,LSMC方法用于资源定价是有效可靠的.研究为煤炭资源价值定价提供了一个完整具有可操作性的工具.     

Numerical solution of a PDE model for a ratchet-cap pricing with BGM interest rate dynamics

In this paper we present a new numerical method to price an interest rate derivative. The financial product consists of a particular ratchet cap contract which contains a set of ratchet caplets. For this purpose, we first pose the PDE pricing model for each ratchet caplet by means of Feynman-Kac theorem. The underlying interest rates are the forward LIBOR rates, the dynamics of which are assumed to follow the recently introduced BGM (LMM) market model. For the set of PDEs associated to the ratchet caplets pricing problems, we propose a second order Crank-Nicolson characteristics time discretization scheme combined with a finite element discretization in the interest rate variables. In order to illustrate the performance of the numerical methods, we present an academic test and a real example of a particular ratchet cap pricing. In the second case, a comparison between the results obtained by Monte Carlo simulation and the proposed method is presented.     

CEV模型和B-P混合驱动模型下亚式期权Monte Carlo模拟

对亚式期权在CEV模型和B-P混合驱动模型限制下进行Monte Carlo模拟定价,建立风险中性测度,模拟出不同弹性因子值下资产价格路径.为了得出优于标准的Monte Carlo模拟,应用方差缩减技术来提高期权定价的精度.最后对亚式期权定价模型进行数值案例分析,得出弹性因子取值、时间步长、模拟次数与期权价值变化的关系.     

An Efficient Monte Carlo Method for Optimal Control Problems with Uncertainty

A general framework is proposed for what we call the sensitivity derivative Monte Carlo (SDMC) solution of optimal control problems with a stochastic parameter. This method employs the residual in the first-order Taylor series expansion of the cost functional in terms of the stochastic parameter rather than the cost functional itself. A rigorous estimate is derived for the variance of the residual, and it is verified by numerical experiments involving the generalized steady-state Burgers equation with a stochastic coefficient of viscosity. Specifically, the numerical results show that for a given number of samples, the present method yields an order of magnitude higher accuracy than a conventional Monte Carlo method. In other words, the proposed variance reduction method based on sensitivity derivatives is shown to accelerate convergence of the Monte Carlo method. As the sensitivity derivatives are computed only at the mean values of the relevant parameters, the related extra cost of the proposed method is a fraction of the total time of the Monte Carlo method.     

Recursive Marginal Quantization of the Euler Scheme of a Diffusion Process

Abstract

We propose a new approach to quantize the marginals of the discrete Euler diffusion process. The method is built recursively and involves the conditional distribution of the marginals of the discrete Euler process. Analytically, the method raises several questions like the analysis of the induced quadratic quantization error between the marginals of the Euler process and the proposed quantizations. We show in particular that at every discretization step t_k of the Euler scheme, this error is bounded by the cumulative quantization errors induced by the Euler operator, from times t₀ = 0 to time t_k. For numerics, we restrict our analysis to the one-dimensional setting and show how to compute the optimal grids using a Newton–Raphson algorithm. We then propose a closed formula for the companion weights and the transition probabilities associated to the proposed quantizations. This allows us to quantize in particular diffusion processes in local volatility models by reducing dramatically the computational complexity of the search of optimal quantizers while increasing their computational precision with respect to the algorithms commonly proposed in this framework. Numerical tests are carried out for the Brownian motion and for the pricing of European options in a local volatility model. A comparison with the Monte Carlo simulations shows that the proposed method may sometimes be more efficient (w.r.t. both computational precision and time complexity) than the Monte Carlo method.     

Arbitrage valuation and bounds for sinking-fund bonds with multiple sinking-fund dates

The paper tackles the problem of pricing, under interest-rate risk, a default-free sinking-fund bond which allows its issuer to recurrently retire part of the issue by (a) a lottery call at par, or (b) an open market repurchase. By directly modelling zero-coupon bonds as diffusions driven by a single-dimensional Brownian motion, a pricing formula is supplied for the sinking-fund bond based on a backward induction procedure which exploits, at each step, the martingale approach to the valuation of contingent-claims. With more than one sinking-fund date, however, the pricing formula is not in closed form, not even for simple parametrizations of the process for zerocoupon bonds, so that a numerical approach is needed. Since the computational complexity increases exponentially with the number of sinking-fund dates, arbitrage-based lower and upper bounds are provided for the sinking-fund bond price. The computation of these bounds is almost effortless when zero-coupon bonds are as described by Cox, Ingersoll and Ross. Numerical comparisons between the price of the sinking-fund bond obtained via Monte Carlo simulation and these lower and upper bounds are illustrated for different choices of parameters.     

Pricing of Ratchet equity-indexed annuities under stochastic interest rates

We consider the valuation of simple and compound Ratchet equity-indexed annuities (EIAs) in the presence of stochastic interest rates. We assume that the equity index follows a geometric Brownian motion and the short rate follows the extended Vasicek model. Under a given forward measure, we obtain an explicit multivariate normal characterization for multiple log-returns on the equity index. Using such a characterization, closed-form price formulas are derived for both simple and compound Ratchet EIAs. An efficient Monte Carlo simulation scheme is also established to overcome the computational difficulties resulting from the evaluation of high-dimensional multivariate normal cumulative distribution functions (CDFs) embedded in the price formulas as well as the consideration of additional complex contract features. Finally, numerical results are provided to illustrate the computational efficiency of our simulation scheme and the effects of various model and contract parameters on pricing.     

光学照明系统设计中若干问题的研究

讨论了光学照明系统设计的各种计算方案 .认为基于蒙特卡洛方法的算法简单 ,适用范围广 ,运行效率高 ,对大范围光场分布可获得一致收敛结果 .数值模拟结果显示 ,对于计算小尺度光源系统在远距离照明的问题具有优势 .围绕 2 0 0 2年数模竞赛 A题 ,给出了若干结果 ,并根据问题的特点进行了分析 .同时讨论了在数值积分中对扩展光源 (线光源 )如何正确处理剖分 :向所有方向漫射的扩展光源符合光学中的朗伯定律 ,剖分时应予考虑 .     

金融中的Levy模型及其仿真

近20年来,金融中Levy模型与蒙特卡洛仿真技术日益受到重视. 在连续时间过程的金融建模中带跳跃的Levy模型相比于连续轨道的布朗运动模型能很好地刻画市场的跳跃,更好地拟合金融数据的统计特征,更准确地对衍生品定价. 但是,相较于经典的Black-Scholes模型,用Levy模型对衍生品定价以及求解对冲策略的计算复杂度大大增加. 蒙特卡洛仿真成为Levy模型计算中最重要的方法之一. 首先详细地介绍了Levy模型引入的背景,并引出仿真方法在其中重要的应用价值. 最后,简要地给出了Levy过程仿真及其梯度估计的基本方法.     

A GARCH option pricing model with α-stable innovations

We develop an option pricing model which is based on a GARCH asset return process with α-stable innovations with truncated tails. The approach utilizes a canonic martingale measure as pricing measure which provides the possibility of a model calibration to market prices. The GARCH-stable option pricing model allows the explanation of some well-known anomalies in empirical data as volatility clustering and heavy tailedness of the return distribution. Finally, the results of Monte Carlo simulations concerning the option price and the implied volatility with respect to different strike and maturity levels are presented.     

Pricing and hedging GLWB in the Heston and in the Black–Scholes with stochastic interest rate models

Valuing Guaranteed Lifelong Withdrawal Benefit (GLWB) has attracted significant attention from both the academic field and real world financial markets. As remarked by Forsyth and Vetzal (2014) the Black and Scholes framework seems to be inappropriate for such a long maturity products. They propose to use a regime switching model. Alternatively, we propose here to use a stochastic volatility model (Heston model) and a Black–Scholes model with stochastic interest rate (Hull–White model). For this purpose we present four numerical methods for pricing GLWB variables annuities: a hybrid tree-finite difference method and a Hybrid Monte Carlo method, an ADI finite difference scheme, and a Standard Monte Carlo method. These methods are used to determine the no-arbitrage fee for the most popular versions of the GLWB contract, and to calculate the Greeks used in hedging. Both constant withdrawal and optimal withdrawal (including lapsation) strategies are considered. Numerical results are presented which demonstrate the sensitivity of the no-arbitrage fee to economic, contractual and longevity assumptions.     

利用随机模拟方法求解第二类积分方程

给出一种求解第二类Fredholm和Volterra积分方程的数值算法,算法在数值积分技术的基础上使用Monte Carlo随机模拟方法求积分方程的近似解.通过数值例子证明了该算法是有效的.