数值分析试题与答案

一. 填空题（本大题共4小题，每小题4分，共16分）

1.设有节点012,,x x x ，其对应的函数()y f x =的值分别为012,,y y y ，则二次拉格朗日插值基函数0()l x 为 。

2.设()2f x x =，则()f x 关于节点0120,1,3x x x ===的二阶向前差分为 。

3.设110111011A -????=--????-??，233x ??

??=??????，则1A ＝ ，1x = 。

4. 1n +个节点的高斯求积公式的代数精确度为 。

二．简答题（本大题共3小题，每小题8分，共24分）

1. 哪种线性方程组可用平方根法求解？为什么说平方根法计算稳定？

2. 什么是不动点迭代法？()x ?满足什么条件才能保证不动点存在和不动点迭代序列收敛于()x ?的不动点？

3. 设n 阶矩阵A 具有n 个特征值且满足123n λλλλ>≥≥≥ ，请简单说明求解矩阵A 的主特征值和特征向量的算法及流程。

三．求一个次数不高于3的多项式()3P x ，满足下列插值条件：

i x 1 2 3 i y 2 4 12 i y '

3

并估计误差。（10分）

四．试用1,2,4n =的牛顿-科特斯求积公式计算定积分1

01

1I dx x

=+?

。（10分） 五．用Newton 法求()cos 0f x x x =-=的近似解。（10分） 六．试用Doolittle 分解法求解方程组：

1232

5610413191963630

x x x -????????????-=??????

??????----?????? （10分）

七．请写出雅可比迭代法求解线性方程组123123123202324812231530

x x x x x x x x x ++=??

++=??-+=? 的迭代格式，并

判断其是否收敛？（10分）

八．就初值问题0(0)y y

y y λ'=??=?考察欧拉显式格式的收敛性。（10分）

《数值分析》（A ）卷标准答案

（2009－2010－1）

一． 填空题（每小题3分，共12分） 1. ()1200102()()

()()

x x x x l x x x x x --=

--； 2.7；3. 3，8；4. 2n+1。

二．简答题（本大题共3小题，每小题8分，共24分）

1. 解：系数矩阵为对称正定的方程组可用平方根法。 （4分）

对于对称正定阵 A ，从2

1

i

ii ik k a l ==

∑

可知对任意k ≤ i 有||ik ii l a ≤。即 L 的元素不会增大，误差可控，不需选主元，所以稳定。 （4分） 2. 解：（1）若()*

*x

x ?=，则称*x 为函数()x ?的不动点。 （2分）

（2）()x ?必须满足下列三个条件，才能保证不动点存在和不动点迭代序列收敛于

()x ?的不动点：

1）()x ?是在其定义域内是连续函数； （2分） 2）()x ?的值域是定义域的子集； （2分） 3）()x ?在其定义域内满足李普希兹条件。 （2分） 3.解：参照幂法求解主特征值的流程 （8分）

步1：输入矩阵A ，初始向量v0，误差限ε，最大迭代次数N; 步2：置k:=1,μ:=0，u0=v0/||v0||∞; 步3：计算vk=Auk-1; 步4：计算

[][]1max ;

k k r i i n

v v ≤≤=

并置mk:=[vk]r, uk:=vk/mk;

步5：若|mk- μ |< ε，计算，输出mk,uk ；否则，转6； 步6：若k

设()2p x 满足 ()()()22212,24,312,p p p ===则()2

2376,p x x x =-+（3分）

再设()()()()()32123p x p x K x x x =+--- （3分） 2K = （1分） ()32329156p x x x x =-+- （1分）

（2）()()()()()()2

4311234!

R x f x x x ξ=

--- （2分） 四．解：应用梯形公式得()()11

012

I I f f ≈=+???? （2分） 0.75= （1分）

应用辛普森公式得：()()21104162I I f f f ?

???

≈=

++ ???????

（2分）

0.69444444= （1分）

应用科特斯公式得：

()()41113703212327190424I I f f f f f ?

???

??

??

≈=

++++ ? ? ???

????????

