Showing papers in "arXiv: Optimization and Control in 2007"

Cooperative Robot Control and Concurrent Synchronization of Lagrangian Systems

Abstract: Concurrent synchronization is a regime where diverse groups of fully synchronized dynamic systems stably coexist. We study global exponential synchronization and concurrent synchronization in the context of Lagrangian systems control. In a network constructed by adding diffusive couplings to robot manipulators or mobile robots, a decentralized tracking control law globally exponentially synchronizes an arbitrary number of robots, and represents a generalization of the average consensus problem. Exact nonlinear stability guarantees and synchronization conditions are derived by contraction analysis. The proposed decentralized strategy is further extended to adaptive synchronization and partial-state coupling.

Approximately bisimilar symbolic models for nonlinear control systems

Abstract: Control systems are usually modeled by differential equations describing how physical phenomena can be influenced by certain control parameters or inputs. Although these models are very powerful when dealing with physical phenomena, they are less suitable to describe software and hardware interfacing the physical world. For this reason there is a growing interest in describing control systems through symbolic models that are abstract descriptions of the continuous dynamics, where each "symbol" corresponds to an "aggregate" of states in the continuous model. Since these symbolic models are of the same nature of the models used in computer science to describe software and hardware, they provide a unified language to study problems of control in which software and hardware interact with the physical world. Furthermore the use of symbolic models enables one to leverage techniques from supervisory control and algorithms from game theory for controller synthesis purposes. In this paper we show that every incrementally globally asymptotically stable nonlinear control system is approximately equivalent (bisimilar) to a symbolic model. The approximation error is a design parameter in the construction of the symbolic model and can be rendered as small as desired. Furthermore if the state space of the control system is bounded the obtained symbolic model is finite. For digital control systems, and under the stronger assumption of incremental input-to-state stability, symbolic models can be constructed through a suitable quantization of the inputs.

Fractional conservation laws in optimal control theory

Abstract: Using the recent formulation of Noether's theorem for the problems of the calculus of variations with fractional derivatives, the Lagrange multiplier technique, and the fractional Euler-Lagrange equations, we prove a Noether-like theorem to the more general context of the fractional optimal control. As a corollary, it follows that in the fractional case the autonomous Hamiltonian does not define anymore a conservation law. Instead, it is proved that the fractional conservation law adds to the Hamiltonian a new term which depends on the fractional-order of differentiation, the generalized momentum, and the fractional derivative of the state variable.

Fractional Optimal Control in the Sense of Caputo and the Fractional Noether's Theorem

Abstract: The study of fractional variational problems with derivatives in the sense of Caputo is a recent subject, the main results being Agrawal’s necessary optimality conditions of Euler-Lagrange and respective transversality conditions. Using Agrawal’s Euler-Lagrange equation and the Lagrange multiplier technique, we obtain here a Noether-like theorem for fractional optimal control problems in the sense of Caputo.

Modeling of Transitional Channel Flow Using Balanced Proper Orthogonal Decomposition

Abstract: We study reduced-order models of three-dimensional perturbations in linearized channel flow using balanced proper orthogonal decomposition (BPOD). The models are obtained from three-dimensional simulations in physical space as opposed to the traditional single-wavenumber approach, and are therefore better able to capture the effects of localized disturbances or localized actuators. In order to assess the performance of the models, we consider the impulse response and frequency response, and variation of the Reynolds number as a model parameter. We show that the BPOD procedure yields models that capture the transient growth well at a low order, whereas standard POD does not capture the growth unless a considerably larger number of modes is included, and even then can be inaccurate. In the case of a localized actuator, we show that POD modes which are not energetically significant can be very important for capturing the energy growth. In addition, a comparison of the subspaces resulting from the two methods suggests that the use of a non-orthogonal projection with adjoint modes is most likely the main reason for the superior performance of BPOD. We also demonstrate that for single-wavenumber perturbations, low-order BPOD models reproduce the dominant eigenvalues of the full system better than POD models of the same order. These features indicate that the simple, yet accurate BPOD models are a good candidate for developing model-based controllers for channel flow.

