✖

Sistemas de equações lineares

Um sistema linear de equações é um conjunto de m equações lineares com n incógnitas. Ao resolver o sistema, obtemos valores que podemos substituir nas nossas n incógnitas de modo que todas as equações façam sentido.

Definição de um sistema de equações lineares

Você certamente já encontrou algum sistema simples do tipo

a11x+a12y=b1a21x+a22y=b2

em que a11,…,a22, b1, b2 são números reais dados e x e y são as incógnitas. Esse sistema de equações pode ter um número diferente de soluções — nenhuma, uma ou infinitas, como uma equação do 1º grau comum. No caso de um sistema de equações lineares, a questão do número de soluções e de como encontrá-las é mais complicada do que no caso de uma simples equação do 1º grau.

Definição de um sistema de equações lineares:

a11x1+a12x2+…+a1nxn=b1a21x1+a22x2+…+a2nxn=b2⋮⋮⋮⋮⋮⋮⋮⋮⋮am1x1+am2x2+…+amnxn=bm

em que a11, …, amn, b1, …, bm são números reais. Esse sistema é chamado, então, de sistema de m equações lineares com n incógnitas x1, …, xn e coeficientes a11, …, amn. Um exemplo concreto de sistema de equações lineares poderia ser assim:

3x1+−2x2+4x3+5x4=15x1+4x2+2x3+−10x4=211x1+8x2+10x3+20x4=4

Mas, se algum aij é negativo, é costume escrever nesta forma:

3x1−2x2+4x3+5x4=15x1+4x2+2x3−10x4=211x1+8x2+10x3+20x4=4

Mas continua valendo que a12 = −2 e a24 = −10.

Notação matricial

Podemos converter esse sistema para a forma matricial. Definimos três matrizes no total: uma para os coeficientes, uma para o lado direito das equações e uma para as próprias incógnitas. As matrizes ficam assim:

A=(a11a12…a1na21a22…a2n⋮⋮⋱⋮am1am2…amn),b=(b1b2⋮bm),x=(x1x2⋮xn)

A matriz A contém os coeficientes das equações originais e a chamamos de “matriz dos coeficientes”, a matriz b contém os lados direitos das equações e a matriz x contém as incógnitas. Então podemos escrever o sistema de equações inteiro como a equação

Ax=b

Essa equação faz sentido, porque no lado esquerdo temos o produto de duas matrizes e no lado direito temos uma matriz, e as ordens delas são compatíveis. Podemos tentar multiplicar o lado esquerdo. Multiplicamos uma matriz de ordem m × n por uma matriz de ordem n × 1, então a matriz resultante será de ordem m × 1, terá m linhas e uma coluna. Por enquanto, está tudo certo.

Agora aplicamos o algoritmo clássico de multiplicação de matrizes, começando com a primeira linha da matriz A e com a primeira coluna da matriz x (não há outra coluna). Obtemos:

a11x1+a12x2+…+a1nxn

Com essa soma, obtivemos o primeiro elemento da nova matriz, ou seja, o valor b1. Substituindo tudo na equação, temos:

a11x1+a12x2+…+a1nxn=b1

Obtivemos exatamente a primeira equação do nosso sistema de equações. Do mesmo modo continuaríamos com as outras linhas. Definimos ainda a matriz ampliada do sistema, que é a matriz dos coeficientes à qual acrescentamos a matriz b assim:

B=(A⋮b)=(a11a12…a1nb1a21a22…a2nb2⋮⋮⋱⋮⋮am1am2…amnbm)

Em seguida, vai nos interessar o posto dessa matriz ampliada do sistema. Representamos o posto da matriz original (não ampliada) dos coeficientes por posto(A), e o da ampliada, por posto(A⋮b). Qual é a relação entre esses postos? Se acrescentamos à matriz A a coluna b, podem ocorrer duas situações: se essa nova coluna é uma combinação linear das colunas da matriz A, o posto não muda. Se não era uma combinação linear, o posto será uma unidade maior.

Teorema de Rouché-Capelli

