Source author record

Bassam Bamieh

Bassam Bamieh appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2024arXiv

Linear-Quadratic Problems in Systems and Controls via Covariance Representations and Linear-Conic Duality: Finite-Horizon Case

Linear-Quadratic (LQ) problems that arise in systems and controls include the classical optimal control problems of the Linear Quadratic Regulator (LQR) in both its deterministic and stochastic forms, as well as $H^\infty$-analysis (the Bounded Real Lemma), the Positive Real Lemma, and general Integral Quadratic Constraints (IQCs) tests. We present a unified treatment of all of these problems using an approach which converts linear-quadratic problems to matrix-valued linear-linear problems with a positivity constraint. This is done through a system representation where the joint state/input covariance (the outer product in the deterministic case) matrix is the fundamental object. LQ problems then become infinite-dimensional semidefinite programs, and the key tool used is that of linear-conic duality. Linear Matrix Inequalities (LMIs) emerge naturally as conal constraints on dual problems. Riccati equations characterize extrema of these special LMIs, and therefore provide solutions to the dual problems. The state-feedback structure of all optimal signals in these problems emerge out of alignment (complementary slackness) conditions between primal and dual problems. Perhaps the new insight gained from this approach is that first LMIs, and then second, Riccati equations arise naturally in dual, rather than primal problems. Furthermore, while traditional LQ problems are set up in $L^2$ spaces of signals, their equivalent covariance-representation problems are most naturally set up in $L^1$ spaces of matrix-valued signals.

preprint2022arXiv

A Short Introduction to the Koopman Representation of Dynamical Systems

The Koopman representation is an infinite dimensional linear representation of linear or nonlinear dynamical systems. It represents the dynamics of output maps (aka observables), which are functions on the state space whose evaluation is interpreted as an output. Conceptually simple derivations and commentary on the Koopman representation are given. We emphasize an important duality between initial conditions and output maps of the original system, and those of the Koopman representation. This duality is an important consideration when this representation is used in data-driven applications such as the Dynamic Mode Decomposition (DMD) and its variants. The adjoint relation between the Koopman representation and the transfer operator of mass transport is also shown.

preprint2022arXiv

A Tutorial on Matrix Perturbation Theory (using compact matrix notation)

Analytic perturbation theory for matrices and operators is an immensely useful mathematical technique. Most elementary introductions to this method have their background in the physics literature, and quantum mechanics in particular. In this note, we give an introduction to this method that is independent of any physics notions, and relies purely on concepts from linear algebra. An additional feature of this presentation is that matrix notation and methods are used throughout. In particular, we formulate the equations for each term of the analytic expansions of eigenvalues and eigenvectors as {\em matrix equations}, namely Sylvester equations in particular. Solvability conditions and explicit expressions for solutions of such matrix equations are given, and expressions for each term in the analytic expansions are given in terms of those solutions. This unified treatment simplifies somewhat the complex notation that is commonly seen in the literature, and in particular, provides relatively compact expressions for the non-Hermitian and degenerate cases, as well as for higher order terms.

preprint2022arXiv

A Tutorial on Solution Properties of State Space Models of Dynamical Systems

The starting point of analysis of state space models is investigating existence, uniqueness and solution properties such as the semigroup property, and various formulas for the solutions. Several concepts such as the state transition matrix, the matrix exponential, the variations of constants formula (the Cauchy formula), the Peano-Baker series, and the Picard iteration are used to characterize solutions. In this note, a tutorial treatment is given where all of these concepts are shown to be various manifestations of a single abstract method, namely solving equations using an operator Neumann series involving the Volterra operator of forward integration. The matrix exponential, the Peano-Baker series, the Picard iteration, and the Cauchy formula can be "discovered" naturally from this Neumann series. The convergence of the series and iterations is a consequence of the key property of asymptotic nilpotence of the Volterra operator. This property is an asymptotic version of the nilpotence property of a strictly-lower-triangular matrix.

preprint2022arXiv

An Explicit Parametrization of Closed Loops for Spatially Distributed Controllers with Sparsity Constraints

In this article, we study the linear time-invariant state-feedback controller design problem for distributed systems. We follow the recently developed system level synthesis (SLS) approach and impose locality structure on the resulting closed-loop mappings; the corresponding controller implementation inherits this prescribed structure. In contrast to existing SLS results, we derive an explicit (rather than implicit) parameterization of all achievable stabilized closed-loops. This admits more efficient IIR representations of the temporal part of the closed-loop dynamics, and it allows for the H2 design problem with closed-loop spatial sparsity constraints to be converted to a standard model matching problem, with the number of transfer function parameters scaling linearly with the closed-loop spatial extent constraint. We illustrate our results with two applications: consensus of first-order subsystems and the vehicular platoons problem. In the case of first-order consensus, we provide analytic solutions and further analyze the architecture of the resulting controller implementation. Results for infinite extent spatially invariant systems are presented to provide insight to the case of a large but finite number of subsystems.

