Exponenta Banner

Решите нелинейную задачу выполнимости, основанную на проблеме

Скрипт Open Live Script

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

Описание задачи

Предположим, что у вас есть следующие ограничения:

(y+x2)2+0.1y2≤1y≤exp(-x)-3y≤x-4.

Сделайте любые точки (x,y) удовлетворить всем ограничениям?

Основанное на проблеме решение

Создайте задачу оптимизации, которая имеет только ограничения, никакую целевую функцию.

x = optimvar('x');
y = optimvar('y');
prob = optimproblem;
cons1 = (y + x^2)^2 + 0.1*y^2 <= 1;
cons2 = y <= exp(-x) - 3;
cons3 = y <= x - 4;
prob.Constraints.cons1 = cons1;
prob.Constraints.cons2 = cons2;
prob.Constraints.cons3 = cons3;
show(prob)
  OptimizationProblem : 

	Solve for:
       x, y

	minimize :

	subject to cons1:
       ((y + x.^2).^2 + (0.1 .* y.^2)) <= 1

	subject to cons2:
       y <= (exp((-x)) - 3)

	subject to cons3:
       y - x <= -4

Создайте псевдослучайную структуру стартовой точки x0 с полями x и y для переменных оптимизации.

rng default
x0.x = randn;
x0.y = randn;

Решите задачу, начинающую с x0.

[sol,~,exitflag,output] = solve(prob,x0)
Solving problem using fmincon.

Local minimum found that satisfies the constraints.

Optimization completed because the objective function is non-decreasing in 
feasible directions, to within the value of the optimality tolerance,
and constraints are satisfied to within the value of the constraint tolerance.

sol = struct with fields:
    x: 1.7903
    y: -3.0102


exitflag = 
    OptimalSolution


output = struct with fields:
              iterations: 6
               funcCount: 9
         constrviolation: 0
                stepsize: 0.2906
               algorithm: 'interior-point'
           firstorderopt: 0
            cgiterations: 0
                 message: '...'
            bestfeasible: [1x1 struct]
     objectivederivative: "closed-form"
    constraintderivative: "forward-AD"
                  solver: 'fmincon'

Решатель находит допустимую точку.

Важность начальной точки

Решатель может не найти решение при запуске с некоторых начальных точек. Установите начальную точку x0.x = –1, x0.y = –4 и затем решите задачу, начинающую с x0.

x0.x = -1;
x0.y = -4;
[sol2,~,exitflag2,output2] = solve(prob,x0)
Solving problem using fmincon.

Converged to an infeasible point.

fmincon stopped because it is unable to find a point locally that satisfies
the constraints within the value of the constraint tolerance.


Consider enabling the interior point method feasibility mode.

sol2 = struct with fields:
    x: -2.1266
    y: -4.6657


exitflag2 = 
    NoFeasiblePointFound


output2 = struct with fields:
              iterations: 121
               funcCount: 279
         constrviolation: 1.4609
                stepsize: 1.9997e-10
               algorithm: 'interior-point'
           firstorderopt: 0
            cgiterations: 261
                 message: '...'
            bestfeasible: []
     objectivederivative: "closed-form"
    constraintderivative: "forward-AD"
                  solver: 'fmincon'

Проверяйте infeasibilities в возвращенной точке.

inf1 = infeasibility(cons1,sol2)
inf1 = 1.1974

inf2 = infeasibility(cons2,sol2)
inf2 = 0

inf3 = infeasibility(cons3,sol2)
inf3 = 1.4609

Оба cons1 и cons3 неосуществимы в решении sol2. Результаты подсвечивают важность использования нескольких стартовых точек, чтобы исследовать и решить задачу выполнимости.

Визуализируйте ограничения

Чтобы визуализировать ограничения, постройте точки, где каждая ограничительная функция является нулем при помощи fimplicit. fimplicit функционируйте передает числовые значения его функциям, тогда как evaluate функция требует структуры. Чтобы связать эти функции, используйте evaluateExpr функция помощника, которая появляется в конце этого примера. Эта функция просто помещает переданные значения в структуру с соответствующими именами.

Примечание: Убедитесь код для evaluateExpr функция помощника включена в конце вашего скрипта или в файле на пути.

Избегайте предупреждения, которое происходит потому что evaluateExpr функция не работает над векторизованными входными параметрами.

s = warning('off','MATLAB:fplot:NotVectorized');
cc1 = (y + x^2)^2 + 0.1*y^2 - 1;
fimplicit(@(a,b)evaluateExpr(cc1,a,b),[-2 2 -4 2],'r')
hold on
cc2 = y - exp(-x) + 3;
fimplicit(@(a,b)evaluateExpr(cc2,a,b),[-2 2 -4 2],'k')
cc3 = y - x + 4;
fimplicit(@(x,y)evaluateExpr(cc3,x,y),[-2 2 -4 2],'b')
hold off

Figure contains an axes object. The axes object contains 3 objects of type implicitfunctionline.

warning(s);

Выполнимая область в красной схеме и ниже в синяках линии. Выполнимая область в нижнем правом углу красной схемы.

Функция помощника

Этот код создает evaluateExpr функция помощника.

function p = evaluateExpr(expr,x,y)
pt.x = x;
pt.y = y;
p = evaluate(expr,pt);
end

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

Похожие темы

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

Поддержка