Рекомендация № МСЭ-R S.1256-0 Международного союза электросвязи
Актуально на 04.11.2026.
Как цитировать: Рекомендация N МСЭ-R S.1256-0 Международного союза электросвязи. МСЭ-R S.1256-0. https://bukva.ai/docs/mse-r-s-1256-0
«Методология определения максимальной суммарной плотности потока мощности на геостационарной спутниковой орбите в полосе 6700 — 7075 мГц от фидерных линий негеостационарных спутниковых систем подвижной спутниковой службы в направлении космос-Земля (Вопрос МСЭ-R 206/4)»
(Принята 28.05.1997 Ассамблеей радиосвязи МСЭ)
МЕЖДУНАРОДНЫЙ СОЮЗ ЭЛЕКТРОСВЯЗИ
РЕКОМЕНДАЦИЯ № МСЭ-R S.1256-0 МЕТОДОЛОГИЯ ОПРЕДЕЛЕНИЯ МАКСИМАЛЬНОЙ СУММАРНОЙ ПЛОТНОСТИ ПОТОКА МОЩНОСТИ НА ГЕОСТАЦИОНАРНОЙ СПУТНИКОВОЙ ОРБИТЕ В ПОЛОСЕ 6700 — 7075 МГЦ ОТ ФИДЕРНЫХ ЛИНИЙ НЕГЕОСТАЦИОНАРНЫХ СПУТНИКОВЫХ СИСТЕМ ПОДВИЖНОЙ СПУТНИКОВОЙ СЛУЖБЫ В НАПРАВЛЕНИИ КОСМОС-ЗЕМЛЯ (ВОПРОС МСЭ-R 206/4)
a) что полоса 6700 — 7075 МГц распределена на первичной основе фиксированной спутниковой службе (ФСС) для передач в направлении космос-Земля для использования фидерными линиями негеостационарных спутниковых сетей подвижной спутниковой службы (ПСС);
b) что полоса 6700 — 7075 МГц распределена на первичной основе также ФСС для передач в направлении Земля-космос, а полоса 6725 — 7025 МГц используется в соответствии с Планом распределения по Приложению 30B Регламента радиосвязи (РР) для геостационарных спутниковых сетей;
c) что в соответствии с положением п. S22.5A Регламента радиосвязи максимальная суммарная плотность потока мощности (п.п.м.), создаваемая в пределах +/- 5° от геостационарной спутниковой орбиты (ГСО) негеостационарной спутниковой системой ФСС, не должна превышать — 168 дБ(Вт/м2) в любой полосе шириной 4 кГц;
d) что Резолюция 115 Всемирной радиоконференции (Женева, 1995 г.) (ВКР-95) предлагает МСЭ-R разработать методологию для определения максимальной суммарной плотности потока мощности на ГСО от негеостационарной спутниковой сети;
e) что параметры орбиты и характеристики передачи негеостационарных спутниковых сетей подвижной спутниковой службы определены в § A.3 vii) Приложения 1 к Резолюции 46 (Пересм. ВКР-95),
1 чтобы для определения максимальной суммарной плотности потока мощности (дБ(Вт/м2)) в любой полосе шириной 4 кГц, в любой точке в пределах углов +/- 5° от ГСО, создаваемой фидерными линиями негеостационарной спутниковой сети, работающей в полосе 6700 — 7075 МГц в направлении космос-Земля, использовалась методология, приведенная в Приложении 1.
Для вычисления суммарной п.п.м., создаваемой спутниковой сетью со спутником, расположенным на негеостационарной орбите, в одиночной испытательной точке на ГСО, требуется компьютерное моделирование всей НГСО спутниковой группировки и испытательной точки на ГСО.
Учитывая, что в обычных условиях ГСО спутник делает полный оборот на геостационарной орбите за период времени около TГСО = 24 ч. и что орбитальный период НГСО спутника (TНГСО) вовсе не обязательно является полным делителем TГСО, для оценки наихудшего возможного сценария, при котором п.п.м. в точках ГСО имела бы максимальный уровень, необходимо выполнять статистическое моделирование, требующее значительных затрат времени.
Для оценки максимальной п.п.м. в любой точке ГСО можно выполнить простое моделирование, требующее намного меньших затрат времени. Вместо реального вращающегося на орбите ГСО спутника рассматривается фиксированная испытательная точка на ГСО, орбитальная позиция которой фиксирована на Картезианской системе координат 0xyz (см. рисунок 1), но не во вращающейся системе координат Земли. Учитывая это, поскольку орбитальный период НГСО спутников равен TНГСО, подразумевается, что позиция НГСО спутников, как она видна из фиксированной точки на ГСО (см. рисунок 1), будет повторяться, по крайней мере, один раз за время орбитального периода TНГСО. Более того, в том случае, когда НГСО спутники в каждой орбитальной плоскости распределены равномерно, одно и то же геометрическое расположение НГСО спутников будет повторяться с периодом, равным TНГСО/Ns (где Ns представляет собой количество НГСО спутников, равномерно распределенных по плоскости). На основе этих положений суммарный уровень п.п.м. (суммированный для всех видимых НГСО спутников) в испытательной точке на ГСО будет повторяться с этим же периодом.
Рисунок 1 Геометрия ГСО/НГСО группировки для вычисления п.п.м.: \Delta\Omega = 0^{\circ}
Для каждого отрезка времени можно рассчитать п.п.м., а максимальная суммарная п.п.м. для выбранной испытательной точки на ГСО можно получить за период моделирования от T0 до T0 + TНГСО / Ns.
Значение, полученное для конкретной точки на ГСО на рисунке 1, не обязательно будет максимальным уровнем п.п.м. Для того чтобы определить максимально возможный уровень суммарной п.п.м., эта процедура может быть выполнена для других точек ГСО при увеличении угла \Delta\Omega (см. рисунок 2) между испытательной точкой на ГСО и линией, соединяющей узлы НГСО. Этот второй итерационный процесс будет выполняться для углов \Delta\Omega от 0° до \Delta\Omega_{max} = 360^{\circ}N_{p} , где Np — количество орбитальных плоскостей НГСО спутника. Если Np — четное (как для ОЗО-F и ОЗО-D), то \Delta\Omega_{max} = 180^{\circ}N_{p} .
Этот метод можно также использовать для любой НГСО группировки, которая не удовлетворяет вышепоказанным орбитальным требованиям (например, распределение спутников на орбите — неравномерное, эллиптические орбиты). В таких случаях моделирование будет выполнено для времени, равного минимальному периоду повторения конфигурации группировки, которое во многих случаях равно периоду вращения группировки TНГСО.
В § 2 приведены все основные уравнения, необходимые для получения суммарного уровня п.п.м., создаваемой данной НГСО сетью в данной испытательной точке на ГСО, а на рисунке 3 показан алгоритм программной реализации описанной здесь методологии.
Рисунок 2 Геометрия ГСО/НГСО группировки для вычисления п.п.м.: \Delta\Omega \neq 0^{\circ}
Рисунок 3 Последовательность операций методологии
Шаг 1: Орбитальная позиция НГСО спутников
Рисунок 4 Негеостационарная орбита и эталонные системы
На рисунке 4 показаны различные параметры, значения которых необходимо знать для полноценной оценки любого из моментов позиции любого НГСО спутника на его орбите. Эти параметры перечислены в § A.3 vii) Приложения 1 к Резолюции 46 (Пересм. ВКР-95):
a: большая полуось, для круговой орбиты большая полуось является постоянной и равна радиусу орбиты
I: наклон орбиты относительно экваториальной плоскости
\Omega_{j} : прямое восхождение восходящего угла для j-той орбитальной плоскости, измеренное против часовой стрелки в экваториальной плоскости от направления точки весеннего равноденствия на точку, где спутник пересекает экваториальную плоскость, переходя из южного полушария в северное (0° \leq \Omega_{j}<360^{\circ})
\omega_{p} : аргумент перигея, для круговой орбиты перигей равен апогею и, следовательно, \omega_{p} можно приравнять 0°
\omega_{i} : первоначальный фазовый угол i-го спутника на его орбитальной плоскости в эталонный момент времени t = 0, измеренный от точки восходящего узла (0° \leq \omega_{i}<360^{\circ})
\theta : истинная аномалия спутника.
Для группировки НГСО спутников, использующих круговые орбиты, a и I будут постоянными, а \omega_{p} будет равен нулю, при этом изменение позиции каждого спутника будет определяться параметрами \Omega и \theta .
Для круговой орбиты угловая скорость спутника постоянна, следовательно, угловая позиция спутника равна его истинной аномалии и вычисляется как:
\theta(t)_{i,j} = \frac{360^{\circ}}{T} t + \omega_{i,j} (1)
для i = 1 — Ns и j = 1 — Np, где Ns — количество спутников в каждой орбитальной плоскости, Np — количество орбитальных плоскостей, а T — орбитальный период в секундах, вычисляемый как:
T = 2\pi \sqrt{a^{3}\mu} , (2)
где \mu — геоцентрическая гравитационная постоянная, равная 3,986 E14(м3с-2).
Различные значения \Omega_{j} будут зависеть от геометрии группировки и будут внесены в множество элементов, содержащихся в § A.3 vii) Приложения 1 к Резолюции 46 (Пересм. ВКР-95). Тот же принцип применяется к значениям \omega_{i,j} .
Зная для каждого спутника истинную аномалию \theta_{i,j}(t) и прямое восхождение восходящего узла \Omega_{j} , можно вычислить его геоцентрические координаты:
x(t)_{i,j} = a[cos\Omega_{j}cos\theta(t)_{i,j} - cosIsin\Omega_{j}sin\theta(t)_{i,j}] (3)
y(t)_{i,j} = a[sin\Omega_{j}cos\theta(t)_{i,j} + cosIcos\Omega_{j}sin\theta(t)_{i,j}] (4)
z(t)_{i,j} = a[sinIsin\theta(t)_{i,j}] . (5)
Положение испытательной точки ГСО относительно линии узлов НГСО группировки определяется параметром \Delta\Omega (см. § 1). Отсюда в уравнениях (3), (4) и (5) \Omega_{j} = \Omega_{j,0} + \Delta\Omega , где \Delta\Omega лежит в пределах от 0 до \Delta\Omega_{max} (см. § 1) и \Omega_{j,0} = \Omega_{j} для \Delta\Omega = 0 .
Шаг 2: Расстояние между НГСО спутником и испытательной точкой на ГСО
xГСО, yГСО и zГСО представляют собой геоцентрические координаты испытательной точки ГСО, вычисленные как:
aГСО: большая полуось геостационарной орбиты (42 164 км)
IГСО: наклон геостационарной орбиты (- 5° <= IГСО <= 5°).
Эти уравнения остаются постоянными во время моделирования, хотя проще изменять \Omega_{j} в уравнениях (3), (4) и (5) путем увеличения сдвига \Delta\Omega .
Теперь расстояние между НГСО спутником и позицией испытательной точки на ГСО может быть вычислено с использованием теоремы Пифагора:
d(t)_{i,j} = \sqrt{(x_{GSO} - x(t)_{i,j})^{2} + y(t)_{i,j}^{2} + (z_{GSO} - z(t)_{i,j})^{2}} . (9)
Шаг 3: Вычисление внеосевого угла излучения НГСО антенны в направлении на испытательную точку ГСО
На рисунке 5 в виде двумерной диаграммы показана геометрия внеосевого излучения НГСО спутника относительно позиции испытательной точки на ГСО.
Рисунок 5 Вычисление \phi_{i,j}
Внеосевой угол излучения НГСО антенны можно определить с использованием теоремы Карно (известной под названием теоремы «косинуса»):
\phi(t)_{i,j} = arccos(\frac{a^{2} + d(t)_{i,j}^{2} - a_{\text{ГСО}}^{2}}{2ad(t)_{i,j}}) . (10)
Шаг 4: Вычисление коэффициента усиления НГСО антенны в боковых лепестках в направлении на испытательную точку ГСО
Зная внеосевой угол излучения, вычисленный в уравнении (10), можно вычислить коэффициенты внеосевого усиления антенны, G(\phi(t)_{i,j}) , для каждого видимого спутника. Однако, как видно из рисунка 5, это необходимо, только если \phi(t)_{i,j} превышает минимальное значение \phi_{min} , полученное из:
\phi_{min} = arcsin(Ra) . (11)
Шаг 5: Вычисление уровня суммарной п.п.м. в направлении на испытательную точку ГСО
Суммарный уровень п.п.м. можно выразить как:
\text{п}.\text{п}.\text{м}.(t) = \frac{P_{\text{пик}},4\text{кГц}}{4\pi} i \cdot j \sum_{= 1 - N} (t)_{v} \frac{G(\phi(t)_{i,j})}{d(t)_{i,j}^{2}} \text{для }\phi(t)_{i,j} \geq \phi_{min} , (12)
Pпик, 4 кГц: пиковая мощность в наихудшей полосе шириной 4 кГц на входе антенны НГСО спутника, величина которой предполагается постоянной и одинаковой для всех НГСО спутников
N(t)v: количество ГСО спутников, видимых из испытательной точки на ГСО в момент времени t.
3 Общее количество шагов моделирования и приращений шага моделирования
Для вычисления максимальной суммарной п.п.м. в направлении на ГСО от НГСО сети, времени шага \Delta t и шага прямого восхождения \delta\Omega требуются два шага моделирования.
Хотя нет непосредственных сонаправленных помех от НГСО спутников (используют ли они антенны с малым коэффициентом усиления или помехи обусловлены боковыми лепестками передающей антенны), различные варианты моделирования (для LEO-D и LEO-F) показывают, что для получения достоверных результатов достаточно использовать приращение угла с шагом не более 0,5°. Следовательно, вычисление должно выполняться как:
\Delta t = \frac{T(s) \times 0,5^{\circ}}{360^{\circ}}
Полное время моделирования для каждой испытательной точки ГСО и общее количество испытательных точек ГСО показано в § 1.