0

我很难创建一个代码来查找函数的根作为定点方法的输入,

在这里,我使用 Newton-Raphson 方法完成了它:

clc,close all
syms x;
fprintf('Newton Raphson\n');
Fun = input('\nType a function \n');
x0 = input('\nType initial value \n');
f = sym(Fun);
df = diff(f,x);
while (1)
   a = subs(f, 'x', x0);
   b = subs(df, 'x', x0);
   x1 = x0 - a/b
   er = (abs((x1 - x0)/x1))*100
   if ( er <= 0.05)
      break;
   end
   x0 = x1;
end

这是一个使用定点方法的代码,但在一个不作为输入的固定函数上:

clc,close all
x0 = 0.5
while (1)
   x1 = (exp(-x0) - sin(x0)) / 5
   er = (abs((x1 - x0)/x1))*100
   if ( er <= 0.05)
      break;
   end
x0 = x1;

结尾

但是我不能使用函数作为输入来做到这一点,因为定点方法在左侧对变量 x 进行了分割。我该怎么做?

4

0 回答 0