Semidefinite Representation of Convex Sets

Abstract: Let $S =\{x\in \re^n: g_1(x)\geq 0, ..., g_m(x)\geq 0\}$ be a semialgebraic set defined by multivariate polynomials $g_i(x)$. Assume $S$ is convex, compact and has nonempty interior. Let $S_i =\{x\in \re^n: g_i(x)\geq 0\}$, and $\bdS$ (resp. $\bdS_i$) be the boundary of $S$ (resp. $S_i$). This paper discusses whether $S$ can be represented as the projection of some LMI representable set. Such $S$ is called semidefinite representable or SDP representable. The contributions of this paper: {\bf (i)} Assume $g_i(x)$ are all concave on $S$. If the positive definite Lagrange Hessian (PDLH) condition holds, i.e., the Hessian of the Lagrange function for optimization problem of minimizing any nonzero linear function $\ell^Tx$ on $S$ is positive definite at the minimizer, then $S$ is SDP representable. {\bf (ii)} If each $g_i(x)$ is either sos-concave ($- abla^2g_i(x)=W(x)^TW(x)$ for some possibly nonsquare matrix polynomial $W(x)$) or strictly quasi-concave on $S$, then $S$ is SDP representable. {\bf (iii)} If each $S_i$ is either sos-convex or poscurv-convex ($S_i$ is compact convex, whose boundary has positive curvature and is nonsingular, i.e. $ abla g_i(x) ot = 0$ on $\bdS_i \cap S$), then $S$ is SDP representable. This also holds for $S_i$ for which $\bdS_i \cap S$ extends smoothly to the boundary of a poscurv-convex set containing $S$. {\bf (iv)} We give the complexity of Schm\"{u}dgen and Putinar's matrix Positivstellensatz, which are critical to the proofs of (i)-(iii).

Guaranteed Minimum-Rank Solutions of Linear Matrix Equations via Nuclear Norm Minimization

Abstract: The affine rank minimization problem consists of finding a matrix of minimum rank that satisfies a given system of linear equality constraints. Such problems have appeared in the literature of a diverse set of fields including system identification and control, Euclidean embedding, and collaborative filtering. Although specific instances can often be solved with specialized algorithms, the general affine rank minimization problem is NP-hard. In this paper, we show that if a certain restricted isometry property holds for the linear transformation defining the constraints, the minimum rank solution can be recovered by solving a convex optimization problem, namely the minimization of the nuclear norm over the given affine space. We present several random ensembles of equations where the restricted isometry property holds with overwhelming probability. The techniques used in our analysis have strong parallels in the compressed sensing framework. We discuss how affine rank minimization generalizes this pre-existing concept and outline a dictionary relating concepts from cardinality minimization to those of rank minimization.

Symbolic Models for Nonlinear Control Systems: Alternating Approximate Bisimulations

Abstract: Symbolic models are abstract descriptions of continuous systems in which symbols represent aggregates of continuous states. In the last few years there has been a growing interest in the use of symbolic models as a tool for mitigating complexity in control design. In fact, symbolic models enable the use of well known algorithms in the context of supervisory control and algorithmic game theory, for controller synthesis. Since the 1990's many researchers faced the problem of identifying classes of dynamical and control systems that admit symbolic models. In this paper we make a further progress along this research line by focusing on control systems affected by disturbances. Our main contribution is to show that incrementally globally asymptotically stable nonlinear control systems with disturbances admit symbolic models. When specializing these results to linear systems, we show that these symbolic models can be easily constructed.

Non-linear Symmetry-preserving Observer on Lie Groups

Abstract: In this paper we give a geometrical framework for the design of observers on finite-dimensional Lie groups for systems which possess some specific symmetries. The design and the error (between true and estimated state) equation are explicit and intrinsic. We consider also a particular case: left-invariant systems on Lie groups with right equivariant output. The theory yields a class of observers such that error equation is autonomous. The observers converge locally around any trajectory, and the global behavior is independent from the trajectory, which reminds of the linear stationary case.

