Mathematica在数值分析中的应用 数学软件与数学实验 教学课件.pptVIP

求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插值多项式。 例．给定数据表 x 0 1 2 3 y=f(x) 1 3 5 12 用Lagrange插值法求三次插值多项式,并给出函数f(x)在x =1.4的近似值。 解： 执行Lagrange插值程序后，在输入的两个窗口中按提示分别输入{0, 1, 2, 3}、{1, 3, 5, 12}，每次输入后用鼠标点击窗口的“OK”按扭，得如下插值函数。 6 + 22 x - 15 x2 + 5 x3 ----------------------- 6 所以得到三次插值多项式L3(x)=1+11 x/3-5 x2/2+5 x3/6 接着键入“lag[1.4]”，则输出3.52，因此f(x)在x =1.4的近似值为3.52，即f(1.4)?3.52. 【例】 多项式插值的误差估计式中可以看到，当插值节点越多时误差会越小，这个结论正确吗？通过实验说明该结论的正确性。 解: 考虑函数f(x) = (1+x2)-1 在区间[-4,4]内选取不同个数的等距插值节点做观察，这里分别选[-4,4] 内的9个和11个的等距节点来做实验，将对应的插值函数图与被插函数f(x) = (1+x2)-1画在一起做观察，为简单起见，这里用Mathematica 命令做实验，对应命令为 u=Table[{x,(1+x^2)^-1},{x,-4,4}] ; (*采取f(x) 在[-4,4] 内的9个插值点 *) g=ListPlot[u, PlotStyle-PointSize[0.04]] (*将散点图图形文件存放在变量g中*) s=InterpolatingPolynomial[u , x] ; (*将插值函数存放在变量s中*) t= Plot[{s, (1+x^2)^-1}, {x,-4,4}, PlotStyle-{{Thickness[0.005]},{Thickness[0.008]}}] (*将插值函数s与f(x)画在一起的图形文件存放在变量t中*) Show[t, g] (*将散点图, 插值函数s, f(x)画在一个坐标系中*) 在[-4,4]中选9个等距节点的插值函数与被插函数图，粗线为被插函数图 u1=Table[{x,(1+x^2)^-1},{x,-4,4,0.8}] ; (*采取f(x) 在[-4,4] 内的11个插值点 *) g1=ListPlot[u, PlotStyle-PointSize[0.04]] s1=InterpolatingPolynomial[u1 , x] ; t1= Plot[{s1, (1+x^2)^-1}, {x,-4,4}, PlotStyle-{{Thickness[0.005]},{Thickness[0.008]}}] Show[t1, g1] 在[-4,4]中选11个等距节点的插值函数与被插函数图 * * 插 值 法 Mathematica 在数值分析 中的应用 Mathematica For Windows 实际问题 期望 试验数据 观测数据 期望 内在规律 期望 函数关系 一、数学的期望 实验数据是否存在内在规律? 实验数据的内在规律是