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

книги из ГПНТБ / Голенко Д.И. Статистические модели в управлении производством

.pdf
Скачиваний:
20
Добавлен:
25.10.2023
Размер:
15.19 Mб
Скачать

Рассмотрим метод (Б) определения оптимального со­ вокупного варианта программы, который целесообразно применять в тех случаях, когда:

а) смешанная сетевая модель содержит очень боль­

шое количество альтернативных вершин

(типов а и а ) ,

л построение дерева исходов D(A, V) является

весьма

трудоемкой задачей;

 

 

 

б) неоднородное дерево исходов D(A,

V)

невелико по

объему, но структурные взаимосвязи

альтернативных

ветвлений различной логической природы

и а) на­

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

Идея предлагаемого метода состоит в использовании

отдельных вычислительных

блоков альфа-алгоритма

[5.24] и

имитации

структуры

развертывания

многовари­

антной

программы

во времени на основе статистическо­

го моделирования

[5.7, 5.38]

альтернативной

сети. При

этом ставится задача последовательного просмотра сете­ вой модели G(Y, U) (либо дерева исходов D(А, V)), на­ чиная от начальной вершины, и отсечения в каждом де­ терминированном ветвлении а «неоптимальных» направ­

лений. Для иллюстрации

изложения

используем

при­

мер

на рис. 5.3.5, имея в виду имитацию принятия

реше­

ний

в ос-вершинах дерева

исходов D(A,

V). Распростра­

нение предлагаемого метода на случай непосредственного

выделения в

сети G(Y, U)

стохастической

подсети

G<r)(Yir\ £ЛГ)),

соответствующей

оптимальному

совокуп­

ному варианту, на наш взгляд, трудностей не вызывает.

Предположим,

что, исследуя

дерево

исходов D(A,

V),

представленное

на рис. 5.3.5,

мы двигались по дуге

(«ь

а 2 )

и достигли

альтернативной вершины первого поколе­

ния,

которая, допустим,

принадлежит

к типу а. В связи

с тем, что аг

отражает

ситуацию

детерминированного

(управляемого) ветвления вариантов, необходимо опре­ делить оптимальное направление (а 2 , «*) и отсечь неоп­ тимальные наборы путей, включающие вершины

ЛЛ

Га 2 / Г а* •

Положим, a,i = a%<=A и зададимся числом N реализа-

дий обхода графа D(A, V) по каждому 1-му направлению; / ^ { 1 , 2, ... , ПІ), щ= | Г а г | . Выбираем одно из щ направ­ лений. Двигаясь вдоль выбранного направления и одно­ временно определяя характеристики {X} пройденных пу­

тей, будем

использовать

следующее правило

имитации

встречающихся

ветвлений:

 

 

 

 

J_)

достигнув

а-вершины

стохастического типа

( а ч є

е Л ) ,

осуществляем

розыгрыш

случайной величины gg ,

равномерно

распределенной

в

интервале [0,1],

и в

соот­

ветствии с

вероятностями

{р}

исходящих дуг (aq,

а) и-

значением

| д определяем

направление дальнейшего пе­

рехода;

 

 

 

 

 

 

 

 

2) направление,_исходящее из вершины детерминиро­

ванного типа ( а 9 є Л ) ,

разыгрываем также с помощью | 3

при допущении, что всем детерминированным

альтерна­

тивам

приписана в е р о я т н о с т ь ( я д = | Г а д | число

альтернатив).

После достижения финального исхода и расчета ха­ рактеристик {X} комплекса операций, связывающего оцениваемую альтернативную вершину аг- с данным ис­

ходом, определяем значение критерия качества

где

п—-номер реализаций; / — номер оцениваемого

направ­

ления (здесь можно использовать любой из рассмотрен­ ных выше критериев Fj оценки полных вариантов). Вы­ полнив N реализаций по данному 1-му направлению, рас­ считаем итоговое значение функции качества 1-го направ­ ления:

(п)

1

N

(га)

 

Ft = M[Fi ] =

2

Fi

