Low-complexity modeling of partially available second-order statistics: theory and an efficient matrix completion algorithm

Abstract

State statistics of linear systems satisfy certain structural constraints that arise from the underlying dynamics and the directionality of input disturbances. In the present paper we study the problem of completing partially known state statistics. Our aim is to develop tools that can be used in the context of control-oriented modeling of large-scale dynamical systems. For the type of applications we have in mind, the dynamical interaction between state variables is known while the directionality and dynamics of input excitation is often uncertain. Thus, the goal of the mathematical problem that we formulate is to identify the dynamics and directionality of input excitation in order to explain and complete observed sample statistics. More specifically, we seek to explain correlation data with the least number of possible input disturbance channels. We formulate this inverse problem as rank minimization, and for its solution, we employ a convex relaxation based on the nuclear norm. The resulting optimization problem is cast as a semidefinite program and can be solved using general-purpose solvers. For problem sizes that these solvers cannot handle, we develop a customized alternating minimization algorithm (AMA). We interpret AMA as a proximal gradient for the dual problem and prove sub-linear convergence for the algorithm with fixed step-size. We conclude with an example that illustrates the utility of our modeling and optimization framework and draw contrast between AMA and the commonly used alternating direction method of multipliers (ADMM) algorithm.

Submitted to IEEE Transactions on Automatic Control

References in corpus (5)

• Color of turbulence
• Model-based design of transverse wall oscillations for turbulent drag reduction
• G-AMA: Sparse Gaussian graphical model estimation via alternating minimization
• Fast cooling for a system of stochastic oscillators
• Optimal steering of a linear stochastic system to a final probability distribution

Cited by in corpus (12)

• From bypass transition to flow control and data-driven turbulence modeling: An input-output viewpoint
• Resolvent-based estimation of space-time flow statistics
• Model-based design of riblets for turbulent drag reduction
• Stochastic dynamical modeling of turbulent flows
• Proximal algorithms for large-scale statistical modeling and sensor/actuator selection
• Stochastic receptivity analysis of boundary layer flow
• Low-rank Optimization with Convex Constraints
• Model-based spectral coherence analysis
• Modeling mode interactions in boundary layer flows via the Parabolized Floquet Equations
• Low-Rank Inducing Norms with Optimality Interpretations
• Input-output analysis of stochastic base flow uncertainty
• Efficient Proximal Mapping Computation for Unitarily Invariant Low-Rank Inducing Norms