mathanalysis302.m

来自「数值分析的实验程序」· M 代码 · 共 18 行

N=10;
i=0:N;
x=-1+2/N*i;
fx=1./(1+25*(x.^2));
figure;
plot(x,fx,'r*');
grid on;
hold on;
xx=-1:0.01:1;
del=0;
ffx=zeros(1,length(xx));
for i=1:N
    del=(fx(i+1)-fx(i))/(x(i+1)-x(i));
    F=find(xx>=x(i)&xx<x(i+1));
    ffx(F)=(xx(F)-x(i))*del+fx(i);    
end
plot(xx,ffx);
title('分段线性插值');