数值分析习题汇总
第一章 引论(习题)
2.证明:x 的相对误差约等于x 的相对误差的1/2.
证明 记 x x f =
)( ,则
)
()(*
**
x x x x x x
x x f E r +-=
-=
)(21**x E x x x x x x
r ≈-?+=
. □
3.设实数a 的t 位β进制浮点机器数表示为)(a fl . 试证明 t
b a b a fl -≤
+*=*12
1||),1/()()(βδδ, 其中的记号*表示+、-、?、/ 中一种运算. 证明: 令:
)
()
()(b a fl b a fl b a **-*=
δ
可估计: 1|)(|-≥*c b a fl β (c 为b a *阶码),
故:
121||--≤
c t c ββδt -=12
1β 于是: )1()()(δ+*=*b a b a fl . □
4.改变下列表达式使计算结果比较精确:
(1) ;1||,
11211<<+--+x x
x
x 对
(2)
;1,11>>-
-+
x x
x x
x 对
(3) 1||,0,cos 1<<≠-x x x
x
对.
解 (1) )21()1(22
x x x ++.
(2)
)
11(2x x x x x
-++.
(3) x
x
x x x x x cos 1sin )cos 1(sin cos 12+≈+=-. □
6.设937.0=a 关于精确数x 有3位有效数字,估计a 的相对误差. 对于x x f -=1)(,估计)(a f 对于)(x f 的误差和相对误差. 解 a 的相对误差:由于 31021|)(|-?≤
-=a x x E . x
a
x x E r -=)(, 221018
1
10921)(--?=?≤
x E r . (1Th ) )(a f 对于)(x f 的误差和相对误差. |11||)(|a x f E ---==
()25
.0210
11321??≤
-+---a
x x a =3
10-
33
104110
|)(|--?=-≤a f E r . □
9.序列}{n y 满足递推关系:1101.100-+-=n n n y y y . 取01.0,110
==y y 及
01.0,
101150=+=-y y ,试分别计算5y ,从而说明该递推公式对于计算是不稳
定的.
解 递推关系: 1101.100-+-=n n n y y y (1) 取初值 10=y , 01.01=y 计算 可得: 110
01.1002
2-?=-y 10001.1-=410-=
6
310-=y , 8
410
-=y , 10
510-=y , …
(2) 取初值 5
0101-+=y , 2
110
-=y ,
记:
n n n y y -=ε,
序列 {}n ε ,满足递推关系,且 5
010--=ε , 01=ε
1101.100-+-=n n n εεε, 于是: 5210-=ε,
531001.100-?=ε, 55241010)01.100(---?=ε,
5
5351002.20010)01.100(--?-?=ε,
可见随着 n ε 的主项 52
10)01.100(--?n 的增长,说明该递推关系式是不稳定的.
第二章 多项式插值 (习 题)
1. 利用Lagrange 插值公式求下列各离散函数的插值多项式(结果要简化):
(1)
(2)
解(2):
方法一. 由 Lagrange 插值公式
)
30)(20)(10()
3)(2)(1()(0x x x x x x x x x x x x x l ------=
,
)()()()()(332211003x l f x l f x l f x l f x L ?+?+?+?=
)1)((31)2)()(1()1)(()(21
23210---=-----=x x x x x x x l ,
))(1(2)
1)()(1()(2122
1
211--=--+=
x x x x x x l ,
x x x x x x l )1()()1()1!()(238
2
1
21232--=-??-+=
, )()1(12)()1()(212
1213-+=??-+=x x x x x x x l . 可得: )21()(2
3-=x x x L
方法二. 令
)()21()(3B Ax x x x L +-= 由 23)1(3-
=-L , 2
1
)1(3=L , 定A ,B (称之为待定系数法) □
2. 设)(,),(),(10x l x l x l n Λ是以n x x x ,,,10Λ为节点的n 次多项式插值问题的基函数. (1)证明
.,,2,1,0,
)(0n k x x l x n
i k i k i Λ==∑=
(2)证明
Λ+----+--+
=)
)(()
)((1)(2010101000x x x x x x x x x x x x x l
)
())(()
())((02010110n n x x x x x x x x x x x x ------+
-ΛΛ.
证明(1) 由于 j i j i x l ,)(δ= 故: =
)(x L n ∑=n
i i k i x l x
)( , 当 j x x = 时有: k j j n x x L =)( ,
n j ,,1,0Λ= )(x L n 也即为 k
x 的插值多项式,由唯一性,有:
∑==n
i k i k i x x l x
)( , n k ,,1,0Λ=
证明(2):
利用Newton 插值多项式
)(],[)()(0100x x x x f x f x N n -+=)()(],,[100---++n n x x x x x x f ΛΛΛ )()
()()
()()(00101x l x x x x x x x x x f n n =----=
ΛΛ
差商表:
f(x) 一阶 二阶 … n 阶差商
0x 1
1x 0
1
01
x x -
M
M 0 )
()(1
1020x x x x --
M
M M O
n x 0 0
)
()(1
010n x x x x --Λ
代入)(*式有:
)
()()()
()(1)(020*******n n n x x x x x x x x x x x x x x x N -----++--+
=-ΛΛΛ.
)(0x l 为n 次代数多项式,由插值多项式的唯一性:
有 )()(0x N x l n ≡. □
4. 设a b b a C x f -<<∈ε0],,[)(3.考虑以b a a ,,ε+为节点的Lagrange 插值公式当
0→ε时的极限.证明成立公式
)()()(x R x p x f +=. 其中
)()
()
2)(()(2
a f a
b a b x x b x p --+-=
),()
()()())((2
2
b f a b a x a f a b x b a x --+'---+ )()()(6
1
)(2ξf b x a x x R '''--=
,),(b a ∈ξ 并计算)(),(),(a p b p a p '.
解 作)(x f 以b a a ,,ε+为节点的Lagrange 插值多项式,有: )()()(22x R x L x f +=, 其中:
)()()
()()()()()()()(2εεεεε+-+--+-----=a f b a b x a x a f b a b x a x x L
)()()()
()(b f a b a b a x a x εε------+
, )()()(!
3)
()(2b x a x a x f x R ----'''=εζ , b a <<ζ 令: 0→ε 有
)()(6)
()()(22b x a x f x R x R --'''=
→ζ, 又: )()()()([
)()(2a f a b a
x a f a b a x x b x L εεεεε----+----= )]()()
()()(a f a b a x a f a b a x -------+
εεεε )()
()()
()(b f a b a b a x a x εε------+
)()
()
2()(2
a f a
b a b x x b --+-→)()()()(a f a b a x x b '---+
)()()
()(2
2
x P b f a b a x =--+ 故当 0→ε 时,成立公式: )()()(x R x P x f +=. □
5. 给出13)(4+-=x x x f 的数值表
解:因为34)(3
'
-=x x f ,2
'
'12)(x x f = )(x f 为凹函数.又从数值表可见:
1) 当]5.0,1.0[∈x 时,)(x f 单调下降.有反函数)(1
y f x -=
)(y f
的Newton 插值多项式:
)
17440.0)(10810.0)(40160.0)(70010.0(01225.0)
10810.0)(40160.0)(70010.0(01531.0)
40160.0)(70010.0(0096436.0)
70010.0(33500.01.0)(4+---+------+--=y y y y y y y y y y y N
.337.0)0(4*≈=N x □
7、若 1)(3
7
++=x x x f ,问:
?]2,,2,2[710=Λf ; ?]2,,2,2[810=Λf .
解 1)(3
7
++=x x x f .有:
=]2,,2,2[7
1
Λf !
7)()7(ξf =1, !8)(]2,,2,2[)8(8
10ηf f =Λ0=. □
9.证明下列关系的正确性:
(1) i i i i i i f g g f g f ?+?=??+1)( (2)
1/)()/(+?-?=?i i i i i i i i g g g f f g g f
(3) )()(/!)1()/1(nh x h x x h n x n
n n ++-=?Λ 证明:
(1) =?-?=??++i i i i i i g f g f g f 11)(i i i i i i i i g f g f g f g f ?-?+?-?++++1111
i i i i f g g f ?+?=+1.
(3) n x
n n )1()1(-=?!)()(nh x h x x h
n
++Λ
此题可利用数学归纳法:
设 k n = 成立,证明 1+=k n 成立.又 1=n 时是成立的. □
10. 利用差分性质证明:
2333]2/)1([21)(+=+++=n n n n g Λ.
[提示:考虑差分)()1()(n g n g n g -+=?,并利用差分和函数值可互相表示] 证明: 记: 2
]2/)1([)(+=n n n f ,3
3
3
21)(n n g +++=Λ 有: 3)1()()1()(+=-+=?n n f n f n f 故: ∑-=?=
10
)()(n k k f n g ∑-=-+=1
)]()1([n k k f k f
2]2/)1([)0()(+=-=n n f n f . □
13、求次数4≤的多项式)(x P ,满足1)1()0()1(===-P p P ,2)1()1(='=-'P P 解 作重节点差商的Newton 插值公式 )1(]1,1[)1()(+--+-=x f f x P
2
2
)1(]1,0,1,1[)1(]0,1,1[+--++--+x x f x f )1()1(]1,1,0,1,1[2
-+--+x x x f 重节点差商表:
i x i f 一阶 二阶 三阶 四阶
10-=x 1
10-=x 1 2
01=x 1 0 -2
12=x 1 0 0 1
12=x 1 2 2 1 0
得 2
2
)1()1(2)1(21)(+++-++=x x x x x P 13
+-=x x . □
17、 设 ]1,0[)(2
C x f ∈ ,并且 0)0(=f ,1)1()2
1(==f f
求证
?
≥''1
212))((dx x f
证: 取 ,00=x 2
11=x , 12=x , 2
1=h 00=f , 11=f , 12=f 记: )(i i x s M ''= , 2,1,0=i
有 h x x M h x x M x S 01
101
)(-+-=''x M x M 102)2
1
(2+-= )2
1(2)1(2)(212
-+-=''x M x M x S 又三弯矩方程为:( 2],,[210-=x x x f )
244210-=++M M M , )24(4
1
201M M M ++-
=. 分段积分:
?
?+''=''?
1
2
12
2
1
)]([)]([dx x s dx x s ?
''1
22
2
1)]([dx x s ?
+-+=2
1
201)]21
([4dx x M x M ?-+-12
1221)]21()1([4dx x M x M
??-+-+-+-=1211
2
1221201)]21()1([4)]1()21([4dx
x M x M dx x M x M
由于 ?
=-12
1224
1)21(dx x ,
?
=-12
1224
1
)1(dx x ,
?
=--12
148
1)1()21(dx x x ,于是:
?++++=''?
1
02
2212110202]2[6
1))((M M M M M M M dx x S 又: )24(4
1
201M M M ++-=
记 =
),(20M M I ?
?
''1
2))((dx x S =)()24(4
1
[6120202
22
0M M M M M M +++-
+ ])24(8
1
220M M +++
由
00
=??M I ,
02=??M I
. 得: ??
?=+-=-0
70
72020M M M M 即当: 020==M M 时, ),(20M M I 达最小
故:
?
=??≥''?
1
2212)24(8
1
61))((dx x S ,由最小模原理: ?
≥''1
212)]([dx x f . □
20.
解 利用三弯矩方法 )(i i x s M ''= , 2,1,0=i 10=x , 22=x , 32=x
???
??-=+=++=+54
23646221
21010M M M M M M M
解得: 70-=M , 201=M , 372-=M
]2,1[∈x 72
431729)(231-+-=
x x x x s
]3,2[∈x
1052
29367219)(232+-+-
=x x x x s . □
第四章 数值积分方法与数值微分 (习 题)
1.直接验证梯形公式(1.2)与中矩形公式(1.3)具有1次代数精度,而辛甫生公式
(1.4)则具有3次代数精度. 解
梯形公式:
?+-≈
b
a
b f a f a
b dx x f )]()([2
)(. 矩形公式: ?+-≈b a b
a f a
b dx x f )2
(
)()(. 以上两求积公式以 ,1)(=x f x 代入公式两边,结果相等,而以2
)(x x f =
代入公式两边,其结果不相等.故梯形公式的代数精度等于1. Simpson 公式:
?
+++-≈
b a
b f b
a f a f a
b dx x f )]()2
(4)([6)(. 容易验证:以2
,,1)(x x x f =分别代入Simpson 公式两边,结果相等。
以3
)(x x f =代入 左边=
)(444
13a b dx x b
a
-=?
右边=[
]
32
322322332
3
3
3
36246b ab b a a a b b b a a a b +++-=???
?????+???
??++- =
).(4
144
a b - Simpson 公式两边,结果相等。而以4
x 代入Simpson 公式两边,其结果 不相等。故Simpson 求积公式的代数精度为3. □
3.对于?
=h dx x f I 30
)(的数值积分公式?
=
h h dx x p I 30
)(,其中)(x p 为对)(x f 在
h h x 2,,0=进行插值的2次多项式.证明:)()0(8
354h O f h I I h +'''?=-.
证明: )(x P 为)(x f 于h h x 2,,0=进行插值的二次多项式,则: )()()(x R x P x f += 其中: )2()(!
3)
()(h x h x x f x R --'''=ξ. 求积分公式误差
?
?-=h
h
dx x P dx x f f E 3030
)()()(
?
--'''=h
dx h x h x x f 30
)2()(!
3)
(ξ ?
--'''=
h
dx h x h x x f 30
)2()(!
3)
0( ?
--+
h
dx h x h x x f 30
)4()2()(!
3)
(ξη, 其中: h 20<<ξ, )(0)0(8
3
)(54h h f f E +?'''=. □
4.证明 中矩形公式的Peano 核误差公式为:
?
?
''=
-h h dt t f t k h hf dx x f 0
)()()2/()(,
其中 ??
?≤≤-≤
≤=,
,
2/)(,
0,
2/)(2
2
2
2h t t h t t t k h h
并由此导出误差形式
?
∈''+?=h h f h h f h dx x f 0
3
],0[),
(24
)2/()(ξξ.
解 已知中矩形公式对于一次多项式精确成立,由Taylor 展开: ?''-+'+=+h
dt t f t x x f f x f 0
)()()0()0()(
)()(1x R x P +=. ?
?-=
h h
f h dx x f f E 0
)2
()()(
又: )()()(1R E R P E f E =+=
])()([)(0?
''-=+h
dt t f t x E R E
??''-=+
h
h dx dt t f t x 0
)()(?''-?-+
h
dt t f t h
h 0
)()2
(
??''---=+
+h h
dt t f t h
h dx t x 00
)(])2
()([
?++-?--=h
t h h dx t x t K 0)2()()(+---=
)2
(2)(2t h h t h
?
?
???≤≤-≤≤=h t h
t h h
t t 2
!2)(2
0,
2
2
2
. □
)(24
)()()()()(3
00
ξξf h dt t K f dt t f t K f E h
h
''=''=''=??
5. 求系数321,A A A 和,使求积公式
?
-+-+-≈11
321),3/1()3/1()1()(f A f A f A dx x f
对于次数2≤的一切多项式都是精确成立的. 解:求积公式
?
-+-+-≈1
1
321)3
1()31()1()(f A f A f A dx x f 是一个插值型求积公式,令 2
,,1)(x x x f =得: ,2321=++A A A
,031
31321=+--A A A 329
1
91321=++A A A ,
解得: 2
11=A , 02=A ,233=A
12. 确定参数a 使求积公式的代数精度尽可能地高
)].()0([)]()0([2
)(20
h f f ah h f f h
dx x f h '-'++≈
?
(*) 解 令: n
x x f =)( , 2≥n 得:
11
12111+++-=+n n n anh h h n , an n -=+2
111, )1(21+-=n n n a
(*)
公式
对1)(=x f 、x 精确成立.
当 2=n 时, 12
1=a , 3=n 时, 12
1
=a ,4=n 时,40
3
=a ,
故:当取 12
1=a 时,(*)具有3次代精确度. □
13 假定求积公式
?
-≈11
002)()(x f A dx x f x
对于1,x 精确成立,试求00,A x 解: 由
?-=1102A dx x , ?-=1
1003
x A dx x
可得: 3231130
==-x A ,?-=112
01dx x A x 0= 故: ?-≈112
)0(3
2)(f dx x f x . □
14. 建立Gauss 型求积公式:)()()(221110
x f A x f A dx x
x f +≈?
.
解: 令:
)()()(211x x x x x f --=, 和 )()()(212x x x x x x f --= 代入 得: 7
6
21=+x x , 35321=?x x
5672731?-=
x , 5
6
72732+=x 653111+
=A , 6
5
3112-=A 。 □
16. 求数值微分公式的余项.
h h x f h x f x f x f 2/))2()(4)(3()(0000+-++-≈'.
解:于 0x ,h x 20+ ,h x +0三点作)(x f 的Lagrange 插值多项式: )(2)
2()()(02
002x f h h x x h x x x L ----=
)()2()(0
200h x f h h x x x x +----+
)2(2)
()(02
00h x f h h x x x x +---+. )(2322)(02
02
x f h
h
x x x L --=' )()
222(02
0h x f h
h x x +---+
)2(2)22(020h x f h h x x +--+. 令 0x x = ,得:
)()(02
0x L x f '≈'h h x f h x f x f 2/))2()(4)(3(000+-++-= 余项:因为
)2)()((!
3)
()()()(000)3(2h x x h x x x x f x L x f x R -----=
-=ξ 有
.3
)()()()(2
)3(02
00h f x L x f x R ξ='-'='
第五章 线性代数方程组的解法 (习 题)
3.设A 为n 阶按行严格对角占优矩阵,经Gauss 消去法一步后A 变为如下形式 ?
????
?=)2(221211)2(0A A a A
试证)
2(22A 是1-n 阶按行严格对角占优矩阵.
证明: 已知: ∑
=≠>n
j i j ij ii a a 1
, ),,2,1(n i ΛΛ=, 证明:
>)
2(ii
a
∑≠=n
j j ij
a 1
2
)
2( , ),,,3,2(n i ΛΛΛ= 右=
∑
≠=?-
n
i j j j i ij a a a a 211
11∑
≠=+≤n
i
j j ij a 2∑≠=?n
i
j j j i a a a 2
1111
∑≠=-=
n
i
j j ij
a
11
i a ?+111a a i ???
? ??-∑=n
j i j a a 2
11 11
11a a a a i i ii +
-<()i a a
111
-
11
111111a a a a a a a a i i ii i i ii ?-≤?-
=)
2(ii a = . (对列也可以证明) □ 4.设A 为实对称非奇异矩阵,且各阶顺序主子式
,,,1,0n k k Λ=≠?
证明:A 可以分解为T LDL A =,其中L 为具有正对角元的下三角阵,D 为对角阵,其对角元1||=ii d .
证明:LDU A = (已知) L —单位下三角阵,U —单位上三角阵,D=diag(nn d d ,,11Λ)
由于: A A T
= , 故 U L T
= ,写:
?????????
???=nn d d d D O 22
11=????
????????212111nn d d O . ???????
?????
????nn nn d d d d 0011
11O ??????
???????2121
11nn d d O ,
21
21
~
E D E =.
T T
L E D E L L D L A )~(212
1==T T L D L LE D LE ~~~)(~)(212
1==
D ~的对角元绝对值等于1.L ~
的对角元即为21E 的对角元>0. □
6. 假定已知n n R A ?∈的三角分解:LU A =,试设计一个算法来计算1-A 的),(j i 元素.
解: 记: ),,,(211
n a a a A Λ=-,其中: j a 为1-A 的第j 列元素
由于: I AA =-1 , 故:j j e a A = , T
j e )0,0,1,,0(ΛΛ=
T
nj j j j a a a a ),,,(21Λ= 矩阵1-A 的),(j i 元素即为ij a
因为 U L A ?= 有 j j e a U L =?)(, 记: (**)
.
b a U j = 有 (*)
.
j e Lb =
解系数矩阵为三角阵的方程组(*)求得b ,然后解系数矩阵为上三角阵的方程
组(**),即可求得j a ,因而得到ij a 的结果。
7.试证对n 维向量x 有
∞∞
≤≤x n x x
1.
证明: ∑==n
k k x x
1
1
∞=≤x n x n k k
max
又 k k
x x max 1≥, 故: ∞∞
≤≤x n x x
1. □
8.设A 为n 阶实矩阵,试证 .
12
F F
A A A n
≤
≤
证明: nxn ij a A )(= ∑
∑===n
i n
j ij
F
a
a
1
1
22)(A A tr T =
又: n T
A A tr λλλ+++=Λ21)(
其中:i λ为A A T 的特征值,由于A A T
为半正定矩阵,有0≥i λ
故: i i
T
A A tr λmax )(≥, 有:()
22
1
max A A
i
i
F
=≥λ.
又: (
)
i i
T
n A A tr λmax )(≤, 故:
21A A n
F
≤. □
9.设‖·‖是向量范数,A 为n n ?实矩阵,x 是n 维向量,证明Ax 是x 的连续函数. 证明 Ay Ax - Ay Ax -≤y x A -?≤ 有: Ay Ax → 当 y x → 时. □
10.设B A ,为n 阶非奇异矩阵,‖·‖ 表示矩阵的任一种从属范数,试证
(1)A A /11≥-,
(2)B A B A B A -??≤-----1111.
证明:
(1) A
A 1
1≥-
事实上: I A A =-1
故: A A I ?≤-1, 而 1=I , 有: A
A 1
1≥-
(2) B A B A B A -??≤-----1111 事实上: )(111
1
-----=-AB I A B
A 11)(---=
B A B A
故: 1111----?-?≤-B A B A B A . □
(习题) 证明:矩阵的∞-范数(行和范数): ∑=∞
=n
j ij i
a A
1
||max 。
证明:因为:
Ax x
Ax A x x 1
max max
=≠==
对任何1=∞
x
,
}||||{max }{|max 1
1
11∑∑==≤≤≤≤∞
?≤=n
j n
j j ij n
i j ij n
i x a x a Ax
∞=????
?
??
=∑x a n j ij i
1||max 从而
∑=∞
≤n
j ij i
a Ax
1
||max ,
∑=∞
=≤∞
n
j ij i
x
a Ax
1
1||max max
现设∑∑===n
j n
j kj ij
i
a a
1
1
||||max
,令
)(kj j a sign y =
显然,T
n y y y y ),,,(21Λ=满足1=∞
y
,并且
∑=≥∞∞=∞
j
j j i i
x
y a Ay Ax ,1max max
∑∑====n
j n
j ij i
kj
a a
11
||max ||,
从而∑=∞
=n j ij i
a A
1
||max
11. 设‖·‖是由向量范数‖·‖诱导的矩阵范数, 证明:若n n R A ?∈非奇异,则
.min 1
1
1Ax A x =--=
证明: 事实上: x
x A x A A
x x 111
1
max
max -≠-=-==θ
,
令: y x A =-1
, Ay x =, 于是: Ay
y A
y max
1
θ
≠-=y
Ay y θ≠=
min 1Ax
x 1
min 1
==
故: Ax A x 1
1
1min =--=. , □
12.设LU A =是n n R A ?∈的三角分解,其中1||≤ij l .并设T i T i u a ,分别表示A 和U 的第i 行, 验证
∑-=-=1
1i j T j
ij T i T i u l a u , 并证明 ∞-∞
≤A U
n 12.
证: 由 U L A ?= 故:
=
????????
?
?????????T n T
i T a a a M M 1??????
???????????????????????????
?T n T T
nn n n ii i i u u u l l l l l l l l l M ΛΛΛΛΛΛ212121222111
0 有: T
i ii T i T i T i u l u l u l a +++=Λ2211, 1=ii l , n i ,,2,1Λ=
故: ∑-=?-=1
1
i j T j ij T
i
T i
u l a u
.
由 T
T
a u 11= 故: ∞≤=A a u T T
1
1
1
1
归纳法,设 ∞-≤A u k T k 11
α
由: =+1
1T
k u
1
1
,11
∑=++?-k
j T j j k T k u
l a
+≤+1
1T
k a
∑=k
j T j
u 1
1
∞-∞++++≤A A k )222(11Λο∞=A k 2,
故: ∞-≤A u n T n 11
2. 而: T j j
u U
m ax =∞
有: ∞-∞
≤A U
n 1
2
. □
13.设A 非奇异,i λ是方阵A A T 的特征值,证明
i i
i i
T A A cond A cond λλmin /max )()(22
2==.
证明:首先证明A A T 为正定矩阵,事实上,由于A 非奇异,若Ax Y =,
0≠x 则 0≠Y ,考虑当0≠x 时有
0),(),(>=Ax Ax x Ax A T
即A A T 为正定矩阵,设它的特征值为:021>≥≥≥n λλλΛ 2
1
2
2)
()(-?=A A A A A A Cond T
T
T
122
)(])[(λρρ===A A A A A A T T T
))((]))[(()
(1212
1
---==A A A A A A T T T
ρρ=
n
λ1
故: i
i
n T
A A Cond λλλλmin max )(12== 又: 2
12
2
2
2)(-?=A A A Cond , 由于A 为对称矩阵,有 2
2
2A A A
T
= , 21
2
1)(--=A A A
T 故: i
i T
A A Cond A Cond λλmin max )()(222==. □
14.设
??
?
???=??
????=11,
98.099.099.01b A . 已知方程组b Ax =的精确解为T x )100,100(*-=. (1)计算条件数∞
-∞?1A A ;
(2)取T T x x )5.99,5.100(~,)0,1(~21-==,
分别计算它的残余向量.本题的结果说明了什么问题? 解
??
?
???=98.099
.099.01A , 99.1=∞A 又: ???
?
??????--=-----4444
11011099.01099.01098.0A
41
1099.1?=∞
-A
故: (1) ∞
-∞
∞?=1)(A A
A cond .3960110)99.1(42=?=
( 2)取T
x )0,1(~
1=,得: ()T
x A b r 01.0,0~11=-=
取T
x )5.99,5.100(~
2-=,得:
???
?
??--=???? ???-?-?--=-=985.0995.0)5.9998.05.10099.0(1)5.9999.05.100(1~22x A b r
1~x 、2~x 均可作为方程组的近似解, 但1~x 和2~x 有很大的差别,说明方程组是病态的.
15.求矩阵Q 的p
Q
,∞=,2,1p ,以及)(Q cond ∞,其中
????
?
????
???------=1111111111111111
Q . 解: ,41==∞
Q
Q
)(2
Q Q Q
T ?=ρ
????
?????
???=4000040000400004Q Q T ,故:4)(=Q Q T ρ..22=Q
又:I Q Q T 4=, 故:T
Q Q 411=-
)(Q cond ∞=4141
=?=?∞
-∞
Q Q
. □
16.若存在正定矩阵P ,使 PH H P B T -=为对称正定阵, 试证迭代法
b Hx x k k +=+)()1(, Λ,2,1,0=k
收敛.
证明: 已知 0),(>x Px , 0),(>x Bx , 证明
1 设λ,u 分别为H 的任意特征值和相对应的特征向量,有: ),)((),(u u PH H P u Bu T -=),(),(Hu PHu u Pu -= ),(),(2 u Pu u Pu λ-= ),()1(2 u Pu λ-=. 由B P ,的正定性,有:012>-λ , 1<λ. □ 17.设有方程组b Ax =,其中 ?? ??? ?????-=220122101A , ???? ? ?????-=3/23/12/1b . 已知它有解T x )0,3/1,2/1(-=.如果右端有小扰动2/106-∞=b δ,试估计由此引起的解的相对误 差. 解: ∞ ∞∞ -∞ ∞ ∞ ?≤-b b A A x x x δ1 ~******* ????? ?????-=220122101A ?? ?? ??????----+-=-11223121111A 5=∞A , 2 91 =∞ -A 故: 3 22102 95~6 -∞ ∞ ? ?≤-x x x 8101356-?=. □ 18.设有迭代格式 g Bx x k k +=-)1()(,Λ,2,1=k 其中 ????? ???? ?-=05.02/15.005 .02/15.00B , ???? ??????--=5.015.0g . 试证该迭代格式收敛.并取T x )0,0,0() 0(=,计算)4(x . 解: 可先验证: 0)(=B ρ. 迭代格式收敛,由于0)(=B ρ,故: ?? ?? ? ?????=01010~J B 即 T J T B 1-=, 03 =B 记: ) (*k k x x -=ε 于是: 0εεk k B = 03=ε 即 *) 3(x x = , )4(x 即为精确解. □ 19.给定方程组 ?? ?? ? ?????=????????????????? ???-11112211 1221321x x x . 特别声明:考试时需带计 算器作辅助计算 1.2015x *=是经四舍五入得到的近似值,则其相对误差* r e ≤-31 104 ?. 2. 01(),(), ,()n l x l x l x 是以01,, ,n x x x 为节点的拉格朗日插值基函数,则 3.设(0)1(1)3(2)4(3)2f =,f =,f =,f =,[0123]f =,,,1 3 - . 4. 利用Simpson 公式求?2 1 2dx x = 7.3 5. 设求积公式1 0()d (),(1)n k k k f x x A f x n ≈≥∑?=是Gauss 型求积公式,则3 n k k k A x == ∑1 .4 6. 数值微分公式(2)(2) ()i i i f x h f x h f x h +≈ --'的截断误差为 2().O h 7. 设1101A ?? = ??? ,则A 的谱半径()A ρ= 1 ,A 的条件数1cond ()A = 4. 8. 用牛顿下山法求解方程3 03 x x -=根的迭代公式是 2 13 3(1),3n n n n x x x x x λ+-=-- 下山条件是 1()().n n f x f x +< 9.对任意初始向量(0)x 及任意向量f ,线性方程组的迭代公式(1)()(0,1,2,)k k k +=+=x Bx f ,迭代序列()k x 收敛于方程组的精确解x *的充分必要条件是()1.ρ 数值分析试题 一、 填空题(2 0×2′) 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值,则x 有 2 位 有效数字。 2. 若f (x )=x 7-x 3+1,则f [20,21,22,23,24,25,26,27]= 1 , f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设,‖A ‖∞=___5 ____,‖X ‖∞=__ 3_____, ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ,则使用该迭代 函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时,若所求节点靠近首节点,应该选用等距节点下牛顿差商 公式的 前插公式 ,若所求节点靠近尾节点,应该选用等距节点下牛顿差商公式的 后插公式 ;如果要估计结果的舍入误差,应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是:=∑=n i i x a 0)( 1 ;所以当 系数a i (x )满足 a i (x )>1 ,计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%,至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ,线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收 敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。 12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的,其中的残差 r i = (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ,(i =0,1,…,n )。 13. 在非线性方程f (x )=0使用各种切线法迭代求解时,若在迭代区间存在唯一解,且f (x ) 第一章典型例题 例3 ln2=0.…,精确到10-3的近似值是多少 解 精确到10-3=,即绝对误差限是=, 故至少要保留小数点后三位才可以。ln2 第二章典型例题 例1 用顺序消去法解线性方程组 ??? ??1 -=4+2+4=+2+31 -=4++2321 321321x x x x x x x x x 解 顺序消元 ?? ?? ??????---???→???????????---????→???????????--=-?+-?+-?+1717005.555.00141 25.025.105.555.001412142141231412]b A [)3()2/1()2/3(231312r r r r r r M 于是有同解方程组 ?? ? ??-==--=++17175.555.0142332321x x x x x x 回代得解 x 3=-1, x 2=1,x 1=1,原线性方程组的解为X =(1,1,-1)T 例2 取初始向量X (0)=(0,0,0)T ,用雅可比迭代法求解线性方程组 ??? ??5 =+2+23=++1=2-2+321 321321x x x x x x x x x 解 建立迭代格式 ???????+--=+--=++-=+++5223122) (2)(1)1(3 ) (3)(1)1(2 ) (3)(2)1(1k k k k k k k k k x x x x x x x x x (k =1,2,3,…) 第1次迭代,k =0 X (0)=0,得到X (1)=(1,3,5)T 第2次迭代,k =1 ???????-=+?-?-=-=+--==+?+?-=3 532123 351515232)2(3) 2(2)2(1x x x X (2)=(5,-3,-3)T 第3次迭代,k =2 ???????=+-?-?-==+---==+-?+-?-=1 5)3(2521 3)3(511)3(2)3(2)2(3) 3(2)3(1x x x X (3)=(1,1,1)T 第4次迭代,k =3 ???????=+?-?-==+--==+?+?-=1 512121 311111212)2(3) 2(2)2(1x x x X (4)=(1,1,1)T 例4 证明例2的线性方程组,雅可比迭代法收敛,而高斯-赛德尔迭代法发散。 证明 例2中线性方程组的系数矩阵为 A =?? ?? ? ?????-122111221 于是 D =?? ?? ??????100010001 D -1=D ??????????=022001000L ~ ????? ?????-=000100220U ~ 雅可比迭代矩阵为 第四章 数值积分与数值微分 1.确定下列求积公式中的特定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度: 101210121 12120 (1)()()(0)(); (2)()()(0)(); (3)()[(1)2()3()]/3; (4)()[(0)()]/2[(0)()]; h h h h h f x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-?? ?? 解: 求解求积公式的代数精度时,应根据代数精度的定义,即求积公式对于次数不超过m 的多项式均能准确地成立,但对于m+1次多项式就不准确成立,进行验证性求解。 (1)若101(1) ()()(0)()h h f x dx A f h A f A f h --≈-++? 令()1f x =,则 1012h A A A -=++ 令()f x x =,则 110A h Ah -=-+ 令2 ()f x x =,则 3 221123 h h A h A -=+ 从而解得 01 1431313A h A h A h -?=?? ?=?? ?=?? 令3 ()f x x =,则 3()0h h h h f x dx x dx --==? ? 101()(0)()0A f h A f A f h --++= 故 101()()(0)()h h f x dx A f h A f A f h --=-++? 成立。 令4 ()f x x =,则 数值分析典型习题 特别声明:考试时需带计 算器作辅助计算 1.2015x *=是经四舍五入得到的近似值,则其相对误差* r e ≤ -31 104 ?. 2. 01(),(),,()n l x l x l x L 是以01,,,n x x x L 为节点的拉格朗日插值基函数,则 3.设(0)1(1)3(2)4(3)2f =,f =,f =,f =,[0123]f =,,,1 3 - . 4. 利用Simpson 公式求?2 1 2dx x = 7.3 5. 设求积公式1 0()d (),(1)n k k k f x x A f x n ≈≥∑?=是Gauss 型求积公式,则3 n k k k A x == ∑1 .4 6. 数值微分公式(2)(2) ()i i i f x h f x h f x h +≈ --'的截断误差为 2().O h 7. 设1101A ?? = ??? ,则A 的谱半径()A ρ= 1 ,A 的条件数1cond ()A = 4. 8. 用牛顿下山法求解方程3 03 x x -=根的迭代公式是 2 13 3(1),3n n n n x x x x x λ+-=-- 下山条件是 1()().n n f x f x +< 9.对任意初始向量(0)x 及任意向量f ,线性方程组的迭代公式(1)()(0,1,2,)k k k +=+=L x Bx f ,迭代序列()k x 收敛于方程组的精确解x *的充分必要条件是()1.ρ 华南理工大学研究生课程考试 《数值分析》试卷C (2015年1月9日) 1. 考前请将密封线内各项信息填写清楚; 所有答案请按要求填写在本试卷上; 课程代码:S0003004; 4. 考试形式:闭卷; 5. 考生类别:硕士研究生; 本试卷共八大题,满分100分,考试时间为150分钟。 一、(12分)解答下列问题: 1)设近似值0x >,x 的相对误差为δ,试证明ln x 的绝对误差近似为δ。 2)利用秦九韶算法求多项式 542()681p x x x x x =-+-+ 在3x =时的值(须写出计算形式),并统计乘法次数。 (12分)解答下列问题: 1)设()235f x x =+,求[]0,1,2f 和[]0,1,2,3f 。 2)利用插值方法推导出恒等式: 33220,0[]j j i i x j i x i j =≠=-=-∑∏ 。 (1)设{}∞ =0)(k k x q 是区间[]1,0上带权1=ρ而最高次项系数为1的正交多项式族,其中1)(0=x q ,求1()q x 和2()q x 。 (2)求形如2y a bx =+的经验公式,使它与下列数据拟合: 四、(14分)对积分()10I f x dx = ?,试 (1)构造一个以012113,,424 x x x ===为节点的插值型求积公式; (2)指出所构造公式的代数精度; (3)用所得数值求积公式计算积分1 203x dx ?的精确值; (4)指出所得公式与一般的Newton-Cotes 型公式在形式上的重要区别。 (1)设?? ????=4321A ,计算1A 、()Cond A ∞和()A ρ。 (2)用列主元Gauss 消去法解方程组: 12312315410030.112x x x ????????????=????????????-?????? 六、(13分)对2阶线性方程组 11112212112222 a x a x b a x a x b +=??+=? (11220a a ≠ ) (1)证明求解此方程组的Jacobi 迭代与Gauss-Seidel 迭代同时收敛或同时发散; (2)当同时收敛时,试比较它们的收敛速度。 《数值分析》 重点考察内容及各章作业答案 学院: 学号: 姓名: 重点考察内容 基本概念(收敛阶,收敛条件,收敛区域等), 简单欧拉法。 第一章基础 掌握:误差的种类,截断误差,舍入误差的来源,有效数字的判断。 了解:误差限,算法及要注意的问题。 第二章插值 掌握:Hermite插值,牛顿插值,差商计算,插值误差估计。 了解:Lagrange插值 第三章数据拟合 掌握:给出几个点求线性拟合曲线。 了解:最小二乘原理 第四章数值积分微分 掌握:梯形公式,Simpson公式,代数精度,Gauss积分,带权Gauss积分公式推导,复化梯形公式推导及算法。 了解:数值微分,积分余项 第五章直接法 掌握:LU分解求线性方程组,运算量 了解:Gauss消去法,LDL,追赶法 第六章迭代法 掌握:Jacobi,Gauss-Seidel迭代格式构造,敛散性分析,向量、矩阵的范数、谱半径 了解:SOR迭代 第七章Nolinear迭代法 掌握:牛顿迭代格式构造,简单迭代法构造、敛散性分析,收敛阶。 了解:二分法,弦截法 第八章ODE解法 掌握:Euler公式构造、收敛阶。 了解:梯形Euler公式、收敛阶,改进Euler公式 题目类型:填空,计算,证明综合题 第一章 误差 1. 科学计算中的误差来源有4个,分别是________,________,________,________。 2. 用Taylor 展开近似计算函数000()()'()()f x f x f x x x ≈+-,这里产生是什么误差? 3. 0.7499作 3 4 的近似值,是______位有效数字,65.380是舍入得到的近似值,有____几位有效数字,相对误差限为_______. 0.0032581是四舍五入得到的近似值,有_______位有效数字. 4. 改变下列表达式,使计算结果比较精确: (1)11,||1121x x x x --++ (2 ||1x (3) 1cos ,0,|| 1.x x x x -≠ (4)sin sin ,αβαβ-≈ 5. 采用下列各式计算61)时,哪个计算效果最好?并说明理由。 (1) (2 )99-3 )6 (3-(4 6. 已知近似数*x 有4位有效数字,求其相对误差限。 上机实验题: 1、利用Taylor 展开公式计算0! k x k x e k ∞ ==∑,编一段小程序,上机用单精度计算x e 的函数 值. 分别取x =1,5,10,20,-1,-5,-10,-15,-20,观察所得结果是否合理,如不合理请分析原因并给出解决方法. 2、已知定积分1 ,0,1,2,,206 n n x I dx n x ==+? ,有如下的递推关系 111 110 0(6)61666 n n n n n x x x x I dx dx I x x n ---+-===++-? ? 可建立两种等价的计算公式 (1) 1016,0.154n n I I I n -= -=取; (2) 12011),0.6n n I nI I n -=-=(取 第四章 数值积分与数值微分 1.确定下列求积公式中的特定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度: 101210121 12120 (1)()()(0)(); (2)()()(0)(); (3)()[(1)2()3()]/3; (4)()[(0)()]/2[(0)()]; h h h h h f x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-?? ?? 解: 求解求积公式的代数精度时,应根据代数精度的定义,即求积公式对于次数不超过m 的多项式均能准确地成立,但对于m+1次多项式就不准确成立,进行验证性求解。 (1)若101(1) ()()(0)()h h f x dx A f h A f A f h --≈-++? 令()1f x =,则 1012h A A A -=++ 令()f x x =,则 110A h Ah -=-+ 令2 ()f x x =,则 3 221123 h h A h A -=+ 从而解得 011431313A h A h A h -?=?? ? =?? ?=?? 令3 ()f x x =,则 3()0h h h h f x dx x dx --==? ? 101()(0)()0A f h A f A f h --++= 令4()f x x =,则 455 1012()5 2 ()(0)()3 h h h h f x dx x dx h A f h A f A f h h ---== -++=? ? 故此时, 101()()(0)()h h f x dx A f h A f A f h --≠-++? 故 101()()(0)()h h f x dx A f h A f A f h --≈-++? 具有3次代数精度。 (2)若 21012()()(0)()h h f x dx A f h A f A f h --≈-++? 令()1f x =,则 1014h A A A -=++ 令()f x x =,则 110A h Ah -=-+ 令2 ()f x x =,则 3 2211163 h h A h A -=+ 从而解得 1143 8383A h A h A h -?=-?? ? =?? ?=?? 令3 ()f x x =,则 22322()0h h h h f x dx x dx --==? ? 101()(0)()0A f h A f A f h --++= 第一章典型例题 例3…,精确到10-3的近似值是多少? 解 精确到10-3=,即绝对误差限是?=, 故至少要保留小数点后三位才 可以。ln2? 第二章典型例题 例1 用顺序消去法解线性方程组 解 顺序消元 于是有同解方程组 回代得解 x 3=-1, x 2=1,x 1=1,原线性方程组的解为X =(1,1,-1)T 例2 取初始向量X (0)=(0,0,0)T ,用雅可比迭代法求解线性方程组 解 建立迭代格式 ??? ????+--=+--=++-=+++5223122)(2)(1)1(3) (3)(1)1(2 )(3)(2)1(1k k k k k k k k k x x x x x x x x x (k =1,2,3,…) 第1次迭代,k =0 X (0)=0,得到X (1)=(1,3,5)T 第2次迭代,k =1 X (2)=(5,-3,-3)T 第3次迭代,k =2 X (3)=(1,1,1)T 第4次迭代,k =3 X (4)=(1,1,1)T 例4 证明例2的线性方程组,雅可比迭代法收敛,而高斯-赛德尔迭 代法发散。 证明 例2中线性方程组的系数矩阵为 A =?? ?? ? ?????-122111221 于是 D =?? ?? ??????100010001 D -1 =D ?? ?? ? ?????=022001000L ~ ?? ?? ? ?????-=000100220U ~ 雅可比迭代矩阵为 B 0=?? ?? ? ?????--=??????????-??????????-=+--022101220022101220100010001)U ~L ~(D 1 得到矩阵B 0的特征根03,2,1=λ,根据迭代基本定理4,雅可比迭代法收敛。 高斯-赛德尔迭代矩阵为 G =-U ~ )L ~D (1-+ =-?? ?? ??????----=??????????-??????????---=??????????-??????????-2003202200001002201200110010001002201220110011 解得特征根为?1=0,?2,3=2。由迭代基本定理4知,高斯-赛德尔迭代发散。 例5 填空选择题: 1. 用高斯列主元消去法解线性方程组 作第1次消元后的第2,3个方程分别为 。 数值分析典型例题 例1 对下列各数写出具有5位有效数字的近似值。236.478, 0.00234711, 9.000024, 9.0000343 10?. 解:按照定义,以上各数具有5位有效数字的近似值分别为:236.478, 0.0023471, 9.0000, 9.0000310?。 注意: *x =9.000024的5位有效数字是9.0000而不是9,因为9 是1位有效数字。 例2 指出下列各数具有几位有效数字。2.0004, -0.00200, -9000, 9310?, 23 10-?。 解:按照定义,以上各数的有效数字位数分别为5, 3, 4,1,1 例3 已测得某物体行程* s 的近似值s=800m ,所需时间* s 的近似值为t=35s ,若已知m s s s t t 5.0||,05.0||**≤-≤-,试求平均速度v 的绝对误差和相对误差限。 解:因为t s v /=,所以)()(1)()()(2t e t s s e t t e t v s e s v v e -=??+??≈ 从 而 05.00469.035 800 5.0351|)(||||)(|1|)(|22≤≈+?≤+≤t e t s s e t v e 同样v v e v e r )()(≈)()()()(t e s e t e v t t v s e v s s v r r r -=??+??= 所以00205.035 05 .08005.0|)(||)(||)(|≈+≤+≤t e s e v e r r r 因此绝对误差限和相对误差限分别为0.05和0.00205。 例4试建立积分20,,1,05 =+=n dx x x I n n 的递推关系,并研究它的误差 传递。 解:151 --= n n I n I ……………………………………………..…...(1) 5ln 6ln 0-=I ,计算出0I 后可通过(1)依次递推计算出1I ,…,20I 。 但是计算0I 时有误差0e ,由此计算出的1I ,…,20I 也有误差,由(1)可 知近似值之间的递推关系为 151 --= n n I n I ……………………………………………….…..(2) (1)-(2)可得 01)5(5e e e n n n -=-=-,由0I 计算n I 时误差被放大了n 5倍。所以(1)不稳 定。 (1) 可以改写为 n I I n n 51 511+ -=- ……………………………………… (3) 如果能先求出20I ,则依次可以求出19I ,…,0I ,计算20I 时有误差,这样根据(3)计算19I ,…,0I 就有误差,误差传播为 n n n e e ?? ? ??-=-511 ,误差依次减少。 例5 用二分法求解方程012)(23=+--=x x x x f 在区间[0,1]内的1个实根,要求有3为有效数字。 解:因为0)1()0( 数值分析习题集 (适合课程《数值方法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 -=…) 计算到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)的值.若开平方用六位函数表,问求对数时误差有多大? 数值分析第8章 数值积分与数值微分 8.1 填空题 (1)n+1个点的插值型数值积分公式∫f(x)dx b a ≈∑A j n j=0f(x j )的代数精度至少是 n ,最高不超过 2n+1 。【注:第1空,见定理8.1】 (2)梯形公式有 1 次代数精度,Simpson 公司有 3 次代数精度。【注:分别见定理8.1,8.3】 (3)求积公式∫f(x)dx h 0≈h 2[f (0)+f (h )]+ah 2[f ′(0)?f ′(h)]中的参数a= 1/12 时,才能保证该求积公式的代数精度达到最高,最高代数精度为 3 。 解:令f(x)=1,x,x 2带入有, { h 2[1+1]+ah 2[0?0]=h h 2[0+h ]+ah 2[1?1]=12 (h 2)h 2[0+h 2]+ah 2[0?2h ]=13 (h 3) //注:x 的导数=1 解之得,a=1/12,此时求积公式至少具有2次代数精度。 ∴ 积分公式为:∫f(x)dx h 0≈h 2[f (0)+f (h )]+h 2 12[f ′(0)?f ′(h)] 令 f(x)= x 3带入求积公式有:h 2 [0 +h 3]+ h 212 [0?3h 2]=14 (h 4),与f(x)= x 4的定积分计算值1 4 (h 4)相等, 所以,此求积公式至少具有3次代数精度。 令f(x)= x 4带入求积公式有,h 2[0+h 4]+h 2 12[0?4h 3]=1 6(h 5),与f(x)= x 5的定积分计算值1 5(h 5)不相等,所以,此求积公式的最高代数精度为3次代数精度。 8.2 确定下列求积公式的求积系数和求积节点,使其代数精度尽量高,并指出其最高代数精度。 解题思路:按照P149 中8.3式进行求解,根据求积公式中未知量n 的数量决定代入多少f(x),当积分公式代入求积节点x n 的计算结果与定积分的计算结果一致,继续代入求积节点X n+1,,若计算结果与对应的定积分计算结果不一致时,求积公式拥有最高n 次的代数精度。 (1)∫f(x)dx 2h 0≈A 0f (0)+A 1f (h )+A 2f(2h) 解:令f(x)=1,x,x 2代入有,【注:本例中需求解A 0、A 1、A 2共3个未知量,故需3个相异求积节点f(x)】 {A 0+A 1+A 2=2h A 1h +A 22h =1 2(2h )2A 1h 2+A 2(2h )2=1 3(2h )3 求解得A 0=13h ,A 1=43h ,A 2=1 3h , ∴求积公式为:∫f(x)dx 2h 0≈13hf (0)+43hf (h )+1 3 hf(2h) ∵该求积公式对3个相异节点1,x,x 2均有余项E (f )=0, //注:参见P149定理8.1 ∴该求积公式至少具有2次代数精度。 令f(x)= x 3,代入求积公式有:4 3hh 3+1 3h (2h )3=4h 4 ∵函数f(x) = x 3的定积分结果为:∫x 3dx 2h 0=1 4(2h )4=4h 4 ,与求积公式计算值相等, ∴该求积公式具有3次代数精度。 第四版 数值分析习题 第一章 绪 论 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.1 求下列各数的具有四位有效数字的近似值, 并指出其绝对误差限和相对误差限 )1.0ln(,121,101 1,1014321== = = x x x x 1.2 下列各数都是对准确值进行四舍五入得到的近似值, 指出它们的绝对误差限、相对误差限和有效数字的位 数。 3 * 5* 4* 3* 2* 1100.5,5000,50.31,3015.0,0315.0?=====x x x x x 1.3 为了使 3 1的近似值的相对误差不超过0.1%, 问应取几位有效数字? 1.4 怎样计算下列各题才能使得结果比较精确? (1) x x sin )sin(-+ε,其中ε充分小 (2) ? ++1 2 1N N x dx ,其中N 是充分大的正数 (3) x x sin cos 1-,其中x 充分小 (4) o 1cos 1- (5) 1001.0-e (6) )11010ln(84-- 1.5 求方程01562=+-x x 的两个根, 使至少具有四位有效数字。 习题二 2.1 证明方程043 =-+x x 在区间[1,2]内有且仅有一个根。如果用二分法求它具有五位有效数字的根,试问需对 分多少次?(不必求根) 2.2 用二分法求方程0134 =+-x x 在[0.3, 0.4]内的一个根, 精度要求2 10 2 1-?= ε。 2.3 找出下列方程的有根区间,选择适当的初始点用二分法求方程的根,精度要求2 10 -=ε。 (1) 02 =--x x ; (2) 06cos 2 =-++-x e x x ; (3) 01tan =--x x ; (4) 0sin 2=--x e x 。 2.4 考虑方程032 =-x e x ,将其改写为3 x e x ± =,取00=x ,用两种迭代公式迭代,分别收敛到1.0和-0.5附 近的两个根(取精度要求3 10-=ε)。 ()()()()()()()()()收敛较慢 代入上式得:将解: 收敛速度次并分析该迭代公式的迭代的根求方程 取试用迭代公式∴≠<<*'*+++-='∴+*+*=*∴=+?+?? ? ??===++= =∴++= ==-++=++=++014.01022220||10 2202613381013202132020 132010212010220. 2.0 20102110220 4.1222 222212012123021x x x x x x x x x x x x x x x x x x x x x x x x k k k k k k k ?????? )))()()()[]()()[])49998.0cos 215.0cos 2 1,022,00cos 2 102 12,0210,2,0.cos 2 10sin 2 11,cos 2 113cos 2 12; 1.0cos 2 12.4120101==== ==->-=<-=-=>+='-===-+x x x x x x x f f x x x f x x f x x x f x x x x k k 则 取上有一个根在所以上在为单调递增函数故则令解: 位有效数字求出这些根,精确到用迭代公式分析该方程有几个根给定方程ππππ 500 .0105.0102.0||3412≈*?=---x x x 所以方程的根 41444444466666.6663.4k k S S S S s +=+=++++++=+故迭代公式为可知: 由解: 动点迭代公式:导出下列连根公式的不Λ ΛΛΛ()()()()()()()()()()()()))()))() )()?得到的是什么迭代公式步迭代时选取第?得到的是什么迭代公式选取使收敛速度快; 选取的单根附近收敛; ,使迭代在选取值写出迭代公式是参数其中的迭代公式 给定方程不收敛 解: 都不收敛于迭代则对任何初值都有数证明:如果对于任何实为一实数设k k k k k k k k k k k k k k k k k k k x f k x f x f x f x x x x x x x G x G x x x x x G x G x x x G x G x G x Gx x x x x x G x x x G x '='==-=∴*-≥≥-≥-=*-∴-≥-∴≥--='*=*≠≥'**=*+++++++++1514302. 1. , 1.45.41 ,,1.,4.40101111111100λλλλλλΛΛ 数值分析试卷1 一、填空题(每空2分,共30分) 1. 近似数231.0=*x 关于真值229.0=x 有____________位有效数字; 2. 设)(x f 可微,求方程)(x f x =根的牛顿迭代格式是_______________________________________________; 3. 对1)(3++=x x x f ,差商=]3,2,1,0[f _________________; =]4,3,2,1,0[f ________; 4. 已知??? ? ??-='-=1223,)3,2(A x ,则=∞||||Ax ________________,=)(1A Cond ______________________ ; 5. 求解线性方程组?????=+=+045 11532121x x x x 的高斯—赛德尔迭代格式为_______________________________________;该迭代格式迭代矩阵的谱半径=)(G ρ_______________; 二、(12分)(1)设LU A =,其中L 为下三角阵,U 为单位上三角阵。已知 ?????? ? ??------=2100121001210012A ,求L ,U 。 (2)设A 为66?矩阵,将A 进行三角分解:LU A =,L 为单位下三角阵,U 为上三角阵,试写出L 中的元素65l 和U 中的元素56u 的计算公式。 三、给定数据表如下 x 0.20.40.60.81 1.2f(x)212523202124 (1) 用三次插值多项式计算f ( 0.7 ) 的近似值; (2) 用二次插值多项式计算f ( 0.95 ) 的近似值: (3) 用分段二次插值计算 f ( x ) )2.12.0(≤≤x 的近似值能保证有几位有 第四版 数值分析习题 第一章绪论 1.设x>0,x的相对误差为δ,求的误差. 2.设x的相对误差为2%,求的相对误差. 3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字: 4.利用公式求下列各近似值的误差限: 其中均为第3题所给的数. 5.计算球体积要使相对误差限为1%,问度量半径R时允许的相对误差限是多少 6.设按递推公式 ( n=1,2,…) 计算到.若取≈(五位有效数字),试问计算将有多大误差 7.求方程的两个根,使它至少具有四位有效数字(≈. 8.当N充分大时,怎样求 9.正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝ 10.设假定g是准确的,而对t的测量有±秒的误差,证明当t增加时S的绝对误差增加,而 相对误差却减小. 11.序列满足递推关系(n=1,2,…),若(三位有效数字),计算到时误差有多大这个计算过程 稳定吗 12.计算,取,利用下列等式计算,哪一个得到的结果最好 13.,求f(30)的值.若开平方用六位函数表,问求对数时误差有多大若改用另一等价公式 计算,求对数时误差有多大 14.试用消元法解方程组假定只用三位数计算,问结果是否可靠 15.已知三角形面积其中c为弧度,,且测量a ,b ,c的误差分别为证明面积的误差满足 第二章插值法 1.根据定义的范德蒙行列式,令 证明是n次多项式,它的根是,且 . 2.当x= 1 , -1 , 2 时, f(x)= 0 , -3 , 4 ,求f(x)的二次插值多项式. 3. 4.给出cos x,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数 字,研究用线性插值求cos x 近似值时的总误差界. 5.设,k=0,1,2,3,求. 6.设为互异节点(j=0,1,…,n),求证: i) ii) 7.设且,求证 8.在上给出的等距节点函数表,若用二次插值求的近似值,要使截断误差不超过,问使用函 数表的步长应取多少 9.若,求及. 10.如果是次多项式,记,证明的阶差分是次多项式,并且为正整数). 11.证明. 12.证明 13.证明 14.若有个不同实根,证明 15.证明阶均差有下列性质: i)若,则; ii)若,则. 16.,求及. 17.证明两点三次埃尔米特插值余项是 并由此求出分段三次埃尔米特插值的误差限. 18.求一个次数不高于4次的多项式,使它满足并由此求出分段三次埃尔米特插值的误差限. 19.试求出一个最高次数不高于4次的函数多项式,以便使它能够满足以下边界条件,,. 20.设,把分为等分,试构造一个台阶形的零次分段插值函数并证明当时,在上一致收敛到. 21.设,在上取,按等距节点求分段线性插值函数,计算各节点间中点处的与的值,并估计误 差. 22.求在上的分段线性插值函数,并估计误差. 23.求在上的分段埃尔米特插值,并估计误差. i) ii) 25.若,是三次样条函数,证明 i); ii)若,式中为插值节点,且,则. 26.编出计算三次样条函数系数及其在插值节点中点的值的程序框图(可用式的表达式). 第三章函数逼近与计算 1.(a)利用区间变换推出区间为的伯恩斯坦多项式. (b)对在上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误 差做比较. 2.求证: (a)当时,. (b)当时,. 3.在次数不超过6的多项式中,求在的最佳一致逼近多项式. 数值分析教案 土建学院 工程力学系 2014年2月 一、课程基本信息 1、课程英文名称:Numerical Analysis 2、课程类别:专业基础课程 3、课程学时:总学时32 4、学分:2 5、先修课程:《高等数学》、《线性代数》、《C 语言》 6、适用专业:工程力学 二、课程的目的与任务: 数值分析是工程力学专业的重要理论基础课程,是现代数学的一个重要分支。其主要任务是介绍进行科学计算的理论方法,即在计算机上对来自科学研究和工程实际中的数学问题进行数值计算和分析的理论和方法。通过本课程的学习,不仅使学生初步掌握数值分析的基本理论知识,而且使学生具备一定的科学计算的能力、分析问题和解决问题的能力,为学习后继课程以及将来从事科学计算、计算机应用和科学研究等工作奠定必要的数学基础。 三、课程的基本要求: 1.掌握数值分析的常用的基本的数值计算方法 2.掌握数值分析的基本理论、分析方法和原理 3.能利用计算机解决科学和工程中的某些数值计算应用问题,增强学生综合运用知识的能力 4.了解科学计算的发展方向和应用前景 四、教学内容、要求及学时分配: (一) 理论教学: 引论(2学时) 第一讲(1-2节) 1.教学内容: 数值分析(计算方法)这门课程的形成背景及主要研究内容、研究方法、主要特点;算法的有关概念及要求;误差的来源、意义、及其有关概念。数值计算中应注意的一些问题。 2.重点难点: 算法设计及其表达法;误差的基本概念。数值计算中应注意的一些问题。3.教学目标: 了解数值分析的基本概念;掌握误差的基本概念:误差、相对误差、误差限、相对误差限、有效数字;理解有效数字与误差的关系。学会选用相对较好的数值计算方法。数值分析典型习题
数值分析试题及答案汇总
数值分析典型例题
数值分析第4章答案
数值分析典型习题资料
华南理工大学数值分析试题-14年下-C
郑州大学数值分析重点考察内容及各章习题
数值分析第四章数值积分与数值微分习题答案
数值分析典型例题
数值分析典型例题
数值分析习题集及答案Word版
《数值分析》杨大地-标准答案(第八章)
数值分析第四版习题及答案
数值分析习题
数值计算第四章课后习题答案
数值分析试题1
数值分析第四版习题和答案解析
数值分析教案