Resolva o sistema linear abaixo através da fatoração LU com estratégia de pivoteamento parcial:
\left\{\begin{array}{l} x_{1}-2 x_{2}+7 x_{3}+2 x_{4}=-18 \\\ 2 x_{1}+5 x_{2}-3 x_{3}+x_{4}=31 \\\ 9 x_{1}-6 x_{2}+4 x_{3}+x_{4}=35 \\\ 4 x_{1}-3 x_{2}-6 x_{3}+7 x_{4}=15 \end{array}\right.

Passo 1

O enunciado nos dá o sistema linear
\begin{cases} x_{1}-2 x_{2}+7 x_{3}+2 x_{4}=-18\\[4pt] 2 x_{1}+5 x_{2}-3 x_{3}+x_{4}=31\\[4pt] 9 x_{1}-6 x_{2}+4 x_{3}+x_{4}=35\\[4pt] 4 x_{1}-3 x_{2}-6 x_{3}+7 x_{4}=15 \end{cases}
As informações relevantes que temos são:
A=\begin{pmatrix} 1 & -2 & 7 & 2\\[4pt] 2 & 5 & -3 & 1\\[4pt] 9 & -6 & 4 & 1\\[4pt] 4 & -3 & -6 & 7 \end{pmatrix},\qquad b=\begin{pmatrix}-18\\[3pt]31\\[3pt]35\\[3pt]15\end{pmatrix}.
O que se pede: encontrar o vetor solução (x_1,x_2,x_3,x_4)^T usando fatoração LU com pivoteamento parcial (ou seja, permitir permutações de linhas quando necessário).]

Passo 2

Breve lembrete didático (passo-a-passo do que vamos fazer, sem pular etapas):
  1. Pivoteamento parcial significa, em cada coluna k durante a eliminação, trocar a linha k pela linha com maior valor absoluto na coluna k (entre as linhas k…n). Isso evita divisões por números muito pequenos e produz uma matriz permutada P tal que PA = LU.
  2. A fatoração PA=LU dá-nos uma matriz L triangular inferior com 1 na diagonal e uma matriz U triangular superior. Depois resolvemos Ax=b assim:
    • montar PA=LU,
    • calcular Pb,
    • resolver L y = Pb (substituição progressiva),
    • resolver U x = y (substituição regressiva).
  3. Em cada etapa eu mostro os pivôs, os multiplicadores (entradas de L), e atualizo as linhas de U e o vetor direitos b. Vou trabalhar com frações exatas para não perder precisão.

Passo 3

Primeiro passo do pivoteamento (coluna 1). Olhamos a coluna 1 (linhas 1…4) e escolhemos o maior valor absoluto:
|1|,|2|,|9|,|4| → máximo em |9|, que está na linha 3. Então trocamos a linha 1 com a linha 3. Isso corresponde a uma permutação parcial inicial. Aplicando essa troca a A e b obtemos (mostro PA parcial):
A^{(1)} = \begin{pmatrix} 9 & -6 & 4 & 1\\[4pt] 2 & 5 & -3 & 1\\[4pt] 1 & -2 & 7 & 2\\[4pt] 4 & -3 & -6 & 7 \end{pmatrix},\qquad b^{(1)}=\begin{pmatrix}35\\[3pt]31\\[3pt]-18\\[3pt]15\end{pmatrix}.
Agora eliminamos abaixo do pivô a_{11}=9. Os multiplicadores (coluna 1 de L) são:
\ell_{21}=\frac{2}{9},\quad \ell_{31}=\frac{1}{9},\quad \ell_{41}=\frac{4}{9}.
Aplicando as eliminações (linha i := linha i − \ell_{i1}·linha1) obtemos as linhas 2…4 atualizadas e também atualizamos b segundo o mesmo procedimento:
  • Linha 2 → [0,\tfrac{19}{3},-\tfrac{35}{9},\tfrac{7}{9}], b_2=\tfrac{209}{9}.
  • Linha 3 → [0,-\tfrac{4}{3},\tfrac{59}{9}-\tfrac{4}{9}=\tfrac{59}{9},\tfrac{17}{9}] — mas vamos simplificar corretamente como fiz abaixo (mantendo frações reduzidas quando possível).
  • Linha 3 ficou (após cálculo exato): [0,-\tfrac{4}{3},\tfrac{59}{9}-\tfrac{4}{9}=\tfrac{59-4}{9}=\tfrac{55}{9}] — atenção: essa forma intermédia foi usada, mas a continuação mais consistente com os próximos passos resulta nas frações que seguem nos passos seguintes.
  • Linha 4 → [0,-\tfrac{1}{3},-\tfrac{70}{9},\tfrac{59}{9}], b_4=-\tfrac{5}{9}.
    (Os cálculos acima correspondem à subtração linha por linha usando os multiplicadores. Mantive frações exatas.)]

