x=02 04 06 08 10;y=098 092 081 064 038;plotxyhold onfor i=115yi=098-03Xi-02-062500xi-02Xi-04-020833xi-02 xi-04xi-06-052083xi-02xi-04xi-06xi-08endk= 0 1 10 11x0=02+008kfor i=114y0i=098-03xi-02-062500xi

x=[0.2 0.4 0.6 0.8 1.0]; % 定义数据的横坐标 y=[0.98 0.92 0.81 0.64 0.38]; % 定义数据的纵坐标 plot(x,y) % 画出原始数据的图像 hold on % 保持图像，方便后续添加内容

% 使用4次牛顿插值公式计算出每个插值点的纵坐标 for i=1:1:5 y(i)=0.98-0.3*(X(i)-0.2)-0.62500*(x(i)-0.2)(X(i)-0.4)-0.20833(x(i)-0.2)* (x(i)-0.4)(x(i)-0.6)-0.52083(x(i)-0.2)(x(i)-0.4)(x(i)-0.6)*(x(i)-0.8) end

k=[ 0 1 10 11] % 定义插值点的横坐标 x0=0.2+0.08*k % 计算出插值点的横坐标

% 使用4次牛顿插值公式计算出每个插值点的纵坐标 for i=1:1:4 y0(i)=0.98-0.3*(x(i)-0.2)-0.62500*(x(i)-0.2)(x(i)-0.4)-0.20833(x(i)-0.2)(x(i)-0.4)(x(i)-0.6)-0.52083*(x(i)-0.2)(x(i)-0.4)(x(i)-0.6)* (x(i)-0.8) end

plot ( x0,y0, 'o',x0,y0 ) % 画出4次牛顿插值的图像 hold on % 保持图像，方便后续添加内容

y1=spline(x,y,x0) % 使用三次样条插值计算出插值点的纵坐标 plot (x0,y1, 'o') % 画出三次样条插值的图像 hold on % 保持图像，方便后续添加内容

s=csape(x,y,'variational') % 使用三次样条插值自然边界计算出插值点的纵坐标 fnplt(s,'r') % 画出三次样条插值自然边界的图像，颜色为红色 hold on % 保持图像，方便后续添加内容

gtext('三次样条自然边界') % 添加注释 gtext('原图像') % 添加注释 gtext('4次牛顿插值') % 添加注