Resolver Sistema de Ecuaciones en R

En este artículo, discutiremos cómo resolver un sistema de ecuaciones en el lenguaje de programación R. 

La función solve() en lenguaje R se usa para resolver la ecuación. Aquí la ecuación es como a*x = b, donde b es un vector o array yx es una variable cuyo valor se va a calcular.

Sintaxis: solve(a, b)

Parámetros:

  • a: coeficientes de la ecuación
  • b: vector o array de la ecuación

Ejemplo 1: Resolución de la ecuación del sistema de tres ecuaciones

Given Equations:
x + 2y + 3z = 20  
2x + 2y + 3z = 100  
3x + 2y + 8z = 200

Matrix A and B for solution using coefficient of equations:
A->
1   2   3
2   2   3
3   2   8
B->
20
100
200

Para resolver esto usando dos arrays en R usamos el siguiente código:

R

# create matrix A and B using given equations
A <- rbind(c(1, 2, 3), 
           c(2, 2, 3), 
           c(3, 2, 8))
B <- c(20, 100, 200)
  
# Solve them using solve function in R
solve(A, B)

Producción:

80 -36 3.99999999999999

lo que significa que x=80, y=-36 y z=4 es la solución para ecuaciones lineales.

Ejemplo 2: Resolución de la ecuación del sistema de tres ecuaciones

Para obtener soluciones en forma de fracciones, usamos la biblioteca MASS en lenguaje R y ajustamos la función de resolución en fracciones.

Given Equations:
19x + 32y + 31z = 1110  
22x + 28y + 13z = 1406  
31x + 12y + 81z = 3040
Matrix A and B for solution using coefficient of equations:
A->
19   32   31
22   28   13
31   12   81
B->
1110
1406
3040

Para resolver esto usando dos arrays en R usamos el siguiente código:

R

# Load package MASS
library(MASS)
  
# create matrix A and B using given equations
A <- rbind(c(19, 32, 31), 
           c(22, 28, 31), 
           c(31, 12, 81))
B <- c(1110, 1406, 3040)
  
# Solve them using solve
# function wrapped in fractions
fractions(solve(A, B))

Producción:

[1] 159950/2243 -92039/4486  29784/2243

lo que significa que x=159950/2243, y=-92039/4486 y z=29784/2243 es la solución para la ecuación lineal anterior.

Ejemplo 3: Resolución de array inversa

R

# create matrix A and B using given equations
A <- matrix(c(4, 7, 3, 6), ncol = 2)
print(A)
  
print("Inverse matrix")
  
# Solve them using solve function in R
print(solve(A))

Producción:

     [,1] [,2]
[1,]    4    3
[2,]    7    6
[1] "Inverse matrix"
          [,1]      [,2]
[1,]  2.000000 -1.000000
[2,] -2.333333  1.333333

Publicación traducida automáticamente

Artículo escrito por mishrapriyank17 y traducido por Barcelona Geeks. The original can be accessed here. Licence: CCBY-SA

Deja una respuesta

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *