Abstract
We incorporate future information in the form of the estimated value offuture gradients in online convex optimization. This is motivated by demandresponse in power systems, where forecasts about the current round, e.g., theweather or the loads' behavior, can be used to improve on predictions made withonly past observations. Specifically, we introduce an additional predictivestep that follows the standard online convex optimization step when certainconditions on the estimated gradient and descent direction are met. We showthat under these conditions and without any assumptions on the predictabilityof the environment, the predictive update strictly improves on the performanceof the standard update. We give two types of predictive update for variousfamily of loss functions. We provide a regret bound for each of our predictiveonline convex optimization algorithms. Finally, we apply our framework to anexample based on demand response which demonstrates its superior performance toa standard online convex optimization algorithm.
Quick Read (beta)
Predictive Online Convex Optimization\thanksreffootnoteinfo
Abstract
We incorporate future information in the form of the estimated value of future gradients in online convex optimization. This is motivated by demand response in power systems, where forecasts about the current round, e.g., the weather or the loads’ behavior, can be used to improve on predictions made with only past observations. Specifically, we introduce an additional predictive step that follows the standard online convex optimization step when certain conditions on the estimated gradient and descent direction are met. We show that under these conditions and without any assumptions on the predictability of the environment, the predictive update strictly improves on the performance of the standard update. We give two types of predictive update for various family of loss functions. We provide a regret bound for each of our predictive online convex optimization algorithms. Finally, we apply our framework to an example based on demand response which demonstrates its superior performance to a standard online convex optimization algorithm.
AND
ucb]Antoine LesageLandry, is]Iman Shames, ut]Joshua A. Taylor
${}^{\mathrm{a}}$Energy & Resources Group, University of California, Berkeley, USA ${}^{}$
${}^{\mathrm{b}}$Department of Electrical and Electronic Engineering, University of Melbourne, Parkville, Australia ${}^{}$
${}^{\mathrm{c}}$The Edward S. Rogers Sr. Department of Electrical and Computer Engineering, University of Toronto, Toronto, Canada ${}^{}$
Key words: Convex optimization; learning algorithms; machine learning; power systems; renewable energy systems; load dispatching
1 Introduction
Online convex optimization (OCO) has found applications in fields like network resource allocation [8, 7, 6] and demand response in power systems [20, 21]. It is used for sequential decisionmaking when contextual information or feedback is only revealed to the decision maker at the end of the current round. Theoretical results showing that OCO algorithms have bounded regret guarantee the performance of these algorithms under mild assumptions.
In many applications, the decision maker has access to both revealed past information and estimates about future rounds. For example, in power systems, weather forecasts or historical load patterns can be used to estimate the future regulation needs [22, 4]. In this work, we present the predictive online convex optimization (POCO) framework. POCO works under the assumption that an estimate of the gradient of the loss function for the next round is available to the decision maker. In POCO, a standard OCO update is first applied using past information to compute the next decision. Then, the decision maker checks the quality of the estimated information available to them. If the estimated gradient is considered accurate enough, the decision maker implements an additional projected gradient step based on the estimated gradient to improve their decision for this round. This last step is referred as the predictive update.
We introduce explicit criteria for determining if the quality of the estimated gradient is high enough to guarantee an improvement over a standard OCO step when the predictive update is applied. A regret bound is obtained for all our algorithms. We conclude this work by presenting numerical examples where a POCO algorithm is used to improve on the performance of demand response with standard OCO. This example is motivated by the fact that a load aggregator often has access to an estimate of the power imbalance they have to counteract for regulation purposes.
Literature review. Recent work in online convex optimization has focused on including prior or future information. Reference [28], which builds on [9], assumes that the problem’s unknown and uncertain parameters follow a predictable process plus some noise [27] for their OCO algorithm. As in our setting, a second update with an estimated gradientlike term follows a mirror descent update. This second update is used by the algorithm in every step regardless of the quality of the estimated gradient. For this reason, the algorithm is referred to as optimistic. Optimistic algorithms were also studied in [31, 23, 34]. No conditions are provided about the estimated gradient in this case except that it comes from past observations and/or side information via an oracle. The authors of [28] show that the optimistic mirror descent can lead to a tighter bound than a standard online mirror descent algorithm if the process is indeed predictable. In [19], the authors provide a dynamic regret bound for the optimistic mirror descent. There is, however, no guarantee that in a given round the optimistic update does not do worse than the standard OCO update. An algorithm similar to [28] is given in [18]. In their work, they make the stronger assumption in which the exact gradient of the next round loss function is available and then provide a static regret bound for their setting. This differs from our setting in that we provide dynamic regretbounded algorithms and use an estimated gradient which entails less restrictive assumptions. Several other authors have studied different ways to incorporate future information in OCO like using state information [17] or the direction of the loss function’s gradient in an online linear optimization setting [10].
The projected gradient descent, inexact gradient descent, and proximal algorithms [1, 29, 2] from conventional convex optimization resemble our setting. These algorithms differ from ours because they aim to minimize the same objective function throughout all descent steps. In OCO, we minimize a sequence of objective functions ${\left\{{f}_{t}\right\}}_{t=1}^{T}$ and at each time $t$ provide a decision to minimize the current loss function. The loss function in a given round is only observed after we have committed to a decision. OCO will be introduced formally in Section 2.
Model predictive control (MPC) [14, 3] is another widelyused sequential decisionmaking framework. In MPC, the decision maker solves to optimality a receding horizon optimization problem that relies on models of future round loss functions. This thus requires significantly more contextual information and computational resources. These limitations are absent in OCO, making it a more suitable tool for realtime decision making with small computational resources.
Because we characterize conditions under which the predictive step improves performance, we guarantee improvement over conventional OCO and require no predictability assumptions. These conditions can be checked at each round of OCO, and if satisfied, the predictive update is implemented. In sum, in this work we make the following contributions:

$\bullet $
We introduce a novel predictive online convex optimization framework and provide conditions for when to use side information.

$\bullet $
We propose a predictive update with a predetermined step size for loss functions that have a Lipschitz gradient. We show that this update leads to a strict improvement over an OCO update when used (Section 4).

$\bullet $
We give a predictive update with backtracking line search that applies to a broader family of problems. We show that it leads to strict improvement over an OCO update (Section 5).

$\bullet $
We obtain sublinear regret bounds in the number of rounds for all algorithms.

