Teoria

Lembra de quando vimos as técnicas de pivoteamento que iriam resolver aquele problema da eliminação de Gauss?

Então, na fatoração L U nós precisamos aplicar a eliminação de Gauss, ou seja, novamente podemos ter o problema de um pivô ser igual ou muito próximo de zero! E é aí que entra a fatoração P L U , que é simplesmente a fatoração L U com o pivoteamento!

A única diferença é que só vamos estudar o pivoteamento parcial (que bom né!). O que acontece agora é que teremos uma terceira matriz P que guardará essa informação da troca de linhas. Daí o nome fatoração P L U .

E como obtemos essa matriz P ? É bem simples na verdade, é só pegarmos a matriz identidade e fazermos as mesmas trocas de linha que fizermos durante a eliminação de Gauss só que de trás pra frente. É como se tivéssemos desfazendo as trocas. Vamos ver um exemplo 3 × 3 :

1 3 4 5 2 7 - 2 1 - 9

Aplicando o pivoteamento parcial, fazemos L 1 ↔ L 2 :

5 2 7 1 3 4 - 2 1 - 9

E a matriz de pivoteamento P será dada por:

0 1 0 1 0 0 0 0 1

Que é exatamente a matriz identidade com as primeiras duas linhas trocadas. Se mais alguma troca tivesse sido feita depois teríamos de aplicar essa última troca primeiro para dar certo!

Então, seja o sistema linear A x = b , temos:

P L U x = b

L U x = P - 1 b

Onde P - 1 é a inversa de P , que é obtida simplesmente fazendo as mudanças de linhas na matriz identidade na mesma ordem em que fizemos no escalonamento. Ou seja, para obtermos a matriz P da fatoração P L U aplicamos as mudanças de linha na ordem inversa. Se quisermos resolver o sistema encontramos P - 1 aplicando as mudanças de linha na ordem certa. Assim, precisamos resolver dois sistemas triangulares:

L y = P - 1 b

U x = y

E assim vamos obter a solução x do nosso sistema.

Obs.: O chato aqui é que se forem feitas trocas de linhas no depois de já termos definido parte da matriz L , esta será alterada. Então eu recomendo ir completando a matriz L junto da eliminação de Gauss, só não pode aplicar o escalonamento nos multiplicadores sem querer! Vai ficar mais claro no exemplo:

Encontre a fatoração P L U do sistema abaixo:

x 1 + x 2 + x 3 = 1 x 1 - x 2 + x 3 = 2 2 x 1 + x 2 - x 3 = 3

Passando para a forma matricial:

1 1 1 1 - 1 1 2 1 - 1 x 1 x 2 x 3 = 1 2 3

Aplicando a eliminação de Gauss com pivoteamento parcial na matriz A :

1 1 1 1 - 1 1 2 1 - 1

L 1 ↔ L 3 :

2 1 - 1 1 - 1 1 1 1 1

L 2 ← L 2 - 1 2 L 1 m 21 = 1 2 :

2 1 - 1 { 0.5 } - 1.5 1.5 1 1 1

Entendeu o que eu quis dizer agora? Nós vamos botando os multiplicadores junto, mas marcando que eles não fazem parte da matriz em si! Aqui eu estou pondo em chaves, mas vocês podem marcar isso de outra forma! Só não pode confundir haha.

L 3 ← L 3 - 1 2 L 1 m 31 = 1 2 :

2 1 - 1 { 0.5 } - 1.5 1.5 { 0.5 } 0.5 1.5

Já está pivoteado, então continuamos a eliminação.

2 1 - 1 { 0.5 } - 1.5 1.5 { 0.5 } 0.5 1.5

L 3 ← L 3 - - 1 3 L 2 m 32 = - 1 3 :

2 1 - 1 { 0.5 } - 1.5 1.5 { 0.5 } { - 0.33 } 2

Então ficamos com:

P = 0 0 1 0 1 0 1 0 0

Que foi formada da matriz identidade 1 0 0 0 1 0 0 0 1 com as linhas 1 e 3 trocadas.

L = 1 0 0 0.5 1 0 0.5 - 1 3 1

