Файл: Керблай, Т. С. О траекториях коротких радиоволн в ионосфере.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 24.10.2024

Просмотров: 62

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Т а б л и ц а 19

СПараметры

(/о^2)ор

РтЩср

( У т Р 2 ) с р

(/о^2)Ср

+ 6

У т + 3

Ѵ о Щ с р

Л т - 3

У т G

Д,

гр а д

? F 2

4 F 2

k 1 со

9

4

9 - 1 0

2

7

 

А , гр а д

Параметры

 

3F 2

4 F 2

ioF2 +

а

 

(^т)ср

3

8

 

 

( У т ) С р

 

 

faF2 — б

 

(^т)ср

/рас >

ѴПЧ

0 / т ) ср

 

 

о межслоевой ионизации. Можно предполагать, что оно значи­ тельно больше вследствие изменении способа распространения.

До сих пор рассматривались эффекты изменчивости ионо­ сферы на сравнительно коротких радиолиниях, в основном одиоскачковых.

В табл. 19 приведены некоторые сведения о влиянии флук­ туаций ионосферы па углы прихода радиоволны фиксированной частоты для радиолинии •— 9000 км, пролегающей в северном полушарии в полосе шпрот 20—60°. Для расчета выбран наи­ более простой случай (зима, ночь, W ~ 80), когда ионосферу можно считать однослойной. При расчете предполагалось, что флуктуации ионосферы имеют один знак па протяжении всей радиолинии, что в общем случае не всегда может выполняться. Из таблицы следует, что иа протяженных радиолиниях эффект флуктуаций ионосферы, выражающийся в изменении А, не пре­ вышает аналогичный эффект на одиоскачковых трассах. Однако следует отметить, что изменение способа распространения, о чем говорилось выше в применении к коротким радиолиниям, на длинных трассах играет еще большую роль. При этом присходит и изменение углов А. Некоторое представлепне о разнооб­ разии способов распространения в дневное время па длинных трассах можно получить из материалов гл. VI и прил. 4.


ПРИЛОЖЕНП Я

1. Карты распределения горизонтальных градиентов критической частоты слоя F 2

На рис. 75 и 76 даны карты, характеризующие распределение медианных горизонтальных градиентов критической частоты слоя

F2, связанных с изменением местного времени. На рис. 77

и 78

приведены карты

градиентов f 0F2, построенные

по данным из­

мерений f aF2 за

один

деиь.

местного

вре­

По оси абсцисс на

рис. 75—78 нанесены часы

мени, по оси ордииат — геомагнитные широты. Сплошные линии соответствуют положительному направлению градиента при дви­ жении с запада на восток, пунктирные — отрицательному.

На рис. 79—82 приведены карты, характеризующие распре­ деление медианных горизонтальных градиептов f 0F2 в зависи­ мости от широты для часов восхода и захода. На картах нанесена географическая сетка: по оси абсцисс отмечены долготы, по оси ордииат — широты. Сплошные линии соответствуют положитель­ ному направлению градиента при движении с севера па юг, пунк­ тирные — отрицательному.

Величины градиентов даны в мегагерцах и отнесены к 100 км.

Рпс. 75. Карта градиентов Январь, 1964 г.

f0F2,

связанных с изменением местного времен]

Ш и р о т а

S

Ш и р о т а

Рис. 76. Карта градиентов f0F2, связанных с изменением местного времени. Июнь 1964 г.


Рис. 7 7 .’ Карта градиентов f0F2, связанных с изменением местного времени 5 января 1958 г.

Рис. 78. Карта градиентов f0F2, связанных с изменением местного времени

28 января 195S г.

Рис. 79. Карты градиентов /0/’2, связанных с изменением широты. Январь 1958 г., период восхода

Рис. 80. Карта градиентов Гf„F2, связанных с изменением широты. Июнь 1958 г.,'период восхода

Рис. 81. Карта градиентов /0F2, связанных с изменением шпроты. Январь 1958 г., период захода

Рис. 82. Карта градиентов f0F2, связанных о изменением широты. Июнь 1958 г., период захода

2.Программа расчета'траекторий луча

визотропной трехмернонеоднородной параболической ионосфере

Программа решения системы обыкновенных дифференциаль­ ных уравнений (4.2) методом Рунге — Кутта четвертого порядка с автоматическим выбором шага интегрирования позволяет рас­ считать траекторию луча, расстояние, перекрываемое лучом, углы прихода в вертикальной и горизонтальной плоскостях,

время распространения луча. Программа

составлена

на

языке

«Алгол-60» для

a-транслятора (машина БЭСМ- ) н

реализована

в ВЦ

АН СССР.

 

 

лучей, вышедших

Расчет может быть произведен для пучка

6

 

— угол

из точки излучения под разными углами Д

и а

 

возвышения, а

— азимутальный угол между направлением луча

и заданной линией). Программа составлена

