Комплексный подход к оценке имущественных отношений между экономическими субъектами Российской Федерации (микро- и мезоэкономические аспекты)
С.В. Лушин
аспирант дневного отделения Центрального экономико-математического института Российской академии наук (ЦЭМИ РАН)
Часть III1
Все проекты зело исправны быть должны, дабы казну зряшно не разорять и Отечеству ущерба не чинить. Кто станет абы как ляпать, того чина лишу и кнутом драть велю.
Петр I
Во многих сферах деятельности инвестиции, направленные на совершенствование орудий длительного пользования, зачастую осуществляются людьми, не являющимися собственниками этих орудий.
Артур Пигу
В части II мы подробно рассмотрели решение двух локальных постановочных задач: ми-нисуммной и максисуммной, для которых был выбран единый способ решения - алгоритм Егервари. На очереди - максиминная и минимаксная постановочные задачи. Напомним читателю их вербальные модели (приведем способ и обобщенный результат решения).
Максиминная постановка задачи
Способ решения - алгоритм Гросса. Макросреда условно-стабильная. Функционал (f): составить оптимальный инвестиционный портфель (ОИП), при котором самый большой объем инвестиций дочернего предприятия окажется минимально возможным, т. е. минимизировать max a[i, p(i)] по всем многовариантным версиям.
Как показывает решение, приведенное далее, во всех многовариантных версиях преобладает почти средний объем объявленных дочерними предприятиями инвестиций (три инвестиционных потока по 8, один - 4 и один - 6 млрд р.), что характеризует все еще осторожное отношение дочерних предприятий к изменению состояния макросреды.
Минимаксная постановка задачи
Способ решения - алгоритм Гросса. Макросреда стабильная. Функционал (f): составить ОИП, при котором самый маленький объем инвестиций дочернего предприятия окажется максимально возможным, т. е. максимизировать min a[i, p(i)] по всем многовариантным версиям.
Как показывает изложенное далее решение, дочерние предприятия стали более уверенными в макросреде, что подтверждается преобладанием в портфеле инвестиционных потоков, которые превышают средний объем объявленных дочерними предприятиями инвестиций (два инвестиционных потока по 11, один - 12, один - 13 и один - 15 млрд р.).
Теперь перейдем к развернутому решению каждой из двух постановочных локальных задач, используя для этой цели алгоритм Гросса (см. приложения).
1 Части I и II см. // Имущественные отношения в Российской Федерации. 2004, № № 3, 4.
ю
о
Таблица 1
Решение максиминной задачи (алгоритм Гоосса)
ю
СП
и
и
ю
о
о
а
з
5Г
<в
&
І
Содержание и иллюстрации в виде одномерных массивов, матриц (таблиц) и графов
с[
О
т
т
а(і, І) -» Ьі(і, і)
А[і, І] -> Вк[і, І
Функционал
Орграфы (0|)
к
з-
л
со
с
со
3"
X
f = В (большое число)
В цикле по I = 1-ь5, по] = 1'+6'
Р(!) = I, Маэ_и = {1,2, 3, 4, 5}. fstart = 14, Т. к. = тах(а[1,1]) (максимальный элемент из главной диагонали квадратной матрицы)_____________________
\
V
12
13
11
8
2'
11
12
3'
4'
10
10
14
12
10
11
5'
14
8
15
8
6'
10
Івіаіі —14
з»
В
ю
О
V
2'
3'
4'
5'
6'
В циклах по I = 1-5-5, по] = 1^6', если а[1, Л < 14, то Ь-|[1, Л = 1, иначе Ь-|[1, ^ = 0
Ві[У]
_1_
2
3_
А_
5
fnew — В-| < В При I = 3 находится новый функционал а[3, 1] = 11
|new
= 11
V
В циклах по \ = 1-н5, по] = 1'+6', если а[1, Д < 11, то Ьг[1, Л = 1, иначе Ьг[1, Л = 0
В2[У]
_3.
_4
5
2'
3'
4'
5'
6'
^е\« _ В2 ^ В1 При I = 1 находится новый функционал а[1, 3] = 10
fnew —10
Вопросы подписки по тел. (095) 129-0001
В цикле по \ = 1-И5, по] = Г-5-6', если а[1, Я < 10, то Ьзр, Л = 1, иначе Ьз[1, Л = 0 і \ V 2' 3' 4' 5' 6' — Вз < В2 (Т) © © © ©
1 0 0 0 0 0 1 При \ = 3 находится і
3 Вз[і,Л 2 0 1 1 1 1 1 новый функционал
3 0 1 0 1 0 1 а[3, 4] = 9
4 1 1 0 1 1 0 ^е\« = 9 © © © © (£) ©
5 1 0 0 0 1 1
В цикле по 1 - 1-5-5, по] = 1'+6', если а[1, Л < 9, то Ь^, Л = 1, иначе Ь^, Л = 0 V 2' 3' 4' 5' 6' © © © © ©
1 0 0 0 0 0 1 При \ = 3 находится
4 В4[і,л 2 0 1 1 1 1 0 новый функционал
3 0 1 0 0 0 1 а[3, 2] = 8
4 1 1 0 1 ї 0 = в(финиш) пч (?) © ®
5 1 0 0 0 ї ї
В цикле по 1 = 1-г-5, по] = Г-5-6', если а[1, Л < 8, то Ь5[1, Л = 1, иначе Ь5[1, Л = 0 \ V 2' 3' 4' 5' 6' Стоп! №з_и={6,3,-1,2,-1} Выходим из циклов, т. к. №з_[3] = -1, №б_[5] = -1 Лучший результат получен на итерации 4 Э Э 0
1 0 0 0 0 0 7
5 в5[і,Л 2 0 1 1 0 0 0
3 0 0 0 0 0 1
4 1 1 0 0 ї 0 © © © © ©®
5 0 0 0 0 0 ї
Результат f = 8 (см. итерацию № 4); Р - {6, 3, -1, 2, -1}. На базе таблицы В4/1, Л возможны комбинации: 1- (6, 3, 2,1, 5); II IV- (6, 3, 2, 5, 1); V- (6, 4, 2, 1, 5); VI- (6, 4, 2, 5, 1); VII- (6, 5, 2, 4, 1) -(6,3, 2,4, 1);///-(6, 3,2,4, 5);
Ю
Таблица 2
Продолжение решения максиминной задачи (формирование рациональных версий)
Двудольный орграф с указанием весов цветных дуг в наибольшем паросочетании Значение функционала (^, структура инвестиционного портфеля (№з и = { }), слагаемые весов наибольшего паросочетания Итоговая сумма инвестицион- ного портфеля
с4(1) (гГ (г; @ © f = 8 Маз_и = {6, 3, 2,1, 5} 6 + 4 + 8 + 5 + 8 = 31 31
с4(2) (у) (?) (у) {£) ®}\^ f = 8 Мае и = {6, 3, 2, 4,1} 6 + 4 + 8 + 8 + 8 = 34 34
с4(3) © (г; © © ф)(£) f = 8 Мае и = {6, 3, 2, 4, 5} 6 + 4 + 8 + 8 + 8 = 34 34
с4(4) f = 8 Мае и = {6, 3, 2, 5,1} 6 + 4 + 8 + 6 + 8 = 32 32
с4(5) 0^ (й/ ® f = 8 Мав.и = {6, 4, 2,1, 5} 6 + 8 + 8 + 5 + 8 = 35 35
с4(6) (2<© 5$ f = 8 Мав_и = {6, 4, 2, 5,1} 6 + 8 + 8 + 6 + 8 = 36 36
с4(7) f = 8 Мав_и = {6, 5, 2, 4,1} 6 + 8 + 8 + 8 + 8 = 38 38
Комментарий к решению максиминной задачи. Функционал f = 8 дает ответ на вопрос: когда при логистическом решении максиминной постановки задачи на основе теории графов, нужно выходить из циклов? Однако функционал далеко не всегда, например как в данном случае, дает однозначный ответ (единственное решение) о структуре инвестиционного портфеля, ибо итерация № 4 (решающая, см. табл. 1) предлагает лицу, принимающему решение (ЛПР), многоваринатные версии (орграфов G4 много - семь, см. табл. 2). Это объясняется особенностями управления инвестиционным портфелем: его оптимизация сложна и отличается субъективностью суждений о тенденциях как состояния микро-, макросреды, так и риска на фондовом рынке, об отдельных финансовых активах, намечаемых к включению в портфель. Эти обстоятельства порождают необходимость проводить дополнительно, уже в рамках имитационного моделирования на основе стратегии и тактики конкретного ЛПР, сложные многовариантные расчеты.
Цель экономиста-математика - в рамках постановочных задач подготовить для ЛПР математически обоснованные многовариантные версии структуры инвестиционного портфеля. В таблице 2 мы это проиллюстрировали, построив всевозможные орграфы G4(n) на базе таблицы В4р, ]] из итерации 4. Построение орграфов G4(n) выполнено по единичным решениям таблицы В4р, ]] «жадным» алгоритмом.
Решение данной максиминной задачи - наглядный пример многоцелевой, динамичной (многоступенчатой и гибкой) экономико-математической модели.
Ответ. Для имитационного моделирования в рамках нашей максиминной постановки задачи предлагается 7 (семь) версий структуры инвестиционного портфеля со следующими итоговыми суммами: 31, 32, 34, 34, 35, 36, 38 млрд рублей (рис. 1).
Развитие поиска однозначного решения для данной максиминной задачи состоит в формировании единственного ОИП на базе семи полученных рациональных версий, а это - предмет имитационного моделирования, на чем подробнее остановимся далее.
Рис. 1. Итоговые суммы рациональных решений (млрд р.)
Максиминная постановка задачи
Двумерное отображение решения в виде матрицы BM(i, \), таблицы B1[i, j] и совмещенного (раскрашенного) двудольного орграфа Grv (обобщение решения)
А(У) =
(\2 11 10 10 14 6
13 3 4 8 8 9
11 8 14 9 15 3
5 7 12 8 6 10
І8 12 10 11 8 6
ВМ (!>.)') =
0 0 0 0 0
0 3 4 8 8
0 5 0 0 0
5 7 0 8 6
8 0 0 0 8
6л
0
3
о
6
Вт[У]
V і \ 1' 2' 3' 4' 5' 6'
1 1
2 1 1 1 1
3 1 1
4 1 1 1 1
5 1 1 1
Примечание. В целях лучшего восприятия рисунка, последовательность в нумерации совмещенного (раскрашенного) двудольного орграфа Grv не соблюдена, а связи с отрицательной полярностью не показаны.
Расшифровка дуг орграфа Grv:
........ f = 8, Nas_u = {6, 3, 2, 1, 5}, cost_n = {6, 4, 8, 5, 8}, S = 6 + 4 + 8 + 5 + 8 = 31;
-------- f = 8, Nas_u = {6, 3, 2, 4, 1}, cost_n = {6, 4, 8, 8, 8}, S = 6 + 4 + 8 + 8 + 8 = 34;
-------- f = 8, Nas_u = {6, 3, 2, 5, 1}, cost_n = {6, 4, 8, 6, 8}, S = 6 + 4 + 8 + 6 + 8 = 32;
-------- f = 8, Nas_u = {6, 3, 2, 4, 5}, cost_n = {6, 4, 8, 8, 8}, S = 6 + 4 + 8 + 8 + 8 = 34;
-------- f = 8, Nas_u = {6, 4, 2, 1, 5}, cost_n = {6, 8, 8, 5, 8}, S = 6 + 8 + 8 + 5 + 8 = 35;
-------- f = 8, Nas_u = {6, 4, 2, 5, 1}, cost_n = {6, 8, 8, 6, 8}, S = 6 + 8 + 8 + 6 + 8 = 36;
-------- f = 8, Nas_u = {6, 5, 2, 4, 1}, cost_n = {6, 8, 8, 8, 8}, S = 6 + 8 + 8 + 8 + 8 = 38,
где f - значение функционала;
Nas_u = {...} - структура инвестиционного портфеля (качественная);
cost_n = {....} - структура инвестиционного портфеля (количественная);
X - общая сумма инвестиций ПО;
Grv - орграф рациональных версий (всего семь).
Вопросы подписки по тел. (095) 129-0001
Решение
Минимаксная постановка задачи
а
а
5Г
<в
&
І
Содержание и иллюстрации в виде одномерных массивов, матриц (таблиц) и графов
§
00
СП
Ґ\2 11 10 10 14 6Л
13 3 4 8 8 9
11 8 14 9 15 3
5 7 12 8 6 10
^8 12 10 11 8 6 у
о;
з
со
со
СЦ
со
3-
I
а(і, І) -> Ьі(і, і)
А[і, Л -> Вк[і, ]]
Функционал (1)
Орграфы (С|)
f = Э (малое число)
В циклах по I = 1-5-5,
по] = Г-5-6',
Рр) = I, Маэ_и = {1,2, 3, 4, 5}.
Там = з, т. к. = гп1п(а[1, I]). Минимальный элемент из главной диагонали квадратной матрицы
А[і. Л
1'
12
13
11
2'
3'
11
10
14
12
12
10
4'
10
8
8
11
5'
14
15
6'
6
10
6
3- Го
О 3
3'
4'
5'
В циклах по і = 1-5-5, по \ = 1'+6', если а[і, Л > 3, то Ь[і, Л = 1, иначе Ь[і, Л = 0
Ві[і, Л
— Э1 > Э.
При I = 2 находится новый функционал: а[2, 3] = 4; рр) = Маэ_и = {1,3,2,4,5}
В циклах по \ = 1-5-5, по] = Г-5-6', если а[1, Л >4, то Ь[\, Д = 1, иначе Ь[1, Л = 0
Вгр.Л
N 1' 2' 3' 4' 5' 6'
1 1 1 1 1 1 1
2 1 0 0 1 1 1
3 1 1 1 1 1 0
4 1 1 1 1 1 1
5 1 1 1 1 1 1
fnew — Эг > Э1.
При I = 2 находится новый функционал: а[2, 4] = 8; рР) = Маз_и = {1,4,3,2,5} = 8
ю
Сп
УПРАВЛЕНИЕ
ю
о
В циклах по \ = 1^-5, по] = Г-6', если а[\, Я > 8, то Ь[\, \] = 1, иначе Ьр, Л = 0 X 1' 2' 3' 4' 5' 6' 1пе\« — бз > Эг-При \ = 2 находится новый функционал: а[2, 6] = 9; р(!) = №8_и = {1,6,4,3,2} 7) © © © ©
1 1 1 1 1 1 0
3 ВзР.Л 2 1 0 0 0 0 1
3 1 0 1 1 1 0
4 0 0 1 0 0 1 (7 ) © © © © ©
5 0 1 1 1 0 0
В циклах по \ = 1н-5, по] = 1'-М5\ если а[\, \] > 9, то Ъ[\, \] = 1, иначе Ь[1, Л = 0 1' 2' 3' 4’ 5' 6' — §4 > Эз. При \ = 4 находится новый функционал: а[4, 6] = 10; р(0 = Маэ_и = {2,1,3,6,4} ^пе\« = 10 7) © © (7) ©
1 1 1 1 1 1 0
4 в4р, Л 2 1 0 0 0 0 0
3 1 0 1 0 1 0
4 0 0 1 0 0 1 (Г ) © © © © ©
5 0 1 1 1 0 0
В циклах по \ = 1н-5, по] = Г-6', если а[1, Л > 10, то Ь[1, Л = 1, иначе Ь[1, ]] = 0 X V 2' 3' 4' 5' 6' 1пе\« = Эб > 84. При 1 = 5 находится новый функционал: аГ5 41 = 11 ■ Р(0 = N88111 = {2.1,5,3,4} ^е\« — 11 7) © © © ©
1 1 1 0 0 1 0
5 в5р,Л 2 1 0 0 0 0 0
3 1 0 1 0 1 0
4 0 0 1 0 0 0 (Г ) © © © © ©
5 0 1 0 1 0 0
В циклах по \ = 1-5-5, по] = 1'+6\ если а[1, Л > 11, то Ь[1, Л = 1, иначе Ь[1, Л = 0 X 1' 2' 3' 4’ 5' 6' N38 и = {5,1,-1,3,2} э©@@©
1 1 0 0 0 1 0 Выходим из циклов, /81ор! / \
6 Вер, Л 2 1 0 0 0 0 0 т. к. пае [3] = -1.
3 0 0 1 0 1 0 Лучший результат
4 0 0 1 0 0 0 найден на итерации5 © ) © © ©©) ©
5 0 1 0 0 0 0
Ответ Г = 11, Р = {5, 1,-1, 3,2}, Го = 11 + 13 + 15+ 12 + 11 = = 62
10
СП
и
и
N
О
о
Переход от моделей теории графов к моделям имитации выбора инвестиционных решений
Имитационная модель выбора инвестиционных решений определяется как модельное описание действительного хода процесса с помощью определенной системы понятий и ограниченного набора показателей.
Отличительная особенность метода имитационного моделирования состоит в том, что при разработке модели заранее не делается предположений о принадлежности ее к какому-либо специальному классу, не используется какая-либо специальная схема, но описание процесса в избранной системе показателей выполняется по возможности точно и близко к реальности, несмотря на импровизацию ЛПР.
Имитационную модель не следует понимать как уже установившуюся, фиксированную форму, она сама служит полезным средством для своего собственного усовершенствования и уточнения, скорее - это саморазвивающаяся система. В самом деле, оказывается, что создание и использование имитационной модели в прикладных целях открывает путь к анализу и оценке уже построенной модели (в нашем случае это подтверждается необходимостью развития решения максиминной постановки задачи), а затем к ее же самосовершенствованию.
Необходимость применения имитационного моделирования к задаче, вербальная модель которой изложена выше, возникает по причине того, что максиминная постановка задачи имеет не одно, а несколько решений: семь версий структуры потенциальных инвесторов с шестью итоговыми суммами инвестиций, а наша конечная цель - построить ОИП в широком смысле, когда количественные и качественные показатели не имеют альтернативы, конечно же, с поправкой на предпочтения ЛПР.
Далее покажем наше видение продолжения решения максиминной постановки задачи с помощью пяти показателей (наших условных предпочтений), оцененных по десятибалльной системе (таблицы 3а и 3б). Согласно последнему столбцу таблицы 3б (Всего...) выбор инвестиционного решения выпадает на одномерный массив в строчке ] = 2', которому по оценке соответствует 211 баллов.
Результат. Из семи решений максиминной постановки задачи с помощью имитационного моделирования на базе предпочтений ЛПР нами определен четвертый ОИП со следующими количественными и качественными характеристиками:
Г = 34 млрд р. = Г0 ^ Nas_u = {6, 3, 2, 4, 1}, Р0 = {6, 4, 8, 8, 8}.
Замечание. Предложенный вариант имитационного моделирования наглядно показывает сложность выбора инвестиционного решения, ибо составленный четвертый ОИП отличается от других версий всего одним или четырьмя баллами (+0,5 %; +1,9 %), за которыми стоит разность в - 4; - 2; - 1; 0; + 2; + 3 млрд р. (- 10,5 %; - 5,6 %; - 2,9 %; 0 %; + 6,3 %; + 9,7 %), т. е. сохраняется возможность для совершенствования имитационной модели в процессе ее эксплуатации.
Таблица 3а
Система предпочтений ЛПР к эмитентам ( = 1'+6') по десятибалльной шкале
№ п/п Содержание предпочтений ЛПР Эмитенты 0)
V 2' 3' 4' 5' 6'
1 Отношения эмитента с государственными органами управления и контроля 7 10 8 6 9 10
2 Прозрачность эмитента 10 8 6 9 7 8
3 Общее экономическое состояние эмитента 10 9 10 9 9 10
4 Рейтинг эмитента в отрасли 8 9 10 9 8 7
5 Социальный климат в коллективе эмитента 9 6 10 8 7 5
6 Итого (баллов) 44 42 44 41 40 40
Таблица 3б
Определение ОИП на базе решений максиминной постановки задачи с помощью имитационного моделирования (табл. 3а)
\ -1 л 1 2 3 4 5 Определение ОИП
Решения максиминной задачи (№з и = { }) Слагаемые предпочтений Всего (баллов)
г 44 Маэ_и = {6, 3, 2, 1, 5} 40 + 44 + 42 + 44 + 40 210
Маэ_и = {6, 3, 2, 5, 1} 40 + 44 + 42 + 40 + 44 210
2' 42 Мав_и = {6, 3, 2,4, 1} 40 + 44 + 42 + 41 + 44 211
3' 44 Маэ_и = {6, 3, 2, 4, 5} 40 + 44 + 42 + 41 + 40 207
4' 41 Маэ_и = {6, 4, 2, 1, 5} 40 + 41 + 42 + 44 + 40 207
5' 40 Мав_и = {6, 4, 2, 5, 1} 40 + 41 + 42 + 40 + 44 207
6' 40 Маэ_и = {6, 5, 2,4, 1} 40 + 40 + 42 + 41 + 44 207
Примечание. В целях лучшего визуального восприятия и сохранения логической цифровой цепочки таблица 3б составлена максимально упрощенно: клетки по горизонтали (кроме диагональных) не заполнены, т. к. значения в них совпадают (по определению) со значением, указанным в диагональной клетке данной горизонтали.
Замечание. Как правило, запись имитационной модели осуществляется на обычном алгебраическом языке. Однако для расчетов на ЭВМ необходимо дать представление этой модели на некотором входном языке ЭВМ. Здесь возникает проблема подбора наилучшего машинного языка, соответствующего данной имитационной модели, и формирования машинной имитационной модели в режиме человеко-машинного диалога2.
Завершим изложение первых 3-х частей исследования иллюстрацией схемы «Последовательность (этапы) формирования ОИП коллективным инвестором в широком экономическом смысле (логистический подход: микро- и мезоэкономические аспекты)». Наблюдательный читатель без труда убедится в эффективности алгоритмов неконтинуальной математики, которые применены нами в единичных решениях локальных постановочных задач, что подтверждает возможность и целесообразность сосуществования единичных и программных (электронных) решений. Такое сосуществование увеличивает не только степень свободы ЛПР, экономистов и практиков, но и пространство их предпочтений при выборе ими средств для реализации своих тактических планов.
Расшифровку этапов формирования ОИП на нашем конкретном материале можно представить как:
1) формирование Ры - базового инвестиционного портфеля на основе оценок эффективности инвестиционных проектов, так называемая «черновая» работа, когда коллективный инвестор в лице его доверенных экономистов-математиков освобождается от предложений эмитентов, не вписывающихся в его инвестиционную стратегию;
2)формирование предпочтений (объявленных инвестиционных потоков) дочерних предприятий - построение матрицы А(1, ]);
3) построение неориентированного графа G <1, J>, который в нашем случае порождает 720 различных комбинаций ориентированных графов;
2 Багриновский К.А., Егорова Н.Е. Имитационные системы в планировании экономических объектов. М.: Наука, 1980.
4) формирование ОИП [как возможный промежуточный этап - формирование рациональных инвестиционных портфелей (см. единичное решение максиминной постановочной задачи)] для конкретного состояния макросреды - Gs <1, J>, при этом соотношение «результат - исходная база данных» измеряется величиной 1/720, что соответствует
0,139 процента;
5) формирование ОИП для всех состояний макросреды: непредсказуемого, условностабильного, стабильного и динамично развивающегося - построение G0 <1, J>. При необходимости алгоритмы теории графов (Гросса и Егервари) могут дополняться моделями и методами имитации выбора из рациональных портфелей оптимального инвестиционного на основе предпочтений конкретного ЛПР;
6) оформление сводной таблицы - количественные (инвестиционные потоки) и качественные (структура коллективного инвестора) показатели ОИП.
Количественный результат логистического подхода к формированию ОИП коллективным инвестором - отбраковано 99,444 процента предпочтений дочерних предприятий (удельный вес четырех ОИП, «отфильтрованных» из 720 различных комбинаций предпочтений коллективного инвестора составил 0,556 процента).
Примечание. Наше единичное решение комплексной задачи не только подтверждает истину: в условиях неопределенности «здравый смысл» нередко приводит к ошибочным заключениям, но и защищает ЛПР от возможных негативных последствий использования компьютеров, о чем предупреждают нас Майкл Мичем и Третий закон Гоида (см. // Имущественные отношения в Российской Федерации. 2004. № 7).
Замечание. Имитационное моделирование возникает всякий раз, когда на поставленный ЛПР вопрос нет однозначного решения. При формировании ОИП коллективным инвестором по предлагаемой нами концепции (см. // Имущественные отношения в Российской Федерации. 2004, № № 3, 4) такая ситуация может возникнуть при решении любой локальной постановочной задачи - все зависит от конкретного материла, т. е. от базового инвестиционного портфеля (Ры), составленного на основе оценок эффективности инвестиционных проектов и предпочтений дочерних предприятий - А(1, ]). Требовательный и взыскательный читатель может в этом убедиться в режиме реального времени. Для этого достаточно изменить некоторые параметры Ры и А(1, ]) и решить самостоятельно все локальные постановочные задачи.
С вербальными результатами нашего исследования, итогами применения многоцелевых, динамичных (многоступенчатых и гибких) экономико-математических моделей, выводами и др. читатель сможет ознакомиться в следующем номере журнала, а с тезисами автора «Методы исследования и моделирования механизма формирования оптимального инвестиционного портфеля (ОИП)» в материалах 5-го Всероссийского симпозиума «Стратегическое планирование и развитие предприятий»3.
3 Материалы 5-го Всероссийского симпозиума «Стратегическое планирование и развитие предприятий» / Под ред. проф. Г.Б. Клейнера. М.: ЦЭМИ РАН, 2004. 178 с.
Схема
Последовательность (этапы) формирования ОИП коллективным инвестором в широком экономическом смысле (логистический подход: микро-и мезоэкономические аспекты)
{5+13; 3+12; 4+14; 8+11; 6+15; 3+10}
А(У) =
V
Ры 1 — (5+ 13;
12 11 ▼ 10 10 14 6
13 3 4 8 8 9
11 8 14 9 15 3
5 7 12 8 6 10
8 12 10 11 8 6
Л
0<1,>
у
Оо<1,Х>
Оптимальный инвестиционный портфель (ОИП) Справочно
Качественная структура Л/ав и = { } Количественная структура Ры = {.....} Общий объем инвестиций в млрд р. (1о) Макросреда Категория постановки задачи
Мав_и = {4, 2, 6, 1,5} Роі = {1 о, 3, 3, 5, 8} 29 Непредсказуемая Минисуммная
Маэ_и = {6, 3, 2, 4, 1} РЗІ = {6, 4, 8, 8, 8} 34 Условно-стабильная Максиминная
Маэ_и = {2, 1,5,3, 4} Роі = {11,13,15,12,11} 62 Стабильная Минимаксная
Мав_и = {5, 1, 3, 6, 2} Роі = {14,13,14,10,12} 63 Динамично развивающаяся Максисуммная
Приложение 1
Исходные данные: матрица A(i, j), граф С<!, J>
где А(1, ]) - матрица предпочтений дочерних предприятий;
i = (1^5) - номер (условный) дочернего предприятия (строка матрицы);
] = (1'^6') - номер (условный) эмитента (столбец матрицы);
а(1, ]) - объем заявленных инвестиций ьго дочернего предприятия в ]-й вид акций инвестиционного портфеля, вес ребра и(1, ]);
G<I, J> - двудольный неориентированный граф, построенный по данным матрицы Ар, ]);
I - множество дочерних предприятий, I = {1, 2, 3, 4, 5};
J - множество финансовых активов (видов акций), J = {1',2',3',4',5',6'}; ир, ]) - ребра, соединяющие множества I и J на основе матрицы Ар, ]); f - функционал (критерий);
Nas_u - одномерный массив, инвестирование I в J, Nas_u = {.......};
f = ?, Nas_u = ?
Постановка задач, способы их решения и рассматриваемые состояния макросреды в обобщенном виде
^^^^ГЛакросреда КПЗ/СР Непредсказуемая (А) Условно- стабильная (В) Стабильная (С) Динамично развивающаяся (0) Производный эффект
Минисуммная^^^^ ^^^^алгоритм Егервари ? ? ? ? ?
Максиминная^^^^^ ^^^^алгоритм Г росса ? ? ? ? ?
Минимаксная^^^^^ ^^^^лгоритм Г росса ? ? ? ? ?
Максисуммная^^^^ ^^^^алгоритм Егервари ? ? ? ? ?
где КПЗ - категория постановки задачи (минисуммная, максиминная, минимаксная, мак сисуммная);
СР - способ решения задач (алгоритмы Егервари, Гросса);
А, В, С, D - итоговые суммы оптимального (рационального) инвестиционного портфеля (А < В < С < D);
? - поиск результатов.
Приложение 44
Математическое описание алгоритма Гросса
5 КАТЕГОРИЯ ПОСТАНОВКИ ЗАДАЧИ
Итераци МИНИМАКСНАЯ МАКСИМИННАЯ
Ввод Ар,)] Ар, \]
Инициализация р[М] - стартовое (начальное) решение, присвоить р: = (1, 2...М) - число на месте показывает номер инвестиционного проекта в портфеле, в котором будет принимать участие не дочернее предприятие ПО; f - значение функционала, присвоить Ь = М!п(а[!, I]) - минимальный элемент на главной диагонали квадратной матрицы, размещенной внутри Ар, Л. f = Э (малое число) р[М], р(!) = I, №8_и = {1, 2...М} Л = Мах(а[1, 1]) - максимальный элемент на главной диагонали квадратной матрицы, размещенной внутри Ар, Л f = В (большое число) С0 - двудольный орграф
>э _ Э. (0 з а § з 1 Породить матрицу Вр, Л: в цикле по \ = 1-5-М, по] = 1-5-М если ар, Л > Э, то Ьр, Л = 1, иначе Ь[1, Л = 0 если ар, Л < В, то Ьр, Л = 1, иначе Ьр, Л = 0
2 Использовать алгоритм нахождения наибольшего паросочетания, ввод матрица В[1, Л, выход двудольный орграф Сч, матрица №з_и Использовать алгоритм нахождения наибольшего паросочетания, ввод матрица Вр, Л, выход двудольный орграф 01, матрица №э_и
3 Проверить, является ли паросочетание наибольшим = Э1 > Э (Эч — большее число), в цикле ПО \ = 1-5-М, если паэ_и[1] ф -1 и ар, пав_и(1)] > Э, то = ар, паэ_и(1)] = иначе выход (паросочетание не является наибольшим и лучший результат был получен ранее) \ (паросочетание наибольшее, обновляем значение функционала и улучшаем результат) в цикле по \ = 1-5-М рр] = паз_ир] = В1 < В (В-| - меньшее число)
Вывод Р^
Ответ Сумма весов лучшего результата - есть объем инвестиций ПО по оптимальному портфелю, а одномерный массив Мае и = { } раскрывает внутреннюю структуру участия дочерних предприятий в инвестициях ПО Сумма весов лучшего результата есть объем инвестиций ПО по оптимальному портфелю, а одномерный массив Мае и = { } раскрывает внутреннюю структуру участия дочерних предприятий в инвестициях ПО
4 Приложения 2 и 3 см. // Имущественные отношения в Российской Федерации. 2004. № 4.
Приложение 5
Комментарий к математическому описанию алгоритма Гросса
Минимаксная постановка. В условиях стабильного состояния макросреды совет директоров объединения разрабатывает инвестиционную стратегию, предусматривающую такую структуру оптимального инвестиционного портфеля, при которой наименьший объем инвестиций дочернего предприятия окажется максимально возможным, т. е. максимизировать min a[i, p(i)] по всем многовариантным версиям.
Максиминная постановка. В условиях условно-стабильного состояния макросреды совет директоров объединения разрабатывает инвестиционную стратегию, предусматривающую такую структуру оптимального инвестиционного портфеля, при которой наибольший объем инвестиций дочернего предприятия окажется минимально возможным, т. е. минимизировать max a[i, p(i)] по всем многовариантным версиям.
В рамках теории графов в качестве модели каждой из задач мы получим двудольный орграф, ребра которого имеют неотрицательные веса. Обе задачи могут быть решены посредством одного алгоритма: дан двудольный граф G<M, N> с взвешенными ребрами, A(i, j) - матрица весов. Найти наибольшее паросочетание M* с тем, чтобы:
а) минимальный вес по всем ребрам, вошедшим в паросочетание M*, был максимальным по всем M*;
б) максимальный вес по всем ребрам, вошедшим в паросочетание M*, был минимальным по всем M*.
Основная идея.
Составить начальное наибольшее паросочетание - стартовая версия структуры инвестиционного портфеля (по главной диагонали базовой матрицы). Согласно целям задачи трансформировать базовую матрицу A(i, j), в матрицу B(i, j), для которой попытаться построить наибольшее паросочетание и, если это удастся сделать, попробовать вновь видоизменить матрицу A(i, j) в B'(i, j), запомнив при этом значение функционала. Это и есть алгоритм Гросса.
Возможные сложности связаны с разрешением следующих вопросов.
1. Как составить первое наибольшее паросочетание?
2. Как видоизменить (трансформировать) базовую матрицу A(i, j), чтобы следующий результат (значение функционала f) был не хуже предыдущего?
3. Когда завершить процесс улучшения значения функционала (f)?
Способы преодоления.
1. Для первого (начального) инвестиционного портфеля выбрать диагональные элементы базовой матрицы A(i, j), т. е. i = 1^M, j = 1^M (<N) и i = j, а это означает, что 1-ое дочернее предприятие инвестирует в 1-й вид ценных бумаг, 2-ое во 2-й и т. д.
2. Порождение матрицы B(i, j) производится по правилу:
b(i, j) = 1, если размер инвестиций i-го инвестора в j-й вид ценных бумаг больше найденного на предыдущем шаге значения функционала (Sn > Sn-1 >...> S1 > S0) и b(i, j) = 0, в противном случае - для минимаксной задачи;
b(i, j) = 1...(Bn < Bn-1 <...< B1 < B0) и b(i, j) = 0, в противном случае - для максиминной задачи.
3. Признаком окончания итерационного процесса является невозможность на определенном (n + 1) этапе построить наибольшее паросочетание. Тогда лучший результат (нижняя или верхняя «планка» инвестиций) и соответствующая структура инвестиций определены на предыдущем (n-м) шаге.
Замечание. Решение по инвестиционному портфелю может быть: а) единственным (как по объему инвестиций, так и по структуре); б) единственным по объему инвестиций и многовариантным по структуре; в) многовариантным (как по объему инвестиций, так и по их структуре, причем вариантов по структуре всегда не меньше, чем вариантов по объему).
Окончание в следующем номере журнала