Создайте проблему оптимизации.

x = optimvar('x',3);
y = optimvar('y',3,3);
prob = optimproblem('Objective',x'*y*x);

Преобразование проблемы в структуру.

problem = prob2struct(prob);

Получение линейных индексов в problem из всех prob переменные.

idx = varindex(prob);
disp(idx.x)

     1     2     3

disp(idx.y)

     4     5     6     7     8     9    10    11    12

Получить y только индексы.

idxy = varindex(prob,'y')

idxy = 1×9

     4     5     6     7     8     9    10    11    12

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

x = optimvar('x',3,1,'LowerBound',1,'UpperBound',1);
y = optimvar('y',3,3,'LowerBound',-1,'UpperBound',1);
prob = optimproblem('Objective',x'*y*x + [2 3 4]*x);
rng default
x0.x = rand(3, 1);
x0.y = rand(3, 3);
[solp,fvalp,exitflagp,outputp] = 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.

Затем преобразуйте проблему в форму на основе решателя с помощью prob2struct. Чтобы иметь fmincon решатель использовать автоматические градиенты в задаче, установить SpecifyObjectiveGradient опция для true.

solverprob = prob2struct(prob,x0);
solverprob.options = optimoptions(solverprob.options,"SpecifyObjectiveGradient",true);

Решить проблему с помощью fmincon.

[sols,fvals,exitflags,outputs] = fmincon(solverprob);

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.

Чтобы преобразовать fmincon решение в форме структуры, возвращенной solve, создайте соответствующие структуры с помощью varindex.

idx = varindex(prob);
sol.x = sols(idx.x);
sol.y = sols(idx.y);

y индекс, который varindex использует является линейным индексом. Изменение формы переменной sol.y чтобы иметь размер x0.y.

sol.y = reshape(sol.y,size(x0.y));

Проверьте, что две структуры решения идентичны.

isequal(sol,solp)

ans = logical
   1

Причина того, что два подхода не являются полностью эквивалентными, заключается в том, что fmincon может возвращать больше аргументов, таких как множители Лагранжа, тогда как solve не может.