我正在尝试为 Matlab 中的函数编写 Taylor 求和,我实际上评估了 McLaurin making x=0
,在此符号之后在此代码中命名为 a :
这是我到目前为止尝试过的代码:
>> a = -100;
b = 100;
n = 20;
vectorx = linspace(a,b,50);
vectory = [];
sumterms = [];
syms x y a;
y = sin(a);
for i = 1:n
t = (diff(y,i-1) / factorial(i-1)) * (x-0)^(i-1);
sumterms = [sumterms;t];
end;
sumterms
for j = 1:length(vectorx)
x_i = vectorx(j);
aux = 0;
for k = 1:length(sumterms)
sumterm = sumterms(k);
aux = aux + subs(sumterm, [a,x], [0,x_i]);
end
vectory = [vectory;aux];
end
length(vectory)
length(vectorx)
plot (vectorx, vectory)
但是我没有得到正确的结果我已经跳过了每个句子,我真的看不出它有什么问题。
这是我对 sin(x) 的绘图结果:
这是指数(x)的图
每个 Sumterms 结果出现在图像捕获中,它们似乎没问题,但是我认为错误是在评估中。