Файл: Чесноков, Н. И. Оптимизация решений при разработке урановых месторождений.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) и

Условимся

далее,

что

і= 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