ВУЗ:
Составители:
Рубрика:
74
Первая формула называется явной разностной схемой , поскольку дает
возможность вычислить значение
1+i
y в следующей точке сетки по явной
формуле
(
)
iiii
yxhfyy ,
1
+
=
+
.
Две другие формулы представляют неявные разностные схемы.
Метод Эйлера
Основу метода Эйлера составляет явная разностная схема на
равномерной разностной сетке. Интервал интегрирования [a,b] обыкно -
венного дифференциального уравнения первого порядка разбивается на n
частей узловыми точками
ihxx
i
+
=
0
, ni ,...,2,1
=
, и значения
i
y
в этих точках
вычисляется по формуле
(
)
iiii
yxhfyy ,
1
+
=
+
, 1,...,2,1,0
−
=
ni .
Метод Эйлера относится к группе одношаговых методов, в которых
для расчета точек
(
)
11
,
++ ii
yx необходима информация о предыдущей
вычисленной точке
(
)
ii
yx , .
На практике для оценки погрешности расчета функции
(
)
xy
используется правило Рунге . Сначала проводится вычисление
h
i
y с шагом
(
)
nabh
−
=
. Затем число интервалов разбиения удваивается
(
)
nn 2
←
, и
проводятся вычисления
2
h
i
y . За оценку погрешности вычислений на сетке с
шагом 2h принимается величина
2
2
11
max
h
i
h
i
ni
yy −
−≤≤
.
Пример 8.4. Найти методом Эйлера численное решение задачи Коши
1−
=
′
x
ey
при 1)1(
=
y в интервале [1;3] с шагом 0,5. При каком числе
разбиений будет обеспечена погрешность 1,0
=
ε
?
Решение . Проведем решение данной задачи в среде MathCAD.
74 Первая формула называется явной разностной схемой, поскольку дает возможность вычислить значение yi +1 в следующей точке сетки по явной формуле yi +1 = yi +hf (xi , yi ) . Две другие формулы представляют неявные разностные схемы. Метод Эйлера Основу метода Эйлера составляет явная разностная схема на равномерной разностной сетке. Интервал интегрирования [a,b] обыкно- венного дифференциального уравнения первого порядка разбивается на n частей узловыми точками xi =x0 +ih , i =1,2,..., n , и значения yi в этих точках вычисляется по формуле yi +1 = yi +hf (xi , yi ) , i =0,1,2,..., n −1 . Метод Эйлера относится к группе одношаговых методов, в которых для расчета точек (xi +1 , yi +1 ) необходима информация о предыдущей вычисленной точке (xi , yi ) . На практике для оценки погрешности расчета функции y(x ) используется правило Рунге. Сначала проводится вычисление y ih с шагом h =(b −a ) n . Затем число интервалов разбиения удваивается (n ← 2n ) , и h проводятся вычисления y . За оценку погрешности вычислений на сетке с i 2 h шагом h 2 принимается величина 1max ≤i ≤n −1 y −y . h i 2 2i Пример 8.4. Найти методом Эйлера численное решение задачи Коши y ′ =e x −1 при y(1) =1 в интервале [1;3] с шагом 0,5. При каком числе разбиений будет обеспечена погрешность ε =0,1 ? Решение. Проведем решение данной задачи в среде MathCAD.
Страницы
- « первая
- ‹ предыдущая
- …
- 72
- 73
- 74
- 75
- 76
- …
- следующая ›
- последняя »