Подскажите пожалуйста, какие шаги (в общем) необходимо сделать для имплементации решения интеграла в Матлаб методом монте карло?
Привет! Схема следующая:
Задаем количество случайных чисел
>>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);
Ответы
Привет! Схема следующая:
Задаем количество случайных чисел>>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);