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 3: State and state construction Lecturer: Paolo Rapisarda


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 3: State and state construction Lecturer: Paolo Rapisarda

slide-3
SLIDE 3

Outline

The axiom of state Discrete-time systems First-order representations State maps The shift-and-cut map Algebraic characterization Continuous-time systems Computation of state-space representations

slide-4
SLIDE 4

Questions

  • Are state representations “natural"?
slide-5
SLIDE 5

Questions

  • Are state representations “natural"?
  • Mechanics ❀ 2nd order;

SYSID, transfer functions ❀ high-order;

slide-6
SLIDE 6

Questions

  • Are state representations “natural"?
  • Mechanics ❀ 2nd order;

SYSID, transfer functions ❀ high-order;

  • First principles and “tearing and zooming" modelling

❀ systems of high-order differential equations

slide-7
SLIDE 7

Questions

  • Are state representations “natural"?
  • Mechanics ❀ 2nd order;

SYSID, transfer functions ❀ high-order;

  • First principles and “tearing and zooming" modelling

❀ systems of high-order differential equations

  • Algebraic constraints among variables
slide-8
SLIDE 8

Questions

  • Are state representations “natural"?
  • Mechanics ❀ 2nd order;

SYSID, transfer functions ❀ high-order;

  • First principles and “tearing and zooming" modelling

❀ systems of high-order differential equations

  • Algebraic constraints among variables
  • What makes a latent variable a “state"?
slide-9
SLIDE 9

Questions

  • Are state representations “natural"?
  • Mechanics ❀ 2nd order;

SYSID, transfer functions ❀ high-order;

  • First principles and “tearing and zooming" modelling

❀ systems of high-order differential equations

  • Algebraic constraints among variables
  • What makes a latent variable a “state"?
  • What does that imply for the equations?
slide-10
SLIDE 10

Questions

  • Are state representations “natural"?
  • Mechanics ❀ 2nd order;

SYSID, transfer functions ❀ high-order;

  • First principles and “tearing and zooming" modelling

❀ systems of high-order differential equations

  • Algebraic constraints among variables
  • What makes a latent variable a “state"?
  • What does that imply for the equations?
  • How to construct a state from the equations?
slide-11
SLIDE 11

Questions

  • Are state representations “natural"?
  • Mechanics ❀ 2nd order;

SYSID, transfer functions ❀ high-order;

  • First principles and “tearing and zooming" modelling

❀ systems of high-order differential equations

  • Algebraic constraints among variables
  • What makes a latent variable a “state"?
  • What does that imply for the equations?
  • How to construct a state from the equations?
  • How to construct a state representation from the

equations?

slide-12
SLIDE 12

The basic idea

It’s the Mariners’ final game in the World Series. You’re late... The current score is what matters...

slide-13
SLIDE 13

The basic idea

  • The state contains all the relevant information

about the future behavior of the system

  • The state is the memory of the system
  • Independence of past and future given the state
slide-14
SLIDE 14

The axiom of state

Σ = (T, W, X, Bfull) is a state system if (w1, x1), (w2, x2) ∈ Bfull and x1(T) = x2(T) ⇓ (w1, x1) ∧

T (w2, x2) ∈ Bfull

T is concatenation at T:

(f1 ∧

T f2)(t) :=

  • f1(t) for t < T

f2(t) for t ≥ T

slide-15
SLIDE 15

Graphically...

(w1, x1), (w2, x2) ∈ Bfull and x1(T) = x2(T) ⇓ (w1, x1) ∧

T (w2, x2) ∈ Bfull

slide-16
SLIDE 16

Graphically...

(w1, x1), (w2, x2) ∈ Bfull and x1(T) = x2(T) ⇓ (w1, x1) ∧

T (w2, x2) ∈ Bfull

(w ,x )

1 1

X W

2

(w ,x )

2

time

slide-17
SLIDE 17

Graphically...

(w1, x1), (w2, x2) ∈ Bfull and x1(T) = x2(T) ⇓ (w1, x1) ∧

T (w2, x2) ∈ Bfull

(w ,x )

1 1

X W

2

(w ,x )

2

time

