Function Intersection using Newton's Method

Question

SB 2012-11-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/54149-function-intersection-using-newton-s-method

编辑： Torsten 2022-11-27

采纳的回答： Matt J

在 MATLAB Online 中打开

Hi everyone, I'm trying to write a function that finds a point at

which two functions f(x) and g(x) intersect, i.e., f(x) = g(x). I'm using

Newton's Method and making and h(x)=f(x)-g(x) and h'(x) as well, but I'm not getting the right x-value. Please help me debug my code!

 % function x = fgIntersect(f, df, g, dg, x0, tol, maxIter)
h=f(x0)-g(x0)
dh=df(x0)-dg(x0)
k=1;
while k<=maxIter
x=x0-h/dh;
if abs(x-x0)<tol*abs(x)
return
end
x0=x;
k=k+1;
end
end

2 个评论
显示无隐藏无

Matt J 2012-11-19

Give us example data that let's us reproduce the failure.

SB 2012-11-19

编辑：SB 2012-11-19

format compact; format long;

f = @(x) exp(x) - 3;

df = @(x) exp(x);

g = @(x) sqrt(x);

dg = @(x) .5*x^(-.5);

x = fgIntersect(f, df, g, dg, 1, 1e-6, 50)

x should equal 1.434542442506692

Another case:

format compact; format long;

p1 = [1 -2 3 -8];

p2 = [1 -3 2 -4];

f = @(x) polyval(p1,x);

df = @(x) polyval(polyder(p1),x);

g = @(x) polyval(p2,x);

dg = @(x) polyval(polyder(p2),x);

x = fgIntersect(f, df, g, dg, 2, 1e-6, 50)

x should be 1.561552842846145

请先登录，再进行评论。

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

Answer 1

Matt J 2012-11-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/54149-function-intersection-using-newton-s-method#answer_65825

You're not updating h and dh within your loop.

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

SB 2012-11-19

Thank you so much! I got it, I can't believe I forgot to update them within my loop.

Torsten 2022-11-26

编辑：Torsten 2022-11-27

在 MATLAB Online 中打开

f = @(x) exp(x) - 3;
df = @(x) exp(x);
g = @(x) sqrt(x);
dg = @(x) .5*x^(-.5);
x = fgIntersect(f, df, g, dg, 1, 1e-6, 50)
x = 1.4345
p1 = [1 -2 3 -8];
p2 = [1 -3 2 -4];
f = @(x) polyval(p1,x);
df = @(x) polyval(polyder(p1),x);
g = @(x) polyval(p2,x);
dg = @(x) polyval(polyder(p2),x);
x = fgIntersect(f, df, g, dg, 2, 1e-6, 50)
x = 1.5616
function x = fgIntersect(f, df, g, dg, x0, tol, maxIter)
  h=f(x0)-g(x0);
  dh=df(x0)-dg(x0);
  k=1;
  while k<=maxIter
    x=x0-h/dh;
    if abs(x-x0)<tol*abs(x)
        return
    end
    x0=x;
    h = f(x0)-g(x0);
    dh = df(x0)-dg(x0);
    k=k+1;
  end
end