Distributed parameter system

In control theory, a distributed-parameter system (as opposed to a lumped-parameter system) is a system whose state space is infinite-dimensional. Such systems are therefore also known as infinite-dimensional systems. Typical examples are systems described by partial differential equations or by delay differential equations.

Linear time-invariant distributed-parameter systems

Abstract evolution equations

Discrete-time

With U, X and Y Hilbert spaces and ${\displaystyle A\,}$ ∈ L(X), ${\displaystyle B\,}$ ∈ L(U, X), ${\displaystyle C\,}$ ∈ L(X, Y) and ${\displaystyle D\,}$ ∈ L(U, Y) the following difference equations determine a discrete-time linear time-invariant system:

${\displaystyle x(k+1)=Ax(k)+Bu(k)\,}$

${\displaystyle y(k)=Cx(k)+Du(k)\,}$

with ${\displaystyle x\,}$ (the state) a sequence with values in X, ${\displaystyle u\,}$ (the input or control) a sequence with values in U and ${\displaystyle y\,}$ (the output) a sequence with values in Y.

Continuous-time

The continuous-time case is similar to the discrete-time case but now one considers differential equations instead of difference equations:

${\displaystyle {\dot {x}}(t)=Ax(t)+Bu(t)\,}$,

${\displaystyle y(t)=Cx(t)+Du(t)\,}$.

An added complication now however is that to include interesting physical examples such as partial differential equations and delay differential equations into this abstract framework, one is forced to consider unbounded operators. Usually A is assumed to generate a strongly continuous semigroup on the state space X. Assuming B, C and D to be bounded operators then already allows for the inclusion of many interesting physical examples,[1] but the inclusion of many other interesting physical examples forces unboundedness of B and C as well.

Example: a partial differential equation

The partial differential equation with ${\displaystyle t>0}$ and ${\displaystyle \xi \in [0,1]}$ given by

${\displaystyle {\frac {\partial }{\partial t}}w(t,\xi )=-{\frac {\partial }{\partial \xi }}w(t,\xi )+u(t),}$

${\displaystyle w(0,\xi )=w_{0}(\xi ),}$

${\displaystyle w(t,0)=0,}$

${\displaystyle y(t)=\int _{0}^{1}w(t,\xi )\,d\xi ,}$

fits into the abstract evolution equation framework described above as follows. The input space U and the output space Y are both chosen to be the set of complex numbers. The state space X is chosen to be L^2(0, 1). The operator A is defined as

${\displaystyle Ax=-x',~~~D(A)=\left\{x\in X:x{\text{ absolutely continuous }},x'\in L^{2}(0,1){\text{ and }}x(0)=0\right\}.}$

It can be shown[2] that A generates a strongly continuous semigroup on X. The bounded operators B, C and D are defined as

${\displaystyle Bu=u,~~~Cx=\int _{0}^{1}x(\xi )\,d\xi ,~~~D=0.}$

Example: a delay differential equation

The delay differential equation

${\displaystyle {\dot {w}}(t)=w(t)+w(t-\tau )+u(t),}$

${\displaystyle y(t)=w(t),}$

fits into the abstract evolution equation framework described above as follows. The input space U and the output space Y are both chosen to be the set of complex numbers. The state space X is chosen to be the product of the complex numbers with L^2(−τ, 0). The operator A is defined as

${\displaystyle A{\begin{pmatrix}r\\f\end{pmatrix}}={\begin{pmatrix}r+f(-\tau )\\f'\end{pmatrix}},~~~D(A)=\left\{{\begin{pmatrix}r\\f\end{pmatrix}}\in X:f{\text{ absolutely continuous }},f'\in L^{2}([-\tau ,0]){\text{ and }}r=f(0)\right\}.}$

It can be shown[3] that A generates a strongly continuous semigroup on X. The bounded operators B, C and D are defined as

${\displaystyle Bu={\begin{pmatrix}u\\0\end{pmatrix}},~~~C{\begin{pmatrix}r\\f\end{pmatrix}}=r,~~~D=0.}$

Transfer functions

As in the finite-dimensional case the transfer function is defined through the Laplace transform (continuous-time) or Z-transform (discrete-time). Whereas in the finite-dimensional case the transfer function is a proper rational function, the infinite-dimensionality of the state space leads to irrational functions (which are however still holomorphic).

