Файл: Автоматизированная система обработки и интерпретации результатов гравиметрических измерений..pdf

ляется одна зона — 1010

кодов.

Вводится пакет 1. Программа ВсПЗ запишется па МЛ. Вво­

дится пакет 2. Последовательно будут введены

и

записаны

на МЛ

все массивы МЗ. Вводится пакет 3. Последовательно будут

введены

и записаны на МЛ все массивы М4.'Вводится

пакет 4.

В зависимости от объема исходных данных в течение

нескольких

минут создается массив

(ж, г/, z, gH}, записывается

на

МЛ и

печа­

тается следующая информация: 1) количество

тетрад х, у, z, gn и

массив этих тетрад; 2) количество пар {ПРПК, g„}, не

имеющих

координат,

п массив

этих

пар; 3) количество тетрад х,

у, z,

ПРПК,

не имеющих соответствующих пар ПРПК.

Каждый

массив

{.г, у, z, gH }, записанный

на

МЛ,

имеет

вид

гг. {х, у, z, gn}, где п — число

тетрад.

данных

для ОсП8.

4°. ВсП8 служит для подготовки исходных

С ее помощью составляется матрица заданных

прямоугольных

раз­

меров из исходных данных, записанных на МЛ в виде

квадратных

матриц 16s X 16s.

Программа использует СП с номерами: 0000; 0042.

Исходные данные для ВсП8 готовятся

в

виде

двух

мас­

сивов.

Массив

M l — информация

о МЛ готовится

так, как

указано

в гл. X V I I , разделе

1.

Массив М2 — информация о счете составляется для ВсП8 и запи­

сывается в

десятичном виде:

ная

матрица

из п строк

и m столбцов, записанная последовательно

по

строкам).

При этом

подразумевается, что m <1000, п < 1000,

номера 0-f-?i—1,

столбцы

0-г-пг—1.

(і =

А л г о р и т м .

В

данном массиве

значений

функции

Frj

= 0 , 1,

. . ., п—1;

/

=

0,1,

. . .,

т—1)

определяются

^ т а х

и

Fmin

и

число

изолиний,

которые

должны

быть проведены в

этом массиве,

по

формуле

Fmax — F min

где

С — сечение

изолинии.

Если

I >

100,

то

происходит

передача

управления

в

ячейку

по А2, а вычисления не производятся.

Затем по программе

определяется величина

максимальной

изо­

линии

по

формуле

/і

=

если

^ т а х ^ 0 ,

[х] — entier (х),

[ С [ % ^ ] >

е с л и і ? т а х < 0 ,

где

т. е.

наименьшее

целое

число,

не

превыша­

ющее

X.

Величина

других

изолинии

определяется

из

соотношения

fk+i =

fk — C-

Если fk

5г

^тіпі то проверяется,

находится

ли данная

изолиния

ордината

изолинии

(d) вычисляется линейной интерполяцией: если

в

строке

матрицы

изолиния

величины

fk

находится

между

Ftj

и

-^ѵ+і (изолиния

пересекает і-ую строку

между /-ым

и /' +

1-ым

столбцами),

то

Если

fk

<Fmin,

происходит

переход

от

строк матрицы

к ее

В и д в ы д а ч и

р е з у л ь т а т о в .

Программа печатает вели­

чину изолиний,

а

затем — ее координаты вместе с

номером той

строки, где эта

изолиния имеет такую

координату.

Печать начи­

и координат изолиний,

находящихся в строках

матрицы,

аналогично

осуществляется печать

по столбцам. Величина

изолинии

печатается

с пятью знаками мантиссы (последние

четыре

цифры — пробелы),

координата'изолинии — также с пятью

знаками

мантиссы, отделен­

строки

(столбца), где эта изолиния имеет

такую

координату.

Например, имеется исходный массив

1,0

2,3

0,9

0,7

3,1

3,5

m =

3, п = 2 и С — 1,0. Тогда будут отпечатаны данные,

приведен­

ные в

табл.

80.

Т а б л и ц а 80

+++

01

3000-

Интерполяция по строкам

-Максимальная

изолиния

++ +

00

95833

001-

-Номер строки

+++

t

-Координата изолинии

01

20000*

-Значение изолинии

+++

{_

-Координата изолинии

00

76923

000*

-Нулевая строка

+++

01

12143

000

Координата изолинп в нулевой строке

+++

00

54167

001

Координата изолинии

в первой

строке

+++

01

10000

000

+++

00

00000

Интерполяция

по столбцам

+++

01

19286

000

+++

00

12500

001

+++

01

30000*

-Максимальная

изолиния

+++

00

87500

ООІн

-Номер столбца-

+++

t _

-Координата. изолинии

00

80769

002

+++

01

20000

002

+++

00

42308

+++

01

10000

000

+++

00

00000

+++

01

38422

002

фровываются соответственно значению

максимальной изолинии.

5. Аналогично

строятся все остальные изолинии (рис. 39).

О б р а щ е н и е

к СП. Программа

оформлена как стандартная

х - 1

ООО

16

X

7501

7610

X

52

т.

№ С П

п

х + 1

52

a

A

ß

Рже.

39. Пример жостроення

карты

по СП-0154.

1

где

п +

я3

[РА] — число строк матрицы (в восьмиричной системе);

m +

п1

[РА] — число

столбцов

матрицы

(в восьмиричной си­

стеме) ;

а +

я х

[РА] — адрес

ячейки,

где находится

шаг изолинии;

куда

ß +

я 3

[РА] — адрес начала матрицы (F00);

А — адрес ячейки,

передается

управление,

если

изолиний

больше 100.

Х а р а к т е р и с т и к а

СП. Длина СП равна

0234 ячейкам;

х - 1

000

16

X

7501

7610

X

00я3

00

ß

0155

Y

х + 1

яіОхСз

00

п

0000

п+1

х + 2

00лз

00

а

0000

со

где ß — ячейка,

в которой хранится дата (год, месяц,

день) в

деся­

тичном виде (І-г-4

разряды мантиссы — год, 5 и 6

разряды

ман­

тиссы — месяц,

7

и 8 разряды мантиссы — день);

у — ячейка,

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

в радианах; (п + 1) — номер ячейки,

в которой записана в двоич­

ном

виде

географическая

долгота точки

наблюдения,

выраженная

в радианах; а — помер ячейки, в которой

записан в двоичном виде

номер часового пояса;

со — ячейка, в которую

помещается резуль­

тат,

т. е.

вычисленная суммарная поправка за влияние лунно-

солнечного

притяжения.

Длпна

СП равна

0142

ячейкам.

Рабочие

ячейки

0001—0003.

V—1

000

16

у.

7501

7610

х

Лі0я3

. 00

а

0156 -

сс+1

и + 1

00л 3

00

0000

00U0

ш

находится дробная

часть долей суток в

восьмеричном виде:

со — ячейка, начиная

с которой помещается

результат в двоичном

х - 1

000

16

X

7501

7610

у.

ят0л3

00

ß

0157

со

где

ß — ячейка,

в которой

записывается

юлианская дата,

выражен­

ная

в днях; ю — ячейка,

начиная

с которой

помещается

результат