一类随机线性互补问题的求法

考虑一类随机互线性补问题的求解方法,目的是通过定义NCP函数来使正则化期望残差最小化.通过拟蒙洛包洛方法产生一系列观察值并且证得离散近似问题最小值解的聚点就是相应随机线性互补问题的期望残差最小值ERM,同时得到利用ERM到解为有界的充分条件.进一步证明ERM法能够得到具有稳定性和最小灵敏度的稳健解.     

The pricing of Asian options under stochastic interest rates

The purpose of this paper is to analyse the effect of stochastic interest rates on the pricing of Asian options. It is shown that a stochastic, in contrast to a deterministic, development of the term structure of interest rates has a significant influence. The price of the underlying asset, e.g. a stock or oil, and the prices of bonds are assumed to follow correlated two-dimensional Itô processes. The averages considered in the Asian options are calculated on a discrete time grid, e.g. all closing prices on Wednesdays during the lifetime of the contract. The value of an Asian option will be obtained through the application of Monte Carlo simulation, and for this purpose the stochastic processes for the basic assets need not be severely restricted. However, to make comparison with published results originating from models with deterministic interest rates, we will stay within the setting of a Gaussian framework.     

Numerical comparison of three stochastic methods for nonlinear PN junction problems

We apply the Monte Carlo, stochastic Galerkin, and stochastic collocation methods to solving the drift-diffusion equations coupled with the Poisson equation arising in semiconductor devices with random rough surfaces. Instead of dividing the rough surface into slices, we use stochastic mapping to transform the original deterministic equations in a random domain into stochastic equations in the corresponding deterministic domain. A finite element discretization with the help of AFEPack is applied to the physical space, and the equations obtained are solved by the approximate Newton iterative method. Comparison of the three stochastic methods through numerical experiment on different PN junctions are given. The numerical results show that, for such a complicated nonlinear problem, the stochastic Galerkin method has no obvious advantages on efficiency except accuracy over the other two methods, and the stochastic collocation method combines the accuracy of the stochastic Galerkin method and the easy implementation of the Monte Carlo method.     

Unbiased multi-index Monte Carlo

We introduce a new class of Monte Carlo-based approximations of expectations of random variables such that their laws are only available via certain discretizations. Sampling from the discretized versions of these laws can typically introduce a bias. In this paper, we show how to remove that bias, by introducing a new version of multi-index Monte Carlo (MIMC) that has the added advantage of reducing the computational effort, relative to i.i.d. sampling from the most precise discretization, for a given level of error. We cover extensions of results regarding variance and optimality criteria for the new approach. We apply the methodology to the problem of computing an unbiased mollified version of the solution of a partial differential equation with random coefficients. A second application concerns the Bayesian inference (the smoothing problem) of an infinite-dimensional signal modeled by the solution of a stochastic partial differential equation that is observed on a discrete space grid and at discrete times. Both applications are complemented by numerical simulations.     

Stochastic Volatility: Option Pricing using a Multinomial Recombining Tree

In this paper, we present a natural mathematical framework to model trader behavior as a continuous time discrete event process, and derive stochastic differential equations for aggregate behavior and price dynamics by passing to diffusion limits. In particular, we model extraneous, value, momentum and hedge traders. Through analysis and numerical simulation we explore some of the effects these trading strategies have on price dynamics.     

带补偿的随机二次规划的近似Lagrange-Newton方法

In this paper, two-stage stochastic quadratic programming problems with equality constraints are considered. By Monte Carlo simulation-based approximations of the objective function and its first (second)derivative,an inexact Lagrange-Newton type method is proposed.It is showed that this method is globally convergent with probability one. In particular, the convergence is local superlinear under an integral approximation error bound condition.Moreover, this method can be easily extended to solve stochastic quadratic programming problems with inequality constraints.     

Revisit of stochastic mesh method for pricing American options

From an importance sampling viewpoint, Broadie and Glasserman [M. Broadie, P. Glasserman, A stochastic mesh method for pricing high-dimensional American options, Journal of Computational Finance 7 (4) (2004) 35–72] proposed a stochastic mesh method to price American options. In this paper, we revisit the method from a conditioning viewpoint, and derive some new weights.     

A methodology for stochastic inventory models based on a zero‐adjusted Birnbaum‐Saunders distribution

The Birnbaum–Saunders (BS) distribution is receiving considerable attention. We propose a methodology for inventory logistics that allows demand data with zeros to be modeled by means of a new discrete–continuous mixture distribution, which is constructed by using a probability mass at zero and a continuous component related to the BS distribution. We obtain some properties of the new mixture distribution and conduct a simulation study to evaluate the performance of the estimators of its parameters. The methodology for stochastic inventory models considers also financial indicators. We illustrate the proposed methodology with two real‐world demand data sets. It shows its potential, highlighting the convenience of using it by improving the contribution margins of a Chilean food industry. Copyright © 2015 John Wiley & Sons, Ltd.     

随机线性互补问题的无约束优化再定式

