数值计算课后问题详解4文档格式.docx

1、2证明： 1由拉格朗日插值定理，以x0,x1,x2,xn为插值节点，对y=f（x）=xk作n次插值，插值多项式为，而yi=xik,同时，插值余项结论得证。2取函数对此函数取节点，如此对应的插值多项式为由余项公式，得令t=x，4、给定数据f（x）（1）试用线性插值计算f（2.3）的近似值，并估计误差；（2）试用二次Newton插值多项式计算f（2.15）的近似值，并估计误差。用线性插值计算f（2.3），取插值节点为2.2和2.4，如此相应的线性插值多项式是用x=代入，得（2）作差商表如下一阶差商二阶差商三阶差商0.35010.34070.6599根据定理2，f（x）f（x0）fx0，x1（x-x

2、0）+fx0，x1，x2（x-x0）（x-x1）+fx0，x1，xn（x-x0）（x-x1）（x-xn1）+fx0，x1，xn，x（x） 。 以表中的上方一斜行中的数为系数，得f（2.15）1.41421（-2.0） （2.15-2.0） （2.15-2.1） 1.663725指出：误差未讨论。5、给定函数表245164688试求各阶差商，并写出牛顿插值多项式和插值余项。作差商表如下四阶差商73032188根据定理2，以表中的上方一斜行中的数为系数，得余项未讨论。5*、给定函数表3试求各阶差分，并求等距节点插值。由条件，显然，x0=0，h=1，x=t。作差分如下一阶差分二阶差分三阶差分四阶差分

3、1421214042142130根据等距节点插值公式，指出:在此题这种情况下,实际上，也就是说,在这样的条件下,t的多项式就是x的多项式，可以直接转换。一般情况下，把t的关系转换为x的关系需要根据x=x0+th，将t用x表示，即将代入得到的多项式。6、给定数据表0.1250.5000.750试用三次牛顿差分插值公式计算f（0.1581）与f（0.636）。所给节点是等距结点：计算差分得五阶差分0125079618-0.022840250077334-0375074371050007041306250656320750060228令,根据等距结点插值公式，得7、设f（x）在-4,4有连续的4阶导

4、数，且（1）试构造一个次数最低的插值多项式p（x），使其满足 （2）给出并证明余项f（x）-p（x）的表达式。（1）由7*可以求出满足的三次埃尔米特插值多项式设，如此p（x）满足由得（2）余项具有如下结构作辅助函数如此显然在点处有6个零点其中0，3是二重零点，即不妨假设由罗尔定理，存在使得再注意到，即有5个互异的零点再次由罗尔定理得，存在第三次应用罗尔定理得，存在第四次应用罗尔定理得，存在第五次应用罗尔定理得，存在注意到中p（t）是4次函数，其5次导数为0。代入余项表达式，有此题是非标准插值问题，比拟简单的求解方法有：求插值问题的根本方法是待定系数法。以此题来说，有5个条件，可以确定一个4次的

5、插值多项式，设为，将条件代入，建立一个5元的线性方程组，求出各参数，就可以求出插值多项式。求插值问题的第二种方法是基函数法，即根据给定条件设定插值多项式的结构和各基函数的结构，根据条件确定基函数即可。具体方法与拉格朗日插值基函数构造和埃尔米特插值基函数构造相似。以标准插值为根底的方法是一种更简单的方法，此题中，首先利用4个条件构造一个埃尔米特插值，在此根底上设定所求插值多项式的一般形式，保证其满足埃尔米特插值条件，代入未利用条件解方程组，求出其中的未知参数，即可求出插值多项式。此题也可以先利用构造一个2次插值多项式，以此为根底构造4次插值多项式的结构是满足再根据列出两个线性方程组成的方程组，求

6、出a、b两个参数，即可求出所求的插值多项式。求插值函数余项的常用方法是：应具有如下形式以此题为例处有6个零点其中0，3是二重零点。反复应用罗尔定理，直到至少有一个，使得此时即有代入余项表达式即可求出。7*、设f（x）在-4,4有连续的4阶导数，且试用两种方法构造三次埃尔米特插值多项式H（x），使其满足解一待定系数法：，如此由插值条件得解二基函数法：因为线性拉格朗日插值基函数为由得同理由得8、设，试作一个二次多项式p（x），使其满足，并导出余项估计式。设此二次式为所以，由条件将其代入，得所以，要求的二次多项式为因为0是2重零点，1是1重零点，因此可以设余项具有如下形式：其中K（x）为待定函数。固

7、定x，作辅助函数显然再次应用罗尔定理，存在中p（t）是2次函数，其3次导数为0。石瑞民数值计算关于余项讨论很清楚。9、给出sinx在0,上的等距结点函数表，用线性插值计算sinx的近似值，使其截断误差为，问该函数表的步长h取多少才能满足要求？为等距结点，步长为h，如此当时，作f（x）的线性插值如此有由此易知因此关于最大值的计算与12题一样。10、求在区间a,b上的分段埃尔米特插值，并估计误差。由分段三次埃尔米特插值多项式的分段埃尔米特插值为其中其余项估计式为11、数据表i10求三次样条插值函数。这是第一类边界条件，要求解方程组将以上数据代入方程组将获得的数据代入到中，得12、设具有二阶连续导数

8、，且f（a）=f（b）=0，证明：以a、b为节点进展插值，得在处取得最大值，故101用两种方法求其二次拟合曲线。解一：设所求的拟合函数为对a、b、c分别求偏导，并令偏导数等于0，得将各数据点的数值代入，得方程组为解之得a=，b=0。42，c=0.0857，所以数据点所反映的函数的近似关系为解二：将数据代入方程得方程组的系数矩阵和右端向量为14、试验数据1925313844190323490733978用最小二乘法求形如的经验公式，并计算均方误差。对a、b分别求偏导，并令偏导数等于0，得将数据代入得化简得第二个方程减去第一个方程乘以1065进一步化简得如此x与y的函数关系是y=1.x2。此时，平方逼近误差为所以,均方误差为均方误差实际上就是按最小二乘法如此确定的残差。15、观测物体的直线运动，得出如下数据：时间t（s）距离s（m）5080110求运动方程。设运动方程为sa+bt如此将上述数据代入方程组得方程组所以，利用统计型计算器，有关中间数据可以简单求出。16、在某化学反响中，由实验得分解物浓度与时间关系如下：时间t152035浓度