З досвіду застосування активного індуктивного моделювання до прогнозування цін

Ознайомлення з циклами абсорбції (поглинання) та редукції (засвоєння). Розгляд прикладу прогнозування цін на молібден за допомогою різних методів. Характеристика різних економічних критеріїв. Аналіз проблеми нестабільності у лінійному моделюванні.

Рубрика Экономико-математическое моделирование
Вид статья
Язык украинский
Дата добавления 02.10.2018
Размер файла 91,2 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

р

Міжнародний центр інформаційних технологій та систем

З досвіду застосування активного індуктивного моделювання до прогнозування цін

УДК 519.237 + 004.896

Ю.В. Дзядик

iurius@i.com.ua

проспект академіка Глушкова, 40, Київ

Анотації

У роботі коротко висвітлюються такі теми:

(1) проблема нестабільності у лінійному моделюванні;

(2) факторний аналіз та стабілізація; метод двох порогів (МДП), або (в, г)-метод;

(3) принцип стабілізації; гіпотеза про сутність МГУА;

(4) економічний критерій моделей прогнозування цін;

(5) активні агентні моделі; цикли абсорбції (всмоктування) та редукції (засвоєння).

Наведено приклади, що демонструють перевагу МДП над відомими методами прогнозування GAME та МГУА.

Ключові слова: індуктивне моделювання, проблема стабілізації, активна агентна модель, факторний аналіз, (в, г) метод, метод двох порогів (МДП), цикли абсорбції та редукції.

This paper shortly considers the next topics:

(1) problem of unstability in linear modelling;

(2) factor analysis and stabilization; method of two thresholds (MTT), or (в, г)-method;

(3) stabilization principle; the hypothesis about essence of GMDH;

(4) economic criteria of forecast models;

(5) active agent models; cycles of absorption and reduction.

On some examples we demonstrate the advantage of MTT over known methods GAME and GMDH.

Keywords: inductive modeling, stabilization problem, active agent models, factor analysis, method of two thresholds (MTT), (в, г) method, economic criterion, cycles of absorption and reduction

В статье кратко рассмотрены следующие вопросы.

(1) Проблема нестабильности в линейном моделировании.

(2) Факторный анализ и стабилизация. Метод двух порогов (МДП), или (в, г)-метод.

(3) Принцип стабилизации. Гипотеза о сущности МГУА.

(4) Экономический критерий моделей прогнозирования цен.

(5) Активные агентные модели. Циклы поглощения и редукции.

Некоторые примеры показывают преимущество МДП перед известными методами прогнозирования GAME и МГУА.

Ключевые слова: индуктивное моделирование, проблема стабилизации, активная агентная модель, факторный анализ, (в, г) метод двух порогов (МДП), цикли абсорбции и редукции.

Вступ

Нехай ми будуємо індуктивну модель y = f (t1, t2, …, t l), яка є лінійною відносно базисних функцій {x1, x2, …, xk}, де xi = ц i (t1, t2, …, t l) - деякі функції [1, 7]. Нехай n - розмірність статистики, що використовується. Ми отримуємо вектор y та k векторів {x1, x2, …, xk} у дійсному n_вимірному просторі R n.

Позначимо

X = (x1- м(x1), x2- м(x2), …, xk - м(xk)), де wRn: м(w) = (1)

Проблема нестабільності у лінійному моделюванні

Ми шукаємо форму a лінійної залежності y = Xa, де n_вектор y та (n, k)_матриця X відомі, k_вектор коефіцієнтів a є шуканим. Відомо, що чисто символічний шлях розв'язання формального рівняння y = Xa:

y = Xa XTy = XTXa (XTX)-1XTy = a(2)

приводить до результату, тотожного з розв'язком, що дає метод найменших квадратів (МНК):

a = (XTX)-1XTy.(3)

Позначимо квадратну (k, k)_матрицю Грама [3] XTX = W.

Нехай {л1, л2,… лk} - впорядкована множина власних значень W, занумерованих так, що л1 л2 ... лk. Нагадаємо, що для будь-якої матриці Грама всі власні значення лi дійсні та невід'ємні: лi 0 i. Нагадаємо, що тоді

