Оптимизация нелинейных функций

Минимизация функций одной переменной

Учитывая математическую функцию единственной переменной, можно использовать функцию fminbnd, чтобы найти локальный минимизатор функции в данном интервале. Например, рассмотрите функцию humps.m, которой предоставляют ^MATLAB®. Следующие данные показывают график humps.

x = -1:.01:2;
y = humps(x);
plot(x,y)
xlabel('x')
ylabel('humps(x)')
grid on

Чтобы найти минимум функции humps в области значений (0.3,1), использовать

x = fminbnd(@humps,0.3,1)

x = 0.6370

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

opts = optimset('Display','iter');
x = fminbnd(@humps,0.3,1,opts)

 
 Func-count     x          f(x)         Procedure
    1       0.567376      12.9098        initial
    2       0.732624      13.7746        golden
    3       0.465248      25.1714        golden
    4       0.644416      11.2693        parabolic
    5         0.6413      11.2583        parabolic
    6       0.637618      11.2529        parabolic
    7       0.636985      11.2528        parabolic
    8       0.637019      11.2528        parabolic
    9       0.637052      11.2528        parabolic
 
Optimization terminated:
 the current x satisfies the termination criteria using OPTIONS.TolX of 1.000000e-04

x = 0.6370

Итеративное отображение показывает текущую стоимость x и значения функции в f(x) каждый раз, когда функциональная оценка происходит. Для fminbnd одна функциональная оценка соответствует одной итерации алгоритма. Последний столбец показывает, какая процедура используется в каждой итерации, или золотой поиск раздела или параболическая интерполяция. Для получения дополнительной информации смотрите Итеративное Отображение.

Минимизация функций нескольких переменных

Функция fminsearch подобна fminbnd за исключением того, что это обрабатывает функции многих переменных. Задайте стартовый вектор x0, а не стартовый интервал. fminsearch пытается возвратить вектор x, который является локальным минимизатором математической функции около этого стартового вектора.

Чтобы попробовать fminsearch, создайте функциональный three_var трех переменных, x, y и z.

function b = three_var(v)
x = v(1);
y = v(2);
z = v(3);
b = x.^2 + 2.5*sin(y) - z^2*x^2*y^2;

Теперь найдите минимум для этой функции с помощью x = -0.6, y = -1.2 и z = 0.135 как начальные значения.

v = [-0.6,-1.2,0.135];
a = fminsearch(@three_var,v)

a =
    0.0000   -1.5708    0.1803

Максимизация функций

fminbnd и решатели fminsearch пытаются минимизировать целевую функцию. Если у вас есть проблема максимизации, то есть, проблема формы

$\underset{}{} maxxf (x),$

затем задайте g (x) = –f (x) и минимизируйте g.

Например, чтобы найти максимум загара (cos (x)) рядом x = 5, оцените:

[x fval] = fminbnd(@(x)-tan(cos(x)),3,8)

x =
    6.2832

fval =
   -1.5574

Максимум 1.5574 (отрицание fval, о котором сообщают) и происходит в x = 6.2832. Этот ответ правилен с тех пор к пяти цифрам, максимум коричнев (1) = 1.5574, который происходит в x = 2π = 6.2832.

`алгоритм fminsearch`

fminsearch использует алгоритм симплекса Nelder-меда, как описано в Lagarias и др. [1]. Этот алгоритм использует симплекс n + 1 точка для n-мерных векторов x. Алгоритм сначала делает симплекс вокруг исходного предположения x0 путем добавления 5% каждого компонента x0 (i) к x0. Алгоритм использует эти n векторы в качестве элементов симплекса в дополнение к x0. (Алгоритм использует 0.00025 в качестве компонента i если x0 (i) = 0.) Затем алгоритм неоднократно изменяет симплекс согласно следующей процедуре.

Примечание

Ключевые слова для fminsearch итеративное отображение появляются полужирным после описания шага.

Позвольте x (i), обозначают список точек в текущем симплексе, i = 1..., n+1.
Закажите точки в симплексе от самого низкого значения функции f (x (1)) к самому высокому f (x (n+1)). На каждом шаге в итерации алгоритм отбрасывает текущую худшую точку x (n+1) и принимает другую точку в симплекс. [Или, в случае шага 7 ниже, это изменяет весь n точек со значениями, больше, чем f (x (1))].
Сгенерируйте отраженную точку
r = 2 м – x (n+1),
где
m = Σx (i)/n, i = 1... n,
и вычислите f (r).
Если f (x (1)) ≤ f (r) <f (x (n)), примите r и отключите эту итерацию. Отразиться
Если f (r) <f (x (1)), вычислите точку s расширения
s = m + 2 (m – x (n+1)),
и вычислите f (s).
1. Если f (s) <f (r), примите s и отключите итерацию. Расширение
2. В противном случае примите r и отключите итерацию. Отразиться
Если f (r) ≥ f (x (n)), выполните сокращение между m и лучше x (n+1) и r:
1. Если f (r) <f (x (n+1)) (то есть, r лучше, чем x (n+1)), вычислить
  c = m + (r – m)/2
  и вычислите f (c). Если f (c) < f (r), примите c и отключите итерацию. Контракт снаружи В противном случае, продолжите Шаг 7 (Уменьшение).
2. Если f (r) ≥ f (x (n+1)), вычислить
  cc = m + (x (n+1) – m)/2
  и вычислите f (cc). Если f (cc) < f (x (n+1)), примите cc и отключите итерацию. Контракт внутри В противном случае, продолжите Шаг 7 (Уменьшение).
Вычислите n точек
v (i) = x (1) + (x (i) – x (1))/2
и вычислите f (v (i)), i = 2..., n+1. Симплекс в следующей итерации является x (1), v (2)..., v (n+1). Уменьшение

Следующие данные показывают точки, что fminsearch может вычислить в процедуре, наряду с каждым возможным новым симплексом. Исходный симплекс имеет полужирную схему. Итерации продолжают, пока они не соответствуют останавливающемуся критерию.

Ссылка

[1] Lagarias, J. C. J. A. Тростники, М. H. Райт и П. E. Райт. “Convergence Properties Симплекс-метода Nelder-меда в Низких Размерностях”. SIAM Journal Оптимизации, Издания 9, Номера 1, 1998, стр 112–147.

Документация