（2分） 0.6931746= （2分） 五．解：由零点定理，cos 0x x -=在(0,

)2

π

内有根。 （2分）

由牛顿迭代格式1cos 0,1,......1sin n n

n n n

x x x x n x +-=-=+ （4分）

取04

x π

=

得， 12340.73936133;

0.739085178

0.7390851330.739085133

x x x x ==== （3分）

故取*

40.739085133x x ≈= （1分） 六．解：对系数矩阵做三角分解：

1112

1321222331323325610

0413********u u u l u u l l u -??????

??????-=????????????---??????

（2分）

125621373414A LU -????

????=-=????????-????

（4分）

若Ly b =，则12310,1,4y y y ==-=； （2分） 若Ux y =，则(3,2,1)T x =。 （2分）

七．解：（1）对于方程组，雅可比方法的迭代矩阵为

00.50.51010.50.50B -????=--??

????

（2分）

其特征多项式为()

2

det() 1.25I B λλλ-=+，且特征值为

1230, 1.25, 1.25i i λλλ===- （2分） 故有() 1.251B ρ=>，因而雅可比迭代法不收敛。 （1分） （2）对于方程组，Gauss-Seidel 迭代法迭代矩阵为

00.50.500.50.5000.5B -????=--????-??

（2分） 其特征值为1230,0.5λλλ=== （2分） 故有()0.51B ρ=<，因而雅可比迭代法收敛。 （1分） 八．证明题（本大题共2小题，每小题7分，共14分） 1. 证：该问题的精确解为0()x

y x y e

λ= （2分）

欧拉公式为1(1)i i i i y y h y h y λλ+=+=+ （2分） 对任意固定的i x x ih ==， 有/1/00(1)[(1)]i i x h

x h i y y h y h λλλλ=+=+， （2分）

则0()i

x i y e

y x λ= （1分）

2.证：牛顿迭代格式为125,0,1,2,66n n n

x a

x n x +=

+= （3分）

因迭代函数为()25,66x a

x x

?=+而()35,63a x x ?'=+又*3x a =， （2分）

则

()

()

3

3

351

062

3a a a

?'

=

+

=

≠。 故此迭代格式是线性收敛的。 （2分）

《数值分析》参考解答

三．计算题（每小题7分，共42分）：

1. 设 x

e x

f =)(, 试构造基函数求)(x f 的2次插值多项式 )(2x P ,满足:

)1()1(),0(')0('),0()0(222f P f P f P ===.

解 设)(2x P 的基函数为)(),(),(010x x x βαα，则它们满足下列关系

（1分）

x 0 1 x 0

)(2x P 1 e )('2x P 1

)

(0x α 1 0

)('

0x α 0

)(1x α 0 1 )('1x α 0

)

