Файл: Падалко Л.П. Математические методы оптимального планирования развития и эксплуатации энергосистем учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 04.07.2024

Просмотров: 118

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Таблица 1-2

 

 

 

Свободные переменные

Базисные

Свободные

 

 

переменные

члены

х,

xs

 

 

4 + 1

P i —

 

aIs—^ais'xlj

Xk+l

а д — X a / l a-lj

als—lals*ij

а1к—Л"lkal і

xj

X

Xk+m

 

F

Ik—Ішу

 

Т а к им образом, при использовании симплексных т а б ­

лиц

нет необходимости к а ж д ы й раз

проделывать

с л о ж ­

ные

алгебраические преобразования .

Стандартный

алго­

ритм перехода от одной таблицы к другой позволяет до ­ статочно просто программировать решение на Э Ц В М . П р и з н а к о м оптимальности решения является неположи ­ тельность всех элементов последней строки таблицы .

§1 . 5 .Отыскание допустимого базисного решения

Пр и изложении вычислительной схемы симплексного метода мы требовали, чтобы система ограничений была приведена к виду (1.2). Иными словами, предваритель ­ ным шагом к решению задачи симплексным методом я в ­ ляется нахождение допустимого базисного решения. В о

многих з а д а ч а х базисное

решение находится

весьма про­

сто, как, в частности, в

используемом нами

иллюстра ­

тивном примере. В других з а д а ч а х нахождение исходно­ го допустимого решения оказывается более сложным . Н е останавливаясь подробно на этом вопросе, мы отметим один из приемов отыскания допустимого базисного реше­ ния, называемый обычно методом искусственного базиса.

Пусть система ограничений записана в общем виде:

anXi

+ a12x2

+...+а1пхп

 

=

Ьг;

 

 

 

&21Х1

OonXv, -\- ... -f- QonXn

=

Ь2 і

 

 

/1 л \

a m l X l

+ am"X1

+

•••

+ amnXn

=

bm- .

 

 

 

Свободные члены д о л ж н ы

быть

неотрицательными.

Если для какого-либо уравнения это не выполняется,

мы

у м н о ж а е м его на —1.

 

 

 

 

 

 

 

 

Введем вспомогательные, или искусственные, пере­

менные

уи у%

 

ут,

связанные

с Х\,

х2,

хп уравне ­

ниями:

 

 

 

 

 

 

 

 

 

 

 

у1 =

Ь1 — (ахх

+

а1 2 л-2

+

... +

а1пхп);

 

 

У-2 =

Ь» — ( « 2 1 * 1

+

й22Х2

+

•••

+ а2„Хп)'<

 

/1

сч

У,п =

Ьт — (атЛхх

+ ат,х,

+

... +

атпхп).

 

 

31


Очевидно, решению системы

(1.4) д о л ж н о

отвечать

решение системы (1.5) с условием

г/і=0,

(/2 = 0,

ут = 0 .

В системе (1.5) переменные

у2,

Ут

образуют

базис. Предположим, что нам удалось, отправляясь от

этого

базиса,

перейти

к

другому,

не

с о д е р ж а щ е м у

уг,

І/2,

Ут- Иными словами,

допустим,

что после

преобра­

зований

мы получим систему:

 

 

 

 

 

 

 

 

* i =

P

i — ( a i i - V m + H -

•••

+ a i n * „

+ Y

i i #

i

+

. . .

Ут)\

x2=$2—{a.llxm,l+

... +аіпхппу1+

 

 

 

...

+y.lmym);

 

Xm = $m—(amlxm + l+

••• +атпХп

+ УтіУі +

- +УттУт)-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1-6)

П о л а г а я здесь tj\,

г/г,

Ут

равными нулю,

получим

систему:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

* і = Рі ( « і Л + і +

-

+ « і п - 0 ;

і

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.7)

Л =

 

Рш

 

( а ш А ! 1 ~Г

•••

" Т " а т ; Л і ) .

j

 

 

 

 

 

 

Э т а система

д о л ж н а быть равносильна

 

системе

(1.4), и

поэтому з а д а ч а нахождения допустимого базисного

ре­

шения

выполнена.

 

 

 

 

 

 

 

 

 

 

 

Д л я

того

чтобы перейти от

системы

 

(1.4)

к

системе

.(1.7),

необходимо решить

задачу

минимизации:

 

 

min F = yl

+ ya + ... + уп

 

 

 

 

 

 

 

(1.8)

при ограничениях (1.5).

 

 

 

 

 

 

 

 

 

 

Если

при

рассмотрении

последней

симплекс-таблицы

о к а ж е т с я ,

что переменные

у\, г/2,

Ут

 

 

будут

свободны­

ми, то это

означает, что

мы пришли к системе

(1.6),

т. е.

з а д а ч а

решена. Ц е л е в а я

функция

(1.8)

приобретает

при

этом нулевое

значение.

 

 

 

 

 

 

 

 

 

 

Если

ж е о к а ж е т с я ,

что min F

>

0, то это означает,

что

система

(1.4)

не имеет неотрицательных

 

решений. Отсю­

да вытекает,

что и задача

линейного

программирования

т а к ж е

не имеет решения.

 

 

 

 

 

 

 

 

 

 

32


§ 1.6. Теория двойственности в линейном

программировании

Р а с с м о т р им следующую

з а д а ч у

об

использовании

сырья. Д л я

производства

двух

видов продукции

исполь­

зуются четыре вида сырья . З а п а с ы сырья

и его расход на

