---
title: 'PIETOOLS: MATLAB Toolbox for PIE Analysis'
url: https://www.emergentmind.com/topics/pietools
type: topic
---

# PIETOOLS: MATLAB Toolbox for PIE Analysis

PIETOOLS is a MATLAB-based software suite for the analysis, control, and robust optimization of infinite-dimensional systems represented as Partial Integral Equations (PIEs). Leveraging an algebra of Partial Integral (PI) operators, PIETOOLS generalizes Linear Matrix Inequality (LMI) methods to the operator level, enabling convex optimization over classes of partial differential equations (PDEs), delay-differential equations (DDEs), and distributed-parameter ODE–PDE/DDE–PDE systems. Its architecture integrates advanced operator-theoretic formulations, semidefinite programming, and symbolic representation, providing a tractable and numerically precise environment for stability analysis, input-output norm estimation, observer/controller synthesis, and robust/μ-analysis of complex dynamical systems [1910.01338][2501.17854][2004.03638][2511.14896].

## 1. Foundational Concepts: Partial Integral Equations and PI Operators

The core mathematical structure in PIETOOLS is the Partial Integral Equation (PIE) formalism, which encodes the dynamics of infinite-dimensional systems using PI operators that act on spaces such as $\mathbb{R}^m \times L_2^n[a,b]$ (or multi-dimensional analogues). A 4-PI operator $\mathcal{P}: \mathbb{R}^m \times L_2^n[a,b] \to \mathbb{R}^p \times L_2^q[a,b]$ is parameterized by a matrix $P$, two function-valued matrices $Q_1$, $Q_2$, and a 3-kernel operator $R$ [1910.01338]:

\[
\mathcal{P} \begin{bmatrix} x \\ y(\cdot) \end{bmatrix} =
\begin{bmatrix}
P\,x + \int_a^b Q_1(\theta)\,y(\theta)\, d\theta \\
Q_2(s) x + R_0(s)\, y(s) + \int_a^s R_1(s,\theta) y(\theta) d\theta + \int_s^b R_2(s,\theta) y(\theta) d\theta
\end{bmatrix}.
\]

These PI operators form a $*$-subalgebra under addition, composition, and adjoint [1910.01338]. In finite dimensions, the Lyapunov and control synthesis problems are posed as LMIs on matrices; in the PIE setting, these generalize to linear PI inequalities (LPIs) where the variable is a PI operator [2004.03638]. This abstraction is essential for developing convex stability and performance certificates for infinite-dimensional, spatially-distributed, and delay systems [1812.06794][2211.05326][1910.03881].

## 2. PIE Representation and System Conversion

PIETOOLS automates the conversion of a wide range of system classes into PIE representation. Any suitably well-posed linear ODE–PDE, DDE, or coupled system can be mapped to an algebraic PIE, where all boundary, continuity, and delay constraints are absorbed into the operator structure [2501.17854][2004.03638][2211.05326][2205.03735][1812.06794]. The “fundamental state”—typically the highest-order spatial derivative—forms an unconstrained state in an appropriate $L_2$ space. The conversion process is critical for transforming implicit constraint-laden models into boundary-free operator equations amenable to convex optimization.

The user specifies the PDE, boundary conditions, distributed delays, or ODE couplings in a symbolic MATLAB format. PIETOOLS applies analytic formulas, often relying on generalized Cauchy integrals or Gohberg–Kaashoek–Lay inversion, to compute the mapping operators (e.g., T, A, B, C, D) [2205.03735][2508.14840][2101.02050]. For multi-dimensional domains, PIETOOLS constructs multivariate PI operators using tensor products of univariate kernels, imposing explicit consistency conditions for domain intersections [2508.14840].

## 3. Convex Operator Inequalities: Stability and Performance Analysis