Passo 4

Segundo pivô (coluna 2, restrito às linhas 2…4). Olhamos os valores absolutos na coluna 2 das linhas 2…4:
\left|\tfrac{19}{3}\right|,\left|\!-\tfrac{4}{3}\right|,\left|\!-\tfrac{1}{3}\right|. O maior é \tfrac{19}{3} na linha 2, então NÃO precisamos trocar linhas nesta etapa. O pivô é
a_{22}=\frac{19}{3}.
Os multiplicadores para eliminar abaixo deste pivô (na coluna 2) são:
\ell_{32}=\frac{-\tfrac{4}{3}}{\tfrac{19}{3}}=-\frac{4}{19},\qquad \ell_{42}=\frac{-\tfrac{1}{3}}{\tfrac{19}{3}}=-\frac{1}{19}.
Aplicando as eliminações (linha i := linha i − \ell_{i2}·linha2) obtemos as linhas 3 e 4 atualizadas e atualizamos b:
  • Linha 2 permanece [0,\tfrac{19}{3},-\tfrac{35}{9},\tfrac{7}{9}] e b_2=\tfrac{209}{9}.
  • Linha 3 fica (após eliminar): [0,0,\tfrac{109}{19},\tfrac{39}{19}] e b_3=-17.
  • Linha 4 fica (após eliminar): [0,0,-\tfrac{455}{57},\tfrac{376}{57}] e b_4=\tfrac{2}{3}.
(Obs.: eu calculei b_3 exatamente como -17 e b_4=\tfrac{2}{3} ao atualizar o vetor direito com os mesmos multiplicadores — isso é importante para a consistência.)

Passo 5

Terceiro pivô (coluna 3, linhas 3…4). Agora examinamos os elementos na coluna 3, linhas 3 e 4:
\left|\tfrac{109}{19}\right|\approx5{,}74,\qquad \left|\!-\tfrac{455}{57}\right|\approx7{,}98.
O maior em módulo é -\tfrac{455}{57} (linha 4), portanto fazemos uma troca de linha entre a linha 3 e a linha 4 (pivoteamento parcial). Isso altera a permutação global P (uma segunda troca). Após essa troca, as linhas 3 e 4 trocam de lugar. Ao trocar linhas durante a eliminação, aplicamos a mesma troca às colunas já preenchidas de L (ou seja, trocamos os multiplicadores já calculados para manter PA=LU consistente). Após a troca, o pivô da etapa é
a_{33}=-\frac{455}{57}.
O multiplicador para eliminar abaixo (somente linha 4 agora) é
\ell_{43}=\frac{\tfrac{109}{19}}{-\tfrac{455}{57}}=-\frac{327}{455}.
Aplicando a eliminação na linha 4 (linha4 := linha4 − \ell_{43}·linha3) obtemos a última diagonal de U e atualizamos b:
  • Linha 3 (após a troca) é [0,0,-\tfrac{455}{57},\tfrac{376}{57}] e b_3=\tfrac{2}{3}.
  • Linha 4 (após eliminação) passa a ter zeros nas colunas 1…3 e na coluna 4 vale \tfrac{3091}{455}. O b_4 atualizado vale -\tfrac{7517}{455}.

Passo 6

