R相当于scipy.integrate.simps()?

问题描述 投票:1回答:1

我想在函数上进行一些数值积分,其中我只有函数的样本,并且不能生成任意新的样本。我理解在python宇宙中这可以通过scipy.integrate.simps()来实现,但是我的工作流程目前在R中。有关R函数/包的任何建议来实现这一点吗?

python r scipy numerical-methods integral
1个回答
1
投票

假设你有一个x值向量和一个y值向量,你可以很容易地应用Simpson的规则。

simpson <- function(x, y)
{
    if(length(x) < 5)
        stop("Must have at least 5 values")
    if(length(x) %% 2 == 0)
        stop("Number of values must be odd")
    ord <- order(x)
    x <- x[ord]
    y <- y[ord]
    diffs <- diff(x)
    delta <- mean(diffs)
    if((max(diffs) - min(diffs))/delta > 1e-6)
        stop("X-values must be equally spaced")
    coefs <- c(1, 4, rep(c(2, 4), (length(x) - 3)/2), 1)
    sum(coefs*y)*delta/3
}

simpson(1:7, (1:7)^2
# [1] 114
© www.soinside.com 2019 - 2024. All rights reserved.