изготовление единицы к а ж д о г о

вида продукции

з а д а ю т с я

табл .

1.3. Коэффициент

аи- показывает,

сколько

требу­

ется

сырья

1-го вида д л я

производства

единицы

 

продук­

ции вида /.

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1.3

 

 

 

 

Продукция

 

 

Виды

сырья

Запас сырья

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

ь.

 

«и

 

 

«12

 

2

 

 

 

 

«22

 

3

 

 

«31

 

 

«32

 

4

 

 

«41

 

 

«42

Доход

 

 

 

 

 

 

 

В

последней строке у к а з а н

доход,

который

получает

предприятие от реализации единицы к а ж д о г о вида про­ дукции.

Математическая

формулировка задачи

такова:

 

max

ххг

+

с2х2);

 

 

 

(1.9)

a i i x i + а и х 2 < Ьх;

 

 

 

 

азххх

+ а32х2

< b3;

 

 

 

« 4 1 * 1 +

« 4 2 * 2 < bA.

 

 

 

Иными

словами,

требуется

максимизировать

доход

при соблюдении ограничений (1.10). Здесь

х \ и х 2

— ко­

личество продукции

соответственно 1-го и

2-го видов.

Рассмотрим теперь другую задачу . Предположим, что

некоторая

организация ж е л а е т

приобрести

сырье,

кото­

рым располагает предприятие. Спрашивается: по какой

цене предприятие

будет продавать

сырье

организации?

Обозначим через

уи у2,

Уз, УІ цену

к а ж д о г о

вида сырья.

Ясно, что предприятие

не

будет продавать сырье, если

выручка от этой

п р о д а ж и

окажется

меньше дохода, ко-

3 Л. П. Падалко

33


торый оно получает при изготовлении продукции. Этому условию отвечают неравенства:

аііУі + а,2у2 + a3„y3 + ai2y4 > c2 .

В то ж е время о б щ а я стоимость всех запасов сырья, приобретаемых организацией,составит

hUi + bzy2 + b3y3 + ЬЛ yt.

Ясно, что организация стремится минимизировать за­ траты на приобретение сырья. Таким образом, мы при­ ходим к следующей з а д а ч е линейного программирова ­ ния:

min

 

+

b„y„ + b3tj3

+ biiji);

 

 

(1.11)

вцУі

+

а»гУі +

апУз

+ СЦ\УІ

>

сі,

)

(1.12)

 

 

 

 

 

 

 

 

 

 

вііУі

+

« 2 2 0 2 +

аз°Уз

+

a^JJi

>

с».

)

 

Сравнивая

данную

задачу

с предыдущей, мы заме ­

чаем следующие

особенности:

 

 

 

 

1) матрица

 

 

 

 

 

 

 

 

w 1 2

а 2 ї а32

а 4 „/

 

 

 

 

 

 

при переменных у\, у2, уз, УІ получается из матрицы при переменных хи х2

^21

^ 2 °

 

а 3 1

а32

 

с 4 1

а 4 2

 

транспонированием, т. е. заменой

строк столбцами;

2) неравенства в системах ограничений направлены в.

разные

стороны;

 

3)

роль свободных членов

в системе ограничений

(1.12) выполняют коэффициенты

максимизируемой функ-

34


ции (1.9). В то ж е время коэффициенты минимизируемой функции являются свободными членами системы ограни­

чений

(1.10);

 

 

 

 

 

 

 

 

 

 

 

4)

в .первой

задач е

целевая

функция

максимизирует­

ся, во второй — минимизируется.

 

 

 

 

 

Д в е

задачи

линейного

программирования,

имеющие

такие связи м е ж д у собой, называются взаимно

 

двойствен­

ными.

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, если п р я м а я

з а д а ч а

линейного про­

граммирования

имеет такой вид:

 

 

 

 

 

max

+ с2 *2 +

... +

спхп);

 

 

 

 

 

 

% Л + а12х„

+ ... + ах„

< V .

 

 

 

 

 

°ПХ1

 

+ ^22*2 +

••• + а2пХП

<

Ь „ \

 

 

 

 

 

а « Л

+ а„.2х2

+

••• + a m n x n <

Ьт;

 

 

 

 

 

* х > 0 , - v 2 > 0 ,

... ,

л - „ > 0 ,

 

 

 

 

 

 

то двойственная к ней будет следующей:

 

 

 

min [Ьш + b2y2 + ... +

bmijm}\

 

 

 

 

 

ОііУі

+ аиУг

+

••• +

 

атт>с\>

 

 

 

 

 

аігУі

+ а-ігУг +

••• +

атіут

> с3 ;

 

 

 

 

 

ащУі

+ а2пУ2

+

••• + атпут

>

ст.

 

 

 

 

 

Основой теории двойственности является

следующая

теорема: если

одна

из двойственных

задач

 

линейного

программирования

 

 

имеет решение,

то и

другая

задача

также

имеет решение,

при

этом

максимум

одной

задачи

равен

минимуму

 

другой.

 

 

 

 

 

 

 

Эту

теорему можно доказать,

прибегнув

к

множите­

л я м Л а г р а н ж а . С этой

целью составим

функцию

Л а г р а н -

ж а для прямой

задачи:

 

 

 

 

 

 

 

 

 

m

п

 

 

п

 

 

 

 

 

 

 

^ і = 2 с Л + 2 yfti

- 2

Ч*Ъ-

 

 

 

 

 

 

 

i=\

/=1

 

1=1

 

 

 

 

 

 

 

3*

35