Студопедия
Случайная страница | ТОМ-1 | ТОМ-2 | ТОМ-3
АвтомобилиАстрономияБиологияГеографияДом и садДругие языкиДругоеИнформатика
ИсторияКультураЛитератураЛогикаМатематикаМедицинаМеталлургияМеханика
ОбразованиеОхрана трудаПедагогикаПолитикаПравоПсихологияРелигияРиторика
СоциологияСпортСтроительствоТехнологияТуризмФизикаФилософияФинансы
ХимияЧерчениеЭкологияЭкономикаЭлектроника

Усовершенствованный метод Эйлера

Читайте также:
  1. I. Определение и проблемы метода
  2. I. ОПРЕДЕЛЕНИЕ И ПРОБЛЕМЫ МЕТОДА
  3. I. ОРГАНИЗАЦИОННО-МЕТОДИЧЕСКИЙ РАЗДЕЛ
  4. I. Экспертные оценочные методы
  5. II МЕТОДИЧЕСКИЕ РЕКОМЕНДАЦИИ ПО ПОДГОТОВКЕ К ПРАКТИЧЕСКОМУ ЗАНЯТИЮ
  6. II. Категории и методы политологии.
  7. II. МЕТОДИЧЕСКИЕ УКАЗАНИЯ ДЛЯ ВЫПОЛНЕНИЯ САМОСТОЯТЕЛЬНОЙ РАБОТЫ

Характерной особенностью метода является использование в качестве направления поиска каждой последующей точки численного решения касательной, определяемой в центре отрезка [ xk, xk +1]. Как видно на рис.4 для k = 0, последовательно выполняются два шага по методу Эйлера.

Первый из них используется для вычисления тангенса угла наклона касательной в средней точке отрезка [ xk, xk +1]

,

,

,

.

После этого делается второй шаг – для вычисления новой точки

,

.

Если тангенсы углов наклона касательных заменить на правую часть дифференциального уравнения f (x, y), то приведённые выше фор­мулы для выполнения одного шага по усовершенствованному методу Эйлера сведутся к следующим

, ,

, .

Этот метод является методом второго порядка точности. Он даёт меньшую погрешность численного решения на шаге h, чем метод Эйлера. Его абсолютная погрешность ε абс(xk +1, h) на каждом шаге пропорциональна величине h 3, а решение совпадает с истинным в случае, когда оно представимо квадратичной функцией y = a 1+ a 2 x + a 3 x 2. Однако это достигается тем, что его трудоёмкость увеличивается примерно в два раза, поскольку для одного шага приходится два раза вычислять значение правой части дифферен­циального уравнения.

Метод Рунге–Кутта (C.Runge 1895, W.Kutta 1901)

Под наименованием «метод Рунге–Кутта» принято подразумевать целое семейство методов численного решения задачи Коши, объединенных одной идеей выбора некоторого усреднённого направления поиска каждой новой точки приближённого решения. Алгоритм одного из этих методов: метода Рунге–Кутта 4-го порядка, строится по следующей схеме. Из начальной точки, как это можно видеть на рис.5 для k = 0, в направлении касательной делается шаг по методу Эйлера на величину h/ 2 и в полученной точке вычисляется тангенс угла наклона касательной

,

, ,

.

Из начальной точки по найденному направлению делается второй шаг на ту же величину и определяется тангенс угла наклона во второй средней точке

,

.

Третий шаг делается опять из начальной точки по последнему найденному направлению (), но на шаг h, где опять находится тангенс угла наклона касательной

,

,

.

Последний шаг из начальной точки на величину h делается по присущему этому методу осреднённому направлению и находится ордината следующей точки приближённого решения

.

Описанные действия, необходимые для получения координат каждой последующей точки приближённого решения, реализуются в виде следующей последовательности вычислений по формулам

, ,

,

, ,

.

Получающееся численное решение имеет на каждом шаге h погрешность ε абс(xk +1, h), которая пропорциональнавеличине h 5. Оно совпадает с истинным решением в случае, когда оно – многочлен четвёртого порядка y = a 1+ a 2 x + a 3 x 2+ a 4 x 3+ a 5 x 4.


Дата добавления: 2015-07-08; просмотров: 112 | Нарушение авторских прав


Читайте в этой же книге: Метод основывается на приведении исходного уравнения к форме | Относительная разница между значениями приближения корня на третьей и четвёртой итерациях составляет | Метод хорд | Контрольные задания | Справочная информация | Метод Гаусса с выбором главного элемента | Метод простых итераций | О выборе метода решения систем уравнений | Контрольные задания | Кусочно-линейная интерполяция |
<== предыдущая страница | следующая страница ==>
Справочная информация| Оценка погрешностей методов

mybiblioteka.su - 2015-2024 год. (0.017 сек.)