Читать онлайн Моделирование канала коротковолновой радиосвязи бесплатно

Информация о научной работе

Реферат: В работе предложена методика расчета отношения сигнал/шум для заданных трассы КВ радиосвязи, времени года и времени суток, с использованием оперативного прогноза станций ИЗМИРАН для определения напряженности поля сигнала в заданной точке приема.

Поскольку напряженность электрического поля помехи в точке приема в соответствии с МСЭ-R P.372-9 определяется для условной антенны симметричный вибратор, находящейся в свободном пространстве, а на практике имеем дело в основном со штатными для приемных центров антеннами ВГД, выведены формулы для пересчета мощности помехи, рассчитываемой в соответствии с МСЭ-R P.372-9 для условной антенны, в мощность помехи, принимаемой реальной антенной ВГД.

Рассмотрены различные способы моделирования помехи с огибающими, распределенными по нормальному, логнормальному и полунормальному законам распределения, проанализирована возможность и адекватность их применения, даны рекомендации по их применению.

Разработана программа KVkanal для Matlab, приведенная в приложении А, в которой вычисляется отношение сигнал/шум для каждой из действующих помех и для суммарной помехи, а также мощность помехи, которую требуется создавать при моделировании КВ радиоканала. Инструкция по работе с программой KVkanal приведена в приложении Б.

Депонировано №122/В217 от 05.10.2017 г. в Федеральном государственном бюджетном учреждении науки Всероссийский институт научной и технической информации Российской академии наук (ВИНИТИ РАН), г. Москва.

Ключевые слова: радиосвязь, коротковолновая радиосвязь, канал связи, напряженность поля, диаграмма направленности, симметричный вибратор, закон распределения, мощность помех, мощность сигнала.

Abstract: In work the design procedure of the relation a signal/noise for set lines КВ of a radio communication, a season and time of days, with use of the operative forecast of stations IZMIRAN for definition of intensity of a field of a signal in the set point of reception is offered.

As intensity of electric field of a hindrance in a reception point, according to

МСЭ-R P.372-9, is defined for the conditional aerial the symmetric vibrator which are in free space, and in practice we deal basically with regular aerials VGD for the reception centres, formulas for recalculation of capacity of the hindrance counted according to МСЭ-R P.372-9 for the conditional aerial, in capacity of the hindrance accepted by real aerial VGD are deduced.

Various ways of modelling of a hindrance with bending around, distributed on normal are considered, логнормальному and to seminormal laws of distribution, possibility and adequacy of their application is analysed, recommendations about their application are given.

Program KVkanal for Matlab, resulted in appendix A in which the relation a signal/noise for each of operating hindrances and for a total hindrance is calculated, and also capacity of a hindrance which is required to be created at modelling КВ of a radio channel is developed. The instruction on work with program KVkanal is resulted in the appendix Б.

Key Words: radio communication, short-wave radio communication, communication channel, intensity of a field, the orientation diagramme, the symmetric vibrator, the distribution law, capacity of hindrances, capacity of a signal.

Введение

Для прогнозирования качества работы канала КВ радиосвязи или для оценки качества работы вновь разработанных средств радиосвязи без проведения натурных испытаний требуется моделировать такой канал радиосвязи с помощью ПЭВМ. В настоящее время такое моделирование можно осуществлять с использованием данных сети ионосферных станций ИЗМИРАН.

Пример оперативного прогноза станций ИЗМИРАН приведен на рис. 1.

Рисунок 1. Пример оперативного прогноза.

Для прогнозирования качества КВ радиосвязи нужно оценить прогнозируемое отношение с/ш на входе приемника, для чего нужно знать напряженности электрического поля сигнала и помехи в точке приема.

