Решение дифференциальных уравнений

Теоретические сведения

Порядок выполнения работы

1. Получить у преподавателя вариант задания, включающий в себя

· дифференциальное уравнение (F(x))

· интервал (а,b)

· шаг (h)

· краевое значение функции (у 0 )

2. Написать подпрограмму для каждого метода (Эйлера, Эйлера-Коши,
Рунге-Кутта)

3. Написать подпрограмму процедуры-функции

4. Написать головной модуль

5. Отладить программу и получить результаты

6. Построить график решений дифференциального уравнения для всех 3-х методов в Excel.

Содержание отчета

1. Содержательная постановка задачи

2. Исходные данные

3. Краткое описание методов

4. Блок схема подпрограмм и блок схема головного (или управляющего) модуля

5. Листинг подпрограмм и управляющего модуля

6. Распечатка полученных результатов

7. Распечатка результатов в Excel.

Дифференциальные уравнения очень часто встречаются при построении

моделей динамики объектов исследования. Они описывают, как правило, изменения параметров объекта во времени. Результатом решения дифференциальных уравнений являются функции, а не числа, как при решении конечных уравнений, вследствие чего методы решения их более трудоемки.

Дифференциальные уравнения описывают также процессом, тепло-массообмен, изменение концентрации вещества, процессы кристаллизации сахара и многие другие. При использовании численных методов решения дифференциальных уравнений:

или y^’= f (x,y) представляется в табличном виде, т.е. получается
dx совокупность значений Y_i и X_i.

Решение носит шаговый характер, т.е. по одной или нескольким начальным точкам (х, у) за один шаг находят следующую точку и т.д. Разница между двумя соседними значениями аргумента h = x_i+1 - x_i называется шагом.

Наибольшее распространение имеют задачи Коши, в которых заданы начальные условия: при x = x₀, y(x₀) = y₀. Имея их, легко начинать процесс решения, т.е. найти при x₁ , y₂ - при х₂ и т.д.

Основная идея получения простейших вычислительных алгоритмов в одношаговых методах сводится к разложению исходного решения у(х) в ряд Тейлора.

Количество оставленных членов ряда определяет порядок и, следовательно, точность метода. По полученному разложению, зная значения у в точке разложения у_i и производную f(x_i, y_i), находят значения у через шаг h:
y_i+1 = y_i + ∆y_i.

Если в разложении удерживается большее число членов, то необходимо рассчитывать f(x_i, y_i) в несколько точках (таким способом избегают необходимости прямого вычисления высших производных, присутствующих в разложении в ряд Тейлора).

Расчётные алгоритмы многошаговых методов базируются на построении интерполяционных или аппроксимирующих функций, от которых берётся интеграл.

Численными методами решаются не только отдельные уравнения, но и системы уравнений (чаще всего первого порядка), причем большинство методов решения одного уравнения легко распространяются на решения систем.

К классу одношаговых методов относятся методы Эйлера,
Рунге – Кутта и Эйлера-Коши.

Функциональное уравнение у¢ = f(x,у), связывающее между собой независимую переменную, искомую функцию у(х) и ее производную у (х), называется дифференциальным уравнением 1-го порядка.

Решением (частным) решением уравнения на интервале (а, b) называется любая функция у = (х), которая, будучи подставлена в это уравнение вместе со своей производной ¢ (x) обращает его в тождество относительно xÎ (а,b). Уравнение Ф. (х,y) = 0, определяющее это решение как неявную функцию, называется интегралом дифференциального уравнения. На плоскости с фиксированной декартовой прямоугольной системой координат уравнение Ф (х,y) =0 определяет некоторую кривую, которая называется интегральной кривой дифференциального уравнения.

Если в дифференциальном уравнении у¢ = f(x,у) функция f(x,у) непрерывна в некоторой области D, плоскости Оху и имеет в этой области ограниченную частную производную (x,y), то для любой точки (x₀,y₀) Î D, в некотором интервале х₀ — h £ х £ х₀ + h, существует и притом единственное решение у (х) этого уравнения, удовлетворяющее начальному условию

у (х_о) - у_о.

Это утверждение известно как теорема Коши о существовании и единственности решения дифференциального уравнения с заданным начальным условием.

Для задач подобного типа, выделенных в целый класс задач Коши, помимо аналитических методов решения разработаны методы численного решения.

Дата добавления: 2015-05-29