Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Фролов, С. А. Кибернетика и инженерная графика

.pdf
Скачиваний:
15
Добавлен:
20.10.2023
Размер:
7.75 Mб
Скачать

мых tg ф = отличаются от спектров базисных прямых лишь тем, что составляющие их растрэлементы распола­ гаются не в одной строке (столбце), а в нескольких. В этом случае количество строк (столбцов) равно величине чис­ лителя (знаменателя) дроби ~ (рис. 96).

Отмеченные выше зависимости установлены для пря­ мых, толщина которых равна линейному размеру растрэлемента —0,2 мм. При увеличении толщины линии эти зависимости сохраняются, причем значения N и 5 из­

меняются пропорционально

ц — коэффициенту неодно-

 

6

родности, который равен t] =

-j-,

где о — толщина линии,

d — диаметр светового пятна анализирующего луча. Стро­ гая закономерность в расположении растрзлементов, обра­ зующих прямую, и простое выражение этой закономер­ ности позволяют использовать свойства спектра прямой для автоматизации решения задачи о рассортировке коор­

динат растрзлементов по

их принадлежности

к

каждой

из заданных линий, когда

они пересекаются.

На

рис. 97

в дискретной форме представлены две пересекающиеся прямые / и II. Цифры показывают, в каком порядке будут записываться в а-массив растрэлементы в процессе обычного слежения. Чтобы суметь пройти пересечение, необходимо определить структуру прямой. Это можно

сделать в процессе слежения за ее

«чистым» участком.

Для составления спектра прямой

нужно выявить все

типы растрзлементов, образующих прямую, и установить существующую между ними внутреннюю связь.

На примере (рис. 97) показано, как с помощью спектра прямой проходится зона пересечения. Осуществление про­ граммы слежения машина начинает с растрэлемента /. При обследовании его окрестности установлено, что в зо­ нах А, В, D имеются свободные растрэлементы. Двигаясь по прямой в направлении, указанном стрелкой, обнару­

живают растрэлементы 3, 13, 16, не имеющие

в зонах

А,

В, D растрзлементов, принадлежащих линии,

за которой

осуществляется слежение. Поэтому растрэлементы /,

3,

13, 16 принадлежат к одному типу. Выяснив, в каких зонах растрэлементы участка прямой / имеют свободные растрэлементы, можно составить предварительный (не­ однородный) спектр прямой, который для рассматривае­ мого случая имеет следующий вид:

140

 

Признаки

 

 

 

Признаки

 

 

 

 

ряда

(свободные

М» растр-

(свободные

№ растр-

растрэлемен-

элементов

ряда

растрэлемен-

элементов

 

ты)

 

 

 

 

 

ты)

 

 

 

 

 

I

А,

В,

D

1, 3,

13,

IV

А,

В,

С,

D

10,

21

 

А,

Н

 

16 . . .

V

Е,

F,

G,

Н

9,18.

 

. .

I I

 

2, 4,

11,

V I

А,

В,

Н

7,

19

• •

 

 

 

 

14 . . .

JVII

A,

 

G, Н

8,

17 . . .

ГП

Е,

G,

II

6, 5,

12,

 

 

 

 

 

 

 

 

 

 

 

 

 

15 . . .

 

 

 

 

 

 

 

 

 

Для получения однородного спектра раскладываем длинные неоднородные ряды в однородные. С этой целью проверяем, не попали ли в один ряд растрэлементы, при­ надлежащие к разным однородным рядам.

•В

 

 

 

 

 

 

 

131ПЛП21

 

38

39

 

 

 

16

14

 

32

33

 

21

19

17

18

 

 

27

щ

22

20

 

 

 

 

 

 

23

 

 

 

40

34

 

 

 

 

 

 

41

35

 

30

31

 

 

 

42

 

 

36

37

 

 

 

• • • и

 

 

43][44]

 

 

 

 

 

 

 

 

 

 

I—It—1

 

 

 

 

 

 

 

Рис. 97

Проверку осуществляем с помощью соседних растрэлементов. Если эти растрэлементы одинаково располо­ жены по отношению к однородному ряду, то и рассматри­ ваемые также находятся в однородном ряду. Следуя таким

путем,

устанавливаем,

что у растрэлемента, координаты

которого

отличаются

от обследуемого

на

величину у

=

=

#обсл

+

1 (зона G — растрзлемент

4),

свободные

об-

 

141

