Вычислительная математика (стр. 17 из 20)

R »

|y

- y

| < e. (6.8)

Для метода Эйлера условие (6.8) примет вид

R » |y

- y

| < e (6.9)

Приближенным решением будут значения y

, i = 0, 1, …, n.

Пример 6.1.

Найдем решение на отрезке [0, 1] следующей задачи Коши:

y' (t) = y –

, (6.10)

y(0) = 1.

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

= 5.

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

y_i+1 = y_i + 0.2

, y₀ = 1, i = 0, 1, 2, 3, 4, 5.

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

Таблица 6.1

i	0	1	2	3	4	5
ti	0	0.2	0.4	0.6	0.8	1.0
yi	1.0000	1.2000	1.3733	1.5294	1. 6786	1.8237

Уравнение (6.10) есть уравнение Бернулли. Его решение можно найти в явном виде:

y =

. (6.11)

Для сравнения точного и приближенного решений представим точное решение (6.11) в виде таблицы 6.2:

Таблица 6.2

i	0	1	2	3	4	5
ti	0	0.2	0.4	0.6	0.8	1.0
y(ti)	1.0000	1.1832	1.3416	1.4832	1. 6124	1.7320

Из таблицы видно, что погрешность составляет R =

| y(t_i) – y_i| = 0.0917.

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

i	ti	yi	f(ti, yi)	t	y	h f
0 1 2 3 4 5	0 0.2 0.4 0.6 0.8 1.0	1 1.1836 1.3426 1.4850 1.6152 1.7362	0.1 0.0850 0.0747 0.0677 0.0625	0.1 0.3 0.5 0.7 0.9	1.1 1.2682 1.4173 1.5527 1.6777	0.1836 0.1590 0.1424 0.1302 0.1210

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