数值分析例题.docxVIP

例2-2 给定函数值表nx 2.3025852.3978952.4849072.564949用二次插值计算 ln(11.25) 的近似值,并估计误差。 在区间[10,12]上lnx 的三阶导数 (2/x3) 的上限 M3=0.002,可得误差估计式注：实际上,ln(11.25)=2.420368, |R2(11.25)|=0.000058 x1.01.41.82.0yi=f(xi)-2.0-0.80.41.2分析：求解如上问题等价于求解x关于y的反函数问题。均差计算表xi ?(xi) 一阶均差二阶均差三阶均差…n阶均差x0x1 x2 x3 xn?(x0) ?(x1) ?(x2) ?(x3) ?(xn)?[x0,x1] ?[x1,x2] ?[x2,x3] ?[xn-1,xn]?[x0,x1,x2] ?[x1,x2,x3] ?[xn-2,xn-1,xn]?[x0,x1,x2,x3] ?[xn-3,xn-2,x2,x3]……………… ?[x0,x1,…,xn]例如 由函数y=(x)的函数表写出均差表.i 0123xi -2-112(xi) 531721解 均差表如下 i xi?(xi) 一阶均差二阶均差三阶均差0123-2-112531721-2743-1-1例2-6 对例如中的 (x),求节点为 x0,x1 的一次插值x0,x1,x2 的 二次插值和 x0,x1,x2,x3 的三次插多项式.解 由差商表知[x0,x1]=-2,[x0,x1,x2]=3, [x0,x1,x2,x3]=-1,于是有N1(x)=5-2(x+2)=1-2xN2(x)=1-2x+3(x+2)(x+1)=3x2+7x+7N3(x)=3x2+7x+7-(x+2)(x+1)(x-1)=-x3+x2+8x+9i xi?(xi) 一阶均差二阶均差三阶均差0123-2-112531721-2743-1-1例2-7 给出 f(x) 的函数表，求4次牛顿插值多项式，并计算f(0.596) 的近似值。xi?(xi) 一阶均差二阶均差三阶均差四阶均差五阶均差0.400.550.650.800.901.050.410750.578150.696750.888111.026521.253821.116001.186001.275731.384101.515330.280000.358930.433480.524930.197330.21300 0.228630.031340.03126-0.00012x012(x)129’(x) 030解法二（用重节点的均差表建立埃尔米特多项式）xi?(xi) 一阶均差二阶均差三阶均差01121229137241例3-3 已知实测数据表如下,求它的拟合曲线xi1 2 3 4 5 yi ωi4 4.5 6 8 8.52 13 1 1 0xy 2 4 68642例3-4 已知实测数据表如下,确定数学模型 y=aebx,用最小二乘法确定a，b。i 0 1 2 3 4 xi yi1.00 1.25 1.50 1.75 2.005.10 5.79 6.53 7.45 8.46 分析：根据给定数据描图也可确定拟合曲线方程，但它不是线性形式。因此首先要将经验曲线线性化。本题可以采取等式两边取对数的形式线性化。数据表中的数值也相应的转化为取对数之后的数值，见下表。i 0 1 2 3 4 xi yi yi1.00 1.25 1.50 1.75 2.005.10 5.79 6.53 7.45 8.461.629 1.756 1.876 2.008 2.135 i 0 1 2 3 4 xi yi 19 25 31 38 4419.0 32.3 49.0 73.3 97.8 xi0 1/8 1/4 3/8 1/2 f (xi)1 0.9973978 0.9896158 0.9767267 0.9588510xi 5/8 3/4 7/81 f (xi)0.9361556 0.9088516 0.8414709 0.8414709例6-3 用SOR迭代