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

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

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

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

Добавлен: 17.07.2024

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

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

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

Определяются значения

р.:

 

 

 

jj-—j—I ^

0,

]1 ^

1;

 

V + З ^ О ,

— .

 

 

 

Г

4

п

3

 

x t

меняет знак, следова­

о точке

|j. = — переменная

тельно Я4 должен выйти из базиса, а в базис вводится — Р2. Текущая таблица преобразуется методом последователь­ ного уточнения оценок.

Базис

X,

Ха

СВ. Ч-

1\

Р2

Р ,

Р 4

 

 

ft

Номер

 

 

 

 

 

 

 

 

 

 

 

 

итерации

/>,

 

4

 

 

 

 

 

2

 

 

1

 

Т

~ 5

 

 

 

 

5

 

_

5

 

 

 

 

 

 

 

 

 

1

П_

 

 

1

1

 

3

12

 

9

2(1)

5

 

 

 

~5

~

5

 

_

5

 

 

 

 

 

-

1

5

— ]

0

0

1

 

1

—6

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

По направлению S решается однопараметрическая за­ дача, фиксируя отрезки и базисы, порождающие функцию

4 > W = f (A(0,+ l ‘ S) на этих отрезках.

В результате получается:

A -" = A « + ^ = ( i , i ) 4 - f ( - 1 . - Ч = ( т ' Y >

Из полученной симплексной таблицы определяется строка, базисные переменные которой в точке ji,nax меняют знак (строка Яа).

Из первых трех столбцов строки Я2 составляется урав­ нение

и находятся промежуточные значения К


Если л1==0, то />2 = - 3 -, а если Х2 = 0, то >., = 3.

В результате получается промежуточная точка Л(2\«= (3,0).

Теперь на этой прямой выбирается новое направление

Sv Так как Ли = (3,0), а Л(,,= ( - L - 1 ) , т0

^ - [ ( з - т Х о Ч ) ] - ^ - - ^

По выбранному направлению снова перейдем к'-одно- параметрической задаче и определим коэффициенты при

т - ^ + ( - т ) ( - - т ) - т .

( Х ) ^ + ( - ^ ) ( - т ) = о ,

4 l + 5( ~ t )= 4 -

Свободные члены будут:

_3_ _3_____L — J_

5 * 4

5 ~ 4 ’

270

Имеющаяся таблица примет вид:

 

Рис.

14.3.

 

З д е с ь

т. е.

----- —, следовательно ис­

ключению подлежит

вектор

Я,, а включению

в базис —

вектор Я4.

 

 

 

Найдем точку А<3), соответствующую р = —

—:

л ^ л » + , 5 = ( Л

- л . ) -

Преобразуем симплексную таблицу:

Базис

X,

h

СВ. ч.

Яд

Я,

Р}

Я4

 

 

Номер

 

 

 

 

 

 

 

 

 

 

итерации

Pi

1

2

— 1

5

 

____L

 

15

5

 

2

2

 

 

8

8

 

 

 

 

 

2

 

 

 

1

 

 

3

 

1

 

 

 

4(3)

_

'2

 

 

2

 

~

 

9

3

 

 

 

 

 

3

з

0

5

0

3

0

27

9

 

Для выбора нового направления, определяется проме­

жуточная

точка Л":

 

 

 

 

 

 

- ~ K + 2 h - i = 0-,

 

 

 

если ^ =

0, то х2 =

- р

а если хг= 0. то

х, =

— 2.

 

Отсюда Л "= (0 ,

-i-j.;; Однако Л,1П= ( — - l

-L'j.

Следа-

вательно: S |2> = [^0j

(

~ ) ) , ( | - ; |

) ] r

( f

■ } ) .

272


Новые значения коэффициентов при [* будут соответст­ венно равны:

------!----- - + 2 •— = 0 ,

1

. „ - - Н - Ц - Т — А . А + 3 . Л „ ± ,

а свободных членов:

 

 

 

- Н -

т )

+ т -

°

-

 

 

 

Таким образом,

текущая

таблица

примет вид:

 

Базис

X,

X,

СВ.

ч.

Р,

Р 2

Р3

Р 4

н-

р

Номер

итерации

Р4

1

2

— 1

5

0

1

1

0

0

 

 

 

 

 

2

_

т

 

 

 

 

Л

1

— 1

0

а

1

1

0

1

0

5(4)

~~2

~ 2

 

~

 

 

 

 

~2

~

 

 

-

 

,3

0

5

0

 

0

 

0

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По выбранному направлению 5 (2), реализуя один шаг метода последовательного уточнения оценок (при р- > 0 базисная переменная х 2 меняет знак), получаем следую­ щую таблицу:

Базис

X,

X,

СВ- ч.

Pj

Р г

Р 3

Р 4 и

?

Номер

итерации

Рл

 

11

—1

0

5

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

Pj

Т

 

0

1

2

1

0

0

 

3

3

~~з~

6(5)

 

2

4

0

0

5

2

 

0

 

 

3

3

3

3

 

 

 

 

 

 

 

 

т

18—644


В последней таблице признак

оптимальности

не накла­

