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

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

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

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

Добавлен: 24.10.2024

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

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

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

jRmin от центра Земли до луна и расстояние, соответствующее пути луча от точки излучения до точки R m-,п. На печать в этом случае выведены рабочая частота, угол возвышения А, расстоя­ ние скачка D, R mm. Параметры ионосферы Д, R m, ут вводятся числовым массивом FOF, Н М , YMF. Остальные данные вводятся

вследующем порядке: рабочая частота в ячейке с номером Рх;

вячейки Р 3Ps записываются нули, в процессе счета в них записываются данные из числового массива: Р 0 — угол возвы­

шения (в градусах), далее по порядку 0, Е, N, М, Тh. Нх, шаг по углу А, А/(, шаг по Д, f rk, /Д. Здесь Е — точность; N — чис­ ло уравнений; М — номер уравнения, по которому происходит

сравнение (в

данной программе М = 0); T k — условный

верх­

ний предел

интегрирования; /Д — первоначальный шаг

интег­

рирования при вычислении траектории от начальной точки от­

ражения;

ІД — первоначальный шаг интегрирования при вы­

числении

второй половинки

траектории (Нх — 1, Н2 =

0,01,

Е — 0,0001, N = 2, М — 0).

Текст программы приведен

ниже.

BEGIN REAL FC, FC2, RM, YM, L, BETA, TAY, X, X2, LG, FR2, В, CO, COT, COT2, OTR, KOR, ZI, Z2, Z3, Fl, H, E, FR, FIO, T, TK, DELTA, TETA, TOTR, DELTA2, TETA2, FIB, RNOL, MINUS, Wl, W2, W3, U4, U5, FIRD1, DG1, UGOL, FIBRD2, DG2, D2, Dl, UGOLPR, D, POLD, D3, FCK, RMK, YMK, FRFCK, KOSIN, COX, RADIUS, R07, DO, N1, N2, N3, R071, KRIT, N5; INTEGER N, M, I, DS1, DS2; REAL ARRAY

P[1 : 21],

Y(1 : 2],

Z[1 : 2],

FOF[0 :

],

HM[0 : ], YMF[0 : ],

RO[0:6];;; INPUT(P);

MINUS : = 1 0 ( -6

16);

BEGIN PROCEDURE

 

6

 

6

F(T, Y, Z); VALUE T; ARRAY Y, Z; REAL T; BEGIN FC : =

=

P[2'] + P [

] X Y[2]; FC2 : = FC X FC; RM : =

P[3] + РІ7ІХ

X

Y[2]

;

YM : =

P[4] + P[ ] X

Y2;

 

L : =

FC X YM X P[ ];

 

 

 

6

X

 

 

TAY : = FC2 x P[7];

X : = (RM -

BETA : =

FC2

 

P[

];

 

8

 

 

 

 

 

 

 

6

-

T)/YM;

X2 : = X8X X; LC : = ABS(LC); FR2 : = FR X FR;

В : =

(1 — FC2 X (1 — X2)/FR2)

x T x T ;

 

CO : = P[5] X SIN

(FIO

X 3,14/180);

COT : =

CO +

Y[l];

COT2 : =

СОТ X COT;

OTR : = В — COT2; IF

OTR <

MINUS

 

THEN

GO

TO Ml;

KOR : -

SQRT(OTR);

ZI : = LC X 6370/(FR2 X

KOR

X YM);

Z2 : = 6370 X

COT/(KOR x T); Z3 : = COT/SQRT(B); TOTR : = T;

COX : =

COT;

 

FI : = (ARCSIN(Z3)) X 180/3,14;

Z[l] : = ZI;

Z[2] : = Z2; IF T >

RM THEN BEGIN U4 : = UGOL;

GO TO

M7 END;

END; INPUT(FOF, HM, YMF); FOR I : = О STEP

1

UNTIL

 

DO

ROII] : =

HMII] -

YMF[1];

I : = O;

R07 : =

=

ROII];

R ll: I : =

I + 1 ;