Central to PIETOOLS is the implementation of LPI-based certificates for system properties. To test exponential stability or input-output performance, PIETOOLS parameterizes self-adjoint positive PI operators (often via sum-of-squares or finite-dimensional basis expansions) and formulates operator inequalities that are direct analogues of LMIs [2004.03638][2508.14840][1812.06794][2203.15257]. For example, the Lyapunov stability for a PIE
\[
T \dot{x}(t) = A x(t)
\]
is certified by the feasibility of
\[
A^* P T + T^* P A \preceq -\delta\, T^* T, \quad P = P^* \succeq \varepsilon I,
\]
where $P$ is a PI operator [1812.06794][2501.17854]. For $H_\infty$/$L_2$-gain and impulse-to-peak norm bounds, block-operator KYP conditions are used [2004.03638][2604.03399].

PIETOOLS translates these infinite-dimensional inequalities into finite semidefinite programs (SDPs) by expanding operator kernels over polynomial or monomial bases and encoding positivity via matrix constraints [2203.15257][2604.03399]. Feasibility guarantees non-conservative, numerically certified stability and performance results for a wide variety of ODE–PDE, DDE, and multivariate PDE models.

## 4. Synthesis: Observer, Controller, and Robust Optimization

PIETOOLS supports the synthesis of optimal state feedback controllers, Luenberger observers, and robust controllers under uncertainty. Synthesis problems are posed as convex LPIs over PI operator variables. For $H_\infty$-optimal state-feedback and observer design, the duality theory for PIEs implies the convexification variable change (e.g., $K = Z P^{-1}$ for controllers, $L = P^{-1} Z$ for observers), mirroring finite-dimensional results [2004.03638][2004.04482][2402.05061][2501.17854].

Robust and $\mu$-optimal synthesis in the presence of structured uncertainty uses Integral Quadratic Constraint (IQC) multipliers, also modeled as PI operators [2511.14896]. The design procedure involves joint optimization over Lyapunov PI operators and IQC multipliers, yielding upper bounds on $\mu$ or worst-case gain with significantly reduced conservatism compared to unstructured analysis.

For nonlinear or polynomial PDE dynamics, PIETOOLS interfaces with distributed sum-of-squares (SOS) representations in tensor-product PI operator algebras, enabling local stability analysis via polynomial kernel expansions [2604.01115].

## 5. Software Architecture, User Workflow, and Example Applications

PIETOOLS is structured as a MATLAB toolbox, introducing the `opvar` (operator variable) and `pio` classes for symbolic and numerical manipulation of PI operators [1910.01338]. Operator arithmetic—addition, composition, adjoint—is overloaded. The toolbox integrates with SOSTOOLS/MULTIPOLY for polynomial manipulations and interfaces with SDP solvers (e.g., SeDuMi, MOSEK, SDPT3) for LPI feasibility [1910.01338][2501.17854][2101.02050].

A typical user workflow [2501.17854][2101.02050]:
1. Define spatial, temporal, and auxiliary variables (e.g., `pvar s, t, theta`).
2. Specify system equations, boundary conditions, and inputs/outputs through symbolic PDE or DDE models.
3. Convert to PIE form using conversion routines (e.g., `convert(PDE,'pie')`).
4. Declare PI-operator decision variables and form LPI constraints via built-in functions (`sosprogram`, `sos_posopvar`, `sos_opineq`).
5. Solve the resulting SDP and extract operator-valued solutions (`sossolve`, `sosgetsol_opvar`).
6. Construct synthesized controllers or observers and simulate open- or closed-loop PIE systems using `PIESIM`.

PIETOOLS ships with end-to-end demo scripts for canonical PDEs (e.g., tip-damped wave equation, 2D heat and plate, reaction–diffusion with delay, boundary ODE–PDE couplings), and supports batch benchmarking for scalability studies [1812.06794][2604.03399][2203.15257].

## 6. Performance, Scalability, and Comparison

Numerical results indicate that PIETOOLS achieves high accuracy—matching analytic or high-fidelity finite-difference results to four significant digits in $H_\infty$ and impulse-peak norms [2203.15257][2604.03399][2004.04482]—and can handle problems with up to tens of PDE states, operator degrees up to four, and moderate polynomial basis sizes in seconds to a few minutes on contemporary desktops [1812.06794][2508.14840][2501.17854]. The direct operator-theoretic approach eliminates the need for spatial discretization or ad hoc state augmentation.