(5.5.29)

 

N

n = i

 

 

Подобная процедура имитации выполняется для всех щ направлений, «сходящих из а ^ є А. В качестве опти­ мального выбирается /*-е направление с первой ветвью (аг-, аг*), для которого:

У> = externum [ Л ] .

(5.5.30)

1= T^ni

 

В дальнейшем исключаем из рассмотрения все «неоптямальные» направления (с номерами Іфі*) и перехо­ дим к а-вершине следующего поколения (аг* ) . Если ока-

зывается, что ш * & 4 , то полагаем ai = ai* и повторяем описанную процедуру розыгрыша направлений, исходя­

щих из детерминированного ветвления а*. В случае,

если

С ' с е Л , необходимо просмотреть все стохастические

на­

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

нальными

исходами. Фиксируем количество

указанных

 

 

Л

_

а-вершин

\{ak}\=K,

для которых аьєГаг* ;

a,k^A\A';

Г _ 1 а ь є Л . Полагая

далее аг- = а& последовательно для де­

терминированных

ветвлений аь а.2,... а.к, повторяем про­

цедуру розыгрыша оптимальных направлений. Оптималь­ ный совокупный вариант считается найденным, когда все определенные направления вида (ai, ai* ) завершаются в финальных исходах, т. е. а;* є А'. Рассмотренный метод обеспечивает довольно быстрое нахождение оптимально­ го совокупного варианта за счет последовательного суже­ ния зоны его поиска.

Вернемся к примеру на рис. 5.3.5. Предположим, что на основе применения описанного метода мы получили следующие результаты:

в вершине аг выбрано направление

2, а 3 )

и отсе-

' чены

все пути, минующие а 3 (т. е. определены

«неопти­

мальные» исходы аіб—ai9 и аю, а п ) ;

_

 

 

— в вершине а7 выбрано направление

(ал, ctis).

Таким образом, найденный оптимальный вариант со-

 

 

*

*

*

ответствует стохастическому дереву исходов DO (АО, V(r>) для г = 2 на рис. 5.3.6.

§ 5. 6. Принятие решений в точках ветвления

смешанных альтернативных сетей

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

шений, отображаемых в сети событиями типа а. Управ­ ление программой осуществляется аппаратом управле­ ния путем корректировки полной смешанной сетевой мо­ дели G(Y, U) с учетом текущего состояния программы и принимаемых решений. Необходимость корректировки сети может возникнуть вследствие: изменения оценок (продолжительности, стоимости и др.) элементарных опе­ раций; изменения вероятности и сравнительной эффек­ тивности конкурирующих частичных вариантов; отказа от запроектированных ранее направлений (вариантов) программы, ставших маловероятными или неэффектив­ ными; появления новых вариантов достижения конечных целей программы.

Кроме использования указанных «пассивных» управ­ ляющих воздействий, вызывающих составление нового плана реализации программы, возникает возможность и необходимость выработки на основе смешанной модели «активных» управляющих воздействий. Последними слу­ жат процедуры принятия решений в альтернативных си­ туациях программы. Естественно, что для выполнения этих задач система управления должна иметь в своем со­ ставе постоянную службу СПУ [5.24], в обязанности ко­ торой входят периодическая корректировка смешанной альтернативной сетевой модели и выполнение необходи­ мых модельных расчетов путем реализации комплекса вычислительных алгоритмов на ЭВМ. В результате служ­ ба СПУ вырабатывает рекомендации по принятию реше­ ний в процессе оперативного управления программой. Рассмотрим методы использования смешанной модели для определения оптимальных решений в возникающих альтернативных ситуациях в процессе реального осуще­ ствления программы.

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

следующем. Исследуется альтернативная сеть, первое ветвящее событие которой отображается наступающей альтернативной ситуацией типа а. Первое ветвление

(обозначим его через аг) рассматривается в

соответствии

с его логической природой как управляемая

ситуация вы­

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

Задача состоит в исследовании стохастических под­ сетей с корневой вершиной в наступающей ситуации щ и первыми сетевыми фрагментами G,j, соответствующими детерминированным частичным вариантам, возникаю­ щим в «г. Таким образом, количество детерминирован­ ных направлений, возникающих в а*, определяет количе­ ство подлежащих оценке и сравнению стохастических сетей (tii= | Г а г | ) . Сформулированная задача выбора оптимального направления в ситуации аг- решается путем имитации на ЭВМ вероятностных процессов, отобража­ емых стохастическими подсетями. С этой целью для каж­ дого из I направлений (/=1,2, ... , пг) и соответствующих подсетей выполняется N реализаций, в результате чего получаются оценки математических ожиданий функции

качества {Fi},

1=1, щ.

Решение в наступающей

ситуации,

моделируемой в сети

событием at-, принимается путем

сравнения Fi, 1= 1, пг- и выбора номера I* такого

направле­

ния

реализации

программы,

для

которого

Fi* =extremum [Fi].

 

 

 

1=1,

ПІ

 

 

 

 

П.

При

приближении моделируемого

смешанной

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

нованный на теории

статистических

решений [5.27,

5.28].

В терминологии этой теории

под «статистическими

реше­

ниями» понимаются

такие

решения,

которые принима-

ются на основе изучения соответствующей статистичес­ кой обработки совокупности случайных величин и слу­ чайных процессов. По тем же соображениям, что и в слу­ чае анализа ситуации а, будем рассматривать последую­ щие альтернативные ситуации программы как стохасти­ ческие. Задача состоит в принятии решения в точках ветв­ ления стохастической сети, первое ветвление которой со­ ответствует исследуемой ситуации а*. Используя альфаалгоритм [5.24] или какую-либо из его модификаций [5.12—5.14], построим для исследуемой сети дерево исхо­ дов с разъединительными путями D (Л, V). Рассмотрим предлагаемый метод (см. также [5.39]). В данном случае, говоря о программе создания сложного комплекса, мы будем иметь в виду ту ее часть, которую еще предстоит выполнить. Таким образом, корневой вершиной дерева исходов/)(Л, V) служит момент моделирования програм­ мы, определяющий текущее состояние, в котором нахо­ дится процесс^ создания сложного комплекса. Пусть де­ рево исходов Ь(А, V) характеризуется выборкой апосте­ риорных данных D, содержащих сведения относительно возможной продолжительности программы, ее трудоем­ кости, стоимости, интенсивности потребления ресурсов и т. д. В общем случае данные, характеризующие моде­ лируемый процесс, могут быть представлены вектором

D = {DU D2,...,Dn}.

Обозначим

множество

вариантов

осуществления

программы через

5 = { 5 Ь

S2,

...,Sn).

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

процесс), либо 5 задается как функция одного

или не­

скольких

параметров

V={vu

v2,...,vi}.

Плотность

веро­

ятностей

а (5) может

быть

получена

на основании

соот­

ношений, характеризующих

каждую

а-вершину

D(A,

V),

и введенного критерия. В зависимости от характера

зада­

чи исследования а (5)

может быть

дискретной

и непре­

рывной функцией. Если 5 рассматривается как случай­ ный процесс, характеризующийся параметрами V, экви­ валентом a(S) может служить плотность вероятностей

a(V). Решение

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

соответствии

с

некоторым

правилом. Каждое

решение

(обозначим

его

через уІ)

можно рассматривать

как

эле­

мент множества

возможных решений у= {уь у 2 , . . . } .

Пра­

вило в этом случае является некоторой функцией у име­

ющихся данных D и определяется решающей

функцией

6(y/D). Таким

образом, компоненты

D образуют

случай­

ную выборку

апостериорных данных

из пространства Г,

на основании

которых принимается

решение

у.

Наряду

с введенными составляющими рассматриваемого процес­ са (S, D, V, у, б) необходимо учесть влияние случайных воздействий, действующих на систему управления созда­ нием сложного комплекса. Через W(N) обозначим плот­ ность вероятностей случайных воздействий, выражаю­ щихся в первую очередь в случайном колебании налич­ ных ресурсов в период дальнейшего выполнения програм­ мы и непредвиденных изменениях условий ее осуществле­ ния. Таким образом, результаты (D*) функционирования системы формируются в зависимости от характера про­ странства возможных способов осуществления програм­ мы (5*) и вида пространства случайных воздействий (Л'"), а также от распределений вероятностей вариантов про­ цесса а ( 5 ) и распределения вероятностей случайных

воздействий W(N), т. е. D*=SN.

Всилу существования самых различных продолжений

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

является

функцией реализации случайного процесса,

С (5, у)

является случайной функцией. На этом основа­

нии оценку риска или потерь следует искать как вероят­ ностную характеристику, например как математическое

ожидание. Выражения условной и средней оценки потерь при выбранном правиле в пространстве решений Д име­ ют вид:

г (S,6) = IF„ (Я/5) dD $ С(S,y) б(y/D)dy,

(5.6.1)

Т

д

R(o,b)=M{r(S,6)}=

lr(S,t>)o(S)dS

или

"

t (5.6.2)

R(a,6) =

lo(S)dS-§Fn(D/S)dD'$

C(S,y)6(y,/D)dy,

(5.6.3)

 

а

г

д

 

 

где Fn(D/S)

—функция распределения

данных

для из­

вестного

процесса S; C(S,

у) —функция

потерь

в стои­

мостном выражении, часто являющаяся функцией пара­

метров

V.

 

 

Как

следует

из соотношений (5.6.1) — (5.6.3),

сред­

ний риск зависит от принятого решающего правила

б и,

следовательно,

может служить для сравнения различных

правил принятия решений. Решающее правило, для кото­ рого средний риск (5. 6.3) оказывается минимальным, называется байесовым решением относительно априор­ ного распределения o(S). Используя соотношения (5.6.2)

и

(5.6.3)

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

Fn(D/S)

и о(S), можно определить

зависимость средне­

го

риска

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

параметра ст. Так, в

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

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

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

Щ

Ю

15

2°(6Cj-aLj)

 

 

Рис. 5. 6. I. Зависимость средного риска от размаха распределения.

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

не значатся в дереве исходов D(A,

V), но как

бы незри­

мо присутствуют в нем.

 

 

Рассмотрим дерево

исходов, множество а-вершин ко­

торого задано, причем

для каждой

а-вершины

известны

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

мер, а-вершину, характеризующуюся

двумя

исходами

( | Г « | = 2 )

с вероятностями

рх и

1рх.

Пусть

исход, ве­

роятность

которого равна

1ри

носит,

в свою очередь,

альтернативный характер и может разветвляться на два исхода, один из которых также ветвится, и т. д. Можно принять схему ветвления, показанную в табл. 5.6.1. При этом правило ветвления основывается на сохранении в

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

а-верши-

на удовлетворяет полной группе исходов.

 

 

 

 

 

 

 

 

Т а б л и ц а

5.6.1.

 

 

Обобщенный

 

 

 

 

 

Априорная вероятность

показатель,

Параметр

Примечание

возможного исхода

характеризую­

ветвления

 

 

 

 

 

 

щий исход

 

 

 

 

 

 

Pi

*1

Нет

Второй

исход

(1—р.) (I—ря )

 

 

 

 

 

 

ветвится

 

( 1

— Р з )

 

( 1 - й )

Третий,

вновь

 

 

 

 

образовавший­

 

 

 

 

ся

исход,

так­

 

 

 

 

же

ветвится

( 1 — Р і ) № ( 1 — Pi)

•ч

( 1 - А )

Четвертый,

 

 

 

 

вновь

образо­

 

 

 

 

вавшийся

ис­

 

 

 

 

ход

также вет­

 

 

 

 

вится

 

 

(1—рі)Р2Рз-Рп-і(\—Рп)

f/i + l

( 1 - А , )

И т.

д.

 

(1—Pi)

РїРі...Рп

Рп

 

 

 

 

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

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