上海交通大学计算方法作业答案【DOC精选】.docxVIP

P50 – 1 %%牛顿插值多项式function [ c, d] = newpoly( x,y )%这里 x为n个节点的横坐标所组成的向量，y为纵坐标所组成的向量。 %c为所求的牛顿插值多项式的系数构成的向量。n=length(x);d=zeros(n, n); d(: , 1)=y;for j=2 : nfor k=j : n d(k, j)=(d(k, j-1) - d(k-1, j-1)) / (x(k)-x(k-j+1));endend c =d(n, n);for k=(n-1) : - 1 : 1 c =conv(c, poly(x(k))); m=length(c); c(m)=c(m)+d(k, k);end x = 0.2:0.2:1 ; y =[ 0.98,0.92,0.81,0.64,0.38] ; c= newpoly(x, y )c =-0.5208 0.8333 -1.1042 0.1917 0.9800%%三次样条插值 x=[0.2,0.4,0.6,0.8,1.0]; y=[0.98,0.92,0.81,0.64,0.38];x0 = [0.2,0.28,1.0,1.08]; pp=csape(x,y,variational);%