Модифицированные методы Эйлера

6.3 Модифицированные методы Эйлера

Первый модифицированный метод Эйлера. Суть этого метода состоит в следующем. Сначала вычисляются вспомогательные значения искомой функции y в точках t = t_i +  с помощью формулы:

y = y_i +  f_i = y_i +f(t_i, y_i).

Затем находится значение правой части уравнения (6.1) в средней точке

f = f(t, y)

и затем полагается

y_i+1 = y_i + h f, i = 0, 1, …, n – 1. (6.12)

Формулы (6.12) являются расчетными формулами первого модифицированного метода Эйлера.

Первый модифицированный метод Эйлера является одношаговым методом со вторым порядком точности

Второй модифицированный метод Эйлера – Коши. Суть этого метода состоит в следующем. Сначала вычисляются вспомогательные значения

 = y_i + h f(t_i, y_i). (6.13)

Затем приближения искомого решения находятся по формуле:

y_i+1 = y_i + [f(t_i, y_i) + f(t_i+1, )],  i = 0, 1, …, n – 1. (6.14)

Формулы (6.14) являются расчетными формулами второго модифицированного метода Эйлера – Коши.

Второй модифицированный метод Эйлера – Коши, так же, как и первый, является одношаговым методом со вторым порядком точности.

Оценка погрешности. Приближенная оценка погрешности модифицированных методов Эйлера осуществляется как и для простого метода Эйлера с использованием правила Рунге (см. предыдущий раздел 6.2). Так как оба модифицированных метода Эйлера имеют второй порядок точности, т. е. p = 2, то оценка погрешности (6.6) примет вид

R  » |y- y|. (6.15)

Используя правило Рунге, можно построить процедуру приближенного вычисления решения задачи Коши модифицированными методами Эйлера с заданной точностью e. Нужно, начав вычисления с некоторого значения шага h, последовательно уменьшать это значение в два раза, каждый раз вычисляя приближенное значение y, i = 0, 1, …, n. Вычисления прекращаются тогда, когда будет выполнено условие:

R  » |y- y| < e. (6.16)

Приближенным решением будут значения y, i = 0, 1, …, n.

Пример 6.2.

Применим первый модифицированный метод Эйлера для решения задачи Коши

y' (t) = y – , y(0) = 1,

рассмотренной ранее в примере 6.1.

Возьмем шаг h = 0.2. Тогда n =  = 5.

В соответствии с (6.3) получим расчетную формулу первого модифицированного метода Эйлера:

y_i+1 = y_i + h f =  y_i + 0.2 f, где

f = f(t, y) = y – ,

t = t_i +  = t_i + 0.1,

y = y_i +f(t_i, y_i) = y_i +0.1,

t₀  = 0, y₀ = 1, i = 0, 1, …, 4.

Решение представим в виде таблицы 6.3:

Таблица 6.3

Третий столбец таблицы 6.3 содержит приближенное решение y_i, i = 0, 1, …, 5.

Сравним полученное приближенное решение с точным решением (6.11), представленном в таблице 6.2. Виднм, что погрешность составляет R = | y(t_i) – y_i| = 0.0042.

Пример 6.3.

Применим второй модифицированный метод Эйлера – Коши для решения задачи Коши

y' (t) = y – , y(0) = 1,

рассмотренной ранее в примерах 6.1 и 6.2. Так же, как и ранее, зададим шаг h = 0.2. Тогда n =  = 5.

В соответствии с (6.14) получим расчетную формулу метода Эйлера – Коши:

y_i+1 = y_i + [f(t_i, y_i) + f(t_i+1, )] = y_i + 0.1[f(t_i, y_i) + f(t_i+1, )],

где

f(t_i, y_i) = y_i –

 = y_i + h f(t_i, y_i) = y_i + 0.1

t₀  = 0, y₀ = 1, i = 0, 1, …, 4.

Решение представим в виде таблицы 6.4:

Таблица 6.4

Таблица 6.4 заполняется последовательно по строкам, сначала первая строка, затем вторая и т. д. Третий столбец таблицы 6.4 содержит приближенное решение y_i, i = 0, 1, …, 5.

Сравним полученное приближенное решение с точным решением (6.11), представленном в таблице 6.2. Видим, что погрешность составляет R = | y(t_i) – y_i| = 0.0222.

Похожие работы

Лекции по вычислительной математике

Скачать

10356

0

0

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

Основные задачи вычислительной математики

Скачать

7721

1

1

... если  - предельная абсолютная погрешность приближённого числа , то  (1.2) отсюда следует, что  (1.3) Значение предельной абсолютной погрешности, обычно, подбирается интуитивно по смыслу задачи. Пример 2: Определить предельную абсолютную погрешность числа , заменяющего число , точное значение которого нам неизвестно. Так как мы знаем, что , ...

Выдающиеся личности в истории вычислительной техники. Августа Ада Лавлейс

Скачать

37333

0

0

... удивили меня…, хоть речь идёт обо мне самой. Они действительно написаны прекрасным стилем, который превосходит стиль самого очерка" /2/. 2.3. Рождение первенца и критическое перенапряжение Августа Ада Лавлейс работает с большим напряжением. В письмах к Бэббиджу она неоднократно жалуется на утомление, болезни, плохое самочувствие. Наконец, 6 августа Бэббидж отсылает Аде свои последние замечания ...

Роль женщин в развитии вычислительной техники

Скачать

54819

0

0

... в Украине, бывшем Советском Союзе и за рубежом научная школа теоретического программирования. В 2001-м году ее не стало... Но не только в научном плане велика роль женщин в развитии вычислительной техники. Со временем образуется огромное количество различных фирм по разработке и продаже программного и аппаратного обеспечения. Следовательно, разыгрываются человеческие трагедии капиталистического ...