U = 2 1 - 1 0 - 1.5 1.5 0 0 2

Prontinho! Se tivéssemos que resolver o sistema era só botarmos na forma L U x = P - 1 b e resolver como de costume. Repare que quando fazemos apenas uma troca de linhas, temos P - 1 = P .

Algoritmo

Novamente podemos dividir esse algoritmo em três fases bem distintas. Fase 1 será a eliminação de Gauss. Só que agora com o pivoteamento! E além das matrizes L e U essa fase também teria como saída a matriz P . A Fase 2 receberia a matriz L vinda da Fase 1 e o vetor P b para resolver o sistema L y = P b , que seria resolvido não por retrosubstituição, mas por substituição simples. A Fase 3 receberia o resultado da fase 2 , o vetor y , e a matriz U da Fase 1, e aplicaríamos nosso algoritmo da retrosubstituição. Vamos lá então:

Algoritmo da Fase 1:

Passo 1. Dados iniciais (o que seu programa vai receber de você):

a. Matriz quadrada A n × n

Passo 2. Definimos P como a matriz identidade

Passo 3. Percorremos as linhas de A com uma variável i (de 1 a n ).

Passo 4. l = i + 1

Passo 5. p o s M a i o r = i

Passo 6. m a i o r = a i i

Passo 7. Se a l i > m a i o r , então m a i o r = a l i e p o s M a i o r = l .

Passo 8. l = l + 1

Passo 9. Enquanto l ≤ n , volte para o Passo 6. Caso contrário, continue.

Passo 10. Troque as linhas i e p o s M a i o r de A e de P .

Passo 11. Escolher pivô da i -ésima linha de A   e armazenar na variável p .

Passo 12. Para cada ‘ i ’, percorremos as linhas de A de i + 1 até n com uma variável k .

Passo 13. Definimos o multiplicador m i k = a k i p .

Passo 14. Guardar m i k na posição linha i e coluna k de A .

Passo 15. Para cada ‘ k ’, percorremos as colunas de A i + 1 até n com uma variável j .

Passo 16. a k j = a k j - m i k ∙ a i j

Passo 17. Depois que acabarem todos os ' i ' do Passo 4, é importante separarmos as matrizes L e U da matriz A . A matriz L será uma matriz triangular inferior com 1 na diagonal e elementos a i j com i > j iguais ao da matriz A . Já a matriz U será uma matriz triangular superior com os elementos a i j com i ≤ j iguais ao da matriz A .

Esse último passo é necessário por causa daquela história de armazenar os multiplicadores na própria matriz que estamos escalonando. Então temos que separar no final. Repare também que como o objetivo aqui é resolver sistemas no algoritmo P é o nosso P - 1 da teoria já que fazemos as mudanças na ordem que elas acontecem.

Pronto, tá feito o nosso algoritmo. Daí vamos tirar nossas matrizes P , U e L . Agora o algoritmo da Fase 2 segue igual ao da fatoração L U , mas recebendo o vetor P b ao invés de b , então não vale a pena escrevê-lo aqui de novo.

Para a Fase 3 também temos exatamente a mesma coisa da fatoração L U .

Algoritmo

Exercícios Resolvidos

Exercício Resolvido #1

Considere as seguintes matrizes. Encontre a matriz de permutação P tal que P A pode ser fatorada no produto L U , onde L é uma matriz triangular inferior com 1 s em sua diagonal e U é triangular superior.

a A = 1 2 - 1 2 4 0 0 1 - 1

b A = 0 1 1 1 - 2 - 1 1 - 1 1

c A = 1 1 - 1 0 1 1 4 3 2 - 1 2 4 2 - 1 2 3

d A = 0 1 1 2 0 1 1 - 1 1 2 - 1 3 1 1 2 0

Passo 1

A primeira coisa que temos que perceber é que não devemos usar a fatoração P L U em si! Ele apenas quer que façamos o pivoteamento quando for necessário para a fatoração L U ! Vamos lá.

Passo 2

a Realizando a Eliminação de Gauss:

1 2 - 1 2 4 0 0 1 - 1

L 2 ← L 2 - 2 L 1 :

