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

Г. МИС-19 (ячейки 2400-2477) - исправление уровня матрицы.

Для работы программы требуется числовой материал на перфо­ ленте: первое восьмеричное число - адрес первой ячейки на НМЛ для переписывания исправленной матрицы, второй десятичный "с запятой" код - значение дТ , которое нужно отнять из значений матрицы. Два этих кода в указанном порядке без адресов набиваются на пер­ фоленте в границах ввода. После передачи управления в ячейку 2400 производятся вводы числового материала, вызывается таблица матри­ цы, программа настраивается на параметры матрицы, вызывает и пере­ рабатывает числовой материал, переписывая его на ;/казанное место НМЛ. Останов СчАК 2410 - не совпали циклические суммы двух вводов числового материала. Если пуск!, вводы повторяются.

Печать программы состоит из двух чисел: первый код - цикли­ ческая сумма новой матрицы, второй код - номер ячейки на ленте, следующий за последним числом новой матрицы. Останов СчАК 2500конец работы программы. Используются СП: 4 - перевод 10—2; 5 - обмен с НМЛ. Для работы программы необходимы ЛПМ 3, а также, если перерабатываемая матрица не размещается на этой ЛПМ, то и другой лентопротяжный механизм по указанию обработчика.

Д. МИС-12 - подключение программы КИ .

Программа КИ предназначена для графического вывода из ЭВМ "Минск-2" (с двумя блоками МОЗУ) числовой информации о рельефе любой функции в виде карты изолиний с постоянным или переменным сечением изолиний, с вычерчиванием рамок карты и выполнением оцифровки изолинии. Числовая информация о рельефе функции пред­ ставляет собой матрицу значений этой функции в узлах прямоуголь­ ной сети, записанных на магнитной ленте [ 8 ] . Для детального озна­ комления с работой программы КИ читатель отсылается к инструкции Л.Г.Перфильева и B.C. Крамака "Вычерчивание карты изолиний на

ЭВМ "Минск-2" с помощью устройства АТЛАС"[8]. Уточним следующие обстоятельства.

1.Номер начального слова матрицы +0000 0301 0000 или тот, который был выбран для исправленной матрицы в МИС-19.

2.Д0ПИНФ к КИ, полученный в результате действий "Матрицы", дополняется границами и перед первым числом пробивается адрес 5000 и код передачи адреса.

111

+ 0 0 0 0 0 1 9 7 1 + 0 0 0 0 0 0 0 ? f + 0 0 0 0 0 0 0 2 1 + 0 0 0 0 0 2 , 0 0 + 0 0 0 0 1 0 , 0 0

+ 0 0 0 0 0 0 0 6 0 + 0 0 0 0 0 0 1 00 + 0 0 0 0 0 0 1 4 0 + 0 0 0 0 0 0 2 0 0 + 0 0 0 0 0 0 3 0 0 + 0 0 0 0 0 0 4 0 0 + 0 0 0 0 0 0 6 0 0

+++++++++

+ 8 9 3 3 9 9 9 + 0 4 + 7 3 4 9 9 9 9 04

+ 1 0 0 0 0 0 0 + 0 5 + 8 0 4 9 9 99 + 04

+ S 4 4 5 0 0 0 + 0 5

ВЫЧЕРЧЕНО УСТРОЙСТВОМ „ Д Т Л Я С " 1 7 / 1 - 7 - 7 2

Рис. 34.

IS-2629

113

114

+ 0 0 0 0 0 197 1 + 0 0 0 0 0 0 0 6 6 + 0 0 0 0 0 0 0 7 1 + 0 0 0 0 0 2 , 0 0

+ оооо10,00

+ 0 0 0 0 0 0 0 0 1 + 00000002 - 1

+ 0 0 0 0 0 0 0 0 0

-ooooooiоо

- 0 0 0 0 0 0 0 6 0

- O O O O D O O I O

+ 8 9 9 9 9 9 9 + 04 + 8 0 4 9 9 9 9 + 0 4

+ 1 0 0 0 0 0 0 + 05 + 8 6 9 9 9 9 9 + 04

+ 5 4 4 5 0 0 0 + 0 5

ВЫЧЕРЧЕНО УСТРОЙСТВОМ „ДТЛЯС 2 7 / 1 - 9 - 7 2

