如何通过我自己的函数制作逆矩阵

时间:2019-05-14 15:32:52

标签: r function for-loop matrix-inverse

我正在编写代码(带有R)来解决带for循环和函数的逆矩阵。 编写for-loop子句存在问题。

我正在用我的函数求解带有协因子的逆矩阵。

##This is my example matrix

    A <- matrix(data = c(1, 3, -4, 5, -2, 1, -1, 3, 3, -2, 2, -1, 4,-1, -3, -1),
            byrow = TRUE,
            nrow = 4, n col = 4)




##This is co factor function
library(functional)

cof = function(A, i, j)
{ stop if not(n col(A) == n row(A))

  det(A[-i,-j])*(-1)^(i+j)
}                         

## this is the problem-writing for loop

n <- nrow(A)
f <- c(n)
f[1] <- 1

for (i in 1:n) {
  for (j in 1:n) {
  f[i,j] <- cof(A, i, j)
  }
}
f

我如何分配 我的辅因子函数(A,i,j)到a(i,j) 从11到44 带有for循环子句。

请帮助

0 个答案:

没有答案