Search

We prove that solution operators of elliptic obstacle-type variational inequalities (or, more generally, locally Lipschitz continuous functions possessing certain pointwise-a.e. convexity properties) are Newton differentiable when considered as maps between suitable Lebesgue spaces and equipped with the strong-weak Bouligand differential as a generalized set-valued derivative. It is shown that this Newton differentiability allows to solve optimal control problems with H1-cost terms and one-sided pointwise control constraints by means of a semismooth Newton method. The superlinear convergence of the resulting algorithm is proved in the infinite-dimensional setting and its mesh independence is demonstrated in numerical experiments. We expect that the findings of this paper are also helpful for the design of numerical solution procedures for quasi-variational inequalities and the optimal control of obstacle-type variational problems., Comment: Minor revisions. Published in SIAM Journal on Control and Optimization, Vol. 61, Iss. 3 (2023)

In this paper, we study optimal control problems of semilinear elliptic and parabolic equations. A tracking cost functional, quadratic in the control and state variables, is considered. No control constraints are imposed. We prove that the corresponding state equations are well posed for controls in L2. However, it is well known that in the L2 framework the mappings involved in the control problem are not Frechet differentiable in general, which makes any analysis of the optimality conditions challenging. Nevertheless, we prove that every L2 optimal control belongs to L∞, and consequently standard optimality conditions are available. The first author was supported by MCIN/AEI/10.13039/501100011033 under research project PID2020-114837GB-I00. The second author was partially supported by the German Research Foundation (DFG) under project grant Wa 3626/3-2.

Given a parameter dependent fixed point equation $x = F(x,u)$, we derive an abstract compactness principle for the fixed point map $u \mapsto x^*(u)$ under the assumptions that (i) the fixed point equation can be solved by the contraction principle and (ii) the map $u \mapsto F(x,u)$ is compact for fixed $x$. This result is applied to infinite-dimensional, semi-linear control systems and their reachable sets. More precisely, we extend a non-controllability result of Ball, Marsden, and Slemrod [1] to semi-linear systems. First we consider $L^p$-controls, $p>1$. Subsequently we analyze the case $p=1$.

Necessary conditions for asymptotic stability and stabilizability of subsets for dynamical and control systems are obtained. The main necessary condition is homotopical and is in turn used to obtain a homological one. A certain extension is ruled out. Questions are posed., Comment: Accepted to SIAM J Control and Optimization

The paper is concerned with a zero-sum Stackelberg stochastic linear-quadratic (LQ, for short) differential game over finite horizons. Under a fairly weak condition, the Stackelberg equilibrium is explicitly obtained by first solving a forward stochastic LQ optimal control problem (SLQ problem, for short) and then a backward SLQ problem. Two Riccati equations are derived in constructing the Stackelberg equilibrium. An interesting finding is that the difference of these two Riccati equations coincides with the Riccati equation associated with the zero-sum Nash stochastic LQ differential game, which implies that the Stackelberg equilibrium and the Nash equilibrium are actually identical. Consequently, the Stackelberg equilibrium admits a linear state feedback representation, and the Nash game can be solved in a leader-follower manner.

We study non-autonomous observation systems \begin{align*} \dot{x}(t) = A(t) x(t),\quad y(t) = C(t) x(t),\quad x(0) = x_0\in X, \end{align*} where $(A(t))$ is a strongly measurable family of closed operators on a Banach space $X$ and $(C(t))$ is a family of bounded observation operators from $X$ to a Banach space $Y$. Based on an abstract uncertainty principle and a dissipation estimate, we prove that the observation system satisfies a final-state observability estimate in $\mathrm{L}^r(E; Y)$ for measurable subsets $E \subseteq [0,T], T > 0$. We present applications of the above result to families $(A(t))$ of uniformly strongly elliptic differential operators as well as non-autonomous Ornstein-Uhlenbeck operators $P(t)$ on $\mathrm{L}^p(\mathbb{R}^d)$ with observation operators $C(t)u = u|_{\Omega(t)}$. In the setting of non-autonomous strongly elliptic operators, we derive necessary and sufficient geometric conditions on the family of sets $(\Omega(t))$ such that the corresponding observation system satisfies a final-state observability estimate., Comment: minor corrections

The objective of this paper is to analyze the existence of equilibria for a class of deterministic mean field games of controls. The interaction between players is due to both a congestion term and a price function which depends on the distributions of the optimal strategies. Moreover, final state and mixed state-control constraints are considered, the dynamics being nonlinear and affine with respect to the control. The existence of equilibria is obtained by Kakutani's theorem, applied to a fixed point formulation of the problem. Finally, uniqueness results are shown under monotonicity assumptions.