O posto da matriz ampliada será importante para nós, porque vale um teorema importante, o teorema de Rouché-Capelli: O sistema de equações Ax = b tem solução se e somente se o posto da matriz dos coeficientes é igual ao posto da matriz ampliada do sistema: posto(A)=posto(A⋮b).

Vamos tentar mostrar que esse teorema vale pelo menos em um sentido — se x = β é solução da equação Ax = b, então posto(A)=posto(A⋮b). Se esses postos são iguais, deve valer que a coluna b é uma combinação linear das colunas da matriz A. Vamos detalhar. A matriz β tem a forma:

β=(β1β2⋮βn)

E afirmamos que, se substituirmos x1 por β1 etc., todas as equações do sistema de equações serão verdadeiras:

a11β1+a12β2+…+a1nβn=b1a21β1+a22β2+…+a2nβn=b2⋮⋮⋮⋮⋮⋮⋮⋮⋮am1β1+am2β2+…+amnβn=bm

Agora reescrevemos um pouco as equações. Se olharmos para a primeira colunazinha, veremos que ali sempre temos os valores a11, a21, …, am1, que multiplicamos pelo mesmo valor β1. Se isolarmos a colunazinha em uma matriz separada, podemos colocar em evidência o valor β1 diante da matriz, obtendo:

(a11β1a21β1⋮am1β1)=β1(a11a21⋮am1)

Agora podemos reescrever assim o sistema inteiro:

β1(a11a21⋮am1)+β2(a12a22⋮am2)+…+βn(a1na2n⋮amn)=(b1b2⋮bm)

Para que essa equação valha, a matriz do lado direito deve ser uma combinação linear das matrizes-coluna do lado esquerdo, e os valores β1, …, βn devem ser os coeficientes dela.

Com isso, demonstramos que, para que a matriz β seja solução do sistema de equações, o lado direito do sistema, a matriz b, deve ser uma combinação linear das colunas da matriz dos coeficientes. E, se a coluna da matriz b é uma combinação linear das colunas da matriz A, então a matriz A deve ter o mesmo posto que a matriz (A⋮b).

Ao mesmo tempo, vale que, se posto(A)=posto(A⋮b)=k e k é igual ao número de incógnitas, ou seja, k = n, o sistema tem exatamente uma solução. Se vale k < n, o sistema tem infinitas soluções e, para escrevê-las, precisaremos de n − k parâmetros.

Operações elementares com as linhas

Além disso, pela equação anterior se vê que, se realizarmos na matriz (A⋮b) operações elementares com as linhas, a solução do sistema de equações não muda. Se somássemos à segunda linha a primeira linha, obteríamos a equação:

β1(a11a21+a11⋮am1)+β2(a12a22+a12⋮am2)+…+βn(a1na2n+a1n⋮amn)=(b1b2+b1⋮bm)

Se detalharmos a segunda linha, obtemos:

β1(a21+a11)+β2(a22+a12)+…+βn(a2n+a1n)=b1+b2

Eliminamos os parênteses:

β1a21+β1a11+β2a22+β2a12+…+βna2n+βna1n=b1+b2

Agora só reorganizamos um pouco:

(β1a11+β2a12+…+βna1n)+(β1a21+β2a22+…+βna2n)=b1+b2

No primeiro parêntese temos, na verdade, a primeira linha e, no segundo parêntese, a segunda linha original. Sabemos que a primeira linha é igual ao valor b1 e a segunda linha é igual ao valor b2. Sabemos isso pelas equações originais. Então, no lugar do primeiro parêntese podemos escrever b1 e, no lugar do segundo, b2:

b1+b2=b1+b2.

Do mesmo modo, podemos multiplicar uma linha inteira por uma constante c ≠ 0 sem alterar a validade do sistema de equações. Atenção, porém: não podemos realizar operações com colunas.

Sistema com uma única solução

Com um exemplo, vamos mostrar como resolver um sistema que tem exatamente uma solução. Considere este sistema de equações:

2x1+3x2+7x3=473x1+8x2+x3=503x2+3x3=27

A matriz A ficará assim:

A=(237381033)

Agora calculamos o posto dessa matriz:

(237381033)∼(69216162033)∼(692107−19033)∼(6921021−5702121)∼(6921021−570078)

