我在使用定点方法创建用于查找函数根作为输入的代码时遇到了麻烦,
我在这里使用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进行处理。 我该怎么办?