Курсовая работа: Решение дифференциальных уравнений. Обзор

начальные условия;

набор точек в которых нужно найти решение;

само дифференциальное уравнение, записанное в некотором специальном виде, который будет описан ниже.

Один из наиболее эффективных алгоритмов интегрирования ОДУ основан на численном методе Рунге-Кутты четвертого порядка. Функция, реализующая этот метод, имеет вид rkfixed (y,x₁ ,x₂ , npoints,D)

Здесь:

y-вектор начальных условий размерности n, где n- порядок дифференциального уравнения или число уравнений в системе (если решается система уравнений);

x₁ , x₂ – граничные точки интервала, на котором ищется решение дифференциального уравнения. Начальные условия ,заданные в векторе y,- это значение решения в точке x₁ ;

npoints- число точек (не считая начальной точки), в которых ищется приближенное решение. При помощи этого аргумента определяется число строк (1+npoints) в матрице, возвращаемой функцией rkfixed;

D(x,y) – функция,возвращающая значение в виде вектора n элементов, содержащих первые производные неизвестных функций.

2.1 Метод Эйлера

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

Численное решение на отрезке [a, b] задачи Коши

y' =f(x, y), y(a) =y₀

состоит в построении таблицы приближенных значений

y₀ ,y₁ , ...,y_i , ...y_N

решенияy(x)в узлах сетки

a=x₀ <x₁ < ... <x_i < ...<x_N =b, y(x_i )@y_i .

Еслиx_i =a+ i h, h=(b-a)/ N,то сетка называетсяравномерной.

Численный метод решения задачи Коши называетсяодношаговым , если для вычисления решения в точкеx₀ +hиспользуется информация о решении только в точкеx₀ .

Простейший одношаговый метод численного решения задачи Коши -метод Эйлера . В методе Эйлера величиныy_i вычисляются по формуле

y_i +1 =y_i +hf(x_i ,y_i ), i= 0, 1

Найдем методом Эйлера на [0, 1] с шагом h=0.2 приближенное решение задачи Коши

Для того чтобы изменить стиль изображения, щелкните дважды по полю графиков и установите соответствующие параметры

Определим правую часть уравнения

Расчетные формулы метода Эйлера для решения этой задачи имеют вид

x0=0, y0= 1, xi+1 = xi + 0.2, yi+1 = yi + 0.2(sinxi - cosyi), i =0, 1, ..., 4.

Изобразим приближенное решение графически.

y' = sin x – cos y, y(0)=1.

Определим диапазон изменения номера точки i=0,1, ..., 4