Calculating A k using Fulmers Method Rasheen Alexander, Katie - - PowerPoint PPT Presentation

calculating a k using fulmer s method
SMART_READER_LITE
LIVE PREVIEW

Calculating A k using Fulmers Method Rasheen Alexander, Katie - - PowerPoint PPT Presentation

Introduction Examples Proving why Fulmers Method works for A k Example of e At from A k Summary Calculating A k using Fulmers Method Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker July 12, 2013 Calculating A k using Fulmers


slide-1
SLIDE 1

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary

Calculating Ak using Fulmer’s Method

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker July 12, 2013

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-2
SLIDE 2

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary

Table of Contents

1

Introduction Why? Definitions

2

Examples Partial Fractions Decomposition Fulmer’s Method

3

Proving why Fulmer’s Method works for Ak Linear Independence Generalizing to an n × n matrix

4

Example of eAt from Ak

5

Summary

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-3
SLIDE 3

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Why find Ak and eAt?

Ak is essential to find the solutions to difference equations. Calculating eAt, the matrix exponential. eAt is used in solving matrix linear differential equations.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-4
SLIDE 4

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Definition Let n be a nonnegative integer. The falling factorial is the sequence kn, with k = 0, 1, 2, . . . given by the following formula.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-5
SLIDE 5

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Definition Let n be a nonnegative integer. The falling factorial is the sequence kn, with k = 0, 1, 2, . . . given by the following formula. kn = k(k − 1)(k − 2) · · · (k − n + 1). If k were allowed to be a real variable then kn could be characterized as the unique monic polynomial of degree n that vanishes at 0, 1, . . . , n − 1. Observe also that kn|k=n = n!.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-6
SLIDE 6

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Definition Let n be a nonnegative integer and a be a complex number. We define a sequence ϕn,a(k) as ϕn,a(k) =

  • ak−nkn

n!

a = 0, δn(k) a = 0.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-7
SLIDE 7

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Definition Let n be a nonnegative integer and a be a complex number. We define a sequence ϕn,a(k) as ϕn,a(k) =

  • ak−nkn

n!

a = 0, δn(k) a = 0. where δn(k) =

  • k = n

1 k = n

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-8
SLIDE 8

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Example With ϕ0,0(k), ϕ1,4(k), and ϕ2,5(k), we have

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-9
SLIDE 9

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Example With ϕ0,0(k), ϕ1,4(k), and ϕ2,5(k), we have ϕ0,0(k) = δ0(k) = (1, 0, 0, 0, 0, . . .)

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-10
SLIDE 10

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Example With ϕ0,0(k), ϕ1,4(k), and ϕ2,5(k), we have ϕ0,0(k) = δ0(k) = (1, 0, 0, 0, 0, . . .) ϕ1,4(k) = 4k−1k = (0, 1, 8, 48, 256, . . .)

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-11
SLIDE 11

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Example With ϕ0,0(k), ϕ1,4(k), and ϕ2,5(k), we have ϕ0,0(k) = δ0(k) = (1, 0, 0, 0, 0, . . .) ϕ1,4(k) = 4k−1k = (0, 1, 8, 48, 256, . . .) ϕ2,5(k) = 5k−2k(k−1)

2

= (0, 0, 1, 15, 150, . . .)

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-12
SLIDE 12

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Definition Let y(k) be a sequence of complex numbers. We define the Z-transform of y(k) to be the function Z{y(k)}(z), where z is a complex variable, by the following formula: Z{y(k)}(z) =

  • k=0

y(k) zk

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-13
SLIDE 13

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Now that we’ve defined the Z-Transform, we can now apply it to ϕn,a. Let a ∈ C and n ∈ N and we obtain the following formula. Z{ϕn,a(k)}(z) = z (z − a)n+1

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-14
SLIDE 14

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Theorem Let A be an n x n matrix with entries in the complex plane. Then Z{Ak}(z) = z(zI − A)−1 where I is the n x n identity matrix.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-15
SLIDE 15

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Why? Definitions

Note that the Z-Transform is one-to-one and linear. Therefore, the Z-Transform has an inverse. Now that we know that the Z-Transform is invertable we obtain the following formula Ak = Z−1{z(zI − A)−1}

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-16
SLIDE 16

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Example Find Ak if A = 2 −1 1

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-17
SLIDE 17

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Recall Z{Ak}(z) = z(zI − A)−1. First, we would compute zI − A. zI − A = z − 2 1 −1 z

  • Next, compute the inverse of zI − A.

(zI − A)−1 = 1 z2 − 2z + 1 z −1 1 z − 2

  • =

1 (z − 1)2 z −1 1 z − 2

  • =
  • z

(z−1)2 −1 (z−1)2 1 (z−1)2 z−2 (z−1)2

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-18
SLIDE 18

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

