当前位置:文档之家› 数值分析第六章课后习题答案

数值分析第六章课后习题答案

数值分析第六章课后习题答案
数值分析第六章课后习题答案

第六章课后习题解答

(1)()()123(1)()213(1)()()312(0

1.21125551154213351010(1,1,1),17( 4.0000186,

2.99999k k k k k k k k k T

x x x x x x x x x x x

+++ì??=---??????=-+í???

??=-++????==-(17)

解:(a )因系数矩阵按行严格对角占优,故雅可比法与高斯-塞德尔均收敛。(b )雅可比法的迭代格式为取迭代到次达到精度要求

(1)()()

123(1)(1)()

213(1)(1)(1)

312

(0)

(8)

15,2.0000012)

2112555

115

4213

3510

10

(1,1,1),8( 4.0000186,2.9999915,2.0000012)

T

k k k k k k k k k T

T

x x x x x x x x x x

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

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

??=-+

+

????==-高斯塞德尔法的迭代格式为x 取迭代到次达到精度要求

1

2

1

2:00.40.4.0.40

0.80.40.80||(0.8)(0.80.32)

()1.09282031,00.4

0.4()

00.160.6400.032

0.672D

L U I B

D L U l l l l

-

-骣--÷

?÷?÷?

÷?÷=+=

--

?÷?÷÷?

÷?÷--÷

?桫

-=-+-

=>-?

--

????=-=-????è

l

J

J

J

S

解(a )雅可比法的迭代矩阵

B

()B

B

故雅可比迭代法不收敛

高斯塞德尔法迭代矩阵

1

3

1

()||||0.810

2210122

0||02

2

02

300

2S

J

B

D

L U I B

D L U

l l

-

-?÷÷÷÷

÷÷÷÷÷÷?÷??<

骣-÷

?÷?÷?

÷?÷=+=

--?÷?÷÷?

÷?÷--?÷

-=骣

-÷?÷?÷?÷

?÷=-=

-?÷?÷÷?÷?÷?桫

l

l

S

J

J S

B

故高斯-塞德尔迭代法收敛。(b )雅可比法的迭代矩阵

B

(), (B )=0〈1故雅可比迭代法收敛。高斯-塞德尔法的迭代矩阵

B

()2

||2S

I B

l l l ÷

-=-

l

S (), (B )=2 〉1故高斯-塞德尔法不收敛。

(1lim lim lim ,(0,...0,1,0,...0)

(1,2,...,)(k k k T

i

k k i i

i n A x a

. ?萎

===(k )

k i j i j K k k K K n

K i 3:证明

必要条件:由A =A ,知a =a ,从而有||A -A ||0(K

).

故对任意的x ,有||A x -A x ||||A -A || ||x ||0 (k )

即 A x A x , A x =A x.

充分条件: 对任意的x R 有A x A x (k ),取x )

()()212()

,,...,)

()

(,,...,)(1,2,...,;1,2,...),lim .

k k T

i n i i T

i i i n i k ji

ji k k

k a a A x k

A x a a a a a j n i n A A A A

=?

=?

故即14.()12(),10.84||0.296a A A ====l V V V J 23解:不一定,因其谱半径B 不一定小于。

对习题对称,又〉0,〉0,〉0,故A 正定,但其雅可比迭代法不收敛

5.解答见例6-4

()()()

()

1123(1)()(1)()

()

22123(1)()(1)(1)

3312

3()1221()

55

511(5)42

313()

10510k k k k k k k k k k k k k k x x w x x x x x w x x x x x w x x x +++++ì??=+----??????=++--í???

??=+-+-????

(k +1)16.解:S O R 迭代格式为

(1,1,1),T

x

=取初始值计算如表.

(8)

(7)

4

||||0.00037710

,

x

x

-¥-=<因

||1|1()|120()

22,0,||1,()1,.

()

w A w A w w B A m l l m b

l -<<

<

<<

<

22

3

23

2

18.:1,2

1d e t 10,d e t ()(1)(12)01.000d e t ()32()(2)

(J

J J a a

a A a a a

A a a

B a a a a

a a I B a

a a a

a a a a B l l l l l l l l -<<骣÷?÷?=->=-+>÷?÷?÷

骣--÷

?÷?÷?÷?÷=--?÷?÷÷?÷?÷

--?桫

÷?÷?÷?

÷?÷-==-+=-+?÷?÷÷?÷?÷?÷桫

l 证明当时由

故是正定的又雅可比法迭代矩阵

故11

11)|2|,,.

22

9.:0

1,1,,()0,00i i

r

i i i n n n r a a G J P G J J n n

a J x G x g x g G l l =

=′=-

<<

÷骣?÷?÷?÷

÷??÷÷

??÷÷?÷?=÷?÷?÷÷??÷÷?÷?÷?÷?÷

÷?桫

÷

?桫

=+==?

O

O

-1

故当时雅可比迭代法收敛证明相似与它的若当标准行J ,即存在可逆阵P ,使 G =P 由于的特征值全为零,故J 一定有如下形式J =方程组等价于(I -G )由于故*

*

*

(1)

*

()

*

()

*

(0)

*

1

()*

()

*

()10,

.()

()

,k k n n n

n

n n G x g x x

G x

g

x x

G x

x x

x G x

x G P J P x

x

o x x

l +-= ==+-=--=-==-==(I -G )=12-从而I -G 非奇异,即(I -G )有唯一解于是与所述迭代格式相减,有故又故即

1

23(1)()2

131.2112

555115k k k k k k x x x x x x ++ì??=---??????=-+í解:(a)因系数矩阵按行严格对角占优,故雅可比法与高斯-塞德尔均收敛。

(b)雅可比法的迭代格式为

数值分析课后答案

1、解:将)(x V n 按最后一行展开,即知)(x V n 是n 次多项式。 由于 n i i i n n n n n i n x x x x x x x x x x V ...1...1... ......... ...... 1 )(21110 20 0---= ,.1,...,1,0-=n i 故知0)(=i n x V ,即110,...,,-n x x x 是)(x V n 的根。又)(x V n 的最高 次幂 n x 的系数为 )(...1...1... ...... .........1),...,,(101 1 21 11 2 2221 02001101j n i j i n n n n n n n n n n n x x x x x x x x x x x x x x V -== ∏-≤<≤-----------。 故知).)...()()(,...,,()(1101101------=n n n n x x x x x x x x x V x V 6、解:(1)设 .)(k x x f =当n k ,...,1,0=时,有.0)()1(=+x f n 对 )(x f 构造Lagrange 插值多项式, ),()(0 x l x x L j n j k j n ∑== 其 0)()! 1() ()()()(1)1(=+=-=++x w n f x L x F x R n n n n ξ, ξ介于j x 之间,.,...,1,0n j = 故 ),()(x L x f n =即 .,...,1,0,)(0 n k x x l x k j n j k j ==∑= 特别地,当0=k 时, 10) (=∑=n j x j l 。 (2) 0)()1(1) ()1()()(0000=-=??? ? ??-??? ? ??-=--=-===∑∑∑∑k j j i j i k j k i i j i i k j n j k i i j k n j j x x x x i k x l x x i k x l x x )利用(。 7、证明:以b a ,为节点进行线性插值,得 )()()(1 b f a b a x a f b a b x x P --+--= 因 0)()(==b f a f ,故0)(1=x P 。而 ))()(("2 1 )()(1b x a x f x P x f --= -ξ,b a <<ξ。 故)("max )(8 122)("max )(max 2 2 x f a b a b x f x f b x a b x a b x a ≤≤≤≤≤≤-=??? ??-≤。 14、解:设 ))...()(()(21n n x x x x x x a x f ---=, k x x g =)(,记)() (1 ∏=-=n j j n x x x w ,则 ),()(x w a x f n n =).()(' j n n j x w a x f = 由差商的性质知 [])! 1()(1,..,,1) (' 1 )(')('1 211 11 -== ==-===∑∑∑ n g a x x x g a x w x a x w a x x f x n n n n n j j n k j n n j j n n k j n j j k j ξ, ξ介于n x x ,...,1之间。 当20-≤≤ n k 时,0)()1(=-ξn g , 当 1-=n k 时,)!1()(1-=-n g n ξ, 故 ???-=-≤≤=-= --=∑1,,20,0)!1()(1) ('1 11 n k a n k n g a x f x n n n n j j k j ξ 16、解:根据差商与微商的关系,有 [] 1! 7! 7!7)(2,...,2,2)7(7 10===ξf f , [ ] 0! 80 !8)(2,...,2,2)8(8 1 ===ξf f 。 ( 13)(47+++=x x x x f 是7次多项式, 故 ,!7)()7(=x f 0)()8(=x f )。 25、解:(1) 右边= [][]dx x S x f x S dx x S x f b a b a ??-+-)(")(")("2)(")("2 = [] d x x S x f x S x S x S x f x f b a ?-++-)("2)(")("2)(")(")("2)(" 222 = [] d x x S x f b a ?-)(")(" 22 = [][]dx x S dx x f b a b a 2 2 )(")("??- =左边。 (2)左边= ? -b a dx x S x f x S ))(")(")(("

数值分析课后题答案

数值分析 第二章 2.当1,1,2x =-时,()0,3,4f x =-,求()f x 的二次插值多项式。 解: 0120121200102021101201220211,1,2, ()0,()3,()4;()()1 ()(1)(2)()()2()()1 ()(1)(2) ()()6 ()()1 ()(1)(1) ()()3 x x x f x f x f x x x x x l x x x x x x x x x x x l x x x x x x x x x x x l x x x x x x x ==-===-=--==-+-----==------= =-+-- 则二次拉格朗日插值多项式为 2 20 ()()k k k L x y l x ==∑ 0223()4() 14 (1)(2)(1)(1)23 537623 l x l x x x x x x x =-+=---+ -+= +- 6.设,0,1,,j x j n =L 为互异节点,求证: (1) 0()n k k j j j x l x x =≡∑ (0,1,,);k n =L (2)0 ()()0n k j j j x x l x =-≡∑ (0,1,,);k n =L 证明 (1) 令()k f x x = 若插值节点为,0,1,,j x j n =L ,则函数()f x 的n 次插值多项式为0 ()()n k n j j j L x x l x == ∑。 插值余项为(1)1() ()()()()(1)! n n n n f R x f x L x x n ξω++=-= + 又,k n ≤Q

(1)()0 ()0 n n f R x ξ+∴=∴= 0()n k k j j j x l x x =∴=∑ (0,1,,);k n =L 0 000 (2)()() (())()()(()) n k j j j n n j i k i k j j j i n n i k i i k j j i j x x l x C x x l x C x x l x =-==-==-=-=-∑∑∑∑∑ 0i n ≤≤Q 又 由上题结论可知 ()n k i j j j x l x x ==∑ ()()0 n i k i i k i k C x x x x -=∴=-=-=∑原式 ∴得证。 7设[]2 (),f x C a b ∈且()()0,f a f b ==求证: 21 max ()()max ().8 a x b a x b f x b a f x ≤≤≤≤''≤- 解:令01,x a x b ==,以此为插值节点,则线性插值多项式为 10 101010 ()() ()x x x x L x f x f x x x x x --=+-- =() () x b x a f a f b a b x a --=+-- 1()()0()0 f a f b L x ==∴=Q 又 插值余项为1011 ()()()()()()2 R x f x L x f x x x x x ''=-= -- 011 ()()()()2 f x f x x x x x ''∴= --

最新第六章习题答案-数值分析

第六章习题解答 2、利用梯形公式和Simpson 公式求积分2 1 ln xdx ? 的近似值,并估计两种方法计算值的最大 误差限。 解:①由梯形公式: 21ln 2 ()[()()][ln1ln 2]0.3466222 b a T f f a f b --= +=+=≈ 最大误差限 3''2 ()111 ()()0.0833******** T b a R f f ηη-=-=≤=≈ 其中,(1,2)η∈ ②由梯形公式: 13()[()4()()][ln14ln()ln 2]0.38586262 b a b a S f f a f f b -+= ++=++≈ 最大误差限 5(4)4()66 ()()0.0021288028802880 S b a R f f ηη-=-=≤≈, 其中,(1,2)η∈。 4、推导中点求积公式 3''()()()()() ()224 b a a b b a f x dx b a f f a b ξξ+-=-+<

数值分析习题集及答案[1].(优选)

数值分析习题集 (适合课程《数值方法A 》和《数值方法B 》) 长沙理工大学 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2%,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出 它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式(3.3)求下列各近似值的误差限: ********12412324(),(),()/,i x x x ii x x x iii x x ++其中**** 1234 ,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1%,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 1n n Y Y -=( n=1,2,…) 计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字27.982). 8. 当N 充分大时,怎样求2 1 1N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对 误差增加,而相对误差却减小. 11. 序列 {}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字), 计算到 10y 时误差有多大?这个计算过程稳定吗? 12. 计算6 1)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若