Напряженность электрического поля сигнала определяется по таблицам 1 и 2 из рисунка 1 для выбранной рабочей частоты, которая должна выбираться как можно ближе к максимально применимой частоте (МПЧ) для обеспечения минимального затухания сигнала. При применении направленных антенн понадобиться также угол прихода сигнала для определения коэффициента усиления антенны. Будем считать, что в точке приема используется направленная на корреспондента стандартная для приемных центров России антенна Надененко ВГДШ c длиной плеча

Напряженность электрического поля шума в точке приема определяется в соответствии с рекомендациями МСЭ-R P.372-9 [3]. Основными видами шума в КВ диапазоне являются промышленные помехи, атмосферные помехи за счет удаленных грозовых разрядов и галактический шум. Тепловые шумы в этом диапазоне вносят незначительный вклад и их можно не учитывать.

По предлагаемой ниже методике рассчитывается отношение сигнал/шум на входе приемника в заданной точке приема для заданного времени года и времени суток. На основе полученного отношения сигнал/шум производится моделирование канала КВ радиосвязи для заданной трассы и заданного времени работы.

Расчет отношения сигнал/шум на входе приемника

Отношение с/ш обычно обозначается , и вычисляется по формуле:

где Р_с – мощность сигнала на входе приемника;

Р_п – мощность шума (помехи) на входе приемника.

Зная напряженность поля сигнала в точке приема Е_с и параметры приемной антенны, по формулам, приведенным в [2], можно определить напряжение сигнала U_с на входе приемника:

(1)

где F(φ) – функция направленности антенны в вертикальной плоскости;

R_a – волновое сопротивление антенны;

R_f – волновое сопротивление фидера;

– действующая длина антенны,

где F_max – максимальное значение функции направленности.

Соответственно мощность сигнала будет определяться по формуле:

Напряжение помехи рассчитывается по формулам, приведенным в [3] для симметричного полуволнового вибратора в свободном пространстве, то есть, без учета влияния земли. Поскольку в качестве приемной антенны у нас используется симметричный вибратор с характеристиками от полуволнового до волнового и высотой подвеса h=λ/2, диаграмма направленности которого отличается от диаграммы направленности полуволнового вибратора в свободном пространстве, нужно определить коэффициент для пересчета напряжения шума. Для этого сравним диаграммы направленности полуволнового и волнового симметричных вибраторов.

Функция направленности симметричного вибратора, находящегося в свободном пространстве в плоскости, проходящей через ось вибратора записывается в виде [2]:

(2)

где k=2π/λ – волновое число;

– длина плеча вибратора в частях длины волны λ;

α – угол в радианах, отсчитываемый от оси вибратора.

Поскольку для полуволнового симметричного вибратора =λ/4, то для него функция направленности запишется:

(3)

Максимальное значение F(α)_λ/2max=1 при α=π/2.

Функцию направленности в плоскости, перпендикулярной оси вибратора можно записать:

(4)

Для волнового симметричного вибратора =λ/2 и функция направленности будет иметь вид:

(5)

Максимальное значение F(α)_λmax=2 при α=π/2.

Функцию направленности в плоскости, перпендикулярной оси вибратора можно записать:

(6)

Влияние земли на диаграмму направленности антенны учитывают с помощью формулы [2]:

(7)

где k_o – коэффициент отражения от земли, примем k_o=1;

h=λ/2 – высота подвеса антенны;

β – сдвиг фаз между антенной и ее зеркальным отражением, для горизонтальных антенн β=180^о;

φ₁ – угол, отсчитываемый от вертикали.

Тогда

(8)

После перехода к дополнительному углу φ=90^о-φ₁, отсчитываемому от поверхности земли, будем иметь:

(9)

Тогда функции направленности полуволнового и волнового вибраторов в вертикальной плоскости можно записать:

(10)

      (11)

Максимальные значения этих функций будут F(φ)_λ/2max=2 и F(φ)_λmax=4, при двух значениях углов φ=30^о и φ=150^о.

