数值分析第五章答案.pdfVIP

数值分析第五章答案--第1页

数值分析第五章答案

【篇一：数值分析第五版计算实习题】

第二章

2-1

程序：

clear;clc;

x1=[0.20.40.60.81.0];

y1=[0.980.920.810.640.38];

n=length(y1);

c=y1(:);

orj=2:n%求差商

fori=n:-1:j

c(i)=(c(i)-c(i-1))/(x1(i)-x1(i-j+1));

end

end

symsxdfd;

df(1)=1;d(1)=y1(1);

fori=2:n%求牛顿差值多项式

df(i)=df(i-1)*(x-x1(i-1));

d(i)=c(i)*df(i);

end

disp(4次牛顿插值多项式);

p4=vpa(collect((sum(d))),5)%p4即为4次牛顿插值多项式,并保留

小数点后5位数pp=csape(x1,y1,variational);%调用三次样条函数

q=pp.coefs;

disp(三次样条函数);

fori=1:4

s=q(i,:)*[(x-x1(i))^3;(x-x1(i))^2;(x-x1(i));1];

s=vpa(collect(s),5)

end

x2=0.2:0.08:1.08;

dot=[121112];

figure

ezplot(p4,[0.2,1.08]);

holdon

y2=fnval(pp,x2);

x=x2(dot);

数值分析第五章答案--第1页

数值分析第五章答案--第2页

y3=eval(p4);

y4=fnval(pp,x2(dot));

plot(x2,y2,r,x2(dot),y3,b*,x2(dot),y4,co);

title(4次牛顿插值及三次样条);

结果如下：

4次牛顿插值多项式

p4=-0.52083*x^4+0.83333*x^3-1.1042*x^2+0.19167*x+

0.98三次样条函数

x∈[0.2,0.4]时，s=-1.3393*x^3+0.80357*x^2-0.40714*x+

1.04x∈[0.4,0.6]时，s=0.44643*x^3-1.3393*x^2+0.45*x+

0.92571x∈[0.6,0.8]时，s=-1.6964*x^3+2.5179*x^2-1.8643*x

+1.3886x∈[0.8,1.0]时，s=2.5893*x^3-7.7679*x^2+6.3643*x-

0.80571输出图如下

2-3（1）

程序：

clear;

clc;

x1=[01491625364964];

y1=[012345678];%插值点

n=length(y1);

a=ones(n,2);

a(:,2)=-x1;

c=1;

fori=1:n

c=conv(c,a(i,:));

end

q=zeros(n,n);

r=zeros(n,n+1);

fori=1:n

[q(i,:),r(i,:)]=deconv(c,a(i,:));%wn+1/(x-xk)

end

dw=zeros(1,n);

fori=1:n

dw(i)=y1(i)/polyval(q(i,:),x1(i));%系数