1 2 - 1 0 0 2 0 1 - 1

Repare que agora, para continuarmos com o processo precisamos fazer a troca das linhas.

L 2 ↔ L 3 :

1 2 - 1 0 1 - 1 0 0 2

Pronto. Repare que aí já terminamos o processo da fatoração. Mas ele não está interessado nas matrizes L e U , apenas na P que é dada pela matriz identidade com as linhas 2 e 3 trocadas.

P = 1 0 0 0 1 0 0 0 1 → L 1 ↔ L 2 1 0 0 0 0 1 0 1 0

Passo 3

b Realizando a Eliminação de Gauss:

0 1 1 1 - 2 - 1 1 - 1 1

L 1 ↔ L 2 :

1 - 2 - 1 0 1 1 1 - 1 1

L 3 ← L 3 - L 1 :

1 - 2 - 1 0 1 1 0 1 2

L 3 ← L 3 - L 2 :

1 - 2 - 1 0 1 1 0 0 1

Trocamos apenas as linhas 1 e 2 , então:

P = 1 0 0 0 1 0 0 0 1 → L 2 ↔ L 3 0 1 0 1 0 0 0 0 1

Passo 4

c Realizando a Eliminação de Gauss:

1 1 - 1 0 1 1 4 3 2 - 1 2 4 2 - 1 2 3

L 2 ← L 2 - L 1 :

1 1 - 1 0 0 0 5 3 2 - 1 2 4 2 - 1 2 3

L 3 ← L 3 - 2 L 1 :

1 1 - 1 0 0 0 5 3 0 - 3 4 4 2 - 1 2 3

L 4 ← L 4 - 2 L 1 :

1 1 - 1 0 0 0 5 3 0 - 3 4 4 0 - 3 4 3

L 2 ↔ L 3 :

1 1 - 1 0 0 - 3 4 4 0 0 5 3 0 - 3 4 3

L 4 ← L 4 - L 2 :

1 1 - 1 0 0 - 3 4 4 0 0 5 3 0 0 0 - 1

Trocamos apenas as linhas 2 e 3 , então:

P = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 → L 2 ↔ L 3 1 0 0 0 0 0 1 0 0 1 0 0 0 0 0 1

Passo 5

d Realizando a Eliminação de Gauss:

0 1 1 2 0 1 1 - 1 1 2 - 1 3 1 1 2 0

L 1 ↔ L 3 :

1 2 - 1 3 0 1 1 - 1 0 1 1 2 1 1 2 0

L 4 ← L 4 - L 1 :

1 2 - 1 3 0 1 1 - 1 0 1 1 2 0 - 1 3 - 3

L 3 ← L 3 - L 2 :

1 2 - 1 3 0 1 1 - 1 0 0 0 3 0 - 1 3 - 3

L 4 ← L 4 + L 2 :

1 2 - 1 3 0 1 1 - 1 0 0 0 3 0 0 4 - 4

L 3 ↔ L 4 :

1 2 - 1 3 0 1 1 - 1 0 0 4 - 4 0 0 0 3

Trocamos as linhas 1 e 3 , e em seguida as linhas 3 e 4 , então:

P = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 → L 3 ↔ L 4 1 0 0 0 0 1 0 0 0 0 0 1 0 0 1 0 → L 1 ↔ L 3 0 0 0 1 0 1 0 0 1 0 0 0 0 0 1 0

Resposta

a   1 0 0 0 0 1 0 1 0

b 0 1 0 1 0 0 0 0 1

c 1 0 0 0 0 0 1 0 0 1 0 0 0 0 0 1

d 0 0 0 1 0 1 0 0 1 0 0 0 0 0 1 0

Exercício Resolvido #2

Obtenha a fatoração na forma A = P L U para as seguintes matrizes.

a A = 0 2 3 1 1 - 1 0 - 1 1

b A = 1 2 - 1 1 2 3 2 - 1 4

c A = 1 - 2 3 0 3 - 6 9 3 2 1 4 1 1 - 2 2 - 2

d A = 1 - 2 3 0 1 - 2 3 1 1 - 2 2 - 2 2 1 3 1

Passo 1