Kalman Filtering with Equality and Inequality State Constraints

Abstract: Both constrained and unconstrained optimization problems regularly appear in recursive tracking problems engineers currently address – however, constraints are rarely exploited for these applications. We define the Kalman Filter and discuss two different approaches to incorporating constraints. Each of these approaches are first applied to equality constraints and then extended to inequality constraints. We discuss methods for dealing with nonlinear constraints and for constraining the state prediction. Finally, some experiments are provided to indicate the usefulness of such methods.

A Semismooth Newton Method for Tikhonov Functionals with Sparsity Constraints

Abstract: Minimization problems in $\ell^2$ for Tikhonov functionals with sparsity constraints are considered. Sparsity of the solution is ensured by a weighted $\ell^1$ penalty term. The necessary and sufficient condition for optimality is shown to be slantly differentiable (Newton differentiable), hence a semismooth Newton method is applicable. Local superlinear convergence of this method is proved. Numerical examples are provided which show that our method compares favorably with existing approaches.

On the Nonexistence of Quadratic Lyapunov Functions for Consensus Algorithms

Abstract: We provide an example proving that there exists no quadratic Lyapunov function for a certain class of linear agreement/consensus algorithms, a fact that had been numerically verified in [5]. We also briefly discuss sufficient conditions for the existence of such a Lyapunov function.

Optimal transportation on non-compact manifolds

Abstract: In this work, we show how to obtain for non-compact manifolds the results that have already been done for Monge Transport Problem for costs coming from Tonelli Lagrangians on compact manifolds. In particular, the already known results for a cost of the type $d^r,r>1$, where $d$ is the Riemannian distance of a complete Riemannian manifold, hold without any curvature restriction.

GloptiPoly 3: moments, optimization and semidefinite programming

Abstract: We describe a major update of our Matlab freeware GloptiPoly for parsing generalized problems of moments and solving them numerically with semidefinite programming.

On Distributed Averaging Algorithms and Quantization Effects

Abstract: We consider distributed iterative algorithms for the averaging problem over time-varying topologies. Our focus is on the convergence time of such algorithms when complete (unquantized) information is available, and on the degradation of performance when only quantized information is available. We study a large and natural class of averaging algorithms, which includes the vast majority of algorithms proposed to date, and provide tight polynomial bounds on their convergence time. We also describe an algorithm within this class whose convergence time is the best among currently available averaging algorithms for time-varying topologies. We then propose and analyze distributed averaging algorithms under the additional constraint that agents can only store and communicate quantized information, so that they can only converge to the average of the initial values of the agents within some error. We establish bounds on the error and tight bounds on the convergence time, as a function of the number of quantization levels.

The Fermat-Torricelli problem in normed planes and spaces

Abstract: We investigate the Fermat-Torricelli problem in d-dimensional real normed spaces or Minkowski spaces, mainly for d=2. Our approach is to study the Fermat-Torricelli locus in a geometric way. We present many new results, as well as give an exposition of known results that are scattered in various sources, with proofs for some of them. Together, these results can be considered to be a minitheory of the Fermat-Torricelli problem in Minkowski spaces and especially in Minkowski planes. This demonstrates that substantial results about locational problems valid for all norms can be found using a geometric approach.

Singular arcs in the generalized Goddard's Problem

Abstract: We investigate variants of Goddard's problems for nonvertical trajectories. The control is the thrust force, and the objective is to maximize a certain final cost, typically, the final mass. In this report, performing an analysis based on the Pontryagin Maximum Principle, we prove that optimal trajectories may involve singular arcs (along which the norm of the thrust is neither zero nor maximal), that are computed and characterized. Numerical simulations are carried out, both with direct and indirect methods, demonstrating the relevance of taking into account singular arcs in the control strategy. The indirect method we use is based on our previous theoretical analysis and consists in combining a shooting method with an homotopic method. The homotopic approach leads to a quadratic regularization of the problem and is a way to tackle with the problem of nonsmoothness of the optimal control.

