- •Глава 7. Дифференциальные уравнения
- •7.1. Общие понятия
- •7.2. Дифференциальные уравнения первого порядка
- •7.3. Задача Коши. Теорема о существовании и единственности решения дифференциального уравнения первого порядка
- •7.4. Дифференциальные уравнения с разделяющимися переменными
- •7. 5. Однородные дифференциальные уравнения (дифференциальные уравнения с однородными функциями)
- •7.6. Линейные дифференциальные уравнения первого порядка
- •7.7. Уравнение Бернулли
- •7.8. Дифференциальные уравнения первого порядка в полных дифференциалах
- •7.9. Дифференциальные уравнения высших порядков
- •7. 10. Дифференциальное уравнение вида
- •7.11. Дифференциальные уравнения второго порядка, приводимые к дифференциальным уравнениям первого порядка
- •7.12. Линейные дифференциальные уравнения n-ого порядка. Свойства их решений
- •7.13. Линейная зависимость функций. Определитель Вронского
- •7.14. Структура общего решения линейного неоднородного дифференциального уравнения n-ого порядка
- •7.15. Комплексные числа и действия над ними
- •Действия над комплексными числами
- •7.16. Показательная функция с комплексным показателем
- •7.17. Показательная форма комплексного числа. Формула Эйлера
- •7.18. Общее решение линейного однородного дифференциального уравнения n-ого порядка с постоянными коэффициентами.
- •7.19. Частное решение линейного неоднородного дифференциального уравнения n-ого порядка с постоянными коэффициентами
- •7.20. Метод Эйлера численного интегрирования дифференциальных уравнений
7.20. Метод Эйлера численного интегрирования дифференциальных уравнений
Метод Эйлера является самым простым из множества методов численного интегрирования дифференциальных уравнений.
Пусть для дифференциального уравнения на отрезке требуется решить задачу Коши, т. е. необходимо найти решение дифференциального уравнения , удовлетворяющее начальным условиям , где некоторое заданное значение, .
Разобьем отрезок на n элементарных отрезков длиной с помощью точек . На каждом элементарном отрезке заменим интегральную кривую касательной к этой кривой (рис. 85) и найдем приближенные значения решения дифференциального уравнения в этих точках .
Рис. 85
На первом элементарном отрезке дано значение искомой функции в граничной точке . Найдем приближенное значение функции в правой граничной точке рассматриваемого отрезка.
Получаем
.
Данное значение является исходным для нахождения значения искомой функции на отрезке , где , .
Получаем
.
Аналогично вычисляются приближенные значения решения на следующих элементарных отрезках по формулам
.
Вычисления продолжаются до тех пор, пока значение х достигнет конечной точки отрезка .
Общая схема расчета состоит в том, что сначала проводят расчет при некотором произвольно выбранном значении n, получают значение искомой функции (частного решения) в конечной точке отрезка интегрирования . Затем увеличивают число элементарных отрезков (обычно в два раза, 2n). Снова проводят расчет, находят . И сравнивают полученные результаты
,
где заданная точность расчета.
Если два последовательных приближения отличаются друг от друга менее, чем на заданное значение , то последнее значение функции (решения) принимается за окончательное. Иначе, число элементарных отрезков увеличивается и расчет продолжается.
Пример 7.28. Дифференциальное уравнение при проинтегрировать на отрезке .
Пусть , . Результаты расчета приведены в таблице.
i |
|||||
0 |
0 |
1 |
1 |
0,1 |
1,1 |
1 |
0,1 |
1,1 |
1,2 |
0,12 |
1,22 |
2 |
0,2 |
1,22 |
1,42 |
0,142 |
1,362 |
3 |
0,3 |
1,362 |
1,662 |
0,1662 |
1,5282 |
4 |
0,4 |
1,5282 |
1,9282 |
0,19282 |
1,72102 |
5 |
0,5 |
1,72102 |
|
|
|
Следовательно .