So we must perform partial fraction decomposition to obtain z (z − 1)2 = A (z − 1) + B (z − 1)2 z = A(z − 1) + B If z = 1, 1 = A(1 − 1) + B 1 = A(0) + B 1 = 0 + B B = 1

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-19
SLIDE 19

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

z (z − 1)2 = A (z − 1) + B (z − 1)2 z = A(z − 1) + B If z = 0, 0 = A(0 − 1) + B 0 = A(−1) + B 0 = −A + B A = B A = 1

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-20
SLIDE 20

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Plug in results: z (z − 1)2 = 1 (z − 1) + 1 (z − 1)2

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-21
SLIDE 21

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

z − 2 (z − 1)2 = A (z − 1) + B (z − 1)2 z − 2 = A(z − 1) + B If z = 1, 1 − 2 = A(1 − 1) + B −1 = A(0) + B −1 = 0 + B B = −1

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-22
SLIDE 22

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

z − 2 (z − 1)2 = A (z − 1) + B (z − 1)2 z − 2 = A(z − 1) + B If z = 0, 0 − 2 = A(0 − 1) + B −2 = A(−1) + B −2 = −A + B 2 + B = A 2 + (−1) = A 1 = A

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-23
SLIDE 23

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Plug in results: z − 2 (z − 1)2 = 1 (z − 1) + −1 (z − 1)2 (zI − A)−1 =

  • 1

(z−1) + 1 (z−1)2 −1 (z−1)2 1 (z−1)2 1 (z−1) + −1 (z−1)2

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-24
SLIDE 24

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Finally, we obtain the following: (zI − A)−1 = 1 z − 1 1 1

  • +

1 (z − 1)2 1 −1 1 −1

  • which separates the original matrix into two matrices that have a

common denominator for each entry. Multiplying z into the equation, we obtain z(zI − A)−1 = z z − 1 1 1

  • +

z (z − 1)2 1 −1 1 −1

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-25
SLIDE 25

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

By previous definition and taking the inverse Z-transform, we find that ϕn,a(k) = Z−1

  • z

(z − a)n+1

  • and since we have fractions that are of this form, we may apply the

inverse Z-transform to obtain Ak = Z−1

  • z

z − 1 1 1

  • +

z (z − 1)2 1 −1 1 −1

  • = ϕ0,1(k)

1 1

  • + ϕ1,1(k)

1 −1 1 −1

  • = 1k−0k 0

0! 1 1

  • + 1k−1k 1

1! 1 −1 1 −1

  • =

1 1

  • + k

1 −1 1 −1

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-26
SLIDE 26

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Therefore, we obtain that Ak = 1 1

  • + k

1 −1 1 −1

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-27
SLIDE 27

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Example Find Ak by using the Fulmer’s Method if A = 2 −1 1

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-28
SLIDE 28

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Recall Ak = Z−1{z(zI − A)−1}. We first want to compute the determinant of z(zI − A)−1. We start by computing zI − A. zI − A = z − 2 1 −1 z

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-29
SLIDE 29

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Recall Ak = Z−1{z(zI − A)−1}. We first want to compute the determinant of z(zI − A)−1. We start by computing zI − A. zI − A = z − 2 1 −1 z

  • Next,

det(zI − A) = (z − 1)2

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-30
SLIDE 30

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Recall Ak = Z−1{z(zI − A)−1}. We first want to compute the determinant of z(zI − A)−1. We start by computing zI − A. zI − A = z − 2 1 −1 z

  • Next,

det(zI − A) = (z − 1)2 this implies det((zI − A)−1) = 1 (z − 1)2

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-31
SLIDE 31

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Recall Ak = Z−1{z(zI − A)−1}. We first want to compute the determinant of z(zI − A)−1. We start by computing zI − A. zI − A = z − 2 1 −1 z

  • Next,

det(zI − A) = (z − 1)2 this implies det((zI − A)−1) = 1 (z − 1)2 Lastly, we must multiply both sides by z z det((zI − A)−1) = z (z − 1)2

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-32
SLIDE 32

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

We can write Ak in the following way Ak = ϕ0,1(k)M + ϕ1,1(k)N. where M and N are our unknown matrices. By the way we defined the ϕ function, we obtain Ak = 1k−0k0 0! M + 1k−1k1 1! N = M + kN

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-33
SLIDE 33

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Let k = 0 I = M

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-34
SLIDE 34

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Let k = 0 I = M Let k = 1 A = M + N

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-35
SLIDE 35

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Let k = 0 I = M Let k = 1 A = M + N Solve for M and N. From the first equation, we get M = 1 1

  • From the second equation, we obtain

N = A − M = 2 −1 1