$\bullet $
We apply our framework to demand response in power systems and find that it outperforms a standard OCO algorithm (Section 6).
2 Background
In OCO, one must make a decision at each round to minimize their cumulative loss [30, 15]. The current round’s loss function and any other rounddependent parameters are not available at the moment when the decision is made. Only information about previous rounds can be used to make the decision. Once the decision has been made, information about the current round is observed.
Let $t$ denote the current round index and $T$ be the time horizon. Let $\mathcal{X}\subset {\mathbb{R}}^{N}$, $N\in \mathbb{N}$, be the decision set, and let ${\mathbf{x}}_{t}\in \mathcal{X}$ be the decision variable at time $t$. We denote the differentiable convex loss function by ${f}_{t}({\mathbf{x}}_{t})$ for $t=1,2,\mathrm{\dots},T$. Let $\parallel \cdot \parallel $ be the Euclidean norm. We denote the projection operator onto the set $\mathcal{Y}$ as ${\mathrm{proj}}_{\mathcal{Y}}(\mathbf{x})\in {\mathrm{arg}\mathrm{min}}_{\mathbf{y}\in \mathcal{Y}}\parallel \mathbf{x}\mathbf{y}\parallel $.
The goal of the decision maker is to sequentially solve the following sequence of problems:
$$\underset{{\mathbf{x}}_{t}\in \mathcal{X}}{\mathrm{min}}{f}_{t}({\mathbf{x}}_{t})$$  (1) 
for $t=1,2,\mathrm{\dots},T$. The decision maker observes the loss function ${f}_{t}$ after choosing ${\mathbf{x}}_{t}$. For this reason, even if the loss function has a simple form, an analytical solution to the round optimization problem (2) is not obtainable. The decision ${\mathbf{x}}_{t}$ is computed using a gradient descentbased [35], mirrored descentbased [13] or Newton stepbased rule [16]. For example, in the online gradient descent (OGD) [35] algorithm, the decision at round $t+1$, ${\mathbf{x}}_{t+1}$, is given by the update:
$${\mathbf{x}}_{t+1}={\mathrm{proj}}_{\mathcal{X}}\left({\mathbf{x}}_{t}\eta \nabla {f}_{t}({\mathbf{x}}_{t})\right),$$  (2) 
where $\eta \propto {T}^{1/2}$ to guarantee a sublinear upper bound on the dynamic regret of OGD [35, Theorem 2].
Assumption 1
The set $\mathrm{X}$ is convex and compact.
The decision set $\mathcal{X}$ represents all constraints on ${\mathbf{x}}_{t}$. In this version of OCO, we only consider timeinvariant constraints.
Assumption 2
The loss function is $B$bounded: $\mathrm{}{f}_{t}\mathit{}\mathrm{(}{\mathrm{x}}_{t}\mathrm{)}\mathrm{}\mathrm{\le}B$ for $t\mathrm{=}\mathrm{1}\mathrm{,}\mathrm{2}\mathrm{,}\mathrm{\dots}\mathrm{,}T$ and $$.
Assumption 3
The gradient of the loss function is $G$bounded: $\mathrm{\parallel}\mathrm{\nabla}\mathit{}{f}_{t}\mathit{}\mathrm{(}{\mathrm{x}}_{t}\mathrm{)}\mathrm{\parallel}\mathrm{\le}G$ for $t\mathrm{=}\mathrm{1}\mathrm{,}\mathrm{2}\mathrm{,}\mathrm{\dots}\mathrm{,}T$ and $$.
As a consequence of Assumption 1, the decision variable is also $X$bounded: $\parallel {\mathbf{x}}_{t}\parallel \le X$ for $t=1,2,\mathrm{\dots},T$. We define the diameter of the compact set $\mathcal{X}$ as $\mathrm{diam}\mathcal{X}=sup\{\parallel \mathbf{x}\mathbf{y}\parallel \mathbf{x},\mathbf{y}\in \mathcal{X}\}$and let $D=\mathrm{diam}\mathcal{X}$, a positive scalar. The remainder of the assumptions will be stated when a specific technical result requires it.
The design tool of OCO algorithms is the regret [30, 15]. In this work, we use the dynamic regret [8, 19, 35, 24]:
$${\text{\U0001d681\U0001d68e\U0001d690}}_{T}^{d}=\sum _{t=1}^{T}{f}_{t}({\mathbf{x}}_{t}){f}_{t}({\mathbf{x}}_{t}^{\ast}),$$  (3) 
where ${\mathbf{x}}_{t}^{\ast}\in {\mathrm{arg}\mathrm{min}}_{\mathbf{x}\in \mathcal{X}}{f}_{t}(\mathbf{x})$. The dynamic regret compares the loss suffered by the decision maker to optimal performance in each round. Other versions of the regret exists, e.g., static regret [30, 15, 35], which is defined in terms of the optimal stationary decision, ${\mathbf{x}}^{\ast}\in {\mathrm{arg}\mathrm{min}}_{\mathbf{x}\in \mathcal{X}}{\sum}_{t=1}^{T}{f}_{t}(\mathbf{x})$ in (2). In this work, we only consider the dynamic regret because it yields a stronger theoretical guarantee. This theoretical guarantee is also more relevant in the context of timevarying optimization. For this reason, we refer to the dynamic regret, ${\text{\U0001d681\U0001d68e\U0001d690}}_{T}^{d}$, simply as the regret. Note that a bounded dynamic regret implies a bounded static regret [8]. The goal when designing an OCO algorithm is to show that the regret is sublinearly bounded above in the number of rounds. An OCO algorithm with a sublinearly bounded dynamic regret in the number of rounds will on average perform as well as the round optimal decision at each round [30, 15, 16].
We conclude this section by defining the quantity ${V}_{T}={\sum}_{t=2}^{T}\parallel {\mathbf{x}}_{t}^{\ast}{\mathbf{x}}_{t1}^{\ast}\parallel $. The term ${V}_{T}$ quantifies the variation of the optimal predictions through all rounds.
3 Predictive OCO
We now introduce our POCO framework. We let ${\overline{\mathbf{x}}}_{t+1}\in \mathcal{X}$ be the decision computed by an OCO algorithm in round $t$. This OCO algorithm can be, for example, the aforementioned OGD. The decision ${\overline{\mathbf{x}}}_{t+1}$ is then given by the update (2). In POCO, we consider an $\u03f5$forecaster introduced in Assumption 4. Let ${\mathbf{g}}_{t}({\overline{\mathbf{x}}}_{t+1})\in {\mathbb{R}}^{N}$ be the estimated gradient of the loss function ${f}_{t+1}$ at ${\overline{\mathbf{x}}}_{t+1}$.
Assumption 4 ($\u03f5$forecaster)
The $\u03f5$forecaster has access to an estimate of the gradient of the next round’s loss function evaluated at ${\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}$, ${\mathrm{g}}_{t}\mathit{}\mathrm{\left(}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{\right)}$, and the maximum estimation error, $\u03f5\mathrm{>}\mathrm{0}$, such that $\mathrm{\parallel}{\mathrm{g}}_{t}\mathit{}\mathrm{\left(}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{\right)}\mathrm{}\mathrm{\nabla}\mathit{}{f}_{t\mathrm{+}\mathrm{1}}\mathit{}\mathrm{(}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{)}\mathrm{\parallel}\mathrm{\le}\u03f5$ for all time $t\mathrm{=}\mathrm{1}\mathrm{,}\mathrm{2}\mathrm{,}\mathrm{\dots}\mathrm{,}T$.
In other words, we consider a forecaster that has access to limited information about the next round in the form of ${\mathbf{g}}_{t}\left({\overline{\mathbf{x}}}_{t+1}\right)$. This could represent a prediction based on historical data, e.g., a weather or demand forecast. The decision maker uses this information to improve on the OCO update. For conciseness, we denote the estimated gradient by ${\mathbf{g}}_{t}$. We omit its dependency on ${\overline{\mathbf{x}}}_{t+1}$ because it is always evaluated at the OCO update output, ${\overline{\mathbf{x}}}_{t+1}$, and no other points. The decision maker can meet Assumption 4 by relying on an exogenous model to estimate the gradient $\nabla {f}_{t+1}\left({\overline{\mathbf{x}}}_{t+1}\right)$. In the context of demand response, historical data of the load’s consumption and generator output’s patterns, weather history and the historical values of the gradient, for example, can be used to build a statistical model to estimate the value of the $\nabla {f}_{t+1}$ at the decision given by OCO update. The parameter $\u03f5$ can then be set according to, for example, a high confidence interval or a worstcase performance parameter. The forecaster would then provide ${\mathbf{g}}_{t}$ using this model.
Then, if certain conditions are met, the following update rule for our proposed POCO algorithm is used.
Definition 1 (Predictive update)
Let ${\beta}_{t}\mathrm{>}\mathrm{0}$ be an appropriately chosen step size. The predictive update is
$${\mathbf{x}}_{t+1}={\mathrm{proj}}_{\mathcal{X}}\left({\overline{\mathbf{x}}}_{t+1}{\beta}_{t}{\mathbf{g}}_{t}\right).$$  (4) 
The predictive update is to be used directly after the OCO update and will lead to a strict improvement over the OCO update under certain conditions. The aforementioned conditions will be discussed in the next sections and depend on the properties of the loss function. If the conditions are not met, ${\overline{\mathbf{x}}}_{t+1}$ is directly used. Let $\delta >0$ be the desired improvement when using the predictive update. We define the counter ${c}_{t}$:
$${c}_{t+1}=\{\begin{array}{cc}{c}_{t}+1\hfill & \text{if}\parallel {\mathbf{x}}_{t+1}{\overline{\mathbf{x}}}_{t+1}\parallel \ge \delta \hfill \\ {c}_{t}\hfill & \text{otherwise}\hfill \end{array}$$ 
with ${c}_{0}=0$. The variable ${c}_{t}$ represents the number of predictive updates as described in Definition 1. Let $\nu ={c}_{T}/T$ be the ratio of rounds using the predictive update to the total number of rounds.
Depending on the loss function, any regretbounded OCO update can be used in the POCO framework. Back to the OGD example, the predictive OGD uses the update (2) and if certain conditions are met, ${\mathbf{x}}_{t+1}$ is provided by (1) and if not, ${\mathbf{x}}_{t+1}={\overline{\mathbf{x}}}_{t+1}$. We write ${\mathbf{x}}_{t+1}({\beta}_{t})={\mathrm{proj}}_{\mathcal{X}}\left({\overline{\mathbf{x}}}_{t+1}{\beta}_{t}{\mathbf{g}}_{t}\right)$ as a function of the step size ${\beta}_{t}>0$ and let ${\mathbf{d}}_{t+1}={\mathbf{x}}_{t+1}({\beta}_{t}){\overline{\mathbf{x}}}_{t+1}$ be the descent direction.
Next, we provide sufficient conditions for the estimated gradient ${\mathbf{g}}_{t}$ to be a feasible descent direction. Later, we consider the step size selection problem. Particularly, two cases are considered where (i) the step sizes are constant and chosen a priori based on a property of the loss functions, or (ii) the step sizes are selected through the application of a backtracking line search that enforces a modified online version of the Armijo condition [2].
The following lemma introduces a sufficient condition for the estimated gradient ${\mathbf{g}}_{t}$ to be a descent direction of the OCO problem (2).
Lemma 1 (Estimated descent direction)
The vector $\mathrm{}{\mathrm{g}}_{t}$ provided by the $\u03f5$forecaster is a descent direction for ${f}_{t\mathrm{+}\mathrm{1}}\mathit{}\mathrm{(}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{)}$ if $\mathrm{\parallel}{\mathrm{g}}_{t}\mathrm{\parallel}\mathrm{>}\u03f5$.
The proof of Lemma 1 is presented in Appendix A. The next lemma is adapted from [2] and ensures that the predictive step follows a feasible descent direction.
Lemma 2 (Feasible estimated descent direction)
For all ${\beta}_{t}\mathrm{>}\mathrm{0}$ and ${\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{\in}\mathrm{X}$, if $\mathrm{\parallel}{\mathrm{g}}_{t}\mathrm{\parallel}\mathrm{>}\u03f5$ and ${\mathrm{x}}_{t\mathrm{+}\mathrm{1}}\mathit{}\mathrm{(}{\beta}_{t}\mathrm{)}\mathrm{\ne}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}$ , then ${\mathrm{x}}_{t\mathrm{+}\mathrm{1}}\mathit{}\mathrm{(}{\beta}_{t}\mathrm{)}\mathrm{}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}$ is a feasible descent direction at ${\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}$ and ${\mathrm{g}}_{t}^{\mathrm{\top}}\mathit{}\mathrm{\left(}{\mathrm{x}}_{t\mathrm{+}\mathrm{1}}\mathit{}\mathrm{(}{\beta}_{t}\mathrm{)}\mathrm{}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{\right)}\mathrm{\le}\mathrm{}\frac{\mathrm{1}}{{\beta}_{t}}\mathit{}{\mathrm{\parallel}{\mathrm{x}}_{t\mathrm{+}\mathrm{1}}\mathit{}\mathrm{(}{\beta}_{t}\mathrm{)}\mathrm{}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{\parallel}}^{\mathrm{2}}\mathrm{.}$
4 POCO with fixed step size
We now present a predictive update where step sizes ${\beta}_{t}$ are fixed and based on a propriety of the sequence of loss functions. We conclude this section by providing regretbounded algorithms using these updates. In this section, we add the following assumption:
Assumption 5
Let $$. The loss function ${f}_{t}\mathit{}\mathrm{(}\mathrm{x}\mathrm{)}$ has an $L$uniformly Lipschitzcontinuous gradient: $\mathrm{\parallel}\mathrm{\nabla}\mathit{}{f}_{t}\mathit{}\mathrm{(}\mathrm{x}\mathrm{)}\mathrm{}\mathrm{\nabla}\mathit{}{f}_{t}\mathit{}\mathrm{(}\mathrm{y}\mathrm{)}\mathrm{\parallel}\mathrm{\le}L\mathit{}\mathrm{\parallel}\mathrm{x}\mathrm{}\mathrm{y}\mathrm{\parallel}$ for all $t\mathrm{=}\mathrm{1}\mathrm{,}\mathrm{2}\mathrm{,}\mathrm{\dots}\mathrm{,}T$ and $\mathrm{x}\mathrm{,}\mathrm{y}\mathrm{\in}\mathrm{X}$.
We propose a predictive update with fixed step size next. We state sufficient conditions that guarantee a strict improvement over an OCO update. These sufficient conditions can be checked at each round to determine if the estimated information is accurate enough, and therefore if the predictive update should be used in the current round.
Lemma 3 (Predictive update with fixed step size)
Suppose that Assumption 5 holds and $\mathrm{\parallel}{\mathrm{g}}_{t}\mathrm{\parallel}\mathrm{>}\u03f5$. If $\beta \mathrm{\le}\frac{\mathrm{1}}{L}$ and $\mathrm{\parallel}{\mathrm{d}}_{t\mathrm{+}\mathrm{1}}\mathrm{\parallel}\mathrm{=}\mathrm{\parallel}{\mathrm{proj}}_{\mathrm{X}}\mathit{}\mathrm{\left(}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{}{\beta}_{t}\mathit{}{\mathrm{g}}_{t}\mathrm{\right)}\mathrm{}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{\parallel}\mathrm{\ge}\frac{\u03f5}{L}\mathrm{+}\sqrt{\frac{{\u03f5}^{\mathrm{2}}}{{L}^{\mathrm{2}}}\mathrm{+}\frac{\mathrm{2}\mathit{}\delta}{L}}$, then the predictive update (1) used by the $\u03f5$forecaster strictly improves on the OCO update and the improvement is bounded below by $\delta \mathrm{>}\mathrm{0}$.
The proof of Lemma 3 is provided in Appendix C. We now present regret bounds for POCO algorithms. This algorithm uses the predictive update with fixed step size to improve the performance of OCO algorithms.
Theorem 1 (POCO regret bound)
Consider an OCO algorithm with a sublinear regret upper bound. Suppose that the forecaster uses the predictive update (1) only at rounds $t$ when the estimated gradient ${\mathrm{g}}_{t}$ and feasible descent direction ${\mathrm{d}}_{t\mathrm{+}\mathrm{1}}\mathrm{=}{\mathrm{proj}}_{\mathrm{X}}\mathit{}\mathrm{\left(}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{}{\beta}_{t}\mathit{}{\mathrm{g}}_{t}\mathrm{\right)}\mathrm{}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}$ satisfy the assumptions of Lemma 3. If the ratio of rounds satisfying these assumptions is greater than $\nu $, then the regret of the POCO algorithm is bounded above by
$${\text{Reg}}_{T}^{d}(POCO)\le {\text{Reg}}_{T}^{d}(OCO)T\nu \delta .$$ 
The proof of Theorem 1 is presented in Appendix D. This theorem leads to the following corollary which provides a regret bound for the OGD with predictive updates (POGD).
Corollary 1 ($O\left(\sqrt{T}\right)$ regret bound for POGD)
Suppose that the ratio $\nu $ of rounds that respects the assumptions of Lemma 3 is $\nu \mathrm{>}\frac{\mathrm{1}}{\sqrt{T}}$. Then the predictive OGD algorithm’s regret is bounded above by
${\text{Reg}}_{T}^{d}(\text{\U0001d67f\U0001d67e\U0001d676\U0001d673})$  $\le {\text{Reg}}_{T}^{d}(\text{\U0001d67e\U0001d676\U0001d673})\delta \sqrt{T},$  
$=\left({\displaystyle \frac{7{X}^{2}}{4}}+{\displaystyle \frac{{G}^{2}}{2}}+X{V}_{T}\delta \right)\sqrt{T},$ 
which is sublinear and tighter than the OGD regret bound.
5 POCO with backtracking line search
In this section, we do not require Assumption 5 to hold. We however use the following proposition:
Proposition 1
The loss function ${f}_{t}\mathit{}\mathrm{(}\mathrm{x}\mathrm{)}$ is $\mathrm{\Delta}$timeLipschitz with $$, that is:
$$\left{f}_{t}(\mathbf{x}){f}_{t+\tau}(\mathbf{x})\right\le {\mathrm{\Delta}}_{t}(\mathbf{x})\tau \le \mathrm{\Delta}\tau $$ 
for all $\tau \mathrm{\in}\mathrm{\left\{}i\mathrm{\in}\mathrm{Z}\mathrm{\right}\mathrm{0}\mathrm{\le}t\mathrm{+}i\mathrm{\le}T\mathrm{\}}$ at all $t\mathrm{=}\mathrm{1}\mathrm{,}\mathrm{2}\mathrm{,}\mathrm{\dots}\mathrm{,}T$, and all $\mathrm{x}\mathrm{\in}\mathrm{X}$.
Proposition 1 always holds because Assumption 2 implies that $\mathrm{\Delta}=2B/\left\tau \right$ is sufficient. Under Proposition 1, we consider functions that are $(t,\mathbf{x})$locally and globally Lipschitz in their time argument, respectively, for the intermediary bound (${\mathrm{\Delta}}_{t}(\mathbf{x})$) and the upper bound ($\mathrm{\Delta}$). This can represent, for example, loss functions like squared tracking error functions, in which the timevarying targets are always contained in a closed set.
In the case of the POCO with backtracking (POCOb), we reexpress the update (1) given in Definition 2. The backtracking line search for predictive update is given in Algorithm 2.
Definition 2 (POCOb update)
Let $\zeta $ be a positive scalar and ${\beta}^{m}$ be determined by a backtracking line search algorithm. The predictive update with backtracking line search is
$${\mathbf{x}}_{t+1}={\overline{\mathbf{x}}}_{t+1}+{\beta}^{m}\left({\mathrm{proj}}_{\mathcal{X}}\left({\overline{\mathbf{x}}}_{t+1}\zeta {\mathbf{g}}_{t}\right){\overline{\mathbf{x}}}_{t+1}\right)$$  (5) 
algorithm[tb] \[email protected]@algorithmic[1] \STATEParameters: Given $\beta \in (0,1)$ and $M\in \mathbb{N}$. \STATEInitialization: Set $\zeta >0$.
${\mathbf{d}}_{t+1}={\mathrm{proj}}_{\mathcal{X}}\left({\overline{\mathbf{x}}}_{t+1}\zeta {\mathbf{g}}_{t}\right){\overline{\mathbf{x}}}_{t+1}$ \STATE$m=0$. \WHILE${f}_{t}\left({\overline{\mathbf{x}}}_{t+1}+{\beta}^{m}{\mathbf{d}}_{t+1}\right)>{f}_{t}({\overline{\mathbf{x}}}_{t+1})+{\beta}^{m}\left({\mathbf{g}}_{t}^{\top}{\mathbf{d}}_{t+1}\u03f5\parallel {\mathbf{d}}_{t+1}\parallel \right)2\mathrm{\Delta}$ and $m\le M$ \STATE$m=m+1$.
$m>M$
$\beta =0$.
The next lemma shows that the backtracking line searchbased predictive update improves on the OCO update. Our claim relies on the modified Armijo condition for gradient projection. This condition ensures a sufficient decrease in the objective when using an estimated gradient projection descent direction [33]. We adapt this condition to the estimated gradient and online setting. The modified Armijo condition for gradient projection [2] on ${f}_{t+1}$ and feasible descent direction ${\mathbf{d}}_{t+1}={\mathrm{proj}}_{\mathcal{X}}\left({\overline{\mathbf{x}}}_{t+1}+\zeta {\mathbf{g}}_{t}\right){\overline{\mathbf{x}}}_{t+1}$ for some $\zeta >0$ with step size ${\beta}^{m}$ is given by:
${f}_{t+1}\left({\overline{\mathbf{x}}}_{t+1}+{\beta}^{m}{\mathbf{d}}_{t+1}\right)$  $\le {f}_{t+1}({\overline{\mathbf{x}}}_{t+1})$  (6)  
$\mathrm{\hspace{1em}}+{\beta}^{m}\nabla {f}_{t+1}{({\overline{\mathbf{x}}}_{t+1})}^{\top}{\mathbf{d}}_{t+1}.$ 
Lemma 4 (Sufficient decrease of POCOb update)
Suppose $\mathrm{\parallel}{\mathrm{g}}_{t}\mathrm{\parallel}\mathrm{>}\u03f5$. If Algorithm 2 terminates to a step size ${\beta}^{m}\mathrm{>}\mathrm{0}$, then the predictive update with backtracking line search (5) used by the $\u03f5$forecaster satisfies the modified Armijo condition (6), and leads to a sufficient decrease in the loss function.
The proof of the lemma can be found in Appendix E.
Remark 1
Algorithm 2 ensures that when $\beta \mathrm{\ne}\mathrm{0}$, ${\beta}^{m}$ satisfied:
${f}_{t}\left({\overline{\mathbf{x}}}_{t+1}+{\beta}^{m}{\mathbf{d}}_{t+1}\right)$  $\le {f}_{t}({\overline{\mathbf{x}}}_{t+1})+{\beta}^{m}{\mathbf{g}}_{t}^{\top}{\mathbf{d}}_{t+1}$  
$\mathrm{\hspace{1em}}{\beta}^{m}\u03f5\parallel {\mathbf{d}}_{t+1}\parallel 2\mathrm{\Delta}$  (7) 
Every element of (7) is available at time $t$, which is not the case in (6). This allows us to use a backtracking line search algorithm to determine ${\beta}_{t}$ in an OCO setting. Algorithm 2 also ensures that the step size is not too small (cf. [33, Section 3.1]).
Note that there is an additional $\u03f5\parallel {\mathbf{d}}_{t+1}\parallel $ term in the modified Armijo condition for estimated gradient projection. This is a consequence of not having access to the exact gradient of ${f}_{t}$. Hence, to ensure that the update is valid, the modified Armijo condition is augmented by a term proportional to the error of the estimated gradient. The second additional term, $2\mathrm{\Delta}$, is due to the timevarying setting of OCO.
We now discuss the existence of step sizes ${\beta}_{t}$ that satisfy (7) at round $t$. Before stating the main result, for a given ${\overline{\mathbf{x}}}_{t+1}$ and ${\mathbf{g}}_{t}$, define the set of step sizes that comply with line 5 in the line search algorithm, which is the modified Armijo condition for online settings (7):
$\mathcal{S}=\{\beta >0$  ${f}_{t}\left({\overline{\mathbf{x}}}_{t+1}+\beta {\mathbf{d}}_{t+1}\right)\le {f}_{t}({\overline{\mathbf{x}}}_{t+1})$  
$+\beta {\mathbf{g}}_{t}^{\top}{\mathbf{d}}_{t+1}\beta \u03f5\parallel {\mathbf{d}}_{t+1}\parallel 2\mathrm{\Delta}\}.$ 
Theorem 2
Suppose ${\mathrm{d}}_{t\mathrm{+}\mathrm{1}}\mathrm{=}{\mathrm{proj}}_{\mathrm{X}}\mathit{}\mathrm{\left(}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{}\zeta \mathit{}{\mathrm{g}}_{t}\mathrm{\right)}\mathrm{}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{\ne}\mathrm{0}$ is a feasible descent direction and ${f}_{t}$ is bounded below for all $t$. Then there exists $\mathrm{x}\mathrm{\in}\mathrm{X}$ such that ${f}_{t}\mathit{}\mathrm{(}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{)}\mathrm{}{f}_{t}\mathit{}\mathrm{(}\mathrm{x}\mathrm{)}\mathrm{>}\mathrm{2}\mathit{}\mathrm{\Delta}$ if and only if $\mathrm{S}\mathrm{\ne}\mathrm{\varnothing}$.
The proof for Theorem 2 is presented in Appendix F. We note that Theorem 2 does not guarantee that the backtracking algorithm, Algorithm 2, will find a nonzero step size. Other techniques like exact line searches, might be required to identify an adequate step size in some problem instances. Using Theorem 2, we can provide a lower bound on the improvement of the predictive update with backtracking line search.
Corollary 2 (POCOb update improvement)
Suppose that the assumptions of Lemma 4 hold and ${\beta}^{m}\mathrm{>}\mathrm{0}$, then the predictive update with backtracking line search improves on the OCO update by a minimum of $\mathrm{2}\mathit{}\mathrm{\Delta}$.
The proof of Corollary 2 is given in Appendix G. We now state a regret bound for the POCOb algorithm.
Theorem 3 (POCOb regret bound)
Consider an OCO algorithm with bounded regret. Suppose that the assumptions of Lemma 4 are met. If the ratio of rounds with $\beta \mathrm{>}\mathrm{0}$ and satisfying these assumptions to $T$ is greater than $\nu $, then the regret of the POCO algorithm with backtracking used by the $\u03f5$forecaster is bounded above by
$${\text{Reg}}_{T}^{d}(POCOb)\le {\text{Reg}}_{T}^{d}(OCO)2T\nu \mathrm{\Delta}$$  (8) 
and thus outperforms the OCO algorithm.
Remark 2
Note that if the locally Lipschitz statement of Proposition 1 is used, then $\mathrm{2}\mathit{}\mathrm{\Delta}$ is replaced by ${\mathrm{\Delta}}_{t\mathrm{,}\mathrm{1}}\mathit{}\mathrm{(}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{+}{\beta}^{m}\mathit{}{\mathrm{d}}_{t\mathrm{+}\mathrm{1}}\mathrm{)}\mathrm{+}{\mathrm{\Delta}}_{t\mathrm{,}\mathrm{1}}\mathit{}\mathrm{(}{\overline{\mathrm{x}}}_{t\mathrm{+}\mathrm{1}}\mathrm{)}$ in the modified Armijo condition for online settings (7), and the bound (8) can be recomputed accordingly.
6 Example
In this section, we apply POCO algorithms to demand response (DR) in power systems [5, 26], specifically regulation and curtailment. At each time step, a DR aggregator sends instructions to their loads to follow a regulation signal, e.g., a power imbalance due to a sudden change in renewable power generation [4, 32]. Each load responds to the signal by adjusting its power consumption. The power consumption is constrained by a storage capacity, which could represent physical storage like a battery or the load’s limits, e.g., thermal constraints. The regulation signal is unknown at the time the DR instructions are sent. This can be due, for example, to a drop in renewable power generation which is only assessed after the generator has committed to some amount of power. The objective of the DR aggregator is, therefore, to predict the DR dispatch at each time instance. This problem can be formulated as POCO, in which an estimate of the regulation signal is available to the load aggregator.
We consider $N$ loads. Let ${\mathbf{x}}_{t}\in {\mathbb{R}}^{N}$ denote the decision variable at round $t$. The variable ${\mathbf{x}}_{t}$ represents the instructions sent to the loads. Let ${r}_{t}\in \mathbb{R}$ be the regulation signal at time $t$. Let $\overline{\mathbf{x}},\underset{\xaf}{\mathbf{x}}\in {\mathbb{R}}^{N}$ be the maximum and minimum power that can be consumed or delivered for all loads. Define the decision set $\mathcal{X}=\left\{\mathbf{x}\in {\mathbb{R}}^{N}\right\underset{\xaf}{\mathbf{x}}\le \mathbf{x}\le \overline{\mathbf{x}}\}$. We let ${\mathbf{s}}_{t}\in {\mathbb{R}}^{N}$ denote the state of charge vectors of the loads at time $t$ and $\mathbf{c}\in {\mathbb{R}}^{N}$ the vector of load energy capacities. The state of charge of a load $i$ at time $t$ is ${s}_{t}(i)={s}_{0}(i)+{\sum}_{n=1}^{t}{x}_{n}(i)$. In the current case, we assume that there is no leakage nor energy losses. The OCO problem takes the following form:
$$\underset{{\mathbf{x}}_{t}\in \mathcal{X}}{\mathrm{min}}{({r}_{t}{\mathrm{\U0001d7cf}}^{\top}{\mathbf{x}}_{t})}^{2}+\sigma \parallel {\mathbf{s}}_{t1}+{\mathbf{x}}_{t}\frac{\mathbf{c}}{2}{\parallel}^{2}.$$  (9) 
The loss function has two terms: (i) a regulation term where the aggregated loads are dispatched to follow a regulation signal ${r}_{t}$ and (ii) a state of charge objective added to keep the loads near half their energy capacity. The loss function given in (9) is $\sigma $strongly convex. For this reason we use the OGD for strongly convex functions ($\sigma $OGD) proposed in [24], which offers tighter regret bound than the standard OGD. The following corollary gives an upper bound on the regret of predictive OGD for strongly convex function ($\sigma $POGD).
Corollary 3 (POGD for strongly convex functions)
Suppose ${f}_{t}$ is $\sigma $strongly convex and satisfies Assumption 5 for all $t$. Consider the $\sigma $OGD update
${\overline{\mathbf{x}}}_{t+1}$  $={\mathbf{x}}_{t}+\eta \left({\mathrm{proj}}_{\mathcal{X}}\left({\mathbf{x}}_{t}{\displaystyle \frac{1}{\gamma}}\nabla {f}_{t}({\mathbf{x}}_{t})\right){\mathbf{x}}_{t}\right)$ 
where $\eta \mathrm{\in}\mathrm{(}\mathrm{0}\mathrm{,}\mathrm{1}\mathrm{]}$ and $$. Then, the $\sigma $POGD with fixed step size, given that the assumptions of Lemma 3 hold for a ratio of the total rounds greater than $\nu $, has a regret bounded above by
${\text{Reg}}_{T}^{d}(\sigma \text{\U0001d67f\U0001d67e\U0001d676\U0001d673})$  $\le {\text{Reg}}_{T}^{d}(\sigma \text{\U0001d67e\U0001d676\U0001d673})T\nu \delta ,$  
$\le O\left({V}_{T}+1\right)T\nu \delta .$ 
The result follows from the proof of Theorem 1 and the $\sigma $OGD regret bound from [24]. We now present simulation results. All optimizations are solved using CVXPY [11] and the ECOS [12] solver.
Fixed step size example.
The load and algorithm parameters for this example are gathered in Table 1. The initial state of charge of each load is set to half its capacity. The regulation signal is ${r}_{t}=0.2\mathrm{sin}\left(\frac{2\pi}{T}t\right)+{w}_{t}$. The parameter ${w}_{t}\sim \mathrm{N}(0,0.01)$ is a Gaussian noise used to model sudden changes. We assume that the aggregator has access to estimated gradient for different level of accuracy $\u03f5$. This represents, for example when $\u03f5=0.01$, a relative error of at least $4\%$ of the actual gradient norm. The parameter $\sigma $ is set to achieve adequate regulation performance without deviating too much from each load’s desired state of charge.
Parameter  Value  Unit 

