《数值分析》第五章答案
习题5
1．导出如下3个求积公式，并给出截断误差的表达式。
(1) 左矩形公式：?-≈b
a a
b a f dx x f ))(()(
(2) 右矩形公式：))(()(a b b f dx x f b
a
-≈?
(3) 中矩形公式：?-+≈b
a
a b b
a f dx x f ))(2
(
)( 解：(1) )()(a f x f ≈， )()()()(a b a f dx a f dx x f b
a b
a -=≈??
????-=-=--b
a b
a b
a b
a dx a f x f dx a f dx x f a
b a f dx x f ))()(()()())(()(
),()(2
1
)()()()(2
ηηξf a b dx a x f dx a x f b
a b a
'-=-'=-'=??
),(,b a ∈ηξ
(2) )()(b f x f ≈，??-=≈b a
b
a
a b a f dx b f dx x f ))(()()(
????-=-≈--b a b a b a b
a dx
b f x f dx b f dx x f a b b f dx x f )]()([)()())(()(
)()(2
1)()()()(2ηηξf a b dx b x f dx b x f b
a b a
'--=-'=-'=??
，),(,b a ∈ηξ
(3) 法1 )2
(
)(b
a f x f +≈ ， ??-+=+≈b
a
b
a a
b b
a f dx
b a f dx x f ))(2
()2()(
?-+-b
a
a b b a f dx x f ))(2(
)(??+-=b a b a dx b a f dx x f )2
()( dx b a f x f b a ???????
+-=)2()( dx b a x f b a x b
a f
b a ???
?
??
?
+-''++-+'=2)2)((21)2)(2(ξ dx b a x f dx b a x b a f b
a b a 2)2
()(21)2()2(??+-''++-+'=η
3))((24
1
a b f -''=
η 法2 可以验证所给公式具有1次代数精度。作一次多项式 )(x H 满足 )2()2(
b a f b a H +=+，)2
()2(b
a f
b a H +'=+'，则有 )2
)(2()2(
)(b
a x
b a f b a f x H +-+'++= 2
)2)((!21)()(b a x f x H x f +-''=
-ξ， ),(b a ∈ξ
))(2
())(2()(a b b
a f a
b b a H dx x H b
a -+=-+=? 于是
dx x H dx x f a b b a f dx x f b a b
a b
a
???-=-+-)()()()2
(
)( []dx b a x f dx x H x f b
a
b
a 2
)2
(!2)()()(+-''=-=??
ξ 32))((24
1
)2(2)(a b f dx b a x f b a -''=+-''=?ηη 2．考察下列求积公式具有几次代数精度：
（1）
?'+
≈1
0)1(2
1
)0()(f f dx x f ； （2）
)3
1()31()(1
1f f dx x f +-
≈?-。 解： （1）当1)(=x f 时，左=1，右=1+0=1，左=右；
当x x f =)(时，左21=，右=2
1
210=+，左=右；
当2
)(x x f =时，左=3
1，右=1，左≠右，代数精度为1。
（2）当1)(=x f 时，左=2，右=2，左=右；
当x x f =)(时，左=0，右=03
1)31(=+-
，左=右； 当2
)(x x f =时，左32=
，右3
2
3131=+=，左=右； 当3
)(x x f =时，左0=，右0)3
1(
)31(33=+-
=，左=右；
当4
)(x x f =时，左52
=
，右9
2)31()31(22=+=，左≠右。代数精度为3。
3．确定下列公式中的待定参数，使其代数精度尽量高，并指出其代数精度的次数。
（1）
?-++-≈1
1)](3)(2)1([31
)(βαf f f dx x f ；
（2）
)]()([)()]()([2
)(2b f a f a b a b f a f a
b dx x f b
b
'-'-++-≈
?-； （3）
)1()0()1()(211
10f a f a f a dx x f ++-≈?-。
解：)1( 当1)(=x f 时，左2=，右2)321(3
1=++=，左=右；
当x x f =)(时，左0=，右)321(3
1βα++-=，
当2
)(x x f =时，左32=
，右)321(3
12
2βα++=； 要使所给求积公式至少具有2次代数精度当且仅当α、β满足
0)321(3
1
=++-βα 3
2)321(3122=++βα 132=+βα
13222=+βα
)21(3
1αβ-=
1)21(3
1
222=-+αα
314462
2
=+-+ααα 024102=--αα 01252=--αα
56156512,1±=±=
α ， 15
6251)61(521312,1 =??????±-=β 求积公式（1）：
??
?
???-+++-≈?-)156251(3)561(2)1(31)(1
1f f f dx x f （A ）
求积公式（2）：
??
?
???++-+-≈?-)156251(3)561(2)1(31)(1
1f f f dx x f （B ）
当3
)(x x f =时，（A ）的左端为1。
（A ） 的右端1)156251(3)561(213133≠??
????-++?+-=
（B ） 的右端1)156251(3)561(213133≠??
????++-?+-=
∴ （A ）和（B ）的代数精度均为2。
（2）
[]?'-'-++-≈
b
a
b f a f a b b f a f a
b dx x f )]()([)()()(2
)(2α 当1)(=x f 时，左a b -=，右a b a
b -=+-=)11(2
当x x f =)(时，左)(2122a b -=
，右)(2
1][222a b b a a b -=+-= 当2
)(x x f =时，左)(3
13
3a b -=，
右)2)(()(2
22
b a a b b a a b --++-=
αα ])(2)(2
1)[(2
22a b a b a b --+-=α
要使求积公式具有2次代数精度，当且仅当
)(3
1
])(2)(21)[(33222a b a b a b a b -=--+-α
)(3
1
)(2)(2122222a ab b a b a b ++=--+α )2(6
1)(222
2
a a
b b a b +-=
-α 121=α
)]()([)(12
1
)]()([2)(2b f a f a b b f x f a b dx x f b
a
'-'-++-≈?
当3
)(x x f =时，左),(4
144
3
a b dx x b
a
-=
=? 右]33[)(12
1][22223
3b a a b b a a b --++-=
])(222[4)(22222a b b ab a a b --+--=
2222222))((41
))((21a b a b b ab a a b ---+--=
)]2(222)[(4
1
222222a ab b b ab a a b +--+--= )(4
1
44a b -= 当4
)(x x f =时，左)(51554
a b dx x b
a
-=
=?，5b 的系数5
1=。 右)44()(12
1
][233244b a a b b a a b --++-=
，
其中5
b 的系数5
1
61)4(12121≠=-?+=。因而 代数精度为3。 5．设函数)(x f 由下表给出：
x 1.6 1.8 2.0 2..2 2.4 2.6 )(x f 4.953 6.050 7.389 9.025 11.023 13.464
x 2.8 3.0 3.2 3.4 3.6 3.8 )(x f 16.445 20.086 24.533 29.964 36.598 44.701
解： x 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4
)(x f 6.050 7.389 9.025 11.023 13.464 16.445 20.086 24.533 29.964
(1) 复化梯形公式
2.0=h ， ih x i +=8.1，8,,2,1,0 =i
])())()((21
[71
808∑=++=i i x f x f x f h T
023.11025.9389.7)964.29050.6(2
1[2.0++++?= ]533.24086.20445.16464.13++++ 9149.23=
（2）4.0=h
)]2.2()0.2(4)8.1([64.04f f f S ++=
)]6.2()4.2(4)2.2([64.0f f f +++ )]4.3()2.3(4)0.3([6
4
.0)]0.3()8.2(4)6.2([64.0f f f f f f ++++++
)]0.3()6.2()2.2([2)4.3()8.1({6
4
.0f f f f f ++?++=
)]}2.3()8.2()4.2()0.2([4f f f f +++?+
]086.20464.13025.9[2964.29050.6{6
4.0++?++=
]}533.24445.16023.11389.7[4+++?+
9149.23=
（3） Romberg 算法
8
4
42
22
1
111T S T C S T R C S T
8112.28)]4.3()8.1([28.14.31=+-=f f T
1768.25)]6.2(6.1[2
1
12=?+=f T T
2328.24))]0.3()2.2((8.0[21
24=+?+=f f T T
))]2.3()8.2()4.2()0.2((4.0[2
1
48f f f f T T +++?+=
9944.23=
9653.2331
34121=-=T T S
9181.2331
34242=-=T T S
9149.2331
34483=-=T T S
91495.2315
1
1516121=-=
S S C 91469.23151
1516242=-=S S C
91469.2363
1
6364121=-=
C C R 7．试用复化梯开公式计算曲线x x f tan )(=在区间[4
,
0π
]上这一段的弧长，取
3102
1-?=ε。
解： x x f tan )(=， x
x f 2
c o s 1)(=
'
dx x
dx x f S ??+
='+=40
4
40
2
cos 11)(1π
π
x
x g 4c o s 1
1)(+
= 2)0(=
g 5)4
(=π
g
43346.1)]4
()0([20
4
1=+-=π
π
g g T 54032
.137255583.28
c o s
1
1)8
(4
==+
=π
π
g 32161.1]54032.14
43346.1[21)]8(4[2112=?+=?+=πππg T T
32π 323π 32
5π 327π 0 16π 8π 16
3π
4π
44246
.1)16
(
=π
g 75848.1)16
3(=π
g ))]16
3()16((8[2124π
ππg g T T ++=
28931.1)]75848.144246.1(832161.1[21=++=
π
01077.03
1
24=-T T
??
?
?????? ??++++=)327()325()323()32(162148ππππg g g g T T T
??
????++++=
)94953.162881.148071.142109.1(1628931.121π 28084.1=
00282.03
1
48=-T T ??? ??++++=
)64
7()645()643()64(3221816πππππg g g g T T ??
?
???++++)6415()6413()6411(
)649(ππππg g g g 21=
[32
28084.1π
+(50746.145925.142986.141592.1+++ )]07792.284462.168747.158033.1++++ 27869.1=
38161072.03
1
-?=-T T [)128
7()1285()1283()128(64211632πππππg g g g T T ++ ??++=
)128
15()12813()12811()1289(π
πππg g g g ++++
)128
23()12821()12819()12817(
π
πππg g g g ++++ ??
?
???++++)12831()12829()12827()12825(
ππππg g g g [(43566.142501.141807.141464.164
27869.121
++++=
π 52307.149338.146936.145031.1+++
+ 72128.165675.160341.155935.1+++
+ )]15281
.201043.289446.179946.1++++
27816.1)27762.127869.1(2
1
=+=
331632102
1
10177.031--?
=82
4ln 1
dx x
计算4ln 时，若采用复化梯形公式，
点才能使其误差绝对值不超过
5102
1
-?。 解： 2=a ， 8=b ， 21)(=
x f ， 21)(x x f -='， 32)(x
x f =''
2)(12
)()(h f a
b f T dx x f b
a
n ξ''--
=-?， )8,2(∈ξ 要使 52102
1
)(1228-?≤''-ξf h 只要
5
3210212
221-?≤?h 5
2210
2-?≤h
522
1026-?≤??
? ??n
1010263
-?≤n
68.9481030010
1033
==?≥
n
取 949=n
答：取950个等距节点，则有
51
1
102
1
)(--?≤-?n T dx x f
方法2 []22222181121)()(121)()(h h b f a f f T f I n ??????-='-'≈- 52102
164141121)()(-?≤??? ??-≈
-h f T f I n
52102
1
6415121-?≤?h 6
2
10644-??≤h
3
1082-??≤h
310166-?≤n 375
1253108
3
3=?=?≥n 10．用Romberg 方法求?
821
dx x ，要求误差不超过5
102
1-?。从所取节点个数与上题结果比较中体会这2种方法的优缺点。
解： 将区间[2，8]作16等分，8
3
1628=- x
x f 1)(=
x 2， 2+8198
3
=
， 822， 825， 828， 831， 8
34
， 837，
)(x f 21， 198 228 258， 288， 318， 34
8， 378
x 840， 843
， 846， 849， 852， 855， 858， 861， 864
)(x f 408, 438, 468, 498, 528, 558, 588, 61
8, 648
[]875.16482126)8()2(2281=???
???+?=+-=
f f T []5375.14086875.121)5(62112=??
?????+=?+=
f T T ????????? ??+?+=
)852()8
28(32124f f T T
428090659.1528288
35375.121=??
??????? ??+?+?=
397126249.1)858()846()834()8
22
(5.12148=????????? ??+++?+=
f f f f T T )843()837()831()825()819
(75.031816f f f f f T T ++++??? ???+=
??
?
???+++)861()855()849(f f f 38903085.1=
875.11=T 425
.11=S 3893956.11=C 38643748.11=R 5375.12=T 39162087.12=S 3864837.12=C 3862979.12=R 428090659
.14=T 386804775.14=S 3863008.14=C 397126246.18=T 386332385.18=S 38903085
.116=T 3)4(3
13412121T T T T S -=-= 15)16(121S S C -= 63)64(121C C R -=
5712102
1
1047010549.5)(2551--?≤?-=-R R 38630
.1≈I 实际上
38630.1386294361.14ln ≈= 12．用3点Gauss-Legendre 公式求dx e I x
?-=10
。 解：
dx e x
?
-10
)1(2
1
t x +=
三点Gauss 公式
)5
3
(95)0(98)53(95)(1
1
g g g dt t g ++-≈?-
????
?
?
???
??
???????
??+++????
?
?
??-?=?253195)21(9825319521)(1
0f f f dx x f
??
???
????
???++?≈+
-----?25312125311095989521e e e dx e x
???
??
??
?++?=--2
6.026.02
1585181e e e 632120255.0=
21．根据下列x x f tan )(=的数值表：
x 1.20 1.24 1.28 1.32 1.36
)(x f 2.572 15 2.911 93 3.341 35 3.903 35 4.673 44
解： x x f tan )(= x x
x f 2
2t a n 1c o s 1)(+==
'
x x x x x x x f 3
2
t a n 2t a n 2)t a n 1(t a n 2)(t a n t a n
2)(+=+?='?='' h
h x f h x f h x D 2)
()(),(000--+=，
)(6
1),()(2
00ξf h h x D x f '''-=-'，),(00h x h x +-∈ξ
1330625
.1316.057215
.267344.420.136.1)20.1()36.1()08.0,28.1(=-=--=
f f D 39275
.1208
.091193
.290335.324.132.1)24.1()32.1()04.0,28.1(=-=--=
f f D )tan 1)(tan 62())(tan tan 62()(222x x x x x f ++='+=''' 16461982.1228.1tan 1)28.1(2=+='f
)
36.1(08.06
1
)(08.061)08.0,28.1()28.1(22f f D f '''??≤'''??=-'ξ )67344.41()67344.462(08.06
1
222+??+??=
241509202.3=
实际误差 96844268
.0)08.0,28.1()28.1(=-'D f )
90335.31()90335.362(04.06
1
)04.0,28.1()28.1(222+??+??≤-'D f 4044611.0=
实际误差 22813018.0)04.0,28.1()28.1(=-'D f
14597917.12)08.0,28.1(3
1
)04.0,28.1(34)04.0,28.1(~
=-=D D D
04.0=h
h
h x f h x f h h x f h x f 4)
2()2(312)()(3400~
0~
0--+-
--+= [])2()(8)(8)2(121
0000h x f h x f h x f h x f h
-+--+++=
018640653.0)04.0,28.1()28.1(~
=-'D f
)24.128.1()20.128.1(!
5)
()04.0,28.1()28.1()5(~
-?-=-'ξf D f
)36.128.1)(32.128.1(--? )1)(t a n 2t a n 15tan 15(8)(22
4)
5(+++=x x x x f
)267344.41567344.415(8!
51
)04.0,28.1()28.1(24~
+?+???≤
-'D f 2
2
2
04.008.0)167344.4(??+? 116713518.0=
数值分析课后题答案
数值分析 第二章 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)的值.若开平方用六位函数表,问求对数时误差有多大?若
数值分析习题集及答案
（适合课程《数值方法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、解：将)(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?当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
第六章习题答案数值分析.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
数值分析第四版习题及答案
第四版 数值分析习题 第一章绪论 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 近似值时得总误差界、
(完整版)数值分析第7章答案
第七章非线性方程求根 一、重点内容提要 （一）问题简介 求单变量函数方程 ()0f x = (7.1) 的根是指求*x （实数或复数），使得(*)0f x =.称*x 为方程(7.1)的根,也称*x 为 函数()f x 的零点.若()f x 可以分解为 ()(*)()m f x x x g x =- 其中m 为正整数,()g x 满足()0g x ≠,则*x 是方程(7.1)的根.当m=1时,称*x 为单根;当m>1时,称*x 为m 重根.若()g x 充分光滑,*x 是方程(7.1)的m 重根,则有 (1)() (*)'(*)...(*)0,(*)0m m f x f x f x f x -====≠ 若()f x 在[a,b]上连续且()()0f a f b <,则方程(7.1)在(a,b)内至少有一个实根,称[a,b]为方程(7.1)的有根区间.有根区间可通过函数作图法或逐次搜索法求得. (二)方程求根的几种常用方法 1.二分法 设()f x 在[a,b]上连续,()()0f a f b <,则()0f x =在(a,b)内有根*x .再设()0f x =在 (a,b)内仅有一个根.令00,a a b b ==,计算0001 ()2x a b =+和0()f x .若0()0f x =则*x x =,结束计算;若00()()0f a f x >,则令10,1a x b b ==,得新的有根区间11[,]a b ;若 00()()0 f a f x <,则令 10,10 a a b x ==,得新的有根区间 11[,]a b .0011[,][,]a b a b ?,11001()2b a b a -=-.再令1111 ()2x a b =+计算1()f x ,同上法得 出新的有根区间22[,] a b ,如此反复进行,可得一有根区间套 1100...[,][,]...[,] n n n n a b a b a b --????
最新数值分析课程第五版课后习题答案(李庆扬等)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. 计算下列矩阵的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 迭代的收敛性；如果都收敛，比较哪种方法收敛更快.
东南大学_数值分析_第七章_偏微分方程数值解法
第七章 偏微分方程数值解法 ——Crank-Nicolson 格式 ****（学号） *****（姓名） 上机题目要求见教材P346,10题。 一、算法原理 本文研究下列定解问题(抛物型方程) 22(,) (0,0)(,0)() (0) (0,)(), (1,)() (0)u u a f x t x l t T t x u x x x l u t t u t t t T ?αβ???-=<<≤≤???? =≤≤??==<≤?? (1) 的有限差分法，其中a 为正常数，,,,f ?αβ为已知函数，且满足边界条件和初始条件。关于式(1)的求解，采用离散化方法，剖分网格，构造差分格式。其中，网格剖分是将区域{}0,0D x l t T =≤≤≤≤用两簇平行直线 (0) (0)i k x x ih i M t t k k N τ==≤≤?? ==≤≤? 分割成矩形网格，其中,l T h M N τ==分别为空间步长和时间步长。将式(1)中的偏导数使用不同的差商代替，将得到不同的差分格式，如古典显格式、古典隐格式、Crank-Nicolson 格式等。其中，Crank-Nicolson 格式具有更高的收敛阶数，应用更广泛，故本文采用Crank-Nicolson 格式求解抛物型方程。 Crank-Nicolson 格式推导：在节点(,)2 i k x t τ +处考虑式(1)，有 22(,)(,)(,)222 i k i k i k u u x t a x t f x t t x τττ??+-+=+?? (2) 对偏导数 (,)2 i k u x t t τ ?+?用中心差分展开 []2311+13 1(,)(,)(,)(,) ()224k k i k i k i k i i k i k u u x t u x t u x t x t t t t ττηητ++??+=--<
数值分析课后习题答案
第一章 题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 =??=?? =??=-?
数值分析第七章非线性方程求根习题答案
第七章非线性方程求根 （一）问题简介 求单变量函数方程 ()0f x = (7.1) 的根是指求*x （实数或复数），使得(*)0f x =.称*x 为方程(7.1)的根,也称*x 为函数() f x 的零点.若()f x 可以分解为 ()(*)()m f x x x g x =- 其中m 为正整数,()g x 满足()0g x ≠,则*x 是方程(7.1)的根.当m=1时,称*x 为单根;当m>1时,称*x 为m 重根.若()g x 充分光滑,*x 是方程(7.1)的m 重根,则有 (1)() (*)'(*)...(*)0,(*)0m m f x f x f x f x -====≠ 若()f x 在[a,b]上连续且()()0f a f b <,则方程(7.1)在(a,b)内至少有一个实根,称[a,b]为方程(7.1)的有根区间.有根区间可通过函数作图法或逐次搜索法求得. (二)方程求根的几种常用方法 1.二分法 设()f x 在[a,b]上连续,()()0f a f b <,则()0f x =在(a,b)内有根*x .再设()0f x =在(a,b)内 仅有一个根.令00,a a b b ==,计算0001()2x a b =+和 0()f x .若0()0f x =则*x x =,结束计算;若 00()()0 f a f x >,则令 10,1a x b b ==,得新的有根区间 11[,] a b ;若 00()()0 f a f x <,则令 10,10a a b x ==,得新的有根区间11[,]a b .0011[,][,]a b a b ?,11001 () 2b a b a -=-.再令1111 ()2x a b =+计算1()f x ,同上法得出新的有根区间22[,]a b ,如此反复进行,可得一有根区 间套 1100...[,][,]...[,] n n n n a b a b a b --???? 且110011 *,0,1,2,...,()...() 22n n n n n n a x b n b a b a b a --<<=-=-==-. 故 1 lim()0,lim lim ()* 2n n n n n n n n b a x a b x →∞→∞→∞-==+=
数值分析课后题答案
数值分析 第二章 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 ''∴=--
第六章非线性方程的数值解法习题解答
第六章非线性方程的数值解法习题解答 填空题： 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 λ<< 时，均收敛于方程的根。
数值计算方法第七章习题 2013
计算方法 第七章 习题 复习与思考题 1．设f ∈C [a , b ]，写出三种常用范数2 1 f f 及∞ f 。 2．f , g ∈C [a , b ]，它们的内积是什么？如何判断函数族{? 0, ? 1, …, ? n }∈C [a , b ]在[a ,b ]上线性无关？ 3．什么是函数f ∈C [a , b ]在区[a , b ]上的n 次最佳一致逼近多项式？ 4．什么是f 在[a , b ] 上的n 次最佳平方逼近多项式？什么是数据{}m i f 0的最小二乘曲 线拟合？ 5．什么是[ a , b ]上带权ρ (x )的正交多项式？什么是[ -1, 1 ]上的勒让德多项式？它有什 么重要性质？ 6．什么是切比雪夫多项式？它有什么重要性质？ 7．用切比雪夫多项式零点做插值得到的插值多项式与拉格朗日插值有何不同？ 8．什么是最小二乘拟合的法方程？用多项式做拟合曲线时，当次数n 较大时为什么不直接求解法方程？ 9．哪种类型函数用三角插值比用多项式插值或分段多项式插值更合适？ 10．判断下列命题是否正确？ （1）任何f (x ) ∈C [a , b ]都能找到n 次多项式P n (x ) ∈ H n ，使| f (x ) - P n (x ) | ≤ ε ( ε 为任给的误差限)。 （2）n n H x P ∈)(* 是f (x )在[ a , b ]上的最佳一致逼近多项式，则)()(lim * x f x P n n =∞ →对 ],[b a x ∈?成立。 （3）f (x ) ∈C [a , b ]在[a , b ]上的最佳平方逼近多项式P n (x ) ∈ H n 则)()(lim x f x P n n =∞ →。 （4）)(P ~ x n 是首项系数为1的勒让德多项式，Q n (x ) ∈ H n 是任一首项系数为1的多项式，则 ? ? --1 1 21 1 2d )(d )](P ~ [x x Q x x n n 。 （5）)(T ~ x n 是[-1 , 1]上首项系数为1的切比雪夫多项式。Q n (x ) ∈ H n 是任一首项系数为1的多项式，则 .)(max )(~ max 1 11 1x Q x T n x n x ≤≤-≤≤-≤ （6）当数据量很大时用最小二乘拟合比用插值好。
数值分析复习题答案
数值分析复习题 一、填空 Chapter1 绪论 近似数x*=0.4231关于真值x=0.4229有 3 位有效数字. 用1000.1近似真值1000时，其有效数字有 4 位， 已知准确值x*与其有t 位有效数字的近似值12 10.10(0)s n x a a a a =?≠的绝对误差为 1 x*-x 102s t -≤ ?。 设 2.40315x * =是真值 2.40194x =的近似值，则x * 有 3 位有效数字。 设一近似数x*=2.5231具有5位有效数字，则其相对误差限是44 11 1010224--?=?? ，其绝对误差限是4 1 102-?。 当x 很大时，为防止损失有效数字，应该使 = 。 Chapter2 插值方法 设642 ()3651f x x x x =+-+，则[3,2,1,0,1,2,3]f ---= 3 。 若 42 f(x)=2x +x -3, 则f[1,2,3,4,5,6]= 0 。 对 32f(x)=x +3x -x+5,差商f[0,1,2,3,4]= 0 。 设 643()35f x x x x =-+-，则差商[0,1,2,3,4,5,6]f = 1 。 已知y=f(x)的均差 021[,,]5f x x x =, 402[,,]9f x x x =, f[x4, x3, x2]=14， f[x0, x3, x2]=8 ,.那么 均差f[x4, x2, x0]= 9 。（交换不变性） 设有数据112 032 x y -则其 2 次 Larange 插值多项式为 32 (1)(2)(1)(1)23x x x x -+-++-，2次拟合多项式为 （最佳平方逼近可求）。？？？ 以n + 1个 整 数 点k ( k =0,1,2,…,n) 为 节 点 的 Lagrange 插 值 基 函 数 为 ()k l x ( k =0,1,2,…,n)，则 n k k=0 kl (x)= ∑ x 。？？（注： k y k =，则有拉格朗日插值公式：