Национальный цифровой ресурс Руконт - межотраслевая электронная библиотека (ЭБС) на базе технологии Контекстум (всего произведений: 634938)
Контекстум
Руконтекст антиплагиат система
Ракетно-космическое приборостроение и информационные системы  / №3 2016

Использование методов линейного программирования в задаче планирования сеансов приема целевой информации с КА орбитальной группировки ДЗЗ

0   0
Первый авторЧернов
Страниц8
ID577058
АннотацияРассмотрена постановка и решение задачи оптимального распределения (планирования) сеансов приема целевой информации (ЦИ), получаемой от КА орбитальной группировки (ОГ) ДЗЗ, по ограниченному количеству приемных комплексов (ПК) НКПОР, с учетом приоритетности ЦИ и возможности одновременного попадания в зону радиовидимости пункта приема информации (ЗРВ ППИ) нескольких КА. Решение рассматриваемой задачи производится в два этапа. На первом этапе решается задача формирования критерия эффективности планирования расписания сеансов приема ЦИ. Критерий эффективности формируется на основе «весовой матрицы» («матрицы приоритетов»), определяющей относительную ценность сеансов приема ЦИ на прямом произведении множеств всех КА ОГ и всех ПК НКПОР. Алгоритм расчета «весовой матрицы» позволяет учесть некоторые априорные требования Оператора ОГ, которые он предъявляет к конечному результату планирования. На втором этапе построенный критерий эффективности используется для постановки и решения собственно задачи оптимального планирования расписания сеансов приема [2]. Задачи и первого и второго этапов ставятся как транспортные задачи целочисленного линейного программирования. Ограничения и особенности задачи планирования расписания сеансов приема позволяют поставить ее как известную в целочисленном линейном программировании «задачу о назначениях». Приведены примеры расчета «весовых матриц» и оптимального расписания сеансов приема, полученные с помощью разработанного программного комплекса.
УДК629.7.054:621.396
Чернов, А.А. Использование методов линейного программирования в задаче планирования сеансов приема целевой информации с КА орбитальной группировки ДЗЗ / А.А. Чернов // Ракетно-космическое приборостроение и информационные системы .— 2016 .— №3 .— С. 46-53 .— doi: 10.17238/issn2409-0239.2016.3.46 .— URL: https://rucont.ru/efd/577058 (дата обращения: 02.05.2024)

Предпросмотр (выдержки из произведения)

РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ 2016, том 3, выпуск 3, c. <...> 46–53 АЭРОКОСМИЧЕСКИЕ МЕТОДЫ ЗОНДИРОВАНИЯ ЗЕМЛИ УДК 629.7.054:621.396 Использование методов линейного программирования в задаче планирования сеансов приема целевой информации с КА орбитальной группировки ДЗЗ А. А. <...> Рассмотрена постановка и решение задачи оптимального распределения (планирования) сеансов приема целевой информации (ЦИ), получаемой от КАорбитальной группировки (ОГ) ДЗЗ, по ограниченному количеству приемных комплексов (ПК) НКПОР, с учетом приоритетности ЦИ и возможности одновременного попадания в зону радиовидимости пункта приема информации (ЗРВ ППИ) нескольких КА. <...> На первом этапе решается задача формирования критерия эффективности планирования расписания сеансов приема ЦИ. <...> Критерий эффективности формируется на основе «весовой матрицы» («матрицы приоритетов»), определяющей относительную ценность сеансов приема ЦИ на прямом произведении множеств всех КАОГ и всех ПК НКПОР. <...> Алгоритм расчета «весовой матрицы» позволяет учесть некоторые априорные требования Оператора ОГ, которые он предъявляет к конечному результату планирования. <...> На втором этапе построенный критерий эффективности используется для постановки и решения собственно задачи оптимального планирования расписания сеансов приема [2]. <...> Задачи и первого и второго этапов ставятся как транспортные задачи целочисленного линейного программирования. <...> Ограничения и особенности задачи планирования расписания сеансов приема позволяют поставить ее как известную в целочисленном линейном программировании «задачу о назначениях». <...> Приведены примеры расчета «весовых матриц» и оптимального расписания сеансов приема, полученные с помощью разработанного программного комплекса. <...> Ключевые слова: орбитальная группировка, оптимальное планирование сеансов приема информации, критерий оптимальности, приоритетность КА Application of Linear Programming Methods in Scheduling of <...>
Использование_методов_линейного_программирования_в_задаче_планирования_сеансов_приема_целевой_информации_с_КА_орбитальной_группировки_ДЗЗ.pdf
РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ 2016, том 3, выпуск 3, c. 46–53 АЭРОКОСМИЧЕСКИЕ МЕТОДЫ ЗОНДИРОВАНИЯ ЗЕМЛИ УДК 629.7.054:621.396 Использование методов линейного программирования в задаче планирования сеансов приема целевой информации с КА орбитальной группировки ДЗЗ А. А.Чернов к. т. н., АО «Российские космические системы» e-mail: a2chernov@yandex.ru Аннотация. Рассмотрена постановка и решение задачи оптимального распределения (планирования) сеансов приема целевой информации (ЦИ), получаемой от КАорбитальной группировки (ОГ) ДЗЗ, по ограниченному количеству приемных комплексов (ПК) НКПОР, с учетом приоритетности ЦИ и возможности одновременного попадания в зону радиовидимости пункта приема информации (ЗРВ ППИ) нескольких КА. Решение рассматриваемой задачи производится в два этапа. На первом этапе решается задача формирования критерия эффективности планирования расписания сеансов приема ЦИ. Критерий эффективности формируется на основе «весовой матрицы» («матрицы приоритетов»), определяющей относительную ценность сеансов приема ЦИ на прямом произведении множеств всех КАОГ и всех ПК НКПОР. Алгоритм расчета «весовой матрицы» позволяет учесть некоторые априорные требования Оператора ОГ, которые он предъявляет к конечному результату планирования. На втором этапе построенный критерий эффективности используется для постановки и решения собственно задачи оптимального планирования расписания сеансов приема [2]. Задачи и первого и второго этапов ставятся как транспортные задачи целочисленного линейного программирования. Ограничения и особенности задачи планирования расписания сеансов приема позволяют поставить ее как известную в целочисленном линейном программировании «задачу о назначениях». Приведены примеры расчета «весовых матриц» и оптимального расписания сеансов приема, полученные с помощью разработанного программного комплекса. Ключевые слова: орбитальная группировка, оптимальное планирование сеансов приема информации, критерий оптимальности, приоритетность КА Application of Linear Programming Methods in Scheduling of Reception Sessions of the Target Information from SC of ERS Orbital Constellation А. А.Chernov candidate of engineering science, Joint Stock Company “Russian Space Systems” e-mail: a2chernov@yandex.ru Abstract. This article considers the problem of optimal allocation (scheduling) of target information (TI) reception sessions from the spacecraft of the ERS orbital constellation (OC) with a limited number of core ground segment (CGS) reception complexes (RC), and suggests its solution, taking into account the priorities of the TI and the possibility of simultaneous presence of several SC in the radio visibility zone (RVZ) of the receiving station (RS). The solution of the problem is carried out in two stages. The first step is to work out a criterion of efficiency of scheduling of the TI reception sessions. The criterion of effectiveness is based on the “weight matrix” (“priority matrix”), which determines the relative value of the TI reception sessions on the direct product of factors of the orbital constellation SC and all the CGS RC. The algorithm for the “weight matrix” calculation accommodates some a priori requirements of the OC operator, which he specifies for the end result of the scheduling. On the second stage the criterion of efficiency is used to set and solve the problem of optimal scheduling of the reception sessions [2]. The problems of the first and second stages are viewed as transportation problems of integer linear programming. Limitations and features of the problem of reception session scheduling make it possible to view it as the “assignment problem”, well known in integer linear programming. Examples of “weight matrix” calculation and optimal scheduling of reception sessions using the developed software are presented. Keywords: orbital constellation, optimal scheduling of information reception sessions, criterion of efficiency, SC priority
Стр.1
ИСПОЛЬЗОВАНИЕ МЕТОДОВ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ Введение Одной из особенностей задачи планирования приема ЦИ является учет возможности одновременного попадания в ЗРВ ППИ (зону радиовидимости пунктов приема информации) нескольких КА, в результате чего прием ЦИ с некоторых КА в это время становится невозможным. Для иллюстрации этого положения на рис. 1 приведены результаты моделирования прохождений КАОГ в ЗРВ ППИ «Отрадное» (ППИ «102») 05.08.2015 11:16 ДМВ (676 мин). Орбитальная группировка включает 9 КА, имена которых приведены на рис. 1. Из рис. 1 следует, что 05.08.2015 11:16 ДМВ в зоне радиовидимости ППИ «Отрадное» одновременно находилось 5 КА: «Метеор-М», «МетеорМ2», «Канопус-В1», «Ресурс-П2», TERRA. На рис. 2 показано развитие ситуации с количеством КАв ЗРВ ППИ в течение указанных суток. В табл. 1 приведены итоговые результаты моделирования прохождений КАОГ в ЗРВ ППИ за 200 сут 2015 г. Таблица 1. nmax в течение суток одновременно попадающих в ЗРВ ППИ «102», и количество суток (в процентах от времени работы ОГ), в которые это событие происходит KA — максимальное количество КА, nmax KA 2 3 4 5 Количество сут, % 19 70 10 1 Из табл.1 следует, чтовкаждыесутки одновременное попадание в ЗРВ ППИ нескольких КА орбитальной группировки надо скорее считать правилом, чем исключением. Более того, для рассматриваемой ОГ из 9 КАсобытие nmax ся примерно в 3,7 раза чаще, чем событие nmax Следовательно, при функционировании данной ОГ неизбежно будут возникать потери целевой информации, если количество ПК на ППИ будет менее nmax KA = 3 реализуетKA = 2. KA = 3. К другим особенностям, которые должны учитываться при постановке задачи, необходимо отнести требования к оперативности получения принимаемой ЦИ; различия в технических характеристиках и составе ПК ППИ; использование на КАцелевой аппаратуры различного разрешения; существование преимуществ для некоторых КАна некото47 рых ПК в определенные периоды времени в обслуживании перед другими КА, а также ряд других особенностей технического характера. Таким образом, сеансы приема ЦИ от каждого КАОГ на каждом ПК НКПОР в условиях возможной потери ЦИ, естественно, приобретают большую или меньшую значимость (приоритетность) по отношению друг к другу и возникает необходимость ранжировать КАотносительно операции приема ЦИ, т. е. для всех КАна всех ПК установить определенные приоритеты pij по приему ЦИ (здесь i = 1, ... ,M —номер КА; j = 1, ... ... ,N — номер ПК), с учетом которых и строить расписание сеансов приема. Приоритеты pij объединяются в матрицу приоритетов W размерности M×N; pij — целые положительные числа, наивысший приоритет равен 1 (как в спорте — первое место), самый низкий (когда i-й КАне обслуживается на j-м ПК) равен большому числу. W =         ... .. . ... ... pM1 pM1 ... pMN p11 p12 p21 p21 ... p1N ... p2N .         Организованная таким образом матрица приоритетов W может быть использована для построения критерия оптимальности F в задаче оптимального планирования расписания сеансов приема ЦИ с КАОГ. Однако отсутствие регулярного алгоритма для получения матрицы W, которая должна обладать определенными свойствами, обусловленными спецификой задачи, вызывает значительные трудности. Формализуем алгоритм построения матрицы W. Установление отношений порядка на множествах КА ОГ иПКНКПОР иполучение матрицы приоритетов W Положим, что оператор при получении матрицы W способен сравнить предпочтительность операции приема ЦИ с каждого КАгруппировки на некотором (одном!) ПК НКПОР и упорядочить множество КАна этом ПК по степени предпочтения, иначе говоря, установить очередность РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 3 вып. 3 2016
Стр.2
48 А.А.ЧЕРНОВ Рис. 1. Прохождения КАОГ в ЗРВ ППИ «Отрадное» Рис. 2. Количество КАОГ, одновременно попадающих в ЗРВ ППИ «Отрадное» в течение суток 05.08.2015 обслуживания каждого КА ОГ на заданном ПК. Отметим, что от оператора не требуется, чтобы он соизмерял «предпочтительность» в каких-то единицах, достаточно, чтобы он был способен ранжировать все КА-группировки по степени предпочтительности их обслуживания на данном ПК. Затем такое установление отношений предпочтения (т. е. очередности обслуживания) среди КАОГ оператор последовательно производит для каждого ПК, задействованного для обслуживания ОГ, составляя из этих отношений столбцы матрицы C (column — столбец) (табл. 2). Если некоторый ПК не задействуется для обслуживания некоторого КА, то соответствующий элемент матрицы C (на пересечении строки «КА» и столбца «ПК») полагается равным большому числу, например 100. Условимся также, что очередность обслуживания строгая, т. е. два КАне могут иметь один и тот же номер в очереди на обслуживание на данном ПК. В результате матрица C приобретает вид, подобный приведенному в табл. 2 (числа 100, которые должны стоять в свободных элементах матрицы, не показаны, чтобы упростить зрительное восприятие данных). Подчеркнем, что информация одного столбца никак не связана с информацией другого: столбцы заполняются совершенно независимо друг от друга. На следующем этапе работы оператор сравнивает предпочтительность операции приема ЦИ РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 3 вып. 3 2016
Стр.3
ИСПОЛЬЗОВАНИЕ МЕТОДОВ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ Та б лица 2. Приоритетность сеансов приема ЦИ с каждого из КАОГ на конкретном ПК НКПОР (в столбцах матрицы C) 49 Та б лица 3. Приоритетность сеансов приема ЦИ с конкретного КАОГ на каждом ПК НКПОР (в строках матрицы R) снекоторогоКАгруппировкинавсехПКНКПОР и упорядочивает множество ПК для рассматриваемого КАпо степени предпочтения, иначе говоря, устанавливает очередь, в которой рассматриваемый КА «выбирал бы» ПК для приема передаваемой ЦИ. Все, что сказано об установлении отношений предпочтения в столбцах матрицы C,действительно теперь по отношению к строкам матрицы R (row — строка). В результате матрица R приобретает вид, подобный приведенному в табл. 3. Подчеркнем, что информация в одной строке никак не связана с информацией в другой: строки заполняются совершенно независимо друг от друга. Понятно, что цифра 100 стоит на одних и тех же местах как в матрице C,такивматрице R, поскольку смысл элементов, равных 100, в обеих матрицах одинаков: j-й ПК с i-го КАинформацию не принимает. Итак, согласно табл. 2 (т. е. матрице C)ранжируются сеансы приема ЦИ с каждого из КАОГ на конкретном ПКj, j = 1, 2, . . . , 7, а согласно табл.3 (т. е.матрице R)—сеансы приемаЦИ на каждом ПК НКПОР с конкретного КАi, i = = 1, 2, . . . , 9. Важно, что опытный оператор устанавливает эти отношения предпочтения, а значит, строит матрицы C и R, без особого труда. Однако использовать любую из этих матриц непосредственно при формировании критерия эффективности оптимальной задачи как «весовую» матрицуW нельзя, так как последняя применяется в алгоритмах оптимизации при одновременном варьировании как КА, так и ПК. А матрицы C и R этого делать не позволяют: матрица C предоставляет корректную информацию для сравнения при переборе ее элементов только в столбцах, а матрица R —только в строках. Более того, матрицы C и R, вообще говоря, могут противоречить друг другу. Поэтому для формирования критерия эффективности задачи оптимального планирования сеансов приема ЦИ необходимо найти матрицу W = f(C,R), удовлетворяющую (если это возможно) одновременно как матрице C, так и матрице R.Будем говорить в этом случае, что матрица W есть результат «слияния» матриц C и R, и используем обозначение W = C R. Представим способ построения матрицы приоритетности W. Для упрощения изложения, но без потери общности, способ получения матрицы W проиллюстрируем на конкретном примере. Пусть матрицы C и R с положительными элементами имеют размерность 4×3(M = 4 — четыре КАи N =3— три ПК). Зададим C и R, а также искомую матрицу W той же размерности в виде C =         14 3 21 4 32 1 43 2         ; R =         132 312 231 321         ; РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 3 вып. 3 2016
Стр.4
50 W = C R =         x11 x12 x13 x21 x22 x23 x31 x32 x33 x41 x42 x43         А.А.ЧЕРНОВ . (1) Для нахождения матрицы W потребуем, чтобы элементы в ее столбцах удовлетворяли тем же отношениям порядка, что и элементы в соответствующих столбцах матрицы C, а также чтобы наименьший по величине элемент в каждом столбце матрицы W был не меньше единицы. В каждом столбце матрицы C сравним очередность обслуживания КАна ПК, записывая неравенства для пар КАв порядке возрастания очередности их обслуживания. Тогда, если матрица W существует, элементы в первом ее столбце должны удовлетворять системе неравенств:          x11 0; y21 > 0; y31 > 0; y41 > 0; y12 > 0; y22 > 0; y32 > 0; y42 > 0; y13 > 0; y23 > 0; y33 > 0; y43 > 0, Аналогично записываются системы неравенств для второго и третьего столбцов матрицы C. Добавим в неравенства (2) неизвестные y11 > образующие матрицу Y , чтобы получить из (2) соответствующие системы линейных уравнений:                           1x11 −1x21 +0x31 +0x41 +1y11 = 0; 0x11 +1x21 −1x31 +0x41 +1y21 = 0; 0x11 +0x21 +1x31 −1x41 +1y31 = 0; 1x11 +0x21 +0x31 +0x41 −1y41 −1 = 0; 0x12 +1x22 −1x32 +0x42 +1y12 = 0; 0x12 +0x22 +1x32 −1x42 +1y22 = 0; −1x12 +0x22 +0x32 +1x42 +1y32 = 0; 0x12 +1x22 +0x32 +0x42 −1y42 −1 = 0; 0x13 +0x23 +1x33 −1x43 +1y13 = 0; −1x13 +0x23 +0x33 +1x43 +1y23 = 0; 1x13 −1x23 +0x33 +0x43 +1y33 = 0; 0x13 +0x23 +1x33 +0x43 −1y43 −1 = 0. Итак, информация, имеющаяся в матрице C, реализована при построении системы линейных (3) уравнений, которые можно использовать для определения основных (в матрице W)ивспомогательных (в матрице Y ) неизвестных. Следуя той же логике, что и при анализе матрицы C, перейдем к анализу матрицы R. Потребуем, чтобы элементы в четырех строках искомой матрицы W удовлетворяли тем же отношениям порядка, что и элементы в соответствующих строках матрицы R, а также чтобы наименьший по величине элемент в каждой строке матрицы W был не меньше единицы. Добавив в получившиеся неравенства неизвестные z11 > 0; z21 > 0; z31 > 0; z12 > 0; z22 > 0; z32 > 0; z13 > 0; z23 > 0; z33 > 0; z14 > 0; z24 > 0; z34 > 0, образующие матрицу Z, получим, что элементы в строках матрицы W должны удовлетворять системам линейных уравнений                   1x11 +0x12 −1x13 +1z11 = 0; 0x11 −1x12 +1x13 +1z21 = 0; 1x11 +0x12 +0x13 −1z31 −1 = 0; 0x21 +1x22 −1x23 +1z12 = 0; −1x21 +0x22 +1x23 +1z22 = 0; 0x21 +1x22 +0x23 −1z32 −1 = 0; −1x31 +0x32 +1x33 +1z13 = 0; 1x31 −1x32 +0x33 +1z23 = 0; 0x13 +0x23 +1x33 −1z33 −1 = 0; 0x41 −1x42 +1x43 +1z14 = 0; −1x41 +1x42 +0x43 +1z24 = 0; 0x41 +0x42 +1x43 −1z34 −1 = 0. Задача, таким образом, сводится к решению системы 2M · N линейных уравнений (3) и (4) с3M · N целыми положительными основными и вспомогательными неизвестными. Несмотря на то, что количество уравнений меньше количества неизвестных, система (а вместе с ней и задача получения матрицы W) может не иметь решений. Если же рассматриваемая система уравнений (3) и (4) совместна, то она имеет не единственное решение. Поэтому будем искать нормальное псевдорешение [1] систем (3) и (4), т. е. такой вектор решения w0 =(x11,x21,x31,x41,x12,. . . ,x33,x43)T , который среди всех решений системы имеет наименьшую норму (имеется в виду первая векторная норма  w1 = M·N РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 3 вып. 3 2016 k=1 |xk|). Таким образом, (4)
Стр.5
ИСПОЛЬЗОВАНИЕ МЕТОДОВ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ математическая формулировка задачи нахождения матрицы приоритетов W имеет вид: найти вектор w решения систем (3) и (4), минимизирующий функционал N вии, что векторы неизвестных w, y и z содержат только целые положительные компоненты. Это задача целочисленного линейного программирования. Решая эту задачу применительно к матриj=1M i=1 xij при услоцам C и R размерности 4 Ч 3, заданных выражениями (1), получаем матрицу приоритетов: W =         19 3 51 4 67 1 98 2         . денных в табл. 2 и табл. 3 соответственно, в результате использования описанного метода была получена матрица приоритетов: W =                   45 3 10 5910 4 5 26 78 3 6 11 15 67 29 8 23 1 6 7 109 711 89 14 34 2 12 8 13 .                   Для матриц C и R размерности 9× 7, приве51 по времени, получим поток заявок на входе в систему массового обслуживания (СМО), которой, собственно, и является совокупность приемных комплексов. Дисциплина обслуживания — «обслуживание с приоритетом», когда некоторые прохождения КАобслуживаются вне очереди. Функционирование СМО моделируется в дискретные моменты модельного времени T с интервалом в одну минуту в течение одних суток. Вкаждыймоментвремени T задается текущее «22-минутное скользящее окно» [2] и определяются все попавшие в него КАОГ, а также ПК, в ЗРВ которых эти КАимеют прохождения. Далее оцениваются состояния этих КАОГ (обслуживается/не обслуживается) и этих ПК (принимает/не принимает) в текущий момент времени T. Таким образом, вмомент T выявляется множество КА, требующих обслуживания (обозначим их количество через mM), и множество ПК, которые могут его предоставить (обозначим их количество через nN). С учетом этих данных формируется матри(5) ца приоритетностиW размерности m×n, строится функционал F = f(W) и применительно к моменту T решается распределительная задача, ответом к которой является назначение ПК для обслуживания КАиз «22-минутного скользящего окна». Функционал задачи, подлежащий минимизаЭлементы в столбцах матрицW целые и положительные (как и требуется от приоритетов) и удовлетворяют отношениям порядка, установленным оператором ОГ в соответствующих столбцах матрицы C и строках матрицы R. Получив матрицу приоритетности W,переходим к решению задачи оптимального планирования расписания сеансов приема ЦИ на ПК НКПОР. Постановка задачи планирования сеансов приема ЦИ на ПК НКПОР Для решения задачи ПСП применим метод математического моделирования, часто используемый при решении задач теории массового обслуживания. Определив прохождения КАОГ в ЗРВ ПК НКПОР за сутки и упорядочив эти прохождения ции в текущий момент T, формируется как F = m n  i=1  j=1 pijxij,(6) где pij — элементы матрицы приоритетовW; xij — элементы матрицы решения X. Элемент матрицы решения xij = 1, если ведется прием информации с i-го КАна j-м ПК, xij = 0—впротивном случае. Минимизация функционала обеспечивает максимальный «суммарный приоритет», получаемый на решении задачи. При поиске решения полагается, что каждый из m необслуженных КАдолжен передать информациюнаодинитольконаодиниз n свободных ПК; каждый свободный ПК должен быть задействован в приеме информации от одного и только одного, необслуженного КА. Таким образом, в каждом столбце и в каждой строке матрицы решения X можетнаходитьсятолькоодин элемент, РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 3 вып. 3 2016
Стр.6
52 А.А.ЧЕРНОВ равный 1, а все остальные их элементы равны нулю. Следовательно, xij могут рассматриваться как неотрицательные целочисленные решения системы m+n уравнений, являющихся ограничениями задачи оптимизации: xi1 +xi2 +.. .+xin = 1 (i = 1, . . . ,m); x1j +x2j +...+xmj = 1 (j = 1, . . . ,n). Сформулированная задача относится к классу транспортных задач целочисленного линейного программирования и известна как «задача о назначениях». Таким образом, при решении задачи планирования сеансов приема применительно к моменту T дискретного времени осуществляется распределение КАпо ПК. Если в отношении первого КА в «22-минутном скользящем окне» в рамках компьютерной программы принято решение о проведении обслуживания, то формируется сигнал «Разрешено обслуживание КА.. . на ПК . . .». Если такого решения нет — формируются сигналы «Перенос обслуживания КА... на ПК ...», «КА... уже обслуживается» или «Отказ обслуживания КА.. . на ПК . . .». Сигнал «Отказ обслуживания» при прохождении КА в ЗРВ ПК формируется в следующих случаях: 1) ПК при прохождении КАзанят; 2) ПК свободен, но в МП стоит запрет приема в «22-минутном скользящем окне» следует более приоритетный КА, который должен быть обслужен в первую очередь. Применительно к остальным КА из «22-минутного скользящего окна» решение, найденное в «задаче о назначениях», остается нереализованным, если момент их входа в ЗРВ ППИ еще не наступил. Описанные действия повторяются для следус i-го КАна j-м ПК (pij = 100); 3) ПК свободен, pij < 100, но за текущим КА ющего момента модельного времени (T + 1)   1440 мин и «задача о назначениях» решается вновь уже применительно к новому «22-минутному скользящему окну». В качестве примера в табл. 4 приведено решеи приемный комплекс, определенный в результате решения «задачи о назначениях» для его обслуживания. Так как приемных комплексов 7, а космических аппаратов 9, то 2 КА(«Метеор-М» и TERRA) ние «задачи о назначениях» для 9 КАи 7 ПК с матрицей приоритетовW, заданной выражением (5): Втабл.4элементы xij = 1указывают КА Таблица 4 (7) остались не обслуженными. Решение обеспечивает максимальный суммарный приоритет и соответствует тем предпочтениям и ограничениям, которые задал оператор при формировании матриц C и R. Следующий пример иллюстрирует получение расписания сеансов приема для рассматриваемой ОГ из 9 КАна 08.08.2015. На рис. 3 приведен фрагмент расписания сеансов приема, когда в качестве НКПОР используются всего два приемных комплекса на ППИ г. Железногорска (на рис. 3 использовано обозначение 6ПК и 7ПК ППИ «111»). ОГ имеет 54 прохождения за сутки в ЗРВ ППИ «111». Матрица приоритетов задана выражением (5). В результате решения задачи получено расписание, содержащее 35 сеансов приема (64,8% от 54). Отказы в обслуживании КА (19 отказов) вызваны запретами в матрице приоритетов (5), а также использованием свободных ПК для приема ЦИ с более приоритетных КА. Например, строки 18, 19 и 20 расписания демонстрируют случай, когда программа «отказывает» в обслуживании КАTERRA и «Метеор-М» на свободном 7ПК, т. к. через 9 мин назначит 7ПК для приема ЦИ с более приоритетного КА«Метеор-М2». Заключение Рассмотрена задача оптимального распределения сеансов приема целевой информации, получаемой от орбитальной группировки КАДЗЗ, по приемным комплексам НКПОР в условиях их недостаточности. Распределение сеансов приема ставилось РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 3 вып. 3 2016
Стр.7
ИСПОЛЬЗОВАНИЕ МЕТОДОВ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ 53 Рис. 3. Фрагмент отчета «Планирование сбросов ЦИ» в зависимость от приоритетности информации, принимаемой на интервале планирования, на множестве всех КАи ПК. Приоритетность информации задавалась «матрицей приоритетов», определяющей относительную важность сеансов приема целевой информации с КАОГ на ПК. В работе предложен автоматизированный алгоритм нахождения «матрицы приоритетов» (т. е. «весовой» матрицы), который позволяет учесть некоторые априорные требования Оператора ОГ, предъявляемые к конечному результату планирования. «Матрица приоритетов» затем используется для построения критерия оптимальности и решения собственно задачи оптимального планирования расписания сеансов приема. Задачи построения «матрицы приоритетов» и оптимального планирования расписания сеансов приема ставятся как транспортные задачи целочисленного линейного программирования. Ограничения и особенности задачи планирования расписания сеансов приема позволяют поставить ее как известную в целочисленном линейном программировании «задачу о назначениях». Приведены примеры расчета «весовых матриц» и оптимального расписания сеансов приема, полученнные с помощью разработанного программного комплекса. Данный алгоритм может быть задействован для построения «весовых» матриц в задачах оптимизации, в которых уместно использование понятия приоритетности при формировании критерия оптимальности. Список литературы 1. Воеводин В. В. Вычислительные основы линейной алгебры. М.: Наука, 1977. 304 с. 2. Гончаров А.К., Чернов А.А. Планирование сеансов приема информации с космических аппаратов орбитальной группировки при ограниченном количестве приемных комплексов // Космонавтика и ракетостроение, 2014, вып. 1 (74). С. 180–189. РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 3 вып. 3 2016
Стр.8