книги из ГПНТБ / Фролов, С. А. Кибернетика и инженерная графика
.pdfмых 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 %.