《数值计算课后习题四答案.pdf》由会员分享，可在线阅读，更多相关《数值计算课后习题四答案.pdf（16页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、习 题 四 解 答1、设，写出的一次插值多项式，并估计插值误差。解：根据已知条件，有x 0 1y 1设插值函数为，由插值条件，建立线性方程组为解之得则因为所以，插值余项为所以02、给定函数表-0.1 0.3 0.7 1.10.995 0.995 0.765 0.454选用合适的三次插值多项式来近似计算f(0.2)和f(0.8)o解：设三次插值多项式为，由插值条件，建立方程组为即解之得则所求的三次多项式为。所以3、设 是n+1个互异节点，证明:(1);(2)。证 明:(1)由拉格朗日插值定理，以x0,x1,x2,xn为插值节点，对y=f(x)=xk作n次插值，插值多项式为而 yi=xik,所以同
2、时，插值余项所以结论得证。(2)取函数对此函数取节点，则对应的插值多项式为由余项公式，得所以令 t=x,4、给定数据()x 2.0 2.1 2.2 2.4f(x)1.414214 1.449138 1.48320 1.54919(1)试用线性插值计算f(2.3)的近似值，并估计误差；(2)试用二次Newton插值多项式计算f(2.15)的近似值，并估计误差。解：用线性插值计算f(2.3),取插值节点为2.2和2.4,则相应的线性插值多项式是用x=2.3代入，得(2)作差商表如下Xf(x)一阶差商2.0 1.4142140.35012.1 1.449138-0.0470.3407 4.10752
3、.2 1.48320 1.5960.65992.4 1.54919二阶差商 三阶差商根据定理2,f(x)=f(x0)+fx0,x1(x-x0)+fx0,x1,x2(x-x0)(x-x1)+fx0,x1,，xn(x-x0)(x-x1)-(x-xn1)+fx0.x1,，xn,x 兀(x)o以表中的上方一斜行中的数为系数，得f(2.15)=1.41421+0.3501 X(2.15-2.0)-0.047 X(2.15-2.0)X(2.15-2.1)=1.663725指出：误差未讨论。5、给定函数表x 0 1 2 4 5y 0 16 46 88 0试求各阶差商，并写出牛顿插值多项式和插值余项。解：作差
4、商表如下x f(x)一阶差商 二阶差商 三阶差商 四阶差商0 0161 16 7302 46-3214 88-8 85 0根据定理2,以表中的上方一斜行中的数为系数，得指出：余项未讨论。5*、给定函数表x 0 1 2 3 4y 0 16 46 88 0试求各阶差分，并求等距节点插值。解：由已知条件，显然，xO=0,h=1,x=t(作差分如下Xf(x)一阶差分二阶差分 三阶差分四阶差分00116161430-224612-1 4 08842-1 4 24-1 3 0-8 850根据等距节点插值公式,指出：在本题这种情况卜，实 际 上，也就是说，在这样的条件下,t 的多项式就是x 的多项式，可以直
5、接转换。一般情况下，把t的关系转换为X的关系需要根据x=xO+th,将t用x表示，即将代入得到的多项式。6、给定数据表x 0.125 0.250 0.375 0.500 0.625 0.750f(x)0.79618 0.77334 0.74371 0.70413 0.656320.60228试用三次牛顿差分插值公式计算f(0.1581)及f(0.636)。解：所给节点是等距结点：O计算差分得X f(x)一阶差分二阶差分 三阶差分 四阶差分五阶差分0.1250.79618-0.022840.2500.77334-0.00679-0.02963-0.003160.3750.74371-0.0099
6、50.00488-0.039580.00172-0.004600.5000.70413-0.008230.00028-0.047810.002000.6250.65632-0.00623-0.054040.7500.60228令，根据等距结点插值公式，得则O7、设f(x)在-4,4有连续的4阶导数，且(1)试构造一个次数最低的插值多项式p(x),使其满足(2)给出并证明余项f(x)-p(x)的表达式。解:(1)由7*可以求出满足的三次埃尔米特插值多项式O设，则 P(x)满足由得所以O(2)余项具有如卜结构作辅助函数则显然在点处有6 个零点(其中0,3 是二重零点)，即不妨假设。由罗尔定理，存在
7、，使得，再注意到，即 有 5 个互异的零点再次由罗尔定理得，存在，使得第三次应用罗尔定理得，存在使得，第四次应用罗尔定理得，存在使得，第五次应用罗尔定理得，存在使得注意到（中p是4 次函数，其 5 次导数为0）。所以代入余项表达式，有指出：本题是非标准插值问题，比较简单的求解方法有：求插值问题的基本方法是待定系数法。以本题来说，有 5 个条件，可以确定一个4次的插值多项式，设为，将条件代入，建立一个5 元的线性方程组，求出各参数，就可以求出插值多项式。求插值问题的第二种方法是基函数法，即根据给定条件设定插值多项式的结构和各基函数的结构，根据条件确定基函数即可。具体方法与拉格朗日插值基函数构造和
8、埃尔米特插值基函数构造相似。以标准插值为基础的方法是一种更简单的方法，本题中，首先利用4 个条件构造一个埃尔米特插值，在此基础上设定所求插值多项式的一般形式，保证其满足埃尔米特插值条件，代入未利用条件解方程（组），求出其中的未知参数，即可求出插值多项式。本题也可以先利用构造一个2 次插值多项式，以此为基础构造4 次插值多项式，的结构是满足再根据列出两个线性方程组成的方程组，求出a、b 两个参数，即可求出所求的插值多项式。求插值函数余项的常用方法是：应具有如下形式（以本题为例）作辅助函数则在点处有6 个零点（其中0,3 是二重零点）。反复应用罗尔定理，直到至少有一个，使 得。此时即有代入余项表达
9、式即可求出。7*、设f(X)在卜4,4 有连续的4 阶导数，且试用两种方法构造三次埃尔米特插值多项式H(x),使其满足0解一(待定系数法)：解：设，则由插值条件得解之得，所以。解二(基函数法)：解：设，因为线性拉格朗日插值基函数为，由得同理由得则O8、设，试作一个二次多项式p(x),使其满足,并导出余项估计式。解：设此二次式为，因为，所以，由已知条件将 其 代 入，得所以，要求的二次多项式为O因 为0是2重零点，1是1重零点，因此可以设余项具有如下形式:其 中K(x)为待定函数。固 定x,作辅助函数显然不 妨 假 设。由罗尔定理，存 在，使 得，再注意到再次由罗尔定理得，存 在，使得再次应用罗
10、尔定理，存在使得O注意到(中p 是2次函数，其3次导数为0)。所以代入余项表达式，有指出:石瑞民 数值计算关于余项讨论很清楚。9、给出sinx在 0,叫上的等距结点函数表，用线性插值计算sinx的近似值，使其截断误差为，问该函数表的步长h 取多少才能满足要求？解：设为等距结点，步长为h,则当时，作f(x)的线性插值则有由此易知因此由，得。指出：关于最大值的计算与12题相同。10、求 在区间 a,b 上的分段埃尔米特插值，并估计误差。解：由分段三次埃尔米特插值多项式则的分段埃尔米特插值为其中其余项估计式为O11、已知数据表i 0 1 22.5 7.5 104.0 7.0 5.00.13-0.13
11、求三次样条插值函数。解：这是第一类边界条件，要求解方程组其中将以上数据代入方程组解之得将获得的数据代入到中，得12、设(具有二阶连续导数)，且f(a)=f(b)=0,证明:证明：以a、b 为节点进行插值，得因为在处取得最大值，故13.给定数据表x-2 -1 0 1 2y-0.1 0.1 0.4 0.9 1 .6用两种方法求其二次拟合曲线。解一：设所求的拟合函数为，则。对 a、b、c 分别求偏导，并令偏导数等于0,得将各数据点的数值代入，得方程组为解之得 a=0.4086,b=0o 42,c=0.0857,所以数据点所反映的函数的近似关系为解二：设所求的拟合函数为，将数据代入方程得方程组的系数矩
12、阵和右端向量为因为所以解之得 a=0.4086,b=0o 42,c=0.0857,所以数据点所反映的函数的近似关系为14、已知试验数据x 19 25 31 38 44y 19.0 32.3 49.0 73.3 97.8用最小二乘法求形如的经验公式，并计算均方误差。解：设则对 a、b 分别求偏导，并令偏导数等于0,得将数据代入得化简得第二个方程减去第一个方程乘以1065进步化简得解之得则x 与y 的函数关系是y=1.01+0.05x2,此时，平方逼近误差为所以，均方误差为。指出：均方误差实际上就是按最小二乘法则确定的残差。15、观测物体的直线运动，得出如下数据：时间 t(s)0 0.9 1.9
13、3.0 3.9 5.0距离 s(m)0 10 30 50 80 110求运动方程。解：设运动方程为s=a+bt则将上述数据代入方程组得方程组解之得所以，。指出：利用统计型计算器，有关中间数据可以简单求出。16、在某化学反应中，由实验得分解物浓度与时间关系如下:时 间 t 0 5 10浓度15 20 25 30 351.27 2.162.86 3.44 3.87 4.15 4.37时 间 t 4045 50 550浓度4.514.58 4.62 4.64用最小二乘法求y=f(t)o解：描草图，观察草图可以发现，该组数据分布近似于指数函数曲线，而且随着t 的增大，y 的增速放缓，故设两边取对数，得
14、令，则拟合函数转化为线性拟合关系。将上述数据代入得解之得所以指出：（1）T=0,该拟合函数不适用。（2）专业的变化规律（经验函数）应当由专业人员给出。仅仅从有限数据的草图得出的规律可能不具普遍性。17、给定数据表x 7.2 2.7 3.5 4.1 4.8y 65 60 53 50 46用最小二乘法求形如的经验公式。解：对两边取对数，得令，则代入数据，建立方程组为解之得所以O18、用最小二乘法求方程组的近似解。分析：这是方程个数多于未知数个数的超定方程组，是矛盾方程组，用最小二乘法求解。解：设方程组中各个方程的一般形式为，则对x、y分别求偏导，并令偏导数等于0,得将数据代入得解之得19、已知数据表X12345678y 15.3 20.5 27.4 36.6 49.1 65.6 87.8 117.6它有形如的拟合函数，试求本问题的最小二乘解。解：令，则拟合函数变形为,原拟合问题转化为线性拟合问题。则。对 a、b 分别求偏导，并令偏导数等于0,得将数据代入，得解之得所以，所求的拟合函数为O20、在平面上给出三个点，它们的坐标是，每个点对应一个函数值，找出一个通过这三个点的平面。解：这实际上是求过三个点的平面方程。由解析儿何知识可知，平面的三点式方程为将三点坐标代入，解此方程就可求出所求平面方程。（以下从略）