так, что сначала для

01

 

01

 

01

 

одной

рабочей01частоты происходит перебор

углов

Д01, затем-

углов а01. Шаг

по углам Д п а0і может быть

произвольным.

Программа

состоит из

двух блоков: внешнего

и

внутрен­

него.

Во внешнем блоке

рассчитывается

траектория

луча от

01

 

 

 

 

 

 

точки излучения до точки отражения. Во внутреннем блоке происходит расчет характеристик траектории от точки отражения до точки выхода луча из слоя. Внешний блок содержит процедуру F счета правых частей уравнений (4. 2):

tgqp =

(соі +

б і)

V 1 +

t g 2 l|) COS2 ф

 

n-R- -

(coi + öl)2 (1 -b tg-opcos Ф)

tg ^

V

dS =

< tö /l

+

4

tg*q>, d2ö ^ ^dQd S ,

tg? >+

tg ф coi • 6 i

dn

d62= ^ - d S ,

dX

dB = tg ф

d%=tgty^',

dt = dSßOOii.

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

,

.

І к п 4 - І * ° . ѵ

о

_ Д

,

дНт д .

дНт ..

/ с

/с о \

QQ ^ " Г

%■>

“ m

-“

mo п

^ “Г

д%

ди

ди

У т = У т 0 + T

e + ^ f Ь

гд е fco, Rmo> Уто — параметры

ионосферы (/с — критическая ча­

стота; R m — высота максимума,

отсчитываемая от центра Земли;

ут — полутолщина); dfJdQ, d R J 50, dyJdQ — градиенты па­

раметров в направлении заданной линии (параметры ионосферы и их градиенты в плоскости линии привязаны к начальной точке траектории в ионосфере); dfjd%, dR J d %, дут/д% — градиенты пара­ метров в поперечном направлении. Расстояние Ѳотсчитывается от точки вхождения луча в слой. Показатель преломления п и его производные вычисляются по формулам § 1 гл. IV.

138



Перед обращением к стандартной процедуре решения системы дифференциальных уравнений URKVH происходит расчет ха­ рактеристик траектории луча от точки излучения до точки вхо­ ждения в слой. По специальной подпрограмме определяется высота нижней границы слоя і?01 (%0> В0) Для луча, направлен­ ного под углом а к заданной линии. Затем рассчитываются:

Яоі

(Хоі) =

Rmo (Xoi) — Ут (x)>

Xoi = Xoi (B i);

 

— проекция

пути

луча

на

горизонтальную

плоскость;

Хоі — отклонение на­

чальной точки

траектории луча в

ионосфере от плоскости за­

 

0

0O1

 

данной линии; угол ср01, характеризующий угол между каса­ тельной к траектории и радиусом R, проведенным из центра Земли в точку вхождения луча в слой (см. рис. 43).

Нижний предел интегрирования равен величине R 01 (Хоі)- Для верхнего предела интегрирования условно выбрана вели­ чина R > і ? т , поскольку в процедуре F интегрирование пре­ кращается или при выполнении условия отражения, или при выполнении условия R > R m (х,Ѳ).

Внутренний блок включает процедуру FIG счета правых частей уравнений, записанных несколько по-иному и обращение

кпроцедуре URKVH. Иная запись уравнений вызвана переходом

кновой переменной интегрирования х = (R m R)/ym. В ре­

зультате нижний

предел

интегрирования

становится

равным

хотр = (Rm R0Tp)lym, где

йотр — высота

точки

отражения.

Верхний предел равен единице.

 

F.

Здесь

Процедура FIG построена аналогично процедуре

также происходит

переход

от скользящей

системы

координат

к системе координат, привязанной к заданной линии, и вычи­

сление параметров

ионосферы в текущей точке траектории по

тем же формулам,

что и на первой половине траектории.

В результате

расчета второй части траектории становятся

известными характеристики траектории в точке выхода луча из слоя: угол ср02; угол ф02; расстояние В, соответстующее пути от точки отражения до точки выхода луча из слоя; %— отклонение точки выхода от дуги, соединяющей точки излучения и приема;

R Q (Ѳ,х) — высота нижней границы слоя

в

точке выхода.

Далее происходит расчет расстояния Ѳ

по известному

углу

ср02и высоте і?02, углов Д2, а2и

(%і2характеризует отклонение

2

линии).02В

 

 

точки прихода луча от заданной

том случае,

когда

%02

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

На печать выдаются следующие параметры: рабочая частота (в мегагерцах); угол возвышения в точке излучения (в градусах); угол прихода А в вертикальной плоскости (в градусах); угол а, характеризующий отклонение направления прихода от дуги, соединяющей точки излучения и прихода (в градусах); откло­ нение точки прихода луча от плоскости дуги заданной линии % (в километрах); расстояние, перекрываемое лучом одним скачком (в километрах); время распространения (в мсек).

139