Шаги для имплементации интегралов методом Монте Карло в Матлаб
22.01.2020
Подскажите пожалуйста, какие шаги (в общем) необходимо сделать для имплементации решения интеграла в Матлаб методом монте карло?
Подскажите пожалуйста, какие шаги (в общем) необходимо сделать для имплементации решения интеграла в Матлаб методом монте карло?
Комментарии
Привет! Схема следующая:
Задаем количество случайных чисел
>>N = 1e6;
Генерируем вектор случайных чисел с равномерным распределением
>>x = a + rand(1,N)*(b-a);
Считаем интеграл
S = (b-a)*sum(exp(-(x-3).^2))/N;
Можно сделать проверку
S1 = integral(@(x) exp(-(x-3).^2),a,b);