det W = л1 л2лk ; trace W = л1 + … + лk; cond W = л1 (л k)-1;(4)

причому c: cond cW = cond W, де cond W = ||W|| ||W-1|| - число обумовленості, || || - будь-яка “розумна” норма, c - довільна константа.

Тепер дослідимо вираз (3): a = W-1XTy. Він не має сенсу, якщо det W = 0 лk = 0 cond W = . Більш того, якщо cond W є близьким до нескінченності, тобто найменше власне значення лk є близьким до нуля, то рішення y = Xa є нічого не вартим для екстраполяції чи прогнозування [2]. Це очевидно з добре відомого частинного випадку: точна поліноміальна інтерполяція, яка є розв'язком деякого матричного рівняння вигляду y = Xa, не має сенсу за межами інтервалу інтерполяції.

Визначимо міру стабільності матриці X як

stab X = k лk Tr (XTX)-1 = k лk (л1 + л2 + … + лk)-1 . (5)

Очевидно, для довільної матриці X та константи c мають місце властивості: абсорбція ціна економічний

0 stab (X) 1; c: stab (cX) = stab (X);

stab (X) cond (XTX) = k л1 (л1 + … + лk) -1 .(6)

Тоді з (6) та нерівностей 1 (л1 + л2 + … + лk) / л1 k отримуємо:

1 stab (X) cond (XTX) k stab X [cond (XTX)] -1. (7)

Тепер ми готові підійти до головної проблеми - як уникнути безглуздості, беззмістовності, нікчемності моделі y = Xa поза межами області, у якій вона була побудована?

1. Факторний аналіз та стабілізація. Метод двох порогів (МДП), або (в, г)_метод

Приведемо матрицю Грама XTX = W ортогональним перетворенням S до такого діагонального вигляду STWS = D, що dii = лi, л1 л2 ... лk. Тоді вектори (стовпці) матриці XS = Z = (z1, z2, …, zk) називаються факторами X.

Зауважимо, що i: м(zi) = 0, отже, {i, y, c}: y + c, zi = y, zi.

Позначимо через L лінійну оболонку вектор-стовпців матриці X:

L = L (x1- м(x1), x2- м(x2), …, xk - м(xk)).(8)

За побудовою, перші ненульові фактори {z1, z2, …, zp}, де

p = dim L = rank X min (k, n),

утворюють ортогональний базис простору L. Отже, довільна лінійна модель y (x1, x2, …, xk) застосуванням перетворення може бути представлена у вигляді

y = y0 + y1 z1 + y2 z2 + … + yp zp,(9)

yi = , 0 < i p.(10)

Далі, визначимо для кожного фактора zi дві характеристики: стабільність stab (zi) та суттєвість essn (y, zi)

stab (zi) = ; (11)

essn (y, zi) = corr2 (y, zi) = ,(12)

де corr (y, zi) = cos (y - y0, zi), тобто

corr (y, zi) = .(13)

Назвемо (в, г)_редукцією моделі y за допомогою (в, г)-методу двох порогів (МДП), таку діагональну проекцію P моделі (9), diag P = (s1, s2, …, sp)

ys = y0 + s1 y1 z1 + s2 y2 z2 + … + sp yp zp,(14)

де si = si (в, г) = 0, якщо stab (zi) < в або essn (y, zi) < г, та si = 1 для всіх інших факторів. Іншими словами, (в, г)_редукція є викреслювання всіх в_нестабільних та г_несуттєвих факторів.

Якщо г = 0, то (в, 0)_редукцію будемо називати в_стабілізацією.

Очевидно, для кожного в > 0 існує таке число j N, що всі в_стабільні фактори - це просто підмножина перших j факторів {z1, z2, …, zj}.

2. Принцип стабілізації у лінійному моделюванні

Нагадаємо, що простір L визначено формулою (8). Нехай V - довільний підпростір L, P - відповідна проективна матриця, так що XP - це проекція X на підпростір V. Ми можемо вище замінити матрицю X матрицею XP та застосувати всі визначення: cond (PTWP), stab (XP) тощо.

