你好同学,不需要每次求导的!!而且符号运算只是为了推导公式,不需要全程出现在关键的程序代码里。
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