matlab牛顿迭代法求根时出错

Question

使用matlab牛顿迭代法求根时出现错误，这段代码如何修改呢？最后的结果要求精确。

Answer 1

你好同学，不需要每次求导的！！而且符号运算只是为了推导公式，不需要全程出现在关键的程序代码里。

clc;clear
% 前面用符号计算导数
syms t 
F = 0.99403+1.671e-4*t+9.7215e-8*t^2-9.5838e-11*t^3+1.952e-14*t^4-1.2;
f = diff(F);
% 这里将符号转化成函数
F = matlabFunction(F);
f = matlabFunction(f);
% 之后跟符号运算没关系了
error = 1e-8;
x = 1; 
n = 0;
dx = F(x)/f(x);
while(abs(dx)>error)
    n=n+1;
    x = x - dx;
    dx = F(x)/f(x);
end
disp(x)
disp(n)
disp(abs(dx))

得到结果：

     1.126009750841875e+03

     5

     2.087453837960381e-13