Файл: Баженов, Ю. М. Перспективы применения математических методов в технологии сборного железобетона.pdf

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

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

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

Добавлен: 29.10.2024

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

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

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

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

На изготовление арматурных каркасов может быть затрачено 11 чел.-дней. Соответственно за это время может быть произведено 200 м3 бетона, а па изготовление конструкций может быть затрачено 210 чел.-дней (включая вспомогательные работы), для прогрева из­ делий может быть использовано 22 г пара.

Возможности

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

определяют

ограничения

задачи.

Обозначим

количество изготовляемых

конструкций соответственно

*2, *з> х4. Примем, что эффективность всех конструкций одина­

ковая, т. е.

Сі =

С2=

Сз =

с4=

 

1.

 

 

 

 

 

 

Математическая формулировка задачи:

 

 

 

 

 

 

*і +

+

*э +

*4

макс.

 

 

 

(V.5)

 

 

 

 

 

 

х > 0

 

 

 

 

(Ѵ.6)

 

0 . 25*! +

0,10л:а +

0,09*3 +

0,15*,, <

11;

 

 

 

1, 6а'і

+

1 ,5л-,

+

1,6хэ

+

1.9*і

< 200;

 

 

 

2,5*і

+ 1 .9*2

+ 3 * з

+

8

<210;

 

(Ѵ-7>

 

0,28*1 +

0,24*0 +

0,2*3

-j

0 ,2.Vj <( 22.

 

 

Возможным планом задачи является вектор

* = (* |.

*2,

*з, *і),

удовлетворяющий

условиям

(Ѵ.6)

и

(Ѵ.7).

Если же

этот

вектор

удовлетворяет еще и условие (Ѵ.5), то он будет оптимальным пла­ ном задачи.

Для нахождения оптимального плана симплексным методом си­ стема неравенств (Ѵ.7) преобразуется в систему линейных уравне­ ний путем введения дополнительных неизвестных. Обозначим их *5, *6- *7. *8 и назовем дополнительными переменными. Каждая допол­

нительная переменная входит в соответствующее уравнение с ко­ эффициентом, равным единице, а в уравнение целевой функции — с нулевым коэффициентом.

Запишем полученные уравнения, поменяв местами правые и ле­ вые части и для удобства дальнейших вычислений умножив каждый

элемент уравнений на 100:

 

 

 

 

1100 =

25*і +

19*2 + 9*з + 15*і + *5.

 

 

20 000 =

160*і +

150*2 +

160*з +

190*і +

*G;

(V .8)

21 000 =

250*і +

190*і + 300*з +

800*4 +

*7;

 

2200 =

28*1 +

24*2 +

20*3 +

20*4 +

*8;

 

 

L =

*і + *з +

*з -}- *4 = макс.

(V.9)

Решение задачи состоит из следующих этапов: составление исходной симплексной таблицы; проверка плана на оптимальность;

последовательное улучшение плана путем пересчета симп­ лексной таблицы.

На приведенной ниже схеме показаны основные элементы симплексной таблицы.

В верхней части таблицы записываются все неизвестные, вхо-

166


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

В центре симплексной таблицы располагается матрица коэффи­ циентов при неизвестных, в которой число столбцов равно числу всех

неизвестных,

а число строк — числу

уравнений

ограничений. Каж­

дый столбец

матрицы представляет

собой вектор-столбец х, и в

дальнейшем будет называться просто

вектор Xj.

Слева располагает­

ся основная часть матрицы, справа — единичная матрица, состав­ ленная коэффициентами при дополнительных переменных.

Внизу таблицы располагается целевая строка, служащая для определения оптимальности плана на каждом этапе вычислений. Элементы этой строки показывают для каждого вектора Xj, насколь­ ко увеличится значение функционала L цели при включении в план

данного вектора.

Столбец Я Л>— содержит неизвестные, входящие в план на дан­ ной итерации к. Число компонентов вектора Pw должно быть давно числу уравнений ограничений т. Входящие в план векторы должны

быть линейно

независимы.

 

 

 

 

Схема составления симплексной таблицы

 

 

 

 

Коэффици­

 

 

 

 

 

енты эффек­

 

 

 

 

 

тивности

 

 

 

 

 

Шапка мат­

 

 

<7

p W

*0

рицы (номера

а

ß

неизвестных

 

 

 

Xj )

 

 

Коэффици­

Неизве­

Итоговый

Матрица

у *пКЛ

Коэффици­

A0f a if

енты эффек­

стные,

столбец

коэффициен­

 

енты для

тивности

входя­

 

тов при

 

пересчета

при неиз­

щие

 

неизвестных

 

элементов

вестных,

в план

 

 

 

строк

входящих

 

 

 

 

матрицы

в план

 

 

 

 

 

 

 

L w

Целевая

 

 

 

 

 

строка

 

 

В исходной таблице в столбец PW записывают дополнительные

неизвестные. Эти неизвестные независимы, они обозначают фиктив­ ную продукцию, другими словами, решение по симплексному спосо­

бу начинается с момента, когда ничего не производится, т. е. lS°}= 0.