Рис. 36.

115

3. В случае, если ДОПИНФ, приготовленный "Матрицей", не под­ ходит (не хватает градаций, нужен другой номер шрифта, необходимо изменить масштаб, жестко задаваемый в системе и т.п.), нетрудно приготовить новый ДОПИНФ вручную.

На рис. 33, 34, 35 и 36 приведены образцы карт графиков и карт изодвнам лТ , вычерченных с использованием программы Ш на АТЛАСЕ. Сеть значений для каждого из образцов в масштабе изобра­ жений равна 2 х 10 мм или I х 10 мм, т.е. отличаются большой ани­ зотропией. Исходные значения поля не подвергались какому-либо сглаживанию, кроме того, который здесь именуется первым уровнем (сглаживания). Между тем карты имеют удовлеторительный рабочий вид, хотя и возможно, нуждаются в некоторой редакции.

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

Е. МИС-ГЗ. Тригонометрическое сглаживание-анализ (ячейки 3500-4170)

Программа выполняет гармонический анализ заданного профиля с выбором параметра сглаживания и наглядным изображением резуль­ татов анализа (описание алгоритма см. в § 8, раздел "В"). В про­ грамме используется ключ 0001 (включен,если анализ ведется с ри­ сованием; выключен, если рисовать не надо), исходный числовой ма­ териал на ЛПМ I . Номер профиля дается с пульта на клавишном набо ре. Перед обращением к МИС-13 ШАПКУ следует пополнить уровнем сглаживания

I . 3500-3513. Вызов с ЛПМ I ДОПИНФА-4, снятие о клавишного набора номера анализируемого профиля. Поиск в ДОПИНФЕ-4 информа­ ции о заданном профиле. Предусмотрен останов СчАК 3513 (с замком), если набран номер, отсутствующий в КАТАЛОГЕ: (СМ) = (клавишный набор). Используется СП 5 - обмен с НМД.

11Б

2.3514-3566. Определенно со (8-17) числа палеток г, умещаю­ щихся в профиле, определение по (8-18) и (8-19) местоположения и числа точек в каждой палетке. Используется СП 10 - перевод целого двоичного числа в 2 с.с. с плавающей запятой.

3.3567-3625. Начало цикла по г (число палеток). Вызов (во­ зобновление) профиля в ГОЗУ, забивание кодов ошибок интерполиро­ ванными значениями. Формирование и выдача на печать значений дТ,

относящихся к очередной палетке. Используются СП: 5 - обмен с НМЛ; 12 - ш1терполяция на ошибочные точки; 7 - перевод 2—10 с плавающей запятой массива чисел.

4.3626-3651. Вычисление по формулам (8-8) и (8-7) линейного фона и д. Последняя величина выдается на печать. Используются СП: 10 - перевод целого двоичного числа в 2 с.с. с плавающей за­ пятой; 7 - перевод 2—10 с плавающей запятой массива чисел.

5.3652-3703. Вычисление по (8-10) значений Ь(К) и выдача их на печать в 10 с.с. с плавающей запятой. Используются СП:

10 - перевод целого двоичного числа в 2 с.с. с плавающей запятой; 16 - stnx(cosx) ; 7 - перевод 2—10 с плавающей запятой массива чисел.

6.3704-3753. Определение из условия минимума (8-16) номера последней оставляемой при синтезе гармоники m , величины R ;

выдача на печать п" (число точек в палетке), m и R . Исполь­ зуются СП: 10 - перевод целого двоичного числа в 2 с.с. с плаваю­ щей запятой; 7 - перевод 2—10 с плавающей запятой массива чисел.

7.3754-4002. Вычисление по ( 8 - I I ) и (8-12) сглаженной кри­ вой. Выдача на печать сглаженной кривой. Используются СП: 16 - sinх(cosх) ; 7 - перевод 2—10 с плавающей запятой массива чисел.

8.4003-4005. Переход по ключу к рисованию на АТЛАСЕ. Цикл

по г (количество палеток). Останов СчАК 4006.

9. 4010-4170. Рисование результатов гармонического анализа на АТЛАСЕ. Используются СП: 14, 15 - "Линия", "Надпись"; 9 - пе­ ревод 2—10 с кодовой запятой и заданной точностью; 3 - 'График". Предусмотрен останов СчАК 4011 - в сумматоре порядковый номер па­ летки: включить АТЛАС, перописец в левый нижний угол планшета, пускI

