Distributionally Robust Covariance Steering
with Optimal Risk Allocation
Venkatraman Renganathan, Joshua Pilipovsky, Panagiotis Tsiotras
Abstract—This article extends the optimal covariance steering
(CS) problem for discrete time linear stochastic systems modeled
using moment-based ambiguity sets. To hedge against the un-
certainty in the state distributions while performing covariance
steering, distributionally robust risk constraints are employed
during the optimal allocation of the risk. Specifically, a distri-
butionally robust iterative risk allocation (DR-IRA) formalism
is used to solve the optimal risk allocation problem for the CS
problem using a two-stage approach. The upper-stage of DR-IRA
is a convex problem that optimizes the risk, while the lower-
stage optimizes the controller with the new distributionally robust
risk constraints. The proposed framework results in solutions
that are robust against arbitrary distributions in the considered
ambiguity set. Finally, we demonstrate our proposed approach
using numerical simulations. Addressing the covariance steering
problem through the lens of distributional robustness marks the
novel contribution of this article.
Index Terms—Covariance Steering, Distributional Robustness,
Iterative Risk Allocation, Stochastic Systems, Moment Based
Ambiguity Sets.
I. INTRODUCTION
Intelligent and adaptive systems of the “smart world” that
work under operational constraints seek to solve some instance
of a constrained optimal control problem for optimizing their
performance. Such constrained optimal control problems can
now be increasingly solved efficiently using several numerical
optimization techniques. For instance, robot path planning in
uncertain environments [1]–[6] has gained the attention of
researchers worldwide as robots are being increasingly de-
ployed to solve many real-world problems. Apart from realistic
constraints, reliability of operation of these systems is often
thwarted by the ineffective handling of system uncertainties,
which can be either deterministic or stochastic.
Control of stochastic systems can be best formulated as a
problem of controlling the distribution of trajectories over time
subject to constraints. Recently, the finite horizon covariance
steering (CS) problem, namely, the problem of steering an
initial distribution to a final distribution at a specific final
time step subject to linear time varying dynamics has been
explored [7]. Specifically, the control problem in the CS
problem setting involves steering the mean and the covariance
to the desired terminal values. When the decision-making
process relies blindly on the functional form of the process
that models the stochastic uncertainty, it is known to result
in potentially severe miscalculation of risk. For instance,
Gaussianity assumptions made in the name of tractability in
V. Renganathan is with the Department of Automatic Control LTH,
Lund University, Lund, Sweden. J. Pilipovsky and P. Tsoitras are with the
School of Aerospace Engineering, Georgia Institute of Technology, Atlanta,
GA 30332-0150 USA. E-mail: venkatraman.renganathan@control.lth.se,
jpilipovsky3@gatech.edu, tsiotras@gatech.edu.
several modeling regimes are actually rarely justifiable, as
the true distribution that governs the uncertain data might
be non-Gaussian. Such shortcomings can be mitigated with
risk-based stochastic optimization where the risk of wrong
decisions can be appropriately handled to result in risk-averse
decision making. One such tool is the Distributionally Robust
Optimization (DRO) advocated in [8], [9] which enables
modelers to explicitly incorporate ambiguity in probability
distributions into the optimization problem.
Control of stochastic systems often involves optimizing
the system’s objective subject to chance constraints, where
one assumes that the system uncertainties follow a known
distribution and enforces that the system constraints hold with
high probability as a function of the decision variables. The
number of constraint violations, called the total risk budget,
is usually a user-defined a priori specification and is a natural
metric to assess risk. Hence, one may consider the problem
of finding a risk allocation procedure that will allocate the
probability of violating each individual chance constraint at
each time step. Given a number of chance constraints across a
finite horizon, the total risk budget has to be allocated for all
chance constraints across all time steps [10]. It is a common
practice to consider a uniform risk allocation, i.e., allocate the
same risk for all constraints and across all time steps. However,
risk allocation can be optimized as in [11], [12] to reduce the
conservatism resulting from a uniform risk allocation. If the
probability distribution of the system uncertainties are known
exactly, then non-uniform risk allocation can be performed
effectively. However, risk allocation optimization with arbi-
trary distribution of the system uncertainties has not yet been
explored till now. This article addresses this shortcoming.
Contributions: Since authors in [13] solved the CS problem
for the Gaussian case, this article extends it with arbitrary
distributions using the theory of distributional robustness (DR).
To the best of our knowledge, this article is the first one to ex-
tend the CS problem using distributionally robust optimization
techniques for both polytopic and convex conic state constraint
sets. Our main contributions in this article are as follows:
1) We extend the covariance steering problem tailored
between Gaussian distributions to arbitrary distributions
modeled using moment based ambiguity sets.
2) We enforce distributionally robust risk constraints for
both polytopic and convex cone state constraint satisfac-
tion, while solving the covariance steering problem, and
obtain the optimal risk allocation through a distribution-
ally robust iterative risk allocation (DR-IRA) algorithm.
3) We demonstrate our approach using simulation examples
and show the effectiveness of the proposed generaliza-
tion for covariance steering problems between arbitrary
distributions in moment-based ambiguity sets.
arXiv:2210.00050v2 [math.OC] 7 Nov 2022