preprint2022arXiv

Discovering Transforms: A Tutorial on Circulant Matrices, Circular Convolution, and the Discrete Fourier Transform

How could the Fourier and other transforms be naturally discovered if one didn't know how to postulate them? In the case of the Discrete Fourier Transform (DFT), we show how it arises naturally out of analysis of circulant matrices. In particular, the DFT can be derived as the change of basis that simultaneously diagonalizes all circulant matrices. In this way, the DFT arises naturally from a linear algebra question about a set of matrices. Rather than thinking of the DFT as a signal transform, it is more natural to think of it as a single change of basis that renders an entire set of mutually-commuting matrices into simple, diagonal forms. The DFT can then be "discovered" by solving the eigenvalue/eigenvector problem for a special element in that set. A brief outline is given of how this line of thinking can be generalized to families of linear operators, leading to the discovery of the other common Fourier-type transforms, as well as its connections with group representations theory.

preprint2022arXiv

On Structured-Closed-Loop versus Structured-Controller Design: the Case of Relative Measurement Feedback

We consider the optimal distributed controller design problem subject to two structural requirements: locality, i.e. available measurements and sub-controllers' interactions are governed by a graph structure, and relative feedback, i.e. only differences of measurements are available to the controller. We formalize controller locality in terms of the controller's transfer function, state-space realization, or resulting closed-loop mapping. We demonstrate that the relative feedback requirement can be written as a convex constraint on the controller and (in special cases) on the resulting closed-loop, and we characterize the allowable structures of relative feedback controllers. We prove that sparse closed-loop design is a convex relaxation of structured controller state-space design, even in the continuous time IIR setting. This formalizes and extends results of the recently developed System Level Synthesis framework. We take a first step toward quantifying the performance gap associated with this convex relaxation by constructing a class of examples (based on relative feedback requirements) for which the difference in performance, measured by an H2 norm, is infinite. The results presented are used to contrast several issues of structural constraints in distributed control design that remain as open problems.

preprint2020arXiv

Diffusion and Consensus in a Weakly Coupled Network of Networks

We study diffusion and consensus dynamics in a Network of Networks model. In this model, there is a collection of sub-networks, connected to one another using a small number of links. We consider a setting where the links between networks have small weights, or are used less frequently than links within each sub-network. Using spectral perturbation theory, we analyze the diffusion rate and convergence rate of the investigated systems. Our analysis shows that the first order approximation of the diffusion and convergence rates is independent of the topologies of the individual graphs; the rates depend only on the number of nodes in each graph and the topology of the connecting edges. The second order analysis shows a relationship between the diffusion and convergence rates and the information centrality of the connecting nodes within each sub-network. We further highlight these theoretical results through numerical examples.

preprint2018arXiv

An Input-Output Approach to Structured Stochastic Uncertainty

We consider linear time invariant systems with exogenous stochastic disturbances, and in feedback with structured stochastic uncertainties. This setting encompasses linear systems with both additive and multiplicative noise. Our concern is to characterize second-order properties such as mean-square stability and performance. A purely input-output treatment of these systems is given without recourse to state space models, and thus the results are applicable to certain classes of distributed systems. We derive necessary and sufficient conditions for mean-square stability in terms of the spectral radius of a linear matrix operator whose dimension is that of the number of uncertainties, rather than the dimension of any underlying state space models. Our condition is applicable to the case of correlated uncertainties, and reproduces earlier results for uncorrelated uncertainties. For cases where state space realizations are given, Linear Matrix Inequality (LMI) equivalents of the input-output conditions are given.

preprint2013arXiv

Consensus and Coherence in Fractal Networks

We consider first and second order consensus algorithms in networks with stochastic disturbances. We quantify the deviation from consensus using the notion of network coherence, which can be expressed as an $H_2$ norm of the stochastic system. We use the setting of fractal networks to investigate the question of whether a purely topological measure, such as the fractal dimension, can capture the asymptotics of coherence in the large system size limit. Our analysis for first-order systems is facilitated by connections between first-order stochastic consensus and the global mean first passage time of random walks. We then show how to apply similar techniques to analyze second-order stochastic consensus systems. Our analysis reveals that two networks with the same fractal dimension can exhibit different asymptotic scalings for network coherence. Thus, this topological characterization of the network does not uniquely determine coherence behavior. The question of whether the performance of stochastic consensus algorithms in large networks can be captured by purely topological measures, such as the spatial dimension, remains open.

