实验报告插值法.doc

实验目的和要求实验内容和步骤endv=v+w*y (k);end 实验目的和要求 实验内容和步骤 end v=v+w*y (k); end 计算机上机实验报告 专业和班级 姓名 成绩 学号 课程名称 数值计算方法 实验名称 插值法 买验H的 1、 掌握用MATLAB计算拉格朗日、分段线性、三次样条三种插 值的方法，改变节点的数目，对三种插值结果进行初步分析。 2、 掌握用MATLAB作线性最小二乘的方法。 3、 通过实例学习如何用插值方法与拟合方法解决实际闷题，注意 二者的联系和区别。 实验的主要内容 1、编制拉格朗日、牛顿插值程序，并运行一个简单的实例。 (1)拉格朗日插值程序： function v=polyinterp(x,y, u) n=length(x); v=zeros(size(u)); for k=l:n w=ones (size (u)); for j=[l:k-l k+1:n] w=(u-x ( j)) ./(x(k)-x (j)) .*w; (2)牛顿插值程序： function y=newinterp(X,Y,x)% 牛顿插值函数 m=length(X); for k=2:m for j=l:k-l Y(k)= (Y(k)- Y(j))/(X(k)-X(j)); end end y=Y(m); for :1 y=y.*(x-X(j))+Y(j)； end 实例:当x=144, 169, 225时，y=12, 13, 15,用牛顿差值 法求根号175。如下： Current Directory - C:\MATLAB7\workAll Fifes L.newinlerp mFile Type Current Directory - C:\MATLAB7\work All Fifes L. newinlerp m File Type M-file Lad McxSticd 2013-4-10 15:41:02 ccmmana winaow ? K=[144 169 225】: ? Y=[12 13 15]; ? x=l?5; ? y=newinterp (3J, Y^x) y = 13. 2302 2JCurrent Orectory | Workspace | 2J Current Orectory | Workspace | Command History X kf=【12 13 15】： ? [140:1:230]: ? y=newiivterp(X, Y,x); ? plot(x,y) ? 2、给定函数，已知： /(2.0) = .1.414214 /(2. l) = .l .449138 /(2.2) = .l .483240 /(2.3) = .1.516575 /(2.4) = .l .549193 用牛顿插值法求4次插值多项式在2.15处的值，以此作为函 Command Window ? X=[2. 0 2. 1 2. 2 2. 3 2.4]; ? Y=[l.414214 1.449138 1.483240 1.516575 1. 549193]; * ? x=2.15; y=newint erp (H x) y = 1.4663 IJ ? x=[1.9:0.006:2.5]; y=newint erp (H x); plot (Xj y) ? 选择函数y=exp（-x2） （-2x2），在n个节点上（n不要太大，如 5?11）用拉格朗円、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m要适中，如50?100）。通过数值和图形输出, 将三种插值结果与精确值进行比较。适当增加n，在作比较，由此 作初步分析。 程序： %不同插值方法是否会出现震荡runge现象 %M文件 function runge 10 ［X，Y］=fenduan（10，l）;%将［-1，］区间分成10等份，返回对应的（x，y） 五组数据 x=linspace（-2,2，l00）;%将［-1，1 ］划分成100等份，以便作出样条插值 多项式的图形。 for i=l :length（x）%绘制原函数曲线图 y（i）=exp（-x（i）A2）; end hold on plot（x，y）; text（0,1 /\leftarrow原函数’）％对曲线添加标注 y=newinterp（X，Y，x）;%多项式插值中的牛顿插值法 hold on plot（x，y）; title（’插值函数中的runge现象，区间等分为10段’）;％添加标题 xlabel（’X 轴?）； ylabel（’Y 轴?）; text（-0.9，l .5，’\leftarrow牛顿插值’）％对曲线添加标注 y=interpl（X，Y，x）; pl