Контрольная работа по «Экономико-математические методы и прикладные модели»

 

Поскольку все элементы матрицы В неотрицательны, то матрица А продуктивна.

 

2. Руководствуясь балансовым методом планирования и экономическим смыслом прямых материальных затрат будем иметь следующую модель межотраслевого баланса:

 

          1,0 х₁ - 0,4х₂ - 0,1х₃ = 160

          -0,4х₁ + 0,9х₂ - 0,0х₃ = 180

          -0,3х₁ + 0,0х₂ + 0,9х₃ = 150

Для решения воспользуемся  пакетом Exсel.

В результате решения будем иметь  следующие объемы валового продукта по предприятиям: Х₁ =325,35;  Х₂ =344,60; Х₃ =275,12.

 

Схема межотраслевого баланса будет  выглядеть следующим образом:

Межотраслевой баланс производства и  распределения продукции
Производящие	Потребляющие структуры			Конечный	Валовой
структуры	1	2	3	продукт	продукт
1	0	137,84	27,51	160	325,35
2	130,14	34,46	0	180	344,60
3	97,61	0	27,51	150	275,12
Условно чистая	98,61	172,30	220,09	490	-
продукция
Валовой продукт	326,35	344,60	275,12	-	946,07

 

 

 

Задача 4

 

В течение  девяти последовательных недель фиксировался спрос Y(t) (млн. р.) на кредитные ресурсы финансовой компании. Временной ряд Y(t) этого показателя (повариантно) приведен ниже в таблице

 

Номер варианта	Номер наблюдения ( t = 1,2,…,9)
	1	2	3	4	5	6	7	8	9
8	8	13	15	19	25	27	33	35	40

 

 

Требуется: 

1) Проверить наличие  аномальных наблюдений. 

2) Построить линейную  модель   , параметры которой оценить МНК ( - расчетные, смоделированные значения временного ряда). 

3) Оценить адекватность построенных моделей, используя свойства независимости остаточной компоненты, случайности и соответствия нормальному закону распределения (при использовании R/S-критерия взять табулированные границы 2,7—3,7).

4) Оценить точность моделей на основе использования средней относительной ошибки аппроксимации.

5) По двум построенным моделям осуществить прогноз спроса на следующие две недели  (доверительный интервал прогноза рассчитать при доверительной вероятности р = 70%).

6) Фактические значения показателя, результаты моделирования и прогнозирования представить графически.

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

 

Решение:

 

1. Введем исходные данные.

Таблица 3.1

t	yt, факт.
1	8
2	13
3	15
4	19
5	25
6	27
7	33
8	35
9	40

 

Для проверки наличия аномальных наблюдений воспользуемся пакетом Excel. В результате решения будем иметь следующие данные:

 

Таблица 3.2

/yt.-yt-1/	с.к.о.	Характеристика Ирвина
5	10,9023	0,45861696
2		0,183446784
4		0,366893568
6		0,550340352
2		0,183446784
6		0,550340352
2		0,183446784
5		0,45861696

 

Чтобы найти  среднее квадратическое отклонение воспользуемся функцией СТАНДОТКЛОН Мастера функций Excel.

Табличное значение Величины Ирвина равно 1,5 , следовательно, в соответствии с методом Ирвина аномальные наблюдения не выявлены.

 

2. Построим линейную однопараметрическую модель регрессии Y от t. Для проведения регрессионного анализа воспользуемся надстройкой Excel Анализ данных. В результате получим следующее:

 

Таблица 3.3

Регрессионная статистика
Множественный R	0,996406256
R-квадрат	0,992825427
Нормированный R-квадрат	0,991800487
Стандартная ошибка	0,987219922
Наблюдения	9

 

Таблица 3.4

Дисперсионный анализ
	df	SS	MS	F	Значимость F
Регрессия	1	944,0666667	944,0666667	968,6677524	9,12922E-09
Остаток	7	6,822222222	0,974603175
Итого	8	950,8888889

 