In this paper, we analyze and discuss the well-posedness of a new variant of the so-called sweeping process, introduced by J.J. Moreau in the early 70's \cite{More71} with motivation in plasticity theory. In this variant, the normal cone to the (mildly non-convex) prox-regular moving set $C(t)$, supposed to have an absolutely continuous variation, is perturbed by a sum of a Carath\'{e}odory mapping and an integral forcing term. The integrand of the forcing term depends on two time-variables, that is, we study a general integro-differential sweeping process of Volterra type. By setting up an appropriate semi-discretization method combined with a new Gronwall-like inequality (differential inequality), we show that the integro-differential sweeping process has one and only one absolutely continuous solution. We also establish the continuity of the solution with respect to the initial value. The results of the paper are applied to the study of nonlinear integro-differential complementarity systems which are combination of Volterra integro-differential equations with nonlinear complementarity constraints. Another application is concerned with non-regular electrical circuits containing time-varying capacitors and nonsmooth electronic device like diodes. Both applications represent an additional novelty of our paper., Comment: 36 pages, 1 figures

Based on the stochastic maximum principle for the partially coupled forward-backward stochastic control system (FBSCS for short), a modified method of successive approximations (MSA for short) is established for stochastic recursive optimal control problems. The second-order adjoint processes are introduced in the augmented Hamiltonian minimization step since the control domain is not necessarily convex. Thanks to the theory of bounded mean oscillation martingales (BMO martingales for short), we give a delicate proof of the error estimate and then prove the convergence of the modified MSA algorithm. In a special case, we obtain a logarithmic convergence rate. When the control domain is convex and compact, a sufficient condition which makes the control returned from the MSA algorithm be a near-optimal control is given for a class of linear FBSCSs., Comment: 30 pages

We extend the method of controlled Lagrangians with kinetic shaping to those mechanical systems on semidirect product Lie groups with broken symmetry, more specifically to the Euler--Poincar\'e equations with advected parameters. We find a matching condition for the controlled Lagrangian for such systems whose configuration manifold is a general semidirect product Lie group $\mathsf{G} \ltimes V$. Our motivating examples are a bottom-heavy underwater vehicle and a top spinning on a movable base. Their configuration space is the special Euclidean group $\mathsf{SE}(3) = \mathsf{SO}(3) \ltimes \mathbb{R}^{3}$, where the $\mathsf{SE}(3)$-symmetry is broken by the gravity. The controls resulting from the matching condition stabilize unstable equilibria of these examples. Furthermore, the matching helps us find additional dissipative controls that asymptotically stabilize those unstable equilibria., Comment: 29 pages, 10 figures

We present several characterizations, via some weak observability inequalities, on the complete stabilizability for a control system $[A,B]$, i.e., $y'(t)=Ay(t)+Bu(t)$, $t\geq 0$, where $A$ generates a $C_0$-semigroup on a Hilbert space $X$ and $B$ is a linear and bounded operator from another Hilbert space $U$ to $X$. We then extend the aforementioned characterizations in two directions: first, the control operator $B$ is unbounded; second, the control system is time-periodic. We also give some sufficient conditions, from the perspective of the spectral projections, to ensure the weak observability inequalities. As applications, we provide several examples, which are not null controllable, but can be verified, via the weak observability inequalities, to be completely stabilizable., 28 pages

This paper studies an infinite horizon optimal control problem for discrete-time linear systems and quadratic criteria, both with random parameters which are independent and identically distributed with respect to time. A classical approach is to solve an algebraic Riccati equation that involves mathematical expectations and requires certain statistical information of the parameters. In this paper, we propose an online iterative algorithm in the spirit of Q-learning for the situation where only one random sample of parameters emerges at each time step. The first theorem proves the equivalence of three properties: the convergence of the learning sequence, the well-posedness of the control problem, and the solvability of the algebraic Riccati equation. The second theorem shows that the adaptive feedback control in terms of the learning sequence stabilizes the system as long as the control problem is well-posed. Numerical examples are presented to illustrate our results., Comment: 24 pages, 3 figures

We consider a pointwise tracking optimal control problem for a semilinear elliptic partial differential equation. We derive the existence of optimal solutions and analyze first and, necessary and sufficient, second order optimality conditions. We devise two strategies of discretization to approximate a solution of the optimal control problem: a semidiscrete scheme where the control variable is not discretized -- the so-called variational discretization approach -- and a fully discrete scheme where the control variable is discretized with piecewise constant functions. For both solution techniques, we analyze convergence properties of discretizations and derive error estimates.