1 1

  • =

1 −1 1 −1

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-36
SLIDE 36

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Partial Fractions Decomposition Fulmer’s Method

Finally, Ak = 1 1

  • + k

1 −1 1 −1

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-37
SLIDE 37

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Theorem The standard basis Bq = {ϕ1λ1, ϕ1λ2, · · · , ϕ2λ1, ϕ2λ2, · · · ϕrλm} is linearly independent.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-38
SLIDE 38

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Write Bq as a linear combination. a11ϕ1λ1 + a12ϕ1λ2 + · · · + a1mϕ1λm+ a21ϕ2λ1 + a22ϕ2λ2 + · · · + a2mϕ2λm+ . . . ar1ϕrλ1 + ar2ϕrλ2 + · · · + armϕrλm = 0

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-39
SLIDE 39

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

a11 z (z − λ1)2 + a12 z (z − λ2)2 + · · · + a1m z (z − λm)2 + a21 z (z − λ1)3 + a22 z (z − λ2)3 + · · · + a2m z (z − λm)3 + . . . ar1 z (z − λ1)r+1 + ar2 z (z − λ2)r+1 + · · · + arm z (z − λm)r+1 = 0

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-40
SLIDE 40

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Regroup to form like denominators a11(z(z − λ1)r−1) + a21(z(z − λ1)r−2) + · · · + ar1(z) (z − λ1)r+1 + a12(z(z − λ2)r−1) + a22(z(z − λ2)r−2) + · · · + ar2(z) (z − λ2)r+1 + . . . a1m(z(z − λm)r−1) + a2m(z(z − λm)r−2) + · · · + arm(z) (z − λm)r+1 = 0

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-41
SLIDE 41

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Let nb(z − λb) be a polynomial. n1(z − λ1) (z − λ1)r+1 + n2(z − λ2) (z − λ2)r+1 + · · · + nm(z − λm) (z − λm)r+1 = 0

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-42
SLIDE 42

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Let nb(z − λb) be a polynomial. n1(z − λ1) (z − λ1)r+1 + n2(z − λ2) (z − λ2)r+1 + · · · + nm(z − λm) (z − λm)r+1 = 0 If n1(z − λ1) = 0, then lim

z→λ1

n1(z − λ1) (z − λ1)r+1 + n2(z − λ2) (z − λ2)r+1 + · · · + nm(z − λm) (z − λm)r+1

  • = ∞+C = 0

where C is a constant. Thus, we get a contradiction.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-43
SLIDE 43

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Therefore, n1(z − λ1) = 0 which implies that n1(z) = 0 and a11, a12, · · · , a1m = 0. You can continue this argument by induction to obtain ∀a′s = 0

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-44
SLIDE 44

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Therefore, n1(z − λ1) = 0 which implies that n1(z) = 0 and a11, a12, · · · , a1m = 0. You can continue this argument by induction to obtain ∀a′s = 0 With this Theorem we may now show that we may apply Fulmer’s method to any n × n matrix.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-45
SLIDE 45

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Definition Let n be a natural number and a(k) to be a sequence with complex

  • terms. We define E as a shift operator for sequences where

E n{a(k)} = a(k + n) Example E(2k) = 2k+1 E 2(2k) = 2k+2 E 3(2k) = 2k+3

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-46
SLIDE 46

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Can this work for any n × n matrix?

From Thai’s paper, we know that Ak =

r

  • a = 1

ma−1

  • n = 0

Mna ϕn,ar (k)

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-47
SLIDE 47

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Can this work for any n × n matrix?

From Thai’s paper, we know that Ak =

r

  • a = 1

ma−1

  • n = 0

Mna ϕn,ar (k) However, because there are a finite amount of terms, we may drop the double subscript in favor of a single subscript.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-48
SLIDE 48

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Can this work for any n × n matrix?

From Thai’s paper, we know that Ak =

r

  • a = 1

ma−1

  • n = 0

Mna ϕn,ar (k) However, because there are a finite amount of terms, we may drop the double subscript in favor of a single subscript. Ak =

R

  • n = 1

Mn ϕn(k)

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-49
SLIDE 49

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Can this work for any n × n matrix?

Knowing what the phi sequences are, we may create a system of equations with the shift operator. Setting k = 0, we have Ak = I = M1ϕ1(0) + M2ϕ2(0) + M3ϕ3(0) + . . . E{Ak} = A = M1ϕ1(1) + M2ϕ2(1) + M3ϕ3(1) + . . . E 2{Ak} = A2 = M1ϕ1(2) + M2ϕ2(2) + M3ϕ3(2) + . . . . . . = . . . = . . . E R−1{Ak} = AR−1 = M1ϕ1(R − 1) + M2ϕ2(R − 1) + . . .

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-50
SLIDE 50

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Can this work for any n × n matrix?

