数值分析次作业及参考答案.docVIP

数值计算方法第一次作业及参考答案 已测得函数的三对数据：（0，1），（－1，5），（2，－1）， （1）用Lagrange插值求二次插值多项式。（2）构造差商表。（3）用Newton插值求二次插值多项式。 解：(1)Lagrange插值基函数为 同理 故 （2）令，则一阶差商、二阶差商为 实际演算中可列一张差商表： 一阶差商 二阶差商 0 1 －1 5 －4 2 －1 －2 1 （3）用对角线上的数据写出插值多项式 在上给出的等距节点函数表，若用二次插值求的近似值，要使截断误差不超过，问使用函数表的步长应取多少？ 解： 求在［a,b］上的分段线性插值函数，并估计误差。 解： 已知单调连续函数的如下数据 －0.11 0.00 1.50 1.80 －1.23 －0.10 1.17 1.58 用插值法计算约为多少时（小数点后至少保留4位） 解：作辅助函数则问题转化为为多少时，此时可作新的关于的函数表。由单调连续知也单调连续，因此可对的数值进行反插。的牛顿型插值多项式为 故 设函数在区间[0,3]上具有四阶连续导数,试用埃尔米特插值法,求一个次数不高于3的多项式，使其满足，，, 。并写出误差估计式。 解：由所给条件可用埃尔米特插值法确定多项式， 由题意可设 为确定待定函数，作辅助函数： 则在[0,3]上存在四阶导数且在[0,3]上至少有5个零点为二重零点），反复应用罗尔定理，知至少有一个零点使，从而得。 故误差估计式为 设函数在节点的函数值均为零，试分别求满足下列边界条件下的三次样条插值函数： （1） （2） 解：（1）取处的一阶导数作为参数，。由于 以及由三转角方程 得 由于从而 解之可得 故 （2）取处的二阶导数作为参数，。由于 以及由三弯矩方程 由于代入方程可得 故 7．编程实现题： 略。 8、试求 最佳一次一致逼近多项式。 解：因为在内不变号，故最佳一次一致逼近多项式为 式中 从而 9、给定，试利用最小零偏差定理，即切比雪夫多项式的最小零偏差性质，在上求的三次最佳一致逼近多项式。 解：令 设为在上的三次最佳一致逼近多项式，由于的首项系数为，故 10、设，分别在上求一函数，使其为的最佳平方逼近，并比较其结果。 解： 由结果知（1）比（2）好。 11、用最小二乘法求一个形如的经验公式，使它与下列数据拟合，并计算均方误差。 19 25 31 38 44 19.0 32.3 49.0 73.3 97.8 解： 12、用格拉姆－施密特方法构造正交多项式求在［0，1］上的二次最佳平方逼近多项式。（参考讲义与参考书） 解： 构造正交多项式 于是 所以，在［0，1］上的二次最佳平方逼近多项式为 13、求在［－1，1］上的三次最佳平方逼近多项式。（参考讲义与参考书，利用Legendre正交多项式） 解 先计算。 ； ； ； ； 又有 ， ， ， 得 均方误差 14、 A、B、C三点连成一条直线，AB长为，BC长为，某人测量的结果为米，米，为控制丈量的准确性，又测量米，试合理地决定和的长度。（小数点后取四位有效数字） 解：令为AB的所求值，为BC的所求值，则 在最小二乘意义下，要达到极小， 即求的极小点。 令 解的。故应取。 15、求函数在区间［－1，1］上的近似3次最佳一致逼近多项式有哪几种方法？选一种方法解本题，并估计误差。（参考讲义与参考书） 解：三种方法，见参考讲义。 截断切比雪夫级数 由富利叶级数系数公式得 ， 它可用数值积分方法计算，得到 由 及的公式得到 拉格朗日插值余项的极小化 由的4个零点 做插值点可求得 ， 台劳级数项数的节约 应用的台劳展开，取，得 作为的近似，其误差为 ， 由于 则 其中 用做的逼近多项式，其误差为 若再用代入可求出 16.编出用正交多项式(格拉姆－施密特)作最小二乘拟合的程序或框图。（参考讲义与参考书） 略。