Файл: Арутюнян А.Г. Применение математических методов и ЭВМ в народном хозяйстве.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 17.07.2024
Просмотров: 176
Скачиваний: 0
чает, что найден отрезок Nv В этом случае при дальней ших переходах отмечаются границы последующих отрез ков Nr и соответствующие выражения целевой функции
|
|
F (}) = } |
|
. |
).£N.г, |
т = 1, 2,... |
|
|
Через |
конечное |
число |
переходов |
получится либо луч |
||||
( — оо, |
Xmi„), на котором |
функция f(>.) |
порождается |
одним |
||||
и тем |
же |
базисом, |
либо |
точка |
левее которой |
область |
N не имеет ни одной точки.
Возможно, что указанный выше случай встретится при наличии в базисе хотя бы одного искусственного вектора. Очевидно, в этом случае N = '0 .
Таким образом, применением метода последовательного улучшения плана, а затем (при движении в сторону умень шения параметра) — метода последовательного уточнения оценок, а также пользуясь формулами определения границ отрезков Nx (см. [18]), за конечное число шагов указанной процедуры строится функция на всей области ее опре деления.
Блок-схема решения задачи (15.1)— (15.3) приведена на стр. 278.
§ 6. Решение задачи параметрического программирования с вектором условий, линейно зависящим от параметра
Предположим, для ясности, что урезанная задача (11.1Г )—(11.13') разрешима. Применяя к ней метод после довательного улучшения плана, через конечное число ша гов получим оптимальный опорный план. Добавлением к этому плану переменной х , с нулевым значением получим опорный план задачи (11.11)—(11.13). Последний будет оптимальным опорным планом для некоторого множества изменения параметра Л. Для определения границ такого множества вектор
P ,(A ) = Q »,+ V k^l
(где Qk,= ( a ku, ak21,..., akmi), £ = 1, 2 ..... К) представляется
277
Блок-схема решения простейшей параметрической задачи
Начало
278
по векторам имеющегося базиса. Пусть лг0и -|- £
_ к=1
коэффициент при базисном векторе Р . (занимающем /-ую
позицию в базисе) в представлении вектора РДЛ) по век торам базиса. Оценка
ДДЛ) = |
sTj(Л) — сг= |
V с1{ |
V jckixXk |
|
|
к |
|
|
— С, ~ |
1 И - \ |
X km+ 1, 1 ).к |
|
|
к=1 |
|
вектора РДЛ) |
является линейной функцией от Л. Отсюда |
||
ясно, что условие г,(Л ) — |
определяет некоторое по |
лупространство. Так как на этом полупространстве значе ние целевой функции не зависит от Л, то оно постоянно.
Для построения функции |
Р (А ) на остальной части ее |
||||
определения рассматривается |
полупространство |
|
|||
|
к |
|
|
|
|
* ’ш+1. 1+ |
Е Л |
+1.1 |
К < 0 . |
|
(16.1) |
|
k. I |
|
|
|
|
Может оказаться, что на этом полупространстве каж |
|||||
дое из выражений |
к |
|
|
|
|
__ |
|
|
|
|
|
x il(A )= x °ii-f -^ |
|
, t — 1 , |
2 ,..., т |
(16.2) |
|
принимает только |
неположительные |
значения. |
Ясно, что |
||
в этом случае указанное |
полупространство не |
содержит |
ни одной точки области N.
Если же в полупространстве, определяемом условием (16.1) , хотя бы одно из выражений (16.2) остается поло жительным (а это означает, что пересечение двух полу пространств
Л - Н .1+ £ Л +1.1'к < о и |
^ 1и>-к> о |
к«1 |
к-1 |
является непустым множеством), то фиксируется некото рая точка на граничной гиперплоскости полупространства (16.1) , выбирается определенное направление и на этом направлении многопараметрическая задача приводится к однопараметрической. В качестве такого направления мож-
279
но выбрать, например, |
вектор внутренней нормали |
к гра |
ничной гиперплоскости полупространства (16.1). |
|
|
Таким образом, ставится задача построения функции |
||
'f(a) = F (A °-f-aS ) на |
полуоси 0 ^ я < ( + оо. Ясно, |
что на |
полуоси — o e < a s s :0 |
имеет место |
|
-р (» ) — ^ ( А “).
В результате решения однопараметрической задачи строится последовательность отрезков, на каждом из ко торых функция <р(я) порождается одним и тем же базисом (вернее, одной и той же группой базисов). Получив гра ничную точку nm„ , выбирается новое направление и по этому направлению снова решается однопараметрическая задача, не рассматривая отрезки, соответствующие базисы которых уже были определены. Так как каждый раз вы бором направления обращаемся к однопараметрической за даче, то следует детально изучить однопараметрическую задачу. На указанной полуоси вектор Р ,(Л ) будет пред ставлен в виде
|
|
|
Р г (-4 = Q"i + |
»<?,- |
|
|
|
Коэффициентами при базисных векторах в разложении |
|||||||
этого |
вектора будут |
|
|
|
|
||
|
х\р = |
x 'n -\-zx'\v / = 1, 2 ,..., т, |
|
(16.3) |
|||
|
1 |
>1 |
Т |
+ К я |
V о |
|
(16.4) |
|
|
|
|
к—1 |
|
|
|
представляет собой значение базисной переменной |
в точке |
||||||
Л = Л |
(являющейся началом выбранного луча aS). |
Коэф |
|||||
фициенты при параметре а в выражении (16.3) |
опреде |
||||||
ляются формулой |
|
к |
|
|
|
|
|
|
х"-п = |
-£kii% / = |
1, 2 ..... т. |
|
|
||
|
|
|
|
||||
|
|
|
к = 1 |
|
|
|
|
Ясно, что оценка вектора РДЛ) тоже будет представ |
|||||||
лена |
в виде (16.3), |
т. |
е. вместо |
полупространства |
(16.1) |
||
имеется полуось |
0 < а < -| -о о , на которой |
|
|
||||
т. е. |
|
Х'т + 1,1 “Ь г -X''m+l, 1 < О, |
|
|
|||
|
|
|
_ |
к |
|
|
|
|
Х и + 1.1 =а Л ш + 1. 1( Л в) « - 0 , |
Х па + х .\ ** V X km+l,l Sy. |
280
В зависимости от значений величин х 'и и х "н рассмат
риваются следующие возможности. Если |
= 0 и *"^*£0 |
для i = l , 2 ,..., т, то на полуоси (0 , -f-oo) |
функция 'р(я) не |
определена, так как целевая функция задачи неограничена сверху на множестве планов задачи. А если jc'ix — 0 для 1 = 1, 2 ,..., т, но хотя бы для одного индекса г, то имеется возможность исключить из базиса -некоторый век тор, вводя вместо него в базис вектор Р1 (Л), и получить таким образом основу, являющуюся базисом на некотором отрезке (0 , а'|.
Аналогичный результат получается и в случае, когда имеется хотя бы одна положительная величина х\г. На этом интервале рассматриваются отношения
■*11(<0
но тем строкам, для которых -*ц(*)1>0. Сравнивая эти отношения, отрезок [0 , а'] можно разделить на несколько кусков, на каждом из которых одно из них принимает ми нимальное значение. Из полученных отрезков выбирается тот, левой границей которого является точка а = 0. На этом отрезке осуществляется преобразование симплексной таблицы. Полученные на последующих шагах симплексные таблицы уже целиком будут зависеть от параметра я. Каж дый элемент таблицы выражается в виде дробно-линейной функции от я:
< » ( • ) - “ *Щ-**Ч, i = l , 2..... |
т + 1; / = 1, 2..... |
п, |
(16.5) |
ау + г
где величина « у + 2 представляет собой значение опре делителя матрицы, состоящей из базисных векторов — столбцов.
Так как
х 1о(0 ) = дг>‘° 5=0, i = l , 2 ,..., т,
то область базисности данной основы определяется при помощи знаков величин у и х'-10. В связи с этим в отдель ности рассматриваются следующие случаи:
я) у = 0 ; б) у ф 0.
251
В случае а) выполняется неравенство
где отношения в правой части неравенства берутся по всем
индексам /, для которых |
•>мо<С.О (т- е ' |
и xS>o имеют |
||||
противоположные знаки). |
|
|
|
|||
В случае б) |
необходимо анализировать знак величины г. |
|||||
При 2 > 0 |
(а |
также |
при |
z < 0 |
в случае |
а)) граница <х' |
определяется |
системой |
неравенств |
(16.6). |
|
Если zy^>0, то правая граница области базисности дан ной основы определяется системой неравенств (16.6). Если
же zy |
0 , |
то для определения |
правой границы области |
базисности |
к системе (16.6) добавляется неравенство |
||
|
|
* '< 7 - |
<16Г> |
Совершенно аналогичными рассуждениями можно опре делить верхнюю границу области оптимальности данной основы, т. е.
a ' s s ---- J I fJJ...(Где |
х 2т 11, j < 0), |
(16.8) |
|||
|
x m+l. J |
|
|
|
|
если y = 0 |
или y z > 0 . А |
в случае y z < 0 |
к системе |
(16.8) |
|
добавляется неравенство вида (16.7). |
|
|
|||
Ясно, |
что величина |
а, = |
ягш(а', а") |
представляет со |
бой верхнюю границу области базисности и оптимальности данной основы. В зависимости от взаимного расположения областей базисности и оптимальности данной основы для перехода к соседней справа области применяется либо ме тод последовательного улучшения плана (при а,=а"), либо метод последовательного уточнения оценок (при е^гввя').
На последующих областях изложенный анализ коэф фициентов выражения (16.5) повторяется со следующим
изменением: |
вместо величин г и |
х 2^ рассматриваются со |
|
ответственно |
величины |
—|—atjc'jj |
и z-f-a ty, где а,— ле |
вая граница изучаемой области (т. е. верхняя граница предыдущей области).
Теперь проанализируем возможности перехода от верх ней границы at.
282