实验2：插值法应用实训.docVIP

河北联合大学实验教案 第2012-2013-1学期 实验2： 插值法应用实训 教学目的与要求：了解插值问题及其适用的场合，理解并掌握常用的插值算法的构造和计算，了解差商概念、Runge现象及样条插值方法，学习用计算机求近似函数的一些科学计算方法和简单的编程技术。 ■ 教学内容： 程序设计基础知识： Interpolation插值函数的一般形式： Interpolation[data，InterpolationOrder-n ] 对数据data进行插值运算，并可设置插值多项式的次数n，默认值为3，Interpolation生成一个InterpolatingFunction[插值范围，]目标（如Out[2]），系统不显示所构造的插值函数，因此，用户直接用生成的插值函数计算函数的近似值。 编写一个多项式插值的Mathematica程序： In[1] : =data=Table[{{1,16},{2,12},{4,8},{5,9}}]； In[2] : =g=Interpolation[data ,InterpolationOrder-3] Out[2]=InterpolatingFunction[[1,5],] In[3] : =g[1.2]（*计算x=1.2处函数的近似值*） Out[3]=15.1307 插值法程序示例：f[x_]:=Exp[x] A=Table[{x,f[x]},{x,0,0.8,0.2}]//N g1=ListPlot[Table[A],Prolog-AbsolutePointSize[18]]; Interpolation[A,InterpolationOrder-3] g2=Plot[%[x],{x,0,0.8}] Show[g1,g2] N[%%%[0.12],20] N[%%%%[0.72],20] N[f[0.12],20] N[f[0.72],20] 2) 求Lagrange插值多项式程序 Clear[lag,xi,x,yi]; xi=Input[xi=] yi=Input[yi=] n=Length[xi]-1; p=Sum[yi[[i]]* (Product[(x-xi[[j]])/(xi[[i]]-xi[[j]]),{j,1,i-1}]* Product[(x-xi[[j]])/(xi[[i]]-xi[[j]]), {j,i+1,n+1}]),{i,1,n+1}]; lag[x_]=Simplify[p] 说明 本程序用于求n次Lagrange插值多项式。程序执行后，按要求通过键盘输入插值基点xi:{x0 , x1, ... , xn }和对应函数值yi:{ y0 , y1 , … , yn }后，程序即可给出对应的n次Lagrange插值多项式lag[x]。 程序中变量说明 xi：存放插值基点{x0 , x1, ... , xn } yi: 存放对应函数值{y0 , y1 , … , yn} lag[x]: 存放求出的n次Lagrange插值多项式Ln(x) 注：语句lag[x_]=Simplify[p]用简化形式给出对应的n次Lagrange插值多项式。）求Newton插值多项式程序 Clear[newt,s,x]; xi=Input[xi=] yi=Input[yi=] n=Length[xi]; (*下面计算差商表*) f=Table[0,{n},{n}]; Do[f[[i,1]]=yi[[i]],{i,1,n}] Do[f[[i,j+1]]=(f[[i,j]]-f[[i+1,j]])/(xi[[i]]-xi[[i+j]]), {j,1,n-1},{i,1,n-j}] Print[差商表] Do[Print[xi[[i]], ,f[[i]]],{i,1,n}] (*下面求Newton插值多项式*) fa=1; s=f[[1,1]]; Do[fa=(x-xi[[k]])*fa; s=s+fa*f[[1,k+1]], {k,1,n-1}] newt[x_]=s Simplify[%] 说明 本程序用于求n次Newton插值多项式。程序执行后，按要求通过键盘输入插值基点xi:{x0 , x1, ... , xn }和对应函数值yi:{ y0 , y1 , … , yn }后，程序依次给出输入的数据表、计算出的差商表、Newton插值多项式、Newton插值多项式的简化形式。 )求三次样条插值函数程序 mar=Input[边界条件类型=1或2(1:第一型,2：第二型)]; f1=Input[输入左边界条件函数值=] f2=Input[输入右边界条件函数值=] xi=Input[输入插值节点xi