0Pricing
R Academy · Aula

Resolvendo sistemas lineares com solve()

Encontre soluções para sistemas Ax = b e calcule inversas de matrizes.

Resolvendo sistemas lineares com solve() é uma aula grátis de R Academy no CoddyKit. Esta é a aula 2 de 4. Você pode ler a aula completa abaixo gratuitamente — depois pratica ao vivo no navegador com um editor de código integrado e um tutor de IA 24/7. Faz parte do caminho de aprendizado de R Academy, e seu progresso é sincronizado entre a web e o app CoddyKit. O curso de R Academy inclui 4 aulas no total.

Sistemas lineares: Ax = b

Um sistema de equações lineares pode ser escrito como Ax = b, em que A é uma matriz de coeficientes, x é o vetor de incógnitas e b é o lado direito. Resolver x analiticamente significa calcular x = A⁻¹b.

# System of equations:
# 2x + y = 5
# x + 3y = 7

# Matrix form: A %*% x = b
A <- matrix(c(2, 1,
              1, 3), nrow = 2, byrow = TRUE)
b <- c(5, 7)

# What are A and b?
print(A)
print(b)
cat('We want to find x such that A %*% x = b')

solve(A, b): a solução direta

solve(A, b) resolve Ax = b para x. Internamente, usa a decomposição LU, que é numericamente mais estável e eficiente do que calcular A⁻¹ explicitamente e depois multiplicá-la por b.

A <- matrix(c(2, 1,
              1, 3), nrow = 2, byrow = TRUE)
b <- c(5, 7)

# Solve Ax = b
x <- solve(A, b)
print(x)  # x[1] = ?, x[2] = ?

# Verify: A %*% x should equal b
residual <- A %*% x - b
print(residual)  # Should be near zero

# Manual check:
# 2*(8/5) + (9/5) = 16/5 + 9/5 = 25/5 = 5 ✓
# 1*(8/5) + 3*(9/5) = 8/5 + 27/5 = 35/5 = 7 ✓

solve(A): a inversa da matriz

Chamar solve(A) com apenas um argumento retorna a inversa de A: A⁻¹, tal que A %*% A⁻¹ = I. Evite usá-la para resolver Ax=b — use solve(A,b) diretamente, pois é mais rápido e mais estável.

A <- matrix(c(4, 3,
              3, 2), nrow = 2, byrow = TRUE)

# Compute inverse
A_inv <- solve(A)
print(A_inv)

# Verify: A %*% A_inv = I
A %*% A_inv  # Should be identity matrix
round(A %*% A_inv, 10)

# Also A_inv %*% A = I
round(A_inv %*% A, 10)

# det(A) != 0 required for invertibility
det(A)  # -1 (nonzero, so invertible)

Verificando sua solução

Sempre verifique a solução calculando A %*% x - b. Devido à aritmética de ponto flutuante, o resíduo não será exatamente zero, mas deverá estar próximo do épsilon da máquina (~1e-15). Use norm() para obter uma única magnitude do resíduo.

A <- matrix(c(3, -1, 2,
              1,  4, 0,
             -2,  1, 5), nrow = 3, byrow = TRUE)
b <- c(1, 2, 3)

# Solve
x <- solve(A, b)
cat('Solution x:\n'); print(x)

# Residual check
residual <- A %*% x - b
cat('Residual vector:\n'); print(residual)

# Residual norm (should be near 0)
resid_norm <- sqrt(sum(residual^2))
cat('Residual norm:', resid_norm, '\n')
# Expected: something like 2e-16

Número de condição: kappa()

O número de condição de A mede o quanto a solução é sensível a perturbações em b. Um número de condição grande significa que pequenos erros em b causam grandes erros em x — o sistema é mal condicionado.

# Well-conditioned matrix
A_good <- matrix(c(2, 1, 1, 3), nrow = 2, byrow = TRUE)
kappa(A_good)  # Small -> good

# Ill-conditioned (nearly singular) matrix
A_bad <- matrix(c(1.000, 1.001,
                  1.001, 1.002), nrow = 2, byrow = TRUE)
kappa(A_bad)  # Very large -> bad!

