数值分析报告插值多项式.docVIP

插值多项式 一、实验目的： 熟悉插值多项式构造，掌握数据插值的思想方法，了解用内部函数和掌握编程方法解决实验问题；了解龙格现象，通过拉格朗日插值、分段线性插值以及三次样条插值进行插值效果的比较和分析；掌握上机编程与调试能力。 二、实验内容： 给定被插函数插值区间为[1，8]，节点为，分别用分段线性和Hermite多项式进行插值。并和7次Lagrange或Newton多项式进行比较。 三、实验要求： 1．分别用内部函数和编程实现。2．可以自拟实验题，要求同上。 四、实验程序及结果： 方法一：内部函数法 函数：Interpolation InterpolatingPolynomial Table[t[i]=i,{i,1,8}]; Table[f[i]=Sin[t[i]]/t[i],{i,1,8}]; Data=Table[{t[i],f[i]},{i,1,8}]; a0=Interpolation[Data,InterpolationOrder?1]; a1=InterpolatingPolynomial[Data,x]//Expand; c0=Plot[a0[x],{x,1,8},PlotStyle?{Black,Thick},Frame?True]; c1=Plot[a1,{x,1,8},PlotStyle?{Red,Thick},Frame?True]; b=Graphics[{PointSize[0.02],Black,Point[Data]}]; d=Table[Plot[Evaluate[Sin[x]/x],{x,t[i],t[i+1]},PlotStyle?{Blue,Thickness[0.001]}],{i,1,7}]; Print[被插函数（蓝）、7次多项式差值（红）与分段线性差值（黑）如下图：] Show[b,c1,c0,d,AspectRatio?2/3,Frame?True] 被插函数（蓝）、7次多项式差值（红）与分段线性差值（黑）如下图： 方法二：编程实现法 h[i_]:=t[i+1]-t[i]; ?[x_,i_]:=(x-t[i])/h[i]; Table[t[i]=i,{i,1,8}]; l[x_,n_,xi_,xj_]:=Product[(x-xi)[[i]],{i,1,n}]/Product[(xj-xi)[[i]],{i,1,n}]; L8[x_]:=l[x,7,{t[2],t[3],t[4],t[5],t[6],t[7],t[8]},t[1]]*f[1]+l[x,7,{t[1],t[3],t[4],t[5],t[6],t[7],t[8]},t[2]]*f[2]+l[x,7,{t[1],t[2],t[4],t[5],t[6],t[7],t[8]},t[3]]*f[3]+l[x,7,{t[1],t[2],t[3],t[5],t[6],t[7],t[8]},t[4]]*f[4]+l[x,7,{t[1],t[2],t[3],t[4],t[6],t[7],t[8]},t[5]]*f[5]+l[x,7,{t[1],t[2],t[3],t[4],t[5],t[7],t[8]},t[6]]*f[6]+l[x,7,{t[1],t[2],t[3],t[4],t[5],t[6],t[8]},t[7]]*f[7]+l[x,7,{t[1],t[2],t[3],t[4],t[5],t[6],t[7]},t[8]]*f[8]; PL[x_,i_]:=f[i](1-?[x,i])+f[i+1]?[x,i]; Table[f[i]=Sin[t[i]]/t[i],{i,1,8}]; f[x_]:=Sin[x]/x;s1=Table[fig[i]=Plot[Evaluate[(PL[x,i]-Sin[x]/x)],{x,t[i],t[i+1]},PlotStyle?{Thickness[0.005],Blue}],{i,1,7}]; s2=Table[fig[i]=Plot[Evaluate[L8[x]-Sin[x]/x],{x,t[i],t[i+1]},PlotStyle?{Thickness[0.005],Red}],{i,1,7}]; k=Table[fig[i]=Plot[Evaluate[PL[x,i]],{x,t[i],t[i+1]},PlotStyle?{Thickness[0.005],Black}],{i,1,7}]; g=Table[fig[i]=Plot[Evaluate[L8[x]],{x,t[i],t[i+1]},PlotStyle?{Thickness[0.005],Red}],{i,1,7}]; s=Table[fig[i]=Plot[Evaluate[Sin[x]/x],{x,t[i],t[i