R中未知值的矩阵代数

Question

有一种方法可以在R中解决矩阵代数，就像在站点中可以做到的那样： https://www.symbolab.com/solver/matrix-calculator

我正试图在R上解决这个问题： image with aexample of what i'm trying to solve

我也试过使用ktsolve包，但我认为它不能做我想要的。

感谢您的时间。

Answer 1

根据您的目标，优化可能是

中最快速，最简单的选择

calculate <- function(x){
    A <- matrix(c(1 ,2 ,x), nrow = 1)
    B <- matrix(c(1, 2, 3, 4, 5, 2, 6, 2, 1), nrow = 3)
    C <- matrix(c(1, 2, x), nrow = 3)
    return(abs(A %*% B %*% C - 1))
}


optimize(calculate, c(-100, 100))