Нехай XPS = U = (u1, u2, …, uw) - фактори матриці XP.

Визначимо стабільність, суттєвість та достатність підпростору V як

stab (V) = min {stab (v): vV} = min {stab (ui): ui U} = stab (uw);(15)

essn (y, V) = min {essn (y, v): vV} = min {essn (y, ui): ui U};(16)

suff (y, V) = .(17)

Будемо говорити, що підпростір V є:

в_стабільним, якщо stab (V) > в;

г_суттєвим, якщо essn (V) > г;

д_достатнім, якщо suff (y, V) > 1- д.

Якщо в, г чи д є фіксованими, ми будемо говорити, що підпростір V є просто стабільним, суттєвим чи достатнім. Як правило, неявні значення є: в = 10_3, г = 10_4, д = 5%.

Нехай у = {i1, i2, …, iу} - деяка підмножина у K = {1, 2, …, k}. Позначимо через V(у, X), V(у, Z) лінійну оболонку векторів-стовпців {xi - м(xi): iу} та {zi: iу} відповідно.

Гіпотеза. Сутність МГУА полягає у проектуванні простору L на такий стабільний та суттєвий підпростір V(у, X), який є найбільш достатнім.

Загальним (в, г)-методом, або принципом стабілізації, ми назвемо пошук розв'язку наступної оптимізаційної задачі. Нехай задано деяку множину 2L, де 2L є множина усіх лінійних підпросторів L, а також пороги в та г. Знайти у множині такий в_стабільний та г_суттєвий підпростір V, який є найбільш достатнім:

{stab (V) > в} & {essn (y, V) > г} & {suff (y, V) = max} & {V }(18)

Приклади множин :

(a) = {V(у, X): у2K };

(b) = {V(у, Z): у2K };

(c) = 2L.

3. Активні агентні моделі. цикли абсорбції (поглинання) та редукції (засвоєння)

Якщо ми не знайшли достатньо показників { t1, t2, …, t l }, від яких залежить прогнозована величина y, то ніякий метод не дасть задовільного прогнозу. Отже, (в, г)-метод двох порогів (МДП) лише відсіє нестабільні та несуттєві (зайві) змінні і забезпечить, що модель y = Xa буде стабільною. Але точність моделі не може бути більшою, ніж максимальна достатність серед стабільних підпросторів у 2L. Другою частиною циклу є абсорбція (поглинання) релевантних змінних (показників, в термінах економіки) у інформаційному просторі (як правило, у базах даних в Інтернет). Ця частина індуктивного моделювання має справу з такими термінами, як побудова інтелектуальних агентів (constructing intelligent agents, CIA) [8], штучний інтелект, мультиагентні системи (МАС) тощо. Більшість з цих тем ще не входять до теорії індуктивного моделювання. Дамо кілька початкових визначень.

Активна модель - модель, яка шукає усі релевантні показники (t1, t2, …, t l), що необхідні для її функціонування (та статистику цих показників) у інформаційному просторі (як правило, в Інтернет).

Активна агентна модель - модель, яка провадить пошук у Інтернет за допомогою інтелектуальних пошукових агентів. Синонім: активна інтелектуальна модель.

Активне прогнозування є прогнозування за допомогою активної інтелектуальної моделі.

4. Деякі результати. Порівняння з іншими методами

Метод двох порогів був успішно застосований до моделювання та прогнозування цін на феромолібден (FeMo) та молібден (Mo). Особливий інтерес викликало прогнозування місячних цін у 2004-07, коли вони надзвичайно швидко змінювались. У 2004 році деякі часописи, які багато десятиліть надавали інформацію про FeMo та Mo ціни, перестали публікувати цю інформацію. Навіть можливість будь-якого прогнозування цін при таких початкових даних ставилася під сумнів.

Нами для місячного прогнозування цін на молібден у 2005-07, у Інтернет (після чотирьох циклів абсорбції та редукції) було вибрано 9 показників:

t1, t2 - ціни на експорт та імпорт молібдену, обчислені на підставі даних у таблицях [4], стовпці “Exports” та “Imports”, ряд “Molybdenum”,

t3 - ціни на імпорт міді, аналогічно t1, t2 обчислені з таблиць [4]

t4, t5, …, t 9 - ціни на сталь, взяті з [5].

Табл. 1. Прогнозування цін на молібден ($/кг) у 2005-07 рр. за допомогою різних методів: МДП (стовпець 4), МГУА критерій Акаіке (5), МГУА регулярний критерій (6), прогноз зміни ціни за МДП (7, 8), та рекомендації плану закупок за прогнозом МДП: на скільки місяців (9) та на яку суму (10)

T

місяць

факт

МДП

FPE Akaike

AR regular

прогноз змін

план

abs

%

міс.

$

1

2

3

4

5

6

7

8

9

10

24

Oct 05

93,51

-23,15

< 0

25

Nov 05

67,45

70,36

87,46

80,81

6,43

10%

2

134,90

26

Dec 05

81,06

73,88

79,16

79,42

-18,66

< 0

27

Jan 06

80,58

62,40

57,21

60,11

-15,08

< 0

28

Feb 06

76,19

65,50

47,80

67,84

0,90

< 0

29

Mar 06

73,15

77,08

69,72

63,21

0,22

< 0

30

Apr 06

68,49

73,37

107,15

76,02

7,01

10%

1

68,49

31

May 06

73,38

75,50

65,89

67,48

3,17

4%

1

73,38

32

Jun 06

66,83

76,55

44,19

63,08

10,63

16%

2

133,66

33

Jul 06

67,55

77,46

34,03

74,84

12,18

18%

2

135,10

34

Aug 06

66,29

79,73

73,04

67,11

3,32

5%

1

66,29

35

Sep 06

68,37

69,60

61,42

68,09

16,86

25%

3

205,11

36

Oct 06

86,76

85,23

100,67

86,95

-3,23

< 0

37

Nov 06

71,01

83,53

97,35

98,46

16,03

23%

2

142,02

38

Dec 06

76,35

87,05

92,08

87,31

8,67

11%

1

76,35

39

Jan 07

68,18

85,02

86,49

82,16

15,84

23%

2

136,36

40

Feb 07

84,67

84,02

79,65

62,20

СКВ

9,62

20,22

12,51

сума

1269,82

17

1171,66

З показників {t1, t2, …, t 9} ми утворюємо залежну змінну y(ф) та k = 12 вхідних змінних {x1, x2, …, x12}:

y (ф) = t1(ф+1);

x1(ф) = t1(ф), x2(ф) = t1(ф-1), x3(ф) = t1(ф-2),

x4(ф) = t2(ф), x5(ф) = t2(ф-1),

xi(ф) = ti-3(ф), i = 6..12.

Для отримання прогнозів використовувалося декілька алгоритмів. Для усіх алгоритмів спільним є прогнозування методом ковзного вікна на місяць наперед: ми вибираємо значення вхідних змінних {x1(t), x2(t), …, xk(t)} для M послідовних значень часу t = q+1, q+2,…, q+M, та обраховуємо прогнозоване значення величини y(t) у (M+1)-й точці t = q+M+1. Величина q приймає значення, виходячи з умов

3 ? q+1 ? t ? q+M+1 ? 39, звідки отримуємо 2 ? q ? 38-M .

У більшості випадків ми покладали M = 21, але для порівняння якості прогнозу вибирали й інші значення.

Отримуємо наступні результати (табл. 1 та рис. 1).

Рис.1.Результати моделювання цін на молібден

Інший приклад базується на медичних даних з Motol hospital у Празі. Ці дані надав P. Kordik для порівняння (в, г)-методу з результатами [6], які були отримані за допомогою методу GAME (Group of Adaptive Models Evolution).

Табл. 2. Порівняння (в, г)-методу (стовпець MTT) з МГУА (стовпці GMDH) та GAME на прикладах прогнозування цін на молібден (ряд Mo) у 2005-07 та прогнозу вмісту CO2 у мозку (ряд CO2): СКВ помилки передбачення

