Файл: Автоматизированная система обработки и интерпретации результатов гравиметрических измерений..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