Abstract:Stochastic approximation (SA) and stochastic gradient descent (SGD) algorithms are work-horses for modern machine learning algorithms. Their constant stepsize variants are preferred in practice due to fast convergence behavior. However, constant step stochastic iterative algorithms do not converge asymptotically to the optimal solution, but instead have a stationary distribution, which in general cannot be analytically characterized. In this work, we study the asymptotic behavior of the appropriately scaled stationary distribution, in the limit when the constant stepsize goes to zero. Specifically, we consider the following three settings: (1) SGD algorithms with smooth and strongly convex objective, (2) linear SA algorithms involving a Hurwitz matrix, and (3) nonlinear SA algorithms involving a contractive operator. When the iterate is scaled by $1/\sqrt{\alpha}$, where $\alpha$ is the constant stepsize, we show that the limiting scaled stationary distribution is a solution of an integral equation. Under a uniqueness assumption (which can be removed in certain settings) on this equation, we further characterize the limiting distribution as a Gaussian distribution whose covariance matrix is the unique solution of a suitable Lyapunov equation. For SA algorithms beyond these cases, our numerical experiments suggest that unlike central limit theorem type results: (1) the scaling factor need not be $1/\sqrt{\alpha}$, and (2) the limiting distribution need not be Gaussian. Based on the numerical study, we come up with a formula to determine the right scaling factor, and make insightful connection to the Euler-Maruyama discretization scheme for approximating stochastic differential equations.

Recursive And Iterative Stochastic Gradient Algorithms Based On Two-Step Update Estimation For Wiener Systems

Online Successive Convex Approximation for Two-Stage Stochastic Nonconvex Optimization

Two‐stage gradient‐based iterative algorithms for the fractional‐order nonlinear systems by using the hierarchical identification principle

Parameter Estimation of Wiener Systems Based on the Particle Swarm Iteration and Gradient Search Principle

Recursive estimation for nonlinear stochastic systems with multi-step transmission delays, multiple packet dropouts and correlated noises

Two‐stage auxiliary model gradient‐based iterative algorithm for the input nonlinear controlled autoregressive system with variable‐gain nonlinearity

Hierarchical Stochastic Gradient Algorithm and its Performance Analysis for a Class of Bilinear-in-Parameter Systems

Particle Swarm Optimization Iterative Identification Algorithm and Gradient Iterative Identification Algorithm for Wiener Systems with Colored Noise

A Successive over Relaxation Implicit Iterative Algorithm for Solving Stochastic Linear Systems with Markov Jumps

Continuous-time stochastic gradient descent for optimizing over the stationary distribution of stochastic differential equations

Two-stage Gradient-based Iterative Estimation Methods for Controlled Autoregressive Systems Using the Measurement Data

Recursive identification methods for general stochastic systems with colored noises by using the hierarchical identification principle and the filtering identification idea

Stationary Behavior of Constant Stepsize SGD Type Algorithms: An Asymptotic Characterization

Two‐stage and three‐stage recursive gradient identification of Hammerstein nonlinear systems based on the key term separation

Stochastic gradient identification of Wiener system with maximum mutual information criterion

Escaping Saddle Points with Stochastically Controlled Stochastic Gradient Methods

Stochastic Modified Equations and Adaptive Stochastic Gradient Algorithms

A Class of Convex Optimization-Based Recursive Algorithms for Identification of Stochastic Systems

Stochastic Point-to-Point Iterative Learning Tracking Without Prior Information on System Matrices

The modified extended Kalman filter based recursive estimation for Wiener nonlinear systems with process noise and measurement noise

Asymptotic error analysis for stochastic gradient optimization schemes with first and second order modified equations