ласти находятся в зонах А и Н. Для следующего растрэлемента первого ряда (растрзлемент 3) в зоне G располо­

жен растрзлемент 8, имеющий свободные области

в зо­

нах A, G, Н. Следовательно, растрэлементы 1 и 3,

хотя

и принадлежат к одному типу, не могут находиться в одном

однородном

ряду, так как

их положение по отношению

к соседним

растрэлементам

прямой различно.

Проверяем растрзлемент 13. В зоне G находится растрэлемент 14, для которого области А и Н свободные, т. е. растрэлементы 4 к 14 относятся к одному однородному ряду. Поступая таким образом, раскладываем неоднород­ ные ряды в однородные и получаем окончательный одно­ родный спектр прямой. Осуществляя слежение за «сво­ бодным» участком прямой, можно определить ее спектр г .

 

Признаки

 

растрэле­

 

Признаки

растрэле­

(свободные №

№.

(свободные

ряда

растрэлемен­

 

ментов

ряда

растрэлемен­

 

ментов

 

 

ты)

 

 

 

 

 

ты)

 

 

I

|

А,

В,

D

1,

13 . . .

V I I

Е,

F,

G, Н

9,

18 . . .

I I

3,

16 . . .

 

 

 

 

 

 

 

 

 

 

V I I I

А,

В,

Н

7,

19 . . .

 

 

 

 

2,

11 . . .

I I I

}

А,

Н

 

 

 

 

 

 

 

 

 

4,

14 . . .

IX

A,

G, Н

8,17 . . .

IV

 

 

 

 

 

 

 

 

6,

12 . . .

V

}

Е,

G,

Н

X

А,

В,

С, D

 

 

V I

5,

15 . . .

10,

21 . . .

 

 

 

 

Зная спектр, определяем, к какому его ряду относится тот или иной растрзлемент прямой. Это даст возможность «пройти» зону пересечения прямых.

Ниже на конкретном примере приводится методика, пользуясь которой, можно установить, к какому однород­ ному ряду спектра принадлежит исследуемый растрзле­ мент (рис. 97).

1. Растрзлемент 20 имеет свободные области в зонах Е, G, Н, поэтому он может принадлежать к V или V I рядам. Чтобы установить, к какому из них относится рассматри­ ваемый растрзлемент, проверяем, в каком ряду находится граничащий с ним растрзлемент, заведомо принадлежащий прямой. Для этого проверку переводим на зону С. В ней

1 Полученный спектр по внешним признакам отличается от спек­ тров, которые мы приводили раньше. Это объясняется тем, что нуме­ рация растрэлементов в данном случае принята иной, чем раньше.

142

расположен растрзлемент 17 (IX ряд). В зоне F от растрэлемента IX ряда лежит растрзлемент V ряда. Следова­ тельно, растрзлемент 20 входит в состав спектра прямой /.

 

2.

Растрзлемент 22 имеет свободные области в зонах А

и

Н.

Поэтому

можно

сразу утверждать,

что

он

входит

в состав спектра. Для

уточнения

принадлежности

к ряду

( I I I

или

IV)

выясняем,

какой

растрзлемент

находится

в

зоне

С.

Находим

растрзлемент 19,

принадлежащий

к

V I I I

ряду.

В зоне

F растрэлементов этого

ряда

лежат

растрэлементы

I I I ряда. Следовательно,

растрзлемент 22

также

входит

в состав

спектра

прямой

и принадлежит

кI I I ряду.

3.Растрзлемент 23 не имеет Е, G, Н, в зоне С от него находится растрзлемент 20, относящийся к V ряду. Сле­ довательно, растрзлемент 23 принадлежит к спектру и лежит в V I ряду.

Начиная с растрэлемента 24, мы вступаем в зону пере­ сечения, в которой картина расположения свободных растрэлементов существенно меняется, но накопленный к этому моменту опыт слежения позволяет безошибочно установить принадлежность растрэлементов этой зоны к спектру прямой, за которой велось слежение. Например, растрзлемент 24 имеет свободную область в зоне В. Та­ кого типа в спектре прямой / нет, но утверждать, что растрэлемент 24 не принадлежит прямой /, может быть прежде­ временно, так как возможно, что зоны А и D (растрзлемент со свободными областями в зонах А, В и D в спектре пря­ мой / есть) заполнены растрэлементами линии, пересекаю­ щей прямую /. Допустим, что растрзлемент 24 принадле­ жит прямой /. Чтобы убедиться в этом, достаточно уста­ новить, что в зоне С находится растрзлемент, не имею­ щий А, В, С, D. Аналогично можно установить, что растрэлементы 25, 26, 28, 29 принадлежат спектру прямой /. Дальнейшее слежение за линией I (выявление растрэле­ ментов 34, 35, 40, 41 и т. д.) не вызывает никаких затруд­ нений и осуществляется по обычной программе.

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

