Численные методы решения дифференциальных уравнений

Сентябрь 3, 2022

Главная
Алгебра
Численные методы решения дифференциальных уравнений

Содержание

2. Рассмотрим дифференциальное уравнение первого порядка: соотношение часто удается записать в виде: Последнее уравнение называется дифференциальным уравнением,
3. Метод Эйлера Дано дифференциальное уравнение y’=f(x,y), удовлетворяющее начальному условию y(x0)=y0 . Требуется найти решение на отрезке
4. Значение функции y1 в точке x1 можно определить как точку пересечения касательной проведенной к функции y=y(x)
5. Пример y1 = -2+0.5*(-(-2/1)) = -1 x1 = 1+0.5 =1.5 y2 = -1+0.5*(-(-1/1.5)) = -0.667 x2
6. Модифицированный метод Эйлера Графическая интерпретация. Определяем точку и вычисляем значение функции в этой точке Значение функции
7. произвольную точку определим где i = 1, 2, 3, …, n Пример y1/2 = -2+0.25*(-(-2/1)) =
9. Скачать презентацию

Слайд 2

Рассмотрим дифференциальное уравнение первого порядка:
соотношение часто удается записать в виде:
Последнее уравнение

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

Общим решением уравнения будет являться семейство функций y=y(x,c1) различающихся значение постоянной c1. Задаем одно начальное условие y(x0)=y0, которое определяет значение c1и конкретное частное решение – задача Коши.

Для простейшего дифференциального уравнения y’=3x2. Общее решение имеет вид y=x3+c, а подставив в общее решение начальное условие x0=1, y0=2 вычислим с=1 и определим частное решение как: y=x3+1

Слайд 3

Метод Эйлера
Дано дифференциальное уравнение y’=f(x,y), удовлетворяющее начальному условию y(x0)=y0 . Требуется

найти решение на отрезке [a,b]. Разобьем отрезок интегрирования на n равных частей: x0=a, x1= a+h, x2=x1+h,…,xi=xi–1+h,…,xn=b, тогда величина шага интегрирования будет равна:

Слайд 4

Значение функции y1 в точке x1 можно определить как точку пересечения

касательной проведенной к функции y=y(x) в точке (x0,y0) с вертикальной прямой проходящей через точку x1.

Тангенс угла наклона касательной есть значение производной в точке (x0,y0) и задается правой частью дифференциального уравнения, т.е. tg( β)=f(x0,y0). С другой стороны из геометрического представления метода можно записать:

Следовательно

Откуда

и т.д.

Решение будет заключаться в последовательном применении формул:

где i = 1, 2, 3, …, n

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

Слайд 5

Пример
y1 = -2+0.5(-(-2/1)) = -1 x1 = 1+0.5 =1.5
y2 = -1+0.5(-(-1/1.5))

= -0.667 x2 = 1.5+0.5 =2
и т.д.

Слайд 6

Модифицированный метод Эйлера
Графическая интерпретация.
Определяем точку
и вычисляем значение функции в

этой точке

Значение функции y1 в точке x1 определяем, как точку пересечения касательной, вычисленной в точке (x1/2,y1/2) и проведенной к функции y=y(x) в точке (x0,y0) , с вертикальной прямой проходящей через точку x1.

Слайд 7

произвольную точку определим
где i = 1, 2, 3, …, n

Пример

y1/2 = -2+0.25*(-(-2/1)) = -1.5; x1/2 = 1+0.25 = 1.25
y1 = -2+0.5*(-(-1.5/1.25)) = -1.4; x1 = 1+0.5 = 1.5
y3/2 = -1.4 +0.25*(-(-1.4/1.5)) = -1.1667; x3/2 = 1.5+0.25 = 1.75
y2 = -1.4 +0.5*(-(-1.1667/1.75)) = -1.0667; y2 = 1.5+0.5 = 2
и т.д.