how to set Fmincon function tolerance?

excuse me, how can set a tolerance? because I want to solve this:( I mean finding "cc" which satisfies below)

* F(cc)<=min(F(c))+(0.5)^k *

and k k is iteration number. in paper, writer said that tolerance should be set on "(0.5)^k". I've written a az below,but it doesn't work.

objfcn = @(cc) fmp_cmp(cc, Uu, Cc, mm, p);
options = optimset('TolFun', delta^k);
A = [ones(1,24)];
b =[1];
cc0=rand(1,24);
cc0=cc0/sum(cc0);
[cc,fval]=fmincon(objfcn,cc0,[],[],A,b,zeros(1,24),ones(1,24),options);

Answer 1

[Matt J]

Your 'options' should be the 10th input argument to fmincon. You have passed it as the 9th. If you have no nonlinear constraints, set that to [].

Answer 2

[Matt J]

There may be a typo in your post. The additive constant (0.5)^k doesn't affect the solution in any way. cc=argmin(F) is a solution for all k. Perhaps you really meant to write an equality,

    F(cc)=min(F)+(0.5)^k

If so, the first step is to minimize F() as I think you have done in other posts. If Fmin is this minimum, then the above reduces to

    Fmin+(0.5)^k -F(cc)=0

You could indeed use fminsearch to solve this nonlinear equality

   cc=fminsearch(@(cc) abs(Fmin+(0.5)^k -F(cc)) ,...)

However, it would be more direct to use fzero

    cc = fzero(@(cc) Fmin+(0.5)^k -F(cc) ,...)

[Maryamm Bodaghi]

Dear Star Strider, I've attached papers. algorithm 2, page 6, step(1) .

the algorithm is explained in file named "SP_10.1023-2FB-3AJOTA.0000037605.19435.63".

Dear Matr, realy I don't have Fmin, because I'm trying to get it, with YALMIP , but after 2 iteration I get NANs. also k is iteration of the main loop, not fmincon.

thanks

[Matt J]

realy I don't have Fmin, because I'm trying to get it,

But why can't you find Fmin with fmincon? And what effect is 0.5^k supposed to have? It just adds a constant to F(c) without changing the location of its minima.

[Maryamm Bodaghi]

thanks Dear Matt;

I'me trying to do your advice, but it doesn't accepr vector.

   it=1;  mm=1;  
    p=1e-9;   delta=0.5;   eta=10^(-6);

    c=sdpvar(1,24)
    const=[0<=c<=1,sum(c)==1];

function F_m_p_cmp=fmp_cmp(F_m_p_c,it,cc, Uu, Cc, mm, p)
for ii=1:mm
    for t=1:24
    ex(ii,t)=1*(cc(t)*((log(1-(1./(Uu(1,t,mm)+Cc(1,t,mm)))))+(log(1-(1./(Uu(2,t,mm)+Cc(2,t,mm)))))+...
    (log(1-(1./(Uu(2,t,mm)+Cc(2,t,mm)))))));
    end
end
expp=p*ex;
eexp=exp(-1*expp);
F_m_p_cmp=(1/p)*log10(sum(sum(eexp)));
fmp_cmp=F_m_p_c+(0.5)^it-F_m_p_cmp
end

    for ii=1:mm
        for t=1:24
        ex(ii,t)=-1*(c(t)*((log(1-(1./(Uu(1,t,mm)+Cc(1,t,mm)))))+(log(1-(1./(Uu(2,t,mm)+Cc(2,t,mm)))))+...
        (log(1-(1./(Uu(2,t,mm)+Cc(2,t,mm)))))));
        end
    end
    expp=p*ex;
    eexp=exp(expp);
    F_m_p_c=(1/p)*log(sum(sum(eexp)));

    solvesdp(const,-F_m_p_c)

    double(F_m_p_c)
    c=double(c)
    objfcn = @(cc) fmp_cmp(F_m_p_c,it,cc, Uu, Cc, mm, p)
    cc = fzero(objfcn,[0 1])

cc should be a 1x24 vector

[Matt J]

I'me trying to do your advice

Not really sure what "my advice" was. I don't think you've finished explaining what the goal is.