Vemos que posto(A) = 3. Disso decorre também que posto(A⋮b)=3, porque uma matriz de ordem 3 × 4 pode ter posto máximo três. O sistema tem três incógnitas, então o sistema inteiro tem exatamente uma solução. Como a encontramos?

Eliminação de Gauss

A eliminação de Gauss consiste em transformar a matriz ampliada do sistema em forma escalonada. Ao transformarmos a matriz em forma escalonada, conseguiremos descobrir o valor de uma incógnita. Depois de encontrarmos o valor de uma incógnita, podemos começar a substituir esse valor nas outras equações.

A matriz (A⋮b) tem a forma:

(A⋮b)=(237473815003327)

Vamos transformá-la em forma escalonada. Faremos as mesmas operações do passo anterior, só que desta vez as aplicaremos também à quarta coluna:

(237473815003327)∼(6921141616210003327)∼(692114107−19−4103327)∼(6921141021−57−12302121189)∼(6921141021−57−1230078312)

Como realizamos operações elementares com as linhas, essa matriz descreve um sistema de equações equivalente — um sistema que tem o mesmo conjunto solução. Vamos reescrever o sistema da notação matricial para a usual, usando os coeficientes da nossa última matriz:

6x1+9x2+21x3=14121x2−57x3=−12378x3=312

Agora partimos da última equação. Ela nos diz que 78x3 = 312. Que valor deve ter x3? É uma simples equação do 1º grau, então

78x3=312x3=31278x3=4

Temos o valor da primeira incógnita, x3. Substituímos esse valor na segunda equação. Obtemos a equação:

21x2−57x3=−12321x2−57⋅4=−12321x2=−123+22821x2=105x2=10521x2=5

E, por fim, substituímos os dois valores calculados na primeira equação:

6x1+9x2+21x3=1416x1+9⋅5+21⋅4=1416x1+45+84=1416x1=12x1=2

O valor da incógnita x1 é igual a dois. Neste momento, temos a solução completa. Podemos escrever que x = (x1, x2, x3) = (2, 5, 4). Como se vê, a eliminação de Gauss é um método simples, mas eficaz, para resolver um sistema de equações lineares.

Infinitas soluções

Se posto(A)=posto(A⋮b)=k e k é menor que o número de incógnitas, k < n, então o sistema tem infinitas soluções. Vamos tomar um exemplo simples:

4x1+x2=512x1+3x2=15

Vemos que a segunda equação é igual ao triplo da primeira equação. Se multiplicarmos a primeira equação por menos três e somarmos à segunda linha, obtemos:

4x1+x2=50x1+0x2=0

Neste momento, o número de incógnitas é igual a dois, n = 2, o posto da matriz dos coeficientes é igual a um, posto(A) = 1, e o posto da matriz ampliada também é igual a um, posto(A⋮b)=1. Vale que k < n, então o sistema tem infinitas soluções. Como encontramos essas soluções?

Olhamos para as equações. A segunda equação é satisfeita para todos os x1, x2, então só nos interessa a primeira equação. Na verdade, todos os pares x1, x2 que satisfazem a primeira equação formarão o conjunto de todas as soluções do sistema de equações dado. Então calculamos a solução da equação 4x1 + x2 = 5. Transformamos a equação na forma x2 = 5 − 4x1.

Em seguida, vamos nos ajudar com um parâmetro. Escolhemos um parâmetro t. Agora vamos tentar expressar o par x1, x2 por meio do parâmetro t, de modo que possamos escrever que, por exemplo, todos os pares da forma (t, t + 2) são solução do sistema. Ou seja, pares como (1, 3) ou (14, 16). O objetivo é expressar o valor de uma incógnita por meio de outra incógnita, de modo que saibamos que, quando o valor de x1 for tal e tal, o valor de x2 será, por exemplo, três vezes maior. Como ponto de partida, escolhemos x1.

Impomos, então, a igualdade t = x1. Se a incógnita x1 tiver o valor t, que valor terá a incógnita x2? Vemos isso pela equação anterior: terá o valor x2 = 5 − 4x1. Ou seja, se x1 = t, então x2 = 5 − 4t. Podemos escrever, então, que todos os pares da forma (t, 5 − 4t) são solução do sistema de equações.

