Занятие № 16. Численные методы решения задачи Коши для обыкновенных дифференциальных уравнений. Классификация методов. Метод Пикара

Будем рассматривать обыкновенное дифференциальное уравнение (ОДУ) первого порядка

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

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

y(х0) = у0, (2)

где f(x) — некоторая заданная, в общем случае, нелинейная функция двух переменных. Будем считать, что для данной задачи (1)-(2), называемой начальной задачей или задачей Коши, выполняются требования, обеспечивающие существование и единственность на отрезке [х0,b] ее решения у=у(х).

Несмотря на внешнюю простоту уравнения (1), решить его аналитически, т.е. найти общее решение у=у(х, С) с тем, чтобы затем выделить из него интегральную кривую у=у(х), проходящую через заданную точку (х00), удается лишь для некоторых специальных типов таких уравнений. Поэтому, как и в родственной для (1)-(2) задаче вычисления интегралов, приходится делать ставку на приближенные способы решения начальных задач для ОДУ, которые можно разделить на три группы:

1)приближенно-аналитические методы;

2)графические или машинно-графические методы;

3)численные методы.

К методам первой группы относят такие, которые позволя­ют находить приближение решения у(х) сразу в виде некоторой «хорошей» функции φ(х).Например, широко известен метод степенных рядов, в одну из реализаций которого заложено представление искомой функции у(х) отрезком ряда Тейлора, где тейлоровские коэффициенты, содержащие производные высших порядков, находят последовательным дифференцирова­нием самого уравнения (1). Другим представи­телем этой группы методов является метод последовательных приближений, суть которого приведена чуть ниже.

Название графические методы говорит о приближенном представлении искомого решения у(х) на промежутке [x0,b] в виде графика, который можно строить по тем или иным прави­лам, связанным с графическим толкованием данной задачи. Фи­зическая или, возможно, точнее будет сказать, электротехниче­ская интерпретация начальных задач для определенных видов уравнений лежит в основе машинно-графических методов при­ближенного решения. Реализуя на физико-техническом уровне заданные электрические процессы, на экране осциллографа на­блюдают поведение решений дифференциальных уравнений, описывающих эти процессы. Изменение параметров уравнения приводит к адекватному изменению поведения решений, что по­ложено в основу специализированных аналоговых вычисли­тельных машин (АВМ).

Наконец, наиболее значимыми в настоящее время, характе­ризуемое бурным развитием и проникновением во все сферы че­ловеческой деятельности цифровой вычислительной техники, являются численные методы решения дифференциальных уравнений, предполагающие получение числовой таблицы приближенных значений yi искомого решения у(х) на некото­рой сетке Занятие № 16. Численные методы решения задачи Коши для обыкновенных дифференциальных уравнений. Классификация методов. Метод Пикара - student2.ru значений аргумента х. Этим способам и будет посвящено дальнейшее изложение. Что делать с получае­мыми численными значениями решения, зависит от прикладной постановки задачи. Если речь идет о нахождении только значе­ния у(b), тогда точка b включается как конечная в систему рас­четных точек хi, и все приближенные значения yi ≈y(xi), кроме последнего, участвуют лишь как промежуточные, т.е. не требуют ни запоминания, ни обработки. Если же нужно иметь прибли­женное решение у(х) в любой точке х, то для этого к получае­мой числовой таблице значений yi можно применить какой-либо из способов аппроксимации табличных функций, рассмотренных ранее, например, интерполяцию или сплайн-интерполя­цию. Возможны и другие использования численных дан­ных о решении.

Коснемся одного приближенно-аналитического способа решения начальной задачи (1)-(2), в котором искомое ре­шение у=у(х) в некоторой правой окрестности точки х0 явля­ется пределом последовательности получаемых определенным образом функций уп(х).

Проинтегрируем левую и правую части уравнения (1) в границах от х0 до х:

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

Отсюда, с учетом того, что одной из первообразных для у'(х) служит у(х), получаем

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

или, с использованием начального условия (2),

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

Таким образом, данное дифференциальное уравнение (1) с на­чальным условием (2) преобразовалось в интегральное урав­нение (неизвестная функция здесь входит под знак интеграла).

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

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

достаточно обстоятельно рассматривавшийся приме­нительно к системам линейных и нелинейных алгебраических и трансцендентных уравнений. Беря в качестве начальной функции y0(х) заданную в (2) по­стоянную y0, по формуле (4) при п=0 находим первое при­ближение

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

Его подстановка в (4) при п=1 дает второе приближение

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

и т.д. Таким образом, этот приближенно-аналитический метод, называемый методом последовательных приближений или методом Пикара определяется формулой

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

где n=0,1, 2,... и у0(х)=y0.

Отметим две характеристики метода последовательных приближений Пикара, которые можно отнести к негативным. Во-первых, в силу известных проблем с эффективным нахождением первообразных, в чистом виде метод (5) редко реализуем. Во-вторых, как видно из вышеприведенного утверждения, этот ме­тод следует считать локальным, пригодным для приближения решения в малой правой окрестности начальной точки. Большее значение метод Пикара имеет для доказательства существования и единственности решения задачи Коши, нежели для его практи­ческого нахождения.

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