(0x β0 0

)('0x β 1

（2分）

(1) 令00200)(c x b x a x ++=α，则有???

??===++===0

)0(0)1(1

)0(0'

0000000b c b a c ααα,

即1,0,1000==-=c b a . 所以1)(20+-=x x α. 或由0)1(0=α，先得))(1()(0l kx x x +-=α.

再由1)0(0=α，得1=-l ，即1-=l . 由1)0(0

='α，得0=-k l ，即1-==l k . 所以1)1)(1()(20+-=+--=x x x x α.

（1分）

(2) 令112

11)(c x b x a x ++=α，则有?????===++===0

)0(1)1(0

)0(1'

1111111b c b a c ααα,

即0,0,1111===c b a . 所以21)(x x =α.

或由0)0()0(1

1='=αα,先得21)(kx x =α. 再由1)1(1=α，得1=k . 所以21)(x x =α.

（1分）

(3) 令22220)(c x b x a x ++=β，则有???

??===++===1

)0(0)1(0

)0(2'

0222020b c b a c βββ，

即 0,1,1222==-=c b a . 所以x x x +-=20)(β 或由0)1()0(00==ββ,先得)1()(0-=x kx x β.

再由1)0(0='

β，得1=-k ，即1-=k . 所以x x x x x +-=--=20)1()(β （1分） 最后得 1)2()()0()()1()()0()(20'102++-=++=x x e x f x f x f x P βαα. （1分）

2. 求 x x x x f ++=2323)( 在区间 [-1,1] 上的２次最佳一致逼近多项式;

解 设所求的2次最佳一致逼近多项式为)(*2x P . 令)]()([3

1

)(*2x P x f x Q -=. （2分）

则)(x Q 的首项系数为1, 并且当)(2

1)]()([31)(32*

2x T x P x f x Q =-=时, )(x Q 与0的偏

差最小, 即)(x f 与)(*

2x P 的偏差最小.

（2分） 因为]1,1[-上的3次切比雪夫Chebyshev 多项式为x x x T 34)(33-=.

（1分） 所以x x x x x x x x T x f x P 4

13

2)493(23)(43)()(23233*

2+=--++=-

=. （2分）

3．利用龙贝格公式计算定积分(计算到1R 即可)：

?

-+7

1

.2dx x

解 ]7,1[,2)(-∈+=x x x f ，16)]7()1([2

)

1(71=+---=

f f T , （1分）

94428.16548)3(28

2112=?+=+=f T T ,

22774.17)73(247214.8)]5()1([24

2124=+?+=++=f f T T ,

)]6()4()2()0([2

2

2148f f f f T T ++++=

30599.178********.8=++++=, （2分）

25904.173134121=-=T T S ,32223.1731

34242=-=T T S ,

33207.173

1

34484=-=T T S , （2分） 32644.1715

11516121=-=S S C , T S C R n=1 16 17.25904 17.32644 17.33283 n=2 16.94428 17.32223 17.33273 n=4 17.22774 17.33207 n=8

17.30599

33273.17151

1516242=-=

S S C , 33283.1763

1

6364121=-=C C R .(2分)

4．利用改进的尤拉方法求解常微分方程初值问题：（要求取步长2.0=h 计算）

解 令x y y x f +=),(，则改进的尤拉公式为：

))],(,(),([2

1n n n n n n n n y x hf y h x f y x f h

y y ++++=+

（2分）

22)2(2)2(12h x h h y h h n n +++??

?

???++=. （2分)

取2.0=h 得，02.022.022.11++=+n n n x y y . （1分)

计算结果如下：

x

y 1 1 1.2 1.46 1.4 2.0652 1.6

2.84754

（2分)

5．用牛顿法求方程 023)(3

=--=x x x f 在 30=x 附近的根（只要求迭代２步）。

解 牛顿迭代公式为：)

()

('

1n n n n x f x f x x -

=+ （2分)

)

1(32

2

----

=n n n n x x x x .

（2分)

取迭代初值为30=x ，则迭代结果如下表所示：

（3分） 6．写出解如下线性方程组的高斯－塞德尔迭代公式，并讨论其收敛性。如果不收敛，则应怎样处理才能得到收敛的高斯－塞德尔迭代公式？

n

n x

0 3

1

2.33333

2 2.05555

??

?=≤≤+='.

