CC BY
57
УДК 330.45
08.00.00 Экономические науки
ПРОГНОСТИЧЕСКОЕ ИССЛЕДОВАНИЕ
ПРИРОДНО-ЭКОНОМИЧЕСКОГО
ПРОЦЕССА
Кумратова Альфира Менлигуловна к.э.н., доцент SPIN-код=2144-8802
Василенко Игорь Иванович к.с-х.н., доцент SPIN-код: 1828-1453
Ксёнз Станислав Юрьевич магистрант
ФГБОУ ВПО «Кубанский государственный аграрный университет» г. Краснодар, Россия
Ратушная Елена Анатольевна к. э. н.
Филиал Ростовского государственного экономического университета, Черкесск, Россия
Увеличение объема обрабатываемых данных и бурное развитие мониторинга окружающей среды, моделирования, предвидения, анализа, визуализации, прогнозирования в современных условиях связано с последовательным ростом уровня их формализации. Базисом этому всему послужили требования существенно изменившейся стохастики природных и экономических процессов. В статье предлагается один из новых методов нелинейной динамики, а именно метод последовательного R/S-aнaлизa. В статье авторами уделено методу фрактального анализа временных рядов. Основоположником фрактального анализа является британский гидролог Х.Е. Херст. Он показал, что природные и естественные явления, такие как стоки рек, осадки, температуры, солнечная активность следуют «смещенному случайному блужданию», т. е. тренду с шумом. Уровень шума и устойчивость тренда оцениваются изменением нормированного размаха уровней временного ряда по истечению времени, или, другими словами, насколько введенная им величина, называемая показателем Херста, превосходит значение 0,5. Для составления прогноза весьма существенную информацию носит циклическая компонента. Тем самым, возникает необходимость дальнейшего исследования природно-экономических процессов на базе новых математических моделей. Эти методы привносят в прогноз новые полезные методологические элементы, отсутствующие в непрерывной методологии, такие понятия, как: «цвет шума», персистентность и антиперсистентность рядов,
UDC 330.47 Economic science
PROGNOSTIC RESEARCH ON THE NATURAL AND ECONOMIC PROCESSES
Kumratova Alfira Menligulovna Cand. Econ. Sci., Assistant professor SPIN-codе=2144-8802
Vasilenko Igor Ivanovich Cand.Agr.Sci., Assistant professor SPIN-codе= 1828-1453
Ksenz Stanislav Yurievich Master of Science
Federal state budget institution of higher professional education "Kuban state agrarian University", Krasnodar, Russia
Ratushnaya Elena Anatolievna Cand.Econ.Sci.
Branch of Rostov State Economic University, Cherkessk, Russia
The increase in volume of processed data and the rapid development of environmental monitoring, modeling, forecasting, analysis, visualization, prediction in modern conditions is connected with the consistent increase in their level of formalization. The bases for all this are requirements of significantly changed stochastics natural and economic processes. A new method of nonlinear dynamics, namely the method of sequential R/S-analysis is proposed. In the article, the authors paid attention to the method of fractal analysis of time series. The founder of fractal analysis is a British hydrologist H.E Hurst. He showed that natural phenomena such as river flows, rainfall, temperature, solar activity is followed by «biased random walk», i.e. trend with noise. The noise level and trend resistance are estimated in change in the normalized amplitude levels of the time series for the expiration time, or, in other words, how they entered a quantity called the Hurst exponent exceeds the value of 0.5. Rather essential information is a cyclical component to forecast. Thus, there is a need for further study of natural and economic processes based on the new mathematical models. These methods bring to forecast new useful methodological elements that are not in continuous methodology, concepts such as «noise color» persistence and anti-persistent series, Hurst, «long-term memory», R/S-trajectory and the trajectory of the Hurst exponent, etc.
показатель Херста, «долговременная память», R/S-траектории и траектории показателя Херста и пр
Ключевые слова: ВРЕМЕННОЙ РЯД, Keywords: TIME SERIES, PRE-FORECASTED
ПРЕДПРОГНОЗНЫЙ АНАЛИЗ, ANALYSIS, SERIAL R/S-ANALYSIS, HURST
ПОСЛЕДОВАТЕЛЬНЫЙ R/S-АНАЛИЗ, ПОКАЗАТЕЛЬ ХЕРСТА
Прогностическое исследование природно-экономического процесса начинается в первую очередь с разработки задания на прогноз, т.е. предпрогнозная ориентация. Другими словами, создание такого документа, в котором определены предмет и объект прогнозирования, порядок разработки, задачи и цели прогнозирования, временной интервал ретроспекции и дальности прогнозных оценок, что является актуальной задачей во многих сферах экономики. В связи с этим уже в самом задании на прогноз содержится основание для разработки прогноза.
Развитие мониторинга окружающей среды, моделирования, предвидения, анализа, визуализации, прогнозирования в современных условиях связано с последовательным ростом уровня их формализации, а также роста объема обрабатываемых данных. Базисом этому всему послужили требования существенно изменившейся стохастики природных и экономических процессов. В статье предлагается один из новых методов нелинейной динамики, а именно метод последовательного R/S-анализа.
Актуальность статьи продиктована следующим обстоятельством. Как отмечено в книге математика-экономиста Э. Петерса [1313], современная экономическая теория вступила в новую фазу своего развития. Это обусловлено несколькими, факторами. В первую очередь, усложнением и глобализацией мировой экономики, а также, вторжением в науку математических методов нелинейной динамики. Вследствие этого, появлением новых компьютерных и информационных технологий, сделавших возможным исследование сложных процессов и явлений на экране дисплея.
Адаптированный авторами [19] метод последовательного R/S-анализа зарегистрирован в Российском агентстве по патентам и товарным знакам (свидетельство № 2003611093).
Актуальность предмета настоящего исследования приобрела особо важное экономическое значение после наводнения, имевшего место по р. Кубань в июне 2002 г. Материальный ущерб этого наводнения оценивается суммой свыше 0,5 млрд. долларов.
Отметим, что в результате наводнения на юге России пострадало свыше 240 населенных пунктов, разрушено более 110 километров газопровода, 269 мостов, 1490 километров автодорог, погибло 102 человека. Общее же число пострадавших в Южном федеральном округе достигло 340 тысяч человек, а материальный ущерб превысил 15 миллиардов рублей. Огромные потери претерпело и сельское хозяйство, ущерб которого составляет около 250 миллионов рублей.
В статье авторами уделено методу фрактального анализа временных рядов (ВР) [6,13]. Отметим, что основателем фрактального анализа является британский гидролог Х.Е. Херст. Он показал, что природные и естественные явления, такие как речные стоки, температуры, осадки, числа Вольфа солнечной активности, следуют «смещенному случайному блужданию», т.е. тренду с шумом. Такие показатели, как устойчивость тренда и уровень шума могут быть оценены тем, как изменяется нормированный размах исследуемого временного ряда со временем, или, другими словами, насколько введенная им величина H е (0; 1),
называемая показателем Херста, превосходит значение 0,5.
Если уровни временного ряда отражают «чисто случайный процесс» (являются независимыми случайными величинами), то в соответствии с классической статистикой для такого ВР значения показателя Херста H»0,5. Исследованные Херстом, а позже Мандельбротом и др. многочисленными природные временные ряды (выпадение осадков, пятна
на солнце, годичные кольца и т.д.) обладают так называемой долговременной памятью [13], в силу чего показатель Херста для каждого из этих ВР принимал значение, превосходящее число 0,7.
Отметим, что к моменту выхода в свет монографий [12,13] уже сложилось общепризнанное мнение о том, что для подавляющего большинства природных систем величина Н значительно отклоняется от 0,5 в область черного шума [13].
В англоязычной литературе исследуемый метод называется как «метод нормированного размаха Херста» (в источнике [14] он назван «Е/Б-анализ: Руководство шаг за шагом»)
В течение более чем полувекового периода Е/Б-анализ природных, экономических и др. ВР исследователи осуществляют на базе опубликованного в работе [18] алгоритма, который.
В качестве модельного временного ряда рассмотрим ряд ежемесячных объемов стока горной реки Кубань (материалы по данным Карачаево-Черкесской Республики) за период с 1926 г. по 2003 г. включительно. Графическое представление части исходного ВР дано на рисунке 1.
1000
800
600
400
200
спспспспспспспспспспспс^спс^спс^с^с^с^с^с^с^с^с^с^с^с^с^с^с^
Рисунок 1 - Графическое изображение ежемесячных объемов стоков р. Кубань
за период 1988 - 2003 гг.
Из рисунка 1 видно, что рассматриваемому ряду присуще свойство
0
010000020100010101020201010102000001010248
периодичности.
Для верификации метода НР Херста на базе ВР, представленного на рисунке 1, строим чисто периодический ВР:
г = , (1)
70=( 2 0
который представляет собой последовательность одного и того же цикла, представленного на рисунке 2. В дальнейшем ВР 70 называем термином «эталонный ВР».
5 4 3 2 1 0
1 2 3 4 5 6 7
Рисунок 2 - Графическое изображение цикла эталонного временного ряда (по оси х - временной интервал, по оси у значение измерения во времени 2? = 1, 20 = 2, 230 = 3, 2 40 = 4, 250 = 3, 260 = 2, 270 = 1)
С учетом длины цикла (т.е. количества составляющих его уровней) Ь = 7 в качестве длины эталонного ВР (1) принимаем значения N = 700, т.е. ВР (1) состоит из 100 следующих один за другим циклов длины Ь = 7 вышеуказанного вида.
Далее, применим к эталонному ВР (1) алгоритм нормированного размаха Херста для последовательности длин отрезков:
щ = 10, П2 = 11,..., пк = к + 9,..., щ341 = 350
(2)
В результате реализации описанных этапов в [13] получим множество точек, определяемых логарифмическими координатами вида:
Хк = 1св(Ик), Ук = 1о§((%)к), к = 1,2,...,т . (3)
Геометрическое расположение этих точек, а также линейного тренда для них представлено на рисунке 3, где для удобства видуализации
множество (3) представлено в виде Л/^-траектории, которая в свою очередь получается путем соединения отрезков каждой пары соседних точек
{хк, Ук), {хк+ъ Ук+\), к =1 т, т = 341.
= 0,0183х + 0,5785
0,8 0,7
СО
^ 0,6
га 1
о
0,5 0,4
1 1,2 1,4 1,6 1,8 2 2,2 2,4 2,6 2,8
|од(п)
Рисунок 3 - Линейный тренд эталонного временного ряда (1) и ее Я/Б-траектория
Визуализация рисунка 3 дает основание для следующего вывода: полученная линия регрессии практически параллельна оси абсцисс, т.е. угол наклона прямой регрессии у = ах + Ь близок к 0. В силу этого получаем низкую оценку показателя Херста Н, значение которой находится в зоне розового шума [13], т.е. существенно меньше 0,5. Последнее является противоречием теории ^/^-анализа [13,14], согласно которой для периодических ВР показатель Херста находится в области черного шума, т.е., как указанно в [12], его значения обычно превосходят 0,7 .
1200 1000 800 600 400 200
* >
* >
* >
* >
* >
* >
* >
* >
Рисунок 4 - Гистограмма временного ряда 2 (ежемесячных объемов стоков р. Кубань)
за период с 01.1988 г.- 12.2003 г.
0
В работе с достаточной полнотой реализованы численные расчеты с целью обосновать верхнюю оценку глубины памяти рассматриваемых ВР (алгоритм а1) [13]. Работа алгоритма а1 подробно описана в работах [6,13].
Результат работы алгоритма а1 исследуемого временного ряда 7 представлен в таблице 1.
Представленные в последней строке таблицы 1 значения /(/) функции принадлежности элементов / нечеткому множеству М(7) пропорциональны числам ё(/), / е ь(7), они получаются путем нормирования значений долей ё(/) так, что /(/)<1 для всякого /е ь(Б). Результат работы этапа 2 для ВР 7 представлен в таблице 1. Значения элементов /(/) последней строки в таблице 1 вычисляются следующим образом. Сначала находим максимальную долю ё * = тах ё (/) (например, в
1еь(г)
таблице 1 значение ё * =0,17) и соответствующую ей глубину / * (ё *(/ ) = / *, в таблице 1 значение / *=16).
Таблица 1 - Глубина памяти для ВР ежемесячных объемов стока р. Кубань за период с 01.1988 г. по 12.2003 г.
Глубина памяти 12 13 14 15 16 17 18 19 20 21
Количество глубины памяти 1 6 12 16 17 17 9 9 8 5
Доля глубины памяти 0,01 0,06 0,12 0,16 0,17 0,17 0,09 0,09 0,08 0,05
Функция принадлежности 0,05 0,3 0,64 0,85 0,9 0,9 0,47 0,47 0,42 0,26
глубины памяти
Далее для этой глубины /* экспертным путем устанавливается значение функции принадлежности /и* = /и(/*) (в таблице 1 значение. т'=м(16) =0,9). После чего для остальных элементов /е ь(7)
соответствующие им значения функции принадлежности /(/)
*
вычисляются по формуле /(/ ) = Ц ё (/). Формирование НМ М (7)
ё *
осуществляется путем попарного объединения элементов первой и
последней строк таблицы 1. Полученная оценка глубины памяти для ВР 7 представляется в виде следующего нечеткого множества:
М (7) = {(12,0,05), (13;0,3), (14;0,64), (15;0,85), (16;0,9), (17;0,9), (5)
(18;0,47), (19;0,47), (20;0,42), (21;0,2б)} (5)
Геометрическое изображение нечеткого множества рассматриваемого ВР 7 представлено на рисунке 6.
Ч 1
0,9 0,8 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0
12 13 14 15 16 17 18 19 20 21
Рисунок 6 - Графическое представление нечеткого множества глубины памяти для временного ряда 7
Исходя из выполненных расчетов получены следующие выводы:
1. Глубина памяти исследуемого временного ряда ежемесячных объемов стока р. Кубань не является постоянным числом. Из таблицы 1 видно, что численные значения глубины памяти меняются в диапазоне натурального ряда 12, 14,..., 21.
2. Для численного представления глубины памяти рассматриваемого ВР (4) 7 наиболее целесообразным является математический аппарат теории нечетких множеств, т.е. оцениваемая глубина представляет собой нечеткое множество М(7) = {(/,/(/))}, /е {/0,/0 +1,...,Ь0 }, где / - численное значение встречающейся глубины памяти, /(/) - значение функции принадлежности для этой глубины.
22 21 20 19 18 17 16 15 14 13 12 11 10
12 месяцев
00°° оо
co°°gco
оооооого^ооооооа>°>та>а>а>а>°>та>а>а>а>°>та>а>а>а>°>у
дата
Рисунок 7 - Графическое представление номеров смены тренда ^/¿'-траекторий временных рядов 2Г из семейства £(Z)
Как отмечено в [14], .R/S-анализ представляет собой инструмент для выявления циклов, как периодических так и непериодических. На рисунке 7 представлено изображение в виде гистограммы динамики точек смены тренда Ä/S-траекторий временных рядов Zr е S(Z) для индексов r, поставленных в соответствие месяцам, начиная с января 1988 по апрель 1996 г. Исходя из этого, можно выносить суждение о новых дополнительных «циклических» возможностях «алгоритма последовательного Ä/S-анализа» по отношению к «алгоритму нормированного размаха Херста».
Примечание. Для изложения выводов, вытекающих из визуализации рисунка 6, необходимо использовать понятие «цикл». На рисунке 7 представлено графическое представление ВР L = {¡г), i = 1,2,..., N.
Согласно определению, под термином «цикл» подразумевается отрезок:
h, W., , (6)
состоящий из двух частей, первая из которых является возрастающей (убывающей), а вторая - убывающей (возрастающей). При этом указанный отрезок (6) является локально-максимальным по своей длине (r +1).
Обнаружение долговременной памяти в ежемесячных объемах стоков горных рек является не самоцелью, а должно послужить объективным обоснованием принципиальной возможности построения прогнозной модели, в процессе работы которой учитываются все существенные факторы (например, наличие циклов вида (6)), которыми обусловлено наличие этой памяти. В контексте проблемы прогнозирования уместно отметить уже сложившееся, т.е. ставшее классическим основное положение декомпозиционного анализа временных рядов [17]. В общем случае временной ряд условно разделяется на следующие составляющие части: тренд; циклическая компонента; сезонное колебание; нерегулярная или остаточная компонента.
При этом циклическая компонента может нести весьма существенную информацию для составления прогноза. Тем самым, возникает необходимость дальнейшего исследования природно-экономических процессов на базе новых математических моделей. Эти методы привносят в прогноз новые полезные методологические элементы, отсутствующие в непрерывной методологии, такие понятия, как: «цвет шума», персистентность и антиперсисистентность рядов, показатель Херста, «долговременная память», Я/8-траектории и траектории показателя Херста и пр.
Литература
1. Кумратова А. М. Выявление свойств прогнозируемости методами классической статистики / А. М. Кумратова // В сборнике: Актуальные проблемы социально-экономических исследований сборник материалов 6-й Междунар. научно-практ. конф. НИЦ «Апробация». - 2014. - С. 99-101.
2. Кумратова А. М. Исследование тренд-сезонных процессов методами классической статистики / А.М. Кумратова // Политематический сетевой электронный научный журнал КубГАУ. - 2014. - № 103. - С. 312-323.
3. Кумратова А. М. Методы искусственного интеллекта для принятия решений и прогнозирования поведения динамических систем / А. М. Кумратова // Политематический сетевой электронный научный журнал КубГАУ. - 2014. - № 103. -С. 324-341.
4. Кумратова А. М. Методы многокритериальной оптимизации и классической статистики для оценки риск-эктремальных значений / А. М. Кумратова, Е. В. Попова, Н. В. Третьякова // Известия КубГУ. Естественные науки. 2014. № 1. С. 55-60.
5. Кумратова А. М. Методы нелинейной динамики как основа построения двухуровневой модели прогноза / А. М. Кумратова // В сборнике: Экономическое прогнозирование: модели и методы материалы X международной научно-практической конференции. Воронеж, 2014. С. 169-174.
6. Кумратова А. М. Оценка и управление рисками: анализ временных рядов методами нелинейной динамики: монография / А. М. Кумратова, Е. В. Попова. -Краснодар: КубГАУ, 2014. - 212 с.
7. Кумратова А. М. Прогноз динамики экономических систем: клеточный автомат / А. М. Кумратова. - Краснодар: КубГАУ, 2015. - 241 с.
8. Кумратова А. М. Прогнозирование и выявление сезонных компонент временного ряда туристского потоками /А. М. Кумратова, Е. В. Попова, М. И. Попова // В сборнике: Актуальные проблемы социально-экономических исследований сборник материалов 6-й Международной научно-практической конференции. НИЦ «Апробация». 2014. С. 89-98.
9. Кумратова А. М. Сопоставительный анализ прогноза урожайности для зон рискового земледелия / А. М. Кумратова // Экономическое прогнозирование: модели и методы: Материалы X междунар. научно-практ. конф. Воронеж. - 2014. - С. 174-179.
10. Кумратова А. М. Точный прогноз как эффективный способ снижения экономического риска агропромышленного комплекса / А. М. Кумратова // Политематический сетевой электронный научный журнал КубГАУ. - 2014. - № 103. -С. 293-311.
11. Кумратова А. М. Экономико-математическое моделирование риска в задачах управления ресурсами здравоохранения / А. М. Кумратова, Е. В. Попова, А. З. Биджиев. - Краснодар: КубГАУ, 2014. - 168 с.
12. Перепелица В.А., Попова Е.В. Математические модели и методы оценки рисков экономических, социальных и аграрных процессов. - Ростов н/Д: Изд-во Рост. Ун-та, 2002. -202с.
13. Петерс Э. Хаос и порядок на рынках капитала. Новый аналитический взгляд на циклы, цены и изменчивость рынка / Э. Петерс. - М.: Мир, 2000. - 333 с.
14. Петерс Э. Фрактальный анализ финансовых рынков: Применение теории Хаоса в инвестициях и экономике. - М.: Интернет-трейдинг, 2004. - 304с.
15. Попова Е. В. Методы моделирования поведения экономических систем на основе анализа временных рядов / Е. В. Попова, А. М. Кумратова, М. И. Попова // Экономическое прогнозирование: модели и методы. Материалы X междунар. научно-практ. конф. - Воронеж. - 2014. - С. 200-206.
16. Попова Е. В. Управление рисками в вопросах безопасности инвестиций в АПК / Е. В. Попова, А. М. Кумратова // Экономическое прогнозирование: модели и методы. Материалы X междунар. научно-практ. конф. - Воронеж, - 2014. - С. 194-200.
17. Сигэл Э. Практическая бизнес-статистика. - М.: Издательский дом «Вильямс», 2004. -1056с.
18. Федер Е. Фракталы. - М.: Мир, 1991. - 260 с.
19. Янгишиева А. М. Моделирование экономических рисков методами нелинейной динамики. Автореферат диссертации на соискание ученой степени кандидата экономических наук / Ставропольский государственный университет. Ставрополь, 2005. - 24 с.
References
1. Kumratova A. M. Vyjavlenie svojstv prognoziruemosti metodami klassicheskoj statistiki / A. M. Kumratova // V sbornike: Aktual'nye problemy social'no-jekonomicheskih issledovanij sbornik materialov 6-j Mezhdunar. nauchno-prakt. konf. NIC «Aprobacija». -2014. - S. 99-101.
2. Kumratova A. M. Issledovanie trend-sezonnyh processov metodami klassicheskoj statistiki / A.M. Kumratova // Politematicheskij setevoj jelektronnyj nauchnyj zhurnal KubGAU. - 2014. - № 103. - S. 312-323.
3. Kumratova A. M. Metody iskusstvennogo intellekta dlja prinjatija reshenij i prognozirovanija povedenija dinamicheskih sistem / A. M. Kumratova // Politematicheskij setevoj jelektronnyj nauchnyj zhurnal KubGAU. - 2014. - № 103. - S. 324-341.
4. Kumratova A. M. Metody mnogokriterial'noj optimizacii i klassicheskoj statistiki dlja ocenki risk-jektremal'nyh znachenij / A. M. Kumratova, E. V. Popova, N. V. Tret'jakova // Izvestija KubGU. Estestvennye nauki. 2014. № 1. S. 55-60.
5. Kumratova A. M. Metody nelinejnoj dinamiki kak osnova postroenija dvuhurovnevoj modeli prognoza / A. M. Kumratova // V sbornike: Jekonomicheskoe prognozirovanie: modeli i metody materialy X mezhdunarodnoj nauchno-prakticheskoj konferencii. Voronezh, 2014. S. 169-174.
6. Kumratova A. M. Ocenka i upravlenie riskami: analiz vremennyh rjadov metodami nelinejnoj dinamiki: monografija / A. M. Kumratova, E. V. Popova. - Krasnodar: KubGAU, 2014. - 212 s.
7. Kumratova A. M. Prognoz dinamiki jekonomicheskih sistem: kletochnyj avtomat / A. M. Kumratova. - Krasnodar: KubGAU, 2015. - 241 s.
8. Kumratova A. M. Prognozirovanie i vyjavlenie sezonnyh komponent vremennogo rjada turistskogo potokami /A. M. Kumratova, E. V. Popova, M. I. Popova // V sbornike: Aktual'nye problemy social'no-jekonomicheskih issledovanij sbornik materialov 6-j Mezhdunarodnoj nauchno-prakticheskoj konferencii. NIC «Aprobacija». 2014. S. 89-98.
9. Kumratova A. M. Sopostavitel'nyj analiz prognoza urozhajnosti dlja zon riskovogo zemledelija / A. M. Kumratova // Jekonomicheskoe prognozirovanie: modeli i metody: Materialy X mezhdunar. nauchno-prakt. konf. Voronezh. - 2014. - S. 174-179.
10. Kumratova A. M. Tochnyj prognoz kak jeffektivnyj sposob snizhenija jekonomicheskogo riska agropromyshlennogo kompleksa / A.M. Kumratova // Politematicheskij setevoj jelektronnyj nauchnyj zhurnal KubGAU. - 2014. - № 103. - S. 293311.
11. Kumratova A. M. Jekonomiko-matematicheskoe modelirovanie riska v zadachah upravlenija resursami zdravoohranenija / A. M. Kumratova, E. V. Popova, A. Z. Bidzhiev. -Krasnodar: KubGAU, 2014. - 168 s.
12. Perepelica V.A., Popova E.V. Matematicheskie modeli i metody ocenki riskov jekonomicheskih, social'nyh i agrarnyh processov. - Rostov n/D: Izd-vo Rost. Un-ta, 2002. -202s.
13. Peters Je. Haos i porjadok na rynkah kapitala. Novyj analiticheskij vzgljad na cikly, ceny i izmenchivost' rynka / Je. Peters. - M.: Mir, 2000. - 333 s.
14. Peters Je. Fraktal'nyj analiz finansovyh rynkov: Primenenie teorii Haosa v investicijah i jekonomike. - M.: Internet-trejding, 2004. - 304s.
15. Popova E. V. Metody modelirovanija povedenija jekonomicheskih sistem na osnove analiza vremennyh rjadov / E. V. Popova, A. M. Kumratova, M. I. Popova // Jekonomicheskoe prognozirovanie: modeli i metody. Materialy X mezhdunar. nauchno-prakt. konf. - Voronezh. - 2014. - S. 200-206.
16. Popova E. V. Upravlenie riskami v voprosah bezopasnosti investicij v APK / E. V. Popova, A. M. Kumratova // Jekonomicheskoe prognozirovanie: modeli i metody. Materialy X mezhdunar. nauchno-prakt. konf. - Voronezh, - 2014. - S. 194-200.
17. Sigjel Je. Prakticheskaja biznes-statistika. - M.: Izdatel'skij dom «Vil'jams», 2004. -1056s.
18. Feder E. Fraktaly. - M.: Mir, 1991. - 260 s.
19. Jangishieva A. M. Modelirovanie jekonomicheskih riskov metodami nelinejnoj dinamiki. Avtoreferat dissertacii na soiskanie uchenoj stepeni kandidata jekonomicheskih nauk / Stavropol'skij gosudarstvennyj universitet. Stavropol', 2005. - 24 s.
