数值分析第三章答案.pdf

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

数值分析第三章答案

【篇一：常州大学数值分析作业第三章】

答：matlab程序

function[a,y]=lagrange(x,y,x0)%检验输入参数

ifnargin2||nargin3

error(incorrectnumberofinputs);end

iflength(x)~=length(y)

error(thelengthofxmustbeequaltoitofy);end

m=length(x);n=m-1;l=zeros(m,m);%计算基本插值多项式的系数

fori=1:n+1c=1;

forj=1:n+1ifi~=j

ifabs(x(i)-x(j))epsabs(x(i)-x(j))eps

error(therearetwotwosamenodes);end

c=conv(c,poly(x(j)))/(x(i)-x(j));

endend

l(i,:)=c;end

%计算lagrange插值多项式的系数a=y*l;

%计算f(x0)的近似值ifnargin==3

y=polyval(a,x0);

工程（专）学号

end

a=fliplr(a);return

[a,y]=lagrange(x,y,x0);p1=vpa(poly2

sym(a),3)y

[a,y]=lagrange(x,y,x0);p2=vpa(poly2sym(a),3)y

p2=x2-0.109x-0.336y=0.5174

[a,y]=lagrange(x,y,x0);p4=vpa(poly2sym(a),3)y

p4=x4+0.00282x3-0.514x2+0.0232x+0.0287y=0.5001

次多项式在2.8处的值。

答：matlab程序function

[t,y0]=aitken(x,y,x0,t0)ifnargin==3t0=[];end

n0=size(t0,1);

m=max(size(x));n=n0+m;

t=zeros(n,n+1);

t(1:n0,1:n0+1)=t0;t(n0+1:n,1)=x;t(n0+1:n,2)=y;ifn0==0i0=2;

else

i0=n0+1;end

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

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

fori=i0:n

forj=3:i+1

t(i,j)=fun(t(j-2,1),t(i,1),t(j-2,j-1),t(i,j-1),x0);endend

y0=t(n,n+1);return

function[y]=fun(x1,x2,y1,y2,x)y=y1+(y2-y1)*(x-x1)/(x2-x1);

return

%选取0、1、3、4四个节点，求三次插值多项式x=[0,1,3,4];

y=[0.5,1.25,3.5,2.75];x0=2.8;

[t,y0]=aitken(x,y,x0)t=

00.50000001.01.25002.6000003.03.50003.29993.23000

4.02.75002.07502.28503.4190y0=3.419000000000000

16、选取适当的函数y=f(x)和插值节点，编写matlab程序,分别利

用lagrange插值方法,newton插值方法确定的插值多项式，并将函

数y=