代数(线性？)的特征函数函数

Question

当k+j-s-i =0时，Ai，j，k= ps ( ps是给定的向量)

这可以由特征函数写成ps*i( k+j-s-i)，或者由Kronecker增量函数写成ps*i(0，k+j-s-i)。在R中有没有什么“内建”函数--我的意思是有没有内建的？

或者我必须自己写？

我认为，当f(x) =0时，构建函数返回1，否则返回0，这将非常有用，至少对于线性f(X)是这样。

Answer 1

我会把它重写为

A[i,j,k] = p[k+j-i] if that exists, otherwise 0

，然后可以将其实现为

p <- c(1,2,3,4,5)

pfun <- function(x) {
  if (x < 1 | x > length(p)) {
    0
  } else {
    p[x]
  }
}

n <- 5
A <- array(0, c(n, n, n)) 

for (i in 1:n) {
  for (j in 1:n) {
    for (k in 1:n) {
      A[i,j,k] <- pfun(k+j-i)
    }
  }
}

可能有一些比三重嵌套for循环更优雅的东西。

至于您询问的函数，就像下面这样简单

as.numeric(f(x)==0)

会起作用的。