MatLab定点法找到函数的根作为输入

时间:2016-03-28 05:26:01

标签: matlab fixed-point

我在使用定点方法创建用于查找函数根作为输入的代码时遇到了麻烦,

我在这里使用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进行处理。 我该怎么办?

0 个答案:

没有答案