$N$  $25$  loads 
$h$  $30$  seconds 
$\overline{\mathbf{x}}/h$  $\mathrm{Uniform}[1,3]$  kW 
$\underset{\xaf}{\mathbf{x}}/h$  $\overline{\mathbf{x}}$  kW 
$\mathbf{c}$  $\mathrm{Uniform}[10,15]$  kWh 
$\u03f5$  $0.1$, $0.05$ & $0.01$  — 
$\delta $  ${10}^{6}$  — 
$\sigma $  $0.005$  — 
$\eta $  $1$  — 
$\gamma $  $L$  — 
$\beta $  $1/L$  — 
We now present the performance of our POCO algorithm with a fixed step size. We implement the OMD from [28] for comparison. This algorithm uses ${\mathbf{g}}_{t}$ without validating the estimated information. Figure 1 shows an instance of the experimental regret for the POCO with three different values of $\u03f5$, the conventional OCO algorithm, their respective regret bounds and OMD’s regret. POCO outperforms its bound, the OCO and the OMD algorithms. We remark that as expected the number of predictive updates increases with the accuracy of the estimated gradient, the performance of the POCO algorithm also improves.
Backtracking line search example.
We now present an example of POCO with backtracking. We consider a curtailment scenario. We let ${p}_{t}$ be the total power to be curtailed by the loads at time $t$ for $t=1,2,\mathrm{\dots},T$. When a contingency occurs in the network, flexible loads are called to curtail their power consumption, e.g., by temporarily shutting down their HVAC system. Contrary to the regulation case, the loads are not contracted to follow a setpoint and no penalties are assessed on loads curtailing more than asked. Similar to the regulation setting, the curtailment signal is unknown until immediately after the current round. This setting can be modeled as POCO where an estimated curtailment signal is available to the aggregator at each round. We use the same notation as the previous examples. Let ${[\cdot ]}^{+}=\mathrm{max}\{0,\cdot \}$. This curtailment scenario is modeled by loss function given below:
$${f}_{t}({\mathbf{x}}_{t})={\left({\left[{p}_{t}{\mathrm{\U0001d7cf}}^{\top}{\mathbf{x}}_{t}\right]}^{+}\right)}^{2}+\sigma {\parallel \alpha {\mathbf{s}}_{t1}+{\mathbf{x}}_{t}\frac{\mathbf{c}}{2}\parallel}^{2}$$  (10) 
where we have added a recovery coefficient to the state of charge objective term used previously. This coefficient models the usual evolution of the load (e.g., ambient temperature heating for a thermostatic load). We let $\alpha =1.001$. This is equivalent to a recovery coefficient of $1.13$ per hour. The function ${f}_{t}$ given in (10) is not gradient Lipschitz and Assumption 5 does not hold. We model the curtailment signal to be quickly increasing at first and then slowly plateauing to represent new level of available generation. This event is assumed to be limited in time, after which the network goes back to its normal state and no curtailment is then required. We let ${p}_{t}=0.04{t}^{0.3}+{w}_{t}$ where ${w}_{t}\sim \mathrm{N}(0,0.01)$ for $t=1,2,\mathrm{\dots},T/4$ and then ${p}_{t}=0.04{\left(T/4\right)}^{0.3}+{w}_{t}^{\prime}$ where ${w}_{t}^{\prime}\sim \mathrm{N}(0,0.001)$ for $t=T/4,T/4+1,\mathrm{\dots},T$. The noise variance is equivalent to approximatively $10\%$ of curtailment signal’s value at first and then about $1\%$.
Parameter  Value 