a Realizando a Eliminação de Gauss com pivoteamento parcial:

0 2 3 1 1 - 1 0 - 1 1

L 1 ↔ L 2 :

1 1 - 1 0 2 3 0 - 1 1

L 3 ← L 3 - - 1 2 L 1 m 32 = - 1 2 :

1 1 - 1 0 2 3 0 - 1 2 2.5

Então ficamos com:

P = 0 1 0 1 0 0 0 0 1 ; L = 1 0 0 0 1 0 0 - 1 2 1 ; U = 1 1 - 1 0 2 3 0 0 2.5

Passo 2

b Realizando a Eliminação de Gauss com pivoteamento parcial:

1 2 - 1 1 2 3 2 - 1 4

L 1 ↔ L 3 :

2 - 1 4 1 2 3 1 2 - 1

L 2 ← L 2 - 1 2 L 1 m 21 = 1 2 :

2 - 1 4 1 2 2.5 1 1 2 - 1

L 3 ← L 3 - 1 2 L 1 m 31 = 1 2 :

2 - 1 4 1 2 2.5 1 1 2 2.5 - 3

L 3 ← L 3 - L 2 m 32 = 1 :

2 - 1 4 1 2 2.5 1 1 2 1 - 4

Então ficamos com:

P = 0 0 1 0 1 0 1 0 0 ; L = 1 0 0 0.5 1 0 0.5 1 1 ; U = 2 - 1 4 0 2.5 1 0 0 - 4

Passo 3

c Realizando a Eliminação de Gauss:

1 - 2 3 0 3 - 6 9 3 2 1 4 1 1 - 2 2 - 2

L 1 ↔ L 2 :

3 - 6 9 3 1 - 2 3 0 2 1 4 1 1 - 2 2 - 2

L 2 ← L 2 - 1 3 L 1 m 21 = 1 3 :

3 - 6 9 3 1 3 0 0 - 1 2 1 4 1 1 - 2 2 - 2

L 3 ← L 3 - 2 3 L 1 m 31 = 2 3 :

3 - 6 9 3 1 3 0 0 - 1 2 3 5 - 2 - 1 1 - 2 2 - 2

L 4 ← L 4 - 1 3 L 1 m 41 = 1 3 :

3 - 6 9 3 1 3 0 0 - 1 2 3 5 - 2 - 1 1 3 0 - 1 - 3

L 2 ↔ L 3 :

3 - 6 9 3 2 3 5 - 2 - 1 1 3 0 0 - 1 1 3 0 - 1 - 3

Reparou o que aconteceu?? Os multiplicadores trocaram de lugar junto com as linhas! É nesses casos que armazenar os multiplicadores junto da matriz que estamos escalonando é útil pois evita problemas no final!

L 3 ↔ L 4 :

3 - 6 9 3 2 3 5 - 2 - 1 1 3 0 - 1 - 3 1 3 0 0 - 1

Eita, várias trocas de linha né? Vamos ver com calma como ficou nossa matriz P . Lembrando que temos que fazer as trocas de trás pra frente:

P = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 → L 3 ↔ L 4 1 0 0 0 0 1 0 0 0 0 0 1 0 0 1 0 → L 2 ↔ L 3 1 0 0 0 0 0 0 1 0 1 0 0 0 0 1 0 → L 1 ↔ L 2 0 0 0 1 1 0 0 0 0 1 0 0 0 0 1 0

Então ficamos com:

P = 0 0 0 1 1 0 0 0 0 1 0 0 0 0 1 0 ; L = 1 0 0 0 2 3 1 0 0 1 3 0 1 0 1 3 0 0 1 ; U = 3 - 6 9 3 0 5 - 2 - 1 0 0 - 1 - 3 0 0 0 - 1

Passo 4

d Realizando a Eliminação de Gauss:

1 - 2 3 0 1 - 2 3 1 1 - 2 2 - 2 2 1 3 1

L 1 ↔ L 4 :

2 1 3 1 1 - 2 3 1 1 - 2 2 - 2 1 - 2 3 0

L 2 ← L 2 - 1 2 L 1   m 21 = 1 2 :