Выдачу на печать МИС-13 открывают два десятичных числа, на­ печатанных в форме "с запятой". Первое число - номер профиля, второе - номер палетки. Далее через интервал печатается п* де­ сятеричных чисел с плавающей запятой, составляющих палетку

117

анализируемой величины.

Графический вывод на АТЛАС (рис. 27) открывается заголовком, в котором после номера профиля указывается номер палетки. Во вто­ рой строке заголовка указываются значения N, М, R , масштаб графиков.

Графики - исходный, сглаженный и коэффициентов Ь(К) вычер­ чиваются с шагом 0,5 см между точками. На первых двух графиках оцифровка означает номера точек (в палетке), на последнем - но­ мера гармоник. На сглаженном графике первая и последняя точки, не входящие в сглаживание, не показываются. Аналогично, нулевая и последняя гармоники, не имеющие смысла, не показываются на графике коэффициентов Ь(К) .

Ж. МИС-14. Тригонометрическое сглаживание числового материала (ячейки 5500-6010)

Программа обеспечивает тригонометрическое сглаживание профи­ лей в пределах участка съемки или фрагмента (по включенному ключу 0001). При сглаживании по фрагменту профили, не имеющие общих то­ чек с полосой y m l n , y m m t , подвергаются первому уровню сгла­ живания (линейная интерполяция на ошибочные точки). Исходный чис­ ловой материал на ЛПМ I,сглаженные профили вместе с Д0ПИНФ0М-4бис переписываются на ЛПМ 2. Перед обращением к МИС-14 ШАПКУ следует пополнить параметром сглаживания R и, если расчеты будут вес­ тись по фрагменту, уточнить координаты его угловых точек.

1. 5500-5505, 5750-5771. Вызов с ЛПМ I ДОПИНФА-4, определе­ ние полосы ymi.„, ymox,включающую фрагмент. Используется СП-5

-обмен с НМЛ.

2.5506-5514, 5772-6010. Начало цикла по профилям. Вызов с М03У профиля с ЛПМ I . Первый уровень сглаживания и определение (при включенном ключе 0001) имеет ли профиль общие точки с поло­

сой y m l n , y m e x . Если имеет - переход к следующему блоку; ина­

116

и числа точек в каждой палетке. Используется СП 10 - перевод це­ лого двоичного числа в двоичную систему счисления с плавающей запятой.

4. 5577-5716, 5745-5747. Цикл по г (число палеток) внутри профиля. Выделение палеток, вычисления по формулам (8-8) и (8-7)-

и (8-12) сглаженной кривой. Композиция с учетом перекрытий сгла­ женного проф;ля. Используется СП: 10 - перевод целого двоичного числа в 2 с.с. с плавающей запятой; 16 - sinx (cosx).

5.5717-5736. Циклическое суммирование сглаженного профиля. Запись профиля на ЛПМ 2. Исправление КАТАЛОГА для Д0ПИНФА-4-бис. Цикл по профилям. Используется СП 5 - обмен с НМЛ.

6.5733-5744. Заключительный блок в программе: циклическое суммирование Д0ПИНФА-4-бис, переписывание его на ЛПМ 2. Исполь­ зуется СП 5 - обмен с НМЛ. СчАК 5745 - останов.

3. МИС-17, 18. Итерационное сглаживание (ячейки 6160-6776, 7671-7776)

Две программы МИС-17 (анализ сглаживания по Савинскому) и МИС-18 (сглаживание числового материала по Савинскому) объединя­ ются здесь в одном описании в силу того, что имеют ряд общих блоков, выполненных в стиле СП, в отличие от более автономных

-31 00 6640 0017 в ячейку 0040 заносится код +п-1 0000 0000, где п - число точек в сглаживаемой кривой, а сама кривая распо­ лагается в ячейках МОЗУ, начиная с 1500 (действует обычное в си­ стеме ограничение - не больше I 02 IJQ точки). На выходе коэффици­

т о ч к а м (ячейки 7671-7712). Блок осуществляет однократное сглаживание по формуле (8-28) числового материала, располагающе­

)