Markovian Foundations for Quasi-Stochastic Approximation with Applications to Extremum Seeking Control

TL;DR

This paper develops a Markovian framework for quasi-stochastic approximation, providing new theoretical insights, error bounds, and stability results applicable to optimization and extremum seeking control, especially under non-Lipschitz conditions.

Contribution

It introduces a novel Markovian foundation for QSA, derives an exact ODE representation, and extends stability and error analysis to non-Lipschitz algorithms in extremum seeking.

Findings

Error bound of order O(α) reduced to O(α^2) with filters

New stability results for non-Lipschitz extremum seeking algorithms

Potential for error bounds better than O(α) with Markovian noise

Abstract

This paper concerns quasi-stochastic approximation (QSA) to solve root finding problems commonly found in applications to optimization and reinforcement learning. The general constant gain algorithm may be expressed as the time-inhomogeneous ODE $ \frac{d}{dt}\Theta_t=\alpha f_t (\Theta_t)$, with state process $\Theta$ evolving on $\mathbb{R}^d$. Theory is based on an almost periodic vector field, so that in particular the time average of $f_t(\theta)$ defines the time-homogeneous mean vector field $\bar{f} \colon \mathbb{R}^d \to \mathbb{R}^d$ with $\bar{f}(\theta^*)=0$. Under smoothness assumptions on the functions involved, the following exact representation is obtained: \[\frac{d}{dt}\Theta_t=\alpha[\bar{f}(\Theta_t)-\alpha\bar\Upsilon_t+\alpha^2\mathcal{W}_t^0+\alpha\frac{d}{dt}\mathcal{W}_t^1+\frac{d^2}{dt^2}\mathcal{W}_t^2]\] along with formulae for the smooth signals $\{\bar \Upsilon_t , \mathcal{W}_t^i : i=0, 1, 2\}$. This new representation, combined with new conditions for ultimate boundedness, has many applications for furthering the theory of QSA and its applications, including the following implications that are developed in this paper: (i) A proof that the estimation error $\|\Theta_t-\theta^*\|$ is of order $O(\alpha)$, but can be reduced to $O(\alpha^2)$ using a second order linear filter. (ii) In application to extremum seeking control, it is found that the results do not apply because the standard algorithms are not Lipschitz continuous. A new approach is presented to ensure that the required Lipschitz bounds hold, and from this we obtain stability, transient bounds, and asymptotic bias of order $O(\alpha^2)$, and asymptotic variance of order $O(\alpha^4)$. (iii) It is in general possible to obtain better than $O(\alpha)$ bounds on error in traditional stochastic approximation when there is Markovian noise.