Одномерные ограничения полу-Бога

Найдите значения x , которые минимизируют

f (x) = (x ₁ – 0.5)² + (x ₂ – 0.5)² + (x ₃ – 0.5)²

где

$\begin{array}{l} K_{1} (x, w_{1}) = \sin (w_{1} x_{1}) \cos (w_{1} x_{2}) - \frac{1}{1000} {(w_{1} - 50)}^{2} - \sin (w_{1} x_{3}) - x_{3} \leq 1, \\ K_{2} (x, w_{2}) = \sin (w_{2} x_{2}) \cos (w_{2} x_{1}) - \frac{1}{1000} {(w_{2} - 50)}^{2} - \sin (w_{2} x_{3}) - x_{3} \leq 1, \end{array}$

для всех значений w ₁ и w ₂ в областях значений

1 ≤ w ₁ ≤ 100,
1 ≤ w ₂ ≤ 100.

Обратите внимание на то, что полубесконечные ограничения одномерны, то есть, векторы. Поскольку ограничения должны быть в форме _Ki (x, _wi) ≤ 0, необходимо вычислить ограничения как

$\begin{array}{l} K_{1} (x, w_{1}) = \sin (w_{1} x_{1}) \cos (w_{1} x_{2}) - \frac{1}{1000} {(w_{1} - 50)}^{2} - \sin (w_{1} x_{3}) - x_{3} - 1 \leq 0, \\ K_{2} (x, w_{2}) = \sin (w_{2} x_{2}) \cos (w_{2} x_{1}) - \frac{1}{1000} {(w_{2} - 50)}^{2} - \sin (w_{2} x_{3}) - x_{3} - 1 \leq 0. \end{array}$

Во-первых, запишите файл, который вычисляет целевую функцию.

function f = myfun(x,s)
% Objective function
f = sum((x-0.5).^2);

Во-вторых, запишите файл mycon.m это вычисляет нелинейное равенство и ограничения неравенства и полубесконечные ограничения.

function [c,ceq,K1,K2,s] = mycon(X,s)
% Initial sampling interval
if isnan(s(1,1)),
   s = [0.2 0; 0.2 0];
end
% Sample set
w1 = 1:s(1,1):100;
w2 = 1:s(2,1):100;

% Semi-infinite constraints 
K1 = sin(w1*X(1)).*cos(w1*X(2)) - 1/1000*(w1-50).^2 -...
       sin(w1*X(3))-X(3)-1;
K2 = sin(w2*X(2)).*cos(w2*X(1)) - 1/1000*(w2-50).^2 -...
       sin(w2*X(3))-X(3)-1;

% No finite nonlinear constraints
c = []; ceq=[];

% Plot a graph of semi-infinite constraints
plot(w1,K1,'-',w2,K2,':')
title('Semi-infinite constraints')
drawnow

Затем вызовите стандартную программу оптимизации.

x0 = [0.5; 0.2; 0.3];      % Starting guess
[x,fval] = fseminf(@myfun,x0,2,@mycon);

После восьми итераций решение

Значение функции и максимальные значения полубесконечных ограничений в решении x

fval
fval =
    0.0771

[c,ceq,K1,K2] = mycon(x,NaN); % Initial sampling interval
max(K1)
ans =
   -0.0077
max(K2)
ans =
   -0.0812

График полубесконечных ограничений производится.

Этот график показывает, как peaks в обоих ограничениях находится на границе ограничений.

Команда plot в mycon.m замедляет расчет. Удалите эту линию, чтобы улучшить скорость.

Смотрите также

fseminf

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

Поддержка

Памятка переводчика

1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.

2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.

3. Сохраняйте структуру оригинального текста - например, не разбивайте одно предложение на два.

4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.

5. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.

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