Файл: Арутюнян А.Г. Применение математических методов и ЭВМ в народном хозяйстве.pdf

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

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

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

Добавлен: 17.07.2024

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

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

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

. . .

lA.p I_________ (Хя")

 

 

 

 

(X'^J

о

 

A

 

P u C

11 1

множеств следует из конечности количества всевозможных основ (следовательно и базисов) задачи (11.1)—(11.3).

Простыми примерами можно показать, что в случае

произвольных выпуклых функций утверждение

1 неверно.

Отметим, наконец, что на множестве АД,

0 функ­

ция Р(Л, М) является непрерывной функцией,

причем на

этом множестве она порождается хотя бы двумя базисами (базисами, порождающими F{ А, М) на ЛД, и на ЛД,).

3°. Постановка задачи параметрического линей­ ного программирования с линейно зависящим о т па­ р а м е тр а вектором условий. На практике часто могут встретиться задачи ПЛП, в которых один из векторов условий линейно зависит от векторного параметра. Не на­ рушая общности, можем предположить, что от параметра зависит первый вектор условий. Задача формулируется

следующим образо_м.

 

 

 

Для каждого Л £ RK необходимо

построить функцию

F (А ) = max f [А (А)1 =*ягах ^

с х (А ),

(П .11)

°1

°А ,«1

 

 

где множество

определяется условиями

 

232

Ф, (А )* , (Л) + Е

вцлгДЛ)— »,, ,'= ,1 , 2....

т\ (11.12)

x , { \ ) s * 0 ,

j = \, 2.....

я.

(11.13)

Функции а,, (Л) предполагаются линейными, т. е.

К

= a°h +

Oki1/-k. 1 = 1, 2,..., rn.

Множество

О= |А£Дк|С_ ^ 0 |

называем множеством допустимости задачи (11.11)—(11,13). Ясно, что область определения функции F (A )

N = |A£D \тах f\ X (А)] < ool

является подмножеством множества D , т. е. N ( D .

Для решения поставленной задачи в первую очередь необходимо изучить свойства множеств D и N.

Множество D может быть пустым или же состоять из единственной точки. Здесь пока будем рассматривать слу­

чай, когда множество D

содержит хотя

бы две

точки

A '£ D и A "£D . Исследуем

в отдельности

каждый

из сле­

дующих возможных (одновременно взаимно-дополняющих) случаев:

а) лгДА'). х, (Л") = 0 , б) JTj (Л'). х ,( А " )^ 0 .

Как видно,

случай

а) означает,

что хотя бы’

одна из

