Лабораторная работа №5 Численное интегрирование ОДУ. Задача Коши
Технические требования:
Среда разработки для написания кода на Python Компилятор Python 3 (+ стандартный набор библиотек)
Требования к знаниям студентов:
Основы синтаксиса Python 3:
-умение работать с переменными и знание типов данных
-умение работать со списками
-умение считывать данные из файла
-умение выводить данные в консоль/файл
-знание условных конструкций if-elif-else
-умение писать циклы for/while
-умение создавать и вызывать функции, возращать значения
-умение генерировать случайные числа
Цель работы:
Знакомство с численными методами для решения задач:
-численное решение дифференциального уравнения
-численное решение задачи Коши
Теоретическое введение:
Блок 1 “Численное интегрирование ОДУ. Задача Коши”
Дифференциальные уравнения — это фундаментальная область математики, связанная с моделированием физических процессов, инженерных систем и многих других явлений. Они позволяют описать изменение величин в зависимости от различных факторов, таких как время, пространство или другие переменные. В отличие от аналитического подхода, который требует нахождения точного решения уравнений, численные методы используются тогда, когда найти точное решение сложно или невозможно. Рассмотрим несколько популярных алгоритмов численного решения обыкновенных дифференциальных уравнений (ОДУ): метод Эйлера, метод Хойна и метод Рунге-Кутта 4 порядка.
Обыкновенное дифференциальное уравнение первого порядка можно записать следующим образом:
где y — искомая функция, зависящая от переменной t, а f(t,y) — известная функция. Задача состоит в том, чтобы найти такую функцию y(t), которая удовлетворяет этому уравнению и начальным условиям.
Метод Эйлера - один из простейших методов численного интегрирования. Его идея заключается в приближенном вычислении значения функции y(t) через малые шаги по времени Δt. Формула метода Эйлера выглядит так:
где h — шаг интегрирования,yn — значение функции в точке tn, f(tn,yn) — производная функции в этой же точке.
Метод Эйлера является явным методом первого порядка точности, однако он имеет большую погрешность и нестабильность при больших значениях шага h.
Метод Хойна улучшает стандартную схему Эйлера путем усреднения двух значений производной на каждом шаге. Формула выглядит так:
где k1=f(tn,yn), k2=f(tn+h,yn+h k1).
Метод Хойна имеет второй порядок точности и обеспечивает лучшее поведение по сравнению с обычным методом Эйлера.
Метод Рунге-Кутта четвертого порядка является одним из наиболее широко используемых методов численного интегрирования благодаря своей высокой точности и устойчивости. Этот метод базируется на оценке нескольких промежуточных значений производной на каждом шаге:
Предложенный метод Рунге-Кутта обладает высоким порядком точности (четвертый порядок) и часто используется там, где важны высокая точность и стабильность расчета.
Практическая часть
Задание 1.1 (основная часть) “Численное решение задачи Коши”
Задача Коши представляет собой одну из ключевых постановок задачи для обыкновенного дифференциального уравнения (или системы). Она формулируется следующим образом:
Пусть дано дифференциальное уравнение первого порядка вида:
и задано начальное условие:
Тогда задача Коши состоит в нахождении такой функции y(x), которая одновременно удовлетворяет данному дифференциальному уравнению и начальному условию.
Реализуйте решение задачи Коши для произвольно выбранных функции f(t, y) и начального значения y(0) всеми 3 изученными методами. Сравните их эффективности (критерий эффективности предлагается выбрать и обосновать самостоятельно).
Заключение:
В качестве результата необходимо представить отчет, включающий в себя листинг программы (текст исходного кода), полученные результаты, а также их анализ.
При сдаче работы студент должен знать изученные численные методы, а также разбираться в коде написанной программы.