数值分析上机实验——插值

《数值分析上机实验——插值》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、庶爾牘N理2厚卖验报告课程名称数值分析实验项目名称插值实验实验类型上机实验学吋2班级20111131学号2011113130姓名张振指导教师沈艳实验室名称理学楼407实验时间2013.10.11实验成绩预习部分实验过程表现实验报告部分总成绩教师签字日期哈尔滨工程大学教务处制实验一插值实验一.插值法的基本思想1.Lagrange插值多项式的基本思想及公式：Lagrange插值多项式的基本思想是把多项式的构造问题转化成n+1个基函数.（X）的构造。其公式如下：L”（x）二外仏⑴.其中ljx）=（7）…（f）（fm）用0,1

2、,…,n.（耳_兀0）…（忑（"一忑+J…（耳-斗）2.Newton插值多项式基本思想及公式：Newton插值多项式是Lagrange插值多项式的改进，是由基函数｛1,（x・x°）,…，（x・xo）.・.（x・x“_］）・｝逐次地推得到的多项式。其公式如下：P“(x)=ao+a】(x・Xo)+art(x-x0)...(x-xw_1).其中=f[Xo，X[,…,Xr],k=0,l,…n.3.Hermite插值的基木思想及公式：Hermite插值不仅满足各插值点上函数值相等，而且各节点上导数值也想等,甚至高阶导数值也相等。

3、在这里只给出两种典型的Hermite插值：①P(x)=f(x())+f[x0,x)](x-x())+f[x(),xl,x2](x-x0)(x-xI)+A(x・x0)(x-x!)(x-x2)其屮系数A二f'(xJ・f[Xo，X

4、]・f[Xo，X

5、,X2](X

6、・Xo)八八(Xj-x0)(x,-x2)余项R(x)=lf(4)(©)(x・x())(x・xJ2(x・X2).4!〜②其中(1+2「x-x{I*xoJ⑴-州丿/、/、21+2—兀一X。I®~x)l坷-“0丿h°(x)h(兀)=/、2【x-xxHo（兀）=（兀一兀

7、o）卫0一州丿H(%)=hQ(x)yQ+/?j(x)y{+(兀)几+H〕(x)y、2W1(x)=(x-x1)余项R3(x)=f(x)-H3(x)=(x-x0)2(x-x()2兵(a,b)4•三次样条函数S(x)wL［a,b］，且在每个小区间抹沖鬥］上是三次多项式，其中a=x0

8、•实验题目及实验目的1.(第2章习题20)给定数据表如下：X/0.250.300.390.450.53y；0.50000.54770.62450.67080.7280试求三次样条插值S(x),并满足条件：(1)S(0.25)=1.0000,S(0.53)=0.6868；(2)S"(0.25)=S(0.53)=0.2.(第2章计算实习题3)下列数据点的插值X01491625364964Y012345678可以得到平方根函数的近似，在区间［0,64］上作图。(1)用这9个点作8次多项式插值L*(x).(2)用三次样条(第一

9、边界条件)程序求S(x).从得到结果看在［0,64］上，哪个插值更精确；在区间［0,1］上，两种插值哪个更精确？三.实验手段：指操作环境和平台:'¥ir)7系统b*MATLABR2009a程序语言：一种类似C语言的程序语言，但比C语言要宽松得多，非常方便。四•程序1.①三次样条插值程序：functiony=scyt(x0,yO,x)n=length(xO)；m=length(x)；y=zeros(1#m)；p(n)=0;q(l)=1;d(l)=ZL.OOOO;d(n)=0.6868;fork=2：n-1,h(1)=xO

10、(k)-xO(k-1)；h(2)=xO(k+1)-xO(k)；p(k)=h(l)/(h(l)+h(2));q(k)=h(2)/(h(l)+h(2));d(k)=6*((yO(k+l)-yO(k))/h(2)-(yO(k)-yO(k-1))/h(l))/(h(l)+h(2))；endb(l)=q(l)/2;fork=2zn-l,b(k)=q(k)/(2-p(k)*b(k-1))；endtemp(1)=d(1)/2；fork=2zn-l,temp(k)=(d(k.)・p(k)*temp(k-1))/(2-p(k)*b(k-

11、1))；endM(l)=0;M(n)=0;fork=1：-1：2,M(k)=temp(k)-b(k)*M(k+l)；endfori=1：m#fork=1：n-l,if(x(i)>=xO(k))&&(x(i)<=xO(k+1)),h=xO(k+1)-xO(k)；y(i)=M(k)*(xO(k+1)-x(i))^3/(6*h)x(i)