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

插 值 法;实际问题;实验数据是否存在内在规律?; 插值与拟合的基本问题;而不知道函数在其他点x上的取值，这时只能用一 个经验函数y=g(x)对真实函数y=f(x)作近似。;;求任一插值点;?;六、插值法的一般定义;; Lagrange 插 值;定理1;令：; 如何确定 ;已知 n+1个节点;;;;;;拉格朗日(Lagrange) 插值多项式;四、拉格朗日插值的余项与误差估计;;例3;;多项式插值法：;求Lagrange插值多项式程序;程序中变量说明 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个等距节点的插值函数与被插函数图