We study the temperature control problem for Langevin diffusions in the context of non-convex optimization. The classical optimal control of such a problem is of the bang-bang type, which is overly sensitive to errors. A remedy is to allow the diffusions to explore other temperature values and hence smooth out the bang-bang control. We accomplish this by a stochastic relaxed control formulation incorporating randomization of the temperature control and regularizing its entropy. We derive a state-dependent, truncated exponential distribution, which can be used to sample temperatures in a Langevin algorithm, in terms of the solution to an HJB partial differential equation. We carry out a numerical experiment on a one-dimensional baseline example, in which the HJB equation can be easily solved, to compare the performance of the algorithm with three other available algorithms in search of a global optimum.

We study controlled filter stability and its effects on the robustness properties of optimal control policies designed for systems with incorrect priors applied to a true system. Filter stability refers to the correction of an incorrectly initialized filter for a partially observed stochastic dynamical system (controlled or control-free) with increasing measurements. This problem has been studied extensively in the control-free context, and except for the standard machinery for linear Gaussian systems involving the Kalman Filter, few studies exist for the controlled setup. One of the main differences between control-free and controlled partially observed Markov chains is that the filter is always Markovian under the former, whereas under a controlled model the filter process may not be Markovian since the control policy may depend on past measurements in an arbitrary (measurable) fashion. This complicates the dependency structure and therefore results from the control-free literature do not directly apply to the controlled setup. In this paper, we study the filter stability problem for controlled stochastic dynamical systems, and provide sufficient conditions for when a falsely initialized filter merges with the correctly initialized filter over time. These stability results are applied to robust stochastic control problems: under filter stability, we bound the difference in the expected cost incurred for implementing an incorrectly designed control policy compared to an optimal policy. A conclusion is that filter stability leads to stronger robustness results to incorrect priors (compared with setups without controlled filter stability). Furthermore, if the optimum cost is that same for each prior, the cost of mismatch between the true prior and the assumed prior is zero.

Motivated by models of epidemic control in large populations, we consider a Stackelberg mean field game model between a principal and a mean field of agents evolving on a finite state space. The agents play a non-cooperative game in which they can control their transition rates between states to minimize an individual cost. The principal can influence the resulting Nash equilibrium through incentives so as to optimize its own objective. We analyze this game using a probabilistic approach. We then propose an application to an epidemic model of SIR type in which the agents control their interaction rate and the principal is a regulator acting with non pharmaceutical interventions. To compute the solutions, we propose an innovative numerical approach based on Monte Carlo simulations and machine learning tools for stochastic optimization. We conclude with numerical experiments by illustrating the impact of the agents' and the regulator's optimal decisions in two models: a basic SIR model with semi-explicit solutions and a more complex model with a larger state space., 31 pages, 10 figures, 3 algorithms, 3 tables

We consider nonlinear filters for diffusion processes when the observation and signal noises are small and of the same order. As the noise intensities approach zero, the nonlinear filter can be approximated by a certain variational problem that is closely related to Mortensen's optimization problem(1968). This approximation result can be made precise through a certain Laplace asymptotic formula. In this work we study probabilities of deviations of true filtering estimates from that obtained by solving the variational problem. Our main result gives a large deviation principle for Laplace functionals whose typical asymptotic behavior is described by Mortensen-type variational problems. Proofs rely on stochastic control representations for positive functionals of Brownian motions and Laplace asymptotics of the Kallianpur-Striebel formula., 28 pages

We consider a general formulation of the random horizon Principal-Agent problem with a continuous payment and a lump-sum payment at termination. In the European version of the problem, the random horizon is chosen solely by the principal with no other possible action from the agent than exerting effort on the dynamics of the output process. We also consider the American version of the contract, which covers the seminal Sannikov's model, where the agent can also quit by optimally choosing the termination time of the contract. Our main result reduces such non-zero-sum stochastic differential games to appropriate stochastic control problems which may be solved by standard methods of stochastic control theory. This reduction is obtained by following Sannikov's approach, further developed by Cvitanic, Possamai, and Touzi. We first introduce an appropriate class of contracts for which the agent's optimal effort is immediately characterized by the standard verification argument in stochastic control theory. We then show that this class of contracts is dense in an appropriate sense so that the optimization over this restricted family of contracts represents no loss of generality. The result is obtained by using the recent well-posedness result of random horizon second-order backward SDE.

