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

книги из ГПНТБ / Кафаров, В. В. Принципы математического моделирования химико-технологических систем (введение в системотехнику химических производств) учеб. пособие

.pdf
Скачиваний:
86
Добавлен:
22.10.2023
Размер:
12.55 Mб
Скачать

Граф называется п о л н ы м, если

любая пара

его вершин

соединена ребром (рис. 1V-5).

т о п о л о г и ч е с к о й

Граф является не геометрической, а

ф и г у р о й. Последней называют такую фигуру,

определенные

свойства которой инвариантны при взаимнонепрерывном и взаимно­

однозначном

пространственном

преобразовании.

 

 

 

Существенные инвариантные свойства графа отражают только

число вершин,

число дуг (ребер)

и характер связей между верши­

 

 

 

 

 

 

нами. Так как граф яв­

 

 

л

 

 

ляется

топологической

 

 

 

 

фигурой, то один и тот же

 

1

т

 

граф может быть изобра­

I

5

жен различными спосо­

г

 

3

4

бами: вершины можно рас­

 

 

Рпс.

ІѴ-5. Полные графы.

 

полагать

в

произвольном

их ребра (дуги)

 

 

порядке,

а

соединяющие

проводить в виде прямых или ломаных линий. Вместе

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

ь

а

6

Рпс. ІѴ-6. Исходный граф (а) и изоморфный ему граф (б).

Примером неориентированных изоморфных графов могут слу­ жить графы, изображенные на рис. ІѴ-6, а, б. Обязательным усло­ вием изоморфности ориентированных графов является одинаковая ориентация всех дуг. Получение графа, изоморфного некоторому исходному графу, можно наглядно представить, изобразив этот исходный граф на упругой поверхности, например на листе резины. Какой бы деформации без разрушения не подвергалась поверхность листа резины, изображенный на ней граф не претерпит топологиче­ ских изменений. Каждый вновь образующийся граф при данной деформации листа резины будет изоморфен исходному графу, хотя геометрические фигуры, изображающие графы, при этом существен­ но отличаются друг от друга.

Граф G', содержащийся в графе G, будем называть п о д г р а ­ ф о м данного графа. Так, если граф G имеет вид, представленный на рис. ІѴ-7, а, то графы G', G', показанные сплошными линиями

118

на рис. ІѴ-7, б, в, будут подграфами графа G. Подграф G' содержит в общем случае лишь часть вершин и ребер (дуг), входящих в граф G, однако не может содержать вершины или ребра (дуги), не име­ ющиеся в G.

Рис. ІѴ-7. Исходный граф (а), подграфы (б, в) и их дополнения (г, д).

Граф, состоящий из вершин и ребер (дуг), которыми исходный

граф G отличается от подграфа

G\ будем называть д о п о л н е ­

н и е м D подграфа G'. Для

подграфов, изображенных на

 

а

Рис. ІѴ-8. К определению степеней вершин неориентированного (а) и ориентированного (б) графов.

рис. ІѴ-7, б, в, дополнения показаны пунктирными линиями и изобра­ жены отдельно на рис. ІѴ-7, г, д. Подграф G' и его дополнение D, взятые вместе, всегда образуют исходный граф G.

Число ребер, инцидентных некоторой вершине і, будем назы­

вать с т е п е н

ь ю

в е р ш и н ы

і и обозначать р (г). Так, степени

вершин графа

на рис. ІѴ-8, а соответственно равны:

p(a) =

p(d) = 3; р(Ь) = 5;

р (с) =4; р(«) = 1

119

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

входящих р'

(і) и выходящих р"

(і) дуг.

Для графа, показанного

на рис. ІѴ-8,

б, имеем:

А

 

 

р '(а )= р '(6 )= 1 ;

Р '(с)=3;

р''(а) =

р " Д )=2; р" (с) = 1

Очевидно,

что в ориентированном графе