数值分析课后题答案

数值分析 2?当x=1,—1,2时,f(x)=O, 一3,4,求f(x)的二次插值多项式。解: X 0 =1,x j = — 1,x 2 = 2, f(X。)= 0, f (xj = -3, f (x2)= 4; l o(x)=(x-xi^~x2\=-1(x 1)(x-2) (x o -X/X o _x2) 2 (x -x0)(x -x2) 1 l i(x) 0 2(x-1)(x-2) (x i ~x0)(x i ~x2) 6 (x—x0)(x—x,) 1 l2(x) 0 1(x-1)(x 1) (X2 -X°)(X2 - X i) 3 则二次拉格朗日插值多项式为 2 L 2(X)= ' y k 1 k ( x) kz0 = -3l°(x) 4l2(x) 1 4 =(x_1)(x—2) 4 (x-1)(x 1) 2 3 5 2 3 7 x x - 6 2 3 6?设Xj, j =0,1,||(,n 为互异节点,求证: n (1 )7 x:l j(x) =x k(k =0,1川,n); j=0 n (2 )7 (X j -x)k l j(x)三0 (k =0,1川,n); j £ 证明 (1)令f(x)=x k

n 若插值节点为X j, j =0,1,|l(, n,则函数f (x)的n次插值多项式为L n(x)八x k l j(x)。 j=0 f (n 十)(?) 插值余项为R n(X)二f(X)-L n(X) n1(X) (n +1)!

.f(n1)( ^0 R n(X)=O n 二瓦x k l j(x) =x k(k =0,1川,n); j :o n ⑵、(X j -x)k l j(x) j卫 n n =為(' C?x j(—x)k_L)l j(x) j =0 i =0 n n i k i i =為C k( -x) (、X j l j(x)) i =0 j=0 又70 _i _n 由上题结论可知 n .原式二''C k(-x)k_L x' i=0 =(X -X)k =0 -得证。 7设f (x) c2 la,b 1且f (a) =f (b)二0,求证: max f(x)兰一(b-a) max a $至小一*丘f (x). 解:令x^a,x^b,以此为插值节点,则线性插值多项式为 L i(x^ f(x o) x x f (xj X o —人x -X o X —X o x-b x-a ==f(a) f(b)- a - b x -a 又T f (a) = f (b)二0 L i(x) = 0 1 插值余项为R(x)二f (x) - L,(x) f (x)(x - X Q)(X - xj 1 f(x) = 2 f (x)(x -X g)(X -xj

数值分析习题与答案

第一章绪论 习题一?1.设x>0,x*的相对误差为δ,求f(x)=ln x的误差限。 解:求lnx的误差极限就是求f(x)=lnx的误差限,由公式(1.2.4)有 已知x*的相对误差满足,而 ,故 即 2.下列各数都是经过四舍五入得到的近似值,试指出它们有几位有效数字,并给出其误差限与相对误差限。 解:直接根据定义和式(1.2.2)(1.2.3)则得?有5位有效数字,其误差限,相对误差限 有2位有效数字, 有5位有效数字, 3.下列公式如何才比较准确? (1)?(2) 解:要使计算较准确,主要是避免两相近数相减,故应变换所给公式。

(1)?(2) 4.近似数x*=0.0310,是 3 位有数数字。 5.计算取,利用 :式计算误差最小。 四个选项: 第二、三章插值与函数逼近 习题二、三 1. 给定的数值表 用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解:仍可使用n=1及n=2的Lagrange插值或Newto n插值,并应用误差估计(5.8)。线性插值时,用0.5及0.6两点,用Newton插值??误差限 ,因,

故? 二次插值时,用0.5,0.6,0.7三点,作二次Newton插值 ?误差限,故? 2. 在-4≤x≤4上给出的等距节点函数表,若用二次插值法求的近似值,要使误差不超过,函数表的步长h应取多少? 解:用误差估计式(5.8), ?令 因?得 3. 若,求和.

解:由均差与导数关系 ?于是 4. 若互异,求 的值,这里p≤n+1. 解:,由均差对称性 可知当有?而当P=n +1时 ?于是得 5. 求证. 解:解:只要按差分定义直接展开得 ? 6. 已知的函数表

数值分析习题集及答案

(适合课程《数值方法A 》和《数值方法B 》) 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2%,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位 有效数字: ***** 123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式(3.3)求下列各近似值的误差限: * * * * * * * * 12412324(),(),()/,i x x x ii x x x iii x x ++其中* * * * 1234,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1%,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 11783 100 n n Y Y -=- ( n=1,2,…) 计算到100Y .若取783≈27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字(783≈27.982). 8. 当N 充分大时,怎样求 2 11N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设2 12S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加, 而相对误差却减小. 11. 序列{}n y 满足递推关系1101 n n y y -=-(n=1,2,…),若02 1.41y =≈(三位有效数字),计算到10 y 时误差有多大?这个计算过程稳定吗? 12. 计算6 (21)f =-,取 2 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 6 3 11,(322), ,9970 2. (21) (322) --++ 13. 2 ()ln(1)f x x x =- -,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等 价公式 2 2 ln(1)ln(1)x x x x - -=-+ + 计算,求对数时误差有多大? 14. 试用消元法解方程组{ 10 10 12121010; 2. x x x x +=+=假定只用三位数计算,问结果是否可靠? 15. 已知三角形面积 1sin , 2 s ab c = 其中c 为弧度, 02c π << ,且测量a ,b ,c 的误差分别为,,.a b c ???证 明面积的误差s ?满足 . s a b c s a b c ????≤ ++ 第二章 插值法 1. 根据( 2.2)定义的范德蒙行列式,令

数值分析课后习题答案

习 题 一 解 答 1.取3.14,3.15, 227,355113 作为π的近似值,求各自的绝对误差,相对误差和有效数字的位数。 分析:求绝对误差的方法是按定义直接计算。求相对误差的一般方法是先求出绝对误差再按定义式计算。注意,不应先求相对误差再求绝对误差。有效数字位数可以根据定义来求,即先由绝对误差确定近似数的绝对误差不超过那一位的半个单位,再确定有效数的末位是哪一位,进一步确定有效数字和有效数位。有了定理2后,可以根据定理2更规范地解答。根据定理2,首先要将数值转化为科学记数形式,然后解答。 解:(1)绝对误差: e(x)=π-3.14=3.14159265…-3.14=0.00159…≈0.0016。 相对误差: 3()0.0016 ()0.51103.14r e x e x x -==≈? 有效数字: 因为π=3.14159265…=0.314159265…×10,3.14=0.314×10,m=1。 而π-3.14=3.14159265…-3.14=0.00159… 所以│π-3.14│=0.00159…≤0.005=0.5×10-2=21311 101022 --?=? 所以,3.14作为π的近似值有3个有效数字。 (2)绝对误差: e(x)=π-3.15=3.14159265…-3.14=-0.008407…≈-0.0085。 相对误差: 2()0.0085 ()0.27103.15r e x e x x --==≈-? 有效数字: 因为π=3.14159265…=0.314159265…×10,3.15=0.315×10,m=1。 而π-3.15=3.14159265…-3.15=-0.008407… 所以│π-3.15│=0.008407……≤0.05=0.5×10-1 =11211101022 --?=? 所以,3.15作为π的近似值有2个有效数字。 (3)绝对误差: 22 () 3.14159265 3.1428571430.0012644930.00137e x π=-=-=-≈- 相对误差:

数值分析第四版习题及答案

第四版 数值分析习题 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2%,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式求下列各近似值的误差限: ********12412324(),(),()/,i x x x ii x x x iii x x ++其中**** 1234 ,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1%,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 1n n Y Y -=…) 计算到100Y .(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字. 8. 当N 充分大时,怎样求 2 11N dx x +∞ +? ? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±秒的误差,证明当t 增加时S 的绝对误 差增加,而相对误差却减小. 11. 序列 {}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字), 计算到 10y 时误差有多大?这个计算过程稳定吗? 12. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式 ln(ln(x x =- 计算,求对数时误差有多大?

数值分析第四版习题及答案

第四版 数值分析习题 第一章绪论 1.设x>0,x得相对误差为δ,求得误差、 2.设x得相对误差为2%,求得相对误差、 3.下列各数都就是经过四舍五入得到得近似数,即误差限不超过最后一位得半个单位,试指 出它们就是几位有效数字: 4.利用公式(3、3)求下列各近似值得误差限: 其中均为第3题所给得数、 5.计算球体积要使相对误差限为1%,问度量半径R时允许得相对误差限就是多少? 6.设按递推公式 ( n=1,2,…) 计算到、若取≈27、982(五位有效数字),试问计算将有多大误差? 7.求方程得两个根,使它至少具有四位有效数字(≈27、982)、 8.当N充分大时,怎样求? 9.正方形得边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝? 10.设假定g就是准确得,而对t得测量有±0、1秒得误差,证明当t增加时S得绝对误差增 加,而相对误差却减小、 11.序列满足递推关系(n=1,2,…),若(三位有效数字),计算到时误差有多大?这个计算过程 稳定吗? 12.计算,取,利用下列等式计算,哪一个得到得结果最好? 13.,求f(30)得值、若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式 计算,求对数时误差有多大? 14.试用消元法解方程组假定只用三位数计算,问结果就是否可靠? 15.已知三角形面积其中c为弧度,,且测量a ,b ,c得误差分别为证明面积得误差满足 第二章插值法 1.根据(2、2)定义得范德蒙行列式,令 证明就是n次多项式,它得根就是,且 、 2.当x= 1 , -1 , 2 时, f(x)= 0 , -3 , 4 ,求f(x)得二次插值多项式、 3. 4., 研究用线性插值求cos x 近似值时得总误差界、

第六章习题答案数值分析.docx

第六章习题解答 2 2、利用梯形公式和 Simpson 公式求积分 ln xdx 的近似值, 并估计两种方法计算值的最大 1 误差限。 解:①由梯形公式: T ( f ) b a [ f (a) f (b)] 2 1 [ln1 ln 2] ln 2 0.3466 2 2 2 最大误差限 R ( f ) (b a)3 f '' ( ) 1 1 1 0.0833 T 12 12 2 12 12 其中, (1,2) ②由梯形公式: b a 4 f ( b a f (b)] 1 4ln( 3 ln 2] 0.3858 S( f ) [ f (a) ) [ln1 ) 6 2 6 2 最大误差限 R S ( f ) (b a)5 f (4) ( ) 6 6 0.0021, 2880 2880 4 2880 其中, (1,2) 。 4、推导中点求积公式 f ( x)dx (b a) f ( a b ) (b a) 3 (a b) b a 2 24 证明: 构造一次函数 P ( x ),使 P a 2 b f a b , P ' ( a b ) f ' ( a b ), P '' ( x) 0 2 2 2 则,易求得 P( x) f ' ( a b )( x a b ) f ( a b ) 2 2 2 且 P(x)dx f ' ( a b )( x a b ) f ( a b ) dx b b a a 2 2 2 f ( a b )dx (b a) f ( a b ) ,令 P(x)dx I ( f ) b b a 2 2 a 现分析截断误差:令 r ( x) f ( x) P(x) f ( x) f ' ( a b )( x a b ) f ( a b ) 2 2 2 由 r ' ( x) f ' (x) f ' ( a b ) 易知 x a 2 b 为 r (x) 的二重零点, 2 a b )2 , 所以可令 r (x) ( x)( x 2

数值分析课后习题答案

第一章 题12 给定节点01x =-,11x =,23x =,34x =,试分别对下列函数导出拉格朗日插值余项: (1) (1) 3 ()432f x x x =-+ (2) (2) 4 3 ()2f x x x =- 解 (1)(4) ()0f x =, 由拉格朗日插值余项得(4)0123() ()()()()()()0 4!f f x p x x x x x x x x x ξ-=----=; (2)(4) ()4!f x = 由拉格朗日插值余项得 01234! ()()()()()() 4! f x p x x x x x x x x x -= ----(1)(1)(3)(4)x x x x =+---. 题15 证明:对于()f x 以0x ,1x 为节点的一次插值多项式()p x ,插值误差 012 10()()()max () 8x x x x x f x p x f x ≤≤-''-≤. 证 由拉格朗日插值余项得 01() ()()()()2!f f x p x x x x x ξ''-= --,其中01x x ξ≤≤, 01 0101max ()()()()()()()() 2!2!x x x f x f f x p x x x x x x x x x ξ≤≤''''-=--≤-- 01210()max () 8x x x x x f x ≤≤-''≤. 题22 采用下列方法构造满足条件(0)(0)0p p '==,(1)(1)1p p '==的插值多项式 ()p x : (1) (1) 用待定系数法; (2) (2) 利用承袭性,先考察插值条件(0)(0)0p p '==,(1)1p =的插值多项式 ()p x . 解 (1)有四个插值条件,故设230123()p x a a x a x a x =+++,2 123()23p x a a x a x '=++, 代入得方程组001231123010231 a a a a a a a a a =? ?+++=?? =? ?++=? 解之,得01230 021 a a a a =??=?? =??=-?

最新数值分析课程第五版课后习题答案(李庆扬等)1

第一章 绪论(12) 1、设0>x ,x 的相对误差为δ,求x ln 的误差。 [解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=* ****1)()(ln )(ln x x x x x , 相对误差为* * ** ln ln ) (ln )(ln x x x x r δ εε= = 。 2、设x 的相对误差为2%,求n x 的相对误差。 [解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而n x 的误差为n n x x n x n x x n x x x ** 1 *** %2%2) ()()()(ln * ?=='=-=εε, 相对误差为%2) () (ln )(ln *** n x x x n r == εε。 3、下列各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出它们是几位有效数字: 1021.1*1=x ,031.0*2=x ,6.385*3=x ,430.56*4=x ,0.17*5 ?=x 。 [解]1021.1*1 =x 有5位有效数字;0031.0* 2=x 有2位有效数字;6.385*3=x 有4位有效数字;430.56* 4 =x 有5位有效数字;0.17*5?=x 有2位有效数字。 4、利用公式(3.3)求下列各近似值的误差限,其中* 4*3*2*1,,,x x x x 均为第3题所给 的数。 (1)* 4*2*1x x x ++; [解]3 334* 4*2*11** *4*2*1*1005.1102 1 10211021)()()()()(----=?=?+?+?=++=? ??? ????=++∑x x x x x f x x x e n k k k εεεε; (2)* 3*2 *1x x x ;

(参考资料)数值分析课后答案1

1第一章 习题解答 1 设x >0,x 的相对误差限为δ,求 ln x 的误差。 解:设 x 的准确值为x *,则有 ( | x – x * | /|x *| ) ≤ δ 所以 e (ln x )=| ln x – ln x * | =| x – x * | ×| (ln x )’|x=ξ·≈ ( | x – x * | / | x *| ) ≤ δ 另解: e (ln x )=| ln x – ln x * | =| ln (x / x *) | = | ln (( x – x * + x *)/ x *) | = | ln (( x – x * )/ x * + 1) |≤( | x – x * | /|x *| ) ≤ δ 2 设 x = – 2.18 和 y = 2.1200 都是由准确值经四舍五入而得到的近似值。求绝对误差限ε( x ) 和 ε( y ) 。 解:| e (x ) | = |e (– 2.18)|≤ 0.005,| e (y ) | = |e ( 2.1200)|≤ 0.00005,所以 ε( x )=0.005, ε( y ) = 0.00005。 3 下近似值的绝对误差限都是 0.005,问各近似值有几位有效数字 x 1=1.38,x 2= –0.0312,x 3= 0.00086 解:根据有效数字定义,绝对误差限不超过末位数半个单位。由题设知,x 1,x 2, x 3有效数末位数均为小数点后第二位。故x 1具有三位有效数字,x 2具有一位有效数字,x 3具有零位有效数字。 4 已知近似数x 有两位有效数字,试求其相对误差限。 解:| e r (x ) | ≤ 5 × 10– 2 。 5 设 y 0 = 28,按递推公式 y n = y n-1 – 783/ 100 ( n = 1,2,…) 计算到y 100。若取≈78327.982 (五位有效数字),试问,计算 y 100 将有多大的误差? 解:由于初值 y 0 = 28 没有误差,误差是由≈78327.982所引起。记 x = 27.982,783?=x δ。则利用理论准确成立的递推式 y n = y n-1 – 783/ 100 和实际计算中递推式 Y n = Y n-1 – x / 100 (Y 0 = y 0) 两式相减,得 e ( Y n ) = Y n – y n = Y n-1 – y n-1 – ( x – 783)/ 100 所以,有 e ( Y n ) = e ( Y n-1) – δ / 100 利用上式求和 δ?=∑∑=?=100111001)()(n n n n Y e Y e 化简,得 e ( Y 100) = e ( Y 0) – δ = δ 所以,计算y 100 的误差界为 4100105001.05.0)(?×=×=≤δεY 6 求方程 x 2 – 56x + 1 = 0的两个根,问要使它们具有四位有效数字,D=ac b 42 ?至少要取几位有效数字? 如果利用韦达定理,D 又应该取几位有效数字? 解:在方程中,a = 1,b = – 56,c = 1,故D=4562?≈55.96427,取七位有效数字。

数值分析 第六章 习题

第六章 习 题 1. 计算下列矩阵的1A ,2A ,A ∞三种范数。 (1)1101A ???=????,(2)312020116A ????=??????? . 2. 用Jacobi 方法和Gauss-Seidel 迭代求解方程组 1231231 238322041133631236x x x x x x x x x ?+=??+?=??++=? 要求取(0)(0,0,0)T x =计算到(5)x ,并分别与精确解(3,2,1)T x =比较。 3. 用Gauss-Seidel 迭代求解 12312312 35163621122x x x x x x x x x ??=??++=???+=?? 以(0)(1,1,1)T x =?为初值,当(1)() 310k k x x +?∞?<时,迭代终止。 4. 已知方程组121122,2,x x b tx x b +=?? +=? (1)写出解方程组的Jacobi 迭代矩阵,并讨论迭代收敛条件。 (2)写出解方程组的Gauss-Seidel 迭代矩阵,并讨论迭代收敛条件. 5. 设有系数矩阵 122111221A ?????=?????? , 211111112B ?????=??????? , 证明:(1)对于系数矩阵A ,Jacobi 迭代收敛,而Gauss-Seidel 迭代不收敛. (2)对于矩阵B ,. 6. 讨论方程组 112233302021212x b x b x b ?????????????=??????????????????? 用Jacobi 迭代和Gauss-Seidel 迭代的收敛性;如果都收敛,比较哪种方法收敛更快.

第六章非线性方程的数值解法习题解答

第六章非线性方程的数值解法习题解答 填空题: 1. 求方程()x f x =根的牛顿迭代格式是__________________。 Ans:1()1()n n n n n x f x x x f x +-=- '- 2.求解方程 在(1, 2)内根的下列迭代法中, (1) (2) (3) (4) 收敛的迭代法是(A ). A .(1)和(2) B. (2)和(3) C. (3)和(4) D. (4)和(1) 3.若0)()(,故迭代发散。 以上三中以第二种迭代格式较好。 2、设方程()0f x =有根,且'0()m f x M <≤≤。试证明由迭代格式1()k k k x x f x λ+=- (0,1,2,)k =L 产生的迭代序列{}0k k x ∞ =对任意的初值0(,)x ∈-∞+∞,当2 0M λ<< 时,均收敛于方程的根。

数值分析第四版习题及答案

第四版 数值分析习题 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2%,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试 指出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式(3.3)求下列各近似值的误差限: ********12412324(),(),()/,i x x x ii x x x iii x x ++其中**** 1234 ,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1%,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 1n n Y Y -= ( n=1,2,…) 计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字27.982). 8. 当N 充分大时,怎样求 211N dx x +∞ +? ? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的 绝对误差增加,而相对误差却减小. 11. 序列{}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y ≈(三位有效数

字),计算到10y 时误差有多大?这个计算过程稳定吗? 12. 计算6 1)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大? 若改用另一等价公式 ln(ln(x x =- 计算,求对数时误差有多大? 14. 试用消元法解方程组 { 101012121010; 2. x x x x +=+=假定只用三位数计算,问结果是否可靠? 15. 已知三角形面积 1sin ,2s ab c = 其中c 为弧度, 02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ???证明面积的误差s ?满足 .s a b c s a b c ????≤++ 第二章 插值法 1. 根据( 2.2)定义的范德蒙行列式,令 20000112111 2 1 ()(,,,,)11 n n n n n n n n n x x x V x V x x x x x x x x x x ----== 证明()n V x 是n 次多项式,它的根是01,,n x x - ,且 101101()(,,,)()()n n n n V x V x x x x x x x ---=-- . 2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式. 3. 给出f (x )=ln x 的数值表用线性插值及二次插值计算ln 0.54 的近似值.

数值分析课后题答案

数值分析 第二章 2.当1,1,2x =-时,()0,3,4f x =-,求()f x 的二次插值多项式。 解: 0120121200102021101201220211,1,2, ()0,()3,()4; ()()1()(1)(2)()()2()()1()(1)(2)()()6()()1()(1)(1)()()3 x x x f x f x f x x x x x l x x x x x x x x x x x l x x x x x x x x x x x l x x x x x x x ==-===-=--= =-+-----= =------==-+-- 则二次拉格朗日插值多项式为 2 20()()k k k L x y l x ==∑ 0223()4() 1 4(1)(2)(1)(1)23 537623 l x l x x x x x x x =-+=---+-+=+- 6.设,0,1,,j x j n =L 为互异节点,求证: (1)0 ()n k k j j j x l x x =≡∑ (0,1,,);k n =L (2) 0()()0n k j j j x x l x =-≡∑ (0,1,,);k n =L 证明 (1) 令()k f x x = 若插值节点为,0,1,,j x j n =L ,则函数()f x 的n 次插值多项式为0()()n k n j j j L x x l x ==∑。 插值余项为(1)1()()()()()(1)! n n n n f R x f x L x x n ξω++=-=+ 又,k n ≤Q

(1)()0()0 n n f R x ξ+∴=∴= 0 ()n k k j j j x l x x =∴=∑ (0,1,,);k n =L 000(2)()() (())()()(())n k j j j n n j i k i k j j j i n n i k i i k j j i j x x l x C x x l x C x x l x =-==-==-=-=-∑∑∑∑∑ 0i n ≤≤Q 又 由上题结论可知 0()n k i j j j x l x x ==∑ 0()()0 n i k i i k i k C x x x x -=∴=-=-=∑原式 ∴得证。 7设[]2 (),f x C a b ∈且()()0,f a f b ==求证: 21max ()()max ().8 a x b a x b f x b a f x ≤≤≤≤''≤- 解:令01,x a x b ==,以此为插值节点,则线性插值多项式为 10101010()() ()x x x x L x f x f x x x x x --=+-- =()()x b x a f a f b a b x a --=+-- 1()()0 ()0 f a f b L x ==∴=Q 又 插值余项为1011()()()()()()2 R x f x L x f x x x x x ''=-=-- 011()()()()2 f x f x x x x x ''∴=--

相关主题
文本预览
相关文档 最新文档