Podemos conferir. Se substituirmos t por um, temos: x1 = 1 e x2 = 5 − 4 · 1 = 1. Substituímos esses valores na equação: 4 · 1 + 1 = 5, a equação confere.

Se substituirmos t por cinco, temos: x1 = 5, x2 = 5 − 4 · 5 = −15. Substituindo na equação, temos: 4 · 5 − 15 = 5. Confere.

Solução geral e solução particular

A solução de um sistema de equações escrita por meio de parâmetros é chamada de solução geral. Uma solução concreta, isto é, quando substituímos os parâmetros por números concretos, é chamada de solução particular.

Então, do exemplo anterior: (x1, x2) = (t, 5 − 4t) é a solução geral do sistema de equações, e (1, 1) e (5, −15) são soluções particulares.

Se temos a matriz dos coeficientes A e a matriz ampliada do sistema (A⋮b) e os postos delas são iguais, posto(A)=posto(A⋮b)=k, então, para expressar a solução geral, precisamos de n − k parâmetros. Repare que, se n = k, precisamos de zero parâmetros. Esse é o caso em que o sistema tem exatamente uma solução — nenhum parâmetro é necessário.

Exemplo

Resolva o seguinte sistema de equações lineares:

3x1−2x2+4x3+5x4=15x1−4x2+2x3+10x4=211x1−8x2+10x3+20x4=4

Primeiro calculamos o posto das matrizes A e (A⋮b). Transformamos a matriz ampliada em forma escalonada, e dela descobrimos o posto das duas matrizes. Vamos às operações. (Multiplicamos a primeira linha por dois, depois somamos a primeira e a segunda linhas e subtraímos esse resultado da terceira linha. Em seguida, subtraímos da segunda linha a primeira linha e, por fim, dividimos a primeira linha de volta por dois.)

(A⋮b)=(3−24515−4210211−810204)∼(6−481025−4210211−810204)∼(6−481025−4210200000)∼(6−48102−10−60000000)∼(3−2451−10−60000000)

Dessa última matriz já vemos que posto(A)=posto(A⋮b)=k=2. O número de incógnitas é n = 4, então precisaremos de n − k = 2 parâmetros. Vamos escrever as equações de acordo com a última matriz:

3x1−2x2+4x3+5x4=1−x1−6x3=0

Da segunda equação obtemos

−x1−6x3=0−x1=6x3x1=−6x3

Escolhemos o primeiro parâmetro, t = x3. Então podemos escrever que x1 = −6t. Podemos substituir isso na primeira equação:

3(−6t)−2x2+4t+5x4=1−18t−2x2+4t+5x4=1−14t−2x2+5x4=1−2x2=1+14t−5x4x2=−12−7t+52x4

Entra em jogo o segundo parâmetro, s = x4. Então podemos escrever que

x2=−12−7t+52s

Com isso obtemos a solução geral do sistema:

(x1,x2,x3,x4)=(−6t,−12−7t+52s,t,s)

Podemos experimentar algumas soluções particulares. Por exemplo, s = t = 0. Então obtemos: (x1,x2,x3,x4)=(0,−12,0,0). Substituindo nas equações originais:

3x1−2x2+4x3+5x4=15x1−4x2+2x3+10x4=211x1−8x2+10x3+20x4=4

obtemos as equações:

0−2⋅(−12)+0+0=10−4⋅(−12)+0+0=20−8⋅(−12)+0+0=4

Depois de eliminar os elementos nulos e de multiplicar as frações, obtemos:

1=12=24=4

Vamos tentar outra solução particular: s = 2, t = 1. Então obtemos:

x1=−6x2=−12−7+5=−52x3=1x4=2

Substituímos nas equações:

−18−2(−52)+4+10=1−30−4(−52)+2+20=2−66−8(−52)+10+40=4

Simplificamos:

−4+5=1−8+10=2−16+20=4

E, por fim, obtemos:

1=12=24=4

Parece que a solução geral obtida está correta.

Capítulos: Sistemas de equações lineares, Regra de Cramer, Sistemas lineares homogêneos