The Behavioral Approach to Systems Theory Paolo Rapisarda, Un. of - - PowerPoint PPT Presentation

the behavioral approach to systems theory
SMART_READER_LITE
LIVE PREVIEW

The Behavioral Approach to Systems Theory Paolo Rapisarda, Un. of - - PowerPoint PPT Presentation

The Behavioral Approach to Systems Theory Paolo Rapisarda, Un. of Southampton, U.K. & Jan C. Willems, K.U. Leuven, Belgium MTNS 2006 Kyoto, Japan, July 2428, 2006 Lecture 2: Representations and annihilators of LTIDS Lecturer: Jan C.


slide-1
SLIDE 1

The Behavioral Approach to Systems Theory

Paolo Rapisarda, Un. of Southampton, U.K. & Jan C. Willems, K.U. Leuven, Belgium MTNS 2006 Kyoto, Japan, July 24–28, 2006

slide-2
SLIDE 2

Lecture 2: Representations and annihilators of LTIDS Lecturer: Jan C. Willems

slide-3
SLIDE 3

Issues

  • What is a linear time-invariant differential system

(LTIDS) ?

slide-4
SLIDE 4

Issues

  • What is a linear time-invariant differential system

(LTIDS) ?

  • How are they represented?
slide-5
SLIDE 5

Issues

  • What is a linear time-invariant differential system

(LTIDS) ?

  • How are they represented?
  • The annihilators
  • Differential annihilators
  • Rational annihilators
slide-6
SLIDE 6

Issues

  • What is a linear time-invariant differential system

(LTIDS) ?

  • How are they represented?
  • The annihilators
  • Differential annihilators
  • Rational annihilators
  • Controllability, transfer functions, and image

representations

slide-7
SLIDE 7

Issues

  • What is a linear time-invariant differential system

(LTIDS) ?

  • How are they represented?
  • The annihilators
  • Differential annihilators
  • Rational annihilators
  • Controllability, transfer functions, and image

representations

  • Representations using proper stable rational functions
slide-8
SLIDE 8

LTIDS

slide-9
SLIDE 9

The class of systems

We discuss the fundamentals of the theory of dynamical systems Σ = (R, Rw, B) that are

slide-10
SLIDE 10

The class of systems

We discuss the fundamentals of the theory of dynamical systems Σ = (R, Rw, B) that are 1. linear, meaning (‘superposition’) [ [(w1, w2 ∈ B) ∧ (α, β ∈ R)] ] ⇒ [ [αw1 + βw2 ∈ B] ]

slide-11
SLIDE 11

The class of systems

We discuss the fundamentals of the theory of dynamical systems Σ = (R, Rw, B) that are 1. linear, meaning (‘superposition’) [ [(w1, w2 ∈ B) ∧ (α, β ∈ R)] ] ⇒ [ [αw1 + βw2 ∈ B] ]

  • 2. time-invariant, meaning

[ [(w ∈ B) ∧ (t′ ∈ R)] ] ⇒ [ [σt′w ∈ B)] ] σt′: backwards t′-shift: σt′w(t) := w(t + t′).

slide-12
SLIDE 12

The class of systems

We discuss the fundamentals of the theory of dynamical systems Σ = (R, Rw, B) that are 1. linear, meaning (‘superposition’) [ [(w1, w2 ∈ B) ∧ (α, β ∈ R)] ] ⇒ [ [αw1 + βw2 ∈ B] ]

  • 2. time-invariant, meaning

[ [(w ∈ B) ∧ (t′ ∈ R)] ] ⇒ [ [σt′w ∈ B)] ] σt′: backwards t′-shift: σt′w(t) := w(t + t′). 3. differential, meaning B consists of the sol’ns of a system of diff. eq’ns.

slide-13
SLIDE 13

The class of systems

w variables: w1, w2, . . . ww, up to n-times differentiated, g equations. ❀

Σw

j=1R0 1,jwj + Σw j=1R1 1,j

d dt wj + · · · + Σw

j=1Rn 1,j

dn dtn wj = Σw

j=1R0 2,jwj + Σw j=1R1 2,j

d dt wj + · · · + Σw

j=1Rn 2,j

