Is there any other alternative fuction that I can use instead of fminunc?.
7 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
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
0 comentarios
Respuestas (3)
Matt Kindig
el 20 de Sept. de 2012
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 comentarios
Walter Roberson
el 23 de Sept. de 2012
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 comentarios
Walter Roberson
el 25 de Sept. de 2012
You appear to already be constructing options, so include 'FunValCheck', 'on' when constructing the options.
Sean de Wolski
el 25 de Sept. de 2012
fmincon's Interior-Point algorithm can recover from inf and nan output.
doc fmincon
8 comentarios
Sean de Wolski
el 1 de Oct. de 2012
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);
Ver también
Categorías
Más información sobre Solver Outputs and Iterative Display en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!