数值分析习题集及答案 下载本文

?1?A?3???43124??2?1??正交相似约化为对称三对角阵。 8. 设A?Rn?n

,且

Tai1,aj1不全为零,

PijPA为使aj1?0的平面旋转阵,试推导计算ij第i行,第j行元

(2)素公式及APij第i列,第j列元素的计算公式。

9. 设An?1是由豪斯荷尔德方法得到的矩阵,又设y是An?1的一个特征向量。 (a)证明矩阵A对应的特征向量是x?P1P2?Pn?2y; (b)对于给出的y应如何计算x? 10. 用带位移的QR方法计算

?1?A?2???02?10??1?3???3?B?1???0120??1?1??11(a) , (b)

全部特征值。

11. 试用初等反射阵A分解为QR,其中Q为正交阵,R为上三角阵,

?1?A?2???21?1?41???1?5??。

数值分析习题简答

(适合课程《数值方法A》和《数值方法B》)

第一章 绪论习题参考答案

?(x)1. ε(lnx)≈

?r(x)?n*xn*??r(x)??*。

?nx*n?1?(x)x*n?(x)n*2.

*x*?n?(x)x**?0.02n。

***3. x1有5位有效数字,x2有2位有效数字,x3有4位有效数字,x4有5位有效数字,x5有

*2位有效数字。

4. ?(x1?x2?x4)??(x1)??(x2)??(x4)?0.5?10*****************?4?0.5?10?3?0.5?10?3?1.05?10?3?(x1x2x3)?x2x3?(x1)?x1x3?(x2)?x1x2?(x3)?0.214790825?(x2x4**)?*1x4*?(x)?3V4?*2x2x4**2?(x4)?8.855668?10136?V2*?6。

?(V)/123?r(R)??r(3)?133V4??1?(V)3V?13?r(V)?0.0033335. 6.

?(Y100)?100?100?12?10?3??10?3。

7. x1?28?783?55.982,

x2?28?783?128?783?155.982?0.01786。

?8.

9.

??N11?xdx?2?2?arctgN

gt?(t)12gt2?(x)??(S)?12S?12?(S)?0.005?r(S)??2?(t)t?0.2t10. ?(S)?gt?(t)?0.1gt,加,相对误差减小。

11. 12.

f2?,故t增加时S的绝对误差增

?(y10)?10?(y0)?61012?108,计算过程不稳定。

f?(2?1)?0.005051,

如果令

32?1.,4则

f1?(2?1)?0.0040961(3?22)36,

1(2?1)6?0.005233,

f3?(3?22)?0.008,

f4??0.005125,

f5?99?702?1,f4的结果最好。

13.

f(30)??4.094622,开平方时用六位函数表计算所得的误差为

??12?10?4,分别代

入等价公式

?(f1)?ln(1??x?2f1(x)?ln(x?)?x?)?x?x?1),f2(x)??ln(x??(x?160x?1)??60?1222x?1)12?10?42中计算可得

?3?x?12?3?10x?1,

?7?(f2)?ln(1??x?x?12?x?12???10?4?8.33?10。

999999998999999999?1.00000014.

方程组的真解为

?ssx1?1000000000999999999?1.000000,x2?,而无论用方

程一还是方程二代入消元均解得x1?1.00,x2?1.00,结果十分可靠。

?tanc?c??bsinc?a?asinc?b?abcosc?cabsinc??aa??bb??cc15.

第二章 插值法习题参考答案

n?1Vn(x)?1.

?(x?xi?0i)?(x0?j?i?n?1i?xj);

Vn?1(x0,x1,?,xn?1)?(x?1)(x?2)(1?1)(1?2)2?(x0?j?i?n?1i?xj).

2.

L2(x)?0??56?(?3)?(x?1)(x?2)(?1?1)(?1?2)?4?(x?1)(x?1)(2?1)(2?1)

x?32x?73.

3. 线性插值:取x0?0.5,x1?0.6,y0??0.693147,y1??0.510826,则

ln0.54?L1(0.54)?y0?y1?y0x1?x0?(0.54?x0)??0.620219;

,则

二次插值:取

x0?0.4,x1?0.5,x2?0.6,y0??0.916291,y1??0.693147,y2??0.510826ln0.54?L2(0.54) ?y0?(0.54?x1)(0.54?x2)(x0?x1)(x0?x2)?y1?(0.54?x0)(0.54?x2)(x1?x0)(x1?x2)?y2?(0.54?x0)(0.54?x1)(x2?x0)(x2?x1)=-0.616707 . 4.

R1(x)?f(x)?L1(x)?|R1(x)|?1212f??(?)(x?x0)(x?x1)

,其中??[x0,x1].

所以总误差界

x0?x?x1max|cos??(x)|?max|(x?x0)(x?x1)|x0?x?x1?12?1?(x1?x0)42???1?8??????1.06?108?60180? . 125. 当

x?x0?4?37?hl2(x)?(x?x0)(x?x1)(x?x3)(x2?x0)(x2?x1)(x2?x3)

时,取得最大值

x0?x?x3max|l2(x)|?10?7727 .

k6. i) 对f(x)?x,(k?0,1,?,n)在x0,x1,?,xn处进行n次拉格朗日插值,则有