example\method

MTT

GMDH, FPE

GMDH, AR

GMDH, AC

GAME

CO2

0,0380

--

--

0,0704

0,0386

Mo

9,62

20,22

12,51

--

--

Дані з табл. 2 та порівняння графіків на рис. 2 демонструють перевагу МДП перед МГУА та GAME [6]. Варто відзначити, що GAME вимагає велетенську навчальну вибірку (декілька сот спостережень) [6], водночас для МДП (як і для МГУА) досить 10_20 спостережень. У таблиці 2 наведено середньоквадратичне відхилення помилки передбачення.

Рис. 2. Графіки фактичних цін на молібден (лінія “fact”) та прогнозних значень, обчислених за МДП (лінія “Praha”) та МГУА (лінії “FPE” та “AR”)

5. Економічні критерії

Нехай y - ціна деякого товару (чи, більш загально, блага). Можна купувати цей товар щомісяця за ринковою ціною yi, тоді видатки будуть дорівнювати Z0 = ni=1 yi. Наприкінці періоду можна визначити найменшу ціну ymin = min ni=1 yi. Очевидно, видатки на придбання товару не могли бути меншими за Zmin = n ymin. Якщо ми будуємо прогноз за методом m (method), ми можемо за деяким правилом w (way) використати цей прогноз для планування закупок. Тоді ми можемо підрахувати відповідні видатки Zm,w :

Zm,w = ,(19)

де pi (m, w) - кількість товару, придбаного у i_му місяці, застосовуючи метод прогнозування m та спосіб використання прогнозу для планування закупок w.

Можна запропонувати кілька різних економічних критеріїв. Наприклад, відносний критерій

.(20)

Покажемо, що для рекомендацій закупок молібдену згідно МДП, наведених у стовпцях 9-10 табл. 1, цей критерій дорівнює біля 68%. Проведемо разрахунки.

pmin = 66,29, отже, Zmin = n pmin = 17 (міс) 66,29 ($/кг/міс) = 1126,93 ($/кг).

Z0 = 1269,82 ($/кг), Zm,w = 1171,66 ($/кг)

= 68.7%

Нехай деякий завод щомісяця використовує Q = 20 тон молібдену. Якщо він застосовує прогноз цін за (в, г)-методом (стовпець 4 таблиці 1) та відповідні рекомендації (стовпці 9, 10) для прийняття рішення про величину закупок для формування стратегічного запасу, то за період з жовтня 2005 по лютий 2007 (17 місяців) абсолютна економія цього заводу порівняно із щомісячними закупками за ринковими цінами (стовпець 3) склала б

(Z0 - Zm,w) Q = (1269,82 - 1171,66) $/кг 20 000 кг = 98,16 20 000 $,

що дорівнює 1 963 200 $, тобто майже два мільйони доларів.

Посилання

1. Ивахненко А. Г., Ивахненко Г. А. Обзор задач, решаемых по алгоритмам Метода Группового Учета Аргументов (МГУА) - http://www.gmdh.net/articles/rus/algorith.zip/algorith.doc, http://www.gmdh.net/articles/rus/algorith.pdf

2. Коппа Ю. В., Степашко В. С. Сравнение прогнозирующих свойств моделей регрессионного типа и МГУА - http://www.gmdh.net/articles/rus/compare.pdf

3. Gramian matrix - http://www.answers.com/topic/gramian-matrix

4. Metals Statistics, U. S. Metals Trade - http://www.ita.doc.gov/td/metals/statindx.html

5. MEPS (International) Ltd. - Independent Steel Industry Analysts, Consultants, Steel Prices, Reports and Publications. World Stainless Steel Product Prices - http://www.meps.co.uk/Stainless Prices.htm

6. Josef Bouљka, Pavel Kordik. Time Series Prediction by means of GMDH Analogues Complexing and GAME - http://www.gmdh.net/articles/iwim/iwim_39.pdf

