В. Давнис Прогнозные модели экспертных предпочтений
.pdfмногомерного прогнозирования. Поэтому логичным продолжением развития идеи комбинированного прогнозирования является созда ние комбинированных моделей для многомерных процессов. Ос новой построения таких моделей может служить адаптивный мат ричный предиктор, различные варианты которого разработаны в предыдущих параграфах настоящей главы. Ниже приводится до статочно подробное алгоритмическое описание построения одно го из вариантов подобной модели.
I.Набор данных, имеющихся на момент построения модели
1.Динамика показателей:
хи |
х12 |
••• |
х1п, |
Х2\ |
-^22 |
' ' ' |
Х2п' |
Xil |
Xl2 |
'" |
ХнГ |
2. Ожидаемые темпы роста показателей:
вариант |
1: qn |
qt2 |
••• |
qm; |
вариант |
2: q\ |
q]2 |
••• |
q]n\ |
вариант |
m: q't" |
q"2 |
••• |
q'"x. |
II. Обозначения
1. Константы:
N — число имитационных экспериментов; п — число моделируемых переменных; М — число экспертов; т — число альтернативных вариантов;
т — число прогнозных периодов; в — число наблюдений для построения начального значения
предиктора; / — число наблюдений, используемых для построения прогноз
ной модели.
2. Оцениваемые параметры и величины:
Q* — диагональная матрица, диагональными элементами ко торой являются ожидаемые темпы роста показателей по k-му ва рианту (к = 1, т);
190
dp s. — математическое ожидание и среднеквадратическая ошибка случайной величины, характеризующей изменения абсо лютных приростов /-го показателя (г = 1, /г);
d, s — математическое ожидание и среднеквадратическая ошибка случайной величины, характеризующей ожидаемые значе
ния элементов весовой матрицы W; |
|
|
|
W — матрица весовых |
коэффициентов w., |
характеризующих |
|
вклад j-ro показателя в прирост /'-го, |
имеющая |
вид |
|
|
1 |
|
|
XV- |
Пп |
|
|
W = |
21 |
|
|
|
л 2 |
1 |
|
|
XV |
|
|
В некоторых случаях, когда число прогнозируемых переменных невелико, элементы матрицы W могут быть настраиваемыми.
В имитационных экспериментах элементы матрицы W форми руются случайным образом.
3. Настраиваемые параметры:
а — параметр адаптации, регулирующий в предикторе соотно шение старых и новых тенденций;
а* — оптимальное значение параметра а; ц — параметр, регулирующий комбинирование прямых и кос
венных темпов прироста. В предикторе используется матрица М,
элементы которой зависят от ц: |
|
||
|
1-/л |
1-И |
1-М |
М = |
|
1-М |
|
|
|
|
|
|
1-/л |
\-ц |
|
М* — матрица М с оптимальным |
значением параметра /и*; |
к — параметр, регулирующий вклад в прогнозную траекторию экстраполяционных тенденций и экспертных ожиданий. В блоч ном предикторе используется диагональная матрица с элементами,
зависящими от X: |
Л |
О |
о |
|
|||
Л, |
О Л2 |
о |
|
|
|
|
|
|
о |
о |
Лг |
191
[Л |
0 • •• |
0" |
!лт0 |
0 |
• •• |
0 |
0 Л • • |
0 |
ят |
••• |
0 |
||
|
|
; |
А Т = |
|
|
|
,0 |
0 • • к ^ 0 |
0 • • ят |
4. Случайные величины:
й/у — элемент /-й строки y'-го столбца матрицы w v , форми руемый датчиком случайных чисел в V-M имитационном экспери менте;
W" — матрица весовых коэффициентов, сформированная из имитируемых датчиком случайных чисел элементов wjj;
5,v+/, — прогнозная ошибка /-го показателя в периоде t+l, ими тируемая датчиком случайных чисел в V-M эксперименте;
v,+/,y — элемент /-й строки y'-го столбца матрицы Vv, рассчи тываемый по имитируемой случайной ошибке 8^+и в V-M экспе рименте для периода t+l;
V" — комбинированная матрица прямых и косвенных темпов прироста, сформированная по имитируемым прогнозным ошиб кам в v-м эксперименте;
Av — корректирующая матрица в v-ом имитационном экспе рименте;
А^+. — блочный предиктор для расчета вектора прогнозных оценок по всем упреждающим периодам в V-M имитационном экс перименте;
Х)'+. — вектор прогнозных оценок для всех упреждающих пе риодов, полученных в v-ом имитационном эксперименте.
5. Процедура моделирования случайной величины:
F — процедура моделирования псевдослучайных чисел, имити рующих элементы весовой матрицы;
Т — процедура моделирования псевдослучайных чисел, имити рующих прогнозные ошибки.
6. Процедура экспертного оценивания:
Е — процедура интерактивного экспертного оценивания веро ятности реализации того или иного прогнозного варианта.
7. Результаты экспертного оценивания:
~Рн-и ~ вероятностная оценка, полученная от r-го эксперта, относительно реальности А:-го варианта /-го показателя для момен та времени / + /;
"p,k+li — групповая экспертная оценка относительно реальности k-ro варианта /-го показателя для момента времени t+l;
192
xt+n ~ усредненное по групповым экспертным оценкам значе ние /-го показателя для момента времени / + /;
~х,+и — наиболее вероятное значение /-го показателя для мо мента времени по результатам экспертного опроса.
8. Моделируемые переменные и величины:
x,i — значение /-го показателя в момент времени /;
хн — расчетное значение /-го показателя в момент времени /; xt+u ~ расчетное значение /-го показателя по к-му варианту в момент времени / + / (с целью упрощения алгоритма для каж
дого / задается фиктивный вариант х,'"!1, принимающий макси мально возможное значение);
X, =(хп, х,2, •••, х,„) — вектор-столбец из значений моделиру емых показателей в момент времени t,
х, — вектор расчетных значений моделируемых показателей для момента времени t,
X, =(х,,х,,...,х,) — вектор-столбец с компонентами из одина ковых вектор-столбцов хг;
X =(х+1, x,+,,...,x,+r) ~~ вектор-столбец с компонентами из вектор-столбцов прогнозных значений моделируемых показателей;
Х^+. =(х^+|,х^+2,...,х^+г) — вектор-столбец значений комбини рованной прогнозной траектории для моментов времени / + 1, / + 2,
•••'r Zт;
Mt+u — индикаторная переменная, принимающая значение О или 1 в зависимости от результата сравнения полученной величи ны /-го показателя в v-м эксперименте для периода / + 1 с соот ветствующим к-м вариантом;
Ае — начальное значение предиктора;
Ve — начальное значение комбинированной матрицы прямых и косвенных темпов прироста;
V,+. — корректирующая матрица прямых и косвенных темпов прироста для моментов времени t + \, t + 2, , t + г,
Vt+uj ~ элемент /-й строки у'-го столбца матрицы V,+., опре деляемый по отклонениям постпрогнозных расчетов от фактичес ких значений для периода / + 1;
Z(a, /л) — критерий настройки параметров адаптации а и fi. А^+. _ скорректированное на прогнозные ошибки значение
предиктора для моментов времени в + 1, в + 2, ..., в + г,
193
Ar + . |
— корректирующая |
матрица для моментов времени |
t + 1, |
|||||||||||||
Г +2,..., |
t + |
г, |
|
|
|
|
|
|
|
|
|
|
|
|
||
А,+ . |
— предиктор для моментов |
времени |
/ + 1, / + 2, ..., t+ г, |
|||||||||||||
рк |
и |
— вероятность |
реальности |
|
значения |
/-го показателя по |
||||||||||
к-му варианту в момент |
времени t + 1; |
|
|
|
|
|
||||||||||
А? + . — предиктор комбинированной модели для моментов вре |
||||||||||||||||
мени |
t |
+ 1, t + 2, ..., |
t + |
т . |
|
|
|
|
|
|
|
|||||
III. Вычислительная |
схема |
|
|
|
|
|
|
|
|
|||||||
1. |
Расчет |
альтернативных |
вариантов: |
|
|
|
|
|
||||||||
|
|
( S-J |
|
\ |
(глк |
|
о |
|
|
|
О |
|
V |
х |
Л |
|
|
|
хк |
|
|
Q, |
|
|
|
|
|
|
|
|
|||
|
|
|
|
О |
|
Q*Q* |
|
|
О |
|
|
|
k = l,m. |
(5.41) |
||
|
|
|
|
|
О |
|
0 |
|
••• |
Qk |
Qk |
\х>; |
|
|
||
2. |
Построение |
начального значения |
предиктора: |
|
||||||||||||
|
|
|
|
|
Д*/ |
= |
ч>; |
|
|
|
/ = 1, |
п; |
|
|
|
(5.42) |
|
|
|
|
|
|
|
0 - 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
vy |
|
, |
i, |
j = |
l,n; |
|
|
(5.43) |
||
|
|
|
|
|
|
|
|
v0 = |
|
|
|
|
|
|
(5.44) |
|
|
|
|
|
|
Ae =(l-M*W*Ve )-l |
|
|
(5.45) |
||||||||
3. Настройка параметров предиктора по постпрогнозным расче |
||||||||||||||||
там: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
fх |
|
\ |
(к |
О |
|
О Vx, > |
|
|
||||||
|
|
|
|
J + 2 |
|
О |
А2е |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(5.46) |
|
|
|
|
|
|
|
О |
О |
|
|
Ag j ухв |
; |
|
|
||
|
|
)+/;/ |
_ хе+« |
-*•'е+к , |
/, j |
= 1, п, |
|
I = 1, т; |
(5.47) |
|||||||
|
|
|
|
|
%+/i |
|
|
|
|
|
|
|
|
|
194
|
v e + . = |
e+lij\ |
|
(5.48) |
|
A0 + . =(l - M*W*V 0 + . ) - 1 ; |
|
(5.49) |
|||
A-£+.=a*Ae+.+(l-a%+.A |
e+. |
(5.50) |
|||
|
|
|
|
|
|
|
A |
- Aad |
|
(5.51) |
|
|
X0+. - A0+#X0; |
|
(5.52) |
||
Z(a, fi) = max -%+// |
-^fl+Zi |
/ = 1, n, / = 1, т; |
(5.53) |
||
|
ve+« |
|
|
|
|
a |
= Argmin Z(a, /i); |
|
(5.54) |
||
|
|
||||
v^ , |
|
|
|
||
|
|
|
|
||
|
V,+ .= |
'f+fty |
|
(5.55) |
|
|
|
|
|
|
(5.56) |
A,+. = a*A0+. + (l - a* )&,+.Ae+. |
(5.57) |
||||
|
|
|
|
|
(5.58) |
4. Имитационные расчеты: |
|
|
|
||
|
|
|
|
|
(5.59) |
|
|
|
|
|
(5.60) |
|
1, если i = j ; |
|
|
||
Щ |
w-j / w\, если i * j ; |
|
(5.61) |
||
|
|
|
|||
|
ГУ |
|
|
|
|
- V |
(5.62) |
|
|
|
(5.63) |
195
~v |
't+li |
i, j = l,n, I = 1, т; |
V,t+lij |
|
|
|
V v = |
/+///• |
|
|
|
|
|
|
|
|
AV =&-M**WV , *VV )"1 ; |
|||||
Al.=a*A„.+(l-a*)kvAt+.; |
|
||||
|
X ; + . = A,+ .X,, |
|
|||
J0, |
если |
х^+1 |
< xt+l, |
|
|
• |
|
к |
„ -v |
^ k+\ |
k = l,m', |
[l, |
если |
x,+l |
< x,+i |
< *,+i |
, |
#,; = ЛГ1 £й#й
v = \
5. Экспертное оценивание
P = E(X, P);
~k |
1 |
M |
~kr |
1 v |
|||
Pt+u = 7 7 L, |
Л+ |
||
~ |
_ Д |
~* |
.* |
-*?+// |
|
2-iPr+uxt+ii; |
(5.64)
(5.65)
(5.66)
(5.67)
(5.68)
(5.69)
(5.70)
(5.71)
(5.72)
(5.73)
|
|
|
|
(5.74) |
6. Расчет значений |
комбинированной прогнозной траекторииш |
|||
|
|
Хг + . = А,+ .Х,, |
(5.75) |
|
_ |
xt+li |
~~ xi+li |
i, j = l,n, l = 1, т ; |
(5.76) |
"/+/// |
|
' |
||
|
|
*/ +// |
|
|
|
|
V ; + . = |
-f+/ |
(5.77) |
|
|
|
|
(5.78) |
A,c+. = ЛА Г + . + (I - Л ) А , + . А , + # ; |
(5.79) |
196
Хгс+.=А,С+.ХГ |
(5.80) |
IV. Описание вычислительной схемы модели
Вычислительная схема модели представляет собой систему ин формационно связанных между собой расчетных блоков. В пер вом из них осуществляется расчет альтернативных вариантов мо делируемых показателей в соответствии с ожидаемыми темпами роста этих показателей, задаваемых матрицей со специально уст роенной блочной структурой диагонали. В следующем блоке оп ределяется начальное значение предиктора по средним темпам прироста каждого показателя. Для этого используется известная схема построения комбинированной матрицы прямых и косвенных темпов прироста. Для удобства в расчетную формулу предиктора введена операция блочного умножения.
Далее с помощью полученного начального значения строится блочный предиктор, позволяющий осуществлять прогнозные рас четы сразу на несколько периодов. Путем последовательных кор ректировок и одновременного соответствующего подбора парамет ров адаптации он трансформируется в текущий предиктор с оп тимальными параметрами.
Текущий предиктор используется в блоке имитационных расче тов, в котором с помощью датчика случайных чисел получают матрицу весовых коэффициентов и комбинированную матрицу прямых и косвенных темпов прироста. Полученные матрицы ис пользуются для построения корректирующей матрицы, которая, в свою очередь, применяется для корректировки текущего предик тора. Скорректированный предиктор, по сути, является случай ным оператором, имитирующим одну из возможных многомерных траекторий развития моделируемых процессов. Многократное об ращение к датчикам случайных чисел позволяет воссоздать доста точно полную картину будущего. Сравнение имитируемых вари антов с альтернативными обеспечивает построение вероятностного распределения, характеризующего степень реальности альтернатив ных вариантов. Именно это вероятностное распределение и явля ется основным результатом комплекса расчетов этого блока, ко торый затем подлежит экспертному оцениванию.
В целом, блок экспертного оценивания предназначен для формирования оценок, отражающих экспертное представление о будущей динамике моделируемых показателей. Ориентиром экс пертных ожиданий является эмпирическое вероятностное распре-
197
деление, полученное в четвертом блоке. Строгих формальных правил получения экспертных ожиданий моделью не предусматри вается. Возможен, например, интерактивный режим формирова ния экспертных оценок.
Последний блок вычислительной схемы реализует главный за мысел модели — комбинирование на принципах адаптации экстраполяционных оценок и экспертных ожиданий.
Ниже приводятся результаты вычислительного эксперимента, иллюстрирующие весь комплекс расчетов, предусмотренных этой моделью. В расчетах использовались фактические данные, харак теризующие экономическое развитие Воронежской области. С целью упрощения для более легкого восприятия логики расчетов в соответствии с вычислительной схемой комбинированной моде ли прогнозирования многомерных процессов из всех макроэконо мических показателей были выбраны только три — валовой регио нальный продукт (ВРП), продукция промышленности и продук ция сельского хозяйства, динамика которых приведена в табл. 5.3.
Для проведения расчетов были сформированы три альтернатив ных варианта возможного развития Воронежской области в бли жайшие два года (2003—2004 гг.). Учитывая, что исходные дан ные содержат инфляционную составляющую (вариант I — 14 %, вариант II — 15 %, вариант III — 16 %), в эти варианты были внесены поправки на ожидаемые уровни инфляции. При разра ботке вариантов для упрощения расчетов принято допущение, что темпы роста и темпы инфляции остаются неизменными на про тяжении 2003—2004 гг. В табл. 5.4 представлены данные о трех вариантах возможной динамики макроэкономических показателей.
Т а б л и ц а 5.3
Основные макроэкономические показатели экономического развития Воронежской области
(1995-1997 гг. - млрд р., 1998-2002 гг. - млн р.)
Годы |
Валовой региональный |
Продукция |
||
продукт |
промышленности |
сельского хозяйства |
||
|
||||
1 |
2 |
3 |
4 |
|
1995 |
14 547,90 |
9598,00 |
4118,80 |
|
1996 |
20 158,30 |
12 909,00 |
5731,00 |
|
1997 |
23393,10 |
15 384,70 |
6902,70 |
|
1998 |
24075,00 |
14 653,30 |
6783,80 |
|
1999 |
40 710,10 |
24 237,10 |
15 500,00 |
198
|
|
О к о н ч а н и е |
т а б л . 5.3 |
1 |
2 |
3 |
4 |
2000 |
52 100,40 |
33151,50 |
18 247,20 |
2001 |
63 866,90 |
41 623,90 |
22 300,20 |
2002 |
75284,00 |
51 573,10 |
25443,70 |
Т а б л и ц а 5.4
Альтернативные варианты экономического развития Воронежской области
|
|
Валовой регио- |
|
Продукция |
|
||
Годы |
№ |
нальный продукт |
промышленности |
сельского хозяйства |
|||
варианта |
Значение, |
Темп |
Значение, |
Темп |
Значение, |
Темп |
|
|
|
млн р. |
роста, % |
млн р. |
роста, % |
млн р. |
роста, % |
|
I |
88 827,59 |
118,0 |
60 439,55 |
117.2 |
29 875,99 |
117,4 |
2003 |
11 |
90 997,28 |
119,4 |
62 277,61 |
119,0 |
30 400,13 |
118,5 |
|
III |
89 909,29 |
120,9 |
61 354,02 |
120,8 |
30 138,06 |
119,5 |
|
I |
104 807,68 |
139,2 |
70 830,31 |
137,3 |
35080,39 |
137,9 |
2004 |
II |
109 990,23 |
142,6 |
75203,95 |
141,5 |
36 322,08 |
140,3 |
|
III |
107 375,81 |
146,1 |
72 989,90 |
145,8 |
35698,54 |
142,8 |
Динамические ряды исходных данных были разделены на две части: первые пять наблюдений были использованы для построе ния начального значения адаптивного матричного предиктора, а три последних — для его обучения (настройки параметра адапта ции) по постпрогнозным расчетам. В рамках процедуры постро ения начального значения рассчитывались средние приросты (5.42) показателей за 1995—1999 гг., которые использовались при фор мировании комбинированной матрицы прямых и косвенных тем пов прироста (5.43)—(5.44)
'0,1607 |
0,2699 |
0,4220^ |
Vfl = 0,0899 |
0,1510 |
0,2361 |
0,0699 |
0,1174 |
0,1836 |
|
'0,15 |
0,85 |
0,85^ |
1 |
0,3 |
0,7 |
М |
0,85 |
0,15 |
0,85 |
W = 0,7 |
1 |
0,3 |
|
0,85 |
0,85 |
0,15 |
0,7 |
0,3 |
1 |
С помощью этих матриц было получено начальное значение предиктора (5.45)
199