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.     

Relations among eigenvalues of Sturm-Liouville problems with different types of leading coefficient functions

For any Sturm-Liouville problem with a separable boundary condition and whose leading coefficient function changes sign (exactly once), we first give a geometric characterization of its eigenvalues λn using the eigenvalues of some corresponding problems with a definite leading coefficient function. Consequences of this characterization include simple proofs of the existence of the λn's, their Prüfer angle characterization, and a way for determining their indices from the zeros of their eigenfunctions. Then, interlacing relations among the λn's and the eigenvalues of the corresponding problems are obtained. Using these relations, a simple proof of asymptotic formulas for the λn's is given.     

The solution existence of general variational inclusion problems

We propose general variational inclusion problems which are slightly different from corresponding problems considered in several recent papers in the literature and show that they are advantageous. Sufficient conditions for the solution existence are established. As applications we derive consequences for several special cases of variational inclusion problems, quasioptimization problems, equilibrium problems and implicit variational inequalities and show that they improve the results of some recent existing papers.     

Fixed-sign solutions for a system of singular focal boundary value problems

We consider a system of focal boundary value problems where the nonlinearities may be singular in the independent variable and may also be singular in the dependent arguments. Using Schauder fixed point theorem, we establish criteria such that the system of boundary value problems has at least one fixed-sign solution.     

Patterns in parabolic problems with nonlinear boundary conditions

We obtain existence of asymptotically stable nonconstant equilibrium solutions for semilinear parabolic equations with nonlinear boundary conditions on small domains connected by thin channels. We prove the convergence of eigenvalues and eigenfunctions of the Laplace operator in such domains. This information is used to show that the asymptotic dynamics of the heat equation in this domain is equivalent to the asymptotic dynamics of a system of two ordinary differential equations diffusively (weakly) coupled. The main tools employed are the invariant manifold theory and a uniform trace theorem.     

Identification of a dissipation coefficient by a variational method

A generalized inverse problem for the identification of the absorption coefficient for a hyperbolic system is considered. The well-posedness of the problem is examined. It is proved that the regular part of the solution is an L ₂ function, which reduces the inverse problem to minimizing the error functional. The gradient of the functional is determined in explicit form from the adjoint problem, and approximate formulas for its calculation are derived. A regularization algorithm for the solution of the inverse problem is considered. Numerical results obtained for various excitation sources are displayed.     

Numerical solution of differential equations using Sinc method based on the interpolation of the highest derivatives

In general, we will use the numerical differentiation when dealing with the differential equations. Thus the differential equations can be transformed into algebraic equations and then we can get the numerical solutions. But as we all have known, the numerical differentiation process is very sensitive to even a small level of errors. In contrast it is expected that on average the numerical integration process is much less sensitive to errors. In this paper, based on the Sinc method we provide a new method using Sinc method incorporated with the double exponential transformation based on the interpolation of the highest derivatives (SIHD) for the differential equations. The error in the approximation of the solution is shown to converge at an exponential rate. The numerical results show that compared with the exiting results, our method is of high accuracy, of good convergence with little computational efforts. It is easy to treat nonhomogeneous mixed boundary condition for our method, which is unlike the traditional Sinc method.     

Multisplitting Iteration Schemes for Solving a Class of Nonlinear Complementarity Problems

We consider several synchronous and asynchronous multisplitting iteration schemes for solving aclass of nonlinear complementarity problems with the system matrix being an H-matrix.We establish theconvergence theorems for the schemes.The numerical experiments show that the schemes are efficient forsolving the class of nonlinear complementarity problems.     

四阶微分方程奇异边值问题的正解

本文利用锥上的不动点指数理论,在更一般的边界条件下讨论了一类四阶奇异边值问题正解的存在性并给出了应用．     

Noncooperative facility location games

A noncooperative game theoretical approach is considered for the multifacility location problem. It turns out that the facility location game is a potential game in the sense of Monderer and Shapley and some properties of the game are studied.