A passivity-based approach to stability of spatially distributed systems with a cyclic interconnection structure

Abstract: A class of distributed systems with a cyclic interconnection structure is considered. These systems arise in several biochemical applications and they can undergo diffusion driven instability which leads to a formation of spatially heterogeneous patterns. In this paper, a class of cyclic systems in which addition of diffusion does not have a destabilizing effect is identified. For these systems global stability results hold if the "secant" criterion is satisfied. In the linear case, it is shown that the secant condition is necessary and sufficient for the existence of a decoupled quadratic Lyapunov function, which extends a recent diagonal stability result to partial differential equations. For reaction-diffusion equations with nondecreasing coupling nonlinearities global asymptotic stability of the origin is established. All of the derived results remain true for both linear and nonlinear positive diffusion terms. Similar results are shown for compartmental systems.

Geometric structure-preserving optimal control of the rigid body

Abstract: In this paper we study a discrete variational optimal control problem for the rigid body. The cost to be minimized is the external torque applied to move the rigid body from an initial condition to a pre-specified terminal condition. Instead of discretizing the equations of motion, we use the discrete equations obtained from the discrete Lagrange--d'Alembert principle, a process that better approximates the equations of motion. Within the discrete-time setting, these two approaches are not equivalent in general. The kinematics are discretized using a natural Lie-algebraic formulation that guarantees that the flow remains on the Lie group SO(3) and its algebra so(3). We use Lagrange's method for constrained problems in the calculus of variations to derive the discrete-time necessary conditions. We give a numerical example for a three-dimensional rigid body maneuver.

Nonlinear optimal control via occupation measures and LMI-relaxations

Abstract: We consider the class of nonlinear optimal control problems (OCP) with polynomial data, i.e., the differential equation, state and control con- straints and cost are all described by polynomials, and more generally for OCPs with smooth data. In addition, state constraints as well as state and/or action constraints are allowed. We provide a simple hierarchy of LMI (lin- ear matrix inequality)-relaxations whose optimal values form a nondecreasing sequence of lower bounds on the optimal value. Under some convexity assump- tions, the sequence converges to the optimal value of the OCP. Preliminary results show that good approximations are obtained with few moments.

Global optimal control of perturbed systems

Abstract: We propose a new numerical method for the computation of the optimal value function of perturbed control systems and associated globally stabilizing optimal feedback controllers. The method is based on a set oriented discretization of state space in combination with a new algorithm for the computation of shortest paths in weighted directed hypergraphs. Using the concept of a multivalued game, we prove convergence of the scheme as the discretization parameter goes to zero.

A higher order model for image restoration: the one dimensional case

Abstract: The one-dimensional version of the higher order total variation-based model for image restoration proposed by Chan, Marquina, and Mulet in [4] is analyzed. A suitable functional framework in which the minimization problem is well posed is being proposed and it is proved analytically that the higher order regularizing term prevents the occurrence of the staircase effect.

Newton's method on Gra{\ss}mann manifolds

Abstract: A general class of Newton algorithms on Gra{\ss}mann and Lagrange-Gra{\ss}mann manifolds is introduced, that depends on an arbitrary pair of local coordinates. Local quadratic convergence of the algorithm is shown under a suitable condition on the choice of coordinate systems. Our result extends and unifies previous convergence results for Newton's method on a manifold. Using special choices of the coordinates, new numerical algorithms are derived for principal component analysis and invariant subspace computations with improved computational complexity properties.

On second order shape optimization methods for electrical impedance tomography

