Supress Warning in a Loop

Question

ercan duzgun 2021-2-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/754969-supress-warning-in-a-loop

评论： ercan duzgun 2021-2-24

I have a matrix table. For each row in the table, I am using Newton-Raphson iteration. And for some variables of the table, the iteration doesn't converge, and it gives warning errors. How can I supress these warning messages?

I receive this kind of message in command window :

Warning: Matrix is singular,
close to singular or badly
scaled. Results may be
inaccurate. RCOND = NaN. 
> In
SP_FK_NewtonRaphson_function
(line 8)
  In fonk_cemberYont_FK (line
  72)

My Newton-Raphson function file is:

  function [X,i]=SP_FK_NewtonRaphson_function(verilen_L,X0)
maxIter=1000; % maximum iteration
tolX=1e-6; % Tolerance for error
X=X0;
Xold=X0;
for i=1:maxIter
     [f,j]=myfunction2(verilen_L,X);
     X=X-inv(j)*f;
     
     err(:,i)=abs(X-Xold);
     Xold=X;
     Tablo(:,i)=X;
     if (err(:,i)<tolX)
         break;
     end
end
  end

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

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2021-2-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/754969-supress-warning-in-a-loop#answer_632169

在 MATLAB Online 中打开

X=X-inv(j)*f;

replace with

if rank(j)<length(j) 
  X = nan;
else
  X = X - j\f;
end

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

ercan duzgun 2021-2-24

在 MATLAB Online 中打开

If I write as you say, I have "error" message, in state of "warning " message.

I receive this eror message:

Index exceeds the number of array elements (1).
Error in myfunction2 (line 4)
ty=X(2);
Error in SP_FK_NewtonRaphson_function (line 8)
     [f,j]=myfunction2(verilen_L,X);
Error in fonk_cemberYont_FK (line 72)
[Tablo_NR_FK_xyzABG(s,:),Tablo_NR_iter(s,:)]=SP_FK_NewtonRaphson_function(Tablo_L(:,1),[G_0(1:3,s);ABG_0(1:3,s)]);

But if I write my previous command, it execute the iteration well, however there are lots of warning messages for several iterations on the screen. Displaying the warning messages each time on sceen seems to be time consuming. Therefore I wanted to supress the warning messages.

@Walter Roberson, If I change it with your suggested code, the program stop with that error message.

And I want to ask isn't it very time consuming for the program to check your "if else" command. Because it is in a function file, and function file is in a loop for 1000 iterations. Checking the condition of X with if-else in a 1000 loop cycle isn't bad for the program? I am trying to learn, and know your comment. Thanks in advance.

Walter Roberson 2021-2-24

在 MATLAB Online 中打开

function [X,i]=SP_FK_NewtonRaphson_function(verilen_L,X0)
    maxIter=1000; % maximum iteration
    tolX=1e-6; % Tolerance for error
    X=X0;
    Xold=X0;
    for i=1:maxIter
        [f,j]=myfunction2(verilen_L,X);
        if rank(j)<length(j) 
            X = nan(size(X0));
        else
            X = X - j\f;
        end     
        err(:,i)=abs(X-Xold);
        Xold=X;
        Tablo(:,i)=X;
        if (err(:,i)<tolX)
         break;
        end
    end
end