<<
>>

7.1. Метод Эйлера (метод Рунге-Кутта 1-го порядка).

Разобьем [a, b] на n равных частей – элементарных отрезков, x0, x1,…,xn будем называть узлами сетки, h = (b-a)/n - шаг сетки.

Очевидно, что , ; , .

Заменим в уравнении (7.1) y’ в точке xi её приближенной оценкой – отношением приращений (это следует из определения производной):

Тогда получаем:

Отсюда формула Эйлера:

(7.4)

, – номер узла

Зная y0 в точке x0 (начальное условие) можно найти y1, затем, используя уже известные значения x1 и y1, вычислить x2 и y2 и так далее.

Рассмотрим геометрическую иллюстрацию метода Эйлера. В координатах (x,y) отобразим известные данные: отрезок [a,b] на оси Х и начальное условие y0 – точка А с координатами (a, y0). Отрезок [a,b] разобьем на n равных частей, получим узлы равномерной сетки a = x0, x1, x2, … , xn = b. Вычислим значения первой производной искомой функции в точке А, используя координату этой точки и исходное уравнение (7.3)

Полученное значение позволяет построить касательную к искомой функции в точке А. Эту касательную можно использовать для вычисления приближенного значения искомой функции в новом узле х1 (кривую y(x) заменяем на отрезком АВ на элементарном отрезке [x0, x1]).

Рис. 7.1. Геометрическая иллюстрация метода Эйлера.

Зная (x1,y1), можно аналогично получить новую точку (x2,y2) и т.д.

Из геометрической иллюстрации следует, что:

1. На каждом шаге есть погрешность (на рисунке это отрезок BD). Погрешность тем больше, чем больше шаг.

2. Ошибка может накапливаться.

Формула Эйлера (7.4) имеет погрешность метода

Для практического выбора h с целью обеспечения заданной точности решения задачи e применяется следующий прием.

Выполняются 2 расчета: с n и 2n узлами. Если полученные значения функции в во всех узлах отличаются не более чем на e, задача считается решенной. Если нет, число узлов вновь удваивают и опять сравнивают полученные значения функций.

Таким образом, расчет продолжается до достижения условия

(7.5)

Значение n может достигать большой величины – более 1000. Чтобы не печатать столько значений функции, в алгоритме решения ОДУ методом Эйлера нужно предусмотреть печать не всех рассчитанных значений, а только части их, например, 10-ти значений, распределенных равномерно по всему отрезку.

Рис. 7.2. Алгоритм расчета новой точки методом Эйлера:

Пример 7.1. Дано уравнение

Найти решение для отрезка [0; 1], если y(0) = 1.

Выберем n = 10, тогда шаг h =(1-0)/10 = 0,1.

Запишем уравнение в каноническом виде

Начальная точка x0 = 0, y0 = 1.

Вычислим первую точку

x1 = x0 + h = 0 + 0,1 = 0,1

Вычислим вторую точку

Аналогично нужно вычислить еще восемь точек (выбрано n=10).

Рис. 7.3. Алгоритм решения ОДУ 1-го порядка методом Эйлера.

<< | >>
Источник: Мухамадеев И.Г.. АЛГОРИТМЫ ВЫЧИСЛИТЕЛЬНОЙ МАТЕМАТИКИ. КУРС ЛЕКЦИЙ. 2007

Еще по теме 7.1. Метод Эйлера (метод Рунге-Кутта 1-го порядка).:

  1. 3.1. Имитационное моделирование развития комплекса предприятий
  2. 1.3.1 Математическое описание алгоритма модели движения НКА
  3. 2.4. Моделирование работы подвижного состава с использованием марковских случайных процессов
  4. 1.3 Исследование динамической системы с использованием пакета Matlab
  5. 1.5 Исследование динамической системы с использованием пакета Mathcad
  6. Метод Рунге – Кутта.
  7. 4.2. СОДЕРЖАНИЕ РАЗДЕЛОВ ДИСЦИПЛИНЫ
  8. 2.1. Рабочая программа (объем дисциплины 150 часов)
  9. 7.4. Метод Рунге – Кутта
  10. ПРИЛОЖЕНИЕ
  11. ОГЛАВЛЕНИЕ
  12. 7.1. Метод Эйлера (метод Рунге-Кутта 1-го порядка).
  13. 7.2. Модифицированный метод Эйлера (метод Рунге-Кутта 2-го порядка).
  14. 7.4. Метод Рунге-Кутта 4 порядка.
  15. Тема 8 Решение систем обыкновенных дифференциальных уравнений.
  16. 3. Явные методы типа Рунге - Кутта.
  17. 5. Практические способы оценки погрешности приближенного решения.
  18. 8. Построение начальных значений.
  19. 1.10.3. Распространение ошибок в начальных данных при решении обыкновенных дифференциальных уравнений.
  20. V. ЗАДАЧА КОШИ ДЛЯ ОДУ.