Открытая модель. Метод потенциалов.
В результате геологических разработок и усовершенствований наших ученых добыча газа на месторождении Янбург увеличилась до 180 млн. м3/сутки
Получилась открытая модель, где Σi ai<Σj bj, для того, чтобы уравновесить модель можно ввести еще одну фиктивную страну, затраты на транспортировку в которую равны 0 (там будут собираться запасы). Объем потребления этой фиктивной страны должен быть равен | Σj bj - Σi ai |
Найдем начальный опорный план любым из методов:
Добыча
Продажа |
М1 |
М2 |
М3 |
|
170 |
160 |
180 |
||
S1 |
150 |
50 |
30 |
3150 |
S2 |
90 |
370 |
420 |
50 |
S3 |
70 |
270 |
20 |
30 |
S4 |
140 |
40 |
2140 |
30 |
S5 |
60 |
030 |
00 |
030 |
F=150*3+70*3+20*4+70*2+140*2+30*0+30*0=1160
Чтобы план был оптимален необходимо:
1) Сумма потенциалов должна быть равна cij
2) Для каждой занятой клетки сумма потенциалов ≤cij
Построим для каждой свободной переменной плана числа γij=cij-(Ui+Vj) у оптимального решения они все должны быть ≥0 т.к. число потенциалов = 8, а система состоит из 7 уравнений, то для нахождения решения какому-либо потенциалу нужно придать произвольное значение, пусть U2=0. Исходя из этого найдем другие потенциалы:
Добыча
Продажа |
М1 |
М2 |
М3 |
|
|
170 |
160 |
180 |
|||
S1 |
150 |
50 |
30 |
3150 |
U1=0 |
S2 |
90 |
3 70+ |
420- |
50 |
U2=0 |
S3 |
70 |
270 |
20 |
30 |
U3=-1 |
S4 |
140 |
40 |
2140 |
30 |
U4=-2 |
S5 |
60 |
030- |
00+ |
030 |
U5=-3 |
|
V1=3 |
V2=4 |
V3=3 |
|
Теперь найдем все γij для свободных переменных.
γ11=5-(0+3)=2≥0
γ12=3-(0+4)=-1<0
γ23=5-(0+3)=2≥0
γ32=2-(4-1)=-1<0
γ33=3-(3-1)=1≥0
γ41=4-(3-2)=3≥0
γ43=3-(3-2)=2≥0
γ52=0-(4-3)=-1<0
Видим, что γ12, γ32, γ52 <0 , значит, опорный план можно улучшить. Т.к. среди γ12, γ32, γ52 минимального нет, то вектор в базис можно вводить в любой из трех клеток. Пусть это будет клетка (5;2) Но, так как добавляя нагрузку на (5;2) нарушаем баланс строк и столбцов распределительной таблицы, то следует изменить объемы поставок в ряде других занятых клеток. А чтобы число базисных переменных осталось прежним, необходимо вывести из базиса одну переменную. Выводим ту, у которой загрузка в цикле минимальна.
Строим цикл замен: начальная клетка (5;2) приписываем этой клетке знак + (увеличиваем нагрузку), и чередуем знаки в клетках (если увеличили нагрузку на одну клетку, то для достижения равновесия системы надо уменьшить нагрузку на другой). (5;2)+ - (5;1)- - (2;1)+ - (2;2)-
Нагрузка увеличивается и уменьшается на одну и ту же величину θ=min(xij), где xij-содержимое клеток со знаком минус.
Получим новый опорный план и проверим его на оптимальность:
Добыча
Продажа |
М1 |
М2 |
М3 |
|
|
170 |
160 |
180 |
|||
S1 |
150 |
50 |
30 |
3150 |
U1=3 |
S2 |
90 |
390 |
40 |
50 |
U2=3 |
S3 |
70 |
270 |
20 |
30 |
U3=2 |
S4 |
140 |
40 |
2140 |
30 |
U4=2 |
S5 |
60 |
010 |
020 |
030 |
U5=0 |
|
V1=0 |
V2=0 |
V3=0 |
|
F=150*3+90*3+70*2+140*2=1140
γ11=5-(0+3)=2≥0
γ12=3-(0+3)=0≥0
γ22=4-(0+3)=1≥0
γ23=5-(0+3)=2≥0
γ32=2-(0+2)=0≥0
γ33=3-(0+2)=1≥0
γ41=4-(0+2)=2≥0
γ43=3-(0+2)=1≥0
Видим, что все γ≥0, а значит план больше улучшить нельзя, получили оптимальное решение! F=1140