Файл: Автоматизированная система обработки и интерпретации результатов гравиметрических измерений..pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 207
Скачиваний: 0
параметры заданы верно п функция удовлетворяет условию задачи, то min e s ' s _ 1 будет найден и вслед за триадами на печать будет выда ваться регуляризированное приближение, соответствующее min es -s _ 1 . Далее управление автоматически передается на счет приближений нового уровня.
Если для некоторого уровня решение отсутствует, то на конец счета указывает останов по команде 77 в ячейке 0512. Возможен «ложный» останов в ячейке 0512, не означающий отсутствия реше ния в следующих случаях:'
1. При продолжении на первый уровень заданный шаг Да ариф метической сетки велик, а локальные минимумы и максимумы функ
ционала |
е находятся на отрезке |
[10 - 1 °, |
Да], на отрезке |
[Да, 103 ] |
||||||
функционал е монотонно убывает. Тогда нужно повторить |
счет,, за |
|||||||||
дав достаточно маленький шаг Да и а 0 |
^ |
5Да, посмотреть на пове |
||||||||
дение |
функционала |
е вблизи a m i n = Ю - 1 |
0 . |
|
|
[10~1 0 , 103 ]. |
||||
2. |
Функционал |
е монотонно |
возрастает на отрезке |
|||||||
Тогда |
вполне возможно, что локальный |
минимум |
лежит вне этого |
|||||||
отрезка справа. В данном случае нужно |
увеличить |
а т |
а х |
— ввести |
||||||
перфокарту 50, с помощью которой а т а х |
увеличится |
на порядок и |
||||||||
счет продолжится. |
|
|
|
|
|
|
|
|
||
В модификации сглаживания автоматический поиск |
min e s ' s _ 1 |
|||||||||
ведется на уровне |
задания исходной функции (гг |
= 0). |
Вслед за |
|||||||
триадами (as , ßs , e s > s - 1 ) выдаются на печать сглаженные |
значения |
|||||||||
исходной |
функции. |
На конец |
счета |
указывает |
останов |
в ячей |
||||
ке 0713. |
|
|
|
|
|
|
|
|
|
Возможен останов по команде 77 в ячейке 0523 при расшифровке компактной записи, когда размеры числовой матрицы исходной функции не ' соответствуют заданным в перфокарте информации о счете.
При работе СП-0154 (вычисление изолиний и координат) в случае,
когда |
число |
изолиний |
превышает |
100, автоматически |
передается |
|
управление на выдачу результатов |
в виде числовых матриц. |
|||||
|
|
б) Расчеты с использованием магнитной ленты |
||||
Составляются два пакета к счету: |
|
|||||
Пакет |
1 |
собирается |
из перфокарт: «вызов ИС-2»; |
чистая; 51; |
||
чистая, |
2 |
|
48 — ОсП8; |
чистая; МЗ. |
|
|
Пакет |
2 |
собирается |
из перфокарт: «Вызов ИС-2»; |
чистая; 52 |
и 53 — вызов ОсП8 с МЛ; чистая; МЗ; чистая; M l .
Вводится пакет 1. Записывается ОсП8 на МЛ и машина останав ливается по команде 77 в ячейке 0013.
С помощью ВсП8 готовится матрица заданных размеров и запи сывается на МБ (см. гл. X I X , раздел 1, п. 4°).
Вводится пакет 2. Вызывается с МЛ ОсП8, настраивается по параметрам, вызывается матрица исходных данных с МБ и пере дается управление на счет. Счет протекает так, как описано выше.
Вид выдачи результата зависит от решаемой задачи:
12 Заказ 76 |
177 |
1. Если функция продолжается в область |
верхнего полупро |
||||||||||
странства (гг < 0 ) , то на печать |
будет выдаваться |
результат |
после |
||||||||
довательно |
для каждого уровня. |
|
|
|
|
|
|||||
2. Если |
функция или ее вычисленная первая |
производная про |
|||||||||
должается в- область нижнего |
полупространства |
( z > 0 ) , |
то для |
||||||||
каждого |
уровня |
на печать |
последовательно |
будут |
выдаваться: |
||||||
а) триады |
чисел as , ßs , es -s _ 1 , |
которые являются |
информацией |
||||||||
о характере |
последовательных |
приближений |
и |
поведении |
нормы |
||||||
es -s _ 1 . Печать ' триад прекращается, когда es >s _ 1 |
достигает |
мини |
|||||||||
мума и начинает расти на восходящей ветви графика; |
|
|
|||||||||
б) регулярпзовашюе приближение, отвечающее min es >s _ 1 , т. е. |
|||||||||||
решение |
задачи для данного |
уровня z,-; |
|
|
|
|
|||||
в) два числа — оптимальные |
параметры а о п , ß o n , при которых |
||||||||||
получен |
минимум |
нормы. |
|
|
|
|
|
|
|
||
Если для какого-нибудь zt минимум нормы отсутствует, то про |
|||||||||||
изойдет «авост» без выдачи приближения. |
|
|
|
|
|||||||
3. При сглаживании функции на уровне zx |
= 0 на печать будет |
||||||||||
выдана |
информация, аналогичная пп. «а» п «б». Вслед за сглажен |
||||||||||
ным полем будут напечатаны четыре |
числа: ссо п , ß o n , Б, а, после чего |
||||||||||
произойдет |
останов. |
|
|
|
|
|
|
|
|||
Норма s есть модуль максимального отклонения |
сглаженной |
||||||||||
функции от исходной: |
|
|
|
|
|
|
|
||||
|
|
|
|
e = max|V(z, |
у)сгл |
— Ѵ(х, j/)„c x |. |
|
|
|
||
|
|
|
|
- V, У |
|
|
|
|
|
|
|
Норма о 2 характеризует среднеквадратпческое уклонение сгла женной и ИСХОДНОЙ функций
M N
j = l i=l
Значения результативной функции в зависимости от содержания седьмой строки M l (информация о счете) могут быть выданы в двух видах:
1. Если ?і =_0, то результат будет напечатан в виде матриц NM значений поля по равномерной сетке s, где построчная последователь ность чисел сохраняется той же, что была в исходной функции. Ман тисса чисел укорочена и состоит из шести цифр. Три последние деся тичные цифры содержат порядковый номер этого числа в строке прямоугольной матрицы.- (Выдача по СП-0161).
Графическое представление такого результативного поля сво
дится к переписыванию матрицы значений на |
миллиметровку |
в масштабе: 1 см соответствует шагу s. При этом |
рамку матрицы |
178
удобно предварительно разметить через 1 см: начиная с левого верх
него угла, вправо пишутся номера по строке от нуля до М~і. |
В двух |
|||
мерном случае |
результат выдается |
в виде значений |
по |
профилю. |
2. Если п = |
Сп, где С„ > 0, то |
результативное |
поле |
выдается |
в виде значений изолиний с сечением Сп и координат этих изолиний.
Порядок печати, рекомендации по нанесению координат |
изолиний |
на миллиметровку и построению поля изолиний подробно |
изложены |
в описании к СП-0154 (см. гл. X I X , раздел 2). Функции двух пере менных выданные по профилю в виде координат изолиний удобно строить в так называемой вертикальной плоскости, для чего исполь
зуется несколько |
уровней |
нижнего |
полупространства, нанесенных |
||
в |
масштабе |
один |
под другим. Координаты равных значений поля |
||
на |
уровнях |
соединяются |
плавными |
линиями. |
В табл. 69 приведены исходные данные к задаче сглаживания.
+ - + |
01 |
521 |
528 |
534 |
|
+ - + |
01 |
539 |
546 |
543 |
|
+ - + |
01 |
551 |
544 |
543 |
|
+ - + |
01 |
539 |
534 |
534 |
|
+ - + |
01 |
537 |
534 |
532 • |
|
+ - + |
01 |
544 |
544 |
542 |
|
+ - + |
01 |
543 |
545 |
545 |
|
+ - + |
01 |
545 |
535 |
535 |
|
+ - + |
01 |
533 |
533 |
533 |
|
+ - + |
01 |
530 |
521 |
533 |
|
+ - + |
01 |
523 |
519 |
525 |
|
+ - + |
01 |
526 |
520 |
523 |
|
+-+ |
01 |
519 |
515 |
505 |
|
+ - + |
01 |
498 |
503 |
497 |
|
+ - + |
01 |
490 |
477 |
470 |
|
+ - + |
01 |
401 |
466 |
459 |
|
+ - + |
01 |
463 |
462 |
452 |
|
+ - + |
01 |
450 |
440 |
439 |
|
+ - + |
01 |
435 |
431 |
414 |
|
+ - + |
01 |
418 |
411 |
400 |
|
+ - + |
01 |
397 |
380 |
382 |
|
+ - + |
01 |
369 |
356 |
.349 |
|
+-+ |
01 |
337 |
336 |
335 |
|
+ - + |
01 |
330 |
336 |
333 |
|
+ - + |
01 |
323 |
310 |
313 |
|
+-+ |
01 |
313 |
297 |
282 |
|
+ — |
+ |
01 |
290 |
275 |
271 |
H |
Ь |
01 |
260 |
251 |
257 |
H—+ |
01 |
267 |
265 |
260 |
|
+ — |
+ |
01 |
247 |
240 |
245 |
H—4- |
01 |
252 |
241 |
242 |
|
H |
h |
01 |
255 |
260 |
259 |
+-+ |
01 |
259 |
261 |
273 |
|
+-+ |
01 |
273 |
281 |
291 |
|
|
|
|
|
|
|
Т а б л и ц а . |
69 |
+ - + |
|
01 |
303 |
305 |
305 |
+ - + |
|
01 |
302 |
302 |
303 |
+ - + |
|
01 |
314 |
307 |
302 |
+ — h |
01 |
306 |
327 |
335 |
|
+ - + |
|
01 |
341 |
342 |
357 |
+ - + |
|
01 |
368 |
368 |
386 |
+—h |
01 |
397 |
423 |
428 |
|
+ - + |
|
01 |
453 |
483 |
511 |
+ — + |
01 |
541 |
580 |
595 |
|
1 — + |
01 |
649 |
655 |
674 |
|
01 |
691 |
708 |
728 |
||
+-+ |
01 |
766 |
785 |
783 |
|
+ - + |
|
01 |
798 |
807 - |
802 |
+ - + |
|
01 |
821 |
828 |
827 |
+-+ |
|
01 |
838 |
829 |
819 |
-1 |
h |
01 |
816 |
816 |
819 |
+ - + |
|
01 |
799 |
791 |
800 |
+ - + |
|
01 |
802 |
7S4 |
781 |
+ - + |
|
01 |
766 |
760 |
750 |
+ - + |
+ |
01 |
741 |
736 |
742 |
+ - |
01 |
722 |
712 |
708 |
|
+ — h |
01 |
695 |
695 |
678 |
|
+j - + |
1_ |
-01 |
665 |
676 |
679 |
01 |
680 |
678 |
688 |
||
-1—+ |
|
01 |
668 |
668 |
665 |
+ — h |
01 |
669 |
677 |
678 |
|
+ - + |
|
01 |
680 |
672 |
668 |
+—V |
01 |
665 |
685 |
678 |
|
' + - + |
|
01 |
665 |
641 |
631 |
+ - + |
|
01 |
629 |
609 |
596 |
+ - + |
|
01 |
590 |
583 |
578 |
_j—_|_ |
01 |
572 |
559 |
543 |
|
+ - + |
|
01 |
532 |
519 |
504 |
П р п ы е ч а н п е . В таблице приведены 67 значений Ѵг по профилю |
|
с одинаковым шагом. |
|
12* |
179 |
|
В табл. 70 приведены результаты расчетов по рассматриваемому профилю по Ос118. 2. Время расчета с момента ввода пакета 2 со ставляет S мин 40 с.
Т а б л и ц а 70
++-
++-
+++
++-I—! -
i i
+ ++
+++
+++
++-
+++
_j i
+- +
+- +
+-+
+-+
+- +
+- 4 -
тг
+1 -+>
тг
- 1 — ] — U
++4 -
++ +
++-
01 |
999999999 |
01 |
004975124 |
01 |
284904195 |
00 |
200000000 |
01 |
009950248 |
01 |
•182885804 |
01 |
895430240 |
01 |
004454876 |
00 |
321935339 |
01 |
984973264 |
01 |
004900364 |
00 |
322334174 |
01 |
535553000 |
01 |
535917001 |
01 |
536847002 |
01 |
379822060 |
01 |
373234061 |
01 |
365833062 |
01 |
550211198 |
01 |
546478199 |
01 |
545132220 |
01 |
895430240 |
01 |
004454876 |
01 |
268354025 |
00. |
253113941 |
Сглаженные значения исходной функции
а0 п
ßon
е
0
На рис. 35 показана часть исходной (табл. 69) и сглаженной кри вых.
В табл. 71 приводятся результаты пересчета Ѵг в нижнее полу пространство на модельном примере (исходные данные являются результатом расчетов по ОсП5 по тесту, рассмотренному в гл. X I X ) . Время расчетов 3 мин.
На рис. 36, а изображено поле аномалий силы тяжести на уровне z = 1,5 км. На рис. 36, б показано поле Ѵгг на уровне z—
=1 км.
Втабл. 72 приведены результаты расчетов Ѵгг в нижнем полу пространстве. Время расчетов 3 мин.
180