内容发布更新时间 : 2024/5/12 23:55:01星期一 下面是文章的全部内容请认真阅读。

2008信息与计算科学专业计算方法习题参考解答 江世宏编

?8(x??/4)(x??/2)?2?82x(x??/2)??2

L()?6?8(?/6??/4)(?/6??/2)82?/6(?/6??/2)?2??2?2?42?0.8508 9绝对误差为：cos?32?4293?4?82?L()????0.0153 662918?L()66L()6cos相对误差为：????93?4?824?82?0.0179 余项为：

r(x)?sin?x(x??/4)(x??/2)，其中，0????/2 3!1x(x??/4)(x??/2) 6其余项的上界为：r(x)?1??????3r()?(?)(?)?4?0.0239 66664626?比较可知，实际计算所得的绝对误差较余项公式所估计出的值要小一些。

6 已知函数值f(0)?6,f(1)?10,f(3)?46,f(4)?82,f(6)?212，求函数的四阶均差

f[0,1,3,4,6]和二阶均差f[4,1,3]。（均差的计算）

解：采用列表法来计算各阶均差，有

x 0 1 3 4 6 y 6 10 46 82 212 一阶均差 4 18 36 65 二阶均差 14/3 6 29/3 三阶均差 1/3 11/15 四阶均差 1/15 从表中可查得：f[0,1,3,4,6]?x 4 1 3 y 82 10 46 1。 15一阶均差 72/3 18 二阶均差 6 故f[4,1,3]?6。其实，根据均差的对称性，f[4,1,3]?f[1,3,4]?6，该值在第一个表中就可以查到。

5

2008信息与计算科学专业计算方法习题参考解答 江世宏编

7 设f(x)?(x?x0)(x?x1)?(x?xn)求f[x0,x1?xp]之值，其中p?n?1，而节点

xi(i?0,1,?n?1)互异。（均差的计算）

解：由均差可以表示成为函数值的线性组合，有

f[x0,x1?xp]??i?0pf(xi)

(xi?x0)(xi?x1)?(xi?xi?1)(xi?xi?1)?(xi?xp?1)(xi?xp)而f(xi)?0 0?i?p，故f[x0,x1?xp]?0。 8 如下函数值表

x f(x) 0 1 1 9 2 23 4 3 建立不超过三次的牛顿插值多项式。（牛顿插值多项式的构造） 解：

先构造均差表 x 0 1 2 4 f(x) 1 9 23 3 一阶均差 8 14 -10 二阶均差 3 -8 三阶均差 -11/4 故 N(x)?1?8x?3x(x?1)?11x(x?1)(x?2)。 49求一个次数小于等于三次多项式p(x)，满足如下插值条件：p(1)?2，p(2)?4，

p?(2)?3，p(3)?12。（插值多项式的构造）

解法一（待定系数法）：设p(x)?ax3?bx2?cx?d，则

p?(x)?3ax2?2bx?c，由插值条件，有

?a?b?c?d?2?8a?4b?2c?d?4? ?12a?4b?c?3???27a?9b?3c?d?12解得：a?2,b??9,c?15,d??6。 故 p(x)?2x3?9x2?15x?6

6

2008信息与计算科学专业计算方法习题参考解答 江世宏编

解法二（带重节点的均差法）：据插值条件，造差商表

x 1 2 2 3 y 2 4 4 12 一阶差商 2 3 8 二阶差商 1 5 三阶差商 2 故 p(x)?2?2(x?1)?(x?1)(x?2)?2(x?1)(x?2)2?2x3?9x2?15x?6 10 构造一个三次多项式H(x)，使它满足条件H(0)?1,H(1)?0,H(2)?1,H?(1)?1（埃尔米特插值）。

解：设H(x)?ax3?bx2?cx?d，H?(x)?3ax2?2bx?c 利用插值条件，有

?d?1?a?b?c?d?0? ??8a?4b?2c?d?1??3a?2b?c?1解得：a??1,b?4,c??4,d?1。

H(x)??x3?4x2?4x?1

11 设f(x)?x,x0?1/4,x1?1,x2?9/4。(1)试求f(x)在?1/4,9/4?上的三次埃尔米特插值多项式H(x)，使得H(xj)?f(xj),j?0,1,2,H?(x1)?f?(x1)，H(x)以升幂形式给出。(2)写出余项R(x)?f(x)?H(x)的表达式。（埃尔米特插值及其余项的计算）。

321192733解：f()?，f(1)?1，f()?，f?(x)?x2，f?(1)?

484822设H(x)?ax?bx?cx?d，H?(x)?3ax?2bx?c

3221111?1a?b?c?d??641648?a?b?c?d?1???72981927 ?64a?16b?4c?d?8??3a?2b?c?3?2?解得：a??

142632331，b?，c?，d??。 225450450257

2008信息与计算科学专业计算方法习题参考解答 江世宏编

故 H(x)??14326322331x?x?x?。 225450450255193?219R(x)??(x?)(x?1)2(x?)，其中，???。 441284412 若f(x)?c2[a,b],f(a)?f(b)?0，试证明：

max|f (x)|?a?x?b1?b?a?2max|f?? (x)|（插值余项的应用）

a?x?b8解：以f(a)?f(b)?0为插值条件，作线性插值多项式，有

L(x)?其余项为

x?bx?a?f(a)??f(b)?0 a?bb?af??(?)(x?a)(x?b) 2!1a?ba?b1?a)(b?)?(b?a)2maxf??(x)。 故 maxf(x)?maxf??(x)?(a?x?ba?x?b2a?x?b228R(x)?f(x)?L(x)?f(x)?13 设f(?2)??1,f(0)?1,f(2)?2,求p(x)使p(xi)?f(xi)(i?0,1,2)； 又设 |f???(x)|?M ，则估计余项r(x)?f(x)?p(x)的大小。（插值误差的估计） 解：由插值条件，有

?4a?2b?c??1? ?c?1?4a?2b?c?2??a??1/8?解得：?b?3/4

?c?1?从而 p(x)??其余项为

123x?x?1 84r(x)?f(x)?p(x)?f???(?)(x?2)x(x?2)??(?2,2) 3!r(x)?M3M1683(x?4x)?3?M 66927 8