2 р'(0=2>'(0

и

2 р(о= 2 р' w + S p' w

Для подсчета общего числа ребер е некоторого графа G, имеющего V вершин, удобно пользоваться следующим соотношением:

V

(іѵ д)

1

Так, для графа G, приведенного на рис. ІѴ-7, а, суммируем сте­ пени всех его вершин и делим сумму пополам, получая в итоге

а Ь

с d е {

_ 2 + 3 + 3 + 3 + 3 + 2

*--------------

2--------------

8

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

В случае ориентированного графа общее число дуг может быть подсчитано по одному из следующих очевидных выражений:

е = 2 Р' (О 1

И Л И

е = І > '( 0 1

Например, для графа, представленного на рис. ІѴ-4, а, легко получить:

 

 

a b

e

d

a b

e

d

 

 

 

е = 1 + 3 + 1 + 1==2 + 0 + 2 + 2 = 6

 

Д е р е в о м

с в я з н о г о

 

г р а ф а

 

G называется связный

граф Т, содержащий все вершины графа G, но не имеющий циклов.

Ребра

(дуги) графа, входящие

в

дерево,

называются в е т в я м и ,

а ребра (дуги), входящие в дополнение

дерева, х о р д а м и .

На

рис. ІѴ-9, а изображен некоторый граф G, а на рис. ІѴ-9, б, в пока­

заны

сплошными

линиями

два

возможных

дерева этого графа

Т х

и Т 2.

Дополнения деревьев обозначены на соответствующих рисун­

ках пунктирными

линиями.

 

 

 

 

 

120

Дерево, все ветви которого имеют общую вершину, называют иногда л а г р а н ж е в ы м д е р е в о м (рис. ІѴ-9, б).

Несвязный граф имеет несколько деревьев — по числу связных компонентов. Совокупность деревьев несвязного графа называют л е с о м F графа. На рис. ІѴ-10, б лес несвязного графа выделен сплошными линиями.

Рис. ІѴ-9. Исходный граф (а), два возможных дерева [одно — лагранжево дерево (б); другое — в] и фундаментальные циклы (г)

графа.

Дерево (лес) и его дополнение играют важную роль в практиче­ ских приложениях теории графов к исследованию ХТС. Рассмотрим несколько теорем, определяющих связь между числом вершин, ветвей и хорд графа.

Рис. ІѴ-10. Несвязный граф (а) и лес несвязного графа (б). К составлению матрицы смежности связного графа (в).

Теорема ІѴ-1. Любой связный граф имеет дерево.

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

подграф не окажется деревом.

Так, для превращения графа G (см. рис. ІѴ-9, а)

в дерево Т2 (см. рис. ІѴ-9, в)

последовательно удалялись ребра 4, 5,

6 и т. п.

Теорема ІѴ-2. Если дерево имеет ѵт вершин, то число его ветвей

ет опре­

деляется соотношением

= Vrp" 1

(1\ ,2)

 

Теорема, очевидно, справедлива, если дерево состоит из одной ветви и имеет две вершины. Добавление еще одной связной ветви увеличит ѵти ет на единицу,

т. е. не нарушит соотношения (IV,2). Это соотношение не нарушится, иока конфигурация с добавленными ветвями остается деревом. Например, любое

121

дерево Т графа G (см. рис.

ІѴ-9, а) содержит ѵт = 4 — 1

= 3 ветви. Из теоремы

ІѴ-2 вытекает ряд важных следствий.

 

Следствие 1. У графа,

содержащего циклы, число вершин ѵ и число ребер е

находится в соотношении

е +> ѵ — 1. Действительно,

добавление к дереву

ребра, образующего цикл, увеличивает е на единицу, но не изменяет ѵ. Следствие 2. Если лес состоит из к компонентов и имеет ѵр вершин, то

число ветвей определяется соотношением

 

eF = vF —K

аѵ ,з)

Действительно, для каждого связного компонента справедливо выражение (IV,2). Суммирование таких выражений для к компонентов леса с учетом оче­

видных

соотношений

ер =

и

vf =

2

ѴТ приводит к выражению (IV,3 ) •

Для графа, изображенного на рис.

ІѴ-10, а,

число ветвей е =

9 — 3 =

6 (вы

делены сплошными линиями).

входящих

в

дополнение дерева, составляет:

Следствие 3.

Число хорд,

 

 

 

egT = e

еТ = е

і-'у + І

 

 

( I V , 4)

Для

графа,

показанного на рис. ІѴ-9, а, число хорд es

= 6

— 4 + 1 =

= 3 (пунктирные линии).

входящих

в дополнение леса, равно:

 

Следствие 4.

Число хорд,

 

 

 

 

eSF = e —eF = e —vf + k

 

 

(1V,5)

Для

графа, представленного на рис.

ІѴ-10, а, число хорд es

= 8

— 9 +

+ 3 = 2

(пунктирные линии).

 

 

 

 

 

 

 

Число хорд называют также д и п л о м а т и ч е с к и м

ч и с-

л о м,

ч и с л о м

Б е т т и

или

д е ф и ц и т о м

г р а ф а .

Это число указывает, сколько ребер графа должно быть удалено, чтобы превратить этот граф в дерево (лес).

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

нулю.

Введем понятие ранга графа, смысл которого выяснится позднее. Ранг R графа G описывается соотношением

Д = (IV,6)

где к —- число связных компонентов несвязного графа. Следова­ тельно, ранг связного графа на единицу меньше числа его вершин.

Информация, содержащаяся в каком-либо графе, может быть представлена в иной форме — в форме м а т р и ц ы . Далее мы будем использовать обе формы представления информации — графическую (в виде графа) и численную (в виде матрицы).

2. ПРЕДСТАВЛЕНИЕ ГРАФОВ С ПОМОЩЬЮ МАТРИЦ

Информация, содержащаяся в некотором графе, может быть представлена в алгебраическом виде посредством матриц. Эта связь графа и матрицы имеет чрезвычайно важное значение при практиче-

122

сном приложении топологических методов к математическому .опи­ санию ХТС, так как позволяет перевести структурные особенности системы на язык чисел, фигурирующих в математических уравне­

ниях.

М а т р и ц е й с м е ж н о с т и , соответствующей некоторому графу G = (X , Y), который состоит из п вершин X t (і = 1, 2,. . ., п), называется матрица [Н] порядка (п X п) с элементами

(

0,

если вершина ж,- не связана дугой с вершиной xj

(I ’ ’7)

hij = (

1,

 

 

 

 

I

если вершина ж/ связана дугой с вершиной ж/

 

Например, матрица [Н] некоторого

графа (см. рис.

ІѴ-10, в)

имеет следующий вид:

 

 

 

 

 

 

 

“0

0

0

1

1'

 

 

 

0

1

0

0

1

 

 

 

1

0

0

0

0

 

 

 

0

1

0

0

1

 

 

 

_0

0

0

1

0.

 

Заметим, что диагональный элемент /г22 = 1, что соответствует

наличию петли при вершине 2.

 

 

М а т р и ц е й

и н ц и д е н ц и й , которая соответствует неко­

торому графу G

=

(к,

q), состоящему из п вершин к(- =

1,

2, . . .,п)

и из т дуг

qj

(j

= 1,

2, . . ., т), называется матрица

[S]

порядка

(п X т) с

элементами

 

 

 

 

 

—1, если дуга q,- выходит из вершины щ

 

 

 

sч ~

+1, если дуга ?/ входит в вершину к(-

 

( I V , 8)

 

 

 

О,

если дуга д/ не инцидентна вершине к,-

 

Матрицу [S] называют также с т р у к т у р н о й м а т р и ц е й графа.

В качестве примера составим матрицы инциденций для графов,

изображенных на рис. ІѴ-4,

а, б.

ІѴ-4, а,

имеем:

Для графа, приведенного

на

рис.

і

І

1

2

3

4

5

6

 

 

 

 

 

 

 

а

"-1

0

0

1

0

-1

Ъ

 

1

1

0

0

1

0

с

 

0

—1

—1

0

0

- 1

d

 

0

0

1

—1

—1

0

123

Д л я н есв язн ого граф а, п редставлен н ого на ри с. ІѴ -4, б, получим :

 

І 1

2

3

і

s

6

7

а

і

 

 

 

 

 

 

- - 1

1

0

0

0

0

<r

Ъ

1

—1

0

0

0

0

0

с

0

0

- 1

0

0

0

1

d

0

0

1

—1

—1

0

0

е

0

0

0

1

0

1

0

/ _ 0

0

0

0

1

- 1 1_

или

 

 

 

О ■

 

 

 

 

 

[ S 2] =

 

 

 

 

 

 

. о

S i .

 

 

 

 

 

 

 

 

 

Здесь S' и S' — подматрицы инциденций, соответствующие связным компонентам несвязного графа.

Свойства матриц инциденций отражают топологические особен­ ности соответствующих графов и могут быть сформулированы в виде трех теорем.

Теорема ІѴ-3. Любой определитель, содержащийся в матрице инциденций [S], равен нулю, + 1 или —1.

Если каждый столбец определителя Д, содержащегося в [S], имеет два ненулевых элемента (+ 1 и —1), то Д = О, так как сумма всех строк дает строку, все элементы которой равны нулю. Если хотя бы один столбец А имеет только один ненулевой элемент, то при разложении Д по элементам этого столбца можно в результате прийти к определителю первого порядка. Таким образом, Д равен 0, + 1 или —1.

Теорема ІѴ-4. Определитель любой подматрицы матрицы инциденций, отвечающей циклу, равен нулю.

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

Теорема ІѴ-5. Определитель (то—1) порядка подматрицы (то—1) ранга матрицы инциденций (то ■— число строк этой матрицы), отличный от нуля, отвечает дереву исходного связного графа.

Обратная теорема: определитель, соответствующий дереву графа, отличен от нуля.

Число ветвей дерева, как было отмечено ранее, на единицу мень­ ше числа вершин, следовательно, определитель дерева должен иметь порядок — 1). Если определитель порядка — 1) равен нулю, то, согласно теореме ІѴ-4, он отвечает некоторому циклу в графе и, значит, не соответствует дереву. Наоборот, если такой определитель не равен нулю, отвечающий ему граф, который имеет пг — 1 дуг, не содержит циклов и, следовательно, является деревом.

Теорема ІѴ-5 может быть распространена и на несвязные, графы. В этом случае определитель, соответствующий лесу графа, имеет порядок — к) и отличен от нуля, и обратно — определитель по­ рядка — к), отличный от нуля, отвечает лесу несвязного графа.

Таким образом, р а н г

м а т р и ц ы и н ц и д е н ц и й всегда

равен числу ветвей дерева

(леса) графа.

124

Применение приведенных теорем можно проиллюстрировать на примере графов, изображенных на рис. ІѴ-4, а, б. Так, легко убе­ диться, что определитель третьего порядка, записанный для эле­ ментов 1, 2, б (см. рис. ІѴ-4, а), равен нулю:

1 2 6

а1 0 г

Ді

1

1

0 = 0

с

о

1

1_

асоответствующая ему конфигурация в графе является циклом. Любой определитель третьего порядка для групп элементов 1, 2,

3; 5, 4, 6; 1, 4, 3 и т. д. отличен от нуля:

II 1С- <

или

Д3 =

 

1

2

3

а

_-1

0

сг

о»

1

1

0

с

- 0

—1

1_

 

1

2

3

ъ "1

1

0"

с

0

-1 - -1 =

d _0

0

1_

Указанные конфигурации элементов представляют собой деревья графа. Подобные результаты легко получить и для несвязного графа, показанного на рис. ІѴ-4, б.

Из приведенных выше теорем следует важное утверждение о ра­ венстве ранга графа, определяемого выражением (IV,6), и ранга

матрицы пнциденций, т. е.

 

rs=Ä

(IV,9)

Например, ранг графа, представленного на рис. ІѴ-4, а, равен:

R = v — 1 = 4—1 = 3

Матрица инциденций этого графа путем элементарных преобра­ зований приводится к следующему виду:

1 0 0 1 0 -1

[Si] =

0 1 0 1 1 —1

0 0 —1 1 1 0

1 о

0

0

0

0

0

т. е. ранг [SJ также rs = 3.

125

Д л я граф а, и зоб р а ж ен н о го на ри с. ІѴ -4, б, имеем:

R = v —k = 6 —2 = 4

Матрица инциденций [S2] этого графа приводится к виду

"1

1

0

0

0

0

0'

0

0

0

0

0

0

0

0

0

1

0

1

0

1

0

0

0

1

1

0

1

0

0

0

0

0

- 1

1

_0

0

0

0

0

0

0.

т. е. ее ранг rs = 4.

Равенство рангов графа и соответствующей матрицы инциденций позволяет в ряде случаев значительно упростить вычисления.

Кроме матрицы инциденций [S],

для каждого графа может быть

записана другая матрица, называемая матрицей циклов

[М], м а т ­

р и ц е й с о в п а д е н и й

или

с о е д и н е н и й .

 

Общий элемент

матрицы [М] будет:

 

 

0,

еслп ребро

не входит в цикл р/

 

т,7 =

1,

если ребро входит в цикл п его ориентация

совпадает с ориентацией цикла

(IV, 10)

 

—1, если ребро входит в цикл и его ориента­

 

ция противоположна ориентации цикла

 

Таким образом,

для составления матрицы [М] в графе выделяют

все возможные

циклы, придавая

им определенную

ориентацию.

Так, в графе, показанном на рис. ІѴ-4, а, содержатся следующие циклы: 1 - 1 , 5 , 4; II - 1, 2, 3, 4; III - 5, 2, 3; IV — 1, 2, 6; V — 4, 3, 6; VI — 4, 5, 2, 6; VII — 1, 5, 3, 6, а в графе, представленном на рис. ІѴ-4, б, содержатся циклы: I — 1, 2; II — 3, 5, 7; III — 4, 43, 5; IV — 3, 4, 6, 7 (принята ориентация циклов по часовой стрелке).

Матрицы циклов для графов, изображенных на рис. ІѴ-4, а, б, имеют следующий вид:

 

I

II

III

IV V VI VII

1

- 1

1

0

1

0

0

г

2

0

-1 —1 - -1 . 0

- 1

0

3

0

1

1

0

—1

0

- 1

4

1

1

0

0

—1 - 1

0

5

—1

0

1

0

0

1

- 1

6

_ 0

0

0 - -1

—1 - 1

-1 _

126

 

I

II

III

IV

 

1 "

1

0

0

0

-

2

1

0

0

0

 

3

0

1

0

1

 

[>Ь] = 4

0

0

1

1

 

5

0

1

—1

0

 

6

0

0

—1

—1

 

7

0

1

0

1

_

пли

- М( 0 '

[М,] =

0 щ .

Здесь М' и М' — подматрицы, соответствующие связным компо­ нентам несвязного графа.

Теорема ІѴ-6. Ганг матрицы циклов равен числу хорд соответствующего графа, т. е.

rM = e —v+ k

(IV,11}

или

Доказательство теоремы можно выполнить в два приема. Как уже отме­ чалось, дерево (лес) не имеет циклов. Если к дереву добавить одну хорду, то* образуется один цикл и соответственно появится столбец в матрице циклов. Убрав первую хорду и добавив новую, получим еще один новый цикл — новый столбец в матрице циклов. Таким образом, каждой хорде будет отвечать неза­ висимый столбец в матрице циклов, а общее число независимых столбцов (ранг матрицы) будет по крайней мере не меньше числа хорд. Можно доказать, что* ранг матрицы не может быть больше числа хорд, т. е. всегда справедливо соот­ ношение (IV,11).

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

Для формального дерева «Г3» в графе (см. рис. ІѴ-9, г) можно выделить следующие фундаментальные циклы (рг — образован хор­ дой q2 и ветвями qr, g4 и т. д.):

{32.91,94}; Р2={Дз> 94. 95};

Гз = (?б> Чи 44, Чъ}

Матрицу циклов, составленную только из векторов-циклов, соот­

ветствующих фундаментальным циклам, называют ц и к

л о м

э т и ­

ч е с к о й м а т р и ц е й

[С]. Порядок матрицы

[С]

равен X

X — р)], а ее ранг

тс = гм = е — р = es .

При

записи

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

12Т

Соседние файлы в папке книги из ГПНТБ