ВСЕРОССИЙСКИЙ ЗАОЧНЫЙ ФИНАНСОВО-ЭКОНОМИЧЕСКИЙ ИНСТИТУТ
Филиал г. Волгоград
КОНТРОЛЬНАЯ РАБОТА
по дисциплине «Экономико-математические методы и прикладные модели»
ВАРИАНТ № 3
Исполнитель:
Специальность «Финансы и кредит»
Группа ФНО
№ зачетной книжки
Руководитель:
Волгоград – 2009
Задача 1
Решить графическим методом типовую задачу оптимизации
Некоторая фирма выпускает два набора удобрений для газонов: обычный и улучшенный. В обычный набор входит 3 кг азотных, 4 кг фосфорных и 1 кг калийных удобрений, а в улучшенный – 2 кг азотных, 6 кг фосфорных и 3 кг калийных удобрений. Известно, что для некоторого газона требуется по меньшей мере 10 кг азотных, 20 кг фосфорных и 7 кг калийных удобрений. Обычный набор стоит 3 ден. Ед., а улучшенный – 4 ден. Ед. Какие и сколько наборов удобрений нужно купить, чтобы обеспечить эффективное питание почвы и минимизировать стоимость?
Построить экономико-математическую модель задачи, дать необходимые комментарии к ее элементам и получить решение графическим методом. Что произойдет, если решать задачу на максимум, и почему?
Решение:
Сформулируем прямую оптимизационную задачу.
Пусть х1 – количество обычных наборов удобрений;
х2 – количество улучшенных наборов удобрений.
Содержание в двух данных наборах азотных удобрений: 3х1 + 2х2
А для некоторого газона требуется по крайней мере 10 кг азотных удобрений, следовательно:
3х1 + 2х2 ? 10
Содержание в двух данных наборах фосфорных удобрений должно быть не менее 20 кг, т.е.:
4х1 + 6х2 ? 20
И содержание в двух данных наборах калийных удобрений должно быть не менее 7 кг, т.е.:
1х1 + 3х2 ? 7
Стоимость необходимых наборов удобрений составит:
3х1 + 4х2.
Таким образом, получим следующую экономико-математическую модель задачи:
min ?(х) = 3х1 + 4х2
3х1 + 2х2 ? 10
4х1 + 6х2 ? 20
х1 + 3х2 ? 7
х1 ? 0, х2 ? 0
Построим область решений системы ограничений. Для этого рассмотрим равенства и построим их графики – прямые.
1) 3х1 + 2х2 ? 10
3х1 + 2х2 = 10
Для нахождения полуплоскости, соответствующей данному неравенству, берем любую точку, не лежащую на граничной прямой, и подставляем ее координаты в неравенство.
Возьмем точку О(0;0):
3*0 + 2*0 ? 10
0 ? 10
Неравенство не выполняется, значит, исходному неравенству соответствует та полуплоскость, которая не содержит точку (0;0).
2) 4х1 + 6х2 ? 20
4х1 + 6х2 = 20
4*0 + 6*0 ³ 20
0 ³ 20
Неравенство не выполняется, значит, исходному неравенству соответствует полуплоскость, не содержащая точку О(0;0).
3) х1 + 3х2 ? 7
х1 + 3х2 = 7
1*0 + 3*0 ³ 7
0 ³ 7
Неравенство не выполняется, значит, исходному неравенству соответствует полуплоскость, не содержащая точку О(0;0).
4) х1 ? 0
х1 = 0 – ось ОХ2.
5) х2 ? 0
х2 = 0 – ось ОХ1.
Следовательно, область решений системы ограничений находится только в первой четверти декартовой системы координат.
Рис.1. Графическое решение ЗЛП
Находим общую часть всех построенных полуплоскостей. Это выпуклая заштрихованная область.
Для нахождения оптимального решения задачи изобразим графически функцию цели:
?(х) = d1x1 + d2x2
?(х) = 3х1 + 4х2
Для этого строим вектор d, начало которого в точке (0;0), а конец в точке (d1;d2).
d = (3; 4).
И строим одну из линий уровня функции цели (это линия, на которой функция цели принимает постоянное значение).
Для определения минимума данной функции, передвигаем линию уровня в направлении, противоположном вектору d, и видим, что она последний раз соприкасается с областью решений в точке В, где и будет достигнут min?(х).
Определим координаты точки В:
3х1 + 2х2 = 10 *(-3)
4х1 + 6х2 = 20
-9х1 – 6х2 = -30
4х1 + 6х2 = 20
Складываем почленно уравнения и получаем:
-5х1 = -10
х1 = 2
EMBED Equation.3
В(2; 2)
max ?(х) = 3*2 + 4*2 = 14 (ден. ед.)
Таким образом, чтобы минимизировать стоимость удобрений, нужно купить 2 обычных набора удобрений и 2 улучшенных набора удобрений. При этом минимальные затраты на покупку удобрений составят 14 денежных единиц.
Если решать данную задачу на максимум, то конечного оптимума не найдем, т.к. функция цели неограниченна, область решений системы ограничений бесконечна.
Задача 2
Использовать аппарат теории двойственности для экономико-математического анализа оптимального плана задачи линейного программирования.
Для изготовления четырех видов продукции используется три вида сырья. Запасы сырья, нормы его расхода и цены реализации единицы каждого вида продукции приведены в таблице.
Требуется:
1. Сформулировать прямую оптимизационную задачу на максимум выручки от реализации готовой продукции, получить оптимальный план выпуска продукции.
2. Сформулировать двойственную задачу и найти ее оптимальный план с помощью теорем двойственности.
3. Пояснить нулевые значения переменных в оптимальном плане.
4. На основе свойств двойственных оценок и теорем двойственности:
- проанализировать использование ресурсов в оптимальном плане исходной задачи;
- определить, как изменятся выручка и план выпуска продукции при увеличении запасов сырья I и II вида на 8 и 10 ед. соответственно и уменьшении на 5 ед. запасов сырья III вида;
- оценить целесообразность включения в план изделия Д ценой 10 единиц, на изготовление которого расходуется по две единицы каждого вида сырья.
Решение:
1. Сформулируем прямую оптимизационную задачу.
Пусть х1 – количество продукции вида А;
х2 – количество продукции вида Б;
х3 – количество продукции вида В;
х4 – количество продукции вида Г.
mах ?(х) = 5x1 + 7x2 + 3x3 + 6x4 – общая стоимость продукции;
2x1 + x2 + 3x3 + 2x4 ? 200
левая часть – расход сырья I типа на производство всей продукции
правая часть (200 ед.) – запас сырья I типа
x1 + 2x2 + 4x3 + 8x4 ? 160
левая часть – расход сырья II типа на производство всей продукции
правая часть (160 ед.) – запас сырья II типа
2x1 + 4x2 + x3 + x4 ? 170
левая часть – расход сырья III типа на пр-во всей продукции
правая часть (170 ед.) – запас сырья III типа
xj ? 0 (j = 1,2,3,4)
Оптимальный план задачи получим с помощью надстройки Excel Поиск решения (см. приложение 1).
Таким образом, получен оптимальный план:
x1 = 80
x2 = 0
x3 = 0
x4 = 10
mах ?(х) = 460 (ден. ед.)
Т.е. максимальный доход в 460 денежных единиц можно получить, если продукции вида А выпустить 80 единиц, вида Г – 10 ед.; вида Б и В не выпускать (т.к. х2 = 0 и х3 = 0); при этом ресурсы II и III типа используются полностью, а ресурс I не полностью, его остаток составляет 20 ед.
2. Сформулируем двойственную задачу.
уi – цена единицы ресурса i-го типа (i = 1,2,3)
min g(у) = 200у1 + 160у2 + 170у3 – общая стоимость ресурсов
2у1 + у2 + 2у3 ? 5
левая часть - стоимость ресурсов, расходуемых на производство
одного изделия вида А (если предприятие решило
продать имеющиеся ресурсы)
правая часть – сумма, которую предприятие может получить при
переработке сырья в готовую продукцию
у1 + 2у2 + 4у3 ? 7
стоимость ресурсов, идущих на пр-во одного изделия вида Б
3у1 + 4у2 + у3 ? 3
стоимость ресурсов, идущих на пр-во одного изделия вида В
2у1 + 8у2 + у3 ? 6
стоимость ресурсов, идущих на пр-во одного изделия вида Г
у1, у2, у3 ? 0
Найдем оптимальное решение двойственной задачи:
2у1 + у2 + 2у3 = 5
2у1 + 8у2 + у3 = 6
у1 = 0
у2 + 2у3 = 5
8у2 + у3 = 6 *(-2)
у2 + 2у3 = 5
-16у2 - 2у3 = -12
-15у2 = -7
у2 = EMBED Equation.3
у3 = EMBED Equation.3
у = (0; EMBED Equation.3 ; EMBED Equation.3 )
g(у) = 200*0 + 160* EMBED Equation.3 + 170* EMBED Equation.3 = 460
max g(у) = min ?(х) = 460
Следовательно, у = (0; EMBED Equation.3 ; EMBED Equation.3 ) – оптимальный план двойственной задачи.
3. Поясним нулевые значения переменных в оптимальном плане.
x2 = 0
x3 = 0
Т.е. продукцию вида Б и В выпускать не следует, т.к. их выпуск не приведет к максимальной выручке от реализации готовой продукции.
4. Проанализируем использование ресурсов в оптимальном плане.
Оценим использование ресурсов:
I типа: 2*80 + 1*0 + 3*0 + 2*10 = 180 < 200 – сырье I типа используется не полностью, остаток составляет 20 ед.; х5 = 20
II типа: 1*80 + 2*0 + 4*0 + 8*10 = 160 = 160 – сырье II типа используется полностью; х6 = 0
III типа: 2*80 + 4*0 + 1*0 + 1*10 = 170 = 170 – сырье III типа используется полностью; х7 = 0
х = (80, 0, 0, 10, 20, 0, 0)
Т.к. у1 = 0, то ресурс I типа используется не полностью.
Т.к. у2 ? 0 и у3 ? 0, то ресурсы II и III типов используются полностью. При чем ресурс III типа является более ценным, более дефицитным, чем ресурс II типа, т.к. у3 ? у2.
Определим, как изменятся выручка от реализации продукции и план ее выпуска при увеличении запасов сырья I и II видов на 8 и 10 единиц соответственно и уменьшении на 5 единиц запасов сырья III вида.
Изменение плана выпуска продукции ?Х находим следующим образом:
?Х = Д * ?В
EMBED Equation.3
EMBED Equation.3
Матрица А* состоит из столбцов матрицы А (канонической формы), которые соответствуют основным переменным.
х = (80, 0, 0, 10, 20, 0, 0)
Определитель матрицы А* равен:
|А*| = 0 + 0 + 1 – (16 + 0 + 0) = -15
Транспонируем матрицу А*:
Находим элементы присоединенной матрицы (Ã*):
EMBED Equation.3
EMBED Equation.3
EMBED Equation.3
EMBED Equation.3
EMBED Equation.3
EMBED Equation.3
EMBED Equation.3
EMBED Equation.3
EMBED Equation.3
Получили присоединенную матрицу:
EMBED Equation.3
EMBED Equation.3
?х1 = EMBED Equation.3 х1* = 80 EMBED Equation.3 = EMBED Equation.3
?х4 = EMBED Equation.3 х4* = 10 + EMBED Equation.3 = EMBED Equation.3
?х5 = EMBED Equation.3 х5* = 20 + EMBED Equation.3 = EMBED Equation.3
Х* = ( EMBED Equation.3 ; 0; 0; EMBED Equation.3 ; EMBED Equation.3 ; 0; 0)
?(Х*) = 5* EMBED Equation.3 + 7*0 + 3*0 + 6* EMBED Equation.3 = EMBED Equation.3
??(Х*) = 5* EMBED Equation.3 + 6* EMBED Equation.3 = EMBED Equation.3
Таким образом, при увеличении запасов сырья I и II вида на 8 и 10 ед. соответственно и уменьшении на 5 ед. запасов сырья III вида выпуск продукции вида А уменьшится на EMBED Equation.3 ед., вида Г – увеличится на EMBED Equation.3 ед. Остаток сырья I вида увеличится на EMBED Equation.3 ед.
Определим целесообразность включения в план изделия Д ценой 10 единиц, на изготовление которого расходуется по две единицы каждого вида сырья.
С5 = 10
а15 = 2; а25 = 2; а35 = 2
?j = ?аijyj – cj
?5 = а15y1 – а25y2 – а35y3 – c5
?5 = EMBED Equation.3
?5 < 0, следовательно изделие вида Д выпускать выгодно.
Так как затраты на ресурсы при выпуске одного данного изделия вида Д будут меньше цены изделия вида Д.
Задача 3
Исследовать динамику экономического показателя на основе анализа одномерного временного ряда
В течение девяти последовательных недель фиксировался спрос Y(t) (млн. руб.) на кредитные ресурсы финансовой компании. Временной ряд Y(t) этого показателя приведен в таблице.
Требуется:
1. Проверить наличие аномальных наблюдений.
2. Построить линейную модель Y(t) = а0 + а1t, параметры которой оценить МНК.
3. Построить адаптивную модель Брауна Y(t) = a0 +a1k с параметром сглаживания а = 0,4 и а = 0,7; выбрать лучшее значение параметра сглаживания
4. Оценить адекватность построенных моделей, используя свойства независимости остаточной компоненты, случайности и соответствия нормальному закону распределения (при использовании R/S-критерия взять табулированные границы 2,7 – 3,7).
5. Оценить точность моделей на основе использования средней относительной ошибки аппроксимации.
6. По двум построенным моделям осуществить прогноз спроса на следующие две недели (доверительный интервал прогноза рассчитать при доверительной вероятности р = 70%).
7. Фактические значения показателя, результаты моделирования и прогнозирования представить графически.
Решение:
1. Проверим наличие аномальных наблюдений.
Для изобразим на графике фактические уровни.
EMBED Excel.Chart.8 \s
Рис.2. Фактические уровни спроса в течение девяти последовательных недель
Резко отличающихся отдельных значений нет, следовательно аномальных наблюдений в данном временном ряду нет.
2. Построим линейную трендовую модель y = а0 + а1t, определив ее параметры с помощью метода наименьших квадратов.
Для этого решим систему нормальных уравнений:
na0 + а1 ?t = ?yt
a0 ?t + a1 ?t2 = ?tyt
Составим расчетную таблицу:
9а0 + 45 а1 = 132 *(-5)
45 а0 + 285 а1 = 822
+ -45 а0 – 225 а1 = -660
45 а0 + 285 а1 = 822
60 а1 = 162
а1 = 2,7
EMBED Equation.3
Уравнение линейного тренда имеет вид: y = 2,7 + 1,2 t
4. Оценим адекватность построенной модели.
а) Проверим случайность значений остатков по критерию пиков (поворотных точек):
EMBED Equation.3
m – количество поворотных точек;
[ … ] – целая часть числа.
Уровень Et считается поворотной точкой, если он меньше (или больше) двух рядом с ним стоящих уровней.
График остатков имеет вид:
EMBED Excel.Chart.8 \s
Рис.3. График остатков
Получили пять поворотных точек, т.е. m = 5.
EMBED Equation.3
m > 2 ? свойство выполняется; значения остатков случайные.
б) Проверим независимость уровней ряда остатков по критерию Дарбина-Уотсона (критические уровни: d1 = 0,82 и d2 = 1,32):
EMBED Equation.3
dрасч. > 2 ? рассчитываем d’ = 4 - dрасч. = 4 – 2,21 = 1,79
d2 ? d’ ? 2 ? свойство выполняется; остатки независимы, автокорреляция отсутствует.
в) Проверим нормальность распределения остатков по R/S-критерию (критические уровни: 2,7 – 3,7):
EMBED Equation.3
Emax = 2,2; Emin = -2,5
EMBED Equation.3
EMBED Equation.3
Т.к. значение R/S-критерия попадает в интервал 2,7 – 3,7, то свойство выполняется; остаточная компонента подчинена нормальному закону распределения.
Таким образом, последовательностью остатков выполняются все свойства по выбранным критериям, следовательно, можно сделать вывод о том, что модель
y = 2,7 + 1,2 t является адекватной.
5. Оценим точность построенной трендовой модели.
Для оценки точности модели найдем среднее квадратическое отклонение от линии тренда Sy и среднюю относительную ошибку аппроксимации S:
EMBED Equation.3
где р – количество параметров модели, р = 2.
EMBED Equation.3
Т.е. в среднем расчетные значения по линейному тренду отличаются от фактических уровней на 8,9%.
Т.к. S ? 5%, то модель является не точной.
6. Построим точечный и интервальный прогнозы спроса на следующие две недели.
y = 1,2 + 2,7 t
Точечный прогноз:
y(10) = 1,2 + 2,7*10 = 28,2 (млн. руб.)
y(11) = 1,2 + 2,7*11 = 30,9 (млн. руб.)
Интервальный прогноз:
yn+k = t? * Sпрогноз.
t? = 1,05
k – период упреждения;
t? * Sпрогноз. – ширина доверительного интервала;
Sпрогноз. – средняя квадратическая ошибка прогноза.
EMBED Equation.3
EMBED Equation.3
tn+k = n + k
1 шаг k = 1
tn+k = 9 + 1 = 10
EMBED Equation.3
28,2 ± 1,05 * 1,786
(26,3 – 30,1) – интервальный прогноз на один шаг.
2 шаг k = 2
tn+k = 9 + 2 = 11
EMBED Equation.3
30,9 ± 1,05 * 1,890
(28,9 – 32,9) – интервальный прогноз на второй шаг.
7. Результаты моделирования и прогнозирования представим на графике.
EMBED Excel.Chart.8 \s
Рис.4. Динамика спроса и прогноз на две надели вперед
3. Построим адаптивную модель Брауна y = а0 + а1k с параметром сглаживания а = 0,4 и а = 0,7.
Начальные оценки параметров получим по первым пяти точкам при помощи метода наименьших квадратов.
EMBED Equation.3
EMBED Equation.3
EMBED Equation.3 EMBED Equation.3
EMBED Equation.3 EMBED Equation.3
EMBED Equation.3
где ? = 1 – ? – коэффициент дисконтирования.
Возьмем k = 1 и ? = 0,4, следовательно, ? = 1 – 0,4 = 0,6.
Подробно покажем расчет на первых двух шагах, остальное представим в таблице.
t = 1 EMBED Equation.3
t = 2 EMBED Equation.3
и т.д.
На последнем шаге получена модель: Yр(N+k) = 24,6 + 2,4 k.
Среднее квадратическое отклонение равно:
EMBED Equation.3
Построим адаптивную модель Брауна с параметром сглаживания ? = 0,7, т.е. ? = 1 – 0,7 = 0,3.
t = 1 EMBED Equation.3
t = 2 EMBED Equation.3
t = 3 EMBED Equation.3
t = 4 EMBED Equation.3
t = 5 EMBED Equation.3
t = 6 EMBED Equation.3
t = 7 EMBED Equation.3
t = 8 EMBED Equation.3
t = 9 EMBED Equation.3
На последнем шаге получена модель: Yр(N+k) = 23,5 + 0,9 k.
Среднее квадратическое отклонение равно:
EMBED Equation.3
Т.к. Sy(1) < Sy(2), то точнее первая модель Yр(N+k) = 24,6 + 2,4 k с параметром сглаживания ? = 0,4.
Оценим качество первой модели Yр(N+k) = 24,6 + 2,4 k.
4. Оценим адекватность построенной модели.
а) Проверим случайность значений остатков по критерию пиков (поворотных точек).
График остатков имеет вид:
EMBED Excel.Chart.8 \s
Рис.5. График остатков
Получили пять поворотных точек, т.е. m = 5.
EMBED Equation.3
m > 2 ? свойство выполняется; значения остатков случайные.
б) Проверим независимость уровней ряда остатков по критерию Дарбина-Уотсона (критические уровни: d1 = 0,82 и d2 = 1,32):
EMBED Equation.3
d2 ? dрасч ? 2 ? свойство выполняется; остатки независимы, автокорреляция отсутствует.
в) Проверим нормальность распределения остатков по R/S-критерию (критические уровни: 2,7 – 3,7):
EMBED Equation.3
Emax = 1,5; Emin = -4,6
EMBED Equation.3
EMBED Equation.3
Т.к. значение R/S-критерия попадает в интервал 2,7 – 3,7, то свойство выполняется; остаточная компонента подчинена нормальному закону распределения.
Таким образом, последовательностью остатков выполняются все свойства по выбранным критериям, следовательно, модель Yр(N+k) = 24,6 + 2,4 k является адекватной.
5. Оценим точность построенной модели.
Для оценки точности модели найдем среднюю относительную ошибку аппроксимации S:
EMBED Equation.3
Т.к. S > 5%, то модель не является точной.
6. Построим точечный и интервальный прогнозы спроса на следующие две недели.
Yр(N+k) = 24,6 + 2,4 k
Точечный прогноз:
y(9+1) = 24,6 + 2,4*1 = 27,0
y(9+2) = 24,6 + 2,4*2 = 29,4
Интервальный прогноз:
t? = 1,05
1 шаг k = 1
tn+k = 9 + 1 = 10
EMBED Equation.3
27,0 ± 1,05 * 2,565
(24,3 – 29,7) – интервальный прогноз на один шаг.
2 шаг k = 2
tn+k = 9 + 2 = 11
EMBED Equation.3
29,4 ± 1,05 * 2,714
(26,6 – 32,3) – интервальный прогноз на второй шаг.
7. Результаты моделирования и прогнозирования представим на графике.
EMBED Excel.Chart.8 \s
Рис.6. Динамика спроса и прогноз на две надели вперед
Задача 4
Провести моделирование и решить специальную задачу линейного программирования
Транспортная задача
Компания, занимающаяся ремонтом автомобильных дорог, в следующем месяце будет проводить ремонтные работы на пяти участках автодорог. Песок на участки ремонтных работ может доставляться из трех карьеров, месячные объемы предложений по карьерам известны. Из планов производства ремонтных работ известны месячные объемы потребностей по участкам работ. Имеются экономические оценки транспортных затрат (в у.е.) на перевозку 1 тонны песку с карьеров на ремонтные участки.
Требуется предложить план перевозок песка на участки ремонта автодорог, который обеспечивает минимальные совокупные транспортные издержки.
Матрица планирования:
Решение:
Для того чтобы решить транспортную задачу, воспользуемся 2я этапами реализации метода потенциалов для транспортной задачи.
Этап 1. Первоначальное закрепление участков работ за карьерами.
Рассмотрим два метода получения начального распределения (начального опорного плана): метод северо-западного угла и метод наименьших стоимостей.
Метод северо-западного угла:
В первую очередь заполняется клетка (из числа невычеркнутых), стоящая в верхнем левом (северо-западном) углу матрица планирования. Результат показан в табл. 1: заполняется клетка (1;1) и вычеркивается первый столбец, заполняется клетка (1;2) и вычеркивается столбец, заполняется клетка (2;2) и вычеркивается столбец, заполняется клетка (2;3) и вычеркивается строка, заполняется клетка (3;3) и вычеркивается столбец, заполняется клетка (3;4) и вычеркивается столбец, заполняется клетка (3;5) и вычеркивается последние строка и столбец. Число занятых клеток равно m + n = 3 + 5 – 1 = 7. Суммарные затраты на реализацию данного плана перевозок составят:
f(x) = 4*40 + 2*20 + 4*30 + 3*60 + 4*90 + 6*50 + 2*0 = 1160.
Таблица 1.
Недостатком данного метода является то, что не учитываются значения элементов EMBED Equation.3 матрицы транспортных расходов, в результате чего полученное этим методом начальное распределение (начальный опорный план перевозок) может быть достаточно далеко от оптимального. Поэтому используем метод наименьших стоимостей.
Метод наименьших стоимостей:
В модификациях метода наименьших стоимостей заполнение клеток матрицы планирования проводится с учетом значений величин EMBED Equation.3 . Так, в модификации «двойного предпочтения» отмечают клетки с наименьшими стоимостями перевозок сначала по каждой строке, а затем по каждому столбцу. Из двух клеток с одинаковой стоимостью перевозок предпочтение отдается клетке, через которую осуществляется больший объем перевозок. Вычеркивание строк и столбцов при заполнении клеток проводится по описанным выше правилам. Результат начального распределения методом наименьших стоимостей представлен в таблице 2.
Таблица 2.
Порядок заполнения клеток: (1;5), (2;1), (1;2), (2;3), (3;3), (3;2), (3;4). Суммарные затраты на перевозки, представленные в таблице 2. составляют
f(x) = 2*10 + 1*50 + 2*40 + 3*50 + 5*30 + 4*90 + 6*20 = 930.
Следовательно, данный план перевозок значительно ближе к оптимальному, чем план, составленный по методу северо-западного угла.
Этап 2. Проверка оптимальности полученного плана перевозок.
Введем специальные показатели EMBED Equation.3 для каждой строки и показатели EMBED Equation.3 для каждого столбца. Эти показатели называются потенциалами. Потенциалы подбираются таким образом, чтобы для заполненной клетки (i, j) выполнялось равенство EMBED Equation.3 (1)
Рассмотрим процесс нахождения потенциалов для базисного начального распределения по методу северо-западного угла представленного в таблице 1. Задав EMBED Equation.3 = 0 и используя формулу (1) для заполненных клеток (1;1) и (1;2), находим EMBED Equation.3 =4, EMBED Equation.3 =2, зная EMBED Equation.3 , по заполненной клетке (2;2) находим EMBED Equation.3 = -2, а зная EMBED Equation.3 , по заполненной клетке (2;3) находим EMBED Equation.3 = 1. Зная EMBED Equation.3 , по заполненной клетке (3;3) находим EMBED Equation.3 = -3, зная EMBED Equation.3 , по заполненным клеткам (3;4), (3;5) находим EMBED Equation.3 = 3 EMBED Equation.3 = -1
Результаты представлены в таблице 3., где потенциалы приведены в последнем столбце и последней строке.
Таблица 3.
Аналогичные результаты для начального распределения по методу наименьших стоимостей, приведенного в таблице 2., представлены в таблице 4.
Таблица 4.
Чтобы оценить оптимальность распределения, для всех клеток (i;j) матрицы планирования определяются их оценки, которые обозначим через EMBED Equation.3 , по формуле: EMBED Equation.3 (2)
Условием оптимальности распределения служит условие неотрицательности оценок свободных клеток матрицы планирования. Оценки клеток по формуле (2) удобно представить в виде матрицы оценок. Для рассматриваемого распределения, полученного методом северо-западного угла, матрица оценок клеток имеет вид:
0 0 2 1 2
EMBED Equation.3 -4 0 0 0 5
-1 0 0 0 0
Наличие отрицательных оценок свободных клеток свидетельствует о том, что данный план перевозок далек от оптимального (суммарные затраты на перевозку по этому плану равны 1160).
Для распределения, полученного методом наименьших стоимостей (таблица 4), матрица оценок имеет вид:
4 0 2 1 0
EMBED Equation.3 0 0 0 0 3
3 0 0 0 0
Так как все оценки неотрицательны, то не имеется возможности улучшить данный план, т.е. он оптимален (суммарные затраты на перевозку по этому плану равны 930). Наличие нулевых оценок свободных клеток в оптимальном плане перевозок, свидетельствует о неединственности оптимального плана.
Приложение 1.
Для того чтобы получить оптимальный производим ввод исходных данных (рис. 1.1)
Рис. 1.1
1) Вводим зависимость для целевой ячейки F4. С помощью «Мастера функции», функция СУММПРОИЗВ вводим данные (рис. 1.2) и получаем введенную функцию в ячейку F4.
Рис. 1.2
2) Копируем значения ячейки F4 в ячейки F8, F9, F10.
3) В ячейке F4 находим значение с помощью Поиска решения, введя данные показанные на рис. 1.3.
рис. 1.3
