N=10; %%定义插值节点的个数
x=-5:10/N:5; %%依据课本xi=-5+i*h(i=0,1,...N),h=10/N
y=x./(1+x.^4); %%依据课本(1)插值公式
xi=-5:0.25:5; %%依据课本xk=-5+0.25k(k=0,1,...40)
m=length(x); %%求出插值节点向量长度
n=length(y);
s=zeros(1,length(xi));
if m ~= n , error('向量x与y的长度必须一致');
%%这里肯定一致,只是为了消除直接选取x,y数据的时候出错而设定的
end;
%%这里定义一个长度为计算节点个数的单位向量,可以加快运行速度
for m=1:length(xi)
for b=1:n
if (xi(m)>x(b))&&(xi(m)<=x(b+1))
%%寻找取值区间
s(m)=(xi(m)-x(b))*(y(b)-y(b+1))/(x(b)-x(b+1))+y(b);
%%两点式直线公式
else
continue
end
end
end
fprintf('\t\t\t插值点X\t\t插值点Y\n');
for i=1:length(x)
fprintf('I(%2.0f)\t\t%5.3f\t\t%5.3f\n',i,x(i),y(i));
end
fprintf('\t\t\t计算点X\t\t计算点Y\n');
for i=1:length(s)
fprintf('I(%2.0f)\t\t%5.3f\t\t%5.3f\n',i,xi(i),s(i));
end
hold on
%以下是设置绘图区域的相关指令
plot(x,y,'r*');
ezplot('x./(1+x.^4)',[-5,5]);
plot(xi,s,'k-');
plot(xi,s,'go');
title('分段线性插值');
legend('插值点','原始函数','分段线性插值函数','计算点');
xlabel('x的取值');
ylabel('y的取值');
axis([-5,5,min([min(y) min(f)])-0.2,max([max(y) max(f)])+0.2]);
set(gca,'xtick',[-5:0.25:5]);
set(gca,'ytick',[-10:0.1:10]);
grid on
hold off