dn dtn wj = . . . . . . . . . Σw

j=1R0 g,jwj + Σw j=1R1 g,j

d dt wj + · · · + Σw

j=1Rn g,j

dn dtn wj =

Coefficients Rk: 3 indices! i = 1, . . . , g : for the i-th differential equation, j = 1, . . . , w : for the variable wj involved, k = 1, . . . , n : for the order dk

dtk of differentiation.

slide-14
SLIDE 14

The class of systems

In vector/matrix notation: w =      w1 w2, . . . ww      , Rk =      Rk

1,1

Rk

1,2

· · · Rk

1,w

Rk

2,1

Rk

2,2

· · · Rk

2,w

. . . . . . · · · . . . Rk

g,1

Rk

g,2

· · · Rk

g,w

     .

slide-15
SLIDE 15

The class of systems

In vector/matrix notation: w =      w1 w2, . . . ww      , Rk =      Rk

1,1

Rk

1,2

· · · Rk

1,w

Rk

2,1

Rk

2,2

· · · Rk

2,w

. . . . . . · · · . . . Rk

g,1

Rk

g,2

· · · Rk

g,w

     . R0w + R1 d dt w + · · · + Rn dn dtn w = 0, with R0, R1, · · · , Rn ∈ Rg×w.

slide-16
SLIDE 16

The class of systems

In vector/matrix notation: w =      w1 w2, . . . ww      , Rk =      Rk

1,1

Rk

1,2

· · · Rk

1,w

Rk

2,1

Rk

2,2

· · · Rk

2,w

. . . . . . · · · . . . Rk

g,1

Rk

g,2

· · · Rk

g,w

     . R0w + R1 d dt w + · · · + Rn dn dtn w = 0, with R0, R1, · · · , Rn ∈ Rg×w. With polynomial matrix R(ξ) = R0 + R1ξ + · · · + Rnξn ∈ R [ξ]g×w we obtain the mercifully short notation

R( d dt )w = 0.

slide-17
SLIDE 17

Definition of the behavior

What shall we mean by the behavior of R( d dt )w = 0 ? Solutions in C∞ (R, Rw)? As many times differentiable as there appear derivatives ap- pear in DE ? Distributional solutions in Lloc(R, Rw)? In Lloc

2 (R, Rw)?

Distributions?

slide-18
SLIDE 18

Definition of the behavior

What shall we mean by the behavior of R( d dt )w = 0 ? Solutions in C∞ (R, Rw)? As many times differentiable as there appear derivatives ap- pear in DE ? Distributional solutions in Lloc(R, Rw)? In Lloc

2 (R, Rw)?

Distributions? The easy way out B := {w ∈ C∞ (R, Rw) | R( d dt )w(t) = 0 ∀ t ∈ R} Notation: B = ker(R( d

dt ))

slide-19
SLIDE 19

Notation

R [ξ] : polynomials with real coeff., indeterminate ξ R [ξ]n×m: polynomial matrices R [ξ]•ו: appropriate number of rows, columns Lw, L•: linear differential systems B ∈ Lw := (R, Rw, B) ∈ Lw; B = ker(R( d

dt ))

R(ξ) : rational f’ns with real coeff., indeterminate ξ R(ξ)n×m: matrices of rat. f’ns R(ξ)•ו: appropriate number of rows, columns

slide-20
SLIDE 20

Rational symbols

We also want to give a meaning to

F( d dt )w = 0

with F ∈ R(ξ)•×w, i.e. a matrix of rational functions. What do we mean by a solution?

slide-21
SLIDE 21

Rational symbols

We also want to give a meaning to

F( d dt )w = 0

with F ∈ R(ξ)•×w, i.e. a matrix of rational functions. What do we mean by a solution? We do this in terms of a left co-prime polynomial factoriza- tion. F(ξ) = P(ξ)−1Q(ξ) with P, Q ∈ R [ξ]•ו , det(P) = 0,

  • P

Q

  • left prime.
slide-22
SLIDE 22

Rational symbols

We also want to give a meaning to

F( d dt )w = 0

