Как использовать метод Монте-Карло в Matlab?

Мне нужно проверить, что pi = интеграл от 0 до 2 от (sqrt (4-x ^ 2)) dx, используя метод Монте-Карло и 10 ^ 7 случайных чисел. Я немного смущен тем, как реализовать это в Matlab. Может ли кто-нибудь помочь мне?


matlab integration montecarlo
person SKLAK    schedule 03.02.2013    source источник
comment
Я пытаюсь сделать то, что вы только что опубликовали прямо сейчас.   -  person SKLAK    schedule 03.02.2013

Ответы (1)


arrow_upward
4
arrow_downward

Интеграл описывает четверть окружности. Чтобы вычислить число пи, выполните следующие действия:

  1. Сгенерировать 10 ^ 7 точек с координатами (x, y) в [0 2] x [0 2]
  2. Подсчитайте точки со свойством x ^ 2 + y ^ 2 ‹= 4. Пусть n будет их количеством.
  3. Вычислите значение числа пи: pi = 4 * n / 10 ^ 7
person kaspersky    schedule 03.02.2013
comment
хм, я попробовал это, и все точки, кажется, оказываются внутри круга, поэтому значение числа пи равно 4. Может быть, я делаю что-то не так. - person SKLAK; 03.02.2013
comment
неважно, я не считал очки в нужной области. Позвольте мне попробовать это снова. - person SKLAK; 03.02.2013