Двойной/тройной интеграл в R

Я пытаюсь решить приведенное ниже (изображение) уравнение в R для двойного интеграла по площади. Я использую следующий код с помощью следующего сообщения: Triple интеграл в R (как указать домен)

Не уверен, если я что-то упустил, как если бы я изменил значение A, это все равно дает мне те же результаты. Правильно ли следующий код описывает уравнение, которое я пытаюсь решить?

library(cubature)

PW <-  0.0144
Cw <- 400/(0.91*70)
A <- 1960                
w <- as.numeric(PW, Cw, A)

lower <- rep(0,3)
upper <- rep(1,3)

fxyz <- function(w) {
    x <- w[1]
    y <- w[2]
    z <- w[3]
    as.numeric(x <= y)*as.numeric(y <= z)
}

Int <- adaptIntegrate(f=fxyz,lowerLimit=lower,upperLimit=upper,doChecking=TRUE,
       maxEval=2000000,absError=10e-5,tol=1e-5)

спасибо, Крина

введите здесь описание изображения


r integral numerical-integration calculus
person Krina M    schedule 27.03.2017    source источник

Ответы (1)


arrow_upward
0
arrow_downward

w, который является локальным для определения fxyz, затеняет глобальный w, который вы определили выше. Изменение A повлияет на глобальный w, но это совершенно не имеет отношения к w, который является локальным для fxyz. Причина, по которой изменение A не влияет на значение интеграла, заключается в том, что ни подынтегральная функция, ни пределы интегрирования не зависят от A.

В ответ на ваш вопрос: «Правильно ли следующий код описывает уравнение, которое я пытаюсь решить?»: Код для самого интеграла кажется правильным, но код для fxyz, по-видимому, нет. Кажется, вы подразумеваете некоторую зависимость между этой функцией и параметрами Pw, Cw, A, но не объяснили, какой должна быть эта зависимость. Возможно, вы могли бы уточнить, либо здесь, либо в другом вопросе.

person John Coleman    schedule 27.03.2017
comment
Спасибо, Джон. Pw — проницаемость, Cw — концентрация, A — площадь поверхности кишечника. Я пытаюсь вычислить двойную интеграцию (Pw*Cw) по домену A, который представляет собой площадь поверхности кишечника. Не уверен, какой будет точная зависимость, основанная на приведенном выше уравнении. Я постараюсь углубиться в теорию поверхностных интегралов. Я не очень знаком с этим. - person Krina M; 28.03.2017