with F ∈ R(ξ)•×w, i.e. a matrix of rational functions. What do we mean by a solution? We do this in terms of a left co-prime polynomial factoriza- tion. F(ξ) = P(ξ)−1Q(ξ) with P, Q ∈ R [ξ]•ו , det(P) = 0,

  • P

Q

  • left prime.

Define the behavior of this ‘diff. eq’n’ to be that of

Q( d dt )w = 0

Whence ∈ L•.

slide-23
SLIDE 23

Elimination

slide-24
SLIDE 24

Problem

Assume (w1, w2) governed by R1( d dt )w1 = R2( d dt )w2 R1, R2 ∈ R [ξ]•ו. Behavior B. Obviously B ∈ L• Define the ‘projection’ B1 := {w1 | ∃ w2 such that (w1, w2) ∈ B}

slide-25
SLIDE 25

Problem

Assume (w1, w2) governed by R1( d dt )w1 = R2( d dt )w2 R1, R2 ∈ R [ξ]•ו. Behavior B. Obviously B ∈ L• Define the ‘projection’ B1 := {w1 | ∃ w2 such that (w1, w2) ∈ B} Does B1 belong to L• ?

slide-26
SLIDE 26

Problem

Assume (w1, w2) governed by R1( d dt )w1 = R2( d dt )w2 R1, R2 ∈ R [ξ]•ו. Behavior B. Obviously B ∈ L• Define the ‘projection’ B1 := {w1 | ∃ w2 such that (w1, w2) ∈ B} Does B1 belong to L• ? Theorem: It does indeed, also with R1, R2 ∈ R (ξ)•ו. Algorithms?

slide-27
SLIDE 27

Examples

The input/output behavior of d dt x = Ax + Bu, y = Cx + Du. Every B ∈ L• admits such a representation w ∼ = u y

  • .

Also representation P( d dt )y = Q( d dt )u, w ∼ = u y

  • .
slide-28
SLIDE 28

Examples

The input/output behavior of d dt x = Ax + Bu, y = Cx + Du. The manifest behavior of R( d dt )w = M( d dt )ℓ, R, M ∈ R (ξ)•ו

slide-29
SLIDE 29

Examples

The input/output behavior of d dt x = Ax + Bu, y = Cx + Du. The manifest behavior of R( d dt )w = M( d dt )ℓ, R, M ∈ R (ξ)•ו Any combination of variables in a signal flow graph with rational t’f f’ns in the edges, is an LTIDS.

slide-30
SLIDE 30

Examples

The input/output behavior of d dt x = Ax + Bu, y = Cx + Du. The manifest behavior of R( d dt )w = M( d dt )ℓ, R, M ∈ R (ξ)•ו Any combination of variables in a signal flow graph with rational t’f f’ns in the edges, is an LTIDS. The port behavior of a circuit with (a finite number) linear re- sistors, capacitors, inductors, transformers, and gyrators.

slide-31
SLIDE 31

Examples

The input/output behavior of d dt x = Ax + Bu, y = Cx + Du. The manifest behavior of R( d dt )w = M( d dt )ℓ, R, M ∈ R (ξ)•ו Any combination of variables in a signal flow graph with rational t’f f’ns in the edges, is an LTIDS. The port behavior of a circuit with (a finite number) linear re- sistors, capacitors, inductors, transformers, and gyrators. Expect this to be a particular situation for LTIDS – but also holds for linear constant coefficient PDE’s.

slide-32
SLIDE 32

The annihilators

slide-33
SLIDE 33

Polynomial annihilators

Let B ∈ Lw, and n ∈ R [ξ]1×w. Call n a polynomial annihilator of B :⇔ n( d dt )w = 0 ∀ w ∈ B, i.e. iff n( d dt )B = 0. Denote the set of annihilators by NR[ξ]

B .

The term consequence is also used.

slide-34
SLIDE 34

Polynomial annihilators

Let B ∈ Lw, and n ∈ R [ξ]1×w. Call n a polynomial annihilator of B :⇔ n( d dt )w = 0 ∀ w ∈ B, i.e. iff n( d dt )B = 0. Denote the set of annihilators by NR[ξ]

B .

Easy: NR[ξ]

B

