研一上学期,学校开了数值分析这门课,而且要求用matlab进.docVIP

研一上学期，学校开了数值分析这门课，而且要求用matlab进行上机试验，经过一段时间的学习，对matlab有了最初步的了解，结合所学课程，废了好大的劲才编了几个小程序，希望对后来学习的人有点帮助。 1．题目/要求： 对Runge函数R(x)=1/1+x2, x∈[-5，5]分别用下列方法做插值函数，并与R(x)的图形进行比较。 （1）取等距节点xi=-5+i (i=0,1,……,10), 绘其10次Newton插值多项式。 （2）取节点(i=0,1,……,20) 绘其20次langrange插值多项式。 （3）取等距节点xi=-5+i (i=0,1,……,10), 绘其分段线性插值函数。 （4）取等距节点xi=-5+i (i=0,1,……,10), 绘其3次自然样条插值函数。 ? 2.算法程序 ① 牛顿插值 function yy=newton(x,y,xx) % yy=newton(x,y,xx)是求牛顿插值，其中x是插值节点，y是对应的向量，xx %是插值点，也可以是向量 disp(牛顿插值); x=-5:0.2:5 y=1./(1+x.^2) plot(x,y,r); hold on; xx=-5+[0:1:10] n=length(x); m=length(y); if m~=n error(向量x和向量y的长度必须一致) end w=ones(1,length(