7. Yuriy V. Dzyadyk. Some Results of the Synthesis of GMDH and Factor Analysis for Inductive Modelling - http://www.gmdh.net/articles/iwim/iwim_19.pdf

8. Joseph P. Bigus, Jennifer Bigus. Constructing Intelligent Agents Using Java™: Professional Developer's Guide. John Wiley & Sons, 2001, second edition. ISBN 0471191353 - http://www.research.ibm.com/people/b/bigus

9. Тимашова Л. А., Дзядик Ю. В., Лещенко В. А., Бондар Л. А. Інтелектуальна система прогнозування цін // Проблеми впровадження інформаційних технологій в економіці. Тези доповідей VI Міжнародної науково-практичної конференції. Ірпінь, 2007 - http://www.i.com.ua/~iurius/dzyadyk/sipro.htm

Размещено на Allbest.ru

...

Подобные документы

  • Методи економічного прогнозування, їх відмінні особливості, оцінка переваг та недоліків. Моделі прогнозування соціально-економічних об’єктів. Принципи вибору моделей та комбінування прогнозів. Прогнозування показників розвитку банківської системи.

    курсовая работа [813,1 K], добавлен 18.02.2011

  • Економіко-математичне моделювання як спосіб вивчення господарської діяльності. Аналіз коефіцієнтів оборотності капіталу. Оцінка факторів, що впливають на ділову активність. Застосування моделей прогнозування для підприємств гірничообробної промисловості.

    курсовая работа [274,5 K], добавлен 06.09.2013

  • Структурна схема ВАТ "Вагоно-ремонтний завод". Аналіз фінансового та економічного стану підприємства. Методики побудови апроксимаційних нелінійних залежностей за допомогою методу Ньютона нелінійного оптимального пошуку. Розробка методики прогнозування.

    дипломная работа [986,3 K], добавлен 08.03.2010

  • Вихідні поняття прогнозування, його сутність, принципи, предмет і об'єкт. Суть адаптивних методів. Прогнозування економічної динаміки на основі трендових моделей. Побудова адаптивної моделі прогнозування прибутку на прикладі стоматологічної поліклініки.

    дипломная работа [1,2 M], добавлен 18.06.2015

  • Введення в міжнародний валютний ринок FOREX, проблема прогнозованості, аналіз математичних методів. Формалізація задачі прогнозування валютних курсів на основі теорії нечітких множин, оцінка адекватності результатів на основі запропонованого методу.

    дипломная работа [985,4 K], добавлен 12.06.2013

  • Економетричні моделі - системи взаємопов'язаних рівнянь і використовуються для кількісних оцінок параметрів економічних процесів та явищ. Прикладні економетричні моделі Франції та США. Макроеконометричні моделі України та прогнозування економіки.

    реферат [20,6 K], добавлен 01.02.2009

  • Теоретичні основи економічного прогнозування: сутність, види і призначення, принципи і методи. Особливості вибору моделей та створення систем державних прогнозів і соціально-економічних програм України. Порядок моделювання динаміки господарської системи.

    курсовая работа [869,6 K], добавлен 16.02.2011

  • Використання методів економетричного моделювання, аналізу і прогнозування на всіх напрямках економічних досліджень: мікро- та макроекономіка, міжнародна економіка, фінансові ринки. Розробка і використання адекватних статистичних (економетричних) моделей.

    контрольная работа [330,4 K], добавлен 25.01.2015

  • Теоретичні аспекти дослідження ID-IS моделей. Попит та пропозиція як економічні категорії. Особливості моделей перехідної економіки. Аналіз підходів щодо моделювання сукупного попиту та пропозиції. Процес досягнення рівноваги та прогнозування ціни.

    курсовая работа [639,7 K], добавлен 15.11.2010

  • Сутність прогнозу та прогнозування. Теоретичні основи наукового передбачення. Класифікація прогнозів і прогнозування за періодичністю проведення та ступенем вірогідності, за формами конкретизації управління. Аналіз процесів і тенденцій у сучасному світі.

    реферат [34,5 K], добавлен 09.12.2013

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.