Максимальные значения функций направленности в плоскости вибратора и в плоскости перпендикулярной оси вибратора должны быть равны. Если максимальное значение функции направленности в плоскости перпендикулярной оси вибратора увеличилось в

раз, то и в плоскости проходящей через ось вибратора и расположенной под углом к горизонту, соответствующему максимальному значению функции направленности в плоскости перпендикулярной оси вибратора, максимальное значение

увеличилось в

раз. Поэтому, функции направленности в плоскости проходящей через ось вибратора и расположенной под углом

φ

=30

^о

к горизонту, то есть плоскости, проходящей через середину одного из двух лепестков диаграммы направленности, нужно пересчитать по формулам:

(12)

                                          (13)

В дальнейшем приведенные выше функции направленности (12) и (13) будем считать функциями диаграммы направленности в горизонтальной плоскости.

Мощность помехи, приходящей с некоторого направления под углами α и φ будет определяться по формуле:

(14)

где U_п – напряжение помехи на входе приемника;

– действующая длина антенны;

R_a – волновое сопротивление антенны;

R_f – волновое сопротивление фидера;

E_п – напряженность поля помехи в точке приема;

– коэффициент пропорциональности;

– функция направленности антенны;

– нормированная функция направленности антенны;

F_max – максимальное значение функции направленности антенны.

Будем считать, что помеха принимается антенной со всех направлений верхней полусферы с одинаковой интенсивностью, фазы случайны и равновероятны. Тогда мощность принимаемых помех будет суммой элементарных мощностей ΔР_п, то есть интегралом по полусфере:

(15)

где

(16)

где

– нормированная функция направленности в горизонтальной плоскости;

– нормированная функция направленности в вертикальной плоскости.

Мощность помех, принимаемых из верхней полусферы эквивалентной антенной (симметричным полуволновым вибратором в свободном пространстве), для которого F_н(φ)=1 и F_max=1, будет вычисляться по формуле:

(17)

Чтобы перейти от мощности помех в эквивалентной антенне к мощности помех в реальной антенне, введем коэффициент пересчета, определяемый как:

(18)

Тогда мощность помехи в реальной антенне будет вычисляться по формуле:

(19)

Подставив в формулу (18) значения Р_п и Р_пэ из (15) и (17), получим выражение для вычисления коэффициента пересчета:

(20)

Вычислим значения коэффициентов пересчета для двух реальных антенн с высотой подвеса над землей h=λ/2, полуволнового симметричного вибратора и волнового симметричного вибратора.

Расчет напряжения промышленных и галактических помех и отношения сигнал/шум на входе приемника

Медианные значения коэффициентов промышленного и галактического шума приведены на рис. 10 в рекомендациях МСЭ-R P.372-9 [3]. Для получения этих значений также можно воспользоваться следующей формулой [3]:

(21)

где f – рабочая частота в МГц;

c и d – коэффициенты, приведенные в таблице 1.

Таблица 1. Значения параметров

В таблице для примера приведены коэффициенты только для тех категорий окружающей среды, в которых размещаются приемные центры. В таблице также приведены децили коэффициента шума относительно медианного значения и соответствующие среднеквадратические отклонения, рассчитанные по формуле (22) для верхней децили.

Верхняя дециль D_u соответствует случайной величине, не превышаемой в течение 90% времени.

Нижняя дециль D_l соответствует случайной величине, не превышаемой в течение 10% времени.

Среднеквадратическое отклонение σ соответствует случайной величине, не превышаемой в течение 84% времени.

Медианные значения коэффициентов промышленного и галактического шума F_ппm и F_пgm рассчитываем по формуле (21), причем для галактического шума только для частот f≥10МГц.

Среднеквадратическое отклонение коэффициента промышленного шума от медианного значения σ определяется по графику интегральной функции распределения или по формуле:

(22)

или можно взять из таблицы 1 уже рассчитанное, а для галактической помехи, в соответствии с таблицей 1, σ=1,56 (дБ).