Lição
A teoria — Eliminação de Gauss269 palavras
A eliminação transforma um sistema noutro equivalente triangular, que pode depois ser solucionado de baixo para cima. O seu valor não está em funcionar. Vários métodos funcionam. Está naquilo que custa, porque é esse custo que determina o método que a máquina vai realmente executar.
O que significa cada símbolo
n- o número de incógnitas, e a única coisa de que o custo depende.
pivot- o coeficiente pelo qual uma linha é dividida. A ordem fixada pelo algoritmo é o que torna o processo livre de decisões. O número que calha nesse lugar é o que assegura a sua sobrevivência em vírgula flutuante.
- Pressupõe
- Um sistema denso e sem qualquer estrutura explorável. Sistemas esparsos, em banda ou simétricos ficam mais baratos. Saber escolher o método especializado adequado resume grande parte da álgebra linear numérica.
- Falha quando
- O custo é cúbico. Cúbico é uma barreira. A eliminação sobre n incógnitas requer cerca de
n³/3multiplicações e igual número de adições. Contando diretamente numa eliminação real, temos 333.300 multiplicações comn = 100face às 333.333 da estimativa. Aquilo que o expoente nos cobra é uma penalização fixa de crescimento: duplicar as incógnitas multiplica o trabalho por exatamente 8. Desse modo, 100 incógnitas requerem 333.300 multiplicações e 200 requerem 2.666.600. Uma máquina que despache mil incógnitas num segundo precisa de oito segundos para duas mil. Levará cerca de duas horas e um quarto para vinte mil. Esse é o formato da barreira. É a razão pela qual a pergunta interessante da álgebra linear numérica quase nunca passa por como resolver um sistema denso, mas sim como evitar sequer tê-lo.
Problemas resolvidos na íntegra
-
Resolver 2x + y − z = 8 e o custo em precisão 7 passos
Resolva 2x + y − z = 8, −3x − y + 2z = −11, −2x + y + 2z = −3 — e depois avalie o que a eliminação lhe custou em exatidão.
-
Escreva como uma matriz aumentada. As letras não fazem nenhum trabalho, por isso descarte-as e mantenha as colunas.
-
Faça o pivotamento pela maior entrada da primeira coluna, que é −3 na linha dois, e desloque-a para cima. Cada troca de linha inverte o sinal do determinante, portanto acompanhe a conta.
-
Zere o restante da primeira coluna subtraindo múltiplos da linha do pivô.
-
Repita o procedimento na segunda coluna. A seguir, leia a última linha: ela contém uma incógnita e o seu valor.
-
Faça a substituição regressiva de baixo para cima. A última linha fornece z, a linha do meio dá y e a do topo revela x.
-
Multiplique os pivôs e aplique o sinal correspondente a duas trocas para obter o determinante. Verifique o resultado substituindo-o nas equações originais em vez das reduzidas.
-
O resíduo aproxima-se de 10⁻¹⁶, o que não é zero. Trata-se do tamanho do arredondamento que um computador não consegue evitar, e é a medida honesta a relatar. Uma resposta anunciada como exata ou provém de aritmética inteira ou não está a prestar atenção. Altere o coeficiente principal para 10⁻¹⁷ com o pivotamento desligado e o mesmo resíduo passará a 1, o que equivale ao mesmo aviso a soar aos gritos.
Resposta
x = 2, y = 3, z = −1, determinante −1. O resíduo perto de 10⁻¹⁶ significa que a aritmética foi tão exata quanto o ponto flutuante permite. Desligue o pivotamento num sistema mal escalonado e esse será o número a avisar que a resposta está errada.
-
-
Um determinante de 0, duas vezes, com dois significados diferentes 7 passos
Esta predefinição não tem solução e a seguinte tem uma infinidade delas, e ambas apresentam um determinante de 0. Determine o que realmente as separa e decida o que mudar o segundo membro alguma vez lhe poderia comprar.
-
A pivotagem parcial percorre a primeira coluna à procura do maior coeficiente. Encontra o 3 na última linha, pelo que o primeiro movimento é uma troca.
-
Anule a coluna: subtraia ⅔ da nova primeira linha à segunda e ⅓ dela à terceira.
-
Na segunda coluna, 4/3 já está acima de 2/3, portanto nada precisa de ser trocado. Subtraia metade da linha 2 à linha 3.
-
A linha 3 fica agora vazia à esquerda e com −½ à direita, o que se lê 0 = −½. Foram colocados apenas dois pivôs, 3 e 4/3, onde três incógnitas exigem três, e um pivô em falta anula o produto.
-
O que o determinante esconde é a característica. A matriz dos coeficientes tem duas linhas independentes. A matriz ampliada tem três, porque aquele −½ que sobra é uma linha que nada consegue anular. Característica 2 contra característica 3 é o que significa «sem solução».
-
A predefinição Equação repetida faz a mesma aritmética e diz outra coisa. Cada linha é um múltiplo da primeira, pelo que ambas as características valem 1 e o determinante volta a ser 0. O conjunto das soluções tem 3 − 1 = 2 direções livres: o plano inteiro x + y + z = 3. A ferramenta diz «uma infinidade» sem dizer em quantas direções.
-
Mude um número e veja qual das grandezas se mexe. Coloque 6 onde está o 7, à direita da linha 2. A matriz dos coeficientes fica intacta, pelo que o determinante nem estremece, mas a característica da ampliada desce para 2 e a resposta passa de nenhuma solução para uma reta de soluções.
Resposta
O determinante limita-se a dizer se a matriz dos coeficientes tem o conjunto completo de pivôs. Nunca olha para o segundo membro, e é por isso que não consegue separar as duas predefinições, ao passo que as duas características conseguem: sem solução quando discordam, e um conjunto de soluções de dimensão 3 − característica quando coincidem. Assim que o determinante é 0, nenhum segundo membro lhe comprará uma resposta única: obtém nenhuma ou toda uma família, e a escolha entre essas duas é tudo o que b decide. O cartão do resíduo apaga-se pela mesma razão por que existe: um resíduo precisa de uma solução para repor nas equações.
-
Percurso de aprendizagem
Isolar o x, a partir do sinal de menos
Referências (1)
- Two centuries of the method being reinvented and misattributed, including how Gauss's name came to be on it: J. F. Grcar, "How ordinary elimination became Gaussian elimination." Historia Mathematica 38:2 (2011), 163–218.