2 1 3 1 0.5 - 2.5 1.5 0.5 1 - 2 2 - 2 1 - 2 3 0

L 3 ← L 3 - 1 2 L 2 m 31 = 1 2 :

2 1 3 1 0.5 - 2.5 1.5 0.5 0.5 - 2 . 5 0.5 - 2.5 1 - 2 3 0

L 4 ← L 4 - 1 2 L 2 m 41 = 1 2 :

2 1 3 1 0.5 - 2.5 1.5 0.5 0.5 - 2.5 0.5 - 2.5 0.5 - 2.5 1.5 - 0.5

L 3 ← L 3 - L 2 m 32 = 1 :

2 1 3 1 0.5 - 2.5 1.5 0.5 0.5 1 - 1 - 3 0.5 - 2.5 1.5 - 0.5

L 4 ← L 4 - L 2 m 42 = 1 :

2 1 3 1 0.5 - 2.5 1.5 0.5 0.5 1 - 1 - 3 0.5 1 0 - 1

Então ficamos com:

P = 0 0 0 1 0 1 0 0 0 0 1 0 1 0 0 0 ; L = 1 0 0 0 0.5 1 0 0 0.5 1 1 0 0.5 1 0 1 ; U = 2 1 3 1 0 - 2.5 1.5 0.5 0 0 - 1 - 3 0 0 0 - 1

Resposta

a   A = 0 1 0 1 0 0 0 0 1 1 0 0 0 1 0 0 - 1 2 1 1 1 - 1 0 2 3 0 0 2.5

b A = 0 0 1 0 1 0 1 0 0 1 0 0 0.5 1 0 0.5 1 1 2 - 1 4 0 2.5 1 0 0 - 4

c A = 0 0 0 1 1 0 0 0 0 1 0 0 0 0 1 0 1 0 0 0 2 3 1 0 0 1 3 0 1 0 1 3 0 0 1 3 - 6 9 3 0 5 - 2 - 1 0 0 - 1 - 3 0 0 0 - 1

d A = 0 0 0 1 0 1 0 0 0 0 1 0 1 0 0 0 1 0 0 0 0.5 1 0 0 0.5 1 1 0 0.5 1 0 1 2 1 3 1 0 - 2.5 1.5 0.5 0 0 - 1 - 3 0 0 0 - 1

Exercício Resolvido #3

Justifique se for verdadeira ou dê um contraexemplo se for falsa a afirmação:

“Dada uma matriz A , n × n , sua fatoração L U , obtida com estratégia de pivoteamento parcial, é tal que todos os elementos da matriz L têm módulo menor ou igual a 1 .”

Passo 1

-

Essa afirmativa é verdadeira! Mas por quê? Bom, quais são os elementos que compõe a matriz L ? De forma geral ela é dada por:

1 0 ⋯ 0 0 m 21 1 ⋯ 0 0 ⋮ ⋮ ⋱ ⋮ ⋮ m n - 1 1 m n - 1 2 ⋯ 1 0 m n 1 m n 2 ⋯ m n n - 1 1

Ou seja, ela é uma matriz triangular inferior cuja diagonal é toda igual a 1 e contém os multiplicadores usados na eliminação de Gauss.

Bom, a única coisa que pode ter módulo maior ou igual a 1 são os multiplicadores certo? E como eles são formados? Assim:

m i j = a i j p ; i > j

onde p é o pivô da linha j . Como estamos usando a estratégia de pivoteamento parcial, o pivô de uma linha, se comparado com todos os números abaixo dele na coluna, será sempre maior em módulo, e isso é garantido pelo pivoteamento.

Sendo assim, a divisão a i j p terá sempre seu módulo menor ou igual a 1 . Logo, a afirmativa é verdadeira.

Resposta

Ei, a resposta está no passo a passo :)

Exercício Resolvido #4

Resolva os sistemas lineares abaixo através da fatoração L U com estratégia de pivoteamento parcial. Discuta a existência ou não de soluções.

a   1 x 1 + 1 x 2 + 1 x 3 = 4 2 x 1 + 1 x 2 - 1 x 3 = - 1 3 x 1 + 2 x 2 + 0 x 3 = 3