Abstract: This paper is devoted to the analysis of a second order method for recovering the \emph{a priori} unknown shape of an inclusion $\omega$ inside a body $\Omega$ from boundary measurement. This inverse problem - known as electrical impedance tomography - has many important practical applications and hence has focussed much attention during the last years. However, to our best knowledge, no work has yet considered a second order approach for this problem. This paper aims to fill that void: we investigate the existence of second order derivative of the state $u$ with respect to perturbations of the shape of the interface $\partial\omega$, then we choose a cost function in order to recover the geometry of $\partial \omega$ and derive the expression of the derivatives needed to implement the corresponding Newton method. We then investigate the stability of the process and explain why this inverse problem is severely ill-posed by proving the compactness of the Hessian at the global minimizer.

Convergence of the Min-Sum Algorithm for Convex Optimization

Abstract: We establish that the min-sum message-passing algorithm and its asynchronous variants converge for a large class of unconstrained convex optimization problems

Maxwell strata in Euler's elastic problem

Abstract: The classical Euler's problem on stationary configurations of elastic rod with fixed endpoints and tangents at the endpoints is considered as a left-invariant optimal control problem on the group of motions of a two-dimensional plane $\E(2)$. The attainable set is described, existence and boundedness of optimal controls are proved. Extremals are parametrized by Jacobi's elliptic functions of natural coordinates induced by the flow of the mathematical pendulum on fibers of the cotangent bundle of $\E(2)$. The group of discrete symmetries of Euler's problem generated by reflections in the phase space of the pendulum is studied. The corresponding Maxwell points are completely described via the study of fixed points of this group. As a consequence, an upper bound on cut points in Euler's problem is obtained.

Conjugate points in Euler's elastic problem

Abstract: For the classical Euler's elastic problem, conjugate points are described. Inflectional elasticae admit the first conjugate point between the first and the third inflection points. All the rest elasticae do not have conjugate points.

Remarks on the calculus of variations on time scales

Abstract: The calculus of variations is a classical subject which has gain throughout the last three hundred years a level of rigor and elegance that only time can give. In this note we show that, contrary to the classical field, available formulations and results on the recent calculus of variations on time scales are still at the heuristic level.

Controlled Synchronization of One Class of Nonlinear Systems under Information Constraints

Abstract: Output feedback controlled synchronization problems for a class of nonlinear unstable systems under information constraints imposed by limited capacity of the communication channel are analyzed. A binary time-varying coder-decoder scheme is described and a theoretical analysis for multi-dimensional master-slave systems represented in Lurie form (linear part plus nonlinearity depending only on measurable outputs) is provided. An output feedback control law is proposed based on the Passification Theorem. It is shown that the synchronization error exponentially tends to zero for sufficiantly high transmission rate (channel capacity). The results obtained for synchronization problem can be extended to tracking problems in a straightforward manner, if the reference signal is described by an {external} ({exogenious}) state space model. The results are applied to controlled synchronization of two chaotic Chua systems via a communication channel with limited capacity.

Discrete Control Systems

Abstract: Discrete control systems, as considered here, refer to the control theory of discrete-time Lagrangian or Hamiltonian systems. These discrete-time models are based on a discrete variational principle, and are part of the broader field of geometric integration. Geometric integrators are numerical integration methods that preserve geometric properties of continuous systems, such as conservation of the symplectic form, momentum, and energy. They also guarantee that the discrete flow remains on the manifold on which the continuous system evolves, an important property in the case of rigid-body dynamics. In nonlinear control, one typically relies on differential geometric and dynamical systems techniques to prove properties such as stability, controllability, and optimality. More generally, the geometric structure of such systems plays a critical role in the nonlinear analysis of the corresponding control problems. Despite the critical role of geometry and mechanics in the analysis of nonlinear control systems, nonlinear control algorithms have typically been implemented using numerical schemes that ignore the underlying geometry. The field of discrete control system aims to address this deficiency by restricting the approximation to choice of a discrete-time model, and developing an associated control theory that does not introduce any additional approximation. In particular, this involves the construction of a control theory for discrete-time models based on geometric integrators that yields numerical implementations of nonlinear and geometric control algorithms that preserve the crucial underlying geometric structure.