Is there any other alternative fuction that I can use instead of fminunc?.

8 次查看（过去 30 天）
salva
salva 2012-9-20
Dear all,
I am using the function fminunc as I want to find the mode and the hessian matrix of a function. I use Matlab 2009 but I get an error related to this function
??? Error using ==> roots at 28
Input to ROOTS must not contain NaN or
Inf. Error in ==>
lineSearch>globalMinimizerOfPolyInInterval
at 304
stationaryPoint = roots([3*coeff(1)
2*coeff(2) coeff(3)]); Error in ==>
lineSearch>pickAlphaWithinInterval at 279
z =
globalMinimizerOfPolyInInterval(zlb,zub,coeff); Error in ==> lineSearch>sectioningPhase at
193
alpha =
pickAlphaWithinInterval(brcktEndpntA,brcktEndpntB,a,b,f_a,fPrime_a,f_b,fPrime_b); Error in ==> lineSearch at 55
[alpha,f_alpha,grad,exitflag,funcCount]
= sectioningPhase(funfcn, ... Error in ==> fminusub at 208
[alpha,f,grad,exitflagLnSrch,funcCountLnSrch]
= ... Error in ==> fminunc at 376
[x,FVAL,GRAD,HESSIAN,EXITFLAG,OUTPUT] =
fminusub(funfcn,x, ... Error in ==> testxa2 at 252
[xxx,fval,exitflag,output,gradient,H]=fminunc('modeA,x00,options,...
Is there any other alternative fuction that I can use instead of fminunc?. The function produces no error in Matlab2012 but only in matlab 2009
thank you

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

回答（3 个）

Matt Kindig
Matt Kindig 2012-9-20
Hi Salva,
You could use fminsearch or fminbnd, two derivative-free optimizers. However, given that the error is thrown when the input to roots() is invalid, it would probably be better to catch that error before entering roots. You could put a check in your globalMinimizerOfPolyInInterval() function prior to calling roots() that checks that coeff(1:3) are all not NaN or Inf. Something like:
if all(isfinite(coeffs(1:3)),
stationaryPoint = roots([3*coeff(1) 2*coeff(2) coeff(3)]);
else
%do something to force the optimizer to reject this solution
end
  2 个评论
salva
salva 2012-9-21
thank you Matt. It seems that the
else
%do something to force the optimizer to reject this solution
is the most important of all.
Any suggestion for
%do something to force the optimizer to reject this solution?
thankss

请先登录，再进行评论。

Walter Roberson
Walter Roberson 2012-9-23
If you use an options structure that turns on function value checking, then the minimization will stop at the point that inf or NaN is produced.
  3 个评论
显示 1更早的评论
Walter Roberson
Walter Roberson 2012-9-25
You appear to already be constructing options, so include 'FunValCheck', 'on' when constructing the options.
salva
salva 2012-9-25
thank you walter. Actually, I am looking the opposite. I want the process to continue even when the point that is being produced is inf or NaN (as in matlab2012)
thanks

请先登录，再进行评论。

Sean de Wolski
Sean de Wolski 2012-9-25
fmincon's Interior-Point algorithm can recover from inf and nan output.
doc fmincon
  8 个评论
显示 6更早的评论
Sean de Wolski
Sean de Wolski 2012-10-1
Well since you don't need constraints you can ignore all of those inputs, just put some lower and upper bounds on your variables and use it just like fminunc
fmincon(foo,x0,[],[],[],[],lb,ub,options);

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Solver Outputs and Iterative Display 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by