Численное решение задачи Коши для обыкновенных дифференциальных уравнений

ТЕОРЕТИЧЕСКИЕ ОСНОВЫ

ОСНОВНЫЕ численные методы

Численное дифференцирование функций

Пусть функция f(x) дифференцируема в точке x. Тогда значение производной определяется следующими пределами [1]

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (1)

Вычисление первой производной

Для вычисления производной необходимо проведение ряда операций. Можно вычислять значения функции и проводить с ними арифметические действия. Но мы не можем вычислять пределы, поскольку это требует бесконечных затрат ресурсов (времени, памяти и т.д.). Получим приближенные формулы:

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (2)

Пусть Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru - шаг разбиения. Введем обозначение Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru и т.д. Тогда (2) можно переписать в виде

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (3)

Первое из этих трех отношений носит название правой разностной производной, второе – левой, третье – центральной разностной производной.

Вычисление второй производной

Для приближенного вычисления второй производной в качестве примера используем формулу [1]

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru , (4)

где Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru определяется по формуле (3).

Отметим, что значения правой и левой разностных производных в точке Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru одновременно являются центральными разностными производными Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru и Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru , рассчитанными соответственно в точках Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru и Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru (см. рис. 1).

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru

Рис 1. Схема численного дифференцирования

Тогда

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (5)

Численное интегрирование функций

Пусть необходимо вычислить определенный интеграл от некоторой непрерывной функции f(x)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (6)

Численное значение интеграла равно площади, заключенной между кривой y= f(x), осью x и вертикальными прямыми x=a и x=b (рис. 2).

Разобьем отрезок интегрирования на n частей. Введем в рассмотрение последовательность узловых точек xjÎ[a,b], xj=a+jh, j=0,...,n. Величина Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru называется шагом разбиения. Обозначим fj=f(xj).

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

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru

Рис. 2. Численное интегрирование

Таким образом, интеграл представляется суммой интегралов

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (7)

Метод левых прямоугольников

Вместо площади криволинейной фигуры вычисляется площадь прямоугольника, высота которого равна значению функции f(x) в крайней левой точке (рис. 3,а)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (8)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (9)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru а Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru б

Рис. 3. Методы левых (а) и правых (б) прямоугольников

Метод правых прямоугольников

Вместо площади криволинейной фигуры вычисляется площадь прямоугольника, высота которого равна значению функции f(x) в крайней правой точке (рис. 3,б)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (10)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (11)

Метод средних прямоугольников

Для вычисления интеграла можно использовать значение функции в середине отрезка (рис. 4,а)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (12)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (13)

Метод трапеций

Теперь рассмотрим численное интегрирование фигуры, полученной путем соединения отрезком прямой двух точек на графике функции. Полученная фигура является трапецией (рис. 4,б). Тогда интеграл приближенно равен ее площади

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru , (14)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru (15)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru а Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru б

Рис. 4. Методы средних прямоугольников (а) и трапеций (б)

Численное решение задачи Коши для обыкновенных дифференциальных уравнений

Постановка задачи

Рассмотрим сначала следующую задачу. Пусть задано обыкновенное дифференциальное уравнение первого порядка

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru (16)

и начальное условие

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (17)

Такая задача называется задачей Коши для обыкновенного дифференциального уравнения.

Известно, что решение обыкновенного дифференциального уравнения первого порядка для непрерывной функции f(x) графически представляется в виде однопараметрического семейства кривых (см. например, рис. 5,а). Задание начального значения позволяет выбрать из этого семейства одну кривую, которая является решением задачи Коши.

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru а Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru б

Рис. 5. Решение задачи Коши для обыкновенного дифференциального уравнения

Для численного решения задачи введем сетку xj, j=0,...,n, x0=a. Верхняя граница xn=b выбирается исходя из условий практического приложения результатов расчета. В случае равномерного разбиения с шагом h узловые точки определяются следующим образом: xj=a+jh. Обозначим yj=y(xj). Тогда

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (18)

Однако, в отличие от задачи интегрирования функций, под знак интеграла входит неизвестная зависимость y(x). Поэтому интегрирование не может быть осуществлено точно. Но существуют способы сделать это приближенно.

Метод Эйлера

Самый простой способ заключается в использовании метода левых прямоугольников

Численное решение задачи Коши для обыкновенных дифференциальных уравнений - student2.ru . (19)

Такой метод решения задачи Коши называют методом Эйлера. Результаты применения этого метода показаны на рис. 5,б. Интегрирование на каждом шаге проводится при постоянном значении f(x,y), а значение этой функции равно производной искомой функции y(x). Отсюда следует, что метод Эйлера сводится к сдвигу вдоль касательной, проведенной к графику функции y(x) из начальной точки. Далее определяется угол наклона касательной в полученной точке и проводится новый шаг.

Наши рекомендации