Файл: Коваль Л.А. Автоматизированная система обработки данных магниторазведки с применением ЭВМ (АСОМ-АМ).pdf

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

оси ох , дс] - шаг в метрах в направлении оси оу. 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