atken.m

来自「数值分析中几种常用的插值算法的matlab源程序。」· M 代码 · 共 26 行

function f = Atken(x,y,x0)
syms t;
if(length(x) == length(y))
    n = length(x);    
else
    disp('x和y的维数不相等！');
    return;
end                                      %检错

y1(1:n) = t;                             %符号函数数组要赋初值
for(i=1:n-1)   
    for(j=i+1:n)
        y1(j) = y(j)*(t-x(i))/(x(j)-x(i))+y(i)*(t-x(j))/(x(i)-x(j));    
    end
    y = y1;
    simplify(y1);
end

if(nargin == 3)
    f = subs(y1(n),'t',x0);            %计算插值点的函数值
else
    simplify(y1(n));                   %化简
    f = collect(y1(n));                %将插值多项式展开
    f = vpa(f,6);                      %将插值多项式的系数化成6位精度的小数
end