Discrete-time

In discrete-time the transfer function is given in terms of the state-space parameters by ${\displaystyle D+\sum _{k=0}^{\infty }CA^{k}Bz^{k}}$ and it is holomorphic in a disc centered at the origin.[4] In case 1/z belongs to the resolvent set of A (which is the case on a possibly smaller disc centered at the origin) the transfer function equals ${\displaystyle D+Cz(I-zA)^{-1}B}$. An interesting fact is that any function that is holomorphic in zero is the transfer function of some discrete-time system.

Continuous-time

If A generates a strongly continuous semigroup and B, C and D are bounded operators, then[5] the transfer function is given in terms of the state space parameters by ${\displaystyle D+C(sI-A)^{-1}B}$ for s with real part larger than the exponential growth bound of the semigroup generated by A. In more general situations this formula as it stands may not even make sense, but an appropriate generalization of this formula still holds.[6] To obtain an easy expression for the transfer function it is often better to take the Laplace transform in the given differential equation than to use the state space formulas as illustrated below on the examples given above.

Transfer function for the partial differential equation example

Setting the initial condition ${\displaystyle w_{0}}$ equal to zero and denoting Laplace transforms with respect to t by capital letters we obtain from the partial differential equation given above

${\displaystyle sW(s,\xi )=-{\frac {d}{d\xi }}W(s,\xi )+U(s),}$

${\displaystyle W(s,0)=0,}$

${\displaystyle Y(s)=\int _{0}^{1}W(s,\xi )\,d\xi .}$

This is an inhomogeneous linear differential equation with ${\displaystyle \xi }$ as the variable, s as a parameter and initial condition zero. The solution is ${\displaystyle W(s,\xi )=U(s)(1-e^{-s\xi })/s}$. Substituting this in the equation for Y and integrating gives ${\displaystyle Y(s)=U(s)(e^{-s}+s-1)/s^{2}}$ so that the transfer function is ${\displaystyle (e^{-s}+s-1)/s^{2}}$.

Transfer function for the delay differential equation example

Proceeding similarly as for the partial differential equation example, the transfer function for the delay equation example is[7] ${\displaystyle 1/(s-1-e^{-s})}$.

Controllability

In the infinite-dimensional case there are several non-equivalent definitions of controllability which for the finite-dimensional case collapse to the one usual notion of controllability. The three most important controllability concepts are:

• Exact controllability,
• Approximate controllability,
• Null controllability.

Controllability in discrete-time

An important role is played by the maps ${\displaystyle \Phi _{n}}$ which map the set of all U valued sequences into X and are given by ${\displaystyle \Phi _{n}u=\sum _{k=0}^{n}A^{k}Bu_{k}}$. The interpretation is that ${\displaystyle \Phi _{n}u}$ is the state that is reached by applying the input sequence u when the initial condition is zero. The system is called

• exactly controllable in time n if the range of ${\displaystyle \Phi _{n}}$ equals X,
• approximately controllable in time n if the range of ${\displaystyle \Phi _{n}}$ is dense in X,
• null controllable in time n if the range of ${\displaystyle \Phi _{n}}$ includes the range of A^n.

Controllability in continuous-time

In controllability of continuous-time systems the map ${\displaystyle \Phi _{t}}$ given by ${\displaystyle \int _{0}^{t}{\rm {e}}^{As}Bu(s)\,ds}$ plays the role that ${\displaystyle \Phi _{n}}$ plays in discrete-time. However, the space of control functions on which this operator acts now influences the definition. The usual choice is L^2(0, ∞;U), the space of (equivalence classes of) U-valued square integrable functions on the interval (0, ∞), but other choices such as L^1(0, ∞;U) are possible. The different controllability notions can be defined once the domain of ${\displaystyle \Phi _{t}}$ is chosen. The system is called[8]

• exactly controllable in time t if the range of ${\displaystyle \Phi _{t}}$ equals X,
• approximately controllable in time t if the range of ${\displaystyle \Phi _{t}}$ is dense in X,
• null controllable in time t if the range of ${\displaystyle \Phi _{t}}$ includes the range of ${\displaystyle {\rm {e}}^{At}}$.

