increasing fsolve accuracy when fval values are different in terms of order of magnitude

Question

Sayed 2015-6-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/225212-increasing-fsolve-accuracy-when-fval-values-are-different-in-terms-of-order-of-magnitude

评论： Walter Roberson 2015-6-29

I have a set of six equations-six unknowns. When I solve it by fsolve, it "stops prematurely" and gives me 6 values for unknowns that are not exactly the answers that I'm expecting (I've tried different initial guesses, but it's the same). Answer set that fval gives is:

[x1=6.242, x2=2.148, x3=2.688, x4=1.563, x5=-3.359, x6=0.996]

and values of "fval" after solving are:

[-9e-9, 1e-10, -2e-6, 0.0054, -0.1208, -0.1615 ]

(The last 3 fvals are comparatively larger). Please let me know what I can do to enhance the results accuracy? Thanks

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

请先登录，再进行评论。

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

Answer 1

John D'Errico 2015-6-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/225212-increasing-fsolve-accuracy-when-fval-values-are-different-in-terms-of-order-of-magnitude#answer_183920

编辑：John D'Errico 2015-6-25

The simple answer is the one you need to use.

SCALE YOUR VARIABLES. Scale them so they are all roughly the same order of magnitude. A factor of 10 is irrelevant. A dynamic range of 1e10 is a problem.

Then after the solve is done, scale them back.

6 个评论
显示 4更早的评论隐藏 4更早的评论

Walter Roberson 2015-6-26

Why shouldn't they be very different? I showed a case where they must be very different. It is what is demanded by the system of equations.

Questions I would consider would include:

was a mistake made in entering the equations?
if the "known" answers are plugged into the equations then do the expected results get produced? This is a crosscheck onto whether the MATLAB expression has been implemented properly
are there multiple solutions to the equations?
when the solutions found are plugged into the MATLAB expression of the equations, do the equations hold to within tolerance? This tests whether MATLAB is doing a reasonable job with the function that it has been given
when the solutions found are plugged into the equations on a different system, do the equations hold to within tolerance? This is a crosscheck on whether the equations have been expressed properly in MATLAB

What cause does fsolve() give when it stops "prematurely" ?

And I do not understand what you mean by drawing a difference between "Answer set that fval gives" and "values of "fval" after solving" ? What is "giving" the "Answer set" ? What is "solving" to give those other numbers?

Sayed 2015-6-29

Dear Walter,

Thanks for commenting.

When we use fsolve, it returns two different things:

1- an answer set: that depending to your equations have different components and of course could be very different in their order of magnitude as well.

2- fval values that is the values of your equations function calculated with the final answer that Matlab has found (F==0). So, fval can't basically have very different values and all its components should be very close to zero; while in my case only the first 3 component are close to zero, but the rest are not [-9e-9, 1e-10, -2e-6, 0.0054, -0.1208, -0.1615 ]

This is fsolve message when stops:

Solver stopped prematurely.

fsolve stopped because it exceeded the function evaluation limit, options.MaxFunEvals = 9000 (the selected value).

Walter Roberson 2015-6-29

Pass in a larger options.MaxFunEvals

请先登录，再进行评论。

increasing fsolve accuracy when fval values are different in terms of order of magnitude

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

回答（1 个）

6 个评论
显示 4更早的评论隐藏 4更早的评论

另请参阅

类别

标签

Community Treasure Hunt

increasing fsolve accuracy when fval values are different in terms of order of magnitude

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

回答（1 个）

6 个评论 显示 4更早的评论隐藏 4更早的评论

另请参阅

类别

标签

Community Treasure Hunt

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

6 个评论
显示 4更早的评论隐藏 4更早的评论