Студопедия

КАТЕГОРИИ:


Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748)

Выравнивание при помощи ряда Фурье.




Решение.

1. Анализ данных.

Строим график на основе данных таблицы.

Период времени
Объем закупок, млн. руб.

Визуальный анализ данных позволяет сделать вывод о возможности использования мультипликативной модели, т.к. размах вариаций фактических значений относительно линии тренда постоянно возрастает:

Y=T´ S´ E

 

2. Расчет сезонной компоненты.

Процедура этого этапа включает расчет скользящей средней с шагом, равным 4, центрированной скользящей средней, оценку сезонной компоненты на основе деления уровня ряда на значение центрированной скользящей средней за соответствующий момент времени. Результаты расчетов приведены в таблице.

 

Период времени Номер квартала Объем закупок, млн. руб. Скользящая средняя за 4 квартала Центрированная скользящая средняя Коэффициент сезонности
Год Квартал
  I     I                
           
II                
           
III       138,25 0,940
      140,5    
IV       140,5 1,011
      140,5    
Период времени Номер квартала Объем закупок, млн. руб. Скользящая средняя за 4 квартала Центрированная скользящая средняя Коэффициент сезонности
  Год   Квартал
    II I       140,5 1,21
      141,5    
II       144,25 0,915
           
III       148,25 0,904
      149,5    
IV       150,25 1,092
      151,0    
  III I       152,25 1,103
      153,25    
II       154,63 0,892
      156,00    
III       178,50 0,909
      161,00    
IV          
           
IV I          

Далее определяются средние значения сезонной компоненты по кварталам, и проводится их корректировка.

Среднее значение для каждого квартала года вычисляется на основе фактически имеющихся данных по одноименным кварталам по рассматриваемым годам. Сумма оценок сезонной компоненты должна равняться 4. Если сумма не равна 4, то производится корректировка сезонной компоненты.

Результаты расчетов приведены в таблице.

  Период времени
Год Квартал
I II III IV  
I 0,940 1,011
II 1,121 0,915 0,904 1,092
III 1,103 0,892 0,909
Итого     2,224 1,807 2,753 2,103
Оценка сезонной компоненты     1,112 0,903 0,918 1,051 Сумма = 3,984
Скорректированная сезонная компонента     1,116 0,907 0,922 1,052 Сумма = 4

 

 

3. Определение тренда.

На этом этапе проводится десезонализация данных путем деления фактических значений уровня ряда на скорректированные коэффициенты сезонности по кварталам. Таким образом, получают значения, содержащие тренд и случайную компоненту:

Номер квартала Объем закупок, млн.руб. Коэффициент сезонности Десезонализированный объем, млн.руб. Y /S=ТЕ
    1,116 125,4
    0,907 145,6
    0,922 141,2
    1,055 134,3
    1,116 141,6
    0,907 145,6
    0,922 145,4
    1,055 155,4
    1,116 150,4
    0,907 152,2
    0,922 156,4
    1,055 164,8
    1,116 168,4

 

Далее осуществляется построение модели тренда методом наименьших квадратов на основе десезонализированных данных объема экспорта. Визуальное сравнение графиков фактических данных и десезонализированных данных даёт возможность принять гипотезу о линейности модели тренда:

Tt=a0+a1t

Номер квартала
Десезонализированный объем, млн.руб.

Модель тренда с численными параметрами будет иметь вид:

Tt= 125,8 + 3,2 t

4. Определение качества модели и расчет ошибок.

О качестве модели можно судить по ошибкам. Рассчитываются они по формуле:

Е =

 

Номер квартала Объем закупок, млн.руб. Сезонная компонента Трендовое значение Т, млн.руб. Ошибка
T S
    1,116 129,0 143,9 0,97
    0,907 132,2 119,9 1,10
    0,922 135,4 124,8 1,04
    1,055 138,6 146,2 0,97
    1,116 141,8 158,2 1,01
    0,907 145,0 131,5 1,00
    0,922 148,2 136,6 0,98
    1,055 151,4 159,7 1,03
    1,116 154,6 172,5 0,97
    0,907 157,8 143,1 0,96
    0,922 161,0 148,4 1,03
    1,055 164,2 173,2 1,00
    1,116 167,4 186,8 1,01

 

Значение ошибок приведены в последней графе таблицы. Значения ошибок невелики и составляют порядка 1%. Это говорит о таком качестве модели, которое делает её вполне пригодной для прогнозирования (1 % - 3 %).

 

5. Построение прогноза с учетом сезонных колебаний.

Здесь осуществляется расчет прогноза на II квартал IV года по модели тренда с последующей корректировкой прогноза на сезонную компоненту.

Рассчитаем прогноз на II квартал по модели тренда. Порядковый номер II квартала IV года равен 14, тогда прогноз значения объема экспорта по уравнению тренда будет следующим:

TII кв.IV год = 125,8 + 3,2х14

Скорректируем это значение с учетом фактора сезонности:

Ŷ скорII кв.IV год = (125,8+3,2х15) 0,907 = 154,7

Однако не следует забывать о том, что чем больше период упреждения, тем меньше степень обоснованности прогноза.

 

 

