Реферат: Прогнозирование емкости и конъюнктуры рынка
Гипероглавление:Раздел 1 «Прогнозирование на основе стационарного временного ряда»
2. Отобразим исходные данные на рис.1:
3. Анализ методом коэффициента Кендэла:
4. Прогнозирование стационарного процесса:
2. Отобразим исходные данные на рис.2:
3. Сглаживание исходных данных методом скользящей средней:
4. Анализ методом Фостера-Стюарта:
5. Анализ методом коэффициента Кендэла:
6. Определение параметров линейного тренда методом усреднения по левой и правой половине ряда:
7. Определение параметров линейного тренда методом наименьших квадратов:
8. Определение параметров нелинейной трендовой модели:
9. Рассчитаем параметры для показательного тренда:
10.Сравним линейный и выбранный нелинейный – показательный тренды:
11. Сравнение трендовых моделей с помощью критерия наименьшей суммы квадратов отклонений:
12. Оценка адекватности выбранной трендовой модели:
13. Прогнозирование на основе трендовой модели:
13.1 Точечный прогноз:
Раздел 3 «Прогнозирование на основе сезонного цикла временного ряда»
2. Отобразим исходные данные на рис.8:
3. Анализ методом коэффициента Кендэла:
4.Сезонное прогнозирование:
Раздел 4 «Прогнозирование с помощью метода экспоненциального сглаживания»
2. Отобразим исходные данные на рис.9:
3. Метод экспоненциального сглаживания:
4.Критерий наименьшей суммы квадратов отклонений.
--PAGE_BREAK-- продолжение
--PAGE_BREAK--3. Анализ методом коэффициента Кендэла:
Оценим наличие тенденции среднего уровня ряда в исходных данных, с помощью коэффициента Кендэла. Отобразим расчет в таблице 2.
Таблица 2
t
Yt
P
t
1
8,7
-
2
10,0
1
3
10,3
2
4
10,1
2
5
10,4
4
6
9,1
1
7
9,4
2
8
9,2
2
9
9,9
4
10
9,7
4
Итого
-
22
Рассчитаем число случаев превышения текущим уровнем ряда предыдущих ему уровней ряда. Первый уровень ряда у1=8,7 не с чем сравнить (нет предыдущих уровней ряда), поэтому в графе 3 поставим прочерк. Второй уровень ряда у2=10,0 больше предыдущего у1=8,7 и к тому же он всего один, поэтому в графе 3 ставим 1. Третий уровень у3=10,3 больше у2=10,0 и у1=8,7, поэтому в графе 3 ставим 2. Четвертый уровень у4=10,1 меньше у3=10,3, но больше у2=10,0 и у1=8,7, поэтому в графе 3 ставим 2… Аналогичным образом определим число таких случаев и для остальных уровней ряда.
Подведя итог по графе Pt, найдем общее число случаев, когда текущий уровень ряда больше предыдущих по формуле :
Р=Σ Рt=1+2+2+4+1+2+2+4+4=22.<img width=«12» height=«23» src=«ref-2_325690324-73.coolpic» v:shapes="_x0000_i1244"> (1)
Определим расчетное значение коэффициента Кендэла по формуле:
<img width=«340» height=«47» src=«ref-2_325690397-685.coolpic» v:shapes="_x0000_i1245"> (2)
Рассчитаем теоретическую дисперсию по формуле:
<img width=«331» height=«47» src=«ref-2_325691082-1335.coolpic» v:shapes="_x0000_i1025"> (3)
Для оценки наличия в ряде тенденции среднего уровня ряда выберем вероятность, равную 0,95 (95%). С учетом выбранной вероятности коэффициент доверия t=1,96.
Сопоставим расчетное и теоретическое значения коэффициента Кендэла. При сопоставлении может возникнуть три варианта:
Первый вариант, когда с вероятностью t во временном ряде нет тренда;
<img width=«252» height=«38» src=«ref-2_325692417-921.coolpic» v:shapes="_x0000_i1026"> (4)
<img width=«353» height=«30» src=«ref-2_325693338-1307.coolpic» v:shapes="_x0000_i1246">
<img width=«250» height=«24» src=«ref-2_325694645-939.coolpic» v:shapes="_x0000_i1247">
— соотношение выполняется
Второй вариант, когда с вероятностью t во временном ряде есть убывающая тенденция среднего уровня ряда;
<img width=«147» height=«40» src=«ref-2_325695584-605.coolpic» v:shapes="_x0000_i1027"> (5)
<img width=«422» height=«25» src=«ref-2_325696189-1099.coolpic» v:shapes="_x0000_i1248">
Третий вариант, когда с вероятностью – t во временном ряде есть возрастающая тенденция среднего уровня ряда.
<img width=«150» height=«43» src=«ref-2_325697288-680.coolpic» v:shapes="_x0000_i1028"> (6)
<img width=«407» height=«25» src=«ref-2_325697968-1029.coolpic» v:shapes="_x0000_i1249">
Из трех вариантов мы выбираем первый, поскольку только в нем выполняется необходимое соотношение расчетного и теоретического значений коэффициента Кендэла
Из установленного соотношения следует, что с вероятностью 95% во временном ряде отсутствует тенденция среднего уровня ряда.
На основе ранее полученных частных выводов, можно сделать обобщенный вывод: с высокой вероятностью 95% во временном ряде отсутствует тенденция среднего уровня ряда. С учетом обобщенного вывода можно считать, что временный ряд является стационарным временным рядом или стационарным процессом.
продолжение
--PAGE_BREAK--4. Прогнозирование стационарного процесса:
Исходя из ранее произведенных анализов, выяснилось, что временной ряд является стационарным. При нахождении прогноза по такому временному ряду вначале определяется точечный прогноз, а затем уже – интервальный.
Определим точечный и интервальный прогнозы оборота овощной палатки на следующий (одиннадцатый) уровень ряда.
Для простоты расчетов воспользуемся таблицей 3.
Таблица 3
t
yt
<img width=«60» height=«25» src=«ref-2_325698997-150.coolpic» v:shapes="_x0000_i1250">
<img width=«68» height=«26» src=«ref-2_325699147-169.coolpic» v:shapes="_x0000_i1251">
1
2
3
4
1
8,7
-0,98
0,9604
2
10,0
0,32
0,1024
3
10,3
0,62
0,3844
4
10,1
0,42
0,1764
5
10,4
0,72
0,5184
6
9,1
-0,58
0,3364
7
9,4
-0,28
0,0784
8
9,2
-0,48
0,2304
9
9,9
,22
0,0484
10
9,7
0,02
0,0004
Итого:
96,8
2,836
Для расчета дисперсии ряда при определении интервального прогноза понадобиться разница между уt
и<img width=«19» height=«34» src=«ref-2_325699316-185.coolpic» v:shapes="_x0000_i1252">. Рассчитаем эту разницу для каждого уровня ряда и занесем в соответствующие строки графы 3. В графе 4 возведем выражение <img width=«60» height=«25» src=«ref-2_325698997-150.coolpic» v:shapes="_x0000_i1253"> каждого ряда в квадрат.
Определим точечный прогноз по формуле:
<img width=«173» height=«43» src=«ref-2_325699651-399.coolpic» v:shapes="_x0000_i1254"> (7)
Перед определением интервального прогноза рассчитаем дисперсию ряда по формуле:
<img width=«233» height=«47» src=«ref-2_325700050-509.coolpic» v:shapes="_x0000_i1255"> (8)
Для того чтобы найти интервальный прогноз нам необходимо узнать значение tγ. Для этого выбираем уровень значимости, равный 0,05, т.е. а=0,05. Отсюда доверительная вероятность γ=1−а γ=1–0,05=0,95. Определим число степеней свободы k=n–1 k=10–1=9. Зная доверительную вероятность и число степеней свободы по приложению 2 (2, стр. 70), найдем табличное значение t
γ. Оно будет равно 2,262.
Найдем интервальный прогноз по формуле:
<img width=«457» height=«55» src=«ref-2_325700559-1010.coolpic» v:shapes="_x0000_i1256"> (9)
Отсюда верхняя граница интервального прогноза 10,082 (9,68+0,402), а нижняя – 9,278 (9,68–0,402).
Таким образом, с вероятностью 95% прогнозный оборот овощной палатки на следующую декаду будет лежать между 10,082 и 9,278.
Раздел 2 «Прогнозирование на основе тренда временного ряда»
Исходные данные для Варианта №23:
1. Таблица с исходными данными:
Таблица 4. Ежедневный оборот магазина «Ткани для дома»
t
yt
1
10,2
2
10,8
3
10,4
4
11,9
5
12,2
6
12,5
7
13,1
8
12,4
9
13,6
10
14,3
11
14,9
12
13,8
продолжение
--PAGE_BREAK--2. Отобразим исходные данные на рис.2:
<img width=«660» height=«455» src=«ref-2_325701569-32737.coolpic» v:shapes="_x0000_i1257">
На основе визуального анализа со средней вероятностью можно сделать вывод: во временном ряду присутствует тенденция в виде тренда, и он имеет линейный характер.
3. Сглаживание исходных данных методом скользящей средней:
Чтобы полученная визуальная оценка была более убедительной и наглядной, осуществим сглаживание временного ряда с помощью метода скользящей средней с интервалом сглаживания, равным трем. Рассчитаем сглаженные уровни ряда по формуле:
<img width=«144» height=«43» src=«ref-2_325734306-485.coolpic» v:shapes="_x0000_i1029"> (10)
Так, первый сглаженный уровень ряда:
<img width=«304» height=«43» src=«ref-2_325734791-630.coolpic» v:shapes="_x0000_i1258">
Аналогично рассчитаем остальные сглаженные уровни ряда:
<img width=«306» height=«43» src=«ref-2_325735421-610.coolpic» v:shapes="_x0000_i1259">
<img width=«297» height=«43» src=«ref-2_325736031-587.coolpic» v:shapes="_x0000_i1260">
<img width=«299» height=«43» src=«ref-2_325736618-619.coolpic» v:shapes="_x0000_i1261">
<img width=«297» height=«43» src=«ref-2_325737237-608.coolpic» v:shapes="_x0000_i1262">
<img width=«308» height=«43» src=«ref-2_325737845-641.coolpic» v:shapes="_x0000_i1263">
<img width=«305» height=«43» src=«ref-2_325738486-596.coolpic» v:shapes="_x0000_i1264">
<img width=«309» height=«43» src=«ref-2_325739082-643.coolpic» v:shapes="_x0000_i1265">
<img width=«314» height=«43» src=«ref-2_325739725-656.coolpic» v:shapes="_x0000_i1266">
<img width=«320» height=«43» src=«ref-2_325740381-675.coolpic» v:shapes="_x0000_i1267">
Результаты расчетов внесем в таблицу 5:
Таблица 5
t
yt
y’t
1
10,2
-
2
10,8
10,47
3
10,4
11,03
4
11,9
11,50
5
12,2
12,20
6
12,5
12,60
7
13,1
12,67
8
12,4
13,03
9
13,6
13,43
10
14,3
14,27
11
14,9
14,33
12
13,8
-
И построим график сглаженного ряда на рис. 3:
<img width=«660» height=«410» src=«ref-2_325741056-27720.coolpic» v:shapes="_x0000_i1268">
На основе визуального анализа сглаженного временного ряда с высокой вероятностью можно сделать вывод: во временном ряду присутствует тенденция в виде тренда, и он имеет линейный характер.
4. Анализ методом Фостера-Стюарта:
Оценим данные, приведенные в таблице 5, с помощью метода Фостера–Стюарта с точки зрения наличия в них тенденции среднего уровня ряда и дисперсии. Расчет проведем с помощью таблицы 6.
Таблица 6
Для реализации этого метода вначале определим ut иlt
по следующим условиям :
<img width=«271» height=«60» src=«ref-2_325768776-1456.coolpic» v:shapes="_x0000_i1030"> в графе 3 определим ut.
Так как дляу1=10,2 нет предыдущего уровня у0, то в графе 3 поставим прочерк. Сравниваем у2=10,8 со всеми предыдущими уровнями ряда. Он всего один − у1=10,2. Поскольку у2>y1, постольку в графе 3 ставим 1. Сравниваем у3=10,4 со всеми предыдущими уровнями ряда (у2=10,8; у1=10,2). Так как у3 меньше хотя бы одного из предыдущих, а именно у2, то в графе 3 ставим 0. Сравниваем у4=11,9 со всеми предыдущими уровнями ряда (у3=10,4; у2=10,8; у1=10,2). Он больше у1, у2 и у3, поэтому в графе 3 ставим 1. Аналогично проводится сравнение и других уровней ряда.
<img width=«271» height=«60» src=«ref-2_325770232-1461.coolpic» v:shapes="_x0000_i1031"> для графы 4 определим lt
.
Расчет проводится так же как и для графы 3, но с обратным условием: текущий уровень ряда уtдолжен быть меньше всех предыдущих уровней:
Для у1 нет предыдущего уровня, значит в графе 3 ставим прочерк;
Для у2 – ставим 0 (у2>у1);
Для у3 – ставим 0 (у3< у2; у3>у1);
Для у4 – ставим 0 (у4>у1, у3, у2);
Для у5 – ставим 0 (у5>у4, у3,у1, у2,);
Для у6 – ставим 0 (у6>у4, у5, у3, у2,у1);
Для у7 – ставим 0 (у7>у1, у2,у3, у4, у5, у6);
Для у8 – ставим 0 (у8>у5, у4, у3,у2,у1; у8< у7, у6,);
Для у9 – ставим 0 (у9>у8, у7, у6, у4, у5, у3, у2,у1);
Для у10 – ставим 0 (у10>у9, у8, у6, у5, у3, у2, у1, у7,у4);
Для у11 – ставим 0 (у11>у10, у9, у8, у7, у6, у5, у4, у3; у2, у1);
Для у12– ставим 0 (у12<у11, у10; у12>у9, у8, у7, у6, у5,у4, у3, у2,у1);
Затем, на основе величин ut иlt
,для графы 5 определим величину Sпо формуле:
S=ΣSt, гдеSt= ut +lt; (11)
Для t=1 в графе 5 поставим прочерк. Рассчитаем величну Sдля t=2:
S
2
=
u
2 +l
2=1+0=1
Дляt=2 в графе 5 поставим 1. Для остальных уровней ряда проводится аналогичные расчеты. Результаты заносятся в графу 5 таблицы 6.
Затем найдем итоговую сумму по графе 5:
S
=
ΣSt
=1+0+1+1+1+1+0+1+1+1+0=8 (12)
Для графы 6 по формуле:
D=ΣDt, гдеDt= ut – lt (13)
Dtдля t=1 в графе 6 ставим прочерк. Найдем значения Dt
дляt
=2:
D
2
=
u
2 – l
2
=1-0=1
в графе 6 ставим 1. Для остальных уровней ряда проводится аналогичные расчеты. Результаты заносятся в графу 6 таблицы 6.
Затем найдем итоговую сумму по графе 6:
D
=
ΣDt
=1+0+1+1+1+1+0+1+1+1+0=8 (14)
Выдвинемнулевую гипотезу: во временном ряде нет тенденции среднего уровня и нет тенденции дисперсии.
Проверим выдвинутую нулевую гипотезу по формулам:
<img width=«76» height=«49» src=«ref-2_325771693-356.coolpic» v:shapes="_x0000_i1032">, <img width=«82» height=«50» src=«ref-2_325772049-375.coolpic» v:shapes="_x0000_i1033">, (15)
Найдем значения μ, σ1,σ2. В приложении 1 (2, стр. 70) приведены данные дляn=10 и для n=15, а нам надо найти данные для n=12.
Найдемμ для n=12 следующим образом. Значение μ для n=10 равно 3,858, для n=15 равно 4,363. Увеличение μ при изменении nна 2 шага найдем следующим образом:
<img width=«328» height=«44» src=«ref-2_325772424-1299.coolpic» v:shapes="_x0000_i1034"> (16)
Отсюда μ(12)=μ(10)+Δμ=3,858+0,311=4,169
Найдем σ1. Значение σ1 для n=10 равно 1,288, для n=15 равно 1,521. Увеличение σ1 при изменении nна 2 шага найдем следующим образом:
<img width=«353» height=«45» src=«ref-2_325773723-1268.coolpic» v:shapes="_x0000_i1035"> (17)
Отсюда σ1 (12)=σ1 (10)+Δσ1=1,288+0,093=1,381.
Найдем σ2. Значение σ2 для n=10 равно 1,964, для n=15 равно 2,153. Увеличение σ2 при изменении nна 2 шага найдем следующим образом:
<img width=«357» height=«45» src=«ref-2_325774991-1276.coolpic» v:shapes="_x0000_i1036"> (18)
Отсюда σ2 (12)=σ2 (10)+Δσ2=1,964+0,076=2,040
Найдем значения t1 и t2:
<img width=«191» height=«45» src=«ref-2_325776267-442.coolpic» v:shapes="_x0000_i1269"> <img width=«172» height=«45» src=«ref-2_325776709-408.coolpic» v:shapes="_x0000_i1270">
Найдем табличное значение tγ. Для этого зададимся уровнем значимости, например а=0,05 (это стандартная величина). Затем определим доверительную вероятность γ=1– а=1– 0,05=0,95 и число степеней свободы k
=
n
– 1=12 –1=11. Относительно найденных значений γ и kпо приложение 2 (2, стр. 70) найдем табличное значение tγ— оно равно 2,201.
Сопоставим значения t1 и t2с tγ:
Поскольку |t1=2,77|>|tγ=2,201|, гипотеза отвергается, следовательно, во временном ряде имеет место тенденция дисперсии.
Поскольку |t2=3,92|>|tγ=2,201|, гипотеза отвергается, следовательно, во временном ряде имеет место тенденция среднего уровня.
На основе сопоставлений с выбранной вероятностью 95% можно утверждать, что во временном ряде присутствуют тенденция дисперсии и тенденция среднего уровня.
продолжение
--PAGE_BREAK--5. Анализ методомкоэффициента Кендэла:
Оценим наличие тенденции среднего уровня ряда в исходных данных, с помощью коэффициента Кендэла. Расчет проведем с помощью данных таблицы 4, а результаты расчета приведем в таблице 7.
Таблица 7
t
yt
P
t
1
10,2
-
2
10,8
1
3
10,4
1
4
11,9
3
5
12,2
4
6
12,5
5
7
13,1
6
8
12,4
5
9
13,6
8
10
14,3
9
11
14,9
10
12
13,8
9
Итого
-
61
Рассчитаем число случаев превышения текущим уровнем ряда предыдущих ему уровней ряда.
Первый уровень ряда у1=10,2 не с чем сравнить (нет предыдущих уровней ряда), поэтому в графе 3 поставим прочерк. Второй уровень ряда у2=10,8 больше предыдущего у1=10,2, поэтому в графе 3 ставим 1. Третий уровень у3=10,4 больше у1=10,2 и меньше у2=10,8), поэтому в графе 3 ставим 1. Четвертый уровень у4=11,9 больше у3=10,4, у1=10,2, у2=10,8, поэтому в графе 3 ставим 3.
Аналогичным образом определим число таких случаев и для остальных уровней ряда.
Подведя итог по графе 3, найдем общее число случаев, когда текущий уровень ряда больше предыдущих по формуле (1):
Р
=ΣР
t
=1+1+3+4+5+6+5+8+9+10+9=61.
Определим расчетное значение коэффициента Кендэла по формуле (2):
<img width=«339» height=«47» src=«ref-2_325777117-699.coolpic» v:shapes="_x0000_i1271">
Рассчитаем теоретическую дисперсию по формуле (3):
<img width=«320» height=«47» src=«ref-2_325777816-814.coolpic» v:shapes="_x0000_i1272">
Для оценки наличия в ряде тенденции среднего уровня ряда выберем вероятность, равную 0,95 (95%). С учетом выбранной вероятности коэффициент доверия t=1,96,
Сопоставим расчетное и теоретическое значения коэффициента Кендэла. При сопоставлении может возникнуть три варианта.
Первый вариант, когда с вероятностью tво временном ряде нет тренда;
<img width=«252» height=«38» src=«ref-2_325692417-921.coolpic» v:shapes="_x0000_i1037">
<img width=«354» height=«30» src=«ref-2_325779551-1432.coolpic» v:shapes="_x0000_i1273">
<img width=«250» height=«24» src=«ref-2_325780983-986.coolpic» v:shapes="_x0000_i1274">соотношение не выполняется
Второй вариант, когда с вероятностью tво временном ряде есть убывающая тенденция среднего уровня ряда;
<img width=«147» height=«40» src=«ref-2_325695584-605.coolpic» v:shapes="_x0000_i1038">
<img width=«174» height=«25» src=«ref-2_325782574-738.coolpic» v:shapes="_x0000_i1275">соотношение не выполняется
Третий вариант, когда с вероятностью – tво временном ряде есть возрастающей тенденция среднего уровня ряда.
<img width=«150» height=«43» src=«ref-2_325697288-680.coolpic» v:shapes="_x0000_i1039">.
<img width=«159» height=«25» src=«ref-2_325783992-728.coolpic» v:shapes="_x0000_i1276"> соотношение выполняется
С выбранной вероятностью 95 % можно сделать вывод, что во временном ряде есть возрастающая тенденции среднего уровня.
Этот вывод согласуется с выводами, полученными нами ранее при визуальном анализе и применении метода Фостера-Стюарта.
продолжение
--PAGE_BREAK--6. Определение параметров линейного тренда методом усреднения по левой и правой половине ряда:
Рассчитаем параметры линейного тренда графическим методом – методом усреднения по левой и правой половине данных. Отобразим в таблице 8.
Таблица 8
t
yt
1
10,2
2
10,8
3
10,4
4
11,9
5
12,2
6
12,5
7
13,1
8
12,4
9
13,6
10
14,3
11
14,9
12
13,8
<img width=«185» height=«3» src=«ref-2_325784720-127.coolpic» v:shapes="_x0000_s5235">Разделим данные таблицы на две части. В первую часть попадут данные с 1-го по 6-й день, а во вторую часть – с 7-го по 12-й день работы. Рассчитаем по каждой половине среднее число дней и средние объемы продаж. По формулам:
<img width=«136» height=«54» src=«ref-2_325784847-480.coolpic» v:shapes="_x0000_i1040"> (19)
<img width=«174» height=«54» src=«ref-2_325785327-575.coolpic» v:shapes="_x0000_i1041"> (20)
найдем искомые значения для первой половины данных таблицы:
<img width=«195» height=«44» src=«ref-2_325785902-394.coolpic» v:shapes="_x0000_i1277">
<img width=«328» height=«42» src=«ref-2_325786296-579.coolpic» v:shapes="_x0000_i1278">
По формулам:
<img width=«272» height=«50» src=«ref-2_325786875-828.coolpic» v:shapes="_x0000_i1042"> (21)
<img width=«238» height=«51» src=«ref-2_325787703-508.coolpic» v:shapes="_x0000_i1043"> (22)
найдем искомые значения для второй половины данных таблицы:
<img width=«221» height=«41» src=«ref-2_325788211-459.coolpic» v:shapes="_x0000_i1279">
<img width=«349» height=«44» src=«ref-2_325788670-645.coolpic» v:shapes="_x0000_i1280">
В результате расчетов мы получили координаты двух точек А(3,5; 11,33) и В(9,5; 13,68). Построим эти точки, через них проведем прямую до пересечения с осью ординат (объем продаж) (рис. 4). Точка пересечения а0=9,95.
рис. 4 Определение параметра а0
Теперь определим значение параметра а1 по формуле:
<img width=«96» height=«48» src=«ref-2_325789315-418.coolpic» v:shapes="_x0000_i1044"> (23)
<img width=«185» height=«48» src=«ref-2_325789733-403.coolpic» v:shapes="_x0000_i1281">
В результате расчетов линейная модель <img width=«104» height=«28» src=«ref-2_325790136-434.coolpic» v:shapes="_x0000_i1045"> будет иметь следующий конкретный вид <img width=«134» height=«27» src=«ref-2_325790570-529.coolpic» v:shapes="_x0000_i1282">
продолжение
--PAGE_BREAK--7. Определение параметровлинейного тренда методом наименьших квадратов:
Чтобы найти параметры линейного тренда <img width=«102» height=«28» src=«ref-2_325791099-434.coolpic» v:shapes="_x0000_i1046">, необходимо решить систему нормальных уравнений:
<img width=«159» height=«65» src=«ref-2_325791533-839.coolpic» v:shapes="_x0000_i1047"> (24)
Для расчета параметров линейного тренда <img width=«104» height=«28» src=«ref-2_325790136-434.coolpic» v:shapes="_x0000_i1048"> методом МНК используем данные таблицы 4. Для проведения промежуточных расчетов построим таблицу 9 и проведем в ней необходимые расчеты.
Таблица 9
Найдем параметры а0и а1 подставив цифры из итоговой строки в формулы:
<img width=«450» height=«48» src=«ref-2_325792890-1041.coolpic» v:shapes="_x0000_i1283"> (25)
<img width=«393» height=«44» src=«ref-2_325793931-947.coolpic» v:shapes="_x0000_i1284"> (26)
В результате расчетов линейный тренд <img width=«109» height=«30» src=«ref-2_325794878-444.coolpic» v:shapes="_x0000_i1050"> примет конкретный вид
<img width=«128» height=«24» src=«ref-2_325795322-254.coolpic» v:shapes="_x0000_i1285">.
Рассчитаем значения линейного тренда для каждого момента времени подставив соответствующие значения t в уравнения
<img width=«128» height=«24» src=«ref-2_325795322-254.coolpic» v:shapes="_x0000_i1286">
и отобразим в таблице 10:
Таблица 10
Построим линейный тренд на графике (рис. 5):
<img width=«673» height=«410» src=«ref-2_325796084-35441.coolpic» v:shapes="_x0000_i1288">
На основе визуального анализа можно сделать вывод: что соответствие линейного тренда с трендом, который может иметь место во временном ряду очевидно. Ранее был сделан вывод, что во временном ряду присутствует тенденция в виде тренда, и он имеет линейный характер. На графике присутствует линейный тренд
Сравним параметры линейного тренда, вычисленные графическим методом – а0=9,95 и а1=0,394 и методом МНК – а0=9,929 и а1=0,397. Они достаточно близки.
Таблица 11
продолжение
--PAGE_BREAK--8. Определение параметровнелинейной трендовой модели:
8.1 Предварительным условием выбора параболы для описания исходных данных является относительное постоянство разностей второго порядка по абсолютной величине. Проверим, можно ли описать изменения данных с помощью параболы второго порядка. Для этого необходимо рассчитать абсолютный цепной прирост уровней временного ряда ∆tили разности первого порядка <img width=«21» height=«28» src=«ref-2_325832024-198.coolpic» v:shapes="_x0000_i1051">, разности второго порядка <img width=«20» height=«25» src=«ref-2_325832222-112.coolpic» v:shapes="_x0000_i1052">. Расчет будет производиться по сглаженным данным, так как их колебание меньше, чем колебание исходных данных.
Расчеты внесем в таблицу 12, с расчетами прироста первого (графа 4) и второго порядков (графа 5):
Таблица 12
Так как по сглаженным данным первый и двенадцатый уровни отсутствуют, а для второго уровня отсутствует предыдущий, рассчитаем абсолютный прирост уровней ряда с третьего по одиннадцатый.
Абсолютный цепной прирост третьего и последующих уровней ряда рассчитаем по формуле:
<img width=«108» height=«29» src=«ref-2_325832557-339.coolpic» v:shapes="_x0000_i1055"> (27)
Итак,
<img width=«266» height=«29» src=«ref-2_325832896-754.coolpic» v:shapes="_x0000_i1291">
<img width=«266» height=«29» src=«ref-2_325833650-735.coolpic» v:shapes="_x0000_i1292">
<img width=«266» height=«29» src=«ref-2_325834385-754.coolpic» v:shapes="_x0000_i1293">
<img width=«266» height=«29» src=«ref-2_325835139-776.coolpic» v:shapes="_x0000_i1294">
<img width=«267» height=«29» src=«ref-2_325835915-771.coolpic» v:shapes="_x0000_i1295">
<img width=«266» height=«29» src=«ref-2_325836686-780.coolpic» v:shapes="_x0000_i1296">
<img width=«264» height=«29» src=«ref-2_325837466-741.coolpic» v:shapes="_x0000_i1297">
<img width=«276» height=«29» src=«ref-2_325838207-771.coolpic» v:shapes="_x0000_i1298">
<img width=«277» height=«29» src=«ref-2_325838978-789.coolpic» v:shapes="_x0000_i1299">
Рассчитанные значения приростов занесем в графу 4 таблицы 12.
Условием выбора параболы для описания сглаженных данных является относительное постоянство разностей второго порядка по абсолютной величине
Определим их по формуле:
<img width=«108» height=«28» src=«ref-2_325839767-236.coolpic» v:shapes="_x0000_i1056"> (28)
Расчет начнем с четвертого уровня, так как он является первым, имеющим предыдущий ряд, а закончим одиннадцатым рядом.
Итак,
<img width=«259» height=«28» src=«ref-2_325840003-790.coolpic» v:shapes="_x0000_i1300">
<img width=«243» height=«28» src=«ref-2_325840793-683.coolpic» v:shapes="_x0000_i1301">
<img width=«225» height=«28» src=«ref-2_325841476-609.coolpic» v:shapes="_x0000_i1302">
<img width=«244» height=«28» src=«ref-2_325842085-689.coolpic» v:shapes="_x0000_i1303">
<img width=«252» height=«28» src=«ref-2_325842774-766.coolpic» v:shapes="_x0000_i1304">
<img width=«243» height=«28» src=«ref-2_325843540-704.coolpic» v:shapes="_x0000_i1305">
<img width=«251» height=«28» src=«ref-2_325844244-699.coolpic» v:shapes="_x0000_i1306">
<img width=«262» height=«28» src=«ref-2_325844943-789.coolpic» v:shapes="_x0000_i1307">
Анализ графы 5 показывает, что разности второго порядка по абсолютной величине относительно постоянны (0,2, 0,3 или 0,4). Это позволяет предположить, что для описания исходных данных может быть использована и парабола. Но эта трендовая модель используется тогда, когда значения экономического показателя, в среднем, вначале возрастают (убывают), достигают максимального (минимального) значения, а затем – убывают (возрастает). У нас же, на лицо, наличие возрастающей тенденции, скорее свойственной для степенного или показательного тренда.
8.2 Степенной тренд используется тогда, когда темпы роста экономического показателя, в среднем, либо постепенно возрастают, либо постепенно убывают. Рассмотрим возможность использования степенного тренда для описания исходных данных в таблице 13, рассчитав цепные темпы роста.
Таблица 13
Рассчитаем цепные темпы роста для каждого уровня ряда, начиная со второго, по формуле:
<img width=«143» height=«55» src=«ref-2_325845874-624.coolpic» v:shapes="_x0000_i1058"> (29)
Итак,
<img width=«267» height=«47» src=«ref-2_325846498-588.coolpic» v:shapes="_x0000_i1308"> <img width=«252» height=«47» src=«ref-2_325847086-572.coolpic» v:shapes="_x0000_i1309">
<img width=«267» height=«47» src=«ref-2_325847658-605.coolpic» v:shapes="_x0000_i1310">% <img width=«267» height=«47» src=«ref-2_325848263-628.coolpic» v:shapes="_x0000_i1311">
<img width=«267» height=«47» src=«ref-2_325848891-599.coolpic» v:shapes="_x0000_i1312"> <img width=«267» height=«47» src=«ref-2_325849490-593.coolpic» v:shapes="_x0000_i1313">
<img width=«260» height=«47» src=«ref-2_325850083-622.coolpic» v:shapes="_x0000_i1314"> <img width=«266» height=«47» src=«ref-2_325850705-638.coolpic» v:shapes="_x0000_i1315">
<img width=«274» height=«47» src=«ref-2_325851343-645.coolpic» v:shapes="_x0000_i1316"> <img width=«273» height=«47» src=«ref-2_325851988-618.coolpic» v:shapes="_x0000_i1317">
<img width=«266» height=«47» src=«ref-2_325852606-643.coolpic» v:shapes="_x0000_i1318">
Анализ третьей графы таблицы 13 показывает, что по мере роста объема сбыта темп сбыта продукции остается относительно постоянным. Это позволяет предположить, что для описания исходных данных использование степенного тренда менее предпочтительно, чем показательный тренд.
8.3 Показательный тренд используется при условии, когда экономический показатель имеет относительно постоянный цепной темп роста, который в среднем равен а1.
Оценим возможность описания данных, приведенных в таблице 13 с помощью показательного тренда. Так как цепные темпы роста, приведенные в третьей графе, относительно постоянны, то исходные данные могут быть описаны с помощью показательного тренда.
продолжение
--PAGE_BREAK--9. Рассчитаем параметры для показательного тренда:
С помощью таблицы 14 рассчитаем необходимые промежуточные данные:
Таблица 14
t
yt
t2
lnyt
tlnyt
1
2
3
4
6
1
10.2
1
2.322388
2.322388
2
10.8.
4
2.379546
4.759092
Продолжение таблицы 14
3
10.4
9
2.341806
7.025418
4
11.9
16
2.476538
9.906152
5
12.2
25
2.501436
12.507180
6
12.5
36
2.525729
15.154374
7
13.1
49
2.572612
18.008284
8
12.4
64
2.517696
20.141568
9
13.6
81
2.610070
23.490630
10
14.3
100
2.660259
26.602590
11
14.9
121
2.701361
29.714971
12
13.8
144
2.624669
31.496028
78
-
650
30.23411
201.128675
На основе данных итоговой строки определим параметры линеаризированной модели показательного тренда lna
и lna1:
<img width=«523» height=«53» src=«ref-2_325853249-1526.coolpic» v:shapes="_x0000_i1319"> (34)
<img width=«488» height=«52» src=«ref-2_325854775-1418.coolpic» v:shapes="_x0000_i1320"> (35)
Найдем значение a
, пропотенцировав lna. Тогда <img width=«256» height=«25» src=«ref-2_325856193-493.coolpic» v:shapes="_x0000_i1321">. Найдем a
1, пропотенцировав величину lna1. Тогда <img width=«245» height=«24» src=«ref-2_325856686-438.coolpic» v:shapes="_x0000_i1322">. В результате расчетов параметров показательный тренд будет иметь следующий конкретный вид:
<img width=«125» height=«25» src=«ref-2_325857124-258.coolpic» v:shapes="_x0000_i1323">
продолжение
--PAGE_BREAK--10.Сравним линейный и выбранный нелинейный – показательный тренды:
В одной системе координат, построим 3 графика: исходный временной ряд, линейную и показательную модели. Для построения графиков используем таблицу 15. Рассчитаем значения линейного и показательного трендов для каждого момента времени, подставив соответствующие значения t в уравнения <img width=«128» height=«24» src=«ref-2_325795322-254.coolpic» v:shapes="_x0000_i1324"> и <img width=«140» height=«29» src=«ref-2_325857636-556.coolpic» v:shapes="_x0000_i1325">
Таблица 15
Нарисуем график по исходным данным на рис.6:
<img width=«647» height=«410» src=«ref-2_325859002-19702.coolpic» v:shapes="_x0000_i1328">
На основе визуального анализа можно сделать вывод: что совпадение выбранного нелинейного тренда с реальным трендом временного ряда очевидно.
11. Сравнение трендовых моделей с помощью критерия наименьшей суммы квадратов отклонений:
Выбранные тренды имеют одинаковое число параметров, поэтому для расчета критерия наименьшей суммы квадратов отклонений выберем формулу
<img width=«172» height=«28» src=«ref-2_325878704-644.coolpic» v:shapes="_x0000_i1059">. (36)
Для проведения промежуточных расчетов используем таблицу 16:
Таблица 16
В графу 5 внесем отклонение фактических данных от линейного тренда: для t=1 <img width=«76» height=«28» src=«ref-2_325880584-316.coolpic» v:shapes="_x0000_i1331">=(10,2–10,326)2=0,016, для t=2 <img width=«77» height=«29» src=«ref-2_325880900-319.coolpic» v:shapes="_x0000_i1332">=(10,8–10,723)2=0,006 и т.д. Аналогичным образом в графу 6 впишем отклонения фактических данных от показательного тренда, и подведем итог.
Сравним значения критерия наименьшей суммы квадратов отклонений для линейного и показательного трендов. Для линейного тренда критерий равен 2,789 (итог графы 5), а для показательного тренда 3,109 (итог графы 6). Коэффициент для линейного тренда меньше, чем для показательного, поэтому линейный тренд лучше аппроксимирует исходные данные. Следовательно, для прогнозирования необходимо взять линейный тренд
.
продолжение
--PAGE_BREAK--12. Оценка адекватности выбранной трендовой модели:
Чтобы оценить адекватность выбранной линейной трендовой модели теоретическому тренду временного ряда, найдем разность е
tмежду исходными данными у
tи нашей трендовой моделью <img width=«19» height=«26» src=«ref-2_325881219-138.coolpic» v:shapes="_x0000_i1062"> по формуле
<img width=«94» height=«29» src=«ref-2_325881357-319.coolpic» v:shapes="_x0000_i1063"> (37)
результаты занесем в таблицу 17.
Таблица 17
Построим график ряда отклонения е
t
по данным графы 4 на рис.7:
<img width=«639» height=«345» src=«ref-2_325882269-26107.coolpic» v:shapes="_x0000_i1335">
Визуальный анализ показывает, что колебание величины
е
t
не содержит элементов тенденции, т.е. носит случайный характер.
Оценим адекватность выбранной модели тренда исходному ряду на основе анализа данных ряда отклонений е
t. Величина е
tдолжна отвечать следующим четырем условиям (требованиям):
Условие 1. Колебание величины е
t
должно носить случайный характер.
Проверим условие 1 с помощью критерия поворотных точек.
Величина е
tсчитается поворотной, если она соответствует одному из двух условий:
е
t-1< еt>еt+1 или е
t-1> еt<еt+1
Для этого на базе данных графы 4 определим поворотные точки и в графе 5 проставим соответствующие им значения. Рассматриваемые точки обозначим Рt. Тем точкам, которые будут поворотными, присвоим значение Рt=1. А тем точкам, которые не будут поворотными, присвоим значение Рt=0. результаты занесем в таблицу 18.
Таблица 18
Определим общее число поворотных точек по формуле
<img width=«303» height=«27» src=«ref-2_325908969-537.coolpic» v:shapes="_x0000_i1338">
Результат занесем в итоговую строку графы 6.
Для проверки выполнения условия 1 выдвинем нулевую гипотезу Н0: колебание величины еt
носит случайный характер
.
Для проверки этой гипотезы определим математическое ожидание числа поворотных точек по формуле:
<img width=«291» height=«44» src=«ref-2_325909506-1062.coolpic» v:shapes="_x0000_i1064">; (38)
и его дисперсию по формуле:
<img width=«298» height=«45» src=«ref-2_325910568-966.coolpic» v:shapes="_x0000_i1065"> (39)
Кроме того, для проверки нулевой гипотезы используем вероятность, равную 95%, которой соответствует коэффициент доверия t=1,96. С помощью формулы
<img width=«303» height=«27» src=«ref-2_325911534-836.coolpic» v:shapes="_x0000_i1066"> (40)
проверим нулевую гипотезу, подставив в нее значения М(Р), D(P), t:
<img width=«306» height=«27» src=«ref-2_325912370-554.coolpic» v:shapes="_x0000_i1339"> или 4,029<7<9,305
Расчет показывает, общее число поворотных точек – 7 находится в требуемом интервале. Это позволяет сделать следующий вывод: с вероятностью 0,95 (95%) колебание величины е
t
носит случайный характер и, следовательно, отвечает данному условию.
Условие 2. Распределение величины е
t
соответствует нормальному распределению.
Проверим распределениееtна соответствие нормальному распределению. Вначале определим среднее квадратическое отклонение по формуле:
<img width=«221» height=«50» src=«ref-2_325912924-610.coolpic» v:shapes="_x0000_i1340"> (41)
для определения расчетного значения критерия RSриз графы 4 таблицы 18 найдем максимальное emax=0,604 и минимальное emin=−0,893 значения. Расчетное значение критерия RSр найдем по формуле:
<img width=«320» height=«53» src=«ref-2_325913534-1292.coolpic» v:shapes="_x0000_i1341"> (42)
Следующим шагом проверки условия 2 является нахождение табличного значения RS-критерия – RSTпо приложению 3 (2, стр. 71). В таблице приводятся нижнее и верхнее значения RS-критерия для n=10 и n=20; а у нас n=12. Для нахождения нижнего и верхнего значений RS-критерия для n=12 используем линейную интерполяцию. Найдем величину RS12н:
· Увеличение RSnн при изменеие n на 2 найдем по формуле:
<img width=«393» height=«45» src=«ref-2_325914826-1447.coolpic» v:shapes="_x0000_i1067">
· Значение RSnн при t=12 найдем по формуле:
<img width=«336» height=«23» src=«ref-2_325916273-944.coolpic» v:shapes="_x0000_i1068">
Найдем величину RSnв:
· Увеличение RSnвпри изменеие n на 2 найдем по формуле:
<img width=«393» height=«45» src=«ref-2_325917217-1479.coolpic» v:shapes="_x0000_i1069">
· Значение RSnвпри n=12 найдем по формуле:
<img width=«343» height=«23» src=«ref-2_325918696-1029.coolpic» v:shapes="_x0000_i1070">
Выдвинем нулевую гипотезу Н0: величина еtсоответствует нормальному распределению. Сопоставим по формуле:
RSnн< RSр< RSnв (43)
расчетное значение критерия RSр с табличным –RSТ… Сопоставление показывает, что RSр попадает в интервал, определяемый нижним и верхним табличными значениями RS-критерия, т.е. 2,772<2,835<3,978
Это позволяет нам сделать следующий вывод: с вероятностью 95% нулевая гипотезе принимается, т.е. величина еt
соответствует нормальному распределению и, следовательно, отвечает условию 2.
Условие 3. Математическое ожидание величины е
t
равно нулю
Для проверки данного условия выдвинем нулевую гипотезу Н0: математическое ожидание et=0.
Для промежуточных расчетов используем таблицу 19:
Таблица 19
Вначале определим среднюю арифметическую величину еt, использовав итог графы 3 таблицы 19;
<img width=«192» height=«43» src=«ref-2_325920266-402.coolpic» v:shapes="_x0000_i1343"> (44)
Затем рассчитаем и внесем в графу 4 табл. 19 квадрат отклонения фактического значения еtот ее среднего значения. Так, дляt=1
<img width=«64» height=«28» src=«ref-2_325920668-291.coolpic» v:shapes="_x0000_i1072">=(−0,126−(-0,0012))2 =0,016 и т.д.
Далее определим среднее квадратическое отклонение, используя итог графы 4 таблицы 19 по формуле:
<img width=«242» height=«51» src=«ref-2_325920959-636.coolpic» v:shapes="_x0000_i1344"> (45)
Теперь найдем расчетное значение величиныtpпо формуле:
<img width=«298» height=«48» src=«ref-2_325921595-619.coolpic» v:shapes="_x0000_i1345"> (46)
Чтобы найти табличное значение величины tT, зададимся уровнем значимости а=0,05, относительно которого определим доверительную вероятность γ=1−0,05=0,95, а также число степеней свободы k
=12–1=11. Теперь, зная γ и k
,определим tT по Стьюденту (приложение 2 (2, стр. 70)); tТ =2,201.
Сопоставим расчетное tp=-0,0083 и табличное tT=2,201 значения:
t
p
<
t
T или -0,0083 < 2,201.
Сопоставление показывает, расчетное значение меньше табличного.
Это позволяет нам сделать следующий вывод: с вероятность 0,95 (95%) нулевая гипотеза принимается и мы может утверждать: математическое ожидание е
t
=0.
Условие 4. Независимость членов ряда друг от друга
Оценим наличие автокорреляции в ряде данных еtс помощью критерия Дарбина-Уотсона. Вначале в графу 3 таблицы 20 внесем квадраты величины еt, а в графе 4 – квадраты разниц между текущим и предыдущим значениями еt. Так, для t=1 мы не можем найти требуемое значение квадрата разницы, так как у нас нет значения е0. А для t=2
(е2–е1)2=(0.077 – (–0.126))2=0.041 и т.д. (47)
Таблица 20
Рассчитаем итоговые значения граф 3 и 4 и по ним определим расчетное значение критерия Дарбина — Уотсонаdpпо формуле:
<img width=«247» height=«53» src=«ref-2_325922913-1017.coolpic» v:shapes="_x0000_i1346"> (48)
Расчетное значение критерия Дарбина–Уотсона оказалось больше двух, (следовательно, коэффициент попал в область отрицательной автокорреляции), поэтому пересчитаем его для области с положительной автокорреляцией по формуле:
<img width=«226» height=«29» src=«ref-2_325923930-645.coolpic» v:shapes="_x0000_i1347"> (49)
Найдем табличное значение критерия Дарбина–Уотсона dTпри n=12, и числе факторов в используемой трендовой модели V=1 (в линейном тренде один фактор –время).
При n=12 и V=1 в приложении 4 (2, стр. 71) находим табличное значение критерия Дарбина–Уотсона dT. Однако у нас n=12, а в таблице наименьшее значение n=15, поэтому возьмем табличное значение критерия Дарбина–Уотсона для n= 15. Его нижнее значение равно d1=1,08, а верхнее d2=1,36:
Сопоставим расчетное (1,786) и табличное (1,08; 1,36) значения критерия Дарбина–Уотсона. При этом могут возникнуть три ситуации:
1) d
p
<
d
1, что будет говорить о наличии в ряде автокорреляции
2) d
p
>
d
2, что будет говорить об отсутствии в ряде автокорреляции;
3) d
1
≤
d
p
≤
d
2, что будет говорить о необходимости дополнительной
проверки наличия в ряде автокорреляции.
Расчетное значение больше верхнего табличного, т.е. возникает вторая ситуация, когда d
p
>
d
2или 1,786>1,36.
С учетом этого мы может сделать вывод: с вероятностью 0,95(95%) в ряде е
t
отсутствует автокорреляция.
Поскольку соблюдаются 4 условия:
1.
Условие 1. Колебание величины еt должно носить случайный характер. Это условие означает, что колебание (изменение) величины еt не содержит элементов тенденции.
2.
Условие 2. Распределение величины еt соответствует нормальному распределению.
3.
Условие 3. Математическое ожидание величины еt равно нулю
4.
Условие 4. Независимость членов ряда друг от друга. Это условие означает отсутствие автокорреляции во временном ряде еt
Можно утверждать, что выбранная трендовая модель: <img width=«143» height=«27» src=«ref-2_325924575-570.coolpic» v:shapes="_x0000_i1348"> адекватна тенденции, имеющей место во временном ряде.
продолжение
--PAGE_BREAK--13. Прогнозирование на основе трендовой модели:
13.1 Точечный прогноз:
Определим точечный прогноз на 13-й день. Из условия задачи вытекает: период основания прогноза n=12, а период упреждения прогноза τ=1. Одновременно определим уровень значимости а=0,05. Рассчитаем точечный прогноз по формуле:
<img width=«406» height=«30» src=«ref-2_325925145-1242.coolpic» v:shapes="_x0000_i1349"> (50)
13.2 Интервальный прогноз:
Для расчета интервального прогноза предварительно определим табличное значение критерия Стьюдента с уровнем значимостиа и числом степеней свободы k
=
n−2. Так как мы выбрали, а=0,05, доверительная вероятность γ=1−а=1−0,05=0,95, а число степеней свободы k=n−2=12−2=10. По приложению 2 (2, стр. 70) при γ=0,95 и k=10 табличное значение критерия Стьюдента tT=2,228. Найдем стандартную ошибку тренда, по формуле:
<img width=«261» height=«54» src=«ref-2_325926387-1171.coolpic» v:shapes="_x0000_i1350"> (51)
Определим интервальный прогноз по формуле:
<img width=«465» height=«115» src=«ref-2_325927558-2818.coolpic» v:shapes="_x0000_i1351"> (52)
Отсюда верхняя граница прогнозного интервала 15,09+1,381=16,471, а нижняя 15,09–1,381=13,709. Таким образом, прибыль от продаж на фабрике на 13-й день с вероятностью γ=0,95 будет расположена в интервале от 13,709 … 16,471.
Для расчета интервального прогноза с использованием формулы
<img width=«148» height=«33» src=«ref-2_325930376-566.coolpic» v:shapes="_x0000_i1076"> (53)
определим К. Согласно исходным данным число уровней ряда n=12, а период упреждения прогноза τ=1, поэтому К= 2,1274 (2, стр. 73). Подставим найденное К в формулу и получим интервальный прогноз
<img width=«418» height=«31» src=«ref-2_325930942-1291.coolpic» v:shapes="_x0000_i1352">
Отсюда верхняя граница прогнозного интервала 15,09+1,123=16,213, а нижняя 15,09–1,123=13,967. Таким образом, прогноз оборота магазина на 13-й день с вероятностью γ=0,9 будет расположен в интервале 13,967 … 16,213.
Верхняя и нижняя границы прогнозного интервала отличаются от полученных ранее. Причиной этого является то, что при расчете по 1-й формуле был использован уровень значимости а=0,05, откуда доверительная вероятность γ=0,95, а при расчете по 2-й формуле была использована величина К, которая в приложении 6 (2, стр. 73) рассчитана относительно уровня значимости а=0,1, откуда доверительная вероятность равна 0,9.
Раздел 3 «Прогнозированиена основе сезонного цикла временного ряда»
Исходные данные для Варианта №23:
1. Таблица с исходными данными:
Таблица 21. Объем реализации продукции фирмы АО «Лен» (усл.ед.)
продолжение
--PAGE_BREAK-- 2. Отобразим исходные данные на рис.8:
Вначале в одной системе координат построим два графика – один по исходным данным, другой график линейной трендовой модели <img width=«144» height=«24» src=«ref-2_325932233-291.coolpic» v:shapes="_x0000_i1077"> (рис. 8) и проведем его визуальный анализ.
Визуальный анализ графика временного ряда показывает, что исходный ряд содержит сезонную компоненту, так как характер колебания ряда стабильно повторяется из года в год и имеет приблизительно одинаковый характер изменения. Можно предположить о наличии тенденции в виде тренда. И он может быть описан линейным трендом <img width=«144» height=«24» src=«ref-2_325932233-291.coolpic» v:shapes="_x0000_i1078">
--PAGE_BREAK--4.Сезонное прогнозирование:
Осуществить прогноз сезонного цикла на основе линейной трендовой модели <img width=«144» height=«24» src=«ref-2_325932233-291.coolpic» v:shapes="_x0000_i1088">. Параметры прогнозирования отображены в таблице 23.
Таблица 23
Месяц
Год
2009
Январь
9979
Февраль
9222
Март
10526
Апрель
11740
Май
12466
Июнь
13801
Июль
14604
Август
14638
Сентябрь
13709
Октябрь
11790
Ноябрь
9849
Таким образом, согласно таблице необходимо осуществить прогнозирование с января по ноябрь 2009 года с использованием мультипликативной модели прогнозирования.
Определим объем реализации продукции на основе линейного тренда в январе 2006, 2007 и 2008 годов:
· в январе 2006 года: <img width=«351» height=«24» src=«ref-2_326011429-1042.coolpic» v:shapes="_x0000_i1089">
· в январе 2007 года: <img width=«359» height=«24» src=«ref-2_326012471-1109.coolpic» v:shapes="_x0000_i1090">
· в январе 2008 года: <img width=«367» height=«24» src=«ref-2_326013580-1119.coolpic» v:shapes="_x0000_i1091">
Определим абсолютное отклонение фактических данных от тренда по формуле:
<img width=«67» height=«58» src=«ref-2_326014699-311.coolpic» v:shapes="_x0000_i1354"> (54)
Рассчитаем эти отклонения:
· в январе 2006 года: <img width=«227» height=«53» src=«ref-2_326015010-1024.coolpic» v:shapes="_x0000_i1355">
· в январе 2007 года: <img width=«230» height=«53» src=«ref-2_326016034-1045.coolpic» v:shapes="_x0000_i1356">
· в январе 2008 года: <img width=«236» height=«53» src=«ref-2_326017079-1092.coolpic» v:shapes="_x0000_i1357">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для января по формуле:
<img width=«318» height=«45» src=«ref-2_326018171-640.coolpic» v:shapes="_x0000_i1358"> (55)
Перед определением сезонного тренда найдем значение тренда в январе 2009 года. Исходный ряд содержит данные за три года (k=3), период упреждения прогноза равен одному году (τ=1), январь имеет номер 1, поэтому январь 2009 года будет иметь номер <img width=«41» height=«25» src=«ref-2_326018811-135.coolpic» v:shapes="_x0000_i1092">. Отсюда объем реализации продукции в январе 2009 года по тренду: <img width=«397» height=«27» src=«ref-2_326018946-615.coolpic» v:shapes="_x0000_i1359">
В итоге рассчитаем точечный прогноз продукции в январе по формуле:
<img width=«316» height=«30» src=«ref-2_326019561-1054.coolpic» v:shapes="_x0000_i1360"> (56)
Отразим спрогнозированное значение на январь на рис 8.
Определим объем реализации продукции на основе линейного тренда в феврале 2006, 2007 и 2008 годов:
· в феврале 2006 года: <img width=«356» height=«24» src=«ref-2_326020615-579.coolpic» v:shapes="_x0000_i1093">
· в феврале 2007 года: <img width=«359» height=«24» src=«ref-2_326021194-592.coolpic» v:shapes="_x0000_i1094">
· в феврале 2008 года: <img width=«368» height=«24» src=«ref-2_326021786-604.coolpic» v:shapes="_x0000_i1095">
Рассчитаем отклонения:
· в феврале 2006 года: <img width=«217» height=«53» src=«ref-2_326022390-969.coolpic» v:shapes="_x0000_i1361">
· в феврале 2007 года: <img width=«219» height=«53» src=«ref-2_326023359-991.coolpic» v:shapes="_x0000_i1362">
· в феврале 2008 года: <img width=«227» height=«53» src=«ref-2_326024350-1069.coolpic» v:shapes="_x0000_i1363">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для февраля по формуле (55):
<img width=«320» height=«45» src=«ref-2_326025419-627.coolpic» v:shapes="_x0000_i1364">
Объем реализации продукции в феврале 2009 года по тренду:
<img width=«426» height=«27» src=«ref-2_326026046-680.coolpic» v:shapes="_x0000_i1096">
В итоге рассчитаем точечный прогноз продукции в феврале по формуле (56):
<img width=«309» height=«30» src=«ref-2_326026726-1019.coolpic» v:shapes="_x0000_i1365">
Отразим спрогнозированное значение за февраль на рис 8.
Определим объем реализации продукции на основе линейного тренда в марте 2006, 2007 и 2008 годов:
· в марте 2006 года: <img width=«359» height=«24» src=«ref-2_326027745-587.coolpic» v:shapes="_x0000_i1097">
· в марте 2007 года: <img width=«376» height=«24» src=«ref-2_326028332-614.coolpic» v:shapes="_x0000_i1098">
· в марте 2008 года: <img width=«376» height=«24» src=«ref-2_326028946-616.coolpic» v:shapes="_x0000_i1099">
Рассчитаем отклонения:
· в марте 2006 года: <img width=«227» height=«53» src=«ref-2_326029562-1027.coolpic» v:shapes="_x0000_i1366">
· в марте 2007 года: <img width=«236» height=«53» src=«ref-2_326030589-1114.coolpic» v:shapes="_x0000_i1367">
· в марте 2008 года: <img width=«236» height=«53» src=«ref-2_326031703-1096.coolpic» v:shapes="_x0000_i1368">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для марта по формуле (55):
· <img width=«318» height=«45» src=«ref-2_326032799-629.coolpic» v:shapes="_x0000_i1369">
Объем реализации продукции в марте 2009 года по тренду:
<img width=«436» height=«27» src=«ref-2_326033428-699.coolpic» v:shapes="_x0000_i1100">
В итоге рассчитаем точечный прогноз продукции в марте по формуле (56):
<img width=«324» height=«30» src=«ref-2_326034127-1071.coolpic» v:shapes="_x0000_i1370">
Отразим спрогнозированное значение за март на рис 8.
Определим объем реализации продукции на основе линейного тренда в апреле 2006, 2007 и 2008 годов:
· в апреле 2006 года: <img width=«351» height=«24» src=«ref-2_326035198-577.coolpic» v:shapes="_x0000_i1101">
· в апреле 2007 года: <img width=«367» height=«24» src=«ref-2_326035775-609.coolpic» v:shapes="_x0000_i1102">
· в апреле 2008 года: <img width=«367» height=«24» src=«ref-2_326036384-601.coolpic» v:shapes="_x0000_i1103">
Рассчитаем отклонения:
· в апреле 2006 года: <img width=«206» height=«53» src=«ref-2_326036985-950.coolpic» v:shapes="_x0000_i1371">
· в апреле 2007 года: <img width=«227» height=«53» src=«ref-2_326037935-1074.coolpic» v:shapes="_x0000_i1372">
· в апреле 2008 года: <img width=«222» height=«53» src=«ref-2_326039009-997.coolpic» v:shapes="_x0000_i1373">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для апреля по формуле (55):
<img width=«297» height=«45» src=«ref-2_326040006-591.coolpic» v:shapes="_x0000_i1374">
Объем реализации продукции в апреле 2009 года по тренду:
<img width=«426» height=«27» src=«ref-2_326040597-679.coolpic» v:shapes="_x0000_i1104">
В итоге рассчитаем точечный прогноз продукции в апреле по формуле (56):
<img width=«307» height=«30» src=«ref-2_326041276-952.coolpic» v:shapes="_x0000_i1375">
Отразим спрогнозированное значение за апрель на рис 8.
Определим объем реализации продукции на основе линейного тренда в мае 2006, 2007 и 2008 годов:
· в мае 2006 года: <img width=«363» height=«24» src=«ref-2_326042228-592.coolpic» v:shapes="_x0000_i1105">
· в мае 2007 года: <img width=«376» height=«24» src=«ref-2_326042820-614.coolpic» v:shapes="_x0000_i1106">
· в мае 2008 года: <img width=«376» height=«24» src=«ref-2_326043434-619.coolpic» v:shapes="_x0000_i1107">
Рассчитаем отклонения:
· в мае 2006 года: <img width=«216» height=«53» src=«ref-2_326044053-980.coolpic» v:shapes="_x0000_i1376">
· в мае 2007 года: <img width=«223» height=«53» src=«ref-2_326045033-1021.coolpic» v:shapes="_x0000_i1377">
· в мае 2008 года: <img width=«225» height=«53» src=«ref-2_326046054-1001.coolpic» v:shapes="_x0000_i1378">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для мая по формуле (55):
<img width=«271» height=«45» src=«ref-2_326047055-547.coolpic» v:shapes="_x0000_i1379">
Объем реализации продукции в мае 2009 года по тренду:
<img width=«433» height=«27» src=«ref-2_326047602-692.coolpic» v:shapes="_x0000_i1108">
В итоге рассчитаем точечный прогноз продукции в мае по формуле (56):
<img width=«305» height=«30» src=«ref-2_326048294-980.coolpic» v:shapes="_x0000_i1380">
Отразим спрогнозированное значение за май на рис 8.
Определим объем реализации продукции на основе линейного тренда в июне 2006, 2007 и 2008 годов:
· в июне 2006 года: <img width=«356» height=«24» src=«ref-2_326049274-584.coolpic» v:shapes="_x0000_i1109">
· в июне 2007 года: <img width=«367» height=«24» src=«ref-2_326049858-599.coolpic» v:shapes="_x0000_i1110">
· в июне 2008 года: <img width=«359» height=«24» src=«ref-2_326050457-587.coolpic» v:shapes="_x0000_i1111">
Рассчитаем отклонения:
· в июне 2006 года: <img width=«207» height=«53» src=«ref-2_326051044-917.coolpic» v:shapes="_x0000_i1381">
· в июне 2007 года: <img width=«216» height=«53» src=«ref-2_326051961-955.coolpic» v:shapes="_x0000_i1382">
· в июне 2008 года: <img width=«202» height=«53» src=«ref-2_326052916-879.coolpic» v:shapes="_x0000_i1383">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для июня по формуле (55):
· <img width=«280» height=«45» src=«ref-2_326053795-547.coolpic» v:shapes="_x0000_i1384">
Объем реализации продукции в июне 2009 года по тренду:
<img width=«426» height=«27» src=«ref-2_326054342-679.coolpic» v:shapes="_x0000_i1112">
В итоге рассчитаем точечный прогноз продукции в июне по формуле (56):
<img width=«305» height=«30» src=«ref-2_326055021-977.coolpic» v:shapes="_x0000_i1385">
Отразим спрогнозированное значение за июнь на рис 8.
Определим объем реализации продукции на основе линейного тренда в июле 2006, 2007 и 2008 годов:
· в июле 2006 года: <img width=«363» height=«24» src=«ref-2_326055998-589.coolpic» v:shapes="_x0000_i1113">
· в июле 2007 года: <img width=«375» height=«24» src=«ref-2_326056587-617.coolpic» v:shapes="_x0000_i1114">
· в июле 2008 года: <img width=«375» height=«24» src=«ref-2_326057204-609.coolpic» v:shapes="_x0000_i1115">
Рассчитаем отклонения:
· в июле 2006 года: <img width=«206» height=«53» src=«ref-2_326057813-996.coolpic» v:shapes="_x0000_i1386">
· в июле 2007 года: <img width=«225» height=«53» src=«ref-2_326058809-1074.coolpic» v:shapes="_x0000_i1387">
· в июле 2008 года: <img width=«225» height=«53» src=«ref-2_326059883-977.coolpic» v:shapes="_x0000_i1388">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для июля по формуле (55):
· <img width=«281» height=«45» src=«ref-2_326060860-558.coolpic» v:shapes="_x0000_i1389">
Объем реализации продукции в июле 2009 года по тренду:
<img width=«435» height=«27» src=«ref-2_326061418-696.coolpic» v:shapes="_x0000_i1116">
В итоге рассчитаем точечный прогноз продукции в июле по формуле (56):
<img width=«317» height=«30» src=«ref-2_326062114-1037.coolpic» v:shapes="_x0000_i1390">
Отразим спрогнозированное значение за июль на рис 8.
Определим объем реализации продукции на основе линейного тренда в августе 2006, 2007 и 2008 годов:
· в августе 2006 года: <img width=«351» height=«24» src=«ref-2_326063151-578.coolpic» v:shapes="_x0000_i1117">
· в августе 2007 года: <img width=«368» height=«24» src=«ref-2_326063729-605.coolpic» v:shapes="_x0000_i1118">
· в августе 2008 года: <img width=«367» height=«24» src=«ref-2_326064334-601.coolpic» v:shapes="_x0000_i1119">
Рассчитаем отклонения:
· в августе 2006 года: <img width=«203» height=«53» src=«ref-2_326064935-979.coolpic» v:shapes="_x0000_i1391">
· в августе 2007 года: <img width=«216» height=«53» src=«ref-2_326065914-962.coolpic» v:shapes="_x0000_i1392">
· в августе 2008 года: <img width=«215» height=«53» src=«ref-2_326066876-912.coolpic» v:shapes="_x0000_i1393">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для августа по формуле (55):
· <img width=«283» height=«45» src=«ref-2_326067788-564.coolpic» v:shapes="_x0000_i1394">
Объем реализации продукции в августе 2009 года по тренду:
<img width=«426» height=«27» src=«ref-2_326068352-680.coolpic» v:shapes="_x0000_i1120">
В итоге рассчитаем точечный прогноз продукции в августе по формуле (56):
<img width=«305» height=«30» src=«ref-2_326069032-992.coolpic» v:shapes="_x0000_i1395">
Отразим спрогнозированное значение за август на рис 8.
Определим объем реализации продукции на основе линейного тренда в сентябре 2006, 2007 и 2008 годов:
· в сентябре 2006 года: <img width=«359» height=«24» src=«ref-2_326070024-589.coolpic» v:shapes="_x0000_i1121">
· в сентябре 2007 года: <img width=«375» height=«24» src=«ref-2_326070613-615.coolpic» v:shapes="_x0000_i1122">
· в сентябре 2008 года: <img width=«376» height=«24» src=«ref-2_326071228-611.coolpic» v:shapes="_x0000_i1123">
Рассчитаем отклонения:
· в сентябре 2006 года: <img width=«215» height=«53» src=«ref-2_326071839-972.coolpic» v:shapes="_x0000_i1396">
· в сентябре 2007 года: <img width=«225» height=«53» src=«ref-2_326072811-1022.coolpic» v:shapes="_x0000_i1397">
· в сентябре 2008 года: <img width=«225» height=«53» src=«ref-2_326073833-1025.coolpic» v:shapes="_x0000_i1398">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для сентября по формуле (55):
· <img width=«271» height=«45» src=«ref-2_326074858-543.coolpic» v:shapes="_x0000_i1399">
Объем реализации продукции в сентябре 2009 года по тренду:
<img width=«433» height=«27» src=«ref-2_326075401-699.coolpic» v:shapes="_x0000_i1124">
В итоге рассчитаем точечный прогноз продукции в сентябре по формуле (56):
<img width=«308» height=«30» src=«ref-2_326076100-966.coolpic» v:shapes="_x0000_i1400">
Отразим спрогнозированное значение за сентябрь на рис 8.
Определим объем реализации продукции на основе линейного тренда в октябре 2006, 2007 и 2008 годов:
· в октябре 2006 года: <img width=«347» height=«24» src=«ref-2_326077066-570.coolpic» v:shapes="_x0000_i1125">
· в октябре 2007 года: <img width=«373» height=«24» src=«ref-2_326077636-609.coolpic» v:shapes="_x0000_i1126">
· в октябре 2008 года: <img width=«368» height=«24» src=«ref-2_326078245-605.coolpic» v:shapes="_x0000_i1127">
Рассчитаем отклонения:
· в октябре 2006 года: <img width=«185» height=«53» src=«ref-2_326078850-863.coolpic» v:shapes="_x0000_i1401">
· в октябре 2007 года: <img width=«219» height=«53» src=«ref-2_326079713-1011.coolpic» v:shapes="_x0000_i1402">
· в октябре 2008 года: <img width=«217» height=«53» src=«ref-2_326080724-1014.coolpic» v:shapes="_x0000_i1403">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для октября по формуле (55):
· <img width=«282» height=«45» src=«ref-2_326081738-569.coolpic» v:shapes="_x0000_i1404">
Объем реализации продукции в октябре 2009 года по тренду:
<img width=«432» height=«27» src=«ref-2_326082307-687.coolpic» v:shapes="_x0000_i1128">
В итоге рассчитаем точечный прогноз продукции в октябре по формуле (56):
<img width=«323» height=«30» src=«ref-2_326082994-1020.coolpic» v:shapes="_x0000_i1405">
Отразим спрогнозированное значение за октябрь на рис 8.
Определим объем реализации продукции на основе линейного тренда в ноябре 2006, 2007 и 2008 годов:
· в ноябре 2006 года: <img width=«367» height=«24» src=«ref-2_326084014-600.coolpic» v:shapes="_x0000_i1129">
· в ноябре 2007 года: <img width=«375» height=«24» src=«ref-2_326084614-616.coolpic» v:shapes="_x0000_i1130">
· в ноябре 2008 года: <img width=«376» height=«24» src=«ref-2_326085230-609.coolpic» v:shapes="_x0000_i1131">
Рассчитаем отклонения:
· в ноябре 2006 года: <img width=«218» height=«53» src=«ref-2_326085839-1022.coolpic» v:shapes="_x0000_i1406">
· в ноябре 2007 года: <img width=«195» height=«45» src=«ref-2_326086861-489.coolpic» v:shapes="_x0000_i1407">
· в ноябре 2008 года: <img width=«227» height=«53» src=«ref-2_326087350-1010.coolpic» v:shapes="_x0000_i1408">
На основе отклонений фактических данных от тренда определим среднее значение абсолютного отклонения, то есть сезонную компоненту для ноября по формуле (55):
· <img width=«281» height=«45» src=«ref-2_326088360-563.coolpic» v:shapes="_x0000_i1409">
Объем реализации продукции в ноябре 2009 года по тренду:
<img width=«441» height=«27» src=«ref-2_326088923-699.coolpic» v:shapes="_x0000_i1132">
В итоге рассчитаем точечный прогноз продукции в ноябре по формуле (56):
<img width=«312» height=«30» src=«ref-2_326089622-1015.coolpic» v:shapes="_x0000_i1410">
Отразим спрогнозированное значение за ноябрь на рис 8.
продолжение
--PAGE_BREAK--Раздел 4 «Прогнозированиес помощью метода экспоненциального сглаживания»
Исходные данные для Варианта №23:
1. Таблица с исходными данными:
Таблица 24. Курсы акций АО «Московская швея»
Таблица 24
2. Отобразим исходные данные на рис.9:
<img width=«644» height=«410» src=«ref-2_326090637-23635.coolpic» v:shapes="_x0000_i1411">
По графику исходных данных интуитивно можно предположить, что курс акций на 13 день составит от 666-675.
3. Метод экспоненциального сглаживания:
Рассчитаем при помощи метода экспоненциального сглаживания прогнозные оценки на 7-, 8-, 9-, 10-, 11-, 12-ый дни
Параметр сглаживания, а=0,2
Рассчитаем прогноз курсов акций с параметрами сглаживания, а=0,1 по формуле
<img width=«175» height=«29» src=«ref-2_326114272-650.coolpic» v:shapes="_x0000_i1133"> (57)
· На 7 день:<img width=«367» height=«25» src=«ref-2_326114922-566.coolpic» v:shapes="_x0000_i1412">
· На 8 день:<img width=«380» height=«25» src=«ref-2_326115488-574.coolpic» v:shapes="_x0000_i1413">
· На 9 день: <img width=«400» height=«25» src=«ref-2_326116062-616.coolpic» v:shapes="_x0000_i1414">
· На 10 день:<img width=«410» height=«25» src=«ref-2_326116678-630.coolpic» v:shapes="_x0000_i1415">
· На 11 день:<img width=«423» height=«25» src=«ref-2_326117308-641.coolpic» v:shapes="_x0000_i1416">
· На 12 день:<img width=«418» height=«24» src=«ref-2_326117949-639.coolpic» v:shapes="_x0000_i1417">
· На 13 день:<img width=«421» height=«24» src=«ref-2_326118588-634.coolpic» v:shapes="_x0000_i1418">
Занесем полученные прогнозное значение в последнюю строчку графы 2 и в предпоследнюю строчку графы 3 таблицы 25 (например, расчет, сделанный седьмого дня на восьмой равен 657,8 – это значение занесем в последнюю строчку графы 3 и в предпоследнюю строчку графы 4).
Параметр сглаживания, а=0,3
Рассчитаем прогноз курсов акций с параметрами сглаживания, а=0,3 по формуле (57):
· На 7день:<img width=«365» height=«25» src=«ref-2_326119222-572.coolpic» v:shapes="_x0000_i1419">
· На 8 день:<img width=«379» height=«25» src=«ref-2_326119794-573.coolpic» v:shapes="_x0000_i1420">
· На 9 день:<img width=«399» height=«25» src=«ref-2_326120367-604.coolpic» v:shapes="_x0000_i1421">
· На 10 день:<img width=«407» height=«25» src=«ref-2_326120971-620.coolpic» v:shapes="_x0000_i1422">
· На 11 день:<img width=«418» height=«25» src=«ref-2_326121591-642.coolpic» v:shapes="_x0000_i1423">
· На 12 день:<img width=«414» height=«24» src=«ref-2_326122233-626.coolpic» v:shapes="_x0000_i1424">
· На 13 день: <img width=«420» height=«24» src=«ref-2_326122859-635.coolpic» v:shapes="_x0000_i1425">
Занесем полученные результаты в таблицу 26.
Параметр сглаживания, а=0,45
Рассчитаем прогноз курсов акций с параметрами сглаживания, а=0,45 по формуле (57):
· На 7 день:<img width=«383» height=«25» src=«ref-2_326123494-592.coolpic» v:shapes="_x0000_i1426">
· На 8 день:<img width=«405» height=«25» src=«ref-2_326124086-606.coolpic» v:shapes="_x0000_i1427">
· На 9 день:<img width=«416» height=«25» src=«ref-2_326124692-624.coolpic» v:shapes="_x0000_i1428">
· На 10 день:<img width=«417» height=«25» src=«ref-2_326125316-629.coolpic» v:shapes="_x0000_i1429">
· На 11 день: <img width=«438» height=«25» src=«ref-2_326125945-671.coolpic» v:shapes="_x0000_i1430">
· На 12 день: <img width=«437» height=«24» src=«ref-2_326126616-659.coolpic» v:shapes="_x0000_i1431">
· На 13 день:<img width=«439» height=«24» src=«ref-2_326127275-670.coolpic» v:shapes="_x0000_i1432">
Занесем полученные результаты в таблицу 27.
--PAGE_BREAK--
еще рефераты
Еще работы по экономике
Реферат по экономике
Технико-экономический анализ деятельности предприятия 3
2 Сентября 2013
Реферат по экономике
Системный подход к исследованию конъюнктуры рынка
2 Сентября 2013
Реферат по экономике
Лизинговая деятельность в Украине развитие, особенности, законодательная база
2 Сентября 2013
Реферат по экономике
Лизинг как форма аренды и вид инвестиционной деятельности
2 Сентября 2013