Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • SDFSDFSD

    Cargado por

    mojonchico
    43 vistas5 páginas
    El documento resume los métodos de las matrices de banda y el método de Crout para resolver sistemas de ecuaciones lineales tridiagonales. El método de Crout factoriza una matriz tridiagonal A en matrices triangulares L y U de forma que también son tridiagonales, reduciendo así el problema.

    Descripción original:

    fsdfsdfsd

    Título original

    (4)sdfsdfsd

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    El documento resume los métodos de las matrices de banda y el método de Crout para resolver sistemas de ecuaciones lineales tridiagonales. El método de Crout factoriza una matriz tridiagonal A en matrices triangulares L y U de forma que también son tridiagonales, reduciendo así el problema.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Marcar por contenido inapropiado
    Descargar como pdf o txt
    43 vistas5 páginas

    SDFSDFSD

    Cargado por

    mojonchico
    El documento resume los métodos de las matrices de banda y el método de Crout para resolver sistemas de ecuaciones lineales tridiagonales. El método de Crout factoriza una matriz tridiagonal A en matrices triangulares L y U de forma que también son tridiagonales, reduciendo así el problema.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Marcar por contenido inapropiado
    Descargar como pdf o txt
    de 5

    Métodos Numéricos, Sesión 4:

    Contenidos.
    a) Matriz de banda.
    b) Método de Crout.

    1. Matrices de banda.

    Denición 1 Una matriz n × n recibe el nombre de matriz de banda si existen los enteros p y q
    con 1 < p, q < n, con la propiedad de que aij = 0 siempre que p ≤ j − i o q ≤ i − j. El ancho de
    banda o banda de una matriz se dene como w = p + q − 1

    Ejemplo 1 La matriz  
    7 2 0
    A= 3 5 −1 
    0 −5 −6
    es una matriz de banda con p = q = 2 y ancho de banda 2 + 2 − 1 = 3

    1.1. Matrices tridiagonales

    Las matrices de ancho de banda 3 se presentan cuando p = q = 2, reciben el nombre de


    tridiagonales porque tienen la forma
     
    a11 a12 0 ... ... ··· ···
    ... ..
    .
     
     a21 a22 a23 

    ... ... .. 

    0 a32 a33 a34 . 
    A= .. ... ... ... ... ... ...
     
    .

     
    .. ... ... ... ...
     
    .
     
     an−1,n 
    0 . . . . . . . . . . . . an,n−1 ann

    2. Método de Crout

    Los algoritmos de factorización se simplican considerablemente en el caso de matrices de


    banda debido al gran número de ceros que aparecen en estas matrices en patrones regulares. Es

    Métodos Numéricos 1
    especialmente interesante observar la forma que el método Crout o Doolittle asume en este caso.
    Para ilustrar la situación, suponga que una matriz tridiagonal A se puede factorizar en las matrices
    triangulares L y U . Entonces A tiene máximo (3n−2) entradas diferentes a cero. Por lo que existen
    solamente (3n − 2) condiciones a aplicar para determinar las entradas de L y U , siempre y cuando,
    por supuesto, también se obtengan las entradas cero de A. Suponga que las matrices L y U también
    tienen forma tridiagonal; es decir,
     ... 
    l11 0 ... ... ···
     l21 l22 . . . . . . .. 
    . 

    L =  0 ... ... ... ... .. 


     
    . 

     .. . . . . . . . . . ... ... 
     
     .
    ...

    0 ... ... ln,n−1 lnn
     
    1 u12 0 ··· ··· 0
    ... ... ... ..
    .
     
     0 1 
     .. ... ... ... ...

    U =

    . 0

    .. ... ... ...

    .
     
     un−1,n 
    ... ...
     
    0 ··· ··· 1
    El método de Crout consiste en resolver el sistema de ecuaciones lineales de n × n cuya matriz de
    coecientes es tridiagonal.
    a11 = l11
    ai,i−1 = li,i−1 , para cada i = 2, 3, . . . , n
    aii = li,i−1 ui−1,i + lii , para cada i = 2, 3, . . . , n
    y
    ai,i+1 = lii ui,i+1 , para cada i = 1, 2, . . . , n − 1
    Ejemplo 2 Determine la factorización de Crout de la matriz tridiagonal simétrica
     
    2 −1 0 0
     −1 2 −1 0 
     
     0 −1 2 −1 
    0 0 −1 2
    y utilice esta factorización para resolver el sistema lineal
    2x1 − x2 =1
    −x1 + 2x2 − x3 =0
    −x2 + 2x3 − x4 =0
    −x3 + 2x4 =1

    Métodos Numéricos 2
    Solución: La factorización LU de A tiene la forma
        
    a11 0 0 0 l11 0 0 0 1 u12 0 0
     a21 a22 a23 0   l21 l22 0 0   0 1 u 23 0 
    A=
     0 a32 a33 a34  
    =  
    0 l32 l33 0   0 0 1 u34 
    0 0 a43 a44 0 0 l43 l44 0 0 0 1
     
    l11 l11 u12
     l21 l22 + l21 u12 122 u23 0 
     
     0 132 l33 + l32 u23 l33 u34 
    0 0 l43 l44 + l43 u34
    Esto nos da la factorización de Crout:
    1 − 12
        
    2 −1 0 0 2 0 0 0 0 0
     −1 3 2
    2 −1 0   −1
     
    2
    0 0  0
      1 −3 0 
    A= = 4
     = LU
     0 −1 2 −1  0 −1 3
    0   0 0 1 − 43 
    0 0 −1 2 0 0 −1 45 0 0 0 1

    Al resolver el sistema
        
    2 0 0 0 z1 1
     −1 3
    2
    0 0   z2
        0 
    Lz = 
     0 −1 4
    = 
    3
    0   z3   0 
    0 0 −1 54 z4 1

    De donde se obtiene:    1 
    z1 2
     z2   1 
      3
     z3  =  1
     

    4
    z4 1
    luego:
    y al resolver
    1 − 21 1
        
    0 0 x1 2
    2 1
     0 1 − 3
    0   x2  
    3

    Ux =   = 
     0 0 1 − 43   x3   1
    4

    0 0 0 1 x4 1
    Así:    
    x1 1
     x2   1 
     x3  =  1
       

    x4 1

    Métodos Numéricos 3
    Teorema 1 Suponga que A = [aij ] es tridiagonal con ai,i−1 ai,i+1 6= 0, para cada i = 2, 3, . . . , n − 1
    Si |a11 | > |a12 | , |aii | ≥ |ai,i−1 | + |ai,i+1 | , para cada i = 2, 3, . . . , n − 1, y |ann | > |an,n−1 | entonces A
    es no singular y los valores de lii descritos en el algoritmo de factorización Crout son diferentes a
    cero para cada i = 1, 2, . . . , n

    3. Actividades para el estudiante

    3.1 Resolver los siguientes sistemas lineales:


    a)
    x1 − x2 = 0
    −2x1 + 4x2 − 2x3 = −1 (1)
    −x2 + 2x3 = 1,5

    b)
    3x1 + x2 = −1
    2x1 + 4x2 + x3 = 7 (2)
    2x2 + 5x3 = 9

    c)
    2x1 − x2 = 3
    −x1 + 2x2 − x3 = −3 (3)
    −x2 + 2x3 = 1

    d)
    0,5x1 + 0,25x2 = 0,35
    0,35x1 + 0,8x2 + 0,4x3 = 0,77
    (4)
    0,25x2 + x3 + 0,5x4 = −0,5
    x3 − 2x4 = −2,25

    Métodos Numéricos 4
    4. Respuestas

    3.1 a) x1 = 0,5, x2 = 0,5, x3 = 1


    b) x1 = −0,9999995, x2 = 1,999999, x3 = 1
    c) x1 = 1, x2 = −1, x3 = 0
    d) x1 = −0,09357798, x2 = 1,587156, x3 = −1,167431, x4 = 0,5412844

    Métodos Numéricos 5

    También podría gustarte