Автор работы: Пользователь скрыл имя, 28 Марта 2011 в 17:24, контрольная работа
Изучение рынка сбыта показало, что суточный спрос на краску I никогда не превышает спроса на краску T более чем на 1 т. Кроме того, установлено, что спрос на краску I никогда не превышает 2 т в сутки. Оптовые цены одной тонны красок равны 3000 ден. ед. для краски T и 2000 ден. ед. для краски I. Какое количество краски каждого вида должны производить фабрика, чтобы доход от реализации продукции был максимальным?
Построить экономико-математическую модель задачи, дать необходимые комментарии к ее элементам и получить решение графическим методом. Что произойдет, если решать задачу на минимум, и почему?
Начальные оценки параметров получим по первым пяти точка при помощи метода наименьших квадратов.
na0 + a1∑x = ∑y
a0∑x
+ a1∑x2 = ∑xy
5a0 + 15a1 = 49
15a0
+ 55a1 = 172
∑y ∙ ∑x2 – ∑xy ∙ ∑x
a0 = ————————
n ∑x2 – ∑x ∙ ∑x
49 ∙ 55 – 172 ∙ 15
a0 = ——————— ≈ 2,30
5 ∙ 55 – 15 ∙ 15
n∑xy – ∑y ∙ ∑x
a1 = ———————
n ∑x2 – ∑x ∙ ∑x
5 ∙ 172 – 49 ∙ 15
a1 = ——————– ≈ 2,50
5 ∙ 55 – 15 ∙ 15
Данные
для расчета возьмем в
Сумма | ||||||
x | 1 | 2 | 3 | 4 | 5 | 15 |
y | 5 | 7 | 10 | 12 | 15 | 49 |
x2 | 1 | 4 | 9 | 16 | 25 | 55 |
xy | 5 | 14 | 30 | 48 | 75 | 172 |
Уравнение линейной регрессии имеет вид: yx = 2,30 + 2,50x.
Получили a0(0) = 2,30, a1(0) = 2,50.
Возьмем α = 0,4, k = 1 и β = 1 – α = 1 – 0,4 = 0,6.
Будем находить последующие значения a0(t) и a1(t) по формулам
a1(t) = a1(t – 1) + (1 – β)2 ∙ (Y(t) – Yp(t)) и a0(t) = a0(t – 1) + a1(t – 1) + (Y(t) – Yp(t)) ∙ (1 – β2),
где Yp(t)
= a0(t – 1) + a1(t – 1)k.
Таблица 3.
Номер | Факт | a0 | a1 | Расчет | Отклонение | ε2 |
2,30 | 2,50 | |||||
1 | 5 | 4,93 | 2,53 | 4,80 | 0,200 | 0,04 |
2 | 7 | 7,17 | 2,46 | 7,46 | -0,460 | 0,21 |
3 | 10 | 9,86 | 2,52 | 9,62 | 0,376 | 0,14 |
4 | 12 | 12,14 | 2,46 | 12,38 | -0,383 | 0,15 |
5 | 15 | 14,85 | 2,52 | 14,60 | 0,405 | 0,16 |
6 | 18 | 17,78 | 2,62 | 17,38 | 0,624 | 0,39 |
7 | 20 | 20,14 | 2,56 | 20,40 | -0,397 | 0,16 |
8 | 23 | 22,89 | 2,61 | 22,70 | 0,299 | 0,09 |
9 | 26 | 25,82 | 2,69 | 25,50 | 0,502 | 0,25 |
10 | 28,51 | 1,59 | ||||
11 | 31,19 |
___________________
u(1) = 1,12 ∙ 0,48 ∙ √1 + 1/9 + (9 + 1 – 5)2/60 ≈ 3,25
k = 1 (t = 10).
Нижняя граница: 28,51 – 0,66 = 27,85
Верхняя граница: 28,51 0,66 = 29,17
___________________
u(2) = 1,12 ∙ 0,48 ∙ √1 + 1/9 + (9 + 2 – 5)2/60 ≈ 3,44
k = 2 (t = 10).
Нижняя граница: 31,19 – 0,70 = 30,49
Верхняя
граница: 31,19 + 0,70 = 31,89
Проверим адекватность модели.
Проверка случайности ряда остатков по критерию пиков дает результат 7 больше 2 (критическое число поворотных точек).
Вычислим
d по формуле
∑(εt – εt-1)2 3,95
d = ————— = —— = 2,48
∑ εt2
1,59
При проверке независимости уровней ряда остатков друг от друга значение d′ = 4 – 2,48 = 1,52 при уровне значимости a = 0,025 больше d2 = 1,36, т.е. модель адекватна.
Соответствие
ряда остатков нормальному распределению
установим с помощью формулы
R/S = (εmax – εmin)/Sn,
_______________ _____
Sn
= √ ∑(εt – εср)2/(n – 1)
= √1,44/8 = 0,42
R/S = (0,62 – (-0,46))/0,42 = 2,55
Для n = 9 и a = 0,05 найдем критический интервал: [2,7; 3,7]. Вычисленное значение 2,55 не попадает между табулированными границами с заданным уровнем вероятности. Значит, закон нормального распределения не выполняется.
Результаты
аппроксимации и
Возьмем
α = 0,7, k = 1 и β = 1 – α = 1 – 0,7 = 0,3.
Будем
находить последующие значения a0(t)
и a1(t) по формулам
a1(t) = a1(t – 1) + (1 – β)2 ∙ (Y(t) – Yp(t)) и a0(t) = a0(t – 1) + a1(t – 1) + (Y(t) – Yp(t)) ∙ (1 – β2),
где Yp(t)
= a0(t – 1) + a1(t – 1)k.
Таблица 4.
Номер | Факт | a0 | a1 | Расчет | Отклонение | ε2 |
2,30 | 2,50 | |||||
1 | 5 | 4,98 | 2,60 | 4,80 | 0,200 | 0,04 |
2 | 7 | 7,05 | 2,31 | 7,58 | -0,580 | 0,34 |
3 | 10 | 9,94 | 2,62 | 9,37 | 0,634 | 0,40 |
4 | 12 | 12,05 | 2,35 | 12,57 | -0,567 | 0,32 |
5 | 15 | 14,95 | 2,64 | 14,40 | 0,602 | 0,36 |
6 | 18 | 17,96 | 2,84 | 17,59 | 0,413 | 0,17 |
7 | 20 | 20,07 | 2,45 | 20,81 | -0,807 | 0,65 |
8 | 23 | 22,96 | 2,68 | 22,52 | 0,479 | 0,23 |
9 | 26 | 25,97 | 2,86 | 25,64 | 0,360 | 0,13 |
10 | 28,83 | 2,64 | ||||
11 | 31,69 |
___________________
u(1) = 1,12 ∙ 0,61 ∙ √1 + 1/9 + (9 + 1 – 5)2/60 ≈ 0,85
k = 1 (t = 10).
Нижняя
граница: 28,83 – 0,85 = 27,98
Верхняя граница: 28,83 + 0,85 = 29,68
___________________
u(2) = 1,12 ∙ 0,61 ∙ √1 + 1/9 + (9 + 2 – 5)2/60 ≈ 0,90
k = 2 (t = 10).
Нижняя граница: 31,69 – 0,90 = 30,79
Верхняя
граница: 31,69 + 0,90 = 32,59
Проверим адекватность модели.
Проверка
случайности ряда остатков по критерию
пиков дает результат 7 больше 2 (критическое
число поворотных точек).
Вычислим
d по формуле
∑(εt – εt-1)2 8,08
d = ————— = —— = 3,06
∑ εt2
2,64
При проверке независимости уровней ряда остатков друг от друга значение d′ = 4 – 3,06 = 0,94 при уровне значимости a = 0,025 меньше d1 = 1,08, т.е. модель неадекватна.
Соответствие
ряда остатков нормальному распределению
установим с помощью формулы
R/S = (εmax – εmin)/Sn,
_______________ _____
Sn
= √ ∑(εt – εср)2/(n – 1)
= √2,58/8 = 0,57
R/S = (0,57 – (-0,81))/0,57 = 2,54
Для
n = 9 и a = 0,05 найдем критический интервал:
[2,7; 3,7]. Вычисленное значение 2,54 не попадает
между табулированными границами с заданным
уровнем вероятности. Значит, закон нормального
распределения не выполняется.
Результаты
аппроксимации и
Очевидно,
что лучше взять α = 0,4.
4. Для
того чтобы оценить параметры
и качество этой модели (адекватность
и точность), а также построить
точечный и интервальный прогнозы, заполним
следующую таблицу: