How do I improve solving 18 second order system of algebraic equations in fsolve ?

1 次查看（过去 30 天）

farzad 2020-3-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/509491-how-do-i-improve-solving-18-second-order-system-of-algebraic-equations-in-fsolve

评论： Walter Roberson 2020-3-8

Hi All

trying to solve 18 nonlinear (2nd order)equations using fsolve I get : possibility of inacuaracy in the results. while I have used the options like : 1e6 iterations plus 1e-5 solution tolerance.

this is the options setting :

options = optimoptions(@fsolve,'Display','iter','MaxFunEvals',1e6,'MaxIter',1e6,'TolFun',1e-6,'DerivativeCheck','on','Diagnostics','on');

and this is the error I get

Equation solved, inaccuracy possible.
The vector of function values is near zero, as measured by the default value
of the function tolerance. However, the last step was ineffective.
<stopping criteria details>

is there a way to improve it ?

10 个评论
显示 8更早的评论隐藏 8更早的评论

Walter Roberson 2020-3-7

Do you have the Symbolic Toolbox? Sometimes vpasolve() is better suited when accuracy is a problem.

Anyhow, you should make the solution tolerence smaller.

In some cases, when the coefficients are large enough, the gradient can end up so steep that you have to switch to absolute tolerences instead of relative tolerances.

Remember that a warning of inaccuracy does not mean that the solution returned is too inaccurate for your purposes -- it just means that when the solution was cross-checked, the residue was larger than a tolerance, so it is warning that the solution might not be so good. But if you are asking for tolerences of 1e-5 on a system with 18 variables, then chances are that you are not all that concerned about exact locations.

In the context of polynomials of degree 2, the warning would be valid for the situation where the polynomial has a minima or maxima that is close to zero rather than a zero crossing. The equation (x+a)^2+delta with positive delta could potentially be determined to have a root at -a if delta is close enouugh to zero even though the mathematical root is complex valued.

farzad 2020-3-7

dear Walter, well with the current inaccuracy, my results are not super good and cause large errors in the next post processing on the results and give bigger deviations later. so I would need a close enough result. I thought 10^-5 is the difference of calculated with respect to the expected. I reduced it but nothing chaged. i need to make an executable later, so I dont know if vpasolve could be useful