roots finding using muller method

Question

shiv gaur 2021-7-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/885404-roots-finding-using-muller-method

评论： Walter Roberson 2021-8-2

only one root is showing suppose f(x)=x^2+1 only -i root is showing not +i

here this is program of muller method

function shik3
 % t2=1e-6;
p0 = 10;
p1 = 20;
p2 = 30;
TOL = 10^-9;
N0 = 100;
format long;
h1 = p1 - p0;
h2 = p2 - p1;
DELTA1 = (f(p1) - f(p0))/h1;
DELTA2 = (f(p2) - f(p1))/h2;
d = (DELTA2 - DELTA1)/(h2 + h1);
i=3;
while i <= N0
      b = DELTA2 + h2*d;
      D = (b.^2 - 4*f(p2).*d).^(1/2);
      if abs(b-D) < abs(b+D)
          E = b + D;
      else
          E = b - D;
      end
      h = -2*f(p2)/E;
      p = p2 + h;
      
     
      if abs(h) < TOL
          p
          
         disp(p) 
          break
      end
      p0 = p1;
      p1 = p2;
      p2 = p;
      h1 = p1 - p0;
     h2 = p2 - p1;
     DELTA1 = (f(p1) - f(p0))/h1;
      DELTA2 = (f(p2) - f(p1))/h2;
     d = (DELTA2 - DELTA1)/(h2 + h1);
     i = i + 1;
     disp(p) 
     
     %disp([imag(p) t2]) 
     %plot(t2,real(p))
  
end
if i > N0
      formatSpec = string('The method failed after N0 iterations,N0= %d \n');
      fprintf(formatSpec,N0);
  end
 
function y=f(x)
   y=x^2+1 ;
 
   
    
    
end
end

2 个评论
显示无隐藏无

Jan 2021-7-25

Please edit your question. Do you see the pile of blank lines? Deleteing them and using the code block improves the readability. Thanks.

shiv gaur 2021-7-25

runing program only one root is showing

请先登录，再进行评论。

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

Answer 1

Chunru 2021-7-25

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/885404-roots-finding-using-muller-method#answer_753459

在 MATLAB Online 中打开

The Muller's method is supposed for finding the real solution. Your code is attempting to find the complex solution. I am wondering if it has any other side-effect.

The Muller's method is based on secant method (or in a similar way) to iteratively find the solution. It cannot guarantee to find all solutions of a functions. It just find a solution (normally real) within the specified interval.

For you code, you may find another solution by change the sign of the following line

% D = (b.^2 - 4*f(p2).*d).^(1/2);
D = -(b.^2 - 4*f(p2).*d).^(1/2);

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

shiv gaur 2021-7-25

one working at a time

请先登录，再进行评论。

Answer 2

Walter Roberson 2021-8-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/885404-roots-finding-using-muller-method#answer_758512

Your program only ever tries to find one root. You have a while loop, but that while loop is only dealing with one root.

You have several options:

As you are solving a polynomial with real-valued coefficients, all imaginary roots must occur as complex conjugates. So after you find one root, take its complex conjugate and test it to see whether it works as well
rewrite the code to try several different initial guesses, and keep trying guesses until you reach a sanity limit or you have found as many unique roots as you think you should have
rewrite your code to be working on a vector of guesses, and keep iterating until you reach a sanity limit or enough unique roots have reach tolerance. Note that if you fail to try enough different elements originally then you are not certain that enough of them will lead to unique roots.
use a completely different algorithm, such a symbolic method