Abstract
Turnpike phenomena of nonlinear port-Hamiltonian descriptor systems under minimal energy supply are studied. Under assumptions on the smoothness of the system nonlinearities, it is shown that the optimal control problem is dissipative with respect to a manifold. Then, under controllability assumptions, it is shown that the optimal control problem exhibits a manifold turnpike property.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
This paper is concerned with turnpike phenomena. These phenomena were first noticed in the context of economics [1, 2] and have since been observed in many different situations, see, e.g., [3,4,5] and the references therein. Turnpike phenomena are studied in optimal control problems, where the goal is the minimization of a cost functional C(u). Here, the function u acts as the control of a system of interest. In many cases, it can be observed that an optimal control \(u^*\) will, for a majority of the time horizon, steer the associated state trajectory \(x^*\) to a point [6,7,8], a set [4, 5, 9] or, as in our case, a manifold [10]. In other words, optimal solutions depend mainly on the underlying system and the optimization objective and are more or less independent on the choice of the time horizon and other data, such as initial or final values. The behaviour is reminiscent of an observation from daily life: when traveling a long distance by car, it is usually faster to take a detour via a turnpike than to choose a more direct way on slower streets. Also, the chosen turnpike usually does not heavily depend on the start and end of the route. If one would start the journey a few blocks away, then the fastest path would remain more or less the same.
Here, we consider a special class of systems called port-Hamiltonian (pH) systems. Parts of the origins of port-Hamiltonian systems date back to the late 1950s [11], and the interested reader is referred to [12] for an overview on the origins of this system class. Despite their long history, they continue to be the focus of active research [13,14,15,16,17,18]. Arguably, the key feature of pH systems is their modeling perspective: they focus on taking energy as the lingua franca between subsystems. As a consequence, the class of pH systems is a promising class for modeling real-world processes [18]. Benefits of port-Hamiltonian models include inherent stability and passivity, the invariance under Galerkin projection and congruence transformation, and the possibility to interconnect multiple pH systems in a structure-preserving manner.
When pH systems are considered in an optimal control setting, the objective of minimizing the supplied energy is quite natural. This results in a cost term of the form \(\int _0^T y(t)^{{\textsf{T}}}u(t) ~\textrm{d}t\), where y is a collocated observation of the system, and renders the corresponding optimal control problem singular. In [19, 20], the authors have considered this objective for linear time invariant port-Hamiltonian (descriptor) systems. They have shown that the optimal control problem has a measure turnpike property with respect to the part of the state space where dissipation is present, given by the kernel of the matrix corresponding to the non-conservative system dynamics. The infinite-dimensional linear case is discussed in [21]. In this paper, we are concerned with the finite-dimensional nonlinear case. We show that, under smoothness assumptions on the nonlinearities and controllability assumptions on the system, nonlinear pH descriptor systems admit a turnpike phenomenon with respect to a submanifold of \({\mathbb {R}}^n\). This submanifold corresponds, as in the linear case, to the energy dissipating part of the state space.
The structure of this paper is as follows. In Sect. 2, we recall the definition of port-Hamiltonian systems and state the optimal control problem that is considered. After that, a short repetition of results on submanifolds of \({\mathbb {R}}^n\) follows in Sect. 3. In Sect. 4, we define manifold dissipativity and manifold turnpikes following [10] and recall that, under weak assumptions, manifold dissipativity implies a manifold turnpike property. Section 5 contains the main results of this work, where the previously established results are applied to finite-dimensional nonlinear port-Hamiltonian descriptor systems. The theoretical results are then illustrated by a numerical example in Sect. 6. Finally, in Sect. 7 a conclusion is drawn and an outlook on future research is given.
The paper is supplemented by an appendix which collects existence results for the appearing differential algebraic equations and optimization problems.
1.1 Notation
For a set \(Z \subseteq {\mathbb {R}}^n\) we define \(Z^\circ \) as the interior of Z. We denote the Euclidean norm by \(\Vert \cdot \Vert \), other p-norms by \(\Vert \cdot \Vert _p\), and define the distance of a point \(x\in {\mathbb {R}}^n\) to the set \({\mathcal {M}}\subseteq {\mathbb {R}}^n\) as \(\mathop {\textrm{dist}}\limits (x,{\mathcal {M}}):= \inf _{p \in {\mathcal {M}}} \Vert x-p \Vert \). We denote the set of all k-times continuously differentiable functions from U to V by \(C^k(U,V)\) and define \(C(U,V):= C^0(U,V)\). When the spaces U and V are clear from context, we say that \(f\in C^k(U,V)\) is of class \(C^k\). The derivative of a function f at point x is denoted by \(Df_x\). Further, for a matrix \(A\in {\mathbb {R}}^{n,n}\) we write \(A \succeq 0\) if \(x^{{\textsf{T}}}A x \ge 0\) for all \(x\in {\mathbb {R}}^{n}\), and \(A \succ 0\) if \(x^{{\textsf{T}}}A x > 0\) for all \(x\in {\mathbb {R}}^{n}{\setminus }\{0\}\), where \(\cdot ^{{\textsf{T}}}\) denotes the transpose. The kernel and range of the matrix A are denoted by \(\ker (A)\) and \(\mathop {\textrm{ran}}\limits (A)\), respectively. The non-negative square-root of a positive semidefinite matrix A is denoted by \(A^{1/2}\). Often, we suppress the time dependency of functions and write z instead of z(t).
2 Preliminaries and problem setting
Following the definition of [18], we consider nonlinear port-Hamiltonian descriptor systems of the form
Here, x, u and y are the state, input and output of the system, respectively. We restrict our analysis to pH systems without feedthrough but note that the discussion can be extended to systems with a feedthrough term as introduced in [18]. We consider the state space \({\mathbb {R}}^n\) and a set of admissible controls \({\mathcal {U}}_{\textrm{ad}}\) and require that
Further, the functions J and R have to satisfy \(J(x) = - J(x)^{{\textsf{T}}}\) and \(R(x) = R(x)^{{\textsf{T}}}\succeq 0\) for all \(x \in {\mathbb {R}}^n\). We assume that the system (pH) is associated with a Hamiltonian \({\mathcal {H}} \in C^1({\mathbb {R}}^n,{\mathbb {R}})\) which is bounded from below along any solution of (pH) and satisfies
for each \(x\in {\mathbb {R}}^n\). As in [18], we may assume that \({\mathcal {H}}\) is nonnegative along any solution of (pH), since otherwise we can add its infimum along solutions of (pH). A brief discussion regarding the existence of solutions of (pH) is given in the appendix.
Remark 1
We only consider pH systems which do not explicitly depend on time. The definition given here can be generalized to include explicit time dependence, but these systems can easily be made autonomous [18, Remark 4.2].
When it comes to the optimal control of port-Hamiltonian systems, the cost functional should take into account that pH systems stem from using energy as the lingua franca. Hence, choosing the supplied energy as the optimization objective is quite natural. For this, usually the impedance supply \(y^{{\textsf{T}}}u\) is considered, which is related to the scattering supply \(\Vert u \Vert ^2 - \Vert y \Vert ^2\) via the Cayley transform [22]. We focus on the former and thus consider the optimal control problem
Here and in the following, we assume \(x_0, x_T \in K\), where \(K \subseteq {\mathbb {R}}^n\) is a compact set. Further, as we are interested in the properties of optimal solutions to (\(\text {pH} \text { OCP}_T\)), throughout the paper we assume that an optimal solution \(u^*\) and a corresponding trajectory \(x^*\) exist. The existence of optimal controls is shortly discussed in the appendix.
It can be shown [18] that the power-balance equation
holds along any solution x of (pH). This allows us to rewrite the cost functional \(C_{\text {pH},T}(u)\) as
This equation is called the energy-balance equation, and we can interpret each of the terms physically [23]. The term \({\mathcal {H}}(x_T) - {\mathcal {H}}(x_0)\) measures the conserved energy, while the integral term corresponds to the dissipated energy. By rearranging and plugging in the definition of \(C_{\text {pH},T}\), we see that
Note that this implies dissipativity in the sense of Willems [24], and as a consequence shows the aforementioned passivity of pH systems. We will use both of these equations in Sect. 5.
3 Submanifolds of \({\mathbb {R}}^n\) and the orthogonal projection
This section repeats mostly well-known results regarding submanifolds of \({\mathbb {R}}^n\), with a focus on manifolds defined as the zero locus of some smooth function. The main result of the section is Lemma 5, which provides an upper bound for the distance of a point to such a manifold.
We begin with recalling the classical definition of submanifolds of \({\mathbb {R}}^n\). Following [25, 26], we distinguish manifolds whose tangent spaces locally satisfy a Lipschitz condition.
Definition 2
(submanifolds of \({\mathbb {R}}^n\))
-
Let \({\mathcal {M}}\) be a subset of \({\mathbb {R}}^n\). We call \({\mathcal {M}}\) an s-dimensional \(C^k\) manifold if for each \(p\in {\mathcal {M}}\) there exists an open neighborhood U of p and a \(C^k\) diffeomorphism \(\phi :U \rightarrow \phi (U)\subseteq {\mathbb {R}}^n\) such that
$$\begin{aligned} {\mathcal {M}}\cap U = \{ x\in U ~|~ \phi _{s+1}(x) = \dots = \phi _{n}(x) =0 \}. \end{aligned}$$The function \(\phi \) is called a local coordinate system of \({\mathcal {M}}\) at p.
-
Let \({\mathcal {M}}\subseteq {\mathbb {R}}^n\) be an s-dimensional manifold, \(p\in {\mathcal {M}}\) and let \(\phi :U \rightarrow {\mathbb {R}}^n\) be a local coordinate system of \({\mathcal {M}}\) at p. We define the tangent space at p relative to \({\mathcal {M}}\) as
$$\begin{aligned} T_p{\mathcal {M}}:= D\phi _{\phi (p)}^{-1}(\left\{ y\in {\mathbb {R}}^n ~\big |~ y_{s+1} = \dots = y_n = 0 \right\} ). \end{aligned}$$The space \(N_p{\mathcal {M}}:= T_p{\mathcal {M}}^\perp \) is called the normal space at p relative to \({\mathcal {M}}\).
-
We call \({\mathcal {M}}\subseteq {\mathbb {R}}^n\) an s-dimensional \(C^{k,1}\) manifold if \({\mathcal {M}}\) is an s-dimensional \(C^k\) manifold and for all \(p \in {\mathcal {M}}\) there exists a set \(V\subseteq {\mathcal {M}}\) that is open relative to \({\mathcal {M}}\) and a positive constant \(L>0\) such that \(p \in V\) and for all \({{\widetilde{p}}} \in V\) it holds that
$$\begin{aligned} d_{\textrm{H}}(T_p{\mathcal {M}}, T_{{\widetilde{p}}}{\mathcal {M}}) \le L \Vert p - {{\widetilde{p}}} \Vert . \end{aligned}$$Here, \(d_{\textrm{H}}\) denotes the Hausdorff distance defined by
$$\begin{aligned} d_{\textrm{H}}(T_1, T_2):= \sup \left\{ \inf \left\{ \Vert t_2-t_1 \Vert ~\big |~ t_2 \in T_2 \cap S \right\} ~\big |~ t_1 \in T_1 \cap S \right\} , \end{aligned}$$where \(S:= \left\{ z \in {\mathbb {R}}^n ~\big |~ \Vert z \Vert = 1 \right\} \) is the unit sphere.
Next, we recall the definition of the orthogonal projection on a manifold from [26]. For this, consider a manifold \({\mathcal {M}}\subseteq {\mathbb {R}}^n\) and define the set of points with the unique nearest points property as
Clearly, for each \(x\in \mathop {\textrm{unpp}}\limits ({\mathcal {M}})\) there exists a unique \(p(x)\in {\mathcal {M}}\) with the property
Definition 3
(orthogonal projection on manifold, [26]) Let \({\mathcal {M}}\subseteq {\mathbb {R}}^n\) be a manifold. The function \(p:\mathop {\textrm{unpp}}\limits ({\mathcal {M}}) \rightarrow {\mathcal {M}}, ~ x \mapsto p(x)\) is called the orthogonal projection on \({\mathcal {M}}\).
Often, we write \(p_x\) instead of p(x). The maximal open set on which the orthogonal projection is defined plays a special role in [26] and also in our setting. We will refer to this set as \({\mathcal {E}}({\mathcal {M}}):= \mathop {\textrm{unpp}}\limits ({\mathcal {M}})^\circ \).
The next proposition collects selected results on submanifolds of \({\mathbb {R}}^n\) defined in a particular manner and will be useful in Sect. 5, where it will allow us to study the optimal control problem (\(\text {pH} \text { OCP}_T\)).
Proposition 4
Suppose \(f:{\mathbb {R}}^n \rightarrow {\mathbb {R}}^n\) is of class \(C^2\), assume that
is nonempty and assume that there exists an open neighborhood \(G\subseteq {\mathbb {R}}^n\) of \({\mathcal {M}}\) and a constant s with \(0< s < n\) such that \(\dim (\ker (Df_x)) = s\) for all \(x \in G\). Then
-
(i)
the set \({\mathcal {M}}\) is an s-dimensional \(C^2\) submanifold of \({\mathbb {R}}^n\),
-
(ii)
the tangent space at \(p \in {\mathcal {M}}\) is given by \(T_p{\mathcal {M}}= \ker (Df_p)\),
-
(iii)
the manifold \({\mathcal {M}}\) is a \(C^{2,1}\) manifold, and
-
(iv)
it holds that \({\mathcal {M}}\subseteq {\mathcal {E}}({\mathcal {M}})\), and if \(x\in {\mathcal {E}}({\mathcal {M}})\setminus {\mathcal {M}}\) then \(x-p_x \perp T_{p_x}{\mathcal {M}}\).
Proof
-
(i)
Let \(p\in {\mathcal {M}}\). It can be shown [27] that there exist open neighborhoods \(U_1\subseteq {\mathbb {R}}^n\) of p and \(U_2\subseteq {\mathbb {R}}^n\) of f(p) and \(C^2\) diffeomorphisms
$$\begin{aligned} \phi :U_1 \rightarrow \phi (U_1)~~~~ \text {and} ~~~~ \psi :U_2 \rightarrow \psi (U_2) \end{aligned}$$such that
$$\begin{aligned} f(U_1) \subseteq U_2 \end{aligned}$$and
$$\begin{aligned} \psi \circ f \circ \phi ^{-1}(y_{1},\dots ,y_n) = (y_1,\dots ,y_{n-s}, 0,\dots ,0) \end{aligned}$$for all \(y\in \phi (U_1)\). Thus, for each \(v \in U_1\) we have
$$\begin{aligned} v \in {\mathcal {M}}&\Longleftrightarrow f(v) = 0 \\&\Longleftrightarrow (\psi \circ f)(v) = 0 \\&\Longleftrightarrow (\psi \circ f \circ \phi ^{-1}\circ \phi )(v) = 0 \\&\Longleftrightarrow \phi _1(v) = \dots = \phi _{n-s}(v) = 0. \end{aligned}$$ -
(ii)
Suppose that \(p\in {\mathcal {M}}\) and that \(\phi :U \rightarrow {\mathbb {R}}^n\) is a local coordinate system of \({\mathcal {M}}\) at p. Since
$$\begin{aligned} (f\circ \phi ^{-1}) \big (\{ y\in {\mathbb {R}}^n ~|~ y_{s+1} = \dots = y_n = 0 \} \cap \phi (U)\big ) = \{0\}, \end{aligned}$$we have
$$\begin{aligned} Df_p(T_p{\mathcal {M}}) = Df_p\Big (D\phi _{\phi (p)}^{-1}\big (\left\{ y\in {\mathbb {R}}^n ~\big |~ y_{s+1} = \dots = y_n = 0 \right\} \big )\Big ) = 0. \end{aligned}$$The claim then follows from the fact that \(\ker (Df_p)\) has dimension s.
-
(iii)
In [25, Equations (3.3) and (3.6)], it is shown that a \(C^k\) manifold with \(k\ge 2\) is also a \(C^{k,1}\) manifold, from which the claim follows.
-
(iv)
This claim was mentioned in [26, Remark 4] and is a consequence of the results of [25].
\(\square \)
We finish this section with Lemma 5 and a corresponding remark. The lemma establishes an upper bound on the distance to the manifold \({\mathcal {M}}\) defined in Proposition 4 in terms of the function f. This result will be the key in our application to port-Hamiltonian systems, as it will allow us to deduce a dissipativity property for (\(\text {pH} \text { OCP}_T\)).
Lemma 5
Suppose \(f:{\mathbb {R}}^n \rightarrow {\mathbb {R}}^n\) is of class \(C^2\), assume that
is nonempty and assume that there exists an open neighborhood \(G\subseteq {\mathbb {R}}^n\) of \({\mathcal {M}}\) such that for all \(x\in G\) it holds that \(\dim (\ker (Df_x)) = s\), where the constant s is independent of x and satisfies \(0< s < n\). Further, assume that for each \(x\in G\) the smallest nonzero singular value of \(Df_x\) is bounded from below by \({\widetilde{c}} > 0\). Then \({\mathcal {M}}\) is a \(C^{2}\) manifold and there exists an open set \(V \subseteq {\mathbb {R}}^n\) and a constant \(c>0\) with \({\mathcal {M}}\subseteq V \subseteq {\mathcal {E}}({\mathcal {M}})\) and
for all \(x\in V\).
Proof
We will first show that (2) is true locally. Fix a point \(p\in {\mathcal {M}}\) and notice that due to \(f\in C^2\), for all \(x \in {\mathbb {R}}^n\) we have
where the remainder \(g_p\) satisfies
To establish (2) locally, our first goal is to show that there exists an open set \(U_p \subseteq {\mathcal {E}}({\mathcal {M}})\) such that for all \(x\in V_p:= U_p \cap N_p{\mathcal {M}}\) we have
For the sake of simplicity, let us set \(A_p:= Df_{p}^{{\textsf{T}}}Df_{p}\). Then, by the definition of \({\mathcal {M}}\) and Proposition 4, we have \( \ker (A_p) = \ker (Df_{p}) = T_{p}{\mathcal {M}}. \) Now, let us decompose \({\mathbb {R}}^n = \ker (A_p) \oplus \mathop {\textrm{ran}}\limits (A_p) = T_{p}{\mathcal {M}}\oplus N_{p}{\mathcal {M}}\) and accordingly also \(A_p = 0 \oplus \widetilde{A_p}\). Since \(\widetilde{A_p}\) is symmetric positive definite, for any \(x\in N_p{\mathcal {M}}\) we obtain
Now, using (4) and (6), we obtain
where we set \(c_p:=\lambda _{\min }(\widetilde{A_p})^{1/2} > 0\). In particular, we have
Thus, choosing \(x \in N_p{\mathcal {M}}\) sufficiently close to p we obtain the estimate (5), and we can deduce that an open set \(U_p\) with the sought-after properties has to exist. Now, since (3) and \(p\in {\mathcal {M}}\) imply
using (5) we obtain
To finish the local argument, notice that by (6) we have
which together with (7) shows that (2) holds for all \(x\in V_p\).
To construct the set V, first notice that the differentiability of f implies that an expression of the form (3) is possible on the set \({\mathcal {E}}({\mathcal {M}})\). In other words, there exists a function \(g_{\cdot }(\cdot )\) such that
for all \(x\in {\mathcal {E}}({\mathcal {M}})\). Since the orthogonal projection \(x\mapsto p_x\) is differentiable [26] and \(f\in C^2\), the map \(x \mapsto g_{p_x}(x-p_x)\) is continuous on \({\mathcal {E}}({\mathcal {M}})\). Define the function
Then h is continuous and hence the preimage of the open set \((-\infty ,\tfrac{1}{2})\subseteq {\mathbb {R}}\) under h is open. Note that \({\mathcal {M}}\) is a subset of this preimage. Define
Then for each \(p\in {\mathcal {M}}\) we have
The previous arguments can then be used to show that for \(c:= \tfrac{{\widetilde{c}}}{2}\) and \(x \in V\) we have
finishing the proof. \(\square \)
Remark 6
The estimate (2) is related to the Łojasiewicz inequality [28, 29], which states that for a real analytic function \(g:U \rightarrow {\mathbb {R}}\) defined on an open set \(U\subseteq {\mathbb {R}}^n\) and a compact set \(K \subseteq U\), the distance of \(x \in K\) to the zero locus \({\mathcal {Z}}:= \{ z \in U ~|~ g(z) = 0 \}\) of g may be estimated by
where \(\alpha \) and C are positive constants.
4 Manifold dissipativity and manifold turnpikes
In this section, we recall the definition of dissipativity with respect to a manifold and the definition of manifold turnpikes as introduced in [10]. Further, a theorem relating the two properties is stated. Here and in the following, the set \({\mathcal {K}}\) is defined as
We consider the general optimal control problem
As before, we assume \(x_0, x_T \in K\), where \(K \subseteq {\mathbb {R}}^n\) is a compact set. Here, the function g defines the dynamics of the system and the function h corresponds to possible algebraic constraints. We refrain from further specification of these functions as (\(\text {OCP}_T\)) is only used for general definitions. Throughout this section, we assume that an optimal control \(u^*\) of (\(\text {OCP}_T\)) and an associated trajectory \(x^*\) exist.
We begin with the definition of manifold dissipativity. The definition is related to Willems’ notion of dissipativity [24] and is also found in [10].
Definition 7
(manifold dissipativity) Consider the optimal control problem (\(\text {OCP}_T\)) together with the manifold \({\mathcal {M}}\subseteq {\mathbb {R}}^n\). We say that (\(\text {OCP}_T\)) is dissipative with respect to the manifold \({\mathcal {M}}\) if there exists a function \({\mathcal {S}}:{\mathbb {R}}^n \rightarrow [0,\infty )\) that is bounded on compact sets and a function \(\alpha \in {\mathcal {K}}\) such that all optimal controls \(u^*\) and associated trajectories \(x^*\) satisfy the dissipation inequality
for all \(T >0\).
The function \({\mathcal {S}}\) from Definition 7 is also called storage function. Note that we require the dissipation inequality (8) only to hold along optimal solutions of (\(\text {OCP}_T\)). This is not a severe restriction when turnpike phenomena are studied, as we are only interested in properties of optimal solutions.
Next, we define a manifold turnpike property, again following [10]. The property is essentially a notion of measure turnpikes, see, e.g., [3,4,5, 30].
Definition 8
(manifold turnpike) Consider the optimal control problem (\(\text {OCP}_T\)) together with the manifold \({\mathcal {M}}\subseteq {\mathbb {R}}^n\). We say that (\(\text {OCP}_T\)) has the manifold turnpike property with respect to the manifold \({\mathcal {M}}\) if for all compact sets \(K\subseteq {\mathbb {R}}^n\) and all \(\varepsilon >0\) there exists a constant \(C_{K,\varepsilon } > 0\) such that for all \(T>0\) all optimal trajectories \(x^*\) of (\(\text {OCP}_T\)) satisfy
for all \(x_0,x_T \in K\). Here, \(\lambda \) denotes the Lebesgue-measure.
The next theorem can be found similarly in [10, 30, 31]. The theorem shows that manifold dissipativity implies a manifold turnpike property.
Theorem 9
(manifold dissipativity implies manifold turnpike) Consider the optimal control problem (\(\text {OCP}_T\)) together with a submanifold \({\mathcal {M}}\subseteq {\mathbb {R}}^n\) and assume that
-
(i)
there exists a constant \(C_\ell (K)>0\) such that for all optimal controls \(u^*\) of (\(\text {OCP}_T\)) and the associated trajectories \(x^*\) we have
$$\begin{aligned} \int _0^T \ell (x^*,u^*) ~\textrm{d}t< C_\ell (K) \end{aligned}$$for all \(T > 0\), and
-
(ii)
the optimal control problem is dissipative with respect to the manifold \({\mathcal {M}}\).
Then the optimal control problem (\(\text {OCP}_T\)) has the manifold turnpike property.
Proof
Rearranging (8), we see that
Since the set K is compact and S is bounded on compact sets, there exists a constant \(C_S(K) > 0\) independent of T such that
Now define
and notice that
Thus, for all \(T > 0\) we have
finishing the proof. \(\square \)
Remark 10
Denote the optimal value function of (\(\text {pH} \text { OCP}_T\)) by \(V_T(x_0)\), that is
Then, assumption (i) of Theorem 9 is equivalent to
In [10], assumption (i) is replaced by the assumption that (\(\text {pH} \text { OCP}_T\)) is cost controllable. The latter takes \({\mathcal {U}}_{\textrm{ad}}= L^\infty ([0,T], {\mathbb {R}}^m)\) and bounds the growth of the optimal value function by
where the growth function \(b_K:[0,\infty ) \rightarrow [0,\infty )\) is assumed to be bounded and increasing and \(\ell ^*(x_0):= \inf _{u\in {\mathbb {R}}^m} \ell (x_0,u)\) is assumed to be real for all \(x_0 \in K\). Provided that \(\ell \) is bounded on compact sets, [10] shows that (10) implies (9) as we can bound \(b_K(T)\,\ell ^*(x_0) \le C_\ell (K)\) in this case.
5 Application to port-Hamiltonian systems
Finally, we are ready to apply the previous results to port-Hamiltonian systems and the optimal control problem (\(\text {pH} \text { OCP}_T\)). First, recall that we can rewrite the cost functional \(C_{\text {pH},T}(u)\) as
and that rearranging gives
Equation (11) hints that any optimal trajectory will have to spend most of the time close to the set
and that \({\mathcal {H}}\) can be used as a storage function to derive dissipativity notions with respect to \({\mathcal {M}}\). Our aim will be to formalize these ideas. The first step will be to ensure that \({\mathcal {M}}\) has the necessary manifold structure. For that, we make the following assumptions. See also Remark 16.
Assumption 11
-
(A1)
The map \(x\mapsto R(x)^{1/2} \eta (x)\) is of class \(C^2\).
-
(A2)
The set \({\mathcal {M}}\) is nonempty and there exists an open neighborhood \(G\subseteq {\mathbb {R}}^n\) of \({\mathcal {M}}\) such that for all \(x\in G\) it holds that \(\dim (\ker (D(R^{1/2}\eta )_x)) = s\), where the constant s is independent of x and satisfies \(0< s < n\).
-
(A3)
For each \(x\in G\), the smallest nonzero singular value of \(D(R^{1/2}\eta )_x\) is bounded from below by a positive constant \({\widetilde{c}} > 0\).
-
(A4)
Let V be the open set from Lemma 5. Any optimal trajectory \(x^*\) of (\(\text {pH} \text { OCP}_T\)) remains in V for all times.
With assumptions (A1) and (A2), Proposition 4 ensures that \({\mathcal {M}}\) is an s-dimensional \(C^{2,1}\) submanifold of \({\mathbb {R}}^n\) and that \({\mathcal {M}}\subseteq {\mathcal {E}}({\mathcal {M}})\).
The next step is to show that the problem is dissipative. As we will see shortly, this is ensured by assumptions (A3) and (A4), which allow us to use Lemma 5 to conclude that the optimal control problem (\(\text {pH} \text { OCP}_T\)) is dissipative with respect to the manifold \({\mathcal {M}}\). Notice that (A4) implies \(x_0,x_T \in V\).
Theorem 12
((\(\text {pH} \text { OCP}_T\)) is dissipative) Under assumption 11, the optimal control problem (\(\text {pH} \text { OCP}_T\)) is dissipative with respect to the manifold \({\mathcal {M}}\) with storage function \({\mathcal {H}}\).
Proof
The proof is essentially an application of Lemma 5. Since all assumptions of Lemma 5 are satisfied under assumption 11, there exists an open set \(V\subseteq {\mathbb {R}}^n\) and constant \(c>0\) such that \({\mathcal {M}}\subseteq V \subseteq {\mathcal {E}}({\mathcal {M}})\) and
holds for all \(x\in V\). In particular, assumption (A4) ensures that the estimate holds along any optimal trajectory \(x^*\) of (\(\text {pH} \text { OCP}_T\)). With this and (12), we see that for any optimal control \(u^*\) and the associated trajectory \(x^*\) and output \(y^*\), we have
where \(\alpha :s \mapsto c^2 s^2 \in {\mathcal {K}}\). Finally, note that the Hamiltonian \({\mathcal {H}}\), acting as a storage function here, is bounded on compact sets since it is differentiable. \(\square \)
Remark 13
Let us emphasize that the estimate (2) from Lemma 5 was the key to conclude the dissipativity property of (\(\text {pH} \text { OCP}_T\)). As we have mentioned in Remark 6, the estimate is related to Łojasiewicz’ inequality [28, 29]. In fact, if the map \(g:x \mapsto \Vert R^{1/2}(x)\eta (x) \Vert ^2\) is real analytic, then we may use the Łojasiewicz inequality to derive dissipativity without Lemma 5, as long as any optimal trajectory stays in some compact set \(D\subseteq {\mathbb {R}}^n\).
Now, an application of Theorem 9 yields the following result, showing that the optimal control problem (\(\text {pH} \text { OCP}_T\)) has the manifold turnpike property with respect to \({\mathcal {M}}\).
Theorem 14
((\(\text {pH} \text { OCP}_T\)) has a manifold turnpike) In addition to assumption 11, assume that
-
(A5)
there exists a trajectory \(x_1\) that starts in \(x_0\) and reaches the manifold \({\mathcal {M}}\) in time \(T_1\ge 0\), and
-
(A6)
for any \(p\in {\mathcal {M}}\), there exists a trajectory \(x_2\) that starts in p and reaches \(x_T\) in time \(T_2 \ge 0\).
Then the optimal control problem (\(\text {pH} \text { OCP}_T\)) has a manifold turnpike at the manifold \({\mathcal {M}}\).
Proof
Since the port-Hamiltonian system (pH) is autonomous and there is no control cost on the manifold \({\mathcal {M}}\), the conditions (A5) and (A6) ensure that the total cost of the optimal control \(u^*\) is bounded by
Notice that the constants \(C_{{\mathcal {H}}}(K), ~C_1(K)\) and \(C_2(K)\) are independent of the final time T. Thus, using the results of Theorem 12, we may apply Theorem 9 to conclude that the optimal control problem (\(\text {pH} \text { OCP}_T\)) has the manifold turnpike property at \({\mathcal {M}}\). \(\square \)
In Theorem 14, in order to show that a turnpike property holds true, we needed to make the controllability assumptions (A5) and (A6). This is a common pattern in turnpike results, similar controllability assumptions are made in [3, 10] and [19].
Remark 15
In [19], the authors considered linear pH systems of the form
where \(J=-J^{{\textsf{T}}}\), \(R=R^{{\textsf{T}}}\succeq 0\) and \(Q=Q^{{\textsf{T}}}\succ 0\). They have shown that in this case the optimal control problem (\(\text {pH} \text { OCP}_T\)) admits a subspace turnpike property with respect to \(\ker (R^{1/2} Q)\). We can interpret Theorem 14 as a generalization of this result to the nonlinear case. In the linear case, assumption (A1) is immediately satisfied. Further, the set \({\mathcal {M}}\) is the kernel of \(R^{1/2} Q\), and if the dimension of \(\ker (R^{1/2} Q)\) is \(0<s<n\), then assumptions (A2) and (A3) are also satisfied. Since the distance estimate (13) can be shown to hold true globally [19, Lemma 13], the set V from Lemma 5 is \(V={\mathbb {R}}^n\) and thus assumption (A4) is also satisfied.
Remark 16
In many real-world examples, such as, for example, port-Hamiltonian formulations of the isothermal Euler equations arising from gas network modeling [32], the function R has additional structural properties. One example is block structure, where the manifold \({\mathcal {M}}\) often is a linear subspace. In this case, assumption (A2) usually does not pose a problem. On the other hand, even in this situation, assumption (A3) can be difficult to ensure.
If \(x_0, x_T \in V\) and the controllability assumptions of Theorem 14 are satisfied, then assumption (A4) is not critical in many cases. Examining (11), we observe that leaving the manifold \({\mathcal {M}}\) is connected to a higher cost. Hence, any optimal trajectory should approach the manifold \({\mathcal {M}}\), and thus, unless the set V is fairly irregular and \(x_0\) and \(x_T\) are chosen poorly, optimal trajectories will stay within V.
We finish this section with a simple example illustrating Theorem 14.
Example 17
Consider the functions \(E,J,R,\eta \) and B defined by
for all \(x\in {\mathbb {R}}^2\), which together with the Hamiltonian
form the port-Hamiltonian system
For the system (pH-1), the function \(R^{1/2}\eta \) is given by
where we take \(x=\begin{bmatrix} x_1&x_2 \end{bmatrix}^{{\textsf{T}}}\). Thus, the derivative \(D(R^{1/2}\eta )\) reads as
and the subspace \(\ker (D(R^{1/2}\eta )(x))\) is one-dimensional for all \(x\in {\mathbb {R}}^2\). A simple calculation shows that the nonzero singular value of \(D(R^{1/2}\eta )(x)\) is given by
The zero locus \({\mathcal {M}}\) of \(R^{1/2}\eta \) is given by
Hence, assumptions (A1), (A2) and (A3) are satisfied for (pH-1). As \({\mathcal {M}}\) is a linear subspace, the orthogonal projection on \({\mathcal {M}}\) is well defined globally and we have \({\mathcal {E}}({\mathcal {M}})={\mathbb {R}}^2\). Further, for the set V from Lemma 5 it holds that \(V={\mathbb {R}}^2\) since
for all \(x\in {\mathbb {R}}^2\). Thus, also assumption (A4) is satisfied for (pH-1).
Now, for \(\xi = \begin{bmatrix} \xi _1&\xi _2&\xi _3 \end{bmatrix}^{{\textsf{T}}}\in {\mathbb {R}}^3\), let us define the functions \({\widetilde{E}},{\widetilde{J}},{\widetilde{R}},{\widetilde{\eta }}\) and \({\widetilde{B}}\) by
It is easy to see that also the system
satisfies assumptions (A1) – (A4). The zero locus of the map \(\xi \mapsto {\widetilde{R}}(\xi )^{1/2}{\widetilde{\eta }}(\xi )\) is
6 Numerical example
As an example, we consider the optimal control problem (\(\text {pH} \text { OCP}_T\)) together with the pH systems (pH-1) and (pH-2) from Example 17. For the implementation, we use the open-source software package CasADi [33].
In order to use CasADi, the optimal control problem (\(\text {pH} \text { OCP}_T\)) is formulated as a minimization problem of the form
We follow a similar procedure as [33, Sect. 5.4]. In our implementation, w contains the values \(x(t_i)\) for the discretization points \(t_i \in [0,T]\), and the values \(u(t_i)\) for the discretization points \(t_i \in [0,T], ~ i \ne 0\). The initial condition and possible control constraints are incorporated in \(w_{\text {lb}}\) and \(w_{\text {ub}}\). The function G is used to enforce the final condition and a continuity condition on x by using an integrator scheme to determine the value \(x(t_{i+1})\) given the values \(x(t_i)\) and \(u(t_i)\). This integrator scheme is also used to calculate the cost J via the quad option in CasADi’s integrator function. For the solution of the nonlinear optimization problem (14), Ipopt [34] is used.
In Fig. 1, the solution of the optimal control problem (\(\text {pH} \text { OCP}_T\)) with the system (pH-1) under the control constraint \(-50 \le u(t) \le 50\) is shown. The turnpike behaviour is clearly visible; the first component \(x_1\) of the optimal trajectory \(x^*\) approaches the manifold \({\mathcal {M}}= \{ x \in {\mathbb {R}}^2 ~|~ x_1 = 0 \}\) very quickly and remains there for the majority of the time horizon. The same observation can be made for larger time horizons, which is not shown in Fig. 1.
In Fig. 2, the solution of the optimal control problem (\(\text {pH} \text { OCP}_T\)) with the system (pH-2) under the control constraint \(-200 \le u(t) \le 200\) is shown. Again, the turnpike phenomenon can be observed.
7 Conclusion
In this paper, we have considered the optimal control of port-Hamiltonian systems under minimal energy supply with fixed initial and final values. We have seen that the map \(R^{1/2}\eta \), corresponding to the energy dissipating portion of the right-hand side, and its zero locus \({\mathcal {M}}= \{x ~|~ R^{1/2}(x)\eta (x) = 0 \}\), which corresponds to the dissipative part of the state space, play an important role. It was shown that under smoothness assumptions on \(R^{1/2}\eta \), the set \({\mathcal {M}}\) forms a \(C^2\) submanifold of \({\mathbb {R}}^n\). In particular, using results from [26], we observed that the orthogonal projection onto \({\mathcal {M}}\) is well-defined in an open set \({\mathcal {E}}({\mathcal {M}})\). Further, we have shown that under these assumptions the distance of a point x to \({\mathcal {M}}\) can essentially be bounded by \(\Vert R^{1/2}(x)\eta (x) \Vert \) from above. This fact allowed us to deduce that the considered optimal control problem is dissipative with respect to the manifold \({\mathcal {M}}\). Our main result was a consequence of this dissipativity property. Under additional controllability assumptions, we have seen that the problem has a manifold turnpike property with respect to \({\mathcal {M}}\). This theoretical observation was confirmed in a simple numerical example.
An open question from a theoretical perspective is whether the problem (\(\text {pH} \text { OCP}_T\)) exhibits hidden regularities as in the linear case [35]. Another open topic is the study of stronger turnpike properties, such as exponential turnpikes [31, 36]. Applications of the theoretical results to specific port-Hamiltonian systems such as gas networks [32] will be studied in future works.
References
Dorfman R, Samuelson P, Solow R (1987) Linear programming and economic analysis. Dover books on advanced mathematics. Dover Publications, New York. https://books.google.de/books?id=k5_vzaCNQP4C
McKenzie L (1963) Turnpike theorems for a generalized Leontief model. Econometrica 31(1/2):165. https://doi.org/10.2307/1910955
Carlson D, Haurie A, Leizarowitz A (1991) Infinite horizon optimal control. Springer, Berlin. https://doi.org/10.1007/978-3-642-76755-5
Zaslavski A (2006) Turnpike properties in the calculus of variations and optimal control. Springer, Heidelberg. https://doi.org/10.1007/0-387-28154-1
Trélat E, Zhang C (2018) Integral and measure-turnpike properties for infinite-dimensional optimal control systems. Math Control Signals Syst 30(1):1. https://doi.org/10.1007/s00498-018-0209-1
Porretta A, Zuazua E (2013) Long time versus steady state optimal control. SIAM J Control Optim 51(6):4242–4273. https://doi.org/10.1137/130907239
Porretta A, Zuazua E (2016) Remarks on long time versus steady state optimal control. In: Mathematical paradigms of climate science. Springer, Basel, pp 67–89. https://doi.org/10.1007/978-3-319-39092-5_5
Trélat E, Zuazua E (2015) The turnpike property in finite-dimensional nonlinear optimal control. J Differ Equ 258(1):81–114. https://doi.org/10.1016/j.jde.2014.09.005
Rapaport A, Cartigny P (2004) Turnpike theorems by a value function approach. ESAIM Control Optim Calc Variat 10(1):123–141. https://doi.org/10.1051/cocv:2003039
Faulwasser T, Flaßkamp K, Ober-Blöbaum S, Schaller M, Worthmann K (2022) Manifold turnpikes, trims, and symmetries. Math Control Signals Syst 34(4):759–788. https://doi.org/10.1007/s00498-022-00321-6
Paynter H (1961) Analysis and Design of Engineering Systems: Class Notes for M.I.T. Course 2,751. M.I.T. Press, Cambridge. https://books.google.de/books?id=FcdhAAAAMAAJ
Van der Schaft A, Jeltsema D (2014) Port-Hamiltonian systems theory: an introductory overview. Found Trends Syst Control 1(2–3):173–378. https://doi.org/10.1561/2600000002
Moser T, Lohmann B (2022) Rosenbrock framework for tangential interpolation of port-Hamiltonian descriptor systems. https://arxiv.org/abs/2210.16071
Lamoline F, Hastir A (2022) On Dirac structure of infinite-dimensional stochastic port-Hamiltonian systems. https://arxiv.org/abs/2210.06358
Van der Schaft A, Mehrmann V (2022) Linear port-Hamiltonian DAE systems revisited. https://arxiv.org/abs/2211.06676
Breiten T, Hinsen D, Unger B (2022) Towards a modeling class for port-Hamiltonian systems with time-delay. https://arxiv.org/abs/2211.10687
Gernandt H, Hinsen D (2022) Stability and passivity for a class of infinite dimensional port-Hamiltonian networks with application to power grids. https://arxiv.org/abs/2212.02792
Mehrmann V, Unger B (2023) Control of port-Hamiltonian differential-algebraic systems and applications. Acta Numerica 32:395–515. https://doi.org/10.1017/S0962492922000083
Schaller M, Philipp F, Faulwasser T, Worthmann K, Maschke B (2021) Control of port-Hamiltonian systems with minimal energy supply. Eur J Control 62:33–40. https://doi.org/10.1016/j.ejcon.2021.06.017
Faulwasser T, Maschke B, Philipp F, Schaller M, Worthmann K (2022) Optimal control of port-Hamiltonian descriptor systems with minimal energy supply. SIAM J Control Optim 60(4):2132–2158. https://doi.org/10.1137/21m1427723
Philipp F, Schaller M, Faulwasser T, Maschke B, Worthmann K (2021) Minimizing the energy supply of infinite-dimensional linear port-Hamiltonian systems. IFAC-PapersOnLine 54(19):155–160. https://doi.org/10.1016/j.ifacol.2021.11.071
Staffans O (2002) Passive and conservative infinite-dimensional impedance and scattering systems (from a personal point of view). IMA Vol Math Appl 314:1. https://doi.org/10.1007/978-0-387-21696-6_14
Ortega R, Van der Schaft A, Maschke B, Escobar G (2002) Interconnection and damping assignment passivity-based control of port-controlled Hamiltonian systems. Automatica 38(4):585–596. https://doi.org/10.1016/s0005-1098(01)00278-3
Willems J (1972) Dissipative dynamical systems part I: general theory. Arch Ration Mech Anal 45(5):321–351. https://doi.org/10.1007/bf00276493
Dudek E, Holly K (1994) Nonlinear orthogonal projection. Annales Polonici Mathematici 59(1):1–31
Leobacher G, Steinicke A (2021) Existence, uniqueness and regularity of the projection onto differentiable manifolds. Ann Glob Anal Geom 60:559–587. https://doi.org/10.1007/s10455-021-09788-z
Walter R (2009) Einführung in die Analysis 3. De Gruyter, Berlin. https://doi.org/10.1515/9783110214840
Łojasiewicz S (1959) Sur le probleme de la division. Studia Mathematica 18(1):87–136
Ji S, Kollár J, Shiffman B (1992) A global Lojasiewicz inequality for algebraic varieties. Trans Am Math Soc 329(2):813. https://doi.org/10.2307/2153965
Faulwasser T, Korda M, Jones C, Bonvin D (2017) On turnpike and dissipativity properties of continuous-time optimal control problems. Automatica 81:297–304. https://doi.org/10.1016/j.automatica.2017.03.012
Grüne L, Guglielmi R (2021) On the relation between turnpike properties and dissipativity for continuous time linear quadratic optimal control problems. Math Control Rel Fields 11(1):169–188. https://doi.org/10.3934/mcrf.2020032
Domschke P, Hiller B, Lang J, Mehrmann V, Morandin R, Tischendorf C (2021) Gas network modeling: an overview. https://opus4.kobv.de/opus4-trr154/frontdoor/index/index/docId/411
Andersson JAE, Gillis J, Horn G, Rawlings JB, Diehl M (2019) CasADi: a software framework for nonlinear optimization and optimal control. Math Program Comput 11:1–36. https://doi.org/10.1007/s12532-018-0139-4
Wächter A, Biegler L (2006) On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming. Math Program 106(1):25–57. https://doi.org/10.1007/s10107-004-0559-y
Faulwasser T, Kirchhoff J, Mehrmann V, Philipp F, Schaller M, Worthmann K (2023) Hidden regularity in singular optimal control of port-Hamiltonian systems. https://doi.org/10.48550/arXiv.2305.03790
Trélat E, Zhang C, Zuazua E (2018) Steady-state and periodic exponential turnpike property for optimal control problems in Hilbert spaces. SIAM J Control Optim 56(2):1222–1252. https://doi.org/10.1137/16m1097638
Kunkel P, Mehrmann V (2006) Differential-algebraic equations—analysis and numerical solution. EMS Textbooks in Mathematics, Zürich, S. https://doi.org/10.4171/017
Macki J, Strauss A (1982) Introduction to optimal control theory. Springer, New York, NY. https://doi.org/10.1007/978-1-4612-5671-7
Kunkel P, Mehrmann V (2008) Optimal control for unstructured nonlinear differential-algebraic equations of arbitrary index. Math Control Signals Syst 20(3):227–269. https://doi.org/10.1007/s00498-008-0032-1
Acknowledgements
The author thanks Tobias Breiten, Bernhard Höveler and Volker Mehrmann for their helpful comments on an early version of this manuscript. Further, the author thanks the Deutsche Forschungsgemeinschaft for their support within the subproject B03 in the Sonderforschungsbereich/Transregio 154 “Mathematical Modelling, Simulation and Optimization using the Example of Gas Networks” (Project https://gepris.dfg.de/gepris/projekt/239904186).
Funding
Open Access funding enabled and organized by Projekt DEAL. The author is supported within the subproject B03 in the Sonderforschungsbereich/Transregio 154 “Mathematical Modelling, Simulation and Optimization using the Example of Gas Networks” (Project https://gepris.dfg.de/gepris/projekt/239904186).
Author information
Authors and Affiliations
Contributions
A.K. prepared the manuscript without further contributions.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no Conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Existence of solutions to (pH)
To study the existence of solutions to (pH) for arbitrary final times T, we consider the time interval \([0,\infty )\) and assume \(u\in {\mathcal {U}}_{\textrm{ad}}\) is fixed. We may then write the system (pH) in the general form
Assuming that F is sufficiently smooth and satisfies [37, Hypothesis 4.2] uniformly in \(t\in [0,\infty )\), we arrive at the reduced differential algebraic equation
Adapting [37, Theorem 4.13], we get the following result.
Proposition 18
Assume that F is sufficiently smooth and satisfies [37, Hypothesis 4.2] with characteristic values \(\mu ,a,d\) and \(\mu +1,a,d\) uniformly in \(t\in [0,\infty )\). Then, for any suitable initial condition, the reduced problem (16) has a unique solution defined on \([0,\infty )\), and this solution locally solves the original problem (15).
Remark 19
We may also combine x and u in the generalized state vector \(z=(x,u)\) and rewrite (15) accordingly to arrive at a version of Proposition 18 that does not depend on one particular control \(u\in {\mathcal {U}}_{\textrm{ad}}\).
Existence of solutions to (\(\text {pH} \text { OCP}_T\))
In order to state sufficient conditions for the existence of optimal controls \(u^*\) and associated trajectories \(x^*\) for (\(\text {pH} \text { OCP}_T\)), we restrict ourselves to the ODE case and assume that E(x) is pointwise invertible with sufficiently smooth inverse and use results of [38].
Let \(T>0\) be a final time, \({\mathbb {U}} \subseteq {\mathbb {R}}^m\) be compact. We consider the problem
subject to
with \(\ell \) and g continuous, and set
Theorem 20
( [38, Theorem IV.2]) Consider the problem (17), (18), assume \({\mathcal {U}}_{\textrm{bvp}}\ne \emptyset \) and that all solutions to the boundary value problem (18) satisfy an a priori bound
where \(\alpha = \alpha (T)\) is a constant depending only on T. If the set
is convex for each \(x \in {\mathbb {R}}^n\), then there exists an optimal control.
In [38], Theorem 20 is stated for the case of free final time \(t_1 \le T\) and \(x(t_1)=0\). The remarks made after [38, Theorems IV.1 and IV.2] justify the modifications made here. Further, as mentioned in [38, Example IV.4], the convexity assumption is satisfied for the dynamics (pH) and a cost term of the form \(y^{{\textsf{T}}}u\). Since the bound (19) is satisfied if
for some constants \(\alpha , \beta \ge 0\) and all \((x,v) \in {\mathbb {R}}^n \times {\mathbb {U}}\), see [38], this leaves us with the following result.
Proposition 21
Consider the optimal control problem (\(\text {pH} \text { OCP}_T\)) with \({\mathcal {U}}_{\textrm{ad}}= L^{\infty }([0,T],{\mathbb {U}})\), where \({\mathbb {U}} \subseteq {\mathbb {R}}^m\) is compact. Assume that E is pointwise invertible with continuous inverse, and that for
either of the bounds in (20) holds. If there exists some \(u\in {\mathcal {U}}_{\textrm{ad}}\) that steers (pH) from \(x_0\) to \(x_T\) in time T, then there exists an optimal control.
In the general differential-algebraic case, the situation is more involved, and sufficient conditions are difficult to obtain. A thorough treatment of necessary conditions can be found in [39], where a boundary value problem of necessary optimality conditions is stated. The result can be adapted to problems that involve end conditions. As mentioned in the reference, the boundary value problem can in general not be formulated in terms of the original data. In particular, the port-Hamiltonian structure is likely lost.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Karsai, A. Manifold turnpikes of nonlinear port-Hamiltonian descriptor systems under minimal energy supply. Math. Control Signals Syst. 36, 707–728 (2024). https://doi.org/10.1007/s00498-024-00384-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00498-024-00384-7