数值分析典型习题
题型一:有效数字
1,
的首位数字x 1,
x *的相对误差不超过0.5×10-5,至少
要保留几位有效数字.(2010-2011)
1*115
12
11||10100.510222
6n n r x n e x n ---=≤
?=?≤??≥=解答:设至少要保留位有效数字,则有解得, n 5.7取位有效数字.
2,
0.5×10-4,至少要保留几位有效数字?(2009-2010)
3,已知21.787654为有效数,确定其绝对误差界与相对误差界.(2007-2008)
*6
*1187
11
||102111||1010102224
n r e e x ----=?=?=?=??解答:
4,已知30.49876为有效数,确定其绝对误差界.(2006-2007B) 5,设有效数x=12.4567,确定x 的绝对误差界.(2004-2005)
题型二:插值多项式
1,已知f(x)的函数值:f(0)=-2, f(1)=1, f(2)=5, 用反插值法求f(x)=0在[0,2]内的近似根x *.(2010-2011)
11111202012012010210122021()()()()()()()
()()()()()()()()()()
(2)(5)(2)(1)
012(12)(15)(52)(51)2991422884
y y y y y y y y y y y y y L y f y f y f y y y y y y y y y y y y y y y y y y ----------=?+?+?
------+-+-=+?
+?
+-+-=+-解答:对y=f(x)的反函数x=f 进行二次插值2*229
(0)42
y x L ≈=
故,
2,已知f(x)的如下函数值及导数值:f(-1)=1, f(0)=2, f ’(0)=3, f(1)=7; (1),建立不超过3次的埃尔米特插值多项式H 3(x);
(2),x ∈[-1,1], 确定用H 3(x)代替f(x)的误差界(已知|f (4)(x)|≤M 4,x ∈[-1,1]).(2010-2011)
32001001201232233)),(0,1,2)()()[,]()[,,]()()1(1)2(1)(0)
232()()(1)(0)(1)232()'(i i H x f x i N x f x f x x x x f x x x x x x x x x x x x H x N x k x x x x x k x x H ===+-+--=++++-=++=++--=+++-L 解答:(1),满足插值条件((的二次插值多项式为:也可用拉格朗日插值法满足题设插值条件的插值多项式为:2323(4)23
44
3)43(31)'(0)'(0)3()232
()(2),(1)(0)(1),(1,1)4!
1||=4!496
x x k x H f H x x x f R x x x M M R ζζ=++-===+++--∈-≤?由得:k=0故:误差(x)=则误差界(x)
3,已知f(x)的函数值:f(0)=2, f(1)=4, f(2)=9, 写出二次拉格朗日插值多项式及余项.(2009-2010)
4,已知f(x)的如下函数值及导数值:f(1)=1, f(2)=2, f ’(1)=3, f(3)=9; (1),建立不超过3次的埃尔米特插值多项式;
(2)计算f(1.6)的近似值;若M 4=0.5,估计f(1.6)的误差界.(已知|f (4)(x)|≤M 4).(2009-2010)
5,写出满足条件H(0)=1, H(1)=0, H ’(1)=1, H(2)=1的三次插值多项式,并给出误差估计式.(2008-2009B)
6,已知一组数据,求函数f(x)=0的根.(2008-2009B)
7,已知f(x)的如下函数值及导数值:f(0)=1, f(1)=3, f ’(1)=1, f(2)=9, (1),建立不超过3次的埃尔米特插值多项式,写出误差估计式;
(2),计算f(1.8)的近似值:若M 4=1,估计f(1.8)的误差界.(已知|f (4)(x)|≤M 4).(2007-2008)
8,已知f(x)的如下函数值及导数值:f(1)=2, f(2)=4, f ’(2)=5, f(3)=8, (1),建立不超过3次的埃尔米特插值多项式;
(2),计算f(2.5)的近似值:若M 4=0.5,估计f(2.5)的误差界.(已知|f (4)(x)|≤M 4).(2006-2007)
9,已知f(x)的如下函数值表
选取合适的插值节点,用二次插值多项式计算f(0.35)的近似值.(2005-2006) 10,已知f(x)=sinx 的如下函数值表
用插值多项式计算sin1.8, 并估计误差界.(2004-2005)
11,用f(x)的关于互异节点集112{}{}n n
i i i i x x -==和
的插值多项式g(x)和h(x)构造出关于节点集1{}n
i i x =的插值多项式.(2005-2006)(课后习题)
-11111121111{}(),()(){}(),()()()()))()
())]()
n n i i i i n n n n n n n n n n n n n n q x q x g x x x x x x x x x g A x x g x ==------=----=
-L 解答:法一:设关于节点集x 的插值多项式为则与有共同插值节点x ,则设:q(x)=g(x)+Aw w f(x (x )
由q(x )=f(x 得,w w 故:q(x)=g(x)+[f(x (x )w 法二:设q(x)=g(x)+1-122311111()
()(){}()()()()(),01
()=()[()()]()[()()]
()()()()()()[()()]=-n n i i n n n n n n n n n n x g x h x B g x h x B x x x x x x B x x x g x h x B
A
x x g x h x B
q x f x h x A
h x g x x x g x h x B
A B -=---=---≠----===+--L Aw 由于和有共同插值节点x ,则存在常数,使得则,w 故:q(x)=g(x)+由得得
1
111()
[()()]
()
n n x x x x h x g x x x ----则:q(x)=g(x)+
12,(1),已知f(x)的如下函数值:f(0)=1,f(1)=3,f(3)=5,写出二次拉格朗日插值多项式L 2(x);
(2),若同时已知:f ’(1)=1,用待定系数法求埃尔米特插值多项式H 3(x);
(3),当(3)(4)
1|()|2|()|4,[0,3]f x f x x ≤≤≤≤∈及3时,x 不取节点,[0,3]x ∈,求
32()()
|
|()()
f x H x f x L x --的上界.(2011-2012)
题型三:最佳平方逼近多项式及最小二乘法
1,已知函数值表:
用二次多项式y=C 0+C 1X+C 2X 2
按最小二乘法拟合改组数据,并求平方逼近误差.(2010-2011)(2005-2006)
()
000102030410111213142021222324012()()()()()11111()()()()()21012()()()()()4101401210,5010010010034T T T T x x x x x A x x x x x x x x x x y A AC A y c c c ??????????????????? ? ?
==-- ? ?
? ?
????
==??? ?? ?? ?????解答:
法一:线性拟合的法方程组为:即()()0122
2*
2
0000100011402583
,0,357583357
0581358||||=(y,y)-Y 01210402023531701(,)0,(,)(T c c y x C x x
x δ???α?α???α???? ?=? ?? ????
===-
=-???? ? ? ? ? ?=-= ? ? ? ?- ? ?????
====解得:c 则平方逼近误差:法二:构造首项系数为的正交多项式:
(x)=1
(x)=x-1112
11021100
002
*
22
022
2
2
0,)0
(,)(,)2,()()2(,)(,)46583
()()0(2)(,)514357
(,)8
||||=(y,y)-(,)35i i i i i i i i i
x x y x x x x y ?????βαα?β?????????δ??======----==++-=-=
∑
∑(x)(x)=x 则,平方逼近误差:
2,求2
1
()1f x x =
+在区间[0,1]上的一次最佳平方逼近多项式及平方逼近误差(去
权函数ρ(x)=x).(2009-2010) 3,通过实验获得以下数据:
请用最小二乘法求形如y=a+bx 2
的经验公式.(2008-2009)
T T A AC A y =解析:
4,利用正交多项式的性质构造首项系数为1的正交多项式1{()}i i g x ∞
=,有下列公
式:010
111()1
()()()()(),(1,2,...)
k k k k k g x g x x g x x g x g x k ααβ+--==-=--=其中:
111(,)
,(0,1,2...)
(,)
(,)
,(1,2...)
(,)
k k k k k k k k k k xg g k g g g g k g g αβ---=
=== (1),求[0,1]上首项系数为1的正交多项式(权函数ρ(x)=1),g 0(x),g 1(x),g 2(x) (2),以上述正交多项式为基,求sinx 在区间[0,1]上的二次最佳平方逼近多项式,并求平方逼近误差.(2008-2009B)(2004-2005)
01
00001
1
0001
2
01111211021102110000*010001(1),()1(,)11,()(,)
2
2
1()(,)121(,)2()2
(,)11,()()()()(,)126(,)(,)(2),()(,)(g x xdx xg g g x x x g g dx x x dx xg g g g x dx g g g x x g x g x x x g g g f g f x g g g g αααβαβ?===
==-=-
-===
-===--=-+
=
+????解答:212
12211
1
2
0020
111222
0002
2
2
*
2
2
0(,)
,)(,)
11()sin ()sin sin 11621()()1126()()260.00746 1.09130.23546(,)||||(,)0.000623.
(,)i i i i
g f g g g g g x x xdx x xdx xdx x x x dx x dx x x dx x x f g f f f g g ?=+-+-=?+?-+?-+--+=-+--=-=??????∑平方逼近误差:
5,以正交多项式为基,求函数2
1()1f x x =
+在区间[0,1]上的二次最佳平方逼近多
项式,并求平方逼近误差.(2007-2008)(权函数ρ(x)=x,(2011-2012))
2012012220
1201()1,(),(),111()2,()1,()22422
111122342
1111345411
1112224
5
61.0656,0.503x x x x x f In f f In C F In c c c In ???π???π=====-=-=????
? ??? ? ? ? ? ??=- ? ? ? ? ? ??? ? ?- ? ?
??
??==-解答:
法一:取解得,,,正规方程组为:H 即:解得:c c 2*222*00001000111110110002,0.07423
() 1.06560.503020.07423=(f,f)-F 0.000029041()1
1
(,)223,()1(,)332
(,)8(,)1,(,)15(,)T n p x x x C g x xg g g x x x g g xg g g g g g g g δαααβ=-=--======-=-
====
c 故二次最佳平方逼近多项式:平方逼近误差:法二:构造首项系数为的正交多项式:221100*2
01220120011222*18
82163
()()()()()()15318510
(,)(,)(,)
()()()() 1.06560.503020.07423(,)(,)(,)=(f,f)-F 0.00002904
T n g x x g x g x x x x x f g f g f g p x g x g x g x x x g g g g g g C αβδ=--=-
--=-+=++=--=则:平方逼近误差:
6,通过实验获得以下数据:
请用最小二乘法求形如v =的经验公式,并求平方误差.(2006-2007)
01
:c c v
=+解答转化
题型四:代数精确度
1,确定参数α,使求积公式20
()[(0)()]['(0)'()]2
h
h
f x dx f f h h f f h α≈++-?的代数精确度
尽可能高,并求其代数精确度.(2010-2011)
23322
4423
20()1,,()1
(),=
121()()(0)(03)2121()()0+)(04)212
()[(0)()]['(0)'()]2
h h h f x x f x f x x h f x x f x dx h h h h f x x f x dx h h h h
f x dx f f h h f f h αα====++-=≠+-≈++-???解答:令显然成立
令得又时:时:(故具有三次代数精确度.
2,确定参数A 1,A 2,使求积公式12()()(0)()3
h
h h
f x dx A f h A f f h -≈-++?的代数精确度尽可能高,并求其代数精确度.(2009-2010)
3,建立高斯型求积公式1
211221()()()x f x dx A f x A f x -≈+?.(2009-2010)
23
12
121
13
11221
1
224
11221
1
335
11221
1212
00
010
00
1
,
2
3
2
5
1
3
()1
(,)
0,()
(,)
x
A A x dx
A x A x x dx
A x A x x dx
A x A x x dx
x A A
g x
xg g
g x x x
g g
αα
α
-
-
-
-
+==
+==
+==
+==
=-===
=
===-=
?
?
?
?
解答:
法一:已知求积公式有3次代数精确度,令f(x)=1,x,x得
解上述方程组得:x
法二:构造二次正交多项式
1111
1100
2
21100
212
11
22
21
12
11
1221
12
1
(,)(,)3
0,
(,)(,)5
3
()()()()
5
()0,
11
,
33
1
()[(
3
xg g g g
g g g g
g x x g x g x x
g x x
x x x x
A x dx A x dx
x x x x
x f x dx f f
β
αβ
ρ
--
-
====
=--=-
==-=
--
=?==?=
--
≈+
??
?
令得高斯点: x
故高斯型求积公式为:
方法三:设[-1,1]上权(x)22
2
1
2
2
1
12
2
1
2
2
12
1122
22
1122
33
1122
2
12121
().
223
()0,+0,
535
2
()0,0,0
5
3
().
5
2
:
3
2
5
()()(),(
g x x ax b
b
x g x dx b
a
x xg x dx a
g x x
A A
A x A x
A x A x
A x A x
x x x x x x c x c x
??
-
-
=++
===-
?===
=-
+=
+=
+=
+=
=--=++
?
?
=x,首项系数为1的二次正交多项式为
则有:即
即
所以剩下步骤同法二.
法四
显然
2
2222
11221111222122112211122212
2
2
3322
11122211221112221122
1
1
2
)()0
()()()()()()()
223
0,
535
()()()()()
2
0,0
5
3
(),
5
x
A x A x A x c x c A x c x c A x A x c A x A x c A A
c
c
A x x A x x A x A x c A x A x c A x A x
c
c
x x
?
??
??
?
==
+=+++++=+++++
=+==-
+=+++++
===
=-剩下步骤同法二.
4,确定求积公式()()(0)()h
h f x dx Af h Bf Cf h -≈-++?中的参数A,B,C ,使其代数精度尽量高,并指出其代数精确度.(2008-2009B) 5,确定求积公式1
0211123
()()()()343234
f x dx f f f ≈
-+?的代数精确度.(2006-2007B) 6,确定下列求积公式中的参数,使求积公式的代数精确度尽可能高,并求出代数精确度1
0120113
()(
)()()424
f x dx A f A f A f ≈++?.(2005-2006) 7,确定下列求积公式中的参数,使求积公式的代数精确度尽可能高,并求出代数精确度101()()(0)()h
h f x dx A f h A f A f h --≈-++?.(2004-2005) 8,已知h>0,建立高斯型求积公式:
21122()()()h
h
x f x dx A f x A f x -≈+?
.(2011-2012)
题型五:求积公式的最少节点数
1,设定积分32
0x e dx -
?,问用复化辛普森(Simpson)求积公式进行计算,要求误差小于10-6
,所需要的最少节点数为多少?(2010-2011)
(4)2
2
4
4(4)46
1
(),()16
301[]||()|1018018016960
17.05
19.
x x
S f x e f x e b a h f h f h b a
h
η---==--=-≤?=<-=解答:复化辛普森公式截断误差:|R 解得:h<0.176,n>故应取个节点 2,设定积分13
0x e dx -
?,问用复化梯形求积公式进行计算,要求误差小于10-6
,所
需要的最少节点数为多少?(2009-2010)
(2)
3
3
2
2(2)26
1(),()9
101[]||()|1012189162
2.8
.
x x T f x e f
x e
b a h f h f h b a
h
η-
--==--=-≤?=<-=解答:复化梯形公式截断误差:|R 解得:h<0.357,n>故应取4个节点
3,给定积分2
0cos2xdx ?,问用复化梯形求积公式和复化辛普森(Simpson)求积公式进行计算,要求误差小于10-6
,所需要的最少节点数各为多少? (注:2(2)4(4)
[](),[](),[,]122880
T S b a b a R f h f R f h f a b ηηη--=-
=-∈)(2008-2009B) 4,给定积分14
0x e dx -
?,问用复化梯形求积公式和复化辛普森(Simpson)求积公式进行计算,要求误差小于10-6
,所需要的最少节点数各为多少?(2007-2008) 5,给定积分2
1Inxdx ?,问用复化梯形求积公式和复化辛普森(Simpson)求积公式进行计算,要求误差小于10-6,所需要的最少节点数各为多少? (已知:2(2)4(4)
1212[](),[](),,(,)12180
T S b a b a R f h f R f h f a b ηηηη--=-=-∈)(2006-2007) 6,用积分8
2
1
22dx In x
=?计算In2,要使所得近似值具有7位有效数字,问用复化辛普森求积公式至少需要取多少个节点?(2005-2006)
4(4)
8(4)5
2(4)-7
44(4)4-7[](),[2,8]
18011122,(),()223
|()|,[2,8]
8
1
7[]102
631
[]||()|10180180880282
0.04472,S S S b a R f h f In dx f x f x x x x
f x x R f b a h R f h f h h n h
ηηη-=-∈===≤∈≤?-=-≤?=≤?-≤≥=?解答:复化辛普森公式截断误差公式:
则使所得的近似值具有位有效数字,即令:|134.2
137故至少需要取个节点.
7,用积分6
21
3dx In x
=?计算In3,要使所得近似值具有5位有效数字,问用复化梯形求积公式至少需要取多少个节点?(2004-2005)
8,对于定积分1
0()I f x dx =?,当M 2=1/8,M 4=1/32,用11点的复化辛普森(Simpson)求积公式求I 的截断误差为R s [f],用n 个节点的复化梯形求积公式求I 的截断误差为R T [f],要使R T [f]≤R s [f],n 至少是多少?
(M 2=max|f ”(x)|,M 4=max|f (4)
(x)|,[0,1]x ∈).(2011-2012)
题型六:Doolittle 分解及方程组求解
1,求矩阵212454635??
? ?
?-??
的Doolittle 分解.(2010-2011)
212100212454210030635321001LU ?????? ? ???
== ? ??? ? ???---??????
解答:A=
2,求矩阵114103241??
?- ?
???
的Doolittle 分解.(2009-2010)
3,设线性方程组1234101351
1415241016
2
1
16x x x x ---?????? ? ? ?- ?
? ?
?= ? ? ?- ? ? ?-????
?? (1),对方程组的系数矩阵A 作Doolittle 分解;
(2),用所得的Doolittle 分解求该线性方程组的解.(2007-2008&2005-2006)
123412341
010
1
3101
311000
1321
14
124100
0132241011
1916
2116
210
001313191,,,)(5,0,11,)
13
,,,)(1,1,1,1).T T
T T A LU LY b y y y UX Y x x x --????--?? ??
?
- ? ???-
?
=== ???
--- ? ??
? ? ???
-- ???
???
???==---==--解答:由得:(y 由得:(x
4,设线性方程组123411
4
151
0131241076
2118x x x x -?????? ? ? ?
-- ?
? ?
?=
? ? ? ? ? ?
-????
?? (1),对方程组的系数矩阵A 作Doolittle 分解;
(2),用所得的Doolittle 分解求该线性方程组的解.(2006-2007)
5,设线性方程组:123123123231
53478113
x x x x x x x x x ++=+-=-++=-
(1),对方程组的系数矩阵A 作Doolittle 分解; (2),利用上述分解结果求解该线性方程组.(2004-2005)
6,用高斯顺序消去法求解线性方程组:
1324123424253
2431737
x x x x x x x x x x +=+=+++=+=.(2010-2011)
4321102051
02051
02050
10130
10130
1013=124317
022312
002160
10
370
10
370
00242,2,1, 1.x x x ?????? ? ? ?
? ? ?→→ ? ? ? ? ? ???????
====解答:增广矩阵回代求解:x
7,用高斯顺序消去法求解线性方程组:123123123347
2212320
x x x x x x x x x -+=-+-=---=.(2009-2010)
题型七:条件数及范数
1,求线性方程组121239107
898
1510
x x x x x --=+==的系数矩阵A 的条件数cond 1(A),并说明其含义.
(2010-2011)
1111191008900015910089010015()||||||||1919361
1A A cond A A A A b ----?? ?
= ?
???
?? ?-- ?= ?
?
?
?
?==?=解答:系数矩阵条件数远大于,这说明当和有小扰动时会引起解的较大误差,即该方程组是病态的.
2,设矩阵150
00910089A ?? ?=-- ?
???
,求cond ∞(A).(2009-2010)
3,设三阶对称矩阵A 的特征值分别为:-2,1,3,求||A||2及cond 2(A).(2007-2008)
2121222||||3
||||()|||||||| 3.
A A cond A A A --========解答:则:
4,若n 元线性方程组Ax=b 为病态的,可以得到关于系数矩阵A 的什么性质.(2006-2007)
5,若111123124A ??
?= ?
???
,求cond 1(A).(2005-2006)求cond ∞(A).(2004-2005)
6,设1231032475A -??
?=-- ?
?-??
,求1||||||||A A ∞与.(2007-2008)
7,若1234A ??
= ???
,求谱半径()A ρ
.(2005-2006) 52
ρ解答:最大特征值:(A)=
题型八:雅可比迭代与高斯-赛德尔迭代
1,写出求解方程组12312312373212
41021534818
x x x x x x x x x -+=--=--=的雅可比迭代公式,并说明其收敛
性.(2010-2011)
(1)()()123(1)()()213(1)()()312(0)1
(3212)
71
(4215)
101(3418)8
7324102348.
k k k k k k k k k J x x x x x x x x +++=-+=--++=--++-?? ?
=-- ? ?--??
解答:雅可比迭代公式为:
x 雅可比迭代法迭代矩阵:B 严格对角占优,
故求解该方程组的雅可比迭代法关于任意初始向量x 收敛 2,设有方程组:132********
211
2212
x x x x x x x -=+=-++=,讨论用雅可比迭代法和高斯-赛德尔迭代法
解此方程组的收敛性.(2010-2011)
112330200030000202100002000121221000200020031
()0021102||0,=0=-J J J L D U B D L U E B B λλλλρ---???????? ? ? ? ?
=++=++ ? ? ? ?
? ? ? ?--????????
?
? ? ?
?
=-+=- ? ?
?- ???-=解答:A=雅可比迭代矩阵:得,()<1,故用雅可比迭代法解答此方程组对任意(0)1123(0)20031
-()00211001211
||0,=012
-S S S B D L U E B B λλλλρ-?
? ? ?
?
=-+=- ? ?
? ?
??-===
初始向量x 都收敛.
高斯赛德尔迭代矩阵:得,()<1,故用高斯赛格尔迭代法解答此方程组对任意初始向量x 都收敛.
3,写出求解方程组:12312312353212
4721535818
x x x x x x x x x -+=--=--=的高斯-赛德尔迭代公式,并说明收敛
性.(2009-2010)
4,用雅可比迭代法求解以313132323A ??
?
= ? ?-??
为系数矩阵的线性方程组时,确定其收
敛性.(2009-2010)
5,设线性方程组1231231232211
6
2222
x x x x x x x x x -+=-+-=--+=-,讨论分别用雅可比迭代法和高斯-赛德尔迭
代法解此线性方程组的收敛性,若收敛,请给出迭代格式.(2008-2009B)
6,设线性方程组:1231231232215
20
2225
x x x x x x x x x +-=-++=++=-
(1),证明求解该方程组的雅可比迭代法关于任意初始向量收敛;相应的高斯-赛德尔迭代法不是关于任意初始向量收敛; (2),取(0)
(0,0,0)T x
=,用雅可比迭代法进行求解,要求
(1)()5||||10k k x x +--<.(2007-2008)
11231123022()101220||0,===0)1022()023002||0,0,2,)1
-J J J S S S D L U E B D L U E λλλλρλλλλρ---?? ?
=-+=-- ?
?--??
-=<-??
?
=-+=- ?
???
-====>解答:
(1):B B 解得:,(B B 解得:(B 所以用雅可比迭代法解此方程组对任意初始向量都收敛,而用高斯赛德尔迭代法解此方程组不是对任意初始向量都收敛.(2):(1)()()123(1)
()
()
2
1
3
(1)()()312(0)(1)(2)(3)(4)221520
2225(0,0,0)(15,20,25)(105,60,35)(205,160,65)(205,160,65)k k k k k k k k k T T T T
T
x x x x
x x x x x x x x x +++=-+-=--+=---==--=--=-=-雅可比迭代公式:x 当时,计算得:(精确解).
7,设线性方程组:1231231238210
27325431111
x x x x x x x x x ++=--++=-+=-
(1),写出求解该方程组的雅可比迭代法的迭代公式和高斯-赛德尔迭代法的迭代公式,并确定其收敛性; (2),取(0)
(0,0,0)T x
=,用高斯-赛德尔迭代法计算x (3).(2006-2007)
8,设线性方程组Ax=b 的系数矩阵232131t A t t ?? ?
= ? ?-??
,其中t<0,问t 取何值时雅可
比迭代法关于任意初始向量都收敛.(2006-2007)
12122223021
()0310422
||()0=0=-,=
)12||<1,t<-2,or t>2
0, 2.
J J J t t D L U t t t t E B t t t
t
t t λλλλλλρ-?
?-- ? ? ?
=-+=-- ? ? ?- ?
??
-=-=<<<-解答:雅可比迭代矩阵B 得,,雅可比迭代法对于任意初始向量都收敛,则(B 即:得又故
9,1),设线性方程组:12123234324
3430424
x x x x x x x +=+-=-+=-
写出求解该方程组的雅可比迭代法的迭代公式,并确定该迭代法的收敛性;
2),设线性方程组:123123123104413
410811481025
x x x x x x x x x ++=++=++=
写出求解该方程组的高斯-赛德尔迭代法的迭代公式,并确定该迭代法的收敛性.(2004-2005)
10,给定方程组:123123123225
1
223
x x x x x x x x x +-=++=++=
(1),用三角分解法解此方程组;
(2),写出解此方程组的雅可比迭代公式,说明收敛性;取初始向量x 0=(0,0,0)T
,
当2
1||||10k k x x -+-<时,求其解.(2011-2012)
11,
设()
2
125
3sin 3
421sin
cos 43tan 5k k k k k k k A
k k k k
k
??
- ?+ ? ?
= ?+???
?
,求()lim k k A →∞.(2007-2008)
()
020lim 021205K k A →∞
??
?= ? ???
解答: 12,若()
()11,lim 1sin sin k k k k k k A
A k k k k →∞??
?
+=
? ?
???
求.(2004-2005) ()01lim 10K k A →∞
??
= ???
解答: 题型九:非线性迭代
1,
设计一个算法求.(2008-2009B)
101125
(),0.2k k k
x x x +=+>解答:牛顿迭代公式:x
2,
给出用牛顿法求围.(2010-2011)
66155
6'5"4"*00001050
517001701170
[5]66()170,()60,()300()()0,.1170
(5)61170
()(5)6k k k k k k
x x x x x x x x f x x f x x f x x x f x f x x x x x x g x x x +=-=-=-=+=-=>=>>?><=+=+
解答:的正根.
由牛顿迭代法得迭代公式:当故此时收敛到当0<
设'611*01850
()(5)0,()0,
6:0,.0.
x g x x g x g x
x x x x x ∈=-<∈>=
>>∈>故故回到前段.所以当迭代公式也收敛到综上:
3,给出用牛顿法求近似值的迭代公式,并给出初值的取值范围.(2009-2010) 解答:方法同上.
4,设φ(x)=x+c(x 2
-5),当c 为何值时,x k+1
=φ(x k ),(k=0,1,2…)产生的序列
{x k }
收敛于c 为何值时收敛最快?
(2010-2011)
2'
'**1**'*5),||<1,||<110,0;.
k k cx x c ?????+-=-<<<<解答:(x)=x+c(x (x)=1+2cx
x (x )收敛,则有(x )即1+2cx 又当(x )=0,即
5,设2
()(3)x x c x ?=+-,应如何选取常数c 才能使迭代1(),(0,1,2)k k x x
k ?+==L
具有局部收敛性?C 取何值时,这个迭代收敛最快?取x 0=2,c =计算()
x ?的不动点,要求当61||10k k x x -+-<时结束迭代.(2004-2005)
数值分析典型习题
特别声明:考试时需带计 算器作辅助计算 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.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.ρ华南理工大学数值分析试题-14年下-C
华南理工大学研究生课程考试 《数值分析》试卷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)当同时收敛时,试比较它们的收敛速度。
数值分析习题集及答案[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)的值.若开平方用六位函数表,问求对数时误差有多大?若
数值分析典型例题
第一章典型例题 例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.y' = y , x [0,1] ,y (0) =1 , h = 0.1。 1求解析解。 2 Eular法 3 R-K法 ○1解析法 在MATLAB命令窗口执行 clear >> x=0:0.1:1; >> y=exp(x); >> c=[y]' c = 1.000000000000000 1.105170918075648 1.221402758160170 1.349858807576003 1.491824697641270 1.648721270700128 1.822118800390509 2.013752707470477 2.225540928492468 2.459603111156950 2.718281828459046 ○2Euler法 在Matlab中建立M文件如下: function [x,y]=euler1(dyfun,xspan,y0,h) x=xspan(1):h:xspan(2);y(1)=y0; for n=1:length(x)-1 y(n+1)=y(n)+h*feval(dyfun,x(n),y(n)); end x=x';y=y' 在MATLAB命令窗口执行
clear >> dyfun=inline('y+0*x'); >> [x,y]=euler1(dyfun,[0,1],1,0.1); >> [x,y] 得到 ans = 0 1.000000000000000 0.100000000000000 1.100000000000000 0.200000000000000 1.210000000000000 0.300000000000000 1.331000000000000 0.400000000000000 1.464100000000000 0.500000000000000 1.610510000000000 0.600000000000000 1.771561000000000 0.700000000000000 1.948717100000000 0.800000000000000 2.143588810000000 0.900000000000000 2.357947691000000 1.000000000000000 2.593742460100000 ○3R-K法(龙格-库塔法) 在本题求解中,采用经典4阶龙格-库塔法 首先在Matlab的M文件窗口对4阶龙格-库塔算法进行编程: function [x,y]=RungKutta41(dyfun,x0,y0,h,N) x=zeros(1,N+1);y=zeros(1,N+1);x(1)=x0;y(1)=y0; for n=1:N x(n+1)=x(n)+h; k1=h*feval(dyfun,x(n),y(n)); k2=h*feval(dyfun,x(n)+h/2,y(n)+1/2*k1); k3=h*feval(dyfun,x(n)+h/2,y(n)+1/2*k2); k4=h*feval(dyfun,x(n+1)+h,y(n)+k3); y(n+1)=y(n)+(k1+2*k2+2*k3+k4)/6; end 在MATLAB命令窗口执行 clear >> dyfun=inline('y','x','y'); >> [x,y]=RungKutta41(dyfun,0,1,0.1,10); >> c=[x;y]' 得到
数值分析典型例题
数值分析典型例题 例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)的值.若开平方用六位函数表,问求对数时误差有多大? ¥ 数值分析思考题1 1、讨论绝对误差(限)、相对误差(限)与有效数字之间的关系。 2、相对误差在什么情况下可以用下式代替 3、查阅何谓问题的“病态性”,并区分与“数值稳定性”的不同点。 4、取 ,计算 ,下列方法中哪种最好为什么(1)(3 3-,(2)(2 7-,(3) ()3 1 3+ ,(4) ()6 1 1 ,(5)99- , 数值实验 数值实验综述:线性代数方程组的解法是一切科学计算的基础与核心问题。求解方法大致可分为直接法和迭代法两大类。直接法——指在没有舍入误差的情况下经过有限次运算可求得方程组的精确解的方法,因此也称为精确法。当系数矩阵是方的、稠密的、无任何特殊结构的中小规模线性方程组时,Gauss消去法是目前最基本和常用的方法。如若系数矩阵具有某种特殊形式,则为了尽可能地减少计算量与存储量,需采用其他专门的方法来求解。 Gauss消去等同于矩阵的三角分解,但它存在潜在的不稳定性,故需要选主元素。对正定对称矩阵,采用平方根方法无需选主元。方程组的性态与方程组的条件数有关,对于病态的方程组必须采用特殊的方法进行求解。 数值计算方法上机题目1 1、实验1. 病态问题 实验目的: 算法有“优”与“劣”之分,问题也有“好”和“坏”之别。所谓坏问题就是问题本身的解对数据变化的比较敏感,反之属于好问题。希望读者通过本实验对此有一个初步的体会。 数值分析的大部分研究课题中,如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决,当然一般要付出一些代价(如耗用更多的机器时间、占用更多的存储空间等)。 $ r e x x e x x ** * ** - == 141 . ≈)61 习题一 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-=ε)。 数值分析试卷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 的近似值能保证有几位有 数值分析教案 土建学院 工程力学系 2014年2月 一、课程基本信息 1、课程英文名称:Numerical Analysis 2、课程类别:专业基础课程 3、课程学时:总学时32 4、学分:2 5、先修课程:《高等数学》、《线性代数》、《C 语言》 6、适用专业:工程力学 二、课程的目的与任务: 数值分析是工程力学专业的重要理论基础课程,是现代数学的一个重要分支。其主要任务是介绍进行科学计算的理论方法,即在计算机上对来自科学研究和工程实际中的数学问题进行数值计算和分析的理论和方法。通过本课程的学习,不仅使学生初步掌握数值分析的基本理论知识,而且使学生具备一定的科学计算的能力、分析问题和解决问题的能力,为学习后继课程以及将来从事科学计算、计算机应用和科学研究等工作奠定必要的数学基础。 三、课程的基本要求: 1.掌握数值分析的常用的基本的数值计算方法 2.掌握数值分析的基本理论、分析方法和原理 3.能利用计算机解决科学和工程中的某些数值计算应用问题,增强学生综合运用知识的能力 4.了解科学计算的发展方向和应用前景 四、教学内容、要求及学时分配: (一) 理论教学: 引论(2学时) 第一讲(1-2节) 1.教学内容: 数值分析(计算方法)这门课程的形成背景及主要研究内容、研究方法、主要特点;算法的有关概念及要求;误差的来源、意义、及其有关概念。数值计算中应注意的一些问题。 2.重点难点: 算法设计及其表达法;误差的基本概念。数值计算中应注意的一些问题。3.教学目标: 了解数值分析的基本概念;掌握误差的基本概念:误差、相对误差、误差限、相对误差限、有效数字;理解有效数字与误差的关系。学会选用相对较好的数值计算方法。 1. 正方形的边长大约为100cm ,应怎样测量才能使面积误差不超过1cm 2? 2. 已测得某场地长l 的值为110=*l m ,宽d 的值为80=*d m ,已知 2.0≤-*l l m, 1.0≤-*d d m, 试求面积ld s =的绝对误差限与相对误差限. 3.为使π的相对误差小于0.001%,至少应取几位有效数字? 4.设x的相对误差界为δ,求n x的相对误差界. 5.设有3个近似数a=2.31,b=1.93,c=2.24,它们都有3位有效数字,试计算 p=a+bc的误差界和相对误差界,并问p的计算结果能有几位有效数字? 6. 已知33348 7.034.0sin ,314567.032.0sin ==,请用线性插值计算3367.0sin 的值,并估计截断误差. 7. 已知sin0.32=0.314567, sin0.34=0.333487, sin0.36= 0.352274,用抛物插值计算sin0.3367的值, 并估计误差. 8. 已知 1 6243sin ,sin π ππ== =请用抛物插值求sin50的值,并估计误差 9. . .6,8,7,4,1)(,5,4,3,2,1求四次牛顿插值多项式时设当==i i x f x 10. 已知4)2(,3)1(,0)1(=-=-=f f f , 求函数)(x f 过这3点的2次牛顿插 值多项式 . 11. 设x x f =)(,并已知483240.1)2.2(,449138.1)1.2(,414214.1)0.2(===f f f , 试用二次牛顿插值多项式计算(2.15)f 的近似值,并讨论其误差 12. 设],[)(b a x f 在上有四阶连续导数,试求满足条件)2,1,0()()(==i x f x P i i 及 )()(11x f x P '='的插值多项式及其余项表达式. 13. 给定3201219(),,1,,44f x x x x x ====试求()f x 在1944?? ???? ,上的三次埃尔米特 《计算方法》期中复习试题 一、填空题: 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ,则用辛普生(辛卜生)公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案:2.367,0.25 2、1)3(,2)2(,1)1(==-=f f f ,则过这三点的二次插值多项式中2 x 的系数为 , 拉格朗日插值多项式为 。 答案:-1, )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字; 4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( ); 答案 )(1)(1n n n n n x f x f x x x '--- =+ 5、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 ); 6、计算方法主要研究( 截断 )误差和( 舍入 )误差; 7、用二分法求非线性方程f (x )=0在区间(a ,b )内的根时,二分n 次后的误差限为 ( 1 2+-n a b ); 8、已知f (1)=2,f (2)=3,f (4)=5.9,则二次Newton 插值多项式中x 2系数为 ( 0.15 ); 11、 两点式高斯型求积公式?1 0d )(x x f ≈(?++-≈1 )] 321 3()3213([21d )(f f x x f ),代数精 度为( 5 ); 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少,应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ,为了减少舍入误差,应将表达式 第四版 数值分析习题 第一章绪论 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的多项式中,求在的最佳一致逼近多项式. 数值分析典型例题 Revised as of 23 November 2020 第一章典型例题 例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 于是有同解方程组 ?? ? ??-==--=++17175.555.01 42332321x 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.确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度: (1)()()()()? --++-≈h h h f A f A h f A dx x f 1010; (2)()()()()? --++-≈h h h f A f A h f A dx x f 221010; (3)()()()()[]3/3211 121?-++-≈x f x f f dx x f ; (4)()()()[]()()[]h f f ah h f f h dx x f h '0'2/020 +++≈? 解:(1)求积公式中含有三个待定参数,即101A A A ,,-,将()21x x x f ,,=分别代入求积公式,并令其左右相等,得 ()()??? ???? =+=+-=++---3 1121 110132 02h A A h A A h h A A A 解得h A h A A 34 31011===-,。 所求公式至少具有2次代数精度。又由于 ()() ()() 4 4 4 3 33 3 3 33h h h h dx x h h h h dx x h h h h ? ?--+ -≠ +-≈ 故()()()()? --++-≈h h h f A f A h f A dx x f 1010具有三次代数精度。 (2)求积公式中含有三个待定系数:101A A A ,,-,故令公式对()2 1x x x f ,,=准确成立,得()()??? ???? =+=+-=++---3 1121110131604h A A h A A h h A A A ,解得h h h A h A h A A 34 316424381011-=- =-===-, 故()()()[]()03 43 822hf h f h f h dx x f h h - +-≈ ? - 因()?-=h h dx x f 220 而 ()() []03 83 3 =+-h h h 又[ ]4 45 5 6224 3 83 165 2h h h h h dx x h h += ≠= ? -数值分析习题集及答案Word版
数值分析作业思考题汇总
数值分析习题
数值分析试题1
数值分析教案
数值分析题库答案
《数值计算方法》精彩试题集及问题详解1-6 2
数值分析第四版习题和答案解析
数值分析典型例题
数值分析习题第四章