分段线性插值

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