1-3分段线性插值.pptVIP

内容小结 内 容 小 结 1. 分段线性插值法； 2. 分段线性插值多项式的构造. * 命令格式 插值多项式 命令InterpolatingPolynomial求一个多项式，使给定的数据是准确的函数值，其调用格式如下： InterpolatingPolynomial[{f1，f2，…}，x] 当自变量x为1，2，3，…时的函数值为f1，f2，…。 InterpolatingPolynomial[{{x1，f1}，{x2，f2}，…}，x] 当自变量值为xi时的函数值为fi。 InterpolatingPolynomial[{{x1，{f1，df1，ddf1，…}}，…}，x] 规定xi处的导数函数值为dfi。 如果给定n+1个数据，则结果是n次多项式。 命令后加//Expand，表示按升幂排列。 * 例1：由下列已知条件求插值多项式。 Xi 1 2 3 4 yi 1 2 -1 1 xi 0 1 0.5 yi 1 当x = 1时，y = 2，y′= 2，y″= 1 解：In[1]：=InterpolatingPolynomial[{1，2，-1，1}，x] Out[1]= In[2]：=InterpolatingPolynomial[{{0，1}，{1，Exp[-1]}， {1/2，Exp[-1/2]}}，x] Out[2]= （1） （2） （3） * In[3]：=InterpolatingPolynomial[{{1，{2，2，1}}}，x] Out[3]= 说明：上例中的输出结果表明，所得答案是准确结果，而不是近似函数。其中In[3]是给出当x=1时的函数值和一、二阶导数值，由3个条件得到一个二次多项式。 * InterPolation[{f1，f2，…}] 当自变量值为1，2，3，…时的函数值为f1，f2，…。 InterPolation[{{x1，f1}，{x2，f2}，…}] 当自变量值为xi时的函数值为fi。 InterPolation[{{x1，{f1，df1，ddf1，…}}，…}] 规定xi处的导数函数值为dfi 。 此外还有可选参数： InterPolationOrder→n 指定插值多项式的次数（默认值为3）。 2. 插值函数 函数InterPolation通过在数据点之间进行多项式插值，构造一个近似函数（分段插值函数） ，该函数通过已知的数据点，其调用格式如下： * 执行后返回的结果是形如：InterPolationFunction[domain，table]的插值函数，虽然形式特殊，但能当做普通函数名（如sin）一样理解： InterPolatingFunction[domain，table][x] 就是求x处的函数值。 不过无法直接键入或拷贝上式，使用时还应另起一个函数名来表示它。 * 生成插值函数及其图形，并求当x = 27.1时的函数值。 Xi 19.1 25 30.1 36 40 15.1 50 yi 76.3 77.8 79.25 80.8 82.35 83.9 85.1 例2. 给出数据表如下： 解：In[1]：= data={{19.1,76.3},{25,77.8},{30.1,79.25},{36,80.8},{40,82.35},{45.1,83.9},{50,85.1}}; * f=Interpolation[data] Out[2]= InterpolatingFunction[{{19.1,50.}}， ] In[3]：= pd=ListPlot[data, DisplayFunction-Identity,PlotStyle-PointSize[0.02]]； fd=Plot[f[x],{x,19.1,50},DisplayFunctio-Identity]； Show[{pd,fd},DisplayFunction-$DisplayFunction] * Out[5]= -Graphics- In[6]：=f[27.1] Out[6]=78.3938 说明：上例中生成的插值函数是用f表示的，Out[2]给出插值函数的定义域是[19.1，50]，自变量可以用任何字符表示。由图13-47可知，给定的数据点在函数曲线上。 * f[x_]:=Exp[x] A=Table[{x,f[x]},{x,0,0.8,0.2}]//N g1=