is an R [ξ]-module. This means that [ [n1, n2 ∈ NR[ξ]

B

and p ∈ R [ξ]] ] ⇒ [ [n1 + n2 ∈ NR[ξ]

B

and pn1 ∈ NR[ξ]

B ]

]

slide-35
SLIDE 35

Polynomial annihilators

Call n a polynomial annihilator of B :⇔ n( d dt )w = 0 ∀ w ∈ B, i.e. iff n( d dt )B = 0. Denote the set of annihilators by NR[ξ]

B .

Easy: NR[ξ]

B

is an R [ξ]-module. Theorem:

  • 1. Let B = ker(R( d

dt )). Then NR[ξ] B

is the R [ξ]-module generated by the rows of R.

  • 2. There is a 1:1 relation between Lw and the submodules
  • f R [ξ]1×w, the correspondence being

B → NR[ξ]

B

submodule → {w | n( d

dt )w = 0 ∀n ∈ submodule}

slide-36
SLIDE 36

Properties of Polynomial Annihilators

Every submodule of R [ξ]1×w is finitely generated. Number of generators ≤ w.

slide-37
SLIDE 37

Properties of Polynomial Annihilators

Every submodule of R [ξ]1×w is finitely generated. Number of generators ≤ w. R1( d

dt )w = 0

and R2( d

dt )w = 0 define the same system iff

∃ F1, F2 such that R2 = F1R1, R1 = F2R2

slide-38
SLIDE 38

Properties of Polynomial Annihilators

Every submodule of R [ξ]1×w is finitely generated. Number of generators ≤ w. R1( d

dt )w = 0

and R2( d

dt )w = 0 define the same system iff

∃ F1, F2 such that R2 = F1R1, R1 = F2R2 R( d

dt )w = 0 has minimal number of rows among all kernel

representations of same behavior iff R has full row rank. R1( d

dt )w = 0

and R2( d

dt )w = 0 are minimal kernel repr. of

the same system iff ∃ unimodular F such that R2 = FR1. ❀ canonical forms, etc. Basically, therefore, polynomial kernel representations are unique up to unimodular pre-multiplication

slide-39
SLIDE 39

Examples

p( d dt )w = 0 p ∈ R [ξ] Polynomial annihilators: q ∈ R [ξ] with p as a factor: R [ξ] p. Canonical form: p monic. There are also non-minimal representations, e.g. p1( d

dt )w = 0

p2( d

dt )w = 0

with GCD(p1, p2)=p. Exercise: What are the consequences of d

dt w = Aw?

slide-40
SLIDE 40

Proof of elimination thm

‘Fundamental principle’ . When is the equation F(x) = y y given , x unknown solvable? In particular, when is F( d dt )x = y solvable?

slide-41
SLIDE 41

Proof of elimination thm

‘Fundamental principle’ . When is the equation F(x) = y y given , x unknown solvable? In particular, when is F( d dt )x = y solvable? Obvious necessary condition: N ◦ F = 0 ⇒ N(y) = 0 Is this also sufficient, for a ‘small’ set of N’s?

slide-42
SLIDE 42

Proof of elimination thm

‘Fundamental principle’ . When is the equation F(x) = y y given , x unknown solvable? In particular, when is F( d dt )x = y solvable? Obvious necessary condition: N ◦ F = 0 ⇒ N(y) = 0 Is this also sufficient, for a ‘small’ set of N’s? For example, for F a matrix. Then easy to see n.a.s.c. for solvability: n ∈ R•, nF = 0 ⇒ ny = 0

slide-43
SLIDE 43

Proof of elimination thm

In particular, when is F( d dt )x = y solvable? N.a.s.c. for linear diff. eq’ns: n( d dt )F( d dt ) = 0 ⇒ n( d dt )y = 0 These n’s form a R [ξ]-module: n(ξ) such that n(ξ)F(ξ) = 0. Computable! For what w’s is R( d

dt )w = M( d dt )ℓ solvable for ℓ?

Iff nM = 0 ⇒ n( d

dt )R( d dt )w = 0.

❀ condition R′( d

dt )w = 0: elim’ion th’m + algorithm.

slide-44
SLIDE 44

Proof of elimination thm

The fundamental principle and the elimination theorem also hold for linear constant coefficient PDE’s!

