Вычислительная математика. Практическое задание 6 Вариант 8. Практическое задание 6. Тема Численные методы решения дифференциальных уравнений и систем Лекция Численные методы решения обыкновенных дифференциальных уравнений и систем первого и второго порядка
![]()
|
Практическое задание 6Тема 6. Численные методы решения дифференциальных уравнений и систем Лекция 6.1. Численные методы решения обыкновенных дифференциальных уравнений и систем первого и второго порядка Номера вариантов задач определяются с помощью табл. 1 по первой букве фамилии студента. Таблица 1
Вариант 8 Задание 6.1. Решить задачу Коши для дифференциального уравнения первого порядка на заданном отрезке: методом Эйлера; модифицированным методом Эйлера; методом Рунге – Кутты.
Решение: 1) Метод Эйлера Метод Эйлера численного решения задачи Коши ![]() ![]() ![]() Расчетная формула метода Эйлера: ![]() Шаг 1: ![]() ![]() Дальнейшие рассчеты представим в таблице:
Ответ: ![]() ![]() 2) Модифицированный метод Эйлера В модифицированном методе Эйлера численного решения задачи Коши ![]() ![]() ![]() После чего вычисляется приближенное решение: ![]() Выполним вычисления по формулам: “грубое” приближение: ![]() приближенное решение: ![]() ![]() Шаг 1: “грубое” приближение: ![]() приближенное решение задачи Коши: ![]() ![]() Дальнейшие рассчеты представим в таблице:
Ответ: ![]() ![]() ![]() 3) метод Рунге-Кутты Решение дифференциального уравнения вида ![]() ![]() Где коэффициенты вычисляются по формулам: ![]() ![]() ![]() ![]() Выполним вычисления по формулам: ![]() ![]() ![]() ![]() Шаг 1: Вычисляем значения коэффициентов: ![]() ![]() ![]() ![]() ![]() ![]() ![]() Дальнейшие рассчеты представим в таблице:
Ответ: ![]() ![]() Задание 6.2 1. Найти решение дифференциального уравнения второго порядка с заданными начальными условиями методом неопределенных коэффициентов. 2. Найти первые пять членов решения дифференциального уравнения второго порядка с заданными начальными условиями.
Решение: Используем метод конечных разностей. Произведем дискретизацию области непрерывного аргумента сеткой с шагом ![]() ![]() ![]() Вычислим решение краевой задачи не в произвольных точках отрезка, а только в узлах сетки. Производные ![]() ![]() ![]() Заменяем в каждом из внутренних узлов ![]() ![]() Потребуем теперь, чтобы значения искомой сеточной функции удовлетворяли во всех внутренних узлах сетки уравнениям, полученных заменой знака приближенного равенства на знак строгого равенства: ![]() Или: ![]() В граничных узлах сеточная функция должна удовлетворять равенствам: ![]() ![]() Граничные условия аппроксимируем правой и левой разностной производными: ![]() ![]() Тогда имеем ![]() ![]() ![]() Приводим уравнения для внутренних узлов ![]() ![]() С учетом того, что ![]() ![]() Выражая последовательно: - из первого уравнения: ![]() - из второго уравнения: ![]() - из четвертого уравнения: ![]() Подставляя в третье уравнение, находим: ![]() ![]() Остальные значения: ![]() ![]() ![]() Ответ: ![]() ![]() ![]() ![]() |