Автор: Пользователь скрыл имя, 25 Декабря 2011 в 20:54, курсовая работа
При изучении численных методов основной задачей являлось преобразование известной математической модели к виду, допускающему эффективное численное решение, т. е. сведение всех математических понятий ( таких, как производная, интеграл, дифференциальное уравнение и т. п.) к последовательности элементарных арифметических операций. Для дальнейшего существенно, что вместе с алгоритмом метода мы всегда получали так называемую гарантированную оценку погрешности метода. Алгоритм численного метода считался состоятельным лишь в том случае, если малые погрешности, содержащиеся в исходных данных и внесенные в в процессе преобразований или вычислений, не влияли приемлемым образом на результат вычислений. Процесс отыскания решения при таком подходе являлся строго детерминированным, т. е. при безошибочном повторении неизменно приводил к одному и тому же результату.
Введение………………………………………………………………………….3
1. Метод Монте-Карло
1.1 Общая схема метода Монте-Карло…………………………………………4
1.2 Оценка погрешности метода Монте-Карло………………………………..4
2. Вычисление интегралов методом Монте-Карло.
2.1 Пример вычисления определенного интеграла методом Монте-Карло..6
2.2 Программа вычисления определенного интеграла методом Монте-Карло……………………………………………………………………………..8
Заключение……………………………………………………………………....9
Список использованной литературы…………………………………………..10
Приложение……………………………
2.2. Программа вычисления определенного интеграла методом Монте-Карло
Sub MetodMonte_Karlo()
Dim a, b, k, i, s, ri, n, otv As Single
a = Worksheets("Лист1"). Range("B3"). Value
b = Worksheets("Лист1"). Range("B4"). Value
k = b - a
i = 0
n = Worksheets("Лист1"). Range ("B5"). Value
Do
i = i + 1
ri = Rnd
x = a + (b - a) * ri
s = s + (1 + x)
Loop While i < n
otv = k * s / n
Worksheets ("Лист1"). Range ("H7"). Value = otv
End Sub
где n –
число испытаний ; ri – случайные числа.
Заключение
Теория метода Монте – Карло изучает способы выбора случайных величин X для решения различных задач, а также способы уменьшения дисперсии случайных величин.
Метод Монте-Карло используется очень часто, порой некритично и неэффективным образом. Он имеет некоторые очевидные преимущества:
а) Он не требует никаких предложений о регулярности, за исключением квадратичной интегрируемости . Это может быть полезным, так как часто очень сложная функция, чьи свойства регулярности трудно установить.
б) Он приводит к выполнимой процедуре даже в многомерном случае, когда численное интегрирование неприменимо, например, при числе измерений, большим 10.
в) Его легко применять при малых ограничениях или без предварительного анализа задачи.
Он обладает, однако, некоторыми недостатками, а именно:
а) Границы ошибки не определены точно, но включают некую случайность. Это, однако, более психологическая, чем реальная, трудность.
б) Статическая погрешность убывает медленно.
в) Необходимость
иметь случайные числа.
Список использованной литературы
1. Бережная
Е. В., Бережной В. И. Математические
методы моделирования
2. Гмурман В.Е. Руководство к решению задач по теории вероятностей и математической статистике: Учеб. пособие для студентов втузов. – 3-е изд., перераб. И доп. – М.: Высш. школа, 1979.- 157 с.
3. Гмурман В. Е. Теория вероятностей и математическая статистика. Учеб. пособие для втузов. Изд. 5-е, перераб. и доп. М., «Высш. школа», 1977.- 52 с.
4. Ермаков С. М. Методы Монте-Карло и смежные вопросы. М.: Наука, 1971.- 83 с.
5. Маделунг
Э. Математический аппарат
6. Математика.
Большой энциклопедический
7. Мюллер П., Нойман П., Шторм Р. Таблицы по математической статистике. - М.: Финансы и статистика, 1982. - 278 с.
8. Пирумов У. Г. Численные методы : Учеб.пособие для студ. вузов.- 3-е изд., испр.- М.: Дрофа, 2004.- 104 с.
9. Севастьянов Б. А. Курс теории вероятностей и математической статистики. – М.:Наука,1982.- 74 с.
10. Теннант-Смит
Дж. Бейсик для статистиков. - М.: Мир, 1988.
- 208 с.
Приложение.
Равномерно распределённые случайные числа
0,65551 | 0,01825 | 0,54442 | 0,20811 | 0,73431 | 0,373 | 0,998077 | 0,420728 | 0,99487 | 0,038575 |
0,2316 | 0,3123 | 0,694113 | 0,36796 | 0,31581 | 0,78228 | 0,298135 | 0,969085 | 0,90768 | 0,916715 |
0,87747 | 0,14457 | 0,056795 | 0,19437 | 0,61342 | 0,99741 | 0,013489 | 0,687094 | 0,32224 | 0,483566 |
0,19523 | 0,78549 | 0,403027 | 0,19367 | 0,96759 | 0,48122 | 0,115421 | 0,143864 | 0,3957 | 0,737419 |
0,59423 | 0,35817 | 0,623096 | 0,97778 | 0,50603 | 0,52275 | 0,032991 | 0,226661 | 0,75503 | 0,493881 |
0,84997 | 0,93027 | 0,143681 | 0,94082 | 0,33647 | 0,99103 | 0,39494 | 0,904233 | 0,24189 | 0,04593 |
0,44182 | 0,55239 | 0,659597 | 0,52065 | 0,39085 | 0,91565 | 0,948363 | 0,4897 | 0,06589 | 0,880581 |
0,87457 | 0,8749 | 0,537736 | 0,22535 | 0,63158 | 0,01294 | 0,54677 | 0,977355 | 0,89181 | 0,123295 |
0,22431 | 0,87399 | 0,548723 | 0,70763 | 0,01892 | 0,20878 | 0,304361 | 0,412275 | 0,23414 | 0,843135 |
0,75881 | 0,38493 | 0,083743 | 0,66527 | 0,15772 | 0,65642 | 0,671072 | 0,278115 | 0,45473 | 0,579424 |
Информация о работе Вычисления определенных интегралов методом Монте-Карло