Imposing constraint on gradient coefficients in NonLinear optimization

Question

Jonathan 2014-4-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/126105-imposing-constraint-on-gradient-coefficients-in-nonlinear-optimization

评论： Sean de Wolski 2014-4-18

I'm trying to solve a non linear 6 dimensional optimization problem, using a gradient which I calculate (6 dimensional).

I would like to influence the gradient coefficients in such a way that Coeff(1) == Coeff(2) == Coeff(3) and Coeff(4)==Coeff(5)==Coeff(6)

Is this possible?

Thanks! Jonathan

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Alan Weiss 2014-4-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/126105-imposing-constraint-on-gradient-coefficients-in-nonlinear-optimization#answer_133776

在 MATLAB Online 中打开

I suppose that you could try to minimize the sum of the squares of the differences between various coefficients:

 (Coeff(1) - Coeff(2))^2 + (Coeff(1) - Coeff(3))^2 + (Coeff(4) - Coeff(5))^2 + (Coeff(4) - Coeff(6))^2

You might want to use lsqnonlin and have your vector be (coeff(1) - coeff(2), coeff(1) - coeff(3), ...)

Alan Weiss

MATLAB mathematical toolbox documentation

2 个评论
显示无隐藏无

Jonathan 2014-4-18

Thanks for the answer, but that wasn't quite what I was asking. I will try rephrasing.

I am trying to solve an optimization problem, minimizing a scalar energy functional, which is a function of 6 parameters - 3 translation parameters and 3 rotation parameters.

I'm trying to feed the 6 dimensional calculated gradient into a gradient based optimization algorithm.

The constrain I want to impose is that in each iteration, all translation gradients have the same step size, and all rotation gradients have a second step size.

If t is the translation vector, r the rotation vector, at iteration i the propagation will be:

[t(1, i), t(2, i), t(3, i), r(1, i), r(2, i), r(3, i)] = [t(1, i-1), t(2, i-1), t(3, i-1), r(1, i-1), r(2, i-1), r(3, i-1)] - [a, a, a, b, b, b].*[grad(1), grad(2), grad(3), grad(4), grad(5), grad(6)]

Is this possible? Thanks Jonathan