BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Matthew M. Peet

Matthew M. Peet contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.OCmath.DSMachine Learningmath.AP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

preprint2022arXiv

Converse Lyapunov Functions and Converging Inner Approximations to Maximal Regions of Attraction of Nonlinear Systems

This paper considers the problem of approximating the "maximal" region of attraction (the set that contains all asymptotically stable sets) of any given set of locally exponentially stable nonlinear Ordinary Differential Equations (ODEs) with a sufficiently smooth vector field. Given a locally exponential stable ODE with a differentiable vector field, we show that there exists a globally Lipschitz continuous converse Lyapunov function whose 1-sublevel set is equal to the maximal region of attraction of the ODE. We then propose a sequence of d-degree Sum-of-Squares (SOS) programming problems that yields a sequence of polynomials that converges to our proposed converse Lyapunov function uniformly from above in the L1 norm. We show that each member of the sequence of 1-sublevel sets of the polynomial solutions to our proposed sequence of SOS programming problems are certifiably contained inside the maximal region of attraction of the ODE, and moreover, we show that this sequence of sublevel sets converges to the maximal region of attraction of the ODE with respect to the volume metric. We provide numerical examples of estimations of the maximal region of attraction for the Van der Pol oscillator and a three dimensional servomechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Convex Parametrization of a New Class of Universal Kernel Functions

The accuracy and complexity of kernel learning algorithms is determined by the set of kernels over which it is able to optimize. An ideal set of kernels should: admit a linear parameterization (tractability); be dense in the set of all kernels (accuracy); and every member should be universal so that the hypothesis space is infinite-dimensional (scalability). Currently, there is no class of kernel that meets all three criteria - e.g. Gaussians are not tractable or accurate; polynomials are not scalable. We propose a new class that meet all three criteria - the Tessellated Kernel (TK) class. Specifically, the TK class: admits a linear parameterization using positive matrices; is dense in all kernels; and every element in the class is universal. This implies that the use of TK kernels for learning the kernel can obviate the need for selecting candidate kernels in algorithms such as SimpleMKL and parameters such as the bandwidth. Numerical testing on soft margin Support Vector Machine (SVM) problems show that algorithms using TK kernels outperform other kernel learning algorithms and neural networks. Furthermore, our results show that when the ratio of the number of training data to features is high, the improvement of TK over MKL increases significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Generalized LMI Formulation for Input-Output Analysis of Linear Systems of ODEs Coupled with PDEs

In this paper, we consider input-output properties of linear systems consisting of PDEs on a finite domain coupled with ODEs through the boundary conditions of the PDE. This framework can be used to represent e.g. a lumped mass fixed to a beam or a system with delay. This work generalizes the sufficiency proof of the KYP Lemma for ODEs to coupled ODE-PDE systems using a recently developed concept of fundamental state and the associated boundary-condition-free representation. The conditions of the generalized KYP are tested using the PQRS positive matrix parameterization of operators resulting in a finite-dimensional LMI, feasibility of which implies prima facie provable passivity or L2-gain of the system. No discretization or approximation is involved at any step and we use numerical examples to demonstrate that the bounds obtained are not conservative in any significant sense and that computational complexity is lower than existing methods involving finite-dimensional projection of PDEs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A New Algorithm for Tessellated Kernel Learning

The accuracy and complexity of machine learning algorithms based on kernel optimization are limited by the set of kernels over which they are able to optimize. An ideal set of kernels should: admit a linear parameterization (for tractability); be dense in the set of all kernels (for robustness); be universal (for accuracy). The recently proposed Tesselated Kernels (TKs) is currently the only known class which meets all three criteria. However, previous algorithms for optimizing TKs were limited to classification and relied on Semidefinite Programming (SDP) - limiting them to relatively small datasets. By contrast, the 2-step algorithm proposed here scales to 10,000 data points and extends to the regression problem. Furthermore, when applied to benchmark data, the algorithm demonstrates significant improvement in performance over Neural Nets and SimpleMKL with similar computation time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Partial Integral Equation (PIE) Representation of Coupled Linear PDEs and Scalable Stability Analysis using LMIs

We present a new Partial Integral Equation (PIE) representation of Partial Differential Equations (PDEs) in which it is possible to use convex optimization to perform stability analysis with little or no conservatism. The first result gives a standardized representation for coupled linear PDEs in a single spatial variable and shows that any such PDE, suitably well-posed, admits an equivalent PIE representation, defined by the given conversion formulae. This leads to a new prima facie representation of the dynamics without the implicit constraints on system state imposed by boundary conditions. The second result is to show that for systems in this PIE representation, convex optimization may be used to verify stability without discretization. The resulting algorithms are implemented in the Matlab toolbox PIETOOLS, tested on several illustrative examples, compared with previous results, and the code has been posted on Code Ocean. Scalability testing indicates the algorithm can analyze systems of up to 40 coupled PDEs on a desktop computer.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Computing Input-Output Properties of Coupled PDE systems

