Определитель матрицы с символьными элементами в R - PullRequest
0 голосов
/ 23 мая 2018

У меня есть следующий код:

x1 <- c("x", "y")
x2 <- c("z", "w")
X <-data.frame(x1,x2)
A=as.matrix(X)

Матрица A выглядит следующим образом:

      x1  x2 
[1,] "x" "z"
[2,] "y" "w"

Как найти определитель A?Я хочу получить результат, подобный xw-yz.

Я пытался det(A), но он не работает.

В общем, я хочу иметь возможность найти определитель любой матрицы, элементы которойтипа символов.

1 Ответ

0 голосов
/ 11 сентября 2018

вам нужна библиотека символьных вычислений:

library(rSymPy)

x1 <- List(Var("x"), Var("y"))
x2 <- List(Var("z"), Var("w"))
A = Matrix(x1, x2)
Det <- function(x) Sym("(", x, ").det()")
Det(A)
# [1] "w*x - y*z"

x <- Var("x") 
y <- Var("y") 
z <- Var("z")
A = Matrix(List(y+z,z,y), List(z,z+x,x), List(y,x,x+y))
Det(A)
# [1] "4*x*y*z"

a <- Var("a") 
b <- Var("b") 
c <- Var("c")
A = Matrix(List(0,b-a,c-a), List(a-b,0,c-b), List(a-c,b-c,0))
Det(A)
# [1] "0"
...