A stochastic algorithm using one sample point per iteration and diminishing stepsizes

A stochastic algorithm for finding stationary points of real-valued functions defined on a Euclidean space is analyzed. It is based on the Robbins-Monro stochastic approximation procedure. Gradient evaluations are done by means of Monte Carlo simulations. At each iteratex _i , one sample point is drawn from an underlying probability space, based on which the gradient is approximated. The descent direction is against the approximation of the gradient, and the stepsize is 1/i. It is shown that, under broad conditions, w.p.1 if the sequence of iteratesx ₁,x ₂,...generated by the algorithm is bounded, then all of its accumulation points are stationary.     

Sensitivity analysis for discrete optimal control problems

We present local sensitivity analysis for discrete optimal control problems with varying endpoints in the case when the customary regularity of boundary conditions can be violated. We study the behavior of the optimal solutions subject to parametric perturbations of the problem.     

Coupled intervals in the discrete calculus of variations: necessity and sufficiency

In this work we study nonnegativity and positivity of a discrete quadratic functional with separately varying endpoints. We introduce a notion of an interval coupled with 0, and hence, extend the notion of conjugate interval to 0 from the case of fixed to variable endpoint(s). We show that the nonnegativity of the discrete quadratic functional is equivalent to each of the following conditions: The nonexistence of intervals coupled with 0, the existence of a solution to Riccati matrix equation and its boundary conditions. Natural strengthening of each of these conditions yields a characterization of the positivity of the discrete quadratic functional. Since the quadratic functional under consideration could be a second variation of a discrete calculus of variations problem with varying endpoints, we apply our results to obtain necessary and sufficient optimality conditions for such problems. This paper generalizes our recent work in [R. Hilscher, V. Zeidan, Comput. Math. Appl., to appear], where the right endpoint is fixed.     

Explicit Formulas for Generalized Action–Angle Variables in a Neighborhood of an Isotropic Torus and Their Application

Different versions of the Darboux–Weinstein theorem guarantee the existence of action–angle-type variables and the harmonic-oscillator variables in a neighborhood of isotropic tori in the phase space. The procedure for constructing these variables is reduced to solving a rather complicated system of partial differential equations. We show that this system can be integrated in quadratures, which permits reducing the problem of constructing these variables to solving a system of quadratic equations. We discuss several applications of this purely geometric fact in problems of classical and quantum mechanics.     

Virtual logarithms

The modification of the number field sieve which Joux and Lercier recently used to compute logarithms in a prime field of a record 120 decimal digits makes use of the notion of a virtual logarithm of a prime ideal in a number ring. We provide necessary and sufficient conditions for their method to succeed and give an explicit formula for the virtual logarithm of an ideal.     

On the counting function of the lattice profile of periodic sequences

The lattice profile analyzes the intrinsic structure of pseudorandom number sequences with applications in Monte Carlo methods and cryptology. In this paper, using the discrete Fourier transform for periodic sequences and the relation between the lattice profile and the linear complexity, we give general formulas for the expected value, variance, and counting function of the lattice profile of periodic sequences with fixed period. Moreover, we determine in a more explicit form the expected value, variance, and counting function of the lattice profile of periodic sequences for special values of the period.     

The method of lines solution of discrete ordinates method for radiative heat transfer in cylindrical enclosures

A radiation code based on method of lines solution of discrete ordinates method for radiative heat transfer in axisymmetric cylindrical enclosures containing absorbing-emitting medium was developed and tested for predictive accuracy by applying it to (i) test problems with black and grey walls (ii) a gas turbine combustor simulator enclosing a non-homogeneous absorbing-emitting medium and benchmarking its steady-state predictions against exact solutions and measurements. Comparisons show that it provides accurate solutions for radiative heat fluxes and can be used with confidence in conjunction with CFD codes based on the same approach.     

多维Lotka-Uolterra周期系统周期解的存在性及其全局吸引性

研究了m-维离散Lotka-Volterra竞争系统，并获得了该系统周期解的存在性及其全局吸引性.     

X波段六腔渡越管振荡器的高频特性研究

 从麦克斯韦方程出发，采用时域有限差分法（FDTD）和离散傅里叶变换（DFFT）相结合的方法，通过数值计算得出了六腔开放式谐振腔中前四个谐振频率和场分布，计算出的谐振频率与实验测量结果基本相同。比较了开放腔和封闭腔谐振频率，验证了TEM波吸收边界条件，并在实际编程计算中得以应用。计算结果为六腔渡越管振荡器的机理研究提供了依据。     

A Method for Establishing Optimality Conditions for a Nonsmooth Vector-Valued Minimax Problem

Using a method of uniform approximations, necessary and sufficient conditions for a nonsmooth constrained vector-valued minimax problem are established in terms of Mordukhovich subdifferentials.