книги из ГПНТБ / Коваль Л.А. Автоматизированная система обработки данных магниторазведки с применением ЭВМ (АСОМ-АМ)

Таблица 4

2.0100-0105 Константы

3.0106-0256 Подпрограммы

4.0257-0262 Заплата к первой части МИС-15

5.0263-0266 Свободнн

6.0267-0635 Подпрограммы

7.0636-0706 Первая часть МИС-15 (исправление и пополнение

ШАПКИ)

8.0707-1077 Подпрограммы

9.II00-I374 Константы

12.1473-1475 Свободны

13.1476-1477 Запрещены для неоперативного использования

14.1500-2371 МИС-И ("Матрица")

15.2372-2377 Свободны

16.2400-2477 МйС-19 (исправление уровня матрицы)

17.2500-3477 Свободны

18.3500-4170 МИС-13 (тригонометрическое сглаживание, анализ)

19.4171-5477 Свободны

20.5500-6010 MHC-I4 (тригонометрическое сглаживание числового

материала) 21. 60II-6I57 Свободнн

22. 6160-6776 МИС-17, 18 (сглаживание по Савинскоыу), первая часть

23.6777-7670 Подпрограммы

24.7671-7776 МИС-17, 18 (сглаживание по Савинскому), вторая

часть.

Кроме того, в системе АСОМ-AM используется самостоятельная программа КИ [8] - МИС-12 (сы. блок-схему). В следующую таблицу сведены номера ЛПМ, используемых программами системы МИС.

100

Исправленная матрица располага­ ется в произвольном месте НМД. по указанию обработчика

101

этот восьмеричный код используется для вызова значений поля в предыдущее и последующей точках. Например, если первая точка про­ филя находится в ячейке 10000, то значение дТ п р е д затем зано­ сится в рабочую ячейку 0044 по команде -10 56 0000 0044, а значе­

Рабочие ячейки подпрограммы-0040-0047, индексные-0016, 0017. Исходные данные по выходе из подпрограммы сохраняется.

12. Интерполяция на ошибочные точки в строке матрицы и на профиле (ячейки 0753-1025, 1032-1077). Подпрограмма производит занесение на место ошибки в профиле значения, получающегося ли­ нейной интерполяцией между бллкайшими верными точками. Коды оши­ бок, находящиеся на краю профиля, заменяются бликайшим верным значением. Ограничение: не должно быть целиком ошибочной строки или профиля.

Перед обращением к СП в ячейку 0040 заносится код

+ n - I 0000 0000, где п - число точек в строке или на профиле. Варианты:

а) профиль находится в ячейках 2 блока М03У, начиная с ячей­ ки 14000.- Обращение -31 00 0760 0013;

б) профиль находится в ячейках I блока МОЗУ, начиная с ячей­ ки 1500 . Обращение -31 00 1032 0013.

В подпрограмме используются индексные ячейки 0013-0017, стандартные рабочие 0040-0047, СП-2 - линейная интерполяция.

13. Константы системы (ячейки II00 - I367) . Укажем особо на коэффициенты сглаживающего фильтра (8-23) расположенные по стро­ кам в ячейках I I 0 0 - I I 6 6 .

автономный блок исправления и пополнения ШАПКИ: а) в Д0ПИНФЕ-4-бис на. 2 ЛПМ;

б) в Д0ШШФЕ-4 на I ЛПМ. Числовой материал к МИС-15 готовят­ ся точно так же, как и к программе ПЗР0-8 (см. § 3, пункт "В").

103

Обратим внимание на пять последних чисел ШАПКИ с их адресами для ввода, которые в системе ПЕРО на использовались.

оси ох , дс] - шаг в метрах в направлении оси оу. 4062 - б - среднеквадратичный уровень для сглаживания .

4063 - Я - параметр сглаживания-лшнимальный для синтеза полупериод гармоники .

4064 - К - число итераций для сглаживания по Савинскому.

Как и прежде ШАПКУ можно вводгть заново всю, тогда числовой материал должен готовиться так, как это делалось для ПЕРО-2, но можно исправлять и пополнять отдельные числа - они набиваются группами со своими адресами.

И с п р а в л е н и е Д 0 П И Н Ф А - 4 - б и с н а 2 Л П М . После установки числового материала у читающего уст­ ройства и передачи управления в ячейку 0636 машина производит дважды ввод числового материала, распечатывает все 128 строк ШАПКИ, исправляет Д0ПИНФ-4-бис. Останов СчАК 0662 - неисправен

