拉格朗日插值 function 图在节点大于 28 时出现偏差

Question

我不知道这是不是问这个问题的正确地方。 但我有以下代码基于节点 x 插入 function y ，它基于termNumber 。

termNumber = 28;
x = (2.*(0:termNumber)./termNumber)-1;
y = exp(3.*x).*cos(3.*x);

sum=0;
for i=1:length(x)
    p=1;
    for j=1:length(x)
        if j~=i
            c = poly(x(j))/(x(i)-x(j));
            p = conv(p,c);
        end
    end
    term = p*y(i);
    sum= sum + term;
end
disp(sum);
K=flip(sum);

for  g=length(x):-1:2 %display the coefficient on polynomial form
    fprintf('+(%.2fx^%d)',K(g),g-1);
end

fprintf('+(%.2f)',K(1));

v=linspace(-1,1,100);
b = polyval(sum,v);
b1 = exp(3.*v).*cos(3.*v);
plot(v,b,'o') %display lagrange polynomial
hold on
plot(v,b1) %display actual function

当termNumber小于 28 时，使用上述方法正确显示图形：见下图：当 termNumber = 28 时的图形，其中红色曲线是实际函数（蓝色是拉格朗日插值多项式）。 但是当termNumber >= 29 时，两条曲线开始分开： termNumber = 29 。 当termNumber变大时，两条曲线不再可重新定义 termNumber = 40 。 由于我根据教科书编写了该方法，所以我似乎无法弄清楚我做错了什么。 任何输入表示赞赏。