Vasicek随机利率模型下欧式期权 定价的Mellin变换法

运用Feynman-Kac公式和偏微分方程法得到Vasicek随机利率模型下的零息债券价格公式.利用△-对冲方法建立该模型下欧式期权价值满足的偏微分方程模型,并用Mellin变换法求解该偏微分方程,最终得到欧式期权定价公式.从数值算例的结果可以看出Mellin变换法的有效性以及不同参数对期权价值的影响.     

Stochastic Stability of Some Mechanical Systems

We discuss the behavior, for large values of time, of two linear stochastic mechanical systems. The systems are similar mathematically in that they contain a white noise in their parameters. The initial data may be random as well but are independent of white noise. The expected energy is calculated in both cases. It is well known that for free nonstochastic mechanical systems with viscous damping, the energy approaches zero as time increases. We check that this behavior takes place for the stochastic systems under consideration in the case when the initial data are random but the parameters are not. When the parameters contain a random noise the expected energy may be infinite, approach zero, remain bounded, or increase with no bound. This regime is similar to but more interesting than the known regime for the solutions of differential equations with time dependent periodic coefficients that describes the behavior of a mechanical system with characteristics that are periodic functions of time. We give necessary and sufficient conditions for stability of both systems in terms of the structure of the set of roots of an auxiliary equation.     

两参数Ito型随机微分方程解的收敛定理

设ωｚ是Ｒ＾２＋上的布朗单,考虑两参数Ｉｔｏ型随机微分方程：ｄｘｚ＝ａ（ｚ,ｘｚ）ｄωｚ＋ｂ（ｚ,ｘｚ）ｄｚ（１）ｄｘ＾＊ｚ＝ａｚ（ｚ,ｘ＾＊ｚ）ｄωｚ＋ｂｚ（ｚ,ｘ＾＊ｚ）ｄｚ（２）则在方程系数满足一定条件下,本证明了方程（２）的解向方程（１）的解收敛。     

Parabolic regularzation of a first order stochastic partial differential equation

In this paper we use regularization methods for proving the existence and uniqueness of smooth solutions of a first order semilinear stochastic partial differential equation. The regularizations are chosen in such a way so that the known theory of stochastich parabolic Ito equations can be applied. The existence of the generalized solutions and, if the the time parameter is the whole real axis, the existence of mean square bounded generalized solutions, is also considered     

模糊随机微分方程解的存在唯一性

将实数空间上的随机微分方程推广到模糊数空间,即为模糊随机微分方程.本文用Picard迭代的方法证明了其解的存在唯一性定理,推广了现有文献的结果,并且给出Picard迭代近似解误差的估计式.     

Limit theorems for BSDE with local time applications to non-linear PDE

Given a d -dimensional Wiener process W , with its natural filtration F t , a F T -measurable random variable ξ in R , a bounded measure x on R , and an adapted process ( s , y , z ) M h ( s , y , z ), we consider the following BSDE: Y t = ξ + Z t T h ( s , Y s , Z s ) d s + Z R ( L T a ( Y ) m L t a ( Y )) x (d a ) m Z t T Z s d W s for 0 h t h T . Here L t a ( Y ) stands for the local time of Y at level a . For h =0, we establish the existence and the uniqueness of the processes ( Y , Z ), and if h is continuous with linear growth we establish the existence of a solution. We prove limit theorems for solutions of backward stochastic differential equations of the above form. Those limit theorems permit us to deduce that any solution of that equation is the limit, in a strong sense, of a sequence of semi-martingales, which are solutions of ordinary BSDEs of the form Y t = ξ + Z t T f ( Y s ) Z s 2 d s m Z t T Z s d W s . A comparison theorem for BSDEs involving measures is discussed. As an application we obtain, with the help of the connection between BSDE and PDE, some corresponding limit theorems for a class of singular non-linear PDEs and a new probabilistic proof of the comparison theorem for PDEs.     

Stochastic partial differential equations in bounded domains with Dirichlet boundary conditions

In this paper we prove the existence of a unique solution for a class of stochastic parabolic partial differential equations in bounded domains, with Dirichlet boundary conditions. The main tool is an equivalence result, provided by the stochastic characteristics method, between the stochastic equations under investigation and a class of deterministic parabolic equations with moving boundaries, depending on random coefficients. We show the existence of the solution to this last problem, thus providing a solution to the former.     

The Bismut-Elworthy formula for backward SDE's and applications to nonlinear Kolmogorov equations and control in infinite dimensional spaces

We consider a forward-backward system of stochastic evolution equations in a Hilbert space. Under nondegeneracy assumptions on the diffusion coefficient (that may be nonconstant) we prove an analogue of the well-known Bismut-Elworthy formula. Next, we consider a nonlinear version of the Kolmogorov equation, i.e. a deterministic quasilinear equation associated to the system according to Pardoux, E and Peng, S. (1992). "Backward stochastic differential equations and quasilinear parabolic partial differential equations". In: Rozowskii, B.L., Sowers, R.B. (Eds.), Stochastic Partial Differential Equations and Their Applications , Lecture Notes in Control Inf. Sci., Vol. 176, pp. 200-217. Springer: Berlin. The Bismut-Elworthy formula is applied to prove smoothing effect, i.e. to prove existence and uniqueness of a solution which is differentiable with respect to the space variable, even if the initial datum and (some) coefficients of the equation are not. The results are then applied to the Hamilton-Jacobi-Bellman equation of stochastic optimal control. This way we are able to characterize optimal controls by feedback laws for a class of infinite-dimensional control systems, including in particular the stochastic heat equation with state-dependent diffusion coefficient.     

