Researcher profile

Sachin Shivakumar

Sachin Shivakumar contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
6topics
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

7 published item(s)

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.

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.

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.

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.

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.

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.