IF I < 5 THEN BEGIN IF

R07 <

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<

ROII] THEN GO TO R ll ELSE BEGIN R07 : = ROII]; GO

145


TO RTl END; END; FR : = P[l]; М2 : W3 : - P[18]; U5 : -

=

P[ 17];

 

W2 : =

0;

 

UGOL : =

P[9];

М3 : R12 : FIO : =

=

(ARCSIN((6370 X COS(UGOL x 3,14/180))/RO7)) X 180/3,14;

DO : = (90 -

FIO -

UGOL)

X 111;

DS1 : = 500; DS2 : = DO;

I : =

DS2 и- DS1;

N1 : =

I;

N2 : =

DO/500; N3 : =

N2 — N1;

R071 : =

RO[I]

+

(RO[I

+

1] -

RO[I]) X

N3;

 

 

KRIT : =

= ABS(R071 -

R07); R07 : = R071; IF KRIT >

2 THEN GO

TO

R12;

P[2] : =

 

FOF[I]

+

(FOF[I

+ 1 ]

-FOFU1) X N3; P[3]:

: =

HM[I]

+(HM [I

+ 1 ]

-

HM[I]) X N3;

P[4] : =

YMF[II +

+

(YMF[I

+

1] -

YMFtI])

X

N3;

P[5] : =

P[3] -

P[4J;

N5 : =

=

N1 X 500 + 1000 -

DO;

 

P[ ] : = (FOF[I + 2] -

 

P[2])/N5;

PI7] : = (HM[I

+ 2 ]

-

P[3])/N5; 6 P[

] : =

(YMF[I

+

2] - P[4])/

/N5;

P[15] : =

P[3]; FIO : =

(ARCSIN((6370 X COS (UGOL X

 

 

8

 

 

 

 

 

X

3,14/180))./P[5])) X 180/3,14; T : =

P[5]; TK : = P[15]; Y[l] : =

=

P[10];

YI2] : =

P['ll];

H : =

P[16];

N : =

P[13];

M : -

=

P[14];

E : = P[12];

URKVH(T,

Y, H, F,

N,

M,

E, TK);

Ml : TETA : =

Y[2J; DELTA : =

Y[l]; BEGIN REAL FC1, RM1,

YM1, FC12, LI, BETA!, TAY1, RAD, RAD2, TIG2, LC1, COl, C012, Bl, PKOR, CNAM1, ZIG1, ZIG2, ZIG3, XOTR, TIG,

TIGK,

HAG,

ARRAY ZIG, YIG[1 : 2]; PROCEDURE FI G(TIG,

YIG, ZIG); VALUE TIG; ARRAY YIG, ZIG; REAL TIG; BEGIN

FC1 : =

FC +

P[ ] X YIG12J;

RM1 : = RM + P[7] X YIG12];

 

 

 

 

 

[

 

] X YI G[2];

FC12 : = FC1 X FC1;

LI : =

YM1 : = YM + P6

 

 

BETA1 : = FC12 X P[ ];

TAY1 : =

= FC1

X

YM1

X P[ ];'

X

 

8

 

 

 

TIG;

RAD2 : =

=

FC12

 

P[7];

 

6

 

 

 

8

 

 

 

=

RAD X RAD;

 

TIG2 : = TIG X TIG;

LC1 : =

(LI

X

(1 -

TIG2) -

TAY1 X TIG +BETA1 X TIG2) x

RAD;

LC1 : =

=

ABS(LCl);

 

COl : =

COX -

YIG[1];

C012 : =

COl

X COl;

B1 : =

(1— FC12 X

(1 — TIG2)/FR2) x

RAD2; PKOR : =

B1 —

C012; IF PKOR <

MINUS THEN W1 : =

W1 + 1

ELSE

GO

TO M55; IF W1 <

 

20 THEN BEGIN TOTR

: = TOTOR -

0,05;

GO TO W22;

END

ELSE GO TO M5; M55 : CNAM1 : = FR2 X

x “SQRT(PKOR);

 

ZIG1 : =

6370 x

LC1/CNAM1;

ZIG2 : =

= 6370 X COl X YM1/(RAD x

SQRT(PKOR));

ZIG3 : =

=

COl/SQRT^Bl);

 