Таблица 3.5

	Коэффициенты	Стандартная ошибка	t-статистика	P-Значение	Нижние 95%	Верхние 95%	Нижние 95,0%	Верхние 95,0%
Y-пересечение	4,06	0,717198646	5,654717251	0,000770691	2,359650245	5,751460866	2,359650245	5,751460866
t	3,97	0,127449544	31,12342771	9,12922E-09	3,665296384	4,268036949	3,665296384	4,268036949

 

 

Вывод остатка

Таблица 3.6

Наблюдение	Предсказанное yt, факт.	Остатки
1	8,02	-0,02
2	11,99	1,01
3	15,96	-0,96
4	19,92	-0,92
5	23,89	1,11
6	27,86	-0,86
7	31,82	1,18
8	35,79	-0,79
9	39,76	0,24
	Ср. знач.	0,00

 

Во втором столбце таблицы 3.5 содержатся коэффициенты уравнения регрессии а₀, а₁,. Уравнение регрессии зависимости Y_t (спрос на кредитные ресурсы) от t₁ (время) имеет вид:

 

Y(t) = 4,06 + 3,97t

 

3. Оценка адекватности модели.

1) С помощью  функции СРЗНАЧ Мастера функций Excel по таблице 3.6 найдем среднее значение остатков.

2) Построим  график остатков. По графику видно,  что P=7 больше 2, следовательно свойство случайности остатков выполняется.

 

Рис. 3.1 График остатков

 

3) С помощью функции КОРРЕЛ Мастера функций Excel по таблице 3.6 найдем коэффициент корреляции.

По след. формуле найдем t_расч.:

 

 

 

Коэффициент корреляции r = -0,62 незначим, поскольку

 

t_расч.=2,09 < t_табл.=2,36

 

Следовательно, свойство независимости остатков выполняется.

 

4) С помощью функции СТАНДОТКЛОН Мастера функций Excel по таблице 3.6 найдем среднее квадратическое отклонение.

 

S_ε = 0,92

 

Определим RS-критерий:

 

 

Расчетное значение попадает между табулированными границами (2,7-3,7) (для п=9 и 5-% уровня значимости), значит, остатки следуют нормальному закону распределения.

 

Модель в целом адекватна.

 

4. Оценим точность модели на основе использования средней относительной ошибки аппроксимации. Для этого рассчитаем в Excel следующую таблицу:

Таблица 3.7

	Точность
yt, факт.	abs остатков	Расчет Еотн
8	0,02	0,002778
13	1,01	0,077778
15	0,96	0,063704
19	0,92	0,048538
25	1,11	0,044444
27	0,86	0,031687
33	1,18	0,03569
35	0,79	0,02254
40	0,24	0,006111
	Еотн =	3,70%

 

Чтобы рассчитать вторую колонку воспользуемся функцией ABS Мастера функций Excel. Затем рассчитаем среднюю относительную ошибку аппроксимации:_.

 

Е_отн = 3,70 %

 

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

 

5. Прогноз спроса на кредитные ресурсы на следующие две недели.

 

1) Рассчитаем среднее значение фактора «время» (t_ср) и сумму квадратов отклонений t от его средней величины ( ∑(t-t_ср)² ).

 

Таблица 3.8

t	yt, факт.	t - tср
1	8	-4
2	13	-3
3	15	-2
4	19	-1
5	25	0
6	27	1
7	33	2
8	35	3
9	40	4
5		60
tср		∑(t-tср)2

 

2) Произведем точечный  и интервальный прогнозы на 2 шага  вперед.

Шаг прогноза к=1

t - статистика

1,12

Se =0,99

 

Y(t) = 4,06 + 3,97t

 

Y(10) = 4,06 + 3,97*10 = 43,76

Y(11) = 4,06 + 3,97*11 = 47,73

 

 

В результате расчетов получим  следующую таблицу:

Таблица 3.9

Время t	Шаг k	Прогноз	Нижняя граница	Верхняя граница
10	1	43,76	42,51	45,01
11	2	47,73	46,47	48,99

По полученным данным построим график подбора: