NUMERICAL-ANALYTICAL TRANSFORMATION METHOD OF ANALYZING NONLINEAR MATHEMATICAL MODELS WITH POLYNOMIAL STRUCTURE
Journal: VESTNIK OF ASTRAKHAN STATE TECHNICAL UNIVERSITY. SERIES: MANAGEMENT, COMPUTER SCIENCE AND INFORMATICS
Volume 2022 № 2
, 2022
Rubrics: MATHEMATICAL MODELING
Authors:
1.
ITMO University
(Candidate of Physics and Mathematics, Assistant Professor; Assistant Professor of the Department of Information and Communication Technologies )
, Russian Federation
, Russian Federation
2.
ITMO University
(Postgraduate Student of the Department of Information and Communication Technologies )
Sankt-Peterburg, Russian Federation
Sankt-Peterburg, Russian Federation
Type:
Article
Pages:
from 97
to 109
Status:
Published
Received:
26.01.2022
Accepted:
08.04.2022
Published:
22.04.2022
Subject area:
UDK 004.021 Алгоритмы
Language:
Russian
Keywords:
method of transformations, nonlinear dynamic systems, mathematical model of polynomial structure, stability of system solution, numerical and analytical methods
Abstract (English):
The article considers a scientific problem of investigating mathematical models with polynomial nonlinearity, which are represented by systems of nonlinear differential equations. There is presented a numericall-analytical transformation method for investigating nonlinear dynamical systems with many degrees of freedom of polynomial structure. As opposed to the analogues, this method allows solving a wide range of problems for nonlinear systems of general polynomial structure while reducing the computational resource intensity. An algorithm of the method of transformations for investigating the nonlinear dynamic systems with m degrees of freedom is given. The research of stability of solutions of nonlinear dynamic systems with m degrees of freedom is carried out. Algorithmic formulas of transformation method for the solution of systems with nonlinearity of the fourth degree are given. Computational experiments on solving systems of differential equations with a small nonlinear part in the form of a multinomial of the sixth degree using the transformation method show fourth order accuracy in computation. A nonlinear mathematical model of the vibration protection system is investigated by the presented method of transformations. The theorem on determination of the stationary state for systems of differential equations of polynomial structure by transformation method is proved. Algorithmic formulas for computation are presented. A general matrix form for vector indices is presented. Formulas for economical computation of right-hand sides of polynomial structure are presented and it is proposed to apply Pan's scheme with coefficient preprocessing. This method allows studying the different modes of nonlinear models dynamics, for example, to determine such extreme modes as resonance, sub-harmonic, and polyharmonic modes. As an example of the method application, the vibration protection problem of the tower from external periodic influences was solved. The transformed solution takes into account all nonlinear polynomial components. The method allows to investigate the dynamics of a wide range of nonlinear systems with the necessary accuracy
The article considers a scientific problem of investigating mathematical models with polynomial nonlinearity, which are represented by systems of nonlinear differential equations. There is presented a numericall-analytical transformation method for investigating nonlinear dynamical systems with many degrees of freedom of polynomial structure. As opposed to the analogues, this method allows solving a wide range of problems for nonlinear systems of general polynomial structure while reducing the computational resource intensity. An algorithm of the method of transformations for investigating the nonlinear dynamic systems with m degrees of freedom is given. The research of stability of solutions of nonlinear dynamic systems with m degrees of freedom is carried out. Algorithmic formulas of transformation method for the solution of systems with nonlinearity of the fourth degree are given. Computational experiments on solving systems of differential equations with a small nonlinear part in the form of a multinomial of the sixth degree using the transformation method show fourth order accuracy in computation. A nonlinear mathematical model of the vibration protection system is investigated by the presented method of transformations. The theorem on determination of the stationary state for systems of differential equations of polynomial structure by transformation method is proved. Algorithmic formulas for computation are presented. A general matrix form for vector indices is presented. Formulas for economical computation of right-hand sides of polynomial structure are presented and it is proposed to apply Pan's scheme with coefficient preprocessing. This method allows studying the different modes of nonlinear models dynamics, for example, to determine such extreme modes as resonance, sub-harmonic, and polyharmonic modes. As an example of the method application, the vibration protection problem of the tower from external periodic influences was solved. The transformed solution takes into account all nonlinear polynomial components. The method allows to investigate the dynamics of a wide range of nonlinear systems with the necessary accuracy
Keywords:
method of transformations, nonlinear dynamic systems, mathematical model of polynomial structure, stability of system solution, numerical and analytical methods
method of transformations, nonlinear dynamic systems, mathematical model of polynomial structure, stability of system solution, numerical and analytical methods
Введение
Рассматривается научная задача исследования математических моделей с полиномиальной нелинейностью, которые представлены системами нелинейных дифференциальных уравнений. Для ис-следования таких нелинейных моделей разрабатывается метод преобразований. Традиционные методы исследования нелинейных систем, например преобразование к линейной форме, не учитывают влияние нелинейностей, приводят к упрощению модели и существенным погрешностям. Классические методы исследования нелинейных математических моделей были предложены в работах
А. Пуанкаре, А. М. Ляпунова, Ван дер Поля,
Н. М. Крылова, Н. Н. Боголюбова, А. А. Андронова и др. Применяются аналитические приближенные методы малого параметра, усреднения, Крылова–Боголюбова, линеаризации, метод Ван дер Поля, метод гармонического баланса, метод воз-мущений Пуанкаре, асимптотические методы растянутых параметров и многих масштабов. Каждый метод предназначен для решения определенного класса задач с введенными ограничениями. Например, метод усреднения не учитывает квадратичные составляющие, методы Ван дер Поля
и линеаризации не учитывают кубические составляющие. Определенная ограниченность области решения является значительным недостатком применения приближений Чебышева при построении решений нелинейных дифференциальных уравнений. С целью снижения ресурсоемкости решений задачи при использовании приближенных методов критически важно правильно определить началь-ное приближение [1]. Побочным эффектом некорректного определения начального приближения является не только ощутимое повышение трудоемкости вычислительного алгоритма, но и невысокая точность решения.
Приведем обзор актуального состояния научных исследований в рассматриваемой области.
В работе [2] нашло отражение использование методики нахождения полиномиальных решений канонических уравнений. Данная методика основана на разложении решений в ряды по базисным функциям. Далее, дополнительные решения могут находиться путем комбинирования решений, полученных ранее на предыдущих этапах. Подобная аппроксимация искомого решения легла в основу численно-аналитического метода, представленного автором [2]. Существенным преимуществом такого подхода является то, что решение представляется «…в форме непрерывной дифференцируемой по области решения, удовлетворяющих по точности аппроксимаций» [2, с. 204].
Авторы [3] посвятили свою работу исследованию решения систем нелинейных обыкновенных дифференциальных уравнений второго порядка. Данные уравнения являются уравнениями с выделенной главной положительно однородной частью. Работа обладает выраженной практической значимостью, т. к. ее результаты позволяют обобщать решения для подобных систем.
Один из методов численно-аналитического описания решения уравнения объекта представлен
в работе [4]. Помимо описания собственно решения в работе также вычисляются значения числовых характеристик на основе некорректных наблюдений. По мнению авторов, внедрение предложенного метода «…позволит алгоритмическими и программными средствами оперативно находить оптимальные оценки числовых характеристик движения объекта» [4, c. 11].
Одна из сфер применения численно-аналитического моделирования описана в работе [5]. Авторами была изучена сфера систем массового обслуживания поликомпонентных потоков. Построенные три математические модели и разнообразные области их применения являются существенным достоинством данной работы, а полученные численные результаты помогают рассчитать необходимое количество устройств для систем массового обслуживания поликомпонентных потоков.
Нахождение базисных функций с помощью численно-аналитического метода изучается в работе [6]. Предложенный авторами подход «…заключается в построении конечного элемента, учитывающего особенности решения задачи Дирихле для уравнения Лапласа в окрестности угловых точек» [6, c. 325]. В результате удалось значительно повысить точность решения из-за сокращения порядка системы уравнений.
В качестве еще одного варианта применения численно-аналитических методов следует упомянуть работу авторов [7], в основу которой лег метод поиска близких к периодическим решений, соответствующих режимам авторотации механической системы. Подход, являющийся родственным методу Андронова–Понтрягина, позволяет строить решения вспомогательных подсистем. В результате получается набор периодических функций при выполнении условия сходимости метода.
Современными проблемами исследования не-линейных моделей занимается коллектив Института прикладной математики РАН им. М. В. Келдыша. В работе [8] представлена процедура построения разложения базиса решений одного четырехчленного рекуррентного соотношения с коэффициентами из кольца. Показано, что вне зоны близких собственных значений получаемое разложение является асимптотическим.
Авторы [9] рассматривают метод решения жестких систем дифференциальных уравнений на основе вычисления матричной экспоненты. В работе представлен метод контроля точности в случае линейных задач, приведены результаты численных расчетов для линейных и нелинейных задач.
Дифференциальное уравнение второго порядка, содержащее большой параметр, рассматривается
в [10]. Авторами описан предложенный ранее способ, являющийся альтернативой методу Пуанкаре (теорема о неявной функции). Мотивом к данному исследованию послужило желание авторов изучить влияние применения метода Пуанкаре в сингулярно возмущенной задаче.
Сравнение методов построения многослойных приближенных решений дифференциальных уравнений на основе классических приближенных методов выполняется в работе [11]. Варианты реализации алгоритмов расчета переходных процессов для динамических систем с применением опера-торного метода рассматриваются в работе [12].
В работе [13] асимптотическими методами исследуются полиномиальные дифференциальные уравнения. Рассматривается разложение по степеням независимой переменной, коэффициенты которого есть ряды Лорана по убывающим степеням логарифма.
В работе [14] рассматриваются численно-аналитические методы поиска периодических решений нелинейных систем дифференциальных уравнений. Представлены алгоритмы отыскания начальных условий, соответствующих периодическому решению. В исследовании [15] приводится доказательство теоремы о сходимости формального решения обыкновенных дифференциальных уравнений при случае, когда «множество показателей степени ряда имеют одну комплексную образующую» [15, с. 2].
В работе [16] доказана теорема о разрешимости нелинейной системы уравнений относительно приближенных значений коэффициентов Чебышева старшей производной, входящей в дифференциальное уравнение. Авторам [17] для определения точных решений нелинейных дифференциальных уравнений используется метод логистической функции.
В работах научного коллектива под руководством профессора Г. И. Мельникова предложен метод дифференциальных неравенств и применение приближений Чебышева для построения решений нелинейных дифференциальных уравнений.
В работе [18] рассматривается система дифференциальных уравнений полиномиальной структуры. Проводится преобразование уравнений и исследование колебаний объекта с постоянными параметрами. Приводятся оценки движения, полученные методом дифференциальных неравенств для поло-жительных функций Ляпунова. В [19] рассматривается уравнение автономной динамической системы с одной степенью свободы. Применяется асимптотический метод преобразования Пуанкаре–Дюлака и приближения Чебышева для многочленов высо-кой степени. Исследуются механические системы
в случае отсутствия внутренних резонансов.
Развитие численно-аналитических методов исследований поведений систем нашло свое отражение в работе [20]. Авторами рассматривается случай, при котором наблюдения содержат отсчеты сингулярной помехи. Предложенное улучшение метода позволяет определять оптимальные несме-щенные и инвариантные оценки. Достоинством такого подхода является то, что он помогает избежать расширения пространства состояний.
В рассмотренных научных работах предложены методы для решения определенного класса задач при выполнении заданных условий. В отличие
от рассмотренных научных работ в статье предложен метод, позволяющий решать широкий класс задач для нелинейных систем общей полиномиальной структуры при сокращении ресурсоемкости вычислений.
Постановка задачи исследования
Рассматривается актуальная задача разработки оптимального алгоритма для метода исследования нелинейных моделей. Для увеличения точности исследования нелинейных моделей разрабатывается метод преобразований. Метод позволяет исследовать различные режимы динамики нелинейных моделей, например, определять такие экстремальные режимы, как резонанс, субгармонический, полигар-монический режим. Так, традиционные численные методы Рунге–Кутта с фиксированным шагом
не определяют эффект мгновенного скачка. Для современных систем управления в реальном времени динамическими системами задача своевременного определения и устранения экстремальных режимов является актуальной. При разработке алгоритма метода также решается важная задача сокращения вычислительных ресурсов по сравнению с анало-гичными классическими методами.
Метод преобразований для исследования нелинейных динамических систем с m степенями свободы
Рассмотрим динамическую систему с m степенями свободы. Для такой системы математическая модель представлена m нелинейными дифференциальными уравнениями второго порядка с постоянными и периодическими коэффициентами.
Представим систему в матричной форме с применением векторного индекса
с целочисленными неотрицательными компонентами.
Обозначим – сумму компонентов векторного индекса. Математическую модель системы запишем в матричном виде:
(1)
где I – единичная матрица m × m; B, C – квадратные матрицы m × m; – вектор-столбец нелинейных коэффициентов, ;
t – время; ω – частота; qi – обобщенные координаты системы; – вектор-столбец искомых функций; – вектор-столбец производных.
Предположим, что для системы дифференциальных уравнений (1) выполнены условия теоремы Пикара о существовании и единственности решения задачи Коши. Правая часть системы определена, непрерывна и удовлетворяет условию Липшица.
Для представления суммы по векторному индексу в работе [21] предложено применять последовательное раскрытие вложенных сумм вида
Основные этапы метода преобразований включают представление нелинейной системы в нормальной форме, линейное и многочленное преобразование, решение преобразованной системы.
В алгоритме введено условие комплексно-сопряженных корней с малыми отрицательными вещественными частями для уравнения .
Вводятся новые переменные
и .
Представим периодические функции через введенные дополнительные переменные: и .
Приведем систему уравнений (1) к нормальной форме, дополнив введенными переменными:
(2)
где ; P – квадратная матрица ; – комплексные коэффициенты нелинейной части, определяемые посредством перегруппировки членов нелинейных частей системы (1).
Выполним линейное преобразование системы:
где Y – вектор-столбец новых обобщенных координат системы; L – квадратная матрица линейного преобразования.
Определяем невырожденную матрицу L.
В результате линейного преобразования представим систему в виде
где y1, y2 – новые обобщенные координаты системы. Здесь матрица Λ диагонального вида
с комплексно-сопряженными корнями
.
Для нахождения нелинейных коэффициентов rν выполним перегруппировку членов в нелинейной части (2) после линейного преобразования.
На следующем этапе выполняется преобразование
,
; , (3)
где ys – обобщенные координаты системы; zs – преобразованные обобщенные координаты системы; – коэффициенты преобразования.
Результатом преобразования (3) является система
где – коэффициенты преобразованной системы.
Для нахождения особых значений индекса решаются уравнения [22]
(4)
Для не особых значений индексов определяются , для особых значений индексов определяются .
В работе [21] представлена итерационная формула для вычисления коэффициентов и :
Введем комплексно-сопряженные переменные
(5)
В результате перехода к переменным (5) получим автономную систему
(6)
Далее получаем решение автономной системы.
Для перехода к автономной системе с вещественными коэффициентами выполним экспоненциальную замену переменных:
. (7)
Представим систему (6) с вещественными коэффициентами в виде
(8)
В нерезонансном случае система (8) имеет более простой вид:
Проинтегрировав систему, найдем амплитуду ρS и частоту θS для преобразованных обобщенных координат системы. Установившийся режим найдем, приравнивая правые части (8) к нулю
и решая систему алгебраических уравнений. Для представления решения в исходных переменных подставим решение в (3), (5) и (7). Решение в ис-комых переменных X найдем по формулам обратной замены: .
Методом преобразований исходная нелинейная система дифференциальных уравнений второго порядка с постоянными и периодическими параметрами преобразуется к автономному виду. Для экономичного вычисления правых частей полиномиальной структуры предложено применить схему Пана с предварительной обработкой коэффициентов. Методы экономичного вычисления полиномов с предварительной обработкой коэффициентов представлены в работах В. Я. Пана [23]. В соответствии с двухэтапной схемой Пана для вычисления полинома n-й степени при n ≥ 5 необходимо умножений и n + 1 сложений. Двухэтапная схема экономичного вычисления Пана с предварительной обработкой коэффициентов позволяет в 2 раза сократить количество операций умножения при вы-числении полиномов, что приводит к значительному увеличению производительности вычисления полиномов.
Для вычислений особых индексов в алгоритме метода преобразований удобно применить матричную форму представления особых индексов.
Введем ленточную матрицу из единиц, размерностью .
.
Каждая строка матрицы представляет векторный индекс.
Для получения особых индексов для каждого
s построим сумму ленточной матрицы из единиц и матрицы IS с s – столбцом из единиц: .
Для получим матрицу для определения особых индексов вида
В нерезонансном случае система (4) для вычисления особых индексов имеет
следующий вид:
Подстановка значений векторных индексов (строк матриц ) в систему уравнений для доказывает правильность построения матричных форм определения особых индексов.
Каждая строка матрицы представляет особый векторный индекс.
С учетом того, что сумма особых индексов определяет степени мономов в преобразованной системе дифференциальных уравнений и равенства , в каждом дифференциальном уравнении преобразованной системы присутствуют только мономы первой и третьей степени от искомых переменных.
Методом преобразований система приводится
к автономному виду. Для автономных дифференциальных уравнений решения асимптотически приближаются к стационарному состоянию. При устойчивом стационарном состоянии малые от-клонения не выводят систему из малой окрестно-сти стационарного состояния. Для устойчивого стационарного состояния необходимо условие отрицательности вещественных частей характеристического уравнения.
Приведем теорему об определении стационарного состояния системы.
Теорема 1. При условии существования устойчивого стационарного состояния для системы m дифференциальных уравнений второго порядка в нерезонансном случае с малыми нелинейными частями в форме многочленов четвертых степеней существует система m линейных алгебраических уравнений с вещественными коэффициентами, полученная в результате многочленных преобразований, которая определяет стационарное состояние системы.
Доказательство. Представим доказательство для системы m нелинейных уравнений второго порядка.
Исследуем систему m уравнений второго по-рядка с малыми нелинейными частями, представ-ленными многочленами четвертой степени.
В случае отсутствия резонансов автономная
преобразованная система (8) имеет следующую форму:
Для определения особых индексов запишем матричные формы при
.
Каждая строка матрицы представляет особый векторный индекс. Для суммы особых индексов выполняются равенства при , при .
С учетом того, что сумма особых индексов определяет степени мономов в преобразованной системе дифференциальных уравнений и равенства , в каждом уравнении преобразованной системы присутствуют только мономы первой и третьей степени от искомых переменных. Для определения стационарного решения разделим на первое уравнение (8), приравняем правые части к нулю, получим систему m алгебраических уравнений
.
При делении на ρS степени мономов понизились на единицу, и в уравнении присутствуют только мономы нулевой (свободные члены) и второй степени от искомых переменных.
Проведя квадратичную замену переменных , получаем линейную систему m алгебраических уравнений с вещественными коэффициентами вида
В линейной системе алгебраических уравнений будут присутствовать только мономы нулевой и первой степени в соответствии с равенствами для особых индексов: при ; при .
Таким образом, теорема доказана и методом преобразований определена система линейных алгебраических уравнений.
Получив стационарные значения ρS, из второго уравнения (8) определим θS:
Таким образом, определено стационарное решение для системы m дифференциальных уравнений второго порядка с нелинейными частями
в форме полиномов четвертых степеней в нерезонансном случае.
Оценка погрешностей решений
Представим оценку погрешностей решений, полученных методом преобразований. В работе
Г. И. Мельникова [22] методом дифференциальных неравенств решается задача оценки погрешности решений, полученных методом преобразований.
Рассматривается метрика для определения погрешности , где – точное
и приближенное решение в окрестности , с одинаковым начальным значением .
В работе Г. И. Мельникова получена оценка погрешности при , которая определена по теореме о дифференциальных неравенствах [22]. Погрешность оценивается функцией вида , где u – аргумент функции оценки; ε – малая величина; , , – коэффи-циенты нелинейной правой части.
Таким образом, определена оценка погрешности .
Вычислительные эксперименты по решению систем дифференциальных уравнений с малой нелинейной частью в форме многочлена шестой степени методом преобразований показывают точность четвертого порядка при вычислении.
Исследование устойчивости решения нелинейных динамических систем
Для автономных систем применимы теоремы
А. М. Ляпунова для решения задачи устойчивости
в линейном приближении [24]. Движение
асимптотически устойчиво, если вещественные части характеристического уравнения отрицательны.
Рассмотрим автономную систему с вещественными коэффициентами (8), полученную в результате многочленных преобразований:
,
где – малая величина ошибки в методе преобразований.
Для преобразованной системы решение определено в окрестности пространства
при .
Исследуем непрерывную неотрицательную
форму , которая оценивается на промежутке сверху:
,
где – ограниченные, непрерывные, неотрицательные возрастающие функции.
Определим производную по времени от V:
(9)
Оценим компоненты в равенстве (9), учитывая :
,
где .
В результате получим дифференциальное неравенство вида
, (10)
где
В случае отсутствия резонансов и неравенство (10) имеет вид
По теореме об устойчивости движения системы [22] пусть в решение V системы (9) удовлетворяет
.
Пусть G определена, однозначна, кусочно-непрерывна и удовлетворяет условию Липшица по и .
Тогда все решения (9) с начальными значениями из при удовлетворяют
при , (11)
где C – решение уравнения сравнения
с начальным условием .
В случае отсутствия резонансов уравнение сравнения имеет вид
(12)
Запишем решение дифференциального уравнения (12):
(13)
Решение определено при выполнении неравенства .
Подставив (13) в (11), получим:
(14)
Таким образом, построена оценка (14) окрестности системы, в которую входит любое решение с начальными значениями из .
Приложение метода преобразования для исследования математической модели виброзащитной системы
Проведено исследование полученной авторами математической модели виброзащитной системы для вышки управления полетами представленным методом. Модель вышки состоит из следующих компонентов: вышка, платформа основания, кабина управления, рабочая площадка, система управ-ления. Исследуем упрощенную модель вышки, включающую кабину управления (m1), рабочую площадку (m2). Она установлена на платформу (m3) для оборудования. Для виброзащиты рассматривается кабина управления. Между каждой массой установлены нелинейные амортизаторы и демпферы. Нелинейные характеристики аппроксимированы полиномами третьей степени. Предположим, что в результате воздействия внешних периодических сил основание совершает малые вертикальные колебания вида . Размерности: Введем обо-значения координат центра масс: z1 – кабины управления, z2 рабочей площадки, z3 – платформы. Запишем приведенную нелинейную модель в форме
где – относительные вертикальные координаты; k, l, p – приведенные к единице массы коэффициенты жесткости; b, c, d – приведенные
к единице массы коэффициенты демпфирования. Здесь размерности для величин уравнений:
Посредством математического пакета программ Mathematica при следующих коэффициентах представленным методом преобразований определен режим колебаний виброзащитной системы:
Обозначим малые колебания основания: .
Методом многочленных преобразований найден установившийся полигармонический режим:
;
.
Решение в абсолютных координатах:
;
;
.
Таким образом, исследованная численно-аналитическим методом преобразований виброзащитная модель эффективно снижает влияние внешних воздействий для вышки управления полетами.
Заключение
В работе решается научная задача разработки высокопроизводительного алгоритма численно-аналитического метода преобразований для исследования режимов эксплуатации нелинейных динамических систем. Нелинейности в исследуемых моделях представлены полиномами. В работе представлен метод преобразований для исследо-вания нелинейных математических моделей. Приведены алгоритмические формулы метода для вычисления. Представлена общая матричная форма для векторных индексов. Для экономично-го вычисления правых частей полиномиальной структуры представлены формулы и предложено применить схему Пана с предварительной обра-боткой коэффициентов. Метод позволяет исследовать различные режимы динамики нелинейных моделей, например, определять такие экстремальные режимы, как резонанс, субгармонический, полигармонический режим. В качестве примера применения метода решена задача виброзащиты вышки управления полетами от внешних периодических воздействий. Решение, построенное методом преобразований, учитывает все нелинейные ком-поненты полиномов. Метод позволяет выполнять исследование динамики широкого круга нелинейных систем с необходимой точностью
1. Mel'nikov V. G. Preobrazovanie dinamicheskih mnogochlennyh sistem s primeneniem approksimacii Chebysheva // Nauch.-tehn. vestn. inform. tehnologiy, mehaniki i optiki. 2012. № 4. S. 85–89.
2. Bodunov N. M. Chislenno-analiticheskoe reshenie prikladnyh zadach mehaniki s primeneniem polinomi-al'nyh bazisnyh funkcii // Yuzh.-Sib. nauch. vestn. 2019. T. 2. № 4 (28). S. 204–208.
3. Naimov A. N., Kobilzoda M. M. O razreshimosti periodicheskoy zadachi dlya nelineynogo obyknovennogo differencial'nogo uravneniya vtorogo poryadka // Izv. vuzov. Matematika. 2021. № 8. S. 56–65.
4. Bulychev Yu. G., Mel'nikov A. V., Kondrashov A. G., Radu P. Yu. Chislenno-analiticheskoe opisanie i ocenivanie parametrov mnogomernogo dinamicheskogo ob'ekta pri nekorrektnyh nablyudeniyah // Radiotehnika. 2019. T. 83. № 10. S. 11–16.
5. Titovcev A. S., Kirpichnikov A. P. Chislenno-analiticheskoe modelirovanie sistem massovogo obsluzhivaniya polikomponentnyh potokov // Matematicheskie metody v tehnike i tehnologiyah: sb. tr. Mezhdunar. nauch. konf.: v 12 t. SPb.: Izd-vo Politehn. un-ta, 2019. T. 2. S. 80–84.
6. Bondarenko A. I., Voronova E. Yu. Primenenie chislenno-analiticheskih metodov dlya postroeniya bazisnyh funkciy konechnyh elementov // Sovrem. priklad. issledovaniya. 2020. S. 322–326.
7. Klimina L. A., Masterova A. A., Samsonov V. A., Selyuckiy Yu. D. Chislenno-analiticheskiy metod poiska avtorotaciy mehanicheskoy sistemy s dvumya vraschatel'nymi stepenyami svobody // Izv. RAN. Mehanika tverdogo tela. 2021. T. 56. № 3. S. 392–403.
8. Aptekarev A. I., Tulyakov D. N. Asimptoticheskiy bazis resheniy q-rekurrentnyh sootnosheniy vne zony blizkih sobstvennyh znacheniy // Prepr. In-ta priklad. matematiki im. M. V. Keldysha RAN. 2018. 24 s.
9. Galanin M. P., Konev S. A. Ob odnom chislennom metode resheniya obyknovennyh differencial'nyh uravneniy // Prepr. In-ta priklad. matematiki im. M. V. Keldysha RAN. 2017. 28 s.
10. Troickaya A. V., Sazonov V. V. Periodicheskie resheniya differencial'nogo uravneniya vtorogo poryadka s bol'shim parametrom // Prepr. In-ta priklad. matematiki im. M. V. Keldysha RAN. 2018. № 71. 16 s.
11. Kartavchenko A. E., Tarhov D. A. Sravnenie metodov postroeniya priblizhennyh analiticheskih resheniy differencial'nyh uravneniy na primere elementarnyh funkciy // Sovrem. informac. tehnologii i IT-obrazovanie. 2017. T. 13. № 3. S. 16–23.
12. Ch'e E. U., Shein A. B. Osobennosti realizacii algoritmov analiza perehodnyh processov v dinamicheskih sistemah operatornym metodom // Informatika i sistemy upravleniya. 2017. № 2. S. 131–135.
13. Bryuno A. D. O slozhnyh razlozheniyah resheniy obyknovennyh differencial'nyh uravneniy // Zhurn. vychisl. matematiki i mat. fiziki. 2018. № 3. S. 346–364.
14. Petrov L. F. Poisk periodicheskih resheniy suschestvenno nelineynyh dinamicheskih sistem // Zhurn. vychisl. matematiki i mat. fiziki. 2018. № 3. S. 403–413.
15. Goryuchkina I. V. Teorema o shodimosti formal'nogo resheniya ODU // Prepr. In-ta priklad. matematiki im. M. V. Keldysha RAN. 2011. 16 s.
16. Arushanyan O. B., Zaletkin S. F. K teorii vychisleniya ortogonal'nogo razlozheniya resheniya zadachi Koshi dlya obyknovennyh differencial'nyh uravneniy vtorogo poryadka // Vychisl. metody i programmirovanie. 2018. T. 19. № 2. S. 178–184.
17. Kudryashov N. A. Metod logisticheskoy funkcii dlya nahozhdeniya analiticheskih resheniy nelineynyh differencial'nyh uravneniy // Modelirovanie i analiz informacionnyh sistem. 2015. T. 22. № 1. S. 23–37.
18. Melnikov G. I., Dudarenko N. A., Malykh K. S., Ivanova L. N., Melnikov V. G. Mathematical Models of Nonlinear Oscillations of Mechanical Systems with Several Degrees of Freedom // Nonlinear Dynamics and Systems Theory. IET – 2017. V. 17. N. 4. P. 369–375.
19. Melnikov V. G., Melnikov G. I., Malykh K. S., Dudarenko N. A. Poincare–Dulac method with Chebyshev economization in autonomous mechanical systems simulation problem // 2015 International Conference on Mechanics – Seventh Polyakhov's Reading (St. Petersburg, 02–06 February 2015 g.). Saint-Petersburg, 2015. P. 7106757.
20. Bulychev Yu. G., Mel'nikov A. V. Chislenno-analiticheskiy metod issledovaniya povedeniya dinami-cheskoy sistemy po rezul'tatam nekorrektnyh nablyudeniy bez rasshireniya prostranstva sostoyaniy // Zhurn. vychisl. matematiki i mat. fiziki. 2019. T. 59. № 6. S. 937–950.
21. Ivanov S. E. Algoritmicheskaya realizaciya metoda issledovaniya nelineynyh dinamicheskih sistem // Nauch.-tehn. vestn. inform. tehnologiy, mehaniki i optiki. 2012. № 4 (80). S. 90–92.
22. Mel'nikov G. I. Dinamika nelineynyh mehani-cheskih i elektromehanicheskih sistem. L.: Mashinostroenie, 1975. 198 s.
23. Pan V. Ya. O sposobah vychisleniya znacheniy mnogochlenov // Uspehi mat. nauk. 1966. T. 21. № 1 (127). S. 103–134.
24. Merkin D. R. Vvedenie v teoriyu ustoychivosti dvizheniya. M.: Nauka, 1987. 304 c.