FIB : = (ARCSIN(ZIG3))

X 180/3,14;

ZIGtll

: = ZIG1;

 

ZIG[2] : =

ZIG2;

RNOL : = RM1 - YM1; END;

TOTR : = TOTR -

0,05; W1 : = O;

W22 : XOTR : =

(RM -

TOTR)/YM;

 

TIG : =

XOTR;

TIGK : =

1; FIAG : =

P[21];

E : = P[121;

YIGIlj : = РИ0];

YIG[2] : = P[lll; URKVFI(TIG,

YIG, HAG FIG, N, M, E, TIGK); DELTA2 : = YIGll]; TETA2 : =

146


=

YIG12];

END; DG1 : = 90 — FIO -

UGOL; D1 : = 111 x'

X DG1;

 

RADIUS : = RNOL X SIN(FIO X 3,14/180);

IF

RADIUS > 6370 THEN

BEGIN

D : =

 

D1 + (90 -

FIB)

x

X 111 +TE TA +TETA2;

 

OUTPUT('E',

 

 

 

FR, UGOL,

D,

RADIUS, FIB, 7'); (GO TO M5 END; FIBRD2 : 1,57 -

ARCSIN

(RNOL

x

(SIN(FIB

X 3,14/180))/6370);

 

DG2 : = 90 -

FIB -

FIBRD2

x

180/3,14;

 

D2 : = 111 X DG2;

 

 

 

UGOLPR : = 90 -

FIB — DG2;.

D : =

Dl =

TETA + TETA2 + D2;

 

POLD : =

=

D/2; D3 : =

 

POLD -

Dl; FCK : =

P[2] + P[ ] X D3; RMK :

: -

P[3]

+

P[7]

X

D3;

YMK : - P[4]

+

P[

 

]

X D3; FRFCK : =

 

 

 

6

 

 

FRFGK,

- FR/FCK; OUTPUT('E', FR, UGOL,

RMK,

YMK,

 

 

8

 

 

 

 

 

 

 

UGOLPR,

'/')■

OUTPUT ('E',

D, 7');

M5 : UGOL : =

UGOL +

-f- P[ 17]; IF UG O L<P[18] THEN

GO

TO М3; M7 : P[17] : =

= O; W2 : =

W2 + 1;

 

P[18] : =

O;

 

IF (UGOL >

2

THEN

UGOL ; =

U4 — 2 + 0,

2 x

W2

ELSE

UGOL : = W2 x 0,2;

IF

W2 <

10 THEN GO TO М3 ELSE

BEGIN

P[17] : = U5;

P[18] : = W3;

 

GO TO

M4 END; M4 : FR : =

FR

+ P[19];

IF

FR < P[20] THEN GO TO М2; END END END

 

 

 

 

 

4.

Примеры способов распространения радиоволн

 

 

на радиолиниях разной протяженности и направления

 

 

Ниже приведены таблицы,

в которых содержатся результаты

расчета числа скачков для фиксированных рабочих частот и углов излучения А. Время суток, указанное в таблицах, отно­ сится к пункту излучения. Приведенные в таблицах величины

вычислены со следующей точностью: А = ±

0,°5, і — ± 0 ,1 мсек.

В табл. I—VIII приведены расчеты для линий меридиональ­

ного и широтного направления с пунктом

излучения в Москве.

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

лета.

Длины радиолиний приняты 4000,

6000 и 10 000 км.

В

табл. I —IV содержатся результаты

расчета линий, нап­

равленных к югу от Москвы, в табл. V—VIII — линий, направ­ ленных к востоку.

Эти данные могут быть использованы и для других анало­ гичных радиолиний, пункт излучения которых расположен в полосе широт 50—60° С.

В табл. IX содержатся расчеты для линии в 10 000 км ши­ ротного направления восток — запад для зимнего сезона. Расчеты выполнены для полуденных условий в пункте излучения. В этой таблице приведено также абсолютное время распростра­ нения сигнала в миллисекундах.

Сезон года, указанный в таблицах, относится к северному подушариіо,

147


 

 

 

 

 

 

 

Т а б л и ц а

I

 

 

 

 

 

 

 

 

Зима,

W =

0, направление М—»ІО

 

 

 

 

jD =

4000

 

 

 

D

~

6000

км

 

 

D =

10 000 км

/ ,

М г ц Д , гр ад С п о с о б

/ ,

М г ц Д , гр ад С п о с о б

/ , М г ц Д , гр ад

 

С п о с о б

 

 

 

 

 

 

 

 

 

Утро

 

 

 

 

 

 

 

5

10

2F 2

 

5

 

 

0

2 F 2

5

2

2F2 +

2Е

 

5

20

3F 2

 

5

 

10

3 F 2

5

5

2 F 2 + 3£

 

 

 

 

 

 

 

 

 

 

 

5

10

k F 2

+

2 Е

 

 

 

 

 

 

 

 

 

День

 

 

 

 

 

 

 

5

2

2Е

 

5

 

 

2

3Е

 

5

2

Ъ Е

 

 

 

5

6

3Е

 

5

 

 

5

 

5

4

6Е

 

 

 

5

10

4 Е

 

5

 

 

7

Ъ Е

 

5

5

1 Е

 

 

 

10

2

2 Е

 

5

 

10

Ъ Е

 

5

7

8Е

 

 

 

10

8

З Е

 

10

 

 

2

3Е

 

10

2

Ъ Е

 

 

 

15

0

i F 2

 

10

 

 

5

4 Е

 

10

4

Ъ Е

 

 

 

15

12

2 F 2

 

15

 

 

4

2 F 2

10

7

7Е

+

2 Е

 

 

 

 

 

 

 

 

 

 

 

15

3

2 F 1

 

 

 

 

 

 

 

 

 

 

 

15

8

I F 2 + 2 F 1 + 2Е

 

 

 

 

 

 

 

 

 

 

 

15

12

2 F 2 + I F 1 + 1 F 2

 

 

 

 

 

 

 

 

 

 

 

20

4

2 F 2

+

i F l

 

 

 

 

 

 

 

 

 

Ночь

 

 

 

 

 

 

 

5

И

2 F 2

 

5

 

 

4

2 F 2

5

3

2 F 2 + 2 Е

 

5

20

3 F 2

 

5

 

10

3F 2

5

5

2 F 2 + Ъ Е

 

10

0

I F 2

 

10

 

 

6

2 F 2

5

7

3F 2 + 2Е

 

 

 

 

 

 

 

 

 

 

 

5

10

S F 2

+

4 Е

 

 

 

 

 

 

 

 

 

 

 

10

1

2 F 2

+

1Е

 

 

 

 

 

 

 

 

 

 

 

10

6

3 F 2

+

17?

 

 

 

 

 

 

 

Т а б л и ц а

II

 

 

 

 

 

 

 

Зима, W =

100, направление М —>ІО

 

 

 

 

D — 4000

Nj\t

 

 

 

 

D =

6000

к .At

 

 

D = 10 000 к л і

f ,

М г ц Д , гр ад

С п о с о б

/ ,

М г ц Д , г р а д

С п о с о б

/ , М г ц Д , гр а д

С п о с о б

 

 

 

 

 

 

 

 

 

Утро

 

 

 

 

 

 

 

5

9

2 F 2

 

 

5

 

0

 

2 F 2

5

4

 

2 F 2 + 37?

 

5

20

3 F 2

 

 

5

 

9

 

3 F 2

5

9

 

Ъ Е 2

 

10

2

i F 2

-

 

5

 

15

 

A F 2

5

12

 

6 F 2

 

 

 

 

 

 

 

 

 

 

 

 

5

16

 

1 F 2

148