计算物理作业.docVIP

计算物理作业1 浦正国 2013202020010 已知函数节点 、计算和拉格朗日插值多项式的值， 、求带有自然边界条件的三次样条函数表达式中的参数计算 、画出的图形（各描41点） 运行程序（Matlab） x0=-1:.2:1;y0=1./(1+25.*(x0.^2));%列出已知点x0.y0 syms x l %求拉格朗日插值基 for j=1:11 ll=1; pp=1; for i=1:11 if i==j continue; end ll=(x-x0(i))*ll; pp=(x0(j)-x0(i))*pp; end l(j)=ll/pp; end l(11) y=0; for i=1:11 y=y0(i)*l(i)+y; %拉格朗日差值函数 end disp(y); x1=-1:.05:1; y1=subs(y,x,x1); figure(1) subplot(1,2,1) x2=-1:.05:1; y2=1./(1+25.*(x2.^2));y3=subs(l(11),x,x1); %用x2,y2表示被插函数 plot(x1,y1,-.