УЗБЕКСКОЕ АГЕНТСТВО СВЯЗИ И
ИНФОРМАТИЗАЦИИ
ТАШКЕНТСКИЙ УНИВЕРСИТЕТ ИНФОРМАЦИОННЫХ
ТЕХНОЛОГИЙ
Факультет радиотехники, радиосвязи и
телерадиовещания
Кафедра
антенно-фидерных устройств
Ликонцев Д.Н.
СПЕЦИАЛЬНЫЕ ВОПРОСЫ КУРСА
«РАСПРОСТРАНЕНИЕ РАДИОВОЛН»
для направления образования
5А522105 - «Мобильные системы связи»
ПРЕДИСЛОВИЕ
Учебное пособие «Специальные
вопросы курса Распространение радиоволн» содержит материал, отсутствующий в
учебниках, по основным темам дисциплин Антенны мобильных и спутниковых систем
связи», «Распространение радиоволн» и предназначен для использования в процессе
самостоятельной работы магистрантов и на аудиторных групповых занятиях. В
данном пособии имеются графические и другие материалы, собранные из различных
литературных источников и приведенные к удобному виду для их использования в
учебном процессе.
1.
ОСОБЕННОСТИ РАСПРОСТРАНЕНИЯ ЗЕМНЫХ РАДИОВОЛН В ГОРОДСКИХ УСЛОВИЯХ
1.1. Влияние отдельных зданий на
распределение уровней поля
В УКВ диапазоне крупные и железобетонные строения практически
радионепрозрачны. Казалось бы, в области тени непосредственно за зданием
уровень поля должен быть еще ниже, однако это не всегда так. Причиной тому -
волны, отраженные от других зданий, которые подсвечивают теневые области и
создают в них напряженность поля, значительно превышающую напряженность
прошедшей через здание волны.
Измерения напряженности поля за одиночно стоящими строениями,
выполненные на частоте 627 МГц, показали, что за кирпичным зданием
напряженность поля на 26 дБ ниже, чем перед ним, а за железобетонным строением
уровень сигнала падает на 35 дБ.
Поверхности зданий размеры которых во много раз превышают длину
волны, создают интенсивные отражения,
придающие процессу распространения радиоволн в городе многолучевой характер.
При одновременном разрешении волн по углу и времени прихода сигнала в точку
наблюдении возможна идентификация отражения от крупных строений [1…3]. Сложная конфигурация отражающих поверхностей
(оконные проемы, балконы, лоджии и т.д.) порождает рассеяние радиоволн в
широком секторе углов при крайне изрезанной по форме индикатрисе рассеяния. Вместе
с тем в направлениях,
соответствующих зеркальному отражению, интенсивность излучения значительно возрастает.
Измерения, выполненные на частоте 9,4 ГГц с высоким временным (50 нс) и угловым (10°) разрешением,
показали что уровень сигнала, отраженного в зеркальном направлении, на 30 дБ
превышает уровень сигнала, рассеянного в других направлениях [4].
Приведем результаты двух экспериментов, в которых напряженность
поля измерялась вблизи одиночно стоящих зданий. В первом эксперименте
исследовалось пространственное распределение поля на частоте 88 МГц вокруг
девятиэтажного крупнопанельного здания, высота которого равнялась
Рис.1.1.
Геометрия экспериментального исследования
Напряженность поля
измерялась через каждые
Рис.1.2. Результаты эксперимента на частоте 88 МГц
Измерения амплитуды
отраженной волны проводились, когда передатчик находился в точке O, а приемник перемещался от точки В до точки D (рис.1.1). Результаты эксперимента приведены на
рис.1.3, где амплитуда Е отраженной волны дана по отношению к
амплитуде прямой волны , измеренной в точке D. В среднем уровень
отраженного сигнала оказался на 10 дБ ниже уровня сигнала, приходящего
непосредственно от источника. Резкие изменения амплитуды отраженной волны при
перемещении приемника, очевидно, обусловлены сложной конфигурацией отражающей
поверхности здания.
Рис.1.3. Зависимость соотношения амплитуд
отраженной и падающей волн от расстояния
Для оценки ослабления поля в тени здания измерения проводились
тогда, когда передатчик перемещался по прямой от точки т до точки п (рис.1.1)
параллельно зданию на расстоянии
Второй
эксперимент проводился одновременно на двух частотах: 88 и 627 МГц. [6] Измерялась напряженность поля вокруг одиночного пятиэтажного
панельного здания при вертикальной поляризации излучения. На рис.1.4 схематично
показаны здание и замкнутый путь, по которому перемещался приемник, а также
указаны области I, II, III,
в которых пространственное распределение поля имело различный характер. Границы
этих областей, в основном, определялись геометрией эксперимента, В области II
интерференция прямой и отраженной волн обусловила периодические изменения
напряженности поля, пространственный масштаб которых прямо пропорционален длине
волны.
Рис.1.4.
Геометрия экспериментального исследования
Область III -
область тени, в которой уровень принимаемого сигнала мал, а пространственное
распределение поля неоднородно. В области I интерференция прямой
волны и рассеянного зданием излучения порождает периодические изменения
напряженности поля, глубина которых значительно меньше, чем в области III,
а пространственный масштаб больше.
На рис.1.5 показано изменение уровня принимаемого сигнала для трех
областей при различных положениях передатчика. На рис.1.5. графики построены для частоты 88 МГц при трех отмеченных
точками на оси абсцисс положениях передатчика. На рис.1.5.б аналогичные графики
приведены для частоты 627 МГц.
Поскольку большинство
крупных городских строений практически непрозрачны, что приводит к образованию
в городе теневых зон и в значительной степени определяет свойства
формирующегося поля. Для выявления общих закономерностей в пространственном
распределении поля целесообразно отвлечься от локальных особенностей,
определяемых конкретной планировкой районов, расположением отдельных зданий, и
использовать статистические методы описания городской среды и возникающих в ней
затенений.
Рис.1.5.
Результаты эксперимента на частотах 88 и 627 МГц
В качестве модели городской застройки примем множество крупных
непрозрачных объектов, случайно расположенных на плоской поверхности земли [7].
Рельеф города в этой модели представляется резко пересеченной поверхностью с
крутыми неровностями, по высотам которых можно выделить приземный слой - слой
городской застройки. При описании затенений важной характеристикой слоя
является вероятность прямой видимости между двумя разнесенными в нем точками. К
ее расчету подойдем следующим образом. Проведем луч от одной точки к другой и,
перемещаясь вдоль него, будем фиксировать пересечения с рельефом застройки, но
лишь те, при которых луч входит в неровности. Предположим, что каждое новое
пересечение возникает независимо от уже зафиксированных. В этом случае
последовательность пересечений является простым и ординарным потоком, в котором
число событий подчиняется закону распределения Пуассона, то есть вероятность
того, что отрезок луча длиной r испытывает п пересечений с рельефом, равна
(1.1)
где gо - плотность потока,
равная среднему числу пересечений, фиксируемых на отрезке единичной длины.
Чтобы выяснить, при каких условиях приемлемо сделанное
предположение в [7] провели статистическую обработку планов нескольких
городских районов. В среднем плотность застройки в выбранных районах составляла
90 зданий на 1 км2, число этажей зданий изменялось от 5 до14.
Результаты показали, что при r
≥
Рис.1.6. Статистическое распределения n-числа пересечений
горизонтального отрезка длиной
заметны отклонения:
частота пересечений оказывается меньше, чем следует из формулы Пуассона. На
рис.1.6 показано статистическое распределение п - числа пересечений
горизонтального отрезка длиной r =
Из (1.1) следует формула для вероятности прямой видимости между
двумя точками, разнесенными на расстояние r у плоской поверхности земли. Эта вероятность
совпадает с Рr(0) и равна [6]
. (1.2)
На рис.1.7 показана
зависимость от r величины . Пунктирная прямая соответствует формуле (1.2), х - значения, полученные при обработке
планов городских районов (общее число испытаний 400).
Рис.1.7.
Зависимость от r величины lnPr(0)
Обозначим через r расстояние вдоль луча от
его начала до первого пересечения с рельефом застройки. Если для потока
пересечений принять распределения (1.1) и (1.2), то плотность вероятности для r является
экспоненциальной
(1.3)
со средним значением , которое
определяет среднюю горизонтальную дальность прямой видимости в слое городской
застройки. Очевидно, gо и должны определяться характеристиками
застройки: плотностью, размерами зданий, их взаимным расположением и т.д.
Расчет этих параметров выполнен для двух упрощенных моделей застройки. В первой
множество зданий заменялось множеством плоских отражающих экранов,
ортогональных к земной поверхности и размещенных на ней статистически
независимо [8]. Для gо получено значение
, (1.4)
где n - плотность застройки,
определяемая как среднее число зданий, размещенных на 1 км2 площади
городского района, - средняя длина экранов
(зданий). Во второй модели рассматривалось множество непрозрачных прямоугольных
параллелепипедов, случайно расположенных на плоскости [9]. В среднем периметр
основания одного параллелепипеда равнялся . Для γо
и в этом случае получено выражение, совпадающее с (1.3). Например, средняя
горизонтальная дальность прямой видимости , рассчитанная по (1.4) для v = 90 км-2 и =
Обратимся вновь к формуле (1.2). Она определяет вероятность прямой
видимости между точками у поверхности земли и не учитывает высоты зданий. По мере
подъема точек исчезают затенения, создаваемые низкими зданиями, и дальность
прямой видимости возрастает. Это обстоятельство учтено в работе [10].
Предположим, что в рассматриваемом районе здания, имеющие высоты H1, H2,...,HN, встречаются с частотой р1, р2,..., рN соответственно. Определим вероятность PH(h) того, что высота здания
превышает h следующей формулой
, (1.5)
где - ступенчатая функция
Хевисайда, равная единице при и нулю при Hi<h. Две точки наблюдения А1
и А2 подняты над поверхностью земли на высоту h1 и h2, расстояние между их
проекциями на горизонтальную плоскость равно . В этом
случае вероятность прямой видимости между А1 и А2
с учетом высот зданий при h1>h2 определяется формулой
, (1.6)
где
учитывает улучшение
видимости при подъеме точек над поверхностью земли. Из (1.6) следует, что
средняя дальность прямой видимости из точки А1 на
точку A2 равна
Предположим, что точка А2
принадлежит отражающей поверхности, которую будем рассматривать как плоский
вертикально установленный на поверхности земли экран конечных размеров. Между
точкой А1 и точкой А2 есть прямая
видимость. Очевидно, что при этом также будет прямая видимость на все точки
отражающей поверхности, расположенные над точкой A2. Если из точки А1
будет виден горизонтальный отрезок длиной l, то будет просматриваться вертикальная полоса шириной l над этим отрезком. В [6, 8] показано, что случайно расположенные
на плоской поверхности экраны не затеняют горизонтального отрезка l<< от точки A1 с вероятностью
,
(1.7)
равной произведению
вероятности (1.6) незатенения точки А1 на условную
вероятность незатенения отрезка l, для которого точка А2
является внутренней точкой. Угол y определяет положение
отражающего экрана относительно луча g12, а e12 также, как g12, учитывает улучшение
видимости при подъеме точек и равно
(1.8)
Из (1.8) следует, что на расстоянии от точки А1
источник освещает участки отражающих экранов с горизонтальными размерами l, распределенными экспоненциально с плотностью вероятности
и средним значением
. (1.9)
Приведенные результаты
открывают возможность сравнительно простого расчета пространственного
распределения тех точек рассеяния волн в слое городской застройки, которые
«освещаются» источником и видны из точки наблюдения. Пусть положение
источника А в системе
координат определяется
вектором , а точки наблюдения . Рассчитаем вероятность того, что точка будет точкой
рассеяния, отвечающей указанным условиям - «светящейся» точкой для наблюдателя
в точке В. Предположим, что через С проходит рассеивающий экран высотой Hs, ориентированный под углом j к , причем j Î(0, a) U (p, p+a), где a - угол между векторами и (рис.1.8),
поскольку только в этом случае возможно рассеяние волны, приходящей из точки А
в точку В. Вероятность пересечения бесконечно малой окрестности
точки С отражающим экраном равна vLdn, где dn - элемент длины в ортогональном к экрану направлении. Учитывая (1.6)
и пренебрегая возможной статистической зависимостью пересечений на отрезках АС
и ВС, примем, что точка С не будет затеняться относительно А и
В одновременно с вероятностью P1S ×P2S.
Рис.1.8.
Геометрия расчета
Если лучи АС и ВС не пересекаются затеняющими
экранами, то из точек А и В одновременно на отражающей
поверхности будет видна вертикальная полоска, содержащая точку С. Ее
ширина на уровне точки С будет не
меньше l с вероятностью [6]
следовательно,
элементарный участок отражающей поверхности площадью dldh будет содержать незатененную точку с вероятностью
.
Умножим этот результат на
вероятность vLdn того, что точка С окажется на отражающем экране, и
обозначим - элемент объема слоя застройки. Усредняя результат по js
и Hs и опуская элемент объема,
получим следующее выражение для плотности распределения точек рассеяния
, (1.10)
где HS описывает высотный профиль слоя
Если здания можно сгруппировать в N групп, характеризуемых высотами Hi и частотами pi , i
= 1, 2, … N, то
Заметим, что название
«распределение точек рассеяния» в данном случае неточно отражает смысл (1.10).
При вычислениях учтены средние размеры освещенных участков на рассеивающих
экранах. Поскольку каждый участок - источник рассеянного излучения, то
плотность (1.10) определяет среднее число таких источников в единице объема
слоя городской застройки. В ряде случаев удобно рассматривать распределение
точек рассеяния при их проекции на горизонтальную плоскость, то есть после
интегрирования плотности (1.10) по толщине слоя застройки. Этот результат для
простоты рассмотрим в частном случае, когда высоты всех зданий одинаковы и
равны H. Будем считать, что в точке А находится
подвижный пункт, антенна которого поднята над землей на небольшую высоту h2<<H, а в точке В над
городской застройкой установлена антенна базовой станции h1>H. Примем проекцию точки В
на горизонтальную плоскость за начало полярной системы координат (r, j), в которой отсчет
азимутального угла j ведется от направления
на точку А. В малоугловом приближении, определяемом неравенствами
плотность точек
рассеяния, спроецированных на плоскость, представляется формулой [48]
(1.11)
где - расстояние от точки А до точки рассеяния. Первое слагаемое в
(1.11) описывает относительно редкие переизлучатели, распределенные по большой
площади городского района. Второе слагаемое характеризует
«светимость» небольшой окрестности подвижного пункта,
размер
которой определяется расстоянием прямой видимости r и зависит от плотности застройки. Число
переизлучателей в ней пропорционально отношению (h1 - H)/r, а плотность переизлучателей достигает максимального значения на
уровне верхних этажей зданий и крыш, вероятность затенения которых относительно
антенны базовой станции мала.
Заменим
переменную в (1.11) на
относительное время t распространения волны от
источника до базовой станции а также нормируем плотность распределения на среднее число
однократно рассеянных волн
|
В результате получим
плотность вероятности w(t, j) совместного
распределения относительных задержек t и углов прихода j однократно рассеянных
волн в точку приема [6]
На рис.1.9 изображен
рельеф плотности вероятности w(t,j) в относительных единицах
(коэффициент 1/4 опущен) при расстоянии r=2 км, плотности
Рис.1.9.
Рельеф плотности вероятности w (t, j)
застройки v = 100 км-2, =
По мере уменьшения высоты подвеса h1 антенны базовой станции максимум плотности
расширяется по углу, но вероятность больших значений t уменьшается. При h1 = H, когда антенна базовой
станции установлена на уровне крыш окружающих зданий, распределение сохраняет
унимодальный характер, его максимум остается хорошо выраженным, но значительно
расширяется, занимая по единичному уровню угловой интервал ( - 80°,
+80°). При этом интервал относительных задержек сокращается до
t Î (1; 1,6).
1.2. Влияние застройки
города на уровень поля
Определением величин ослабления УКВ в условиях
города занимались многие исследователи в различных странах
мира и, к настоящему времени накоплен обширный экспериментальный материал
[11…19,…]. Во многих работах подчеркивается крайняя нерегулярность
пространственного распределения электромагнитного поля. Обширные затенения
[11…14, 20], создаваемые строениями, практически исключают возможность прямого
прохождения сигнала между передающей антенной и пунктом приема, а рассеянные и
отраженные волны создают сложную интерференционную структуру поля, в которой
минимумы и максимумы напряженности чередуются с периодом от половины длины
волны до нескольких длин волн. В результате сигналы, передаваемые между
центральной станцией и движущимся в городе объектом, подвержены глубоким и
резким замираниям, повторяющимся с частотой до 100 Гц и достигающим 40 дБ
относительно среднего уровня.
Нерегулярность пространственного распределения напряженности поля,
частые глубокие замирания и относительно редкие, но большой интенсивности
выбросы амплитуды сигнала привели к тому, что внимание исследователей было
обращено на медианное значение поля как более устойчивое в этих условиях по
сравнению со средним значением [12, 21].
Результаты экспериментальных исследований показывают, что
ослабление поля с расстоянием в условиях города нарастает значительно быстрее,
чем в свободном пространстве [14], и является примерно одинаковым для
однотипных городов. По данным [22, 23], ослабление поля на частоте 836 МГц в
среднем на 15…25 дБ больше, чем в свободном пространстве; в [23] для частоты
11,2 ГГц приводятся данные об ослаблении порядка 25…56 дБ для Нью-Йорка и 26…40
дБ для загородной местности в штате Нью-Джерси.
В ряде работ [6, 22…25] отмечается значительное влияние планировки
улиц городского района на величину ослабления радиоволн. На продольных (вдоль
трассы) улицах [14, 26] сигнал лучше распространяется вдоль широких улиц, чем
вдоль узких; ослабление на поперечных улицах на 3…8 дБ больше, чем на
продольных [14, 24, 27].
Наиболее полные экспериментальные результаты, отличающиеся высокой
степенью систематизации и позволяющие установить ряд существенно важных
закономерностей для медианного значения напряженности поля УКВ в условиях
города, получены Окамурой [12], который выполнил обширные экспериментальные
исследования в районе Токио на частотах 453, 922, 1430 и 1920 МГц при высотах
подвеса передающих антенн от 45 до
С увеличением частоты ослабление сигнала в условиях города
возрастает [12, 27, …30]. Если мощность сигнала считать пропорциональной
некоторой степени длины волны, то на расстояниях до 10…15 км эта степень на
частоте 100 МГц оказывается близкой к нулю, а на частотах 1…2 ГГц - к единице
[12, 14]. Этот вывод согласуется с результатами других исследователей.
Например, в [13] приводятся данные о том, что увеличение потерь в городе по
сравнению с потерями над плоской землей на частотах 86, 167 и 441 МГц
составляет 16, 19 и 33 дБ соответственно; на частоте 11,2 ГГц ослабление на 15
дБ больше, чем на частоте 450 МГц [23]. Установленная закономерность
подтверждается [13] в ряде других работ.
На больших расстояниях (20…100 км) частотная зависимость
ослабления сигнала возрастает и средняя интенсивность поля оказывается
пропорциональной квадрату длины волны. В условиях города на ослабление сигнала
различным образом влияют высоты подвеса передающих и приемных антенн. Для
расстояний менее
Установленные эмпирические закономерности отражены в серии
графиков, рекомендованных МККР [31]. Они предназначены для оценки медианного
значения уровня поля УКВ в квазиплоских городских районах при связи подвижного
объекта с базовой станцией. Метод оценки по этим графикам получил название
«метода Окамуры». К его недостаткам можно отнести то, что расчет по серии графиков
включает в себя вычисление довольно большого числа промежуточных параметров и
требует учета корректирующих характеристик. При этом остается неясным,
насколько результаты, справедливые для Токио и его окрестностей, применимы для
городов с другими типами застройки и планировки.
Для расчета среднего ослабления поля УКВ в нижней части диапазона
(ниже 100 МГц) в одной из наиболее ранних работ, предлагалось влияние городских
строений на распространение волн от высоко поднятого передатчика к мобильному
приемному пункту описывать как поглощение в некотором абсорбирующем слое.
Однако эти представления не получили распространения, поскольку результаты
экспериментальных исследований указывают на существенно многолучевой характер
поля в условиях города, на определяющую роль отражений и рассеяния волн на
крупных строениях, а не на поглощения в них энергии радиоволн [2, 4, 14, 15].
Не оправдали себя и попытки расчета среднего ослабления по скорректированной
формуле Введенского, предпринимаемые до настоящего времени [6]. Действительно,
квадратичная формула Введенского справедлива в условиях, когда поле в точке
приема формируется при интерференции прямой и отраженной от подстилающей
поверхности волн. Однако в городских условиях неровности для УКВ велики,
затенения и многолучевость, порождаемые отражениями от строений, становятся
основными факторами, определяющими пространственное распределение поля. Это
обстоятельство и пытаются учесть с помощью корректирующих множителей. При этом
в большинстве случаев полностью изменяется характер зависимости среднего
ослабления от расстояния, частоты и высот подвеса антенн.
Анализ экспериментальных и расчетных данных показывает, что
наиболее полное объяснение установленной экспериментально зависимости
интенсивности принимаемого сигнала от расстояния, высот подвеса антенн, частоты
излучения и других параметров дает расчет, учитывающий статистику затенений и
отражений УКВ в городской среде [6, 32…36]. Он основывается на статистическом
модельном описании городской среды как множества крупных случайно размещенных
на поверхности земли объектов, вертикальные поверхности которых отражают
радиоволны и создают обширные затенения. Считается, что стены зданий
электрически неоднородны и неровны, их свойства статистически независимы.
Случайный комплексный коэффициент отражения падающего на стену излучения имеет
масштабы корреляции по отражающей поверхности lВ и lг в вертикальном и
горизонтальном направлениях, квадрат его модуля в среднем равен Г, а аргумент с
равной вероятностью принимает значения от 0 до 2p. В приближении Кирхгофа при статистическом усреднении по ансамблю
реализаций застройки городского района получено среднее дифференциальное
сечение рассеяния <s1S2> элемента городской
застройки, находящегося в точке , при условии, что источник излучения находится в точке , а приемник - в точке
(1.12)
Здесь - волновое
число, l - длина волны;
; ;
- расстояние от
источника до точки рассеяния, a - от точки рассеяния
до приемника. Формула (1.12) позволяет записать в приближении однократного
рассеяния среднюю интенсивность поля УКВ, создаваемого в городе сосредоточенным
ненаправленным излучателем, находящимся в точке
(1.13)
где интегрирование
ведется по объему V
слоя городской застройки.
Как известно, для расчета поля земной волны в зонах дифракционной
тени и полутени на сферической поверхности используется дифракционный множитель
ослабления F [37]. Введем его в подынтегральное выражение
(1.13), учитывая ослабление падающей на слой волны, и вычислим асимптотическое
значение интеграла при γ0r>>1. В результате получим [6, 35, 38]
(1.14)
при этом высота подвеса
передающей антенны в дифракционном множителе должна отсчитываться от уровня h = H, второй же точкой должна
быть не точка а наиболее вероятная
точка однократного рассеяния с координатами (х1, у1, H). Расчет средней интенсивности (1.14) проведен для слоя городской
застройки высотой H. Реально город является
статистически неоднородным по высоте застройки, поэтому необходимо рассмотреть,
как повлияет на среднюю интенсивность наличие в городе зданий различной высоты.
В работе [38] показано, что введение РH(h) в интегральное соотношение (1.13) и анализ полученных результатов
для широкого класса моделей высотной застройки городов подтверждают
правомерность аппроксимации неоднородного по высоте слоя городской застройки
однородным слоем со средней высотой H.
На рис.1.10. пунктирными линиями 1 - 4 представлены зависимости
медианного значения мощности сигнала от расстояния для частоты 453 МГц при h2 =
Рис.1.10. Медианные значения мощности на входе
прёмного
устройства на частоте 453 МГц
Верхняя пунктирная прямая
показывает уровень сигнала в свободном пространстве. Расчетные зависимости
показаны сплошными линиями 1 - 4. Расхождение расчетных и измеренных
значений становится заметным при h1 =
Рис.1.11. Зависимость ослабления в городе по
отношению к свободному пространству от частоты и расстояния
Сплошные линии представляют расчетные зависимости (1.14),
пунктирные - экспериментальные [12, 14].
Проведенное сопоставление показывает, что предположение о
превалирующей роли отражений в формировании поля УКВ и статистическое описание
затенений в городской среде приводят к расчетным результатам, объясняющим
зависимость r–3 средней интенсивности на
расстояниях, меньших дальности видимого горизонта , где а -
радиус Земли. Учет явления дифракции позволяет продолжить расчетные кривые в
зоны полутени и тени, обеспечивая плавное нарастание ослабления, хорошо
согласующееся с экспериментальными данными. Из (1.14) следует, что характер
частотной зависимости JP изменяется от l-1 в нижней части УКВ
диапазона до l на высоких частотах с
переходом через l° в окрестности 100 МГц,
объясняя экспериментально установленную закономерность.
Отметим еще один интересный эффект. В экспериментах обнаружено
[1], что для небольших городов при подъеме приемной антенны (но при h2<<H) уровень сигнала быстрее возрастает на высоких частотах. Для
крупных городов усиление незначительно и практически не зависит от l. Это явление нельзя объяснить дифракцией
волн на крышах зданий, в то время как в (1.14) зависимость <JР> от h2 определяется отношением (H - h2)/l, что дает эффект, аналогичный отмеченному в эксперименте.
Вместе с тем, при внимательном рассмотрении обнаруживаются
недостатки формулы (1.14). В основном они связаны с геометрооптическим
приближением при описании затенений и проявляются в расчетах при малых высотах
подвеса передающей антенны над городской застройкой. Завышенная оценка
ослабления сигнала по формуле (1.14) обнаруживается при высотах h1 - H ≤ 10…15 м, а при
расположении антенны на уровне крыш окружающих зданий h1=H расчетное значение
средней интенсивности обращается в нуль, что совершенно неприемлемо. В [6]
показано, что приближение физической оптики, если его использовать для описания
поля над городом, позволяет скорректировать высотную зависимость и получить
результаты, согласующиеся с экспериментом при малом подъеме антенны над
городом. Это достигается тем, что в расчетах учитываются дифракционное проникновение
радиоволн в область геометрической тени и связанное с этим перераспределение
энергии по высоте. Для этого используется приближение Гюйгенса - Кирхгофа.
В качестве поверхности виртуальных источников выберем поверхность S (рис.1.12),
перпендикулярную слою городской застройки, и поверхность бесконечной полусферы SR, опирающуюся на S и замыкающую объем,
внутри которого находится источник. При этом вкладом виртуальных источников,
расположенных на полусфере SR, можно пренебречь, так как
при стремлении ее радиуса к бесконечности он стремится к нулю.
Рис.1.12.
К расчету поля виртуальных источников
Поле в точке приема запишем в виде
,
где - поле на
поверхности S, рассчитанное в приближении
однократного рассеяния [35],
где – нормаль к поверхности S.
Для средней интенсивности рассеянного поля получим выражение [6]
(1.15)
где - корреляционная
функция поля на поверхности S при условии, что
источник расположен в точке .
Асимптотическая оценка (1.15) при условии r >> r приводит к следующему результату для средней интенсивности
рассеянного поля
(1.16)
обобщающему (1.14) на
весь интервал высот h1 ≥ H.
При формула (1.16) переходит в (1.14), в то же
время у верхней границы слоя городской застройки формула (1.16)
определяет зависимость средней интенсивности от дальности как степенную <JP> ~
r-n
с
показателем степени n = (2,5…3).
При связи между пунктами, когда обе антенны расположены на высоте
3…5 м от поверхности земли, экранирующее влияние строений проявляется в полной
мере. На расстоянии
Рис.1.13.
Ослабление сигнала в центральной части города
близко к 4. По другим
источникам диапазон возможных значений n колеблется от 2 до 6 [41]; в [42] приводятся данные об
экспоненциальном характере ослабления сигнала; но во всех случаях речь идет о
расстояниях не более 3…5 км, при этом результаты разных экспериментов
согласуются на расстояниях до
Анализ
теоретических и экспериментальных данных показывает [47], что при удалении
пунктов возрастает роль многократных отражений, а интенсивность сигнала в
среднем убывает экспоненциально. На расстояниях, превышающих 1,5…2 км, связь
обеспечивается распространением над городом боковой волны, энергия которой,
перераспределяясь по высоте за счет эффектов дифракции, проникает в слой
городской застройки. На частотах выше 100 МГц поле над городом у верхней
границы слоя городской застройки формируется, в основном, отраженными от зданий
волнами, на частотах около 30 МГц заметным становится вклад прямой волны,
испытавшей дифракцию на крышах ближайших зданий. К низко расположенному приемнику радиоволны
приходят, испытав вторичное отражение или дифракцию в его -окрестности (рис.1.14). Предполагая, что для
дециметровых и сантиметровых волн основной вклад в формирование поля в точке
приема вносят однократные отражения от зданий в -окрестностях передающей и приемной антенн и учитывая
явления дифракции в поле боковой волны над городом в приближении
Гюйгенса-Кирхгофа, представим среднюю интенсивность в виде двукратного
интеграла по поверхности S виртуальных источников [6]
Асимптотическая оценка интеграла
при приводит к следующему
результату
, (1.17)
Рис.1.14. К определению поля в точке
приема
который подтверждает степенной
характер зависимости <JР> от дальности r при связи между двумя подвижными
пунктами. Экспериментальные исследования ослабления сигнала проводились для
этого случая с помощью автоматизированного измерительного комплекса [48] на
частоте 410 МГц в условиях сильных затенений, исключающих возможность прямой
видимости между антеннами. Приемник, антенна которого располагалась на
высоте h2 =
Результаты эксперимента представлены на рис.1.15, где точками
нанесены значения уровня сигнала, усредненного по участку измерений, а кривой -
результаты расчета по формуле (1.17) при
Г = 0,1 и lв=1 м. За
0 дБ
Рис.1.15.
Результаты эксперимента на частоте 410 МГц
взято ослабление поля в
свободном пространстве на расстоянии
Радиальные или почти радиальные по отношению к передающей антенне
улицы оказывают сильное влияние на формирование поля УКВ в городе. Опыты
Реудинга в Нью-Йорке [23] на частоте 11,2 ГГц показали, что сигналы,
направления распространения которых параллельны направлениям улиц, обычно имеют
на 10…20 дБ большую мощность, чем сигналы, достигающие приемной антенны под
другими углами. Столь высокий уровень сигнала для радиальных улиц объясняется
наличием прямой видимости между корреспондентами, в результате чего сигнал в
точке приема формируется за счет сложения рассеянной компоненты поля и приемной
волны. В диапазоне дециметровых волн поперечный размер первой зоны Френеля для
прямой волны во многих случаях не превышает ширины современных городских
магистралей, поэтому в [34] предлагается расчет ослабления проводить по обычной
интерференционной формуле. Очевидно, что даже в этих случаях необходима оценка
той роли, которую могут играть отражения от зданий. В плотно застроенных
районах с высокими домами и узкими улицами, где наблюдается наибольшее
ослабление сигнала, такой подход не применим.
Статистическая оценка влияния отражений и затенений на ослабление
сигнала на радиальных улицах дается в [6]. При связи объекта с базовой
станцией, антенна которой поднята над городской застройкой, прямая и отраженная
от поверхности земли волны ослабляются при частичном затенении зданиями первой
зоны Френеля. Их интерференция определяет характерный высотный профиль уровня
сигнала с чередующимися максимумами и минимумами интенсивности, число и высоты
которых определяются высотами подъема антенн над поверхностью земли и длиной
волны излучения. Средняя интенсивность отражений оценивается формулой
(1.18)
при условии, что ширина
улицы l значительно меньше а высота подвеса
передающей антенны над крышами окружающих зданий превышает радиус первой зоны
Френеля . Первое слагаемое в (1.18) учитывает отражения от зданий,
расположенных вдоль улицы. На широких улицах их освещенность возрастает, а на
узких падает. Второе слагаемое учитывает отражения от зданий в - окрестности приемной антенны, при оно практически
совпадает с (1.14).
На рис.1.16 приведены результаты расчета ослабления поля на
радиальной улице шириной
Если приемная антенна
оказывается на площади или другом открытом месте, уровень сигнала возрастает
даже в отсутствие прямой видимости на передающую антенну [14, 39]. Это можно
объяснить двумя обстоятельствами. Во-первых, уменьшается затененность зданий,
расположенных за подвижным
пунктом по периметру площади, в результате возрастают число и интенсивность
отражений. Во-вторых, благодаря дифракции радиоволн на крышах зданий
уменьшается угол укрытия приемной антенны затеняющими зданиями, что повышает
уровень сигнала, приходящего к пункту от передающей антенны.
F, дБ
Рис.1.16. Результаты расчета ослабления поля на
радиальной улице
На условия связи в городе
существенное влияние оказывает рельеф местности [14, 15]. Оценка влияния холмистой и
наклонной местности на характер ослабления проводится на основе рекомендаций
МККР [31] либо на основе вычислений по компьютерным программам, учитывающим
результаты обработки экспериментальных данных, полученных при распространении
радиоволн над неровной поверхностью [49]. Отметим, что механизм формирования
поля в условиях города отличен от механизма для плоской поверхности земли,
поэтому перенос закономерностей, полученных на приземных трассах, на город со
сложным рельефом местности не совсем обоснован и требует дополнительного
исследования, учитывающего свойства городской среды распространения радиоволн.
По условиям
распространения радиоволн можно выделить три наиболее типичные ситуации.
Первый случай. Рельеф
местности обеспечивает условия незатенения r окрестности приемной антенны
относительно передающей антенны. Реально это могут быть плавные уклоны либо
подъемы местности, пологие неровности со средним уклоном, либо подъемом,
склоны, обрывы, холмы, впадины и т.д.; требуется только, чтобы выполнялось
условие незатенения окрестности
приемной антенны. В этом случае для расчета <J> можно
пользоваться результатами, полученными для «квазигладких»
районов, пересчитав рельеф местности в высоту подвеса передающей антенны. Один
из примеров приведен на рис.1.17.
Второй
случай. За счет рельефа местности выделяется группа зданий, расположенных на
возвышенности, которые не затеняются ни относительно передающей антенны, ни
относительно приёмной антенны. Поле,
Рис.1.17.
Пример рельефа местности
Рис.1.18.
Примеры рельефа местности
Рис.1.19.
Пример рельефа местности
рассеянное этими
зданиями, необходимо рассчитывать по формуле (1.12) в каждом конкретном случае
для заданной формы возвышенности и плотности ее застройки. Отметим, что при
расчете реальных городских трасс распространения УКВ этот случай представляется
наименее вероятным, т.к. его реализация возможна в двух ситуациях:
а) возвышенность имеет достаточно большие размеры, поэтому здания
ближайшего окружения приёмной антенны не перекрывают линию прямой видимости
между приемной антенной и зданиями, находящимися на возвышенности (рис.1.18.). При этом должно выполняться условие
.
Если даже принять, что
ближайшее к приемной антенне здание находится на расстоянии r =, то при H
=
б) вблизи приемной антенны нет зданий, т.е. он
находится на площади или другом открытом месте, и возвышенность просматривается
достаточно хорошо (рис.1.18.б). При этом основной вклад в поле в точке приема
будет давать поле прямой волны, а не поле, рассеянное зданиями на
возвышенности. Действительно, даже если пренебречь ослаблением, связанным с
тем, что отраженные волны проходят большее расстояние, их потери при отражении
составляют в среднем 10…15 дБ (по сравнению с прямой волной). С учетом же
ослабления за счет расстояния вклад рассеянного поля будет еще менее заметным.
Третий случай. Подвижный пункт находится в зоне геометрической
тени, создаваемой неровностями рельефа, расположенными вдоль трассы,
соединяющей приемную антенну с передающей антенной. В этом случае поле в точке
приема можно рассматривать как поле, создаваемое вторичными переизлучателями в -окрестности подвижного пункта,
которое дифрагирует на неровности рельефа, и записать для него выражение вида
<Jg> = <JP>/F, где <JP> – средняя
интенсивность рассеянной компоненты поля для «квазигладкой» поверхности,
определяемая по формуле (1.16),
F – дифракционный множитель ослабления.
Если препятствие возможно аппроксимировать клином (рис.1.19), то
расчет множителя F можно выполнить по
известной методике [15].
Более реальная ситуация, – когда неровность рельефа, на которой
происходит дифракция, является плавной. В этом случае методику оценки влияния
рельефа местности на уровень поля в городе можно разбить на следующие этапы:
1) аппроксимация реального
рельефа, полученного по топографичес-кому плану района, полиномом вида
2) вычисление радиуса
кривизны аппроксимирующей кривой по формуле
;
3) расчет эффективной высоты подвеса передающей антенны
относительно аппроксимирующей поверхности
где x2 и h1 - коэффициенты
передающей антенны, х' может быть найдено из уравнения нормали к
аппроксимирующей кривой, которой принадлежит точка (x2 и h1), имеющей вид
а h' определяется из
аппроксимационного полинома при х = х';
4) вычисление средней
интенсивности рассеянного поля для плоской поверхности
при h1 = h1эф;
5) расчет дифракционного
множителя ослабления Ғ для
высот
расположения точек: (h1эф - H) и 0 на сфере радиусом aэ.
Экспериментальные исследования влияния рельефа на поле УКВ в
городе проводились с помощью автоматизированного измерительного комплекса [50],
антенна которого располагалась на высоте h2 =
Рис.1.20.
Средний высотный профиль рельефа
Эксперимент на частоте 166 МГц проводился в меньшем объеме, чем на
частоте 410 МГц, поскольку многое в анализируемых статистических
закономерностях является общим для двух частот. Целью эксперимента на частоте
166 МГц являлось подтверждение установленных закономерностей на другой частоте
и при другом высотном расположении приёмной антенны (h2 =
Экспериментальные (точки)
и расчетные (кривые 1, 2) результаты приведены на рис.1.21.б. Отметим, что
подъем приемной антенны на
Рис.1.21. Сравнение результатов
расчета с экспериментом
В целом можно отметить,
что расчеты, проведенные по предложенной методике, учитывающей параметры и
характер застройки, а также рельеф местности, согласуются с результатами
экспериментальных исследований.
Разработанная методика позволяет
перейти к построению радиокарт городских районов [47, 51]. На рис.1.22 представлена
радиокарта одного из городских районов, построенная для частоты 410 МГц.
Сплошными изолиниями показаны результаты расчета по описанной методике,
экспериментально полученное среднее ослабление интенсивности поля по ансамблю в
2300 точек нанесено цифрами в кружках. За 0 дБ взято ослабление поля в свободном
пространстве на расстоянии
Рис.1.22. изолинии уровней поля на частоте 410
МГц
Условия распространения УКВ в пригороде лучше, чем условия
крупного города, так как пригороды застроены домами меньшей этажности, в основном
- это одно - и двухэтажные здания, и плотность застройки меньше, чем в крупных
городах. При этом чаще реализуются условия, близкие к условиям распространения
радиоволн над «неровной» поверхностью земли. Однако работ, посвященных
исследованию ослабления сигнала в пригороде, мало; имеющиеся результаты
фрагментарны и трудно сопоставимы. Так, в работе [12] приведены значения
поправочного “коэффициента усиления” сигнала в пригороде относительно уровня
сигнала в городе и отмечено, что его величина не зависит от высоты подвеса
антенны базовой станции и расстояния между антеннами, условия распространения
улучшаются с увеличением частоты сигнала: значения “коэффициента усиления”
изменяются от 6 дБ на частоте 100 МГц до 14 дБ для 3 ГГц [14].
Однако Реудинг [23] приводит для частоты 11,2 ГГц величину
порядка 10 дБ.
Обработка топографических карт показывает, что средние параметры
застройки для пригородных районов можно примерно оценить значениями , т.е. средняя горизонтальная дальность прямой видимости
достигает км, что в пять-десять раз превышает величину для городских районов. При таких характеристиках застройки
велика вероятность приема сигнала непосредственно от источника, определяемая
формулой
1.3.
Пространственно-временные изменения уровней поля
Первые экспериментальные
исследования пространственных и временных изменений уровня сигнала в условиях
города привели к выводу, что городской канал распространения радиоволн является
локально стационарным [11], а пространственные замирания уровня сигнала имеют
двойную природу [14]. Быстрые замирания являются результатом сложения полей
отдельных волн многолучевого поля и имеют интерференционный характер, их
пространственный масштаб колеблется от половины длины волны до двух-трех длин
волн [6, 14]. На рис.1.23 показаны записи сигнала, принимаемого движущимся вдоль
улицы абонентом на частоте 410 МГц при высоте подвеса антенны
Рис.1.23. Пример записи уровней
сигнала на частоте 410 МГц
в котором определяют две пространственно разнесенные точки
наблюдения, выбранные на одинаковой высоте у поверхности земли, а ломаные
скобки < > означают статистическое усреднение по ансамблю реализаций
пространственного распределения поля в различных городских районах. Приведенные
зависимости получены на основе обработки экспериментальных данных на частоте 88
МГц. Чтобы определить медленные изменения уровней поля, интерференционные
замирания сглаживались усреднением по скользящим интервалам протяженностью 6…8
Рис.1.24.
Графики структурных функций
длин волн. Графики
структурных функций имеют два четко выраженных максимума. Первый - при разносе
точек на 15…20 м - можно соотнести со средним размером «освещенных» зон,
обусловленных просветами между зданиями. Второй - при разносе 80…100 м -
отражает характерное для городских условий чередование «освещенных» и «теневых»
зон и соотносится со средней длиной зданий в районе города.
В условиях многолучевого распространения волн в городе
естественным является предположение о статистической независимости отраженных
волн. Напряженность поля, создаваемая каждой волной в точке наблюдения, зависит
от множества случайных факторов: формы и электрических свойств отражающей
поверхности, ее неоднородности, ориентации в пространстве. Это позволяет
рассматривать напряженность поля как комплексную случайную величину,
распределение которой близко к нормальному, а для амплитуды сигнала
предположить рэлеевский или райсовский характер распределения [14, 39].
Экспериментальные исследования флуктуаций амплитуды принимаемого
сигнала, выполненные в разнообразных условиях по типам городских застроек и
рельефу местности, охватывают широкий диапазон частот от 50 МГц до 20 ГГц. При
движении приемного объекта в районах с примерно постоянной плотностью застройки
50% замираний оказываются на 9…10 дБ ниже среднего уровня сигнала, глубокие
замирания достигают 24…35 дБ, при этом глубина замираний практически не зависит
от частоты излучения. Исследование распределения амплитуды сигнала,
принимаемого движущимся объектом в условиях города, осложняется статистической
неоднородностью пространственных вариаций поля, причина которой, во-первых, в
ослаблении поля при увеличении расстояния между антеннами и, во-вторых, в том,
что затенения, порождаемые зданиями, вызывают медленные изменения уровня сигнала
и характера интерференционных замираний. В [6] измерения проводились
по нескольким десяткам участков длиной
Рис.1.25. Гистограмма
распределения измеренных значений амплитуды сигнала на частоте 410 МГц
Медленные пространственные замирания и средние значения амплитуд
по участкам протяженностью 20…30 м подчиняются логарифмически-нормальному
распределению со стандартным отклонением, зависимым от рельефа местности и типа
городской застройки. Стандартное отклонение медленных замираний составляет 4,5
дБ для центральных районов города и 0,5…1 дБ для пригорода. Для городов с
пересеченным рельефом местности стандартное отклонение больше, чем для городов
с равнинной местностью, и может достигать 10 дБ [29, 54].
Длительность интерференционных замираний сигнала пропорциональна его частоте,
но при изменении направления движения приемного объекта от перпендикулярного к
трассе до радиального меняется примерно в 2 раза [14, 23].
Распределение амплитуд на различных временных задержках также, в
основном, близко к рэлеевскому, хотя в некоторых случаях наблюдались
отклонения, размах флуктуаций в разных реализациях менялся от 18 до 30 дБ [55, 56].
Временные флуктуации сигнала, принимаемого неподвижным объектом,
достигают 20…25 дБ при прохождении «возмущающего» объекта (транспорт, люди и
т.д.) в “ближней зоне” приемной антенны и не превышают 3 дБ даже в случае
интенсивных транспортных потоков, находящихся в дальней зоне [56].
Рассеяние волн на деревьях вызывает быстрые флуктуации сигнала, глубина и
частота которых зависят от скорости и направления ветра и планировки городского
района. С увеличением частоты излучения (сантиметровый диапазон) глубина
временных флуктуаций возрастает и в зависимости от метеорологических условий
может достигать 10…13 дБ.
Для статистического описания вариаций поля в условиях города в
рамках корреляционной теории необходимо знать функцию пространственной
корреляции и соответствующую ей спектральную плотность. Эти вопросы
рассматривались во многих работах как на основе анализа экспериментальных
данных, так и с применением расчетных методов [6, 14].
В условиях города пространственное распределение напряженности поля практически
однозначно определяет свойства временных флуктуаций сигнала, принимаемого
движущимся объектом, так что масштабы временной и пространственной корреляции
связаны через скорость и время движения. Например, время корреляции быстрых
интерференционных замираний зависит от частоты излучения и в диапазоне от 50
МГц до 11 ГГц при скорости перемещения 20…40 км/ч изменяется от 40…60 до 3…5
мс, медленные вариации при тех же условиях имеют интервал корреляции от 0,5 до
3…5 с, почти неизменный во всем диапазоне частот.
В ряде работ выполнен расчет спектральной плотности сигнала на
выходе приемной антенны подвижного пункта в предположении, что число
интерферирующих волн велико и их направления прихода равновероятны [14].
Полученные спектры симметричны относительно несущей частоты и занимают полосу,
ширина которой равна удвоенной частоте Доплера.
На рис.1.26 показана форма энергетического спектра флуктуаций
вертикальной электрической компоненты поля при условии, что передаваемый сигнал
вертикально поляризован. Здесь wo
- несущая частота, u -
скорость
движения пункта, wm = wou/c определяет максимальное
доплеровское смещение. Соответствующий коэффициент корреляции представлен
функцией Бесселя Iо(wmt) и показан на рис.1.27.
Полученные на этой основе спектральные плотности огибающей сигнала,
принимаемого движущимся объектом, в целом согласуются с экспериментальными
данными. Отметим несколько моментов. Во-первых, очевидно, что модель призвана
описывать лишь быстрые интерференционные замирания; медленные, порождаемые
затенениями, она не учитывает, рассматривая их как отклонения от
стационарности, ограничивающие область применимости модели. Во-вторых, число
волн, приходящих в точку приема, скорее всего, невелико [1, 2, 47].
По-видимому, оно может рассматриваться как случайная пуассоновская величина со
средним значением 3…5. Направления прихода волн в точку приема также не
равновероятны. На это указывает как статистический анализ многолучевой
структуры поля, так и непосредственные измерения реализации углового
энергетического спектра приходящего излучения, выполненные при использовании
антенн с высоким угловым разрешением [6, 57].
Рис.1.26. Форма
энергетического Рис.1.27. Коэффициент корреляции
спектра флуктуаций для
вертикально
вертикальной составляющей поляризованного
поля
поля
В [58, 59] выполнен расчет корреляционной
функции и спектральной плотности флуктуаций сигнала в системе радиосвязи с
подвижными объектами в предположении, что распределение углов прихода волн
имеет сосредоточенный унимодальный характер. Это, естественно, привело к
появлению ряда новых свойств корреляционной функции и спектра сигнала, в
частности к зависимости от направления и скорости перемещения подвижного
объекта. Однако адекватность модели реальному каналу внутригородской связи
осталась неясной, поскольку сделанные предположения не связаны с анализом
процесса многолучевого распространения и не учитывают каких-либо характеристик
городской среды. Наконец, в целом лучевая картина - это локальная
характеристика поля, изменяющаяся при перемещении подвижного объекта.
Статистическое описание затенений в городской среде позволяет более
последовательно выполнить расчет пространственной корреляционной функции.
Предположим, что антенна стационарного передающего пункта
расположена в точке . Наземный приемный объект, двигаясь со скоростью под углом jо к направлению на
передатчик, за время t переместился из точки в точку , l = ut. В приближении однократного рассеяния функция пространственной
корреляции поля может быть представлена выражением [48]
(1.19)
аналогичным формуле
(1.13) для средней интенсивности и записанным в тех же обозначениях. Два новых
множителя в подынтегральной функции имеют ясную физическую интерпретацию: - вероятность того,
что направленный отрезок полностью освещается волной, рассеянной в точке ; - множитель,
описывающий изменение фазового состояния рассеянной волны вдоль - волновой вектор.
Если расстояние между точками А
и В значительно превышает , то асимптотическая оценка интеграла позволяет представить
корреляционную функцию в виде произведения [8]
(1.20)
в котором множитель
(1.21)
описывает медленные
пространственные изменения напряженности поля, порождаемые затенениями, а
функция
(1.22)
представлена интегралом
Фурье от нормированного углового энергетического спектра излучения, приходящего
к приемному пункту, и является коэффициентом корреляции быстрых
интерференционных замираний.
Функцию можно рассматривать как коэффициент
пространственной когерентности поля, поскольку ее значение определяет
вероятность того, что при перемещении точки наблюдения в направлении jо на расстояние l останутся неизменными
число и направление прихода волн. Как следует из (1.21), средний интервал
когерентности в направлении поперек трассы (j0 = 90°, 270°) имеет наименьшее значение
. (1.23)
При разносе точек вдоль трассы интервал когерентности определяется
не функцией (1.21), а экспоненциально убывающей вероятностью прямой видимости,
так что при jо = 0°, 180° масштаб достигает максимального значения, не превышающего . Оценка (1.23) среднего интервала когерентности согласуется
с экспериментальными данными. Первый максимум структурной функции, график
которой приведен на рис.1.24, соответствует разносу точек 18…20 м. При
плотности v = 90 км-2 на
расстоянии
Рассмотрим функцию , определяемую интегралом (1.22).
Её вид зависит от структуры многолучевого поля, складывающейся в -окрестности подвижного приемного объекта. В [48] приводятся
данные о том, что при расположении антенны стационарного объекта на уровне или
ниже крыш окружающих зданий (h1 ≤ H) максимум интенсивности
излучения, приходящего к приемному объекту, совпадает с направлением на
передатчик. По мере подъема антенны стационарного объекта возрастает
освещенность зданий, расположенных за подвижным приемником, и происходит
перераспределение интенсивности отраженных от зданий волн от направления на
передатчик к противоположному ему. Подъем антенны передатчика на 0,1 h над городской застройкой приводит к тому, что интенсивности с
указанных направлений в среднем угловом энергетическом спектре практически
сравниваются, а при подъеме на 0,5 h максимальная мощность
рассеянного излучения в среднем регистрируется с направления, противоположного
направлению на передатчик. В первом случае (h1 ≤ H) при асимптотическая оценка приводит к интегралу Фурье вида [8]
, (1.24)
где . Разложение (1.24) по
функциям Бесселя
позволяет получить ряд
члены которого при и быстро убывают с
ростом номера п, что делает приведенную формулу удобной для расчетов.
На рис.1.28 приведен пример экспериментальной зависимости модуля
коэффициента корреляции R2(l, jо) от величины kl на расстоянии
Рис.1.28. Зависимость модуля коэффициента корреляции R2(l, j0) от величины kl
Вертикальными линиями
указан 90%-ный доверительный интервал (общий объем выборки ~300 значений).
Пунктиром указана расчетная зависимость для j0 = 75° и параметра городской застройки gо
= 6
км-1,
что соответствует условиям эксперимента. На рис.1.29 приведены
экспериментальная и расчетная зависимости аргумента функции корреляции для
указанных условий.
Рассмотрим второй случай h1>H, когда антенна стационарного пункта поднята над городской
застройкой. На расстояниях при подъеме антенны
передатчика над крышами окружающих зданий ослабевает зависимость R2 (l, jо) как от r - расстояния между
объектами, так и от высоты подвеса h1 антенны стационарной
станции, и выражение для него приобретает вид [6]
. (1.25)
Рис.1.28. Сравнение зависимостей аргумента функции корреляции
Расчеты показывают, что
формула (1.25) удовлетворительно описывает интерференционные замирания на
дальностях от двух километров и более при h1 > l,2H в городских районах с
плотностью застройки v
= 90 км-2. На
рис.1.29. показана зависимость модуля |R2(l, jо)| коэффициента
корреляции (1.25) от kl при разносе точек под
углом j0 = p/2 (кривая 1), p/3 (2),
p/4 (3) и p/6 (4)
к трассе. Графики зависимости Arg R2(l, jо) от величины kl при тех же условиях приведены на рис.1.29.б.
Необходимо отметить некоторое увеличение отношения масштаба
когерентности к интервалу корреляции интерференционных замираний, но эти
изменения не являются существенными.
Наконец, отметим еще один результат, который несложно получить из
(1.24) и (1.25). Рассчитаем средний по направлениям движения наземного объекта
интервал корреляции быстрых интерференционных замираний, считая любые значения jо Î (0, 2p) равновероятными. После усреднения по j0 получим
(1.26)
то есть коэффициент корреляции полностью определяется разносом
точек наблюдения и частотой
излучения и не
зависит ни от расстояния
и
Рис.1.29. Зависимость модуля (а)
и аргумента (б) коэффициента
корреляции при разных значениях j0
высоты подъема
стационарного объекта, ни от плотности и этажности городской застройки.
Спектр пространственных частот быстрых интерференционных
замираний, нормированный на среднюю интенсивность <I> и определенный как
преобразование Фурье коэффициента корреляции R2(l, jо)
отличен от нуля в
интервале . При он определяется выражением [6]
(1.27)
где , . На рис.1.30
приведены графики спектров пространственных частот быстрых
интерференционных замираний, рассчитанные по (1.27). Кривые 1, 2, 3 и 4
соответствуют ориентации баз j0 = p/2, p/3, p/4 и p/6. По горизонтальной оси указаны
значения безразмерного отношения u/k. Асимметрия спектра, проявляющаяся при отклонении угла ориентации
базового участка от p/2, объясняется тем, что
направления прихода волн многолучевого поля не являются равновероятными:
угловой энергетический спектр унимодален и симметричен относительно моды,
совпадающей с направлением на излучатель.
При h1 > H,
когда антенна передатчика поднята над городской застройкой, в результате чего
возрастает освещенность -окрестности наземного приемного пункта, нормированный спектр
пространственных частот определяется выражением [6]
(1.28)
Его вид показан на
рис.1.31 для той же ориентации баз, что и в предыдущем случае. Для него также
характерна асимметрия, однако максимум внутри интервала хÎ(-1, +1) отсутствует.
В условиях города пространственное распределение напряженности
поля практически полностью определяет свойства временных флуктуаций сигнала,
принимаемого подвижным приемным объектом. Если скорость его движения в
направлении передатчика равна u, то энергетический спектр
флуктуаций принимаемого сигнала связан со спектром
пространственных частот (, соотношением
(1.29)
Излучение в условиях города проходит от источника к приемнику
несколькими различными путями. В результате формируется многолучевая структура
поля. При движении наземного пункта связи сигнал в каждом луче приобретает
доплеровское смещение по частоте, определяемое частотой излучения w0, скоростью движения и направлением
движения относительно той точки, в которой сигнал отразился. Отсюда следует,
что спектр временных флуктуаций сигнала, принимаемого подвижным объектом,
является спектром доплеровских частот, который отличен от нуля в интервале , ограниченном максимальными доплеровскими смещениями частоты
.
Внутри этого интервала выполняется соотношение (1.29),
которое в двух наиболее типичных ситуациях может быть записано с учетом (1.27)
и (1.28) при замене на отношение w/wдм.
Экспериментальные данные по исследованию спектральных свойств
сигналов подвижных УКВ-источников в условиях города практически полностью
подтверждают теоретические результаты [6]. На рис.1.32 представлен спектр
доплеровских частот сигнала, принятого от источника излучения с частотой 410
МГц при его движении со скоростью u =
На рис.1.33 показан спектр доплеровских частот для случая, когда
приемная антенна располагалась на крыше пятиэтажного здания, а антенна
передатчика - на высоте
Рис.1.30. Спектр пространственных частот быстрых
интерференционных замираний
Рис.1.31.
Нормированный спектр пространственных частот
Рис.1.32.
Спектр доплеровских частот сигнала
Рис.1.33.
Спектр доплеровских частот сигнала
Движение транспорта вблизи наземного пункта и рассеяние волн на
движущихся объектах могут расширять спектр доплеровских частот [60]. По этой
причине в исключительных случаях спектр может расширяться в два раза и более.
Характеристики распространения импульсных сигналов метрового
диапазона, приведены в [61,62].
1.4. Модель
сигнала при многолучевом
распространении УКВ в
городе
Многолучевой канал городской радиосвязи достаточно полно
характеризуется передаточной функцией, которая при движении одного из корреспондентов
является случайной функцией времени. Моделирование ее реализаций на
фиксированной частоте предоставляет возможность численного исследования
вопросов прохождения по каналу сложных сигналов с произвольными по ширине и
форме спектрами. Рассмотрим предпосылки, необходимые для моделирования
реализаций передаточной функции при связи на УКВ базовой станции с движущимся в
городе абонентом.
Основными факторами, порождающими случайную модуляцию сигнала,
излучаемого движущимся в городе передатчиком, являются многолучевость и
затенения. В каждый момент времени в точку приема приходит несколько лучей. Их
число М случайно и меняется во времени при движении передатчика. Время
жизни каждого луча определяется временем незатенения отражающего объекта
относительно передатчика и приемника и также является случайным. Различные
условия отражения и пути распространения сигнала в лучах придают случайный
характер значениям амплитуды и фазы, а перемещение передатчика в определенном
направлении вызывает в лучах различные доплеровские смещения частоты, зависимые
от того, насколько направления на отражающие объекты отличаются от направления
движения.
Учитывая все эти обстоятельства, в случае непрерывного
монохроматического излучения на частоте w представим сигнал на входе приёмника суммой сигнала в М лучах [6]
(1.30)
где aт – комплексная амплитуда, определяемая условиями отражения, tm характеризует задержку
сигнала в m-м луче, Ωm – доплеровское смещение частоты
(1.31)
определяемое скоростью u движения передатчика и
углом jm между
вектором скорости и направлением на
отражающий объект. Функция (t) характеризует время
существования m-го луча и имеет вид [6]
где и - моменты его появления и исчезновения.
Число лучей М, реализовавшихся за время наблюдения,
является пуассоновской случайной величиной, закон распределения которой имеет
вид
(1.32)
с параметром h, равным среднему числу
лучей, наблюдаемому в единицу времени. Из результатов экспериментальных и
теоретических исследований, приведенных в предыдущих разделах можно получить
оценку среднего числа лучей N, приходящих одновременно в
точку наблюдения, а также среднего расстояния , при перемещении на которое лучевая структура сигнала
полностью меняется. При скорости движения передатчика u пройденный за время
наблюдения путь равен uT, следовательно, лучевая
структура сигнала сменится в среднем раз, а так как одновременно в среднем наблюдается лучей, то за время
Т среднее число лучей окажется равным , а за единицу времени
Моменты появления , образуют за время наблюдения простой и ординарный поток, в котором интервалы между моментами
появления двух лучей независимы и распределены экспоненциально с плотностью
вероятности
(1.33)
при этом среднее значение
.
Причины, ограничивающие случайное «время жизни»
каждого луча те же, что и
причины, определяющие интервал времени между
моментами появления двух последовательно возникающих лучей - это движение
передатчика и затенения. Следовательно, tm является также экспоненциально распределенной случайной
величиной, среднее значение которой
. (1.34)
Рассмотрим характеристики сигнала в m-м луче. «Время жизни»
луча ограничено, поэтому прошедший по нему сигнал будет принят в точке
наблюдения как высокочастотный импульс длительностью tm, комплексная амплитуда
которого равна aт, а частота смещена
относительно несущей на Ωm. По условиям отражения комплексная амплитуда в луче может
рассматриваться как нормальная случайная величина, в среднем по ансамблю лучей
равная нулю. Ее действительную и мнимую части - квадратурные составляющие - в
этом случае можно определить как независимые нормальные случайные величины с
нулевыми средними значениями и равными дисперсиями.
Доплеровское смещение частоты сигнала в m-м луче - случайная
величина Ωm, для которой плотность
вероятности определяется распределением точек отражения в окрестности
движущегося передатчика и его скоростью. В случае, когда антенна базовой
станции приподнята над городской застройкой (z>0,2), плотность вероятности , характеризующая угловое распределение отражателей, имеет
вид [6], аналогичный
(1.35)
и практически не зависит
от каких-либо параметров, кроме угла jо, определяющего направление движения
передатчика. Учитывая (1.35), получим при указанных условиях следующее
выражение для плотности вероятности доплеровского смещения Ω
,
(1.36)
где - максимальное
доплеровское смещение частоты, возникающее в том случае, когда точка отражения
находится на траектории движения передатчика.
Если антенна базовой станции установлена на уровне крыш или ниже,
плотность вероятности описывается как
плотность распределения Коши на окружности [6]
,
параметр b которого при равен . В этом случае плотность вероятности
доплеровского смещения имеет вид, аналогичный (1.27)
(1.37)
где
Формула (1.30) определяет случайную передаточную функцию Х(t, w) городского канала связи, а (1.32)…(1.37) - распределения ее
случайных параметров при связи базовой станции с движущимся абонентом. Если
передатчик излучает модулированный сигнал S(t), спектр которого S(w), то в точке приема
(1.38)
где Ωm для каждой спектральной составляющей
пропорционально ее частоте w. Предположим, что сигнал
S(t) - узкополосный, занимает
полосу частот причем . В этом случае различием доплеровского смещения спектральных
составляющих можно пренебречь, что сигнал в
m-м луче в целом
приобретает доплеровское смещение, определяемое центральной (несущей) частотой wo.
Кроме того, в условиях города дальность связи с подвижным объектом, как
правило, невелика, так что для узкополосного излучения в диапазоне УКВ можно
считать и тем более и пренебречь
соответствующими фазовыми изменениями. Заметим также, что случайная задержка в
любом m-м луче на много порядков
превышает период несущего колебания и вызывает равновероятное на интервале [0,
2p] смещение фазы волны. С
учетом сказанного (1.38) можно преобразовать к виду
,
откуда следует, что то есть влияние канала радиосвязи на передаваемый
узкополосный сигнал можно рассматривать как случайную модуляцию, описываемую
передаточной функцией на частоте несущего колебания.
1.5. Деполяризация
излучения
Сложным интерференционным полям, формирующимся в городе в условиях
многолучевого распространения УКВ, свойственны крайне нерегулярные
пространственные изменения поляризации. Независимо от поляризации излучаемой
волны поле в каждой точке пространства оказывается эллиптически поляризованным,
при этом параметры эллипса и его ориентация существенно изменяются на
расстоянии, определяемом масштабом интерференционных замираний [63, 64]. В этом
случае направление нормали к поляризационному эллипсу не связано с
направлением
распространения какой-либо волны и имеет стохастический характер, так как
зависит от числа волн, приходящих в точку наблюдения, их направлений прихода,
амплитуд и фаз. Эти обстоятельства затрудняют поляризационные измерения и
фактически делают непригодными традиционные методы, требующие совмещения
фазовых плоскостей сигнала и антенны поляриметра. Более общий метод измерения
локальных поляризационных характеристик предложен в [64]. Устройство, его
реализующее, включает в себя три линейных ортогональных вибратора с
совмещенными фазовыми центрами и обеспечивает регистрацию как
амплитуд
трех ортогональных составляющих напряженности электрического поля
, так и разностей фаз между ними [6]. Антенная система поляриметра
определяет прямоугольную декартову систему координат, характеризуемую
тройкой
единичных векторов , относительно которой
фиксируется
ориентация поляризационного эллипса.
Общий случай расчета параметров эллиптической поляризации поля,
основанный на представлении напряженности E(t) в данной точке пространства тремя взаимно ортогональными компонентами
Ei(t),
рассматривался
в [6, 64]. Пусть на частоте w компоненты поля
представляются в виде
(1.39)
Определим соответствующие
квадратурные составляющие
и, принимая их за координаты
квадратурных векторов и запишем равенство
Векторы и не зависят от времени
и определяют плоскость, в которой вращается вектор описывая своим концом за период колебания
плоскую кривую - эллипс. Векторное произведение определяет нормаль к плоскости эллипса, ориентацию которой относительно
выбранной координатной системы характеризуют азимутальный угол
(1.40)
и угол места q Î (0, p),
(1.41)
Полуоси эллипса равны:
где - интенсивность поля,
а коэффициент эллиптичности R,
определяемый отношением малой полуоси к большой равен
(1.42)
Очевидно,
случай N = 0 соответствует линейной
поляризации волны, a 2N = J -
круговой поляризации.
Угол β наклона большой полуоси
отсчитывается в плоскости эллипса от линии ее пересечения с горизонтальной
плоскостью, для которой нормаль совпадает с ортом .
Чтобы
записать расчетную формулу для β,
определим - аналог вектора , который можно рассматривать как нормаль к плоскости эллипса
поля , отличающегося от фазовым смещением на p/2 третьей компоненты
Через компоненты вектора и углы ориентации
плоскости эллипса j и q определяется угол наклона большой полуоси
(1.43)
где
Геометрические параметры эллипса R, j, q, β и интенсивность J
полностью определяют поляризацию поля в точке наблюдения, а их связь
с векторами и предоставляет
возможность расчета этих параметров по данным сравнительно простых измерений.
Действительно, определим высокочастотные сигналы
получаемые попарным
сложением ортогональных компонент Ei(t) при их дополнительном фазовом смещении d, кратном p/2. Амплитуды таких сигналов обозначим
(1.44)
где индексы выбраны так, что их расположение i, k, l образуется круговой
перестановкой чисел 1, 2, 3. Несложно убедиться в том, что компоненты вектора N
равны
Аналогичные выражения
определяют две компоненты вектора
,
третья же компонента совпадает с . Если к тому же учесть, что
а
равенство
позволяет получить
интенсивность любой компоненты поля, то оказывается, что амплитуды (1.44)
суммарных сигналов полностью определяют поляризацию поля в точке наблюдения.
Ненаправленный источник УКВ-излучения порождает в условиях города
многолучевое поле, которому свойственны унимодальность и симметричность
относительно моды статистического распределения углов прихода волн в точку
наблюдения. Такими же свойствами обладает средний по ансамблю трасс угловой
энергетический спектр регистрируемого в точке наблюдения излучения. Анализ показывает
[6, 14], что три ортогональные компоненты напряженности электрического поля оказываются
некоррелированными в точке наблюдения (при усреднении по ансамблю трасс), если
одну из координатных осей ориентировать вертикально, а другую направить на
передатчик. С другой стороны, если учесть, что на значения амплитуды и фазы
каждой из трех компонент поля влияет в условиях многолучевого распространения
большое число случайных факторов, то можно рассматривать каждую компоненту в
фиксированный момент времени как комплексную нормальную случайную величину. Но
тогда ортогональные компоненты в любой фиксированный
момент времени являются статистически независимыми нормальными случайными
величинами с нулевыми средними значениями и различными дисперсиями. Это
обстоятельство позволяет записать плотность совместного распределения шести
компонент квадратурных векторов и в виде произведения шести
одномерных нормальных плотностей вероятности. Переход от квадратурных
составляющих к параметрам характеризующим
поляризационное состояние поля в точке наблюдения, позволяет получить плотность
вероятности их совместного распределения, а затем и одномерные плотности вероятности
для каждого параметра [6].
Теоретические и экспериментальные исследования [6, 63] показали,
что поляризационные признаки, свойственные первичной волне, в разной степени
проявляются и на открытых, и на закрытых трассах, причем при вертикально поляризованной
первичной волне эти признаки проявляются несколько более определенно, чем при
горизонтальной поляризации. Обозначим дисперсию вертикальной компоненты , а горизонтальной и будем рассматривать
отношение при вертикальной
поляризации излучаемой источником волны и обратное отношение при горизонтальной
поляризации первичной волны. Обработка данных измерений по 109 участкам с общим
числом точек более 15000 на частоте 410 МГц показала [63], что в 98% случаев
отношение дисперсии кроссполяризованной компоненты к дисперсии основной не
превышало единицы как на открытых, так и на закрытых трассах. Соответствующие
гистограммы распределений показаны на рис.1.35: ) гистограмма распределения для открытых базовых участков -
радиальные по отношению к передатчику или пересекающие радиальные улицы, при
вертикальной поляризации излучения, б) то же для закрытых базовых участков, в)
гистограмма распределения для открытых базовых
участков при горизонтальной поляризации излучения, г) то же для закрытых
базовых участков. Определенной зависимости отношения от расстояния между пунктами на закрытых трассах установить
не удалось, хотя на открытых (радиальных) участках при
удалении от
Рис.1.35.
Гистограмма отношения дисперсий
источника отношение
дисперсии кроссполяризованной компоненты к основной постепенно в среднем возрастает. Можно
предположить, что на закрытых трассах степень деполяризации излучения
определяется, в основном, характером застройки -окрестности подвижного пункта, а не расстоянием между
пунктами.
На рис.1.36. (1 - 9) приведены в качестве примера гистограммы
распределений коэффициента эллиптичности R (1,4,7), угла b наклона большой полуоси
(2, 5, 8), угла q ориентации нормали к
плоскости поляризационного эллипса (3, 6, 9) при вертикальной линейной
поляризации источника излучения. Эти данные относятся к трем наиболее вероятным
подгруппам: и 0,6…0,8 (7 – 9).
Объем выборки для каждой из групп составляет N = 1750, 730, 1050 соответственно. Измерения
проводились на всех участках с интервалом
Рис.1.36. Гистограммы для вертикального
поляризованного поля
Аналогичные гистограммы для горизонтально поляризованного
излучения приведены на рис.1.37 (1 - 9). В этом случае близким к равномерному
оказалось распределение угла q, в то время как
распределения углов β и j приобрели характерную
для горизонтально поляризованного излучения форму. На передающем пункте
излучающий полуволновый вибратор был ориентирован под углом 60°…90° к радиотрассе. Это определило
положение максимумов распределения угла j, которые сформировались
при j = 90o и j = 270o.
Рис.1.37.
Гистограммы для горизонтального
поляризованного поля
На участках, для которых отношение оказалось равным
0,6…0,8, распределения поляризационных параметров β и j приобрели менее
регулярную форму, при этом в них наметились максимумы, смещенные
относительно основных на 90°.
На рис.1.38.(1 - 5) приведены данные измерений на открытых
радиальных участках при вертикальной (1 - 3) и горизонтальной (4 - 5)
поляризациях излучения. Форма распределения коэффициента эллиптичности R не зависит от ориентации эллипса в пространстве, поэтому
рис.1.38.1 в равной степени дает представление об аналогичных ему
распределениях коэффициента эллиптичности в условиях горизонтальной поляризации
источника. Рис.1.38.(2 - 3) по форме аналогичны рис.1.37.(2 - 3), отличаясь от
них лишь несколько большей сосредоточенностью распределений углов β и q. Столь же близко по
форме распределение β на
рис.1.38.4 к распределению на рис.1.37.2 в случае горизонтальной поляризации
излучаемого сигнала.
Рис.1.38.
Гистограммы результатов измерений на
открытых участках
Распределение j на рис.1.38.5 отличается
от распределения на рис.1.37.3 положением максимумов. На радиальных участках
они сместились на 90° относительно максимумов распределений на закрытых
трассах. Это вполне объяснимо, если учесть, что радиальные участки
располагались вдоль радиотрассы, в то время как закрытые участки в основном
были ориентированы поперек трассы. Слабо выраженный максимум при j = 180° распределения на
рис.1.38.5 указывает на преимущественно правое направление вращения вектора
напряженности электрического поля. Это можно объяснить особенностями планировки
городского района, в котором получена характеризуемая распределением рис.1.38.5
выборка.
Переотраженные и рассеянные волны в городе приводят к появлению
компоненты поля, ортогональной излученному полю. Степень деполяризации зависит
от типа трассы и ослабления на трассе. Максимальная развязка (отношение
напряженности поля первичной поляризации к напряженности поля ортогональной
поляризации) наблюдалась на трассах прямой видимости при расположении приемника
на уровне крыш и составляла в среднем 18 дБ с превышением в 90 и 10% мест
приема 9 и 29 дБ соответственно [65]. При опускании приемной антенны с 10 до
Таблица
1.1.
Зависимость
величины поляризационной развязки от условий распространения [66]
Тип
трассы (в пригородной зоне) |
Среднее
ослабление на трассе относительно поля свободного пространства, дБ |
Поляризационная
развязка в 90% случаев, дБ |
Прямая видимость Слабая дифракция Средняя дифракция Сильная дифракция |
7 26 40 27 |
18 13 4 2 |
Ранее исследованием
деполяризации сигналов в Узбекистане в условиях крупного города занимался
Ликонцев Д.Н. [67…69] Им установлено следующее:
– значения
коэффициентов деполяризации на радиальных улицах больше, чем на поперечных;
– значения
коэффициентов деполяризации в районах с большой плотностью застройки ниже, чем в
районах с малой плотностью застройки;
– наблюдается
увеличение значений коэффициентов деполяризации при увеличении высоты подвеса приёмной антенны;
– значения
горизонтальной и вертикальной составляющих напряженности электрического поля в г.
Ташкенте слабо коррелированы, в районах с большой плотностью застройки и на поперечных
улицах, что можно
объяснить большой интенсивностью переотраженных волн в этих районах и на поперечных
улицах, и не коррелированы в районах с малой плотностью застройки и на
радиальных улицах.
Однако результаты этих
измерений не давали полной картины о соотношении этих составляющих поля на
небольших участках, прилегающих к точкам измерений и в сельской местности.
В связи с этим было
решено провести дополнительное исследования распределения уровней вертикальной
и горизонтальной составляющих напряженности поля вблизи точек измерений. Для
этого необходимо провести измерения уровней вертикальной и горизонтальной
составляющих поля на участках вокруг точек измерений диаметром порядка
(1,0…1,5)l, определить средние значения этих
составляющих и и затем уже определить
значение «усредненного» коэффициента деполяризации , равного
, дБ. (1.45)
Найденные значения коэффициента
помогут оценить
возможность приема сигналов антенной с ортогональной поляризацией.
Для измерений был выбран
селективный микровольтметр STV-401 в комплекс с измерительной антенной (симметричным
вибратором). Высота подвеса измерительной антенны была выбрана равной
Анализ результатов
измерений составляющих напряженности поля и определения значений «усредненного»
коэффициента деполяризации показал:
– практически везде
наблюдалась ортогональная к основной составляющая напряженности поля;
– не наблюдалась явная
зависимость значения «усредненного» коэффициента деполяризации от плотности
застройки и расстояния;
– не наблюдалась
зависимость значений «усредненного» коэффициента деполяризации от расстояния в
сельской местности.
Следовательно, для приёма
сигналов метрового диапазона волн можно использовать антенны практически с
любой поляризацией.
1.6. Особенности распространения
радиоволн внутри зданий
Поле в жилых и производственных зданиях от каких-либо внешних
источников излучения представляет собой совокупность
дифракционно-интерференционного механизма распространения с учетом поглощения
электромагнитной энергии в различных конструкциях зданий.
При проникновении УКВ внутрь зданий напряженность поля заметно
уменьшается. Кроме того, на нижних этажах зданий величина напряженности поля
гораздо меньше, чем на верхних. Поле на седьмом этаже здания может составлять
от 6 до 40% поля на крыше, а в аналогичных условиях на первом этаже - от 3 до
7%. Измерения в двухэтажном здании показали, что на чердаке напряженность поля
на 9 дБ больше, чем на первом этаже, но на 10 дБ меньше, чем на крыше [65].
Проникновение поля в здание зависит от частоты. На частоте 570 МГц
разность величин напряженности поля внутри и вне здания на 7 дБ больше, чем на
частоте 55 МГц. Эта разность зависит также от строительного материала и
плотности застройки. Для бетонных зданий она на 10 дБ больше, чем для
деревянных, и в плотно застроенных районах она на 5 дБ больше, чем в редко
застроенных.
Среднеквадратичное отклонение при проникновении поля в здание при
логарифмически-нормальном распределении составляет 14 дБ.
По данным
[6], излучение на частоте 3,26 ГГц ослабляется на 12 дБ при прохождении сквозь
сухую кирпичную стену толщиной
Ниже приведены результаты
исследований распространения радиоволн внутри зданий, которые получены с
использованием документов ITU-R и данных экспериментов, проведенных в НИИ
Радио, в СНГ и за рубежом.
В реальных условиях
(внутри помещений) мощность принимаемого сигнала отличается от мощности сигнала
в свободном пространстве, где, потери растут в соответствии с законом L~rn, n = 2.
Внутри зданий показатель
степени n может быть больше, а иногда меньше
2. Из-за многолучевости сигнала в месте приема уровень сигнала от точки к точке
изменяется случайным образом. Построив статистическое распределение уровней в
различных точках поперек трассы, можно определить медианное значение потерь на
заданном расстоянии [70]
, (1.46)
где - реальная степень
возрастания потерь; Lg (m) – потери
сигнала в стенах и перекрытиях зданий, дБ; m
- число препятствий; f - частота сигнала, МГц.
Параметры и Lg (m) определяются
по результатам многочисленных экспериментов.
Распределение уровней
сигнала по местоположению в помещении, близко к логарифмически-нормальному
закону со стандартным отклонением σ, дБ. Поэтому потери сигнала,
превышаемые в М% мест площади помещения,
, (1.47)
где К - коэффициент, определяемый из статистического распределения
уровней, нормированного к стандартному отклонению. В частности, при М = 50%, К = 0; при М = 90%, К = 1,3; при М = 99%, К = 2,35.
Необходимо учесть
ослабление из-за поглощения излучаемой мощности и экранирования антенны телом
оператора (абонента). Это ослабление на частоте 460 МГц составляет 22 дБ, а на
частоте 930 МГц - 19 дБ [71] при
расположении антенны в непосредственной близости к абоненту.
Величина различна для разных
помещений и изменяется от частоты. Это положение подтверждают данные табл. 1.5
документа ITU-R [71].
Таблица 1.5
Значения степени возрастания потерь
Частота, ГГц |
Значения |
||
Жилое помещение |
Офис |
Помещение для торговли |
|
0,9 |
- |
3,3 |
2,0 |
1,2…1,3 |
- |
3,2 |
2,2 |
1,8…2,0 |
2,8 |
3,0 |
2,2 |
Примечание. В тех
случаях, когда значение не дано, оно может
быть взято таким же, как для офиса.
Однако эксперименты [70]
показали, что значение показателя степени определяется, прежде всего, видом
трассы внутри помещения. Условно их целесообразно подразделять следующим
образом:
-трассы с прямой
оптической видимостью: а) в пустом помещении;
б) при наличии в помещении большого
числа мешающих металлических предметов (шкафы, станки и т.д.);
-трассы с прямой
радиовидимостью;
-закрытые трассы.
На трассах с прямой
оптической видимостью линия визирования
не перекрывается никакими препятствиями и не касается их. На трассах с
прямой радиовидимостью сигнал принимается в главном направлении излучения
передающей антенны, но за "радиопрозрачной" (деревянной, каменной,
железобетонной) стеной. На закрытых трассах сигнал принимается в стороне от
главного направления излучения передающей антенны вследствие переизлучения
(дифракции, отражения, рассеяния) радиоволн от стен, потолка и пола.
Значения , полученные экспериментально, представлены в табл. 1.6, где m - число стен, через которые проходит
сигнал.
Таблица 1.6
Значения степени возрастания потерь
Частота, ГГц |
Значения для трассы |
||||
с прямой оптической видимостью |
с прямой радиовидимостью |
закрытой |
|||
а) |
б) |
m = 1 |
m = 2 |
|
|
0,05 |
3,0 |
- |
3,0 |
- |
- |
0,7 |
1,5 |
2,0 |
2,3 |
3,6 |
- |
2,7 |
1,5 |
1,7…2,0 |
4,0 |
5,0 |
10…13 |
Анализ результатов
экспериментов [70] показал, что после прохождения сигнала через стену, стоящую
на пути распространения радиоволны, значение возрастает. Иначе
говоря, уровень сигнала за препятствием начинает спадать быстрее. Установлено,
что поляризация волны и коэффициент усиления приемной антенны не оказывает
заметного влияния на скорость спадания. На трассах с прямой оптической видимостью,
например, в пустом коридоре, потери с расстоянием убывают медленнее, чем в
свободном пространстве (<2)
В общем случае значение меняется в широких
пределах. Так, по данным [72] внутри здания, в диапазоне частот 1000...2000 МГц
меняется от 3,8 до
6,5; на частоте около 1300 МГц - от 1,8 до 2,81. В случае прямой видимости = 1,49. Эти данные не
противоречат приведенным в табл. 1.6.
В документе ITU-R [71] приводится
таблица (см. табл. 1.7) для определения потерь сигнала Lg(m),
проникающего через m перекрытий
(этажей) здания, в зависимости от его частоты.
Точное определение потерь
сигнала, проходящего через стены или перекрытия, не так просто, как может
показаться на первый взгляд. Это связано с тем, что в месте приема наблюдаются
значительные пространственные флуктуации (максимумы и минимумы) сигнала, а
также тем, что сигнал частично проходит окольными путями вокруг препятствия.
Таблица 1.7
Частота, ГГц |
Потери сигнала при прохождении
через m перекрытий, дБ |
||
Жилое помещение |
Офис (контора, ведомство) |
Помещение для торговли |
|
0,9 |
|
9 (m = 1) |
|
19 (m = 2) |
|||
24 (m = 3) |
|||
1,8…2 |
4m |
15+4 (m-1) |
6+3(m-1) |
Примечание: Имеются в
виду перекрытия в современных железобетонных зданиях.
Анализ результатов
экспериментального исследования показывает, что наименьшие потери сигнала
наблюдаются в сухих стенах и перекрытиях. Однако при увлажнении таких
препятствий потери резко возрастают. Кроме того, потери зависят от материала,
из которого сделано препятствие (дерево, кирпич, железобетон), и в некоторых
случаях от типа поляризации радиоволн. Частота сигнала и толщина препятствия
незначительно влияют на величину потерь. Измерения в сухой стене, [70]
сделанной из горизонтально положенных бревен, показали, что потери при
вертикальной поляризации (ВП) меньше, чем при горизонтальной (ГП) (см. табл.
1.8).
В мокрой деревянной стене
потери возрастали, как минимум, на 5 дБ, в сильно увлажненной стене на частоте
2,7 ГГц -достигали примерно 13 дБ при обеих видах поляризации.
Таблица 1.8
Частота, ГГц |
Потери в сухой стене, дБ |
||
Деревянной (бревна) ВП/ГП |
Кирпичной ВП и ГП |
Железобетонной ВП и ГП |
|
0,06 |
- |
- |
8,5…10 |
0,7 |
2/4 |
- |
6…12 |
1,0 |
- |
4,5…6 |
- |
2,7 |
3,5/8 |
5…8 |
6…12 |
Результаты измерения
потерь сигнала в сухих кирпичных (толщина 17…40 см) и железобетонных (толщина
15…40 см) стенах при ВП и ГП показаны в табл. 1.8. После полива кирпичной стены
потери возрастают примерно в 3…4 раза.
Результаты измерения
потерь в межэтажном перекрытии приведены в табл. 1.9. Установлено, что
поляризация радиоволн практически не влияет на величину потерь.
Таблица 1.9
Частота, ГГц |
Потери в межэтажном перекрытии, дБ |
Примечание |
0,05…0,1 |
10…15 |
Многоэтажное железобетонное здание
(офис) |
0,3 |
8 |
|
2,7 |
12 |
|
2,7 |
7,5 |
Кирпичное 3-этажное здание (офис) |
Внутри помещений уровень
сигнала изменяется случайным образом от точки к точке пространства. В документе
ITU-R [71] указывается, что вид статистического распределения уровней сигнала
близок к логарифмически-нормальному закону со стандартным отклонением σ.
Значения s для частот сигнала 1,8...2 ГГц из указанного
документа приводятся в табл. 1.10.
Эксперименты, проведенные
авторами в 18-этажном современном здании института (офиса) на частотах 0,7 и
2,7 ГГц, показали примерно такое же значение стандартного отклонения [70].
Таблица 1.10
Частота, ГГц |
Значения , дБ |
||
Жилое помещение |
Офис |
Помещение для торговли |
|
1,8…2 |
8 |
10 |
10 |
Кроме пространственных
флуктуации уровня наблюдаются изменения уровня сигнала во времени из-за
отражения, рассеяния и поглощения волн движущимися объектами (людей, дверей,
лифта и т.д.).
При распространении
радиоволн внутри помещения происходит деполяризация излучаемой радиоволны. У
окна, когда условия приема были близки к условиям приема в свободном
пространстве, коэффициент деполяризации составлял порядка 20дБ а при удалении
точки приема от окна внутрь помещения наблюдалась уменьшение значений
коэффициента деполяризации.
В работе [73] представлены результаты измерений напряженности поля
в помещениях зданий административного типа при размещении передатчика с
излучателем в этих же зданиях. Пространственная амплитудно-поляризационная
структура поля в объеме комнаты измерялась на частотах 0,2 ГГц, 0,7 ГГц, 1,2
ГГц путем проведения локальных разрезов зондирующим полуволновым диполем,
последовательно ориентированным в трех ортогональных направлениях. Анализ
экспериментальных данных показали наличие сложного интерференционного поля во
всех обследованных помещениях зданий. Пространственные флуктуации уровня
напряженности ортогональных составляющих поля внутри помещения и среднем
укладывались в пределы 20 дБ, причем положения экстремумов этих составляющих не
совпадали, и результирующий вектор поляризации радиосигнала в любой точке
пространства оказывался ориентированным случайным
образом. Почти во всех комнатах отмечена сильная деполяризация
радиосигналов, так что средний уровень основной составляющей ненамного
превосходил кроссполяризованные.
На рис.1.39 приведены интегральные распределения коэффициента деполяризации определяемого
отношением кроссполяризованных составляющих к основной, полученные в помещениях
административного здания из сборного железобетона при излучении поля
вертикальной поляризации для трех рабочих частот. Как видно из рис.1.39 для
всех частот интегральные распределения близки к логарифмически-нормальному
закону, в то же время наблюдается обратная зависимость коэффициента деполяризации от частоты.
Рис.1.39. Интегральные распределения коэффициента деполяризации
Идентичные результаты получены и в кирпичном административном
здании. Отметим, что такого рода структура поля наблюдалась во всех внутренних
помещениях здания независимо от их размеров, степени заполнения оборудованием и
удаленности от излучателя. Анализ измерений средних уровней напряженности поля
в комнатах по всему зданию показал их сильную и неоднозначную зависимость от
расстояния до излучателя и расположения по этажам. На рис.1.40 представлены
данные измерений среднего ослабления напряженности поля в комнатах,
расположенных на одном этаже с излучателем и на различных расстояниях до него.
Рис.1.40. Среднее ослабление напряженности поля в комнатах
Из рис.1.40 следует, что уровень поля в здании быстро и нелинейно
убывает с расстоянием. Разброс значений ослабления радиосигнала в комнатах,
равноотстоящих от излучателя,
может достигать 10 дБ и более. Также, отмечается частотная зависимость
ослабления поля в структуре здания. Измерения на других этажах по отношению к
излучателю показали скачкообразное уменьшение напряженности поля. Величина
скачка возрастает с увеличением номера этажа и заметно проявляется (до 15…20
дБ) при малых и средних расстояниях. При больших расстояниях величина
поэтажного скачка ослабления существенно уменьшается (6…8 дБ) до почти полного
исчезновения в зависимости от частоты радиосигнала. Вместе с тем
экспоненциальный характер убывания уровня поля с расстоянием сохраняется на
любом этаже. Такая закономерность может свидетельствовать о разных путях и
механизмах прохождения радиоволн в структурных образованиях зданий и по разному
проявляющихся с изменением расстояния между корреспондентами. На рис.1.41
представлены усредненные экспериментальные зависимости ослабления радиосигнала
частотой 1,2 ГГц по этажам здания из сборного железобетона при установке
излучателя на 2-м этаже.
Рис. 1.41. Зависимость ослабления
сигнала по этажам здания на частоте
1,2 ГГц
Подобные зависимости
построены и для других частот 0,2 ГГц и 0,7 ГГц. Многолучевый характер
прохождения радиоволн в структурах зданий весьма усложняет моделирование
процессов ослабления поля, в связи с чем была рассмотрена эмпирическая модель
убывания напряженности поля на основе аппроксимации экспериментальных данных
логарифмической функцией [73]
В этом выражении первый и
второй члены учитывают зависимость ослабления поля от длины радиоволны λ и относительного расстояния r между корреспондентами, λ0 - единичная длина волны. Третий член определяет
между этажный скачок изменения напряженности поля. В этой формуле np - разность номеров этажей
расположения приемника и передатчика. Это выражение применимо в частотной
области 0,1 ГГц £ f
£ 10 ГГц для зданий административного
типа. В других типах зданий жилых, промышленных цехах, складских помещениях и
т.д. закономерности изменений напряжённости поля могут существенно отличаться.
Расчетные кривые для случая n = 0
нанесены пунктирными линиями на экспериментальные результаты, приведенные на
рис.1.40, и показывают удовлетворительное согласие. На рис.1.41 также пунктиром
приведены расчетные зависимости ослабления поля на частоте 1,41 ГГц с учетом
междуэтажных скачков напряженности при np = 0, 1, 2, 3. Эти зависимости
удовлетворительно согласуются с экспериментальными кривыми.
Для сравнения достижимых
уровней дифракционной и волноводной составляющих волн со «сквозной» компонентой
в работе [74] приведены результаты измерений относительного ослабления сигналов
при их распространении из одной комнаты многоэтажного панельного здания в
другую. Схема расположения и нумерация комнат этажа, в которых располагались
приемные антенны, приведена на рис.1.42. В первом цикле измерений передающая
антенна (АНТ на рис.1.42) размещалась в середине крайней комнаты на том же
этаже, что и приемные антенны, во втором цикле - в такой же комнате, но на один
этаж ниже. Измерения проведены на 3 частотах, выбор которых был продиктован
приближением к действующим стандартам мобильных систем радиосвязи: f = 150, 460, 930 МГц. В качестве
передающих и приемных антенн использовались полуволновые симметричные (для f = 460, 960 МГц) или несимметричные
(для f = 150 МГц) вертикальные
вибраторы, поднятые над полом на высоту 1...1,5 м, приемником служил
микровольтметр SMV-8,5.
Результаты всех
проведенных измерений напряжения на выходе приемника нормировались на некоторое
опорное напряжение, в качестве которого было выбрано максимально достижимое
напряжение при размещении приемной антенны в соседней с передатчиком комнате
(на рис.1.42 - комната №1).
Рис. 1.42. Схема
расположения комнат
Как показали дополнительные
эксперименты, приведенные с использованием приемных рупорных антенн, которые за
счет изменения ориентации оси рупоров позволяли определять главное направление
приема, основной вклад в опорное напряжение вносила «сквозная» (через стену)
составляющая поля.
Обработанные результаты
измерений представлены на рис.1.43 и 1.44 в виде графиков нормированного
напряжения Uотн (дБ) в зависимости от условного
номера комнаты, в которой проводились замеры, при размещении передающей и
приемной антенн на одном и том же этаже (рис.1.43 и на соседних этажах
(рис.1.44) здания. Случай относится к измерениям
в комнатах № 1…6, расположенных по туже сторону коридора, что и комната, в
которой находилась передающая антенна, случай б – к измерениям в комнатах №
7…11, расположенных на противоположной стороне коридора. Цифрами на рисунках
обозначены графики, относящиеся к максимальным (1) и минимальным (2) значениям Uотн, зафиксированным в разных местах
одной и той же комнаты.
Анализ полученных результатов показал
следующее:
а) при расположении передатчика
и приемника на одном этаже. Максимальные значения напряжения в комнатах № 2 и
3, расположенных недалеко от передатчика, достигались у ближайшей к передатчику
стенки при наименьшем геометрическом расстоянии между передающей и приемной
антенной, что являлось доказательством превалирующего вклада «сквозной»
составляющей в формирование поля в этих комнатах.
Рис. 1.43.
Зависимость величины Uотн от номера комнаты
В более отдаленных от передатчика
комнатах максимальные значения наблюдались в различных местах комнат, что,
по-видимому, объясняется доминирующем проявлением дифракционной и волноводной
компонент и их интерференцией.
Примечательно, что
уровень относительного ослабления, характеризуемый величиной Uотн, был примерно одним и тем же для всех
используемых частот. При этом в комнатах, расположенных на противоположной
стороне коридора относительно комнаты размещения передатчика, уровень
относительного ослабления был больше в среднем на 10 дБ, что как выяснилось
впоследствии, связано с подавлением проникающей через окно дифракционной
компоненты на этой стороне здания. Вместе с тем, как следует из уточняющих
экспериментов, проведенных в комнате № 11 с использованием на частоте 930 МГц
приемной рупорной антенны, именно дифракционная составляющая давала основной
вклад в результирующее поле не только вблизи окна, но и даже вблизи двери
комнаты. Таким образом, в данном случае волноводная компонента, зависящая от
уровня возбуждения волны в коридоре, была пренебрежимо мала. Однако она резко
возрастала и в итоге становилась доминирующей при передвижении передающей
антенны с середины комнаты к двери, что вызвано существенным повышением уровня
возбуждения волны в коридоре. При исходной позиции передающей антенны в
середине крайней комнаты (рис.1.43) разница между максимальным и минимальным
уровнями относительного ослабления, измеренных в пределах одной комнаты,
достигала значений 20...35 дБ.
Рис. 1.44.
Зависимость величины Uотн от номера комнаты
На основе обработки
результатов измерений получена следующая аппроксимационная формула для
ориентировочной оценки средних значений Uотн а,б(дБ) = -97N +
, где N = 5...10 -
число промежуточных комнат между комнатой расположения передатчика и комнатой
расположения приемника (случай а) или комнатой на противоположной стороне
коридора (случай б), - дополнительное
ослабление, возникающее вследствие условий приема = 0 дБ, = - 10 дБ;
б) при расположении
передатчика и приемника на соседних этажах. В комнатах, достаточно близко
расположенных к передатчику, максимальные значения напряжения достигались при
размещении приемной антенны у окна, что свидетельствовало о превалирующем
вкладе дифракционной компоненты. В более отдаленных комнатах максимальные
значения наблюдались в различных местах комнат, что, по-видимому, вызвано
конкурирующим проявлением волноводной компоненты и ее интерференцией с
дифракционной. Относительное ослабление сигнала имело выраженную частотную
зависимость, наибольшее ослабление зафиксировано для частоты 930 МГц,
наименьшее - для частоты 460 МГц. При этом в случае б по сравнением со случаем а
расположения приемной антенны уровень ослабления был больше в среднем на 27
дБ. Разница между максимальным и минимальным уровнями относительного
ослабления, измеренных в пределах одной комнаты, в случаe а составляла 30...50 дБ, а в случае б существенно
меньше - 20...35 дБ.
Ожидаемые средние значения Uотн можно оценить по формуле
Uотн
а,б = -8,5N
+ (f), дБ, где N = 5...10 - число
промежуточных комнат между комнатой, расположенной над (или под) комнатой
размещения передатчика на соседнем этаже, и комнатой приемника (случай ) или комнатой на противоположной стороне коридора oт
комнаты приемника (случай б); (150МГц) = -13 дБ, (460МГц) = -5 дБ, (930МГц) =-20 дБ, б(f)
= (f)
= -27 дБ.
Аналогичные результаты приведены также в [75].
Полученные результаты измерений и оценки уровней ослабления
радиоволн, а также результаты анализа влияния «сквозной», дифракционной и
волноводной компонент на формирование результирующего поля в различных комнатах
здания могут быть использованы при корректировке существующих
физико-математических моделей радиоканалов передачи - приема сигналов внутри
здания.
2. Анализ
существующих моделей
распределения уровней поля в городских условиях
2.1. Обзор существующих моделей распределения уровней поля
Существуют разные модели
распространения радиоволн в которых авторами предлагаются различные методы
расчета напряженности поля в условиях города. Одни авторы вводят в расчетные
формулы общий эмпирический коэффициент для всего города без учета его застройки,
высоты подвеса антенн, длины волны передатчика. Другие авторы пытаются ввести
эмпирические коэффициенты в имеющиеся теоретические зависимости. Третьи –
пытаются решить эту задачу с помощью теории дифракции и т.д. Ниже приведен
обзор и краткий анализ методов расчета напряженности поля в городских условиях
[1,2].
2.1.1. Расчет напряженности поля эмпирическими методами
Один из доступных
способов оценки всего комплекса механизмов распространения – использование
эмпирических моделей, которые получены на основе усреднения большого объема
экспериментальных данных. Каждая модель имеет собственный набор эмпирических
коэффициентов, а также свои условия применимости.
Одной из первых работ, в
которой на основе измерений выводятся эмпирические коэффициенты для расчета,
является работа Джонса [3]. Он использовал модель, при которой распространение
радиоволн в городе происходит также, как в свободном пространстве, но в
присутствии некоторого абсорбирующего слоя (см. рис.2.1).
Рис.2.1.
К пояснению модели Джонса
Эффект поглощения Джонс
учитывал введением в формулу экспоненциального множителя с эмпирическими
коэффициентами
, В/м, (2.1)
где Pэ – эффективная мощность излучения (P1·D1), Вт;
P1 – мощность передатчика на входе антенны, Вт;
D1 – коэффициент направленного действия передающей антенны;
l – длина волны передатчика, м;
α, K, х –
эмпирические коэффициенты;
r – расстояние между передатчиком и
приемником, м;
h1, h2 – высоты
подвеса передающей и приемных антенн, м;
H – высота поглощающего
слоя, равная высоте
городской застройки Н3,
м;
S –
расстояние, проходимое волной в
эквивалентном поглощающем слое, м
(2.2)
На практике обычно h1<< r и h2<<r, поэтому
. (2.3)
Постоянная ослабления определялась
эмпирически с помощью формулы
(2.4)
в которую подставлялись результаты
экспериментальных исследований.
По результатам измерений
напряженности электрического поля на длинах волн l =
Другая формула для
расчета напряженности поля в городе была предложена в
, мВ/м, (2.5)
где величины – в кВт, a r – в
км.
Следует отметить, что эта
формула пригодна для расчетов при r ≥ 18 h1h2/l.
В
, мВ/м, (2.6)
где
c и m –
эмпирические коэффициенты, определяемые по результатам измерений напряженности
поля.
и (2.7)
Однако в эти выражения не
входят высоты подвеса передающей h1 и приёмной h2 антенн.
В [6] Трифонов П.М. на
основании обработки экспериментальных данных вывел две следующие эмпирические
формулы при высоте подвеса передающей антенны выше уровня крыш:
для частоты 150 МГц
, мкВ/м (2.8)
и для частоты 300 МГц
, мкВ/м; (2.9)
при высоте подвеса
передающей антенны ниже уровня крыш:
для частоты 150 МГц
, мкВ/м, (2.10)
для частоты 300 МГц
, мкВ/м, (2.11)
где величина Рэ в Вт.
К недостатку этих формул
следует отнести отсутствие зависимости уровня поля от высоты подвеса приёмной
антенны.
В работе В.Ф. Чепуры,
основанной на большом количестве экспериментальных данных, была предложена
следующая эмпирическая формула, приведенная также в [5, 7]
, мВ/м, (2.12)
где коэффициенты В и m такие же
как и в формуле (2.7), определяются путем обработки экспериментальных данных. В
этой формуле и далее значения Рэ
приведены в кВт, а r в км. С целью удовлетворения теореме
взаимности, формула (2.8) была несколько изменена в работе [6].
, мВ/м. (2.13)
В отличие от предыдущих
работ В.Ф. Чепура, учитывая, что значения напряженности электрического поля в
некоторой конкретной точке могут отличаться от средних рассчитанных значений
напряженности поля, ввел в формулу так называемый процент "обеспечения", понимая под ним
процент количества пунктов приема от общего числа предполагаемых пунктов, в
которых значения напряженности поля равны или превосходят значения,
рассчитанные по формуле
, мВ/м, (2.14)
где p – процент "обеспечения".
Данная формула
справедлива при процентах "обеспечения" 10 % ≤ p ≤ 95 %.
Важное место занимает
работа Н.И.Бардина и Н.Д. Дымовича [8] , в которой авторы попытались теоретически
обобщить полученные экспериментальные результаты измерений напряженности поля
в г. Ленинграде.
Авторами рассмотрено несколько частных случаев: прием выше уровня крыш, прием
ниже уровня крыш на радиальных и поперечных улицах, прием внутри зданий.
Случай приема ниже уровня
крыш на радиальных улицах (рис.2.2). Если ширина улицы превосходит размеры
первой зоны Френеля, что может иметь место в дециметровом диапазоне волн на
небольшом расстоянии от передатчика, то напряженность поля можно рассчитать по
интерференционной формуле [6, 7]
, мВ/м, (2.15)
где R – модуль коэффициента отражения от
подстилающей поверхности;
q – угол потери фазы при отражении от
подстилающей поверхности.
Или по формуле Введенского
при r ≥ 18h1h2/l
, мВ/м. (2.16)
При ширине улицы В меньше размеров первой зоны Френеля,
расчет должен производиться по формуле Кирхгофа, но пренебрегая влиянием поля
на соседних улицах, перпендикулярных к той, где находится точка приема, и
положив получается расчетная
формула [8]
, мВ/м, (2.17)
где = – расстояние от уровня
крыш зданий до
приемной антенны, м;
– высота застройки
(уровень крыш зданий), м;
– половина ширины улицы, м.
Рис.2.2.
К расчету поля на радиальной улице
Прием на поперечных
улицах. Эта задача была решена в предположении, что стены зданий являлись
поглощающими. Задача сводилась к определению дифракционного поля
на широкой щели,
прорезанной в плоском экране
(рис.2.3), и поле для горизонтально поляризованных волн определялось по формуле
[8]
, мВ/м, (2.18)
где V – дополнительный множитель,
определяемый из графика на рис.2.4;
– высота подвеса
передающей антенны, отсчитываемая от уровня крыш зданий в месте приема, м.
Формула получена в
предположении, что kz >> 1 и kb >>1, где k = 2π/l.
В работе [9] в случае приема выше уровня крыш
рекомендуется использовать для расчета напряженности поля формулу Введенского
(2.16), в которой высоты подвеса передающих и приемных антенн отсчитываются от
уровня крыш зданий.
Рис.2.3.
К расчету поля на поперечной улице
В рекомендациях
Государственного специализированного проектного института радио и телевидения
(ГСПИ) по расчету напряженности поля УКВ приводится формула Г.З. Рубина [10],
полученная им в результате усреднения кривых МККР для систем связи с подвижными
объектами
, мкВ/м (2.19)
или в дБ относительно 1
мкВ/м
(2.19. а)
В работе [11] Д.Н.
Ликонцев немного модернизировал формулу Рубина, введя эмпирические коэффициенты
С1 (вместо 89 в числителе) и m (в степени расстояния вместо 2,7), и
на основе полученных данных измерений напряжённости электрического поля
телевизионных сигналов в г. Алма-Ата и Душанбе получил новые значения этих
коэффициентов для разных частот.
В [12] приводится
методика расчета напряженности поля, использующая кривые [13] для городской
(рис.2.5) и сельской (рис.2.6) местности. Напряженность поля в точке приема
рассчитывается по формуле
, дБ, (2.20)
где Е(50,50) – напряженность поля, дБ,
соответствующая эффективно излучаемой мощности полуволновым вибратором при Рэф = 1 кВт для 50 % пунктов
и времени приема (значения напряженности поля в дБ берутся относительно 1 мкВ/м);
Bp.н.
– поправка, дБ, учитывающая номинальную мощность Рн передатчика Bp.н. = 10 lg (1000/ Рн,);
Bh2 – поправка на высоту подвеса приемной антенны:
, дБ при использовании рис.2.5;
, дБ при использовании рис.2.6;
Dу – коэффициент усиления передающей антенны, дБ;
ВФ – потери в антенно-фидерном тракте, дБ;
∆h – поправка, дБ, при использовании
кривых МККР для сухой и влажной почвы и снижении высоты подвеса приемной
антенны с
Таблица 2.1
Значения
величин поправки
Характер местности |
Величина
поправки ∆E, дБ, на частотах МГц |
||
30 … 150 |
150 … 250 |
450 … 1000 |
|
городская сельская |
9 9 |
11 7 |
14 14 |
В [12] также приводится
методика расчета напряженности поля на улицах и площадях города с помощью
номограмм (рис.2.7 и 2.8). Номограммы были получены на основе дифракционных
формул работы [8] .
В зависимости от условий
приема пользуются номограммами рис.2.7, позволяющими определить напряженность
поля на поперечных улицах города, а также на площадях, открытой местности и
крышах зданий, или номограммами на рис.2.8, позволяющими определить
напряженность поля на радиальных улицах города.
Номограмма на рис.2.7
состоит из четырех сопряженных между собой семейств, характеристик и графика
для определения поправки на пространственные замирания. В первом квадрате по
оси абсцисс отложены значения, м, z = Н3 -
h2 , где Н3 – высота зданий, которыми застроена улица, м, h2 – высота подвеса приемной антенны над уровнем земли, м.
Параметром семейства
кривых является величина b м, – половина ширины улицы, на
которой размещена приемная антенна. При определении ширины улицы следует
учесть, что формулы выведены для случая сплошной застройки, когда соседние
здания примыкают друг к другу. Если между соседними зданиями имеются просветы,
величину следует увеличить в 1,5…2 раза [8] в зависимости от величины
просветов.
В первом квадранте по оси
ординат отложены значения коэффициента K1′, учитывающего местные условия приема на поперечных улицах
города, дБ
(2.21)
где F – функция ослабления, определенная
экспериментально в [8] .
Во втором квадранте
отложены значения коэффициента K1 учитывающего влияние длины волны l, м, на коэффициент K1′
(2.22)
а также значения коэффициента,
учитывающего расстояние r км между пунктами приема и передачи, дБ
(2.23)
В третьем квадранте
приведено семейство кривых, показывающих зависимость коэффициента K4 от длины волны l и
высоты подвеса приемной антенны h2 над поверхностью земли или крышей здания.
В четвертом квадранте
представлены зависимости между коэффициентом K2, учитывающем характеристики передающей системы, и
высотой подвеса передающей антенны h1, причем параметром семейства кривых является произведение Рэф×D у в кВт,
, (2.24)
где h1эф – высота подвеса передающей антенны
относительно местности пункта приема, м.
Квадрантами первым,
вторым и четвертым пользуются для определения напряженности поля на улицах
города, а вторым, третьим и четвертым - для определения напряженности поля на
открытых площадях и крышах зданий.
На рис.2.7 построен также
график, который позволяет определить поправку ∆Eп в дБ, учитывающую пространственные замирания сигнала в
зависимости от заданного
процента пунктов приема в зоне.
Рис.2.7. Номограммы для расчета напряженности поля на
поперечных улицах, площадях, крышах зданий
Рис.2.8. Номограммы для расчета напряженности поля на
радиальных улицах городов и поселков городского типа
Значение напряженности
поля на поперечной улице
в городе Eсп в дБ определяется по формуле
, (2.25)
а для открытой площади или на крыше
зданий Есо в дБ по формуле
. (2.26)
Номограмма на рис.2.8
состоит из четырех сопряженных между собой семейств характеристик и графика для
определения поправки на пространственные замирания сигнала. Семейства кривых,
построенные в квадрантах первом, втором и четвертом, аналогичны соответствующим
кривым рис.2.7, а коэффициенты соответственно равны:
(2.27)
(2.28)
(2.29)
В третьем квадранте
размещено семейство кривых, отражающих зависимость коэффициента М3 от расстояния r и высоты подвеса приемной антенны h2
(2.30)
Значение напряженности
поля на реальной улице Eср в
дБ, определяется по формуле
(2.31)
В [14] предлагается
методика расчета напряженности в городе, пригодная также для расчета
напряженности поля на сухопутной трассе. Напряженность поля, дБ, в городе на
расстоянии r от передающей станции, превышаемая в
L %
мест приема и в T % времени
определяется по формуле
, дБ, (2.32)
где Е(50,50) – медианное значение
напряженности поля (по 50 %) мест и времени при высоте подвеса приемных антенн h2=10 м, РS.=0 дБ (1
кВт);
F(∆h), F(h2) – поправочные коэффициенты, учитывающие степень неровности местности и высоту
подвеса приемных антенн, дБ;
∆E(L) и ∆Е(T) – отклонения
значений напряженности поля от медианного значения в заданных процентах мест L и времени Т
приема, дБ;
РS – эффективная излучаемая мощность, дБкВт;
F(S) –
дополнительное ослабление в условиях города, дБ.
Предлагаемый в [14] метод
расчета полуэмпирический и базируется на документах МККР и ОИРТ [15…18] и
содержит ряд дополнений и уточнений.
Степень неровности
местности оценивается параметром ∆h, который определяется как разница
высот (отметок) местности, превышаемых на 10 и 90 % на определенном расстоянии.
При помощи параметра ∆h вводится
условная классификация типов местности:
равнинная или водная
поверхность ……………………...…..(0...25) м;
равнинно-холмистая
(среднепересеченная) ………………..(25...75) м;
холмистая (
сильнопересеченная ) ……………….....…. (75 …
150) м;
гористая
…………………………………………………..(150 … 400) м.
На рис.2.9 приведены
зависимость медианного значения напряженности поля от расстояния на
равнинно-холмистой местности E(50,50) ( —
линии для метрового диапазона, - - - для дециметрового диапазона при высоте
подвеса приемной антенны h2 =
Зависимости
напряженности поля от
расстояния, приведенные на
рис.2.9 … 2.11, можно аппроксимировать с помощью формулы, дБ
, (2.33)
где Б0 , Б1, Б2
– коэффициенты, значения которых приведены в табл. 2.2 и 2.3 для метрового и
дециметрового диапазонов.
Таблица 2.2
Значения коэффициентов Б0 , Б1, Б2
для метрового диапазона
h1 эф, м |
r = 2 ... |
||
Б0, дБ |
Б1, дБ |
Б2, дБ |
|
20 |
96,14 |
-35,39 |
-3,64 |
37,5 |
100,81 |
-32,66 |
-5,13 |
50 |
102,44 |
-31,16 |
-5,96 |
75 |
104,99 |
-28,73 |
-7,10 |
100 |
105,73 |
-24,45 |
-9,05 |
150 |
105,62 |
-15,92 |
-12,53 |
200 |
102,60 |
-8,16 |
-15,2 |
250 |
99,29 |
-0,25 |
-17,86 |
300 |
95,97 |
7,66 |
-20,51 |
350 |
94,03 |
11,02 |
-21,34 |
400 |
92,08 |
14,38 |
-22,17 |
450 |
90,96 |
16,65 |
-22,64 |
500 |
89,68 |
19,09 |
-23,17 |
Таблица 2.3
Значения коэффициентов Б0, Б1, Б2
для дециметрового диапазона
h1 эф, м |
r = 2 … |
||
Б0, дБ |
Б1, дБ |
Б2, дБ |
|
20 |
103,53 |
-40,07 |
-3,74 |
37,5 |
103,23 |
-40,3 |
-3,57 |
50 |
105,02 |
-36,07 |
-5,41 |
75 |
104,97 |
-28,83 |
-8,83 |
100 |
104,78 |
-21,48 |
-11,93 |
150 |
112,10 |
-22,52 |
-12,46 |
200 |
108,81 |
-13,51 |
-15,67 |
250 |
105,51 |
-4,49 |
-18,88 |
300 |
102,21 |
4,52 |
-22,09 |
350 |
94,89 |
15,80 |
-25,67 |
400 |
88,57 |
25,39 |
-28,52 |
450 |
81,96 |
35,45 |
-31,53 |
500 |
75,92 |
44,76 |
-34,33 |
Эффективная высота
подвеса передающей антенны h1эф, м, для равнинной и равнинно-холмистой (городской)
местности определяется как высота электрического центра антенны над усредненным
уровнем участка земной поверхности 3…15 км в направлении от передающей антенны
к точкам приема [13]
(2.34)
где – высота подвеса антенны над уровнем моря,
Zср – средняя отметка участка 3…15 км, определяемая как среднее
арифметическое отметок всех впадин и возвышенностей.
Рис.2.9. Зависимости медианного
значения напряженности поля от расстояния на равнинно-холмистой местности
Рис.2.10. Зависимости медианного
значения напряженности поля от расстояния на равнинно-холмистой местности для I – III ТВ
диапазонов
Рис.2.11 Зависимости медианного
значения напряженности поля от расстояния на равнинно-холмистой местности для IV–V ТВ диапазонов
Рис.2.12.К
определению значения коэффициента S
Рис.2.13.
К определению значения коэффициента F(S)
Зависимости E=f(r), приведенные на рисунках 2.9…2.11,
можно также аппроксимировать при помощи формул:
(2.35)
если , дБ, (2.36)
где (2.37)
Y, Х –
коэффициенты, которые для частот 30…250 МГц равны Х = 23, Y = 0,15; а в диапазоне частот 450…1000 МГц Х = 27, Y = 0,165.
Значения поправочного
коэффициента F(∆h) можно
определить по следующей формуле
(2.38)
Значения коэффициентов и определяются из табл.
2.4.
Таблица 2.4
Значения
коэффициентов и
Номер телевиз. диапазона параметр |
1 |
2 |
3 |
4 |
5 |
|
18,6 |
20,6 |
25,7 |
38,5 |
45 |
|
0,15 |
0,15 |
0,15 |
0,25 |
0,25 |
Значение поправочного коэффициента F(h2) для равнинно-холмистой местности определяется по формуле
, дБ, (2.39)
где c = 2,6 и 6 для метрового и дециметрового диапазонов волн соответственно.
В городе
при высоте подвеса приемных антенн
Плотность застройки города (или квартала) определяют по плану как отношение застроенной части к общей площади (рис.2.12)
Из рис.2.13 видно, что дополнительное ослабление F(S) для дециметрового диапазона волн [19] изменяется в больших пределах. В частности, для большого города его значение составляет 10…15 дБ. В диапазоне метровых волн оно получается меньше, примерно на 5 дБ. Дополнительное ослабление уменьшается по мере удаления точек приема от передающей станции (рис.2.14) для дециметрового диапазона волн.
Отмечается,
что для крупных городов уменьшение высоты подвеса приемных антенн с 10 до
Поправочный
коэффициент, дБ, для медиан напряженности поля при h2>
(2.40)
где w = 40…50 в дециметровом диапазоне волн.
Окамурой и Омори в [20] в
результате обработки результатов многочисленных экспериментов была предложена
методика расчета напряженности поля Е
с учетом многих факторов, влияющих на распространение радиоволн.
Ими предложено следующее
выражение
E = E0 Am(r,f)H1(h1,r)H2(h2,f)К1К2К3К4, (2.41)
или в дБ (относительно 1 мкВ/м)
E = 20 lgE0 – Am(r,f)+H1(h1,r)+H2(h2,f)+К1+К2+К3+К4,
мкВ/м, (2.42)
где E0 – напряженность электрического поля в свободном пространстве, мкВ/м;
Am(r, f) – ослабление
в городе относительно свободного пространства, дБ, определяется из графика на
рис.2.15;
H1(h1, r) – “высотный
коэффициент усиления” передающей антенны, дБ, определяется из графика на
рис.2.16;
H2(h2, f) – “высотный
коэффициент усиления” приемной антенны, дБ, определяется из графика на
рис.2.17;
K1 – поправочный коэффициент для пригородных и открытых
площадей, дБ, определяемый из графика на рис.2.18;
K2 – поправочный коэффициент для холмистой местности, дБ,
определяемый из графика на рис.2.19;
K3 – поправочный коэффициент для гористой местности, дБ,
определяемый из графика на рис.2.20;
K4 – поправочный коэффициент для преград типа "суша-море", дБ,
определяемый из графика на рис.2.21.
Реудинк в [21] предложил
рассчитывать мощность на входе приемного устройства Pc по формуле, аналогичной формуле
Окамуры и Омори
(2.43)
или в дБ ( относительно 1 мВт)
, (2.43.a)
где P0 – значение принимаемой мощности в
свободном пространстве, мВт.
Все остальные
коэффициенты остались без изменения.
Акеямой, Такаутсу и
Ебиной в своей работе [22] предложена следующая расчетная формула
, дБ, (2.44)
где Emu – величина
стандартного среднего ослабления, дБ;
РETP – значение поправки на эффективную
излучающую мощность, дБ;
Htu – значение поправки на эффективную
высоту подвеса антенны, дБ;
K – значение поправки на окружающую
обстановку, дБ
; (2.45)
S – поправка, учитывающая процент
(α) покрытия области зданиями, дБ;
Кh – поправка на холмы, дБ;
Khf – общая поправка на холмистость, дБ;
Kiм – поправка на отдельные горы, дБ;
КSP – поправка на наклон земной
поверхности, дБ;
КS – поправка на преграды типа
"суша-море", дБ;