Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 34 vistas

    Ejercicio-10 Literal B

    Cargado por

    Cinthya Narvaez
    El documento describe cómo calcular la exponencial de una matriz usando dos métodos: 1) La matriz A es diagonizable, por lo que existe una matriz P tal que P-1AP es una matriz diagonal D. 2) El algoritmo de Putzer, el cual calcula la exponencial de A como una suma de polinomios de A evaluados en los valores propios de A.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd

    Ejercicio-10 Literal B

    Cargado por

    Cinthya Narvaez
    34 vistas4 páginas
    El documento describe cómo calcular la exponencial de una matriz usando dos métodos: 1) La matriz A es diagonizable, por lo que existe una matriz P tal que P-1AP es una matriz diagonal D. 2) El algoritmo de Putzer, el cual calcula la exponencial de A como una suma de polinomios de A evaluados en los valores propios de A.

    Título original

    ejercicio-10 literal b

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    DOCX, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    El documento describe cómo calcular la exponencial de una matriz usando dos métodos: 1) La matriz A es diagonizable, por lo que existe una matriz P tal que P-1AP es una matriz diagonal D. 2) El algoritmo de Putzer, el cual calcula la exponencial de A como una suma de polinomios de A evaluados en los valores propios de A.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    34 vistas4 páginas

    Ejercicio-10 Literal B

    Cargado por

    Cinthya Narvaez
    El documento describe cómo calcular la exponencial de una matriz usando dos métodos: 1) La matriz A es diagonizable, por lo que existe una matriz P tal que P-1AP es una matriz diagonal D. 2) El algoritmo de Putzer, el cual calcula la exponencial de A como una suma de polinomios de A evaluados en los valores propios de A.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    Está en la página 1de 4

    b.

    Hallar la exponencial de la matriz A usando


    i) El hecho de que la matriz A es diagonizable; es decir que existe la matriz P tal que
    P−1 AP=D.
    ii) El algoritmo de Putzer

    Solución

    i) La matriz A es diagonizable. Es decir que existe una matriz P es singular la que P−1 AP , es decir
    P−1 AP=D. La matriz P tiene como vectores columna a los vectores propios de la matriz A
    calculados en (a), es decir:

    1 −1 −1
    P= 2 4
    1 1( 1
    1 )
    Así la matriz inversa de P es:

    1 1
    0

    ( )
    2 2
    −1 −1 −1
    P−1=
    6 3 2
    −1 −1
    1
    3 3

    1 1
    0

    ( )
    2 2
    1 −1 4 1 −1 −1
    −1
    P AP=
    −1
    6
    −1
    ( −1
    3
    −1
    −1
    2
    1
    3 2 −1 2 4
    2 1 −1 1 1 )(
    1
    1 )
    3 3

    3 0 0
    −1
    P AP= 0 1 0
    0 0 −2 ( )
    Así también se puede verificar que:

    A=PD P−1
    ii) El algoritmo de Putzer

    Del ítem (a) se tiene que los valores propios son λ 1=3 , λ2=1 , λ3=−2 aplicando el algoritmo de
    Putzer se tiene que:

    Po ( A )=I
    −2 −1 4
    P1 ( A )=II 1
    m=1 ( A−λmI )= A− λ1 I = 3 −1 −1
    2 1 −4 ( )
    P2 ( A )=II 2m=1 ( A−λmI )=( A−λ 1 I )( A−λ2 I )

    −2 −1 4 0 −1 4
    (
    ¿ 3 −1 −1 3 1 −1
    2 1 −4 2 1 −2 )( )
    5 5 −15
    (
    ¿ −5 −5 15
    2 1 15 )
    Por otra parte:

    r '1 ( t ) =λ1 r 1 ( t ) ; r 1 ( o )=1

    r '1 ( t )
    =3 ⇒d ( ln ( r 1 ( t ) ) )=3
    r 1 (t)

    ⇒ ln ( r 1 ( t ) )=3 t

    ⇒ r 1 ( t )=e 3 t

    ⇒ r 1 ( 0 )=1

    r '2 ( t ) =λ1 r 2 ( t ) +r 1 ( t ) ; r 1 ( o )=0

    r '2 ( t ) =r 2 ( t ) +e 3 t

    r '2 ( t ) −r 2 ( t )=e3 t

    e−t r '2 (t )¿ e−t r 2 ( t ) =e 2t


    d −t
    [ e r 2 ( t ) ] =e 2t
    dt

    Integrando:

    1
    e−t r 2 ( t ) = e2 t +C
    2
    1
    ⇒ r 2 ( t )= e3 t +C e t
    2
    1 −1
    ⇒ r 2 ( 0 ) = +C ⇒ C=
    2 2

    e3 t −e t
    ⇒ r 2 ( t )=
    2

    r '3 ( t ) =λ3 r 3 (t ) +r 2 ( t ) ; r 3 ( 0 )=0

    ' e3 t −e t
    r 3 ( t ) +2 r 3 ( t ) =
    2

    2t ' 2t e 5 t −e3 t
    e r 3 (t) ¿ 2 e r 3 ( t )=
    2

    d 2t ( ) e 5 t−e3 t
    dt
    [ e r 3 t ]= 2

    Integrando:

    2t e 5 t e3 t
    e r 3 ( t )= − +C
    10 6
    5t 3t
    e e
    ⇒ r 3 ( t )= − +C
    10 6

    e3 t e t −3 t
    ⇒ r3 ( 0 ) = − +C e
    10 6
    1 1 1
    ⇒ r3 ( 0 ) = − +C ⇒ C=
    10 6 15

    e3 t e t 1 −3 t
    ⇒ r 3 ( t )= − + e
    10 6 15
    Luego:
    2
    e tA =∑ r k=1 ( t ) P ( A ) =r 1 ( t ) P 0 ( A ) +r 2 ( t ) P1 ( A ) +r 3 (t ) P 2 ( A )
    k=0

    1 0 0 −1 −1 4 5 5 −15
    e3 t −e t e 3 t et 1 −3 t
    tA
    e =e 3t

    ( 0 1 0 +
    0 0 1 )2
    3 −1 −1 +
    2 1 −4 ( − + e
    10 6 15 )(
    −5 −5 15
    −5 −5 15
    )( )

    También podría gustarte