36

Штраф, заданный обратной функцией, имеет вид

Ω = R [1/g(x)].

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

жении внутрь допустимой области. На самой границе значение P(x,R) не опреде-

лено, как и в предыдущем случае возможно появление недопустимых точек.

Штраф типа квадрата срезки имеет вид

W = R g(x) 2 ,

ìa, если a £ 0,

ï

где a = í

ï0, если a > 0.

î

Этот штраф является внешним и недопустимые точки не создают проблем по сравнению с допустимыми. Различие заключается в том, что в допустимых точ-

ках штраф равен нулю. Этот вид штрафа удобен тем, что Р(х,R) непрерывна и определена всюду. Параметр R положителен и увеличивается от итерации к ите-

рации.

2.3.3. Алгоритм метода

Шаг 1. Задать начальные данные N, J, K, ε1, ε2, ε3, x(0), R(0), где

ε1 – параметр окончания одномерного поиска (если таковой используется в процедуре безусловной оптимизации);

ε2 - параметр окончания процедуры безусловной оптимизации;

ε3 - параметр окончания работы алгоритма; x(0) - начальная точка;

R(0) - начальный вектор штрафных параметров.

Шаг 2. Построить штрафную функцию

P(x, R) = f(x) + Ω (R, g(x), h(x)).

Шаг 3. Найти х(t+1), доставляющий экстремум Р(х(t+1), R(t)) при фиксирован-