Observability

As in the finite-dimensional case, observability is the dual notion of controllability. In the infinite-dimensional case there are several different notions of observability which in the finite-dimensional case coincide. The three most important ones are:

• Exact observability (also known as continuous observability),
• Approximate observability,
• Final state observability.

Observability in discrete-time

An important role is played by the maps ${\displaystyle \Psi _{n}}$ which map X into the space of all Y valued sequences and are given by ${\displaystyle (\Psi _{n}x)_{k}=CA^{k}x}$ if k ≤ n and zero if k > n. The interpretation is that ${\displaystyle \Psi _{n}x}$ is the truncated output with initial condition x and control zero. The system is called

• exactly observable in time n if there exists a k_n > 0 such that ${\displaystyle \|\Psi _{n}x\|\geq k_{n}\|x\|}$ for all x ∈ X,
• approximately observable in time n if ${\displaystyle \Psi _{n}}$ is injective,
• final state observable in time n if there exists a k_n > 0 such that ${\displaystyle \|\Psi _{n}x\|\geq k_{n}\|A^{n}x\|}$ for all x ∈ X.

Observability in continuous-time

In observability of continuous-time systems the map ${\displaystyle \Psi _{t}}$ given by ${\displaystyle (\Psi _{t})(s)=C{\rm {e}}^{As}x}$ for s∈[0,t] and zero for s>t plays the role that ${\displaystyle \Psi _{n}}$ plays in discrete-time. However, the space of functions to which this operator maps now influences the definition. The usual choice is L^2(0, ∞, Y), the space of (equivalence classes of) Y-valued square integrable functions on the interval (0,∞), but other choices such as L^1(0, ∞, Y) are possible. The different observability notions can be defined once the co-domain of ${\displaystyle \Psi _{t}}$ is chosen. The system is called[9]

• exactly observable in time t if there exists a k_t > 0 such that ${\displaystyle \|\Psi _{t}x\|\geq k_{t}\|x\|}$ for all x ∈ X,
• approximately observable in time t if ${\displaystyle \Psi _{t}}$ is injective,
• final state observable in time t if there exists a k_t > 0 such that ${\displaystyle \|\Psi _{t}x\|\geq k_{t}\|{\rm {e}}^{At}x\|}$ for all x ∈ X.

Duality between controllability and observability

As in the finite-dimensional case, controllability and observability are dual concepts (at least when for the domain of ${\displaystyle \Phi }$ and the co-domain of ${\displaystyle \Psi }$ the usual L^2 choice is made). The correspondence under duality of the different concepts is:[10]

• Exact controllability ↔ Exact observability,
• Approximate controllability ↔ Approximate observability,
• Null controllability ↔ Final state observability.

See also

• Control theory
• State space (controls)

Notes

1. Curtain and Zwart
2. Curtain and Zwart Example 2.2.4
3. Curtain and Zwart Theorem 2.4.6
4. This is the mathematical convention, engineers seem to prefer transfer functions to be holomorphic at infinity; this is achieved by replacing z by 1/z
5. Curtain and Zwart Lemma 4.3.6
6. Staffans Theorem 4.6.7
7. Curtain and Zwart Example 4.3.13
8. Tucsnak Definition 11.1.1
9. Tucsnak Definition 6.1.1
10. Tucsnak Theorem 11.2.1

References

• Curtain, Ruth; Zwart, Hans (1995), An Introduction to Infinite-Dimensional Linear Systems theory, Springer
• Tucsnak, Marius; Weiss, George (2009), Observation and Control for Operator Semigroups, Birkhauser
• Staffans, Olof (2005), Well-posed linear systems, Cambridge University Press
• Luo, Zheng-Hua; Guo, Bao-Zhu; Morgul, Omer (1999), Stability and Stabilization of Infinite Dimensional Systems with Applications, Springer
• Lasiecka, Irena; Triggiani, Roberto (2000), Control Theory for Partial Differential Equations, Cambridge University Press
• Bensoussan, Alain; Da Prato, Giuseppe; Delfour, Michel; Mitter, Sanjoy (2007), Representation and Control of Infinite Dimensional Systems (second ed.), Birkhauser