(w ,x ) (w ,x )

W

1 2

^

1 2

X

time

slide-18
SLIDE 18

Example 1: discrete-time system

Σ = (Z, Rw, Rl, Bfull), with Bfull := {(w, ℓ) | F ◦(σℓ, ℓ, w) = 0} where σ : (Rl)Z → (Rl)Z (σ(ℓ))(k) := ℓ(k + 1)

slide-19
SLIDE 19

Example 1: discrete-time system

Σ = (Z, Rw, Rl, Bfull), with Bfull := {(w, ℓ) | F ◦(σℓ, ℓ, w) = 0} Special case: input-state-output equations σx = f(x, u) y = h(x, u) w = (u, y)

slide-20
SLIDE 20

Example 2: continuous-time system

Σ = (R, Rw, Rl, Bfull), with Bfull := {(w, ℓ) | F ◦ ( d

dt ℓ, ℓ, w) = 0}

Special case: input-state-output equations d dt x = f(x, u) y = h(x, u) w = (u, y)

slide-21
SLIDE 21

Outline

The axiom of state Discrete-time systems First-order representations State maps The shift-and-cut map Algebraic characterization Continuous-time systems Computation of state-space representations

slide-22
SLIDE 22

First-order discrete-time representations

Theorem: A ‘complete’ latent variable system Σ = (Z, Rw, Rx, Bfull) is a state system iff Bfull can be described by F ◦ (σx, x, w) = 0

slide-23
SLIDE 23

First-order discrete-time representations

Theorem: A ‘complete’ latent variable system Σ = (Z, Rw, Rx, Bfull) is a state system iff Bfull can be described by F ◦ (σx, x, w) = 0 0-th order in w, 1st order in x

slide-24
SLIDE 24

First-order discrete-time representations

Theorem: A ‘complete’ latent variable system Σ = (Z, Rw, Rx, Bfull) is a state system iff Bfull can be described by F ◦ (σx, x, w) = 0 Linear case: Eσx+Fx+Gw = 0

slide-25
SLIDE 25

First-order discrete-time representations

Theorem: A ‘complete’ latent variable system Σ = (Z, Rw, Rx, Bfull) is a state system iff Bfull can be described by F ◦ (σx, x, w) = 0 Linear case: Eσx+Fx+Gw = 0 1st order in x is equivalent to state property!

slide-26
SLIDE 26

Proof (linear case)

V :=      a b c   | ∃(x, w) ∈ Bfull s. t.   x(1) x(0) w(0)   =   a b c      V linear ⇒ ∃ E, F, G s.t. V = ker( E F G )

slide-27
SLIDE 27

Proof (linear case)

V :=      a b c   | ∃(x, w) ∈ Bfull s. t.   x(1) x(0) w(0)   =   a b c      V linear ⇒ ∃ E, F, G s.t. V = ker( E F G ) ⇓ [(x, w) ∈ Bfull = ⇒ Eσx + Fx + Gw = 0]

slide-28
SLIDE 28

Proof (linear case)

V :=      a b c   | ∃(x, w) ∈ Bfull s. t.   x(1) x(0) w(0)   =   a b c      V linear ⇒ ∃ E, F, G s.t. V = ker( E F G ) Converse by induction, using axiom of state: Eσx + Fx + Gw = 0 on [0, k] = ⇒ (w, x)|[0,k] ∈ Bfull|[0,k] Then apply completeness of B

slide-29
SLIDE 29

State construction: basic idea

Problem: Given kernel or hybrid description, find a state representation Eσx + Fx + Gw = 0

slide-30
SLIDE 30

State construction: basic idea

Problem: Given kernel or hybrid description, find a state representation Eσx + Fx + Gw = 0 First compute polynomial operator in the shift acting

  • n system variables, inducing a state variable:

X(σ)w = x X(σ)

  • w

  • = x
slide-31
SLIDE 31

State construction: basic idea

Problem: Given kernel or hybrid description, find a state representation Eσx + Fx + Gw = 0 First compute polynomial operator in the shift acting

  • n system variables, inducing a state variable:

X(σ)w = x X(σ)

  • w

  • = x

