книги из ГПНТБ / Автоматизированная система обработки и интерпретации результатов гравиметрических измерений

Т а б л и ц а 79

№ МЛ

0000

номер

зоны

0000

№ МЛ

0000

номер

зоны

0000

№ МЛ

0000

номер

зоны

0000

№ МЛ

0000

номер

зоны

0000

1-ая

строка — число

массивов МЗ {x, у, z, ПРПК}.

2-ая

строка — число

массивов М4 {ПРПК, gn}.

3-я

строка — контрольная сумма М2.

Массив МЗ * состоит из нескольких числовых массивов, каждый

из

которых содержит не более чем 85

перфокарт (255 тетрад х,

у,

z, ПРПК) и имеет свою контрольную

сумму. Первое число мас­

Составляются четыре пакета к счету:

Пакет

1 — перфокарты: «вызов ИС-2», чистая;

1 — запись ВсПЗ

на |МЛ; чистая; M l ; чистая; 2-^-28 — ВсПЗ.

Пакет

2 — перфокарты: «вызов ИС-2»;

чистая;

29, 30 — запись

МЗ; чистая; M l ; чистая; М2; чистая;

первый

массив

МЗ; чистая;

второй массив МЗ; чистая;

. . .; чистая; последний массив МЗ.

Пакет

3 — перфокарты:

«вызов .ИС-2»;

чистая;

31 — запись

М4; чистая; M l ; чистая; М2; чистая;

первый

массив

М4; чистая;

второй массив М4; чистая;

. . .; чистая;

последний

массив М4.

Пакет

4 — перфокарты:

«вызов

ИС-2»;

чистая;

33 — счет;

чистая; M l ; чистая; М2.

Размечается магнитная лента (или две) на такое

количество зон

по 2000 кодов, какое необходимо для записи

МЗ и М4. Дополни­

тельно нужно разметить столько зон по 2000 кодов

для результа-

ляется одна зона — 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

со

где

ß — ячейка,

в которой

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

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

выражен­

ная

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

начиная

с которой

помещается

результат

тров N„,

которые "должны быть расположены

в МОЗУ подряд,

начиная

с ячейки Ак

и до

Ак+п^.

Я

КОП

Л'і

0000

0000

Л

КОП

0000

оооо

я.

КОН

0000

0000

45—4-3

42—40

39-37

36—25

24—13

12—1

(разряды)

сх

с2

с3

аг

а2

п3

(содержание

Например в ячейках 0127 и 0150 нужно

сформировать

первый,

а в ячейке 0137 — третий

адрес по

параметру 0060, в

ячейке

0575 — второй

и третий

адреса, в ячейке 0666 — первый

адрес по

параметру 0211.

Содержимое

этих

ячеек

следующее:

0127

1

12

0000

0121

0001

0137

6

02

3000

4000

5000

0150

1

12

0000

. 0116

0001

0575

3

06

0077

2000

2000

0666

1

12

0000

0645

0001

0100

4

14

0127

0137

0150

0101

0

00

0000

0000

0000

0102

3

40

0575

0666

0000

0110

0

00

0060

0000

0000

0111

0

00

0211

0000

0000

Если обращение

к

программе имеет вид

х

0

16

х + 1

7501

7610

х + 1 •

0

02

0100

0160

ОНО

будут сформированы

так:

0127

1

12

0060

0121

0001

0137

6

06

3000

4000

5060

0150

1

12

0060

0166

0001

0575

3

06

0077

2211

2211

0666

1

12

0211

0645

0001

О б р а щ е н и е к СП. Программа оформлена как

стандартная

в системе ИС-2. Для обращения к ней используются

две подряд

X

0

16

х + 1

7501

7610

х + 1

я

0160

А

где п — число параметров; А — ячейка

первого

параметра;

А Ш К —

адрес ячейки шкалы первого параметра.

Если в ячейке

х + 1 я = 0, происходит

формирование

по

команде 13, если я =

2, — по

команде 33, т. е. вычитание

N{.

Ни

шкалы,

можно

сформировать

зависимость от любого Nt параметра;

в этом

случае

меняется лишь вторая

команда обращения.

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

СП. Длина СП равна 0040 ячейкам.

Длина

счетной

части — 0021 ячейка.

Количество нестандартных

подразумевается, что длина

строки матрицы N < 1000.

Двоичная

запись материала

в МОЗУ

сохраняется.

О б р а щ е н и е

к СП. Программа оформлена как

стандарт­

X

000

10

х + 1

7501

7610

х + 1

я

72

ан

0161

ак

х + 2

я

72

0000

0000

N

(в

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

системе); N + я 3 [РА] — длина строки матрицы

(в

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

системе).

к

000

16

х + 1

7501

7610

х + 1

000

72

0100

0161

0105

х + 2

000

72

0000

0000

0003

2 = ^4(1 +cos 2яг/5),

где

А == 0,37 км;

г2 = (х — х0)2

+ (у — y0f;

х, у , z — текущие

координаты;

х0, у 0

— координаты

центра

площади.

Исходные

данные выдаются по площади 15,5 X 5,5 км с шагом

0,5 км

(31 X 11 точек).

Одновременно по этой программе вычисляется значение на пло­

скости

относимости zQ — 0,75 км

(21 X 7

точек). Положение

изо­

линий

рельефа

относительно

исходного

поля

показано

на

рис.

31.

В процессе счета используются СП с номерами: 0010; 0016; 0042;

0071;

0154; 0161.

Для

выдачи результата программу ТП1 вводят в машину в виде

пакета

перфокарт: «вызов ИС-2»; чистая; 1 — 12. Проводятся вычис­

ления

в течение 1 мин и печатаются: 1)

значения

и координаты

изолиний {Vz (s0 )};

2)

матрица {Vz (z0 )};

3)

значения

и

координаты

изолиний

{Ѵг

(z)};

4)

матрица

{Vz (z)};

5)

значения

и

координаты

изолиний {z (ж, у)};

6) матрица {z (ж, у)};

затем перфорируется массив

7і {x,

у, z, У2 (z)};

ТП1 приведен в табл. 81.

Пример

выдачи

Т а б л и ц а 81

1

1

1

02

100

00

Зпачеипя

п координаты,

нзолшшіі

Мат­

~!

1 1

01

15$

27

001

рицы ( 7 z ( z 0 ) }

+++

01

239

79

001

+++

01

497

013

400

Матрица

{Vz{zQ)}

+++

01

727

143

001

__1

1

01

'

470

502

030

Последнее

лиело в матрице

{z (.г,

у)}

2°. Вторая и третья тестовые

программы (ТП2

и ТПЗ)

предна­

значены для выдачи исходных тестовых данных для

ОсП7 и

ОсП5.

В качестве

общего числового

теста используется гравитацион­

ное поле Ѵг на

квадратном участке размерами 24

X 24 км,

созда­

(20;.3,5; 2 км),

(9; 12,5; 2,5 км), (14; 16,7;

3,5 км), (-15; - 1 0 ;

60

км),

с плотностями

а и радиусами R(a; R)

соответственно: (0,544 г/см3 ;

0,5 км), (0,116

г/см3 ; 1 км), (0,263 г/см3 ;

1 км) и (0,286 г/см3 ;

10 км).

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

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

{Vz}