ввод. Если пуск! - вводы повторяются. Останов СчАК 0707 - конец работы программы.

И с п р а в л е н и е Д О П И Н Ф А - 4 н а I Л П М . После установки числового материала у читающего устройства и пе­ редачи управления в ячейку 1420 производится дважды ввод числово­ го материала, распечатывается ШАПКА и исправляется ДОПИНФ-4. Ос­ танов СчАК 1444 - неисправен ввод. Если пуск! - вводы повторяются. Останов СчАК 1473 - конец работы программы. В МЙС-15 используются СП: 5 - обмен с НМЛ, 4 .- групповой перевод 10—2.

В. ЫИС-П "Матрица" (ячейки 1500-2371)

Программа обеспечивает построение прямоугольной матрицы зна­ чений дТ по сглаженным (в результате работы программ ПЕРО-10 или МИС-14, 18) значениям поля, заданным по профилям. В програм­ ме используется ключ 0001 (включен, если матрица строится по , фрагменту, выключен, если по всему участку). Используются: ЛПМ 2- исходный числовой материал, ДОПИНФ-4-бис и профили, в нулевую и первую зоны переписывается МИС; ЛПМ 3 - на этой ленте по. строкам, начиная с нулевой ячейки второй зоны,помещается матрица, в нуле­ вую зону заносится некоторый табличный материал. Ограничения на

104

размеры матрицы соответствуют ограничениям ггрограмш КИ [ 8 ] : р <

^ 2000I 0 , p * q $ i n f (600001 0 , 2000 * ( n - 2 ) ) . В последнем огра­

0100-7776. Останов (с замком) СчАК 1510. Сумма не совпала с эта­ лонной, хранящейся в ячейке 0100. (PI) = 0100, (См) = новая К Е . СчАК 1512, сумма совпала - обязательный останов.

2.I 5 I 2 - I 5 I 5 . Переписывание системы на ЛПМ 2 в 0 и I зоны. СчАК 1516 - санитарный останов, на случай использования только начала программы "Матрица". Используется СП 5 - обмен с НМЛ.

