Задача: Предприятие производит продукцию двух видов Р1 и Р2 из сырья трех видов S1, S2, S3. Запасы сырья равны соответственно b1, b2, b3. Расход i -го вида сырья Si на единицу j - го вида продукции Pj равен aij. Доход, получаемыйпредприятием от реализации единицы j - го вида продукции, равен сj.
Найти план производства, обеспечивающий предприятию максимум дохода (значения всех указанных параметров приведены в таблице).
Вид сырья | Вид продукции | Ресурсы сырья bj (в у.е.) | |
Р1 | Р2 | ||
S1 | a11=4 | a12=5 | b1=6 |
S2 | a21=1 | а22=6 | b2=8 |
S3 | а31=1 | а32=4 | b3=4 |
Доход сj (в у.е.) |
1 шаг алгоритма: пусть х1 и х2 - объем выпускаемой продукции вида Р1иР2 (соответственно). Построим математическую модель(З.Л.П.) приведенной экономической задачи:
f (x) = 7 х1 + 5 х2 ® max
х1³0, х2³0
2 шаг алгоритма: сведем систему ограничений З.Л.П. к виду системы линейных уравнений. Для этого прибавим к левой части каждого неравенства дополнительные переменные. После этого сформулированная нами З.Л.П. примет вид: f (x) = 7 х1 + 5 х2 ® max
х1³0, х2 ³0, х3 ³0, х4 ³0, х5 ³0
|
|
3 шаг алгоритма: путем элементарных преобразований необходимопривести систему ограничений к виду системы линейных уравнений с выделенными переменными. В данном случае мы уже имеем систему с выделенными переменными х3, х4, х5. Так как:
- в первом уравнении системы коэффициент перед переменной х3 равен единице, а в других уравнениях – нулю;
- во втором уравнении коэффициент перед переменной х4 равен единице, в других уравнениях данная переменная отсутствует;
- в третьем уравнении коэффициент перед переменной х5 равен единице, а в остальных уравнениях – нулю.
4 шаг алгоритма: найдем базисное решение, приравняв для этого небазисные переменные х1 и х2 к нулю, а базисные х3, х4, х5 – к свободным членам: (х1, х2, х3, х4, х5) = (0, 0, 6, 8, 4).
5 шаг алгоритма: пользуясь Теоремой 1,проверим на оптимальность полученное базисное решение.
Проведем анализ коэффициентов функции: f (x) = 7 х1 + 5 х2. В формуле целевой функции отсутствуют слагаемые с базисными переменными х3, х4, х5, т.е. коэффициенты при них равны нулю. Коэффициенты при небазисных переменных х1 и х2 положительные (7>0 и 5>0 соответственно).Следовательно, найденное базисное решение (х1, х2, х3, х4, х5) = (0, 0, 6, 8, 4) не оптимальное.
Далее по схеме 1 произведем смену одной базисной переменной. Для этого для всех коэффициентов a11, a21, a31 (так как все они положительные!), стоящих перед переменной x1 в системе ограничений, вычислим разрешающие коэффициенты по формуле:
ri = , где bi - свободный член в i-ом уравнении.
Получим, что r1 = = 1,5; r2 = = 8; r3 = = 4.
Минимальный разрешающий коэффициент r1 = min(r1,r2,r3 )=1,5. Следовательно, в 1 -ом уравнении необходимо преобразовать переменную x1 в базисную путем эквивалентных преобразований.
|
|
1. Умножим обе части первого уравнения на действительное число λ=1/4 (чтобы коэффициент перед x1стал равным единице).
Получили:
Далее, чтобы преобразовать переменную x1 в базисную переменную, необходимо исключить ее из 2-го и 3-го уравнений. Для этого выразим в первом уравнении перемнную x1 через небазисные переменные:
x1 =
и заменим ее на полученное выражение в других уравнениях. Получим:
Приведем подобные члены и получим следующий вид системы ограничений:
Таким образом, по определению переменная x1будет являться базисной. Получили новое базисное решение: (х1, х2, х3, х4, х5) = (6/4, 0, 0, 4, 5/2), где х2, х3 – небазисные переменные (приравниваем к нулю) и х1, х4, х5 - базисные переменные (приравниваем к свободным членам).
Новое базисное решение необходимо проверить на оптимальность. При этом из целевой функции исключим новую базисную переменную x1, выразив ее через небазисные переменные в 1 -ом уравнении.
В процессе решения мы получили, что x1 = .
Следовательно,
f (x) =7 х1 + 5 х2 =7() + 5 х2 = .
Проанализируем коэффициенты перед переменными в целевой функции:
- коэффициенты перед базисными переменными х1, х4, х5 равны нулю (так как отсутствуют слагаемые с этими переменными);
- коэффициенты перед не базисными переменными х2, х3 - отрицательные (-15/4 и –1/4 соответственно).
Следовательно, критерий оптимальности базисного решения выполнен, т.е. решение (х1, х2, х3, х4, х5) = (6/4, 0, 0, 4, 5/2) есть оптимальное решение и max(f(x)) = .
6 шаг алгоритма: интерпретируем полученные результаты.
1. План производства, обеспечивающий предприятию максимум дохода, состоит в выпуске продукции вида Р1 объемом в 6/4 у.е. (так как х1 = 6/4 ) и продукции вида Р2 объемом в 0 у.е. (так как х2 = 0).
2. Максимальная прибыль производства будет равна 10,5 условных денежных единиц (так как max(f(x))=10,5 ).
Замечание*:при увеличении числа неизвестных запись приведенного выше решения З.Л.П. становится громоздкой. Поэтому часто решение З.Л.П. симплекс - методом оформляют в виде таблицы. Такая таблица называется симплекс – таблицей. |
Замечание**: при использовании симплекс-таблиц все слагаемые целевой функции, содержащие переменные, необходимо перенести в левую часть уравнения с изменением знака:
f (x) = с1x1+с2x2 +с3x3+…+сnxn ® f (x) -с1x1-с2x2 -с3x3-…-сnxn = 0.
Замечание***: изменение вида записи целевой функции сопряжено с некоторыми изменениями в условиях основных теорем (Теоремы 1, 2, 3). Учитывая данный факт, сформулируем следующие теоремы.
Теорема 1* (Критерий оптимальности неотрицательного базисного решения в задачах З.Л.П. на нахождение максимума целевой функции): Пусть в З.Л.П. с выделенными переменными все коэффициенты при небазисных переменных в целевой функции неотрицательные, а при базисных равны нулю. Тогда базисное решение, соответствующее данному виду З.Л.П., является оптимальным.
Теорема 2* (Критерий неограниченности значений целевой функции): Пусть в З.Л.П. (на максимум) с выделенными переменными некоторый коэффициент в целевой функции с к < 0 при переменной хк, и пусть в системе ограничений перед переменной хк все коэффициенты а’iк £ 0. Тогда значения целевой функции не ограничены на максимум (f max (x) = ¥).
Теорема 3 (О возможности перехода от одного базисного решения к другому с не меньшим значением целевой функции): Пусть в З.Л.П. (на максимум) с выделенными переменными в целевой функции некоторый коэффициент с к < 0 при переменной хк, и в системе ограничений среди коэффициентов перед этой переменной есть положительные (а’iк >0). Тогда можно перейти от данного базисного решения к другому, с не меньшим значением целевой функции.