New efficient numerical procedures for solving stochastic variational problems with a priori maximum pointwise error estimates

In a previous paper we gave a new formulation and derived the Euler equations and other necessary conditions to solve strong, pathwise, stochastic variational problems with trajectories driven by Brownian motion. Thus, unlike current methods which minimize the control over deterministic functionals (the expected value), we find the control which gives the critical point solution of random functionals of a Brownian path and then, if we choose, find the expected value.This increase in information is balanced by the fact that our methods are anticipative while current methods are not. However, our methods are more directly connected to the theory and meaningful examples of deterministic variational theory and provide better means of solution for free and constrained problems. In addition, examples indicate that there are methods to obtain nonanticipative solutions from our equations although the anticipative optimal cost function has smaller expected value.In this paper we give new, efficient numerical methods to find the solution of these problems in the quadratic case. Of interest is that our numerical solution has a maximal, a priori, pointwise error of O(h3/2) where h is the node size. We believe our results are unique for any theory of stochastic control and that our methods of proof involve new and sophisticated ideas for strong solutions which extend previous deterministic results by the first author where the error was O(h²).We note that, although our solutions are given in terms of stochastic differential equations, we are not using the now standard numerical methods for stochastic differential equations. Instead we find an approximation to the critical point solution of the variational problem using relations derived from setting to zero the directional derivative of the cost functional in the direction of simple test functions.Our results are even more significant than they first appear because we can reformulate stochastic control problems or constrained calculus of variations problems in the unconstrained, stochastic calculus of variations formulation of this paper. This will allow us to find efficient and accurate numerical solutions for general constrained, stochastic optimization problems. This is not yet being done, even in the deterministic case, except by the first author.     

动态随机弹性在期权定价中的应用

给出动态随机弹性的概念及运算性质,讨论了动态随机弹性在期权定价模型中的应用.主要结果有:(1)在波动率为常数时,期权价格对的弹性,得到了动态随机弹性服从运动,并给出了相应的经济解释;(2)由于波动率一般不是常数,也是随机过程,因此本文进一步研究了期权价格对波动率的弹性,就股票价格的波动情况给出了数学描述和金融意义上的解释.     

用参考圆法演示简谐振动

对受迫振动仪进行了改进，通过参考圆法，演示了作镁速圆周运动的质点在直径方向上的投影所作的运动为简谐振动。     

Nonlinear control of systems with multiple equilibria and unknown sinusoidal disturbance

In this paper, nonlinear systems having multiple equilibrium points and low order dynamics are investigated. Roll motions of ships are studied by means of modern nonlinear techniques to exemplify the behavior of such nonlinear systems in the case when they are under the influence of external sinusoidal disturbances with unknown amplitudes. The main objective is to analyze the performance of this system at different operating conditions, including those giving rise to chaos, and to design a controller with an overparameterized structure to stabilize the system at the origin. A nonlinear recursive backstepping controller is proposed and the transient performance is investigated. Lyapunov-based techniques are used to force systematic following of a reference model while introducing a nonlinear parameter estimator to guarantee adaptivity. Robustness problems as well as ways to tune the controller parameters are examined. Simulation results are submitted for the uncontrolled and controlled cases, verifying the effectiveness of the proposed controller. Finally, a discussion and conclusions are given with possible future extensions.     

Runge–Kutta Methods for Itô Stochastic Differential Equations with Scalar Noise

A general class of stochastic Runge–Kutta methods for Itô stochastic differential equation systems w.r.t. a one-dimensional Wiener process is introduced. The colored rooted tree analysis is applied to derive conditions for the coefficients of the stochastic Runge–Kutta method assuring convergence in the weak sense with a prescribed order. Some coefficients for new stochastic Runge–Kutta schemes of order two are calculated explicitly and a simulation study reveals their good performance.     

Second order Runge–Kutta methods for Stratonovich stochastic differential equations

The weak approximation of the solution of a system of Stratonovich stochastic differential equations with a m–dimensional Wiener process is studied. Therefore, a new class of stochastic Runge–Kutta methods is introduced. As the main novelty, the number of stages does not depend on the dimension m of the driving Wiener process which reduces the computational effort significantly. The colored rooted tree analysis due to the author is applied to determine order conditions for the new stochastic Runge–Kutta methods assuring convergence with order two in the weak sense. Further, some coefficients for second order stochastic Runge–Kutta schemes are calculated explicitly. AMS subject classification (2000)  65C30, 65L06, 60H35, 60H10     

A Distributional Approach to Multiple Stochastic Integrals and Transformations of the Poisson Measure

We study a class of ‘nonpoissonian’ transformations of the configuration space and the corresponding transformations of the Poisson measure. For some class of Poisson measures we find conditions which are sufficient for the transformed measure (which in general is nonpoissonian) to be absolutely continuous with respect to the initial Poisson measure and get the expression for the corresponding Radon–Nikodym derivative. To solve this problem we use a distributional approach to Poisson multiple stochastic integrals.     

阻尼落体运动的分析力学研究

应用分析力学理论和方法,研究了两种情况下的阻尼落体运动:1)阻力大小与速度成正比;2)阻力大小与速度平方成正比.对两种运动分别给出了等效的Lagrange函数和Hamilton函数,并应用第一积分法、点变换法、正则变换法和Ham-ilton-Jacobi方程法等不同的求解方法进行了求解.