Methodology for managing the flows of target information in the remote sensing space system Part 2. Interrelated mathematical models systems formation
- Authors: Starkov A.V.1, Emelyanov A.А.2, Grishantseva L.A.2, Zhukovskaya K.I.2, Morozov A.A.2, Trishin A.A.1
-
Affiliations:
- Moscow Aviation Institute (National Research University)
- Russian Space Systems, Research Center for Earth Operative Monitoring
- Issue: Vol 22, No 2 (2021)
- Pages: 148-161
- Section: Articles
- URL: https://journals.rudn.ru/engineering-researches/article/view/27547
- DOI: https://doi.org/10.22363/2312-8143-2021-22-2-148-161
- ID: 27547
Cite item
Full Text
Abstract
In the second part of the series of articles, the issues of the systemic organization of mathematical models for solving the problem of controlling the flows of target information in the Earth remote sensing space system are considered. A description of the interrelated mathematical models of the orbital constellation as components of the information system, the main task of which is to survey ground objects and the formation of the initial volume of information for its further processing, is presented. To calculate the time of servicing the request by the space segment, the following methods of formation are proposed: a model of the evolution of the Earth remote sensing spacecraft orbit; model for forecasting possible spacecraft correction intervals to maintain nominal orbital parameters; model for forecasting possible time intervals for on / off cycles of observation equipment; model for forecasting possible time intervals for dumping the received information to the information reception points. When calculating the cost of servicing a single request from the orbital complex, both the cost of servicing one spacecraft per unit of time and the cost of processing a single request from the ground complex were taken into account. In conclusion, a generalized form of representation of the target information flow model of the Earth remote sensing space system is proposed as an interconnected sequence of functions for changing the amount of information when an appropriate processing process (traffic change functions) is applied to it. General approaches to solving the optimization problem are considered.
Full Text
Введение В первой части данной серии статей были представлены результаты формализации актуальной задачи управления потоками целевой информации при функционировании космической системы дистанционного зондирования Земли (ДЗЗ), сформулирована техническая постановка задачи и предложен общий методический подход к решению поставленной технической задачи управления, который включает в себя: формулировку постановки оптимизационной задачи, выбор общего подхода к ее решению и определение условий применимости [1-9]. Математическая постановка задачи заключается в оптимизации маршрута прохождения заявки по узлам графа, в котором вершиной является заявка, а формализованные данные заявки передаются последовательно на элементы космического и наземного сегментов: с учетом возможности вычисления представленных выше системных показателей, или их свертки, при наличии текущих ограничений Характеристики узлов и каналов связи между ними считаются известными и по ним мы можем определить время и стоимость передачи информации. Прохождение трафика по графу приводит к вычислению системных показателей. Представленная математическая постановка определяет необходимость в разработке согласованных между собой математических моделей более низкого уровня [2-4], декомпозиция которых изображена на рис. 1. 1. Общие подходы к формированию отдельных элементов системы взаимосвязанных моделей К ним относятся (см. рис. 1): модель движения в задаче определения и прогнозирования состояния космического аппарата (КА), а также проведения орбитальных коррекций; модели описания земных объектов модель проведения съемки, включая модели целевой аппаратуры; модели проведения сеансов связи и передачи целевой информации между КА и пунктом приема и передачи информации (ППИ) [2]. Результатом применения моделей являются: - множество интервалов времени вывода КА на техобслуживание, связанное с необходимостью проведения коррекций орбиты [10-18]; - последовательность возможных интервалов времени проведения наблюдений заданного района - последовательность возможных интервалов времени проведения сеансов связи КА ↔ ППИ В качестве модели движения можно принять систему дифференциальных уравнений, описывающих пространственное движение космического аппарата в оскулирующих элементах [10-15]: Рис. 1. Система взаимосвязанных моделей Figure 1. System of interconnected models с учетом влияния нецентральности гравитационного поля Земли, возмущений гравитационной природы от Луны и Солнца, давления солнечного света, торможения атмосферой и работы корректирующей двигательной установки. Основными возмущениями типовых орбит КА орбитальной группировки (ОГ) ДЗЗ являются долгопериодические возмущения гринвичской долготы восходящего узла трасс, вековые возмущения эксцентриситета, аргумента перигея, наклонения и абсолютной долготы восходящего узла. Для типовых орбит эти задачи в основном сводятся [10-18] к формированию (поддержанию) гринвичских долгот восходящих узлов орбит в заданном ограниченном диапазоне, который определяет возможные колебания трасс полета. Поэтому при построении космических систем необходимо решать задачу приведения и стабилизации трасс полета КА. Классическое решение предполагает проведение орбитальных коррекций с целью изменения двух параметров - высоты перигея орбиты и периода обращения КА [10-18]. Задача об изменении высоты перигея формулируется следующим образом. Задана исходная орбита. Полагается, что она определяется апогейным и перигейным расстояниями соответственно . Требуется определить величину, направление и точку приложения управляющего импульса скорости, обеспечивающего изменение перигейного расстояния с на при условии сохранения апогейного расстояния орбиты, т. е. при условии . Оптимальный маневр при изменении большой полуоси орбиты характеризуется известным трансверсальным импульсом, приложенным в перигее орбиты. Стохастический подход к решению задачи управления подразумевает линеаризацию исходных уравнений движения в окрестности средней долготы [13-14]. Это позволяет синтезировать алгоритм управления для проведения динамических операций в окрестности опорной орбиты. Так как главными возмущениями, изменяющими драконический период обращения КА по орбите, являются возмущения вследствие резонанса с долготными членами разложения гравитационного потенциала и доминирующими являются компоненты возмущений с индексами (2, 2, 1, 1) в разложении геопотенциала, то в уравнениях движения можно выделить уравнения, описывающие эволюцию средней долготы и периода обращения (скорости дрейфа). С учетом ошибок реализации управляющего ускорения при корректировании средней долготы и возмущений от коррекции наклонения уравнение движения долготы восходящего узла представимо в дискретном виде (см. [13-14]). Для определения последовательности возможных интервалов времени проведения наблюдений отдельным КА заданного района можно использовать представленную выше модель прогноза с учетом запретов на проведение съемки, связанных с техобслуживанием КА и проведением коррекций. Иллюстрация возможных моделей объектов съемки и проверки возможности наблюдения и сеанса связи приведена на рис. 2-4 [2]. При этом должны учитываться возможности многопроходной съемки, выполняемой в надир, многопроходной съемки, выполняемой с переориентацией оси визирования камеры, съемки в азимутальном режиме, реализуемом за счет синхронизированного программного разворота по тангажу и крену, одновитковой двухполосной съемки со скачкообразной переориентацией оси визирования камеры по тангажу и крену после первого прохода. В модели переориентации КА участкам изменения углов визирования земных объектов бортовой съемочной аппаратурой КА соответствуют участки кусочно-постоянной функции скорости их переориентации, а сам процесс переориентации бортовой съемочной аппаратуры КА характеризуется предельным углом визирования земного объекта и предельной скоростью. Figure 2. Types of shooting objects Рис. 3. Проверка возможности наблюдения Figure 3. Checking the possibility of observation Рис. 4. Проверка возможности проведения сеанса связи КА↔ППИ Figure 4. Checking the possibility of conducting a communication session Spacecraft ↔ information reception points Для расчета ограничений на объем данных отдельно взятого КА можно принять модель динамики бортового запоминающего устройства, в которой занятый объем БЗУ увеличивается/уменьшается линейно по времени, скорость роста зависит от числа одновременно работающих камер, скорость сброса определяется производительностью передающей аппаратуры, а приращение занятого объема бортового запоминающего устройства (БЗУ) за счет съемки одного объекта определяется суммарным объемом снимков одновременно задействованных камер. Как было отмечено ранее обобщенная форма представления модели потока целевой информации космической системы ДЗЗ может быть представлена как взаимосвязанная последовательность функций изменения объема информации при применении к нему соответствующего процесса обработки (функций изменения трафика). Формально каждый элемент рассматриваемой системы представляется как функция преобразования объема входной информации в объем выходной информации инвариантно к ее тематическому содержанию, при одновременном вычислении показателей эффективности: времени, необходимого для выполнения процесса; условной стоимости обработки. И проверке ограничений: на максимальное количество одновременно выполняющихся однотипных процессов; на максимальный общий объем памяти. В рамках такой формализации для космических аппаратов, наземных пунктов и линий необходимо сформировать модель расчета объема выходной информации; модель расчета времени выполнения процесса; модель расчета стоимости. Для оценки объема сбрасываемой с отдельного КА на ППИ информации (начального трафика) можно принять, что для каждого типа аппаратуры наблюдения имеется прямая связь между общей площадью района съемки, заданного в заявке, и объемом информации, передаваемым с каждого конкретного КА на ППИ. Эта связь может быть задана как в виде аналитической функции, так и таблично. Время выполнения процесса определяется как время, необходимое для передачи циклограммы работы на борт КА, для проведения съемки и передачи снимка на ППИ. Ограничением является объем БЗУ. Все эти параметры вычисляются с использованием математического аппарата, методика формирования которого представлена выше. Для оценки затрат на обслуживание единичной заявки со стороны орбитального комплекса можно считать, что стоимость обслуживания одного КА в единицу времени Cуд зависит от сложности КА - QКА и затрат на оплату труда операторов смены. Тогда оценку стоимости обслуживания заявки CiqКА одним КА можно рассчитать по формуле [2-4] где Cоп - стоимость чел.-ч в конкретном ППИ; Nоп - число операторов в смене. Параметр , входящий в эмпирическую формулу, существенно зависит от уровня развития техники и нуждается в периодическом уточнении на основании обработки статистических данных по различным КА. Также возможно более упрощенное, табличное представление уравнения для стоимости обслуживания заявки. Для оценки изменения объема информации при применении к ней процесса обработки будем считать, что имеется прямая связь между объемом входной и выходной информации. Эта связь может быть задана как в виде аналитической функции, так и таблично. Тогда общее время выполнения процесса рассчитывается по формуле [2-4] (1) где Tожидания - время ожидания до момента освобождения ресурсов в следующем узле, когда для каждого момента времени выполнения процесса выполняются условия: количество одновременно запущенных однотипных процессов на каждом наземном пункте строго меньше заданного Vвх+Vвых < Vсвоб; объем входной и выходной информации не превышает уровня объема свободной памяти Nпроц ≤ Nпроц max; Tэталон.проц i - эталонное время выполнения процесса в эталонном наземном пункте, зависящее от объема входной информации; Kппи - коэффициент производительности наземного пункта, в котором выполняется процесс обработки данных. Возможно также более упрощенное представление расчета времени в виде таблицы. Для оценки затрат на обслуживание процесса при обработке единичной заявки со стороны наземного пункта можно принять, что стоимость обслуживания процесса в единицу времени Cпроц.уд зависит от его сложности - Qпроц и затрат на оплату труда задействованных операторов. Тогда для стоимости обслуживания заявки Ciqпроц одним ППИ будем иметь где Cоп - стоимость чел.-ч в конкретном наземном пункте; Nоп - число операторов, задействованных в процессе обработки; Cппи - коэффициент себестоимости наземного пункта, в котором выполняется процесс обработки данных. Возможно более упрощенное представление расчета стоимости в виде таблицы. Представленный подход к формированию математических моделей позволяет определить способы их дальнейшей программной реализации. При этом должна быть учтена необходимость обеспечения возможности оперативного внесения изменений как в форму представления, так и в параметры модели. Механизмы, которые обеспечат такую возможность, формируются одним из ни жеперечисленных способов: 1. Текстовой формулой в синтаксисе языка программирования, например: . В этом случае в параметры правой части формулы помимо VВХОД могут входить как длительность процесса, так и общие константы. 2. В виде полинома, например: , где коэффициенты a, b, c задаются экспертами. 3. Таблично. 4. Внешней вычислительной библиотекой, если расчет показателей связан с привлечением сложной вычислительной процедуры. 2. Описание системы взаимосвязанных математических моделей Модель процесса обработки информации космическим сегментом представляется в виде функции генерации трафика, входом которой являются параметры заявки потребителя. Результатом применения модели является величина информации, которая после проведения сеанса связи попадает для дальнейшей обработки в НКПОР. В процессе вычислений определяются показатели эффективности - время и стоимость [3]. Для расчета времени и стоимости обслуживания заявки космическим сегментом должны быть сформированы: - модель эволюции орбиты КА ДЗЗ; - модель прогноза возможных интервалов коррекции КА для поддержания номинальных параметров орбит (стратегии коррекций); - модель прогноза возможных интервалов времени циклов включения/выключения аппаратуры наблюдения; - модель прогноза возможных интервалов времени для сброса полученной информации на ППИ. Для модели съемки земных объектов должны быть сформированы [2]: - модель съемочного плана, включающая в себя: модели объектов съемки, модели для проверки выполнения условий съемки и модель временных последовательностей выполнения съемочных операций; - модель съемки точечных объектов; - модель съемки площадных объектов, сформированная на основе аппроксимации объекта съемки как в виде замкнутого полигона, так и в виде совокупности точечных объектов. Модель передачи целевой информации между КА и ППИ в общем случае служит для определения интервалов взаимной видимости ППИ→КА и КА→ППИ. Совокупность моделей целевого функционирования аппаратуры бортового информационного комплекса КА составляют модель динамики состояния БЗУ и модели переориентации бортовой съемочной аппаратуры. Для модели изменения памяти в БЗУ можно считать, что в процессе съемки занятый объем БЗУ увеличивается линейно по времени, в процессе сброса - уменьшается также линейно по времени. Скорость роста зависит от числа одновременно работающих камер, скорость сброса определяется производительностью передающей аппаратуры. Приращение занятого объема БЗУ за счет съемки одного объекта определяется суммарным объемом снимков одновременно задействованных камер. При выполнении съемки земных объектов бортовая съемочная аппаратура КА изменяет свою ориентацию. Переориентация аппаратуры совершается равномерно, скорость переориентации представляет собой кусочно-постоянную функцию, участки которой определяются парами последовательных моментов съемки земных объектов. В результате время, необходимое для проведения съемки по заявке, определяется как совокупность суммарного времени проведения съемки, времени, необходимого для переориентации аппаратуры наблюдения, и времени, необходимого для того, чтобы КА Рис. 5. Граф возможных вариантов выполнения заявки Figure 5. Graph of possible options for executing an order оказался в зоне взаимной радиовидимости с одним из ППИ, при условии, что: - КА в это время не совершает корректирующие маневры; - объема памяти БЗУ достаточно для сохранения снимка объекта съемки на протяжении всего интервала наблюдений. В конечном итоге каждый элемент рассматриваемой системы формально представляется как функция преобразования объема входной информации в объем выходной информации инвариантно к ее тематическому содержанию, при одновременном вычислении показателей эффективности: - времени, необходимого для выполнения процесса; - условной стоимости обработки и проверке ограничений: - на максимальное количество одновременно выполняющихся однотипных процессов; - на максимальный общий объем памяти. Таким образом, возникает полностью определенная задача оптимизации распределения процессов обработки с учетом возможности вычисления представленных выше системных показателей или их свертки, при наличии текущих ограничений. Эта задача может быть интерпретирована как задача оптимизации маршрута на взвешенном графе (рис. 5) [4]. 3. Общий подход к решению оптимизационной задачи Как следует из постановки, задача оптимального распределения целевой информации с разнотипных КА ДЗЗ по составным частям системы является двухкритериальной. В связи с тем, что число критериев невелико, для ее решения можно воспользоваться классическим математическим аппаратом многокритериальной оптимизации [19-20], предварительно адаптировав их к условиям функционирования. Это позволит определить и выработать наиболее предпочтительные варианты распределения целевой информации с разнотипных КА ДЗЗ по составным частям системы. Таким образом, целью является принятие решений в условиях многокритериальной неопределенности, позволяющей - рассматривать множество всех допустимых способов учета неопределенности; - формировать комплексную оценку решений в условиях неопределенности, используя любые известные и неизвестные свертки; - рассматривать как количественные, так и качественные критерии; - использовать интервальные оценки частных критериев; - не задавать весовые коэффициенты (они считаются неопределенными); - устанавливать предпочтения между отдельными критериями. В случае, если число процессов и доступных ППИ для обработки заявки невелико, задача оптимального распределения может быть решена методом перебора. Кроме того, возможно применение одного из классических методов: - метод сведения к однокритериальной задаче с использованием линейной свертки критериев; - метод поиска Парето-оптимальных решений; - метод последовательных назначений; - метод случайной оптимизации. В нашем случае, когда число критериев минимально и равном двум, по-видимому, наиболее эффективным вариантом служит представление всего множества Парето в виде графика, что позволяет проанализировать последствия всех возможных эффективных решений. Это не обеспечивает полного решения задачи, так как по-прежнему остается множество альтернатив, для которых необходимо выбрать единственное наилучшее. Один из современных путей развития Парето-оптимальных методов многокритериальной оптимизации состоит в том, чтобы предоставить возможность опираться на все множество допустимых способов учета неопределенности, что, однако, требует больших вычислительных ресурсов. Такой подход может быть реализован с использованием метода уверенных суждений [19-20]. Изучение частных методов и алгоритмов решения задачи позволяет выработать единый методический подход к их решению. Во-первых, можно заметить, что предлагаемые для решения численные методы и алгоритмы являются стандартными при минимальной адаптации к рассматриваемой предметной области, что позволяет достаточно легко выполнить их программную реализацию. Во-вторых, предложенная формализация моделей процессов обработки заявки также иллюстрирует общность в способах их формирования и программной реализации. Наконец, в-третьих, реализация последовательности «модель-алгоритм-оптимизация» также свидетельствует о возможности выработать единый методический подход, в основе которого предлагается последовательное выполнение четырех основных этапов вычислений (расчетов), укрупненная блок-схема которых представлена на рис. 6. Первый этап - моделирование динамики орбитальной группировки КА ДЗЗ. Целью является исследование взаимного расположения космических аппаратов и ППИ на всем интервале планирования. Определяется множество интервалов времени, доступных для проведения сеансов съемок и приема-передачи информации с учетом эволюции орбит и необходимости вывода КА на техническое обслуживание. Рассчитываются параметры для управления бортовой целевой аппаратурой. Рис. 6. Схема формирования и использования математических моделей Figure 6. Scheme of the formation and use of mathematical models Исходными данными являются: - параметры орбиты КА; - характеристики целевой аппаратуры; - циклограммы работы отдельно взятых КА ОГ ДЗЗ; - информация о НКПОР (географическое расположение, зоны радиовидимости ППИ и др.). К внешним условиям относятся: - условия освещенности и наличие теней в районе съемки; - облачность в районе съемки; - совместимость радиочастотных диапазонов приема-передачи информации и др. Необходимо отметить, что для проведения исследований необходимо использовать: - упрощенные модели, которые позволяют интегрировать уравнения движения как без учета внешних управляющих воздействий, так и при использовании наиболее значимых; - полные модели, учитывающие весь спектр возмущений, для расчетов орбитального движения КА на длительных интервалах времени вплоть до периода их активного существования. Результатом первого этапа являются: - возможные интервалы времени проведения съемки и их продолжительность; - возможные интервалы проведения сеансов связи с наземными ППИ, их продолжительность и объем информации; - условная стоимость реализации процесса съемки как параметров заявки потребителя. Второй этап - формирование модели КА и элементов НКПОР как взаимосвязанных составных частей функционирования космической системы ДЗЗ, представляющей собой модель планирования с учетом специфики сценария ее целевого функционирования и постановки оптимизационной задачи. На этом этапе используются и обобщаются результаты выполнения предыдущего этапа в виде параметров математической модели. Такая модель служит основой для реализации процедуры формирования (поиска) плана целевого функционирования. Элементы модели, по существу, выражают значения показателей эффективности выполнения целевых операций КА системы (см. п. 1.4), а также возможность или отсутствие возможности потенциальных переходов между двумя последовательными процессами с учетом имеющихся ограничений. Третий этап - решение задачи оптимизации распределения целевой информации по элементам (узлам) системы. Проводится синтез временной последовательности выполнения операций на основе модели первого и второго этапа. Следует отметить, что на данном этапе допускается возможность существования нескольких (или множества) планов, что особенно характерно для задач многокритериальной оптимизации. Реализовать процедуру формирования оптимального плана можно с применением одного из классических, например, рассмотренных выше методов либо специально разработанных методов оптимизации. Результатом выполнения рассматриваемого этапа методического подхода является один (или набор) оптимальных (близких к оптимальному) планов распределения целевой информации по узлам системы. Четвертый этап - многократное решение задачи распределения потоков информации в целях проведения системных исследований. На этом этапе осуществляется многократное решение задачи управления при различных вариантах исходных данных, в качестве которых могут выступать параметры орбит КА ОГ ДЗЗ, характеристики бортовой аппаратуры КА, параметры и расположение элементов НКПОР, требования потребителей к результатам тематической обработки и т. д. Результаты многократного решения задачи позволяют выявить зависимости показателей эффективности целевого функционирования при реализации предложенного программно-математического подхода, т. е. провести параметрический анализ целевого функционирования системы в целом. Заключение В статье на основе математической постановки оптимизационной задачи определен перечень взаимосвязанных математических моделей, требующих дальнейшей разработки. В совокупности эти модели описывают процесс обработки, в том числе определяют, что происходит с объемом информации, и вычисляют системные показатели эффективности - время, требуемое для реализации процесса и его условную стоимость. Это позволяет окончательно представить модель обработки информации как взаимосвязанную последовательность функций изменения условного объема информации при применении к нему соответствующего процесса обработки (функция изменения трафика), а сама математическая постановка задачи может быть сведена к оптимизации маршрута прохождения заявки по графу. Таким образом, можно сказать, что предлагаемый подход к формированию и использованию математических моделей при решении прикладных задач оптимизации распределения потоков целевой информации является системным, а для доведения ее до практического использования необходимо формирование соответствующих численных методов, алгоритмов и программного обеспечения.
About the authors
Alexander V. Starkov
Moscow Aviation Institute (National Research University)
Author for correspondence.
Email: starkov@goldstar.ru
SPIN-code: 5242-3413
Professor of the Department of System Analysis and Management, Doctor of Technical Sciences
4 Volokolamskoe Shosse, 125993, Moscow, Russian FederationAndrey А. Emelyanov
Russian Space Systems, Research Center for Earth Operative Monitoring
Email: ntsomz@ntsomz.ru
SPIN-code: 4484-1479
Head of the Research Center for Earth Operative Monitoring, Russian Space Systems, Candidate of Technical Sciences
b. 51, h.25, Decabristov St, 127490, Moscow, Russian FederationLyubov A. Grishantseva
Russian Space Systems, Research Center for Earth Operative Monitoring
Email: grishantseva_la@ntsomz.ru
SPIN-code: 9940-8756
Head of the Sector of the Research Center for Earth Operative Monitoring, Russian Space Systems, Candidate of Physical and Mathematical Sciences
b. 51, h.25, Decabristov St, 127490, Moscow, Russian FederationKsenia I. Zhukovskaya
Russian Space Systems, Research Center for Earth Operative Monitoring
Email: zubkova.k@ntsomz.ru
SPIN-code: 4805-5960
Research Engineer of the 1St Category of the Research Center for Earth Operative Monitoring
b. 51, h.25, Decabristov St, 127490, Moscow, Russian FederationAlexander A. Morozov
Russian Space Systems, Research Center for Earth Operative Monitoring
Email: aamorozko@mail.ru
Research Engineer of the 3rd category of the Research Center for Earth Operative Monitoring
b. 51, h.25, Decabristov St, 127490, Moscow, Russian FederationAlexey A. Trishin
Moscow Aviation Institute (National Research University)
Email: trishin0202@mail.ru
Student of the Department of Information and Control Systems of Aircraft
4 Volokolamskoe Shosse, 125993, Moscow, Russian FederationReferences
- Selin VA, Emelyanov AA, Sizov OS, Emelyanov KS, Borisov AV. Medium-resolution optical range space images: consumer expectations. Izvestiya. Atmospheric and Oceanic Physics. 2020;56(9):1182—1189. doi: 10.1134/ S0001433820090224
- Emel’yanov AA, Malyshev VV, Smol’yaninov YuA, Starkov AV. Upravlenie potokami celevoj informacii pri funkcionirovanii kosmicheskoj sistemy distancionnogo zondirovaniya Zemli [Control of target information flows in the functioning of the space system of remote sensing of the Earth: Monograph]. Мoscow: MAI-PRINT Publ.; 2020. (In Russ.)
- Emelyanov AA, Malyshev VV, Nguyen VHN, Starkov AV. Mathematical model of functioning of the ground segment in distributed remote sensing data processing. Nauchnotekhnicheskij vestnik Povolzh’ya [Scientific and Technical Bulletin of the Volga Region]. 2018;2:74—79. (In Russ.)
- Еmelianov АА, Grishantseva LA, Zubkova KI, Malyshev VV, Nguyen VHN, Starkov АV, Zay YW. Mathematical model of ERS data processing ground segment operation in terms of processing distribution. Advances in the Astronautical Sciences (vol. 170). CA, USA: Univelt Inc.; 2020. p. 495—504.
- Leun EV, Leun VI, Sysoev VK, Zanin KA, Shulepov AV, Vyatlev PA. The active control devices of the size of products based on sapphire measuring tips with three degrees of freedom. IOP Conf. Series: Journal of Physics: Conf. Series. 2018;944:012073. doi: 10.1088/1742-6596/944/1/012073
- Zanin KA, Moskatinev IV. Improvement of Methods for Evaluating the Resolving Power of a Space Synthetic Aperture Radar. Solar System Research. 2018;52(7):666-672. doi: 10.1134/S0038094618070213
- Zay Yar Win, Malyshev VV, Bobronnikov VT, Starkov AV. The joint solution of problem of evasion and keeping in a neighborhood reference orbit. Advances in the Astronautical Sciences. CA, USA: Univelt Inc. 2020;170:433—442.
- Malyshev VV, Starkov AV, Fedorov AV. Orbital Corrections of Space Vehicles while Performing Dynamic Operations. J. Comput. Syst. Sci. Int. 2013.52(2):313—325. doi: 10.1134/S1064230713010085
- Malyshev VV, Starkov AV, Zay Yar Win. The Decision of Problems of Evasion When Holding the Geostationary Satellites in the Neighborhood of The Reference Orbit. Journal of Advanced Research in Dynamical and Control Systems. 13-Special Issue, 2018;10:53—58.
- Razoumny Y, Razoumny V, Kozlov P, Baranov A, Varatharajoo R. Method of optimization of the servicing spacebased system orbits and detached units maneuveres parameters in the problem of on-orbit-servicing of the given multisatellite space infrastructure. Proceedings of the International Astronautical Congress, IAC. 2016.
- Baranov AA, Razoumny VY, Razoumny YN, Malyshev VV. Low orbit spacecraft service planning. Proceedings of the International Astronautical Congress, IAC68. 2017:835—844.
- Baranov AA, Grishko DA, Mayorova VI. The features of constellations’ formation and replenishment at near circular orbits in non-central gravity fields. Acta Astronautica. November–December 2015;116:307—317.
- Krasil’shchikov MN, Malyshev VV, Fedorov AV. Autonomous implementation of dynamic operations in a geostationary orbit. I. Formalization of control problem. J. Comput. Syst. Sci. Int. 2015;54(6):916—930. doi: 10.1134/ S1064230715060118
- Voiskovskii AP, Krasil’shchikov MN, Malyshev VV, Fedorov AV. Autonomous implementation of dynamic operations in a geostationary orbit. II. Synthesis of control algorithms. J. Comput. Syst. Sci. Int. 2016;55(6):948—968. doi: 10.1134/ S1064230716060113
- Petukhov VG. Application of the Angular Independent Variable and Its Regularizing Transformation in the Problems of Optimizing Low-Thrust Trajectories. Cosmic Research. 2019; 57(5);351—363. doi: 10.1134/S001095251905006X
- Ivanyukhin AV, Petukhov VG. Low-Energy SubOptimal Low-Thrust Trajectories to Libration Points and HaloOrbits. Cosmic Research. 2019;57(5):378-388. doi: 10.1134/ S0010952519050022
- Petukhov VG, Ivanyukhin AV, Sang Wook W. Joint Optimization of Control and Main Trajectory and Design Parameters of an Interplanetary Spacecraft with an Electric Propulsion System. Cosmic Research. 2019;57(3):188-203. doi: 10.1134/S0010952519030079
- Grechkoseev AK, Krasil’shchikov MN, Kruzhkov DM, Mararescul TA. Refining the Earth Orientation Parameters Onboard Spacecraft: Concept and Information Technologies. J. Comput. Syst. Sci. Int. 2020;59(4):598—608. doi: 10.1134/ S1064230720040061
- Golubev SI, Malyshev VV, Piyavskii SA, Sypalo KI. Decision making in multicriteria problems at the image design stage of aviation rocket technique. J. Comput. Syst. Sci. Int. J. Comput. Syst. Sci. Int. 2020;59(1):83—94. doi: 10.1134/ 2020;59(2):223—231. doi: 10.1134/ S1064230720020057
- Brusov VS, Korchagin PO, Malyshev VV, Piyavsky SA. Advanced «Confident Judgments» Method when Choosing Multicriteria Solutions in a Multipurpose Approach. J. Comput. Syst. Sci. Int. 2020. Vol. 59 (1). P. 83-94. doi: 10.1134/S1064230720010049
Supplementary files