1)1()6.11(y x x y y

??

?=+=-.

123,

0322121x x x x

解 2332A -??=????,2002D ??=????,0030L ??=-????,0300U -??

=-????,??

????=10b . （1分）

则 ()

1

1

2020132324D L --????

-==????-????

, （1分）

得 ()1

2003061132000944G D L U -??????

=-=

=??????--??????

, （1分）

()

1

200011321244f D L b -??????

=-==??????-??????

,

（1分）

(

)

()()

1600419024k k k X GX f X +?

???????=+=+??????

-??????

为高斯-塞德尔迭代公式.

（1分）

这时G 的2个特征值为129

,04

λλ=-=,故()1G ρ>，迭代法不收敛. （1分）

若原方程 1212230321x x x x -=??

+=?改写成为 1212321230

x x x x +=??-=?, 这时3223A ??

=??

-??是严格对角 优势矩阵，则由此得到的迭代法必收敛.

（1分）

四．证明题(每小题9分,共18分)：

1. 证明本试卷第三大题（即计算题）第１小题的插值余项：

)1(6)()()(2

22-=-=x x e x P x f x R ξ)10(<<ξ, 并有误差估计.81|)(|2≤x R

证：方法一：因为()()22()R x f x P x =-，则0,1是()2R x 的零点且0为二重的,

(1分) 于是可设()()2

2()1R x k x x x =??-，令[]2

2()()()()(1

),0,1t f t P t k x t t t φ=---∈ (2分) 则()t φ有4个零点:0,0,1,x ,连续使用三次罗尔定理,则()0,1,ξ?∈使'''

()0φξ=,

(2

分)

即''''''

()()()3!0()3!6

f e f k x k x ξ

ξξ-?=?=

=, 得()()2216e R x x x ξ=??-. (2

分)

方法二: 设2

222

()()()()()(1)(1)

f x P x t f t P t t t x x φ-=----, 则()t φ有3个零点0,1,x , (1分)

'()t φ有2+1个零点,。'''()t φ有一个零点，所以''''''22()()

0()()3!(1)

f x P x f x x φξξ-==-

- (2

分)

'''22()()

()3!(1)

f x P x f x x ξ-=

-

(2

分)

''2

2211()()()(1)(1)66f x P x f x x x e x x ξ-=-=-， 即()()2216

e R x x x ξ=??-.

(2

分)

最后()()()2

2223(1)1||11166628

e e x x R x x x x x x ξ+-????=-≤-≤??

分)

2．证明: 求积公式

)5

3

(95)0(98)53(95)(1

1

f f f dx x f ++-≈

?

-恰有5次代数精度. 证：当()1f x =时，

1

1

()f x dx -?

11

1()2d x -==?，

53853585

()(0)()1(0)1295995999

f f f f -++=?++?=;

(1分)

当()f x x =时， 1

21

1

11

1

()()02x f x dx xd x ---??

===??????,

5385353853

()(0)()()(0)()09599595995

f f f f -++=-++=; (1分)

当2

()f x x =时，1

31

1

211

12()()33

x f x dx x d x ---??===??????,

2253853538532

()(0)()()(0)()95995959953

f f f f -++=-++=; (1分)

当3

()f x x =时，

1

1

31

1

()()0f x dx x d x --==?

?,

335385353853

()(0)()()(0)()09599595995

f f f f -++=-++=; (1分)

当4

()f x x =时， 1

51

1

411

12()()55

x f x dx x d x ---??===??????,

4453853538532()(0)()()(0)()95995959955

f f f f -++=-++=; (1

分)

当5()f x x =时，

1

1

51

1

()()0f x dx x d x --==?

?,

555385353853

()(0)()()(0)()09599595995

f f f f -++=-++=. (1

分)

即求积公式对次数不超过5的多项式准确成立, 但当6

()f x x =时，

1

71

1

611

1

2()()77x f x dx x d x ---??===??????, 6653853538536

()(0)()()(0)()959959599525

f f f f -++=-++=, 不成立. (2

分)

综之，求积公式具有5次代数精度.

(1

分)

数值分析试题1

1. 已知325413

.0,325413*

2*

1==X X 都有6位有效数字，求绝对误差限。（4分） 解：

由已知可知,n=6

5.0102

1

,0,6,10325413.0016*1=?=

=-=?=ε绝对误差限n k k X 2分

620*2102

1

,6,0,10325413.0-?=

-=-=?=ε绝对误差限n k k X 2分

2. 已知?????=001A 220

- ????

?440求21,,A A A ∞ （6分） 解：

{},88,4,1max 1==A 1分 {},66,6,1max ==∞A 1分

()A A A

T max 2

λ= 1分

?????=001A A T 420 ??

???

-420?????001 220

- ?????440=?????001 080 ????

?3200 2分 {}3232,8,1max )(max ==A A T λ 1分

24322==A

3. 设3

2

)()(a x x f -= （6分） ① 写出f(x)=0解的Newton 迭代格式

② 当a 为何值时，)(1k k x x ?=+ （k=0,1……）产生的序列{}k x 收敛于2

解：

①Newton 迭代格式为：

x

a x x x a

x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2

2

32

1

+=

+=---=-=+? 3分

②时迭代收敛即当222,112

10)2(',665)('2<<-<-=-=a a x a x ?? 3分

4. 给定线性方程组Ax=b ，其中：???=13A ???22，??

?

???-=13b 用迭代公式

)()()()1(k k k Ax b x x -+=+α（k=0,1……）求解Ax=b ，问取什么实数α，可使迭代收

敛 （8分）

解：

所给迭代公式的迭代矩阵为?

?

?

--?

??--=-=ααααα21231A I B 2分 其特征方程为

0)

21(2)31(=----=

-αλα

ααλλB I 2分

即，解得αλαλ41,121-=-= 2分 要使其满足题意，须使1)(

5. 设方程Ax=b ，其中?????=211A 212 ????

?-112，??

????????=765b 试讨论解此方程的Jacobi 迭代法的收

敛性，并建立Gauss-Seidel 迭代格式 （9分）

解：

U D L A ++=

??

???--=+-=-210

)(1U L D B J 202-- ?

????-012 3分 0,03213=====-λλλλλJ B I 2分

即10)(<=J B ρ，由此可知Jacobi 迭代收敛 1分 Gauss-Seidel 迭代格式：

?????--=--=+-=++++++)

1(2

)1(1)1(3

)

(3

)1(1)

1(2)

(3)(2)1(12276225k k k k k k k k k x x x x x x x x x （k=0,1,2,3……） 3分

6. 用Doolittle 分解计算下列3个线性代数方程组：i i b Ax =（i=1,2,3）其中

?????=222A 331 ?????421，23121,,974x b x b b ==????

?

?????= （12分）

解：

①11b Ax =

?????222 331 ?????421????

??????=9741x A=?????111 110 ?????100?????002 021 ????

?

211=LU 3分 由Ly=b1，即?????111 110 ?????100y=??????????974 得y=??

????????234 1分 由Ux1=y ，即?????002 021 ?????211x1=??????????234 得x1=????

?

?????111 2分 ②22b Ax =

?????222 331 ?????421x2=????

?

?????111 由Ly=b2=x1，即?????111 110 ?????100y=??????????111 得y=??

???

?????001 1分 由Ux2=y ，即?????002 021 ?????211x2=??????????001 得x2=??

??

?

?????005.0 2分

③33b Ax =

?????222 331 ?????421x3=????

?

?????005.0 由Ly=b3=x2，即?????111 110 ?????100y=??????????005.0 得y=??

?

??

?????-05.05.0 1分

由Ux3=y ，即?????002 021 ?????211x3=??????????-05.05.0 得x3=???

?

?

?????-025.0375.0 2分

7. 已知函数y=f(x)有关数据如下：

要求一次数不超过3的H 插值多项式，使

'11'

33)(,)(y x H y x H i i == （6分）

解：

作重点的差分表，如下：

3分

2

1021101011001003))(](,,,[))(](,,[)](,[][)(x x x x x x x x f x x x x x x x f x x x x f x f x H --+--+-+= =-1+(x+1)-x(x+1)+2x.x(x+1)

=2

3

2x x + 3分

8. 有如下函数表：

试计算此列表函数的差分表，并利用Newton 前插公式给出它的插值多项式 （7分）

解：

由已知条件可作差分表，

3分

i ih x x i =+=0 （i=0,1,2,3）为等距插值节点，则Newton 向前插值公式为： 03

3

2100

22100003!3))()((!2))((!1)()(f h x x x x x x f h x x x x f h x x f x N ?---+?--+?-+

= =4+5x+x(x-1)

=442

++x x 4分

9. 求f(x)=x 在[-1,1]上的二次最佳平方逼近多项式)(2x P ，并求出平方误差 （8分）

解：

令2

2102)(x a x a a x P ++= 2分

取m=1, n=x, k=2

x ,计算得： (m,m)=dx ?

-1

11=0 (m,n)=

dx x ?-11

=1 (m,k)= dx x

?-1

1

2

=0

(n,k)= dx x ?-1

13=0.5 (k,k)= dx x ?

-1

1

4=0 (m,y)= dx x ?-1

1

=1

(n,y)=

dx x ?

-1

1

2=0 (k,y)=

dx x ?

-1

1

3=0.5

得方程组：??

?

??==+=5

.05.005.011201a a a a 3分

解之得c a a c a 2,1,210-=== （c 为任意实数，且不为零）

即二次最佳平方逼近多项式2

22)(cx x c x P -+= 1分 平方误差：3

2

),(2

2222

2

22

=

-=-=∑=i i i y a f

p f ?δ

2分

10. 已知如下数据：用复合梯形公式，复合Simpson 公式计算?+=

1

0214

dx x π的近似值（保

留小数点后三位） （8分）

解：

用复合梯形公式： )}1()]8

7

()43()85()21()83()41()81([2)0({1618f f f f f f f f f T ++++++++=

=3.139 4分

用复合Simpson 公式： )}1()]4

3

()21()41([2)]87()85()83()81([4)0({2414f f f f f f f f f S ++++++++=

=3.142 4分

11. 计算积分?

=

20

sin π

xdx I ，若用复合Simpson 公式要使误差不超过5102

1

-?，问区间

]2,0[π要分为多少等分?若改用复合梯形公式达到同样精确度，区间]2

,0[π

应分为多少等

分? （10分）

解： ①由Simpson 公式余项及x x f x x f sin )(,sin )()4(==得

544)4(2

04102

1)1()4(360)(max )4(1802)(-≤≤?≤=≤n x f n f R x n ππππ

π 2分

即08.5,6654≥≥n n ，取n=6 2分

即区间]2,0[π分为12等分可使误差不超过5102

1

-? 1分

②对梯形公式同样1)(''max 2

0≤≤

≤x f x π

，由余项公式得

5102

1

)2(122)(-?≤≤n f R n ππ

2分

即255,2.254=≥n n 取 2分

即区间]2,0[π分为510等分可使误差不超过5102

1

-? 1分

12. 用改进Euler 格式求解初值问题：?

?

?==++1)1(0

sin 2'y x y y y 要求取步长h 为0.1，计算y （1.1）的近似值 （保留小数点后三位）[提示：sin1=0.84,sin1.1=0.89] （6分）

解：

改进Euler 格式为：

??

???++=+=+-++-

+)]

,(),([2)

,(1111n n n n n n n n n n y x f y x f h

y y y x hf y y 2分 于是有

?????+++-=+-=+-++-

+-+)

sin sin (05.0)

sin (1.0121121

2

1n n n n n n n n n n n n n x y y x y y y y x y y y y （n=0,1,2……） 2分 由y(1)=0y =1，计算得

????

?=≈=+-=-

838

.0)1.1(816

.0)1sin 11(1.01121y y y 2分 即y(1.1)的近似值为0.838

13. ]

[],[],,[lim ],[),,(],,[)(0'

000000'

x f x x f x x f x x f b a x b a C x f x x ==∈∈→证明：定义：设（4分）

证明：

]

['],[],[],[lim ]

[][lim

]['00000000000

x f x x f x x f x x f x x x f x f x f x x x x ===--=→→故可证出 4分

14. 证明：设n

n R

A ?∈，?为任意矩阵范数，则A A ≤)(ρ （6分）

证明：

设λ为A 的按模最大特征值，x 为相对应的特征向量，则有Ax=λx 1分

且λρ=)(A ，若λ是实数，则x 也是实数，得Ax x =λ 1分

而

x x ?=λλ x A x ,?≤??≤λ故x A Ax 2分

由于A x 0x ≤≠λ得到，两边除以 1分 故A A ≤)(ρ 1分 当λ是复数时，一般来说x 也是复数，上述结论依旧成立

数值分析试题2

1、（本题5分）试确定7

22

作为π的近似值具有几位有效数字，并确定其相对误差限。 解 因为 7

22

=3.142857…=1103142857

.0-?

π=3.141592… 所以

312102

11021005.0001264.0722--?=?=<=-

π （2分） 这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7

22

作为π的近似值具有3位有效数字。 （1分） 而相对误差限

3102

1

0005.00004138.0001264.07

22-?=

<≈=

-=

π

π

πε

r （2分） 2、（本题6分）用改进平方根法解方程组：???

?? ??=????? ??????? ??--654131*********x x x ；

解 设???

?

?

??????? ?

?????? ??===????? ??--11111

1

13132111232312132

1

32

3121l l l d d d l l l LDL A T

由矩阵乘法得：

5

7,21,215

27

,25,2323121321-

==-==

-==l l l d d d （3分）

由y D x L b Ly T

1

,-==解得

T T x y )9

23

,97,910(,)563,

7,4(== （3分） 3、（本题6分）给定线性方程组???????=++-=+-+=-+-=-+17

7222382311387

510432143213

21431x x x x x x x x x x x x x x

1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式；

2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性； 解 1）Jacoib 迭代格式为

???????-+-=----=+-=+--=++++7

)2217()8()2323(8)311(10

)57()(3)(2)(1)1(4

)(4)(2)(1)

1(3

)

(3)(1)1(2

)

(4)(3)1(1k k k k k k k k k k k k k k x x x x x x x x x x x x x x （2分）

Gauss-Seidel 迭代格式为

???????-+-=----=+-=+--=++++++++++7

)2217()

8()2323(8

)311(10

)57()1(3

)1(2

)1(1

)1(4

)(4)1(2)

1(1)1(3

)

(3)1(1)1(2)(4)(3)1(1k k k k k k k k k k k k k k x

x

x

x x x

x x x x x x x x （2分）

2）由于所给线性方程组的系数矩阵

????

??

? ??----=72211

8230381510

10A 是严格对角占优的，所以Jacoib 迭代格式和Gauss-Seidel 迭代格式均是收敛的。（2分）

4、（本题6分）已知方程

08.023=--x x

在5.10=x 附近有一个根。将此方程改写成如下2个等价形式：

8.0,8.0332-=

+=x x x x

构造如下两个迭代格式：

1） ,2,1,0,8.032

1=+=+k x x k k 2） ,2,1,0,8.03

1=-=

+k x x k k

判断这两个迭代格式是否收敛；

解 1）记328.0)(x x +=?，则3

22)

8.0(32)('x x

x +=

?， 14755.005.31

)5.18.0(1)5.18.0(35.12)5.1('3

2322322<==+=+?=

? （2分）

所以该迭代格式是局部收敛的。 （1分） 2）记8.0)(3

-=

x x ?，则8

.023)('3

2-=

x x x ?，

1103.28

.05.125.13)5.1('3

2>=-?=

? （2分）

所以该迭代格式是发散的 （1分） 5、（本题6分）设23)()(a x x f -= （1）写出解0)(=x f 的牛顿迭代格式； （2）证明此迭代格式是线性收敛的。

解 （1）因23)()(a x x f -=，故)(6)('32a x x x f -=，由牛顿迭代公式 )

(')

(1n n k k x f x f x x -

=+， ,1,0=k （1分）

得k

k k k k k k x a

x a x x a x x x 665)(6)(32231

+

=---=+， ,1,0=k （2分） （2）因迭代函数2665)(x a

x x +=?， 3365)('x

a

x -=?， （1分）

3*

a x = 故02

1

)(365)('33*

≠=-=

a a x ? 此牛顿迭代格式是线性收敛的。 （2分）

6、（本题9分）给定数据

（1） 写出)(x f 的3次Lagrange 插值多项式)(3x L ； （2） 写出)(x f 的3次Newton 插值多项式)(3x N ； 解 （1）由题意知5,3,2,03210====x x x x

2)(,4)(,3)(,1)(3210=-=-==x f x f x f x f +------=))()(()

)()(()

()(30201032103x x x x x x x x x x x x x f x L

+------)

)()(()

)()(()

(3121013201x x x x x x x x x x x x x f

x 0 2 3 5 f(x) 1 -3 -4 2