Одним из путей повышения точностных характеристик системы ГЛОНАСС считается создание глобальной сети станций сбора измерений (ССИ). Математическая постановка задачи развертывания этой сети и возможные пути решения задачи рассмотрены в [1]. Следует обратить внимание на ряд моментов. Как правило, ССИ должны располагаться в населенных пунктах с необходимой инфраструктурой или вблизи от них. Исключения могут составлять важные с географической точки зрения необитаемые районы Земного шара, например, в Антарктиде. Это обстоятельство предопределяет комбинаторный характер рассматриваемой оптимизационной задачи с большим количеством альтернативных решений. Например, если анализировать возможность размещения 20 станций в 100 населенных пунктах, то множество альтернатив будет включать в себя 1018 вариантов, полный перебор которых на современных персональных компьютерах потребует 170 миллионов лет.
В этих условиях необходимо на первых этапах проектирования оценить ориентировочный количественный состав сети ССИ и ее предпочтительную топологию (диапазон расстояний между ближайшими станциями, направления и т.д.). Такие оценки могут быть получены на основе решения задачи комбинаторной геометрии о покрытии сферы кругами минимального радиуса. Эта задача известна давно, и хотя ею занимались крупные математики (Л.Ф. Тот [4], Б.Л. Ван дер Варден, К. Шютте), она решена лишь для немногих частных случаев. То же относится и к родственной задаче о редчайшем расположении точек на сфере, которой интересовались еще Кеплер и Ньютон [3].
Точная нижняя оценка размера зоны покрытия станции при заданном их количестве в системе может быть получена из решения задачи покрытия сферы кругами одинакового радиуса [2,3]. Ниже будет рассмотрен метод численного поиска локально-оптимального решения, а также некоторые вспомогательные задачи, которые используются при анализе многогранника, образованного точками расположения станций, для расчета показателей качества функционирования глобальной сети станций.
Имеется K кругов одинакового радиуса. Требуется найти на сфере такое расположение этих кругов, чтобы сфера была ими покрыта, а их радиус при этом был бы минимальным.
Пусть координаты центра i-го круга определяются двумя угловыми координатами . Тогда целевая функция – необходимый размер зоны обзора при данном расположении центров кругов, определяемый парами , определится выражением:
,
где – угловое расстояние от центра i-го круга до точки на сфере с координатами .
Требуется определить координаты центров кругов, обеспечивающих минимальное значение целевой функции
В основе решения задачи лежат следующие три достаточно очевидных утверждения.
Если предварительно соединить все ближайшие центры кругов дугами большого круга, то в результате получается некоторый сферический многогранник, состоящий из сферических граней и сферических ребер. Если у обычного многогранника ребрами являются отрезки прямых, а гранями – треугольники или многоугольники, то у рассматриваемого сферического многогранника – части дуг и сферические треугольники или многоугольники соответственно. Понятно, что между размерами этих элементов существует простое геометрическое соотношение. Кроме того, очевидно, что объединение всех сферических граней образует сферу. Отсюда вытекают следующие утверждения.
Утверждение 1. Для того, чтобы сфера была покрыта К кругами, необходимо и достаточно, чтобы ими была покрыта каждая из граней многогранника.
Утверждение 2. Если радиус каждого из К кругов не уступает радиусу окружности, описанной вокруг грани многогранника, то эта грань будет покрыта данными кругами.
Здесь следует иметь в виду два случая:
а) центр описанной вокруг грани окружности принадлежит этой грани. Тогда сформулированное условие является необходимым и достаточным;
б) центр описанной вокруг грани окружности не принадлежит этой грани. Тогда для покрытия грани радиус окружности должен быть не менее половины большего из ребер, образующих данную грань, что несколько меньше, чем радиус описанной окружности. Однако это обстоятельство не имеет смысла учитывать, поскольку в этом случае покрытие центра описанной окружности должно быть учтено при рассмотрении смежной грани. А радиус соответствующей описанной окружности будет заведомо больше, чем исходной. Поэтому в дальнейшем, без ущерба для конечного результата, будем считать, что справедливо следующие утверждение.
Утверждение 3. Требуемый минимальный радиус каждого из К кругов определяется максимальным из радиусов окружностей, описанных вокруг каждой из граней:
,
где J – множество граней многогранника из K вершин, j – номер грани.
Определение радиуса описанной вокруг грани окружности проводится следующим образом. Если грань является многоугольником (не треугольником), то ее можно представить в виде совокупности треугольников. Описанные вокруг них окружности совпадают и являются окружностью, описанной вокруг всей грани. Поэтому, не теряя общности, будем находить радиус окружности, описанной вокруг одного из треугольников, который либо является гранью, либо ее частью.
Из геометрических построений вытекает, что искомый радиус описанной вокруг j-й грани определяется выражением
; ; ,
где R=1 – радиус сферы, hj – расстояние от центра сферы до j-й грани.
Отсюда следует, что максимум целевой функции достигается при минимуме функции hj. Поэтому в дальнейшем будет исследована функция hj, определяющая расстояние от центра сферы до каждой из граней.
Пусть уравнение плоскости, содержащей j-ю грань, в прямоугольной системе координат имеет вид
.
Здесь и далее индекс j опускается и считается, что номера вершин, участвующих в образовании j-й грани, равны 1, 2 и 3. Расстояние от центра сферы до грани определяется соотношением
. (1)
Из аналитической геометрии известно, что уравнение плоскости, проходящей через три точки, может быть представлено в виде
, (2)
где координаты вершины определяются по формулам:
; ; ,
где λ – угол, отсчитываемый против часовой стрелки в плоскости ОХУ от оси ОХ до плоскости, перпендикулярной ОХУ и содержащей одну из вершин, ; φ – угол, отсчитываемый в указанной плоскости против часовой стрелки от ОХУ до самой вершины.
Из выражения (2) можно получить
Отсюда коэффициенты уравнения плоскости определяются выражениями:
; ; ;
Последние соотношения можно преобразовать к виду:
; ; ; . (3)
Метод варьирования положения центров кругов на сфере базируется на традиционной схеме методов синтеза: выбор начального приближения, затем итерационное изменение (варьирование) параметров оптимизации, последовательно улучшающее значение целевой функции до достижения требуемой точности. Одно из центральных мест в этом процессе принадлежит расчету производных и последующему их анализу. Поэтому требуется найти частные производные от целевой функции по начальному положению одного из центров кругов (для определенности – первого круга), воспользовавшись (1):
;
(4)
.
Производные, входящие в правые части этих выражений, определятся соотношениями:
; ;
; ;
; ;
; . (5)
Найденные частные производные по соотношениям (4) с учетом (5) позволяют определять приращения расстояния от центра сферы до каждой из граней следующим образом:
, (6)
где j – номер грани; i – номер вершины, причем i≠1, поскольку первая вершина принимается как базовая, ее положение совмещается с началом системы координат, т.е. , , и в дальнейшем эта вершина не смещается.
Частные производные вычисляются по соотношениям (4) и (5), если i -я вершина принадлежит j -й грани; в противном случае они принимаются равными нулю. Для численного поиска решения используется метод Ньютона.
В таблице 1 представлены параметры, характеризующие топологию многогранников, широты параллелей, на которых располагаются вершины каждого слоя, а также основные характеристики соответствующей сети станций: ρпокр – минимальный радиус кругов, который задает потребный размер зоны радиовидимости, а также δ – угол места, определяющий условия наблюдения космического аппарата (КА) над местным горизонтом. Структура многогранников, представленных в таблице 1, имеет определенную топологию:
– все К вершин (первый столбец таблицы), за исключением одной или двух, равномерно распределены вдоль m параллелей (слоев) по n в каждом слое;
– в каждом из полюсов может располагаться одна из вершин (второй столбец), таким образом количество вершин определяется соотношением К=m*n+Kполюс, где Kполюс – количество вершин в полюсах (величина KПОЛЮС может принимать значение 0, 1 и 2);
– все точки одного слоя вдоль параллели размещены равномерно, на угловом расстоянии друг от друга;
– в очередном слое вершины смещены относительно предыдущего на угол .
Анализ представленных результатов позволяет сделать ряд важных выводов:
1. В идеализированных условиях для однократного покрытия орбит КА системы ГЛОНАСС достаточно четырех станций, при этом угол наблюдения космических аппаратов будет больше величины δ=5 °.
Таблица 1
Параметры многогранника глобальной сети станций сбора измерений
Количество вершин |
Распределение вершин |
ρпокр, град |
Широты параллелей (слоев) |
δ, град |
||||||
Кол-во вершин в полюсах |
Кол-во слоев |
Кол-во в слое |
1 |
2 |
3 |
4 |
5 |
|||
4 |
1 |
1 |
3 |
70,53 |
-19,47 |
|
|
|
|
5,08 |
8 |
- |
2 |
4 |
49,93 |
40,06 |
-40,06 |
|
|
|
27,24 |
1 |
1 |
7 |
61,70 |
-28,29 |
|
|
|
|
14,3 |
|
2 |
2 |
3 |
53,13 |
26,57 |
-26,57 |
|
|
|
23,64 |
|
9 |
- |
3 |
3 |
45,88 |
44,12 |
0 |
-44,12 |
|
|
31,86 |
1 |
2 |
4 |
46,38 |
16,85 |
-43,62 |
|
|
|
31,3 |
|
10 |
- |
3 |
3 |
45,88 |
44,12 |
0 |
-44,12 |
|
|
31,86 |
1 |
2 |
4 |
46,38 |
16,85 |
-43,62 |
|
|
|
31,3 |
|
- |
2 |
5 |
48,03 |
41,97 |
-41,97 |
|
|
|
29,4 |
|
2 |
2 |
4 |
42,31 |
24,47 |
-24,47 |
|
|
|
36,04 |
|
12 |
- |
2 |
6 |
47,06 |
42,94 |
-42,94 |
|
|
|
30,5 |
- |
3 |
4 |
45 |
45 |
0 |
-45 |
|
|
32,88 |
|
- |
4 |
3 |
37,38 |
52,63 |
10,81 |
-10,81 |
-52,63 |
|
41,92 |
|
2 |
2 |
5 |
37,38 |
26,57 |
-26,57 |
|
|
|
41,92 |
|
14 |
- |
2 |
7 |
46,49 |
43,51 |
-43,51 |
|
|
|
31,18 |
2 |
2 |
6 |
34,94 |
27,65 |
-27,65 |
|
|
|
44,86 |
|
2 |
3 |
4 |
36,21 |
35,26 |
0 |
-35,26 |
|
|
43,3 |
|
2 |
4 |
3 |
37,22 |
49,91 |
13,58 |
-13,58 |
-49,91 |
|
42,1 |
|
15 |
- |
3 |
5 |
44,16 |
45,84 |
0 |
-45,84 |
|
|
33,88 |
- |
5 |
3 |
34,04 |
55,96 |
22,83 |
0 |
-22,83 |
-55,96 |
45,96 |
|
1 |
2 |
7 |
39,03 |
17,71 |
-50,97 |
|
|
|
39,92 |
|
16 |
- |
2 |
8 |
46,13 |
43,87 |
-43,87 |
|
|
|
31,6 |
- |
4 |
4 |
33,78 |
56,23 |
5,59 |
-5,59 |
-56,23 |
|
46,28 |
|
2 |
2 |
7 |
33,55 |
28,29 |
-28,29 |
|
|
|
46,62 |
|
17 |
1 |
2 |
8 |
38,29 |
17,79 |
-51,71 |
|
|
|
40,82 |
2 |
3 |
5 |
33,59 |
33,51 |
0 |
-33,51 |
|
|
46,52 |
|
2 |
5 |
3 |
34,21 |
56,90 |
23,29 |
0 |
-23,29 |
-56,90 |
45,76 |
|
18 |
- |
2 |
9 |
45,89 |
44,11 |
-44,11 |
|
|
|
31,86 |
3 |
3 |
|
35,26 |
54,74 |
0 |
-54,74 |
|
|
44,48 |
|
2 |
2 |
8 |
32,67 |
28,70 |
-28,70 |
|
|
|
47,64 |
|
2 |
4 |
4 |
31,43 |
44,87 |
13,59 |
-13,59 |
-44,87 |
|
49,18 |
|
20 |
- |
2 |
10 |
45,72 |
44,28 |
-44,28 |
|
|
|
32,06 |
- |
4 |
5 |
32,32 |
57,68 |
3,46 |
-3,46 |
-57,68 |
|
48,08 |
|
- |
5 |
4 |
30,76 |
59,25 |
27,55 |
0 |
-27,55 |
-59,25 |
50,0 |
|
2 |
2 |
9 |
32,09 |
28,98 |
-28,98 |
|
|
|
48,36 |
|
2 |
3 |
6 |
30,46 |
36,01 |
0 |
-36,01 |
|
|
50,38 |
2. Учет повышенных требований по точности измерений, возможного расположения станций в холмистой местности или вблизи высотных зданий приводит к необходимости расширения сети станций до 8-9 (при этом δ=25 °-30 °).
3. Невозможность расположения станций на территории океанов обуславливает необходимость расширения сети до 15-20 станций.
Анализ геометрических характеристик многогранников, представленных в таблице 1, с учетом их реализуемости с точки зрения географического фактора позволяет обосновать следующие рекомендации по размещению глобальной сети станций:
1. Наиболее предпочтительными географическими широтами расположения станций являются приэкваториальные (±6 ° и 20 °) и средние (56 °-59 °).
2. Расстояние между соседними станциями могут изменятся в пределах от пяти до семи тысяч километров.
Таким образом, предложенная методика позволяет обосновывать рекомендации по выбору облика глобальной сети ССИ системы ГЛОНАСС, включающие приблизительный количественный состав станций, предпочтительные географические широты размещения станций и расстояния между местами расположения станций. Эти рекомендации позволяют значительно сократить множество альтернатив на заключительном этапе обоснования облика глобальной сети станций сбора измерений системы ГЛОНАСС.
Рецензенты:
Петров Г.Д., д.т.н., профессор, начальник кафедры организации эксплуатации и технического обеспечения вооружения, военной и специальной техники ФГКВОУ ВПО «Военно-космическая академия имени А.Ф. Можайского» Министерства обороны Российской Федерации, г. Санкт-Петербург.
Басыров А.Г., д.т.н., профессор, начальник кафедры информационно-вычислительных систем и сетей, ФГКВОУ ВПО «Военно-космическая академия имени А.Ф. Можайского» Министерства обороны Российской Федерации, г. Санкт-Петербург.
Библиографическая ссылка
Власов В.А., Горбулин В.И., Каргу Д.Л., Паршин А.В. МЕТОДИКА ВЫБОРА ВАРИАНТОВ ОБЛИКА ГЛОБАЛЬНОЙ СЕТИ СТАНЦИЙ СБОРА ИЗМЕРЕНИЙ СИСТЕМЫ ГЛОНАСС // Современные проблемы науки и образования. – 2015. – № 1-1. ;URL: https://science-education.ru/ru/article/view?id=18877 (дата обращения: 08.12.2024).