Knowing what the phi sequences are, we may create a system of equations with the shift operator. Setting k = 0, we have Ak = I = M1ϕ1(0) + M2ϕ2(0) + M3ϕ3(0) + . . . E{Ak} = A = M1ϕ1(1) + M2ϕ2(1) + M3ϕ3(1) + . . . E 2{Ak} = A2 = M1ϕ1(2) + M2ϕ2(2) + M3ϕ3(2) + . . . . . . = . . . = . . . E R−1{Ak} = AR−1 = M1ϕ1(R − 1) + M2ϕ2(R − 1) + . . .

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-51
SLIDE 51

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Can this work for any n × n matrix?

Representing as a matrix equation, we have        I A A2 . . . AR−1        =        ϕ1(0) ϕ2(0) ϕ3(0) . . . ϕR(0) ϕ1(1) ϕ2(1) ϕ3(1) . . . ϕR(1) ϕ1(2) ϕ2(2) ϕ3(2) . . . ϕR(2) . . . . . . . . . . . . . . . ϕ1(R − 1) ϕ2(R − 1) ϕ3(3) . . . ϕR(R − 1)               M1 M2 M3 . . . MR        Where        M1 M2 M3 . . . MR        is our unknown matrix to solve for.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-52
SLIDE 52

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Can this work for any n × n matrix?

We will let B equal to        ϕ1(0) ϕ2(0) ϕ3(0) . . . ϕR(0) ϕ1(1) ϕ2(1) ϕ3(1) . . . ϕR(1) ϕ1(2) ϕ2(2) ϕ3(2) . . . ϕR(2) . . . . . . . . . . . . . . . ϕ1(R − 1) ϕ2(R − 1) ϕ3(3) . . . ϕR(R − 1)        B is a special matrix called the Matrix of Casorati, where the matrix is made from a set of functions and their E shift. Its determinant is called the Casoratian If we want to find a unique answer for our coefficient matrices, B must be invertible or the determinant of B must be non-zero.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-53
SLIDE 53

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Can this work for any n × n matrix?

Theorem Let u1(t), u2(t), . . . , un(t) be a set of functions. If these functions are linearly independent, then the Casoratian is non-zero for all t.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-54
SLIDE 54

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary Linear Independence Generalizing to an n × n matrix

Can this work for any n × n matrix?

Knowing that the phi functions are linearly independent, we know we find that the determinant of B is non-zero. Therefore, B is invertible and we may find the set of unknown coefficient matrices. Therefore, we may use Fulmer’s Method for any n × n matrix. B−1        I A A2 . . . AR−1        =        M1 M2 M3 . . . MR       

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-55
SLIDE 55

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary

From Davidson’s class, we learned how to solve eAt by the Laplace Transform. A = 2 −1 1

  • eAt = et

1 1

  • + tet

1 −1 1 −1

  • Now, we will show a quicker way to find eAt by knowing Ak from

the Z-Transform.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-56
SLIDE 56

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary

Find eAt if A = 2 −1 1

  • From a previous example, we found that

Ak = 1 1

  • + k

1 −1 1 −1

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-57
SLIDE 57

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary

Next, we can manipulate this equation by adding a summation, multiplying by tk, and dividing by k! to both sides. We get the following:

  • k=0

tkAk k! =

  • k=0

tk 1 1

  • k!

+

  • k=0

ktk 1 −1 1 −1

  • k!

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-58
SLIDE 58

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary

From calculus, we know we can rewrite this as: eAt =

  • k=0

tkAk k! =

  • k=0

tk 1 1

  • k!

+

  • k=0

ktk 1 −1 1 −1

  • k!

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-59
SLIDE 59

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary

Since the matrices are not dependent on k we can move them out

  • f the summation.

eAt = 1 1 ∞

  • k=0

tk k! + 1 −1 1 −1 ∞

  • k=0

ktk k! = 1 1 ∞

  • k=0

tk k! + 1 −1 1 −1 ∞

  • k=1

tk (k − 1)!

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method

slide-60
SLIDE 60

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary

Again, from calculus we know we can rewrite the summations as a Taylor Series, giving us the following equation: eAt = et 1 1

  • + tet

1 −1 1 −1

  • Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker

Calculating Ak using Fulmer’s Method

slide-61
SLIDE 61

Introduction Examples Proving why Fulmer’s Method works for Ak Example of eAt from Ak Summary

Summary

Using the ϕ functions and the Z-transform, we found that Fulmer’s method may be used to calculate Ak. Using the formula for Ak, we may calculate eAt.

Rasheen Alexander, Katie Huston, Thomas Le, Camera Whicker Calculating Ak using Fulmer’s Method