Фазовый анализ для оценки цикличности временного ряда урожайности зерновых культур
Исследование алгоритма фазового анализа как одного из методов нелинейной динамики для исследования временных рядов на цикличность. Разработка алгоритма формирования нечеткого множества длин квазициклов. Обоснование адекватности предпрогнозного анализа.
Рубрика | Экономико-математическое моделирование |
Вид | статья |
Язык | русский |
Дата добавления | 22.05.2017 |
Размер файла | 111,8 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru
Финансовый университет при Правительстве Российской Федерации
УДК 519.21
08.00.00 Экономические науки
ФАЗОВЫЙ АНАЛИЗ ДЛЯ ОЦЕНКИ ЦИКЛИЧНОСТИ ВРЕМЕННОГО РЯДА УРОЖАЙНОСТИ ЗЕРНОВЫХ КУЛЬТУР
Темиров Астемир Алиевич
аспирант
E-mail: t-ast@mail.ru
Москва, Россия
В статье рассматривается алгоритм фазового анализа как один из методов нелинейной динамики для исследования временных рядов на цикличность. В настоящее время большой интерес для науки представляют методы нелинейной динамики. Особенно стало актуальным применение их в анализе временных рядов. Существующие в классической эконометрике методы оценки цикличности разработаны для случайных систем, динамика которых подчиняется нормальному закону распределения. Известно, что в природе также существуют и неслучайные системы, характеризующиеся трендами, периодическими и непериодическими циклами, называемыми квазициклами. В статье представляется вычислительный процесс выявления квазициклов, который иллюстрируется на примере временного ряда урожайности «зерновых всего» в России за последние 119 лет. Фазовый портрет рассматриваемого временного ряда представляется в пространстве размерности два и состоит из 22-х неустойчивых, в плане периодичности, квазициклов совокупность которых образует странный аттрактор. Статистика квазициклов представлена табличной и геометрической моделями, а так же в виде нечеткого множества. Квазициклы имеют количественную (длина), и качественную (конфигурация), характеристики. Совокупность этих двух характеристик определяют собой такое важное свойство временного ряда как «трендоустойчивость». В статье задействован математический аппарат теории нечеткого множества, на базе которой разработан и представлен алгоритм формирования нечеткого множества длин квазициклов. Фазовый анализ является сильной формой анализа временных рядов на цикличность и служит достаточно адекватным инструментом для предпрогнозного анализа
Ключевые слова: НЕЛИНЕЙНАЯ ДИНАМИКА, ВРЕМЕННОЙ РЯД, ФАЗОВЫЙ АНАЛИЗ, ФАЗОВЫЙ ПОРТРЕТ, АТТРАКТОР, НЕЧЕТКОЕ МНОЖЕСТВО, КВАЗИЦИКЛИЧНОСТЬ
В экономической науке анализ временных рядов (ВР) занимает особое место. Исследование временной структуры данных в реальных экономических процессах позволяет адекватно отразить их в математических и экономических моделях. Исследуя временные ряды, по сути, происходит оценка экономических и социальных связей между переменными, что немаловажно для процесса принятия управленческих решений. Понимание и осознание этого факта всегда приводит к развитию различных методов анализа временных рядов.
В настоящее время большой интерес для науки представляют методы нелинейной динамики. Особенно стало актуальным применение их в анализе временных рядов. Это такие методы как: фрактальный анализ, фазовый анализ с применением теории нечеткого множества и теории детерминированного хаоса. Также стали популярными прогнозные модели на базе интеллектуальных систем: нейронные сети, клеточные автоматы, генетические алгоритмы, дающие хорошие прогнозы с достаточно приемлемой погрешностью [4,5,6,7,8]. Развитию методов нелинейной динамики предшествовало также наметившийся прогресс в компьютерных технологиях, когда стало возможным визуализировать на экране дисплея результаты моделирования достаточно сложных процессов и систем в виде графиков и диаграмм.
Анализ многочисленных публикаций [4,5,6,8] показывает, что свойство нелинейности характерно для многих социальных, экономических, природных и других процессов. Особенно большое количество научной литературы и статей посвящено исследованиям финансовых временных рядов методами нелинейной динамики [6]. Одним из преимуществ этих методов является то, что применение их в анализе временных рядов не требует подчинения исследуемых систем нормальному закону распределения, как это принято в классических моделях анализа и прогнозирования. цикличность временной квазицикл предпрогнозный
При построении прогнозной модели очень важно знать предпрогнозные характеристики исследуемой системы. Опираясь на них можно осуществить выбор адекватной прогнозной модели и вместе с тем и принятие правильного управленческого решения. Одной из таких задач является задача исследования экономического временного ряда на цикличность.
Проблема цикличности в экономических системах всегда волновала умы великих ученых из многих областей науки. Однако среди экономистов, признающих цикличность, нет единого мнения, относительно природы данного явления. Традиционно выделяют экзогенный и эндогенный подходы к объяснению этого явления. Основателем теории внешних факторов является английский экономист У.С. Джевонс, связавший экономический цикл с 11-летним циклом солнечной активности. Он связывал цикличность солнечной активности в основном с сельским хозяйством и торговлей, а его последователи распространили влияние солнечного цикла на всю экономику.
Существует классические методы выявления циклической компоненты в экономических временных рядах. Это известные в эконометрике автокорреляции различных порядков, также некоторые модели из серии ARIMA. Основная особенность этих моделей заключается в том, что поведение исследуемых временных рядов должно подчиняться нормальному закону распределения. Последнее не всегда верно для многих временных рядов [4,6,8].
Фазовый анализ (ФА) как один из инструментальных методов нелинейной динамики представляет интерес в исследовании цикличности в экономических временных рядах [5,6]. Анализ публикаций по теме исследования показывает, что существуют так называемые графические тесты хаоса Гилмора [1,7], которые выявляют неустойчивые квазипериодические конфигурации орбит, заключенные в странном аттракторе. Для определения таких орбит наиболее адекватным представляется метод «построения фазового портрета и его разложение на квазициклы» [5,6,7].
В настоящей статье предлагается следующий алгоритм фазового анализа, состоящий из 4-х этапов.
На первом этапе выбирается размерность фазового пространства . Второй этап заключается в построении фазового портрета изучаемой системы методом соединения соседних точек , либо отрезками, либо кривой. На третьем этапе фазовый портрет раскладывается на «квазициклы» и представляется в виде нечеткого множества.
Выше перечисленные этапы алгоритма фазового анализа проиллюстрируем на примере временного ряда урожайности «зерновых всего» в России за период с 1896 г. по 2014 г., графическое представление которого визуализируется на рис. 1.
Рисунок 1. Временной ряд урожайности зерновых культур в России за период 1896-2007гг
Введем обозначение этого ВР
, . (1)
Число наблюдений для ВР (1) составляет при этом каждое наблюдение приравнивается к одному году. Численные значения наблюдений определяют средний объем выхода урожайности зерновых на текущий год и измеряется в ц/га. Для различных экономических временных рядов достаточным является построение фазового портрета в фазовом пространстве размерности два, т.е. . В литературе [5] отмечено, что для конструирования истинного фазового пространства, необходимо знать все переменные, релевантные изучаемой системе. В случае временных рядов, мы имеем дело с одной лишь известной динамической переменной. Для восстановления фазового пространства по одной динамической переменной Паккард и Рюэль [4] предложили простой метод, который наполняет другие размерности посредством запаздывающих значений одной наблюдаемой переменной. Есть исходный временной ряд, уровни которого откладываются по оси абсцисс, и есть второй временной ряд, значения уровней которого откладываются по оси ординат, и являющийся той же реализацией, но с отставанием на один период. Математическое объяснение Паккарда такому механизму построения фазового портрета состоит в том, что нелинейные динамические системы являются детерминированными, и являются внутренне зависимыми системами. Это означает, что текущие величины каждой переменной есть результат воздействия на них прошлых величин [4,5].
Рассмотрим фазовый портрет исследуемого ВР (1), построенный по механизму Паккарда [4,5] и представленный на рис. 2 в виде траектории точек, в которой каждая соседняя пара соединена кривой. Визуализация этого фазового портрета свидетельствует о циклической природе рассматриваемого ВР . Для получения числовых и качественных характеристик этой цикличности разложим его на «квазициклы». Термин «квазицикл» в переводе с греческого языка означает «как бы цикл» и он близок к определению «цикла». Отличие квазицикла от цикла состоит в следующем:
- во-первых, квазицикл допускает точного несовпадения начала цикла с его концом;
- во-вторых, для квазицикла достаточно вхождение конечной его точки в область начальной точки;
- в-третьих, допускает самопересечение начального и конечного звеньев квазицикла, если это приводит к наименьшему расстоянию между точками начала и конца его.
Рисунок 2. Фазовый портрет временного ряда урожайности зерновых культур (1)
В целом фазовый портрет исследуемого ВР разложился на 22 последовательных неустойчивых в плане периодичности квазицикла , , совокупность которых демонстрирует странный аттрактор [7,8]. Квазициклы имеют количественную (длина), и качественную (конфигурация) характеристики. Число точек в квазицикле называется его длиной и обозначается через . Типичные квазициклы различных длин и конфигураций для временного ряда представлены на рис. 3.
а) б) в)
Рисунок 3. Типичные квазициклы в фазовом пространстве для ВР :
а, б - квазициклы ,, не имеющие точку самопересечения;
в - квазицикл , имеющий точку самопересечения
Типичный квазицикл , длины , представленного на рис. 3а, например, означает, что цикл определяется 5-тью уровнями исходного ВР и при этом ВР обладает свойством 5-летней цикличности. Статистика длин квазициклов исследуемого ВР (1) представлена в табл. 1.
Таблица 1. Статистика длин квазициклов ВР (1)
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
21 |
22 |
||
5 |
7 |
4 |
8 |
6 |
5 |
4 |
4 |
4 |
5 |
6 |
3 |
4 |
3 |
4 |
5 |
3 |
5 |
5 |
4 |
3 |
4 |
Таким образом, представленный на рис. 1. временной ряд урожайности «зерновых всего» в России фактически состоит из 22 завершенных квазициклов, которые в совокупности включают в себя 103 годовых уровня исходного ВР. Отсюда получаем среднее значение длины квазициклов и составляет лет. Характерной особенностью квазициклов рассматриваемого временного ряда урожайности зерновых является то, что чаще всего появляются квазициклы длины 4 и 5 в области значений длин , откуда вытекает 5-летнее значение средней длины годичных квазициклов.
При реальном моделировании сложных систем значительная часть информации, необходимой для их математического описания, существует в виде нечетких представлений или пожеланий экспертов. Для преодоления трудностей представления неточных понятий, анализа и моделирования систем, в которых участвует человек, американским математиком Лотфи Заде была предложена так называемая теория нечетких множеств [3,4]. Представим в виде нечеткого множества (НМ) статистику длин квазициклов. Для этих целей данные частотного распределения квазициклов представим в виде табличной модели (см.табл. 2).
Таблица 2. Табличная модель представления нечеткого множества длин квазициклов
Длина квазициклов, |
3 |
4 |
5 |
6 |
7 |
8 |
|
Частота появления, соответствующих длин квазициклов, |
4 |
6 |
7 |
3 |
1 |
1 |
|
Доля квазицикла в общем числе, |
0,18 |
0,27 |
0,32 |
0,14 |
0,05 |
0,05 |
|
Функция принадлежности, |
0,56 |
0,84 |
0,98 |
0,42 |
0,14 |
0,14 |
Предлагается следующий алгоритм формирования нечеткого множества длин квазициклов.
На первом шаге экспертным путем [3,4], с учетом того, что функция принадлежности может принимать значения в интервале от 0 до 1, присваиваем наиболее часто встречающемуся квазициклу, а это , заведомо высокое значение функции принадлежности, и пусть это будет , т.е. . Значение функции принадлежности равное единице соответствует истинности высказывания, а - соответствует ее отсутствию [3,4].
На втором шаге, для остальных квазициклов вычисляем значения функции принадлежности по формуле:
,
где и заполняем результатами вычисления последнюю строку табл.2.
На третьем шаге попарно объединяем элементы первой и последней строки табл.2 и представляем запись нечеткого множества так, как это принято в литературах по теории нечеткого множества [3]. Полученная оценка длин квазициклов для ВР представляется в виде нечеткого множества вида (2):
. (2)
Для наглядности на рис.4 представлено графическое изображение НМ длин квазициклов для ВР .
Рисунок 4. Графическое представление нечеткого множества длин квазициклов для ВР
Из рис.4 видно, что наиболее значимыми для ВР оказались квазициклы длины 3, 4 и 5. Об этом свидетельствуют не только численные значения их соответствующих функций принадлежности , и , но и попадание их в область выше [3]. Следует отметить при этом, что менее выраженным из них оказался трехлетний цикл.
Центр тяжести нечеткого множества , вычисленный по формуле
,
выдал результат, вполне согласующийся с суждением о приблизительной 5-годичной цикличности урожайности зерновых.
По результатам фазового анализа можно сделать заключение о том что, рассматриваемый временной ряд урожайности «зерновых всего» в России содержит компоненту цикл, которая имеет непериодический характер и тем самым порождает «квазициклы» различных конфигураций.
Свойство цикличности, присущее динамике урожайности зерновых, обязательно должно учитываться при планировании, прогнозировании и управлении сельскохозяйственным производством.
Результаты фазового анализа позволяют сделать обоснованные выводы о том, что в динамике ВР урожайности зерновых присутствует свойство «квазицикличности», которое обычно присуще нелинейным процессам и системам [5, 6, 7]. Последнее дает основание утверждать о неадекватности классических прогнозных моделей для их прогнозирования. Желательно для качественного их прогнозирования применить методы, которые базируются на таких интеллектуальных системах, как: нейронные сети, клеточные автоматы, генетические алгоритмы, зарекомендовавшие себя сильными прогнозными моделями.
Литература
1. Gilmore C.G. A new test for chaos / C.G. Gilmore //Journal of economic behavior and organization, №22, 1993. - P. 209-237.
2. Бессонов В.А. Введение в анализ российской макроэкономической динамики переходного периода/ В.А. Бессонов. - М.: ЦЭМИ РАН, 2003. - 151 с.
3. Кофман А., Хил Алуха Х. Введение теории нечетких множеств в управлении предприятиями / А. Кофман, Х. Хил Алуха.- Минск: Высшая школа, 1992.
4. Перепелица В.А. Структурирование данных методами нелинейной динамики для двухуровневого моделирования /В.А. Перепелица, Ф.Б. Тебуева, Л.Г. Темирова - Ставрополь: Ставропольское книжное издательство, 2006. - 284 с.
5. Петерс Э. Фрактальный анализ финансовых рынков: Применение теории Хаоса в инвестициях и экономике / Э. Петерс - М.: Интернет-Трейдинг, 2004. - 304 с.
6. Петерс Э. Хаос и порядок на рынках капитала. Новый аналитический взгляд на циклы, цены и изменчивость рынка / Э. Петерс - М.: Мир, 2000. - 333 с.
7. Сергеева Л.Н. Моделирование поведения экономических систем методами нелинейной динамики (теории хаоса) / Л.Н. Сергеева - Запорожье: ЗГУ, 2002. - 227 с.
8. Шустер Г. Детерминированный хаос: Введение / Г. Шустер - М.: Мир, 1988. - 240 с.
Размещено на Allbest.ru
...Подобные документы
Предпрогнозное исследование рядов урожайности с применением фрактального и R/S-анализа, бинарной кодировки. Расчет коэффициента Херста природных и экономических рядов. Оценка соотношения "детерминированность-стохастичность" для разных областей Украины.
курсовая работа [2,2 M], добавлен 18.09.2010Основные элементы эконометрического анализа временных рядов. Задачи анализа и их первоначальная обработка. Решение задач кратко- и среднесрочного прогноза значений временного ряда. Методы нахождения параметров уравнения тренда. Метод наименьших квадратов.
контрольная работа [37,6 K], добавлен 03.06.2009Анализ временных рядов с помощью статистического пакета "Minitab". Механизм изменения уровней ряда. Trend Analysis – анализ линии тренда с аппроксимирующими кривыми (линейная, квадратическая, экспоненциальная, логистическая). Декомпозиция временного ряда.
методичка [1,2 M], добавлен 21.01.2011Теория и анализ временных рядов. Построение линии тренда и прогнозирование развития случайного процесса на основе временного ряда. Сглаживание временного ряда, задача выделения тренда, определение вида тенденции. Выделение тригонометрической составляющей.
курсовая работа [722,6 K], добавлен 09.07.2019Классические подходы к анализу финансовых рынков, алгоритмы машинного обучения. Модель ансамблей классификационных деревьев для прогнозирования динамики финансовых временных рядов. Выбор алгоритма для анализа данных. Практическая реализация модели.
дипломная работа [1,5 M], добавлен 21.09.2016Построение эконометрических моделей и адекватная оценка их параметров для принятия обоснованных экономических решений. Проведение анализа и краткосрочного прогнозирования урожайности зерновых культур в Нижнем Поволжье методом многократного выравнивания.
реферат [51,4 K], добавлен 25.02.2011Статистические методы анализа одномерных временных рядов, решение задач по анализу и прогнозированию, построение графика исследуемого показателя. Критерии выявления компонент рядов, проверка гипотезы о случайности ряда и значения стандартных ошибок.
контрольная работа [325,2 K], добавлен 13.08.2010Динамика средней урожайности озимой пшеницы для областей Украины. Циклические изменения объемов урожая. Составление прогнозной модели урожайности зерновых. Методика оценки рисков зернопроизводства на основе связи между урожайностью и рентабельностью.
контрольная работа [1,7 M], добавлен 18.07.2010Изучение понятия имитационного моделирования. Имитационная модель временного ряда. Анализ показателей динамики развития экономических процессов. Аномальные уровни ряда. Автокорреляция и временной лаг. Оценка адекватности и точности трендовых моделей.
курсовая работа [148,3 K], добавлен 26.12.2014Теоретические выкладки в области теории хаоса. Методы, которые используются в математике, для прогнозирования стохастических рядов. Анализ финансовых рядов и рядов Twitter, связь между сентиметными графиками и поведением временного финансового ряда.
курсовая работа [388,9 K], добавлен 01.07.2017Временные ряды и их характеристики. Факторы, влияющие на значения временного ряда. Тренд и сезонные составляющие. Декомпозиция временных рядов. Метод экспоненциального сглаживания. Построение регрессионной модели. Числовые характеристики переменных.
контрольная работа [1,6 M], добавлен 18.06.2012Описание алгоритма культурного обмена и проведение экспериментального исследования средней трудоемкости алгоритма случайного поиска. Основные идеи алгоритма и эффективность итерационных методов решения. Зависимость функции качества от длины генотипа.
курсовая работа [373,3 K], добавлен 24.06.2012Эконометрическое моделирование стоимости квартир в Московской области. Исследование динамики экономического показателя на основе анализа одномерного временного ряда. Параметры линейной парной регрессии. Оценка адекватности модели, осуществление прогноза.
контрольная работа [925,5 K], добавлен 07.09.2011Решение задачи изучения изменения анализируемых показателей во времени при помощи построения и анализа рядов динамики. Элементы ряда динамики: уровни динамического ряда и период времени, за который они представлены. Понятие переменной и постоянной базы.
методичка [43,0 K], добавлен 15.11.2010Анализ упорядоченных данных, полученных последовательно (во времени). Модели компонентов детерминированной составляющей временного ряда. Свободные от закона распределения критерии проверки ряда на случайность. Теоретический анализ системы линейного вида.
учебное пособие [459,3 K], добавлен 19.03.2011Методика отбора факторов, влияющих на выходной показатель в статистике. Выравнивание динамических рядов. Показатели анализа ряда динамики. Множественное уравнение регрессии. Проверка адекватности регрессионной модели. Осуществление прогнозных расчетов.
курсовая работа [2,0 M], добавлен 23.01.2012Эконометрическое моделирование стоимости квартир в московской области. Матрица парных коэффициентов корреляции. Расчет параметров линейной парной регрессии. Исследование динамики экономического показателя на основе анализа одномерного временного ряда.
контрольная работа [298,2 K], добавлен 19.01.2011Влияние девальвации национальной валюты на цены активов и процентных ставок на фондовый рынок. Анализ отраслевых взаимосвязей и закономерностей в динамике биржевых индикаторов и множества других временных рядов. Оценка моделей методом "rolling window".
дипломная работа [1,7 M], добавлен 06.11.2015Структурные компоненты детерминированной составляющей. Основная цель статистического анализа временных рядов. Экстраполяционное прогнозирование экономических процессов. Выявление аномальных наблюдений, а также построение моделей временных рядов.
курсовая работа [126,0 K], добавлен 11.03.2014Понятие о рядах динамики, их роль. Показатели анализа ряда динамики. Средние показатели по рядам динамики. Статистическое изучение сезонных колебаний. Методы анализа основной тенденции в рядах динамики. Экстраполяция тенденции как метод прогнозирования.
курсовая работа [106,6 K], добавлен 14.10.2008