In optimal control theory, infimum gap means that there is a gap between the infimum values of a given minimum problem and an extended problem, obtained by enlarging the set of original solutions and controls. The gap phenomenon is somewhat "dual" to the problem of the controllability of the original control system to an extended solution. In this paper we present sufficient conditions for the absence of an infimum gap and for controllability for a wide class of optimal control problems subject to endpoint and state constraints. These conditions are based on a nondegenerate version of the nonsmooth constrained maximum principle, expressed in terms of subdifferentials. In particular, under some new constraint qualification conditions, we prove that: (i) if an extended minimizer is a nondegenerate normal extremal, then no gap shows up; (ii) given an extended solution verifying the constraints, either it is a nondegenerate abnormal extremal, or the original system is controllable to it. An application to the impulsive extension of a free end-time, non-convex optimization problem with control-polynomial dynamics illustrates the results.

This paper presents the analysis of the stability properties of PID controllers for dynamical systems with multiple state delays, focusing on the mathematical characterization of the potential sensitivity of stability with respect to infinitesimal parametric perturbations. These perturbations originate for instance from neglecting feedback delay, a finite difference approximation of the derivative action, or neglecting fast dynamics. The analysis of these potential sensitivity problems leads us to the introduction of a `robustified' notion of stability called \emph{strong stability}, inspired by the corresponding notion for neutral functional differential equations. We prove that strong stability can be achieved by adding a low-pass filter with a sufficiently large cut-off frequency to the control loop, on the condition that the filter itself does not destabilize the nominal closed-loop system. Throughout the paper, the theoretical results are illustrated by examples that can be analyzed analytically, including, among others, a third-order unstable system where both proportional and derivative control action are necessary for achieving stability, while the regions in the gain parameter-space for stability and strong stability are not identical. Besides the analysis of strong stability, a computational procedure is provided for designing strongly stabilizing PID controllers. Computational case-studies illustrating this design procedure complete the presentation.

In this paper we study the problem of the optimal distribution of two materials on smooth submanifolds $M$ of dimension $d-1$ in $\mathbf R^d$ without boundary by means of the topological derivative. We consider a class of shape optimisation problems which are constrained by a linear partial differential equation on the surface. We examine the singular perturbation of the differential operator and material coefficients and derive the topological derivative. Finally, we show how the topological derivative in conjunction with a level set method on the surface can be used to solve the topology optimisation problem numerically.

This paper studies system theoretic properties of the class of difference inclusions of convex processes. We will develop a framework considering eigenvalues and eigenvectors, weakly and strongly invariant cones, and a decomposition of convex processes. This framework will lead to characterizations of reachability, stabilizability, and (null-)controllability for nonstrict convex processes in terms of spectral properties. These characterizations generalize all previously known results regarding, for instance, linear processes and specific classes of nonstrict convex processes.

We investigate the evolution of epidemics over dynamical networks when nodes choose to interact with others in a selfish and decentralized manner. Specifically, we analyze the susceptible-asymptomatic-infected-recovered (SAIR) epidemic in the framework of activity-driven networks with heterogeneous node degrees and time-varying activation rates, and derive both individual and degree-based mean-field approximations of the exact state evolution. We then present a game-theoretic model where nodes choose their activation probabilities in a strategic manner using current state information as feedback, and characterize the quantal response equilibrium (QRE) of the proposed setting. We then consider the activity-driven susceptible-infected-susceptible (SIS) epidemic model, characterize equilibrium activation probabilities and analyze epidemic evolution in closed-loop. Our numerical results provide compelling insights into epidemic evolution under game-theoretic activation. Specifically, for the SAIR epidemic, we show that under suitable conditions, the epidemic can persist, as any decrease in infected proportion is counteracted by an increase in activity rates by the nodes. For the SIS epidemic, we show that in regimes where there is an endemic state, the infected proportion could be significantly smaller under game-theoretic activation if the loss upon infection is sufficiently high.

This paper studies a continuous-review backlogged inventory model considered by Helmes et al. (2015) but with discontinuous quantity-dependent setup cost for each order. In particular, the setup cost is characterized by a two-step function and a higher cost would be charged once the order quantity exceeds a threshold $Q$. Unlike the optimality of $(s,S)$-type policy obtained by Helmes et al. (2015) for continuous setup cost with the discounted cost criterion, we find that, in our model, although some $(s,S)$-type policy is indeed optimal in some cases, the $(s,S)$-type policy can not always be optimal. In particular, we show that there exist cases in which an $(s,S)$ policy is optimal for some initial levels but it is strictly worse than a generalized $(s,\{S(x):x\leq s\})$ policy for the other initial levels. Under $(s,\{S(x):x\leq s\})$ policy, it orders nothing for $x>s$ and orders up to level $S(x)$ for $x\leq s$, where $S(x)$ is a non-constant function of $x$. We further prove the optimality of such $(s,\{S(x):x\leq s\})$ policy in a large subset of admissible policies for those initial levels. Moreover, the optimality is obtained through establishing a more general lower bound theorem which will also be applicable in solving some other optimization problems by the common lower bound approach., 50 pages

This paper examines mean field linear-quadratic-Gaussian (LQG) social optimum control with volatility-uncertain common noise. The diffusion terms in the dynamics of agents contain an unknown volatility process driven by a common noise. We apply a robust optimization approach in which all agents view volatility uncertainty as an adversarial player. Based on the principle of person-by-person optimality and a two-step-duality technique for stochastic variational analysis, we construct an auxiliary optimal control problem for a representative agent. Through solving this problem combined with a consistent mean field approximation, we design a set of decentralized strategies, which are further shown to be asymptotically social optimal by perturbation analysis.

The paper studies the open-loop saddle point and the open-loop lower and upper values, as well as their relationship for two-person zero-sum stochastic linear-quadratic (LQ, for short) differential games with deterministic coefficients. It derives a necessary condition for the finiteness of the open-loop lower and upper values and a sufficient condition for the existence of an open-loop saddle point. It turns out that under the sufficient condition, a strongly regular solution to the associated Riccati equation uniquely exists, in terms of which a closed-loop representation is further established for the open-loop saddle point. Examples are presented to show that the finiteness of the open-loop lower and upper values does not ensure the existence of an open-loop saddle point in general. But for the classical deterministic LQ game, these two issues are equivalent and both imply the solvability of the Riccati equation, for which an explicit representation of the solution is obtained., Comment: 26 pages

We introduce new parametrized classes of shape admissible domains in R^n , n $\ge$ 2, and prove that they are compact with respect to the convergence in the sense of characteristic functions, the Hausdorff sense, the sense of compacts and the weak convergence of their boundary volumes. The domains in these classes are bounded ($\epsilon$, $\infty$)-domains with possibly fractal boundaries that can have parts of any non-uniform Hausdorff dimension greater or equal to n -- 1 and less than n. We prove the existence of optimal shapes in such classes for maximum energy dissipation in the framework of linear acous-tics. A by-product of our proof is the result that the class of bounded ($\epsilon$, $\infty$)-domains with fixed $\epsilon$ is stable under Hausdorff convergence. An additional and related result is the Mosco convergence of Robin-type energy functionals on converging domains.

We discuss a class of explicitly solvable mean field type control problems/mean field games with a clear economic interpretation. More precisely, we consider long term average impulse control problems with underlying general one-dimensional diffusion processes motivated by optimal harvesting problems in natural resource management. We extend the classical stochastic Faustmann models by allowing the prices to depend on the state of the market using a mean field structure. In a competitive market model, we prove that, under natural conditions, there exists an equilibrium strategy of threshold-type and furthermore characterize the threshold explicitly. If the agents cooperate with each other, we are faced with the mean field type control problem. Using a Lagrange-type argument, we prove that the optimizer of this non-standard impulse control problem is of threshold-type as well and characterize the optimal threshold. Furthermore, we compare the solutions and illustrate the findings in an example.

In statistical decision theory involving a single decision-maker, an information structure is said to be better than another one if for any cost function involving a hidden state variable and an action variable which is restricted to be conditionally independent from the state given some measurement, the solution value under the former is not worse than that under the latter. For finite spaces, a theorem due to Blackwell leads to a complete characterization on when one information structure is better than another. For stochastic games, in general, such an ordering is not possible since additional information can lead to equilibria perturbations with positive or negative values to a player. However, for zero-sum games in a finite probability space, P\k{e}ski introduced a complete characterization of ordering of information structures. In this paper, we obtain an infinite dimensional (standard Borel) generalization of P\k{e}ski's result. A corollary is that more information cannot hurt a decision maker taking part in a zero-sum game. We establish two supporting results which are essential and explicit though modest improvements on prior literature: (i) a partial converse to Blackwell's ordering in the standard Borel setup and (ii) an existence result for equilibria in zero-sum games with incomplete information., Comment: To appear in SIAM Journal on Control and Optimization. A short conference version appeared at the 2020 IEEE International Symposium on Information Theory

An optimal ergodic control problem (EC problem, for short) is investigated for a linear stochastic differential equation with quadratic cost functional. Constant nonhomogeneous terms, not all zero, appear in the state equation, which lead to the asymptotic limit of the state non-zero. Under the stabilizability condition, for any (admissible) closed-loop strategy, an invariant measure is proved to exist, which makes the ergodic cost functional well-defined and the EC problem well-formulated. Sufficient conditions, including those allowing the weighting matrices of cost functional to be indefinite, are introduced for finiteness and solvability for the EC problem. Some comparisons are made between the solvability of EC problem and the closed-loop solvability of stochastic linear quadratic optimal control problem in the infinite horizon. Regularized EC problem is introduced to be used to obtain the optimal value of the EC problem.

Problems involving control of large ensmebles of structurally identical dynamical systems, called \emph{ensemble control}, arise in numerous scientific areas from quantum control and robotics to brain medicine. In many of such applications, control can only be implemented at the population level, i.e., through broadcasting an input signal to all the systems in the population, and this new control paradigm challenges the classical systems theory. In recent years, considerable efforts have been made to investigate controllability properties of ensemble systems, and most works emphasized on linear and some forms of bilinear and nonlinear ensemble systems. In this paper, we study controllability of a broad class of bilinear ensemble systems defined on semisimple Lie groups, for which we define the notion of ensemble controllability through a Riemannian structure of the state space Lie group. Leveraging the Cartan decomposition of semisimple Lie algebras in representation theory, we develop a \emph{covering method} that decomposes the state space Lie group into a collection of Lie subgroups generating the Lie group, which enables the determination of ensemble controllability by controllability of the subsystems evolving on these Lie subgroups. Using the covering method, we show the equivalence between ensemble and classical controllability, i.e., controllability of each individual system in the ensemble implies ensemble controllability, for bilinear ensemble systems evolving on semisimple Lie groups. This equivalence makes the examination of controllability for infinite-dimensional ensemble systems as tractable as for a finite-dimensional single system., Comment: keywords: Ensemble control, Semisimple Lie groups, Approximation theory, Controllability, Infinite-dimensional Systems

We introduce a mean field game model for pedestrians moving in a given domain and choosing their trajectories so as to minimize a cost including a penalization on the difference between their own velocity and that of the other agents they meet. We prove existence of an equilibrium in a Lagrangian setting by using its variational structure, and then study its properties and regularity.

International audience; We consider a continuous time stochastic optimal control problem under both equality and inequality constraints on the expectation of some functionals of the controlled process. Under a qualification condition, we show that the problem is in duality with an optimization problem involving the Lagrange multiplier associated with the constraints. Then by convex analysis techniques, we provide a general existence result and some a priori estimation of the dual optimizers. We further provide a necessary and sufficient optimality condition for the initial constrained control problem. The same results are also obtained for a discrete time constrained control problem. Moreover, under additional regularity conditions, it is proved that the discrete time control problem converges to the continuous time problem, possibly with a convergence rate. This convergence result can be used to obtain numerical algorithms to approximate the continuous time control problem, which we illustrate by two simple numerical examples.

In this paper, we provide a system theoretic treatment of a new class of multilinear time-invariant (MLTI) systems in which the states, inputs and outputs are tensors, and the system evolution is governed by multilinear operators. The MLTI system representation is based on the Einstein product and even-order paired tensors. There is a particular tensor unfolding which gives rise to an isomorphism from this tensor space to the general linear group, i.e. the group of invertible matrices. By leveraging this unfolding operation, one can extend classical linear time-invariant (LTI) system notions including stability, reachability and observability to MLTI systems. While the unfolding based formulation is a powerful theoretical construct, the computational advantages of MLTI systems can only be fully realized while working with the tensor form, where hidden patterns/structures can be exploited for efficient representations and computations. Along these lines, we establish new results which enable one to express tensor unfolding based stability, reachability and observability criteria in terms of more standard notions of tensor ranks/decompositions. In addition, we develop a generalized CANDECOMP/PARAFAC decomposition and tensor train decomposition based model reduction framework, which can significantly reduce the number of MLTI system parameters. We demonstrate our framework with numerical examples., 27 pages, 2 figures, 3 tables, SIAM Journal on Control and Optimization, accepted to appear. arXiv admin note: text overlap with arXiv:1905.07427

This paper studies the graph-theoretic conditions for stability of positive monotone systems. Using concepts from input-to-state stability and network small-gain theory, we first establish necessary and sufficient conditions for the stability of linear positive systems described by Metzler matrices. Specifically, we derive two sets of stability conditions based on two forms of input-to-state stability gains for Metzler systems, namely max-interconnection gains and sum-interconnection gains. Based on the max-interconnection gains, we show that the cyclic small-gain theorem becomes necessary and sufficient for the stability of Metzler systems; based on the sum-interconnection gains, we obtain novel graph-theoretic conditions for the stability of Metzler systems. All these conditions highlight the role of cycles in the interconnection graph and unveil how the structural properties of the graph affect stability. Finally, we extend our results to the nonlinear monotone system and obtain similar sufficient conditions for global asymptotic stability.

Controlling the shape and position of moving and pinned droplets on a solid surface is an important feature often found in microfluidic applications. In this work, we consider a well investigated phase field model including contact line dynamics as the state system for an (open-loop) optimal control problem. Here the spatially and temporally changeable contact angles between droplet and solid are considered as the control variables. We consider a suitable, energy stable, time discrete version of the state equation in our optimal control problem. We discuss regularity of the solution to the time discrete state equation and its continuity and differentiability properties. Furthermore, we show existence of solutions and state first order optimality conditions to the optimal control problem. We illustrate our results by actively pushing a droplet uphill against gravity in an optimal way.

In this paper, we propose a second-order continuous primal-dual dynamical system with time-dependent positive damping terms for a separable convex optimization problem with linear equality constraints. By the Lyapunov function approach, we investigate asymptotic properties of the proposed dynamical system as the time $t\to+\infty$. The convergence rates are derived for different choices of the damping coefficients. We also show that the obtained results are robust under external perturbations.

A linear control system with quadratic cost functional over infinite time horizon is considered without assuming controllability/stabilizability condition and the global integrability condition for the nonhomogeneous term of the state equation and the weight functions in the linear terms in the running cost rate function. Classical approaches do not apply for such kind of problems. Existence and non-existence of overtaking optimal controls in various cases are established. Some concrete examples are presented. These results show that the overtaking optimality approach can be used to solve some of the above-mentioned problems and at the same time, the limitation of this approach is also revealed.

We prove a small-gain theorem for interconnections of $n$ nonlinear heterogeneous input-to-state stable (ISS) control systems of a general nature, covering partial, delay and ordinary differential equations. Furthermore, for the same class of control systems, we derive small-gain theorems for asymptotic gain, uniform global stability and weak input-to-state stability properties. We show that our technique is applicable for different formulations of ISS property (summation, maximum, semimaximum) and discuss tightness of achieved small-gain theorems. Finally, we introduce variations of uniform asymptotic gain and uniform limit properties, which are particularly useful for small-gain arguments and characterize ISS in terms of these notions., Accepted to SICON, 2021

Suppose that a continuous-time linear infinite-dimensional system with a static state-feedback controller is strongly stable. We address the following question: If we convert the continuous-time controller to a sampled-data controller by applying an idealized sampler and a zero-order hold, will the resulting sampled-data system be strongly stable for all sufficiently small sampling periods? In this paper, we restrict our attention to the situation where the generator of the open-loop system is a Riesz-spectral operator and its point spectrum has a limit point at the origin. We present conditions under which the answer to the above question is affirmative. In the robustness analysis, we show that the sufficient condition for strong stability obtained in the Arendt-Batty-Lyubich-V\~u theorem is preserved between the original continuous-time system and the sampled-data system under fast sampling., Comment: 26 pages, 2 figures

In this paper, we study a generalized Kalman-Bucy filtering problem under uncertainty. The drift uncertainty for both signal process and observation process is considered and the attitude to uncertainty is characterized by a convex operator (convex risk measure). The optimal filter or the minimum mean square estimator (MMSE) is calculated by solving the minimum mean square estimation problem under a convex operator. In the first part of this paper, this estimation problem is studied under g-expectation which is a special convex operator. For this case, we prove that there exists a worst-case prior. Based on this worst-case prior we obtained the Kalman-Bucy filtering equation under g-expectation. In the second part of this paper, we study the minimum mean square estimation problem under general convex operators. The existence and uniqueness results of the MMSE are deduced., Comment: 24 pages

We discuss a topology optimization problem for an elastoplastic medium. The distribution of material in a region is optimized with respect to a given target functional taking into account compliance. The incremental elastoplastic problem serves as state constraint. We prove that the topology optimization problem admits a solution. First-order optimality conditions are obtained by considering a regularized problem and passing to the limit.

In this paper, we study an optimal control problem for a macroscopic mechanical tumour model based on the phase field approach. The model couples a Cahn--Hilliard type equation to a system of linear elasticity and a reaction-diffusion equation for a nutrient concentration. By taking advantage of previous analytical well-posedness results established by the authors, we seek optimal controls in the form of a boundary nutrient supply, as well as concentrations of cytotoxic and antiangiogenic drugs that minimise a cost functional involving mechanical stresses. Special attention is given to sparsity effects, where with the inclusion of convex non-differentiable regularisation terms to the cost functional, we can infer from the first-order optimality conditions that the optimal drug concentrations can vanish on certain time intervals., 26 pages

In this paper, we consider a stochastic recursive optimal control problem under model uncertainty. In this framework, the cost function is described by solutions of a family of backward stochastic differential equations. With the help of the linearization techniques and weak convergence methods, we derive the corresponding stochastic maximum principle. Moreover, a linear quadratic robust control problem is also studied., 28 pages

In this paper we establish an observability inequality for the heat equation with bounded potentials on the whole space. Roughly speaking, such a kind of inequality says that the total energy of solutions can be controlled by the energy localized in a subdomain, which is equidistributed over the whole space. The proof of this inequality is mainly adapted from the parabolic frequency function method, which plays an important role in proving the unique continuation property for solutions of parabolic equations. As an immediate application, we show that the null controllability holds for the heat equation with bounded potentials on the whole space., Comment: Welcome for any comment

We address the variational formulation of the risk-sensitive reward problem for non-degenerate diffusions on $\mathbb{R}^d$ controlled through the drift. We establish a variational formula on the whole space and also show that the risk-sensitive value equals the generalized principal eigenvalue of the semilinear operator. This can be viewed as a controlled version of the variational formulas for principal eigenvalues of diffusion operators arising in large deviations. We also revisit the average risk-sensitive minimization problem and by employing a gradient estimate developed in this paper, we extend earlier results to unbounded drifts and running costs., Comment: 29 pages

Let $X,Y$ be Banach spaces, $(S_t)_{t \geq 0}$ a $C_0$-semigroup on $X$, $-A$ the corresponding infinitesimal generator on $X$, $C$ a bounded linear operator from $X$ to $Y$, and $T > 0$. We consider the system \[ \dot{x}(t) = -Ax(t), \quad y(t) = Cx(t) \quad t\in (0,T], \quad x(0) = x_0 \in X. \] We provide sufficient conditions such that this system satisfies a final state observability estimate in $L_r ((0,T) ; Y)$, $r \in [1,\infty]$. These sufficient conditions are given by an uncertainty relation and a dissipation estimate. Our approach unifies and generalizes the respective advantages from earlier results obtained in the context of Hilbert spaces. As an application we consider the example where $A$ is an elliptic operator in $L_p(\mathbb{R}^d)$ for $1

In this article, the shape optimization of a linear elastic body subject to frictional (Tresca) contact is investigated. Due to the projection operators involved in the formulation of the contact problem, the solution is not shape differentiable in general. Moreover, shape optimization of the contact zone requires the computation of the gap between the bodies in contact, as well as its shape derivative. Working with directional derivatives, sufficient conditions for shape differentiability are derived. %The problem is addressed in the general framework of two bodies with smooth boundaries. Then, some numerical results, obtained with a gradient descent algorithm based on those shape derivatives, are presented., Comment: Preprint

We study the controlled dynamics of the {\it ensembles of points} of a Riemannian manifold $M$. Parameterized ensemble of points of $M$ is the image of a continuous map $\gamma:\Theta \to M$, where $\Theta$ is a compact set of parameters. The dynamics of ensembles is defined by the action $\gamma(\theta) \mapsto P_t(\gamma(\theta))$ of the semigroup of diffeomorphisms $P_t:M \to M, \ t \in \mathbb{R}$, generated by the controlled equation $\dot{x}=f(x,u(t))$ on $M$. Therefore any control system on $M$ defines a control system on (generally infinite-dimensional) space $\mathcal{E}_\Theta(M)$ of the ensembles of points. We wish to establish criteria of controllability for such control systems. As in our previous work ([1]) we seek to adapt the Lie-algebraic approach of geometric control theory to the infinite-dimensional setting. We study the case of finite ensembles and prove genericity of exact controllability property for them. We also find sufficient approximate controllability criterion for continual ensembles and prove a result on motion planning in the space of flows on $M$. We discuss the relation of the obtained controllability criteria to various versions of Rashevsky-Chow theorem for finite- and infinite-dimensional manifolds., Comment: 24 pages

For a class of quasi-variational inequalities (QVIs) of obstacle-type the stability of its solution set and associated optimal control problems are considered. These optimal control problems are non-standard in the sense that they involve an objective with set-valued arguments. The approach to study the solution stability is based on perturbations of minimal and maximal elements of the solution set of the QVI with respect to {monotone} perturbations of the forcing term. It is shown that different assumptions are required for studying decreasing and increasing perturbations and that the optimization problem of interest is well-posed., Comment: 29 pages