Методика отыскания алгебраических интегралов дифференциальных уравнений первого порядка

Полный текст

1. Введение При решении системы дифференциальных уравнений чрезвычайно важно найти все алгебраические интегралы. Авторы XIX века видели в их отыскании первый шаг к построению решения в конечном виде, однако даже в том случае, когда система не допускает решения в символьном виде и поэтому будет решаться приближённо по методу конечных разностей, всегда стремятся строить схемы консервативные, то есть сохраняющие алгебраические интегралы [1]. Применительно к уравнению 1-го порядка её можно сформулировать в следующем виде. Задача 1 (Дебон). Для заданного дифференциального уравнения �d� + �d� = 0, �, � ∈ Q[�, �] (1) выяснить, имеет ли оно интеграл вида �(�, �) = const, � ∈ C(�, �), и в случае утвердительного ответа предъявить этот интеграл. Эта задача, пусть и в других терминах, была предложена Дебоном (Florimond de Beaune, 1601-1652) Декарту ещё в 1640-х годах [2], однако она до сих пор в такой постановке не решена: пользователь любого пакета, ищущего рациональный интеграл, должен задавать верхнюю границу для порядка искомого интеграла [3]. Возникающее тут затруднение можно проиллюстрировать очень простым примером. Пример 1. Пусть дано уравнение d� = �d�. Вслед за Декартом зададимся большим числом � и будем искать его интеграл вида ∑︁ �+� � ��,����� по методу неопределённых коэффициентов. Увеличивая �, скажем, от 2 до 10, мы получим несовместную систему для коэффициентов, однако из этого не следует, что так будет проходить при всех �. Несовместность систем при всех � означает, что �� не является алгебраической функцией. Мы, в отличие от Декарта, можем доказать отсутствие алгебраического интеграла очень просто, сославшись на теорему Сохоцкого-Вейерштрасса или любой другой вариант описания существенной особенности при � = ∞. Однако это доказательство не является алгебраическим и существенно использует тот факт, что уравнение d� = �d� допускает решение в целых функциях. Аналитики XIX века ожидали, что анализ особенностей должен указать на существование или несуществование алгебраического интеграла. Однако в силу теоремы Пенлеве все подвижные особенности являются алгебраическими [4], и поэтому по поведению решения в подвижных особых точках нельзя ничего сказать о существовании интеграла. Исследование же поведения решения в неподвижной особой точке представляет само по себе большие трудности. По этой причине мы обратимся к более простой задаче. Задача 2. Для заданного дифференциального уравнения (1) и числа � выяснить, имеет ли оно рациональный интеграл, порядок которого не превосходит �, и в случае утвердительного ответа предъявить этот интеграл. Эта задача может быть решена методом неопределённых коэффициентов, однако на практике такое решение приводит к большим системам нелинейных уравнений на коэффициенты, решение которой приводит к большим вычислительным затруднениям. Накануне Первой мировой войны М. Н. Лагутинский [5-7] разработал метод решения «ограниченной» задачи Дебона. Его метод был надолго забыт и открыт вновь в недавних работах по компьютерной алгебре [8]. Современное изложение метода Лагутинского представлено в работах [3, 8, 9]. Метод Лагутинского хорошо ложится на язык теории колец. Здесь мы хотим обсудить вопросы применения его реализации - пакета Lagutinski для Sage [10] - к решению наиболее часто встречающихся задач. Согласно основной теореме Лагутинского дифференциальное уравнение (1) имеет рациональный интеграл порядка � тогда и только тогда, когда определитель Лагутинского соответствующего порядка равен нулю, а сам интеграл можно записать как отношение миноров этого определителя. Эта теорема сводит решение ограниченной задачи Дебона к вычислению определителей Лагутинского. Замечание 1. Теорема Лагутинского позволяет менять числовые поля в исходной задаче. Если дифференциальное уравнение допускает интеграл в поле C(�, �), то его определитель Лагутинского достаточно большого порядка равен нулю, а поэтому это уравнение допускает интеграл, который представляет собой отношение миноров названного определителя. По предположению коэффициенты дифференциального уравнения - числа целые, поэтому эти миноры принадлежат Q(�, �), а следовательно, этому полю принадлежит интеграл. Поэтому всегда можно ограничивать рассмотрение интегралами из поля Q(�, �). Это замечание увязывает решение задачи Дебона с вопросами теории чисел. Если дифференциальное уравнение допускает решение, которое при целом � принимает трансцендентное значение, то оно Малых М. Д., Юй Ин Методика отыскания алгебраических интегралов . . . 287 не допускает алгебраических интегралов. Например, из трансцендентности числа � сразу следует, что уравнение d� = �d� не допускает интеграл в поле C(�, �). На практике обычно имеют дело с большими выборками дифференциальных уравнений, будь то задачи из того или иного задачника, или возможные модели для одного и того же явления. Требуется выявить из них те, которые допускают рациональный интеграл. Число �, конечно, не дано, однако очень большие значения � не интересны, поскольку с интегралом 100-го порядка едва ли работать легче, чем со степенными рядами. Вообще говоря, не ясно, как провести эту черту, и вероятно, она зависит от самой задачи. Для дифференциальных уравнений, предлагаемых студентам в курсе «Дифференциальных уравнений», в статье [10] была предложена следующая процедура отыскания алгебраических интегральных кривых: 1. Вычислить ∆55 в случайно выбранной точке. Если ∆55 ̸= 0, то интегральные кривые имеют 10-й порядок и более или вовсе являются трансцендентными. Если ∆55 = 0, то интегральные кривые вероятно имеют порядок, не превышающий 9. 2. Если ∆55 = 0, вычислить ∆10 как функцию � и �. Если ∆10 тождественно равен нулю, то интегральные кривые имеют порядок, не превышающий 3, и они находятся по методу Лагутинского без существенных затрат. 3. Если ∆10 не равно тождественно нулю, пытаться вычислить ∆� , постепенно повышая � до тех пор, пока не получится нуль или исчерпаются ресурсы. Поскольку вычисление определителя 55-го порядка может занять весьма много времени, на первом шаге его значение вычисляется только в одной точке. Если получается ненулевое значение, то в силу теоремы Лагутинского рационального интеграла 9-го порядка точно не существует. Если же получается нулевое значения, нельзя быть вполне уверенным в том, что определитель равен нулю и во всех других точках, однако вероятность того, что при случайном выборе точки вероятность отыскания нуля многочлена весьма мала. Поэтому в этом случае вероятно имеется интеграл 9-го порядка. Разумеется, на первом шаге теряются дифференциальные уравнения, интегральные кривые которых имеют 10-й порядок и более. На третьем шаге могут возникнуть трудности с вычислением определителя Лагутинского, обычно при � > 10 работа с ним становиться очень ресурсоёмкой. Поэтому представляется важным понять, работает ли предложенный метод в наиболее часто встречающихся в учебном процессе дифференциальных уравнений. В настоящей работе в качестве собрания таких уравнений использован сборник задач, предлагаемых студентам в курсе «Дифференциальных уравнений». Ниже будет дан отчёт об использовании этой процедуры для интегрирования дифференциальных уравнений 1-го порядка в символьном виде и тестирования названного в его честь пакета на уравнениях, взятых из задачника А. Ф. Филиппова [11]. 1. Результаты апробации На основе задач №№ 301-420 из задачника А. Ф. Филиппова [11] был выделен список из 50 дифференциальных уравнений вида �(�, �)d� + �(�, �)d� = 0, �, � ∈ Q[�, �], 20 из которых имеют рациональные интегралы, а остальные интегрируются в элементарных функциях. Прежде всего, следует заметить, что вычисление ∆55 в случайной точке позволило быстро выделить все 20 номеров, вероятно допускающих рациональные интегралы. При этом про оставшиеся 30 номеров можно с уверенностью сказать, что их интегральные кривые или являются трансцендентными, или имеют порядок, больший 9-го. Во-вторых, вычисление ∆� в случайной точке позволило быстро подобрать � во всех номерах. 288 Вестник РУДН. Серия МИФ. Т. 26, № 3, 2018. С. 285-291 Вычисление определителя как многочлена относительно �, � заняло более часа в двух примерах: №№ 395 и 418. Вычисление интеграла заняло более часа в задаче № 361 и было остановлено. Во всех остальных 17 случаях вычисление интеграла было успешным во всех случаях, ответы совпадали с ответами из задачника после упрощений. Рассмотрим проблемные номера. Пример 2. Интегральные кривые уравнения № 395 (�3 - 2��2)d� + 3�2�d� = �d� - �d� образуют семейство кривых �(�2 + �2)3 - 2/5�5 - 4/3�2�3 - 2�4� + ��5 = 0 6-го порядка. Вычисление определителей в случайной точке однозначно указали на то, что вероятно ∆36 = 0, однако вычислить этот определитель за час не удалось. Пример 3. Интегральные кривые уравнения № 418 (3�� + � + �)�d� + (4�� + � + 2�)�d� = 0 образуют семейство кривых 6�3�4 + 2�3�3 + 3�2�4 = � 7-го порядка. Вычисление определителей в случайной точке однозначно указали на то, что вероятно ∆31 = 0, однако вычислить этот определитель за час не удалось. Пример 4. Интегральные кривые уравнения № 361 (2�2� - 3�2)�′ = 6�2 - 2�� + 1 образуют семейство 2�3 - �2�2 + �3 + � = � 4-го порядка. Мы вычислили ∆13 = 0, ∆12 ̸= 0, поэтому точно существует рациональный интеграл порядка � = 13. Вычислить его не удалось. Глядя в ответы, полученные, разумеется, другими методами, нетрудно понять, что предложенный подход позволил легко и быстро выявить все уравнения, допускающие алгебраические интегралы, но не для всех из них позволил найти явное выражение для интеграла. Наиболее проблемным, таким образом, является третий шаг. 2. Заключение Проделанные численные эксперименты свидетельствуют, что метод Лагутинского позволяет на практике быстро и без особых затрат ресурсов и времени идентифицировать наличие рационального интеграла, однако является весьма затратным как метод вычисления этого интеграла. Следует заметить, что проблема отыскания границы для порядка интеграла, всегда обсуждаемая в теории, на практике оказалась не важной, в задачнике нет дифференциальных уравнений, интегральные кривые которых имели бы 10 порядок или более. Вычисление определителей в случайных точках позволяет быстро выяснить, какие мономы должны фигурировать в выражении для интеграла, поэтому естественное направление оптимизации - предварительный подбор базиса по данным в случайных точках. В частности, наиболее интересны те интегралы, которые можно записать как отношение малочленов. В частности, в пакете Lagutinski [10] можно делать перебор по всевозможным тройкам, составленных из первых � мономов кольца Q[�, �]. Двигаясь в этом направлении, можно пытаться увеличить число элементов в кортеже, заменив затратное вычисление определителей Лагутинского вычислением их значения в случайной точке.

Об авторах

Михаил Дмитриевич Малых

Российский университет дружбы народов

Автор, ответственный за переписку.
Email: malykh_md@rudn.university

кандидат физико-математических наук, доцент кафедры прикладной информатики и теории вероятностей РУДН

ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198

Ин Юй

Российский университет дружбы народов

Email: yingy6165@gmail.com

аспирант кафедры прикладной информатики и теории вероятностей РУДН

ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198

Список литературы