О приближенном решении дифференциальных уравнений, общее решение которых зависит от константы алгебраически Текст научной статьи по специальности «Математика»

Математика и механика

УДК 517.9

О приближенном решении дифференциальных уравнений, общее решение которых зависит от константы алгебраически

М. Д. Малых

Факультет наук о материалах Московский государственный университет им. М.В. Ломоносова Ленинские Горы, Корпус «Б», Москва, Россия, 119991 Кафедра прикладной информатики и теории вероятностей Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198

Методы аналитической теории обыкновенных дифференциальных уравнений (ОДУ) основаны на анализе особенностей, но самый популярный метод для численного решения, а именно метод конечных разностей, не работает вблизи особенностей. Однако Пенлеве дал алгебраический метод для решения в конечном виде дифференциальных уравнений, общие решения которых зависят от константы интегрирования алгебраически. Этот подход, который был представлен как своеобразная теория Галуа, напротив, может быть хорошо увязан с методом конечных разностей.

Как известно, обыкновенное дифференциальное уравнение вида у' = f (х,у), обладающее этом свойством, может быть преобразовано алгебраически заменой к уравнению Риккати. Схема Эйлера уп+1 = уп + /(хп,уп)Ах всегда задаёт (1, ^-соответствие между соседними слоями. В то же время точное решение уравнения Риккати задаёт (1,1)-соответствие между любыми слоями и поэтому мы можем написать схему, задающую (1, 1)-соответствие между соседними слоями. В этом случае ангармоническое отношение четырёх точек не меняется от слоя до слоя не только для точного, но также и для приблизительного решения. Таким образом, если у точного решения имеется полюс, то приближенное решение проходит через бесконечность без накопления ошибки. В представленной статье это свойство (1,1)-схем будет проиллюстрировано двумя примерами: с и без решения в элементарных функциях. Таким образом, причина разрушения приближенного решения около полюса спрятана в саму схему Эйлера. В более общем случае, когда точное решение обыкновенного дифференциального уравнения зависит от постоянной интегрирования алгебраически, мы можем написать схему, которая задаёт (I, ^-корреспонденция между соседними слоями. Приближенное решение, найденное на этом пути, проходит через подвижные алгебраические особенности без накопления ошибки.

Ключевые слова: уравнение Риккати, метод конечных разностей, алгебраические соответствия, sage, sagemath.

1. Введение

Открытию трансцендентных уравнений Пенлеве предшествовали работы Пенлеве, в которых развивался алгебраический подход к интегрированию дифференциальных уравнений. Отправной точкой для Пенлеве стало следующее наблюдение: если обыкновенное дифференциальное уравнение = /(х,у) с рациональ-

-х

ной по у правой частью допускает общее решение, зависящее алгебраически от константы интегрирования, то это уравнение сводится алгебраической заменой к уравнению Риккати, см. [1, С. 46-47.]. В XX веке эти работы Пенлеве воспринимали лишь как способ обосновать неэлементарность трансцендентных уравнений Пенлеве, к тому же не вполне строгий. На самом же деле это — версия теории Га-луа для нелинейных дифференциальных уравнений, своеобразие которой состоит

Статья поступила в редакцию 10 июня 2015 г.

6

Вестник РУДН. Серия Математика. Информатика. Физика. №3,2015. С. 5-9

в том, что в ней изначально не фиксируется список допустимых трансцендентных операций [2,3].

Для вычисления этих трансцендентных функций во времена Пенлеве использовали метод разложения в степенной ряд, однако алгебраические конструкции Пенлеве не связаны с анализом особенностей и поэтому прекрасно сочетаются с методом конечных разностей. На рубеже XIX-XX веков, когда разностные методы решения дифференциальных уравнений только начали теснить применение разложений в степенной ряд, методы Эйлера рассматривали исключительно как способ отыскания частных решений (см., например, [4, §89]), на самом же деле ничто не мешает рассматривать разностную схему вида F(yn,yn+i,xn) = 0 как уравнение, задающее алгебраическое соответствие между п-м и (п + 1)-м слоями.

Замечание 1. Мы рассматриваем уп как координаты точки на проективный прямой, которые называем слоем. Уравнение, описывающее шаг разностной схемы, в таком случае задаёт алгебраическое соответствие между двумя слоями. Употребление термина «слой» в данном случае оправдано не только обычной практикой теории разностных схем для уравнений в 2-мерных областях. Кадзуо Окамото [5] рассматривал расслоения, на которых дифференциальные уравнения задают однозначные слоения. Для уравнения Риккати слоем такого расслоения будет как раз проективная прямая, именуемая там пространством начальных данных.

Если дифференциальное уравнение допускает общее решение, которое зависит алгебраически от константы интегрирования, то оно допускает рациональный интеграл г(х, у) = const, коэффициенты которого являются какими-то, вообще говоря, трансцендентными функциями х. Если обозначить степень г относительно у как I, то уравнение г(хп,уп) = г(хт,ут) позволяет найти I значений ут по заданному уп и I значений уп по заданному ут, то есть точное решение начальной задачи задаёт (I, ¿)-значное алгебраическое соответствие между слоями. Чем примечательны разностные схемы, сохраняющие этой свойство?

