2008
Доклады БГУИР
№ 3 (33)
УДК 621.396:621391.82
ВЕРОЯТНОСТНЫЕ ХАРАКТЕРИСТИКИ ЭЛЕКТРОМАГНИТНОЙ ОБСТАНОВКИ, ОБРАЗУЕМОЙ ЭЛЕКТРОМАГНИТНЫМИ ИЗЛУЧЕНИЯМИ НАЗЕМНЫХ РАДИОЭЛЕКТРОННЫХ СРЕДСТВ, НАД ПОВЕРХНОСТЬЮ ЗЕМЛИ
В.И. МОРДАЧЁВ
Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 21 мая 2008
Приведены результаты моделирования распределений вероятности сигналов радиоэлектронных средств размещенных на поверхности Земли и определяющих электромагнитную обстановку над ее поверхностью, по энергетическому параметру и по углу визирования направления прихода в точке наблюдения, размещенной над поверхностью на определенной высоте.
Ключевые слова: электромагнитная совместимость, электромагнитная обстановка, угол визирования, плотность потока мощности, распределение вероятности.
Введение
Внешние условия эксплуатации радиоэлектронных средств (РЭС), размещаемых на летательных аппаратах различных типов (самолеты, вертолеты, ракеты и т.п.) и на искусственных спутниках Земли (ИСЗ) обладают рядом особенностей, которые следует учитывать при обеспечении электромагнитной совместимости (ЭМС) этих РЭС по отношению к РЭС, размещаемым на поверхности Земли. С увеличением высоты точки размещения РЭС над земной поверхностью уменьшается вероятность присутствия вблизи РЭС мощных источников электромагнитных излучений (ЭМИ), что влечет снижение динамического диапазона сигналов в рассматриваемой точке. Вместе с тем с увеличением высоты подъема РЭС значительно увеличивается зона прямой видимости и соответственно число посторонних радиосигналов на входе радиоприемного устройства (РПУ) РЭС, по уровню превышающих порог его восприимчивости.
В основу подхода к статистическому моделированию электромагнитной обстановки (ЭМО) может быть положена методика [1], предполагающая случайное равномерное размещение источников ЭМИ равной эквивалентной изотропно излучаемой мощности (ЭИИМ) на плоской поверхности, и использующая модель распространения радиоволн (РРВ) в виде гиперболической зависимости плотности потока мощности электромагнитного поля П от расстояния R до его источника:
C P
П = —Cv = const, (1)
r v 5 v > V '
где Petr — эквивалентная изотропно излучаемая мощность (ЭИИМ); Cv — константа; v — параметр, определяющий "скорость" затухания электромагнитного поля по мере увеличения рас-
стояния до источника излучения. В частности, v=2 при РРВ в свободном пространстве; v=4 при РРВ с интерференцией прямого и отраженного лучей в дальней зоне [2-7].
Для плоской модели земной поверхности в [1, 2, 4] определены следующие модели: Модель пространственной области потенциального мешающего взаимодействия РЭС в виде круговой области радиуса
Rmax =(CvPJПт1п ) (2)
вокруг точки размещения на поверхности РЭС-рецептора помех, порог восприимчивости которого равен Птт.
Плотность распределения вероятности удаленности РЭС-источника ЭМИ от РЭС-рецептора помехи:
w(R) = RR; R < Rmax. (3)
max
Плотность распределения вероятности w(T) плотности потока мощности сигналов в точке размещения РЭС-рецептора помех на плоской поверхности:
2n2/v n2/v C P
w (П) = (П2!v ПЧЛТТ(2+v)/v ; П G [ Пmax ] Птах = , (4)
Т max Т min JП mm
где Rmin — минимальное расстояние до ближайшего РЭС-источника ЭМИ — дополнительно введенное в [1] искусственное ограничение на минимально возможное расстояние до ближайшего источника ЭМИ, обеспечивающего существование начальных и центральных моментов распределения (4); при отсутствии этого ограничения
^nrnin
w (ПЬтж/г; П * Пmm . (5)
Плотность распределения вероятности угла визирования источников ЭМИ, размещенных на плоской поверхности, из точки над поверхностью [4]:
* ч 2 cos 8 Г .„1 w(e) = —2--—, 8Е[8min, п/2]. (6)
Ctg 8min S1- 8
Непосредственное использование этих моделей для описания статистических характеристик ЭМО на значительной высоте над поверхностью Земли, когда кривизной земной поверхности и высотой точки размещения РЭС-рецептора помех пренебречь нельзя, без должного обоснования представляется некорректным. Данная работа посвящена обоснованию вида моделей (3)-(6), учитывающих кривизну земной поверхности в зоне прямой видимости и высоту точки наблюдения над поверхностью.
Плотность распределения вероятности сигналов над сферической поверхностью Земли
по энергетическому параметру
Рассмотрим случай, когда точка наблюдения ТН, соответствующая точке размещения РЭС, находится на определенной высоте Н над поверхностью, а также учтем сферичность земной поверхности (рис. 1). Будем полагать, что на поверхности в зоне прямой видимости случайно с плотностью р(а,Л) (единиц/км2) размещены РЭС-источники ЭМИ; здесь и на рис. 1 а, е — горизонтальный (азимут) и вертикальный углы визирования РЭС-источников ЭМИ из точки наблюдения.
Рис. 1. Схема расположения точки наблюдения на высоте H над земной поверхностью: LH — линия горизонта; ен — угол визирования горизонта; RЗ — радиус Земли (с учетом рефракции); RH — удаленность линии горизонта от точки наблюдения; ЯН>Ятах
Площадь 5Шс шарового сегмента высоты h, ограничиваемого удаленностью Я от точки, поднятой на высоту Н над сферической поверхностью радиуса ЯЗ, равна:
S^ =Пз (R2-H2)) R +H).
Отсюда среднее число источников ЭМИ в ограниченном элементарными приращениями da,dR элементе шарового слоя, наблюдаемого на расстоянии R и по азимуту а из точки размещения РЭС, будет определяться соотношением:
R
dN (a, R) = р(а, R) dSшс (R, da, dR) =-Rp(a, R) da dR.
R3 + H
В результате получаем общее выражение для совместной плотности распределения w(a,R) вероятностей азимута и удаленности источников ЭМИ по отношению к точке наблюдения и для среднего числа NScp источников ЭМИ в области потенциального мешающего взаимодействия РЭС, ограниченной удаленностью Rmax:
. 1 dN (a, R) R3 _ , m
w(a, R) = —--/р „ч Rp(a, R), (7)
nscp da dR nscp (R + H)
R Rmax
NScp = R-+H J J R P(a, R) da dR;
h 0
Я е [ Н, Ятах], Ятах < Ян 2 Яз Н + Н2 ае [0,2п]. (8)
При равномерном пространственном распределении источников в пределах поверхности сферического сегмента, ограниченной удаленностью Ятах, имеем р(а,Я)=р=соп81:, и распределение (7) приобретает следующий вид:
^(а, Я) = ^(а)^(К)
Я 1 2Я
- н2 )2п ятах -н
(9)
2 Я
Н2 "
(10)
тах
Таким образом, при р(а,Я)=р=соп81 вид ^(а,Я) инвариантен к кривизне земной поверхности и зависит исключительно от вида функции р(а,Я) средней территориальной плотности источников (от их координат). Для случая равномерного распределения источников по территории плотности распределения (3) и (10) отличаются только областями определения величины Я и, соответственно, нормирующими множителями.
Отсюда следует принципиально важный вывод о том, что и вид плотности распределения вероятности w(П) плотности потока мощности сигналов в точке размещения РЭС-рецептора помех над плоской поверхностью при свободном РРВ инвариантен к кривизне земной поверхности и к высоте подъема точки наблюдения над поверхностью и соответствует модели (4) с областью определения, правая граница которой определяется высотой подъема над земной поверхностью:
В этом нетрудно убедиться, используя модель (10) для определения по известным правилам закона распределения вероятности w(П) с использованием функционального преобразования (1).
Этот вывод существенно расширяет области применимости базовой модели (4) статистической теории ЭМС [1], где эта модель использована для описания статистических характеристик ЭМО в точке на плоской поверхности от размещенных на этой поверхности РЭС-источников ЭМИ. Введенное в [1] искусственное ограничение на минимально допустимое расстояние до ближайшего РЭС-источника ЭМИ, существенно ослабляющее адекватность модели для случаев совпадения (пересечения) областей возможного пространственного размещения точки наблюдения и РЭС-источников ЭМИ, в рассматриваемом случае приобретает естественный характер: минимальное расстояние до ближайшего РЭС-источника ЭМИ может быть принято равным высоте подвеса точки наблюдения над поверхностью с размещенными на ней источниками.
В заключение следует отметить также ряд важных особенностей данных моделей.
1. Модель (7) совместной плотности распределения w(а,Я) вероятностей азимута и удаленности источников ЭМИ по отношению к точке наблюдения позволяет определить статистические характеристики ансамбля сигналов, образующих ЭМО в точке наблюдения, при произвольном виде зависимости средней пространственной плотности р(а,Я) РЭС-источников сигналов от координат. Это имеет существенный практический интерес, поскольку с увеличением высоты точки наблюдения над поверхностью и соответствующим расширением границ области прямой видимости и/или области потенциального мешающего взаимодействия РЭС возможности использования допущения о близком к равномерному случайном пространственном размещении РЭС-источников ЭМИ (р(а,Я)=р=соп81;) в ряде случаев становятся менее очевидными. С другой стороны, использование ГИС технологий позволяет определить функцию р(а,Я) для значительных территорий с учетом границ населенных пунктов и промышленных зон, транспортных коммуникаций и других объектов и территориальных образований, насыщенных радиооборудованием.
2. Очевидно, что инвариантность модели (7) к кривизне земной поверхности и к высоте подъема точки наблюдения над поверхностью означает ее инвариантность к высоте точки наблюдения над поверхностью и в случае, когда кривизна земной поверхности в расчет не принимается, что согласуется с предложениями [4] по стандартизации моделей ЭМО.
3. Все приведенные модели и соображения в силу присущего рассматриваемой задаче дуализма справедливы и для "обратной" ситуации, когда в точке наблюдения над поверхностью
Птах = О,Рег/НV , У=2.
земли размещается РЭС-источник ЭМИ, а в области потенциального мешающего взаимодействия радиуса Ятах<ЯН случайно с плотностью р(а,Я) либо случайно равномерно располагаются РЭС-рецепторы помех, имеющие порог восприимчивости Пт,„. Применительно к этой ситуации модели (7)-(10) описывают плотность распределения w(a,Я) вероятностей азимута и удаленности РЭС-рецепторов ЭМИ по отношению к точке наблюдения и для среднего числа Ы5ср рецепторов ЭМИ в области потенциального мешающего взаимодействия РЭС, а модель (4) описывает плотность распределения вероятности плотности потока мощности сигнала в точке размещения на поверхности произвольно выбранного РЭС-рецептора помех от поднятого над поверхностью РЭС-источника ЭМИ.
Плотность распределения вероятности сигналов над сферической поверхностью Земли
по углу визирования направления прихода
При размещении точки наблюдения над плоской поверхностью с размещенными на ней РЭС-источниками ЭМИ величины R и в (угол визирования источника радиосигнала) связаны; для этого случая R=H / cos в. Отсюда для случая p(a,R)=p=const нетрудно получить модель (6) плотности распределения вероятности угла визирования источника радиосигнала [4].
В случае учета сферичности земной поверхности, вид w(b) существенно усложняется. Если ограничиться рассмотрением углов визирования, соответствующих направлениям прихода радиосигнала между перпендикуляром и касательной к сферической земной поверхности, распределение w(b) имеет следующий вид:
w(b) = 4(в) + А2(в) - 4(в); (11)
ве[вн, п/2], Вн = arccos [Rj(R3 + H)]; (( + H) cos в JR2 -(R3 + H)2cos2 в
4(в) = ■
А2(в) =
Я3Н
(Яз + Н )%т в Б1П28 2ЯзН^Я^ -(Яз + Н )2 соб2 в '
4(в) =(яз + Н^ ;
3 ЯзН
область определения этой модели ограничена углом вН визирования линии горизонта.
Во многих случаях вместо достаточно громоздкой модели (11) может использоваться существенно более простая модель (6). Эта замена особенно эффективна для высот полета самолетов, т.е. для точки размещения РЭС-рецептора помех не более 10-20 км над земной поверхностью, а также при выборе следующей области определения модели (6):
Вфтт, П/2], 8тШ =8н/2- (12)
Для такой области определения модели (6) она обеспечивает очень хорошее совпадение с моделью (11) в области углов визирования [1,5 еН , п/2]. В области углов визирования вблизи горизонта [еН, 1,5еН] следует пользоваться точной моделью (11); часть [еН, еН/2] области (12) введена из технических соображений и не имеет практического смысла.
На рис. 2-5 приведены расчетные кривые распределений (6), (11) для различных значений высоты Н РЭС-рецептора помех над земной поверхностью при Я3=8500 км (с учетом рефракции) и для различных значений левой границы области определения аргумента в (6); на этих
рисунках сплошная линия соответствует модели (6), обозначенной как ^>(в); прерывистая линия соответствует модели (11), обозначенной как ^^в). На правых частях рисунков приведены изображения ^>(в), ^5(в) в области малых углов визирования, в которой различия между моделями (6) и (11) могут быть существенными. На рис. 6 приведены те же расчетные кривые для высоты Н=10 км и для области определения [ен, л/2] модели (6). Сравнение рис. 6 и рис. 3 свидетельствует о практической невозможности замены модели (11) моделью (6) без корректировки (12).
100 10 1
ю-1
10"2
ю-3
ю-4
ю-5
1 1 1 1 -
- \
1 1 1 Гч
20 40 60 80
Рис. 2. Расчетные кривые распределений (6), (11): высота Н РЭС-рецептора помех над земной поверхностью 1 км при ЯЗ=8500 км; еН~0,9°, ет,„=еН/2
100 10 1
ю-1
10"2
ю-3
ю-4
ю-5
1 1 1 1 -
| 1 1 I 1
20 40 60
Рис. 3. Расчетные кривые распределений (6), (11): высота Н РЭС-рецептора помех над земной поверхностью 10 км при ЯЗ=8500 км; еН~2,8°, ет1И=еН/2
100 10 1
ю-1 10"2 ю-3 ю-4 ю-5
1 1 1 1
-Чч
| 1 1 1
20 40 60 80
Рис. 4. Расчетные кривые распределений (6), (11): высота Н РЭС-рецептора помех над земной поверхностью 100 км при ЯЗ=8500 км; еН~8,8°, ет,„=еН/2
100 10 1
- 1 1 1 III-
-
1 1 1 1 1 1
10
'. 1 1 1 1
-
-
1 1 1 1
30 40 50 60 70 80 8 23 30 32 34
Рис. 5. Расчетные кривые распределений (6), (11): высота Н РЭС-рецептора помех над земной поверхностью 1000 км при ЯЗ=8500 км; еН~26,5°, етп=еН/2
Рис. 6. Расчетные кривые распределений (6), (11): высота Н РЭС-рецептора помех над земной поверхностью 10 км при ЯЗ=8500 км; еН~2,8°, етПп=еН
Сравнительный анализ моделей (6) и (11) на основе данных рис.2-6 позволяет сделать следующие выводы:
Применение модели (6) при етп=еН/2 обеспечивает очень хорошее совпадение с моделью (11) для широкого диапазона высот точки наблюдения при углах визирования более 1,5 еН.
Применение модели (6) при высотах точки наблюдения, соответствующих диапазону высот установившегося полета современных пассажирских и транспортных самолетов, и при етп=еН обеспечивает существенно завышенную плотность вероятности при углах визирования менее 3,5°-4,0°, т.е. при превышении углом визирования менее чем на 1,5° угла визирования горизонта.
Применение модели (6) при етп=еН обеспечивает существенно завышенные значения плотности вероятности сигналов РЭС-источников ЭМИ в точке наблюдения для углов визирования е>(1,2.. ,1,5)еН по сравнению с моделью (11).
Как и в предыдущем случае, приведенные соображения в силу дуализма рассматриваемой задачи справедливы и для "обратной" ситуации, когда в точке наблюдения над поверхностью земли размещается РЭС-источник ЭМИ, а в области потенциального мешающего взаимодействия радиуса Ятсх<Ян случайно равномерно располагаются РЭС-рецепторы помех, имеющие порог восприимчивости Птп. Применительно к этой ситуации модели (6), (11) соответственно приближенно и точно описывают плотность распределения вероятностей угла визирования произвольно выбранного РЭС-рецептора помех по отношению к точке наблюдения, в которой расположен поднятый над поверхностью РЭС-источник ЭМИ.
Заключение
На основании изложенного можно сделать вывод, что в важном с практической точки зрения случае, характеризуемом близким к равномерному случайным распределением источников ЭМИ, равной ЭИИМ по земной поверхности, вид плотности распределения вероятностей этих сигналов по энергетическому параметру в точках над поверхностью соответствует модели (4) и инвариантен к кривизне земной поверхности и к высоте точки наблюдения над поверхностью. Доказано, что при условии (12) в качестве плотности распределения вероятности сигналов по углу визирования направления прихода в точку наблюдения над поверхностью может использоваться упрощенная модель (6), полученная без учета кривизны земной поверхности. Полученные результаты существенно расширяют возможности применения основных положений статистической теории ЭМС [1], обеспечивая возможность системного анализа электромагнитной совместимости РЭС, функционирующих над поверхностью Земли, по отношению к РЭС на ее поверхности
PROBABILISTIC CHARACTERISTICS OF ELECTROMAGNETIC ENVIRONMENT ABOVE GROUND LEVEL PRODUCED BY EARTH-BASED RADIO STATIONS
V.I. MORDACHEV Abstract
Results of modeling of power flux density and viewing angle probability distributions of radio signals generated by earth based radio stations and forming an electromagnetic environment above ground level on spatial points elevated on some altitude above the spherical ground surface are presented.
Литература
1. Апорович А.Ф. Статистическая теория электромагнитной совместимости радиоэлектронных средств. Минск, 1984.
2. Мордачев В.И. Статистические характеристики электромагнитной обстановки в рассредоточенных группировках РЭС. Научно-технический отчет инв. № 02880034534. Минск, 1987.
3. Мордачев В.И. // Тр. IX Междунар. Вроцлавского симпозиума по электромагнитной совместимости, июнь 1988 г. Вроцлав, 1988. С. 571-576.
4. Мордачев В.И. // Труды X Международного Вроцлавского симпозиума по электромагнитной совместимости, июнь 1990 г. Вроцлав, 1990. С. 409-414.
5. Mordachev V. // 15th Intern. Wroclaw Symp. on EMC. Poland, Wroclaw. June 27-30. 2000. P. 331-335.
6. Mordachev V. // The IEEE Semi Annual VTC Fall 2000. Boston, Sept. 24-28, 2000.
7. Rappaport T.S. Wireless Communications: Principles and Practice. Prentice Hall, Upper Saddle River, 2002.