原创力文档插值法
下列数据点的插值
x 0 1 4 9 16 25 36 49 64
y 0 1 2 3 4 5 6 7 8
可以得到平方根函数的近似，在区间[0,64]上作图.
(1)用这9个点作8次多项式插值Ls(x).
(2)用三次样条(第一边界条件)程序求S(x).
从得到结果看在[0,64]上，哪个插值更精确；在区间[0,1]上，两种插值哪个更精确？
解：(1)拉格朗日插值多项式，求解程序如下
syms x l;
x1=[0 1 4 9 16 25 36 49 64];
y1=[0 1 2 3 4 5 6 7 8];
n=length(x1);
Ls=sym(0);
for i=1:n
l=sym(y1(i));
for k=1:i-1
l=l*(x-x1(k))/(x1(i)-x1(k));
end
for k=i+1:n
l=l*(x-x1(k))/(x1(i)-x1(k));
end
Ls=Ls+l;
end
Ls=simplify(Ls) %为所求插值多项式Ls(x).
输出结果为
Ls =
63504000*x^2+95549/72072*x-1/3048192000*x^8-2168879/4
文档评论（0）