Построение экономико-математических моделей производства продукции
Нахождение оптимального значения целевой функции, позволяющей минимизировать себестоимость произведенной продукции. Оптимизационные задачи на максимум выручки от реализации готовой продукции. Экономико-математическая модель технологической матрицы.
Рубрика | Экономико-математическое моделирование |
Вид | контрольная работа |
Язык | русский |
Дата добавления | 25.10.2013 |
Размер файла | 248,8 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru
Размещено на http://www.allbest.ru
Министерство образования и науки Российской Федерации
Федеральное агентство по образованию ГОУ ВПО
Всероссийский заочный финансово - экономический институт
Контрольная работа
по экономико-математическим методам и прикладным моделям
Преподаватель: Князева И.В.
Калуга - 2008
Задача 1
Завод - производитель высококачественных элементов для автомобилей выпускает два различных типа деталей - X и Y. Завод располагает фондом рабочего времени в 4000 чел.-ч. в неделю. Для производства одной детали типа X требуется 1 чел.-ч., а для производства одной детали типа Y - 2 чел.-ч. Производственные мощности завода позволяют выпускать максимум 2250 деталей типа X и 1750 деталей типа Y в неделю. Каждая деталь типа X требует 2 кг металлических стержней и 5 кг листового металла, а для производства одной детали типа Y необходимо 5 кг металлических стержней и 2 кг листового материала. Уровень запасов каждого вида металла составляет 10 000 кг в неделю. Кроме того, еженедельно завод поставляет 600 деталей типа X своему постоянному заказчику. Существует также профсоюзное соглашение, в соответствии с которым общее число производимых в течение одной недели деталей должно составлять не менее 1500 штук.
Сколько деталей каждого типа следует производить, чтобы максимизировать общий доход за неделю, если доход от производства одной детали типа X составляет 30 ден. ед., а от производства одной детали типа Y - 40 ден. ед.?
Построить экономико-математическую модель задачи, дать необходимые комментарии к ее элементам и получить решение графическим методом. Что произойдет, если решать задачу на минимум, и почему?
Условие задачи кратко можно представить в таблице.
Показатель |
Типы деталей |
Необходимое кол-во сырья |
||
X |
Y |
|||
Время производства одной детали, чел./ч. |
1 |
2 |
4000 |
|
Необходимое количество металлических стержней для производства 1 детали, кг |
2 |
5 |
10000 |
|
Необходимое количество листового металла для производства 1 детали, кг |
5 |
2 |
10000 |
|
Доход с продажи 1 детали, ден. ед. |
30 |
40 |
Составим экономико-математическую модель.
Обозначим через X, Y объем производства соответствующего типа продукции.
Целевая функция - это математическая запись критерия оптимальности, т.е. выражение, которое необходимо максимизировать f = =30X +40Y.
Ограничения по ресурсам:
X +2Y ? 4000,
2X +5Y ? 10000,
5X +2Y ? 10000,
X ? 2250,
Y ? 1750,
5X +2Y ? 10000,
X ? 600,
X +Y ? 1500,
X, Y ?0.
Приравняем левые и правые части неравенств и построим графики функций:
X +2Y = 4000,
2X +5Y = 10000,
5X +2Y = 10000,
X = 2250,
Y = 1750,
5X +2Y = 10000,
X = 600,
X +Y = 1500,
X, Y = 0.
Учитывая знаки неравенств, получили многоугольник решений - АBCDE.
Найдем значения целевой функции в каждой вершине многоугольника:
1) А (600; 900), f (A) = 30*600 + 40*900 = 54000 ден. ед.
2) В (600; 1700), f (B) = 30*600 + 40*1700 = 86000 ден. ед.
3) С (1500; 1250), f (C) = 30*1500 + 40*1250 = 95000 ден. ед. - max
4) D (2000; 0), f (D) = 30*2000 + 40*0 = 60000 ден. ед.
5) Е (1500; 0), f (E) = 30*1500 + 40*0 = 45000 ден. ед. - min
Таким образом получаем результат: максимальное значение целевой функции равно 95000 ден. ед., что соответствует точке С многоугольника решений. Такой доход будет достигнут при выпуске продукции типа X в количестве 1500 штук и продукции Y в количестве 1250 штук.
Решая данную задачу на минимум, получим противоположное, т.е. минимальное значение целевой функции (45000 ден. ед. в точке Е, при выпуске продукции типа X в количестве 1500 штук и отсутствии выпуска продукции Y).
Решая задачу на минимум мы находим оптимальное значение целевой функции, позволяющей минимизировать себестоимость произведенной продукции. Решая обратную задачу (на максимум) мы находим оптимальное значение целевой функции, позволяющее максимизировать доход от продажи выпущенной продукции.
Задача 2
Предприятие выпускает четыре вида продукции и использует три вида оборудования: токарное, фрезерное, шлифовальное. Общий фонд рабочего времени оборудования каждого вида, нормы расхода и цены реализации единицы каждого вида продукции приведены в таблице.
Тип оборудования |
Нормы расхода ресурса на одно изделие |
Фонд рабочего времени, ч |
||||
А |
Б |
В |
Г |
|||
Токарное |
2 |
1 |
1 |
3 |
300 |
|
Фрезерное |
1 |
0 |
2 |
1 |
70 |
|
Шлифовальное |
1 |
2 |
1 |
0 |
340 |
|
Цена изделия |
8 |
3 |
2 |
1 |
Требуется:
1. Сформулировать прямую оптимизационную задачу на максимум выручки от реализации готовой продукции, получить оптимальный план выпуска продукции.
2. Сформулировать двойственную задачу и найти ее оптимальный план с помощью теорем двойственности.
3. Пояснить нулевые значения переменных в оптимальном плане.
4. На основе свойств двойственных оценок и теорем двойственности:
· проанализировать использование ресурсов в оптимальном плане исходной задачи;
· определить, как изменятся выручка и план выпуска продукции, если фонд рабочего времени шлифовального оборудования увеличить на 24 часа;
· оценить целесообразность включения в план изделия Д ценой 11 единиц, если нормы затрат оборудования 8, 2 и 2 единицы соответственно.
Решение:
1. Сформулируем прямую оптимизационную задачу на максимум выручки от реализации готовой продукции.
Целевая функция будет иметь вид max f(x) = 8X1 + 3X2 + 2X3 + X4.
Ограничения целевой функции будут иметь вид:
2X1 + X2 + X3 + 3X4 ? 300,
X1 + 2X3 + 1X4 ? 70,
X1 + 2X2 + X3 ? 340,
X1, X2, X3, X4 ?0.
Найдем оптимальный план выпуска продукции с использованием программного обеспечение Microsoft Excel. Получаем оптимальный план: X = (70; 135; 0; 0).
2. Сформулируем двойственную задачу:
· если в исходной задаче неравенства имеют вид « ? », то в двойственной задаче неравенства будут иметь вид « ? »;
· число неравенства в двойственной задаче равно числу переменных в целевой функции исходной задачи;
· если исходная задача на «max» двойственная будет на «min»;
· коэффициентами целевой функции двойственной задачи являются правые части неравенств исходной задачи;
· матрица коэффициентов в ограничениях двойственной задачи транспонирована к матрице коэффициентов исходной задачи;
· свободные члены ограничений в двойственной задаче равны коэффициентам целевой функции исходной задачи.
Таким образом, двойственная задача будет иметь вид:
min g(y) = 300Y1 + 70Y2 + 340Y3
Ограничения двойственной функции будут иметь вид:
2Y1 + Y2 + Y3 ? 8,
Y1 + 2Y3 ? 3,
Y1 + 2Y2 + Y3 ? 2,
3Y1 + Y2 ? 1,
Y1, Y2, Y3?0.
Найдем оптимальный план двойственной задачи. Используя теоремы двойственности.
Подставим значения оптимального плана исходной целевой функции в исходные ограничения, получаем: X = (70; 135; 0; 0).
· 2*70 + 1*135 + 1*0 + 3*0 = 275 < 300 (удовлетворяет условию неравенства);
· 1*70 + 0*135 + 2*0 + 1*0 = 70 =70 (удовлетворяет условию равенства);
· 1*70 + 2*135 + 1*0 + 0*0 = 340 = 340 (удовлетворяет условию равенства);
· max f(x) = 8*70 +3*135 +2*0 + 1*0 = 965 ден. ед.
Вышеизложенные расчеты означают, что при оптимальном плане II и III виды оборудования (фрезерное и шлифовальное оборудование) используются полностью, т.е. по максимуму их производственных возможностей. I тип оборудования (токарное) используется не полностью, т.е. существует некоторое время простоя этого оборудования.
Из второй теоремы двойственности следует, что Y1 = 0, т.к. первое ограничение в исходной задаче - строгое неравенство.
Второе и третье ограничения в исходной задаче являются строгими равенствами, таким образом приведем неравенства двойственной задачи к равенствам и решим систему уравнений.
2Y1 + Y2 + Y3 = 8,Y1 + 2Y3 = 3; |
2*0 + Y2 + Y3 = 8,0 + 2Y3 = 3; |
||
Y2 = 13/2,Y3 = 3/2. |
Оптимальный план двойственной задачи равен: Y = (0; 13/2; 3/2).
Т.е. min g (y) = 300*0 + 70*13/2 + 340*3/2 = 965 ден. ед.
3. Поясним нулевые значения переменных в оптимальном плане.
В исходной задаче X3 и X4 равны нулю. Это означает, что для максимизации выручки от продажи продукции целесообразно не выпускать продукцию вида В и Г. Только в этом случаю будет достигнут максимум целевой исходной функции. Из исходных данных видно, что для производства этих двух видов продукции необходимо большой расход мощностей оборудования, а цена очень низкая. Именно поэтому не следует включать в план производства продукты вида В и Г.
4. Поясним нулевые значения двойственной задачи.
В двойственной задаче Y1 = 0. Это означает, что в нашем случае первый вид оборудования (токарное) недефицитно, т.е. предприятие обладает достаточными мощностями по токарных станков для производства продукции.
Второй вид оборудования (фрезерное) более дефицитно, чем третий вид (шлифовальное оборудование), т.к. 13/2 > 3/2.
5. Определим, как изменятся выручка и план выпуска продукции, если фонд рабочего времени шлифовального оборудования увеличить на 24 часа.
Используем теорему об оценках.
,
y3 = 3/2, Дb3 = 24, b3 = 364,
ден. ед.
max f(x) = 965 + 36 = 1001 ден. ед.
Таким образом, общая выручка продукции увеличится на 36 ден. ед.
Найдем новый оптимальный план:
Подставим значения оптимального плана исходной целевой функции в исходные ограничения, получаем: X = (70; 135; 0; 0).
2*X1 + 1*X2 + 1*0 + 3*0 = 300;
1*X1 + 0*X2 + 2*0 + 1*0 = 70;
1*X1 + 2*X2 + 1*0 + 0*0 = 364.
X1 = 70,
X2 = 147.
Новый оптимальный план составит: X = (70; 147; 0; 0).
6. Оценим целесообразность включения в план изделия Д ценой 11 единиц, если нормы затрат оборудования 8, 2 и 2 единицы соответственно.
Исходные данные запишем в таблице.
Тип оборудования |
Объективно обусловленные оценки |
Нормы расхода ресурса на одно изделие |
|
Токарное |
0 |
8 |
|
Фрезерное |
13/2 |
2 |
|
Шлифовальное |
3/2 |
2 |
|
Цена изделия |
11 |
,
ден. ед.
Продукцию «Д» целесообразно включить в общий план выпуска продукции, т.к. выручка от реализации превысит общую себестоимость и тем самым будет получена прибыль в размере 5 ден. ед. за 1 единицу продукции.
Задача 3
Промышленная группа предприятий (холдинг) выпускает продукцию трех видов, при этом каждое из трех предприятий специализируется на выпуске продукции одного вида: первое предприятие специализируется на выпуске продукции первого вида, второе предприятие - продукции второго вида, третье предприятие - продукции третьего вида. Часть выпускаемой продукции потребляется предприятиями холдинга (идет на внутреннее потребление), остальная часть поставляется за его пределы (внешним потребителям, является конечным продуктом). Специалистами управляющей компании получены экономические оценки (I = 1, 2, 3; j = 1, 2, 3) элементов технологической матрицы А (норм расхода, коэффициентов прямых материальных затрат) и элементов вектора конечной продукции Y.
Требуется:
1. Проверить продуктивность технологической матрицы А = () (матрицы коэффициентов прямых материальных затрат).
2. Построить баланс (заполнить таблицу) производства и распределения продукции предприятия холдинга.
экономическая математическая модель матрица
Предприятия (виды продукции) |
Коэффициенты прямых затрат |
Конечный продукт Y |
|||
1 |
2 |
3 |
|||
1 |
0,1 |
0,2 |
0,4 |
100 |
|
2 |
0,0 |
0,4 |
0,1 |
200 |
|
3 |
0,1 |
0,3 |
0,4 |
100 |
Решение
1. Проверим продуктивность матрицы А (матрицы коэффициентов прямых материальных затрат).
Запишем матрицу А:
А = |
0,1 |
0,2 |
0,4 |
|
0,0 |
0,4 |
0,1 |
||
0,1 |
0,3 |
0,4 |
Найдем значение матрицы (Е-А):
(Е-А) = |
1 |
0 |
0 |
- |
0,1 |
0,2 |
0,4 |
= |
0,9 |
-0,2 |
-0,4 |
|
0 |
1 |
0 |
0,0 |
0,4 |
0,1 |
0,0 |
0,6 |
-0,1 |
||||
0 |
0 |
1 |
0,1 |
0,3 |
0,4 |
-0,1 |
-0,3 |
0,6 |
Найдем значение матрицы |E-A|:
|Е-А| = |
0,9 |
-0,2 |
-0,4 |
||
0,0 |
0,6 |
-0,1 |
= 0,9*0,6*0,6 + (-0,1)*(-0,2)*(-0,1) + (-0,4)*0,0*(-0,3) - |
||
-0,1 |
-0,3 |
0,6 |
|||
- (-0,1)*0,6*(-0,4) - 0,6*0,0*(-0,2) - 0,9*(-0,3)*(-0,1) = 0,271 |
Найдем значение матрицы (Е-А)':
(Е-А)' = |
0,9 |
0,0 |
-0,1 |
|
-0,2 |
0,6 |
-0,3 |
||
-0,4 |
-0,1 |
0,6 |
Определим значения каждого элемента матрицы :
А11 = (-1)2 = |
0,6 |
-0,3 |
= 0,33; |
|
-0,1 |
0,6 |
|||
А12 = (-1)3 = |
-0,2 |
-0,3 |
= 0,24; |
|
-0,4 |
0,6 |
|||
А13 = (-1)4 = |
-0,2 |
0,6 |
= 0,26; |
|
-0,4 |
-0,1 |
|||
А21 = (-1)3 = |
0,0 |
-0,1 |
= 0,01; |
|
-0,1 |
0,6 |
|||
А22 = (-1)4 = |
0,9 |
-0,1 |
= 0,5; |
|
-0,4 |
0,6 |
|||
А23 = (-1)5 = |
0,9 |
0,0 |
= 0,09; |
|
-0,4 |
-0,1 |
|||
А31 = (-1)4 = |
0,0 |
-0,1 |
= 0,06; |
|
0,6 |
-0,3 |
|||
А32 = (-1)5 = |
0,9 |
-0,1 |
= 0,29; |
|
-0,2 |
-0,3 |
|||
А33 = (-1)6 = |
0,9 |
0,0 |
= 0,54. |
|
-0,2 |
0,6 |
Таким образом получили матрицу :
= |
0,33 |
0,24 |
0,26 |
|
0,01 |
0,5 |
0,09 |
||
0,06 |
0,29 |
0,54 |
Найдем матрицу В:
0,33 |
0,24 |
0,26 |
= |
1,218 |
0,886 |
0,959 |
||
0,01 |
0,5 |
0,09 |
0,037 |
1,845 |
0,332 |
|||
0,06 |
0,29 |
0,54 |
0,221 |
1,070 |
1,993 |
На основании вышеизложенного можно сделать вывод, о том, что матрица А - продуктивна, т.к. матрица (Е-А)-1 существует о больше нуля.
2. Построим баланс производства и распределения продукции предприятия холдинга. Найдем значение матрицы X:
X = В * Y = |
1,22 |
0,89 |
0,96 |
= |
100 |
= |
396 |
|
0,04 |
1,85 |
0,33 |
200 |
407 |
||||
0,22 |
1,07 |
1,99 |
100 |
435 |
Найдем элементы первого квадранта по формуле:
Xij = |
0,1*396 |
0,2*407 |
0,4*435 |
= |
40 |
81 |
174 |
|
0,0*396 |
0,4*407 |
0,1*435 |
0 |
163 |
44 |
|||
0,1*396 |
0,3*407 |
0,4*435 |
40 |
122 |
174 |
Условно чистую продукцию найдем по формуле:
Z1 = 396 - (40 + 0 + 40) = 316,
Z2 = 407 - (81 + 163 = 122) = 41,
Z3 = 435 - (174 + 44 + 174) = 43.
Таким образом получаем баланс производства и распределения продукции предприятия холдинга.
Производ-ственные отрасли |
Потребляющие отрасли |
Конечная продукция |
Валовая продукция |
|||||
1 |
2 |
3 |
||||||
1 |
I |
40 |
81 |
174 |
II |
100 |
396 |
|
2 |
0 |
163 |
44 |
200 |
407 |
|||
3 |
40 |
122 |
174 |
100 |
435 |
|||
Условно чистая продукция |
III |
316 |
41 |
43 |
IV |
400 |
||
Валовая продукция |
396 |
407 |
435 |
1238 |
Задача 4
В течение девяти последовательных недель фиксировался спрос Y (t) (млн. руб.) на кредитные ресурсы финансовой компании. Временной ряд Y (t) этого показателя (повариантного) приведен ниже в таблице.
T |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
|
Y |
20 |
27 |
30 |
41 |
45 |
51 |
51 |
55 |
61 |
Требуется:
1. Проверить наличие аномальных наблюдений.
2. Построить линейную модель , параметры которой оценить МНК ( - расчетные, смоделированные значения временного ряда).
3. Построить адаптивную модель Брауна
4. Оценить адекватность построенной модели, используя свойства независимости остаточной компоненты, случайности и соответствия нормальному закону распределения (при использовании R/S - критерия взять табулированные границы 2,7 - 3,7).
5. Оценить точность модели на основе использования средней ошибки аппроксимации.
6. По построенной модели осуществить прогноз спроса на следующие две недели (доверительный интервал прогноза рассчитать при доверительной вероятности р = 70%).
Фактические значения показателя, результаты моделирования и прогнозирования представить графически.
Решение:
1. Проверим наличие аномальных наблюдений.
Для диагностики аномальных наблюдений разработаны различные критерии. Например, для всех или только подозреваемых в аномальности наблюдений вычисляем среднее значение и среднеквадратическое отклонение двух (или четырех) соседних с ними значений:
Ycp(t) = (Y(t-1) + Y(t+1))/2, t=2,3...N-1;
Sy(t) =
Далее вычисляем величину h(t):
h(t) = (Y(t) - Y(t-1)) / Sy(t)
t |
Yt |
Yср |
Yt - Yср |
(Yt - Yср)2 |
Yt+1 - Yср |
Yt-1 - Yср |
St |
Ht |
|
1 |
20 |
||||||||
2 |
27 |
25 |
-5 |
25 |
5 |
25 |
5 |
1.4 |
|
3 |
30 |
34 |
-7 |
49 |
7 |
49 |
7 |
0.43 |
|
4 |
41 |
37.5 |
-7.5 |
56.25 |
7.5 |
56.25 |
7.5 |
1.47 |
|
5 |
45 |
46 |
-5 |
25 |
5 |
25 |
5 |
0.8 |
|
6 |
51 |
48 |
-3 |
9 |
3 |
9 |
3 |
2 |
|
7 |
51 |
53 |
-2 |
4 |
2 |
4 |
2 |
0 |
|
8 |
55 |
56 |
-5 |
25 |
5 |
25 |
5 |
0.8 |
|
9 |
61 |
Если рассчитанная величина h(t) превышает табличный уровень (для трех наблюдений он равен 2,1), уровень Y(t) считается аномальным. Как видно из таблицы все значение h(t) ниже табличного, поэтому можно сделать вывод об отсутствии аномальных наблюдений.
2. Построим линейную модель
где а0 и a1 - параметры модели, а t =1, 2,..., n.
В результате получим так называемую систему нормальных уравнений:
Решая эту систему двух линейных уравнений с двумя неизвестными, получим:
где , и -- среднее значения, соответственно, моментов наблюдения и уровней ряда.
t |
Yt |
||||||
1 |
20 |
-4 |
16 |
-22,3 |
89,3 |
22,3 |
|
2 |
27 |
-3 |
9 |
-15,3 |
46,0 |
27,3 |
|
3 |
30 |
-2 |
4 |
-12,3 |
24,7 |
32,3 |
|
4 |
41 |
-1 |
1 |
-1,3 |
1,3 |
37,3 |
|
5 |
45 |
0 |
0 |
2,7 |
0,0 |
42,3 |
|
6 |
51 |
1 |
1 |
8,7 |
8,7 |
47,3 |
|
7 |
51 |
2 |
4 |
8,7 |
17,3 |
52,3 |
|
8 |
55 |
3 |
9 |
12,7 |
38,0 |
57,3 |
|
9 |
61 |
4 |
16 |
18,7 |
74,7 |
62,3 |
|
У 45 |
У 381 |
У 0 |
У 60 |
У 0,0 |
У 300,0 |
У 381,0 |
Параметры линейной модели будут соответственно равны:
Таким образом получена линейная модель . Параметр а1=5 показывает, что при возрастании t на единицу, показатель Yt = в среднем увеличивается на 2,6.
3. Построим модель Брауна.
Расчетное значение в момент времени получается по формуле
Yp(t) = a0(t - 1) + a1(t - 1)k (5.3.1)
где k - количество шагов прогнозирования (обычно k = 1).
Это значение сравнивается с фактическим уровнем, и полученная ошибка прогноза E(t) = Y(f) - Yp(t) используется для корректировки модели. Корректировка параметров осуществляется по формулам:
где в - коэффициент дисконтирования данных, отражающий большую степень доверия к более поздним данным. Его значение должно быть в интервале от 0 до 1 в зависимости от значения б (б = 0,4 или б = 0,7).
Рассмотрим модель Брауна со значением б = 0,4.
t |
Yt |
а0 |
а1 |
Yp(t) |
Е(t) |
| Е(t)| |
|
0 |
17.3 |
5 |
|||||
1 |
20 |
20.8 |
4.6 |
22.3 |
-2.3 |
2.3 |
|
2 |
27 |
26.5 |
4.9 |
25.5 |
1.5 |
1.5 |
|
3 |
30 |
27.3 |
3.9 |
31.4 |
-1.4 |
1.4 |
|
4 |
41 |
32.5 |
4.2 |
31.2 |
9.8 |
9.8 |
|
5 |
45 |
33.9 |
3.5 |
36.7 |
8.3 |
8.3 |
|
6 |
51 |
37.2 |
3.4 |
37.5 |
13.5 |
13.5 |
|
7 |
51 |
36.4 |
2.4 |
40.7 |
10.3 |
10.3 |
|
8 |
55 |
40.0 |
2.7 |
38.8 |
16.2 |
16.2 |
|
9 |
61 |
42.3 |
2.6 |
42.7 |
18.3 |
18.3 |
|
У 45 |
У 381 |
У 306.6 |
У 74.4 |
У 81.6 |
Рассмотрим модель Брауна со значением б = 0,7
t |
Yt |
a0 |
а1 |
Yp(t) |
Е(t) |
| Е(t)| |
|
0 |
17.3 |
5 |
|||||
1 |
20 |
20.2 |
3.9 |
22.3 |
-2.3 |
2.3 |
|
2 |
27 |
26.7 |
5.3 |
24.1 |
2.9 |
2.9 |
|
3 |
30 |
25.4 |
1.7 |
32.0 |
-2.0 |
-2.0 |
|
4 |
41 |
36.6 |
6.9 |
27.1 |
13.9 |
13.9 |
|
5 |
45 |
26.0 |
-2.5 |
43.5 |
1.5 |
1.5 |
|
6 |
51 |
57.6 |
15.8 |
23.5 |
27.5 |
27.5 |
|
7 |
51 |
-18.8 |
-33.9 |
73.4 |
-22.4 |
22.4 |
|
8 |
55 |
130.3 |
64.7 |
-52.7 |
107.7 |
107.7 |
|
9 |
61 |
-339.1 |
-222.9 |
195.0 |
-134.0 |
134.0 |
|
У 45 |
У 381 |
У 388.2 |
У 314.2 |
Как видно из вышеизложенных расчетов модель Брауна со значением б = 0,4 лучше приближает опытные данные, так как сумма отклонений расчетных данных от реальных величин гораздо меньше. Проиллюстрируем на графике вышеуказанные расчеты.
Как видно из графика, выводы сделанные по результатам построения двух моделей Брауна подтверждены графически.
4. Оценим адекватность построенной модели, используя свойства независимости остаточной компоненты, случайности и соответствия нормальному закону распределения.
Построим расчетную таблицу.
t |
Yt |
пов. точ. |
еt2 |
еt - еt-1 |
(еt - еt-1)2 |
||||
1 |
20 |
22,3 |
-2,3 |
- |
5,29 |
- |
0,115 |
||
2 |
27 |
27,3 |
-0,3 |
1 |
0,09 |
2 |
4 |
0,011 |
|
3 |
30 |
32,3 |
-2,3 |
1 |
5,29 |
-2 |
4 |
0,077 |
|
4 |
41 |
37,3 |
3,7 |
1 |
13,69 |
6 |
36 |
0,090 |
|
5 |
45 |
42,3 |
2,7 |
1 |
7,29 |
-1 |
1 |
0,060 |
|
6 |
51 |
47,3 |
3,7 |
1 |
13,69 |
1 |
1 |
0,073 |
|
7 |
51 |
52,3 |
-1,3 |
0 |
1,69 |
-5 |
25 |
0,025 |
|
8 |
55 |
57,3 |
-2,3 |
1 |
5,29 |
-1 |
1 |
0,042 |
|
9 |
61 |
62,3 |
-1,3 |
- |
1,69 |
1 |
1 |
0,021 |
|
У 45 |
У 381 |
У 380,7 |
У 0,0 |
У 6 |
У 54,0 |
У 73 |
У 0,514 |
Для проверки условия случайности возникновении отдельных отклонений от тренда часто используется критерий, основанный на поворотных точках. Значение случайной переменной считается поворотной точкой, если оно одновременно больше соседних с ним элементов или, наоборот, меньше значений предыдущего и последующего за ним члена, т.е.
Если еt-1 < еt > еt+1, если еt-1 >еt < еt+1, тогда р=1.
Если еt-1 < еt < еt+1, если еt-1 >еt > еt+1, тогда р=0.
Существует определенная зависимость между средней арифметической , дисперсией ; количества поворотных точек р и числом членов исходного ряда наблюдений n. В случайной выборке средняя арифметическая (математическое ожидание) числа поворотных точек равна , а их дисперсия вычисляется по формуле . Учитывая эти соотношения, критерий случайности отклонений от тренда при уровне вероятности 0,95 можно представить, как
где р -- фактическое количество поворотных точек в случайном ряду.
В нашем случае n=9, р=6. Подставляем данные значения в формулу и получаем 6 > 2 - неравенство выполняется, таким образом свойство случайности уровня остатков ряда выполняется.
Наличие (отсутствие) автокорреляции проще всего проверить с помощью критерия Дарвина--Уотсона. С этой целью строится статистика Дарвина-- Уотсона (d-статистика), в основе которой лежит расчетная формула:
.
Сравним полученное значение с верхним и нижним критическим значением (d1 = 1,08, d2 = 1,36). Таким образом получили d1 < d < d2 -- область неопределенности, когда нет оснований ни принять, ни отвергнуть гипотезу о существовании автокорреляции.
Соответствие рада остатков нормальному закону распределения проверим с помощью R/S-критерия:
R/S = (еmax - еmin)/S
где еmax и еmin - максимальный и минимальный уровни рада остатков соответственно;
S -- среднеквадратическое отклонение.
.
R/S = (3,7 - (-2,3)) / 2,6 = 2,3 < 2,7, т.к.
расчетное значение R/S критерия не попадает между табличными границами (2,7; 3,7), то уровни ряда остатков не подчиняются закону нормального распределения.
Среднюю относительную ошибку аппроксимации найдем по формуле:
Средняя относительная ошибка аппроксимации не превосходит 15% - точность модели считается приемлемой.
Из вышеизложенного можно сделать вывод о том, что так как не выполняются все условия адекватности, выбранная трендовая модель не является адекватной реальному раду экономической динамики.
Построим точечный прогноз по формуле:
Построим интервальный прогноз по формуле (, т.к. Р=70%):
Таким образом, с вероятностью 70% можно утверждать, что:
- прогнозируемая величина Y10 будет находится в интервале от 63,93 до 70,67;
- прогнозируемая величина Y11 будет находится в интервале от 68,73 до 75,87.
Изобразим графически выполненные расчеты.
Из графика видно, что построенная линейная модель достаточно хорошо приближает опытные данные.
Размещено на Allbest.ru
Подобные документы
Графический метод решения и построение экономико-математической модели производства. Определение выручки от реализации готовой продукции и расчет оптимального плана выпуска продукции. Баланс производства проверка продуктивность технологической матрицы.
задача [203,4 K], добавлен 03.05.2009Исследование методом Жордана-Гаусса системы линейных уравнений. Решение графическим и симплексным методом задач линейного программирования. Экономико-математическая модель задачи на максимум прибыли и нахождение оптимального плана выпуска продукции.
контрольная работа [177,8 K], добавлен 02.02.2010Экономико-математическая модель оптимального плана выпуска продукции. Оптимальная организация рекламной компании. Решение транспортной задачи: нахождение суммарных затрат на перевозку. Задача об оптимальном назначении (линейного программирования).
контрольная работа [812,0 K], добавлен 29.09.2010Расчет связи пунктов отправления и назначения. Обеспечение вывоза всех грузов из пункта отправления и ввоза в места назначения необходимых объемов. Экономико-математическая модель задачи на максимум прибыли, расчет оптимального плана выпуска продукции.
курсовая работа [49,1 K], добавлен 29.07.2011Объявление торгов администрацией штата на определенное количество строительных подрядов для определенного количества фирм. Экономико-математическая модели для минимизации затрат. Определение количества песцов и лисиц для получения максимальной прибыли.
контрольная работа [18,2 K], добавлен 05.03.2010Построение экономико-математической модели задачи, комментарии к ней и получение решения графическим методом. Использование аппарата теории двойственности для экономико-математического анализа оптимального плана задачи линейного программирования.
контрольная работа [2,2 M], добавлен 27.03.2008Задачи, функции и этапы построения экономико-математических моделей. Аналитические, анионные, численные и алгоритмические модели. Экономическая модель спортивных сооружений. Модели временных рядов: тенденции и сезонности. Теории массового обслуживания.
реферат [167,6 K], добавлен 22.07.2009Технико-экономические показатели производства продукции и потребления материальных ресурсов. Производительность и годовые фонды реакторов. Технологические способы изготовления эмалей. Составление экономико-математической модели задачи, анализ результатов.
контрольная работа [32,6 K], добавлен 06.01.2011Исследование методики построения модели и решения на ЭВМ с ее помощью оптимизационных экономико-математических задач. Характеристика программных средств, позволяющих решать такие задачи на ЭВМ. Определение оптимального варианта производства продукции.
лабораторная работа [79,3 K], добавлен 07.12.2013Особенности формирования и способы решения оптимизационной задачи. Сущность экономико-математической модели транспортной задачи. Характеристика и методика расчета балансовых и игровых экономико-математических моделей. Свойства и признаки сетевых моделей.
практическая работа [322,7 K], добавлен 21.01.2010