This paper develops a new conservative functional model for operator-valued Schur functions on the right-half plane, extending classical disk models to continuous-time systems with unbounded operators.
Contribution
It introduces a de Branges-Rovnyak model directly in the right-half plane, incorporating unbounded connecting operators and non-invertible intertwinements for enhanced uniqueness.
Findings
01
Model exhibits structure absent in disk setting
02
Handles unbounded connecting operators in continuous-time systems
03
Strengthens classical uniqueness results
Abstract
We present a solution of the operator-valued Schur-function realization problem on the right-half plane by developing the corresponding de Branges-Rovnyak canonical conservative simple functional model. This model corresponds to the closely connected unitary model in the disk setting, but we work the theory out directly in the right-half plane, which allows us to exhibit structure which is absent in the disk case. A main feature of the study is that the connecting operator is unbounded, and so we need to make use of the theory of well-posed continuous-time systems. In order to strengthen the classical uniqueness result (which states uniqueness up to unitary similarity), we introduce non-invertible intertwinements of system nodes.
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsHolomorphic and Operator Theory · Spectral Theory in Mathematical Physics · Differential Equations and Boundary Problems
Full text
A conservative de Branges-Rovnyak func- tional model for operator Schur functions on C+
We present a solution of the operator-valued Schur-function realization problem on the right-half plane by developing the corresponding de Branges-Rovnyak canonical conservative simple functional model. This model corresponds to the closely connected unitary model in the disk setting, but we work the theory out directly in the right-half plane, which allows us to exhibit structure which is absent in the disk case. A main feature of the study is that the connecting operator is unbounded, and so we need to make use of the theory of well-posed continuous-time systems. In order to strengthen the classical uniqueness result (which states uniqueness up to unitary similarity), we introduce non-invertible intertwinements of system nodes.
Key words and phrases:
Schur function; continuous time; linear system; right half-plane; functional model; de Branges-Rovnyak; realization; reproducing kernel
The classical unitary realization result of de Branges and Rovnyak for Schur functions on the complex unit disk is the following: Let U and Y be separable Hilbert spaces and let ϕ be an operator Schur function on D, i.e., ϕ is analytic with ϕ(z)∈L(U;Y) a contraction for all z∈D. Then the following kernel function on D×D, whose values are bounded linear operators on [YU], is positive:
[TABLE]
Denoting its reproducing kernel Hilbert space (RKHS) by Hs, we obtain that the following operator [AsCsBsDs]:[HsU]→[HsY] is unitary:
[TABLE]
This is the classical de Branges-Rovnyak closely connected unitary functional model for ϕ. Indeed, it is unitary and closely connected (the disk version of the concept of simplicity defined in Def. 2.2 below), its transfer function Gs(z)=zCs(1−zAs)−1Bs+Ds coincides with ϕ on D, and conversely, every closely connected unitary realization of ϕ is obtained from (1.2) by means of a unitary change of state variable.
In this paper, we develop a version of this result in the right-half-plane setting, which requires that we use well-posed systems theory in continuous time. In particular, the analogue of [AsCsBsDs] in (1.2) is unbounded. Our main results can be summarized in the following simplified form which is completely analogous to the above disk case if one restricts to μ∗=μ and λ∗=λ:
Theorem 1.1**.**
Let φ be an operator Schur function on the complex right-half plane C+. Then the kernel function
[TABLE]
μ,μ∗,λ,λ∗∈C+, is positive; denote its associated RKHS by Hs.
The following unbounded linear operator from [HsU] into [HsY] is a simple, conservative system node (definitions in §I.3):
[TABLE]
with domain consisting of all [xu] for which this makes sense:
[TABLE]
The transfer function of [A&BC&D]s is φ, and conversely every simple, conservative realization of φ is unitarily similar to [A&BC&D]s; see Def. 2.7 below.
This paper is a direct continuation of [BKSZ15] published earlier in this journal. We refer to that paper as “Part I” and assume that the reader is familiar with it. In Part I, the research is placed in its context and detailed background on passive system nodes is presented. Results from Part I will be referenced using a capital ‘I’; e.g. Thm. I.5.1.3 refers to item 3 of [BKSZ15, Thm. 5.1]. In a certain sense, the conservative model is a coupling of the two semi-conservative models in Part I, but working with the conservative model is easier than working with those in Part I. Indeed, the conservative model has the same structure as its adjoint, and hence it combines all the good properties of the semi-conservative realizations.
Investigations closely related to that reported here have been undertaken before, starting from the work [dBR66a, dBR66b] of de Branges and Rovnyak; see [Bro78] for a nice historic overview of work on the disk case up to that point. For a good monograph on the disk case, see [ADRdS97]. The first results in the right-half-plane setting are in [AN96]; here Arov and Nudelman used a linear fractional transformation to reduce the half-plane case to the disk case. Most of the more recent publications on half-plane functional models also employ this so-called Cayley transformation, but in the present paper we work the details out directly in the half-plane setting, in order to expose detail that is invisible in the disk setting.
Adamjan and Arov [AA66] showed how to embed the Sz.-Nagy-Foiaş model into a suitably more general version of the Lax-Phillips scattering picture (for the discrete-time setting); much later
Nikolski-Vasyunin [NV86, NV89, NV98] refined this analysis by doing
such an embedding also for a Pavlov model and a suitably modified
version of the de Branges-Rovnyak model. Continuous-time versions of this analysis are also of interest, and we plan to investigate this in a forthcoming publication.
In [BS06] an (implicit) “lurking isometry” argument and Cayley transformations were used to obtain the existence of a conservative realization for any operator Schur function on the disk or the right-half plane. The realization that we describe in the present paper is a more explicit alternative to the realization constructed in [BS06]. The results of [BS06] have been extended to a multi-variable case in [BKV15], and we expect also the present results to have natural extensions to various multi-variable settings.
The continuous-time conservative realization has been studied in the state/signal framework developed by Arov and Staffans, too, in [AKS11]. Here the central idea is to consider in H2(C+;W) the graph of the Toeplitz operator Tφ with symbol φ, where W is a Kreĭn space, without assuming any particular partition W=U∔Y into an input space U and output space Y. In this setting, the action of the realization is a pure shift on the appropriate state space and projections onto input and output components are avoided, which leads to cleaner formulas and intuitively more transparent results; see [AKS11, AS09, AS10] for details. Again, in the present work the objective is to obtain as explicit formulas as possible for the input/state/output setting.
Finally, we mention that a closely related realization of a Nevanlinna family (corresponding to an impedance-passive setting rather than to the present scattering-passive setting) in terms of a boundary relation has been worked out in [BHdS08] (or see [BHdS09] for a more elaborate version).
The paper is laid out as follows: In §2 we briefly present some additional background on conservative and simple system nodes that is needed in the present paper, with the auxiliary proofs on non-invertible intertwinements postponed to Appendix A. The conservative model is introduced in §4 after its state space has been constructed in §3. In §5, we present an explicit identification of the extrapolation space and calculate the (unbounded) control operator of the conservative model. The paper is concluded in §6, where we exhibit the relationship with the classical de Branges-Rovnyak model (1.1)–(1.2).
2. More on passive system nodes
We sharpen the uniqueness results in Part I and also recall a few additional concepts from continuous-time systems theory that are needed in the present paper. The discussion that follows uses the definitions of a passive system node, its main operator A, control operator B, observation operator C, and transfer function given in §I.3.
We recall from Def. I.3.6 that [A&BC&D]:[XU]⊃dom([A&BC&D])→[XY] is called scattering dissipative if for all [xu]∈dom([A&BC&D]):
[TABLE]
By Def. I.3.7 every passive system node is scattering dissipative. By the following rather obvious consequence of [Sta13, Thm. 2.5], every passive system node is even maximal scattering dissipative:
Lemma 2.1**.**
A passive system node has no scattering-dissipative proper extension.
Next, we define the concept of simplicity of a continuous-time system. In the following definition, A−1 denotes the unique extension of A to a closed operator on the extrapolation space induced by A, and A−1d denotes the analogous extension of A∗ to its extrapolation space – this latter operator was denoted by A^{*}\big{|}_{\mathcal{X}} in Part I.
Definition 2.2**.**
A passive system node [A&BC&D] with state space X is simple if
[TABLE]
Comparing simplicity to the notions of controllability and observability in §I.3, one observes that every controllable and every observable passive system is simple; take either γ=0 or ν=0 in (2.1).
The equation (I.3.6), which is valid for every system node, plays an important role in the theory of de Branges-Rovnyak models on C+, e.g., in the proof of Thm. I.4.3 (or its further development Thm. 2.8 below). For conservative systems, we have the additional equality (2.3) below:
Lemma 2.3**.**
Let [A&BC&D] be a conservative system node. The (in general unbounded) adjoint of [A&BC&D] is the system node
[TABLE]
For every λ∈C+ and γ∈y,
[TABLE]
Sometimes we write (2.2) in the “time-flow inverse” form
2.1. The passive input, output and past/future maps
The exposition and terminology of this section loosely follow [AS09, §§5 – 6].
Theorem 2.4**.**
Let [A&BC&D] be a passive system node and denote by Ho and Hc the Hilbert spaces with reproducing kernels Ko and Kc in (I.1.17), respectively.
The following (passive frequency-domain) output map is a contraction from X into Ho:
[TABLE]
Moreover, the mapping
[TABLE]
extends by linearity and operator closure to a contraction mapping Hc into X.
The theorem can be seen as a consequence of [Sta05, Thm. 11.1.6], but the connection requires some explanations, and so we include a proof formulated in the present setup for reading convenience. At the end of the proof we need the following notation which is familiar from Part I:
[TABLE]
Also, we introduce the notation Tφ for the usual Toeplitz operator Tφ with symbol φ∈L∞(iR;U,Y):
[TABLE]
where Mφ:L2(iR;U)→L2(iR;Y) is the operator multiplying by φ and P+ is the orthogonal projection of L2(iR;Y) onto H2(C+;Y). In our case, we always have φ∈S(C+;U,Y)⊂H∞(C+;U,Y), so that Tφu=Mφu for all u∈H2(C+;U,Y); indeed, in Part I, we used the somewhat less precise notation Mφ rather than Tφ.
See [Sta05, §11.1] for background and more details on this proof.
Let (u,x,y) be a stable classical trajectory of the system node [A&BC&D], i.e., u∈L2(R+;U)∩C(R+;U), x∈C1(C+;X), y∈C(R+;Y), and
[TABLE]
By Defs. I.3.6–7, we then have for all t≥0 that
[TABLE]
and integrating this from [math] to T≥0, we obtain
[TABLE]
Letting T→+∞, we obtain that y∈L2(R+;Y) and ∥x(0)∥2≥∥y∥L2(R+;Y)2−∥u∥L2(R+;U)2; moreover
[TABLE]
and so ∥x(T)∥2 is bounded. Thus we may take Laplace transforms, obtaining that ∥x(0)∥2≥∥y∥H2(C+;Y)2−∥u∥H2(C+;U)2 and
[TABLE]
Hence, y=Cx(0)+Tφu and the operator [CTφ] is a contraction from
[TABLE]
(as a subset of [XH2(C+;U)]) into H2(C+;Y). By [Sta05, 4.6.11], the set (2.7) contains [dom(A)H01(R+;U)], where H01(R+;U) is the first order Sobolev space of U-valued functions u with the additional restriction u(0)=0, and [dom(A)H01(R+;U)] is dense in [XH2(C+;U)] because A generates a contraction semigroup on X and the Laplace transformation is unitary.
Hence, [CTφ][C⋄Tφ∗]≤1 on H2(C+;Y), where C⋄ denotes the adjoint of C calculated with respect to the inner product in H2(C+;Y) rather than with respect to the inner product in Ho (in which case we would have written C∗). Thus CC⋄≤1−TφTφ∗ and by Douglas’ lemma there exists some contraction C:X→ker(1−TφTφ∗)⊥ such that C=(1−TφTφ∗)1/2C. This implies that C is a contraction from X into Ho, because for every x∈X:
[TABLE]
Consider now the output map Cd of the passive system node [A&BC&D]∗. From (I.1.17) it follows that Ho constructed with φ is the same as Hc constructed with φ. Thus Cd=B∗ is a contraction from X into Hc, and so B is a contraction from Hc into X.
∎
From now on we let B be the extension of (2.6) by linearity and continuity and we call it the (passive frequency-domain) input map. Please note that the input and output maps of the dual system [A&BC&D]∗ are Cd∈L(X;Hc) and Bd∈L(Ho;X) given by
[TABLE]
respectively; see Prop. I.3.10. We shall later make extensive use of the unobservable subspaceU and the approximately reachable subspaceR of [A&BC&D], which are given by
[TABLE]
(closure in X). Also note that U⊥ is the approximately reachable subspace and R⊥ the unobservable subspace of the dual system node [A&BC&D]∗. Therefore, we denote U⊥=:R† and R⊥=:U†.
Definition 2.5**.**
For a passive system node [A&BC&D], we call the contraction Γ:=CB:Hc→Ho the (frequency domain) past/future-map of [A&BC&D].
We first give the action of the past/future map on kernel functions.
Proposition 2.6**.**
The transfer function φ uniquely determines Γ via
[TABLE]
The adjoint of Γ has the following action on kernel functions:
[TABLE]
Denoting the past/future map of φ by Γφ:Hc→Ho, and similar for Γφ:Ho→Hc, we have
[TABLE]
As a consequence, all passive realizations of the same transfer function have the same past/future map.
The equation (2.10) follows from the following computation:
[TABLE]
Then (2.11) follows from (valid for all λ∗,λ∈C+, ν∈U, and γ∈Y):
[TABLE]
Finally, by (2.10) and (2.11), Γφ∗eo(λ)∗γ=Γφeo(λ)∗γ for all λ∈C+ and γ∈Y. Considering linear combinations of kernel functions and extending by continuity, we obtain (2.12).
∎
This operator played a very important role in the more explicit representation of the energy-preserving controllable model [A&BC&D]c and its extrapolation space in Part I. In the present paper, the operator Γ plays an even more crucial role, already in the proof that the reproducing kernel defining the state space of the simple conservative model is positive.
Example*.*
From (2.15) and (2.5), we obtain that the controllable energy-preserving model [A&BC&D]c in §I.4 has Cc=Γ. Moreover, comparing (2.6) with (I.4.7), it becomes evident that Bc=1. Similarly, (I.5.8) implies that the observable co-energy-preserving model [A&BC&D]o has input map Bo=Γ and output map Co=1Ho.
Let e+(μ) and e−(μ∗) denote point-evaluation of functions in H2(C+), at μ∈C+, and H2(C−), at μ∗∈C−, respectively. The right-hand side of (2.10) equals the action of the Hankel operator \mathrm{H}_{\varphi}:=P_{H^{2}({{\mathbb{C}}^{+}};\mathcal{Y})}\,M_{\varphi}\big{|}_{H^{2}({{\mathbb{C}}^{-}};\mathcal{U})} on a kernel function of H2(C−;U), namely
[TABLE]
Indeed, for all fixed parameters λ∗∈C+, for all fixed vectors ν∈U, and for almost all values of the variable μ∈iR:
[TABLE]
where the first term is in H2(C+;Y) and the second term is in
H2(C−;Y); from here we deduce that
[TABLE]
where e−(−λ∗)∗ν=k−(⋅,−λ∗)ν. Then, using the injections ιc:Hc→H2(C+;U) and ιo:Ho→H2(C+;Y), we have for all μ,λ∗∈C+ that
[TABLE]
see Thm. I.2.4. Taking linear combinations and closing, we get
[TABLE]
with the reflection being \reflectboxR:e+(λ∗)∗ν↦e−(−λ∗)∗ν, λ∗∈C+, ν∈U, extended by linearity and continuity to a unitary operator H2(C+;U)→H2(C−;U).
The contractions Bιc∗:H2(C+;U)→X and ιoC:X→H2(C+;Y) factorizing the Hankel operator are also sometimes referred to as (frequency-domain) input and output maps, but here we refer to B and C by these names. From a systems-theory point of view, it would perhaps be more natural to take the state space Hc of the controllable energy-preserving model to be a subspace of H2(C−;U) rather than a subspace of H2(C+;U), since one in time domain often considers input signals in past time R− rather than in future time R+. In particular, the reflection is then absent in (2.17).
2.2. A generalized uniqueness result
We see that B in the previous subsection is precisely the unitary similarity operator Δ in Thm. I.4.3 and C is the adjoint of Δ in Thm. I.5.2. We now proceed to obtain improvements on these uniqueness results in Part I. First we need to relax the notion of unitary similarity from Thm. I.4.3. Please note that this subsection first considers general system nodes, not only passive ones.
Definition 2.7**.**
Let [A&BC&D]0 and [A&BC&D]1 be two system nodes with state spaces X0 and X1, respectively, and the same input spaces U and output spaces Y. Let E map X0 linearly and boundedly into X1.
We say that Eintertwines [A&BC&D]0 with [A&BC&D]1 if
[TABLE]
If E is a contraction (an isometry), then we call the intertwinement contractive (isometric), and if E is unitary then we say that [A&BC&D]0 and [A&BC&D]1 are unitarily similar.
Alternative characterizations of intertwinement and more detail can be found in Appendix A. Paralleling [AKS11, Thms. 8.4 and 9.5], we have the following uniqueness result which is stronger than Thms. I.4.3 and I.5.2:
Theorem 2.8**.**
The following statements hold for every passive system node [A&BC&D] with transfer function φ, input map B, and output map C:
(1)
The input map B intertwines the energy-preserving model [A&BC&D]c for φ contractively with [A&BC&D]. This intertwinement is isometric if [A&BC&D] is energy preserving. The intertwinement B has range dense in X if and only if [A&BC&D] is controllable. Moreover, B is unitary if and only if [A&BC&D] is controllable and energy-preserving.
2. (2)
The output map C intertwines [A&BC&D] contractively with the co-energy-preserving observable model [A&BC&D]o for φ. This intertwinement is co-isometric if [A&BC&D] is co-energy preserving and C is injective if and only if [A&BC&D] is observable. Furthermore, C is unitary if and only if [A&BC&D] is observable and co-energy preserving.
3. (3)
In particular, Bo=Γ intertwines [A&BC&D]c with [A&BC&D]o. Hence, for all [xu]∈dom([A&BC&D]c) and [zy]=[A&BC&D]c[xu]:
[TABLE]
Proof.
We begin with statement one. In Thm. I.4.3, Δ=B and the intertwinement part of the proof goes through even if this operator is only continuous. The proof that B is isometric if [A&BC&D]c is energy preserving is also the same as in Thm. I.4.3. The connection between controllability and dense range is immediate from (2.9). That [A&BC&D] is energy preserving in the unitary case follows from Thm. A.3.3.
We obtain statement two by duality: The input map of [A&BC&D]∗ is C∗ by (2.8), and this operator intertwines the energy-preserving model [A&BC&D]o∗ for φ(μ)=φ(μ)∗ contractively with [A&BC&D]∗; see the introduction to §I.5. Using Lemma A.2, we obtain that C intertwines [A&BC&D] with [A&BC&D]o. The rest of the claim is immediate from the definitions of co-isometry and co-energy-preserving system node.
By Ex. Example, Bo=Γ, and then assertion one with Def. 2.7 gives
The first step in the development is to construct a positive 2×2-block kernel function using only φ, whose reproducing kernel Hilbert space will be the state space of the realization [A&BC&D]s. As in [AS07, §8], we develop the theory using a four-variable kernel rather than the standard two-variable kernel, hoping to make visible how the observable co-energy-preserving and controllable energy-preserving functional models are combined into the conservative simple model.
We begin with a general result on how a RKHS can arise as the range of a multiplication operator:
Lemma 3.1**.**
Let Ω be a point set, X, Y be
Hilbert spaces and let H:Ω→L(X,Y) be an
operator-valued function. Define a subspace HH of the linear space of Y-valued functions on Ω by
[TABLE]
where MH is the multiplication operator MH:x↦H(⋅)x from X to HH.
Then HH is a reproducing kernel Hilbert space
with reproducing kernel KH(z,w):=H(z)H(w)∗ and MH maps ker(MH)⊥ unitarily onto HH.
Proof.
By Thm. I.1.1, K is the reproducing kernel of some uniquely determined Hilbert space of functions. Let H(⋅)x∈HH and let w∈Ω.
For y∈Y we note that
[TABLE]
Noting that ker(MH)={ξ∈X:H(z)ξ=0 for all z∈Ω}, we get H(w)∗y∈ker(MH)⊥ and
⟨x,H(w)∗y⟩X=⟨x′,H(w)∗y⟩X
where x′=P(KerMH)⊥x. By
construction MH is an isometry from ker(MH)⊥ onto HH. Hence the above calculation
continues as
[TABLE]
and we conclude that K(z,w)=H(z)H(w)∗ works as the RK of HH.
∎
Taking Ω=C+ and either H(μ)=C(μ−A)−1 or H(μ)=B∗(μ−A∗)−1, we get the following interesting special cases (recall Thm. 2.4 and the text around (2.8)):
Corollary 3.2**.**
*Assume that [A&BC&D] is a passive system node with input/state/
output spaces (U,X,Y), input map B, and output map C.*
If {\rm im}\bigl{(}\mathfrak{C}\bigr{)} is equipped with the lifted norm \|{\mathfrak{C}}x\|_{{\rm im}\bigl{(}\mathfrak{C}\bigr{)}}=\|P_{\mathfrak{R}^{\dagger}}x\|_{{\mathcal{X}}} then C maps R† unitarily onto {\rm im}\bigl{(}\mathfrak{C}\bigr{)} which is a RKHS HC with reproducing kernel
[TABLE]
Similarly, B∗ is a unitary identification of R with HB∗, where
[TABLE]
is the RKHS with reproducing kernel
[TABLE]
The following result which draws some inspiration from [ADRdS97, Thm. 2.1.2] determines the kernel function (1.3) needed to define the state space of the conservative realization.
Proposition 3.3**.**
Let [A&BC&D] be an arbitrary system node with transfer function φ.
(1)
If [A&BC&D] is co-energy preserving then the reproducing kernel Ko defining the state space for the observable model [A&BC&D]o factorizes as
[TABLE]
and the output map C maps R† unitarily onto Ho.
2. (2)
Set as before φ(μ):=φ(μ)∗, μ∈C+. If [A&BC&D] is energy preserving then Kc associated to the controllable model [A&BC&D]c factorizes as
[TABLE]
μ∗,λ∗∈C+, and the input map B maps Hc unitarily onto R.
3. (3)
Define
[TABLE]
If [A&BC&D] is conservative, then Ks defined in (1.3) factorizes into
[TABLE]
Clearly, H(μ)=eo(μ)C and G(μ∗)=ec(μ∗)B∗. The previous result also extends (I.4.73) and the first formula in Prop. I.5.8, since the realization of φ is arbitrary – only suitable energy properties are assumed.
Assuming that φ∈S(C+;U,Y), the kernel Ks has removable singularities at μ∗=λ and μ=λ∗. When we remove these singularities by continuity, the kernel becomes holomorphic with its values being bounded operators on [YU]. In the sequel we ignore removable singularities, assuming that they have been removed.
Proof.
Let (U,X,Y) denote the input/state/output spaces of [A&BC&D]. We begin by proving (3.2). By (I.3.6), every system node satisfies
[TABLE]
For [A&BC&D] energy preserving, (I.3.14) gives that for all λ∗,μ∗∈C+, ν,v∈U:
[TABLE]
This implies that
[TABLE]
i.e., that (3.2) holds. That (3.1) holds for a co-energy preservation system node follows by applying (3.2) to the energy-preserving system node [A&BC&D]∗; recall that the transfer function of this dual system is φ and that (φ)=φ.
A conservative system is by definition both energy-preserving and co-energy preserving, and so (3.1) and (3.2) both hold. Moreover, by (2.13) every system node satisfies
[TABLE]
and this implies
[TABLE]
To establish the unitary of C from R† onto Ho in assertion 1, we observe that KC=Ko in the co-energy preserving case. This implies that HC=Ho and then unitarity follows from Cor. 3.1. Analogously, B∗ maps R unitarily onto Hc, which implies that B is an isometry into X with range R.
∎
Alternatively, (3.1) can be inferred from Thm. I.5.2 and (3.2) can also be seen as a consequence of Thm. I.4.3. The existence of a conservative realization of an arbitrary operator Schur function on C+ has been proved in, e.g., [AN96, BS06]. Formula (3.4) provides a Kolmogorov factorization of Ks, which proves that Ks is positive, hence the reproducing kernel of a Hilbert space. In order to keep the present article (together with Part I) self-contained, we provide a short direct proof of the positivity in the style of Part I and [AS10, pp. 3321–3323].
Lemma 3.4**.**
For φ∈S(C+;U,Y), the function Ks in (1.3) is a positive kernel on (C+×C+)2, i.e., for all gk∈Y, vk∈U, and ωk,ζk∈C+, k=1,…,N, we have
[TABLE]
Let Γ be the past/future map determined by φ in (2.10) and let e(μ,μ∗)=[e+(μ)00e+(μ∗)] be point-evaluation of functions in [H2(C+;Y)H2(C+;U)]. Then the kernel can be factorized as
[TABLE]
where I=[ιo00ιc]:[HoHc]→[H2(C+;Y)H2(C+;U)] is the injection and [1Γ∗Γ1] is positive semidefinite on [HoHc].
Proof.
By Thm. I.2.4, we have I∗=[1−TφTφ∗001−TφTφ∗], and then Prop. 2.6 and Lem. I.2.1.3 give
[TABLE]
which proves (3.5). Furthermore, [1Γ∗Γ1] is positive semidefinite on [HoHc] due to the contractivity of Γ. Now the positivity follows upon observing that
[TABLE]
∎
In order to fit into standard reproducing kernel Hilbert space (RKHS) theory, we can alternatively interpret Ks(μ,μ∗,λ,λ∗) as a kernel function of two variables μ:=(μ,μ∗) and λ:=(λ,λ∗), both in C+×C+. Then our positive kernel function has the special 2×2-block form
[TABLE]
Elements of the corresponding RKHS, which we denote by Hs, are densely spanned by the kernel functions
Ks(⋅,λ)[yu], where λ sweeps C+×C+ and [yu] sweeps Y⊕U. Note that
each such function is a column [fg](μ) of the form
[TABLE]
where f and g are analytic on C+; therefore this property continues to hold for
all elements of Hs. By standard RKHS theory, the reproducing property is
[TABLE]
Taking ν=0 and setting (the first column of Ks(μ,λ))
[TABLE]
then gives
[TABLE]
Similarly, taking γ=0 gives
[TABLE]
and the general case (3.7) is the superposition of these two.
Theorem 3.5**.**
Let W denote the positive semidefinite square root of [1Γ∗Γ1]∈L([HoHc]). Then:
(1)
We have \mathcal{H}_{s}={\rm im}\bigl{(}W\bigr{)}\subset\left[\begin{smallmatrix}\mathcal{H}_{o}\\
\mathcal{H}_{c}\end{smallmatrix}\right] with the lifted norm
[TABLE]
2. (2)
The operators [Γ1]:Hc→Hs and [1Γ∗]:Ho→Hs are isometric.
3. (3)
Setting Rs†:=[1Γ∗]Ho and Rs:=[Γ1]Hc, both with the norm of Hs, we obtain
[TABLE]
with the first embedding dense and the second continuous. The adjoint of the injection ι:Hs→[HoHc] is
[TABLE]
and the reproducing kernel Ks of Hs has the representation
[TABLE]
4. (4)
The subspaces Rs† and Rs are closed, both in Hs and in [HoHc].
5. (5)
The following maps are co-isometries from Hs onto Ho and Hc, respectively:
[TABLE]
The initial subspace of Π1 is Rs† and the initial subspace of Π2 is Rs. The operators \pi_{1}:=\Pi_{1}\big{|}_{\mathfrak{R}_{s}^{\dagger}}:\mathfrak{R}_{s}^{\dagger}\to\mathcal{H}_{o} and \pi_{2}:=\Pi_{2}\big{|}_{\mathfrak{R}_{s}}:\mathfrak{R}_{s}\to\mathcal{H}_{c} are unitary.
6. (6)
Denote by PRs and PRs† the orthogonal projections in Hs onto Rs and Rs†, respectively, and let Us:=(Rs†)⊥ and Us†:=(Rs)⊥. Then
[TABLE]
The spaces Us and Rs will turn out to be the unobservable and approximately reachable subspaces of the conservative simple model, respectively.
Modifying the proof of Thm. I.2.4 slightly, we obtain (3.8) and (3.10). Due to (3.8), [Γ1] is isometric:
[TABLE]
and the isometricity of [1Γ∗] is proved the same way.
Eq. (3.11) is only a restatement of (3.6). Trivially, W2[HoHc]⊂W[HoHc]⊂[HoHc], and this establishes (3.9), where the first embedding is dense, because the reproducing kernels of Hs lie in Rs†+Rs by (3.11). Moreover, the norm of [1Γ∗Γ1] as an operator on [HoHc] is at most 2: Using that Γ is a contraction, Cauchy-Schwarz, and completion of squares, we obtain
[TABLE]
For all [fg]=Wh with h∈[HoHc]⊖ker(W) it then holds that
[TABLE]
i.e., ι is continuous with norm at most 2. Claims one to three are proved.
Since Γ is bounded with a closed domain, it follows immediately that Rs is closed in [HoHc], and by the isometricity of [Γ1], Rs is closed also in Hs. Furthermore, it follows from Rs†+Rs⊂Hs that Π2 and π2 are onto Hc:
[TABLE]
An analogous argument shows that Rs† is closed and Π1Hs=Ho.
The operator \Pi_{2}\big{|}_{\mathfrak{R}_{s}}=\pi_{2}:\mathfrak{R}_{s}\to\mathcal{H}_{c} is unitary by (3.12) and (3.14). Due to (3.11) and the isometricity of [Γ1], the space Rs is the closed linear span of the kernel functions
[TABLE]
(λ∈C+ is insignificant) and this implies that ker(Π2)=Hs⊖Rs:
[TABLE]
Splitting Hs=[RsRs⊥], we thus obtain Π2=[π20], and furthermore, by the unitarity of π2:
[TABLE]
Hence, the operator Π2 is a co-isometry with initial space Rs (and final space Hc).
The unitarity of π2 implies that π2∗=π2−1 and this operator equals [Γ1], because π2[Γ1]x=x for all x∈Hc; premultiply by π2∗. Then the formula Γ=Π1π2∗ trivially follows. Moreover, PRs=π2∗Π2, because (π2∗Π2)2=π2∗Π2, {\rm im}\bigl{(}\pi_{2}^{*}\,\Pi_{2}\bigr{)}=\pi_{2}^{*}\mathcal{H}_{c}=\mathfrak{R}_{s}, and ker(π2∗Π2)=ker(Π2)=Hs⊖Rs. Then
[TABLE]
The claims on Π1, π1, Γ∗, PRs†, and PUs are proved in the same way.
∎
We will need the following extension of Prop. I.2.6:
Corollary 3.6**.**
Every [fg]∈Hs satisfies f(μ)→0 in Y as Reμ→+∞ and g(μ∗)→0 in U as Reμ∗→+∞. More precisely, for all [fg]∈Hs and μ,μ∗∈C+:
[TABLE]
Proof.
By Prop. I.2.6 and Thm. I.2.4.2, for all f∈Ho and g∈Hc:
[TABLE]
Restricting to Hs and combining this with (3.13) completes the argument for f, and g is handled the same way.
∎
We end the section with an analogue of Thm. 2.4. It is needed for our uniqueness result for the conservative simple model, which is a variation on Thm. 2.8. Inspired by [AKS11, §10], we define the (frequency-domain) bilateral input map of a passive system [A&BC&D] with state space X as the mapping
[TABLE]
as a first step defined on the dense subspace [1Γ∗Γ1][HoHc] of Hs with range in X. We shall in a moment prove that this mapping can be extended to a contraction Hs→X, and its contractive adjoint is called the (frequency-domain) bilateral output map
[TABLE]
Theorem 3.7**.**
Let [A&BC&D] be a passive system node with transfer function φ.
(1)
The bilateral input map Bbil in (3.18) is a contraction Hs→X and ran(Bbil)=R+R†, where the closure is in X. Hence, [A&BC&D] is simple if and only if Bbil has dense range which holds if and only if Cbil is injective.
2. (2)
Denoting the inverse of the injection ι:Hs→[HoHc] with domain ιH by ι−1, the bilateral output map can be written more explicitly as
[TABLE]
i.e, Cbil is [CB∗] with range space Hs rather than [HoHc].
3. (3)
The bilateral input map has the following action on the kernel functions of Hs:
[TABLE]
λ,λ∗∈C+,γ∈Y,ν∈U, and in the notation of Prop. 3.3.3:
[TABLE]
4. (4)
The operator Cbil maps R+R† unitarily onto HCbil, the RKHS with reproducing kernel
[TABLE]
We have the alternative characterization
[TABLE]
5. (5)
Assume that [A&BC&D] is conservative. Then HCbil=Hs and Bbil maps Hs isometrically into X, unitarily onto R+R†. Moreover, \mathfrak{C}^{*}\Gamma\mathfrak{B}^{*}\big{|}_{\mathfrak{R}}=P_{\mathfrak{R}^{\dagger}}\big{|}_{\mathfrak{R}}
Proof.
For every h∈[HoHc], we have (using (3.18), Γ=CB with C and B contractive, and (3.8)):
[TABLE]
hence Bbil is contractive on {\rm im}\bigl{(}\left[\begin{smallmatrix}1&\Gamma\\
\Gamma^{*}&1\end{smallmatrix}\right]\bigr{)} which is dense in Hs. By (3.10) and (3.18) it holds that Bbilι∗=[C∗B]:[HoHc]→X, and this implies that
compare this to Def. 2.2 to obtain the characterizations of simplicity. Finally, (2.8) and (3.3) give (3.21); the factor ι−1 emphasizes the fact that Cbil maps into Hs rather than [HoHc]. This completes the proof of assertions one to three.
Assertion four follows from Lem. 3.1 upon observing that (2.5), (2.8), and (3.3) imply that
[TABLE]
and that by the above, ker(Cbil)⊥=R+R†. For the rest of the proof, we assume that [A&BC&D] is conservative. Then Prop. 3.3.3 gives that KCbil=Ks and by assertion 4, Cbil maps R+R† (which is isometrically contained in X) unitarily onto Hs; then Bbil=Cbil∗ maps Hs isometrically into X, with {\rm im}\bigl{(}\mathfrak{B}_{bil}\bigr{)}=\overline{\mathfrak{R}+\mathfrak{R}^{\dagger}} by the above; now Bbil has closed range because it is isometric with a closed domain. By Prop. 3.3, C∗ and B are both isometric into X; hence C∗C=PR† and BB∗=PR. Combining this with the Def. 2.5 of Γ gives \mathfrak{C}^{*}\Gamma\mathfrak{B}^{*}\big{|}_{\mathfrak{R}}=P_{\mathfrak{R}^{\dagger}}\big{|}_{\mathfrak{R}}.
∎
We remark that \mathfrak{C}^{*}\Gamma\mathfrak{B}^{*}\big{|}_{\mathfrak{R}}=P_{\mathfrak{R}^{\dagger}}\big{|}_{\mathfrak{R}} for a conservative system node means that \Gamma=P_{\mathfrak{R}^{\dagger}}\big{|}_{\mathfrak{R}} if we make the unitary identification of y∈Ho with C∗y and, similarly, we identify x∈R with B∗x.
4. The conservative simple functional model
We construct the conservative simple realization in the following way; cf. Lem. I.4.1:
Proposition 4.1**.**
Let φ∈S(C+;U,Y) and let Hs be the Hilbert space with reproducing kernel (1.3). The mapping
[TABLE]
extends via linearity and operator closure to define a scattering-isometric closed linear operator [A&BC&D]s:[XU]⊃dom([A&BC&D]s)→[XY].
Proof.
From es(μ,μ∗)es(λ,λ∗)∗=Ks(μ,μ∗,λ,λ∗) and (1.3), we obtain the following:
[TABLE]
i.e., for all ν,η∈U, γ,ξ∈Y, and μ,μ∗,λ,λ∗∈C+:
[TABLE]
Proceeding along the lines of the proof of Lemma I.4.1, one obtains from (4.2) that the extension by linearity and operator closure of the mapping (4.1) is a scattering-isometric, well-defined single-valued operator; please note that [z1z2]⊥x2 for all [x2u2]∈dom([A&BC&D]s) implies
[TABLE]
for all λ,λ∗∈C+, ξ∈Y, and η∈U, i.e., z1=0 and z2=0.
∎
From now on [A&BC&D]s always denotes the extension of the mapping (4.1) by linearity and operator closure.
Theorem 4.2**.**
For all φ∈S(C+;U,Y), the following claims are true:
(1)
The operator [A&BC&D]s is a simple scattering-conservative system node with input/state/output spaces (U,X,Y) and transfer function φ.
2. (2)
The adjoint [A&BC&D]s∗:[XY]⊃dom([A&BC&D]s∗)→[XU] is the extension by linearity and operator closure of
[TABLE]
3. (3)
The (unilateral) input map of [A&BC&D]s is Bs=[Γ1] and the approximately reachable subspace is \mathfrak{R}_{s}={\rm im}\bigl{(}\mathfrak{B}_{s}\bigr{)}=\left[\begin{smallmatrix}\Gamma\\
1\end{smallmatrix}\right]\mathcal{H}_{c}. The adjoint of the (unilateral) output map of [A&BC&D]s is Cs∗=[1Γ∗] and the orthogonal complement of the unobservable subspace is \mathfrak{R}_{s}^{\dagger}={\rm im}\bigl{(}\mathfrak{C}_{s}^{*}\bigr{)}=\left[\begin{smallmatrix}1\\
\Gamma^{*}\end{smallmatrix}\right]\mathcal{H}_{o}. The bilateral input and output maps of [A&BC&D]s are both equal to the identity operator on Hs.
Assertion three can be written more explicitly as (for α∈C+,u∈U,[x1x2]∈Hs)
We obtain that [A&BC&D]s is an energy-preserving system node by generalizing the proof of Thm. I.4.2: If
[TABLE]
then in particular for all λ,λ∗∈C+, γ∈Y, and ν∈U:
[TABLE]
Restricting (4.6) to the case λ∗=1 and γ=0, we obtain that u=0. Keeping γ=0 but taking λ∗=1, we get x2=0. Finally, letting γ run over Y and λ over C+, we obtain that x1=0. Combining this with Prop. 4.1 and the proof of Thm. I.4.2, we obtain that [A&BC&D]s is an energy-preserving system node. In the same way we see that the range of [[10]+[As&Bs][Cs&Ds]] is dense in [HsU]; then [A&BC&D]s is a conservative system node by Thm. I.3.12. Claim two now follows immediately from (2.4).
From (I.3.5) and (4.1) we then have (for λ∗∈C+,ν∈U):
[TABLE]
From (4.7) and (3.15) we have that Bs=[Γ1] and by definition Rs=ran(Bs). However, since Bs:Hc→Hs is isometric by Thm. 3.5.2, {\rm im}\bigl{(}\mathfrak{B}_{s}\bigr{)} is closed. By (2.5) we have eo(λ)C=C(λ−A)−1 and carrying out the calculation (4.7) for [A&BC&D]s∗ in (4.3), we obtain Cs∗=[1Γ∗]:
[TABLE]
where we also used (3.11). Combining this and Bs=[Γ1] with (3.18), we obtain for all [fg]∈[HoHc] that
[TABLE]
this shows that Bs,full acts as the identity on the dense subspace {\rm im}\bigl{(}\left[\begin{smallmatrix}1&\Gamma\\
\Gamma^{*}&1\end{smallmatrix}\right]\bigr{)} of Hs; by Thm. 3.7.4, [A&BC&D]c is simple. Using (3.19) we obtain that also Cs,full=1.
∎
We have the following variant of Thm. 2.8 regarding intertwinement with [A&BC&D]s:
Theorem 4.3**.**
Let [A&BC&D] be a conservative system with state space X and transfer function φ. Then the bilateral input map Bbil of [A&BC&D] in Thm. 3.7 intertwines [A&BC&D]s isometrically with [A&BC&D]. Additionally, [A&BC&D]s is simple if and only if Bbil is unitary.
Proof.
The isometricity of Bbil and the fact that {\rm im}\bigl{(}\mathfrak{B}_{bil}\bigr{)} is dense if and only of [A&BC&D] is simple were shown in Thm. 3.7.
Using (3.20), Lemma 2.3, (I.3.6), and (4.1), for all λ,λ∗∈C+, γ∈Y, and ν∈U:
[TABLE]
[TABLE]
Taking linear combinations of elements [λ∗(λ∗−A−1)−1Bνφ(λ∗)ν] and closing in dom([A&BC&D]s) (equipped with the graph norm), we obtain both (2.18) and (2.19).
∎
As an immediate consequence of the theorem, any two simple conservative realizations with the same transfer function are unitarily similar.
The above formulas (4.1) and (4.3) for [A&BC&D]s and its adjoint only give the action on special, kernel-like elements. Using (2.4), we can obtain explicit formulas for the action of [A&BC&D]s on generic elements of its domain:
Theorem 4.4**.**
The model [A&BC&D]s has the explicit representation given in Thm. 1.1, where the vector y can alternatively be defined as the unique y∈Y for which the function [z1z2] in (1.5) is an element of Hs.
Proof.
By (4.3), for all [xu]∈[HsU], μ,μ∗∈C+, γ∈Y, and ν∈U:
[TABLE]
First assume that [xu]∈dom([A&BC&D]s) and define [zy]:=[A&BC&D]s[xu]∈[XY], so that the left-hand side of (4.10) equals
[TABLE]
for all [γν]. Then (1.5) holds, and moreover (1.4) holds by Cor. 3.6. Thus, the action of [A&BC&D]s is correct and dom([A&BC&D]s) is contained in the set on the right-hand side of (1.6).
Now drop the assumption [xu]∈dom([A&BC&D]s) and instead assume that y∈Y is such that [z1z2] defined by (1.5) is in Hs. Then y satisfies (1.4) and by the definition of [z1z2], (4.10) equals
[TABLE]
for all [es(μ,μ∗)∗[1φ(μ∗)]][γν] which by the definition of [A&BC&D]s span a dense subspace of dom([A&BC&D]s). Thus
[TABLE]
for all [x′y′]∈dom([A&BC&D]s∗), so that \left[\begin{smallmatrix}x\\
u\end{smallmatrix}\right]\in{{\rm dom}}\left(\big{(}\left[\begin{smallmatrix}{A\&B}\cr{C\&D}\end{smallmatrix}\right]_{s}^{*}\big{)}^{*}\right).
∎
By (I.3.7), a system node [A&BC&D] with state space X can be reconstructed from its component operators A, B, C, and its transfer function D(α), for an arbitrary α∈res(A). In the following result, we describe these operators for [A&BC&D]s. In order to state the result, we define a linear operator Rα, α∈C+, on the space of analytic functions C+→[YU] by
[TABLE]
Proposition 4.5**.**
The main operator As of [A&BC&D]s is
[TABLE]
with domain consisting of those [x1x2]∈Hs for which the limit
[TABLE]
exists in Y and the function in (4.12) lies in Hs.
The observation operator is
[TABLE]
For all [x1x2]∈dom(As) it holds (with y as above) that
[TABLE]
The resolvent of As is
[TABLE]
and moreover, for all α∈C+ and [x1x2]∈Hs:
[TABLE]
μ,μ∗∈C+.
Proof.
The claims on As and Cs (including (4.14)) follow by comparing Defs. I.3.1 – 2 to (1.4) – (1.6). For α∈C+ and [w1w2]∈dom(As) arbitrary, set [x1x2]:=(α−As)[w1w2]. From (4.12)–(4.13) we see that
[TABLE]
We conclude that Cs[w1w2]=x1(α), which gives an alternative proof of (4.5), and solving for [w1w2], we get (4.15). Then (4.16) follows from (4.15) and the identity As(α−As)−1=α(α−As)−1−1.
∎
The following are easy to see (for all α∈C+):
[TABLE]
From the former of these formulas, it seems reasonable that As,−1[x1x2]=(μ,μ∗)↦[μx1(μ)−μ∗x2(μ∗)] while the latter hints that the control operator Bs of [A&BC&D]s could be Bsu=[φ(⋅)1]u. In order to properly decouple [As&Bs] into [As,−1Bs] and prove these conjectures, we shall next interpret [As&Bs] as an operator that maps into the extrapolation space Hs,−1 of Hs.
5. The extrapolation space and its reproducing kernel
The formula (4.15) for the resolvent of As suggests a way to concretely identify the (−1)-scaled rigged space Hs,−1 defined abstractly as the completion of the space Hs in the norm
[TABLE]
where β is the fixed rigging parameter. Indeed, we should have [z1z2]∈Hs,−1 if and only if Rβ[z1z2]∈Hs, see (4.11), and in this case
[TABLE]
It is straightforward to verify that Rβ[z1z2]=0⟺[z1z2]∈Zs, where
[TABLE]
in particular, RβZs={0}⊂Hs. Hence, ∥⋅∥Hs,−1 is a norm on the quotient space
[TABLE]
The norm on Hs,−1, and the corresponding inner product, depend on the choice of β∈C+, but different choices of β give equivalent norms.
Theorem 5.1**.**
The space Hs,−1 in (5.3) is a Hilbert space with the norm (5.1).
(1)
The map ι:[x1x2]↦[x1x2]+Zs embeds Hs continuously into Hs,−1 as a dense subspace.
A given element [z1z2]+Zs∈Hs,−1
is of the form ι[x1x2] for some [x1x2]∈Hs if and only if the function Rα[z1z2] is not only in Hs but in fact is in dom(As)⊂Hs, for some, or equivalently for all, α∈C+.
When [z1z2]+Zs=ι[x1x2], the representative [x1x2]∈Hs for the equivalence class [z1z2]+Zs is uniquely determined by the decay of the first component at infinity, i.e., by the condition limReη→+∞x1(η)=0.
2. (2)
When Hs is identified as a linear sub-manifold of Hs,−1 as in assertion one, for all [x1x2]∈Hs,
[TABLE]
is the unique extension of As to a closed operator on Hs,−1. Moreover,
[TABLE]
and (β−As,−1)−1 is a unitary operator from Hs,−1 onto Hs.
3. (3)
With Hs,−1 identified concretely as in (5.3) and ιHs identified with Hs, the control operator Bs:U→Hs,−1 is
[TABLE]
Proof.
The argument follows the proof of Thm. I.4.7, but we provide a more polished formulation. In order to establish that Hs,−1 is complete, take a Cauchy sequence [z1,nz2,n]+Zs in Hs,−1. Then the Cauchy sequence Rβ[z1,nz2,n] converges to some [x1x2] in Hs. Solving [x1x2]=Rβ[z1z2] for [z1z2], we obtain that [z1,nz2,n] converges in Hs,−1 to
[TABLE]
Thus, Hs,−1 is a Hilbert space and R_{\beta}\big{|}_{\mathcal{H}_{s,-1}} clearly maps Hs,−1 unitarily onto Hs.
We next prove assertion one. Combining (5.1) with (4.15), we see that ι is continuous: for all x∈Hs it holds that
[TABLE]
As R_{\beta}\big{|}_{\mathcal{H}_{s,-1}} is unitary from Hs,−1 onto Hs and dom(As)=RβιHs is dense in Hs=RβHs,−1, it follows that ιHs is dense in Hs,−1.
For all x∈Hs and α∈C+,
[TABLE]
hence [z1z2]+Zs=ι[x1x2] only if Rα[z1z2]∈dom(As), for all α∈C+. Conversely, if there exists some α∈C+ such that
so that ι(α−As)[w1w2]=[z1z2]+Zs; please observe that (α−As)[w1w2]∈Hs.
The decay condition picks out the unique representative in Hs due to Cor. 3.6, and by writing ι−1 below, we mean the inverse of ι with domain ιHs.
We next prove assertion two, and for this we temporarily denote the mapping in (5.4) by As. In the beginning of the proof, we showed that
[TABLE]
and x lies in Hs for z+Zs∈Hs,−1, and so x=Rβz if and only if x+Zs=ιRβz. In particular, β−As is injective:
[TABLE]
and we get ιRβ=(β−As)−1. Comparing to (4.15), we see that (β−As)−1 is the unique extension of the densely defined ι(β−As)−1ι−1 to an operator in L(Hs,−1); this implies that As is closed. Furthermore, (β−As)−1=(β−As,−1)−1, and inverting, we get As,−1=As. The above imply that
[TABLE]
and the identification ιHs=Hs means that we may remove ι and ι−1 from the above formulas. Since β is arbitrary in C+ the correctness of (5.5) follows.
It remains only to prove assertion three. By (5.1) and (4.19), the operator Bs in (5.6) maps into Hs,−1. Then (5.5) and (4.19) give Bs=Bs.
∎
Strictly speaking, Hs,−1 is not a RKHS, because its elements are equivalence classes of functions rather than functions. However, if we agree to represent an equivalence class in Hs,−1 by the function whose first component vanishes at the rigging point β then we can obtain a reproducing kernel for Hs,−1.
Proposition 5.2**.**
With β the parameter used in the rigging, the space
[TABLE]
with the norm (5.1) is a Hilbert space of [YU]-valued analytic functions on C+×C+. This space has the reproducing kernel
[TABLE]
The condition x1(β)=0 picks out a unique representative of every equivalence class in the extrapolation space. We call Hs,−1βthe β-normalized extrapolation space of [A&BC&D]s. By (5.8) the natural embedding of the state space Hs into Hs,−1β is
[TABLE]
Proof.
For [x1x2]=Rβ[z1z2]∈Hs with z1(β)=0, we have from (5.4):
[TABLE]
thus (5.9) works as the (unique) reproducing kernel of Hs,−1β.
∎
6. Recovering the unitary model for the disk
As in §I.6 we use the following function for mapping D one-to-one onto C+:
[TABLE]
where the parameter α∈C+ is arbitrary but fixed. The inverse of mα is
[TABLE]
Recalling from (I.6.5) that the Cayley transform with parameter α in (I.6.1) of a passive system node [A&BC&D] has the transfer function \phi_{\alpha}(z):=\varphi\big{(}m_{\alpha}(z)\big{)}, we define a positive kernel function Ks(z,w) on D×D as in (1.1) with ϕ replaced by ϕα:
[TABLE]
z,w∈D. The induced RKHS is denoted by Hs,α and evaluation at z∈D in this space is denoted by es,α(z).
for all w,z∈D, γ,y∈Y, and ν,u∈U. Thus (6.3) can be extended by linearity and operator closure into an isometry from Hs,α into Hs and by (6.4) the range of this isometry contains
[TABLE]
Take [fg]∈H perpendicular to this linear span. Setting ν=0 and observing that mα(D)=C+, we get f=0; then also g=0. Thus Ξs,α is unitary from Hs,α to Hs.
∎
The preceding lemma is in agreement with statements 2 of Propositions I.6.2 and I.6.3. Further, we have the following result:
Proposition 6.2**.**
For all α∈C+, the following claims are true:
(1)
The Cayley transform (I.6.1) with parameter α∈C+ of [A&BC&D]s is the unitary operator [As,αCs,αBs,αDs,α]:[HsU]→[HsY] given by
[TABLE]
2. (2)
The operator (6.3) implements a unitary similarity between [AsCsBsDs] in (6.5) and [AsCsBsDs] in (1.2) built for \phi(z):=\varphi\big{(}m_{\alpha}(z)\big{)}, z∈D:
[TABLE]
Proof.
Eq. (I.6.1) gives Ds,α=φ(α) and for As,α, the formula in (6.5) is immediate from Prop. 4.5. The formulas for Bs,α and Cs,α are (4.4) and (4.5) renormalized. For the intertwinement, we obtain:
[TABLE]
where again z=mα−1(μ)=mα−1(μ∗) and mα−1(α)=0.
∎
We note that (6.5) combines the Cayley transforms (I.6.6) and (I.6.18) of [A&BC&D]o and [A&BC&D]c in a way similar to how [A&BC&D]s in Thm. 4.4 combines (I.5.4–6) and (I.4.43). Also, h(α) in (6.5) plays the role of τc,αx=(Γx)(α) in (I.6.18).
Appendix A Non-invertible intertwinement
This section contains results on intertwinements, which are not part of the main story. In this section, no assumptions on passivity are made.
The standard transfer function of a system node [A&BC&D] with input space U, state space X, and output space Y only considers the input/output behavior of the system. We now extend the concept of transfer function in a way which also provides information on the state trajectory. Namely, we extend it into the mapping [x0u(λ)]↦[x(λ)y(λ)], where x0 is the initial state of a Laplace-transformable trajectory (u,x,y) of the system and the hats denote the right-sided Laplace transforms.
Definition A.1**.**
By the input/state/output (i/s/o) resolvent of a system node with i/s/o spaces (U,X,Y), we mean the following family of bounded linear operators from [XU] into [XY]:
[TABLE]
It is immediate from Prop. I.3.10 that the i/s/o resolvent of the dual system [A&BC&D]∗ at λ∈res(A∗) is Sd(λ)=S(λ)∗. Furthermore, a system node is uniquely determined by its i/s/o transfer function at any single point α∈dom(S): if S(α) is determined then the following operator is also determined:
[TABLE]
where the last operator maps [XU]ontodom([A&BC&D]); see p. 740 in Part I.
In the present paper, the four component operators of S in fact play a more important role than the system-node components A, B, C themselves, and much of the theory could be written in terms of these operators. However, here we choose a more explicit exposition which is more in line with the notation in, e.g., [ADRdS97].
Lemma A.2**.**
The following conditions are equivalent for two system nodes [A&BC&D]0 and [A&BC&D]1 with state spaces X0 and X1, respectively, and a bounded operator E:X0→X1:
(1)
The operator E intertwines [A&BC&D]0 with [A&BC&D]1.
2. (2)
The operator E∗ intertwines [A&BC&D]1∗ with [A&BC&D]0∗.
3. (3)
The following operator inclusion holds:
[TABLE]
where [A&BC&D]1[E001] is defined on its maximal domain
[TABLE]
4. (4)
For all λ∈res(A0)∩res(A1) (which contains some right-half plane of C), the operator E intertwines the i/s/o resolvent of [A&BC&D]0 with that of [A&BC&D]1:
[TABLE]
5. (5)
There exists one λ∈res(A0)∩res(A1) such that (A.3) holds.
Proof.
By the definition of operator inclusion, (2.18)–(2.19) are equivalent to (A.2), i.e., claim one holds if and only if claim three holds. The following calculation shows that claim three implies claim two:
[TABLE]
where the first inclusion holds for all (unbounded) operators, the second inclusion follows from (A.2), and the equality holds because [E001] is bounded; see [Rud73, Thm. 13.2]. This proves that statement one implies statement two, and applying this implication with [A&BC&D]1−k∗ in place of [A&BC&D]k and E∗ in place of E, we obtain also that claim two implies claim one, since the (closed) system nodes and E are all equal to their double adjoints.
In order to prove that statement one implies statement four, fix λ∈res(A0)∩res(A1) arbitrarily and assume (2.18)–(2.19). Then it is easy to see that also the following two identities hold:
[TABLE]
and
[TABLE]
The latter of these implies that
[TABLE]
since [[λ0]−[Ak&Bk][01]]−1 maps [XkU] into dom([A&BC&D]k); see the first three pages of §I.3, by which we can also write (A.1) as
[TABLE]
Multiplying (A.4) from the right by [[λ0]−[A0&B0][01]]−1, and using (A.5) and (A.6), we get (A.3) for every λ∈res(A0)∩res(A1).
Statement four implies statement five because res(A0)∩res(A1) is nonempty. In order to prove that statement five implies statement one, we assume that λ∈res(A0)∩res(A1) is such that (A.3) holds. Then we claim that
[TABLE]
Indeed, the top half follows from (A.3) and the bottom half is trivial. Multiplying (A.7) by [[λ0]−[A0&B0][01]] from the right, we obtain
[TABLE]
and further multiplying by [[λ0]−[A1&B1][01]] from the left, we get
[TABLE]
Hence, (2.18) and the top half of (2.19) hold. Finally, we multiply the bottom half of the identity [E001]S0(λ)=S1(λ)[E001], i.e,
[TABLE]
from the right by [[λ0]−[A0&B0][01]] and use (A.8), which gives us the bottom half of (2.19).
∎
We have the following consequences:
Theorem A.3**.**
Let E intertwine [A&BC&D]0 with [A&BC&D]1. Then:
(1)
If E is surjective then (2.18) and (A.2) hold with equality.
2. (2)
If E is unitary then [A&BC&D]0 is energy preserving, co-energy preserving, or conservative if and only if [A&BC&D]1 has the same property.
3. (3)
Defining E−1α:=(α−A1,−1)E(α−A0,−1)−1 for α∈res(A0)∩res(A1), we get
[TABLE]
The extrapolated intertwinement E−1α is surjective if E is surjective.
Moreover, if α=β (the rigging parameter) then E−1:=E−1β inherits the following properties from E: isometricity, co-isometricity, and unitarity.
Proof.
Assertion one follows from (A.7) and the surjectivity of E:
[TABLE]
Now assume that E is unitary, so that (A.2) holds with equality. Then, for every [xu]∈dom([A&BC&D]1) and [zy]=[A&BC&D]1[xu], we have [E∗xu]∈dom([A&BC&D]0) and [E∗zy]=[A&BC&D]0[E∗xu]. For [A&BC&D]0 energy preserving, we have
[TABLE]
thus [A&BC&D]1 inherits energy preservation from [A&BC&D]0. The converse implication is obtained by swapping the roles of [A&BC&D]0 and [A&BC&D]1 and using E∗ for E. The same argument for [E∗001][A&BC&D]1∗=[A&BC&D]0∗[E∗001] gives that [A&BC&D]0 is co-energy preserving if and only if [A&BC&D]1 is co-energy preserving. Hence, [A&BC&D]0 is conservative if and only if [A&BC&D]1 is conservative.
The second line of (A.9) is trivial by the definition of E−1α, and taking α=0 in (α−A1,−1)E=E−1α(α−A0,−1), we get A1,−1E=E−1αA0,−1. Then, for α=0, we on the other hand get E_{-1}^{\alpha}\big{|}_{\mathcal{X}_{0}}=E. Line 2 of (A.9) combined with the upper-right corner of (A.3) gives E−1αB0=B1. The rest of the assertion is immediate from the definition of E−1α and the unitarity of β−Ak,−1:Xk→Xk,−1.
∎
Bibliography21
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[AA 66] V. M. Adamjan and D. Z. Arov, Unitary couplings of semi-unitary operators , Mat. Issled. 1 (1966), no. vyp. 2, 3–64, (Russian) English translation available in Am. Math. Soc. Trans I. Ser. 2, 95 (1970), 75-129.
2[AD Rd S 97] Daniel Alpay, Aad Dijksma, James Rovnyak, and Henrik de Snoo, Schur functions, operator colligations, and reproducing kernel Hilbert spaces , Operator Theory: Advances and Applications, vol. 96, Birkhäuser-Verlag, Basel Boston Berlin, 1997.
3[AKS 11] Damir Z. Arov, Mikael Kurula, and Olof J. Staffans, Canonical state/signal shift realizations of passive continuous time behaviors , Complex Anal. Oper. Theory 5 (2011), no. 2, 331–402.
4[AN 96] Damir Z. Arov and Mark A. Nudelman, Passive linear stationary dynamical scattering systems with continuous time , Integral Equations Operator Theory 24 (1996), 1–45.
5[AS 07] Damir Z. Arov and Olof J. Staffans, State/signal linear time-invariant systems theory. Part IV: Affine representations of discrete time systems , Complex Anal. Oper. Theory 1 (2007), 457–521.
6[AS 09] by same author, Two canonical passive state/signal shift realizations of passive discrete time behaviors , J. Funct. Anal. 257 (2009), 2573–2634.
7[AS 10] by same author, Canonical conservative state/signal shift realizations of passive discrete time behaviors , J. Funct. Anal. 259 (2010), no. 12, 3265–3327.
8[B Hd S 08] Jussi Behrndt, Seppo Hassi, and Henk de Snoo, Functional models for Nevanlinna families , Opuscula Math. 28 (2008), no. 3, 233–245.