Compared to ad hoc SOS techniques, PIETOOLS significantly extends scalability, handles higher-order and multi-dimensional PDEs, supports robust IQC/μ-analysis, and provides closed-form controller/observer recovery via analytic inversion of PI operators (in particular, for structured PI operators in the separable or rank-one-integral class) [2402.05061][2511.14896][2205.03735]. For delay-coupled PDEs or spatially distributed uncertainties, the reduction in conservatism relative to unstructured $L_2$-gain methods can reach 30–40% [2511.14896].

## 7. Limitations, Extensions, and Future Directions

PIETOOLS, as of its latest documented versions, is restricted to linear, time-invariant systems with polynomial or $L_\infty$ coefficients, and to positive operators parameterized via polynomial bases [2501.17854][2211.05326][2205.03735]. Extensions under development address nonlinear PDEs (via SOS-in-tensor-algebras), periodic boundary conditions (with quotient-space representations), time-varying and switching PIE systems, and multivariate domains beyond hyper-rectangles [2604.01115][2508.14840][2211.05326].

Theoretical limitations include the need for prior model well-posedness checks (the software does not certify Lumer–Phillips conditions or boundary-input regularity), and computational scaling in the polynomial degree and state dimension. Nonlinear generalizations, trajectory-level performance analysis (e.g., input–state–output $L_2$–$L_\infty$ norms), and scalable policy-synthesis for semi-infinite constraints remain active research topics.

---

### References

- [PIETOOLS: A Matlab Toolbox for Manipulation and Optimization of Partial Integral Operators](https://arxiv.org/abs/1910.01338) [1910.01338].
- [PIETOOLS 2024: User Manual](https://arxiv.org/abs/2501.17854) [2501.17854].
- [Duality and $H_{\infty}$-Optimal Control Of Coupled ODE-PDE Systems](https://arxiv.org/abs/2004.03638) [2004.03638].
- [A Partial Integral Equation (PIE) Representation of Coupled Linear PDEs and Scalable Stability Analysis using LMIs](https://arxiv.org/abs/1812.06794) [1812.06794].
- [H\infty-Optimal Observer Design for Linear Systems with Delays in States, Outputs and Disturbances](https://arxiv.org/abs/2004.04482) [2004.04482].
- [Representation of PDE Systems with Delay and Stability Analysis using Convex Optimization](https://arxiv.org/abs/2211.05326) [2211.05326].
- [A State-Space Representation of Coupled Linear Multivariate PDEs and Stability Analysis using SDP](https://arxiv.org/abs/2508.14840) [2508.14840].
- [Extension of the Partial Integral Equation Representation to GPDE Input-Output Systems](https://arxiv.org/abs/2205.03735) [2205.03735].
- [$L_2$-Gain Analysis of Coupled Linear 2D PDEs using Linear PI Inequalities](https://arxiv.org/abs/2203.15257) [2203.15257].
- [A $μ$-Analysis and Synthesis Framework for Partial Integral Equations using IQCs](https://arxiv.org/abs/2511.14896) [2511.14896].
- [PIETOOLS 2022: User Manual](https://arxiv.org/abs/2101.02050) [2101.02050].
- [Impulse-to-Peak-Output Norm Optimal State-Feedback Control of Linear PDEs](https://arxiv.org/abs/2604.03399) [2604.03399].
- [A Distributed SOS Program For Local Stability Analysis of Polynomial PDEs in the PIE Representation](https://arxiv.org/abs/2604.01115) [2604.01115].
- [Representation of Networks and Systems with Delay: DDEs, DDFs, ODE-PDEs and PIEs](https://arxiv.org/abs/1910.03881) [1910.03881].
- [$H_{\infty}$-Optimal Estimator Synthesis for Coupled Linear 2D PDEs using Convex Optimization](https://arxiv.org/abs/2402.05061) [2402.05061].

Source: https://www.emergentmind.com/topics/pietools