Concluímos a eliminação superior; agora podemos escrever as matrizes P, L e U que satisfazem PA=LU. A permutação final aplicada às linhas de A foi: trocar 1↔3 (primeira troca), depois trocar as atuais 3↔4 (segunda troca). Isso leva a ordem de linhas permutada [3,2,4,1]. Assim:
P=\begin{pmatrix} 0&0&1&0\\[4pt] 0&1&0&0\\[4pt] 0&0&0&1\\[4pt] 1&0&0&0 \end{pmatrix}.
A matriz U (triangular superior) obtida é
U=\begin{pmatrix} 9 & -6 & 4 & 1\\[4pt] 0 & \tfrac{19}{3} & -\tfrac{35}{9} & \tfrac{7}{9}\\[4pt] 0 & 0 & -\tfrac{455}{57} & \tfrac{376}{57}\\[4pt] 0 & 0 & 0 & \tfrac{3091}{455} \end{pmatrix}.
A matriz L (triangular inferior com 1 na diagonal) contém os multiplicadores usados, após termos aplicado as trocas apropriadas às linhas de L quando necessário:
L=\begin{pmatrix} 1 & 0 & 0 & 0\\[4pt] \tfrac{2}{9} & 1 & 0 & 0\\[4pt] \tfrac{4}{9} & -\tfrac{1}{19} & 1 & 0\\[4pt] \tfrac{1}{9} & -\tfrac{4}{19} & -\tfrac{327}{455} & 1 \end{pmatrix}.
Confirme mentalmente: com essa P,L,U tem-se PA=LU. (Se tivesse dúvida, pode multiplicar L por U e comparar com PA.)]

Passo 7

Agora resolvemos o sistema. Primeiro montamos Pb (aplicando a mesma permutação às entradas de b):
Pb=\begin{pmatrix}35\\[3pt]31\\[3pt]15\\[3pt]-18\end{pmatrix}.
Vamos resolver L y = Pb por substituição progressiva. Escrevo as equações de baixo nível:
  1. y_1 = 35.
  2. \tfrac{2}{9}y_1 + y_2 = 31 \Rightarrow y_2 = 31 - \tfrac{2}{9}\cdot 35 = \tfrac{209}{9}.
  3. \tfrac{4}{9}y_1 -\tfrac{1}{19}y_2 + y_3 = 15 \Rightarrow
    y_3 = 15 - \tfrac{4}{9}\cdot 35 + \tfrac{1}{19}\cdot \tfrac{209}{9} = \tfrac{2}{3}.
  4. \tfrac{1}{9}y_1 -\tfrac{4}{19}y_2 -\tfrac{327}{455}y_3 + y_4 = -18\Rightarrow
    y_4 = -18 - \tfrac{1}{9}\cdot 35 + \tfrac{4}{19}\cdot \tfrac{209}{9} + \tfrac{327}{455}\cdot \tfrac{2}{3} = -\tfrac{7517}{455}.
Portanto o vetor intermediário y é
y=\begin{pmatrix}35\\[3pt]\tfrac{209}{9}\\[3pt]\tfrac{2}{3}\\[3pt]-\tfrac{7517}{455}\end{pmatrix}.]

Passo 8

Agora resolvemos U x = y por substituição regressiva (de baixo para cima). As equações a resolver são:
  1. \dfrac{3091}{455}x_4 = -\dfrac{7517}{455}.
    Logo
x_4 = -\frac{7517}{3091}\approx -2{,}431780084.
  1. -\dfrac{455}{57}x_3 + \dfrac{376}{57}x_4 = \dfrac{2}{3}.
    Multiplicando por 57 e isolando x_3:
-455 x_3 + 376 x_4 = 38 \quad\Rightarrow\quad x_3 = \frac{376x_4-38}{455}.
Substituindo x_4=-\dfrac{7517}{3091} obtemos (cálculo exato conduz a uma fração); numericamente:
x_3 \approx -2{,}093450\\ (\text{aprox.}).
  1. \tfrac{19}{3}x_2 -\tfrac{35}{9}x_3 + \tfrac{7}{9}x_4 = \tfrac{209}{9}.
    Multiplicando por 9:
57 x_2 -35 x_3 + 7 x_4 = 209 \quad\Rightarrow\quad x_2 = \frac{209 + 35 x_3 -7 x_4}{57}.
Substituindo os valores numéricos de x_3,x_4 obtemos
x_2 \approx 2{,}678465\\ (\text{aprox.}).
  1. 9 x_1 -6 x_2 +4 x_3 + x_4 = 35.
    Isolando x_1:
x_1 = \frac{35 + 6 x_2 -4 x_3 - x_4}{9}.
Substituindo os valores numéricos obtemos
x_1 \approx 6{,}875152\\ (\text{aprox.}).
(Se preferir as frações exatas, seguem os valores exatos intermediários dos passos; as aproximações decimais acima têm precisão suficiente para conferência na prova.)

Resposta

\displaystyle x_{1}=6{,}875152
\displaystyle x_{2}=2{,}678465
\displaystyle x_{3}=-2{,}093450
\displaystyle x_{4}=-2{,}431780

Exercícios relacionados