Then use original eqs. and X to obtain 1st order representation.

slide-32
SLIDE 32

State maps for kernel representations

X ∈ R•×w[ξ] induces a state map X(σ) for ker(R(σ)) if the behavior Bfull with latent variable x, described by R(σ)w = X(σ)w = x satisfies the axiom of state.

slide-33
SLIDE 33

Example

B = {w | r(σ)w = 0} where r ∈ R[ξ], deg(r) = n. (Minimal) state map induced by     1 ξ . . . ξn−1     ❀     w σw . . . σn−1w    

slide-34
SLIDE 34

The axiom of state revisited

A linear system Σ = (T, W, X, Bfull) with latent variable x is a state system if (w, x) ∈ Bfull and x(T) = 0 ⇓ (0, 0) ∧

T (w, x) ∈ Bfull

slide-35
SLIDE 35

The axiom of state revisited

A linear system Σ = (T, W, X, Bfull) with latent variable x is a state system if (w, x) ∈ Bfull and x(T) = 0 ⇓ (0, 0) ∧

T (w, x) ∈ Bfull

  • Time-invariance =

⇒ can choose T = 0;

slide-36
SLIDE 36

The axiom of state revisited

A linear system Σ = (T, W, X, Bfull) with latent variable x is a state system if (w, x) ∈ Bfull and x(T) = 0 ⇓ (0, 0) ∧

T (w, x) ∈ Bfull

  • Time-invariance =

⇒ can choose T = 0;

  • Concatenability with zero trajectory is key.
slide-37
SLIDE 37

When is w ∈ B concatenable with zero?

R0w + R1σw + . . . + RLσLw = 0

. . . w(0) w(1) w(2) w(3) . . . . . .

k = −2 k = −1 k = 0 k = 1 k = 2 k = 3

. . .

slide-38
SLIDE 38

When is w ∈ B concatenable with zero?

R0w + R1σw + . . . + RLσLw = 0

. . . R0 R1 R2 R3 . . . . . . w(0) w(1) w(2) w(3) . . . . . .

k = −2 k = −1 k = 0 k = 1 k = 2 k = 3

. . . R0w(0) + R1w(1) + . . . + RLw(L) = 0

slide-39
SLIDE 39

When is w ∈ B concatenable with zero?

R0w + R1σw + . . . + RLσLw = 0

. . . R0 R1 R2 R3 R4 . . . . . . w(0) w(1) w(2) w(3) . . . . . .

k = −2 k = −1 k = 0 k = 1 k = 2 k = 3

. . . R0w(0) + R1w(1) + . . . + RLw(L) = 0 R1w(0) + R2w(1) + . . . + RLw(L − 1) = 0

slide-40
SLIDE 40

When is w ∈ B concatenable with zero?

R0w + R1σw + . . . + RLσLw = 0

. . . R0 R1 R2 R3 R4 R5 . . . . . . w(0) w(1) w(2) w(3) . . . . . .

k = −2 k = −1 k = 0 k = 1 k = 2 k = 3

. . . R0w(0) + R1w(1) + . . . + RLw(L) = 0 R1w(0) + R2w(1) + . . . + RLw(L − 1) = 0 R2w(0) + R3w(1) + . . . + RLw(L − 2) = 0

slide-41
SLIDE 41

When is w ∈ B concatenable with zero?

R0w + R1σw + . . . + RLσLw = 0

. . . w(0) w(1) w(2) w(3) . . . . . .

k = −2 k = −1 k = 0 k = 1 k = 2 k = 3

. . . R0w(0) + R1w(1) + . . . + RLw(L) = 0 R1w(0) + R2w(1) + . . . + RLw(L − 1) = 0 R2w(0) + R3w(1) + . . . + RLw(L − 2) = 0 . . . = . . .

slide-42
SLIDE 42

When is w ∈ B concatenable with zero?

R0w + R1σw + . . . + RLσLw = 0

. . . RL−3 RL−2 RL−1 RL . . . . . . w(0) w(1) w(2) w(3) . . . . . .

k = −2 k = −1 k = 0 k = 1 k = 2 k = 3

