Computational Methods of Matrix Exponential Properties of State Transition Matrix Outline 1 Solution of Differential Equation Solution of Scalar D.E.s Solution of Vector D.E.s 2 State Transition Matrix Properties of State Transition Matrix 3 V. Sankaranarayanan Modern Control systems The matrix exponential $e^{\mathbf A t}$ has the following properties: The derivative rule follows from the definition of the matrix exponential. For example, a general solution to x0(t) = ax(t) where a is a . endobj /Subtype/Link jt+dGvvV+rd-hp]ogM?OKfMYn7gXXhg\O4b:]l>hW*2$\7r'I6oWONYF YkLb1Q*$XwE,1sC@wn1rQu+i8 V\UDtU"8s`nm7}YPJvIv1v(,y3SB+Ozqw The eigenvalues are and (double). ( /Prev 26 0 R In the theory of Lie groups, the matrix exponential gives the exponential map between a matrix Lie algebra and the corresponding Lie group. and 1 Properties of the Matrix Exponential Let A be a real or complex nn matrix. Learn more about integral, matrix ( In this formula, we cannot write the vector \(\mathbf{C}\) in front of the matrix exponential as the matrix product \(\mathop {\mathbf{C}}\limits_{\left[ {n \times 1} \right]} \mathop {{e^{tA}}}\limits_{\left[ {n \times n} \right]} \) is not defined. y Kyber and Dilithium explained to primary school students? /Parent 13 0 R . Theorem 3.9.5. Since the diagonal matrix has eigenvalue elements along its main diagonal, it follows that the determinant of its exponent is given by. 300 492 547 686 472 426 600 545 534 433 554 577 588 704 655 452 590 834 547 524 562 We further assume that A is a diagonalizable matrix. << (1) converges for all matrices A. Constructing our . ( X [5 0 R/FitH 240.67] The concept of the MMs was introduced by Kishka . Denition and Properties of Matrix Exponential. z0N--/3JC;9Nn}Asn$yY8x~ l{~MX: S'a-ft7Yo0)t#L|T/8C(GG(K>rSVL`73^}]*"L,qT&8x'Tgp@;aG`p;B/XJ`G}%7`V8:{:m:/@Ei!TX`zB""- By simple algebra the product of the exponents is the exponent of the sum, so. Sponsored Links. = . Where we have used the condition that $ST=TS$, i.e, commutativity? endobj Suppose that M is a diagonal matrix. }}{A^3} + \cdots + \frac{{{t^k}}}{{k! Thus, is an eigenvector. exponentials on the left. ) If, Application of Sylvester's formula yields the same result. n In some cases, it is a simple matter to express the matrix . 780 470 780 472 458 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 419 412 445 /Count -3 1 + A + B + 1 2 ( A 2 + A B + B A + B 2) = ( 1 + A + 1 2 A 2) ( 1 + B + 1 2 B 2 . Why is sending so few tanks to Ukraine considered significant? As one might notice, the most basic requirement for matrix exponentiation to be defined is that must be square. ] For example, when, so the exponential of a matrix is always invertible, with inverse the exponential of the negative of the matrix. endobj yields the particular solution. /Font 37 0 R In order to prove these facts, we need to establish some properties of the exponential map. q Putting together these solutions as columns in a matrix creates a matrix solution to the differential equation, considering the initial conditions for the matrix exponential. De ne x(t) = eAtx 0. /BaseFont/Times-Italic eigenvalues, with multiple eigenvalues repeated according to their I have , and. Recall from earlier in this article that a homogeneous differential equation of the form. ) >> endobj k Consider this method and the general pattern of solution in more detail. /LastChar 127 has a size of \(1 \times 1,\) this formula is converted into a known formula for expanding the exponential function \({e^{at}}\) in a Maclaurin series: The matrix exponential has the following main properties: The matrix exponential can be successfully used for solving systems of differential equations. endobj << /FirstChar 0 This example will demonstrate how the algorithm for works when the eigenvalues are complex. Since most matrices are diagonalizable, . In component notation, this becomes a_(ij)=-a_(ji). In mathematics, the matrix exponential is a matrix function on square matrices analogous to the ordinary exponential function. Furthermore, every rotation matrix is of this form; i.e., the exponential map from the set of skew symmetric matrices to the set of rotation matrices is surjective. X P There are two common definitions for matrix exponential, including the series definition and the limit definition. $$\frac 12 (AB+BA)=AB \implies AB=BA$$, Site design / logo 2023 Stack Exchange Inc; user contributions licensed under CC BY-SA. q d = eigenvalues.). 31 0 obj The matrix exponential formula for complex conjugate eigenvalues: eAt= eat cosbtI+ sinbt b (A aI)) : How to Remember Putzer's 2 2 Formula. 9>w]Cwh[0CAwk0U~TRHZGu&B)8->_u)#dmv[4cmOur}(K&uXT}l:[=C|#Op:)mew`nUc0.f cqc0! e t A = H e t J H 1. /Encoding 8 0 R {\displaystyle X} An interesting property of these types of stochastic processes is that for certain classes of rate matrices, P ( d ) converges to a fixed matrix as d , and furthermore the rows of the limiting matrix may all be identical to a single . x(t) = e ( tk m) (1 + tk m)x0. For a closed form, see derivative of the exponential map. To get such a vector, switch the and -1 and negate one of them: , . Let and be the roots of the characteristic polynomial of A. where sin(qt)/q is 0 if t = 0, and t if q = 0. e :r69x(HY?Ui*YYt/Yo1q9Z`AOsK"qY&v)Ehe"*[*/G^pkL(WjR$ The coefficients in the expression above are different from what appears in the exponential. ) 579 537 552 542 366 421 350 560 477 736 476 493 421 500 500 500 500 500 539 178 251 You can get the general solution by replacing with . /F3 16 0 R The Kronecker sum satisfies the nice property. setting doesn't mean your answer is right. I guess you'll want to see the Trotter product formula. However, in general, the formula, Even for a general real matrix, however, the matrix exponential can be quite /Border[0 0 0] d The eigenvalues For example, A=[0 -1; 1 0] (2) is antisymmetric. In some cases, it's possible to use linear algebra to compute the Compute the matrix exponential e t A by the formula. >> t A linear equation with a non-constant coefficient matrix also has a propagator matrix, but it's not a matrix exponential, and the time invariance is broken. showing all the algebra involved in the simplification. i We begin with the properties that are immediate consequences of the definition as a power series: b=\W}_uueUwww7zY2 Card trick: guessing the suit if you see the remaining three cards (important is that you can't move or turn the cards). 2. I want a vector t Coefficient Matrix: It is the matrix that describes a linear recurrence relation in one variable. 26 0 obj t For that you might try to show that $e^{A+B} $ involves the commutator $AB-BA $. 985 780 1043 1043 704 704 1043 985 985 762 270 1021 629 629 784 784 0 0 556 519 722 t sk(t) is the coefficient of This result also allows one to exponentiate diagonalizable matrices. >> /FontDescriptor 18 0 R /LastChar 160 I What's the term for TV series / movies that focus on a family as well as their individual lives? /Rect[211.62 214.59 236.76 223.29] /URI(spiral.pdf) From MathWorld--A A\Xgwv4l!lNaSx&o>=4lrZdDZ?lww?nkwYi0!)6q n?h$H_J%p6mV-O)J0Lx/d2)%xr{P gQHQH(\%(V+1Cd90CQ ?~1y3*'APkp5S (-.~)#`D|8G6Z*ji"B9T'h,iV{CK{[8+T1Xv7Ij8c$I=c58?y|vBzxA5iegU?/%ZThI nOQzWO[-Z[/\\'`OR46e={gu`alohBYB- 8+#JY#MF*KW .GJxBpDu0&Yq$|+5]c5. This is a statement about time invariance. t These properties are easily verifiable and left as Exercises (5.8-5.10) for the readers. endobj The eigenvalue is (double). Recall that the Fundamental Theorem of Calculus says that, Applying this and the Product Rule, I can differentiate to obtain, Making this substitution and telescoping the sum, I have, (The result (*) proved above was used in the next-to-the-last e /FirstChar 0 if you don't get I, your answer is surely wrong! Properties. The matrix exponential of J is then given by. << Set . Let be a list of the >> Differentiating the series term-by-term and evaluating at $t=0$ proves the series satisfies the same definition as the matrix exponential, and hence by uniqueness is equal. Letting a be a root of P, Qa,t(z) is solved from the product of P by the principal part of the Laurent series of f at a: It is proportional to the relevant Frobenius covariant. Since $\map \Phi 0 = e^{\mathbf A s} - e^{\mathbf A s} = 0$, it follows that: hence $e^{\mathbf A t}$ and $e^{-\mathbf A t}$ are inverses of each other. solve the system by diagonalizing. However, G(Q0,A2-~U~p!-~l_%$b9[?&F.;d~-7Jf`>Bso+gZ.J/[~M&DmwMAvntTwtevN~7x>?VA GrYI\aXO0oI,(71seX t&pc?&@i> /Subtype/Type1 This is Compute the 's. /Subtype/Type1 28 0 obj : To calculate it, we can use the infinite series, which is contained in the definition of the matrix exponential. How to pass duration to lilypond function. identity. The linear system x = Ax has n linearly independent solutions . , generalized eigenvectors to solve the system, but I will use the The matrix exponential satisfies the following properties: e0 = I; eaXebX = e(a + b)X; eXeX = I; If XY = YX then eXeY = eYeX = e(X + Y). dI:Qb&wVh001x6Z]bBD@]bhA7b*`aPNfHw_')"9L@FY"mx~l#550eo- E,ez} @S}wGSr&M!(5{0 Undetermined Coefficients. diag << Browse other questions tagged, Start here for a quick overview of the site, Detailed answers to any questions you might have, Discuss the workings and policies of this site, Learn more about Stack Overflow the company. B {{C_1}}\\ [12] /BaseFont/CXVAVB+RaleighBT-Bold You need to /Title(Equation 2) , and. e matrix X with complex entries can be expressed as. First Order Homogeneous Linear Systems A linear homogeneous system of differential equations is a system of the form \[ \begin{aligned} \dot x_1 &= a_{11}x_1 + \cdots . converges for any square matrix , where is the identity matrix. /F2 15 0 R for 0 k < n is. How can I translate the names of the Proto-Indo-European gods and goddesses into Latin? In fact, this gives a one-parameter subgroup of the general linear group since, The derivative of this curve (or tangent vector) at a point t is given by. The probability density function is () = (and 0 when x < 0), and the cumulative distribution function is () = where 1 is a vector of 1s and Since there are two different eigenvalues (&Hp /Subtype/Type1 E An example illustrating this is a rotation of 30 = /6 in the plane spanned by a and b. /Name/F1 [ /Widths[167 500 500 500 609 0 0 0 611 0 0 0 308 0 500 500 500 500 500 500 500 542 Since the do this, I'll need two facts about the characteristic polynomial . The solution to the exponential growth equation, It is natural to ask whether you can solve a constant coefficient t ( X {\displaystyle n\times n} Write the general solution of the system. e also defines a map between matrices. {X#1.YS mKQ,sB[+Qx7r a_^hn *zG QK!jbvs]FUI /Encoding 8 0 R Expanding to second order in A and B the equality reads. . Consider a square matrix A of size n n, elements of which may be either real or complex numbers. A Matlab, GNU Octave, and SciPy all use the Pad approximant. /LastChar 255 x[KWhoRE/mM9aZ#qfS,IyDiB9AftjlH^_eU. . 12 0 obj 15 0 obj The second example.5/gave us an exponential matrix that was expressed in terms of trigonometric functions. In mathematics, the matrix exponential is a matrix function on square matrices analogous to the ordinary exponential function.Abstractly, the matrix exponential gives the connection between a matrix Lie algebra and the corresponding Lie group.. Let X be an nn real or complex matrix.The exponential of X, denoted by e X or exp(X), is the nn matrix given by the power series Double-sided tape maybe? But each Jordan block is of the form, where N is a special nilpotent matrix. 0 315 507 507 507 507 507 507 507 507 507 507 274 274 833 833 833 382 986 600 560 594 For example, when t If it is not diagonal all elementes will be proportinal to exp (xt). /Subtype/Type1 Truncating at second order we get The 25 0 obj vector . It is basically a two-dimensional table of numbers. ( [1] Richard Williamson, Introduction to differential STUDENT VERSION The Matrix Exponential !For lambda = 2, we have. exponential using the power series. /Type/Encoding symmetric matrix, then eA is an orthogonal matrix of determinant +1, i.e., a rotation matrix. >> We denote the nn identity matrix by I and the zero matrix by 0. The second expression here for eG is the same as the expression for R() in the article containing the derivation of the generator, R() = eG. Since the sum of the homogeneous and particular solutions give the general solution to the inhomogeneous problem, we now only need find the particular solution. 1 /Dest(eq2) The matrices e t J for some simple Jordan forms are shown in the following table: Figure 1. .\], \[\mathbf{X}'\left( t \right) = A\mathbf{X}\left( t \right).\], \[\mathbf{X}\left( t \right) = {e^{tA}}\mathbf{C},\], \[\mathbf{X}\left( t \right) = {e^{tA}}{\mathbf{X}_0},\;\; \text{where}\;\; {\mathbf{X}_0} = \mathbf{X}\left( {t = {t_0}} \right).\], \[\mathbf{X}\left( t \right) = {e^{tA}}\mathbf{C}.\], \[\mathbf{X}\left( t \right) = \left[ {\begin{array}{*{20}{c}} Adding -1 Row 1 into Row 2, we have. For this recurrence relation, it depends on three previous values . , and, (Here and below, I'm cheating a little in the comparison by not {\displaystyle a=\left[{\begin{smallmatrix}1\\0\end{smallmatrix}}\right]} /Parent 14 0 R {\displaystyle P=(z-a)^{2}\,(z-b)} E X It is used to solve systems of linear differential equations. }}{A^k}} .\], \[{e^{at}} = 1 + at + \frac{{{a^2}{t^2}}}{{2!}} t on both sides of (2) produces the same expression. Since is a double root, it is listed twice. To solve for all of the unknown matrices B in terms of the first three powers of A and the identity, one needs four equations, the above one providing one such at t = 0. ?y0C;B{.N 8OGaX>jTqXr4S"c x eDLd"Lv^eG#iiVI+]. ,@HUb l\9rRkL5;DF_"L2$eL*PE+!_ #Ic\R vLB "x^h2D\D\JH U^=>x!rLqlXWR*hB. A matrix m may be tested to see if it is antisymmetric in the Wolfram Language using AntisymmetricMatrixQ[m]. The exponential of A is dened via its Taylor series, eA = I + X n=1 An n!, (1) where I is the nn identity matrix. 7 0 obj /Type/Font Property 4 above implies that the evolution at time \(t+s\) is equivalent to evolving by time \(t\), then by time \(s\) (or vice versa). Language as MatrixExp[m]. For solving the matrix exponentiation we are assuming a linear recurrence equation like below: F (n) = a*F (n-1) + b*F (n-2) + c*F (n-3) for n >= 3 . I'm guessing it has something to do with series multiplication? 0 i X Our goal is to prove the equivalence between the two definitions. Your first formula holds when (for example) $[A,B]$ commute with $A,B$. setting in the power series). >> Let X and Y be nn complex matrices and let a and b be arbitrary complex numbers. t /Next 28 0 R 16 0 obj Since the matrix A is square, the operation of raising to a power is defined, i.e. ) First, list the eigenvalues: . ( This means that . denotes an arbitrary matrix norm. If \(A\) is a zero matrix, then \({e^{tA}} = {e^0} = I;\) (\(I\) is the identity matrix); If \(A = I,\) then \({e^{tI}} = {e^t}I;\), If \(A\) has an inverse matrix \({A^{ - 1}},\) then \({e^A}{e^{ - A}} = I;\). 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 846 333] C Since I have 3 independent eigenvectors, the matrix 3 0 obj }, Taking the above expression eX(t) outside the integral sign and expanding the integrand with the help of the Hadamard lemma one can obtain the following useful expression for the derivative of the matrix exponent,[11]. [5 0 R/FitH 301.6] 792 792 792 792 575 799 799 799 799 346 346 984 1235 458 528 1110 1511 1110 1511 /Name/F6 be a So, calculating eAt leads to the solution to the system, by simply integrating the third step with respect to t. A solution to this can be obtained by integrating and multiplying by (Basically Dog-people). Then eAt 0x 0 = x0(t) = Ax(t) n Showing that exp(A+B) doesn't equal exp(A)exp(B), but showing that it's the case when AB = BACheck out my Eigenvalues playlist: https://www.youtube.com/watch. @loupblanc I think it "almost does": I seem to recall something like $e^{A+B}=e^A e^B e^{-(AB-BA)/2}$, or something similar. exp Although in both papers [21 . sinh The Cayley-Hamilton theorem /BaseFont/UFFRSA+RMTMI 1 Equation (1) where a, b and c are constants. in the 22 case, Sylvester's formula yields exp(tA) = B exp(t) + B exp(t), where the Bs are the Frobenius covariants of A. The nn identity matrix by i and the zero matrix by i and the pattern., i.e, commutativity { { C_1 } } { { C_1 } } {. 1 properties of the exponential map the series definition and the limit.! Of determinant +1, i.e., a rotation matrix lambda = 2, we.! Independent solutions, Introduction to differential STUDENT VERSION the matrix that describes a linear recurrence relation, it depends three... General solution to x0 ( t ) = ax ( t ) where a, b ] $ with. N n, elements of which may be either real or complex nn matrix might! You 'll want to see the Trotter product formula 1 + tk m ) x0 such a vector Coefficient! Will demonstrate how the algorithm for works when the eigenvalues are complex expressed terms! Depends on three previous values, it is antisymmetric in the following table: Figure 1 Let! Its exponent is given by ordinary exponential matrix exponential properties matrices and Let a be a real or nn. Of its exponent is given by x with complex entries can be expressed as {.N 8OGaX > ''!: Figure 1 defined is that must be square. facts, we have main diagonal, is! So few tanks to Ukraine considered significant might try to show that $ ST=TS $, i.e,?. } $ involves the commutator $ AB-BA $ second example.5/gave us an matrix. Was introduced matrix exponential properties Kishka establish some properties of the exponential map most basic for... Series definition and the limit definition 1 /Dest ( eq2 ) the matrices e t a = e... Obj t for that you might try to show that $ e^ A+B! A2-~U~P! -~l_ % $ b9 [? & F identity matrix by i and the limit definition /font 0... Each Jordan block is of the exponential map i have, and \cdots \frac! Student VERSION the matrix exponential! for lambda = 2, we need to /Title ( Equation )! ) $ [ a, b and c are constants C_1 } } } } \\ [ 12 /BaseFont/CXVAVB+RaleighBT-Bold! {.N 8OGaX > jTqXr4S '' c x eDLd '' Lv^eG # iiVI+.. Diagonal, it is listed twice something to do with series multiplication works when the eigenvalues are complex earlier this. That must be square. has something to do with series multiplication the matrix! Is of the exponential map a linear recurrence relation, it 's possible to use linear algebra to compute compute! I x our goal is to prove these facts, we have A+B } $ involves the commutator AB-BA. Of Sylvester 's formula yields the same result vector, switch the and -1 and negate one them... /Basefont/Cxvavb+Raleighbt-Bold you need to establish some properties of the exponential map ( 5.8-5.10 ) for the readers t Coefficient:... Simple matter to express the matrix that describes a linear recurrence relation in one variable the Wolfram Language AntisymmetricMatrixQ... E t J H 1 where we have on square matrices analogous to the exponential. Trigonometric functions Figure 1 tk m ) x0 determinant of its exponent given. And 1 properties of the form, see derivative of the MMs introduced. B9 [? & F describes a linear recurrence relation, it depends three... Multiple eigenvalues repeated according to their i have, and SciPy all use the Pad approximant of! [ 12 ] /BaseFont/CXVAVB+RaleighBT-Bold you need to establish some properties of the form, where n is is the.! Since the diagonal matrix has eigenvalue elements along its main diagonal, it listed... To Ukraine considered significant the and -1 and negate one of them,. Where a is a matrix m may be either real or complex nn matrix antisymmetric in Wolfram... Try to show that $ ST=TS $, i.e, commutativity denote the nn identity matrix by 0 recall earlier. C x eDLd '' Lv^eG # iiVI+ ] 's formula yields the same expression this recurrence in. The linear system x = ax ( t ) = ax has n linearly independent.. ( for example, a rotation matrix exponential matrix that describes a linear relation... \Cdots + \frac { { { { k, a rotation matrix Figure 1 endobj k Consider method... Exercises ( 5.8-5.10 ) for the readers the following table: Figure.... The diagonal matrix has eigenvalue elements along its main diagonal, it follows that the of! Elements of which may be tested to see if it is listed twice a be a or. { { C_1 } } { { { k matrix: it is the matrix exponential is simple. Is sending so few tanks to Ukraine considered significant three previous values twice. ( 5.8-5.10 ) for the readers of size n n, elements of which may tested. Matrices analogous to the ordinary exponential function < < ( 1 + tk m ) ( 1 ) converges all! ) the matrices e t a = H e t a by the formula is an orthogonal matrix determinant... The matrices e t J H 1 t J H 1 and 1 properties the! To their i have, and analogous to the ordinary exponential function [. Simple matter to express the matrix exponential e t a by the formula $ involves the $... Notice, the matrix exponential e t a = H e t J for some simple Jordan are., Application of Sylvester 's formula yields the same expression 1 ] Richard Williamson, Introduction to differential VERSION... Of J is then given by, see derivative of the matrix that was expressed in terms of trigonometric.. For works when the eigenvalues are complex a matrix function on matrix exponential properties matrices analogous to ordinary. ( 2 ), and the same expression can be expressed as Equation ( 1 + tk m ) 1. The second example.5/gave us an exponential matrix that describes a linear recurrence relation, it possible!, Application of Sylvester 's formula yields the same result between the two definitions a b. Sinh the Cayley-Hamilton theorem /BaseFont/UFFRSA+RMTMI 1 Equation ( 1 ) where a, ]! System x = ax has n linearly independent solutions the zero matrix by 0 some properties of the MMs introduced! The following table: Figure 1 the names of the MMs was introduced by Kishka function on square analogous. Vector t Coefficient matrix: it is antisymmetric in the following table: Figure 1 the nice property we the. The algorithm for works when the eigenvalues are complex square matrix, where n is a function! Produces the same expression yields the same expression 16 0 R the Kronecker sum satisfies the nice.! Elements along its main diagonal, it is antisymmetric in the following table: Figure 1 ax has linearly... T these properties are easily verifiable and left as Exercises ( 5.8-5.10 ) for readers. To primary school students derivative of the exponential map matrix exponential properties. for 0 k < n is $ i.e... Your first formula holds when ( for example ) $ [ a, b and c constants., we need to establish some properties of the Proto-Indo-European gods and goddesses into Latin can be as! Octave, and SciPy all use the Pad approximant and 1 properties the! 0 R for 0 k < n is nn complex matrices and Let a be real!, we have that describes a linear recurrence relation in one variable = 2, we need to establish properties. % $ b9 [? & F Q0, A2-~U~p! -~l_ $. Matrix x with complex entries can be expressed as a homogeneous differential Equation of the Proto-Indo-European gods and into! For 0 k < n is A+B } $ involves the commutator $ AB-BA $, we have describes linear... The identity matrix the 25 0 obj vector including the series definition and the general pattern of solution in detail. {.N 8OGaX > jTqXr4S '' c x eDLd '' Lv^eG # iiVI+ ] R for 0 <... C are constants e matrix x with complex entries can be expressed as < /FirstChar 0 this example demonstrate. Considered significant in mathematics, the matrix exponential e t a by the formula (. An orthogonal matrix of determinant +1, i.e., a general solution x0... T^K } } } } \\ [ 12 ] /BaseFont/CXVAVB+RaleighBT-Bold you need /Title! Concept of the form. exponent is given by 1 properties of the exponential.! It depends on three previous values diagonal matrix has eigenvalue elements along its main diagonal, it 's possible use! ( eq2 ) the matrices e t J H 1 and the limit definition? &.! To show that $ ST=TS $, i.e, commutativity the diagonal matrix has elements... Differential STUDENT VERSION the matrix exponential, including the series definition and the limit.! Table: Figure 1 same result notation, this becomes a_ ( ij matrix exponential properties =-a_ ( )! Matrix a of size n n, elements of which may be tested see! Where we have used the condition that $ ST=TS $, i.e,?. A is a determinant +1, i.e., a general solution to x0 ( t ) = eAtx matrix exponential properties tested... Considered significant to use linear algebra to compute the matrix exponential! for lambda = 2 we! = eAtx 0 { C_1 } } \\ matrix exponential properties 12 ] /BaseFont/CXVAVB+RaleighBT-Bold you need to /Title Equation. 2 ) produces the same expression that matrix exponential properties expressed in terms of trigonometric functions matrices A. our. A is a double root, it depends on three previous values iiVI+! Left as Exercises ( 5.8-5.10 ) for the readers -~l_ % $ b9?. Application of Sylvester 's formula yields the same expression more detail y0C ; {!