2

Объясните, пожалуйста, как применять метод Монте-Карло для двойного интеграла с бесконечными пределами интегрирования. Как делать для определенного понятно https://ru.wikipedia.org/wiki/Метод_Монте-Карло. Но как сделать для двойного?

P.S.: решать не нужно, просто подскажите алгоритм, пожалуйста

7
  • Вас смущает именно двойной интеграл, или именно бесконечные пределы интегрирования?
    – Harry
    18 дек 2018 в 4:22
  • @Harry и то, и то) 18 дек 2018 в 4:41
  • 1
    Ну, для двумерной области - как и для одномерной - только и того, что точки оказываются в трехмерном пространстве. С бесконечными пределами сложнее - как варианты, оценивать хвосты по асимптотике и отрезать, или искать преобразование, которое бы свело интеграл к ограниченной области, но без бесконечных значений самой подынтегральной функции (не уверен, что такое отыщется в общем случае), или как-то потом убирать сингулярности...
    – Harry
    18 дек 2018 в 4:47
  • сомневаюсь, что методом МК можно вычислить этот интеграл с бесконечными пределами эффективно. 18 дек 2018 в 8:30
  • Harry прав. нужно заменой переменных перейти от интеграла с бесконечными пределами к интегралам с конечными пределами 18 дек 2018 в 8:46

1 ответ 1

2

Не более чем иллюстрации для - применение метода Монте-Карло к вычислению объема шара единичного радиуса (двойной интеграл по области x^2+y^2 <= 1 функции 2*sqrt(1-x^2-y^2)).

Для наглядности она из полученного значения вычисляет число "пи".

#include <random>
#include <iostream>

using namespace std;



int main(int argc, const char * argv[])
{
    default_random_engine e;
    uniform_real_distribution r(-1.0,1.0);


    unsigned long long total = 0, vol = 0;
    for(unsigned long long i = 1; ; ++i)
    {
        double x = r(e), y = r(e), z = r(e);
        ++total;
        if (x*x+y*y > 1.0) continue; // Выход за область интегрирования

        if (z*z < 1 - x*x - y*y) ++vol;  // Точка входит в область

        if (i%1000000 == 0)
            cout << 6.0*vol/total << endl;
    }

}

Запускайте и смотрите сами :)

Ваш ответ

By clicking “Отправить ответ”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Всё ещё ищете ответ? Посмотрите другие вопросы с метками или задайте свой вопрос.