xk?Pn(x)?Rn(x)n

1f(n?1)? 由于f(n?1)?i?0lj(x)xj?nk(n?1)!kj(?)(x?x0)?(x?xn)

(?)?0,故有i?0?l(x)xj?xk.

n次拉格朗日插值,有

kn ii) 构造函数

g(x)?(x?t),在x0,x1,?,xn处进行

kLn(x)?k?(xi?0j?t)lj(x).

插值余项为

(x?t)?Ln(x)?g(n?1)(?)n(n?1)!?(x?xj?0j),

由于

g(n?1)(?)?0,(k?1,2,?,n).故有

n(x?t)nk?Ln(x)??(xi?0j?t)lj(x)k

令t?x,即得

?(xi?0j?t)lj(x)?0k.

7. 以a, b两点为插值节点作f(x)的一次插值多项式

L1(x)?f(a)?f(x)?L1(x)?12f(b)?f(a)b?a(x?a),

据余项定理,

f??(?)(x?a)(x?b),??[a,b],

18由于f(a)?f(b)?0,故

|f(x)?L1(x)|?|f(x)|?R2(x)?16?12max|f??(x)|max|(x?a)(x?b)|?a?x?ba?x?b2(b?a)max|f??(x)|.a?x?b

8. 截断误差 其中

e(x?x0)(x?x1)(x?x2),??[?4,4].

x0?x1?h,x2?x1?h,

?4?x?4则

x?x1?3329,h时取得最大值

3?h3max|(x?x0)(x?x1)(x?x2)|?16e?(4 .

由题意, 9.

?yn?2|R2(x)|?293?h)?103?6

n?1所以,h?0.006.

n?1?2,

n?1/2n

?yn?(22n?2?2n?1)?(2?2)?2,

2nn?1n?1nn则可得

?yn??(?yn)?2.

42?yn?2n?1/2?2, ?yn?(22n?1?2)?(2?2422nn)?2,则可得

?yn??(?yn)?2n?210. 数学归纳法证

当k?1时,?f(x)?f(x?h)?f(x)为m-1次多项式;

.

k假设 ?f(x)(0?k?m)是m-k 次多项式,设为g(x),则

?k?1f(x)?g(x?h)?g(x)为

m-(k+1)次多项式,得证。

11. 右?fk(gk?1?gk)?gk?1(fk?1?fk)?fk?1gk?1?fkgk?左

n?112.

?k?0n?1k?0fk?gk?f0g1?f0g0?f1g2?f1g1???fn?1gn?fn?1gn?1,?fk?f1g1?f0g1?f2g2?f1g2???fngn?fn?1gn.2

?gk?1

n?113.

??j?0yj

?(y2?y1)?(y1?y0)?(y3?y2)?(y2?y1)???(yn?1?yn)?(yn?yn?1)

?(yn?1?yn)?(y1?y0)??yn??y0 .

14. 由于x1,x2,?,xn是f(x)的n个互异的零点,所以

f(x)?a0(x?x1)(x?x2)?(x?xn)

nn?a0?(x?xi)?a0(x?xj)?(x?xi),i?1i?1i?j

对f(x)求导得

?n?nf?(x)?a0??(x?xi)?(x?xj)(?(x?xi))???i?0?i?1??i?j?i?j?,

nf?(xj)?a0?(xj?xi)则

ni?1i?j,

?1a0n?j?1xjkf?(xj)?j?1xjnk.?xi)?(xi?1i?jj

记gk(x)?x,则 由以上两式得

kg(n?1)

n?0,0?k?n?2,(x)???(n?1)!,k?n?1.

?j?1xjkf?(xj)?1a0n?j?1gk(xj)n?1a0gk[x1,x2,?,xn]?i?1i?j(xj?xi)

?1gka0n(n?1)(?) 15. i)

F[x0,x1,?,xn]?n(n?1)!?(xj?0j?0,0?k?n?2,???1?a0,k?n?1.

F(xj)j?x0)?(xj?xj?1)(xj?xj?1)?(xj?xn)

c?f(xj)?c?f[x0,x1,?,xn]? ii) 证明同上。 16.

?(xj?00?x0)?(xj?xj?1)(xj?xj?1)?(xj?xn)f(7).

f[2,2,?,2]?0117(?)7!8?f7!7!(8)?1;(?)

?0.f[2,2,?,2]?7!

17.

R3(xj)?f(xj)?p(xj)?0,

?R3(xj)?f?(xj)?p?(xj)?0,j?k,k?1.即xk,xk?1均为R3(x)的二重零点。因而有形式:

R3(x)?K(x)(x?xk)(x?xk?1).22

作辅助函数?(t)?f(t)?p(t)?K(x)(t?xk)(t?xk?1). 则 ?(xk)?0,?(x)?0,?(xk?1)?0,??(xk)?0,??(xk?1)?0. 由罗尔定理,存在?1?(xk,x),?2?(x,xk?1),使得

??(?1)?0,??(?2)?0.

22