b 1 x 1 + 1 x 2 + 1 x 3 = 4 2 x 1 + 1 x 2 - 1 x 3 = - 1 3 x 1 + 2 x 2 + 0 x 3 = 6  

Passo 1

Se liga numa coisa, a matriz A do sistema matricial associado A x = b é a mesma pros dois casos né? Então vamos achar a decomposição P L U dela:

A = 1 1 1 2 1 - 1 3 2 0

Aplicando a eliminação de Gauss com pivoteamento parcial:

L 1 ↔ L 3 :

3 2 0 2 1 - 1 1 1 1

L 2 ← L 2 - 2 3 L 1 m 21 = 2 3 :

3 2 0 2 3 - 1 3 - 1 1 1 1

L 3 ← L 3 - 1 3 L 1 m 31 = 1 3 :

3 2 0 2 3 - 1 3 - 1 1 3 1 3 1

L 3 ← L 3 - - 1 L 2 m 32 = - 1 :

3 2 0 2 3 - 1 3 - 1 1 3 - 1 0

Então ficamos com:

P = 0 0 1 0 1 0 1 0 0 ; L = 1 0 0 2 3 1 0 1 3 - 1 1 ; U = 3 2 0 0 - 1 3   - 1 0 0 0

Passo 2

a Beleza, agora podemos resolver os sistemas!

Resolvendo L y = P - 1 b :

1 0 0 2 3 1 0 1 3 - 1 1 y 1 y 2 y 3 = 0 0 1 0 1 0 1 0 0 4 - 1 3

Lembrando que como só houve uma troca de linhas, P = P - 1 !

1 0 0 2 3 1 0 1 3 - 1 1 y 1 y 2 y 3 = 3 - 1 4

y 1 = 3

2 + y 2 = - 1

y 2 = - 3

1 + 3 + y 3 = 4

y 3 = 0

Resolvendo U x = y :

3 2 0 0 - 1 3   - 1 0 0 0 x 1 x 2 x 3 = 3 - 3 0

3 x 1 + 2 x 2 = 3

x 1 = 1 - 2 x 2 3

- x 2 3 - x 3 = - 3

x 3 = 3 - x 2 3

Então a solução é x 1 , x 2 , x 3 = 1 - 2 x 2 3 , x 2 , 3 - x 2 3

Como podemos ver, esse sistema tem infinitas soluções ;)

Passo 3

Como a matriz A é a mesma não tem necessidade de fatorarmos ela novamente. É só alterar o vetor b :

Resolvendo L y = P - 1 b :

1 0 0 2 3 1 0 1 3 - 1 1 y 1 y 2 y 3 = 0 0 1 0 1 0 1 0 0 4 - 1 6

Lembrando que como só houve uma troca de linhas, P = P - 1 !

1 0 0 2 3 1 0 1 3 - 1 1 y 1 y 2 y 3 = 6 - 1 4

y 1 = 6

4 + y 2 = - 1

y 2 = - 5

1 + 5 + y 3 = 4

y 3 = - 2

Resolvendo U x = y :

3 2 0 0 - 1 3   - 1 0 0 0 x 1 x 2 x 3 = 6 - 5 - 2

Opa, olha a última equação! Temos 0 = - 2 ! Ou seja, esse sistema não tem solução.

Resposta

a x 1 , x 2 , x 3 = 1 - 2 x 2 3 , x 2 , 3 - x 2 3

b Não possui solução.

Exercício Resolvido #5

Resolva o sistema linear abaixo através da fatoração L U com estratégia de pivoteamento parcial:

1 x 1 - 2 x 2 + 7 x 3 + 2 x 4 = - 18 2 x 1 + 5 x 2 - 3 x 3 + 1 x 4 = 31 9 x 1 - 6 x 2 + 4 x 3 + 1 x 4 = 35 4 x 1 - 3 x 2 - 6 x 3 + 7 x 4 = 15

Passo 1

Primeiro precisamos saber quem é A :

A = 1 - 2 7 2 2 5 - 3 1 9 - 6 4 1 4 - 3 - 6 7

Aplicando a eliminação de Gauss com pivoteamento parcial:

L 1 ↔ L 3 :

9 - 6 4 1 2 5 - 3 1 1 - 2 7 2 4 - 3 - 6 7

L 2 ← L 2 - 0.222 L 1 m 21 = 0.222 :

9 - 6 4 1 0.222 6.333 - 3.889 0.778 1 - 2 7 2 4 - 3 - 6 7

L 3 ← L 3 - 0.111 L 1 m 31 = 0.111 :

9 - 6 4 1 0.222 6.333 - 3.889 0.778 0.111 - 1.333 6.556 1.889 4 - 3 - 6 7

L 4 ← L 4 - 0.444 L 1 m 41 = 0.444 :

9 - 6 4 1 0.222 6.333 - 3.889 0.778 0.111 - 1.333 6.556 1.889 0.444 - 0.333 - 7.778 6.556

L 3 ← L 3 - - 0.211 L 2 m 32 = - 0.211 :

9 - 6 4 1 0.222 6.333 - 3.889 0.778 0.111 - 0.211 5.737 2.053 0.444 - 0.333 - 7.778 6.556

L 4 ← L 4 - - 0.053 L 2 m 42 = - 0.053 :

9 - 6 4 1 0.222 6.333 - 3.889 0.778 0.111 - 0.211 5.737 2.053 0.444 - 0.053 - 7.982 6.596

L 3 ↔ L 4 :

9 - 6 4 1 0.222 6.333 - 3.889 0.778 0 . 444 - 0 . 053 - 7.982 6.596 0 . 111 - 0.211 5.737 2.053

L 4 ← L 4 - - 0.719 L 3 m 43 = - 0.719 :

9 - 6 4 1 0.222 6.333 - 3.889 0.778 0.444 - 0 . 053 - 7.982 6.596 0 . 111 - 0.211 - 0.719 6.793

A P - 1 vai ficar:

P - 1 = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 → L 1 ↔ L 3 0 0 1 0 0 1 0 0 1 0 0 0 0 0 0 1 → L 3 ↔ L 4 0 0 1 0 0 1 0 0 0 0 0 1 1 0 0 0

Então ficamos com:

P - 1 = 0 0 1 0 0 1 0 0 0 0 0 1 1 0 0 0 ; L = 1 0 0 0 0.222 1 0 0 0.444 - 0.053 1 0 0.111 - 0.211 - 0.719 1 ;  

U = 9 - 6 4 1 0 6.333 - 3.889 0.778 0 0 - 7.982 6.596 0 0 0 6.793

Passo 2

Agora é só resolver os sistemas!

Resolvendo L y = P - 1 b :

1 0 0 0 0.222 1 0 0 0.444 - 0.053 1 0 0.111 - 0.211 - 0.719 1 y 1 y 2 y 3 y 4 = 0 0 1 0 0 1 0 0 0 0 0 1 1 0 0 0 - 18 31 35 15

1 0 0 0 0.222 1 0 0 0.444 - 0.053 1 0 0.111 - 0.211 - 0.719 1 y 1 y 2 y 3 y 4 = 35 31 15 - 18

y 1 = 35

7.778 + y 2 = 31

y 2 = 23.222

15.556 - 1.231 + y 3 = 15

y 3 = 0.675

3.889 - 4.9 - 0.485 + y 4 = - 18

y 4 = - 16.504

Resolvendo U x = y :

9 - 6 4 1 0 6.333 - 3.889 0.778 0 0 - 7.982 6.596 0 0 0 6.793 x 1 x 2 x 3 x 4 = 35 23.222 0.675 - 16.504

x 4 = - 2.429

- 7.982 x 3 - 16.022 = 0.675

x 3 = - 2.092

6.333 x 2 + 8.136 - 1.9 = 23.222

x 2 = 2.68

9 x 1 - 16.08 - 8.368 - 2.429 = 35

x 1 = 6.876

Então a solução é x 1 , x 2 , x 3 , x 4 = 6.876,2 . 68 , - 2.092 , - 2.429

Resposta

x 1 , x 2 , x 3 , x 4 = 6.876,2.68 , - 2.092 , - 2.429