Solving parabolic Wick-stochastic boundary value problems using a finite element method

A class of linear parabolic stochastic boundary value problems of Wick-type is studied. The equations are understood in a weak sense on a suitable stochastic distribution space, and existence and uniqueness results are provided. The paper continues to discuss a numerical method for this type of problem, based on a Galerkin type of approximation. Estimates showing linear convergence in time and space are derived, and rate of convergence results for the stochastic dimension are reported.     

A stochastic lagrangean picture for the three dimensional navier-stokes equation +

We derive a representation formula for the solutions of the navier-Stokes flow in three dimensions in absence of boundary in terms of stochastic current lines. This pictures generalizes an analogous one given in two dimensions and on the other hand the classical Lafrangian picture for the Euler flow. We prove a Vanishing viscosity limit for the whole structure     

A hybrid method for pricing European options based on multiple assets with transaction costs

The problem of pricing European options based on multiple assets with transaction costs is considered. These options include, for example, quality options and options on the minimum of two or more risky assets. The value of these options is the solution of a nonlinear parabolic partial differential equation subject to a final condition given by the payoff function associated with the option. A computationally efficient method to solve this final-value problem is proposed. This method is based on an asymptotic expansion of the required solution with respect to the parameters related to the transaction costs followed by the numerical solution of the linear partial differential equations obtained at each order in perturbation theory. The numerical solution of these linear problems involves an implicit finite-difference scheme for the parabolic equation and the use of the fast Fourier sine transform to solve the resulting elliptic problems. Numerical results obtained on test problems with the method proposed here are shown and discussed.     

An efficient iterative method for solving the matrix equation AXB + CYD = E

This paper presents an iterative method for solving the matrix equation AXB + CYD = E with real matrices X and Y. By this iterative method, the solvability of the matrix equation can be determined automatically. And when the matrix equation is consistent, then, for any initial matrix pair [X₀, Y₀], a solution pair can be obtained within finite iteration steps in the absence of round‐off errors, and the least norm solution pair can be obtained by choosing a special kind of initial matrix pair. Furthermore, the optimal approximation solution pair to a given matrix pair [X?, ?] in a Frobenius norm can be obtained by finding the least norm solution pair of a new matrix equation AX?B + C?D = ?, where ? = E ? AX?B ? C?D. The given numerical examples show that the iterative method is efficient. Copyright © 2005 John Wiley & Sons, Ltd.     

中立型方程x(t)+cx(t-ㄒ)+ax(t)+bx(t-ㄒ)=0零解渐近稳定的代数判据

本文建立了方程ｘ（ｔ）＋ｃｔ（ｔ－ｔ）＋ａｘ（ｔ）＋ｂｘ（ｔ－ｔ）＝０零解渐近稳定的充要条件,给出了其零解渐近稳定的代数判据,同时纠正了文［２］中出现的错误．     

Comments on ‘An efficient iterative method for solving the matrix equation AXB+CYD=E’

In this note, a technical error is pointed out in the proof of a lemma in the above paper. A correct proof of this lemma is given. In addition, a further result on the algorithm in the above paper is also given. Copyright © 2009 John Wiley & Sons, Ltd.     

A method for the resolution of the Jacobi equation Y ″ + RY = 0 on the manifold Sp (2)/ SU (2)

In this paper a method for the resolution of the differential equation of the Jacobi vector fields in the manifold V ₁ = Sp(2)/SU(2) is exposed. These results are applied to determine areas and volumes of geodesic spheres and balls.     

An analysis for a high‐order difference scheme for numerical solution to utt = A(x,t)uxx + F(x,t, u,ut, ux)

This article is concerned with a high‐order implicit difference scheme presented by Mohanty, Jain, and George for the nonlinear hyperbolic equation u_tt = A(x, t)u_xx + F(x, t, u, u_t, u_x) with Dirichlet boundary conditions. Some prior estimates of the difference solution are obtained by the energy methods. The solvability of the difference scheme is proved by the energy method and Brower's fixed point theorem. Similarly, the uniqueness, the convergence in L_∞‐norm and the stability of the difference solution are obtained. A numerical example is provided to demonstrate the validity of the theoretical results. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 484–498, 2007     

A method for the resolution of the Jacobi equation <Emphasis Type="Italic">Y</Emphasis>″ + <Emphasis Type="Italic">RY</Emphasis> = 0 on the manifold <Emphasis Type="Italic">Sp</Emphasis>(2)/<Emphasis Type="Italic">SU</Emphasis>(2)

In this paper a method for the resolution of the differential equation of the Jacobi vector fields in the manifold V ₁ = Sp(2)/SU(2) is exposed. These results are applied to determine areas and volumes of geodesic spheres and balls. Work partially supported by DGI (Spain) and FEDER Projects MTM 2004-06015-C02-01 and MTM 2007-65852 (first author) and by Research Project PGIDIT05PXIB16601PR (second author). Authors’ addresses: A. M. Naveira, Departamento de Geometría y Topología. Facultad de Matemáticas, Avda. Andrés Estellés, N1, 46100 – Burjassot, Valencia, Spain; A. D. Tarrío Tobar, E. U. Arquitectura Técnica, Campus A Zapateira. Universidad de A Coru?a, 15192 – A Coru?a, Spain