2025-02-14 | | Total: 29
In this paper we address the speed planning problem for a vehicle over an assigned path with the aim of minimizing a weighted sum of travel time and energy consumption under suitable constraints (maximum allowed speed, maximum traction or braking force, maximum power consumption). The resulting mathematical model is a non--convex optimization problem. We prove that, under some mild assumptions, a convex reformulation of the non--convex problem is exact. In particular, the convex reformulation is a Second Order Cone Programming (SOCP) problem, for which efficient solvers exist. Through the numerical experiments we confirm that the convex relaxation can be solved very efficiently and, moreover, we also provide the Pareto front of the trade-off between the two objectives (travel time and energy consumption).
Achieving optimality in controlling physical systems is a profound challenge across diverse scientific and engineering fields, spanning neuromechanics, biochemistry, autonomous systems, economics, and beyond. Traditional solutions, relying on time-consuming offline iterative algorithms, often yield limited insights into fundamental natural processes. In this work, we introduce a novel, causally deterministic approach, presenting the closed-form optimal tracking controller (OTC) that inherently solves pseudoconvex optimization problems in various fields. Through rigorous analysis and comprehensive numerical examples, we demonstrate OTC's capability of achieving both high accuracy and rapid response, even when facing high-dimensional and high-dynamical real-world problems. Notably, our OTC outperforms state-of-the-art methods by, e.g., solving a 1304-dimensional neuromechanics problem 1311 times faster or with 113 times higher accuracy. Most importantly, OTC embodies a causally deterministic system interpretation of optimality principles, providing a new and fundamental perspective of optimization in natural and artificial processes. We anticipate our work to be an important step towards establishing a general causally deterministic optimization theory for a broader spectrum of system and problem classes, promising advances in understanding optimality principles in complex systems.
In "chemostat"-type population models that incorporate substrate (nutrient) dynamics, the dependence of the birth (or growth) rate on the substrate concentration introduces nonlinear coupling that creates a challenge for stabilization that is global, namely, for all positive concentrations of the biomass and nutrients. This challenge for global stabilization has been overcome in the literature using relatively simple feedback when natural mortality of the biomass is absent. However, under natural mortality, it takes fortified, more complex feedback, outside of the existing nonlinear control design toolbox, to avoid biomass extinction from nutrient-depleted initial conditions. Such fortified feedback, the associated control Laypunov function design, and Lyapunov analysis of global stability are provided in this paper. We achieve global stabilization for two different chemostat models: (i) a lumped model, with two state variables, and (ii) a three-state model derived from an age-structured infinite-dimensional model. The proposed feedback stabilizers are explicit, applicable to both the lumped and the age-structured models, and coincide with simple feedback laws proposed in the literature when the mortality rate is zero. Global stabilization means subject to constraints: all positive biomass and nutrient concentrations are within the region of attraction of the desired equilibrium, and, additionally, this is achieved with a dilution input that is guaranteed to remain positive. For the lumped case with Haldane kinetics, we show that the reproduction rate dominating the mortality (excluding the reproduction and mortality being in balance) is not only sufficient but also necessary for global stabilization. The obtained results are illustrated with simple examples.
Electric vehicle (EV) coordination can provide significant benefits through vehicle-to-everything (V2X) by interacting with the grid, buildings, and other EVs. This work aims to develop a V2X value-stacking framework, including vehicle-to-building (V2B), vehicle-to-grid (V2G), and energy trading, to maximize economic benefits for residential communities while maintaining distribution voltage. This work also seeks to quantify the impact of prediction errors related to building load, renewable energy, and EV arrivals. A dynamic rolling-horizon optimization (RHO) method is employed to leverage multiple revenue streams and maximize the potential of EV coordination. To address energy uncertainties, including hourly local building load, local photovoltaic (PV) generation, and EV arrivals, this work develops a Transformer-based forecasting model named Gated Recurrent Units-Encoder-Temporal Fusion Decoder (GRU-EN-TFD). The simulation results, using real data from Australia's National Electricity Market, and the Independent System Operators in New England and New York in the US, reveal that V2X value stacking can significantly reduce energy costs. The proposed GRU-EN-TFD model outperforms the benchmark forecast model. Uncertainties in EV arrivals have a more substantial impact on value-stacking performance, highlighting the significance of its accurate forecast. This work provides new insights into the dynamic interactions among residential communities, unlocking the full potential of EV batteries.
This paper, in the setting at infinity, presents some relationships between the modulus of metric regularity and the radius of (strong) metric regularity that gives a measure of the extent to which a set-valued mapping can be perturbed before (strong) metric regularity is lost. The results given here can be viewed as versions at infinity of [2, Theorem 1.5] and [3, Theorem 4.6].
The Gromov-Wasserstein (GW) problem is a variant of the classical optimal transport problem that allows one to compute meaningful transportation plans between incomparable spaces. At an intuitive level, it seeks plans that minimize the discrepancy between metric evaluations of pairs of points. The GW problem is typically cast as an instance of a non-convex quadratic program that is, unfortunately, intractable to solve. In this paper, we describe tractable semidefinite relaxations of the GW problem based on the Sum-of-Squares (SOS) hierarchy. We describe how the Putinar-type and the Schmüdgen-type moment hierarchies can be simplified using marginal constraints, and we prove convergence rates for these hierarchies towards computing global optimal solutions to the GW problem. The proposed SOS hierarchies naturally induce a distance measure analogous to the distortion metrics, and we show that these are genuine distances in that they satisfy the triangle inequality. In particular, the proposed SOS hierarchies provide computationally tractable proxies of the GW distance and the associated distortion distances (over metric measure spaces) that are otherwise intractable to compute.
Morse parametric qualification condition is a new condition that generalizes uniform strong convexity. Generic semi-algebraic functions are Morse parametric in a piecewise sense, implying that generic semi-algebraic bilevel problems reduce to mixed-integer programming. In this new framework, we study bilevel gradient algorithms with two strategies: the single-step multi-step strategy, which involves a sequence of steps on the lower-level problems followed by one step on the upper-level problem, and a differentiable programming strategy that optimizes a smooth approximation of the bilevel problem. While the first is shown to be a biased gradient method on the problem with rich properties, the second, inspired by meta-learning applications, is less stable but offers simplicity and ease of implementation.
In this paper, we investigate the multi-objective optimal control problem of ordinary differential equations on Riemannian manifolds. We first obtain the second-order necessary conditions for weak Pareto optimal solutions for multi-objective optimal control problems with fixed terminal time, and then extend these results to multi-objective optimal control problems with free terminal time, deriving the corresponding second-order necessary conditions for weak Pareto optimal solutions. Our main results (i.e., Theorem 2.2 and Theorem 2.4) show that weak Pareto optimal solutions depend on the curvature tensor of the Riemannian manifold. Finally, we provide an example (i.e., Example 2.1) as an application of our main results, illustrating how our findings differ from existing related results (see Remark 2.2).
This paper introduces a robust two-timescale compressed primal-dual (TiCoPD) algorithm tailored for decentralized optimization under bandwidth-limited and unreliable channels. By integrating the majorization-minimization approach with the primal-dual optimization framework, the TiCoPD algorithm strategically compresses the difference term shared among agents to enhance communication efficiency and robustness against noisy channels without compromising convergence stability. The method incorporates a mirror sequence for agent consensus on nonlinearly compressed terms updated on a fast timescale, together with a slow timescale primal-dual recursion for optimizing the objective function. Our analysis demonstrates that the proposed algorithm converges to a stationary solution when the objective function is smooth but possibly non-convex. Numerical experiments corroborate the conclusions of this paper.
We study the asymptotic behavior of solutions to linear-quadratic mean field stochastic optimal control problems. By formulating an ergodic control framework, we characterize the convergence between the finite time horizon control problem and its ergodic counterpart. Leveraging these convergence results, we establish the turnpike property for the optimal pairs, demonstrating that solutions to the finite time horizon control problem remain exponentially close to the ergodic equilibrium except near the temporal boundaries. This result reveals the intrinsic connection between long-term dynamics and their asymptotic behavior in mean field control systems.
This paper presents a comprehensive theoretical analysis of six distinct Mixed-Integer Programming (MIP) formulations for preventive Generator Maintenance Scheduling (GMS), a critical problem for ensuring the reliability and efficiency of power systems. By comparing the tightness of their linear relaxations, we identify which formulations offer superior dual bound and, thus, better computational performance. Our analysis includes establishing relationships between the formulations through definitions, lemmas, and propositions, demonstrating that some formulations provide tighter relaxations that lead to more efficient optimization outcomes. These findings offer valuable insights for practitioners and researchers in selecting the most effective models to enhance the scheduling process of preventive generator maintenance.
In this work, we initiate the research on controlling nonlinear waves propagating on lattices from a completely new perspective. We consider nonlinear waves on a lattice as a system of interacting particles and study their collective flocking behavior. By designing suitable feedback controls, we show that any admissible flock can be reached within a finite amount of time. Finally, we highlight the connection between our flocking problem and a minimal-time problem in the framework of nonlinear Hamilton-Jacobi equations and optimal control theory.
We propose a new bundle-based augmented Lagrangian framework for solving constrained convex problems. Unlike the classical (inexact) augmented Lagrangian method (ALM) that has a nested double-loop structure, our framework features a $\textit{single-loop}$ process. Motivated by the proximal bundle method (PBM), we use a $\textit{bundle}$ of past iterates to approximate the subproblem in ALM to get a computationally efficient update at each iteration. We establish sub-linear convergences for primal feasibility, primal cost values, and dual iterates under mild assumptions. With further regularity conditions, such as quadratic growth, our algorithm enjoys $\textit{linear}$ convergences. Importantly, this linear convergence can happen for a class of conic optimization problems, including semidefinite programs. Our proof techniques leverage deep connections with inexact ALM and primal-dual principles with PBM.
We propose the first analytical stochastic model for optimizing the configuration and implementation policies of fare-free transit. The model focuses on a transportation corridor with two transportation modes: automobiles and buses. The corridor is divided into two sections, an inner one with fare-free transit service and an outer one with fare-based transit service. Under the static version of the model, the optimized length and frequency of the fare-free transit zone can be determined by maximizing total social welfare. The findings indicate that implementing fare-free transit can increase transit ridership and reduce automobile use within the fare-free zone while social equity among the demand groups can be enhanced by lengthening the fare-free zone. Notably, the optimal zone length increases when both social welfare and equity are considered jointly, compared to only prioritizing social welfare. The dynamic model, framed within a market entry and exit real options approach, solves the fare policy switching problem, establishing optimal timing policies for activating or terminating fare-free service. The results from dynamic models reveal earlier implementation and extended durations of fare-free transit in the social welfare-aware regime, driven by lower thresholds compared to the social equity-aware regime.
Efficient management of transportation corridors is critical for sustaining urban mobility, directly influencing transportation efficiency. Two prominent strategies for enhancing public transit services and alleviating congestion, Exclusive Bus Lane (EBL) and High Occupancy Vehicle Lane (HOVL), are gaining increasing attention. EBLs prioritize bus transit by providing dedicated lanes for faster travel times, while HOVLs encourage carpooling by reserving lanes for high-occupancy vehicles. However, static implementations of these policies may underutilize road resources and disrupt general-purpose lanes. Dynamic control of these policies, based on real-time demand, can potentially maximize road efficiency and minimize negative impacts. This study develops cost functions for Mixed Traffic Policy (MTP), Exclusive Bus Lane Policy (EBLP), and High Occupancy Vehicle Lane Policy (HOVLP), incorporating optimized bus frequency and demand split under equilibrium condition. Switching thresholds for policy selection are derived to identify optimal periods for implementing each policy based on dynamic demand simulated using an Ornstein-Uhlenbeck (O-U) process. Results reveal significant reductions in total system costs with the proposed dynamic policy integration. Compared to static implementations, the combined policy achieves cost reductions of 12.0%, 5.3% and 42.5% relative to MTP-only, EBLP-only, and HOVLP-only scenarios, respectively. Additionally, in two real case studies of existing EBL and HOVL operations, the proposed dynamic policy reduces total costs by 32.2% and 27.9%, respectively. The findings provide valuable insights for policymakers and transit planners, offering a robust framework for dynamically scheduling and integrating EBL and HOVL policies to optimize urban corridor efficiency and reduce overall system costs.
Traffic safety is a critical concern in transportation engineering and urban planning. Traditional traffic safety analysis requires trained observers to collect data in the field, which is time-consuming, labor-intensive, and sometimes inaccurate. In recent years, microscopic traffic simulation, which simulates individual vehicles' movements within a transportation network, have been utilized to study traffic safety. However, microscopic traffic simulation only focuses on traffic-related factors, such as traffic volume, traffic signals, and lane configurations, neglecting vehicle dynamics and environment-related factors like weather and lighting conditions, which can significantly impact traffic safety. In light of this, this paper explores the application of digital twin technology in traffic safety analysis, integrating vehicle simulators, which consider vehicle dynamics and environmental factors, and microscopic traffic simulators, which simulate the operations of traffic flow, for enhanced safety evaluations. Various scenarios, including different weather conditions and visibility levels, are simulated using a digital twin of a road segment in Tuscaloosa, Alabama. The simulations employ Surrogate Safety Measures (SSMs) like Time to Collision (TTC) and Deceleration Rate to Avoid a Crash (DRAC) to assess safety under varying conditions. The results demonstrate that traffic digital twin can identify potential safety issues that traditional microscopic simulation cannot, providing insights for improving traffic control strategies and transportation infrastructure to enhance traffic safety.
A novel strategy aimed at cooperatively differentiating a signal among multiple interacting agents is introduced, where none of the agents needs to know which agent is the leader, i.e. the one producing the signal to be differentiated. Every agent communicates only a scalar variable to its neighbors; except for the leader, all agents execute the same algorithm. The proposed strategy can effectively obtain derivatives up to arbitrary $m$-th order in a finite time under the assumption that the $(m+1)$-th derivative is bounded. The strategy borrows some of its structure from the celebrated homogeneous robust exact differentiator by A. Levant, inheriting its exact differentiation capability and robustness to measurement noise. Hence, the proposed strategy can be said to perform robust exact distributed differentiation. In addition, and for the first time in the distributed leader-observer literature, sampled-data communication and bounded measurement noise are considered, and corresponding steady-state worst-case accuracy bounds are derived. The effectiveness of the proposed strategy is verified numerically for second- and fourth-order systems, i.e., for estimating derivatives of up to first and third order, respectively.
The Public Bus Transit Network Design Problem and Operations Planning (PBTNDP&OP) remains a core research area within transportation, in particular, because of the emergence of new transit technologies and services. Analytical approaches and mathematical programming are the most commonly applied methodologies to study this problem. Many studies utilize either of these two methods, often viewed as competing due to the unique benefits each provides that the other does not. This two-part paper systematically reviews the application of analytical approaches and mathematical programming in PBTNDP&OP, analyzing publications from 1968 to 2021. It begins by comparing analytical methods and mathematical programming through various statistical analyses, including the number of published papers, the most active journals and authors, keyword frequencies, keyword co-occurrence maps, and co-authorship maps. Subsequent analysis of the identified papers includes examinations from multiple perspectives: the problems investigated, modeling methods used, decision variables considered, network structures, and key findings. This is followed by a critical review of selected papers. The paper concludes by discussing the advantages and disadvantages of each approach and suggests potential extensions for future research based on identified gaps in existing studies.
This paper investigates the problem of certifying optimality for sparse generalized linear models (GLMs), where sparsity is enforced through an $\ell_0$ cardinality constraint. While branch-and-bound (BnB) frameworks can certify optimality by pruning nodes using dual bounds, existing methods for computing these bounds are either computationally intensive or exhibit slow convergence, limiting their scalability to large-scale problems. To address this challenge, we propose a first-order proximal gradient algorithm designed to solve the perspective relaxation of the problem within a BnB framework. Specifically, we formulate the relaxed problem as a composite optimization problem and demonstrate that the proximal operator of the non-smooth component can be computed exactly in log-linear time complexity, eliminating the need to solve a computationally expensive second-order cone program. Furthermore, we introduce a simple restart strategy that enhances convergence speed while maintaining low per-iteration complexity. Extensive experiments on synthetic and real-world datasets show that our approach significantly accelerates dual bound computations and is highly effective in providing optimality certificates for large-scale problems.
We introduce a computationally efficient method for the automation of inverse design in science and engineering. Based on simple least-square regression, the underlying dynamic mode decomposition algorithm can be used to construct a low-rank subspace spanning multiple experiments in parameter space. The proposed inverse design dynamic mode composition (ID-DMD) algorithm leverages the computed low-dimensional subspace to enable fast digital design and optimization on laptop-level computing, including the potential to prescribe the dynamics themselves. Moreover, the method is robust to noise, physically interpretable, and can provide uncertainty quantification metrics. The architecture can also efficiently scale to large-scale design problems using randomized algorithms in the ID-DMD. The simplicity of the method and its implementation are highly attractive in practice, and the ID-DMD has been demonstrated to be an order of magnitude more accurate than competing methods while simultaneously being 3-5 orders faster on challenging engineering design problems ranging from structural vibrations to fluid dynamics. Due to its speed, robustness, interpretability, and ease-of-use, ID-DMD in comparison with other leading machine learning methods represents a significant advancement in data-driven methods for inverse design and optimization, promising a paradigm shift in how to approach inverse design in practice.
Design of adversarial attacks for deep neural networks, as well as methods of adversarial training against them, are subject of intense research. In this paper, we propose methods to train against distributional attack threats, extending the TRADES method used for pointwise attacks. Our approach leverages recent contributions and relies on sensitivity analysis for Wasserstein distributionally robust optimization problems. We introduce an efficient fine-tuning method which can be deployed on a previously trained model. We test our methods on a range of pre-trained models on RobustBench. These experimental results demonstrate the additional training enhances Wasserstein distributional robustness, while maintaining original levels of pointwise robustness, even for already very successful networks. The improvements are less marked for models pre-trained using huge synthetic datasets of 20-100M images. However, remarkably, sometimes our methods are still able to improve their performance even when trained using only the original training dataset (50k images).
This paper presents a rigorous finite element framework for solving an optimal control problem governed by the steady Navier-Stokes-Brinkman equations, focusing on identifying a scalar permeability parameter $\gamma$ from local velocity observations. Three different finite element discretization schemes are proposed, and a priori error estimates are proven under appropriate regularity assumptions for each one. A key contribution of this paper is the development of residual-based a posteriori error estimators for both fully discrete and semi-discrete schemes, guiding adaptive mesh refinement to achieve comparable accuracy with fewer degrees of freedom. The method of manufactured solutions is used for numerical experiments to validate the theoretical findings, to demonstrate optimal convergence rates and the effectivity index is evaluated to measure their reliability. The framework offers insights into flow control mechanisms and paving the way for extensions to time-dependent, stochastic, or multiphysics problems.
We consider a random dynamical system on $\mathbb{R}^d$, whose dynamics is defined by a stochastic differential equation. The annealed transfer operator associated with such systems is a kernel operator. Given a set of feasible infinitesimal perturbations $P$ to this kernel, with support in a certain compact set, and a specified observable function $\phi: \mathbb{R}^d \to \mathbb{R}$, we study which infinitesimal perturbation in $P$ produces the greatest change in expectation of $\phi$. We establish conditions under which the optimal perturbation uniquely exists and present a numerical method to approximate the optimal infinitesimal kernel perturbation. Finally, we numerically illustrate our findings with concrete examples.
Recently, the fundamental lemma by Willems et. al has been extended towards stochastic LTI systems subject to process disturbances. Using this lemma requires previously recorded data of inputs, outputs, and disturbances. In this paper, we exploit causality concepts of stochastic control to propose a variant of the stochastic fundamental lemma that does not require past disturbance data in the Hankel matrices. Our developments rely on polynomial chaos expansions and on the knowledge of the disturbance distribution. Similar to our previous results, the proposed variant of the fundamental lemma allows to predict future input-output trajectories of stochastic LTI systems. We draw upon a numerical example to illustrate the proposed variant in data-driven control context.
The purpose of this note is to provide an optimal rate of convergence in the vanishing viscosity regime for first-order Hamilton-Jacobi equations with purely quadratic Hamiltonian. We show that for a globally Lipschitz-continuous terminal condition the rate is of order O($\epsilon$ log $\epsilon$), and we provide an example to show that this rate cannot be sharpened. This improves on the previously known rate of convergence O( $\sqrt$ $\epsilon$), which was widely believed to be optimal. Our proof combines techniques involving regularization by sup-convolution with entropy estimates for the flow of a suitable version of the adjoint linearized equation. The key technical point is an integrated estimate of the Laplacian of the solution against this flow. Moreover, we exploit the semiconcavity generated by the equation.