Свойства спектра прямой позволяют решить чрезвы­ чайно важную задачу «прохождения» области пересечения двух линий, причем логическая схема программы полу­ чается достаточно простой, но реализовать этот способ на существующих ЭЦВМ невозможно из-за ограниченного

143

объема их оперативной памяти. Для решения этой задачи необходимо, чтобы объем оперативной памяти ЭЦВМ в де­ сятки раз превышал существующий.

§ 15. О ТОЧНОСТИ МАШИННОГО РЕШЕНИЯ ЗАДАЧ

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

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

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

ряда

элементарных

приемов:

1 — отмерить по масштабу заданную величину цирку­

лем;

 

 

2 — поставить

ножку циркуля в точку пересечения

двух

прямых;

 

3 — провести окружность.

Выполнение каждого из этих приемов неизбежно свя­ зано с появлением ошибки (по данным Д. И. Каргина, величина ошибки составляет для приема 1—0,03—0,05 мм;

для

приема

2 — а = 0,1-^-0,14 мм; b = 0,2-^-0,028 мм;

для

приема

3—0,02 мм *).

При машинном решении выполнение этой операции можно осуществить с высокой степенью точности, так как положение центра окружности (прием 2) находится ана­ литически, как точка пересечения двух прямых, и может быть определено с точностью до 40-разрядного двоичного

1 Здесь а и Ь — малая и большая полуоси среднего эллипса оши­

бок.

144

числа. При машинном решении приемы 1 и 2 вообще не содержат ошибок (в пределах того порядка чисел, которые могут быть записаны в ячейке НФ). В общем случае при геометрографическом подходе определения точности того или иного способа машинного решения необходимо опре­ делять ошибку, которая появляется при выполнении каж­ дого стандартного оператора (§ 6, табл. 1). Этот путь не приведет к желаемому результату, так как конечной целью должно быть определение ошибки полученного ре­ зультата, являющегося следствием многократного выпол­ нения (в различной последовательности) стандартных опе­ раторов, т. е. суммарной ошибки.

Для определения величины суммарной ошибки при­ ходится пользоваться громоздкими формулами:

если ошибка зависит от двух операций — выражением

 

+ СО -\-со

 

 

 

 

P = g

\

J

е-2

^ Т - <а р '

R dPl

dp2;

(70)

 

—со —со

 

 

 

 

 

при ошибке, зависящей от трех

операций, — выра­

жением

 

 

 

 

 

 

 

р = ЬЫ*

(

[ f

e-vx'+itxu+vy+l&xz+ez

+ yyz)d x d y

fa

( 7 1 )

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

В § 9 на рис. 61 приведены результаты машинного ре­ шения задачи по определению линии сечения поверхности параболического цилиндра плоскостью. Выясним степень точности полученного решения путем сравнения его с ана­ литическим определением линии пересечения.

10 С. А. Фролов

145

Уравнение параболического цилиндра в координатной системе O'x'y'z' (рис. 98)

У' = 2рх'\

уравнение плоскости, перпендикулярной к оси г',

г' = 0.

Связь между координатными системами O'x'y'z' и Oxyz запишем с помощью следующих формул преобразования:

 

х' = 1хх + ГП\У + tixz + а;

 

у' = 12х + т2у + я 2 г + Ь\

 

г' = /3 x + т3у + п3 г + с,

где коэффициенты /, m, п выражают­

ся через углы Эйлера:

lx

= cos ijj cos ф — sin ^ cos 6 sin ф ;

m1 = sin i|) cos ф — cos xp cos 9 З 1 п ф ;

 

nx = sin 0 sin ф ;

/ 2

=—cos -ф sin ф —• sin ij; cos 9 cos ф;

m2 = —sin ф sin ф — cos ij? cos 9 cos ф;

n 2

= sin 9 cos ф ;

l3

= sin i|) sin 9;

m3 = —cos i|) sin 9; n 3 = cos 9;

a, b, с — координаты точки О' в системе Oxyz.

Тогда уравнение параболического цилиндра с произ­ вольным направлением прямой — образующей примет вид

