Некоторые задачи требуют, чтобы вы нашли точку, которая удовлетворяет всем ограничениям, без целевой функции, чтобы минимизировать. Например, предположим, что у вас есть следующие ограничения:
Делайте все точки удовлетворить ограничениям? Чтобы выяснить, запишите функцию, которая возвращает ограничения в поле структуры Ineq
. Запишите ограничения в терминах двухэлементного вектора вместо . Запишите каждое неравенство как функцию , что означает неравенство , путем вычитания правой стороны каждого неравенства с обеих сторон. Чтобы включить построение графика, запишите функцию векторизированным способом, где каждая строка представляет одну точку. Код для этой вспомогательной функции с именем objconstr
, появляется в конце этого примера.
Постройте график точек, где три функции удовлетворяют равенствам для и , и указать неравенства путем построения графиков линий для значений функции, равных -1/2.
[XX,YY] = meshgrid(-2:0.1:2,-4:0.1:2); ZZ = objconstr([XX(:),YY(:)]).Ineq; ZZ = reshape(ZZ,[size(XX),3]); h = figure; ax = gca; contour(ax,XX,YY,ZZ(:,:,1),[-1/2 0],'r','ShowText','on'); hold on contour(ax,XX,YY,ZZ(:,:,2),[-1/2 0],'k','ShowText','on'); contour(ax,XX,YY,ZZ(:,:,3),[-1/2 0],'b','ShowText','on'); hold off
График показывает, что допустимые точки существуют около [1,75, -3].
Установите нижние границы -5 и верхние границы 3 и решите задачу используя surrogateopt
.
rng(1) % For reproducibility
lb = [-5,-5];
ub = [3,3];
[x,fval,exitflag,output,trials] = surrogateopt(@objconstr,lb,ub)
surrogateopt stopped because it exceeded the function evaluation limit set by 'options.MaxFunctionEvaluations'.
x = 1×2
1.7312 -3.0312
fval = 1x0 empty double row vector
exitflag = 0
output = struct with fields:
elapsedtime: 31.7004
funccount: 200
constrviolation: -0.0800
ineq: [-0.0800 -0.2083 -0.7625]
rngstate: [1x1 struct]
message: 'surrogateopt stopped because it exceeded the function evaluation limit set by ...'
trials = struct with fields:
X: [200x2 double]
Ineq: [200x3 double]
Проверьте допустимость в возвращенном решении x
.
disp(output.ineq)
-0.0800 -0.2083 -0.7625
Эквивалентно вычислите функцию objconstr
в возвращённом решении x
.
disp(objconstr(x).Ineq)
-0.0800 -0.2083 -0.7625
Эквивалентно исследуйте Ineq
поле в trials
структура для решения x
. Сначала найдите индекс x
в trials.X
поле.
indx = ismember(trials.X,x,'rows');
disp(trials.Ineq(indx,:))
-0.0800 -0.2083 -0.7625
Все значения ограничительных функций отрицательны, что указывает на то, что точка x
является допустимым.
Просмотрите допустимые точки, оцениваемые surrogateopt
.
opts = optimoptions("surrogateopt"); indx = max(trials.Ineq,[],2) <= opts.ConstraintTolerance; % Indices of feasible points figure(h); hold on plot(trials.X(indx,1),trials.X(indx,2),'*') xlim([1 2]) ylim([-3.5 -2.5]) hold off
Этот код создает objconstr
вспомогательная функция.
function f = objconstr(x) c(:,1) = (x(:,2) + x(:,1).^2).^2 + 0.1*x(:,2).^2 - 1; c(:,2) = x(:,2) - exp(-x(:,1)) + 3; c(:,3) = x(:,2) - x(:,1) + 4; f.Ineq = c; end