Help with a minimize problem

Question

Roberto López 2017-3-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/328123-help-with-a-minimize-problem

编辑： Walter Roberson 2017-3-6

f

unction b = two_var(v)
x = v(1);
y = v(2);
b = 242*log(1-x-y)+120*log(2*x-x^2-2*x*y)+79*log(2*y-y^2-2*x*y)+33*log(2*x*y);
end
v = [0.5,0.5]
a = fminsearch(@two_var,v);

I'm trying to minimize this function with the next code. The answer is x=0.247 and y=0.173

I leave the link where the problem comes from. Page 4

https://ocw.mit.edu/courses/mathematics/18-443-statistics-for-applications-fall-2006/lecture-notes/lecture12.pdf

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

John D'Errico 2017-3-5

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/328123-help-with-a-minimize-problem#answer_257350

编辑：John D'Errico 2017-3-5

在 MATLAB Online 中打开

Time to think about your function. ALWAYS do that. Plot it if possible.

fun = @(x,y) 242*log(1-x-y)+120*log(2*x-x^2-2*x*y)+79*log(2*y-y^2-2*x*y)+33*log(2*x*y);
fun(.5,.5)
ans =
  -Inf
fun(-1,-1)
ans =
           609.016175390547 +      625.176938064369i
fun(.247, .173)
ans =
       -455.717889710442

So some values of x and y generate -inf. Some generate complex numbers. That is completely expected, since the log function does nasty stuff at 0 or negative values. Well, nasty in terms of what fminsearch will expect.

Basic rule: fminsearch expects a continuous, real valued function of the inputs.

If that presumption fails, then expect all hell to break loose in the eyes of fminsearch.

ezsurf(fun)

So, what do we see? First, the function is certainly unbounded from below, going to -inf along the line

x + y = 1

Next, it appears to have a MAXIMUM at the location you describe.

So honestly, I think you are confused. Are you trying to MAXIMIZE this function?

fminseach is a MINIMIZATION tool. You can make it maximize by negating the function as returned. It still minimizes, but the negative of your original function, so a maximum.

fun2 = @(xy) -fun(xy(1),xy(2));
[xymax,fmax] = fminsearch(fun2,[.2 .2])
xymax =
         0.246457832567394          0.17315985403955
fmax =
          455.717447610379

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Roberto López 2017-3-6

Thank you sincerly, I was completely confused. It's the solution I was looking for.

请先登录，再进行评论。

Answer 2

kowshik Thopalli 2017-3-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/328123-help-with-a-minimize-problem#answer_257312

在 MATLAB Online 中打开

I dont understand.

    two_var =@(x) 242*log(1-x(1)-x(2))+120*log(2*x(1)-x(1)^2-2*x(1)*x(2))+79*log(2*x(2)-x(2)^2-2*x(1)*x(2))+33*log(2*x(1)*x(2));
v = [0.5,0.5];
x = fminsearch(fun,v)

this gives me x =

0.525000000000000 0.500000000000000

You forgot to attach the link. Changing v would give a different answer. What is it that you exactly want to do?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Roberto López 2017-3-5

https://ocw.mit.edu/courses/mathematics/18-443-statistics-for-applications-fall-2006/lecture-notes/lecture12.pdf

Sorry this is the link. Is for finding numerically the Maximum likelihood estimator for a chi square test. The sum of p+q+r= 1, r is described by p and q as r=1-p-q. The solution =.5 and 0.52 is incorrect because they must sum less than 1. With other values for example [.3, .2] The function does not work even if I increase MaxFunEvals option as Matlab suggest. Do you know how can I solve it?

请先登录，再进行评论。

Answer 3

Roger Stafford 2017-3-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/328123-help-with-a-minimize-problem#answer_257321

It is quite possible for functions such as yours to have more than one point of local minimum. If that is the case, the one that fminsearch will arrive at will depend on the given initial estimate. You need to try a different initial estimate if you are to arrive at the value you refer to in the MIT article.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Roberto López 2017-3-5

Yeah, I've tried with values close to that given by the article but message appears about MaxFunEvals must be increase even If I do it. It's still not working, maybe the problem looks to be singular but I don't know how these values in the article were obtain.

请先登录，再进行评论。