. . . R0w(0) + R1w(1) + . . . + RLw(L) = 0 R1w(0) + R2w(1) + . . . + RLw(L − 1) = 0 R2w(0) + R3w(1) + . . . + RLw(L − 2) = 0 . . . = . . . RL−1w(0) + RLw(1) = 0

slide-43
SLIDE 43

When is w ∈ B concatenable with zero?

R0w + R1σw + . . . + RLσLw = 0

. . . RL−2 RL−1 RL . . . . . . w(0) w(1) w(2) w(3) . . . . . .

k = −2 k = −1 k = 0 k = 1 k = 2 k = 3

. . . R0w(0) + R1w(1) + . . . + RLw(L) = 0 R1w(0) + R2w(1) + . . . + RLw(L − 1) = 0 R2w(0) + R3w(1) + . . . + RLw(L − 2) = 0 . . . = . . . RL−1w(0) + RLw(1) = 0 RLw(0) = 0

slide-44
SLIDE 44

The shift-and-cut map

σ+ : R[ξ] → R[ξ] σ+(n

i=0 piξi) := n−1 i=0 pi+1ξi

“Divide by ξ and take polynomial part" Extended componentwise to vectors and matrices

slide-45
SLIDE 45

Example

R(ξ) = R0 + R1ξ + . . . + RL−1ξL−1 + RLξL

slide-46
SLIDE 46

Example

R(ξ) = R0 + R1ξ + . . . + RL−1ξL−1 + RLξL σ+(R(ξ)) = R1+. . .+RL−1ξL−2+RLξL−1

slide-47
SLIDE 47

Example

R(ξ) = R0 + R1ξ + . . . + RL−1ξL−1 + RLξL σ+(R(ξ)) = R1+. . .+RL−1ξL−2+RLξL−1 σ2

+(R(ξ)) = R2+. . .+RL−1ξL−3+RLξL−2

slide-48
SLIDE 48

Example

R(ξ) = R0 + R1ξ + . . . + RL−1ξL−1 + RLξL σ+(R(ξ)) = R1+. . .+RL−1ξL−2+RLξL−1 σ2

+(R(ξ)) = R2+. . .+RL−1ξL−3+RLξL−2

. . . = . . .

slide-49
SLIDE 49

Example

R(ξ) = R0 + R1ξ + . . . + RL−1ξL−1 + RLξL σ+(R(ξ)) = R1+. . .+RL−1ξL−2+RLξL−1 σ2

+(R(ξ)) = R2+. . .+RL−1ξL−3+RLξL−2

. . . = . . . σL

+(R(ξ)) = RL

slide-50
SLIDE 50

Shift-and-cut and concatenability with zero

w is concatenable with zero ⇔ (σ+(R)(σ)w)(0) = 0 (σ2

+(R)(σ)w)(0)

= 0 . . . = . . . (σL

+(R)(σ)w)(0)

