Can you help me with the output of my optimizer?

I have this program from

par0 = [0.1,0.0001,100];% initial values
history = [];
fun_objetivo = @(par,pfrac)FunObjetivo(par);
options = optimset('MaxIter',3000,'MaxFunEvals',3000,'Display','iter','OutputFcn', @myoutput);
%We use here the Nelder-Mead method

par_optimos = fminsearch(fun_objetivo, par0, options);

 function stop = myoutput(par_optimos,optimvalues,state);
        stop = false;
        if isequal(state,'iter')
          history = [history;par_optimos];
        end
    end

but I can´t make it give me the value of my vector par_optimos for each iteration, can someone help me?

Answer 1

[Sulaymon Eshkabilov]

You'd need define a math formulation of fun_objetivo = @(par,pfrac)FunObjetivo(par);

Once it is done, you can run via loop iteration, i.e.: for ... end

Answer 2

[Alan Weiss]

Use this example as a guide: Example of a Nested Output Function.

Alan Weiss

MATLAB mathematical toolbox documentation

[Jorge Armando Vazquez]

That is what I used

fun_objetivo = @(par,pfrac)FunObjetivo(par);
history=[];
options = optimset('MaxIter',3000,'MaxFunEvals',3000,'Display','iter','OutputFcn', @myoutput);;%,'PlotFcns',@optimplotfval);
%We use here the Nelder-Mead method
par_optimos = fminsearch(fun_objetivo, par0, options);

   function stop = myoutput(par_optimos,optimvalues,state);
        stop = false;
        if isequal(state,'iter')
          history = [history;par_optimos];
          
        end
    end 

but there is error with

Error in Main_optim_15_05_21>myoutput (line 15)
          history = [history;par_optimos];

[Walter Roberson]

You did not nest the function like suggested.

history = do_the_work;

function hist = do_the_work()
    fun_objetivo = @(par,pfrac)FunObjetivo(par);
    history=[];
    options = optimset('MaxIter',3000,'MaxFunEvals',3000,'Display','iter','OutputFcn', @myoutput);;%,'PlotFcns',@optimplotfval);
    %We use here the Nelder-Mead method
    par_optimos = fminsearch(fun_objetivo, par0, options);
    hist = history;
    
   function stop = myoutput(par_optimos,optimvalues,state);
        stop = false;
        if isequal(state,'iter')
          history = [history;par_optimos];
          
        end
   end
end