Can I supply a gradient to fmincon when my objective function is anonymous?

Question

DavidZ 2014-8-13

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/147457-can-i-supply-a-gradient-to-fmincon-when-my-objective-function-is-anonymous

评论： Rustem Devletov 2019-4-26

I am using fmincon and would like to supply a gradient for the objective function to speed up convergence. The trouble is I use an anonymous objective function so I can pass parameters to my real objective function. I believe the anonymous function only picks up the first argument of my real objective function, which is the function value and not the gradient. The basic structure of my code is:

if true
%%%%%%%%%%%%%%%%%%%
[f,gradf]=realobjective(parameters,x]
f=blah;
gradf=blah;
end
%%%%%%%%%%%%%%%%%%%
options=optimoptions(@fmincon,'Algorithm','active-set','MaxFunEvals',1E5,'MaxIter',1E5,'TolFun',1E-6,'TolX',1E-6,'GradConstr','on','GradObj','on');
xguess=0;
func=@(x)realobjective(parameters,x)
argmin=fmincon(func,xguess,[],[],[],[],0,1000,nonlconstr,options)
end

So I think fmincon is not using the gradient I have defined in the real objective function. I'd be incredibly grateful for any advice.

David

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Rustem Devletov 2019-4-26

Hey, friend! I'm a fourth year student in Russia, currently writing my thesis. My supervisor asked me to provide gradient to fmincon

I have triend smth like this

function F = rosenbrockwithgrad(x)

F = 100*(x(2) - x(1)^2)^2 + (1-x(1))^2;

if nargout > 1 % gradient required

g = [-400*(x(2)-x(1)^2)*x(1)-2*(1-x(1));

200*(x(2)-x(1)^2)];

end

But he said that he wants it to be found in another file. How can I implement this? How do I find gradient in one file and pass it to another one?

请先登录，再进行评论。

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

Answer 1

Matt J 2014-8-13

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/147457-can-i-supply-a-gradient-to-fmincon-when-my-objective-function-is-anonymous#answer_148484

编辑：Matt J 2014-8-13

What you've posted looks fine. Call nargout inside realobjective (e.g., at a breakpoint), to test whether fmincon is calling it with two output arguments.

5 个评论
显示 3更早的评论隐藏 3更早的评论

Matt J 2014-8-13

Did you try running with 'DerivativeCheck' set to 'on' to validate your gradient calculation? Sounds like your constraint gradient calculation is working okay, but the objective function gradient calculation might have errors.

Rustem Devletov 2019-4-26

Hey, friend! I'm a fourth year student in Russia, currently writing my thesis. My supervisor asked me to provide gradient to fmincon

I have triend smth like this

function F = rosenbrockwithgrad(x)

F = 100*(x(2) - x(1)^2)^2 + (1-x(1))^2;

if nargout > 1 % gradient required

g = [-400*(x(2)-x(1)^2)*x(1)-2*(1-x(1));

200*(x(2)-x(1)^2)];

end

But he said that he wants it to be found in another file. How can I implement this? How do I find gradient in one file and pass it to another one?

请先登录，再进行评论。

Answer 2

DavidZ 2014-8-13

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/147457-can-i-supply-a-gradient-to-fmincon-when-my-objective-function-is-anonymous#answer_148501

Thanks both. It's very mysterious; when I turn the gradient calls off fmincon does converge (slowly) to the right answer (I know it's right as for certain parameter values I can work out an analytical solution to compare to the numerical solution fmincon gives) but when I turn the gradient calls on it never converges. I tried the debugging tip Matt helpfully suggested and fmincon is calling withing nargout 2 some of the time, but not all of it. Any thoughts? Thanks again for your help.