Palamodov Malgrange

slide-45
SLIDE 45

Rational Annihilators

Let B ∈ Lw, and n ∈ R (ξ)1×w. Call n a rational annihilator of B :⇔ n( d dt )w = 0 ∀w ∈ B, i.e. iff n( d dt )B = 0. Note what this means: n = p−1 q1 q2 · · · qw

  • ;

p, q1, q2, . . . , qw co-prime :⇔ q1( d

dt )w1 + q2( d dt )w2 + · · · + qw( d dt )ww = 0 ∀w ∈ B.

Denote the set of rational annihilators by NR(ξ)

B .

slide-46
SLIDE 46

Rational Annihilators

Let B ∈ Lw, and n ∈ R (ξ)1×w. Call n a rational annihilator of B :⇔ n( d dt )w = 0 ∀w ∈ B, i.e. iff n( d dt )B = 0. Denote the set of rational annihilators by NR(ξ)

B .

It is easy to see that NR(ξ)

B

is R [ξ]-module. (Prove!) But, now, a sub-module of R (ξ)1×w viewed as a R [ξ]- module.

slide-47
SLIDE 47

Rational Annihilators

Call n a rational annihilator of B :⇔ n( d dt )w = 0 ∀w ∈ B, i.e. iff n( d dt )B = 0. Denote the set of rational annihilators by NR(ξ)

B .

Theorem:

  • 1. Let B = ker(R( d

dt )). Then NR(ξ) B

is the R [ξ]-module generated by the rows of R.

  • 2. There is a 1:1 relation between Lw and the R [ξ]

submodules of R (ξ)1×w, the correspondence being

B → NR(ξ)

B

submodule → {w | n( d

dt )w = 0 ∀n ∈ submodule}

Not a nice thm: refers to submodules of a vector space!

slide-48
SLIDE 48

Examples

p( d dt )w = 0 p ∈ R [ξ] Rational annihilators: n1 n2 ∈ R (ξ) with n1, n2 co-prime, and with p a factor of n1.

slide-49
SLIDE 49

Examples

p( d dt )w1 = q( d dt )w2 p, q ∈ R [ξ] Rational annihilators: n1 n2

  • p

−q

  • ,

with n1, n2 ∈ R [ξ], co-prime, and with n2, p, q co-prime. In the special case that p, q are co-prime, this is actually the R (ξ)-vector space generated by

  • p

−q ∼ =

  • 1

−q p

  • !

Why do we get a subspace instead of just a module?

slide-50
SLIDE 50

Controllability & Stabilizability

slide-51
SLIDE 51

Controllability

2 1

w w W time

2

T

1

w w ! w

T

W time W

slide-52
SLIDE 52

Stabilizability

Stabilizability :⇔ legal trajectories can be steered to a desired point.

w’ w W time

slide-53
SLIDE 53

Tests

Theorem: B = ker(R( d

dt )), R ∈ R [ξ]•ו is controllable ⇔

R(λ) has the same rank for all λ ∈ C Same result for rational symbols, but care should be taken in defining rank drop in situations where the symbol has zeros and poles in common points of the complex plane.

slide-54
SLIDE 54

Tests

Theorem: B = ker(R( d

dt )), R ∈ R [ξ]•ו is controllable ⇔

R(λ) has the same rank for all λ ∈ C Same result for rational symbols, but care should be taken in defining rank drop in situations where the symbol has zeros and poles in common points of the complex plane. Example 1:

d dt x = Ax + Bu, dim(x) = n is controllable iff

rank(

  • λIn − A

B

  • ) = n for all λ ∈ C.

Example 2: y = G( d

dt )u, w =

u y

  • is always controllable.
slide-55
SLIDE 55

Tests

Theorem: B = ker(R( d

dt )) ∈ Lw, R ∈ R [ξ]•ו is stabilizable ⇔

R(λ) has the same rank for all λ with Re(λ) ≥ 0 Same result for rational symbols, but care should be taken in defining rank drop in situations where the symbol has zeros and poles in common points of the complex plane.

slide-56
SLIDE 56

Tests

Theorem: B = ker(R( d

dt )) ∈ Lw, R ∈ R [ξ]•ו is stabilizable ⇔

