数据插值，拉格朗日，牛顿，埃尔米特插值，等距节点插值三次样条插值

数据插值，拉格朗日，牛顿，埃尔米特插值，等距节点插值三次样条插值三次样条插值与拉格朗日的比较

function f=Atken（xyx0）
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（f6）;  %将插值多项式的系数化成6位精度的小数
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-06-07 15:04  4\
     文件         527  2016-06-22 16:44  4\Atken.m
     文件         809  2016-06-22 16:43  4\Hermite.m
     文件         927  2016-06-22 16:36  4\Language.m
     文件         905  2016-06-22 16:08  4\M4_1.m
     文件         511  2016-06-22 16:33  4\M4_10.m
     文件         332  2016-06-22 16:33  4\M4_11.m
     文件         161  2016-06-22 16:34  4\M4_12.m
     文件         256  2016-06-22 16:35  4\M4_13.m
     文件         706  2016-06-22 16:38  4\M4_15.m
     文件         335  2016-06-22 16:42  4\M4_18.m
     文件         254  2016-06-22 16:43  4\M4_19.m
     文件         526  2016-06-22 16:08  4\M4_2.m
     文件         862  2016-06-22 16:44  4\M4_21.m
     文件        1729  2016-06-22 16:45  4\M4_22.m
     文件         282  2016-06-22 16:09  4\M4_3.m
     文件         278  2016-06-22 16:10  4\M4_4.m
     文件        1200  2016-06-22 16:11  4\M4_5.m
     文件         711  2016-06-22 16:12  4\M4_6.m
     文件         341  2016-06-22 16:12  4\M4_7.m
     文件         663  2016-06-22 16:14  4\M4_8.m
     文件         248  2016-06-22 16:32  4\M4_9.m
     文件         677  2016-06-22 16:37  4\Newton.m
     文件         812  2016-06-22 16:40  4\Newtonback.m
     文件         822  2016-06-22 16:39  4\Newtonforward.m