Two-dimensional optimization (especially using GLPK)
显示 更早的评论
Suppose I want to perform an optimization where the solution matrix X is a two-dimensional matrix rather than a vector. The constraints and objective functions can therefore operate on various rows/columns of the solution matrix; for example, the sum of the elements in each column of the matrix must be as close as possible to a certain value, but the sum of n consecutive elements in each row must be less than or equal to 1.
Are there any functions or combinations of functions in MATLAB that allow me to do this, and is it also possible to do this using the GLPK package?
回答（1 个）
The examples you describe can be handled by any solver just be reshaping X into a vector first. For example the problem you describe is equivalent to
[m,n]=size(A);
ec=ones(m,1);
er=ones(n,1);
S=kron(speye(n),ec.');
A=kron(er.',speye(m));
x=X(:);
min norm(S*x -certainvalue).^2
s.t.
A*x<=ec
Also, in MATLAB Optimization Toolbox functions, no constraints or objective functions that are expressed via function handles (instead of matrix data) ever need to process X as a vector. Same with ub,lb data, I think.
6 个评论
Mike Vukovich
2013-6-21
No, in both linprog and lsqlin, the objective function and constraints are specified using matrices, A,b,Aeq,beq etc... and not using function handles. So, for those, you would need to reformulate the problem in terms of the vector X(:), using techniques in my example above.
But in lsqnonlin, for example, the objective is expressed using a function handle, so you could write that function to accept a 2D matrix argument X as input. The ub,lb data could also be expressed in 2D form.
In fmincon, meanwhile, you have a hybrid situation. The objective is expressed as a function handle and so can accept a matrix argument X. However, fmincon uses matrices A,Aeq, to specify linear constraints. The fmincon algorithms will internally perform operations like A*X(:) and Aeq*X(:) and so A, Aeq have to be constructed accordingly.
Matt J
2013-6-21
See also the documentation
Mike Vukovich
2013-6-21
编辑：Mike Vukovich
2013-6-21
Matt J
2013-6-21
I have no experience with GLPK/Matlab, I'm afraid. I think if it allowed you to formulate the argument space in matrix form, that would be very obvious in the documentation. I can't see much benefit in having this for linear solvers, though, so I would be surprised if it were supported. Same thing with "Cx-d".
类别
在 帮助中心 和 File Exchange 中查找有关 Solver Outputs and Iterative Display 的更多信息
2013-6-21
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