preprint2013arXiv

On sensor fusion for airborne wind energy systems

A study on filtering aspects of airborne wind energy generators is presented. This class of renewable energy systems aims to convert the aerodynamic forces generated by tethered wings, flying in closed paths transverse to the wind flow, into electricity. The accurate reconstruction of the wing's position, velocity and heading is of fundamental importance for the automatic control of these kinds of systems. The difficulty of the estimation problem arises from the nonlinear dynamics, wide speed range, large accelerations and fast changes of direction that the wing experiences during operation. It is shown that the overall nonlinear system has a specific structure allowing its partitioning into sub-systems, hence leading to a series of simpler filtering problems. Different sensor setups are then considered, and the related sensor fusion algorithms are presented. The results of experimental tests carried out with a small-scale prototype and wings of different sizes are discussed. The designed filtering algorithms rely purely on kinematic laws, hence they are independent from features like wing area, aerodynamic efficiency, mass, etc. Therefore, the presented results are representative also of systems with larger size and different wing design, different number of tethers and/or rigid wings.

preprint2012arXiv

The Price of Synchrony: Resistive Losses due to Phase Synchronization in Power Networks

We investigate the total resistive losses incurred in returning a power network of identical generators to a synchronous state following a transient stability event or in maintaining this state in the presence of persistent stochastic disturbances. We formulate this cost as the input-output $H^2$ norm of a linear dynamical system with distributed disturbances. We derive an expression for the total resistive losses that scales with the size of the network as well as properties of the generators and power lines, but is independent of the network topology. This topologically invariant scaling of what we term the price of synchrony is in contrast to typical power system stability notions like rate of convergence or the region of attraction for rotor-angle stability. Our result indicates that highly connected power networks, whilst desirable for higher phase synchrony, do not offer an advantage in terms of the total resistive power losses needed to achieve this synchrony. Furthermore, if power flow is the mechanism used to achieve synchrony in highly-distributed-generation networks, the cost increases unboundedly with the number of generators.

preprint2011arXiv

Coherence in Large-Scale Networks: Dimension-Dependent Limitations of Local Feedback

We consider distributed consensus and vehicular formation control problems. Specifically we address the question of whether local feedback is sufficient to maintain coherence in large-scale networks subject to stochastic disturbances. We define macroscopic performance measures which are global quantities that capture the notion of coherence; a notion of global order that quantifies how closely the formation resembles a solid object. We consider how these measures scale asymptotically with network size in the topologies of regular lattices in 1, 2 and higher dimensions, with vehicular platoons corresponding to the 1 dimensional case. A common phenomenon appears where a higher spatial dimension implies a more favorable scaling of coherence measures, with a dimensions of 3 being necessary to achieve coherence in consensus and vehicular formations under certain conditions. In particular, we show that it is impossible to have large coherent one dimensional vehicular platoons with only local feedback. We analyze these effects in terms of the underlying energetic modes of motion, showing that they take the form of large temporal and spatial scales resulting in an accordion-like motion of formations. A conclusion can be drawn that in low spatial dimensions, local feedback is unable to regulate large-scale disturbances, but it can in higher spatial dimensions. This phenomenon is distinct from, and unrelated to string instability issues which are commonly encountered in control problems for automated highways.

preprint2010arXiv

Amplification and Nonlinear Mechanisms in Plane Couette Flow

We study the input-output response of a streamwise constant projection of the Navier-Stokes equations for plane Couette flow, the so-called 2D/3C model. Study of a streamwise constant model is motivated by numerical and experimental observations that suggest the prevalence and importance of streamwise and quasi-streamwise elongated structures. Periodic spanwise/wall-normal (z-y) plane stream functions are used as input to develop a forced 2D/3C streamwise velocity field that is qualitatively similar to a fully turbulent spatial field of DNS data. The input-output response associated with the 2D/3C nonlinear coupling is used to estimate the energy optimal spanwise wavelength over a range of Reynolds numbers. The results of the input-output analysis agree with previous studies of the linearized Navier-Stokes equations. The optimal energy corresponds to minimal nonlinear coupling. On the other hand, the nature of the forced 2D/3C streamwise velocity field provides evidence that the nonlinear coupling in the 2D/3C model is responsible for creating the well known characteristic "S" shaped turbulent velocity profile. This indicates that there is an important tradeoff between energy amplification, which is primarily linear and the seemingly nonlinear momentum transfer mechanism that produces a turbulent-like mean profile.