Для систем обыкновенных дифференциальных уравнений (ОДУ) с невырожденной линейной частью в общем и гамильтоновом случаях ставится задача отыскания инвариантных координатных подпространств в координатах ее нормальной формы, вычисленной вблизи положения равновесия. Приведены условия существования таких инвариантных подпространств в терминах резонансных соотношений между собственными числами линейной части системы. Дан алгоритм поиска резонансных соотношений между собственными числами без их явного вычисления, который существенно использует методы компьютерной алгебры и q-аналог субрезультантов многочлена. Обсуждается его реализация в трех распространенных системах компьютерной алгебры – Mathematica, Maple и SymPy. Приведены содержательные модельные примеры.
Идентификаторы и классификаторы
Подход А. Пуанкаре к исследованию систем нелинейных обыкновенных дифференциальных уравнений (ОДУ) заключался в том, чтобы посредством аналитических обратимых преобразований максимально упростить правые части системы. Этот подход привел к теории нормальныхформ систем ОДУ, разработанной в работах А. Дюляка и А.Д. Брюно для общих систем (подробнее см. [1]) и в работах Дж. Биркгофа, Т.М. Черри, Ф.Г. Густавсона, К.Л. Зигеля, Ю. Мозера, А.Д. Брюно и других для систем Гамильтона (подробнее см. [2, Гл. I, II]).
Хотя нормальная форма системы ОДУ в окрестности инвариантного многообразия (положения равновесия, периодического решения, k-мерного тора) является формальным объектом, т.е. переход к нормальной форме обычно является расходящимся преобразованием, она может быть эффективно использована для исследования устойчивости соответствующего инвариантного многообразия [3], [4, Part II], локальной интегрируемости системы в его окрестности [5], поиска периодических решений, первых интегралов [6, 7], а также асимптотического интегрирования исходной системы [8].
Список литературы
- Брюно А.Д. Аналитическая форма дифференциальных уравнений (I) // Тр. ММО. 1971. Т. 25. С. 119-262.
- Брюно А.Д. Ограниченная задача трех тел: Плоские периодические орбиты. М.: Наука, 1990. 296 с.
- Маркеев А.П. Точки либрации в небесной механике и космодинамике. М.: Гл. ред. физ.-мат. литер. изд-ва “Наука”, 1978. 352 с. EDN: TTDUOJ
- Bruno A.D. Local Methods in Nonlinear Differential Equations. Berlin-Heidelberg-New York-London-Paris-Tokyo: Springer-Verlag, 1989. 350 p.
- Bruno A.D., Enderal V.F., Romanovski V.G. On new integrals of the Algaba-Gamero-Garcia system // Computer Algebra in Scientific Computing / Ed. by V.P. Gerdt, W. Koepf, E.W. Mayr, E.V. Vorozhtsov. Berlin Heidelberg: Springer, 2017. V. 10490 of Lecture Notes in Computer Science. P. 40-50.
- Брюно А.Д. Нормальная форма системы Гамильтона с периодическим возмущением // ЖВММФ. 2020. Т. 60. № 1. С. 36-52. EDN: UYZIHM
- Брюно А.Д. Нормализация периодической системы Гамильтона // Программирование. 2020. Т. 46. № 2. С. 76-83. EDN: TIHBWJ
- Журавлев В.Ф., Петров А.Г., Шундерюк М.М. Избранные задачи гамильтоновой механики. М.: ЛЕНАНД, 2015. 304 с.
- Wolfram S. The Mathematica Book. Wolfram Media, Inc., 2003. 1488 p.
-
Thompson I. Understanding Maple. Cambridge University Press, 2016. 228 p.
-
SymPy: symbolic computing in Python / Aaron Meurer, Christopher P. Smith, Mateusz Paprocki et al. // Peer J. Computer Science. 2017. V. 3. P. e103. DOI: 10.7717/peerjcs.103
-
Брюно А.Д. Аналитическая форма дифференциальных уравнений (II) // Тр. ММО. 1972. Т. 26. С. 199-239.
-
Биркгоф Д.Д. Динамические системы. Ижевск: Изд. дом "Удмуртский университет", 1999. 408 с.
-
Lyapunov A.M. Problme gnral de la stabilit du monvement // Ann. Fac. Sci. Toulouse Math. 1892. V. 9 (1907). № 2. P. 204474. Ann Math. Studies (1947). Stability of Motion. New York, London: Academic Press, 1966 (English).
-
Зигель К., Мозер Ю. Лекции по небесной механике. Ижевск: НИЦ "Регулярная и хаотическая динамика", 2001. 384 с. EDN: RXGNIR
-
Калинина Е.А., Утешев А.Ю. Теория исключения: Учеб. пособие. СПб.: Изд-во НИИ химии СПбГУ, 2002. 72 с.
-
Basu S., Pollack R., Roy M.-F. Algorithms in Real Algebraic Geometry. Algorithms and Computations in Mathematics 10. Berlin Heidelberg New York: Springer-Verlag, 2006. ix+662 p.
-
Сушкевич А.К. Основы высшей алгебры. 4-е изд. М.-Л.: ОГИЗ ГИТТЛ, 1941. 460 с.
-
Gathen J. von zur, Leking T. Subresultants revisited // Theoretical Computer Science. 2003. V. 297. P. 199-239. EDN: LSVFOF
-
Батхин А.Б. Параметризация дискриминантного множества вещественного многочлена // Программирование. 2016. Т. 42. № 2. С. 8-21. EDN: VWHVXP
-
Кац В.Г., Чен П. Квантовый анализ. М.: МЦНМО, 2005. 128 с.
-
Батхин А.Б. Параметризация множества, определяемого обобщенным дискриминантом многочлена // Программирование. 2018. № 2. С. 5-17. EDN: YWIHFZ
-
Батхин А.Б. Вычисление резонансного множества многочлена при ограничениях на коэффициенты // Программирование. 2019. № 2. С. 6-15. EDN: HTCHFU
-
Батхин А.Б. Вычисление обобщенного дискриминанта вещественного многочлена // Препринты ИПМ им. М.В. Келдыша. 2017. № 88. 40 с. http://www.keldysh.ru/papers/2016/prep2017_88.pdf. EDN: WTMHWL
-
Джури Э. Инноры и устойчивость динамических систем. М., 1979. 304 с.
-
Акритас А.Г. Основы компьютерной алгебры с приложениями. М.: Мир, 1994. 544 с.
-
Батхин А.Б., Брюно А.Д., Варин В.П. Множества устойчивости многопараметрических гамильтоновых систем // Прикл. мат. мех. 2012. Т. 76. № 1. С. 80-133. EDN: OWOHAR
-
Маркеев А.П. Линейные гамильтоновы системы и некоторые задачи об устойчивости движения спутника относительно центра масс. М.-Ижевск: НИЦ "Регулярная и хаотическая динамика", 2009. 396 с. EDN: TTDUHB
-
Маркеев А.П. Нелинейные колебания симпатических маятников // Нелинейная динамика. 2010. Т. 6. № 3. С. 605-621. EDN: MVDGTN
Выпуск
Другие статьи выпуска
Предлагается пакет символьного построения экспоненциально-логарифмических решений таких линейных обыкновенных дифференциальных уравнений, которые могут иметь неполностью заданные (усеченные) коэффициенты – степенные ряды, для которых известны только начальные члены. Входящие в решения ряды также представляются конечным числом начальных членов. Для каждого такого ряда вычисляется максимально возможное число начальных членов, полностью определенных известными начальными фрагментами коэффициентов уравнения. При этом степень усечения каждого из этих рядов не может превосходить заданной пользователем величины. Последнее обеспечивает окончание вычисления и тогда, когда по известным фрагментам коэффициентов уравнения может быть определено любое число членов рядов, входящих в решения.
На основе рекуррентных формул Ньютона приведен алгоритм построения результанта двух целых функций, что является одним из методов исключения неизвестных для систем неалгебраических уравнений. Алгоритм реализован в системе компьютерной алгебры Maple. Приведены примеры, демонстрирующие работу данного алгоритма.
В последнее время на место основного языка научных и инженерных расчетов выдвигается язык Julia. У ряда пользователей возникает желание работать полностью внутри “экосистемы” Julia, подобно тому, как происходит работа в “экосистеме” Python. Для Julia существуют библиотеки, покрывающие большинство потребностей научно-инженерных расчетов. Перед авторами возникла необходимость использовать символьные вычисления для задач математического моделирования. Поскольку основным языком реализации численных алгоритмов мы выбрали язык Julia, то и задачи компьютерной алгебры хотелось бы решать на этом же языке. Авторы выделили основные функциональные области, задающие разные варианты применения систем компьютерной алгебры. В каждой из областей нами выделены наиболее характерные представители систем компьютерной алгебры на Julia. В результате авторы делают вывод, что в рамках “экосистемы” Julia возможно (и даже удобно) использовать системы компьютерной алгебры.
Рассматривается поведение квантовой запутанности в процессе унитарной эволюции в конструктивных моделях многокомпонентных квантовых систем. Описываются группы симметрий квантовых систем, допускающих возникновение геометрических структур, ассоциированных с квантовой запутанностью. Алгоритмы моделирования динамики квантовой запутанности основаны на методах компьютерной алгебры и вычислительной теории групп. Приводятся примеры конкретных вычислений.
В статье предложены два наиболее простых метода определения положений равновесия спутника, движущегося в центральном ньютоновом силовом поле по круговой орбите под действием гравитационного момента. В первом методе применялись подходы линейной алгебры, во втором алгоритмы компьютерной алгебры. Положения равновесия спутника в орбитальной системе координат при заданных значениях главных центральных моментов инерции определяются корнями системы нелинейных алгебраических уравнений. Для определения равновесных решений проводилась декомпозиция системы алгебраических уравнений с применением методов линейной алгебры и алгоритмов построения базисов Гребнера. Положения равновесия спутника определялись путем исследования числа действительных корней алгебраических уравнений из полученных базисов Гребнера. С использованием предложенного подхода показано, что спутник с неравными главными центральными моментами инерции имеет на круговой орбите 24 положения равновесия.
Факторизация полиномов – классическая алгоритмическая проблема алгебры, имеюшая широкий спектр приложений. Особый интерес представляет факторизация над конечными полями, среди которых поле порядка два является, вероятно, наиболее важным в связи с представлением булевых функций полиномами Жегалкина. В частности, факторизация булевых полиномов соответствует конъюнктивной декомпозиции булевых функций, заданных в алгебраической нормальной форме. Кроме того, факторизация дает решение проблемы декомпозиции функций, заданных в СДНФ и позитивных ДНФ, а также декартовой декомпозиции реляционых данных. Эти приложения демонстрируют важность разработки быстрых алгоритмов факторизации. В статье мы рассматриваем некоторые недавно предложенные алгоритмы факторизации полиномиальной сложности и описываем параллельную MIMD-реализацию, которая использует как параллелизм уровня задачи, так и параллелизм уровня данных. Мы представляем эксперименты, выполненные на бенчмарках логического синтеза и на синтетических (случайных) полиномах, которые показывают значительное ускорение факторизации. В заключение представлены результаты тестирования параллельной реализации алгоритма на массивнопараллельной многоядерной архитектуре (Redefine).
Издательство
- Издательство
- ИЗДАТЕЛЬСТВО НАУКА
- Регион
- Россия, Москва
- Почтовый адрес
- 121099 г. Москва, Шубинский пер., 6, стр. 1
- Юр. адрес
- 121099 г. Москва, Шубинский пер., 6, стр. 1
- ФИО
- Николай Николаевич Федосеенков (Директор)
- E-mail адрес
- info@naukapublishers.ru
- Контактный телефон
- +7 (495) 2767735