Error while trying to code the Newton Method for a System of Equations

Question

Roberto 2012-7-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/43564-error-while-trying-to-code-the-newton-method-for-a-system-of-equations

Hi there. This is my first post in the forum. My name is Roberto and I am from Spain. Since English is not my home language, please do not be too hard to me, thanks.

Well, I coded the Newton Method for a System of Equations. It looks like:

function [x,fx,iters]=sistemanewton_sept(func,jacobi,x0,TOL,maxit)
    %x0 vector
    mitol=2*TOL; iters=1;
    x0=x0(:); %Lo ponemos como vector columna
    while (mitol>TOL && iters<maxit)
       J=jacobi(x0); f=func(x0); en=-J\f; %¿sale vector columna directamente?
       x=x0+en; x0=x;
       mitol=abs(func(x0));
       iters=iters+1;
    end
    fx=func(x0);
    if mitol>TOL
        disp('ADVERTENCIA: No se pudo alcanzar la tolerancia deseada dentro del numero de iteraciones permitidas');
    end
end

I don't know how to color the code. The error I get is:

??? Operands to the || and && operators must be convertible to logical scalar values.
Error in ==> sistemanewton_sept at 10
    end

The functions I used in my code are: func and jacobi. They look like:

function f=func(x)
x1=x(1); x2=x(2);
f=[sqrt(x2*x2*x2/5), 0.25*(sin(x1)+cos(x2))]';

And:

function J=jacobi(x)
x1=x(1); x2=x(2);
J=[0, sqrt(0.45*x2); 0.25*cos(x1), -0.25*sin(x2)];

That error is driving me crazy, since I never got it before. As I could read in some other post, that error is related to the way I use the "&&" operator, which must be only used to compare scalars. I am already comparing Scalars so...

Thank you so much.

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

Star Strider 2012-7-13

编辑：Star Strider 2012-7-13

在 MATLAB Online 中打开

Your code looks good to me, but I may be missing something. To start troubleshooting, I suggest you put before your ‘while’ loop something like:

TestTol = mitol>TOL
TestItr = iters<maxit
TestClass = class([TestTol  TestItr])

The variables ‘TestTol’ and ‘TestItr’ should both evaluate as ‘logical’. If they do not, then you have to find and correct the problem. If they are both ‘logical’ and you continue to get the error, I suggest adding parentheses to your ‘while’ statement:

while (mitol>TOL) && (iters<maxit)

or:

while ((mitol>TOL) && (iters<maxit))

Since ‘func’ and ‘jacobi’ do not use the logical operators, they are probably not the problem.

请先登录，再进行评论。

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

Answer 1

Kye Taylor 2012-7-13

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/43564-error-while-trying-to-code-the-newton-method-for-a-system-of-equations#answer_53450

The first time you enter the while loop, you are comparing scalars (assuming the input TOL is a scalar). However, when you overwrite mitol with the absolute value of the output from func, mitol becomes nonscalar since func outputs a 2-by-1 vector.

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

Roberto 2012-7-13

编辑：Roberto 2012-7-13

Oh... you are right, I forgot f(x0) was a vector... I replaced that line with this new one: mitol=norm(func(x0));

Thank you Kye Taylor, very helpful.

请先登录，再进行评论。

Error while trying to code the Newton Method for a System of Equations

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

回答（1 个）

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

另请参阅

类别

标签

Community Treasure Hunt

Error while trying to code the Newton Method for a System of Equations

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

回答（1 个）

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

另请参阅

类别

标签

Community Treasure Hunt

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

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