величин X j(A '),

.*j(A")

равна нулю.

Пусть,_для

ясности

л "](А ')= 0 . Нетрудно убедиться, что вектор X (А ) = Х (А ') будет планом задачи (11.11)—(11.13) в любой точке отрез­ ка, соединяющего точки А' и А". Более того, этот вектор

будет планом задачи в любой точке A

Значит, в этом

случае D = /?к.

 

Рассмотрим, теперь, случай б). В промежуточной точке

Л = аА ' + (1 - с ) А " (O sSasSl)

(11.14)

определим вектор

 

х (л )==р Л ( ло + ( 1 - р) Х ( л "),

(11.15)

233


где величина ^ определяется формулой

(11.16)

оХ\(Л")

Очевидно, величина ja удовлетворяет условиям:

О ^ |х 1; {г/»_о = 0; ti/e=i =■ 1 -

Непосредственной подстановкой можно убедиться, что вектор X (Л) является планом задачи (11.11)—(11.13) в соответствующей точке Л. Значит, множеству D выпукло.

Подставляя выражение (11.15) вектора X (А ) в целевую

функцию, получаем:

 

 

 

 

f [ X ( A ) } ^ i, f [ X { \ /) \ + ( \ - ^ ) f [ X { T V ' ) \ .

(11.17)

Из формулы (11.15)

видно, что

при

Х1(Л,/)и »Х 1(Л/)

£1= 3. Легко показать,

что в этом

случае

функция F (A )

будет вогнутой (выпуклостью вверх) функцией на

отрезке,

соединяющем точки Л' и А".

Действительно, если

Х (Л ') и

X (Л") — оптимальные

планы

задачи (11.11)—(11.13) соот­

ветственно в точках Л'

и Л", то из равенства (11.17) имеем:

/[X (A )l= «f(A ') K I - . ] F ( ? ) .

Но вектор -А"(Л), вообще говоря, может быть и не­ оптимальным планом задачи в точке Л. Отсюда и вытекает неравенство, определяющее вогнутость функции.

При разработке алгоритма решения задачи (11.11)— (11.13) будем пользоваться результатами следующего утверждения.

Утверждение 6. Пусть:

а )

X «£D ;

 

 

 

б) для любого A^A^^G .-» л:1(Л0)^=0:

 

 

в)

функция /(АГ(Л 0)]

неограничена

сверху на

мно­

жестве

о.

 

 

 

Тогда A=f=0 .

 

 

 

Доказательство. Из

предположения

б) вытекает, что

множество D совпадает с

пространством

и любой

век­

тор A'(A°)^G-o является планом задачи на всем пространст­ ве изменения Л.

234


Согласно предположению в), из множества

G- о можно

выделить такую

последовательность

векторов Л 1(Л 0), что

 

 

 

Пт / [ X 1(Л*)| = t -(- оо .

 

 

Так как каждый из векторов ^ ( А 0)

является планом

задачи

в

любой точке

пространства,

т. е.

А', (Л °)^0^

(/= 1 ,

2,...)

для А £ Я К, то

утверждение

верно.

 

Утверждение 7.

Если Л’ £ D \ N l

и Л"£ D \ ,V , то

/VI) (Л £ R K/Л =

a7V +

(1 - а) А ";

0 ^

а ^ 1 ] = 0 .

Доказательство. Согласно предположению утверж­

дения

можно

построить

последовательности

векторов

X ' { A % G r

(/= 1 ,2 ,...) н X l { A " ) £ G r

(/ =

1, 2,...)

такие, что

Пт f \ X l { \ ') } = -| -оо и / i " i / [ A ' ( A " ) ] = + oo .

Пользуясь формулой (11.15), в промежуточной точке (11.14) можно построить соответствующую последователь­ ность векторов Х ' { А )£ G i. Учитывая равенство (11.17) и свойства величины |i, получаем:

И Г Л f \ X ‘ ( Л ) | жав -|- ОО.

А это означает справедливость утверждения.

Вместе с задачей (11.11)—(11.13) рассматриваем сле­ дующую задачу, которую условно назовем урезанной за­

дачей.

 

 

 

 

 

 

 

Найти

максимальное значение линейной

функции

 

/

' W

£

CjXj

 

( i i . l i ' )

при условиях

 

 

) = 2

 

 

 

 

 

 

 

 

v

au Jfj —

1=1,

2.....

т,

(11.12')

1=2

 

 

 

 

 

 

 

JC,X), j

=

2,

3,...,

п.

(11.13)

Как видно, задача (11.11')—(11.13') отличается от за­ дачи (11.11)—(11.13) только тем, что здесь не участвует

1 D \ N — означает дополнение подмножества N до D, т. е. iVUD\<V= D.

235


первый

вектор условий. Задача (11.И ')— (I I .13') является

обычной

задачей линейного программирования.

Если Х ' — (хг, хи) — план задачи (11.1 Г )—(11.13'), то X — (0, л:,,..., х п) — план задачи (11Л 1)—(11.13). Отсюда ясно, что для построения функции F (А) на множестве ее определения можно использовать задачу (I I .И ')—(11.13').

§ 2. Краткое описание алгоритма решения обшей задачи

Добавлением переменных отклонения, систему ограни­ чений (11.2) приведем к каноническому виду и обеспечим неотрицательность правых частей этих ограничений при достаточно больших значениях параметров При этом, если для больших значений ).к все функции остаются неот­ рицательными, то векторы, соответствующие переменным отклонения, уже образуют исходный допустимый базис для решения поставленной задачи методом последовательного улучшения плана. Если же при таких значениях парамет­ ров Хк хотя бы одна из функций ф| (А) отрицательна, то вос­ пользуемся расширенной задачей (/И—задачей), для реше­ ния которой также применяем метод последовательного улучшения плана2. Предполагая достаточно большими зна­ чения параметров ).к и рг, через конечное число шагов по­ лучим оптимальный план расширенной задачи или же об­ наружим неограниченность целевой функции на множестве планов задачи. Ясно, что во втором случае текущий план расширенной задачи будет планом и исходной задачи, так как, до появления признака неограниченности целевой функ­ ции, в текущей симплексной таблице векторы, соответст­ вующие искусственным переменным, заранее исключаются из базиса. В таких случаях можно сразу перейти к опре­ делению области изменения параметров p.j, на которой дан-1

1 Нетрудно заметить, что при сделанных предположениях относи­ тельно функций i]/j (Л) они сохраняют свой знак при достаточно боль­ ших значениях параметров >.к. Аналогичным свойством обладают функ­ ции

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

236


мый опорный план является оптимальным, т. е. определить область Л',, на которой функция F (\ , Ы) порождается од­ ним и тем же базисом. После этого осуществляется пере­

ход

к одной из соседних областей.

_

Если вне +

£)-мерного параллелепипеда Л ^ Л тах,

 

получено

решение расширенной задачи, то осу­

ществляем переход к такой соседней области, на которой в оптимальном базисе участвует меньшее количество ис­ кусственных векторов. Очевидно, этого можно добиваться фиксацией области изменения параметров ц, и изменением параметров /.* {т. е. переход реализуется только в под­ пространстве изменения параметров >.к). При переходе к соседним областям, выбором определенного направления, каждый раз .многопараметрическую задачу приводим к однопараметрической. Решая по выбранному направлению S однопараметрическую задачу, через конечное число ша­ гов метода последовательного уточнения оценок либо избавимся от искусственных векторов, либо обнаружим, что вдоль выбранного луча исходная задача неразрешима. В последнем случае выбирается другое направление, по которому можно исключить из базиса хотя бы один искусст­ венный вектор.

Приведение многопараметрической задачи к однопа­ раметрической оправдывается тем, что свойство выпуклости функции многих переменных эквивалентно свойству выпук­ лости такой функции одной переменной, которая получает­ ся от данной функции многих переменных фиксированием произвольной точки и произвольного направления, т. е. имеет место следующее утверждение (см. [34], § 2.4, тео­ рема 1):

функция ф(Л) выпукла тогда и только тогда, когда функция ф' (<х) = ф (A w-f- a S) выпукла для произвольных А 0 и S.

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

Может оказаться, что в процессе решения задачи (11.1)—(11.3) вне вышеуказанного параллелепипеда, на

237