l2x + т2у + n2z + Ь = 2р (lxx + tnxy +

+

пхг

+

а)2 ;

(72)

уравнение плоскости, перпендикулярной

к оси г',

13х + т3у

+

« 3 2

+ с = 0.

(73)

Значение координат точек линии сечения параболиче­ ского цилиндра плоскостью, перпендикулярной к обра­ зующей поверхности, определяется совместным решением системы уравнений (72) и (73).

Проекции линии сечения определяются следующими уравнениями:

146

проекция

на

плоскость

/7Х

(хОу)

 

 

 

 

 

г _

 

13х + тэу + с ,

 

 

 

 

 

 

 

 

" з

 

 

 

 

12х + /п2г/ — -£-2- (13х + т3 г/ + с) + * =

 

 

 

 

"3

 

 

 

 

 

 

= 2р Г/^ + т х у

(/З х +

т3 г/ +

с) +

а ! 2 ;

(74)

 

L

 

 

»з

 

 

 

J

 

проекция

на

плоскость Я 2 (хОг)

 

 

 

 

 

 

_

+ п3г + с _

 

 

 

12Х +

« 2 г

~

(^3*

+ П 3 2

+ С ) +

^

=

 

= 2p\l1x

(l3x +

п3 г +

с) +

щг +

й]2 .

(75)

Принимаем углы наклона проекций образующих (а сле­ довательно, и оси г') к осям проекции равными 45°, зна­ чения а = 15, a b = 5, коэффициенты при х, у, z в урав­ нении плоскости равными 1/25, а р = 0,5.

Тогда уравнения (72) и (73) примут вид

 

х —

z =

а — (у — z —

Ь)2;

(76)

 

 

х +

у + z = 25.

 

(77)

Соответственно, уравнения проекций

будут:

на

плоскости

П1

 

 

 

 

2х +

у — 40 = —(2у +

х —

30)3 ;

 

 

г =

25 — х у,

 

(78)

на

плоскости

Я 2

 

 

 

 

х —

z

15 = —(2г +

х

20)2 ;

 

 

г/ = 25 — х — z.

(79)

Результаты решений сведены в табл. 4.

В

табл. 4 (графа

4) приведены

значения координат

тех же точек, полученных при машинном решении (см. под­ черкнутые числа на рис. 61, § 9). Для сравнения (в графе 5) даны координаты тех же точек В', Б 2 , Вй, . . ., Б 6 , опре­ деленные по чертежу, выполненному вручную (см. рис. 99).

10*

147

 

 

 

 

Таблица 4

Точка _

Значение

Аналитиче-

Машинное

Ручное

 

на чертеже

координат

решение

решение

решение

 

X

733,00

733,187500

730,0

 

У

833,50

833,500000

835,0

А

X

816,50

816,562500

815,0

У

816,70

816,750000

818,0

 

X

833,20

833,250000

832,0

 

У

833,20

833,250000

833,0

 

X

783,20

783,250000

780,0

 

У

883,20

883,250000

885,0

R5

X

666,50

666,562500

665,0

В \

У

966,70

966,750000

965,0

Я6

X

591,20

591,218750

585,0

В1

У

102,00

101,725000

1020,0

Приведенные данные подтверждают предположение о том, что решение задачи на ЭЦВМ позволяет получить более высокую точность, чем при «ручном» решении.

Из табл. 4 видно, что результаты машинного решения почти не отличаются от аналитического. В данной задаче

148

определение координат каждой точки не зависит от зна­ чения координат, полученных для предшествующих то­ чек, поэтому ошибки не накапливаются.

Больший интерес представляет сопоставление результа­ тов решения задачи, у которой значение каждой после­ дующей величины зависит от результатов, полученных для предыдущей величины. Таким примером может служить задача по определению положения центра тяжести пла­

стинки

(см. § И, стр. 104). В этой задаче положение центра

тяжести

Z 2 полностью зависит

от Z x ; Z 3 от Z 2 ; Z4

от

Z 3

и т. д. Поэтому ошибка при

определении Zx влияет

на

точность определения Z 2 . В свою очередь, ошибка в опре­

делении Z 2

(совместно с ошибкой при подсчете Zr)

сказы­

вается на

точности нахождения Z 3 и т. д. Сравнение

ре­

зультатов машинного и аналитического решения этой за­ дачи дает хорошее совпадение. Ошибка в определении центра тяжести пластины (Z0 ) составила всего лишь 0,01 %.