Файл: Арутюнян А.Г. Применение математических методов и ЭВМ в народном хозяйстве.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