2. Схемы типа (1,1)

Уравнение Риккати у' = р(х)у2 + q(x)y + г(х) задаёт бирациональное соответствие между слоями. Схема Эйлера уп+\ = Уп + (РпУп + ЯпУп + rn)dx задаёт (1, 2)-соответствие между соседними слоями, и следовательно, не сохраняет это свойство. Однако минимальные поправки в схеме позволят сохранить тип зависимости от уо,:

Уп+1 - Уп = (Рп + ЯпУп + rnynyn+i)Ax. (1)

задаёт (1,1)-соответствие между соседними, а следовательно, и любыми слоями. С точки же зрения аппроксимации уравнения Риккати эта схема не лучше и не хуже схемы Эйлера [6].

Поскольку все известные свойства решения уравнения Риккати выведены из того, что его общее решение задаёт бирациональное отображение между начальными и конечными данными, все они переносятся и на решения указанной разностной схемы. Напр., ангармоническое отношение четырёх частных решений уравнения (1) не меняется от слоя к слою:

(Уп,У'п ,Уп,Уп ) = (Ут, Ут, Ут, Ут), (2)

коль скоро оно вообще сохраняется при дробно-линейных подстановках. Особо подчеркнём, что ангармоническое отношение четырёх частных приближенных решений сохраняется точно, в то время как в схеме Эйлера эта величина сохраняется приближённо, с точностью до неизвестных членов порядка Ах, причём коэффициенты при Ах зависят от уп, и с ростом уп ошибка становится все заметнее.

В частности, если точное решение имеет полюс, то приближенное решение, найденное по схеме (1), уходит на бесконечность и возвращается без заметного накопления ошибки. В самом деле, пусть частное решение уравнения Риккати, равное уо при хо, имеет подвижный полюс при х = а. Любые три другие частные решения у', у", у'" не имеют полюса при х = а, и любая схема 1-го порядка позволяет вычислить их значения до и после х = а без заметного накопления ошибки. Тогда искомое точное решение можно вычислить с хорошей точностью из условия (у,у',у",ут) = (у0,Уо,Уо,Уо'); если при х = хп заменить здесь точные частные решения у(х), у'(х), у''(х) и у'''(х) на их приближения, найденные по схеме (1), получится верное равенство. Значит, близость у', у" , у"' с их приближениями влечёт близость у(х) и его приближения.

Для сравнения на рис. 1 изображены точное решение задачи

у' = 1 + У2, у\х=о = 0 ^ у = 1ап х (3)

и два приближенных её решения, найденные по схеме Эйлера и по схеме (1).

Точное решение имеет подвижный полюс при х = 2, вычисления по схеме Эйлера мы вынуждены остановить на подходе к этой точке, а вычисления по второй схеме можно благополучно продолжить и после особой точки.

На рис. 2 представлены приближенные решения задачи

у' = 1 + х2 + ху + х2у2, у\х=о = 0, (4)

найденные по схеме Рунге-Кутты и по схеме (1).

Рис. 1. Решения задачи (3): точное Рис. 2. Решения задачи (3): схема (пунктир), схема Эйлера (отдельные Рунге-Кутты (пунктир) и

точки) и (1,1)-схема (сплошная) (1,1)-схема (сплошная)

У

У

4

2

X

1

Как видно, решение имеет подвижный полюс где-то у х = 1.2, вычисления по явной схеме мы вынуждены остановить, а вычисления по второй схеме дают на удивление правдоподобный результат и после прохождения через особую точку.

Любопытно, что А.Н. Крылов в [4, С. 394] проиллюстрировал недостатки метода Эйлера вслед за Эйлером примером уравнения Риккати, отмечая «невыгодное обстоятельство», которое имеет место там, где f (х,у) или уничтожается, или становится бесконечно большой. Однако в случае уравнения Риккати все эти обстоятельства происходят от того, что схема Эйлера не сохраняет ангармоническое отношение четырёх частных решений и без труда устраняются без каких-либо выделений особенностей.

s

Вестник РУДН. Серия Математика. Информатика. Физика. №3,2015. С. 5-9

3. Схемы типа (2,2)

Обратимся теперь к уравнениям, когда точное решение начальной задачи задаёт (2, 2)-значное алгебраическое соответствие между слоями. Для примера общее решение уравнения Гамбургера = — ^ даётся формулой х2 + у2 = С; соотношение ж^ + у2 = + задаёт (2, 2)-соответствие между слоями. Схема Эйлера задаёт (1, 2)-соответствие между слоями, как видно на рис. 3, решение, найденное по этой схеме, разваливается возле особой точки и начинает случайным образом «прыгать» с одной интегральной кривой на другую.

Однако нетрудно составить схему, задающую (2, 2)-соответствие между слоями: = у2 — 2ж„Дж; решение найденное по этой схеме, см. рис. 4, легко проходит через алгебраическую особую точку.

Рис. З. Решение уравнения Гамбургера, найденное по методу Эйлера

Рис. 4. Вещественная (сплошная) и мнимая части (точки) решения уравнения Гамбургера, найденного по (2,2)-схеме

У

X

X

4. Заключение

Все сказанное можно подытожить следующим образом: если общее решение дифференциального уравнения l-го порядка с рациональной правой частью зависит от константы алгебраически, то это уравнение не только алгебраической заменой сводится к такому, которое обладает свойством Пенлеве, но и можно составить такую разностную схему, аппроксимирующую это уравнение, которая задаёт (I, ¿)-значное алгебраическое соответствие между любыми двумя слоями; эта схема замечательна тем, что она легко перешагивает через подвижные особые точки.

Приведённые в статье рисунки и вычисление были выполнены при помощи Sage Mathematics Software [7].

Литература

1. Painlevé P. Leçons sur la theorie analytique des equations differentielles. — Paris, l897. — Перепечатаны в первом томе Трудов Пенлеве, l97l.

2. Umemura H. Birational Automorphism Groups and Differential Equations // Nagoya Math. J. — 1990. — Vol. 119. — Pp. 1-80.

3. Малых М. Д. Об интегралах систем обыкновенных дифференциальных уравнений, представимых в конечном виде // Вестник РУДН. Серия «Математика. Информатика. Физика». — 20l4. — № 3. — С. ll-l6.

4. Крылов А. Н. Лекции о приближенных вычислениях. — Л.: АН СССР, l933.

5. Okamoto K. Sur les feuilletages associes aux equation du second odre apoints critiques fixes de P. Painleve // Japan. J. Math. — 1979. — T. 5, n0 1. — Pp. 1-79.

6. Калиткин Н. Н. Численные методы. — СПб.: БХВ-Петербург, 2011.

7. Stein W. A. et al., 2015. — Sage Mathematics Software (Version 6.7). — The Sage Development Team. — http://www.sagemath.org.

UDC 517.9

On the Approximate Solving of the Differential Equations which General Solutions Depend on a Constant of Integration

Algebraically

M. D. Malykh

Faculty of Materials Sciences Lomonosov Moscow State University GSP-1, Leninskie Gory, Moscow, Russia, 119991 Department of Applied Probability and Informatics Peoples' Friendship University of Russia 6, Miklukho-Maklaya str., Moscow, Russian Federation, 117198

Methods of analytically theory of ordinary differential equations are based on the analysis of singularities, but the most popular method for the numerical solving, that is the method of finite differences does not work in neighborhood of singularities. However Painleve gave an algebraic method for the solution of the differential equations in finite terms and general solutions of this equations depend algebraically on a constant of integration. This approach which was presented as Galois theory on the contrary can be well combined with method of finite differences.

It is well known, the ordinary differential equation of form y' = f (x, y) with this property can be algebraically transformed by substitution to Riccati equation. Euler scheme yn+i = yn + f (xn,yn)Ax always determines (1, ¿^-correspondence between neighboring layers. But exact solution of Riccati equation determines (1,1)-correspondence between any layers and thus we can write a scheme which determines (1,1)-correspondence between neighboring layers. In this case anharmonic ratio of 4 points does not change from layer to layer not only for exact solution but also for approximate solution. Thus if an exact solution has a pole then the approximate solution passes through infinity without accumulation of an error. In the presented article this property of (1,1)-scheme will be illustrated by two examples: with and without solution in elementary functions. So the cause of destruction of the approximate solution near a pole is put in Euler scheme itself. In more general case when exact solution of ordinary differential equation depends algebraically on an integration constant we can write a scheme which determines (I, ^-correspondence between neighboring layers. Approximate solution which is found on this way passes through movable algebraic singularities without accumulation of an error.

Key words and phrases: Riccati equation, finite differences method, algebraic correspondence, sage, sagemath.

References

1. P. Painlevé, Leçons sur la theorie analytique des equations différentielles, Paris, 1897, reprinted in the first volume of Penleve's work, 1971.

2. H. Umemura, Birational Automorphism Groups and Differential Equations, Nagoya Math. J. 119 (1990) 1-80.

3. M. D. Malykh, On Integrals of Ordinary Differential Equations Systems which are Representable in Finite Terms, Bulletin of Peoples' Friendship University of Russia. Series "Mathematics. Information Sciences. Physics" (3) (2014) 11-16, in Russian.

4. A. N. Krylov, Lectures on Numerical Calculations, AN USSR, Leningrad, 1933, in Russian.

5. K. Okamoto, Sur les feuilletages associes aux equation du second odre a points critiques fixes de P. Painleve, Japan. J. Math. 5 (1) (1979) 1-79.

6. N. N. Kalitkin, Numerical Methods, BHV-SPb, SPb., 2011, in Russian.

7. W. A. Stein, et al., Sage Mathematics Software (Version 6.7), The Sage Development Team (2015).

URL http://www.sagemath.org