r - R中包含矩阵时的函数积分

Question

我想集成一个看起来像的功能

f <- function(x) 1.96 * sqrt(t(c(1,x)) %*% m %*% c(1,x))

其中 m 是

m <- matrix(c(3.855, -0.206, -0.206, 0.01), nrow = 2, ncol = 2, byrow = T)

由于我的函数中的矩阵乘法会产生一个标量，因此对于 的任何值x，这只是 的一维积分f(x)。我怎样才能顺利解决这个问题？

标签： rmatrixvectormatrix-multiplicationnumerical-integration

Answer 1

只需使用integrateand Vectorize：

integrate(Vectorize(f), lower = 0, upper = 1)