Одним из приемов аналитического выравнивания является выравнивание при помощи ряда Фурье. По тригонометрической системе ряд Фурье определяется для каждой функции f, интегрируемой на отрезке [0, 2 ].

Это ряд:

(1.1)

с коэффициентами:

(1.2)

Аналогично строятся ряды Фурье для функций от многих переменных.

С помощью выражения (1.2) строится ряд Фурье (1.1) для каждой функции интегрируемой на отрезке [0, 2 ].

Простейшую периодическую функцию вида, приведенную ниже, принято называть гармоникой.

Это функция хорошо описывает многие колебательные процессы.

A -называют амплитудой, ω - частотой, φ - начальной фазой, - периодом колебания. Функции вида , называют соответственно второй, третьей гармониками относительно основной гармоники.

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

, (1.3)

Описанные выше функции используются в гармоническом анализе в разделе математики, который позволяет сводить решение некоторых задач к вопросам гармонического анализа.

Довольно часто при исследовании социально-экономических процессов используют следующее уравнение:

Это уравнение позволяет выравнивать процессы с характерной периодичной составляющей, где k определяет в данном случае гармонику ряда.

При изучении процессов с использованием ряда Фурье обычно рассчитывают несколько гармоник (чаще всего четыре), а затем устанавливают, какая из них лучше отражает периодичность изменения.

Оценки параметров уравнения вычисляют по методу наименьших квадратов:

Фактор времени t может быть выражен в градусах или в радианной мере. Исчисляется t от 0 с приростом, равным , где n - число уровней временного ряда.

Например, для n = 7 t будет записано так:

 

Значения синусов и косинусов для различных гармоник и необходимые расчеты для удобства проведения вычислений располагают в таблице по образцу (таблица представлена для двух гармоник).

Период времени t y y cos t y sin t y cos 2 t y sin2 t
                 

и -это теоретические значения показателя для k =1 и k =2, соответственно.

За прогнозное уравнение принимают то, для которого теоретические значения по временным периодам ближе к эмпирическим.

Следует заметить, что результаты по ряду Фурье дает наибольший эффект в тех случаях, когда временной ряд насчитывает большое число уровней.

Пример:

Выравнить по ряду Фурье временной ряд, описывающий представленные ниже сезонные изменения показателей цен на товар:

Период времени (месяцы)                        
Цена (тыс.руб) 3,7 4,0 4,4 5,2 4,6 7,0 6,0 4,8 4,6 3,8 3,6 3,5

Решение:

Определим для n = 12 моменты времени:

Построим по вышеприведенному образцу вспомогательные таблицы, в которые внесем результаты всех проведенных необходимых расчетов:

Период времени t cos t sin t cos 2t sin 2t y
    1,000   1,0   3,70
  π/6 0,866 0,500 0,5 0,866 4,00
  π/3 0,500 0,866 - 0,5 0,866 4,40
  π/2   1,000 - 1,0   5,20
  2π/3 - 0,500 0,866 - 0,5 - 0,866 4,60
  5π/6 - 0,866 0,500 0,5 - 0,866 7,00
  π - 1,000   1,0   6,00
  7π/6 - 0,866 - 0,500 0,5 0,866 4,80
  4π/3 - 0,500 - 0,866 - 0,5 0,866 4,60
  3π/2   - 1,000 - 1,0   3,80
  5π/3 0,500 - 0,866 - 0,5 - 0,866 3,60
  11π/6 0,866 - 0,500 0,5 - 0,866 3,50
            55,20

 

Период времени y cos t y sin t y cos 2t y sin 2t
  3,7 0,000 3,7   3,50 3,79
  3,464 2,000 2,0 3,464 3,93 3,96
  2,2 3,810 - 2,2 3,81 4,55 4,30
    5,200 - 5,2   5,17 4,38
  - 2,3 3,984 - 2,3 - 3,984 5,65 5,62
  - 6,062 3,5 3,5 - 6,06 5,84 6,10
  - 6,0   6,0   5,70 5,99
  - 4,157 - 2,4 2,4 4,157 5,27 5,30
  - 2,3 - 3,984 - 2,3 3,984 4,65 4,40
    - 3,8 - 3,8   3,93 3,64
  1,8 - 3,117 - 1,8 - 3,117 3,55 3,52
  3,031 - 1,75 1,75 - 3,031 3,36 3,62
  - 6,624 3,443 1,75 - 0,778 55,10 55,12

 

 

Итак, уравнение для k = 1 имеет вид:

Для к = 2 имеем:

Значения и представлены в двух последних графах одной из вспомогательных таблиц.

Из таблицы видно, что значения цены, рассчитанные для двух гармоник

(k = 2) ближе к эмпирическим уровням ряда, чем для k = 1. Это означает, что выбор падает на уравнение, полученное для двух гармоники. На основе него следует далее строить прогноз.

 




Поделиться с друзьями:


Дата добавления: 2017-02-01; Просмотров: 126; Нарушение авторских прав?; Мы поможем в написании вашей работы!


Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет



studopediasu.com - Студопедия (2013 - 2026) год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! Последнее добавление




Генерация страницы за: 0.013 сек.