R(λ) has the same rank for all λ with Re(λ) ≥ 0 Same result for rational symbols, but care should be taken in defining rank drop in situations where the symbol has zeros and poles in common points of the complex plane. Example 1:

d dt x = Ax + Bu, dim(x) = n is stabilizable iff

rank(

  • λIn − A

B

  • ) = n for all λ with Re(λ) ≥ 0.

Example 2: y = G( d

dt )u, w =

u y

  • is always controllable,

and hence stabilizable.

slide-57
SLIDE 57

Subspaces of annihilators

Characterization of controllability in terms of the structure

  • f rational annihilators:

Theorem:

  • 1. B ∈ Lw is controllable iff its rational annihilators NR(ξ)

B

form an R (ξ)-subspace of R (ξ)1×w.

  • 2. There is a one-to-one relation between the controllable

systems in Lw and the R (ξ)-subspaces of R (ξ)1×w.

slide-58
SLIDE 58

Subspaces of annihilators

Characterization of controllability in terms of the structure

  • f rational annihilators:

Theorem:

  • 1. B ∈ Lw is controllable iff its rational annihilators NR(ξ)

B

form an R (ξ)-subspace of R (ξ)1×w.

  • 2. There is a one-to-one relation between the controllable

systems in Lw and the R (ξ)-subspaces of R (ξ)1×w. The system P( d dt )y = Q( d dt )u is equal to y = G( d dt )u with G = P−1Q iff controllable (i.e., P, Q left co-prime:

  • P

Q

  • left prime.

Transfer functions deal with controllable systems (only) .

slide-59
SLIDE 59

Kernels and images

Each element of L• is by definition the kernel of a linear constant coefficient differential operator, i.e. [ [ B ∈ L• ] ] :⇔ [ [ ∃R ∈ R [ξ]•ו such that B = ker(R( d dt )) ] ] Consider the manifest behavior of w = M( d dt )ℓ, i.e. B = im(M( d dt )) By the elimination theorem im(M( d

dt )) ∈ L•.

Easy: ∃B ∈ L• that do not admit image representation. What system theoretic property characterizes image repr.?

slide-60
SLIDE 60

Kernels and images

Each element of L• is by definition the kernel of a linear constant coefficient differential operator, i.e. [ [ B ∈ L• ] ] :⇔ [ [ ∃R ∈ R [ξ]•ו such that B = ker(R( d dt )) ] ] Consider the manifest behavior of w = M( d dt )ℓ, i.e. B = im(M( d dt )) By the elimination theorem im(M( d

dt )) ∈ L•.

Easy: ∃B ∈ L• that do not admit image representation. What system theoretic property characterizes image repr.? Controllability !!

slide-61
SLIDE 61

Image Representation

Theorem: The following are equivalent for B ∈ Lw:

  • 1. it is controllable
  • 2. ∃ M ∈ R [ξ]•ו such that B is the manifest behavior of

w = M( d dt )ℓ

  • 3. ∃ M ∈ R (ξ)•ו such that B is the manifest behavior of

w = M( d dt )ℓ Controllable iff ∃ image representation. B = im(M( d

dt )).

But be careful to interpret this in the rational case: M( d

dt )

is then a one-to-many ‘map’.

slide-62
SLIDE 62

Image Representation

Theorem: The following are equivalent for B ∈ Lw:

  • 1. it is controllable
  • 2. ∃ M ∈ R [ξ]•ו such that B is the manifest behavior of

w = M( d dt )ℓ

  • 3. ∃ M ∈ R (ξ)•ו such that B is the manifest behavior of

w = M( d dt )ℓ Controllable iff ∃ image representation. B = im(M( d

dt )).

But be careful to interpret this in the rational case: M( d

dt )

is then a one-to-many ‘map’. We may assume WLOG these image repr. observable .

slide-63
SLIDE 63

Controllable part

The controllable part of B ∈ Lw is defined as the largest controllable system B′ ∈ Lw with B′ ⊆ B.

slide-64
SLIDE 64

Controllable part

The controllable part of B ∈ Lw is defined as the largest controllable system B′ ∈ Lw with B′ ⊆ B. Two i/o systems have the same t’f f’n iff they have same controllable part. Transfer functions deal with controllable parts only.

slide-65
SLIDE 65

Controllable part

The controllable part of B ∈ Lw is defined as the largest controllable system B′ ∈ Lw with B′ ⊆ B. Two i/o systems have the same t’f f’n iff they have same controllable part. Transfer functions deal with controllable parts only. The R (ξ)-span of the rows of R in R( d

dt )w = 0 define the

rational annihilators of the controllable part.

slide-66
SLIDE 66

Prime representations

slide-67
SLIDE 67

Primes in rings

A ring is closed under addition and multiplication. Matrices, uni-modularity, etc. Let R be a ring. A matrix M ∈ R•ו is left prime if M = FM′ ⇒ F is unimodular. The matrices M1, M2, . . . , Mn, ∈ Rmו are said to be left coprime if

  • M1

M2 · · · Mn

  • is left prime.

There is an enormous zoology of rings with all sorts of properties...

slide-68
SLIDE 68

Other rings

Consider

  • 1. R [ξ]: polynomials
  • 2. R (ξ): rational functions
  • 3. R (ξ)proper: proper rational
  • 4. R (ξ)proper/stable: proper (Hurwitz) stable rational

These are all rings, with R (ξ) as field of fractions. R (ξ)proper/stable is an Euclidean domain ⇒ Bézout.

  • Matrices. Primeness, unimodularity, factorization, etc.
slide-69
SLIDE 69

Other rings

Consider

  • 1. R [ξ]: polynomials
  • 2. R (ξ): rational functions
  • 3. R (ξ)proper: proper rational
  • 4. R (ξ)proper/stable: proper (Hurwitz) stable rational

These are all rings, with R (ξ) as field of fractions. R (ξ)proper/stable is an Euclidean domain ⇒ Bézout.

  • Matrices. Primeness, unimodularity, factorization, etc.

Every B ∈ Lw admits by definition a ‘kernel repr.’ over R [ξ] i.e., ∃ R ∈ R [ξ]•ו such that B = ker(R( d

dt )).

How about the other rings? Should we care?

slide-70
SLIDE 70

Other rings

Consider

  • 1. R [ξ]: polynomials
  • 2. R (ξ): rational functions
  • 3. R (ξ)proper: proper rational
  • 4. R (ξ)proper/stable: proper (Hurwitz) stable rational

These are all rings, with R (ξ) as field of fractions. R (ξ)proper/stable is an Euclidean domain ⇒ Bézout.

  • Matrices. Primeness, unimodularity, factorization, etc.

Every B ∈ Lw admits by definition a ‘kernel repr.’ over R [ξ] i.e., ∃ R ∈ R [ξ]•ו such that B = ker(R( d

dt )).

How about the other rings? Should we care? Yes! Youla parametrization, dist. between systems, robustness, etc.

slide-71
SLIDE 71

Ring representations

Relation between system properties and prime representability over various rings. Theorem: Refers to ‘kernel repr.’ with rational symbols.

  • 1. B ∈ L• iff it admits a kernel representation with R in

and left prime over R (ξ)proper.

  • 2. B ∈ L• is stabilizable iff it admits a kernel repr. with R

in and left prime over R (ξ)proper/stable.

  • 3. B ∈ L• is controllable iff it admits a kernel

representation with R ∈ R [ξ]•ו left prime over R [ξ].

  • M. Vidyasagar
slide-72
SLIDE 72

Unitary Representation

To close this lecture, a result on unitary representations. Consider B ∈ Lw, controllable. Define B2 = B∩L2(R, Rw). B2 is a closed linear subspace of L2(R, Rw). Are there kernel or image representations that are adapted to this Hilbert space structure?

slide-73
SLIDE 73

Unitary Representation

G ∈ R (ξ)•ו, and consider the system f2 = G( d dt )f1, with f1, f2 ∈ L2(R, R•). Is this a map f1 → f2? If G is proper, no poles on the imaginary axis, then f2 = G( d

dt )f1 defines a bounded linear operator from

f1 ∈ L2(R, R•) → f2 ∈ L2(R, R•). Norm preserving (:⇔ ||f1||2 = ||f2||2) iff G⊤(−iω)G(iω) = I ∀ω ∈ R.

slide-74
SLIDE 74

Unitary Representation

B (controllable) admits a rational kernel representation R( d dt )w = 0 with R proper stable, left prime, and norm preserving. B (controllable) also admits a rational image representation w = M( d dt )ℓ with M proper stable, right prime, and norm preserving.

slide-75
SLIDE 75

Unitary Representation

B (controllable) admits a rational kernel representation R( d dt )w = 0 with R proper stable, left prime, and norm preserving. Idea of proof: start with minimal pol. repr. R( d

dt )w = 0.

Consider the polynomial matrix factorization equation R⊤(−ξ)R(ξ) = F ⊤(−ξ)F(ξ). Take Hurwitz sol’n H. Define the rational kernel repr. G( d dt )w = 0 with G = RH−1

slide-76
SLIDE 76

Summary

slide-77
SLIDE 77

Summary

  • LTIDS = linear, time-invariant, differential. Behavior

defined as sol’ns of constant coeff. diff. eqn’s. Or with rational symbols.

slide-78
SLIDE 78

Summary

  • LTIDS = linear, time-invariant, differential. Behavior

defined as sol’ns of constant coeff. diff. eqn’s. Or with rational symbols.

  • Closed under +, ∩, projection (elimination) , rational
  • perators, etc.
slide-79
SLIDE 79

Summary

  • LTIDS = linear, time-invariant, differential. Behavior

defined as sol’ns of constant coeff. diff. eqn’s. Or with rational symbols.

  • Closed under +, ∩, projection (elimination) , rational
  • perators, etc.
  • Annihilators: polynomial and rational.
slide-80
SLIDE 80

Summary

  • LTIDS = linear, time-invariant, differential. Behavior

defined as sol’ns of constant coeff. diff. eqn’s. Or with rational symbols.

  • Closed under +, ∩, projection (elimination) , rational
  • perators, etc.
  • Annihilators: polynomial and rational.
  • Controllability ⇔ image representation.
slide-81
SLIDE 81

Summary

  • LTIDS = linear, time-invariant, differential. Behavior

defined as sol’ns of constant coeff. diff. eqn’s. Or with rational symbols.

  • Closed under +, ∩, projection (elimination) , rational
  • perators, etc.
  • Annihilators: polynomial and rational.
  • Controllability ⇔ image representation.
  • Math. characterization of Lw:
  • 1:1 relation between Lw and R [ξ]-submodules
  • 1:1 relation between Lw

controllable and R (ξ)-subspaces

slide-82
SLIDE 82

Summary

  • LTIDS = linear, time-invariant, differential. Behavior

defined as sol’ns of constant coeff. diff. eqn’s. Or with rational symbols.

  • Closed under +, ∩, projection (elimination) , rational
  • perators, etc.
  • Annihilators: polynomial and rational.
  • Controllability ⇔ image representation.
  • Math. characterization of Lw:
  • 1:1 relation between Lw and R [ξ]-submodules
  • 1:1 relation between Lw

controllable and R (ξ)-subspaces

various more refined rational representations

slide-83
SLIDE 83

Discrete time systems

What changes for discrete time systems?? Ring

  • for T = N also R [ξ]
  • for T = Z instead R(ξ, ξ−1). This implies some

differences. All major thms remain valid, mutatis mutandis.

slide-84
SLIDE 84

Discrete time systems

What changes for discrete time systems?? There is a nice, ‘higher level’, definition of a linear time- invariant discrete time system.

slide-85
SLIDE 85

Discrete time systems

What changes for discrete time systems?? There is a nice, ‘higher level’, definition of a linear time- invariant discrete time system. Take T = N. The following are equivalent.

  • B linear, shift-inv., closed (pointwise conv.)
  • B linear, time-inv., complete (‘prefix determined’)

:= [ [ w ∈ B ] ] ⇔ [ [ w[t0,t1] ∈ B[t0,t1] ∀t0, t1 ∈ N ] ]

  • ∃ R ∈ R [ξ]•ו (or ∈ R (ξ)•ו) such that:

B = {w : N → R• | R(σ)w = 0 }

  • and the many more traditional representations