针对随机线性互补问题,提出等价的无约束优化再定式模型,即由D-间隙函数定义的确定性的无约束期望残差极小化问题.通过拟Monte Carlo方法,将样本进行了推广,得到了相关的离散近似问题.在适当的条件下,提出了最优解存在的充分条件,以及探究了离散近似问题的最优解及稳定点的收敛性.另外,在针对一类带有常系数矩阵的随机互补线性问题,研究了解存在的充要条件.     

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.     

非仿射随机波动率的欧式障碍期权定价

研究非仿射随机波动率模型的欧式障碍期权定价问题时,首先介绍了非仿射随机波动率模型,其次利用投资组合和It^o引理,得到了该模型下扩展的Black-Schole偏微分方程.由于这个方程没有显示解,因此采用对偶蒙特卡罗模拟法计算欧式障碍期权的价格.最后,通过数值实例验证了算法的可行性和准确性.     

Monte Carlo methods in applied mathematics and computational aerodynamics

A survey of the Monte Carlo methods developed in the computational aerodynamics of rarefied gases is given, and application of these methods in unconventional fields is described. A short history of these methods is presented, and their advantages and drawbacks are discussed. A relationship of the direct statistical simulation of aerodynamical processes with the solution of kinetic equations is established; it is shown that the modern stage of the development of computational methods is impossible without the use of the complex approach to the development of algorithms with regard for all the specific features of the problem to be solved (its physical nature, mathematical model, the theory of computational mathematics, and stochastic processes). Possible directions of the development of the statistical simulation methods are discussed.     

Monte Carlo and quasi-Monte Carlo sampling methods for a class of stochastic mathematical programs with equilibrium constraints

In this paper, we consider a class of stochastic mathematical programs with equilibrium constraints introduced by Birbil et al. (Math Oper Res 31:739–760, 2006). Firstly, by means of a Monte Carlo method, we obtain a nonsmooth discrete approximation of the original problem. Then, we propose a smoothing method together with a penalty technique to get a standard nonlinear programming problem. Some convergence results are established. Moreover, since quasi-Monte Carlo methods are generally faster than Monte Carlo methods, we discuss a quasi-Monte Carlo sampling approach as well. Furthermore, we give an example in economics to illustrate the model and show some numerical results with this example. The first author’s work was supported in part by the Scientific Research Grant-in-Aid from Japan Society for the Promotion of Science and SRF for ROCS, SEM. The second author’s work was supported in part by the United Kingdom Engineering and Physical Sciences Research Council grant. The third author’s work was supported in part by the Scientific Research Grant-in-Aid from Japan Society for the Promotion of Science.     

Weak Approximation of Stochastic Differential Equations and Application to Derivative Pricing

A new, simple algorithm of order 2 is presented to approximate weakly stochastic differential equations. It is then applied to the problem of pricing Asian options under the Heston stochastic volatility model.

2000 Mathematics Subject Classification, 65C30, 65C05.     

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.     

A Monte Carlo study of old and new frontier methods for efficiency measurement

This study presents the results of an extensive Monte Carlo experiment to compare different methods of efficiency analysis. In addition to traditional parametric–stochastic and nonparametric–deterministic methods recently developed robust nonparametric–stochastic methods are considered. The experimental design comprises a wide variety of situations with different returns-to-scale regimes, substitution elasticities and outlying observations. As the results show, the new robust nonparametric–stochastic methods should not be used without cross-checking by other methods like stochastic frontier analysis or data envelopment analysis. These latter methods appear quite robust in the experiments.     

Approximation for non-smooth functionals of stochastic differential equations with irregular drift

This paper aims at developing a systematic study for the weak rate of convergence of the Euler–Maruyama scheme for stochastic differential equations with very irregular drift and constant diffusion coefficients. We apply our method to obtain the rates of approximation for the expectation of various non-smooth functionals of both stochastic differential equations and killed diffusion. We also apply our method to the study of the weak approximation of reflected stochastic differential equations whose drift is Hölder continuous.     

Solving multistage asset investment problems by the sample average approximation method

The vast size of real world stochastic programming instances requires sampling to make them practically solvable. In this paper we extend the understanding of how sampling affects the solution quality of multistage stochastic programming problems. We present a new heuristic for determining good feasible solutions for a multistage decision problem. For power and log-utility functions we address the question of how tree structures, number of stages, number of outcomes and number of assets affect the solution quality. We also present a new method for evaluating the quality of first stage decisions.     

Comparison of techniques for assessing effects of loading uncertainty upon a long term phosphorous model

This paper compares the feasibility of applying three stochastic techniques to a linear water quality model. The Monte Carlo, first order analysis, and generation of moment equation techniques are applied to a long term phosphorus model of Lake Washington. The effect of uncertainty of the phosphorus loading term on simulated phosphorus levels is analysed. The simulated concentrations of phosphorus in the water column are very responsive to uncertainty in annual phosphorus loading, but the sediment concentrations are relatively insensitive. All three stochastic techniques produced identical results, but the level of preparatory and computational effort required varies considerably. The Monte Carlo technique requires the most computation time of the three stochastic techniques examined. The first order analysis and generation of moment equation techniques are shown to be precise and efficient methods of stochastic analysis. In this application they required less than one-thousandth the computer time of the Monte Carlo technique     

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

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