3.I5I6-I747. Вспомогательные расчеты, переработка каталож­ ного материала и др. Используется ключ 0001. Перед пуском (СчАК

Матрица числовых значений поля будет использована для пост­ роения карты изодинам д Т , поэтому резонно автоматизировать под­ готовку Д0ПИНФА к программе "Карта изолиний" [ 8 ] . В блоке перфо­ рируется (и одновременно распечатывается) таблица в десятичной системе с "запятой": условное число (год съемки) p,q, др и &Ц ( в мм с точностью до второго знака после запятой), число 4 (но­ мер шрифта), I - число изодинам (по ограничениям ШАПКИ не более 2 I J Q ) , число 0, далее, последовательно, градации изодинам. За по­

ym i „ , xm a „ , y m o x матрицы (исходная система координат - обратный переход осуществляется по формулам(5-17)). После-чим пя­ тым числом в этой группе за разделающим рядом минусов печатается число Т о м - уровень нормального поля, от которого исчисляется матрица. СчАК 1736 - останов для проомотра распечатки, отключения перфоратора (можно выключить и ключ А I X Используются СП: 5 - обмен с ТС1Е; 10 - перевод целого двоичного числа в 2 с.с. с пла­ вающей запятой; 9 - перевод 2—10, с кодовой запятой и заданной точностью; 7 - перевод 2—10 с плавающей запятой массива чисел.

4. 1750-2012. Сокращение профилей из КАТАЛОГА для дальней-

105

шего поиска (в цикле по строкам матрицы). Программа осуществляет определение значений поля в узлах матрицы по строкам ("горизон­ тальным рядам"). Для заданного у ряда и известных координат Хрип и хтак из дальнейшего поиска исключаются профили, у которых:

1)|хг-х,1 < др ;

i-nf (ivy|,iv,-v!)> 2 u q-

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

5.2013-2243, 2333-2371. Определение искомого значения поля

вточке матрицы связано с поиском профилей исходного числового материала по правилам, которые были изложены выше. Работа с маг­ нитной лентой накладывает свои ограничения на программу с учетом технических возможностей и экономичности работы. Программа состав­ лена таким образом, что после нахождения профиля, с которым точка совпадает по топографическому допуску, этот профиль вызывается в МОЗУ и хранится там до работы со следующей точкой в строке. По­ следняя сразу проверяется на совпадение с ранее найденным профи­ лем и при благоприятном исходе проверки нового обращения к ленте не происходит.

Если в числовом материале отсутствует профиль, проходящий через точку, то, как было сказано выше, отыскиваются два ближай­ ших с обеих сторон профиля, которые для работы с ниш также, вы­ зываются в М03У. При переходе к следующей точке строки поиск новых пар повторяется и затем производится сравнение характерис­ тик старых и новых профилей. Если профили совпадают, то нового обращения к магнитной ленте не производится.

10S

Случай с одним профилем, находящимся на расстоянии меньше 0,5 Щ от искомой точки, является промежуточным между двумя разобранными выше, поиск профилей осуществляется как во втором варианте, присвоение же значения как в первом.

5-1. 2013-2037, 2333-2345. Проверка: имеются ли среди най­ денных прежде профилей такие, расстояние которых до строки матри­ цы меньше топографического допуска, находится ли искомая точка внутри такого профиля. При положительном ответе переход к присвое­ нию одного значения. При отрицательном ответе - переход к следую­ щему пункту. Используется СП 2 - линейная интерполяция.

5-2. 2040-2104. Поиск наилучшей пары профилей по каталогу для данной точки. Используется СП 2 - линейная интерполяция.

5-3. 2105-2142. Анализ найденной пары профилей - последова­ тельно:

1.Не найдено ни одного профиля с обеих сторон - переход к присвоению кода ошибки; иначе -

2.Среди пары профилей имеется профиль, который включает в себя искомую точку - переход к присвоению одного значения; иначе -

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

снаблюденных профилей на искомую точку; иначе -

4.Профили проверяются: нет ли среди них такого, который со­ ответствует условию:

| л у | < 0 . 5 - A q ,

если есть, то переход к присвоению одного значения, иначе пере­ ход к присвоению кода ошибки.

ной проверкой - вызывался ли на предыдущем шаге найденный профиль с магнитной ленты. Если вызывался, то новое обращение к ленте обходится. Используются СП: 5 - обмен с НМЛ; I I - определение точек на профиле для интерполяции; 2 - линейная интерполяция.

5-6. 2173-2237, 2355-2356. Интерполяция значений, снятых с двух профилей (с уточнением на дТ0 ), на искомую точку с пред­ варительной проверкой: вызывались ли на предыдущем шаге оба про­ филя с магнитной ленты. Если вызывались, то новое обращение к

107

ленте обходится. Используются СП: 5 - обмен с НШ; 12 - определе­ ние точек на профиле для интерполяции; 2 - линейная интерполяция.

5-7. 2240-2243. Цикл по точкам в ряду матрицы.

6. 2244-2306. Цикл по строкам матрицы. Циклическое суммирова­ ние исчисленной строки. Проверка: не получилась ли целиком ошибоч­ ная строка. Бели да, то останов СчАК 2254 (с замком), в сумматоре код, который указывает на номер строки матрицы, оказавшейся среди чрезмерно разреженных наблюдений (например, когда пропущены три и более профиля съемки). Код в сумматоре расшифровывается так: если

Перед переписыванием строки на ее место на ленте коды ошибок заменяются на интерполированные (по строке) значения. Напомним, что первый элемент первой строки матрицы записывается в нулевую ячейку 2 зоны 3 ЛПМ. Номер этой ячейки равняется +0000 0301 0000, номер первой свободной на этой ленте ячейки будет

+0000 0301 0000 +p»q . Кроме того, в описываемой части программы проводятся вычисления для дТт „ , дТт 1 п и дТс р ( ,д н в г 0 по участку. Используются СП: 13 - "забивание" и интерполяция ошибок; 5 - об­ мен с НМЛ.

7. 2307-2332. Заключительная часть программы. Печать в деся­ тичной системе с плавающей запятой дТ^,, дТт 1 п и д Т с р и двух восьмеричных кодов; циклической суммы всей матрицы и номера пер­ вой свободной ячейки за матрицей на 3 ЛПМ.

В ячейки 0000-0537 нулевой зоны ЛПМ 3 переписывается табли­ ца, содержание которой может быть полезным для расширений, сле­ дующих за программой "Матрица" (помимо КИ). Содержание таблицы:

108

щий за последней ячейкой матрицы .

0507 - +0000 0000 р, восьмеричное число точек по "тесовой"

0537 _ f несущественно

Числа в ячейках 0521-0532 в 2 с.с. с плавающей запятой. СчАК 2333 - останов.

Перед работой программы "Матрица" требуется задание коорди­

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

109