Файл: Чесноков, Н. И. Оптимизация решений при разработке урановых месторождений.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.10.2024
Просмотров: 168
Скачиваний: 0
Третьим известным методом решения задач нелиней ного программирования может быть назван градиент ный метод.
Наконец, последний метод — классическая оптимиза ция, заключающаяся в применении дифференциального исчисления. Однако практически этот метод приме няется очень редко из-за трудоемкости вычислений. Классические методы оптимизации применяются в ос новном при теоретических исследованиях. Три первых метода обычно применяют для решения практических задач ( о н и будут рассмотрены ниже).
Преобразование нелинейных задач в приближенные линейные
Преобразование нелинейных задач в линейные является одним нз способов приведения нелинейных за дач к такому виду, чтобы можно было для 11X решения
применить симплекс-метод.
Некоторые нелинейные функции могут быть пред ставлены в виде суммы линейных функций, например
У — f (-'Т. -V.,, . . . , .ѵ„) = /, (лу) і- /о (л-2) |
+ /„ (.ѵ„). |
|
(3.69) |
Нелинейные функции, которые могут быть приведены к линейному виду указанным способом, называются се парабельными.
В нелинейных задачах нелинейными могут быть как целевая функция, так и ограничения.
Рассмотрим пример преобразования нелинейной за дачи в линейную. Представим себе, что па одном из рудников возможно применение двух систем разработ ки— камерной с массовым обрушением целиков и гори зонтальными слоями с полной закладкой выработанно го пространства твердеющими смесями п селективной выемкой. Первая система позволяет достигнуть высокой производительности труда и низкой себестоимости до бычи 1 т руды, но приводит к большому разубоживаншо. Вторая система позволяет добыть руду с меньшим разубоживанпем, но очень трудоемка, и себестоимость
добычи 1 |
т руды при |
ее использовании значительно |
||
выше. |
|
|
объем добычи руды |
камерной |
Обозначим через |
ад |
|||
системой |
разработки |
с |
массовым обрушением |
целиков |
176
(млн. r/год), через х2 объем добычи руды горизонталь ными слоями с полной закладкой выработанного прост ранства твердеющей смесыо и селективной выемкой (млн. т/год). На основании исследований па данном руднике были получены следующие соотношения: себе стоимость 1 кг металла в добытой руде (С|)
Ci — 2Xl + X.,;
себестоимость добычи и переработки 1 кг металла (С2)
С2 = 2xi -f Зхг;
содержание металла в добытой руде (Z, %)
Z 1 12лу — x^ -j-- Хо.
Задача состоит в том, чтобы при себестоимости до бычи 1 кг металла в руде не выше 4 руб. и себестоимо сти добычи и переработки не более 6 руб. за 1 кг металла определить максимальное содержание металла в добытой руде и соотношение систем разработки, обес печивающих соблюдение указанных условий.
Иначе данную задачу можно переписать следующим
образом: найти |
|
|
max Z = 2Xl — xf -j- x3 |
|
|
при условиях |
|
|
2Хі + |
х2 <£4, |
(3.70) |
2хг + |
Зх2 < 6. |
|
Как видно из задачи, целевая функция нелинейная, а ограничения линейные.
Попробуем целевую функцию представить в виде, удобном для применения симплекс-метода, т. е. преобра зовать ее в линейную. Перепишем целевую функцию в виде суммы двух функций:
max Z = /у (хх) + f2(х2), |
|
где fi(xi) =2xi—x j , f2(x2) = а'2 ( Х і > 0 , х 2 > 0 ) . |
Поскольку |
fz(х2) = х2— линейная функция, мы не будем |
ее преоб |
разовывать и выражать через другие переменные, пре образующие нелинейные зависимости в кусочно-линей ные.
Для преобразования /у(хі) разобьем область изме нения Хі на интервалы Ах,(1-= 0,25. При этом следует
12 Н. И. Чесноков в др. |
177 |
иметь в виду, что ли не может быть более 2,0, как это видно из функции М л-,) 11 ограничений, т. е. 0 ^ л',^2.
В |
общ ем виде Дл-Л,■ = л |
* |
,■ при |
/ г = ( 1 , |
/г.), |
Д/л; = |
|||
= /лг— / а- і . і - Е сли |
г ^ л у - г ^ л д ;, то |
|
|
|
|
||||
|
|
+■ = АѴ - і . ( + |
( Л л Д ) 6 * ; > |
|
|
( 3 . 7 1 ) |
|||
где |
б,., = Х‘ |
&хкі |
и |
0 < 6іи < 1. |
|
|
|
|
|
Тогда на основании изложенного можно написать |
|
||||||||
|
|
|
|
|
|
|
|
|
(3.72) |
где |
fi(Xi)— кусочно-линейное |
приближение |
функции |
||||||
Ы + ). |
|
|
|
|
|
|
k— 1) и |
||
Условимся |
далее, |
что |
6«і= 1 при |
«= (1, |
|||||
6Лі- = 0. В этом случае получим |
|
|
|
|
|
||||
|
|
^k—і,і = |
к |
|
|
|
|
(3.73) |
|
|
|
2 |
|
|
|
|
|||
|
|
|
|
|
|
|
|||
|
|
h - , ,■= и—1 |
+ |
|
|
|
(3.74) |
||
|
|
|
|
а |
также |
||||
Итак, допустив, что 0<Дб< 1, |
61U= 0 |
и |
и>/г, |
||||||
Shi > 0 и 6иі = 1 («=(1, |
k— 1)], |
можно |
написать |
|
|
||||
|
|
|
пі |
|
|
|
|
(3.75) |
|
|
|
Хі = |
(Л%) |
|
|
|
|||
|
|
I! № |
|
+ ""’о-’ |
|
|
(3.76) |
||
|
|
|
|
|
|
|
|
|
|
Вернемся к рассматриваемому примеру. Для условия |
|||||||||
Ах/И; = 0,25 результаты |
расчетов по формулам |
(3.70) — |
(3.76) приведены в табл. 42. |
|
|
|||||
С учетом |
изложенного условия рассматриваемой за |
||||||
дачи |
можно |
переписать в |
следующем виде: |
0,5^ 6,,+ |
|||
+ З л-2 + л' з = |
6, |
0 , 5 ^ 6 ft+ х 2 + |
х « = 4 - ( ° < ö h < l ; |
/г = |
1 - 2 > - ’ 8 ’ |
||
х2, |
Л'з, |
х . ^ 0 ) . |
Найти |
max 21=0,43756,+0,312502 + |
+ 0,18756з + 0,062564 — 0,06255s — 0,187566 — 0,312567 —
—0,437508 + л'г. Здесь |
Z — аналитическое выражение ло |
||
маной, |
аппроксимирующей целевую функцию |
Z; ö,t — |
|
новые |
переменные; |
л'з, ли — вспомогательные |
перемен- |
178
|
|
|
Т а б л и ц а 42 |
*1 |
Л (.ѵ,)=2.ѵ, - л-2 |
Щ і |
foi |
0 |
0 |
0 |
0 |
0,25 |
0,4375 |
0,4375 |
0 |
0,50 |
0,7500 |
0,3125 |
0 |
0,75 |
0,9375 |
0,1875 |
0 |
1,00 |
1,0000 |
0,0625 |
0 |
1,25 |
0,9375 |
—0,0625 |
0 |
1,5 |
0,7500 |
—0,1875 |
0 |
1,75 |
0,4375 |
—0,3125 |
0 |
2,00 |
0 |
—0,4375 |
0 |
8
ныв; 2 öft= 0,2561 + 0,2502 0,256з + 0,2564 + 0,256s + 0,256g-Ь
”Ь 0,2507 ~Ь 0,2568-
Таким образом, нелинейная задача преобразована в линейную и может быть решена симплекс-методом.
Применив для решения уже преобразованной задачи
симплекс-метод, |
получим |
|
at= 0 , 7 5 , |
А г = 1 , 5 , max Z— |
= 2,4374. Таким |
образом, |
максимальное содержание |
||
металла в добытой руде 2 |
= 2 ,4 3 7 4 % , |
что обеспечивается |
за счет добычи руды камерной системой с обрушением
целиков в объеме 0,75 млн. т/год |
и системой |
горизон |
|||
тальных |
слоев |
с |
полной за кладкой выработанного |
||
пространства |
и |
селективной |
выемкой в |
объеме |
|
1,5 млн. |
т/год. |
Это |
обеспечивает |
себестоимость |
добычи |
1 кг металла в руде не выше 4 руб. и себестоимость до
бычи и переработки 1 кг металла не более 6 руб.
В связи с большой трудоемкостью расчеты здесь не приводятся. В работах [26, 35, 38, 55] можно более под робно познакомиться и с другими методами линейной аппроксимации нелинейных задач, включая случаи не линейности как целевой функции, так и ограничений и квадратичного программирования, также позволяющего применять для решения симплекс-метод.
Градиентный метод решения нелинейных задач
Градиентный метод является также приближен ным и заключается в последовательном приближении к решению задачи. Каждый шаг решения при использова
12* 179
нии данного метода постепенно приближает к точке максимума пли минимума в зависимости от поставлен ной задачи. Граднентный метод применяется в основном при решении нелинейных задач на быстродействующих электронно-вычислительных машинах. Он не требует сепарабельности функции, однако необходимым усло вием является дифференцируемость и непрерывность как самих функций, так и их первых и вторых произ водных.
Метод назван градиентным потому, что при решении нелинейных задач используется вектор частных произ водных функций, называемый градиентом функции.
Если, например, частные производные вычислены в ка кой-то определенной точке у \ = { х п , *12, •••,*іл), то гра
диентом будет в данном случае либо вектор, либо точка в /г-мерном пространстве. При этом вектор из этой точ ки направлен в сторону максимальной скорости измене ния значения функции.
|
При использовании градиентного метода решения |
|||
нелинейных задач сначала выбирают |
начальную (ис |
|||
ходную) точку. Пусть |
начальной точкой будет |
г/і=(х'ц, |
||
*і2, |
-tin). Затем определяют значение частной произ |
|||
водной в данной точке (А/): |
|
|
||
|
М ( У і ) — |
( * i i > * 1 2 ) • - • > * і п ) > |
( 3 . 7 7 ) |
|
где |
Z\ —первая частная производная |
целевой |
функции, |
; = ( 1, /г). После этого находят значения компонент вто рой точки у2 по формуле r/2j= і/іj+ ctiAf(t/ij); аналогич
но далее
( 3 . 7 8 )
где а — коэффициент, порядок определения которого бу дет рассмотрен в числовом примере.
В результате выполнения указанных расчетов полу
чается последовательность точек ij\, |
у2, у3, ...,г/„, которые |
с каждым шагом приближают к |
искомому оптимуму |
(точке максимума или минимума в пределах области, фиксируемой ограничениями).
При определении значения а следует иметь в виду, что при слишком большом а число шагов будет меньше,
180