# Rule of thumb: kappa > 1/machine_epsilon is trouble
.Machine$double.eps  # ~2.2e-16

# For A_bad: you lose about log10(kappa) digits of precision
cat('Digits lost:', log10(kappa(A_bad)), '\n')

backsolve() para matrizes triangulares superiores

backsolve(R, b) resolve Rx = b quando R é triangular superior, usando substituição regressiva. É muito mais rápido do que solve() geral para sistemas triangulares — O(n²) contra O(n³).

# Upper triangular system: Rx = b
# 2x + 3y + z = 14
#      5y + 2z = 13
#           4z = 8
R <- matrix(c(2, 3, 1,
              0, 5, 2,
              0, 0, 4), nrow = 3, byrow = TRUE)
b <- c(14, 13, 8)

# Solve using back-substitution
x <- backsolve(R, b)
print(x)  # z=2, y=(13-4)/5=9/5, x=(14-3*9/5-2)/2

# Verify
all.equal(as.vector(R %*% x), b)  # TRUE

# Compare with general solve
x_general <- solve(R, b)
all.equal(x, x_general)  # TRUE (same result)

forwardsolve() para matrizes triangulares inferiores

forwardsolve(L, b) resolve Lx = b quando L é triangular inferior, usando substituição progressiva. Complementa backsolve() e, juntas, essas funções formam a base dos solucionadores baseados na decomposição LU.

# Lower triangular system: Lx = b
# 3x            = 6
# 2x + 4y       = 10
# x  + 2y + 5z  = 16
L <- matrix(c(3, 0, 0,
              2, 4, 0,
              1, 2, 5), nrow = 3, byrow = TRUE)
b <- c(6, 10, 16)

# Solve using forward-substitution
x <- forwardsolve(L, b)
print(x)  # x=2, y=(10-4)/4=1.5, z=(16-2-3)/5=2.2

# Verify
all.equal(as.vector(L %*% x), b)  # TRUE

# Use case: solving L*U*x = b
# forwardsolve(L, b) -> y, then backsolve(U, y) -> x

Vários lados direitos

solve(A, B), quando B é uma matriz, resolve AX = B simultaneamente para todas as colunas de B. Isso é mais eficiente do que chamar solve(A, b) separadamente para cada coluna.

A <- matrix(c(2, 1,
              1, 3), nrow = 2, byrow = TRUE)

# Solve for two right-hand sides simultaneously
B <- matrix(c(5, 7,   # first system
              3, 1),  # second system
            nrow = 2, byrow = TRUE)

# X[:,1] solves Ax = B[:,1]
# X[:,2] solves Ax = B[:,2]
X <- solve(A, B)
print(X)

# Verify both solutions
A %*% X  # Should equal B
all.equal(A %*% X, B)  # TRUE

Detectando matrizes singulares

Chamar solve(A) para uma matriz singular gera um erro. Verifique det(A) ou rcond(A) (recíproco do número de condição) antes de resolver. Use tryCatch() para criar código robusto.

# Singular matrix (rows are linearly dependent)
S <- matrix(c(1, 2,
              2, 4), nrow = 2, byrow = TRUE)
det(S)     # 0 -> singular
kappa(S)   # Inf

# Safe solve with tryCatch
safe_solve <- function(A, b) {
  tryCatch(
    solve(A, b),
    error = function(e) {
      cat('Matrix is singular or nearly so:\n')
      cat(e$message, '\n')
      return(NULL)
    }
  )
}

result <- safe_solve(S, c(1, 2))
print(result)  # NULL

Mínimos quadrados com solve()

Para sistemas sobredeterminados (com mais equações do que incógnitas), não há uma solução exata. A solução de mínimos quadrados minimiza ||Ax - b||². Ela resolve as equações normais: A'Ax = A'b.

# Overdetermined: 4 equations, 2 unknowns (y = a + b*x)
set.seed(1)
x_vals <- c(1, 2, 3, 4)
y_vals <- c(2.1, 4.0, 5.9, 8.2)  # approx y = 0 + 2x
A <- cbind(1, x_vals)  # Design matrix (4x2)
b <- y_vals

