实验报告—数值分析.docx

百度文库- 百度文库-让每个人平等地提升自我 PAGE PAGE # 《数值分析》实验报告 姓 名： 学 号： 专 业： 指导教师： 刘建生 教授 日 期： 2015 年12月25日 百度文库- 百度文库-让每个人平等地提升自我 PAGE PAGE # 实验一 Lagra nge/n ewt on 插值 一：对于给定的一元函数 y f(x)的n+1个节点值yj f(Xj),j 0,1,卅，n。试用Lagra nge 公式求其插值多项式或分段二次 Lagrange插值多项式。 数据如下： 求 五Xj 欠 yj 计算 f (0.596) , f (0.99)的值(提示：结果为 f (0.596) 0.625732 f (0.99) 1.05423 ) Xj 1 2 3 4 5 6 7 yj 试构造Lagrange多项式Lg(x)，计算的f (1.8), f (6.15)值。(提示：结果为 f (1.8) 0.164762 , f (6.15) 0.001266 ) 二:实验程序及注释 MATLAB 程序：function f=lagrange(x0,y0,x ) n=len gth(xO); m=le ngth(yO); format long s=; for k=1: n P=; for j=1: n if j~=k p=p*(x-x0(j))/(x0(k)-x0(j)); end end s=s+y0(k)*p; End f=s; end 百度文库- 百度文库-让每个人平等地提升自我 PAGE PAGE # 结果运行: [O. 4 0- bts 0. eb 0* SO 0. 1. Ob ]: yO= L0- 41075 0. 57filG 0. S067B 0. 90 1. 00 t. 253S2 J : x=0? 596: 】^Er^nepCvt), yip x) ana ■ 0. 6257323?7352S6a ? xQ= L0, 4 Q, 55 0, 65 Q, SO 0- 95 1,Q5 J ; jrO= TO?41075 0? 57815 0? 69675 0. 90 L- 00 I. 25332 1 : )r=C. 99 : lagr ari￡e tarO, yO, vt) an 5 — I、 0642277CS127TE ? xO=ri 2 3 4 5 S ；]： ? y(J= [0.363 C_ 135 0_ OSO 0. OLE 0.007 fl. 002 0.001 ]: ? x=l.S: ? Lagrangeyflj x) ans = 0. 1S47?1894400000 lacr arse CxO^ yfl, 6^ 1 5) ana = 0.001265S2550039L 结果与提示值完全吻合，说明 Lagra nge插值多项式的精度是很 高的； f(x) (x x1)(x x2)(x X3)(X X4)(X x5) (x0 x1)(x0 x2)(x0 x3)(x0 x4)(x0 x5) / (x x0)(x x1)(x x2)(x x3)(x x4) (x5 x0)(x5 x1)(x5 x2)(x5 x3)(x5 x4) 同时，若采用三点插值和两点插值的方法， 用三点插值的精度更 高。若同时采用两点插值，选取的节点距离 x越近，精度越高。 百度文库- 百度文库-让每个人平等地提升自我 PAGE PAGE # 三：采用newton插值进行计算 算法程序如下： format long ; x0=[ ]； yo=[ ]; x=； n=max(size(xO)); y=yo(i); %disp(y); s=1; / dx=yO; for i=1: n-1 dxO=dx; for j=1:n-i dx(j)=(dx0a+1)-dx0(j))/(x0(i+j)-x0(j)); end df=dx(1); s=s*(x-xO(i)); y=y+s*df; %%disp(y); end disp(y) %计算 运行结果： df=dxl ； s=s *(x-xO Ci)); y=yH-￡+df: S%dispy); end dispCy) 0.025732334079595 %计算 百度文库- 百度文库-让每个人平等地提升自我 PAGE PAGE # 绘制出曲线图: o o o Okoan 書 A 启 j.一卜-Newton mterpoistion 0 o o o O koan 書 A 启 j. 一卜-Newton mterpoistion 0 AccLrate value —*— Piecewise interpotatran 与结果相吻合。 所以new