Файл: Шепелев, И. Г. Математические методы планирования и управления в строительстве конспект лекций.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.10.2024
Просмотров: 61
Скачиваний: 0
при округлении пользоваться правилами арифметики, то мож но получить целочисленные значения Х\ = 4 и х2 = 3. Точка 4-— —3 также лежит вне зоны допустимых значений переменных. Це лочисленные оптимальные значения переменных надо искать среди всех возможных целочисленных значений, лежащих внут ри многоугольника, построенного ограничениями. На рис.. 13 такие точки выделены черным цветом. Если линию X i х2 передвигать параллельно самой себе, то внутри зоны возмож ных значений максимум этой линии пройдет через целочислен ные точки Х\ — 2; х 2— 3 и Х\ = 3; х 2 — 2, х { + х 2= 5 (птах).
Это и будет целочисленным решением задачи. К таким же ре зультатам можно прийти, решая задачу в соответствии с пер вым алгоритмом Гомори.
§ 6.5. Метод ветвей и границ
Метод относится к классу комбинаторных методов целочис ленного программирования и характерен тем, что из множества возможных планов рассматриваются лишь некоторые. Часть возможных планов отбрасывается, как заведомо не представ ляющая интереса. Однако это не эвристическое, интуитивное отбрасывание, а целесообразное, определяемое строго матема тически.
Смысл метода заключается в том, что определяется нижняя
граница |
целевой функции |
и все множество |
возможных пла |
нов G разбивается на ряд |
подмножеств G'. Это разделение про |
||
изводится |
в зависимости |
от того, имеется ли |
по какому-либо |
переменному нецелочисленное решение в базисном оптимальном плане. При дальнейшем разбивании получается дерево ветвле ний.
Схема ветвлений показана на рис. 14. На каждом шаге вет вления любым известным способом линейного программирова ния определяется оптимальный план. Если оптимальный план отвечает условию целочисленности, то ветвление прекращается. Ветвление прекращается и по той ветви, где задача оказывает ся неразрешимой. Ветвление продолжается по нецелочислен ным в оптимальных планах переменным, имеющим ограничения по целочисленности до тех пор, пока не будет получено цело численное решение.
Одним из первых среди комбинаторных методов разработан алгоритм Лэнд и Дойг, смысл которого заключается в следую
щем. Допустим, поставлена |
частично |
целочисленная |
задача |
z = f (л) = |
П |
min |
(6.5.1) |
^ Cj X j |
|||
|
j |
|
|
88
п р и о г р а н и ч е н и я х |
|
|
|
|
||
|
|
2 |
Ьй |
I = |
1, 2, . . . , m |
(6.5.2) |
|
|
j |
|
|
|
|
|
. |
0 < jc lj< d j; |
j = |
1,2, .... n |
(6.5.3) |
|
|
x L— целые |
числа; |
/ == 1 , 2 , . . . , rij |
(6 .5 .4 ) |
||
ni ^ |
n — условие частичной целочисленное™. |
|
||||
Порядок решения (алгоритм): |
|
задан |
||||
1. |
Решается |
задача линейного программирования, |
||||
ная моделью (6.5.1—6.5.3). |
|
|
|
|
||
Если полученный при этом оптимальный план Х° удовлетво |
||||||
ряет условию целочисленности |
(6.5.4), то задача решена. Если |
оптимальный план не удовлетворяет условию целочисленности (6.5.4), то значение целевой функции Z0 дает нижнюю границу
для искомого оптимума Z*. |
|
|
||
2. Пусть |
некоторая |
переменная -*1о(1 < i0 < щ) не получила |
||
в оптимальном плане Х° целочисленного значения, |
тогда для |
|||
выполнения |
условия |
целочисленности |
необходимо, |
чтобы в |
дальнейшем эта переменная была |
|
|
||
|
|
\х1о] или лу>.] А-1о[; |
(6.5.5) |
|
заметим, что ] x-h [ = [-xio] + 1. Решают |
две задачи |
линейного |
программирования, в которые вводится дополнительное ограни чение (6.5.5). Это и есть ветвление.
3. Если одно из полученных оптимальных решений X' отве чает условию целочисленности, то это будет решением задачи. Если же одно из полученных решений будет пустым, т. е. зада ча неразрешима, то этому решению приписывается оценка + с/х Если одно из решений оптимально, но не отвечает условию це лочисленности, то это решение Z' будет нижней границей для последующих решений, а задача разбивается на две задачи при дополнительных ограничениях
(6.5.6)
(6.5.7)
" Таким образом, решение задачи сводится к последователь ному решению задач линейного программирования при все бо лее сужающихся условиях функционирования переменных х.
Кроме этого метода для решения комбинаторных задач цело численного программирования применяются алгоритм Литтла, Мурти, Суини и Кэрел, аддитивный алгоритм Баллаша, и неко торые другие методы. Большинство «з них рассмотрено в моно-
|
|
|
|
ЛИТЕРАТУРА К ГЛАВАМ V—VI |
|||||
1. |
А р и с Р. Дискретное динамическое программирование. М., «Мир», 1969. |
||||||||
2. |
Б е к Н. Н., Г о л е н к о Д. И. |
Статистические |
методы оптимизации в |
||||||
экономических исследованиях. М., «Статистика», 1971. |
|
||||||||
3. |
Б е л л м а н Р., Д р е й ф у с С . |
Прикладные задачи динамического про |
|||||||
граммирования. М., «Наука», 1965. |
программирование. М., «Экономика», |
||||||||
4. |
Б и р м а н |
И. |
Я. |
Оптимальное |
|||||
1968. |
В е и т ц е л ь |
Е. |
С. |
Элементы |
динамического |
программирования, М., |
|||
5. |
|||||||||
«Наука», 1964. |
|
|
Л. |
В., Г о р с т к о А. |
Б. Математическое оптимальное |
||||
6. |
К а н т о р о в и ч |
||||||||
программирование в экономике, М., «Знание», |
1968. |
Дискретное программиро |
|||||||
7. |
К о р б у т А. А., |
Ф и и к е л ь ш т е й н Ю. Ю.. |
|||||||
вание, М., «Наука», 1969. |
|
|
|
«Прогресс», 1967. |
|||||
8. |
Л а н г е |
О. П. Оптимальные решения. М-> |
|||||||
9. |
М а й м и н а с |
Е. |
3. |
Процессы |
планирования в экономике. Изд. 2. М., |
||||
«Экономика», 1971. |
Методы поиска экстремума. М., |
«Наука», 1967. |
|||||||
10. У а й л д |
Д. |
||||||||
11. С ы р ц о в а |
Е. |
Д. |
Математические методы в планировании и управ |
||||||
лении строительным производством. М., «Высшая школа», 1972. |
|||||||||
12. Ш е п е л е в |
И. |
Г. |
О правомерности применения корреляционных мо |
делей в задачах математического программирования.— В сб. «Техника и техно логия месторождении полезных ископаемых. Вып. VII, М., «Недра», 1969.
4
Г Л А В А V I I
МЕТОД СТАТИСТИЧЕСКИХ ИСПЫТАНИЙ (МОНТЕ-КАРЛО)
§ 7.1. Статистическое моделирование
Метод статистических испытаний (Монте-Карло) представ ляет собой совокупность приемов и методов, позволяющих ими тировать (воспроизводить) случайные процессы и приближенно решать любые задачи на экстремум с применением при этом аппарата случайных чисел. Этот же метод позволяет неодно кратно производить математический эксперимент, т. е. проиг рывать на модели ситуацию несколько раз, причем ряд факто ров выбирается случайно, но в соответствии с законами их рас пределения. Другие параметры (неслучайные) при этом изме няются целенаправленно. В результате можно выбрать такие параметры, которые соответствуют оптимальному значению экономического показателя или оптимальной структуре.
В экономике и в технологии не всегда возможно провести некоторые эксперименты, так как всякий эксперимент связан с определенным риском получить отрицательный результат. Для проведения эксперимента в таких условиях можно использовать экономико-статистические модели. Результаты такого математи ческого экспериментирования только тогда заслуживают внима ния, когда эксперимент проводится на адекватной модели, и если моделирование осуществляется совершенными методами.
Метод статистических испытаний находит широкое приме нение при математическом имитационном моделировании вслед ствие того, что при помощи этого метода можно искусственно воспроизвести любой случайный процесс с заданными законами распределения. Строительный процесс [5], ввиду действия боль шого количества факторов, полностью учесть которые не пред ставляется возможным, также может быть отнесен к сложным динамическим стохастическим процессам. Поэтому на некото рых этапах планирования для нахождения параметров опти мального функционирования строительной организации целесо образно применить метод Монте-Карло.
Примером такого имитационного моделирования может быть решение задачи стохастического программирования, допустим
91
в постановке |
(5.5.8) — (5.5.11). |
Эту задачу нельзя решить ника |
ким другим |
методом из-за |
наличия ограничений (5.5.10) — |
— (5.5.11). При решении задачи задаются разные варианты де терминированных параметров, случайные величины моделируют ся в соответствии с законами их распределения. Решение такой задачи имитирует процесс производства, при котором достига ется максимум (5.5.8).
Но метод Монте-Карло может применяться и как чисто вы числительный, например, для вычисления определенного инте грала [2 ], для нахождения экстремума некоторых нелинейных функций (случайный поиск), для определения параметров кор реляционных формул и т. д. Основой метода является аппарат случайных чисел.
§ 7.2. Равномерная случайная последовательность чисел
Случайные процессы с известными параметрами могут быть реализованы в виде случайной последовательности чисел. Слу
чайная последовательность чисел — такая |
последовательность, |
||||
при которой нельзя наперед предсказать, какое число |
следу |
||||
ет за случайным числом gi_i. |
Наиболее простой и распростра |
||||
ненной является |
равномерная |
случайная |
последовательность |
||
чисел с пределами 0 и |
1. В дальнейшем эту последовательность |
||||
будем называть РСП |
(0,1). Принадлежность числа к РСП (0,1) |
||||
означает, что это |
число ^.больше или равно нулю и |
меньше |
|||
или равно единице, и вероятность появления такого числа |
|||||
Р, -= Р,> = ■• • Pi = Pi hi = Pn = |
— , |
(7.2.1) |
|||
|
|
|
|
n |
|
где п — количество чисел в последовательности, зависящее от числа знаков после запятой. Если после запятой идет один знак, т. е. g = 0; 0,1,..., 0,9; 1, то п = 11 и, следовательно, веро
ятность Р = |
если после запятой два знака, то п = 1 0 1 и |
11 ’
Р = _1_ и т. д.
101
Дифференциальный и интегральный законы распределения РСП (0,1) приведены на рис. 15.
Естественно, что
2 Р , = nP = 1, |
(7.2.2) |
i-i
что видно из определения вероятности Р (7.2.1) и графика 156).
92