(matlab)分段线性插值代码

%y=1/(1+x^2)分段线性xi=-5+10*i/n(i=1,1,...,n)。比较发现，随着n的增大，两者吻合得越来越好，龙格现象并未发生clearall%清除命令空间中所有变量%输入插值区间的等分数disp('给出插值区间的等分数n')n=input('n=');%作被插函数图象u=-5:(10/200):5;v=1./(1+u.^2);plot(u,v,'r');holdon%固化图形屏幕%给出插值条件x=-5:(10/n):5;y=1./(1+x.^2);%取需用分段线性插值函数计算其函数值的点zz=-5:10/(2*n):5;m=max(size(z));%计算分段线性插值函数在z处的值，并作图fork=1:mi=1;whilei=nifz(k)=x(i)&z(k)=x(i+1)Ih(k)=y(i)*(z(k)-x(i+1))/(x(i)-x(i+1))+y(i+1)*(z(k)-x(i))/(x(i+1)-x(i));breakelsei=i+1;endendendplot(z,Ih,'b')holdoff%释放固化的图形屏幕