Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Vse_shpory_po_KMMM.docx
Скачиваний:
29
Добавлен:
21.04.2019
Размер:
1.91 Mб
Скачать

17.1. Классификация методов решения оду.

Рассмотрим наши методы относительно к так называемым шаговым (пошаговым) методам интегрирования ОДУ.

Шаг интегрирования может быть постоянным или переменным. Мы рассмотрим такие методы, где для получения решения в очередной момент времени нам требовалось знать решение только в предыдущий момент времени. Такие методы называются одношаговые.

Существуют методы, в которых для построения решения в данный момент времени нужно знать решение задачи в нескольких предшествующих моментах.

Многошаговый метод:

Многошаговый метод не может начать работу по известному лишь в одной точке начальному условию. То есть для его разгона требуется предварительно запустить какой-то одношаговый метод и с его помощью подготовить решение в нескольких точках. Эти методы еще называются несамостартующими или несаморазгоняющими.

В рассмотренных нами методах значение неизвестного решения выражалось в явном виде через решение в предшествующих точках. Поэтому такие методы называются явными.

Существуют и другие методы, где yi+1 – неизвестное решение – является решением какого-то уравнения, чаще всего нелинейного, связанного с видом правой части ОДУ, т.е. неизвестное решение неявно входит в какую-то зависимость. Такие методы называются неявными. В них дополнительно дополнительно придется решать дополнительное уравнение, в общем случае трансцендентное. Обычно такие уравнения решаются методом Ньютона. Такие методы имеют повышенную точность и соответственно позволяют значительно увеличить шаг по независимой переменной.

18.1. Понятие о методах типа Монте-Карло.

В 1943г. в исследовательских лабораториях Лос- Анджелесе при разработки ядерной бомбы возникла задача об определении глубины проникновения электронов в заданное вещество. Решить её не удалось. Тогда Станислав Улом и Ждон фон Нейманом предложили подход основную стратегию, которую используют игроки при игре в кости. Эта стратегия была стратегия была основана на поведении случайных величин. По имени города- это метод получил название Монте – Карло. В дальнейшим по традиции многие методы стали называться метод Монте – Карло.

Изобразим схему метода.

Заметим, что этот подход может использоваться, как для моделирования явлений имеющих в своей основе поведение случайных величин (такие величины называются - стохастическими) так и для явлений процессов, где случайные величины не присутствуют (детерминированные).

19.1. Задача Бюффона как пример использования случайных величин при решении детерминированной задачи.

В 18в. француз Граф Бюффон прадложил для вычисления с произволбным количеством знаков использовать наблюдение за следующим процессом.

Тонкая игла длинной 2l бросается случайным образом на поверхность разлинованную с шагом 2d.

Тогда ситуация пересечение иглы с линией будет соответствовать тому, что точка с координатой ( ) будет лежать в заштрихованной области.

Рассмотрим углы от 0 до в силу симметричности.

Таким образом рассмотрим прямоугольную область .

И рассмотрим область лежащая ниже , то вероятность возникновения пересечения будет равна отношению . В математическом плане если выбираем без предпочтения координату от 0 до и координату y из интервала , то это соответствует произвольной точки , и если выполняется условие , то точка будет находиться в точке .

Величины:

=

=

Тогда .

Как известно при увеличении следует что: , где есть частота появления события.

При увеличении n получаем значение с большим количеством знаков.

Реализация в Mathcad: runif(a,b,N) – выдаёт n случайных величин с равномерном законом распределения. На отрезке от 0 до W – rnd(W).

Алгоритм:

1.Задать число экспериментов n и ограничить цикл от 0 до n(что будет соответствовать n опытов).

2.Обращаемся к функции runif(0, ,N).

3.Генерировать случайные величины от 0 до t.

4.Проверить выполнения неравенства . Если выполняется то n+1.

5.После проверки используем функцию

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]