Тогда разностная схема четвертого порядка метода Адамса запишется в виде
.Сравнивая метод Адамса с методом Рунге — Кутта той же точности, отмечаем его экономичность, поскольку он требует вычисления лишь одного значения правой части на каждом шаге. Но метод Адамса неудобен тем, что невозможно начать счет по одному лишь известному значению
. Расчет может быть начат лишь с узла . Значения необходимые для вычисления , нужно получить каким-либо другим способом , что существенно усложняет алгоритм. Кроме того, метод Адамса не позволяет изменить шаг в процессе счета; этого недостатка лишены одношаговые методы.Метод Милна
Пусть на отрезке [a, b] требуется найти численное решение дифференциального уравнения
с начальным условием . Разобьем отрезок [a, b] на nравных частей точками , где h=(b-a)/n – шаг интегрирования. Используя начальные данные, находим каким-либо способом последовательные значения искомой функции y(x). Таким образом, становится известным . Приближения и для следующих значений последовательно находятся по формулам Милна – где .Абсолютная погрешность значения
приближенно равна .Пример. Дано дифференциальное уравнение y’=y-x, удовлетворяющие начальному условию x0=0, y(x0)=1,5. Вычислить с точность до 0,01 значение решения этого уравнения при x=1,5.
Решение. Выберем начальный шаг вычисления. Из условия h4<0,01 получим h=0,25 Составим таблицу
i | xi | yi | y’i=f(xi, yi)=yi-xi | y'i= f(xi, yi)=yi-xi | εi | ||
0 1 2 3 4 5 6 | 0 0,25 0,50 0,75 1,00 1,25 1,50 | 1,5000 1,8920 2,3243 2,8084 | 1,5000 1,6420 1,8243 2,0584 | 3,3588 3,9947 4,7402 | 2,3588 2,7447 3,2402 | 3,3590 3,9950 4,7406 | 7*10-5 10-5 1,4*10-5 |
Получаем ответ y=(1,5)=4,74.
Краевые задачи для обыкновенных дифференциальных уравнений
На практике приходится часто решать задачи, когда условия задаются при двух значениях независимой переменной (на концах рассматриваемого отрезка). Такие задачи, называемые краевыми, получаются при решении уравнений высших порядков или систем уравнений. Стандартная постановка краевой задачи для обыкновенных дифференциальных уравнений выглядит следующим образом
, а дополнительные условия ставятся более, чем в одной точке отрезка интегрирования уравнений (в этом случае порядок системы не может быть меньше второго): , , .Общая классификация методов решения краевых задач: существуют точные, приближенные и численные методы.
6. Приближенные методы решения дифференциальных уравнений с частными производными
Кроме обычных дифференциальных уравнений существуют так называемые дифференциальные уравнения с частными производными. Далее они будут рассмотрены более подробно.
Классификация дифференциальных уравнений второго порядка
Рассмотрим уравнение второго порядка
, где - функции и . Говорят, что указанное уравнение в области принадлежит гиперболическому типу, если в этой области . Если , то уравнение в области принадлежит параболическому типу. Если , то уравнение принадлежит эллиптическому типу.Уравнение
называется каноническим уравнением гиперболического типа.Уравнение
называется каноническим уравнением параболического типа.Уравнение
называется каноническим уравнением эллиптического типа.Дифференциальное уравнение
называется уравнением характеристик уравнения .Если последнее уравнение гиперболического типа, то уравнение характеристик имеет два интеграла:
т.е. существуют два семейства вещественных характеристик.С помощью замены переменных
, дифференциальное уравнение приводится к каноническому виду: . Для уравнения параболического типа оба семейства характеристик совпадают, т.е. уравнение характеристик дает лишь один интеграл .В этом случае осуществляем замену переменных
, , где — какая-нибудь функция, для которой . После замены переменных получаем уравнение . Для уравнения эллиптического типа интегралы уравнения характеристик имеют вид , где и — вещественные функции.Полагая
и , приводим уравнение к виду .Постановка краевых задач
Классическим решением краевой задачи называются всяка функция, удовлетворяющая дифференциальному уравнению в каждой точке внутри области задания этого уравнения и непрерывная в рассматриваемой области, включая границу. Соответствующую постановку краевой задачи называют классической. Существует несколько таких задач:
уравнения колебания струны и уравнения теплопроводности.