= 0 col((σi

+(R))i=1,...,L(σ) is a state map!

slide-51
SLIDE 51

Shift-and-cut and concatenability with zero

w is concatenable with zero ⇔ (σ+(R)(σ)w)(0) = 0 (σ2

+(R)(σ)w)(0)

= 0 . . . = . . . (σL

+(R)(σ)w)(0)

= 0 col((σi

+(R))i=1,...,L(σ) is a state map!

Other equations equivalent to shift-and-cut ones = ⇒ different state maps are possible!

slide-52
SLIDE 52

Shift-and-cut and concatenability with zero

w is concatenable with zero ⇔ (σ+(R)(σ)w)(0) = 0 (σ2

+(R)(σ)w)(0)

= 0 . . . = . . . (σL

+(R)(σ)w)(0)

= 0 col((σi

+(R))i=1,...,L(σ) is a state map!

Other equations equivalent to shift-and-cut ones = ⇒ different state maps are possible!

slide-53
SLIDE 53

Example: scalar systems

r0w + r1σw + . . . + σnw = 0

slide-54
SLIDE 54

Example: scalar systems

r0w + r1σw + . . . + σnw = 0 Observe w concatenable with zero iff w = 0. Indeed, σn

+(r)(σ)w

= w σn−1

+

(r)(σ)w = rn−1w + σw . . . = . . . σ+(r)(σ)w = r1w + . . . + σn−1w

slide-55
SLIDE 55

Example: scalar systems

r0w + r1σw + . . . + σnw = 0 Observe w concatenable with zero iff w = 0. Indeed, σn

+(r)(σ)w

= w σn−1

+

(r)(σ)w = rn−1w + σw . . . = . . . σ+(r)(σ)w = r1w + . . . + σn−1w Zero at t = 0 iff (σkw)(0) = 0 for k = 0, . . . , n − 1.

slide-56
SLIDE 56

From kernel representation to state map

Denote col((σi

+(R)))i=1,...,L =: ΣR .

Theorem: Let B = ker(R(σ)). Then R(σ)w = ΣR(σ)w = x is a state representation of B with state variable x.

slide-57
SLIDE 57

Algebraic characterization

Theorem: Let B = ker(R(σ)), and define ΣR as

  • above. Then

ΞR := {f ∈ R1×w[ξ] | ∃ g ∈ R1ו[ξ], α ∈ R1ו s.t. f = αΣR + gR} is a vector space over R.

slide-58
SLIDE 58

Algebraic characterization

Theorem: Let B = ker(R(σ)), and define ΣR as

  • above. Then

ΞR := {f ∈ R1×w[ξ] | ∃ g ∈ R1ו[ξ], α ∈ R1ו s.t. f = αΣR + gR} is a vector space over R. X ∈ R•×w[ξ] is state map for B iff row span(X) = ΞR

slide-59
SLIDE 59

Algebraic characterization

Theorem: Let B = ker(R(σ)), and define ΣR as

  • above. Then

ΞR := {f ∈ R1×w[ξ] | ∃ g ∈ R1ו[ξ], α ∈ R1ו s.t. f = αΣR + gR} is a vector space over R. X ∈ R•×w[ξ] is state map for B iff row span(X) = ΞR X is minimal if and only if its rows are a basis for ΞR.

slide-60
SLIDE 60

Example

(σ2 + 2σ + 3)y = (σ + 3)u

  • ξ2 + 2ξ + 3

| −ξ − 3

slide-61
SLIDE 61

Example

(σ2 + 2σ + 3)y = (σ + 3)u

  • ξ2 + 2ξ + 3

| −ξ − 3

  • σ+ ❀

ξ + 2 | −1 ❀ σ + 2 | −1

slide-62
SLIDE 62

Example

(σ2 + 2σ + 3)y = (σ + 3)u

  • ξ2 + 2ξ + 3

| −ξ − 3

  • σ+ ❀

ξ + 2 | −1 ❀ σ + 2 | −1 If (y, u) ∈ B, then for all g ∈ R[ξ]

σ + 2 | −1 y u

  • =

σ + 2 | −1 y u

  • +

g(σ)

  • σ2 + 2σ + 3

| −σ − 3

  • =0 on B
  • y

u

slide-63
SLIDE 63

Example

(σ2 + 2σ + 3)y = (σ + 3)u

  • ξ2 + 2ξ + 3

| −ξ − 3

  • σ+ ❀

ξ + 2 | −1 ❀ σ + 2 | −1 If (y, u) ∈ B, then for all g ∈ R[ξ]

σ + 2 | −1 y u

  • =

σ + 2 | −1 y u

  • +

g(σ)

  • σ2 + 2σ + 3

| −σ − 3

  • =0 on B
  • y

u

  • ‘equivalence modulo R’
slide-64
SLIDE 64

Example

(σ2 + 2σ + 3)y = (σ + 3)u

  • ξ2 + 2ξ + 3

| −ξ − 3

  • σ+ ❀

ξ + 2 | −1 ❀ σ + 2 | −1 σ2

+ ❀

1 | ❀ 1 |

slide-65
SLIDE 65

Example

(σ2 + 2σ + 3)y = (σ + 3)u

  • ξ2 + 2ξ + 3

| −ξ − 3

  • σ+ ❀

ξ + 2 | −1 ❀ σ + 2 | −1 σ2

+ ❀

1 | ❀ 1 | If (y, u) ∈ B, then for all g ∈ R[ξ]

1 | 0 y u

  • =

1 | 0 y u

  • +

g(σ)

  • σ2 + 2σ + 3

| −σ − 3

  • =0 on B

y u

slide-66
SLIDE 66

Example

(σ2 + 2σ + 3)y = (σ + 3)u

  • ξ2 + 2ξ + 3

| −ξ − 3

  • σ+ ❀

ξ + 2 | −1 ❀ σ + 2 | −1 σ2

+ ❀

1 | ❀ 1 | ΞR = {α ξ + 2 | −1 + g(ξ)

  • ξ2 + 2ξ + 3

| −ξ − 3

  • ,

β 1 | + f(ξ)

  • ξ2 + 2ξ + 3

| −ξ − 3

  • α, β ∈ R, f, g ∈ R[ξ]}
slide-67
SLIDE 67

Example

(σ2 + 2σ + 3)y = (σ + 3)u

  • ξ2 + 2ξ + 3

| −ξ − 3

  • σ+ ❀

ξ + 2 | −1 ❀ σ + 2 | −1 σ2

+ ❀

1 | ❀ 1 | ΞR = {α ξ + 2 | −1 + g(ξ)

  • ξ2 + 2ξ + 3

| −ξ − 3

  • ,

β 1 | + f(ξ)

  • ξ2 + 2ξ + 3

| −ξ − 3

  • α, β ∈ R, f, g ∈ R[ξ]}

Any set of generators of ΞR ❀ a state map

slide-68
SLIDE 68

Outline

The axiom of state Discrete-time systems First-order representations State maps The shift-and-cut map Algebraic characterization Continuous-time systems Computation of state-space representations

slide-69
SLIDE 69

On the space of solutions

C∞-solutions to R( d

dt )w = 0 too small ❀ Lloc 1

Equality in the sense of distributions: R( d

dt )w = 0

⇔ +∞

−∞ w(t)⊤(R(− d dt )⊤f)(t)dt = 0

for all testing functions f.

slide-70
SLIDE 70

The axiom of state revisited

Σ = (T, W, X, Bfull) is a state system if (w1, x1), (w2, x2) ∈ Bfull and x1(T) = x2(T) and x1, x2 continuous at T ⇓ (w1, x1) ∧

T (w2, x2) ∈ Bfull

‘State map’ X( d

dt )

slide-71
SLIDE 71

From kernel representation to state map

Denote col((σi

+(R)))i=1,...,L =: ΣR .

Theorem: Let B = ker(R( d

dt )). Then

R( d dt )w = ΣR( d dt )w = x is a state representation of B with state variable x. ¿How to prove it?

slide-72
SLIDE 72

When is w ∈ B concatenable with zero?

0 ∧

0 w ∈ B

⇐ ⇒ +∞

−∞

(0 ∧

0 w)(t)⊤(R(− d

dt )⊤f)(t)dt = 0 ⇐ ⇒ +∞ w(t)⊤(R(− d dt )⊤f)(t)dt = 0 for all testing functions f

slide-73
SLIDE 73

When is w ∈ B concatenable with zero?

0 ∧

0 w ∈ B

⇐ ⇒ +∞

−∞

(0 ∧

0 w)(t)⊤(R(− d

dt )⊤f)(t)dt = 0 ⇐ ⇒ +∞ w(t)⊤(R(− d dt )⊤f)(t)dt = 0 for all testing functions f Integrating repeatedly by parts on f yields: deg(R)

k=1

deg(R)

j=k

(−1)k−1( dj−k

dtj−k w)(0)⊤R⊤ j ( dk−1 dtk−1f)(0)

+ +∞ (R( d

dt )w)(t)⊤f(t)dt = 0

slide-74
SLIDE 74

When is w ∈ B concatenable with zero?

0 ∧

0 w ∈ B

⇐ ⇒ +∞

−∞

(0 ∧

0 w)(t)⊤(R(− d

dt )⊤f)(t)dt = 0 ⇐ ⇒ +∞ w(t)⊤(R(− d dt )⊤f)(t)dt = 0 for all testing functions f Integrating repeatedly by parts on f yields: deg(R)

k=1

deg(R)

j=k

(−1)k−1( dj−k

dtj−k w)(0)⊤R⊤ j ( dk−1 dtk−1f)(0)

+ +∞ (R( d

dt )w)(t)⊤f(t)dt = 0

slide-75
SLIDE 75

w ∈ B concatenable with zero if and only if...

deg(R)

k=1

deg(R)

j=k

(−1)k−1( dj−k

dtj−k w)(0)⊤R⊤ j ( dk−1 dtk−1f)(0) = 0

    f(0) ( d

dt f)(0)

. . . (−1)deg(R)−1( ddeg(R)−1

dtdeg(R)−1f)(0)

    

(ΣR( d

dt )w)(0) = 0

  • (ΣR( d

dt )w)(0) = 0

The shift-and-cut state map!

slide-76
SLIDE 76

Outline

The axiom of state Discrete-time systems First-order representations State maps The shift-and-cut map Algebraic characterization Continuous-time systems Computation of state-space representations

slide-77
SLIDE 77

From kernel representation to state representation

R ∈ Rg×w[ξ] ❀ state map X ∈ Rn×w[ξ] Find: E, F ∈ R(n+g)×n, G ∈ R(n+g)×w T ∈ R(n+g)×g[ξ] with rank(T(λ)) = g ∀λ ∈ C satisfying EξX(ξ) + FX(ξ) + G = T(ξ)R(ξ) Linear equations, Gröbner bases computations!

slide-78
SLIDE 78

From I/O representation to I/O/S representation

I/O representation R = P −Q ❀ state map Xy Xu

  • Find:

A ∈ Rn×n, B ∈ Rn×m, C ∈ Rp×p, D ∈ Rp×m T ∈ R(n+p)×p[ξ] with rank(T(λ)) = g ∀λ ∈ C satisfying

  • ξXy(ξ)

ξXu(ξ) Ip

  • =
  • A

B C D Xy(ξ) Xu(ξ) Im

  • + T(ξ)R(ξ)
slide-79
SLIDE 79

On the choice of state map

State map + system equations ❀ state-space equations

slide-80
SLIDE 80

On the choice of state map

State map + system equations ❀ state-space equations

slide-81
SLIDE 81

On the choice of state map

State map + system equations ❀ state-space equations ( d2

dt2 + 2 d dt + 3)y = ( d dt + 3)u

  • ξ2 + 2ξ + 3

−ξ − 3

slide-82
SLIDE 82

On the choice of state map

State map + system equations ❀ state-space equations ( d2

dt2 + 2 d dt + 3)y = ( d dt + 3)u

  • ξ2 + 2ξ + 3

−ξ − 3

  • Take X(ξ)

=

  • 1

ξ + 2 −1

  • (‘reverse shift-and-cut’).

Then A =

  • −2

1 −3

  • B =
  • −1

−3

  • C =

1 D = ‘observer canonical form’

slide-83
SLIDE 83

On the choice of state map

State map + system equations ❀ state-space equations ( d2

dt2 + 2 d dt + 3)y = ( d dt + 3)u

  • ξ2 + 2ξ + 3

−ξ − 3

  • Take X(ξ) =
  • 1

ξ −1

  • . Then

A =

  • 1

−3 −2

  • B =
  • 1

1

  • C =

1 D = ‘observable canonical form’

slide-84
SLIDE 84

Summary

slide-85
SLIDE 85

Summary

  • The state is constructed!
slide-86
SLIDE 86

Summary

  • The state is constructed!
  • Axiom of state
slide-87
SLIDE 87

Summary

  • The state is constructed!
  • Axiom of state
  • Concatenability with zero
slide-88
SLIDE 88

Summary

  • The state is constructed!
  • Axiom of state
  • Concatenability with zero
  • State maps
slide-89
SLIDE 89

Summary

  • The state is constructed!
  • Axiom of state
  • Concatenability with zero
  • State maps
  • State maps ❀ state-space equations
slide-90
SLIDE 90

Summary

  • The state is constructed!
  • Axiom of state
  • Concatenability with zero
  • State maps
  • State maps ❀ state-space equations
  • Algorithms!