Затем с помощью последовательных этапов решения (итераций) оп­ ределяют нужный объем производства каждого продукта при за­ данной целевой функции и положенных ограничениях.

В столбце С) записываются коэффициенты при неизвестных, вхо­

167


дящих в план. В исходной таблице здесь записываются коэффициен­ ты при дополнительных неизвестных, т. е. 0.

В столбец Л'о записываются свободные члены уравнений ограни­ чении. В конце процесса вычислении в столбце Я(Л) будут записаны неизвестные оптимального плана, в столбце х0— значение этих пе­

ременных, т. е. итог решения, а в /+ > — максимальное (минималь­ ное) значение функции цели.

Последние два столбца таблицы (а и ß) имеют вспомогатель­ ное значение. Заполнение их удобно для выполнения проверки хода

решения.

 

 

 

 

 

 

 

 

 

 

 

 

 

Составление исходной симплексной таблицы производится в со­

ответствии с системой уравнении (V.8).

В нашем примере число

столбцов

матрицы коэффициентов

равно

8,

число

строк равно 4

(табл. V.2).

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

V.2. Исходная симплексная таблица

 

 

 

 

 

1

1

1

1

0

0

0

0

а

 

С° 1 р (°> |

 

Л'о

Х і

Л*» ■ѵ3

Л'і

*5

* 0

x-,

А*а

Р

 

 

 

0

Хц

1

1 0 0

25

1 0

9

15

1

0

0

0

44

0

Л'О 2 0

0 0 0

160

150

160

190

0

1

0

0

125

6,4

0

*7

2 1

0 0 0

250

190

300

800

0

0

1

0

ПО

1 0

0

ха

2 2 0 0

28 24

2 0

2 0

0

0

0

1

78

1 , 1 2

 

 

 

0

1

— 1

1

1

0

0

0

0

 

—0,04

Значение L<°> и элементов целевой строки вычисляется следую­ щим образом:

 

 

т

 

 

L<°>=

S с?*,0;

(V.IO)

 

 

і=і

 

 

 

т

 

(V .ll)

 

= С1 =

s

ci aa ~~ ci-

 

 

Отсюда (0 • 27+0 • 166+0 ■232+0-28) — (—1) = 1-

 

В нашем

примере значение

 

L<°> = 0 -1100+0-20 000+0• 21 000+

+ 0 -2200= 0,

а значение элемента

столбца х {:

(0-25+0-160+0-250+

+ 0 -28)—1 = —1.

168


На основе сказанного можно установить правила, позволяющие формально производить заполнение симплексной таблицы:

а)

в

столбец Я<°> записываются дополнительные переменные;

б)

в

столбец Хо заносятся значения свободных членов уравне­

ний ограничений; в) элементы целевой строки в столбцах основных переменных

равны коэффициентам эффективности при этих переменных, взятых

с обратным

знаком;

 

 

г) элементы целевой строки в столбцах дополнительных пере­

менных

равны нулю;

нулю.

 

д)

значение L равно

 

Анализ плана на оптимальность производится путем просмотра

целевой

строки. Если все элементы этой строки

г —с3- ^ 0 для ос­

новных

векторов, то план, записанный в столбце

Рік\ является оп­

тимальный

и значение

линейной формы этого

плана равно

В других случаях необходимо провести улучшение плана, в частно­ сти после составления исходной таблицы.

Улучшение плана производится введением в план (столбец Я(*>) вектора той неизвестной, которая может больше, чем другие, увели­ чить (уменьшить) значение целевой функции задачи. А так как в плане может быть лишь т неизвестных, то из плана нужно ис­

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

Алгоритм симплексного метода определяет формальные прави­ ла выбора вектора, который необходимо ввести в план, н вектора,

исключаемого из

плана,

и правила

связанных

с этим

расчетов.

1. Определяем

вектор,

вводимый

в план. Для

этого

среди эле­

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

В нашем примере можно вводить любой из векторов *і, х2, х3 или х4, так как коэффициенты у них одинаковые. Введем вектор Х \ . Столбец х х— направляющий. Этот столбец в табл. Ѵ.2 обведен жир­

ной линией.

2. Определяем вектор, подлежащий исключению нз плана:

Для этого находим для

каждого

положительного элемента

столбца к X ik> 0 отношение

т. е.

отношение свободных чле-

xik

нов уравнений к соответствующим коэффициентам направляющего столбца. Значения этих отношений записываем в столбец а табл. Ѵ.2.

Минимальное из этих отношений определяет исключаемый нз плана вектор. Выбранная строка называется направляющей строкой и обо­

значается

индексом

/.

 

 

 

 

Для нашего примера находим:

 

 

 

и =

Ц 0 ? _

44.

*20

20 000

*з0

21 000

 

А'.п

25

*21

160

’ *31

190

01

*4о __ 2200 = 78

*4і 28

Выбираем минимальное из этих отношений, т. е. 44.

169