In this paper, we propose an LMI-based approach to analyze input-output properties of coupled linear PDE systems. This work expands on a newly developed state-space theory for coupled PDEs and extends the positive-real and bounded-real lemmas to infinite dimensional systems. We show that conditions for passivity and bounded L2 gain can be expressed as linear operator inequalities on RxL2. A method to convert these operator inequalities to LMIs by using parameterization of the operator variables is proposed. This method does not rely on discretization and as such, the properties obtained are prima facie provable. We use numerical examples to demonstrate that the bounds obtained are not conservative in any significant sense and that the bounds are computable on desktop computers for systems consisting of up to 20 coupled PDEs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Duality and $H_{\infty}$-Optimal Control Of Coupled ODE-PDE Systems

In this paper, we present a convex formulation of $H_{\infty}$-optimal control problem for coupled linear ODE-PDE systems with one spatial dimension. First, we reformulate the coupled ODE-PDE system as a Partial Integral Equation (PIE) system and show that stability and $H_{\infty}$ performance of the PIE system implies that of the ODE-PDE system. We then construct a dual PIE system and show that asymptotic stability and $H_{\infty}$ performance of the dual system is equivalent to that of the primal PIE system. Next, we pose a convex dual formulation of the stability and $H_{\infty}$-performance problems using the Linear PI Inequality (LPI) framework. LPIs are a generalization of LMIs to Partial Integral (PI) operators and can be solved using PIETOOLS, a MATLAB toolbox. Next, we use our duality results to formulate the stabilization and $H_{\infty}$-optimal state-feedback control problems as LPIs. Finally, we illustrate the accuracy and scalability of the algorithms by constructing controllers for several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Extensions of the Dynamic Programming Framework: Battery Scheduling, Demand Charges, and Renewable Integration

We consider a general class of Dynamic Programming (DP) problems with non-separable objective functions. We show that for any problem in this class, there exists an augmented-state DP problem which satisfies the Principle of Optimality and the solutions to which yield solutions to the original problem. Furthermore, we identify a subclass of DP problems with Naturally Forward Separable (NFS) objective functions for which this state-augmentation scheme is tractable. We extend this framework to stochastic DP problems, proposing a suitable definition of the Principle of Optimality. We then apply the resulting algorithms to the problem of optimal battery scheduling with demand charges using a data-based stochastic model for electricity usage and solar generation by the consumer.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

H\infty-Optimal Observer Design for Linear Systems with Delays in States, Outputs and Disturbances

This paper considers the H\infty-optimal estimation problem for linear systems with multiple delays in states, output, and disturbances. First, we formulate the H\infty-optimal estimation problem in the Delay-Differential Equation (DDE) framework. Next, we construct an equivalent Partial Integral Equation (PIE) representation of the optimal estimator design framework. We then show that in the PIE framework, the H\infty-optimal estimator synthesis problem can be posed as a Linear PI Inequality (LPI). LPIs are a generalization of LMIs to the algebra of Partial Integral (PI) operators and can be solved using the PIETOOLS toolbox. Finally, we convert the PIE representation of the optimal estimator back into an ODE-PDE representation - a form similar to a DDE, but with corrections to estimates of the infinite-dimensional state (the time-history). Numerical examples show that the synthesis condition we propose produces an estimator with provable H\infty-gain bound which is accurate to 4 decimal places when compared with results obtained using Pade-based discretization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

PIETOOLS: A Matlab Toolbox for Manipulation and Optimization of Partial Integral Operators

In this paper, we present PIETOOLS, a MATLAB toolbox for the construction and handling of Partial Integral (PI) operators. The toolbox introduces a new class of MATLAB object, opvar, for which standard MATLAB matrix operation syntax (e.g. +, *, ' e tc.) is defined. PI operators are a generalization of bounded linear operators on infinite-dimensional spaces that form a *-subalgebra with two binary operations (addition and composition) on the space RxL2. These operators frequently appear in analysis and control of infinite-dimensional systems such as Partial Differential equations (PDE) and Time-delay systems (TDS). Furthermore, PIETOOLS can: declare opvar decision variables, add operator positivity constraints, declare an objective function, and solve the resulting optimization problem using a syntax similar to the sdpvar class in YALMIP. Use of the resulting Linear Operator Inequalities (LOIs) are demonstrated on several examples, including stability analysis of a PDE, bounding operator norms, and verifying integral inequalities. The result is that PIETOOLS, packaged with SOSTOOLS and MULTIPOLY, offers a scalable, user-friendly and computationally efficient toolbox for parsing, performing algebraic operations, setting up and solving convex optimization problems on PI operators.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Relaxing The Hamilton Jacobi Bellman Equation To Construct Inner And Outer Bounds On Reachable Sets

We consider the problem of overbounding and underbounding both the backward and forward reachable set for a given polynomial vector field, nonlinear in both state and input, with a given semialgebriac set of initial conditions and with inputs constrained pointwise to lie in a semialgebraic set. Specifically, we represent the forward reachable set using the value function which gives the optimal cost to go of an optimal control problems and if smooth satisfies the Hamilton-Jacobi- Bellman PDE. We then show that there exist polynomial upper and lower bounds to this value function and furthermore, these polynomial sub-value and super-value functions provide provable upper and lower bounds to the forward reachable set. Finally, by minimizing the distance between these sub-value and super-value functions in the L1-norm, we are able to construct inner and outer bounds for the reachable set and show numerically on several examples that for relatively small degree, the Hausdorff distance between these bounds is negligible.

0 citations0 reviewsNo code linkedOpen access