Is there any other alternative fuction that I can use instead of fminunc?.
13 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
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 Commenti
Risposte (3)
Matt Kindig
il 20 Set 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
Walter Roberson
il 23 Set 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 Commenti
Walter Roberson
il 25 Set 2012
You appear to already be constructing options, so include 'FunValCheck', 'on' when constructing the options.
Sean de Wolski
il 25 Set 2012
fmincon's Interior-Point algorithm can recover from inf and nan output.
doc fmincon
8 Commenti
Sean de Wolski
il 1 Ott 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);
Vedere anche
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!