Экономико-математические методы и модели в отрасли связи
Федеральное агентство связи
Сибирский
Государственный Университет
Межрегиональный
центр переподготовки
специалистов
Контрольная
работа
По
дисциплине:
Экономико-математические методы и модели
в отрасли связи
Выполнил: студент
Группа: ЭПТ-92
Проверил: Батый Ада Рамазановна
Новосибирск,
2011 г.
ЗАДАЧА 1.
На территории города имеется три телефонных станции А, Б и В. Незадействованные емкости станций составляют на станции А - QА, Б - QБ, В - QВ номеров (таблица 1.1). Потребности новых районов застройки города в телефонах составляют: 1 - q1, 2 - q2, 3 - q3, 4 - q4 номеров (таблица 1.2).
Необходимо
составить экономико-
Исходные данные
Таблица 1 - Незадействованные ёмкости телефонных станций.
| Возможности станций, номеров | |
| QА | 600 |
| QБ | 400 |
| QВ | 700 |
Таблица 2 - Спрос на установку телефонов
| Спрос районов, номеров | |
| q1 | 350 |
| q2 | 400 |
| q3 | 500 |
| q4 | 450 |
Таблица 3 - Среднее расстояние от станции до районов застройки, км (для всех вариантов)
| Станции | Районы | |||
| 1 | 2 | 3 | 4 | |
| А | 4 | 5 | 6 | 4 |
| Б | 3 | 2 | 1 | 4 |
| В | 6 | 7 | 5 | 2 |
РЕШЕНИЕ
Стоимость доставки единицы груза из каждого пункта отправления в соответствующие пункты назначения задана матрицей тарифов
Распределительный метод является одним из вариантов базового симплексного метода. Поэтому идея распределительного метода (как и симплексного) содержит такие же три существенных момента.
Прежде всего, отыскивается какое-то решение задачи — исходный опорный план. Затем посредством специальных показателей опорный план проверяется на оптимальность. Если план оказывается не оптимальным, переходят к другому плану. При этом второй и последующие планы должны быть лучше предыдущего. Так за несколько последовательных переходов от не оптимального плана приходят к оптимальному.
| Станции | 1 | 2 | 3 | 4 | Районы |
| 1 | 4 | 5 | 6 | 4 | 600 |
| 2 | 3 | 2 | 1 | 4 | 400 |
| 3 | 6 | 7 | 5 | 2 | 700 |
| 350 | 400 | 500 | 450 |
Проверим
необходимое и достаточное
∑ a = 600 + 400 + 700 = 1700
∑ b = 350 + 400 + 500 + 450 = 1700
Условие баланса соблюдается. Запасы равны потребностям. Следовательно, модель транспортной задачи является закрытой.
Занесем исходные данные в распределительную таблицу.
| Станции | 1 | 2 | 3 | 4 | Районы |
| 1 | 4 | 5 | 6 | 4 | 600 |
| 2 | 3 | 2 | 1 | 4 | 400 |
| 3 | 6 | 7 | 5 | 2 | 700 |
| 350 | 400 | 500 | 450 |
Первая итерация заключается в определении исходного опорного плана и проверке его на оптимальность.
Определение исходного опорного плана. Первый опорный план может быть найден посредством различных способов: по правилу северо-западного угла, приоритету ближайших пунктов, способу минимального элемента С=(cij), способу Фогеля и по способу Лебедева-Тихомирова.
Этап I. Поиск первого опорного плана.
1. Используя метод северо-западного угла, построим первый опорный план транспортной задачи.
План начинается заполняться с верхнего левого угла.
Искомый элемент равен 4
Для этого элемента запасы равны 600, потребности 350. Поскольку минимальным является 350, то вычитаем его.
x11
= min(600,350) = 350.
| 4 | 5 | 6 | 4 | 600 - 350 = 250 |
| x | 2 | 1 | 4 | 400 |
| x | 7 | 5 | 2 | 700 |
| 350 - 350 = 0 | 400 | 500 | 450 | 0 |
Искомый элемент равен 5
Для этого элемента запасы равны 250, потребности 400. Поскольку минимальным является 250, то вычитаем его.
x12 = min(250,400) = 250.
| 4 | 5 | x | x | 250 - 250 = 0 |
| x | 2 | 1 | 4 | 400 |
| x | 7 | 5 | 2 | 700 |
| 0 | 400 - 250 = 150 | 500 | 450 | 0 |
Искомый элемент равен 2
Для этого элемента запасы равны 400, потребности 150. Поскольку минимальным является 150, то вычитаем его.
x22 = min(400,150) = 150.
| 4 | 5 | x | x | 0 |
| x | 2 | 1 | 4 | 400 - 150 = 250 |
| x | x | 5 | 2 | 700 |
| 0 | 150 - 150 = 0 | 500 | 450 | 0 |
Искомый элемент равен 1
Для этого элемента запасы равны 250, потребности 500. Поскольку минимальным является 250, то вычитаем его.
x23 = min(250,500) = 250.
| 4 | 5 | x | x | 0 |
| x | 2 | 1 | x | 250 - 250 = 0 |
| x | x | 5 | 2 | 700 |
| 0 | 0 | 500 - 250 = 250 | 450 | 0 |
Искомый элемент равен 5
Для этого элемента запасы равны 700, потребности 250. Поскольку минимальным является 250, то вычитаем его.
x33 = min(700,250) = 250.
| 4 | 5 | x | x | 0 |
| x | 2 | 1 | x | 0 |
| x | x | 5 | 2 | 700 - 250 = 450 |
| 0 | 0 | 250 - 250 = 0 | 450 | 0 |
Искомый элемент равен 2
Для этого элемента запасы равны 450, потребности 450. Поскольку минимальным является 450, то вычитаем его.
x34 = min(450,450) = 450.
| 4 | 5 | x | x | 0 |
| x | 2 | 1 | x | 0 |
| x | x | 5 | 2 | 450 - 450 = 0 |
| 0 | 0 | 0 | 450 - 450 = 0 | 0 |
| Станции | 1 | 2 | 3 | 4 | Районы |
| 1 | 4[350] | 5[250] | 6 | 4 | 600 |
| 2 | 3 | 2[150] | 1[250] | 4 | 400 |
| 3 | 6 | 7 | 5[250] | 2[450] | 700 |
| 350 | 400 | 500 | 450 |
В
результате получен первый опорный
план, который является допустимым,
так как все грузы из баз
вывезены, потребность магазинов
удовлетворена, а план соответствует
системе ограничений
2.
Подсчитаем число занятых
Значение целевой функции для этого опорного плана равно:
4*350 + 5*250 + 2*150 + 1*250 + 5*250 + 2*450 = 5350
Этап II. Улучшение опорного плана.
Проверка опорного плана на оптимальность. Чтобы установить является ли опорный план оптимальным, надо проверить, как повлияет на величину целевой функции любое возможное перераспределение поставок.
План распределения поставок будет оптимальным лишь в том случае, когда целевая функция имеет минимальное значение, т.е. когда дальнейшее уменьшение затрат на поставку будет невозможно.
Проверим возможность уменьшения суммарных затрат на поставку продукции. С этой целью для каждой свободной от поставки клетки определяется величина Δij, характеризующая изменение суммарных затрат на поставку (в расчете на единицу перераспределяемой продукции), при условии включения в план единичной поставки хij=1 от поставщика Аi к потребителю Вj.
При этом должно быть произведено такое изменение остальных поставок, чтобы получившаяся совокупность поставок не нарушала баланса спроса и поставок транспортной задачи.
Величина Δij называется оценкой свободной клетки (или характеристика).
В исходном решении задачи имеются клетки свободные от поставок.
Необходимо вычислить значение оценок Δij для этих свободных от поставок клеток. С этой целью для каждой свободной клетки составляется означенный цикл перерасчета (или замкнутая цепь, круг, кольцо, контур и т.д.).
Под циклом пересчета (цепью) понимается замкнутая ломаная линия. Вершинами цикла (цепи) являются клетки таблицы, проще – вершины лежат в клетках таблицы.
Причем одна из вершин находится в свободной от поставки клетке, в той, для которой определяется оценка Δij. Все другие вершины находятся в базисных клетках, т.е. клетках, занятых поставками.
Вершины,
в которых поставки при перераспределении
увеличиваются, отмечаются плюсом и
называются положительными вершинами
и, наоборот, вершины, в которых поставки
при перераспределении
В цикле знаки по вершинам расставляют начиная с вершины, лежащей в свободной клетке, для которой определяется Δij. В нее записывают знак плюс, затем знаки по вершинам чередуются: минус, плюс , минус, плюс и т. д., независимо от того, расставляют ли их по часовой стрелке или в обратном направлении. Таким образом, в цикле всегда насчитывается одинаковое число положительных и отрицательных вершин.
Следующий этап решения транспортной задачи заключается в улучшении опорного плана.
Если при каком-то опорном плане оказывается несколько свободных клеток с отрицательными оценками Δij, то за один переход к лучшему плану можно занять поставкой только одну клетку – ту, которая обеспечивает наибольшее снижение целевой функции.
Шаг 1. Определяем оценку для каждой свободной клетки.
(1;3): В свободную клетку (1;3) поставим знак «+», а в остальных вершинах многоугольника чередующиеся знаки «-», «+», «-».
| Станции | 1 | 2 | 3 | 4 | Районы |
| 1 | 4[350] | 5[250][-] | 6[+] | 4 | 600 |
| 2 | 3 | 2[150][+] | 1[250][-] | 4 | 400 |
| 3 | 6 | 7 | 5[250] | 2[450] | 700 |
| 350 | 400 | 500 | 450 |
Цикл приведен в таблице (1,3; 1,2; 2,2; 2,3; ).
Оценка свободной клетки равна Δ13 = (6) - (5) + (2) - (1) = 2.
(1;4): В свободную клетку (1;4) поставим знак «+», а в остальных вершинах многоугольника чередующиеся знаки «-», «+», «-».
| Станции | 1 | 2 | 3 | 4 | Районы |
| 1 | 4[350] | 5[250][-] | 6 | 4[+] | 600 |
| 2 | 3 | 2[150][+] | 1[250][-] | 4 | 400 |
| 3 | 6 | 7 | 5[250][+] | 2[450][-] | 700 |
| 350 | 400 | 500 | 450 |
Цикл приведен в таблице (1,4; 1,2; 2,2; 2,3; 3,3; 3,4; ).
Оценка свободной клетки равна Δ14 = (4) - (5) + (2) - (1) + (5) - (2) = 3.
(2;1): В свободную клетку (2;1) поставим знак «+», а в остальных вершинах многоугольника чередующиеся знаки «-», «+», «-».
| Станции | 1 | 2 | 3 | 4 | Районы |
| 1 | 4[350][-] | 5[250][+] | 6 | 4 | 600 |
| 2 | 3[+] | 2[150][-] | 1[250] | 4 | 400 |
| 3 | 6 | 7 | 5[250] | 2[450] | 700 |
| 350 | 400 | 500 | 450 |
Цикл приведен в таблице (2,1; 2,2; 1,2; 1,1; ).
Оценка свободной клетки равна Δ21 = (3) - (2) + (5) - (4) = 2.
(2;4): В свободную клетку (2;4) поставим знак «+», а в остальных вершинах многоугольника чередующиеся знаки «-», «+», «-».
| Станции | 1 | 2 | 3 | 4 | Районы |
| 1 | 4[350] | 5[250] | 6 | 4 | 600 |
| 2 | 3 | 2[150] | 1[250][-] | 4[+] | 400 |
| 3 | 6 | 7 | 5[250][+] | 2[450][-] | 700 |
| 350 | 400 | 500 | 450 |
Цикл приведен в таблице (2,4; 2,3; 3,3; 3,4; ).
Оценка свободной клетки равна Δ24 = (4) - (1) + (5) - (2) = 6.
(3;1): В свободную клетку (3;1) поставим знак «+», а в остальных вершинах многоугольника чередующиеся знаки «-», «+», «-».
| Станции | 1 | 2 | 3 | 4 | Районы |
| 1 | 4[350][-] | 5[250][+] | 6 | 4 | 600 |
| 2 | 3 | 2[150][-] | 1[250][+] | 4 | 400 |
| 3 | 6[+] | 7 | 5[250][-] | 2[450] | 700 |
| 350 | 400 | 500 | 450 |
Цикл приведен в таблице (3,1; 3,3; 2,3; 2,2; 1,2; 1,1; ).
Оценка свободной клетки равна Δ31 = (6) - (5) + (1) - (2) + (5) - (4) = 1.
(3;2): В свободную клетку (3;2) поставим знак «+», а в остальных вершинах многоугольника чередующиеся знаки «-», «+», «-».
| Станции | 1 | 2 | 3 | 4 | Районы |
| 1 | 4[350] | 5[250] | 6 | 4 | 600 |
| 2 | 3 | 2[150][-] | 1[250][+] | 4 | 400 |
| 3 | 6 | 7[+] | 5[250][-] | 2[450] | 700 |
| 350 | 400 | 500 | 450 |
Цикл приведен в таблице (3,2; 3,3; 2,3; 2,2; ).
Оценка свободной клетки равна Δ32 = (7) - (5) + (1) - (2) = 1.
Из приведенного расчета видно, что ни одна свободная клетка не имеет отрицательной оценки, следовательно, дальнейшее снижение целевой функции Fx невозможно, поскольку она достигла минимального значения.
Таким образом, последний опорный план является оптимальным.
Минимальные затраты составят:
4*350 + 5*250 + 2*150 + 1*250 + 5*250 + 2*450 = 5350
Если в оптимальном решении задачи имеется несколько оценок равных нулю, то это является свидетельством того, что среди бесчисленного множества решений этой задачи существуют еще решения, являющиеся также оптимальными, поскольку значение целевой функции остается одинаковым — минимальным. Их принято называть альтернативными.
Примечание. Основной алгоритм распределительного метода является не лучшим методом решения транспортных задач, так как на каждой итерации для проверки опорного плана на оптимальность приходилось строить [mп—(m+n—1)] циклов пересчета, что при больших размерах матрицы оказывается очень громоздким и трудоемким делом. Так, для расчетов по матрице 10х10 на каждой итерации надо строить 81 цикл, а по матрице 20x20 — 361 цикл.