# Normal equations: (A'A) beta = A'b
AtA <- crossprod(A)
Atb <- crossprod(A, b)
beta_ols <- solve(AtA, Atb)
print(beta_ols)  # Intercept ~0.1, slope ~2.0

# Residual sum of squares
y_hat <- A %*% beta_ols
rss <- sum((b - y_hat)^2)
cat('RSS:', rss)

Usando qr.solve() para obter estabilidade

qr.solve(A, b) é numericamente mais estável do que solve() para sistemas mal condicionados ou sobredeterminados. Ele usa a decomposição QR em vez da LU. lm() usa esse recurso internamente.

# For overdetermined system, qr.solve is preferred
x_vals <- c(1, 2, 3, 4, 5)
y_vals <- c(1.9, 4.1, 6.0, 7.8, 10.1)
A <- cbind(1, x_vals)
b <- y_vals

# qr.solve handles overdetermined systems directly
beta_qr <- qr.solve(A, b)
print(beta_qr)  # intercept, slope

# Equivalent to:
beta_lm <- coef(lm(y_vals ~ x_vals))
all.equal(beta_qr, beta_lm, check.names = FALSE)  # TRUE

# For well-determined square systems, solve() is fine
# For overdetermined or ill-conditioned: use qr.solve()

Verificação rápida

Teste sua compreensão da resolução de sistemas lineares em R.

Recapitulação: resolução de sistemas lineares

Principais conclusões: solve(A, b) resolve Ax=b diretamente (é a opção preferida). solve(A) calcula A⁻¹ (evite usá-lo para resolver sistemas). backsolve(R,b) e forwardsolve(L,b) são rápidos para sistemas triangulares. kappa(A) mede o condicionamento — valores grandes significam soluções instáveis. Para sistemas sobredeterminados, use qr.solve(). Sempre verifique com A %*% x - b.

# Summary of solve() functions:
A <- matrix(c(3, 1, 1, 2), nrow = 2)
b <- c(9, 8)

# Solve Ax = b
x <- solve(A, b)
print(x)  # c(2, 3)

# Check conditioning
kappa(A)  # Small -> well-conditioned

# Verify
max(abs(A %*% x - b))  # Near zero

# For triangular systems:
R <- matrix(c(2, 3, 0, 4), nrow = 2, byrow = TRUE)
backsolve(R, c(8, 4))  # x=c(1, 1)

# For least squares (overdetermined):
# qr.solve(design_matrix, y)

Perguntas Frequentes

A aula “Resolvendo sistemas lineares com solve()” é grátis?

Sim — o texto completo de “Resolvendo sistemas lineares com solve()” é grátis para ler aqui na web. Para praticá-la interativamente (um editor de código integrado e um tutor de IA 24/7) e desbloquear o restante do curso de R Academy, atualize para CoddyKit PRO. O curso de R Academy inclui 4 aulas no total.

O que vou aprender em “Resolvendo sistemas lineares com solve()”?

Encontre soluções para sistemas Ax = b e calcule inversas de matrizes. Você pratica R Academy com código prático que executa diretamente no navegador, e um tutor de IA 24/7 responde suas dúvidas enquanto trabalha na aula.

Preciso ter experiência prévia para começar R Academy?

Nenhuma experiência prévia é necessária. R Academy no CoddyKit é estruturado para alunos iniciantes até avançados, então você pode começar aqui ou desde o início e aprender no seu ritmo. Esta é a aula 2 de 4.

Quanto tempo leva a aula “Resolvendo sistemas lineares com solve()”?

A maioria das aulas CoddyKit leva cerca de 5–10 minutos. Cada uma é compacta e interativa, então você faz progresso constante e retoma exatamente de onde parou entre web e app.

Posso escrever e executar código nesta aula de R Academy?

Sim. Cada aula de R Academy inclui um editor de código integrado, então você escreve e executa código real direto no navegador e recebe feedback de IA instantaneamente — nenhuma configuração local necessária.

Todas as aulas deste curso

  1. Multiplicação de matrizes e determinantes
  2. Resolvendo sistemas lineares com solve()
  3. Autovalores e autovetores
  4. Decomposições SVD, QR e de Cholesky
← Voltar para R Academy