Source author record

Matthew Peet

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

ResearcherUnclaimed source record

math.OC math.NA Mathematical Software Numerical Analysis

Catalog footprint

What is connected

6works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2022arXiv

Computational stability analysis of PDEs with integral terms using the PIE framework

The Partial Integral Equation (PIE) framework was developed to computationally analyze linear Partial Differential Equations (PDEs) where the PDE is first converted to a PIE and then the analysis problem is solved by solving operator-valued optimization problems. Previous works on the PIE framework focused on the analysis of PDEs with spatial derivatives up to $2^{nd}$-order. In this paper, we extend the class of PDEs by including integral terms and performing stability analysis using the PIE framework. More specifically, we show that PDEs with the integral terms where the integration is with respect to the spatial variable and the kernel of the integral operator is matrix-valued polynomials can be converted to PIEs if the boundary conditions satisfy certain criteria. The conversion is performed by using a change of variable where every PDE state is substituted in terms of its highest derivative and boundary values to obtain a new equation (a PIE) in a variable that does not have any continuity requirements. Later, we show that this change of variable can be represented using explicit maps from the parameters of the PDE to the parameters of the PIE and the stability test can be posed as an optimization problem involving these parameters. Lastly, we present numerical examples to demonstrate the simplicity and application of this method.

preprint2022arXiv

Efficient Data Structures for Exploiting Sparsity and Structure in Representation of Polynomial Optimization Problems: Implementation in SOSTOOLS

We present a new data structure for representation of polynomial variables in the parsing of sum-of-squares (SOS) programs. In SOS programs, the variables $s(x;Q)$ are polynomial in the independent variables $x$, but linear in the decision variables $Q$. Current SOS parsers, however, fail to exploit the semi-linear structure of the polynomial variables, treating the decision variables as independent variables in their representation. This results in unnecessary overhead in storage and manipulation of the polynomial variables, prohibiting the parser from addressing larger-scale optimization problems. To eliminate this computational overhead, we introduce a new representation of polynomial variables, the "dpvar" structure, that is affine in the decision variables. We show that the complexity of operations on variables in the dpvar representation scales favorably with the number of decision variables. We further show that the required memory for storing polynomial variables is relatively small using the dpvar structure, particularly when exploiting the MATLAB sparse storage structure. Finally, we incorporate the dpvar data structure into SOSTOOLS 4.00, and test the performance of the parser for several polynomial optimization problems.

preprint2015arXiv

Polynomial Optimization with Applications to Stability Analysis and Control - Alternatives to Sum of Squares

In this paper, we explore the merits of various algorithms for polynomial optimization problems, focusing on alternatives to sum of squares programming. While we refer to advantages and disadvantages of Quantifier Elimination, Reformulation Linear Techniques, Blossoming and Groebner basis methods, our main focus is on algorithms defined by Polya's theorem, Bernstein's theorem and Handelman's theorem. We first formulate polynomial optimization problems as verifying the feasibility of semi-algebraic sets. Then, we discuss how Polya's algorithm, Bernstein's algorithm and Handelman's algorithm reduce the intractable problem of feasibility of semi-algebraic sets to linear and/or semi-definite programming. We apply these algorithms to different problems in robust stability analysis and stability of nonlinear dynamical systems. As one contribution of this paper, we apply Polya's algorithm to the problem of H_infinity control of systems with parametric uncertainty. Numerical examples are provided to compare the accuracy of these algorithms with other polynomial optimization algorithms in the literature.

preprint2014arXiv

Constructing Piecewise Polynomial Lyapunov Functions for Local Stability of Nonlinear Systems Using Handelman's Theorem

In this paper, we propose a new convex approach to stability analysis of nonlinear systems with polynomial vector fields. First, we consider an arbitrary convex polytope that contains the equilibrium in its interior. Then, we decompose the polytope into several convex sub-polytopes with a common vertex at the equilibrium. Then, by using Handelman's theorem, we derive a new set of affine feasibility conditions -solvable by linear programming- on each sub-polytope. Any solution to this feasibility problem yields a piecewise polynomial Lyapunov function on the entire polytope. This is the first result which utilizes Handelman's theorem and decomposition to construct piecewise polynomial Lyapunov functions on arbitrary polytopes. In a computational complexity analysis, we show that for large number of states and large degrees of the Lyapunov function, the complexity of the proposed feasibility problem is less than the complexity of certain semi-definite programs associated with alternative methods based on Sum-of-Squares or Polya's theorem. Using different types of convex polytopes, we assess the accuracy of the algorithm in estimating the region of attraction of the equilibrium point of the reverse-time Van Der Pol oscillator.

preprint2014arXiv

Using SOS for Analysis of Zeno Stability in Hybrid systems with Nonlinearity and Uncertainy

Hybrid systems exhibit phenomena which do not occur in systems with continuous vector fields. One such phenomenon - Zeno executions - is characterized by an infinite number of discrete events or transitions occurring over a finite interval of time. This phenomenon is not necessarily undesirable and may indeed be used to capture physical phenomena. In this paper, we examine the problem of proving the existence and stability of zero executions. Our approach is to develop a polynomial-time algorithm - based on the sum-of-squares methodology - for verifying the stability of a Zeno execution. We begin by stating Lyapunov-like theorems for local Zeno stability based on existing results. Then, for hybrid systems with polynomial vector fields, we use polynomial Lyapunov functions and semialgebraic geometry (Positivstellensatz results) to reduce the local Lyapunov-like conditions to a convex feasibility problem in polynomial variables. The feasibility problem is then tested using an algorithm for sum-of-squares programming - SOSTOOLS. We also extend these results to hybrid system with parametric uncertainty, where the uncertain parameters lie in a semialgebriac set. We also provide several examples illustrating the use of our technique.

preprint2013arXiv

A Sum-of-Squares Approach to the Analysis of Zeno Stability in Polynomial Hybrid Systems

Hybrid dynamical systems can exhibit many unique phenomena, such as Zeno behavior. Zeno behavior is the occurrence of infinite discrete transitions in finite time. Zeno behavior has been likened to a form of finite-time asymptotic stability, and corresponding Lyapunov theorems have been developed. In this paper, we propose a method to construct Lyapunov functions to prove Zeno stability of compact sets in cyclic hybrid systems with parametric uncertainties in the vector fields, domains and guard sets, and reset maps utilizing sum-of-squares programming. This technique can easily be applied to cyclic hybrid systems without parametric uncertainties as well. Examples illustrating the use of the proposed technique are also provided.