Scribd
Carregar
10 visualizações2 páginas

Fatoração LU e Resolução de Sistema Linear

O documento descreve a resolução de um sistema linear utilizando a fatoração LU, onde a matriz A é decomposta em uma matriz triangular inferior L e uma matriz triangular superior U. O processo inclui a eliminação gaussiana para calcular os multiplicadores e resolver o sistema em duas etapas: substituição direta e substituição reversa. A solução final do sistema é apresentada como um vetor d com valores aproximados.

Enviado por

Caio Vinicius Pessoa Freires
Direitos autorais
© All Rights Reserved
Levamos muito a sério os direitos de conteúdo. Se você suspeita que este conteúdo é seu, reivindique-o aqui.
Formatos disponíveis
Baixe no formato PDF, TXT ou leia on-line no Scribd
10 visualizações2 páginas

Fatoração LU e Resolução de Sistema Linear

O documento descreve a resolução de um sistema linear utilizando a fatoração LU, onde a matriz A é decomposta em uma matriz triangular inferior L e uma matriz triangular superior U. O processo inclui a eliminação gaussiana para calcular os multiplicadores e resolver o sistema em duas etapas: substituição direta e substituição reversa. A solução final do sistema é apresentada como um vetor d com valores aproximados.

Enviado por

Caio Vinicius Pessoa Freires
Direitos autorais
© All Rights Reserved
Levamos muito a sério os direitos de conteúdo. Se você suspeita que este conteúdo é seu, reivindique-o aqui.
Formatos disponíveis
Baixe no formato PDF, TXT ou leia on-line no Scribd

Sistema linear inicial: Ad = b

   
3 −2 1 8
A = 1 −3 4 , b=6 (1)
9 4 −5 11

Passo 1: Fatoração A = LU (Eliminação Gaussiana)


Matriz inicial:  
3 −2 1
 1 −3 4  (2)
9 4 −5
Cálculo dos multiplicadores da 1ª etapa:
1 1
L′2 = L2 − L1 =⇒ l21 =
3 3 (3)
L′3 = L3 − 3L1 =⇒ l31 = 3

Após a 1ª etapa de eliminação:


 
3 −2 1
0 −7/3 11/3 (4)
0 10 −8

Cálculo do multiplicador da 2ª etapa:


 
′′ ′ 10 30
L3 = L3 − L′2 =⇒ l32 = − (5)
−7/3 7

Matriz final U (Triangular Superior):


 
3 −2 1
U = 0 −7/3 11/3 (6)
0 0 54/7

Matriz L (Triangular Inferior com os multiplicadores):


   
1 0 0 1 0 0
L = l21 1 0 = 1/3 1 0 (7)
l31 l32 1 3 −30/7 1

Passo 2: Resolução de LU d = b em duas etapas


2a. Resolver Ly = b (Substituição Direta)
    
1 0 0 y1 8
1/3 1 0 y2  =  6  (8)
3 −30/7 1 y3 11

1
Cálculo de y:

y1 = 8
 
1 10
8 + y2 = 6 =⇒ y2 =
3 3 (9)
 
30 10 9
3(8) − + y3 = 11 =⇒ y3 =
7 3 7

2b. Resolver U d = y (Substituição Reversa)


    
3 −2 1 d1 8
0 −7/3 11/3 d2  = 10/3 (10)
0 0 54/7 d3 9/7

Cálculo de d:
54 9 1
d3 = =⇒ d3 =
7
  7 6
7 11 1 10 7
− d2 + = =⇒ d2 = − (11)
3 3 6 3 6
 
7 1 11
3d1 − 2 − + = 8 =⇒ d1 =
6 6 6

Solução Final:    
11/6 1.833
d = −7/6 ≈ −1.167 (12)
1/6 0.167

Você também pode gostar