Тема: Численные методы решения дифференциальных уравнений.
Цель: Научиться решать дифференциальные уравнения методами Пикара, Эйлера и Рунге-Кутта.
Источники:
1. Лапчик М.П. Численные методы, М., 2004
2. Бажин И.И. Информационные системы менеджмента, М: ГУ-ВШЭ, 2000
Оборудование: ПК, ПО ЭТ Excel
Ход и содержание работы:
Часть 1:
Задание 1. Найдите методом последовательных приближений решение дифференциального уравнения , удовлетворяющее начальному условию у(0) = 2.
Вычислить y1, y2, y3 и изобразить на графике точное решение и найденную последовательность y1, y2, y3 в интервале [0;0.3]
Задание 2. Найти приближенное решение дифференциального уравнения у' = cos у + 3х методом Пикара с начальным условием у(0) = 1,3 (расчет вести до третьего приближения).
Задание 3. Решить методом Эйлера дифференциальное уравнение у' = cos у + Зх с начальным значением у(0) = 1,3 на отрезке [0; 1], приняв шаг h = 0,2.
Задание 4. Методом Эйлера найти приближенное решение дифференциального уравнения , удовлетворяющее начальному условию у(0) = 1 на отрезке [0, 2] двумя способами: а) с шагом h = 0,2; б) с с шагом h = 0,1. Сопоставить точность полученных значений функции.
|
|
Задание 5. Составить программу решения дифференциальных уравнений и у' = cos у + Зх методом Эйлера по приведенному на рисунке алфавиту и с заданием пользователя начального условия и интервала.
Часть 2:
Задание 1. Решить дифференциальное уравнение у' = y(1-х) на отрезке [0; 0,5] с начальным условием у(0) = 1 и шагом h = 0,05 методом Эйлера и Рунге-Кутта 4-го порядка.
Задание 2. Проведите сравнение методов Эйлера и Рунге-Кутта по результатам работы соответствующих программ (алгоритмы метода Рунге-Кутта и таблица сравнения представлена на рисунке)
Задание3. Пользуясь программой для метода Рунге-Кутта, решить дифференциальное уравнение у’ = у2 + Зх с начальным условием у(1) = 3 на отрезке [1; 2] сначала с шагом h = 0,1, а затем с шагом h = 0,05. Сопоставить и прокомментировать полученные решения.
Часть 3:
Задание 1. Решить задачу Коши для дифференциального уравнения своего варианта у = f(x,y) на отрезке [a; b] при заданном начальном условии и шаге интегрирования методами Эйлера и Рунге —Кутта с помощью Excel b программы для компьютера с шагом h.
Контрольные вопросы
1. В какой форме получается приближенное решение дифференциального уравнения по методу Пикара?
2. В какой форме получается приближенное решение дифференциального уравнения по методу Эйлера?
3. В чем различие одношаговых методов Эйлера и Рунге —Кутта? Как это различие можно охарактеризовать с графической точки зрения?