问R中的图积分

EN

这里发生了几件事。

1. 看起来您想要在r上进行集成，而不是x。integrate(...)在第一个参数的值上集成了f，因此您需要将参数的顺序切换到f。
2. integrate(...)返回一个包含多个元素的列表。其中之一是value，积分的值。因此，要获得这个结果，您需要参考integrate(...)$value。读文献资料。
3. 如果要在函数中使用plot(...)，则需要对函数进行矢量化。也就是说，它需要以一个向量作为参数，并返回一个长度相同的向量。integrate(...)函数是R中少数没有矢量化的函数之一。幸运的是，有一个特殊的函数，Vectorize(...)，旨在将标量值函数转换为向量化函数。

把这一切都卷起来：

# intergrand f
f <- function(r,x) x*exp(-r)   # order of arguments reversed
# integral
h <- function(x) integrate(f, lower=0, upper=Inf, x=x)$value
g <- Vectorize(h)
x <- seq(-2,4,.1)
plot(x,g(x), xlim=c(-2,4), xlab="x", ylab="y", col="red")

​

​

正如在另一个答案中所指出的，将x*exp(-r)与r集成在[0,Inf]上会产生x，所以我真的不认为这样做有什么意义。然而，这就是代码的工作方式。

您的h函数返回一个包含许多变量的列表，而不是一个可以绘制的点。

不清楚您在做什么，因为h(x)只是返回x，因为它是在R上进行集成，而R是集成到1的。

修改后的代码如下所示：

h<-function(x){
  x
}
plot(h, xlim=c(-2,4), xlab="x", ylab="y", col="red")

请注意，xlim和函数都已更改。如果您想要积分x而不是R，请注意积分是

.5*exp(-r)*(x^2)

这也许是个更好的方法。