$\alpha $  $1.001$ 
$\u03f5$  $0.1$, $0.01$ & $0.001$ 
$M$  $100$ 
$\sigma $  $5\times {10}^{5}$ 
$\zeta $  $0.5$ 
$\beta $  $0.9$ 
$\eta $  $1/10\sqrt{T}$ 
We use the same parameters as in the previous section, except for those in Table 2. The POCOb experimental regret shown in Figure 2 is sublinear in the numbers of rounds and outperforms the OCO’s regret. While the performance is not as strong as POCO with fixed step size, this algorithm can be applied to a broader family of functions because it does not require the loss function to be gradient Lipschitz continuous. The difference in performances between the two POCO updates is explained by the fact that the sufficient conditions for the POCOb are rarely satisfied in this simulation. Improved estimation accuracy, $\u03f5$, and a loss function with a lower maximum temporal change, $\mathrm{\Delta}$, could increase the number of times the backtracking line search is used. Nevertheless, the POCOb achieves a regret reduction of 29% when $\u03f5=1\%$ over a standard OCO algorithm. Lastly, we note that POCOb performs better for larger variations in ${p}_{t}$ and smaller values of $\u03f5$.
7 Conclusion
In this work, we have presented the predictive online convex optimization framework. In POCO, a second update is used after the OCO update to improve performance using an estimated gradient. We have presented two versions of the predictive update that can be used under different assumptions. We have shown a regret upper bound for all of our POCO algorithms. We have applied POCO to demand response in electric power systems and found that they outperform conventional OCO using commonly available forecast information. In the case of fixed step size update, we observed an improvement of $95\%$ in the final regret and of $29\%$ in the backtracking case when having access to a $(\u03f5=0.01)$forecaster.
Acknowledgements
A.L.L. thanks P. Mancarella for his support and for cohosting him at The University of Melbourne.
A Proof of Lemma 1
Define ${\mathbf{e}}_{t}\in {\mathbb{R}}^{n}$ as ${\mathbf{e}}_{t}={\mathbf{g}}_{t}\nabla {f}_{t+1}({\overline{\mathbf{x}}}_{t+1})$ where $\parallel {\mathbf{e}}_{t}\parallel \le \u03f5$ by the definition of the $\u03f5$forecaster. A vector $\mathbf{v}$ is a descent direction if $$. Thus, we have
$0$  $$  
$={\mathbf{g}}_{t}^{\top}({\mathbf{g}}_{t}{\mathbf{e}}_{t})$ 
Equivalently, we have ${\mathbf{g}}_{t}^{\top}{\mathbf{g}}_{t}>{\mathbf{g}}_{t}^{\top}{\mathbf{e}}_{t}$. Taking the norm of both sides and dividing by the norm of ${\mathbf{g}}_{t}$ gives
$$\parallel {\mathbf{g}}_{t}\parallel >\parallel {\mathbf{e}}_{t}\parallel \mathrm{cos}{\theta}_{{\mathbf{g}}_{t},{\mathbf{e}}_{t}},$$  (A.1) 
where ${\theta}_{{\mathbf{g}}_{t},\mathbf{e}}$ is the angle between ${\mathbf{g}}_{t}$ and ${\mathbf{e}}_{t}$. By assumption, $\parallel {\mathbf{g}}_{t}\parallel >\u03f5$ and $\parallel {\mathbf{e}}_{t}\parallel \le \u03f5$. Therefore (A.1) always holds and we have proved the lemma. $\mathrm{\square}$
B Proof of Lemma 2
The identity follows from [2, Proposition 6.1.1] with ${\mathbf{g}}_{t}$ instead of the gradient of the loss function. It then follows from Lemma 1 that ${\beta}_{t}{\mathbf{g}}_{t}$ with ${\beta}_{t}>0$ is a descent direction at ${\overline{\mathbf{x}}}_{t+1}$. Thus, ${\mathbf{d}}_{t+1}={\mathbf{x}}_{t+1}({\beta}_{t}){\overline{\mathbf{x}}}_{t+1}$ is a feasible descent direction because ${\mathbf{d}}_{t+1}\in \mathcal{X}$ and $$ for all $t$ and ${\overline{\mathbf{x}}}_{t+1}\in \mathcal{X}$. $\mathrm{\square}$
C Proof of Lemma 3
By Assumption 5, ${f}_{t+1}$ has an $L$Lipschitz gradient. We use the following inequality from [25, Theorem 2.1.5]
${f}_{t+1}(\mathbf{y})$  $\le {f}_{t+1}(\mathbf{x})+\nabla {f}_{t+1}{(\mathbf{x})}^{\top}(\mathbf{y}\mathbf{x})$  (C.1)  
$+{\displaystyle \frac{L}{2}}{\parallel \mathbf{x}\mathbf{y}\parallel}^{2}$ 
for all $\mathbf{x},\mathbf{y}\in \mathcal{X}$. We substitute $\mathbf{y}={\mathbf{x}}_{t+1}(\beta )$ and $\mathbf{x}={\overline{\mathbf{x}}}_{t+1}$ into (C.1) to obtain
${f}_{t+1}({\mathbf{x}}_{t+1}(\beta ))$  $\le {f}_{t+1}({\overline{\mathbf{x}}}_{t+1})+{\displaystyle \frac{L}{2}}{\parallel {\mathbf{x}}_{t+1}(\beta ){\overline{\mathbf{x}}}_{t+1}\parallel}^{2}$  
$\mathrm{\hspace{1em}}+\nabla {f}_{t+1}{({\overline{\mathbf{x}}}_{t+1})}^{\top}({\mathbf{x}}_{t+1}(\beta ){\overline{\mathbf{x}}}_{t+1}).$ 
For the reminder of the proof, we use ${\mathbf{d}}_{t+1}={\mathbf{x}}_{t+1}(\beta ){\overline{\mathbf{x}}}_{t+1}$ to simplify the notation. We rewrite the gradient in term of the estimated gradient, which yields
${f}_{t+1}({\mathbf{x}}_{t+1}(\beta ))$  $\le {f}_{t+1}({\overline{\mathbf{x}}}_{t+1})+{\mathbf{g}}_{t}^{\top}{\mathbf{d}}_{t+1}{\mathbf{e}}_{t}^{\top}{\mathbf{d}}_{t+1}$  
$\mathrm{\hspace{1em}}+{\displaystyle \frac{L}{2}}{\parallel {\mathbf{d}}_{t+1}\parallel}^{2}.$  (C.2) 
By assumption, ${\mathbf{x}}_{t+1}(\beta )\ne {\overline{\mathbf{x}}}_{t+1}$, which ensures that Lemma 2 holds. We use Lemma 2 to upper bound the second term of the righthand side of (C.2). We then have: ${f}_{t+1}({\mathbf{x}}_{t+1}(\beta ))\le {f}_{t+1}({\overline{\mathbf{x}}}_{t+1})\left(\frac{1}{\beta}\frac{L}{2}\right){\parallel {\mathbf{d}}_{t+1}\parallel}^{2}+\u03f5\parallel {\mathbf{d}}_{t+1}\parallel $. Therefore, the predictive update with fixed step size will improve on the OCO update by a minimum of $\delta >0$ if the following condition is satisfied:
$$\frac{1}{\beta}{\parallel {\mathbf{d}}_{t+1}\parallel}^{2}\frac{L}{2}{\parallel {\mathbf{d}}_{t+1}\parallel}^{2}\u03f5\parallel {\mathbf{d}}_{t+1}\parallel \ge \delta .$$  (C.3) 
Assuming $$, then $\frac{1}{\beta}\ge L$, and if $\frac{L}{2}{\parallel {\mathbf{d}}_{t+1}\parallel}^{2}\u03f5\parallel {\mathbf{d}}_{t+1}\parallel \ge \delta $, then (C.3) also holds for any $\beta \in ]0,\frac{1}{L}]$. Solving for the norm of the feasible descent direction $\parallel {\mathbf{d}}_{t+1}\parallel $, we have
$$\parallel {\mathbf{d}}_{t+1}\parallel =\parallel {\mathbf{x}}_{t+1}(\beta ){\overline{\mathbf{x}}}_{t+1}\parallel \ge \frac{\u03f5}{L}+\sqrt{\frac{{\u03f5}^{2}}{{L}^{2}}+\frac{2\delta}{L}}.$$  (C.4) 
Thus, by setting $$ and satisfying (C.4), we obtain ${f}_{t+1}({\mathbf{x}}_{t+1}(\beta ))\le {f}_{t+1}({\overline{\mathbf{x}}}_{t+1})\delta ,$ where $\delta >0$. This implies that the predictive update strictly improves over the OCO update when the feasible descent direction satisfies the condition (C.4). The improvement is bounded below by $\delta $. $\mathrm{\square}$
D Proof of Theorem 1
Let ${\widehat{\mathbf{x}}}_{t}$ denote the decision variable with ${\beta}_{t}=0$ for all $t$. In other words, ${\widehat{\mathbf{x}}}_{t}$ represents the decision variable computed without the predictive algorithm. Denote the set of assumptions of Lemma 3 at round $t$ by ${\mathcal{A}}_{t}$. Let ${\mathbb{I}}_{{\mathcal{A}}_{t}}$ be the indicator function where ${\mathbb{I}}_{{\mathcal{A}}_{t}}=1$ if the assumptions are satisfied and $0$ otherwise. Observe that the improvement, ${i}_{t}$, is given by
$${i}_{t}{\mathbb{I}}_{{\mathcal{A}}_{t}}={f}_{t}({\widehat{\mathbf{x}}}_{t}){f}_{t}({\mathbf{x}}_{t}),$$  (D.1) 
where ${i}_{t}$ is the improvement when ${\mathbb{I}}_{{\mathcal{A}}_{t}}=1$. The regret of the POCO algorithm is
$${\text{\U0001d681\U0001d68e\U0001d690}}_{T}^{d}(POCO)=\sum _{t=1}^{T}{f}_{t}({\mathbf{x}}_{t}){f}_{t}({\mathbf{x}}_{t}^{\ast}).$$  (D.2) 
Using (D.1), we reexpress ${f}_{t}({\mathbf{x}}_{t})$ in (D.2):
${\text{\U0001d681\U0001d68e\U0001d690}}_{T}^{d}(POCO)$  $={\displaystyle \sum _{t=1}^{T}}{f}_{t}({\widehat{\mathbf{x}}}_{t}){f}_{t}({\mathbf{x}}_{t}^{\ast}){\displaystyle \sum _{t=1}^{T}}{i}_{t}{\mathbb{I}}_{{\mathcal{A}}_{t}}$  
$={\text{\U0001d681\U0001d68e\U0001d690}}_{T}^{d}(OCO){\displaystyle \sum _{t=1}^{T}}{i}_{t}{\mathbb{I}}_{{\mathcal{A}}_{t}}$  (D.3) 
By Lemma 3, the improvement ${i}_{t}$ is bounded below by $\delta $. We rewrite (D.3) as ${\text{\U0001d681\U0001d68e\U0001d690}}_{T}^{d}(POCO)\le {\text{\U0001d681\U0001d68e\U0001d690}}_{T}^{d}(OCO){\sum}_{t=1}^{T}\delta {\mathbb{I}}_{{\mathcal{A}}_{t}}.$ A minimum of $T\nu $ rounds satisfy ${\mathcal{A}}_{t}$ and hence ${\text{\U0001d681\U0001d68e\U0001d690}}_{T}^{d}(POCO)\le {\text{\U0001d681\U0001d68e\U0001d690}}_{T}^{d}(OCO)T\nu \delta $. $\mathrm{\square}$
E Proof of Lemma 4
We show that for some step size ${\beta}^{m}$, the estimated gradient descent projection leads to a sufficient decrease thus outperforming the OCO update. We show that if ${\beta}^{m}$ satisfies (7), then it also satisfies satisfies (6), ensuring a sufficient decrease in the objective function. Note that (7) is the condition under which the backtracking algorithm, Algorithm 2, is used. We can see from the lefthand side of the condition (7) that the update improves over the OCO update because the three last terms are bounded above by $0$, i.e., $$ by Lemma 2. Thus all three terms are less or equal to zero. By assumption, $\beta >0$, and these terms are also bounded away from zero since ${\mathbf{x}}_{t+1}\ne {\overline{\mathbf{x}}}_{t+1}$.
We start from (7) and shows it implies (6). By assumption, $\parallel {\mathbf{e}}_{t}\parallel \le \u03f5$ for all $t$ and hence (7) implies,
${f}_{t}\left({\overline{\mathbf{x}}}_{t+1}+{\beta}^{m}{\mathbf{d}}_{t+1}\right)$  $\le {f}_{t}({\overline{\mathbf{x}}}_{t+1})+{\beta}^{m}{\left({\mathbf{g}}_{t}{\mathbf{e}}_{t}\right)}^{\top}{\mathbf{d}}_{t+1}$  
$\mathrm{\hspace{1em}}2\mathrm{\Delta}.$ 
Rearranging the terms, we have
${f}_{t}\left({\overline{\mathbf{x}}}_{t+1}+{\beta}^{m}{\mathbf{d}}_{t+1}\right)+\mathrm{\Delta}$  $\le {f}_{t}({\overline{\mathbf{x}}}_{t+1})\mathrm{\Delta}$  (E.1)  
$\mathrm{\hspace{1em}}+{\beta}^{m}\nabla {f}_{t+1}{({\overline{\mathbf{x}}}_{t+1})}^{\top}{\mathbf{d}}_{t+1}.$ 
By assumption, ${f}_{t+1}$ is timeLipschitz with constant $\mathrm{\Delta}>0$ for all $\mathbf{x}\in \mathcal{X}$ and all $t$. We can therefore bound below and above respectively the lefthand and righthand side of (E.1). This leads to
${f}_{t+1}\left({\overline{\mathbf{x}}}_{t+1}+{\beta}^{m}{\mathbf{d}}_{t+1}\right)$  $\le {f}_{t+1}({\overline{\mathbf{x}}}_{t+1})$  
$\mathrm{\hspace{1em}}+{\beta}^{m}\nabla {f}_{t+1}{({\overline{\mathbf{x}}}_{t+1})}^{\top}{\mathbf{d}}_{t+1},$ 
the modified Armijo condition (6). $\mathrm{\square}$
F Proof of Theorem 2
Assume ${f}_{t}({\overline{\mathbf{x}}}_{t+1}){f}_{t}(\mathbf{x})>2\mathrm{\Delta}$. This assumption implies that ${f}_{t+1}({\overline{\mathbf{x}}}_{t+1}){f}_{t+1}(\mathbf{x})>0$ by Proposition 1. Thus, ${\overline{\mathbf{x}}}_{t+1}$ is not the minimum point of ${f}_{t+1}$. It follows that $\nabla {f}_{t+1}({\overline{\mathbf{x}}}_{t+1})\ne 0$. By assumption, ${\mathbf{d}}_{t+1}\ne \mathrm{\U0001d7ce}$ is a feasible descent direction and we have
$$  (F.1) 
Let $a\in (0,1)$. Subtracting $\nabla {f}_{t}{\left({\overline{\mathbf{x}}}_{t+1}+a\beta {\mathbf{d}}_{t+1}\right)}^{\top}{\mathbf{d}}_{t+1}$ on both side of (F.1) we obtain,
$(\nabla {f}_{t+1}({\overline{\mathbf{x}}}_{t+1})\nabla $  $$  (F.2)  
$\nabla {f}_{t}{\left({\overline{\mathbf{x}}}_{t+1}+a\beta {\mathbf{d}}_{t+1}\right)}^{\top}{\mathbf{d}}_{t+1}.$ 
If the following condition holds, then (F.2) also holds:
$\parallel \nabla {f}_{t+1}({\overline{\mathbf{x}}}_{t+1})$  $$  (F.3)  
$\nabla {f}_{t}{\left({\overline{\mathbf{x}}}_{t+1}+a\beta {\mathbf{d}}_{t+1}\right)}^{\top}{\mathbf{d}}_{t+1}.$ 
Under Assumption 3, for all $\mathbf{x},\mathbf{z}\in \mathcal{X}$ we have
$\parallel \nabla {f}_{t+1}(\mathbf{x})\nabla {f}_{t}(\mathbf{z})\parallel $  $\le \parallel \nabla {f}_{t+1}(\mathbf{x})\parallel +\parallel \nabla {f}_{t}(\mathbf{z})\parallel \le 2G$ 
and by Assumption 1, we have $\parallel {\mathbf{d}}_{t+1}\parallel \le D$. Then, if
$$  (F.4) 
holds, so does (F.3). We rewrite (F.4) as
$$  (F.5) 
Recalling Taylor’s Theorem [33, Theorem 2.1]:
$${f}_{t}(\mathbf{y}+\mathbf{p})={f}_{t}(\mathbf{y})+\nabla {f}_{t}{(\mathbf{y}+a\mathbf{p})}^{\top}\mathbf{p},$$ 
where $\mathbf{y},\mathbf{x}\in \mathcal{X}$, $\mathbf{p}\in {\mathbb{R}}^{n}$ and for some $a\in (0,1)$. We let $\mathbf{y}={\overline{\mathbf{x}}}_{t+1}$ and $\mathbf{p}={\beta}^{m}{\mathbf{d}}_{t+1}$. We have,
${f}_{t}({\overline{\mathbf{x}}}_{t+1}+\beta {\mathbf{d}}_{t+1})$  $={f}_{t}({\overline{\mathbf{x}}}_{t+1})$  (F.6)  
$\mathrm{\hspace{0.17em}}+\beta \nabla {f}_{t}{({\overline{\mathbf{x}}}_{t+1}+a\beta {\mathbf{d}}_{t+1})}^{\top}{\mathbf{d}}_{t+1}.$ 
We bound above the the last term of (F.6) using (F.5) and obtain
$$  (F.7) 
By setting $\beta \le \frac{\mathrm{\Delta}}{GD}$ in (F.7), we have $$. This shows that there always exists at least one point which satisfies the assumption on the existence of $\mathbf{x}\in \mathcal{X}$ such that ${f}_{t}({\overline{\mathbf{x}}}_{t+1}){f}_{t}(\mathbf{x})>2\mathrm{\Delta}$ that is along the feasible descent direction ${\mathbf{d}}_{t+1}$ from ${\overline{\mathbf{x}}}_{t+1}$.
Next, adapting the proof of [33, Lemma 3.1] for the modified Armijo condition for online settings (7), it follows that there exists $\overline{\beta}\le \frac{\mathrm{\Delta}}{GD}$ such that
${f}_{t}\left({\overline{\mathbf{x}}}_{t+1}+\overline{\beta}{\mathbf{d}}_{t+1}\right)$  $\le {f}_{t}({\overline{\mathbf{x}}}_{t+1})+\overline{\beta}{\mathbf{g}}_{t}^{\top}{\mathbf{d}}_{t+1}\overline{\beta}\u03f5\parallel {\mathbf{d}}_{t+1}\parallel $  
$\mathrm{\hspace{1em}}2\mathrm{\Delta}.$ 
The set $\mathcal{S}$ is therefore nonempty if there exists $\mathbf{x}\in \mathcal{X}$ such that ${f}_{t}({\overline{\mathbf{x}}}_{t+1}){f}_{t}(\mathbf{x})>2\mathrm{\Delta}$.
We now show the converse. Assuming $\mathcal{S}\ne \mathrm{\varnothing}$, then there exists $\underset{\xaf}{\beta}\in \mathcal{S}$ and
${f}_{t}\left({\overline{\mathbf{x}}}_{t+1}+\underset{\xaf}{\beta}{\mathbf{d}}_{t+1}\right)$  $$  (F.8) 
holds since $$ by Lemma 2 and $\u03f5>0$. Thus, (F.8) implies that there exists $\mathbf{x}\in \mathcal{X}$ such that ${f}_{t}({\overline{\mathbf{x}}}_{t+1}){f}_{t}(\mathbf{x})>2\mathrm{\Delta}$ and one of such point is $\mathbf{x}={\overline{\mathbf{x}}}_{t+1}+{\beta}_{2}{\mathbf{d}}_{t+1}$. This completes the proof. $\mathrm{\square}$
G Proof of Corollary 2
Since $\beta >0$, then $\mathcal{S}\ne \mathrm{\varnothing}$. By the converse of Theorem 2, we have ${f}_{t}({\overline{\mathbf{x}}}_{t+1}){f}_{t}(\mathbf{x})>2\mathrm{\Delta},$ where $\mathbf{x}={\overline{\mathbf{x}}}_{t+1}+\beta {\mathbf{d}}_{t+1}$, the decision played by the predictive update (5). The predictive update hence improves on the OCO update by at least $2\mathrm{\Delta}$. $\mathrm{\square}$
H Proof of Theorem 3
Let ${\mathbb{I}}_{{\mathcal{A}}_{t}^{\prime}}$ be the indicator function where ${\mathbb{I}}_{{\mathcal{A}}_{t}^{\prime}}=1$ if at round $t$, $\beta >0$ and $\parallel {\mathbf{g}}_{t}\parallel >\u03f5$ or $0$ otherwise. Using the same approach as in Theorem 1’s proof with Corollary 2, we obtain the regret bound. The last term of (8) is strictly positive and thus the POCOb regret is always bounded above by the OCO algorithm regret. $\mathrm{\square}$
References
 [1] D. P. Bertsekas. Nonlinear programming. Journal of the Operational Research Society, 48(3):334–334, 1997.
 [2] D. P. Bertsekas. Convex optimization algorithms. Athena Scientific Belmont, 2015.
 [3] F. Borrelli, A. Bemporad, and M. Morari. Predictive control for linear and hybrid systems. Cambridge University Press, 2017.
 [4] D. S Callaway. Tapping the energy storage potential in electric loads to deliver load following and regulation, with application to wind energy. Energy Conversion and Management, 50(5):1389–1400, 2009.
 [5] D.S. Callaway and I. A. Hiskens. Achieving controllability of electric loads. Proceedings of the IEEE, 99(1):184–199, 2011.
 [6] X. Cao, J. Zhang, and H. V. Poor. A virtualqueue based algorithm for constrained online convex optimization with applications to data center resource allocation. IEEE Journal of Selected Topics in Signal Processing, 2018.
 [7] T. Chen and G. B. Giannakis. Bandit convex optimization for scalable and dynamic iot management. IEEE Internet of Things Journal, 2018.
 [8] T. Chen, Q. Ling, and G. B. Giannakis. An online convex optimization approach to proactive network resource allocation. IEEE Transactions on Signal Processing, 65(24):6350–6364, 2017.
 [9] C.K. Chiang, T. Yang, C.J. Lee, M. Mahdavi, C.J. Lu, R. Jin, and S. Zhu. Online optimization with gradual variations. In Conference on Learning Theory, pages 6–1, 2012.
 [10] O. Dekel, A. Flajolet, N. Haghtalab, and P. Jaillet. Online learning with a hint. In Advances in Neural Information Processing Systems, pages 5305–5314, 2017.
 [11] S. Diamond and S. Boyd. CVXPY: A Pythonembedded modeling language for convex optimization. Journal of Machine Learning Research, 17(83):1–5, 2016.
 [12] A. Domahidi, E. Chu, and S. Boyd. ECOS: An SOCP solver for embedded systems. In European Control Conference (ECC), pages 3071–3076, 2013.
 [13] J. C. Duchi, S. ShalevShwartz, Y. Singer, and A. Tewari. Composite objective mirror descent. In COLT, pages 14–26, 2010.
 [14] C. E. Garcia, D. M. Prett, and M. Morari. Model predictive control: theory and practice—a survey. Automatica, 25(3):335–348, 1989.
 [15] E. Hazan. Introduction to online convex optimization. Foundations and Trends® in Optimization, 2(34):157–325, 2016.
 [16] E. Hazan, A. Agarwal, and S. Kale. Logarithmic regret algorithms for online convex optimization. Machine Learning, 69(23):169–192, 2007.
 [17] E. Hazan and N. Megiddo. Online learning with prior knowledge. In International Conference on Computational Learning Theory, pages 499–513. Springer, 2007.
 [18] N. HoNguyen and F. KılınçKarzan. Accelerating optimization under uncertainty via online convex optimization. Technical report, August, 2016.
 [19] A. Jadbabaie, A. Rakhlin, S. Shahrampour, and K. Sridharan. Online optimization: Competing with dynamic comparators. In Artificial Intelligence and Statistics, pages 398–406, 2015.
 [20] S.J. Kim and G. B. Giannakis. An online convex optimization approach to realtime energy pricing for demand response. IEEE Transactions on Smart Grid, 8(6):2784–2793, 2017.
 [21] A.e LesageLandry and J.A. Taylor. Setpoint tracking with partially observed loads. IEEE Transactions on Power Systems, 33(5):5615 – 5627, 2018.
 [22] J. L. Mathieu, S. Koch, and D. S. Callaway. State estimation and control of electric loads to manage realtime energy imbalance. IEEE Transactions on Power Systems, 28(1):430–440, 2013.
 [23] M. Mohri and S. Yang. Accelerating online convex optimization via adaptive prediction. In Artificial Intelligence and Statistics, pages 848–856, 2016.
 [24] A. Mokhtari, S. Shahrampour, A. Jadbabaie, and A. Ribeiro. Online optimization in dynamic environments: Improved regret rates for strongly convex problems. In Decision and Control (CDC), 2016 IEEE 55th Conference on, pages 7195–7201. IEEE, 2016.
 [25] Y. Nesterov. Introductory lectures on convex programming volume i: Basic course. Lecture notes, 1998.
 [26] P. Palensky and D. Dietrich. Demand side management: Demand response, intelligent energy systems, and smart loads. IEEE transactions on industrial informatics, 7(3):381–388, 2011.
 [27] A. Rakhlin and K. Sridharan. Online learning with predictable sequences. In Proceedings of the 26th Annual Conference on Learning Theory (COLT), pages 1–27, 2013.
 [28] S. Rakhlin and K. Sridharan. Optimization, learning, and games with predictable sequences. In Advances in Neural Information Processing Systems, pages 3066–3074, 2013.
 [29] M. Schmidt, N. L. Roux, and F. R. Bach. Convergence rates of inexact proximalgradient methods for convex optimization. In Advances in neural information processing systems, pages 1458–1466, 2011.
 [30] S. ShalevShwartz. Online learning and online convex optimization. Foundations and Trends® in Machine Learning, 4(2):107–194, 2012.
 [31] J. Steinhardt and P. Liang. Adaptivity and optimism: An improved exponentiated gradient algorithm. In International Conference on Machine Learning, pages 1593–1601, 2014.
 [32] J. A. Taylor, S. V. Dhople, and D. S. Callaway. Power systems without fuel. Renewable and Sustainable Energy Reviews, 57:1322–1336, 2016.
 [33] S. Wright and J. Nocedal. Numerical optimization. Springer Science, 35(6768), 1999.
 [34] S. Yang and M. Mohri. Optimistic bandit convex optimization. In Advances in Neural Information Processing Systems, pages 2297–2305, 2016.
 [35] M. Zinkevich. Online convex programming and generalized infinitesimal gradient ascent. In Proceedings of the 20th International Conference on Machine Learning (ICML03), pages 928–936, 2003.