дывает нового ограничения на fj-

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

име­

лось в самой постановке однопараметрической задачи),

значит на полуоси

однопараметрическая задача ре­

шена.

 

Возвращаясь к двумерному пространству (обычной плоскости), можно убедиться, что уже завершен процесс решения исходной задачи (см. рис. 14.5). Для части плос­ кости AOjC оптимальным базисом является (Р2, Р4), для

СОгВ (Plt Р2), для АОгВ (выше прямой AB) - (PV Р4).

Иными словами, в подобласти А 0 4С целевая функция по­ рождается одним и тем же базисом (Р2, Р4), в подобласти COjS — базисом {Pv Р2), в подобласти АОгВ (выше прямой АВ) — базисом (Pj, Р4).

§ 5. Решение п р о с т е й ш е й однопараметрической задачи

Однопараметрическая простейшая задача в случае изме­ нения численного параметра на конечном интервале была детально изучена С. Гассом и Т. Саати (см. 120]) при пред­ положении невырожденности задачи. Она является частным видом задачи (14.1)— (14.3), когда

Как видно из предыдущих параграфов, при решении многопараметрнческих задач мы неоднократно обращаемся к однопараметрическим. В частности, при решении задачи (14.1)—(14.3) по данному направлению решается простей­

шая однопараметрическая задача, которая

формулируется

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

 

 

Для каждого

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

F (X) = max f [ X (X)) = т а х £ с ^ ( Х ) ,

(15.1)

 

I=1

 

где область Ж\ определяется условиями:

1]

(X) = 6 i+

X6 'i, (=el,

2,...,

т,

(15.2 )

i-i

 

 

 

 

 

 

 

 

*ДХ)гз=0,

/ =

1,

2.....

га.

 

(15.3)

Здесь сь а.ц,Ьь Ь\ (г=1,

2,...,

/га; /'= 1 ,

2

га)— заданные

действительные числа.

 

 

 

 

 

 

В [20]

приведен алгоритм

построения

функции F (X)

на конечном интервале. Но так как по этому алгоритму мы двигаемся по оси X только в одном направлении, то сле­ дует либо запомнить матрицу, соответствующую начальной точке Х0 (чтобы в дальнейшем построить f(X ) на допол­ нении уже рассмотренной полуоси), либо после получения одной границы области N двигаться в обратном направ­ лении по рассмотренным подобластям Nт до получения другой границы. Следовательно, решение поставленной за­ дачи на ЭВМ в первом случае приводит к потерям объема памяти, а во втором случае — машинного времени. Из до­ казанных в предыдущих параграфах утверждений сле­ дует, что:

а) области М, являются конечными отрезками, кроме, может быть, двух крайних, которые могут быть лучами;

б) если в точке Хт получены все базисы, порождающие F(X,), и ни один из этих базисов не остается допустимым левее (правее) точки Хт, то X- является левой (правой) гра­ ницей области Nil

в) на каждом N. функция F(X) является линейной, а на общих границах двух соседних Мт—непрерывной функцией;

г) функция F(X) вогнута.

Для решения задачи (15.1) — (15.3) пользуемся соот­ ветствующей расширенной задачей, обеспечивая заранее неотрицательность правых частей условий (15.2) для доста-

275


точно больших значений параметра ).{Ь'|3 »0 , / — 1,

2 ,..., т

и, если b'i = Q для данного /,

то

0 ).

 

Ясно, что для достаточно

больших значений X

основа,

состоящая из искусственных векторов, будет базисом рас­ ширенной задачи.

Для расширенной задачи заполняется исходная сим­ плексная таблица, отличающаяся от симплексной таблицы обычной числовой расширенной задачи тем, что здесь для значений базисных переменных отводятся два столбца (вместо одного): в первом столбце заполняются коэффи­ циенты Ь\ при X, а во втором -свободные члены Ь>.

К заполненной симплексной таблице применяется алго­ ритм метода последовательного улучшения плана и через конечное число шагов либо обнаруживается неограничен­ ность целевой функции, либо обеспечивается признак опти­ мальности (неотрицательность оценок внебазисных векто­ ров). В первом случае заключаем, что jV = 0 . Во втором случае имеется базис, являющийся оптимальным базисом расширенной задачи на некотором луче (X ',,+ ©о). Как видно, нахождение точки X't эквивалентно (по объему вы­ числений) решению одной числовой задачи линейного про­ граммирования. Если в полученном оптимальном базисе не

участвует

ни

один

из искусственных векторов с положи­

тельным

значением

базисной неизвестной, то >•',=«Xj и

Nl = (t.v -\-oo),

т. е.

функция f {).) уже определена на этом

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

В [20] приведены формулы преобразования симплексной таблицы и определения отрезка, на котором функция /■’ (/.) порождается одним н тем же баз'исом. Очевидно, эта ме­ тодика приемлема и для построения соответствующей функции расширенной задачи. Через конечное число пере­ ходов строятся последовательность соседних отрезков и выражения целевой функции. Если при таких переходах попадаем на отрезок, где оптимизирующий базис уже не

•содержит ни одного искусственного вектора, то это озн&г

276