д) в разделе Distance measure выберите Squared Euclidean distances (квадрат евклидовой метрики) и нажмите ОК.
В появившемся окне нажмите на кнопку Vertical icicle plot. На экране появится дендрограмма (рис. 18), показывающая объединение объектов, расстояние между которыми является наименьшим, в кластеры.
Рис. 18. Дендрограмма при методе одиночной связи
На вертикальной оси дендрограммы откладываются расстояния между объектами и между объектами и кластерами. Так, расстояние между объектами С, и С2 равно 5 (см. матрицу расстояний D2 , в примере 1). Эти объекты объединяются в
один кластер на первом шаге.
Расстояние между этим кластером и объектом С3 равно 8 (см. матрицу расстояний D2 ). Объект С3 объединяется с
кластером (С1, С2) на втором шаге. Наконец, расстояние между объектом С4 и кластером (С1, С2, С3) равно 13 (см. матрицу расстояний D3 ).
Таким образом, горизонтальные отрезки дендрограммы проводятся на уровнях, соответствующих пороговым значениям расстояний, выбираемым для данного шага кластеризации.
Кластеризация методом одиночной связи (ближайшего соседа) приводит к образованию одного кластера (пороговое
29
расстояние равно 13).
Далее последовательно нажмите Continue... и Cancel и в окне установок процедуры в разделе Amalgamation... выберите Complete Linkage. После выполнения процедуры появится следующая дендрограмма (рис. 19).
Рис. 19. Дендрограмма при методе полной связи
Вторая дендрогамма показывает, что кластеризация методом полной связи (дальнего соседа) при таком же пороговом расстоянии равным 13, приводит к образованию двух кластеров. Сравните полученные дендрограммы с результатами решения примера 1.
Задания для самостоятельного решения
1.Используя иерархические алгоритмы, проведите кластеризацию данных из примера 2. Сравните результаты с результатами в примере 2.
2.На предприятии существует 16 научнопроизводственных отделов, занятых выпуском различной продукции, работ, услуг. Поскольку виды деятельности, количество работающих, рентабельность отделов, существенно различаются между собой, было решено сгруппировать отделы
внесколько однородных групп, а затем для каждой группы разработать свою систему премирования.
30
После тщательного анализа выбрали четыре признака, с помощью которых описывались важные (для указанной цели)
параметры каждого отдела: X1 |
– стоимость активной части ос- |
||||||||||
новных |
производственных |
фондов, тыс. руб.; X 2 |
– средне- |
||||||||
месячный объем работ отдела, тыс. руб.; X 3 |
– удельный вес |
||||||||||
работ/услуг отдела по внутрифирменной кооперации, %; |
X 4 – |
||||||||||
среднемесячная прибыль отдела, тыс. руб. |
|
|
|
||||||||
Исходные данные по отделам приведены ниже (табл. 9). |
|||||||||||
|
|
|
|
|
|
|
|
|
|
Таблица 9 |
|
|
|
|
Исходные данные по отделам |
|
|
|
|||||
№ |
|
Значения признаков |
|
№ |
Значения признаков |
||||||
отдела |
|
|
|
|
|
|
отдела |
|
|
|
|
|
X1 |
X 2 |
X 3 |
X 4 |
|
X1 |
X 2 |
X 3 |
X 4 |
||
|
|
|
|
||||||||
1 |
|
699 |
190 |
53 |
11 |
|
9 |
293 |
391 |
16 |
66 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
532 |
211 |
19 |
42 |
|
10 |
300 |
396 |
29 |
87 |
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
650 |
152 |
46 |
14 |
|
11 |
73 |
160 |
0 |
22 |
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
768 |
216 |
67 |
17 |
|
12 |
862 |
199 |
51 |
22 |
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
67 |
106 |
0 |
32 |
|
13 |
112 |
136 |
0 |
29 |
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
322 |
397 |
26 |
52 |
|
14 |
289 |
388 |
31 |
74 |
|
|
|
|
|
|
|
|
|
|
|
|
7 |
|
736 |
180 |
49 |
18 |
|
15 |
512 |
195 |
6 |
58 |
|
|
|
|
|
|
|
|
|
|
|
|
8 |
|
501 |
239 |
11 |
60 |
|
16 |
490 |
201 |
9 |
65 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Проведите |
кластеризацию |
отделов |
используя |
|||||||
иерархические алгоритмы (Joining): |
|
|
|
|
||||||||
|
|
а) |
|
используя исходные данные; |
|
|
||||||
|
|
б) |
|
используя стандартизованные данные, т. е. данные, |
||||||||
преобразованные по формуле |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X ij |
X j |
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
|
S j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где X ij |
i |
е значение j го признака, i |
1,2,...,16; j |
1,2,3,4; |
||||||||
|
|
|
1 |
16 |
оценка среднего для |
j го признака; |
||||||
|
|
|
||||||||||
|
X j |
|
|
|
xij |
|||||||
|
|
|
|
|||||||||
|
|
|
16 i 1 |
|
|
|
|
|
|
|||
31
|
1 |
16 |
|
|
|
2 |
оценка |
среднего квадратического |
||
S j |
xij |
x j |
||||||||
|
|
|
|
|||||||
15 i 1 |
|
|||||||||
отклонения для |
j |
го признака. |
|
|||||||
|
Процедуру стандартизации |
данных можно выполнить |
||||||||
непосредственно в таблице, используя следующую последовательность действий: курсор на имени переменной 
нажать правую кнопку мыши |
в выпадающем меню выбрать |
||
File/Standardize Block |
Standardize Columns |
ОК. |
|
Сравните результаты |
кластеризации. |
По результатам |
|
кластеризации определите число кластеров и их состав. Найдите статистические характеристики каждого кластера.
Проведите кластеризацию используя метод К-средних (число кластеров задайте равным 4). Сравните результаты (составы кластеров).
3.4. Общие положения разомкнутых сетей массового обслуживания
1. Постановка задачи по разомкнутой сети Джексона
Первой работой, посвященной разомкнутым сетям массового обслуживания (РСеМО), является работа Дж. Р. Джексона, опубликованная в 1957 г. С тех пор произошло существенное расширение спектра изучаемых СеМО как достаточно адекватных моделей реальных сложных информационных, производственных, вычислительных, связевых и т.п. систем и сетей, используемых в самых разнообразных областях человеческой деятельности. При этом учитывались многие особенности функционирования реальных систем и сетей. Большинство результатов получено для стационарного режима.
Для теории и практики непосредственный интерес представляет получение результатов для нестационарного режима. Для этого режима точные результаты по РСеМО получены только для частного случая РСеМО Джексона. Для лучшего понимания последующих результатов для нестационарного режима марковских СеМО и их места в общей теории СеМО целесообразно привести с
32
доказательствами точные результаты для РСеМО Джексона по [3]. Сначала приведем постановку задачи.
Постановка задачи. Сеть состоит из т перенумерованных узлов. На вход сети поступает пуассоновский поток однородных требований интенсивности . Поступившее требование с вероятностью P0i поступает на обслуживание в i -
m
ый узел, P0i 1. Длительность обслуживания в i -ом узле
i 1
является случайной величиной, имеющей экспоненциальное распределение с пара-метром, зависящим от числа требований в узле. Более точно, если в момент t в i -ом узле имеется l требований, то в интервале t,t dt , где dt мало, закончится обслуживание одного требования с вероятностью
i l dt o dt ,0
i l
,l 1,2,...,
или закончится время пребывания одного требования в узле с вероятностью
|
|
vi |
l dt |
o dt ,0 |
vi l |
,l |
1,2,... |
|
||
|
Далее и в том, и в другом случае это требование |
|||||||||
мгновенно поступает на обслуживание в |
j -ый узел с |
|||||||||
вероятностью |
Pij |
и с |
вероятностью |
Pi0 |
покидает |
сеть, |
||||
m |
|
|
|
|
|
|
|
|
|
|
Pij |
Pio 1,i |
1,...,m . Матрица |
Pij |
m m |
неразложима. |
|
||||
j 1 |
|
|
|
|
|
|
|
|
|
|
|
Вводится случайный процесс |
|
|
|
|
|||||
|
|
|
|
v t |
v1 t ,...,vm t , |
|
|
|||
где vi |
t — число требований в i -ом узле в момент t, i |
1,...,m . |
||||||||
В силу предпосылок постановки задачи он является марковским случайным процессом. Обозначим
P k1,...,km ,t |
P v1 t k1,...,vm t |
km , |
Требуется определить |
нестационарное |
распределение |
вероятностей состояний этой сети P k1,...,km ,t .
В практических задачах используются, как правило, не нестационарные распределения вероятностей состояний сети, а
33